Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica��Exchange Molecular Dynamics Simulations

The folding process of the 20 residue Trp-cage mini-protein was investigated using standard temperature replica exchange molecular dynamics (T-RexMD) simulation and a biasing potential RexMD (BP-RexMD) method. In contrast to several conventional molecular dynamics simulations, both RexMD methods sampled conformations close to the native structure after 10–20 ns simulation time as the dominant conformational states. In contrast, to T-RexMD involving 16 replicas the BP-RexMD method achieved very similar sampling results with only five replicas. The result indicates that the BP-RexMD method is well suited to study folding processes of proteins at a significantly smaller computational cost, compared to T-RexMD. Both RexMD methods sampled not only similar final states but also agreed on the sampling of intermediate conformations during Trp-cage folding. The analysis of the sampled potential energy contributions indicated that Trp-cage folding is favored by both van der Waals and to a lesser degree electrostatic contributions. Folding does not introduce any significant sterical strain as reflected by similar energy distributions of bonded energy terms (bond length, bond angle and dihedral angle) of folded and unfolded Trp-cage structures.     

From Side Chains Rattling on Picoseconds to Ensemble Simulations of Protein Folding

Simulations of biological macromolecules have evolved tremendously since the discoveries of the 1970s. The field has moved from simple simulations in vacuo on picosecond scales to milliseconds of accurate sampling of large proteins, and it has become a standard tool in biochemistry and biophysics, rather than a dedicated theoretical one. This is partly due to increasing computational power, but it would not have been possible without huge research efforts invested in new algorithms and software. Here, we illustrate some of this development, both past and future challenges, and in particular, discuss how the recent introduction of modern ensemble methods is breaking the trend of ever-longer simulations to instead focus on throughput and sampling. This has not only helped simulations become much more accurate, but it provides statistical error estimates, which are critical, as simulations are increasingly used to predict properties that have not yet been measured experimentally.     

Combining Coarse-Grained Protein Models with Replica-Exchange All-Atom Molecular Dynamics

We describe a combination of all-atom simulations with CABS, a well-established coarse-grained protein modeling tool, into a single multiscale protocol. The simulation method has been tested on the C-terminal beta hairpin of protein G, a model system of protein folding. After reconstructing atomistic details, conformations derived from the CABS simulation were subjected to replica-exchange molecular dynamics simulations with OPLS-AA and AMBER99sb force fields in explicit solvent. Such a combination accelerates system convergence several times in comparison with all-atom simulations starting from the extended chain conformation, demonstrated by the analysis of melting curves, the number of native-like conformations as a function of time and secondary structure propagation. The results strongly suggest that the proposed multiscale method could be an efficient and accurate tool for high-resolution studies of protein folding dynamics in larger systems.     

A Hamiltonian Replica Exchange Molecular Dynamics (MD) Method for the Study of Folding,Based on the Analysis of the Stabilization Determinants of Proteins

Herein, we present a novel Hamiltonian replica exchange protocol for classical molecular dynamics simulations of protein folding/unfolding. The scheme starts from the analysis of the energy-networks responsible for the stabilization of the folded conformation, by means of the energy-decomposition approach. In this framework, the compact energetic map of the native state is generated by a preliminary short molecular dynamics (MD) simulation of the protein in explicit solvent. This map is simplified by means of an eigenvalue decomposition. The highest components of the eigenvector associated with the lowest eigenvalue indicate which sites, named “hot spots”, are likely to be responsible for the stability and correct folding of the protein. In the Hamiltonian replica exchange protocol, we use modified force-field parameters to treat the interparticle non-bonded potentials of the hot spots within the protein and between protein and solvent atoms, leaving unperturbed those relative to all other residues, as well as solvent-solvent interactions. We show that it is possible to reversibly simulate the folding/unfolding behavior of two test proteins, namely Villin HeadPiece HP35 (35 residues) and Protein A (62 residues), using a limited number of replicas. We next discuss possible implications for the study of folding mechanisms via all atom simulations.     

指数-6流体的自扩散系数的分子动力学模拟

