Vesicles from mixtures of bipolar archaebacterial lipids with egg phosphatidylcholine

Bipolar lipids from the membranes of archaebacterium Caldariella acidophila can form small unilamellar liposomes, when sonicated from lipid mixtures containing at least 25 mol% egg phosphatidylcholine. With increasing contents of archaebacterial lipid the inner radius of highly sonicated vesicles increases (from approx. 90 Å to approx. 160 Å) concomitant with an enhanced asymmetric distribution of the phosphatidylcholine molecules towards the outer face of the lipid bilayer membranes.     

Role of hydrophobic core on the thermal stability of proteins-molecular dynamics simulations on a single point mutant of sso7d

Abstract The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Ab initio and molecular dynamics studies of solid β-HMX: effects of hydrostatic pressure and high temperature

Using first-principles density functional theory and classical molecular dynamics (MD), the structural, electronic and mechanical properties of the energetic material β-HMX have been studied. The crystal structure optimised by the local density approximation calculations compares reasonably with the experimental data. Electronic band structure and density of states indicate that β-HMX is an insulator with a band gap of 3.059 eV. The pressure effect on the crystal structure and physical properties has been investigated in the range of 0–40 GPa. The crystal structure and electronic properties change slightly as the pressure increases from 0 to 2.5 GPa; when the pressure is above 2.5 GPa, further increment of the pressure results in significant changes in crystal structure. There is a larger compression along the b-axis than along the a- and c-axes. Isothermal–isobaric MD simulations on β-HMX were performed in the temperature range of 5–400 K. Phase transition at 360 K, corresponding to a volume interrupt, was found. The computed thermal expansion coefficients show anisotropic behaviour with a slightly larger expansion along the b- and c-axes than along the a-axis. In the temperature range of 5–360 K, β-HMX possesses good plasticity and its stiffness decreases with increasing the temperature.     

Role of hydrophobic core on the thermal stability of proteins - molecular dynamics simulations on a single point mutant of Sso7d abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 s long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Role of Hydrophobic Core on the Thermal Stability of Proteins—Molecular Dynamics Simulations on a Single Point Mutant of Sso7d

Abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Molecular dynamics and quantum chemical studies of solvent effects on cyclo glycylglycine and glycylalanine dipeptides

Hydrogen bond (H-bond) interactions between the two cyclo dipeptides, cyclo(glycyl-glycine) (CGG) and cyclo(glycyl-alanine) (CGA), and water have been studied using molecular dynamics (MD) and quantum chemical methods. The MD studies have been carried out on CGG and CGA in water using fixed charge force field (AMBER ff03) for over 10 ns with a MD time step of 2 fs. The results of this study show that the solvation pattern influences the conformations of the cyclo dipeptides. Following molecular simulations, post Hartree–Fock and density functional theory methods have been used to explore the molecular properties of the cyclo dipeptides in gaseous and aqueous phase environments. The self-consistent reaction field theory has been used to optimise the cyclopeptides in diethyl ether (? = 4.3) and water (? = 78.5), and the solvent effects have been analysed. A cluster of eight water molecules leads to the formation of first solvation shell of CGG and CGA and the strong H-bonding mainly contributes to the interaction energies. The H-bond interactions have been analysed by the calculation of electron density ρ(r) and its Laplacian ▽²ρ(r) at bond critical points using atoms in molecules theory. The natural bond orbital analysis was carried out to reveal the nature of H-bond interactions. In the solvated complexes, the keto carbons registered the maximum NMR chemical shifts.     

Extracting Curvature Preferences of Lipids Assembled in Flat Bilayers Shows Possible Kinetic Windows for Genesis of Bilayer Asymmetry and Domain Formation in Biological Membranes

