Aggregation mechanism investigation of the GIFQINS cross-β amyloid fibril

Catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase

Methane is among the most potent of the greenhouse gases, which plays a key role in global climate change. As an excellent carbon and energy source, methane can be utilized by anaerobic methane oxidizing archaea and aerobic methane oxidizing bacteria. The previous work shows that an anaerobic thermophilic enrichment culture composed of dense consortia of archaea and bacteria apparently uses partly similar pathways to oxidize the C4 hydrocarbon butane. However, the catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase is still unknown. Therefore, molecular dynamics (MD) simulation was used to investigate the dynamics differences of catalytic mechanism between methane coenzyme M reductase (MCR) and alkyl-coenzyme M reductase (ACR). At first, the binding pocket of ACR is larger than that of MCR. Then, the complex of butane and ACR is more stable than that of methane and ACR. Protein conformation cloud suggests that the position of methane is dynamics and methane escapes from the binding pocket of ACR during most of the simulation time, while butane tightly binds in the pocket of ACR. The hydrophobic interactions between butane and ACR are more and stronger than those between methane and ACR. At the same time, the binding free energy between butane and ACR is significantly lower than that between methane and ACR. The dynamics correlation network indicates that the transformation of information flow for ACR-butane is smoother than that for ACR-methane. The shortest pathway for ACR-butane is from Gln144, Ala141, Hie135, Ile133, Ala160, Arg206, Asp97, Met94, Tyr347 to Phe345 with synergistic effect for two butane molecules. This study can insight into the catalytic mechanism for butane/ACR complex.

Selectivity Mechanism of ATP-Competitive Inhibitors for PKB and PKA

Protein kinase B (PKB) acts as a central node on the PI3K kinase pathway. Constitutive activation and overexpression of PKB have been identified to involve in various cancers. However, protein kinase A (PKA) sharing high homology with PKB is essential for metabolic regulation. Therefore, specific targeting on PKB is crucial strategy in drug design and development for antitumor. Here, we had revealed the selectivity mechanism for PKB inhibitors with molecular dynamics simulation and 3D-QSAR methods. Selective inhibitors of PKB could form more hydrogen bonds and hydrophobic contacts with PKB than those with PKA. This could explain that selective inhibitor M128 is more potent to PKB than to PKA. Then, 3D-QSAR models were constructed for these selective inhibitors and evaluated by test set compounds. 3D-QSAR model comparison of PKB inhibitors and PKA inhibitors reveals possible methods to improve the selectivity of inhibitors. These models can be used to design new chemical entities and make quantitative prediction of the specific selective inhibitors before resorting to in vitro and in vivo experiment.

Interplay of sequence, topology and termini charge in determining the stability of the aggregates of GNNQQNY mutants: a molecular dynamics study

This study explores the stabilities of single sheet parallel systems of three sequence variants of ¹GNNQQNY⁷, N2D, N2S and N6D, with variations in aggregate size (5–8) and termini charge (charged or neutral). The aggregates were simulated at 300 and 330 K. These mutations decrease amyloid formation in the yeast prion protein Sup35. The present study finds that these mutations cause instability even in the peptide context. The protonation status of termini is found to be a key determinant of stabilities; other determinants are sequence, position of mutation and aggregate size. All systems with charged termini are unstable, whereas both stable and unstable systems are found when the termini are neutral. When termini are charged, the largest stable aggregate for the N2S and N6D systems has 3 to 4 peptides whereas N2D mutation supports oligomers of larger size (5-and 6-mers) as well. Mutation at 2^nd position (N2S and N2D) results in fewer H-bonds at the mutated as well as neighboring (Gly1/Gln4) positions. However, no such effect is found if mutation is at 6^th position (N6D). The effect of Asn→Asp mutation depends on the position and termini charge: it is more destabilizing at the 2^nd position than at the 6^th in case of neutral termini, however, the opposite is true in case of charged termini. Appearance of twist in stable systems and in smaller aggregates formed in unstable systems suggests that twist is integral to amyloid arrangement. Disorder, dissociation or rearrangement of peptides, disintegration or collapse of aggregates and formation of amorphous aggregates observed in these simulations are likely to occur during the early stages of aggregation also. The smaller aggregates formed due to such events have a variety of arrangements of peptides. This suggests polymorphic nature of oligomers and presence of a heterogeneous mixture of oligomers during early stages of aggregation.

