Salt Bridge Research Articles

Molecular Docking Studies and In Vitro Activity of Pancreatic Lipase Inhibitors from Yak Milk Cheese

Pancreatic lipase serves as a primary trigger for hyperlipidemia and is also a crucial target in the inhibition of hypercholesterolemia. By synthesizing anti-hypercholesterolemic drugs such as atorvastatin, which are used to treat hypercholesterolemia, there were some side effects associated with the long-term use of statins. Based on this idea, in the present study, we identified peptides that inhibited PL by virtual screening and in vitro activity assays. In addition, to delve into the underlying mechanisms, we undertook a dual investigative approach involving both molecular docking analyses and molecular dynamics simulations. The results showed that peptides RK7, KQ7, and TL9, all with molecular weights of <1000 Da and a high proportion of hydrophobic amino acids, inhibited PL well. Molecular docking and molecular dynamics showed that peptides RK7, KQ7, and TL9 bound to important amino acid residues of PL, such as Pro and Leu, through hydrogen bonding, hydrophobic interactions, salt bridges, and π-π stacking to occupy the substrate-binding site, which inhibited PL and identified them as potential PL inhibitors. In vitro tests showed that the IC50 of RK7 and KQ7 on PL were 0.690 mg/mL and 0.593 mg/mL, respectively, and the inhibitory effects of RK7 and KQ7 on PL were significantly enhanced after simulated gastrointestinal digestion. Our results suggested that peptides RK7 and KQ7 from yak milk cheese can be identified as a novel class of potential PL inhibitors.

Unbinding of Alpha Chain of Hemoglobin in Sickle and Normal Structures

Abstract Sickle cell disease, a genetic disorder, is caused by a mutation of glutamic acid into valine in the β chain of hemoglobin at the sixth residue, resulting in structural change of the entire hemoglobin molecule into a sickle shape. We investigated the atomic level interaction between the α chain (chain A) and the remaining three chains to identify the structural modification in sickle hemoglobin using the molecular dynamics simulations. H-bond, solvent accessible surface area (SASA), hydrophobic interactions, salt bridges of sickle and normal hemoglobin have been estimated. The estimated parameters from sickle hemoglobin are compared to normal hemoglobin structure. Steered molecular dynamics (SMD) is utilized to estimate the force required in breaking hydrogen bonds in given chains. The SMD simulations at different pulling velocities show that the decoupling forces depend on values of pulling force. This relation is linear, 6780 pN to 12345 pN with pulling velocities of 0.00020 nm/ps to 0.00040 nm/ps in sickle hemoglobin. Much higher force of 8738 pN to 16557 pN in normal is required in normal hemoglobin with same spring constants values from k = 500 to 1100 kcal mol-1 nm-2 and same pulling velocities.

Dissecting the Binding Affinity of Anti-SARS-CoV-2 Compounds to Human Transmembrane Protease, Serine 2: A Computational Study

The human transmembrane protease, serine 2 (TMPRSS2), essential for SARS-CoV-2 entry, is a key antiviral target. Here, we computationally profiled the TMPRSS2-binding affinities of 15 antiviral compounds. Molecular dynamics (MD) simulations for the docked complexes revealed that three compounds exited the substrate-binding cavity (SBC), suggesting noncompetitive inhibition. Of the remaining compounds, five charged ones exhibited reduced binding stability due to competing electrostatic interactions and increased solvent exposure, while seven neutral compounds showed stronger binding affinity driven by van der Waals (vdW) interactions compensating for unfavorable electrostatic effects (including electrostatic interactions and desolvation penalties). Positive and negative hotspot residues were identified as uncharged and charged, respectively, both lining the SBC. Despite forming diverse interactions with compounds, the burial of positive hotspots led to strong vdW interactions that overcompensated for unfavorable electrostatic effects, whereas negative hotspots incurred high desolvation penalties, negating any favorable contributions. Charged residues at the SBC’s outer rim can reduce binding affinity significantly when forming hydrogen bonds or salt bridges. These findings underscore the importance of enhancing vdW interactions with uncharged residues and minimizing the unfavorable electrostatic effects of charged residues, providing valuable insights for designing effective TMPRSS2 inhibitors.

Design, synthesis and cytotoxic activity of novel lipophilic cationic derivatives of diosgenin and sarsasapogenin.

