Molecular dynamics simulation of liquid water under the influence of an external electric field

Molecular dynamics simulations of liquid water were performed at 258K and a density of 1.0?g/cm³ under various applied external electric field, ranging 0～10¹⁰?V/m. The influence of external field on structural and dynamical properties of water was investigated. The simple point charge (SPC) model is used for water molecules. An enhancement of the water hydrogen bond structure with increasing strength of the electric field has been deduced from the radial distribution functions and the analysis of hydrogen bonds structure. With increasing field strength, water system has a more perfect structure, which is similar to ice structure. However, the electrofreezing phenomenon of liquid water has not been detected since the self-diffusion coefficient was very large. The self-diffusion coefficient decreases remarkably with increasing strength of electric field and the self-diffusion coefficient is anisotropic.     

Examination of structural stability and phase transitions in constant-pressure first-principles molecular dynamics simulations

Constant-pressure first-principles molecular dynamics (FPMD) simulation is a powerful tool for investigations of structures in crystals. However, it needs enourmous computations so that highly accurate calculations for electronic states cannot be employed at present. In this report, we examined the reliability and applicability of constant-pressure FPMD in the study of structural properties under this limitation. Crystalline silicon was employed as a benchmark to perform constant-pressure FPMD simulations (with a deformable simulation cell). It is found that, in high pressure (metallic) phases, crystalline symmetry is broken with the present simulation conditions. Several structural transformations were realized by compression and decompression, but they are not entirely consistent with experiment. We discuss this discrepancy and conclude that the number of k point sampling in the Brillouin zone is crucial. It is recommended that constant-pressure FPMD is employed to explore candidate structures for unknown solid phases at present computational resources.     

Molecular Dynamics Study of High Temperature Phase-Separation in a H2O/N2 Mixture with Exp-6 Interactions

Abstract

Equimolar H₂O/N₂ fluid mixture was studied by molecular dynamics simulations for NVT ensemble. Calculations were performed with the modified Buckingham (exp-6) potentials at T = 2000 K. Particular attention was given to the phase separation at very high pressures relevant to a detonation environment. Calculations of pair correlation functions and local mole fractions clearly indicated the occurrence of the fluid separation into N₂-rich and H₂O-rich phase. The density at the phase boundary between homogeneous and inhomogeneous phase-separated state was determined to be p = 1.35 g/cm³ on the basis of the static cross correlation factor which is defined by the sum of the local mole fractions. The ratio of the self-diffusion coefficients of N₂ and H₂O at p < 1.35 g/cm³ was found to be approximately equal to the value predicted by the kinetic theory of the ideal gas, whereas the ratio was close to unity at the phase-separated state (p > 1.35 g/cm³). In addition, two distinctive behaviors of the system could be observed for the relaxation from the initial uniform mixture to the phase-separated fluid: at lower densities (1.35 < p < 2.0 g/cm³) the fluid mixture began to relax into the phase-separated system without obvious incubation time, while clear incubation period was associated for the separation at higher densities. During this incubation period, discontinuous jumps in the mean square displacements were found.     

Molecular dynamics simulation of the Al2O3 film structure during atomic layer deposition

Based on an understanding of atomic layer deposition (ALD) from prior experimental and computational results, all-atom molecular dynamics (MD) simulations are used to model the Al₂O₃ film structure and composition during ALD processing. By separating the large time-scale surface reactions from the small time-scale structural relaxation, we have focused on the growth dynamics of amorphous Al₂O₃ films at the atomic scale. The simulations are able to reproduce some important properties and growth mechanisms of Al₂O₃ ALD films, and hence provide a bridge between atomic-level information and experimental measurements. Information about the evolution of the microscopic structures of the Al₂O₃ films is generated, and the influence of operation parameters on the Al₂O₃ ALD process. The simulations predict a strong influence of the initial surface composition and process temperature on the surface roughness, growth rate and growth mode of the deposited films.     

Exploring the structural determinants of novel xanthine derivatives as A2B adenosine receptor antagonists: a computational study

