A Molecular Dynamics Study for Lithium Metasilicate: Liquid and Quenched Supercooled States

A molecular dynamics (MD) simulation of the lithium metasilicate (Li₂SiO₃) system has been performed to study the dynamics and geometrical structure of a liquid state at 1673 K and a glassy state at 700 K. The long term dynamics of the glassy state is shown for the self-part of the van Hove correlation functions. In the glassy state, diffusion of lithium ion was found to occur mainly through jump motions among equilibrium sites surrounded by oxygen atoms. A geometric structural analysis with polyhedra made of oxygen atoms around lithium was applied to characterize the structure and dynamical behavior.     

Effect of Coordinated Ions on Structure and Flexibility of Parallel G-quadruplexes: A Molecular Dynamics Study

Abstract

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G₇) quadruplex structures with Na⁺ or K⁺ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Å from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na⁺ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na⁺ ions, within the quadruplex cavity, are more mobile than coordinated K⁺ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.     

Molecular Dynamics Simulation of Self Diffusion in Microclusters

Abstract

Molecular dynamics simulation is carried out to study the mechanism of self diffusion which is characteristic of solid-like microclusters. A two-dimensional system with the Lennard-Jones potential is employed and the temperatures near the triple point of the two-dimensional bulk system are adopted for the simulation. The results show that: a) microclusters consist of two regions, i.e., solid-like cores and liquid-like surface regions, b) the size dependence of the diffusion coefficient for microclusters is weak, and the value of the solid-like core region is not much different from that of the bulk liquid, c) the activation energy of diffusion for microclusters is twenty to thirty times larger than that for the bulk liquid, d) the diffusion mechanism in the solid-like region involves the collective motion of small domains containing ten to twenty atoms which results in the formation of low density regions, sometimes even vacancy clusters, between them, and atoms in the low density regions change their positions to cause diffusion.     

Inhibition of Ice Nucleus Growth in Water by Alanine Dipeptide

A molecular dynamics simulation has been carried out for the mixture of an ice nucleus, supercooled water and a molecule of alanine dipeptide (AD). The dipeptide molecule has been allocated near the nucleus surface which corresponds to the prism plane of ice crystal. The molecule is found to approach the ice surface so that the two hydrophilic sites on one side of the molecule (Oc2 and Hn1) are closest to the surface. The hydrogen bond between Hn1 site and the oxygen atom on the prism plane of the ice nucleus is expected. The perturbations of two hydrophilic sites (Oc1 and Hn2), which are surrounded by hydrophobic sites and are pointing away from the surface, attenuate the approach of water molecules to these sites. Thus, these water molecules diffuse. The hydrogen bond between the oxygen atoms on the prism plane and the hydrogen atoms of water molecules is attenuated by the diffusion.     

On the Approximation of Solvent Effects on the Conformation and Dynamics of Cyclosporin A by Stochastic Dynamics Simulation Techniques

Abstract

The molecular simulation technique of stochastic dynamics (SD) is tested by application to the immunosuppressive drug cyclosporin A (CPA). Two stochastic dynamics simulations are performed, one (SD_CCl4 ) with atomic friction coefficients proportional to the viscosity of the nonpolar solvent CCl₄, and one (SD_H2O) with atomic friction coefficients corresponding to an aqueous solution. The atomic friction coefficients are also taken proportional to an approximate expression for the atomic accessible surface area. The properties of both stochastic dynamics simulations are compared to those of two full molecular dynamics (MD) simulations of cyclosporin A, one in a box with 591 CCl₄ molecules, and one in a box with 632 H₂O molecules.

The properties of cyclosporin A as found in the molecular dynamics simulation in CCl₄ are well reproduced by the SD_CCl4 simulation. This indicates that the neglect of a mean force reresenting the average solvent effects on the solute is justified in the case of nonpolar solvents. For polar solvents, like water, this mean force may not be neglected. The SD_H2O simulation of cyclosporin A clearly fails to reproduce the amount of hydrogen bonding found in the molecular dynamics stimulation of cyclosporin A in water.