Molecular dynamics simulations of amyloid fibrils: an in silico approach

Amyloid fibrils play causal roles in the pathogenesis of amyloid-related degenerative diseases such as Alzheimer's disease, type II diabetes mellitus, and the prion-related transmissible spongiform encephalopathies. The mechanism of fibril formation and protein aggregation is still hotly debated and remains an important open question in order to develop therapeutic method of these diseases. However, traditional molecular biological and crystallographic experiments could hardly observe atomic details and aggregation process. Molecular dynamics (MD) simulations could provide explanations for experimental results and detailed pathway of protein aggregation. In this review, we focus on the applications of MD simulations on several amyloidogenic protein systems. Furthermore, MD simulations could help us to understand the mechanism of amyloid aggregation and how to design the inhibitors.

Atomistic mechanism of microRNA translation upregulation via molecular dynamics simulations

MicroRNAs are endogenous 23–25 nt RNAs that play important gene-regulatory roles in animals and plants. Recently, miR369-3 was found to upregulate translation of TNFα mRNA in quiescent (G0) mammalian cell lines. Knock down and immunofluorescence experiments suggest that microRNA-protein complexes (with FXR1 and AGO2) are necessary for the translation upregulation. However the molecular mechanism of microRNA translation activation is poorly understood. In this study we constructed the microRNA-mRNA-AGO2-FXR1 quadruple complex by bioinformatics and molecular modeling, followed with all atom molecular dynamics simulations in explicit solvent to investigate the interaction mechanisms for the complex. A combined analysis of experimental and computational data suggests that AGO2-FXR1 complex relocalize microRNA:mRNA duplex to polysomes in G0. The two strands of dsRNA are then separated upon binding of AGO2 and FXR1. Finally, polysomes may improve the translation efficiency of mRNA. The mutation research confirms the stability of microRNA-mRNA-FXR1 and illustrates importance of key residue of Ile304. This possible mechanism can shed more light on the microRNA-dependent upregulation of translation.

Insight into the stability of cross-β amyloid fibril from VEALYL short peptide with molecular dynamics simulation

Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides.

Revealing interaction mode between HIV-1 protease and mannitol analog inhibitor

HIV protease is a key enzyme to play a key role in the HIV-1 replication cycle and control the maturation from HIV viruses to an infectious virion. HIV-1 protease has become an important target for anti-HIV-1 drug development. Here, we used molecular dynamics simulation to study the binding mode between mannitol derivatives and HIV-1 protease. The results suggest that the most active compound (M35) has more stable hydrogen bonds and stable native contacts than the less active one (M17). These mannitol derivatives might have similar interaction mode with HIV-1 protease. Then, 3D-QSAR was used to construct quantitative structure-activity models. The cross-validated q(2) values are found as 0.728 and 0.611 for CoMFA and CoMSIA, respectively. And the non-cross-validated r(2) values are 0.973 and 0.950. Nine test set compounds validate the model. The results show that this model possesses better prediction ability than the previous work. This model can be used to design new chemical entities and make quantitative prediction of the bioactivities for HIV-1 protease inhibitors before resorting to in vitro and in vivo experiment.

Drug resistant mechanism of diaryltriazine analog inhibitors of HIV-1 reverse transcriptase using molecular dynamics simulation and 3D-QSAR