Adenosine is a ubiquitous endogenous nucleoside that controls numerous physiological functions via interacting with its specific G-coupled receptors. Activation of adenosine receptors (AdoRs), particularly A_2B AdoRs promotes the release of inflammatory cytokines; reduces vascular permeabilization and induces angiogenesis, thereby making A_2B AdoR becomes a potentially pharmacological target for drug development. Presently, for investigating the structural determinants of 164 xanthine derivatives as A_2B AdoR antagonists, we performed an in silico study integrating with 3D-QSAR, docking and molecular dynamics (MD) simulation. The obtained optimal model shows strong predictability (Q²?=?0.647, R²_ncv?=?0.955, and R²_pred?=?0.848). Additionally, to explore the binding mode of the ligand with A_2B AdoR and to understand their binding mechanism, docking analysis, MD simulations (20?ns), and the calculation of binding free energy were also carried out. Finally, the structural determinants of these xanthine derivatives were identified and a total of 20 novel A_2B AdoR antagonists with improved potency were computationally designed, and their synthetic feasibility and selectivity were also evaluated. The information derived from the present study offers a better appreciation for exploring the interaction mechanism of the ligand with A_2B AdoR, which could be helpful for designing novel potent A_2B AdoR antagonists.

Communicated by Ramaswamy H. Sarma     

The solvation of bromide anion in acetonitrile: a structural study based on the combination of theoretical calculations and X-ray absorption spectroscopy

This work studies the solvation of bromide in acetonitrile by combining quantum mechanics, computer simulations and X-ray absorption near edge structure (XANES) spectroscopy. Three different sets of interaction potentials were tested, one of them derived from literature and the other two are simple modifications of the previous one to include specificities of the bromide–acetonitrile interactions. Results for microsolvation of bromide were obtained by quantum mechanical optimization and classical minimization of small clusters [Br(ACN) _n ]^?  (n = 9, 20). Analysis of molecular dynamics (MD) simulations has provided structural, dynamic and energetic aspects of the solvation phenomenon. The theoretical computation of Br K-edge XANES spectrum in solution using the structural information obtained from the different simulations has allowed the comparison among the three different potentials, as well as the examination of the main structural and dynamic factors determining the shape of the experimental spectrum.     

Molecular dynamics simulations of nafion and sulfonated poly ether sulfone membranes II. Dynamic properties of water and hydronium

We measured the self-diffusion coefficients of water in a Nafion membrane and two sulfonated polyethersulfone (SPES) membranes with varying ion-exchange capacities (IEC) in terms of relative humidity using the pulse field gradient NMR (PFG-NMR) technique. The self-diffusion coefficients were plotted against the number of water molecules per sulfonic acid group, λ, and compare these values with the results of molecular dynamics (MD) simulations. Classical MD simulations for all membranes were carried out using a consistent force field at λ = 3, 6, 9, 12, and 15. The dynamic properties of water (H₂O) and hydronium (H₃O⁺) on a molecular level were estimated as self-diffusion coefficients and residence times around a sulfonate group ( \textSO₃^- {\text{SO}}_3^{-} ). The diffusion coefficients of H₂O and H₃O⁺ followed the order, Nafion > SPES with IEC = 1.4 > SPES with IEC = 1.0 > SPES with IEC = 0.75, which agreed with the experimental data. The residence time distribution of H₂O around \textSO₃^- {\text{SO}}_3^{-} in Nafion was in the range of 1–6 ps, whereas H₂O in the SPES exhibited a residence time of greater than 20 ps.     

Structural Assignment of Spectra by Characterization of Conformational Substates in Bound MbCO