A comparison with a molecular dynamics simulation of cyclosporin A in vacuo shows that both the SD_CCl4 and the SD_H2O simulation come closer to the properties of the molecular dynamics simulations in CCl₄ and in H₂O than a molecular dynamics simulation in vacuo.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

An Interatomic Potential Model for H2O: Applications to Water and Ice Polymorphs

Abstract

A new interatomic potential model for H₂O which consists of 2-body central (O-H, O-O and H-H) and 3-body teams and does not contain artificial constraints on the motions of oxygen and hydrogen atoms is proposed. The interatomic potential function parameters were determined empirically so as to reproduce the fundamental and essential features of water and ice Ih using molecular dynamics (MD) methods.

We carried out the MD simulations of water, and find the Ice Ih, Ice II and Ice IX using this potential model, in structures and physical properties are in agreement with the experimental results except for the compressibilities of both water and ice Ih. We expect that, by refining this model, we can apply this model to problems involving the reaction of water molecules with other components.     

Combination of Neutron Scattering and Molecular Dynamics to Determine Internal Motions in Biomolecules

Abstract

The forms and frequencies of atomic dynamics on the pico- and nanosecond timescales are accessible experimentally using incoherent neutron scattering. Molecular dynamics simulations cover the same space and time domains and neutron scattering intensities can be calculated from the simulations for direct comparison with experiment. To illustrate the complementarity of neutron scattering and molecular dynamics we examine measured and simulation-derived elastic incoherent scattering profiles from myoglobin and from the crystalline alanine dipeptide. Elastic incoherent scattering gives information on the geometry of the volume accessible to the atoms in the samples. The simulation-derived dipeptide elastic scattering profiles are in reasonable accord with experiment, deviations being due to the sampling limitations in the simulations and experimental detector normalisation procedures. The simulated dynamics is decomposed, revealing characteristic profiles due to rotational diffusional and translational vibrational motions of the methyl groups. In myoglobin, for which the timescale of the simulation matches more closely that accessible to the experiment, good agreement is seen for the elastic incoherent structure factor. This indicates that the space sampled by the hydrogen atoms in the protein on the timescale <100 ps is well represented by the simulation. Part of the helix atom fluctuations can be described in terms of rigid helix motions.     

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

A molecular dynamics simulation of the flavin mononucleotide–RNA aptamer complex

We report on an unrestrained molecular dynamics simulation of the flavin mononucleotide (FMN)–RNA aptamer. The simulated average structure maintains both cross‐strand and intermolecular FMN–RNA nuclear Overhauser effects from the nmr experiments and has all qualitative features of the nmr structure including the G10–U12–A25 base triple and the A13–G24, A8–G28, and G9–G27 mismatches. However, the relative orientation of the hairpin loop to the remaining part of the molecule differs from the nmr structure. The simulation predicts that the flexible phosphoglycerol part of FMN moves toward G27 and forms hydrogen bonds. There are structurally long‐lived water molecules in the FMN binding pocket forming hydrogen bonds within FMN and between FMN and RNA. In addition, long‐lived water is found bridging primarily RNA backbone atoms. A general feature of the environment of long‐lived “structural” water is at least two and in most cases three or four potential acceptor atoms. The 2′‐OH group of RNA usually acts as an acceptor in interactions with the solvent. There are almost no intrastrand O2′H(n)⋮O4′(n + 1) hydrogen bonds within the RNA backbone. In the standard case the preferred orientation of the 2′‐OH hydrogen atoms is approximately toward O3′ of the same nucleotide. However, a relatively large number of conformations with the backbone torsional angle γ in the trans orientation is found. A survey of all experimental RNA x‐ray structures shows that this backbone conformation occurs but is less frequent than found in the simulation. Experimental nmr RNA aptamer structures have a higher fraction of this conformation as compared to the x‐ray structures. The backbone conformation of nucleotide n + 1 with the torsional angle γ in the trans orientation leads to a relatively short distance between 2′‐OH(n) and O5′(n + 1), enabling hydrogen‐bond formation. In this case the preferred orientation of the 2′‐OH hydrogen atom is approximately toward O5′(n + 1). We find two relatively short and dynamically stable types of backbone–backbone next‐neighbor contacts, namely C2′(H)(n)⋮O4′(n + 1) and C5′(H)(n + 1)⋮O2′(n). These interactions may affect both backbone rigidity and thermodynamic stability of RNA helical structures. © 1999 John Wiley & Sons, Inc. Biopoly 50: 287–302, 1999     

