Learning probabilistic models of hydrogen bond stability from molecular dynamics simulation trajectories
MetadataShow full item record
AbstractBackground: Hydrogen bonds (H-bonds) play a key role in both the formation and stabilization of protein structures. They form and break while a protein deforms, for instance during the transition from a non-functional to a functional state. The intrinsic strength of an individual H-bond has been studied from an energetic viewpoint, but energy alone may not be a very good predictor.Methods: This paper describes inductive learning methods to train protein-independent probabilistic models of H-bond stability from molecular dynamics (MD) simulation trajectories of various proteins. The training data contains 32 input attributes (predictors) that describe an H-bond and its local environment in a conformation c and the output attribute is the probability that the H-bond will be present in an arbitrary conformation of this protein achievable from c within a time duration ?. We model dependence of the output variable on the predictors by a regression tree.Results: Several models are built using 6 MD simulation trajectories containing over 4000 distinct H-bonds (millions of occurrences). Experimental results demonstrate that such models can predict H-bond stability quite well. They perform roughly 20% better than models based on H-bond energy alone. In addition, they can accurately identify a large fraction of the least stable H-bonds in a conformation. In most tests, about 80% of the 10% H-bonds predicted as the least stable are actually among the 10% truly least stable. The important attributes identified during the tree construction are consistent with previous findings.Conclusions: We use inductive learning methods to build protein-independent probabilistic models to study H-bond stability, and demonstrate that the models perform better than H-bond energy alone. 2011 Chikalov et al; licensee BioMed Central Ltd.
CitationChikalov I, Yao P, Moshkov M, Latombe J-C (2011) Learning probabilistic models of hydrogen bond stability from molecular dynamics simulation trajectories. BMC Bioinformatics 12: S34. doi:10.1186/1471-2105-12-S1-S34.
PubMed Central IDPMC3044290
The following license files are associated with this item:
Except where otherwise noted, this item's license is described as This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
- Statistical and molecular dynamics studies of buried waters in globular proteins.
- Authors: Park S, Saven JG
- Issue date: 2005 Aug 15
- Intra-protein hydrogen bonding is dynamically stabilized by electronic polarization.
- Authors: Duan LL, Mei Y, Zhang QG, Zhang JZ
- Issue date: 2009 Mar 21
- Insights on hydrogen-bond lifetimes in liquid and supercooled water.
- Authors: Martiniano HF, Galamba N
- Issue date: 2013 Dec 19
- Molecular dynamics of hydrogen bonds in protein-D2O: the solvent isotope effect.
- Authors: Sheu SY, Schlag EW, Selzle HL, Yang DY
- Issue date: 2008 Feb 7
- Dynamics of noncovalent interactions in all-α and all-β class proteins: implications for the stability of amyloid aggregates.
- Authors: Jain A, Sankararamakrishnan R
- Issue date: 2011 Dec 27
Showing items related by title, author, creator and subject.
All-Atom Molecular Dynamics Simulation of Protein Translocation through an α-Hemolysin NanoporeDi Marino, Daniele; Bonome, Emma Letizia; Tramontano, Anna; Chinappi, Mauro (American Chemical Society (ACS), 2015-08-06)© 2015 American Chemical Society. Nanopore sensing is attracting the attention of a large and varied scientific community. One of the main issues in nanopore sensing is how to associate the measured current signals to specific features of the molecule under investigation. This is particularly relevant when the translocating molecule is a protein and the pore is sufficiently narrow to necessarily involve unfolding of the translocating protein. Recent experimental results characterized the cotranslocational unfolding of Thioredoxin (Trx) passing through an α-hemolisin pore, providing evidence for the existence of a multistep process. In this study we report the results of all-atom molecular dynamics simulations of the same system. Our data indicate that Trx translocation involves two main barriers. The first one is an unfolding barrier associated with a translocation intermediate where the N-terminal region of Trx is stuck at the pore entrance in a conformation that strongly resembles the native one. After the abrupt unfolding of the N-terminal region, the Trx enters the α-hemolisin vestibule. During this stage, the constriction is occupied not only by the translocating residue but also by a hairpin-like structure forming a tangle in the constriction. The second barrier is associated with the disentangling of this region.
Recruitment of EB1, a Master Regulator of Microtubule Dynamics, to the Surface of the Theileria annulata SchizontWoods, Kerry L.; Theiler, Romina; Mühlemann, Marcus; Segiser, Adrian; Huber, Sandra; Ansari, Hifzur R.; Pain, Arnab; Dobbelaere, Dirk A. E. (Public Library of Science (PLoS), 2013-05-09)The apicomplexan parasite Theileria annulata transforms infected host cells, inducing uncontrolled proliferation and clonal expansion of the parasitized cell population. Shortly after sporozoite entry into the target cell, the surrounding host cell membrane is dissolved and an array of host cell microtubules (MTs) surrounds the parasite, which develops into the transforming schizont. The latter does not egress to invade and transform other cells. Instead, it remains tethered to host cell MTs and, during mitosis and cytokinesis, engages the cell's astral and central spindle MTs to secure its distribution between the two daughter cells. The molecular mechanism by which the schizont recruits and stabilizes host cell MTs is not known. MT minus ends are mostly anchored in the MT organizing center, while the plus ends explore the cellular space, switching constantly between phases of growth and shrinkage (called dynamic instability). Assuming the plus ends of growing MTs provide the first point of contact with the parasite, we focused on the complex protein machinery associated with these structures. We now report how the schizont recruits end-binding protein 1 (EB1), a central component of the MT plus end protein interaction network and key regulator of host cell MT dynamics. Using a range of in vitro experiments, we demonstrate that T. annulata p104, a polymorphic antigen expressed on the schizont surface, functions as a genuine EB1-binding protein and can recruit EB1 in the absence of any other parasite proteins. Binding strictly depends on a consensus SxIP motif located in a highly disordered C-terminal region of p104. We further show that parasite interaction with host cell EB1 is cell cycle regulated. This is the first description of a pathogen-encoded protein to interact with EB1 via a bona-fide SxIP motif. Our findings provide important new insight into the mode of interaction between Theileria and the host cell cytoskeleton. 2013 Woods et al.
Dynamic regulation of genome-wide pre-mRNA splicing and stress tolerance by the Sm-like protein LSm5 in ArabidopsisCui, Peng; Zhang, ShouDong; Ding, Feng; Ali, Shahjahan; Xiong, Liming (Springer Nature, 2014-01-07)Background: Sm-like proteins are highly conserved proteins that form the core of the U6 ribonucleoprotein and function in several mRNA metabolism processes, including pre-mRNA splicing. Despite their wide occurrence in all eukaryotes, little is known about the roles of Sm-like proteins in the regulation of splicing.Results: Here, through comprehensive transcriptome analyses, we demonstrate that depletion of the Arabidopsis supersensitive to abscisic acid and drought 1 gene (SAD1), which encodes Sm-like protein 5 (LSm5), promotes an inaccurate selection of splice sites that leads to a genome-wide increase in alternative splicing. In contrast, overexpression of SAD1 strengthens the precision of splice-site recognition and globally inhibits alternative splicing. Further, SAD1 modulates the splicing of stress-responsive genes, particularly under salt-stress conditions. Finally, we find that overexpression of SAD1 in Arabidopsis improves salt tolerance in transgenic plants, which correlates with an increase in splicing accuracy and efficiency for stress-responsive genes.Conclusions: We conclude that SAD1 dynamically controls splicing efficiency and splice-site recognition in Arabidopsis, and propose that this may contribute to SAD1-mediated stress tolerance through the metabolism of transcripts expressed from stress-responsive genes. Our study not only provides novel insights into the function of Sm-like proteins in splicing, but also uncovers new means to improve splicing efficiency and to enhance stress tolerance in a higher eukaryote. 2014 Cui et al.; licensee BioMed Central Ltd.