共查询到20条相似文献,搜索用时 15 毫秒
1.
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. 相似文献
2.
Pan Chen Yoshiharu Nishiyama Jean-Luc Putaux Karim Mazeau 《Cellulose (London, England)》2014,21(2):897-908
We have performed molecular dynamics calculations using a revised version of the Gromos56Acarbo 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. 相似文献
3.
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. 相似文献
4.
5.
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 Rgx2,Rgx2,Rgz2 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 Rgx2(or Rgy2) 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 Ntr 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. 相似文献
6.
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. 相似文献
7.
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. 相似文献
8.
Lindorff-Larsen K Trbovic N Maragakis P Piana S Shaw DE 《Journal of the American Chemical Society》2012,134(8):3787-3791
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. 相似文献
9.
Muniz-Miranda F Pagliai M Cardini G Righini R 《The journal of physical chemistry. A》2012,116(9):2147-2153
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. 相似文献
10.
11.
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. 相似文献
12.
Ramos J Cruz VL Martínez-Salazar J Campillo NE Páez JA 《Physical chemistry chemical physics : PCCP》2011,13(9):3660-3668
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. 相似文献
13.
A molecular orbital explanation of bond distance variation caused by hydrogen bond formation 总被引:4,自引:0,他引:4
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)2, HCONH2 (formamide) crystal and BF3·2H2O crystal. 相似文献
14.
Ben Hanson 《European Polymer Journal》2010,46(12):2310-2320
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. 相似文献
15.
16.
Bandyopadhyay S Chakraborty S Bagchi B 《Journal of the American Chemical Society》2005,127(47):16660-16667
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. 相似文献
17.
18.
Liwen Huang Pui‐Kin So Zhong‐Ping Yao 《Rapid communications in mass spectrometry : RCM》2019,33(Z3):83-89
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. 相似文献
19.
Hisatomi Y Katagiri D Neya S Hara M Hoshino T 《The journal of physical chemistry. B》2008,112(29):8672-8680
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. 相似文献
20.
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. 相似文献