Modeling of Platinum Clusters in H-Mordenite

The size, location and structure of Pt clusters in H-mordenite have been investigated by molecular mechanics energy minimization and molecular dynamics simulation techniques using the Catalysis software of Molecular Simulations (MSI). Lattice energy minimizations are performed to study the effects of the specific framework aluminum positions on the location and stability of monoatomic Pt sites in H-mordenite. The lattice energies relative to the siliceous platinum-aluminosilicate structure reveal that the stability of a single Pt atom in H-mordenite is remarkably influenced by the specific location of the Al atoms in the lattice. At the studied Si/Al ratio of two Al ions per unit cell, a stabilization of the H-mordenite lattice upon Pt deposition is obtained. Moreover, lattice energy calculations on Pt/aluminosilicate mordenites of different metal contents per unit cell have been performed. An optimum size for the aggregate confined to the 12-ring main channel that is almost independent of the Pt content per mordenite unit cell has been found. The structural features of the resulting clusters at the end of molecular dynamics simulations on Pt/alumina-mordenites reflect a strong metal-zeolite interaction. The present results are consistent with a previous molecular dynamics simulation study on the structure of platinum deposited on SiO₂ surfaces.     

Simulation of a Hard-Sphere Fluid in Bicontinuous Random Media

Abstract

The influence of solid-phase connectivity on size-exclusion partitioning and on diffusion of a dilute hard-sphere fluid in overlapping and nonoverlapping spheres models of porous media is investigated using molecular dynamics and Monte Carlo simulation techniques. Four models are examined, two of which are subject to constrained bicontinuity of the pore and solid phases and two in which the solid spheres in the assemblies are randomly distributed in space. It is shown that at moderate to high porosities, connected (bicontinuous) structures lead to a significant increase in the partition and diffusion coefficients when the particles of the pore fluid are of finite size. The consequences of solid phase connectivity are also clearly illustrated in the long-time decay of the velocity autocorrelation function (VACF) of the diffusing particles, particularly in the vicinity of the percolation threshold. Under these conditions the power law exponents on the long-time tail of the VACF are generally found to be higher in connected models than in random systems and the importance of this result is demonstrated using one of the scaling rules of percolation theory. The simulation results are also compared with the predictions of current theories of partitioning and diffusion in random sphere assemblies and, with reference to experimental data available from the literature, it is shown that bicontinuous models are better representations of real porous media.     

Structural Transformations of Ice at High Pressures Via Molecular Dynamics Simulations

Abstract

A simple classical model is used for the study of the structural transformations of ice under high pressures, such as ice VIII to VII and X, via classical molecular dynamics (MD) simulation. In the present MD simulation, pair potentials of a simple form between pair of atoms and a thee-body potential representing the H-O-H angle dependence, originally developed by Kawamura et al., were used. Starting with a stable ice VIII at low pressure and low temperature, we have carried out two different MD runs, one with increasing pressure keeping the temperature constant (simulation I) and the other with increasing temperature under constant pressure (simulation II). From these MD simulations we have obtained the structural transformations from ice VIII to VII for both simulations; the former was finally transformed into ice X for the simulation I. The present results are compatible with recent experiments on high pressure ices.     

Dynamical Properties of Condensed Charged Polymer Melts

Abstract

