Drude2019IDPC polarizable force field reveals structure-function relationship of insulin.

Intrinsically disordered proteins (IDPs) lack stable tertiary structures under physiological conditions, yet play key roles in biological processes and associated with human complex diseases. Their conformational characteristics and high content of charged residues make the use of polarizable force fields an advantageous for simulating IDPs. The Drude2019IDP polarizable force field, previously introduced, has demonstrated comprehensive enhancements and improvements in dipeptides, short peptides, and IDPs, achieving a balanced sampling between IDPs and structured proteins. However, the performance in simulating 5 dipeptides was found to be underestimate. Therefore, we individually performed reweighting and grid-based energy correction map (CMAP) optimization for these 5 dipeptides, resulting in the enhanced Drude2019IDPC force field. The performance of Drude2019IDPC was evaluated with 5 dipeptides, 5 disordered short peptides, and a representative IDP. The results demonstrated a marked improvement comparing with original Drude2019IDP. To further substantiate the capabilities of Drude2019IDPC, MD simulation and Markov state model (MSM) were applied to wild type and mutant for insulin, to elucidate the difference of conformational characteristics and transition path. The findings reveal that mutation can maintain the monomorphic characteristics, providing insights for engineered insulin development. These results indicate that Drude2019IDPC could be used to reveal the structure-function relationship for other proteins.

Early Aggregation Mechanism of SOD128-38 Based on Force Field Parameter of 5-Cyano-Tryptophan.

The aggregation of superoxide dismutase 1 (SOD1) results in amyloid deposition and is involved in familial amyotrophic lateral sclerosis, a fatal motor neuron disease. There have been extensive studies of its aggregation mechanism. Noncanonical amino acid 5-cyano-tryptophan (5-CN-Trp), which has been incorporated into the amyloid segments of SOD1 as infrared probes to increase the structural sensitivity of IR spectroscopy, is found to accelerate the overall aggregation rate and potentially modulate the aggregation process. Despite these observations, the underlying mechanism remains elusive. Here, we optimized the force field parameters of 5-CN-Trp and then used molecular dynamics simulation along with the Markov state model on the SOD128-38 dimer to explore the kinetics of key intermediates in the presence and absence of 5-CN-Trp. Our findings indicate a significantly increased probability of protein aggregate formation in 5CN-Trp-modified ensembles compared to wildtype. Dimeric ß-sheets of different natures were observed exclusively in the 5CN-Trp-modified peptides, contrasting with wildtype simulations. Free-energy calculations and detailed analyses of the dimer structure revealed augmented interstrand interactions attributed to 5-CN-Trp, which contributed more to peptide affinity than any other residues. These results explored the key events critical for the early nucleation of amyloid-prone proteins and also shed light on the practice of using noncanonical derivatives to study the aggregation mechanism.

Substrate DNA Promoting Binding of Mycobacterium tuberculosis MtrA by Facilitating Dimerization and Interpretation of Affinity by Minor Groove Width.

In order to deepen the understanding of the role and regulation mechanisms of prokaryotic global transcription regulators in complex processes, including virulence, the associations between the affinity and binding sequences of Mycobacterium tuberculosis MtrA have been explored extensively. Analysis of MtrA 294 diversified 26 bp binding sequences revealed that the sequence similarity of fragments was not simply associated with affinity. The unique variation patterns of GC content and periodical and sequential fluctuation of affinity contribution curves were observed along the sequence in this study. Furthermore, docking analysis demonstrated that the structure of the dimer MtrA-DNA (high affinity) was generally consistent with other OmpR family members, while Arg 219 and Gly 220 of the wing domain interacted with the minor groove. The results of the binding box replacement experiment proved that box 2 was essential for binding, which implied the differential roles of the two boxes in the binding process. Furthermore, the results of the substitution of the nucleotide at the 20th and/or 21st positions indicated that the affinity was negatively associated with the value of minor groove width precisely at the 21st position. The dimerization of the unphosphorylated MtrA facilitated by a low-affinity DNA fragment was observed for the first time. However, the proportion of the dimer was associated with the affinity of substrate DNA, which further suggested that the affinity was actually one characteristic of the stability of dimers. Based on the finding of 17 inter-molecule hydrogen bonds identified in the interface of the MtrA dimer, including 8 symmetric complementary ones in the conserved α4-ß5-α5 face, we propose that hydrogen bonds should be considered just as important as salt bridges and the hydrophobic patch in the dimerization. Our comprehensive study on a large number of binding fragments with quantitative affinity values provided new insight into the molecular mechanism of dimerization, binding specificity and affinity determination of MtrA and clues for solving the puzzle of how global transcription factors regulate a large quantity of target genes.

