IR spectra of N-methylacetamide in water predicted by combined quantum mechanical/molecular mechanical molecular dynamics simulations

We applied the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulation method in assessing IR spectra of N-methylacetamide and its deuterated form in aqueous solutions. The model peptide is treated at the Austin Model 1 (AM1) level and the induced dipole effects by the solvent are incorporated in fluctuating solute dipole moments, which are calculated using partial charges from Mulliken population analyses without resorting to any available high-level ab initio dipole moment data. Fourier transform of the solute dipole autocorrelation function produces in silico IR spectra, in which the relative peak intensities and bandwidths of major amide bands are quantitatively compatible with experimental results only when both geometric and electronic polarizations of the peptide by the solvent are dealt with at the same quantum-mechanical level. We cast light on the importance of addressing dynamic charge fluctuations of the solute in calculating IR spectra by comparing classical and QM/MM MD simulation results. We propose the adjustable scaling factors for each amide mode to be directly compared with experimental data.     

Phase transitions of methane using molecular dynamics simulations

Using a short ranged Lennard-Jones interaction and a long ranged electrostatic potential, CH4 under high pressure was modeled. Molecular dynamics simulations on small clusters (108 and 256 molecules) were used to explore the phase diagram. Regarding phase transitions at different temperatures, our numerical findings are consistent with experimental results to a great degree. In addition, the hysteresis effect is displayed in our results.     

Molecular dynamics study of methane hydrate formation at a water/methane interface

We present molecular dynamics simulation results of a liquid water/methane interface, with and without an oligomer of poly(methylaminoethylmethacrylate), PMAEMA. PMAEMA is an active component of a commercial low dosage hydrate inhibitor (LDHI). Simulations were performed in the constant NPT ensemble at temperatures of 220, 235, 240, 245, and 250 K and a pressure of 300 bar. The simulations show the onset of methane hydrate growth within 30 ns for temperatures below 245 K in the methane/water systems; at 240 K there is an induction period of ca. 20 ns, but at lower temperatures growth commences immediately. The simulations were analyzed to calculate hydrate content, the propensity for hydrogen bond formation, and how these were affected by both temperature and the presence of the LDHI. As expected, both the hydrogen bond number and hydrate content decreased with increasing temperature, though little difference was observed between the lowest two temperatures considered. In the presence of PMAEMA, the temperature below which sustained hydrate growth occurred was observed to decrease. Some of the implications for the role of PMAEMA in LDHIs are discussed.     

Nonspecific interaction forces at water-membrane interface by forced molecular dynamics simulations

Nonspecific interactions are the main driving forces for the behavior of molecules with great affinity for biologic membranes. To investigate not only the molecular details of these interactions but to estimate their magnitude as well, the theoretical method of Forced Molecular Dynamics Simulations, based on the Atomic Force Spectroscopy experimental technique, was applied. In this approach, an additional one-dimensional elastic force, representing the cantilever probe, was incorporated to the force field of a Molecular Dynamics computational program. This force represents a spring fixed on one end to a selected atom of the molecule; the other end of the spring is displaced at constant velocity to pull the molecule out of the membrane. The force experimented by the molecule due to the spring, is proportional to the spring elongation relative to its equilibrium position. This value is registered during the entire simulation, and its maximum value will determine the molecule-membrane interaction force. Nonexplicit medium simulations were carried out. Polar and apolar media were considered according to their polarizability degree and a specific dielectric constant value was assigned. In this approach, the membrane was considered as the apolar region limited by two flat surfaces with a polar aqueous medium. The potential energy discontinuity at the interfaces was smoothed by considering the polarization-induced effects using the image method. The results of this methodology are presented using a small system, a single Alanine amino acid model, which enables extended simulations in a microsecond time scale. The confinement of this amino acid at the interface reduces its degrees of freedom and forces it to adopt one of the six defined conformations. A correlation between these stable structures at the water-membrane interface and the interaction force value was determined.     

An extensible interface for QM/MM molecular dynamics simulations with AMBER

We present an extensible interface between the AMBER molecular dynamics (MD) software package and electronic structure software packages for quantum mechanical (QM) and mixed QM and classical molecular mechanical (MM) MD simulations within both mechanical and electronic embedding schemes. With this interface, ab initio wave function theory and density functional theory methods, as available in the supported electronic structure software packages, become available for QM/MM MD simulations with AMBER. The interface has been written in a modular fashion that allows straight forward extensions to support additional QM software packages and can easily be ported to other MD software. Data exchange between the MD and QM software is implemented by means of files and system calls or the message passing interface standard. Based on extensive tests, default settings for the supported QM packages are provided such that energy is conserved for typical QM/MM MD simulations in the microcanonical ensemble. Results for the free energy of binding of calcium ions to aspartate in aqueous solution comparing semiempirical and density functional Hamiltonians are shown to demonstrate features of this interface. © 2013 Wiley Periodicals, Inc.     

