Hydration-dependent protein dynamics revealed by molecular dynamics simulation of crystalline staphylococcal nuclease

Molecular dynamics simulations of crystalline Staphylococcal nuclease in full and minimal hydration states were performed to study hydration effects on protein dynamics at temperatures ranging from 100 to 300 K. In a full hydration state (hydration ratio in weight, h=0.49), gaps are fully filled with water molecules, whereas only crystal waters are included in a minimal hydration state (h=0.09). The inflection of the atomic mean-square fluctuation of protein as a function of temperature, known as the glass-like transition, is observed at approximately 220 K in both cases, which is more significant in the full hydration state. By examining the temperature dependence of residual fluctuation, we found that the increase of fluctuations in the loop and terminal regions, which are exposed to water, is much greater than that in other regions in the full hydration state, but the mobilities of the corresponding regions are relatively restricted in the minimal hydration state by intermolecular contact. The atomic mean-square fluctuation of water molecules in the full hydration state at 300 K is 1 order of magnitude greater than that in the minimal hydration state. Above the transition temperature, most water molecules in the full hydration state behave like bulk water and act as a lubricant for protein dynamics. In contrast, water molecules in the minimal hydration state tend to form more hydrogen bonds with the protein, restricting the fluctuation of these water molecules to the level of the protein. Thus, intermolecular interaction and solvent mobility are important to understand the glass-like transition in proteins.     

Diversity of potential hydrogen bonds in cellulose I revealed by molecular dynamics simulation

We have performed molecular dynamics calculations using a revised version of the Gromos56A_carbo force field to understand the consequences of the different potential hydrogen bonding patterns on the structural stability and thermal behavior of the I_α and I_β forms of native cellulose. For each allomorph, we considered three patterns of hydrogen bonds: two patterns obtained from neutron diffraction data refinement and a regular mixture of the two. Upon annealing, the hydrogen bonding schemes of cellulose I_β, irrespective of the starting structure, re-arranged into the main hydrogen bond pattern experimentally observed (pattern A). On the other hand, the I_α structures, irrespective of the starting hydrogen bonding pattern, converged to a non-experimental structure where the adjacent chains are shifted along the chain direction by 0.12 nm in the hydrogen-bonded plane, and the hydroxymethyl group conformation alternates between gt and tg along the chain. The exotic structure in I_α might be a consequence of a deficiency in force field parameters and/or potential molecular arrangement in less crystalline cellulose.     

Hydrophobic attraction as revealed by AFM force measurements and molecular dynamics simulation

Spherical calcium dioleate particles ( approximately 10 mum in diameter) were used as AFM (atomic force microscope) probes to measure interaction forces of the collector colloid with calcite and fluorite surfaces. The attractive AFM force between the calcium dioleate sphere and the fluorite surface is strong and has a longer range than the DLVO (Derjaguin-Landau-Verwey-Overbeek) prediction. The AFM force between the calcium dioleate sphere and the mineral surfaces does not agree with the DLVO prediction. Consideration of non-DLVO forces, including the attractive hydrophobic force and the repulsive hydration force, was necessary to explain the experimental results. The non-DLVO interactions considered were justified by the different interfacial water structures at calcite- and fluorite-water interfaces as revealed by the numerical computation experiments with molecular dynamics simulation.     

兔子朊蛋白结构稳定性的分子动力学研究(英文)

朊蛋白病是一种能够对人类和动物带来致命影响,并具有高度传染性的神经退行性疾病.兔子是目前已经报道的哺乳类动物中对朊蛋白病免疫的少数几个物种之一.我们将分子动力学和操控式分子动力学模拟相结合,研究了兔子正常朊蛋白的结构稳定性;同时讨论了蛋白结构的收敛性及刚性分布,并揭示了兔子朊蛋白中关键二级结构的动力学以及受力各向异性特征,证实了兔子朊蛋白结构的稳定性特征.     

A study on the structural transition of a single polymer chain by parallel tempering molecular dynamics simulation