Studies on the assembly of pure lipid components allow mechanistic insights toward understanding the structural and functional aspects of biological membranes. Molecular dynamic (MD) simulations on membrane systems provide molecular details on membrane dynamics that are difficult to obtain experimentally. A large number of MD studies have remained somewhat disconnected from a key concept of amphipathic assembly resulting in membrane structures—shape parameters of lipid molecules in those structures in aqueous environments. This is because most of the MD studies have been done on flat lipid membranes. With the above in view, we analyzed MD simulations of 26 pure lipid patches as a function of (1) lipid type(s) and (2) time of MD simulations along with 35–40 ns trajectories of five pure lipids. We report, for the first time, extraction of curvature preferences of lipids from MD simulations done on flat bilayers. Our results may lead to mechanistic insights into the possible origins of bilayer asymmetries and domain formation in biological membranes.     

Low-temperature molecular dynamics simulations of horse heart cytochrome c and comparison with inelastic neutron scattering data

Molecular dynamics (MD) simulation combined with inelastic neutron scattering can provide information about the thermal dynamics of proteins, especially the low-frequency vibrational modes responsible for large movement of some parts of protein molecules. We performed several 30-ns MD simulations of cytochrome c (Cyt c) in a water box for temperatures ranging from 110 to 300 K and compared the results with those from experimental inelastic neutron scattering. The low-frequency vibrational modes were obtained via dynamic structure factors, S(Q, ω), obtained both from inelastic neutron scattering experiments and calculated from MD simulations for Cyt c in the same range of temperatures. The well known thermal transition in structural movements of Cyt c is clearly seen in MD simulations; it is, however, confined to unstructured fragments of loops Ω₁ and Ω₂; movement of structured loop Ω₃ and both helical ends of the protein is resistant to thermal disturbance. Calculated and experimental S(Q, ω) plots are in qualitative agreement for low temperatures whereas above 200 K a boson peak vanishes from the calculated plots. This may be a result of loss of crystal structure by the protein–water system compared with the protein crystal.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.     

Structure prediction and functional analyses of a thermostable lipase obtained from Shewanella putrefaciens

Previous experimental studies on thermostable lipase from Shewanella putrefaciens suggested the maximum activity at higher temperatures, but with little information on its conformational profile. In this study, the three-dimensional structure of lipase was predicted and a 60 ns molecular dynamics (MD) simulation was carried out at temperatures ranging from 300 to 400 K to better understand its thermostable nature at the molecular level. MD simulations were performed in order to predict the optimal activity of thermostable lipase. The results suggested strong conformational temperature dependence. The thermostable lipase maintained its bio-active conformation at 350 K during the 60 ns MD simulations.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.     

Simulated dynamics of Ne@C60 aggregates beyond dissociation

Molecular dynamics (MD) computer simulations are utilized to better understand the dynamics of small (N = 5) endohedral Ne@C₆₀ aggregates. Multiple runs at various temperatures are used to increase the reliability of our statistics. The aggregate holds together until somewhere between T = 1150 and 1200 K, where it dissociates, showing no intermediate sign of melting or fullerene disintegration. When the temperature is increased to around T = 4000 K, the encapsulated neon atoms begin to leave the aggregate, with the fullerene molecules still remaining intact. At temperatures near T = 4400 K, thermal disintegration of the fullerenes preempts the aggregate dissociation. Above this temperature neon atoms are more quickly released and the fullerenes form a larger connected structure, with bonding taking place in atom pairs from different original fullerene molecules. Escape constants and half lives are calculated for the temperature range 4000 K ≤ T ≤ 5000 K. The agreements and disagreements of results of this work with experiments suggest that classical MD simulations are useful in describing fullerene systems at low temperatures and near disintegration, but require development of new techniques before it is possible to accurately model windowing at temperatures below T = 3000 K.     

Non-equilibrium molecular dynamics study of nanoscale thermal contact resistance

Interfaces play an important role in microscale and nanoscale heat transfer processes with molecular dynamics (MD) simulations often used to study these interfacial phenomena. In this study, two models were used to simulate thermal conduction across micro contact points and the thermal contact resistance using non-equilibrium molecular dynamics simulations with consideration of the near field radiation. When the ratio of the length of the micro contact to the length of the conduction region is less than 0.125, the influence of the near field radiation should be considered; but when the ratio is larger than 0.2, it can be neglected. When the computational domain sizes are 8.50 × 10.62 × 8.50 nm and 10.62 × 10.62 × 10.62 nm, the MD results show that the thermal contact resistance exponentially increases with decreasing area of the micro contact point and increases with increasing micro contact layer thickness. The MD thermal contact resistances in nanoscale are much larger than that of the classical thermal analysis since the material thermal conductivity reduction is ignored in the classical model. The results also show that material defects increase the thermal resistance.     

