scholarly journals Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

2018 ◽  
Author(s):  
Yong Wang ◽  
Pengfei Tian ◽  
Wouter Boomsma ◽  
Kresten Lindorff-Larsen

AbstractEnergy landscape theory suggests that native interactions are a major determinant of the folding mechanism of a protein. Thus, structure-based (Gō) models have, aided by coarse-graining techniques, shown great success in capturing the mechanisms of protein folding and conformational changes. In certain cases, however, non-native interactions and atomic details are also essential to describe the protein dynamics, prompting the development of a variety of structure-based models which include non-native interactions, and differentiate between different types of attractive potentials. Here, we describe an all-protein-atom hybrid model, termed ProfasiGo, that integrates an implicit solvent all-atom physics-based model (called Profasi) and a structure-based Gō potential, and its implementation in two software packages (PHAISTOS and ProFASi) that are developed for Monte Carlo sampling of protein molecules. We apply the ProfasiGo model to study the folding free energy landscapes of four topologically similar proteins, one of which can be folded by the simplified potential Profasi, and two that have been folded by explicit solvent, all-atom molecular dynamics simulations with the CHARMM22∗ force field. Our results reveal that the hybrid ProfasiGo model is able to capture many of the details present in the physics-based potentials, while retaining the advantages of Gō models for sampling and guiding to the native state. We expect that the model will be widely applicable to study the folding of more complex proteins, or to study conformational dynamics and integration with experimental data.

2018 ◽  
Vol 122 (49) ◽  
pp. 11174-11185 ◽  
Author(s):  
Yong Wang ◽  
Pengfei Tian ◽  
Wouter Boomsma ◽  
Kresten Lindorff-Larsen

Author(s):  
Balaji Selvam ◽  
Ya-Chi Yu ◽  
Liqing Chen ◽  
Diwakar Shukla

<p>The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane. However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as 1 for the glucose the transport mechanism. SWEETs undergoes structural transition to outward-facing (OF), Occluded (OC) and inward-facing (IF) and strongly support alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly same for OF, OC and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and act as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.</p>


2008 ◽  
Vol 2 ◽  
pp. BBI.S460 ◽  
Author(s):  
Lee-Wei Yang ◽  
Choon-Peng Chng

In this review, we summarize the progress on coarse-grained elastic network models (CG-ENMs) in the past decade. Theories were formulated to allow study of conformational dynamics in time/space frames of biological interest. Several highlighted models and their underlined hypotheses are introduced in physical depth. Important ENM offshoots, motivated to reproduce experimental data as well as to address the slow-mode-encoded configurational transitions, are also introduced. With the theoretical developments, computational cost is significantly reduced due to simplified potentials and coarse-grained schemes. Accumulating wealth of data suggest that ENMs agree equally well with experiment in describing equilibrium dynamics despite their distinct potentials and levels of coarse-graining. They however do differ in the slowest motional components that are essential to address large conformational changes of functional significance. The difference stems from the dissimilar curvatures of the harmonic energy wells described for each model. We also provide our views on the predictability of ‘open to close’ (open→close) transitions of biomolecules on the basis of conformational selection theory. Lastly, we address the limitations of the ENM formalism which are partially alleviated by the complementary CG-MD approach, to be introduced in the second paper of this two-part series.


2019 ◽  
Author(s):  
Balaji Selvam ◽  
Ya-Chi Yu ◽  
Liqing Chen ◽  
Diwakar Shukla

<p>The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane. However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as 1 for the glucose the transport mechanism. SWEETs undergoes structural transition to outward-facing (OF), Occluded (OC) and inward-facing (IF) and strongly support alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly same for OF, OC and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and act as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.</p>


2021 ◽  
Author(s):  
Yulian Gavrilov ◽  
Felix Kümmerer ◽  
Simone Orioli ◽  
Andreas Prestel ◽  
Kresten Lindorff-Larsen ◽  
...  

The conformational heterogeneity of a folded protein can affect both its function but also stability and folding. We recently discovered and characterized a stabilized double mutant (L49I/I57V) of the protein CI2 and showed that state-of-the-art prediction methods could not predict the increased stability relative to the wild-type protein. Here we have examined whether changed native state dynamics, and resulting entropy changes, can explain the stability changes in the double mutant protein, as well as the two single mutant forms. We have combined NMR relaxation measurements of the ps-ns dynamics of amide groups in the backbone and the methyl groups in the side-chains with molecular dynamics simulations to quantify the native state dynamics. The NMR experiments reveal that the mutations have different effects on the conformational flexibility of CI2: A reduction in conformational dynamics (and entropy) of the native state of L49I variant correlates with its decreased stability, while increased dynamics of the I57V and L49I/I57V variants correlates with their increased stability. These findings suggest that explicitly accounting for changes in native state entropy might be needed to improve the predictions of the effect of mutations on protein stability.