Phosphorylation Regulation Mechanism of ß2 Integrin for the Binding of Filamin Revealed by Markov State Model.

Leukocyte adhesion deficiency-1 (LAD-1) disorder is a severe immunodeficiency syndrome caused by deficiency or mutation of ß2 integrin. The phosphorylation on threonine 758 of ß2 integrin acts as a molecular switch inhibiting the binding of filamin. However, the switch mechanism of site-specific phosphorylation at the atom level is still poorly understood. To resolve the regulation mechanism, all-atom molecular dynamics simulation and Markov state model were used to study the dynamic regulation pathway of phosphorylation. Wild type system possessed lower binding free energy and fewer number of states than the phosphorylated system. Both systems underwent local disorder-to-order conformation conversion when achieving steady states. To reach steady states, wild type adopted less number of transition paths/shortest path according to the transition path theory than the phosphorylated system. The underlying phosphorylated regulation pathway was from P1 to P0 and then P4 state, and the main driving force should be hydrogen bond and hydrophobic interaction disturbing the secondary structure of phosphorylated states. These studies will shed light on the pathogenesis of LAD-1 disease and lay a foundation for drug development.

Structural Insight of KSIII (ß-Ketoacyl-ACP Synthase)-like Acyltransferase ChlB3 in the Biosynthesis of Chlorothricin.

Chlorothricin (CHL) belongs to a spirotetronate antibiotic family produced by Streptomyces antibioticus that inhibits pyruvate carboxylase and malate dehydrogenase. For the biosynthesis of CHL, ChlB3 plays a crucial role by introducing the 6-methylsalicylic acid (6MSA) moiety to ChlB2, an acyl carrier protein (ACP). However, the structural insight and catalytic mechanism of ChlB3 was unclear. In the current study, the crystal structure of ChlB3 was solved at 3.1 Å-resolution and a catalytic mechanism was proposed on the basis of conserved residues of structurally related enzymes. ChlB3 is a dimer having the same active sites as CerJ (a structural homologous enzyme) and uses a KSIII-like fold to work as an acyltransferase. The relaxed substrate specificity of ChlB3 was defined by its catalytic efficiencies (kcat/Km) for non-ACP tethered synthetic substrates such as 6MSA-SNAC, acetyl-SNAC, and cyclohexonyl-SNAC. ChlB3 successfully detached the 6MSA moiety from 6MSA-SNAC substrate and this hydrolytic activity demonstrated that ChlB3 has the potential to catalyze non-ACP tethered substrates. Structural comparison indicated that ChlB3 belongs to FabH family and showed 0.6-2.5 Å root mean square deviation (RMSD) with structural homologous enzymes. Molecular docking and dynamics simulations were implemented to understand substrate active site and structural behavior such as the open and closed conformation of the ChlB3 protein. The resultant catalytic and substrate recognition mechanism suggested that ChlB3 has the potential to use non-native substrates and minimize the labor of expressing ACP protein. This versatile acyltransferase activity may pave the way for manufacturing CHL variants and may help to hydrolyze several thioester-based compounds.

Early aggregation mechanism of Aß16-22 revealed by Markov state models.

