Multicanonical Molecular Dynamics Simulations of the N-terminal Domain of Protein L9.
Ontology highlight
ABSTRACT:
SUBMITTER: Yasar F
PROVIDER: S-EPMC4169893 | biostudies-literature | 2014 Feb
REPOSITORIES: biostudies-literature
Publications
Multicanonical Molecular Dynamics Simulations of the N-terminal Domain of Protein L9.
Yaşar Fatih F Jiang Ping P Hansmann Ulrich H E UH
Europhysics letters 20140201 3
We describe multicanonical molecular dynamic simulations of the N-terminal domain of the protein L9. Analyzing free energy landscapes and thermal ordering, we propose a possible folding mechanism for the protein. By comparing our results with that of molecular dynamics runs of the protein at constant temperature, we find that multicanonical molecular dynamics leads to orders of magnitude higher sampling of folding transitions. ...[more]
Similar Datasets
Project description:We compare the efficiency of multicanonical and replica exchange molecular dynamics for the sampling of folding/unfolding events in simulations of proteins with end-to-end β-sheet. In Go-model simulations of the 75-residue MNK6, we observe improvement factors of 30 in the number of folding/unfolding events of multicanonical molecular dynamics over replica exchange molecular dynamics. As an application, we use this enhanced sampling to study the folding landscape of the 36-residue DS119 with an all-atom physical force field and implicit solvent. Here, we find that the rate-limiting step is the formation of the central helix that then provides a scaffold for the parallel β-sheet formed by the two chain ends.
| S-EPMC3815463 | biostudies-literature
Project description:The N-terminal domain of L9 (NTL9) is a 56-residue mixed α-β protein that lacks disulfides, does not bind cofactors, and folds reversibly. NTL9 has been widely used as a model system for experimental and computational studies of protein folding and for investigations of the unfolded state. The role of side-chain interactions in the folding of NTL9 is probed by mutational analysis. ϕ-values, which represent the ratio of the change in the log of the folding rate upon mutation to the change in the log of the equilibrium constant for folding, are reported for 25 point mutations and 15 double mutants. All ϕ-values are small, with an average over all sites probed of only 0.19 and a largest value of 0.4. The effect of modulating unfolded-state interactions is studied by measuring ϕ-values in second- site mutants and under solvent conditions that perturb unfolded-state energetics in a defined way. Neither of these alterations significantly affects the distribution of ϕ-values. The results, combined with those of earlier studies that probe the role of hydrogen-bond formation in folding and the burial of surface area, reveal that the transition state for folding contains extensive backbone structure and buries a significant fraction of hydrophobic surface area, but lacks well developed side-chain-side-chain interactions. The folding transition state for NTL9 does not contain a specific "nucleus" consisting of a few key residues; rather, it involves extensive backbone hydrogen bonding and partially formed structure delocalized over almost the entire domain. The potential generality of these observations is discussed.
| S-EPMC5425357 | biostudies-literature
Project description:High molecular weight glutenin subunits (HMW-GS) are of a particular interest because of their biomechanical properties, which are important in many food systems such as breadmaking. Using fold-recognition techniques, we identified a fold compatible with the N-terminal domain of HMW-GS Dy10. This fold corresponds to the one adopted by proteins belonging to the cereal inhibitor family. Starting from three known protein structures of this family as templates, we built three models for the N-terminal domain of HMW-GS Dy10. We analyzed these models, and we propose a number of hypotheses regarding the N-terminal domain properties that can be tested experimentally. In particular, we discuss two possible ways of interaction between the N-terminal domains of the y-type HMW glutenin subunits. The first way consists in the creation of interchain disulfide bridges. According to our models, we propose two plausible scenarios: (1) the existence of an intrachain disulfide bridge between cysteines 22 and 44, leaving the three other cysteines free of engaging in intermolecular bonds; and (2) the creation of two intrachain disulfide bridges (involving cysteines 22-44 and cysteines 10-55), leaving a single cysteine (45) for creating an intermolecular disulfide bridge. We discuss these scenarios in relation to contradictory experimental results. The second way, although less likely, is nevertheless worth considering. There might exist a possibility for the N-terminal domain of Dy10, Nt-Dy10, to create oligomers, because homologous cereal inhibitor proteins are known to exist as monomers, homodimers, and heterooligomers. We also discuss, in relation to the function of the cereal inhibitor proteins, the possibility that this N-terminal domain has retained similar inhibitory functions.
| S-EPMC2312395 | biostudies-literature
Project description:We have studied room-temperature structural and dynamic properties of the p53 DNA-binding domain in both DNA-bound and DNA-free states. A cumulative 55 ns of explicit solvent molecular dynamics simulations with the particle mesh Ewald treatment of electrostatics was performed. It was found that the mean structures in the production portions of the trajectories agree well with the crystal structure: backbone root-mean-square deviations are in the range of 1.6 and 2.0 A. In both simulations, noticeable backbone deviations from the crystal structure are observed only in loop L6, due to the lack of crystal packing in the simulations. More deviations are observed in the DNA-free simulation, apparently due to the absence of DNA. Computed backbone B-factor is also in qualitative agreement with the crystal structure. Interestingly, little backbone structural change is observed between the mean simulated DNA-bound and DNA-free structures. A notable difference is observed only at the DNA-binding interface. The correlation between native contacts and inactivation mechanisms of tumor mutations is also discussed. In the H2 region, tumor mutations at sites D281, R282, E285, and E286 may weaken five key interactions that stabilize H2, indicating that their inactivation mechanisms may be related to the loss of local structure around H2, which in turn may reduce the overall stability to a measurable amount. In the L2 region, tumor mutations at sites Y163, K164, E171, V173, L194, R249, I251, and E271 are likely to be responsible for the loss of stability in the protein. In addition to apparent DNA contacts that are related to DNA binding, interactions R175/S183, S183/R196, and E198/N235 are highly occupied only in the DNA-bound form, indicating that they are more likely to be responsible for DNA binding.
| S-EPMC2522240 | biostudies-literature
Project description:There are multiple drugs for the treatment of type 2 diabetes, including traditional sulfonylureas biguanides, glinides, thiazolidinediones, α-glucosidase inhibitors, glucagon-like peptide-1 (GLP-1) receptor agonists, dipeptidyl peptidase IV (DPP-4) inhibitors, and sodium-glucose cotransporter 2 (SGLT2) inhibitors. α-Glucosidase inhibitors have been used to control postprandial glucose levels caused by type 2 diabetes since 1990. α-Glucosidases are rather crucial in the human metabolic system and are principally found in families 13 and 31. Maltase-glucoamylase (MGAM) belongs to glycoside hydrolase family 31. The main function of MGAM is to digest terminal starch products left after the enzymatic action of α-amylase; hence, MGAM becomes an efficient drug target for insulin resistance. In order to explore the conformational changes in the active pocket and unbinding pathway for NtMGAM, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed for two NtMGAM-inhibitor [de-O-sulfonated kotalanol (DSK) and acarbose] complexes. MD simulations indicated that DSK bound to NtMGAM may influence two domains (inserted loop 1 and inserted loop 2) by interfering with the spiralization of residue 497-499. The flexibility of inserted loop 1 and inserted loop 2 can influence the volume of the active pocket of NtMGAM, which can affect the binding progress for DSK to NtMGAM. ASMD simulations showed that compared to acarbose, DSK escaped from NtMGAM easily with lower energy. Asp542 is an important residue on the bottleneck of the active pocket of NtMGAM and could generate hydrogen bonds with DSK continuously. Our theoretical results may provide some useful clues for designing new α-glucosidase inhibitors to treat type 2 diabetes.
| S-EPMC8435576 | biostudies-literature
Project description:The phosphorylated kinase-inducible activation domain (pKID) adopts a helix-loop-helix structure upon binding to its partner KIX, although it is unstructured in the unbound state. The N-terminal and C-terminal regions of pKID, which adopt helices in the complex, are called, respectively, ?A and ?B. We performed all-atom multicanonical molecular dynamics simulations of pKID with and without KIX in explicit solvents to generate conformational ensembles. Although the unbound pKID was disordered overall, ?A and ?B exhibited a nascent helix propensity; the propensity of ?A was stronger than that of ?B, which agrees with experimental results. In the bound state, the free-energy landscape of ?B involved two low free-energy fractions: native-like and non-native fractions. This result suggests that ?B folds according to the induced-fit mechanism. The ?B-helix direction was well aligned as in the NMR complex structure, although the ?A helix exhibited high flexibility. These results also agree quantitatively with experimental observations. We have detected that the ?B helix can bind to another site of KIX, to which another protein MLL also binds with the adopting helix. Consequently, MLL can facilitate pKID binding to the pKID-binding site by blocking the MLL-binding site. This also supports experimentally obtained results.
| S-EPMC4030872 | biostudies-literature
Project description:STAR proteins are evolutionary conserved mRNA-binding proteins that post-transcriptionally regulate gene expression at all stages of RNA metabolism. These proteins possess conserved STAR domain that recognizes identical RNA regulatory elements as YUAAY. Recently reported crystal structures show that STAR domain is composed of N-terminal QUA1, K-homology domain (KH) and C-terminal QUA2, and mRNA binding is mediated by KH-QUA2 domain. Here, we present simulation studies done to investigate binding of mRNA to STAR protein, mammalian Quaking protein (QKI). We carried out conventional MD simulations of STAR domain in presence and absence of mRNA, and studied the impact of mRNA on the stability, dynamics and underlying allosteric mechanism of STAR domain. Our unbiased simulations results show that presence of mRNA stabilizes the overall STAR domain by reducing the structural deviations, correlating the 'within-domain' motions, and maintaining the native contacts information. Absence of mRNA not only influenced the essential modes of motion of STAR domain, but also affected the connectivity of networks within STAR domain. We further explored the dissociation of mRNA from STAR domain using umbrella sampling simulations, and the results suggest that mRNA binding to STAR domain occurs in multi-step: first conformational selection of mRNA backbone conformations, followed by induced fit mechanism as nucleobases interact with STAR domain.
| S-EPMC5626755 | biostudies-literature
Project description:The exploration of intrinsically disordered proteins in isolation is a crucial step to understand their complex dynamical behavior. In particular, the emergence of partially ordered states has not been explored in depth. The experimental characterization of such partially ordered states remains elusive due to their transient nature. Molecular dynamics mitigates this limitation thanks to its capability to explore biologically relevant timescales while retaining atomistic resolution. Here, millisecond unbiased molecular dynamics simulations were performed in the exemplar N-terminal region of p53. In combination with state-of-the-art Markov state models, simulations revealed the existence of several partially ordered states accounting for [Formula: see text] 40% of the equilibrium population. Some of the most relevant states feature helical conformations similar to the bound structure of p53 to Mdm2, as well as novel [Formula: see text]-sheet elements. This highlights the potential complexity underlying the energy surface of intrinsically disordered proteins.
| S-EPMC7382488 | biostudies-literature
Project description:Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native NMR or X-ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second-order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein-folding studies.
| S-EPMC4487363 | biostudies-literature
Project description:Interest in the structural and dynamic properties of unfolded proteins has increased in recent years owing to continued interest in protein folding and misfolding. Knowledge of the unfolded state under native conditions is particularly important for obtaining a complete picture of the protein folding process. The C-terminal domain of protein L9 is a globular alpha, beta protein with an unusual mixed parallel and antiparallel beta-strand structure. The folding kinetics and equilibrium unfolding of CTL9 strongly depend on pH, and follow a simple two state model. Both the native and the unfolded state can be significantly populated at pH 3.8 in the absence of denaturant, allowing the native state and the unfolded state to be characterized under identical conditions. Backbone (15)N, (13)C, (1)H and side-chain (13)C(beta), (1)H(beta) chemical shifts, amide proton NOEs, and (15)N R(2) relaxation rates were obtained for the two conformational states at pH 3.8. All the data indicate that the pH 3.8 native state is well folded and is similar to the native state at neutral pH. There is significant residual structure in the pH 3.8 unfolded state. The regions corresponding to the two native state alpha-helices show strong preference to populate helical phi and psi angles. The segment that connects alpha-helix 2 and beta-strand 2 has a significant tendency to form non-native alpha-helical structure. Comparison with the pH 2.0 unfolded state and the urea unfolded state indicates that the tendency to adopt both native and non-native helical structure is stronger at pH 3.8, demonstrating that the unfolded state of CTL9 under native-like conditions is more structured. The implications for the folding of CTL9 are discussed.
| S-EPMC3319011 | biostudies-literature