大鼠心肌线粒体内、外膜磷脂动态结构的研究

我们以DPH为荧光探针.用毫微秒荧光分光光度计测定了大鼠心肌线粒体及线粒体内、外膜的动态微细结构;用HPLC分析了磷脂组成.实验结果提示.完整线粒体膜流动性主要反映了线粒体外膜的运动状态.线粒体内膜微粘度及磷脂分子摇动角大于外膜,扩散速率小于外膜.除去了蛋白质的线粒体内、外膜磷脂脂质体膜流动性无明显差异.提示线粒体内膜的高微粘度与膜中所含有的多量蛋白有关.     

Insight into the impact of environments on structure of chimera C3 of human β-defensins 2 and 3 from molecular dynamics simulations

C3 is a chimera from human β-defensins 2 and 3 and possesses higher antimicrobial activity compared with its parental molecules, so it is an attractive candidate for clinical application of antimicrobial peptides. In continuation with the previous studies, molecular dynamics (MD) simulations were carried out for further investigating the effect of ambient environments (temperature and bacterial membrane) on C3 dynamics. Our results reveal that C3 has higher flexibility, larger intensity of motion, and more relevant secondary structural changes at 363 K to adapt the high temperature and maintain its antimicrobial activity, comparison with it at 293 K; when C3 molecule associates with the bacterial membrane, it slightly fluctuates and undergoes local conformational changes; in summary, C3 molecule demonstrates stable conformations under these environments. Furthermore, MD results analysis show that the hydrophobic contacts, the hydrogen bonds, and disulfide bonds in the peptide are responsible for maintaining its stable conformation. In addition, our simulation shows that C3 peptides can make anionic lipids clustered in the bacterial membrane; it means that positive charges and pronounced regional cationic charge density of C3 are most key factors for its antimicrobial activity.     

Effect of periodic box size on aqueous molecular dynamics simulation of a DNA dodecamer with particle-mesh Ewald method

The particle-mesh Ewald (PME) method is considered to be both efficient and accurate for the evaluation of long-range electrostatic interactions in large macromolecular systems being studied by molecular dynamics simulations. This method assumes "infinite" periodic boundary conditions resembling the symmetry of a crystal environment. Can such a "solid-state" method accurately portray a macromolecular solute such as DNA in solution? To address this issue, we have performed three 1500-ps PME molecular dynamics (MD) simulations, each with a different box size, on the d(CGCGA6CG)-(CGT6CGCG) DNA dodecamer. The smallest box had the DNA solvated by a layer of water molecules of at least 5 A along each orthogonal direction. The intermediate size box and the largest box had the DNA solvated by a layer of water molecules of at least 10 A and 15 A, respectively, along each orthogonal direction. The intermediate size box in the present study is similar to the box size currently chosen by most workers in the field. Based on a comparison of RMSDs and curvature for this single DNA dodecamer sequence, the larger two box sizes do not appear to afford any extra benefit over the smallest box. The implications of this finding are briefly discussed.     

CD4 binding partially locks the bridging sheet in gp120 but leaves the beta2/3 strands flexible

The structure of the free form HIV gp120, critical for therapeutic agent development, is unavailable due to its high flexibility. Previous thermodynamic data, structural analysis and simulation results have suggested a large conformational change in the core domain upon CD4 binding. The bridging sheet, which consists of four beta-strands with beta20/21 nestling against the inner/outer domains and beta2/3 facing outward, more exposed to the solvent, was proposed to be unfolded in the native state. In order to test this proposition and to characterize the native conformations, we performed potential mean force (PMF) molecular dynamics (MD) simulations on the CD4-bound crystal structure. We pushed the bridging sheet away from the inner and outer domain to explore the accessible conformational space for the bridging sheet. In addition, we performed conventional MD simulations on structures with the bridging sheet partially unfolded to investigate the stability of the association between the inner and outer domains. Based on the free energy profiles, we find that the whole bridging sheet is unlikely to unfold without other concurrent conformational changes. On the other hand, the partial bridging sheet, beta strands 2/3, can switch its conformation from the folded to the unfolded state. Furthermore, relaxation of conformation with partially unfolded bridging sheet through MD simulations leads to a conformation with beta strands 20/21 quickly re-anchoring against the inner and outer domains. Such a conformation, although lacking some of the hydrophobic interactions present in the CD4-bound structure, displayed high stability as further indicated by other restrained MD simulations. The relevance of this conformation to the free form structure and the pathway for conformational change from the free form to the CD4-bound structure is discussed in detail in light of the available unliganded SIV gp120 crystal structure.     