Aß16-22 is believed to have critical role in early aggregation of full length amyloids that are associated with the Alzheimer's disease and can aggregate to form amyloid fibrils. However, the early aggregation mechanism is still unsolved. Here, multiple long-term molecular dynamics simulations combining with Markov state model were used to probe the early oligomerization mechanism of Aß16-22 peptides. The identified dimeric form adopted either globular random-coil or extended ß-strand like conformations. The observed dimers of these variants shared many overall conformational characteristics but differed in several aspects at detailed level. In all cases, the most common type of secondary structure was intermolecular antiparallel ß-sheets. The inter-state transitions were very frequent ranges from few to hundred nanoseconds. More strikingly, those states which contain fraction of ß secondary structure and significant amount of extended coiled structures, therefore exposed to the solvent, were majorly participated in aggregation. The assembly of low-energy dimers, in which the peptides form antiparallel ß sheets, occurred by multiple pathways with the formation of an obligatory intermediates. We proposed that these states might facilitate the Aß16-22 aggregation through a significant component of the conformational selection mechanism, because they might increase the aggregates population by promoting the inter-chain hydrophobic and the hydrogen bond contacts. The formation of early stage antiparallel ß sheet structures is critical for oligomerization, and at the same time provided a flat geometry to seed the ordered ß-strand packing of the fibrils. Our findings hint at reorganization of this part of the molecule as a potentially critical step in Aß aggregation and will insight into early oligomerization for large ß amyloids.

Disaggregation mechanism of prion amyloid for tweezer inhibitor.

The aggregation of amyloid has been an important event in the pathology of amyloidogenicity. A number of small molecules have been designed for Amyloidosis treatment. Molecular tweezer CLR01, a potential drug for misfolded ß-amyloids inhibition, was reportedly bind directly to Lysine residues and interrupt oligomerization. However, the disaggregation mechanism of amyloid for this inhibitor is unclear. Here we used long timescale of molecular dynamic simulation to reveal the mechanism of disaggregation for pentamer prion amyloid. Molecular docking and molecular dynamics simulation demonstrate that CLR01 is attached with Lysine222 nitrogen by π-cation interaction of its nine aromatic rings and formation of salt bridge/hydrogen bond of one of the two rotatable peripheral anionic phosphate groups. Upon CLR01 binding, we found a major shifting occurs in initial conformation of the oligomer and stretch out the N-terminal chain A from the rest of the amyloid which seems to be the first stage of disaggregated the fibrils slowly yet efficiently. Moreover, the CLR01 remodelled the pentamer Prion220-272 into a compact structure which might be the resistant conformation for further oligomerization. Our work will contribute to better understand the interaction and deterioration mechanism of molecular tweezer for prions and similar amyloids, and offer significant insights into therapeutic development for Amyloidosis treatment.

Comparison and Evaluation of Force Fields for Intrinsically Disordered Proteins.

Molecular dynamics (MD) simulations of six upgraded empirical force fields were compared and evaluated with short peptides, intrinsically disordered proteins, and folded proteins using trajectories of 1, 1.5, 5, or 10 µs (five replicates of 200 ns, 300 ns, 1 µs, or 2 µs) for each system. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. Here, the MD simulations were done in an explicit solvent with RS-peptide, HEWL19, HIV-rev, ß amyloid (Aß)-40, Aß-42, phosphodiesterase-γ, CspTm, and ubiquitin using ff99IDPs, ff14IDPs, ff14IDPSFF, ff03w, CHARMM36m, and CHARMM22* force fields. The IDP ensembles generated by six all-atom empirical force fields were compared against NMR data. Despite using identical starting structures and simulation parameters, ensembles obtained with different force fields exhibit significant differences in NMR RMDs, secondary structure contents, and global properties such as the radius of gyration. The intrinsically disordered protein (IDP)-specific force fields could substantially reproduce the experimental observables in force field comparison: they have the lowest error in chemical shifts and J-couplings for short peptides/proteins, reasonably well for large IDPs and reasonably well with the radius of gyration. A high population of disorderness was observed in the IDP-specific force field for the IDP ensemble with a fraction of ß sheets for ß-amyloids. CHARMM22* performs better for many observables; however, it still has a preference toward the helicity for short peptides. The results of ß-amyloid 42 starting from two different initial structures (Aß421Z0Q and Aß42model) were also compared with DSSP and NMR data. The results obtained with IDP-specific force fields within 2 µs simulation time are similar, even though starting from different structures. The current force fields perform equally well for folded proteins. The results of currently developed or modified force fields for IDPs are capable of enlightening the overall performance of the force field for disordered as well as folded proteins, thereby contributing to force field development.