Residue motions of the distal heme pocket and bound CO ligand of carbonmonoxy Myoglobin are studied using a combination of molecular dynamics simulations and quantum chemical methods. Using mixed quantum mechanics/molecular mechanics calculations together with sampling from molecular dynamics simulations (QM/MM(MD)), the experimentally observed spectroscopic A₀ and A₁ substates of the bound CO ligand are assigned to the open and closed conformation of His⁶⁴ and the His_ɛ⁶⁴ tautomer, respectively. Several previously proposed origins of the A₃ substate, including rotamers of the doubly protonated His⁶⁴H⁺ side chain, His⁶⁴H⁺ inside the distal pocket, and cooperative motions with Arg⁴⁵, are investigated with QM/MM(MD). However, the signatures of the calculated infrared spectra do not agree with the experimentally observed ones. For additional insight on this, extensive molecular dynamics simulations are used together with improved electrostatics for the bound ligand. A CO fluctuating charge model is developed to describe the ab initio dipole and quadrupole moments of the bound ligand. CO absorption spectra are then obtained directly from the dynamics simulations. Finally, the electrostatics of the heme pocket is examined in detail in an attempt to determine the structural origins of the observed spectroscopic A-states from MD simulations. However, contrary to related simulations for unbound CO in myoglobin, the shifts and splittings for carbonmonoxy Myoglobin are generally small and difficult to relate to structural change. This suggests that coupling of the CO motion to other degrees of freedom, such as the Fe-CO stretching and bending, is important to correctly describe the dynamics of bound CO in myoglobin.     

Dynamics of Small Molecules in a Dense Polymer Matrix: Molecular Dynamics Studies

Abstract

The anomalous diffusion regime appearing in the self-diffusion of small molecules in bulk amorphous polymers has been extensively studied by molecular dynamics simulations. A rather long simulation of duration ～ 10 ^?8 s is performed on a polyethylene-like simple polymer model containing either oxygen molecules or helium atoms as a diffusant. Dynamic properties evaluated for these diffusants are the mean-square displacement, the van Hove self correlation function, and the self part of the density autocorrelation. It is first confirmed that the anomalous diffusion regime appears in a few hundred picoseconds for oxygen molecule, while the Einstein relation adopted beyond this regime results in the self-diffusion coefficient of the order of ～ 10^?5 cm²/s. This anomaly is still observed for helium that diffuses much faster than oxygen. In the anomalous diffusion regime, it is found that the correlation functions for the two diffusants show characteristic features and become essentially the same as time is scaled appropriately. These features allow the estimation of the two characteristic spatial scales which are probably dominated by the microstructure of the polymer matrix, namely, the cage size and the distance between adjacent cages. The time dependence of the mean-square displacements of the two diffusants can be well interpreted by these characteristic spatial scales as time is scaled with the self-diffusion coefficients. It is shown that the anomalous diffusion regime arises from the inhomogeneous microstructure of the polymer matrix.     

Molecular Dynamics Simulations for 1:1 Solvent Primitive Model Electrolyte Solutions

Abstract

Molecular dynamics (MD) simulations at constant temperature have been carried out for systems of 1:1 solvent primitive model (SPM) electrolyte solutions. Equilibrium thermodynamics, mean cluster size, self-diffusion coefficients, and collision frequencies were computed to examine the electrostatic effects on the structural and dynamical properties. Coherent ionic cluster motion was deduced from a cluster analysis and from the dependence of the velocity and force autocorrelation functions (FACFs). The resulting MD data for the collision frequencies and self-diffusivities of both ions and hard-spheres were shown to be in good agreement with the theoretical predictions.     

Analysis of Sub-τc and Supra-τc Motions in Protein Gβ1 Using Molecular Dynamics Simulations

The functions of proteins depend on the dynamical behavior of their native states on a wide range of timescales. To investigate these dynamics in the case of the small protein Gβ1, we analyzed molecular dynamics simulations with the model-free approach of nuclear magnetic relaxation. We found amplitudes of fast timescale motions (sub-τ_c, where τ_c is the rotational correlation time) consistent with S² obtained from spin relaxation measurements as well as amplitudes of slow timescale motions (supra-τ_c) in quantitative agreement with S² order parameters derived from residual dipolar coupling measurements. The slow timescale motions are associated with the large variations of the ³J couplings that follow transitions between different conformational substates. These results provide further characterization of the large structural fluctuations in the native states of proteins that occur on timescales longer than the rotational correlation time.     

Mode‐coupling Smoluchowski dynamics of a double‐stranded DNA oligomer