Self-diffusion coefficients of exponential-six fluids are studied using equilibrium molecular dynamics simulation technique. Mean-square displacements and velocity autocorrelation functions are used to calculate self-diffusion coefficients through Einstein equation and Green-Kubo formula. It has been found that simulation results are in good agreement with experimental data for liquid argon which is taken as exponential-six fluid. The effects of density, temperature and steepness factor for repulsive part of exponential-six potential on self-diffusion coefficients are also investigated. The simulation results indicate that the self-diffusion coefficient of exponential-six fluid increases as temperature increases and density decreases. In addition, the larger self-diffusion coefficients are obtained as the steepness factor increases at the same temperature and density condition.     

Folding and Intrinsic Disorder of the Receptor Tyrosine Kinase KIT Insert Domain Seen by Conventional Molecular Dynamics Simulations

The kinase insert domain (KID) of RTK KIT is the key recruitment region for downstream signalling proteins. KID, studied by molecular dynamics simulations as a cleaved polypeptide and as a native domain fused to KIT, showed intrinsic disorder represented by a set of heterogeneous conformations. The accurate atomistic models showed that the helical fold of KID is mainly sequence dependent. However, the reduced fold of the native KID suggests that its folding is allosterically controlled by the kinase domain. The tertiary structure of KID represents a compact array of highly variable α- and 3₁₀-helices linked by flexible loops playing a principal role in the conformational diversity. The helically folded KID retains a collapsed globule-like shape due to non-covalent interactions associated in a ternary hydrophobic core. The free energy landscapes constructed from first principles—the size, the measure of the average distance between the conformations, the amount of helices and the solvent-accessible surface area—describe the KID disorder through a collection of minima (wells), providing a direct evaluation of conformational ensembles. We found that the cleaved KID simulated with restricted N- and C-ends better reproduces the native KID than the isolated polypeptide. We suggest that a cyclic, generic KID would be best suited for future studies of KID f post-transduction effects.     

Exploring the Protein Folding Dynamics of Beta3s with Two-Dimensional Ultraviolet (2DUV) Spectroscopy

Ultraviolet (UV) spectra of proteins originate from electronic excitations of their backbone chromophore and aromatic side chains and provide a sensitive probe of the secondary structure. Recently developed femtosecond lasers allow multidimensional spectroscopy to be extended into the UV regime. Two-dimensional UV (2DUV) techniques, with short pulses, provide a promising tool to study the structures and dynamics of proteins. We combined 2DUV spectroscopy and molecular dynamics generated free energy profiles to simulate the protein electronic transitions and UV photon echo signals to monitor the protein folding process of the small protein Beta3s. Near-ultraviolet (NUV) and far-ultraviolet (FUV) signals illustrate the variation of the 2D correlation plots when the protein evolves along the underlying free energy landscape. Chiral polarization configurations of the NUV and FUV pulses are sensitive to protein structural evolution. This work provides a protocol for applying multidimensional UV spectroscopy to study protein folding.     

Comparative Molecular Dynamics Simulations of Mitogen-Activated Protein Kinase-Activated Protein Kinase 5

The mitogen-activated protein kinase-activated protein kinase MK5 is a substrate of the mitogen-activated protein kinases p38, ERK3 and ERK4. Cell culture and animal studies have demonstrated that MK5 is involved in tumour suppression and promotion, embryogenesis, anxiety, cell motility and cell cycle regulation. In the present study, homology models of MK5 were used for molecular dynamics (MD) simulations of: (1) MK5 alone; (2) MK5 in complex with an inhibitor; and (3) MK5 in complex with the interaction partner p38α. The calculations showed that the inhibitor occupied the active site and disrupted the intramolecular network of amino acids. However, intramolecular interactions consistent with an inactive protein kinase fold were not formed. MD with p38α showed that not only the p38 docking region, but also amino acids in the activation segment, αH helix, P-loop, regulatory phosphorylation region and the C-terminal of MK5 may be involved in forming a very stable MK5-p38α complex, and that p38α binding decreases the residual fluctuation of the MK5 model. Electrostatic Potential Surface (EPS) calculations of MK5 and p38α showed that electrostatic interactions are important for recognition and binding.     

Study on the Characteristics of Gas Molecular Mean Free Path in Nanopores by Molecular Dynamics Simulations