Structure and dynamics of K channel pore-lining helices: a comparative simulation study

Isolated pore-lining helices derived from three types of K-channel have been analyzed in terms of their structural and dynamic features in nanosecond molecular dynamics (MD) simulations while spanning a lipid bilayer. The helices were 1) M1 and M2 from the bacterial channel KcsA (Streptomyces lividans), 2) S5 and S6 from the voltage-gated (Kv) channel Shaker (Drosophila melanogaster), and 3) M1 and M2 from the inward rectifier channel Kir6.2 (human). In the case of the Kv and Kir channels, for which x-ray structures are not known, both short and long models of each helix were considered. Each helix was incorporated into a lipid bilayer containing 127 palmitoyloleoylphosphatidylcholine molecules, which was solvated with approximately 4000 water molecules, yielding approximately 20, 000 atoms in each system. Nanosecond MD simulations were used to aid the definition of optimal lengths for the helix models from Kv and Kir. Thus the study corresponds to a total simulation time of 10 ns. The inner pore-lining helices (M2 in KcsA and Kir, S6 in Shaker) appear to be slightly more flexible than the outer pore-lining helices. In particular, the Pro-Val-Pro motif of S6 results in flexibility about a molecular hinge, as was suggested by previous in vacuo simulations (, Biopolymers. 39:503-515). Such flexibility may be related to gating in the corresponding intact channel protein molecules. Analysis of H-bonds revealed interactions with both water and lipid molecules in the water/bilayer interfacial region. Such H-bonding interactions may lock the helices in place in the bilayer during the folding of the channel protein (as is implicit in the two-stage model of membrane protein folding). Aromatic residues at the extremities of the helices underwent complex motions on both short (<10 ps) and long (>100 ps) time scales.     

Molecular dynamics and IR spectroscopy in investigation of phase transitions in molecular crystal 4,4'-bis(11-hydroxy-1-undecyloxy)biphenyl

Molecular dynamics (MD) simulations combined with temperature-dependent IR spectroscopic measurements were used to study phase transitions in molecular crystals of the mesogenic diol 4,4'-bis(11-hydroxy-1-undecyloxy)biphenyl. DSC measurements revealed four phase transitions in this molecular crystal at approximately 327.1 K, 389.8 K, 419.1 K and 431.9 K. Analysis of the dynamic trajectories at temperatures of 300 K, 360 K, 400 K and 480 K revealed changes in conformation of the mesogenic diol molecules and consequently changes in crystal packing and crystal structure in the temperature range 300-480 K and enabled us to understand the mechanism of the phase transitions.     

Electrodissociation of clathrate-like structures

In the present work, we develop molecular dynamics (MD) simulations in the NPT (isobaric–isothermic) ensemble to analyse the effect of an external electrostatic field over a cubic methane hydrate crystallite. The amplitude of the field is in the range 0.5–3.0 V/nm. For the simulations, we used the SPC/E rigid water model and a single-site model for methane at a temperature of 248 K and a pressure of 20 bar. When the external electrostatic field is applied, the water dipoles are oriented in such a way that the methane molecules can diffuse far away from the water cages, hence the clathrate dissociation takes place. This last phenomenon was observed for intensities above 1.5 V/nm. Taking the final configuration of each run as input, we develop a new set of MD simulations, and we observe that the stable clathrate is not recovered immediately when the external electrostatic field is turned off due to limitations in the simulation time.