The structural transition of a single polymer chain with chain length of 100,200 and 300 beads was investigated by parallel tempering MD simulation.Our simulation results can capture the structural change from random coil to orientationally ordered structure with decreasing temperature.The clear transition was observed on the curves of radius of gyration and global orientational order parameter P as the function of temperature,which demonstrated structural formation of a single polymer chain.The linear relationships between three components of square radius of gyration R_gx²,R_gx²,R_gz² and global orientational order P can be obtained under the structurally transformational process.The slope of the linear relationship between x（or y-axis） component R_gx²(or R_gy²) and P is negative,while that of RL as the function of P is positive.The absolute value of slope is proportional to the chain length.Once the single polymer chain takes the random coil or ordered configuration,the linear relationship is invalid.The conformational change was also analyzed on microscopic scale.The polymer chain can be treated as the construction of rigid stems connecting by flexible loops.The deviation from exponentially decreased behavior of stem length distribution becomes prominent,indicating a stiffening of the chain arises leading to more and more segments ending up in the trans state with decreasing temperature.The stem length N_tr is about 21 bonds indicating the polymer chain is ordered with the specific fold length.So,the simulation results,which show the prototype of a liquid-crystalline polymer chain,are helpful to understand the crystallization process of crystalline polymers.     

The mechanism of vesicle fusion as revealed by molecular dynamics simulations

We describe molecular dynamics simulations elucidating the molecular details of the process of fusion for small lipid vesicles. The simulations are based on a coarse grained (CG) lipid model that accurately represents the lamellar state of a variety of phospholipids and enables us to observe intermediate stages during fusion at near atomic detail. Simulations were conducted on a variety of systems containing common phospholipids such as phosphatidylcholine (PC), phosphatidylethanolamine (PE), lysoPC, and mixtures of the above. The fusion intermediates found are in general agreement with the stalk-pore mechanism. Transient pores sometimes form adjacent to the stalk, however, resulting in the mixing of lipids from the outer and inner monolayers. The speed of stalk formation and the opening of the fusion pore can be modulated by altering the lipid composition in qualitative agreement with experimental observations.     

Proton transport in a binary biomimetic solution revealed by molecular dynamics simulation

We report the simulation results of the proton transport in a binary mixture of amphiphilic tetramethylurea (TMU) molecules and water. We identify different mechanisms that either facilitate or retard the proton transport. The efficiency of these mechanisms depends on the TMU concentration. The overall picture is more complicated than a recent suggestion that the presence of amphiphilic molecules suppresses the proton mobility by slowing down the reorientation of the surrounding water molecules. It has also been suggested that the hydronium ion induces local water orientational order, which results in an ordered region that has to move along with the proton potentially slowing down the proton transport as suggested by experiment. We find that water-wire like structures formed at low amphiphile concentrations facilitate proton transfer, and reduction of the hydrogen bond connectivity induced at high concentrations retards it.     

Structure and dynamics of an unfolded protein examined by molecular dynamics simulation

The accurate characterization of the structure and dynamics of proteins in disordered states is a difficult problem at the frontier of structural biology whose solution promises to further our understanding of protein folding and intrinsically disordered proteins. Molecular dynamics (MD) simulations have added considerably to our understanding of folded proteins, but the accuracy with which the force fields used in such simulations can describe disordered proteins is unclear. In this work, using a modern force field, we performed a 200 μs unrestrained MD simulation of the acid-unfolded state of an experimentally well-characterized protein, ACBP, to explore the extent to which state-of-the-art simulation can describe the structural and dynamical features of a disordered protein. By comparing the simulation results with the results of NMR experiments, we demonstrate that the simulation successfully captures important aspects of both the local and global structure. Our simulation was ~2 orders of magnitude longer than those in previous studies of unfolded proteins, a length sufficient to observe repeated formation and breaking of helical structure, which we found to occur on a multimicrosecond time scale. We observed one structural feature that formed but did not break during the simulation, highlighting the difficulty in sampling disordered states. Overall, however, our simulation results are in reasonable agreement with the experimental data, demonstrating that MD simulations can already be useful in describing disordered proteins. Finally, our direct calculation of certain NMR observables from the simulation provides new insight into the general relationship between structural features of disordered proteins and experimental NMR relaxation properties.     

Bifurcated hydrogen bond in lithium nitrate trihydrate probed by ab initio molecular dynamics