This paper presents studies on the characteristics of gas molecular mean free path in nanopores by molecular dynamics simulation. Our study results indicate that the mean free path of all molecules in nanopores depend on both the radius of the nanopore and the gas-solid interaction strength. Besides mean free path of all molecules in the nanopore, this paper highlights the gas molecular mean free path at different positions of the nanopore and the anisotropy of the gas molecular mean free path at nanopores. The molecular mean free path varies with the molecule’s distance from the center of the nanopore. The least value of the mean free path occurs at the wall surface of the nanopore. The present paper found that the gas molecular mean free path is anisotropic when gas is confined in nanopores. The radial gas molecular mean free path is much smaller than the mean free path including all molecular collisions occuring in three directions. Our study results also indicate that when gas is confined in nanopores the gas molecule number density does not affect the gas molecular mean free path in the same way as it does for the gas in unbounded space. These study results may bring new insights into understanding the gas flow’s characteristic at nanoscale.     

蛋白质折叠与淀粉样沉积的产生

蛋白质是生物体内一切功能的执行者,而蛋白质正确折叠则是发挥其生理功能的一个重要前提条件。近年来研究发现,蛋白质的错误折叠可以形成淀粉样沉积进而导致一些疾病的产生。探讨了蛋白质折叠和淀粉样纤维产生的机制。     

Significance of Molecular Dynamics Simulations for Life Sciences

This article illustrates by examples the limited acceptance by biologists of predictions made with molecular dynamics simulations of biomolecules. Its purpose is to increase the awareness of biologists of the contribution that simulations can make to our understanding of biomolecule function.     

Molecular Dynamics Simulations of Acylpeptide Hydrolase Bound to Chlorpyrifosmethyl Oxon and Dichlorvos

Acylpeptide hydrolases (APHs) catalyze the removal of N-acylated amino acids from blocked peptides. Like other prolyloligopeptidase (POP) family members, APHs are believed to be important targets for drug design. To date, the binding pose of organophosphorus (OP) compounds of APH, as well as the different OP compounds binding and inducing conformational changes in two domains, namely, α/β hydrolase and β-propeller, remain poorly understood. We report a computational study of APH bound to chlorpyrifosmethyl oxon and dichlorvos. In our docking study, Val471 and Gly368 are important residues for chlorpyrifosmethyl oxon and dichlorvos binding. Molecular dynamics simulations were also performed to explore the conformational changes between the chlorpyrifosmethyl oxon and dichlorvos bound to APH, which indicated that the structural feature of chlorpyrifosmethyl oxon binding in APH permitted partial opening of the β-propeller fold and allowed the chlorpyrifosmethyl oxon to easily enter the catalytic site. These results may facilitate the design of APH-targeting drugs with improved efficacy.     

A Generic Force Field for Protein Coarse-Grained Molecular Dynamics Simulation

Coarse-grained (CG) force fields have become promising tools for studies of protein behavior, but the balance of speed and accuracy is still a challenge in the research of protein coarse graining methodology. In this work, 20 CG beads have been designed based on the structures of amino acid residues, with which an amino acid can be represented by one or two beads, and a CG solvent model with five water molecules was adopted to ensure the consistence with the protein CG beads. The internal interactions in protein were classified according to the types of the interacting CG beads, and adequate potential functions were chosen and systematically parameterized to fit the energy distributions. The proposed CG force field has been tested on eight proteins, and each protein was simulated for 1000 ns. Even without any extra structure knowledge of the simulated proteins, the Cα root mean square deviations (RMSDs) with respect to their experimental structures are close to those of relatively short time all atom molecular dynamics simulations. However, our coarse grained force field will require further refinement to improve agreement with and persistence of native-like structures. In addition, the root mean square fluctuations (RMSFs) relative to the average structures derived from the simulations show that the conformational fluctuations of the proteins can be sampled.     

Exploring the Cold-Adaptation Mechanism of Serine Hydroxymethyltransferase by Comparative Molecular Dynamics Simulations