Novel lipophilic cationic derivatives including quaternary ammonium salt and triphenylphosphine series were designed and synthesized using diosgenin (1) and sarsasapogenin (2) as substrates to improve the cytotoxicity and selectivity. Most of the derivatives showed higher cytotoxicity against all cancer cell lines tested, compound 13 exhibited the most superior activity against A549 cells with an IC50 value of 0.95μM, which was 34-fold of diosgenin. Preliminary cellular mechanism studies elucidated that compound 13 might arrest cell cycle at G0/G1 phase, trigger apoptosis via up-regulating the expression of Bax, down-regulating the expression of Bcl-2 and caspase-3, and induce an increase in the generation of intracellular reactive oxygen species (ROS) in A549 cells. In addition, molecular docking analysis revealed that compound 13 could occupy the active site of p38α-MAPK well and interact to the surrounding amino acids by salt bridge and conjugation. These results suggested that compound 13 had the potential to serve as an antitumor lead agent, probably exert antitumor effect through mitochondrial pathway and p38α MAPK pathway.

Studies on inhibition of α-glucosidase using debittered formulation of Bacopa monnieri juice: Enzyme inhibition kinetics, interaction strategy, and molecular docking approach

Bacopa monnieri juice (BMJ) is traditionally used, reported, and scientifically validated for memory enhancement. However, its efficacy against diabetes is less explored. The extreme bitterness of BMJ restricts its commercial applications. This study investigates the reduction of bitterness of BMJ followed by evaluation for its α-glucosidase inhibitory activity. Initially, debittering of 30 % (v/v) BMJ using ZnSO4 (15 mM) was optimized by time-intensity analysis and molecular docking of ZnSO4 as well as bacoside A3, the main active compound in BMJ, with TAS2R14 taste receptor. The study indicated 5 hydrogen bonds to be involved in binding with bacoside A3 with binding energy of −11.82 Kcal/mol, while hydrogen bond, salt bridges and metal complexes were involved in binding of ZnSO4 with binding energy of −6.65 Kcal/mol. Subsequently, BMJ, ZnSO4 and BMJ + ZnSO4 (debittered juice) were also found to be potent inhibitors of α-glucosidase in dose-dependent manner. These inhibitors showed parabolic mixed inhibition of α-glucosidase, altered the secondary structure, and quenching of fluorescence. In silico studies revealed hydrogen bonding and hydrophobic interactions between inhibitors and α-glucosidase with lowest binding energy of −15.53 and −7.54 Kcal/mol being recorded for bacoside A3 and ZnSO4, respectively. Molecular docking of other bioactive compounds in BMJ such as apigenin, luteolin, quercetin and bacopasaponin C also showed lower binding energy than the standard drug, acarbose (−5.84). This study inferred the binding of bacoside A3 at the active site of α-glucosidase and of ZnSO4 with other sites on the protein. The study proposes a debittered BMJ formulation to control hyperglycemia.

Thermal performance of DispersedInorganic magnetic hybrid nanomaterials into mixed convective flow through flexible porous disks

AbstractThis study applies advanced AI techniques, including machine learning algorithms, to explore the numerically unsteady laminar flow of viscous and incompressible fluids in coaxially swirled porous disks, with applications in engineering sciences. Our focus encompasses the effects of magnetic hybrid nanomaterials and the dynamic behaviors associated with expanding/contracting and injection/suction. Utilizing single‐phase simulations, we address nonlinear coupled ordinary differential equations set against appropriate boundary conditions. Key parameters of our study include permeability and relaxation, as well as the influence of chemical reactions and mixed convection on fluid behaviors. The thermophysical properties of Al2O3/Cu nanoparticles have been by varying their morphological aspects. For the hybrid nanofluids flow, the aggregation of nanoparticle volume fraction has been designed critically in conjunction with an energy and mass transfer equation. Because dimensionless ordinary differential equations are employed, the obtained expression is transmuted using the obliging transformation technique. The desired nonlinear system of ODEs is implemented using an accurate numerical method. Our findings reveal significant impacts of chemical reaction parameters on the Sherwood number and a marked increase in the skin friction coefficient, Nusselt number, and Sherwood number as nanoparticle volume fraction rises from 2% to 7%.

The structural organisation of pentraxin-3 and its interactions with heavy chains of inter-α-inhibitor regulate crosslinking of the hyaluronan matrix.