Diaryltriazine inhibitors have highly potent and effective bioactivities for the wild type of HIV-1 reverse transcription. To design new drug of antimutant HIV-1 reverse transcriptase, the mechanism of drug resistance for four types of mutants was revealed. Molecular dynamics simulations suggest that Lys101, Leu100, Lys103, Tyr181, and Tyr188 are key residues. Different mutants of key residues may have different interaction modes and lead to different drug resistances. Then, CoMFA and CoMSIA methods were employed to construct 3D quantitative structure-activity relationship models. These models were evaluated by test set compounds. These models can be used to make quantitative prediction of their bioactivities for lead compounds before resorting to in vitro and in vivo experimentation.

Induced fit or conformational selection for RNA/U1A folding

The hairpin II of U1 snRNA can bind U1A protein with high affinity and specificity. NMR spectra suggest that the loop region of apo-RNA is largely unstructured and undergoes a transition from unstructured to well-folded upon U1Abinding. However, the mechanism that RNA folding coupled protein binding is poorly understood. To get an insight into the mechanism, we have performed explicit-solvent molecular dynamics (MD) to study the folding kinetics of bound RNA and apo-RNA. Room-temperature MD simulations suggest that the conformation of bound RNA has significant adjustment and becomes more stable upon U1A binding. Kinetic analysis of high-temperature MD simulations shows that bound RNA and apo-RNA unfold via a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound RNA folds in the order of RNA contracting, U1A binding, and tertiary folding. The predicted Phi-values suggest that A8, C10, A11, and G16 are key bases for bound RNA folding. Mutant Arg52Gln analysis shows that electrostatic interaction and hydrogen bonds between RNA and U1A (Arg52Gln) decrease. These results are in qualitative agreement with experiments. Furthermore, this method could be used in other studies about biomolecule folding upon receptor binding.

Optimal molecular structures of prion AGAAAAGA amyloid fibrils formatted by simulated annealing

To date, there is little structural data available on the AGAAAAGA palindrome in the hydrophobic region (113-120) of prion proteins, although many experimental studies have shown that this region has amyloid fibril forming properties. This region belongs to the N-terminal unstructured region (1-123) of prions, the structure of which has proved hard to determine using NMR or X-ray crystallography. This paper reports the successful construction of three amyloid fibril models for this region. The models were formatted by standard simulated annealing using suitable templates from the Protein Data Bank, and were refined using several traditional optimization methods within AMBER. Because the NMR or X-ray structure of the hydrophobic region AGAAAAGA of prion proteins has not yet been determined, these models can be used as a reference for experimental studies on this region. The results presented here confirm standard simulated annealing as an effective tool in molecular modeling. The three constructed models for amyloid fibrils may be useful in furthering the goals of medicinal chemistry in this field.

Molecular dynamics simulation of phosphorylated KID post-translational modification

Background

Kinase-inducible domain (KID) as transcriptional activator can stimulate target gene expression in signal transduction by associating with KID interacting domain (KIX). NMR spectra suggest that apo-KID is an unstructured protein. After post-translational modification by phosphorylation, KID undergoes a transition from disordered to well folded protein upon binding to KIX. However, the mechanism of folding coupled to binding is poorly understood.

Methodology

To get an insight into the mechanism, we have performed ten trajectories of explicit-solvent molecular dynamics (MD) for both bound and apo phosphorylated KID (pKID). Ten MD simulations are sufficient to capture the average properties in the protein folding and unfolding.

Conclusions

Room-temperature MD simulations suggest that pKID becomes more rigid and stable upon the KIX-binding. Kinetic analysis of high-temperature MD simulations shows that bound pKID and apo-pKID unfold via a three-state and a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound pKID folds in the order of KIX access, initiation of pKID tertiary folding, folding of helix α_B, folding of helix α_A, completion of pKID tertiary folding, and finalization of pKID-KIX binding. Our data show that the folding pathways of apo-pKID are different from the bound state: the foldings of helices α_A and α_B are swapped. Here we also show that Asn139, Asp140 and Leu141 with large Φ-values are key residues in the folding of bound pKID. Our results are in good agreement with NMR experimental observations and provide significant insight into the general mechanisms of binding induced protein folding and other conformational adjustment in post-translational modification.