Dynamical properties of condensed charged polymer melts are studied with a two-dimensional model and molecular dynamics simulation. Screened Coulombic interactions are assumed for the interactions between the monomer charges as well as the counterions which were introduced to neutralize the total monomer charges of polymer chains. Through molecular dynamics calculations, we have obtained the radial distribution function and velocity auto-correlation function, and their density dependences. As a structural characteristics in condensed charged polymer melts, we find that the monomers tend to form triangular structures locally. The radial distribution function for the center monomers implies that the polymer chains on our two-dimensional model are not entangled. The diffusion properties for both the counter-ions and innermost monomers are studied in detail. In this paper, we have also computed neutral polymer melts to study the effects of the presence of the long-range Coulombic interactions on the dynamical properties of polymer melts. We find that the Coulombic interactions significantly reduce the self diffusion. Snapshots analyses obtained from the molecular dynamics simulation suggest that the reptation model is not a proper model for two-dimensional polymer-chain melts, in which there are little entanglement effects in agreement with the result of the analyses for the radial distribution function.     

Oxygen Adsorption on Zr(0001): An ab Initio Study

Abstract

The main goal of this paper is to study the oxygen adsorption on the Zr(0001) surface using first-principles total-energy calculations based on density-functional theory. We present preliminary results which concern the atomic oxygen adsorption on the Zr(0001) surface. We first report a static study where we calculate the atomic structure and adsorption energy for oxygen occupying various surface and subsurface sites at three different coverages: Θ = 1/4, 1/2, and 1 ML. We find that oxygen atoms are preferentially adsorbed into the octahedral holes between the 2nd and 3rd metallic layers. We secondly perform ab initio molecular dynamics calculations for the Zr(0001) – (3 × 3) – O system to show how the oxygen can penetrate through the surface and how it finally reaches its equilibrium position, trapped between the 1st and 2nd zirconium layers.     

Study on properties of liquid ammonia via molecular dynamics simulation based on ABEEMσπ polarisable force field

Abstract

The molecular dynamics simulation has been performed to investigate the charge distribution, structural and dynamical properties of liquid ammonia at 273 K using a polarisable force field of the atom-bond electronegativity equalisation method (ABEEMσπ). One ammonia molecule in this model has eight charge sites, one N atomic site, three H atomic sites, three N–H bond sites and one lone-pair electron site. ABEEMσπ model can present the quantitative site charges of molecular ammonias in liquid and their changing in response to their surroundings. The radial distribution functions and dynamical properties are in fair agreement with the available experimental data. The first peak of g_NN(r) appears at N–N distance of ~3.50 ± 0.05 Å where most hydrogen bonds are formed. The average coordination number of the first shell is 13.0 ± 0.1 among which a central ammonia molecule intimately connects 3 ~ 4 ammonia molecules by hydrogen bonds. The power spectrum shows the vibrations of hydrogen bonds. For a reference, a simple estimation of the average hydrogen bonding energy in liquid ammonia is 6.5 ± 0.1 kcal/mol larger than 3.8 ± 0.3 kcal/mol in dimer ammonia. Our simulation results provide more detailed information about liquid ammonia.     

Incoherent elastic and quasi-elastic neutron scattering investigation of hemoglobin dynamics

In this work we investigate the dynamic properties of hemoglobin in glycerolD(8)/D(2)O solution using incoherent elastic (ENS) and quasi-elastic (QENS) neutron scattering. Taking advantage of complementary energy resolutions of backscattering spectrometers at ILL (Grenoble), we explore motions in a large space-time window, up to 1 ns and 14 A; moreover, in order to cover the harmonic and anharmonic protein dynamics regimes, the elastic experiments have been performed over the wide temperature interval of 20-300 K. To study the dependence of the measured dynamics upon the protein quaternary structure, both deoxyhemoglobin (in T quaternary conformation) and carbonmonoxyhemoglobin (in R quaternary conformation) have been investigated. From the ENS data the mean square displacements of the non-exchangeable hydrogen atoms of the protein and their temperature dependence are obtained. In agreement with previous results on hydrated powders, a dynamical transition at about 220 K is detected. The results show interesting differences between the two hemoglobin quaternary conformations, the T-state protein appearing more rigid and performing faster motions than the R-state one; however, these differences involve motions occurring in the nanosecond time scale and are not detected when only faster atomic motions in the time scale up to 100 ps are investigated. The QENS results put in evidence a relevant Lorentzian quasi-elastic contribution. Analysis of the dependence of the Elastic Incoherent Structure Factor (EISF) and of the Lorentzian halfwidth upon the momentum transfer suggests that the above quasi-elastic contribution arises from the diffusion inside a confined space, values of confinement radius and local diffusion coefficient being compatible with motions of hydrogen atoms of the amino acid side chains. When averaged over the whole range of momentum transfer the QENS data put in evidence differences between deoxy and carbonmonoxy hemoglobin and confirm the quaternary structure dependence of the protein dynamics in the nanosecond time scale.     