2019 ◽  
Vol 6 (5) ◽  
pp. 190219 ◽  
Author(s):  
An Luo ◽  
Xinbo Li ◽  
Xuecheng Zhang ◽  
Huadong Zhan ◽  
Hewei Du ◽  
...  

Heat-shock protein of 90 kDa (Hsp90) is a key molecular chaperone involved in folding the synthesized protein and controlling protein quality. Conformational dynamics coupled to ATPase activity in N-terminal domain is essential for Hsp90's function. However, the relevant process is still largely unknown in plant Hsp90s, especially those required for plant embryogenesis which is inextricably tied up with human survival. Here, AtHsp90.6, a member of Hsp90 family in Arabidopsis , was firstly identified as a protein essential for embryogenesis. Thus we modelled AtHsp90.6 in its functionally closed ‘lid-down’ and open ‘lid-up’ states, exploring the nucleotide binding mechanism in these two states. Free energy landscape and electrostatic potential analysis revealed the switching mechanism between these two states. Collectively, this study quantitatively analysed the conformational changes of AtHsp90.6 bound to ATP or ADP. This result may help us understand the mechanism of action of AtHsp90.6 in future.


2018 ◽  
Author(s):  
D. R. Kattnig ◽  
C. Nielsen ◽  
I. A. Solov’yov

AbstractBirds appear to be equipped with a light-dependent, radical-pair-based magnetic compass that relies on truly quantum processes. While the identity of the sensory protein has remained speculative, cryptochrome 4 has recently been identified as the most auspicious candidate. Here, we report on allatom molecular dynamics (MD) simulations addressing the structural reorganisations that accompany the photoreduction of the flavin cofactor in a model of the European robin cryptochrome 4 (ErCry4). Extensive MD simulations reveal that the photo-activation of ErCry4 induces large-scale conformational changes on short (hundreds of nanoseconds) timescales. Specifically, the photo-reduction is accompanied with the release of the C-terminal tail, structural rearrangements in the vicinity of the FAD-binding site, and the noteworthy formation of an α-helical segment at the N-terminal part. Some of these rearrangements appear to expose potential phosphorylation sites. We describe the conformational dynamics of the protein using a graph-based approach that is informed by the adjacency of residues and the correlation of their local motions. This approach reveals densely coupled reorganisation communities, which facilitate an efficient signal transduction due to a high density of hubs. These communities are interconnected by a small number of highly important residues characterized by high betweenness centrality. The network approach clearly identifies the sites restructuring upon photoactivation, which appear as protrusions or delicate bridges in the reorganisation network. We also find that, unlike in the homologous cryptochrome from D. melanogaster, the release of the C-terminal domain does not appear to be correlated with the transposition of a histidine residue close to the FAD cofactor.


eLife ◽  
2017 ◽  
Vol 6 ◽  
Author(s):  
Ximena Steinberg ◽  
Marina A Kasimova ◽  
Deny Cabezas-Bratesco ◽  
Jason D Galpin ◽  
Ernesto Ladron-de-Guevara ◽  
...  

TRPV1 channels support the detection of noxious and nociceptive input. Currently available functional and structural data suggest that TRPV1 channels have two gates within their permeation pathway: one formed by a ′bundle-crossing′ at the intracellular entrance and a second constriction at the selectivity filter. To describe conformational changes associated with channel gating, the fluorescent non-canonical amino acid coumarin-tyrosine was genetically encoded at Y671, a residue proximal to the selectivity filter. Total internal reflection fluorescence microscopy was performed to image the conformational dynamics of the channels in live cells. Photon counts and optical fluctuations from coumarin encoded within TRPV1 tetramers correlates with channel activation by capsaicin, providing an optical marker of conformational dynamics at the selectivity filter. In agreement with the fluorescence data, molecular dynamics simulations display alternating solvent exposure of Y671 in the closed and open states. Overall, the data point to a dynamic selectivity filter that may serve as a gate for permeation.


Sign in / Sign up

Export Citation Format

Share Document