Pentraxin-3 (PTX3) is an octameric protein, comprised of eight identical protomers, that has diverse functions in reproductive biology, innate immunity and cancer. PTX3 interacts with the large polysaccharide hyaluronan (HA) to which heavy chains (HCs) of the inter-α-inhibitor (IαI) family of proteoglycans are covalently attached, playing a key role in the (non-covalent) crosslinking of HC•HA complexes. These interactions stabilise the cumulus matrix, essential for ovulation and fertilisation in mammals, and are also implicated in the formation of pathogenic matrices in the context of viral lung infections. To better understand the physiological and pathological roles of PTX3 we have analysed how its quaternary structure underpins HA crosslinking via its interactions with HCs. A combination of X-ray crystallography, cryo-electron microscopy (cryo-EM) and AlphaFold predictive modelling revealed that the C-terminal pentraxin domains of the PTX3 octamer are arranged in a central cube, with two long extensions on either side, each formed from four protomers assembled into tetrameric coiled-coil regions, essentially as described by (Noone et al., 2022; doi:10.1073/pnas.2208144119). From crystallography and cryo-EM data, we identified a network of inter-protomer salt bridges that facilitate the assembly of the octamer. Small angle X-ray scattering (SAXS) validated our model for the octameric protein, including the analysis of two PTX3 constructs: a tetrameric 'Half-PTX3' and a construct missing the 24 N-terminal residues (Δ1-24-PTX3). SAXS determined a length of ∼520 Å for PTX3 and, combined with 3D variability analysis of cryo-EM data, defined the flexibility of the N-terminal extensions. Biophysical analyses revealed that the prototypical heavy chain HC1 does not interact with PTX3 at pH 7.4, consistent with our previous studies showing that, at this pH, PTX3 only associates with HC•HA complexes if they are formed in its presence. However, PTX3 binds to HC1 at acidic pH, and can also be incorporated into pre-formed HC•HA complexes under these conditions. This provides a novel mechanism for the regulation of PTX3-mediated HA crosslinking (e.g., during inflammation), likely mediated by a pH-dependent conformational change in HC1. The PTX3 octamer was found to associate simultaneously with up to eight HC1 molecules and, thus, has the potential to form a major crosslinking node within HC•HA matrices, i.e., where the physical and biochemical properties of resulting matrices could be tuned by the HC/PTX3 composition.

Influence of Viscosity and Thermal Conductivity in Boger Nanofluid Flow through Porous Disk: Finite Difference Analysis

Influence of Viscosity and Thermal Conductivity in Boger Nanofluid Flow through Porous Disk: Finite Difference Analysis

Functional Roles of the Charged Residues of the C- and M-Gates in the Yeast Mitochondrial NAD+ Transporter Ndt1p.

Mitochondrial carriers transport organic acids, amino acids, nucleotides and cofactors across the mitochondrial inner membrane. These transporters consist of a three-fold symmetric bundle of six transmembrane α-helices that encircle a pore with a central substrate binding site, whose alternating access is controlled by a cytoplasmic and a matrix gate (C- and M-gates). The C- and M-gates close by forming two different salt-bridge networks involving the conserved motifs [YF][DE]XX[KR] on the even-numbered and PX[DE]XX[KR] on the odd-numbered transmembrane α-helices, respectively. We have investigated the effects on transport of mutating the C-gate charged residues of the yeast NAD+ transporter Ndt1p and performed molecular docking with NAD+ and other substrates into structural models of Ndt1p. Double-cysteine substitutions and swapping the positions of the C-gate charged-pair residues showed that all of them contribute to the high transport rate of wild-type Ndt1p, although no single salt bridge is essential for activity. The in silico docking results strongly suggest that both the C-gate motif mutations and our previously reported M-gate mutations affect gate closing, whereas those of the M-gate also affect substrate binding, which is further supported by molecular dynamics. In particular, NAD+ most likely interferes with the cation-π interaction between R303-W198, which has been proposed to exist in the Ndt1p M-gate in the place of one of the salt bridges. These findings contribute to understanding the roles of the charged C- and M-gate residues in the transport mechanism of Ndt1p.

New Results of Differential Subordination for a Specific Subclass of p-Valent Meromorphic Functions Involving a New Operator