Cold-adapted enzymes feature a lower thermostability and higher catalytic activity compared to their warm-active homologues, which are considered as a consequence of increased flexibility of their molecular structures. The complexity of the (thermo)stability-flexibility-activity relationship makes it difficult to define the strategies and formulate a general theory for enzyme cold adaptation. Here, the psychrophilic serine hydroxymethyltransferase (pSHMT) from Psychromonas ingrahamii and its mesophilic counterpart, mSHMT from Escherichia coli, were subjected to μs-scale multiple-replica molecular dynamics (MD) simulations to explore the cold-adaptation mechanism of the dimeric SHMT. The comparative analyses of MD trajectories reveal that pSHMT exhibits larger structural fluctuations and inter-monomer positional movements, a higher global flexibility, and considerably enhanced local flexibility involving the surface loops and active sites. The largest-amplitude motion mode of pSHMT describes the trends of inter-monomer dissociation and enlargement of the active-site cavity, whereas that of mSHMT characterizes the opposite trends. Based on the comparison of the calculated structural parameters and constructed free energy landscapes (FELs) between the two enzymes, we discuss in-depth the physicochemical principles underlying the stability-flexibility-activity relationships and conclude that (i) pSHMT adopts the global-flexibility mechanism to adapt to the cold environment and, (ii) optimizing the protein-solvent interactions and loosening the inter-monomer association are the main strategies for pSHMT to enhance its flexibility.     

气体在水中的分子动力学模拟

采用分子动力学（ＭＤ）模拟的方法在常温及工业应用背景条件下对ＣＨ４、ＮＨ３、ＣＯ２、Ｏ２这些气体在水中的结构及扩散情形进行了研究。ＭＤ模拟可以为这些涉及到气体在水中的工业应用情形的机理提供分子水平的解释,同时ＭＤ模拟还可为一些不易实验测定扩散性质的体系提供工程初步设计和过程开发所需的数据。     

KEAP1 Cancer Mutants: A Large-Scale Molecular Dynamics Study of Protein Stability

We have performed 280 μs of unbiased molecular dynamics (MD) simulations to investigate the effects of 12 different cancer mutations on Kelch-like ECH-associated protein 1 (KEAP1) (G333C, G350S, G364C, G379D, R413L, R415G, A427V, G430C, R470C, R470H, R470S and G476R), one of the frequently mutated proteins in lung cancer. The aim was to provide structural insight into the effects of these mutants, including a new class of ANCHOR (additionally NRF2-complexed hypomorph) mutant variants. Our work provides additional insight into the structural dynamics of mutants that could not be analyzed experimentally, painting a more complete picture of their mutagenic effects. Notably, blade-wise analysis of the Kelch domain points to stability as a possible target of cancer in KEAP1. Interestingly, structural analysis of the R470C ANCHOR mutant, the most prevalent missense mutation in KEAP1, revealed no significant change in structural stability or NRF2 binding site dynamics, possibly indicating an covalent modification as this mutant’s mode of action.     

Studies of Conformational Changes of Tubulin Induced by Interaction with Kinesin Using Atomistic Molecular Dynamics Simulations

The transition between strong and weak interactions of the kinesin head with the microtubule, which is regulated by the change of the nucleotide state of the head, is indispensable for the processive motion of the kinesin molecular motor on the microtubule. Here, using all-atom molecular dynamics simulations, the interactions between the kinesin head and tubulin are studied on the basis of the available high-resolution structural data. We found that the strong interaction can induce rapid large conformational changes of the tubulin, whereas the weak interaction cannot. Furthermore, we found that the large conformational changes of the tubulin have a significant effect on the interaction of the tubulin with the head in the weak-microtubule-binding ADP state. The calculated binding energy of the ADP-bound head to the tubulin with the large conformational changes is only about half that of the tubulin without the conformational changes.     

Molecular Dynamics Simulations of the Elusive Matrix-Open State of Mitochondrial ADP/ATP Carrier

We studied and compared in detail an elusive matrix open (m-state) and cytoplasmic open (c-state) state of ADP/ATP carrier (AAC) protein embedded in the DOPC bilayer by microsecond molecular dynamics (MD) simulations. We analyzed both states with and without cardiolipin (CDL) molecules, with a special emphasis on the recently obtained crystallographic structure of the AAC m-state. The obtained results show that both states of the protein are stable in the DOPC bilayer and impermeable to water. In comparison with the c-state of AAC, the m-state is more dynamic, but at the same time possesses a larger occluded area in the protein cavity. Both states of the protein are less flexible in simulations when CDL molecules are present, which is especially visible for the m-state. Finally, the analysis of the protein conformational changes during MD simulations shows that protein parts at the protein/lipid boundaries are prone to larger conformational changes in contrast to central region of the protein embedded in the bilayer core, thus further supporting the cycling mechanism suggested for ADP/ATP exchange by AAC.