Decoding allosteric communication pathways in protein lysine acetyltransferase.

In bacteria, protein lysine acetylation circuits can control core processes such as carbon metabolism. In E. coli, cyclic adenosine monophosphate (cAMP) controls the transcription level and activity of protein lysine acetyltransferase (PAT). The M. tuberculosis PatA (Mt-PatA) resides in two different conformations; the activated state and autoinhibited state. However, the mechanism of cAMP allosteric regulation of Mt-PatA remains mysterious. Here, we performed extensive all-atom molecular dynamics (MD) simulations (three independent run for each system) and built a residue-residue dynamic correlation network to show how cAMP mediates allosteric activation. cAMP binds at the regulatory site in the regulatory domain, which is 32 Å away from the catalytic site. An extensive conformational restructuring relieves autoinhibition caused by a molecular Lid (residues 161-203) that shelters the substrate-binding surface. In the activated state, the regulatory domain rotates (~40°) around Ser144, which links both domains. Rotation removes the C-terminus from the cAMP site and relieves the autoinhibited state. Also, the molecular Lid refolds and creates an activator binding site. A conserved residue, His173, was mutated into Lys in the Lid, and during an MD trajectory of the activated state, positioned itself near an acetyl donor molecule in the catalytic domain, suggesting a direct mechanism for acetylation. This study describes the allosteric framework for Mt-PatA and prerequisite intermediate states that permit long-distance signal transmission.

Allosteric Modulation of Intrinsically Disordered Proteins.

The allosteric property of globular proteins is applauded as their intrinsic ability to regulate distant sites, and this property further plays a critical role in a wide variety of cellular regulatory mechanisms. Recent advancements and studies have revealed the manifestation of allostery in intrinsically disordered proteins or regions as allosteric sites present within or mediated by IDP/IDRs facilitates the signaling interactions for various biological mechanisms which would otherwise be impossible for globular proteins to regulate. This thematic review has highlighted the biological outcomes that can be achieved by the mechanism of allosteric regulation of intrinsically disordered proteins or regions. The similar mechanism has been implemented on Adenovirus 5 early region 1A and tumor apoptosis protein p53 in correspondence with other partners in binary and ternary complexes, which are the subject of the current review. Both these proteins regulate once they bind to their partners, consequently, forming either a binary or a ternary complex. Allosteric regulation by IDPs is currently a subject undergoing intense study, and the ongoing research work will ensure a better understanding of precision and efficiency of cellular regulation by them. Allosteric regulation mechanism can also be researched by intrinsically disordered protein-specific force field.

Gain-of-Function SHP2 E76Q Mutant Rescuing Autoinhibition Mechanism Associated with Juvenile Myelomonocytic Leukemia.

Juvenile myelomonocytic leukemia (JMML) is an invasive myeloproliferative neoplasm and is a childhood disease with very high clinical lethality. The SHP2 is encoded by the PTPN11 gene, which is a nonreceptor (pY)-phosphatase and mutation causes JMML. The structural hierarchy of SHP2 includes protein tyrosine phosphatase domain (PTP) and Src-homology 2 domain (N-SH2 and C-SH2). Somatic mutation (E76Q) in the interface of SH2-PTP domain is the most commonly identified mutation found in up to 35% of patients with JMML. The mechanism of this mutant associated with JMML is poorly understood. Here, molecular dynamics simulation was performed on wild-type and mutant (E76Q) of SHP2 to explore the precise impact of gain-of-function on PTP's activity. Consequently, such impact rescues the SHP2 protein from autoinhibition state through losing the interface interactions of Q256/F7 and S502/Q76 or weakening interactions of Q256/R4, Q510/G60, and Q506/A72 between N-SH2 and PTP domains. The consequences of these interactions further relieve the D'E loop away from the PTP catalytic site. The following study would provide a mechanistic insight for better understanding of how individual SHP2 mutations alter the PTP's activity at the atomic level.

