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.     

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

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

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.     

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.     

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.     

A single CH/pi weak hydrogen bond governs stability and the photocycle of the photoactive yellow protein

Importance of the CH/pi interaction on the structure and function of the photoactive yellow protein (PYP) was substantiated. Focusing on the phenyl ring of Phe6 adjacent to the alkyl chain of Lys123, the mutants for these amino acid residues were characterized. The results demonstrated that the mutants lacking the pi-electron at position 6 or the alkyl chain at position 123 show substantial malfunction. This is a clear example that single CH/pi weak interaction plays a crucial role in the normal action of the protein.     

Crystallization of helical oligomers with chirality selection. I. A molecular dynamics simulation for bare helix

Helical polymers often exhibit pronounced chirality recognition during crystallization. By molecular dynamics simulation, we have already shown that the helical polymers crystallize with or without marked chirality selection depending on structural details of the polymer molecules. We have there classified the helical polymers into two categories: the bare helices made of only backbone atoms which show rather tolerant chirality selection, and the general helices with large side groups showing strict chirality recognition. Polymer crystallization is in general largely hampered and retarded by slow dynamics of the entangled chains, and therefore short helical oligomers are very suitable models for studying the chiral crystallization. We here report on molecular simulations of crystallization in the bare helical oligomer molecules by the use of Monte Carlo and molecular dynamics simulations. First we confirm the low temperature chiral crystal phase and the reversible order-disorder transition. We also observe frequent inversions of the helical sense, and the helix reversal defects propagating along the chains. Then we investigate crystallization from the melt into the chiral crystal phase. We find that the crystallization rate depends very sensitively on the degree of undercooling. The crystallization is found to be the first order transition that conforms well to the traditional picture of crystal growth in small molecules. Even when the crystallization directly into the chiral crystal phase is conducted, marked chirality selections are not observed at the early stage of crystallization; the chains adhere to the crystal surfaces selecting their helical senses rather at random resulting in racemic crystallites. The isothermal crystallization for a sufficiently long time, however, yields lamellar crystals composed of well-developed chiral domains, the growth of which seems to be accomplished through the transition back into the ordered chiral crystal phase.     

Photoionization-induced water migration in the hydrated trans-formanilide cluster cation revealed by gas-phase spectroscopy and ab initio molecular dynamics simulation

Photoionization-induced water migration in the trans-formanilide-water 1:1 cluster, FA-(H(2)O)(1), has been investigated by using IR-dip spectroscopy, quantum chemical calculations, and ab initio molecular dynamics simulations. In the S(0) state, FA-(H(2)O)(1) has two structural isomers, FA(NH)-(H(2)O)(1) and FA(CO)-(H(2)O)(1), where a water molecule is hydrogen-bonded (H-bonded) to the NH group and the CO group, respectively. In addition, the S(1)-S(0) origin transition of FA(CO)-(H(2)O)(2), where a water dimer is H-bonded to the CO group, was observed only in the [FA-(H(2)O)(1)](+) mass channel, indicating that one of the water molecules evaporates completely in the D(0) state. These results are consistent with a previous report [Robertson, E. G. Chem. Phys. Lett., 2000, 325, 299]. In the D(0) state, however, [FA-(H(2)O)(1)](+) produced by photoionization via the S(1)-S(0) origin transitions of FA(NH)-(H(2)O)(1) and FA(CO)-(H(2)O)(1) shows essentially the same IR spectra. Compared with the theoretical calculations, [FA-(H(2)O)(1)](+) can be assigned to [FA(NH)-(H(2)O)(1)](+). This means that the water molecule in [FA-(H(2)O)(1)](+) migrates from the CO group to the NH group when [FA-(H(2)O)(1)](+) is produced by photoionization of FA(CO)-(H(2)O)(1). [FA-(H(2)O)(1)](+) produced by photoionization of FA(CO)-(H(2)O)(2) also shows the IR spectrum corresponding to [FA(NH)-(H(2)O)(1)](+). In this case, the water migration from the CO group to the NH group occurs with the evaporation of a water molecule. Ab initio molecular dynamics simulations revealed the water migration pathway in [FA-(H(2)O)(1)](+). The calculations of classical electrostatic interactions show that charge-dipole interaction between FA(+) and H(2)O induces an initial structural change in [FA-(H(2)O)(1)](+). An exchange repulsion between the lone pairs of the CO group and H(2)O in [FA-(H(2)O)(1)](+) also affects the initial direction of the water migration. These two factors play important roles in determining the initial water migration pathway.