Dynamic properties of water around a protein-DNA complex from molecular dynamics simulations

Formation of protein-DNA complex is an important step in regulation of genes in living organisms. One important issue in this problem is the role played by water in mediating the protein-DNA interactions. In this work, we have carried out atomistic molecular dynamics simulations to explore the heterogeneous dynamics of water molecules present in different regions around a complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. It is demonstrated that such heterogeneous water motions around the complex are correlated with the relaxation time scales of hydrogen bonds formed by those water molecules with the protein and DNA. The calculations reveal the existence of a fraction of extraordinarily restricted water molecules forming a highly rigid thin layer in between the binding motifs of the protein and DNA. It is further proved that higher rigidity of water layers around the complex originates from more frequent reformations of broken water-water hydrogen bonds. Importantly, it is found that the formation of the complex affects the transverse and longitudinal degrees of freedom of surrounding water molecules in a nonuniform manner.     

Dissociative water potential for molecular dynamics simulations

A new interatomic potential for dissociative water was developed for use in molecular dynamics simulations. The simulations use a multibody potential, with both pair and three-body terms, and the Wolf summation method for the long-range Coulomb interactions. A major feature in the potential is the change in the short-range O-H repulsive interaction as a function of temperature and/or pressure in order to reproduce the density-temperature curve between 273 K and 373 at 1 atm, as well as high-pressure data at various temperatures. Using only the change in this one parameter, the simulations also reproduce room-temperature properties of water, such as the structure, cohesive energy, diffusion constant, and vibrational spectrum, as well as the liquid-vapor coexistence curve. Although the water molecules could dissociate, no dissociation is observed at room temperature. However, behavior of the hydronium ion was studied by introduction of an extra H+ into a cluster of water molecules. Both Eigen and Zundel configurations, as well as more complex configurations, are observed in the migration of the hydronium.     

Influence of electrostatic interactions on the properties of cyanobiphenyl liquid crystals predicted from atomistic molecular dynamics simulations

The influence of force field details in all-atom molecular dynamics (MD) simulations on the predicted thermodynamic, structural, and dynamic properties of bulk 4-cyano-4?-pentylbiphenyl (5CB) systems have been investigated in the 292–368 K temperature range. The effect of the molecular dipole moment and the details of dihedral potential for biphenyl unit were investigated using both polarisable (POL) and non-polarisable (NP) versions of the quantum chemistry-based force field. The predicted densities for the nematic and isotropic phases of bulk 5CB were found to be in excellent agreement with available experimental data. The nematic-isotropic transition temperature (T_NI) showed strong sensitivity to the force field details, MD simulations with partial atomic charge distributions and molecular dipole moment corresponding to high-level quantum chemistry calculations predicted an overestimation of the T_NI by about 30 K. Rescaling the charges to allow the molecular dipole to be closer to experimentally reported values of 5CB dipole in condensed phases, significantly improved the prediction of T_NI as well as other thermodynamic and dynamic properties of 5CB. We also discuss how the structural, thermodynamic, and dynamic properties of bulk 5CB are affected by the flexibility of the central biphenyl dihedral and the inclusion of induced polarisation effects.     

Enthalpies of mixing predicted using molecular dynamics simulations and OPLS force field

Enthalpy of mixing (EOM) is one of the most basic thermodynamic properties of mixtures. To assess feasibility of predicting EOM using force field simulation methods, fifteen (15) representative binary mixtures were investigated using MD simulations based on OPLS and TIP4P force fields. The simulation conditions and errors were carefully examined. The precision level of 0.04 kJ/mol was obtained for calculated EOM data. However, the predictions, measured by deviations from experimental data, were only qualitatively correct. Among various factors influencing the accuracy of predictions, force field quality representing interactions among different molecules plays the most significant role. Using methanol/benzene and ethanol/benzene as examples, we demonstrated that non-additive interaction terms between polarizable atoms can be used to significantly improve the quality of predictions. In addition, it appears that charge-dependent LJ parameters are required in order to represent the polarization effects accurately.     

Properties of inhibitors of methane hydrate formation via molecular dynamics simulations