The hydrogen-bond dynamics of lithium nitrate trihydrate has been studied by a combined approach based on ab initio molecular dynamics simulations and wavelet analysis. The simultaneous bifurcated interaction between one hydrogen atom of water molecules and two oxygen atoms of nitrate ions is the pivotal feature of the crystal structure: this bifurcated interaction has deep effects on the O-H stretching region of the vibrational spectrum. The structural, dynamic, spectroscopic, and electronic properties of the bifurcated hydrogen bond have been investigated computationally, elucidating at the molecular level the differences with weak and strong hydrogen bonds present in the crystal. These studies corroborate the very recent IR experiments performed on the lithium nitrate trihydrate crystal, offering new perspectives to interpreting the vibrational spectra. In fact, this approach allows obtaining two-dimensional plots, which summarize the essential features of both the hydrogen-bond network and IR spectra, resulting in a peculiar "signature" of the bifurcated interaction.     

Capturing coacervate formation and protein partition by molecular dynamics simulation

Biomolecules localize and function in microenvironments where their local concentration, spatial organization, and biochemical reactivity are regulated. To compartmentalize and control the local properties of the native microenvironment, cellular mimics and artificial bioreactors have been developed to study the properties of membraneless organelles or mimic the bio-environment for life origin. Here, we carried out molecular dynamics simulation with the Martini 3.0 model to reproduce the experimental salt concentration and pH dependency of different complex coacervates. We showed that coacervates inside vesicles are able to change their shape. In addition, we used these coacervate systems to explore the partitioning of the ubiquitous cytoskeletal protein actin and found that actin spontaneously partitions to all the coacervate peripheries. Therefore, we believe that our study can provide a better understanding of the versatile coacervate platform, where biomolecules partition and gather to fulfill their biological duties.

We carried out molecular dynamics simulation with the Martini 3.0 model to reproduce the experimental salt concentration and pH dependency of different complex coacervates. We showed that coacervates inside vesicles are able to change their shape.     

Dissimilar interaction of CB1/CB2 with lipid bilayers as revealed by molecular dynamics simulation

Cannabinoid receptors CB1 and CB2 are a striking class of transmembrane proteins involved in a high number of important biological processes. In spite of the inherent similarity (40% in aminoacid sequence) these receptors are found in different cell environments. In addition to this, CB1 activity has been intimately associated with lipid rafts whereas CB2 has not. In this work we have performed a 50 nanosecond molecular dynamics simulation of the inactive conformations of both receptors inserted in a POPC lipid bilayer. Although in both cases the overall protein structure is maintained along the entire simulation we have found important differences in the protein-lipid interaction. While CB1 tends to distort the lipid bilayer regularity, especially in the extracellular moiety, CB2 has a minor influence on the lipid distribution along the plane of the bilayer. This observation is consistent with some experimental facts observed in these cannabinoid receptors with regard to lipid/protein interaction.     

A molecular orbital explanation of bond distance variation caused by hydrogen bond formation

For hydrogen bond systems X–D–HA–Y, a simple molecular orbital model is proposed to understand the mechanism of the bond distance variations caused by the hydrogen bond formation. This model explains the bond distance variations for X–D and A–Y as follows. Electrostatic potential that the electrons in a molecule receive from other molecules causes the changes in atomic orbital energy differences between the bonded atoms. Then, the changes in the orbital energy differences make the bond orders larger or smaller and consequently the bond distances vary. The validity of this model has been confirmed by the effective fragment potential method, using the test systems of (HCOOH)₂, HCONH₂ (formamide) crystal and BF₃·2H₂O crystal.     

The effect of hydrogen bonding on oligoleucine structure in water: A molecular dynamics simulation study

Molecular dynamics simulations were conducted in order to improve our understanding of the forces that determine polyleucine chains conformations and govern polyleucine self-assembly in aqueous solutions. Simulations of 10 repeat unit oligoleucine in aqueous solution were performed using the optimized potential for liquid simulations (OPLS) - all atom force field using the canonical ensemble for a minimum of 1.3 ns. These simulations provided information on conformations, chain collapse and intermolecular aggregation. Simulations indicate that single isotactic oligoleucine chains in dilute solution assume tightly packed, regular hairpin conformations while atactic oligoleucine assumes a much less regular and less compact structure. The regular, compact collapsed isotactic chain exhibited a greater degree of intramolecular hydrogen bonding and an increased level of hydrophobic t-butyl functional group aggregation compared to the atactic chain. This occurs at the expense of reduced leucine-water hydrogen bonding.     