The present article aims to significantly improve geometric function theory by making an important contribution to p-valent meromorphic and analytic functions. It focuses on subordination, which describes the relationships of analytic functions. In order to achieve this, we utilize a technique that is based on the properties of differential subordination. This approach, which is one of the most recent developments in this field, may obtain a number of conclusions about differential subordination for p-valent meromorphic functions described by the new operator IHp,q,s j,pν1,n,α,lJ(ζ) within the porous unit disk Δ. Numerous mathematical and practical issues involving orthogonal polynomials, such as system identification, signal processing, fluid dynamics, antenna technology, and approximation theory, can benefit from the results presented in this article. The knowledge and comprehension of the unit’s analytical functions and its interacting higher relations are also greatly expanded by this text.

Novel Tripeptides as Tyrosinase Inhibitors: In Silico and In Vitro Approaches.

Tyrosinase is a key enzyme responsible for the formation of melanin (a natural skin pigment with ultraviolet-protection properties). However, some people experience melanin overproduction, so new, safe, and biocompatible enzyme inhibitors are sought. New tripeptide tyrosinase inhibitors were developed using molecular modeling. A combinatorial library of tripeptides was prepared and docked to the mushroom tyrosinase crystal structure and investigated with molecular dynamics. Based on the results of calculations and expert knowledge, the three potentially most active peptides (CSF, CSN, CVL) were selected. Their in vitro properties were examined, and they achieved half-maximal inhibitory concentration (IC50) values of 136.04, 177.74, and 261.79 µM, respectively. These compounds attach to the binding pocket of tyrosinase mainly through hydrogen bonds and salt bridges. Molecular dynamics simulations demonstrated the stability of the peptid-tyrosinase complexes and highlighted the persistence of key interactions throughout the simulation period. The ability of these peptides to complex copper ions was also confirmed. The CSF peptide showed the highest chelating activity with copper. The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay confirmed that none of the test tripeptides showed cytotoxicity toward the reconstructed human epidermis. Our results indicated that the developed tripeptides were non-toxic and effective tyrosinase inhibitors. They could be applied as raw materials in skin-brightening or anti-aging cosmetic products.

Exploring Binding Sites in Chagas Disease Protein TcP21 Using Integrated Mixed Solvent Molecular Dynamics Approaches.

Chagas disease, caused by the protozoan Trypanosoma cruzi, remains a significant global health burden, particularly in Latin America, where millions are at risk. This disease predominantly affects socioeconomically vulnerable populations, aggravating economic inequality, marginalization, and low political visibility. Despite extensive research, effective treatments are still lacking, partly due to the complex biology of the parasite and its infection mechanisms. This study focuses on TcP21, a novel 21 kDa protein secreted by extracellular amastigotes, which has been implicated in T. cruzi infection via an alternative infective pathway. Although the potential of TcP21 for understanding Chagas disease is promising, further exploration is necessary, particularly in identifying potential binding sites on its surface. Computational tools offer a versatile and effective strategy for preliminary binding site assessment, facilitating a more cost-efficient allocation of experimental resources. In this study, we employed three independent computational approaches─mixed solvent molecular dynamics simulations (MSMD), fragment-based molecular docking, and pharmacophore model docking coupled with molecular dynamics simulations─to identify potential binding sites and provide comprehensive insights into TcP21. The three methodologies converged on a common site located on the external surface of the protein, characterized by key residues such as GLU55, ASP52, VAL70, ILE62, and TRP77. The protonated amino, acetamido, and phenyl groups of the pharmacophore probe were consistently observed to interact with the site via a network of salt bridges, hydrogen bonds, charge-charge interactions, and alkyl-π interactions, suggesting these groups play a significant role in ligand binding. This study does not aim to propose specific therapeutic hits but to highlight a still unknown and unexplored protein involved in T. cruzi cell invasion. In this regard, given the strong correlation between the three distinct approaches used for mapping, we consider this study offers valuable insights for further research into P21 and its role in Chagas disease.

Computational revealing Infigratinib resistance to the FGFR2 N549H and N549K mutations

