scholarly journals Coupled molecular dynamics mediate long- and short-range epistasis between mutations that affect stability and aggregation kinetics

2018 ◽  
Vol 115 (47) ◽  
pp. E11043-E11052 ◽  
Author(s):  
Haoran Yu ◽  
Paul A. Dalby

Multiple mutations are typically required to significantly improve protein stability or aggregation kinetics. However, when several substitutions are made in a single protein, the mutations can potentially interact in a nonadditive manner, resulting in epistatic effects, which can hamper protein-engineering strategies to improve thermostability or aggregation kinetics. Here, we have examined the role of protein dynamics in mediating epistasis between pairs of mutations. With Escherichia coli transketolase (TK) as a model, we explored the epistatic interactions between two single variants H192P and A282P, and also between the double-mutant H192P/A282P and two single variants, I365L or G506A. Epistasis was determined for several measures of protein stability, including the following: the free-energy barrier to kinetic inactivation, ∆∆G‡; thermal transition midpoint temperatures, Tm; and aggregation onset temperatures, Tagg. Nonadditive epistasis was observed between neighboring mutations as expected, but also for distant mutations located in the surface and core regions of different domains. Surprisingly, the epistatic behaviors for each measure of stability were often different for any given pairwise recombination, highlighting that kinetic and thermodynamic stabilities do not always depend on the same structural features. Molecular-dynamics simulations and a pairwise cross-correlation analysis revealed that mutations influence the dynamics of their local environment, but also in some cases the dynamics of regions distant in the structure. This effect was found to mediate epistatic interactions between distant mutations and could therefore be exploited in future protein-engineering strategies.

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>


2005 ◽  
Vol 502 ◽  
pp. 51-56 ◽  
Author(s):  
Sakir Erkoc

The structural and electronic properties of isolated neutral ZnmCdn clusters for m+n £ 3 have been investigated by performing density functional theory calculations at B3LYP level. The optimum geometries, vibrational frequencies, electronic structures, and the possible dissosiation channels of the clusters considered have been obtained. An empirical many-body potential energy function (PEF), which comprices two- and three-body atomic interactions, has been developed to investigate the structural features and energetics of ZnmCdn (m+n=3,4) microclusters. The most stable structures were found to be triangular for the three-atom clusters and tetrahedral for the four-atom clusters. On the other hand, the structural features and energetics of Znn-mCdm (n=7,8) microclusters, and Zn50, Cd50, Zn25Cd25, Zn12Cd38, and Zn38Cd12 nanoparticles have been investigated by performing molecular-dynamics computer simulations using the developed PEF. The most stable structures were found to be compact and three-dimensional for all elemental and mixed clusters. An interesting structural feature of the mixed clusters is that Zn and Cd atoms do not mix in mixed clusters, they come together almost without mixing. Surface and bulk properties of Zn, Cd, and ZnCd systems have been investigated too by performing molecular-dynamics simulations using the developed PEF. Surface reconstruction and multilayer relaxation on clean surfaces, adatom on surface, substitutional atom on surface and bulk materials, and vacancy on surface and bulk materials have been studied extensively.


2021 ◽  
Author(s):  
Hugo Arellano-Santoyo ◽  
Rogelio A Hernandez-Lopez ◽  
Emma Stokasimov ◽  
Ray YR Wang ◽  
David Pellman ◽  
...  

The microtubule (MT) cytoskeleton is central to cellular processes including axonal growth, intracellular transport, and cell division, all of which rely on precise spatiotemporal control of MT organization. Kinesin-8s play a key role in regulating MT length by combining highly processive directional motility with MT-end disassembly. However, how kinesin-8 switches between these two apparently opposing activities remains unclear. Here, we define the structural features underlying this molecular switch through cryo-EM analysis of the yeast kinesin-8, Kip3 bound to MTs, and molecular dynamics simulations to approximate the complex of Kip3 with the curved tubulin state found at the MT plus-end. By integrating biochemical and single-molecule biophysical assays, we identified specific intra- and intermolecular interactions that modulate processive motility and MT disassembly. Our findings suggest that Kip3 undergoes conformational changes in response to tubulin curvature that underlie its unique ability to interact differently with the MT lattice than with the MT-end.


2018 ◽  
Vol 20 (5) ◽  
pp. 3523-3530 ◽  
Author(s):  
Zhi Yue ◽  
Jana Shen

Constant pH molecular dynamics simulations of BBL reveals negligible folding free energy barrier that is pH dependent and a sparsely populated dry molten globule state.


Molecules ◽  
2019 ◽  
Vol 24 (18) ◽  
pp. 3261
Author(s):  
Sousa ◽  
Laurent ◽  
Quéméner ◽  
Mortier ◽  
Questel

Interleukin 15 (IL-15), a four-helix bundle cytokine, is involved in a plethora of different cellular functions and, particularly, plays a key role in the development and activation of immune responses. IL-15 forms receptor complexes by binding with IL-2Rβ- and common γ(γc)-signaling subunits, which are shared with other members of the cytokines family (IL-2 for IL-2Rβ- and all other γc- cytokines for γc). The specificity of IL-15 is brought by the non-signaling α-subunit, IL-15Rα. Here we present the results of molecular dynamics simulations carried out on four relevant forms of IL-15: its monomer, IL-15 interacting individually with IL-15Rα (IL-15/IL-15Rα), with IL-2Rβ/γc subunits (IL-15/IL-2Rβ/γc) or with its three receptors simultaneously (IL-15/IL-15Rα/IL-2Rβ/γc). Through the analyses of the various trajectories, new insights on the structural features of the interfaces are highlighted, according to the considered form. The comparison of the results with the experimental data, available from X-ray crystallography, allows, in particular, the rationalization of the importance of IL-15 key residues (e.g. Asp8, Lys10, Glu64). Furthermore, the pivotal role of water molecules in the stabilization of the various protein-protein interfaces and their H-bonds networks are underlined for each of the considered complexes.


Polymers ◽  
2020 ◽  
Vol 12 (12) ◽  
pp. 3067
Author(s):  
Arantxa Arbe ◽  
Fernando Alvarez ◽  
Juan Colmenero

Combining neutron scattering and fully atomistic molecular dynamics simulations allows unraveling structural and dynamical features of polymer melts at different length scales, mainly in the intermolecular and monomeric range. Here we present the methodology developed by us and the results of its application during the last years in a variety of polymers. This methodology is based on two pillars: (i) both techniques cover approximately the same length and time scales and (ii) the classical van Hove formalism allows easily calculating the magnitudes measured by neutron scattering from the simulated atomic trajectories. By direct comparison with experimental results, the simulated cell is validated. Thereafter, the information of the simulations can be exploited, calculating magnitudes that are experimentally inaccessible or extending the parameters range beyond the experimental capabilities. We show how detailed microscopic insight on structural features and dynamical processes of various kinds has been gained in polymeric systems with different degrees of complexity, and how intriguing questions as the collective behavior at intermediate length scales have been faced.


Sign in / Sign up

Export Citation Format

Share Document