Investigation of molecular dissolution mechanism of ketoprofen binary and ternary solid dispersions by molecular dynamics simulations

Abstract

The research aimed to investigate the molecular dissolution mechanism of both binary and ternary solid dispersions by molecular dynamics simulations. The simulation results indicated that the drug molecules were much easier to be released from surfactant-containing ternary systems than from binary ones. Moreover, sodium dodecyl sulfate as an additive in ternary systems had better effects than Tween 60. The simulation results were in well agreement with the experimental results. This research presented a reasonable explanation of molecular dissolution mechanism for both binary and ternary solid dispersions, which may benefit the future development of solid dispersion formulations.     

Transient hydrogen bonding in uniformly ¹³C,¹⁵N-labeled carbohydrates in water

We report NMR studies of transient hydrogen bonding in a polysaccharide (PS) dissolved in water without cosolvent at ambient temperature. The PS portion of the Escherichia coli O142 lipopolysaccharide is comprised of repeating pentasaccharide units of GalNAc (N-acetyl galactosamine), GlcNAc (N-acetyl glucosamine), and rhamnose in a 3:1:1 ratio, respectively. A 105-ns molecular dynamics (MD) simulation on one pentasaccharide repeat unit predicts transient inter-residue hydrogen bonds from the GalNAc NH groups in the PS. To investigate these predictions experimentally, the PS was uniformly ¹³C,¹⁵N enriched and the NH, carbonyl, C2, C4, and methyl resonances of the GalNAc and GlcNAc residues assigned using through-bond triple-resonance NMR experiments. Temperature dependence of amide NH chemical shifts and one-bond NH J couplings support that NH groups on two of the GalNAc residues are donors in transient hydrogen bonds. The remaining GalNAc and GlcNAc NHs do not appear to be donors from either temperature-dependent chemical shifts or one-bond NH J couplings. These results substantiate the presence of weak or partial hydrogen bonds in carbohydrates, and that MD simulations of repeating units in PSs provide insight into overall PS structure and dynamics. Published 2011 Wiley Periodicals, Inc. Biopolymers 97: 145–154, 2012.     

Molecular dynamics simulation on the decomposition of type SII hydrogen hydrate and the performance of tetrahydrofuran as a stabiliser

Molecular dynamics simulation is used to study the decomposition and stability of SII hydrogen and hydrogen/tetrahydrofuran (THF) hydrates at 150 K, 220 K and 100 bar. The modelling of the microscopic decomposition process of hydrogen hydrate indicates that the decomposition of hydrogen hydrate is led by the diffusive behaviour of H₂ molecules. The hydrogen/THF hydrate presents higher stability, by comparing the distributions of the tetrahedral angle of H₂O molecules, radial distribution functions of H₂O molecules and mean square displacements or diffusion coefficients of H₂O and H₂ molecules in hydrogen hydrate with those in hydrogen/THF hydrate. It is also found that the resistance of the diffusion behaviour of H₂O and H₂ molecules can be enhanced by encaging THF molecules in the (5¹²6⁴) cavities. Additionally, the motion of THF molecules is restricted due to its high interaction energy barrier. Accordingly, THF, as a stabiliser, is helpful in increasing the stability of hydrogen hydrate.