The local dynamics of a double‐stranded DNA d(TpCpGpCpG)₂ is obtained to second order in the mode‐coupling expansion of the Smoluchowski diffusion theory. The time correlation functions of bond variables are derived and the ¹³C‐nmr spin–lattice relaxation times T₁ of different ¹³C along the chains are calculated and compared to experimental data from the literature at three frequencies. The DNA is considered as a fluctuating three‐dimensional structure undergoing rotational diffusion. The fluctuations are evaluated using molecular dynamics simulations, with the ensemble averages approximated by time averages along a trajectory of length 1 ns. Any technique for sampling the configurational space can be used as an alternative. For a fluctuating three‐dimensional (3D) structure using the three first‐order vector modes of lower rates, higher order basis sets of second‐rank tensor are built to give the required mode coupling dynamics. Second‐ and even first‐order theories are found to be in close agreement with the experimental results, especially at high frequency, where the differences in T₁ for ¹³C in the base pairs, sugar, and backbone are well described. These atomistic calculations are of general application for studying, on a molecular basis, the local dynamics of fluctuating 3D structures such as double‐helix DNA fragments, proteins, and protein–DNA complexes. © 1999 John Wiley & Sons, Inc. Biopoly 50: 613–629, 1999     

Molecular dynamics force probe simulations of antibody/antigen unbinding: entropic control and nonadditivity of unbinding forces

Unbinding of a spin-labeled dinitrophenyl (DNP) hapten from the monoclonal antibody AN02 F(ab) fragment has been studied by force probe molecular dynamics (FPMD) simulations. In our nanosecond simulations, unbinding was enforced by pulling the hapten molecule out of the binding pocket. Detailed inspection of the FPMD trajectories revealed a large heterogeneity of enforced unbinding pathways and a correspondingly large flexibility of the binding pocket region, which exhibited induced fit motions. Principal component analyses were used to estimate the resulting entropic contribution of approximately 6 kcal/mol to the AN02/DNP-hapten bond. This large contribution may explain the surprisingly large effect on binding kinetics found for mutation sites that are not directly involved in binding. We propose that such "entropic control" optimizes the binding kinetics of antibodies. Additional FPMD simulations of two point mutants in the light chain, Y33F and I96K, provided further support for a large flexibility of the binding pocket. Unbinding forces were found to be unchanged for these two mutants. Structural analysis of the FPMD simulations suggests that, in contrast to free energies of unbinding, the effect of mutations on unbinding forces is generally nonadditive.     

Probing the interaction between cHAVc3 peptide and the EC1 domain of E-cadherin using NMR and molecular dynamics simulations

The goal of this work is to probe the interaction between cyclic cHAVc3 peptide and the EC1 domain of human E-cadherin protein. Cyclic cHAVc3 peptide (cyclo(1,6)Ac-CSHAVC-NH₂) binds to the EC1 domain as shown by chemical shift perturbations in the 2D ¹H,-¹⁵N-HSQC NMR spectrum. The molecular dynamics (MD) simulations of the EC1 domain showed folding of the C-terminal tail region into the main head region of the EC1 domain. For cHAVc3 peptide, replica exchange molecular dynamics (REMD) simulations generated five structural clusters of cHAVc3 peptide. Representative structures of cHAVc3 and the EC1 structure from MD simulations were used in molecular docking experiments with NMR constraints to determine the binding site of the peptide on EC1. The results suggest that cHAVc3 binds to EC1 around residues Y36, S37, I38, I53, F77, S78, H79, and I94. The dissociation constants (K_d values) of cHAVc3 peptide to EC1 were estimated using the NMR chemical shifts data and the estimated K_ds are in the range of .5 × 10^?5–7.0 × 10^?5 M.     

MD — Simulation of Diffusion of Methane in Zeolites of Type LTA

Using molecular dynamical computer simulations (MD) the dynamics of kinetic processes in zeolites will be discussed on a molecular level. Small changes in lattice parameters can cause dramatical changes in the diffusion coefficient. The presence of cations Na⁺, Ca²⁺ also strongly influences the diffusion. Changes of the self-diffusivities will be discussed that appear if a vibrating lattice instead of a rigid one is used. Nonequilibrium simulations show the correlation between transport-diffusion and self-diffusion in zeolites.     