Within the framework of a proposed two-step mechanism for hydrate inhibition, the energy of binding of four inhibitor molecules (PEO, PVP, PVCap, and VIMA) to a hydrate surface is estimated with molecular dynamic simulations. One key feature of this proposed mechanism is that the binding of an inhibitor molecule to the surface of an ensuing hydrate crystal disrupts growth and therein crystallization. It is found through the molecular dynamic simulations that inhibitor molecules that experimentally exhibit better inhibition strength also have higher free energies of binding, an indirect confirmation of our proposed mechanism. Inhibitors increasing in effectiveness, PEO < PVP < PVCap < VIMA, have increasingly negative (exothermic) binding energies of -0.2 < -20.6 < -37.5 < -45.8 kcal/mol and binding free energies of increasing favorability (+0.4 approximately = +0.5 < -9.4 < -15.1 kcal/mol). Furthermore, the effect of an inhibitor molecule on the local liquid water structure under hydrate-forming conditions was examined and correlated to the experimental effectiveness of the inhibitors. Two molecular characteristics that lead to strongly binding inhibitors were found: (1) a charge distribution on the edge of the inhibitor that mimics the charge separation in the water molecules on the surface of the hydrate and (2) the congruence of the size of the inhibitor with respect to the available space at the hydrate-surface binding site. Equipped with this molecular-level understanding of the process of hydrate inhibition via low-dosage kinetic hydrate inhibitors we can design new, more effective inhibitor molecules.     

Ejection of solvated ions from electrosprayed methanol/water nanodroplets studied by molecular dynamics simulations

The ejection of solvated small ions from nanometer-sized droplets plays a central role during electrospray ionization (ESI). Molecular dynamics (MD) simulations can provide insights into the nanodroplet behavior. Earlier MD studies have largely focused on aqueous systems, whereas most practical ESI applications involve the use of organic cosolvents. We conduct simulations on mixed water/methanol droplets that carry excess NH(4)(+) ions. Methanol is found to compromise the H-bonding network, resulting in greatly increased rates of ion ejection and solvent evaporation. Considerable differences in the water and methanol escape rates cause time-dependent changes in droplet composition. Segregation occurs at low methanol concentration, such that layered droplets with a methanol-enriched periphery are formed. This phenomenon will enhance the partitioning of analyte molecules, with possible implications for their ESI efficiencies. Solvated ions are ejected from the tip of surface protrusions. Solvent bridging prior to ion secession is more extensive for methanol/water droplets than for purely aqueous systems. The ejection of solvated NH(4)(+) is visualized as diffusion-mediated escape from a metastable basin. The process involves thermally activated crossing of a ~30 kJ mol(-1) free energy barrier, in close agreement with the predictions of the classical ion evaporation model.     

Investigation of interfacial and structural properties of CTAB at the oil/water interface using dissipative particle dynamics simulations

We have used dissipative particle dynamics (DPD) to simulate the system of cetyltrimethylammonium bromide (CTAB) monolayer at the oil/water interface. The interfacial properties (interfacial density, interfacial thickness, and interfacial tension), structural properties (area compressibility modulus, end to end distance, and order parameter), and their dependence on the oil/water ratio and the surfactant concentration were investigated. Three different microstructures, spherical oil in water (o/w), interfacial phase, and water in oil (w/o), can be clearly observed with the oil/water ratio increasing. Both the snapshots and the density profiles of the simulation show that a well defined interface exists between the oil and water phases. The interface thickens with CTAB concentration and oil/water ratio. The area compressibility modulus decreases with an increase in the oil/water ratio. The CTAB molecules are more highly packed at the interface and more upright with both concentration and oil/water ratio. The root mean square end-to-end distance and order parameter have a very weak dependence on the oil/water ratio. But both of them show an increase with CTAB concentration, indicating that the surfactant molecules at the interface become more stretched and more ordered at high concentration. As CTAB concentration increases further, the order parameter decreases instead because the bending of the interface. At the same time, it is shown that CTAB has a high interfacial efficiency at the oil/water interface.     

Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

Molecular dynamics simulations of amphiphilic graft copolymer molecules at a water/air interface

Fully atomistic molecular dynamics simulations of amphiphilic graft copolymer molecules have been performed at a range of surface concentrations at a water/air interface. These simulations are compared to experimental results from a corresponding system over a similar range of surface concentrations. Neutron reflectivity data calculated from the simulation trajectories agrees well with experimentally acquired profiles. In particular, excellent agreement in neutron reflectivity is found for lower surface concentration simulations. A simulation of a poly(ethylene oxide) (PEO) chain in aqueous solution has also been performed. This simulation allows the conformational behavior of the free PEO chain and those tethered to the interface in the previous simulations to be compared.     

Ab initio molecular dynamics simulations of the liquid/vapor interface of sulfuric acid solutions