Fibroblast growth factor receptor 2 (FGFR2) is a receptor tyrosine kinase that is involved in many human cancers such as intrahepatic cholangiocarcinomas and hepatocellular carcinomas. The clinically acquired FGFR2 N549H/K mutations make the treatment of Infigratinib ineffective. Here, molecular docking, multiple-replica molecular dynamics (MD) simulations, and binding free energy calculations were used to reveal the mechanism of Infigratinib resistance to the FGFR2 N549H/K mutants. MD simulations indicated that both N549H/K mutations disrupt the local hydrogen bonds or salt bridges formed by His544, Asn549, Glu565, and Lys641. Binding free energy calculation results suggested that both the decrease of van der Waals interactions and electrostatic interaction are responsible for weakening the binding affinity of Infigratinib to the FGFR2 N549H/K mutants. Residue-based energy decomposition analysis further showed that the interactions of Infigratinib with key residues Leu487, Glu565, Ala567, Gly570, and Leu633 become significantly weaker in the N549H/K mutants compared to the wild-type kinase. We anticipate that our results will provide a useful guide to design potential inhibitors against the FGFR2 N549H/K mutants.

Understanding Acute Hemolytic Anemia Severity Through Computational Analysis of G6PDChatham Variant: Designing a New Activator as a Potential Drug

Glucose-6-phosphate dehydrogenase deficiency (G6PDD) is a major enzymatic disease affecting human red blood cells (RBCs), causing hemolytic anemia due to the diminish of the nicotinamide adenine dinucleotide phosphate hydrogen (NADPH) synthesis and altered redox balance within erythrocytes. This study sought to correct the defect in G6PDChatham (Ala355Thr) caused by the loss of interactions (hydrogen bonds and salt bridges) by docking the AG1 molecule at the dimer interface, thus restoring these lost interactions. The enzyme conformation was then analyzed before and after AG1 binding using molecular dynamics simulation (MDS). The reasons behind the severity of acute hemolytic anemia (AHA) were explained using several parameters, such as root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), hydrogen bonds, salt bridges, radius of gyration (Rg), solvent accessible surface area (SASA), and covariance matrix analysis. Structural alterations in G6PDChatham, including the absence of interactions in a key region of the variant structure, can significantly impact protein stability and function, subsequently contributing to disease severity. Upon AG1 binding, these missing interactions were resorted to correct the structural defect of the variant. This restoration improves dimer stability and restores G6PD function. To develop new G6PD activators, several new analogues (SY7, SY8, SY9, and SY10) were rationally developed by substituting the linker region of the AG1 structure with other functional groups using the Avogadro software. These compounds were successfully synthesized and docked with G6PDChatham where the best binding affinity ranged between -8.0 and -9.1 kcal/mol. SY8, a promising G6PD activator, is predicted to be easily metabolized and excreted, making it less likely to cause toxicity. Its high drug score, drug-likeness, and favorable safety profile make it a strong candidate for synthesis and cellular testing. The toxicity risk assessment supported the overall drug score, increasing confidence in finding additional small-molecule activators for G6PDD disorder. Amidst the absence of effective treatments, such discovery hopes to improve the lives of those with AHA by assisting the development of appropriate pharmaceuticals for G6PDD.

Abstract A015: Design and deciphering of precision peptide inhibitors for cancer stemness using generative deep learning and molecular dynamics simulations