Secondary structure sensitivity of hydrogen bond lifetime dynamics in the protein hydration layer

The heterogeneous nature of a protein surface plays an essential role in its biological activity and molecular recognition, and this role is mediated at least partly through the surrounding water molecules. We have performed atomistic molecular dynamics simulations of an aqueous solution of HP-36 to investigate the correlation between the dynamics of the hydration layer water molecules and the lifetimes of protein-water hydrogen bonds. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the quasi-bound states in the surface and quantification of the dynamic equilibrium between quasi-bound and free water molecules in terms of time-dependent rate of interconversion. It is noticed that, irrespective of the structural heterogeneity of different segments of the protein, namely the three alpha-helices, the positively charged amino acid residues form longer-lived hydrogen bonds with water. The overall relaxation behavior of protein-water hydrogen bonds is found to differ significantly among the three helices of the protein. Study of water number density fluctuation reveals that the hydration layer of helix-3 is much less rigid, which can be correlated with faster structural relaxation of the hydrogen bonds between its residues and water. This also agrees excellently with faster translational and rotational motions of water near helix-3, and hence the lower rigidity of its hydration layer. The lower rigidity of the helix-3 hydration layer also correlates well with the biological activity of the protein, as several of the active-site residues of HP-36 are located in helix-3.     

Protein dynamics revealed by hydrogen/deuterium exchange mass spectrometry: Correlation between experiments and simulation

Hydrogen deuterium exchange mass spectrometry (HDX‐MS) is a powerful technique for studying protein dynamics, which is an important factor governing protein functions. However, the process of hydrogen/deuterium exchange (HDX) of proteins is highly complex and the underlying mechanism has not yet been fully elucidated. Meanwhile, molecular dynamics (MD) simulation is a computational technique that can be used to elucidate HDX behaviour on proteins and facilitate interpretation of HDX‐MS data. This article aims to summarize the current understandings on the mechanism of HDX and its correlation with MD simulation, to discuss the recent developments in the techniques of HDX‐MS and MD simulation and to extend the perspectives of these two techniques in protein dynamics study.     

Analysis of the unfolding process of green fluorescent protein by molecular dynamics simulation

Molecular dynamics simulation of the enforced stretching of circularly permuted green fluorescent protein (cpGFP) was performed to observe the detailed process of unfolding of beta-sheets in cpGFP and to clarify the structural change arising from the force. The simulation using the generalized Born method with original force field parameters enabled us to observe the unfolding process of the entire region of the protein and to clarify atom motion of the individual domain during the stretching. The force required for the stretching of cpGFP was estimated from the differential of the computed potential energy. A prominent rise in force appeared three times during the stretching. The amplitude and the position of these three peaks were consistent with the observation in atomic force microscopy (AFM) experiments. Further, the movements of atoms involved in each peak were shown to be closely related to the dissociation of hydrogen bonds. Additional simulations for the unfolding process of titin and spectrin also gave satisfactory interpretation of the results of previous AFM experiments. The difference in the enforced stretching process between cpGFP and wild-type GFP was further discussed through the MD simulation.     

Car-Parrinello and path integral molecular dynamics study of the hydrogen bond in the chloroacetic acid dimer system

We have studied the double proton transfer (DPT) reaction in the cyclic dimer of chloroacetic acid using both classical and path integral Car-Parrinello molecular dynamics. We also attempt to quantify the errors in the potential energy surface that arise from the use of a pure density functional. In the classical dynamics a clear reaction mechanism can be identified, where asynchronized DPT arises due to coupling between the O-H stretching oscillator and several low energy intermolecular vibrational modes. This mechanism is considerably altered when quantum tunneling is permitted in the simulation. The introduction of path integrals leads to considerable changes in the thermally averaged molecular geometry, leading to shorter and more centered hydrogen bond linkages.