Toward a high-resolution structure of IP3R channel

The ability of cells to maintain low levels of Ca²⁺ under resting conditions and to create rapid and transient increases in Ca²⁺ upon stimulation is a fundamental property of cellular Ca²⁺ signaling mechanism. An increase of cytosolic Ca²⁺ level in response to diverse stimuli is largely accounted for by the inositol 1,4,5-trisphosphate receptor (IP₃R) present in the endoplasmic reticulum membranes of virtually all eukaryotic cells. Extensive information is currently available on the function of IP₃Rs and their interaction with modulators. Very little, however, is known about their molecular architecture and therefore most critical issues surrounding gating of IP₃R channels are still ambiguous, including the central question of how opening of the IP₃R pore is initiated by IP₃ and Ca²⁺. Membrane proteins such as IP₃R channels have proven to be exceptionally difficult targets for structural analysis due to their large size, their location in the membrane environment, and their dynamic nature. To date, a 3D structure of complete IP₃R channel is determined by single-particle cryo-EM at intermediate resolution, and the best crystal structures of IP₃R are limited to a soluble portion of the cytoplasmic region representing ∼15% of the entire channel protein. Together these efforts provide the important structural information for this class of ion channels and serve as the basis for further studies aiming at understanding of the IP₃R function.     

Interrelationship between water film thicknesses and contact angles and a model for CO2 adhesion

Thermodynamics of contact angle phenomena is strongly affected by the presence of thin liquid films. However, at present, studies for CO₂/brine/mineral systems only consider the films apart from contact angles. In this paper, molecular dynamics (MD) simulations have been performed to simultaneously investigate the interrelationship between water film thicknesses and water contact angles. Two types of contact angles were considered namely Young’s contact angle (no water film is present) and contact angle with film (a stable film is present). The results showed that as Young’s contact angle increased, film thickness decreased which leading to increasing of contact angle with film. The effects of CO₂-mineral pre-contact have also been investigated and it has been found that on mediate hydrophilic surfaces (Q³), water films were present when CO₂ droplets were placed above the surfaces, however, water films were absent when CO₂ droplets directly contact with the surfaces. This phenomenon implies that water films on mineral surfaces have a possibility to rupture and a film rupture mechanism for CO₂ adhesion on hydrated mineral surfaces was proposed. These results may provide new information on interactions among CO₂, water/brine and mineral to better understand the behaviour of CO₂ during geologic sequestration.     

Molecular dynamics simulations reveal initial structural and dynamic features for the A2AR as a result of ligand binding

G-protein-coupled receptors (GPCRs) are membrane proteins that have a wide variety of physiological roles. Adenosine receptors belong to the GPCR family. Adenosine receptors are implicated in many physiological disorders, such as Parkinson's disease, Huntington's disease, inflammatory and immune's disease and many others. Interestingly, crystal structures of the active and inactive conformations of the A₂-subtype adenosine receptor (A₂AR) have been solved. These two structures could be used to get insights about the conformational changes that occur during the process of activation/inactivation processes of this receptor. Therefore, two ligand-free simulations of the native active (PDB code: 3QAK) and inactive (PDB code: 3EML) conformations of the A₂AR and two halo-simulations were carried out to observe the initial conformational changes induced by coupling adenosine to the inactive conformation and caffeine to the active conformation. Furthermore, we constructed an A₂AR model that contained four thermostabilising mutations, L48A, T65A, Q89A and A54L, which had previously been determined to stabilise the bound conformation of the agonist, and we ran molecular dynamics simulations of this mutant to investigate how these point mutations might affect the inactive conformation of this receptor. This study provides insights about the initial structural and dynamic features that occur as a result of the binding of caffeine and adenosine in the active and inactive A₂AR structures, respectively, as well as the introduction of some mutations on the inactive structure of the A₂AR. Moreover, we provide useful and detailed information regarding structural features such as toggle switch and ionic lock during the activation/inactivation processes of this receptor.