Abstract The use of artificial intelligence (AI) and machine learning (ML) for exploring vast amino acid sequence spaces has recently gained traction in the discovery of antibiotics and the design of biomaterials with a wide array of algorithms. However, despite some advancements, designing peptide inhibitors to specifically modulate protein-protein interactions remains a significant challenge. In this contribution, we explore the use of a Long Short-Term Memory (LSTM) network - a type of recurrent neural network - to model peptide sequences, given its ability to process sequential data and capture long-term dependencies, critical for peptide design. Our research focuses on developing precision peptides that target the interaction between the Notch intracellular domain (NICD) and CBF1/RBPJ transcription factors, key regulators of the Notch signaling pathway implicated in breast and pancreatic cancer stemness. We inferred hydrophobic, hydrophilic, van der Waals, and salt bridge interactions from experimentally determined three-dimensional protein complex structures, which were used for feature engineering via one-hot encoding. Peptides, each 20 amino acids in length, were generated using temperature scaling in the LSTM model. These peptides were then structurally optimized and subjected to molecular dynamics (MD) simulations, followed by molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) analysis to assess their interactions and binding affinities with the Notch receptor. The MD simulations provide valuable molecular-level insights into the peptide-Notch interactions, helping to evaluate their binding strength. Further biological testing is underway to validate the efficacy of these lead peptide inhibitors and elucidate their molecular mechanisms in targeting cancer stem cells associated with breast cancer. Citation Format: Gurudeeban Selvaraj, Satyavani Kaliamurthi, Gilles H. Peslherbe. Design and deciphering of precision peptide inhibitors for cancer stemness using generative deep learning and molecular dynamics simulations. [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Optimizing Therapeutic Efficacy and Tolerability through Cancer Chemistry; 2024 Dec 9-11; Toronto, Ontario, Canada. Philadelphia (PA): AACR; Mol Cancer Ther 2024;23(12_Suppl):Abstract nr A015

Binding Interaction Between Two Mutant Myocilin Olfactomedin Domain Monomers in a Homodimer.

In myocilin-associated glaucoma, pathogenic missense mutations accumulate mainly in the olfactomedin domain (mOLF) of myocilin. This makes the protein susceptible to aggregation, where mOLF-mOLF dimerization is possibly an initial stage. Nevertheless, there are no molecular level studies that have probed the nature of interactions occurring between two mOLF domains and the key characteristics of the resulting dimer complex. In this work, we used AlphaFold2 to obtain an I477N mutant mOLF structure with high quality followed by a stable I477N mOLF-mOLF homodimer model using molecular docking combined with molecular dynamics simulations. Moreover, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) methods coupled with per-residue energy decomposition studies are carried out to identify the key residues involved in the binding interaction. Based on these results, we provide insights into the molecular level understanding of the intermolecular interaction between two mOLF domains in an I477N homodimer. Hydrogen bonds, salt bridges, and favorable van der Waals interactions are observed in the binding interface of the homodimer. Additionally, our results suggest that I477N mutant mOLF aggregation could be a multistep process, beginning with an initial mOLF-mOLF dimerization mainly mediated by residues such as Asp395 and Arg681. Also, the peptides P1 (residues 326-337) and P3 (residues 426-442) of the mOLF domain, previously identified as pertinent for myocilin aggregation, could potentially contribute to a subsequent stage of myocilin aggregation, the first step being mOLF-mOLF dimerization.

Probing the Conformational Ensemble of the Amyloid Beta 16-22 Fragment with Parallel-Bias Metadynamics.

Aβ(16-22) is a segment of the Alzheimer's-related β-amyloid peptide that plays a crucial role in its aggregation. This study applies well-tempered parallel-bias metadynamics to investigate the impact of several denaturants and osmolytes on the conformational ensembles of both termini-capped and uncapped Aβ(16-22) monomers. Comparison of the different sets of collective variables in the metadynamics bias shows that using the set of backbone torsional angles results in better and faster convergence of simulations than employing more general structural characteristics of the short peptide. The equilibrium conformational ensembles of the peptides are characterized in pure water and in the presence of TMAO, urea, guanidinium chloride, and trifluoroethanol. In particular, trifluoroethanol and TMAO are found to increase the population of compact peptide conformations, whereas urea and guanidinium chloride favor extended structures. The analysis of the free energy surfaces in the presence of various substances with a comparison of the behavior of the capped and uncapped peptide forms reveals the role of different types of intrapeptide interactions such as salt bridges, hydrophobic contacts, and hydrogen bonds in stabilization of the compact or extended structures. As compounds reducing the abundance of the compact states of Aβ(16-22) and other disordered peptides are likely to suppress their amyloid fibril formation, simulations in the systems with this short peptide may be useful for the virtual screening of such compounds.

Mechanistic insights into allosteric regulation of the reductase component of p-hydroxyphenylacetate 3-hydroxylase by p-hydroxyphenylacetate: a model for effector-controlled activity of redox enzymes.