The Landscape of Protein Tyrosine Phosphatase (Shp2) and Cancer.

Role of Shp2: The dysregulation of cell signaling cascades associated with the cell differentiation and growth, due to the deletion, insertion or point mutation in specific amino acids which alters the intrinsic conformation of the protein, can ultimately lead to a fatal cancer disease. The protein tyrosine phosphatase has been recognized as a key regulator of extracellular stimuli such as cytokine receptor and receptor tyrosine kinase signaling. In the last era, the PTPN11 gene (encode a Shp2 protein) and its association with acute myeloid, juvenile myelomonocytic, and B-cell acute lymphoblastic leukemia, Noonan syndrome, and myelodysplastic have been recognized as the cause of such deadly disease due to the occurrence of germline mutations in the interface of PTP and SH2 domain. Conclusion: The current study was designed to focus on the allosteric regulation (autoinhibition) of the of Shp2 protein. Subsequently, it will cover the last 10-year recap of Shp2 protein, their role in cancer, and regulation in numerous ways (allosteric regulation).

Inhibitory mechanism of 5-bromo-3-indoleacetic acid for non-structural-3 helicase hepatitis C virus with dynamics correlation network analysis.

Hepatitis C virus, Nonstructural 3 helicase (NS3h) protein is a well-studied segment of Non-structural 3/4 A helicase-protease protein that is crucial for the RNA duplex unwinding and RNA translocation during the process of HCV replication. Similar to other HCV nonstructural proteins, helicase is a potential target for antiviral drugs and several antiviral molecules have been used to target the RNA-binding cleft, despite the fact that none of those helicase antivirals have advanced the clinical trials. Compound 2t9 (5-bromo-1H-indol-3-yl acetic acid) has been identified through the integrated strategies and considered as a potential lead compound for the inactivation of HCV helicase. This inhibitor bind to the 3´-terminal RNA-binding cleft, and reduced the RNA binding and unwinding activity of the targeted protein. In the current study, using all-atom molecular dynamic simulation and correlations network strategy, we scrutinized the inhibitory mechanism of compound 2t9 that needs to be elucidated for the improvement of indole-based and similar HCV helicase inhibitors. Consequently, by comparing the structural dynamics of free (NS3hWT) and bound (NS3h/2t9WT) protein, we identified that the inhibitor-bound protein achieved a conformation resemblance to the open conformation, where the RNA is displaced results in destabilization of RNA-binding cleft, disruption in ATP/ADP binding site and alter the inter-domain communication. The results were evaluated by using the W501 L mutated system. The information based on detailed dynamic aspects of the drug targeted protein will facilitate the researchers in the development of HCV antiviral drugs.

Allosteric mechanism of cyclopropylindolobenzazepine inhibitors for HCV NS5B RdRp via dynamic correlation network analysis.

HCV RNA dependent RNA polymerase (RdRp) nonstructural protein 5B (NS5B) is a major target against hepatitis C virus (HCV) for antiviral therapy. Recently discovered cyclopropylindolobenzazepine derivatives have been considered as the most potent for their ability to bind the thumb site 1 domain and allosterically inhibit HCV NS5B RdRp activity. However, the allosteric mechanism for these derivatives has not been clarified at the molecular level. In this study, fluctuation correlation networks were constructed based on all-atom molecular dynamics simulations to elucidate the allosteric mechanism. The fluctuation correlation networks between free and M2 bound NS5B are significantly different. Information can better transfer from the allosteric site to the catalytic site for bound NS5B than for free NS5B. Thus, the hypothesis of "binding induced allosteric regulation" is proposed to link the enzyme activation and inhibitor binding and then confirmed by the mutant network. Finally, one possible allosteric pathway was identified with the shortest path and evaluated by the perturbation of the network. These methods will be helpful to identify the allosteric pathway of other proteins and to design new drugs targeting the pathway.