Ab initio molecular dynamics simulations of the liquid-vapor interface are presented for thin slabs of 72 water molecules containing a single molecule of sulfuric acid. Trajectories in the 306-330 K range are calculated for two functionals with double- and triple-ζ quality basis sets. Comparisons are made between BLYP and HCTH/120 results for the slab simulations and for bulk simulations of one H(2)SO(4) in a periodic box with 63 waters. Good agreement is found with the available experimental data and the results of other relevant AIMD studies with respect to ionization of the acid, size of the coordination shells, partitioning of the ions with the hydronium exhibiting a surface preference and the anions in the interior, and the orientational distributions for the hydronium ions and for the surface/subsurface water molecules. The major differences in the performance of the two functionals are attributable to the greater basicity of the anion oxygen atoms with the HCTH functional and the more structured aqueous solution with BLYP. The enhanced basicity results in larger aqueous coordination shells for the anion oxygens. The structuring of the BLYP aqueous solution is observed in the corrugation of the water density profile, the higher first peak in g(OO)(r), and a smaller water self-diffusion constant. This structuring with the BLYP functional yields anion hydrogen bonds that endure longer and where the dissociated ions more rapidly and directly segregate in the slab. The simulations indicate that aqueous surfaces containing ionizable diprotic acids can be modeled with rather modest sized systems and be informative.     

Kinetics of NaCl nucleation in supercritical water investigated by molecular dynamics simulations

At system pressures between 17 MPa and 25 MPa the nucleation and growth of NaCl nanoparticles in water at supercritical conditions was investigated by molecular dynamics simulations at different system temperatures and system densities. Our results show that particle formation takes place within a few hundred picoseconds after the jump from ambient to supercritical conditions. After nucleation a phase of growth by adding monomers is followed by growth via cluster-cluster collisions. We present results on the time development of distributions of cluster sizes, cluster compositions, and cluster temperatures as well as radial distribution functions and nucleation rates.     

Structure and dynamics of water at water|Pt interface as seen by molecular dynamics computer simulation

Five molecular dynamics computer simulations were performed to study the structural and dynamical properties of water next to uncharged and charged Pt surfaces. The results show that the structure of a water layer adsorbed on the metal surface is very sensitive to the details of the water–metal potential. While patches of short-living hexagonal ice-like structure are observed in the adsorbed water layer next to the uncharged Pt(111) surface, a square lattice solid-like structure is seen for the layer on top of the uncharged Pt(100) surface. The orientational ordering for the following two layers of water next to uncharged Pt is displaying a preference towards the orientations that are characteristic of hexagonal ice-I, while water is liquid-like in these layers. In the presence of a high value external electric field water reorients and undergoes a layering transition.     

Partition coefficients of organic molecules in squalane and water/ethanol mixtures by molecular dynamics simulations

Accurate partition coefficient data of migrants between a polymer and a solvent are of paramount importance for estimating the migration of the migrant over time, including the concentration of the migrant at infinite time in the two solvents. In this article it is shown how this partition coefficient can be estimated for both a small hydrophilic and a hydrophobic organic molecules between squalane (used here to mimic low density poly ethylene) and water/ethanol solutes using thermodynamic integration to calculate the free energy of solvation. Molecular dynamics simulations are performed, using the GROMACS software, by slowly decoupling of firstly the electrostatic and then the Lennard-Jones interactions between molecules in the simulation box. These calculations depend very much on the choice of force field. Two force fields have been tested in this work, the TraPPE-UA (united-atom) and the OPLS-AA (all-atom). The computational cheaper TraPPE-UA force field showed to be more accurate over the whole range of systems compared to the OPLS-AA force field. Moreover, some of the calculations were done with five different water models to investigate the influence of the specific water model on the calculations. It was found that the combination of the TraPPE-UA force field and the TIP4p water model gave the best results. Based on the methodology proposed in this article, it is possible to obtain good partition coefficients only knowing the chemical structure of the molecules in the system.     

Electronic properties of liquid water by sequential molecular dynamics/density functional theory

Electronic properties of liquid water were analysed by a sequential molecular dynamics (MD)/density functional theory approach. MD simulations are based on a polarisable model for water. Emphasis was placed on the prediction of the water dipole moment, liquid state polarisability, ionisation potential (IP), and vertical electron affinity. The dipole moment of the water molecule in liquid water is not dependent on the number of molecules included in the quantum mechanical calculations. The polarisability of the water molecule in liquid water is 4% lower than its gas phase value. The IP of liquid water (9.7 ± 0.06 eV) is in good agreement with recent experimental data.