Understanding how an enzyme regulates its function through substrate or allosteric regulation is crucial for controlling metabolic pathways. Some flavin-dependent monooxygenases (FDMOs) have evolved an allosteric mechanism to produce reduced flavin while minimizing the use of NADH and the production of harmful hydrogen peroxide (H2O2). In this work, we investigated in-depth mechanisms of how the reductase component (C1) of p-hydroxyphenylacetate (HPA) 3-hydroxylase (HPAH) from Acinetobacter baumanii is allosterically controlled by the binding of HPA, which is a substrate of its monooxygenase counterpart (C2). The C1 structure can be divided into three regions: the N-terminal domain (flavin reductase); a flexible loop; and the C-terminal domain, which is homologous to NadR, a repressor protein having HPA as an effector. The binding of HPA to NadR induces a conformational change in the recognition helix, causing it to disengage from the NadA gene. The HPA binding site of C1 is located at the C-terminal domain, which can be divided into five helices. Molecular dynamics simulations performed on HPA-docked C1 elucidated the allosteric mechanism. The carboxylate group of HPA maintains the salt bridge between helix 2 and the flexible loop. This maintenance shortens the loop between helices 2 and 3, causing helix 3 to disengage from the N-terminal domain. The aromatic ring of HPA induces a conformational change in helices 1 and 5, pulling helix 4, analogous to the recognition helix in NadR, away from the N-terminal domain. A Y189A mutation, obtained from site-saturation mutagenesis, confirms that HPA with an unsuitable conformation cannot induce the conformational change of C1. Additionally, an HPA-independent effect is observed, in which Arg20, an NADH binding residue on the N-terminal domain, occasionally disengages from helix 4. This model provides valuable insights into the allosteric regulation of two-component FDMOs and aromatic effector systems.

Structural comparison of homologous protein-RNA interfaces reveals widespread overall conservation contrasted with versatility in polar contacts.

Protein-RNA interactions play a critical role in many cellular processes and pathologies. However, experimental determination of protein-RNA structures is still challenging, therefore computational tools are needed for the prediction of protein-RNA interfaces. Although evolutionary pressures can be exploited for structural prediction of protein-protein interfaces, and recent deep learning methods using protein multiple sequence alignments have radically improved the performance of protein-protein interface structural prediction, protein-RNA structural prediction is lagging behind, due to the scarcity of structural data and the flexibility involved in these complexes. To study the evolution of protein-RNA interface structures, we first identified a large and diverse dataset of 2,022 pairs of structurally homologous interfaces (termed structural interologs). We leveraged this unique dataset to analyze the conservation of interface contacts among structural interologs based on the properties of involved amino acids and nucleotides. We uncovered that 73% of distance-based contacts and 68% of apolar contacts are conserved on average, and the strong conservation of these contacts occurs even in distant homologs with sequence identity below 20%. Distance-based contacts are also much more conserved compared to what we had found in a previous study of homologous protein-protein interfaces. In contrast, hydrogen bonds, salt bridges, and π-stacking interactions are very versatile in pairs of protein-RNA interologs, even for close homologs with high interface sequence identity. We found that almost half of the non-conserved distance-based contacts are linked to a small proportion of interface residues that no longer make interface contacts in the interolog, a phenomenon we term "interface switching out". We also examined possible recovery mechanisms for non-conserved hydrogen bonds and salt bridges, uncovering diverse scenarios of switching out, change in amino acid chemical nature, intermolecular and intramolecular compensations. Our findings provide insights for integrating evolutionary signals into predictive protein-RNA structural modeling methods.

Characterization of molecular interactions between HDAC7 and MEF2A

Interactions of transcriptional corepressors such as histone deacetylase 7 (HDAC7), a class IIa HDAC, with myocyte enhancer factor-2 (MEF2) regulate MEF2 activity. Despite previous investigations exploring interactions between HDAC7 and MEF2, a detailed characterization of the HDAC7-MEF2 functional complex is still lacking. Herein, we first modeled the structure of the HDAC7-MEF2A complex and investigated the inter-protein interactions using all-atom molecular dynamics (MD) simulations. We identified specific amino acids within HDAC7 and MEF2A that participate in interactions such as salt bridges, hydrogen bonds, and hydrophobic interactions. Our results reveal a salt bridge formed between LYS96(HDAC7) and ASP63(MEF2A). Our analysis also predicted formations of reliable hydrogen bonds between SER82(HDAC7) and ASP63(MEF2A) as well as LYS96(HDAC7) and ASP63(MEF2A). In addition, clustering of hydrophobic residues at the interface contributes in stabilizing the HDAC7-MEF2A complex. Results from multiple sequence alignment show that most of the HDAC7 residues that are predicted to associate with MEF2A are conserved in at least three class IIa HDACs and all predicted residues in MEF2A are conserved in MEF2s. We also found that the association of DNA to MEF2A has no significant effect on HDAC7-MEF2A interactions. Our results may also provide useful insights into the interactions between other class IIa HDACs and MEF2s.