On the Quest of Dioxygen by Monomeric Sarcosine Oxidase. A Molecular Dynamics Investigation

It is reported here on random acceleration molecular dynamics (RAMD) simulations with the 2GF3 bacterial monomeric sarcosine oxidase (MSOX), O₂, and furoic acid in place of sarcosine, solvated by TIP3 H₂O in a periodic box. An external tiny force, acting randomly on O₂, accelerated its relocation, from the center of activation between residue K265 and the si face of the flavin ring of the flavin adenine dinucleotide cofactor, to the surrounding solvent. Only three of the four O₂ gates previously described for this system along a composite method technique were identified, while two more major O₂ gates were found. The RAMD simulations also revealed that the same gate can be reached by O₂ along different pathways, often involving traps for O₂. Both the residence time of O₂ in the traps, and the total trajectory time for O₂ getting to the solvent, could be evaluated. The new quick pathways discovered here suggest that O₂ exploits all nearby interstices created by the thermal fluctuations of the protein, not having necessarily to look for the permanent large channel used for uptake of the FADH cofactor. To this regard, MSOX resembles closely KijD3 N‐oxygenase. These observations solicit experimental substantiation, in a long term aim at discovering whether gates and pathways for the small gaseous ligands inside the proteins are under Darwinian functional evolution or merely stochastic control operates.     

A Leap-frog Algorithm for Stochastic Dynamics

Abstract

A third-order algorithm for stochastic dynamics (SD) simulations is proposed, identical to the powerful molecular dynamics leap-frog algorithm in the limit of infinitely small friction coefficient γ. It belongs to the class of SD algorithms, in which the integration time step Δt is not limited by the condition Δt ≤ γ^?1, but only by the properties of the systematic force. It is shown how constraints, such as bond length or bond angle constraints, can be incorporated in the computational scheme. It is argued that the third-order Verlet-type SD algorithm proposed earlier may be simplified without loosing its third-order accuracy. The leap-frog SD algorithm is proven to be equivalent to the verlet-type SD algorithm. Both these SD algorithms are slightly more economical on computer storage than the Beeman-type SD algorithm.     

MD Simulations of Liquids with Td,Oh Molecular Symmetry

Abstract

A new method is proposed for the calculation of intermolecular interactions in Molecular Dynamics simulations of liquids with T_d, O_h molecular symmetry. The new algorithm is based on the separation of the pair potential into a short-range and a long-range contribution described by a site-site and a spherical centre-centre potential model respectively using an additional cutoff distance. Test calculations for the Lennard-Jones fluids CCl₄ and SF₆ show significant savings in CPU time. We compare thermodynamic properties, pair correlation functions and a few dynamic autocorrelation functions obtained with the novel strategy with results of the commonly used algorithm for systems containing 864 molecules. Since no significant differences appear the new algorithm may be suggested as a useful contribution to the area of Molecular Dynamics simulation of liquids with these rather high molecular symmetries.     

On the Atomic Velocities in Molecular and Langevin Dynamics Simulations of Soft-Sphere Systems

Abstract

Time dependent probability distributions of the changes of direction of atomic velocities are considered in order to examine in detail the shape of the trajectories obtained through molecular simulations. We have analysed the atomic motions obtained from molecular dynamics simulations of soft-sphere systems at three very different states, i.e. a dilute fluid, a liquid at high density, and a solid. The methodology has also been used to check the reliability of the velocity evolution obtained when it is assumed that a single particle obeys the generalized Langevin equation and the effect of the other particles is represented by friction and random forces.     

Polymers of intrinsic microporosity for gas permeation: a molecular simulation study

We report a molecular simulation study for gas permeation in two membranes constructed from polymers of intrinsic microporosity (PIM-1 and PIM-7). With rigid ladder polymer chains, the membranes posses approximately 47.7 and 46.6% fractional free volumes (FFVs) in PIM-1 and PIM-7, respectively. The voids in the membranes have a diameter up to 9 Å and are largely interconnected. The sorption and diffusion of four gases (H₂, O₂, CH₄ and CO₂) were calculated by Monte Carlo and molecular dynamics simulations. The solubility coefficients increase in the order of H₂ < O₂ < CH₄ < CO₂, while the diffusion coefficients increase in the following order: CH₄ < CO₂ < O₂ < H₂. The simulation results agree well with experimental data, particularly for the solubility coefficients. The solubility and diffusion coefficients correlate well separately with the critical temperatures and effective diameters of gases. These molecular-based correlations can be used in the prediction for other gases. As attributed to the microporous structure, PIM-1 and PIM-7 outperform most glassy polymeric membranes in sorption and diffusion. PIM-1 has larger solubility and diffusion coefficients than PIM-7 because the cyano groups in PIM-1 lead to a stronger affinity and a larger FFV. The simulated solubility, diffusivity and permeation selectivities of CO₂/H₂, CO₂/O₂ and CO₂/CH₄ are consistent with experimental data. The quantitative microscopic understanding of gas permeation in the PIM membranes is useful for the new development of high-performance membranes.     

Molecular Dynamics Study of BaB2O4 in Crystalline and Molten States

Abstract

Structural aspects of BaB₂O₄ liquids have been investigated by the molecular dynamics simulation including the determination on the parameters of the interatomic potential applicable to BaB₂O₄ in both crystalline and molten states. The structure and physical properties of BaB₂O₄ crystals were successfully reproduced by the MD simulation for both α and β phases. The simulated interference function in the liquid state was also in good agreement with the experimental one. Several interesting features on the relaxation phenomena just after melting were reproduced by the simulation that the structure factors of simulated liquid maintain the characteristic features of the original crystal structure for more than 40ps after melting, and the variation of the number of rings formed by B-O bondings was found to increase after melting.     

Molecular Dynamics Simulation of CaF2 Structure in Crystalline,Superionic and Molten Phases

Abstract

We describe here a number of molecular dynamics simulations on calcium fluoride over a range of temperatures spanning the transitions to the superionic and molten state. The simulation temperatures are 1400, 1590, 1800, 2000, 2200 K. By using the bond spherical harmonics method with equal neighbor number, we have studied the structure and bond orientation of cation sublattice and anion sublattice in superionic conductor CaF₂. The bond order parameters Q₁ have been calculated both for the computer generated instataneous configurations from the simulation system and for the standard configurations from the normal distribution model of bond orientation. The comparison of Q₁ between the molecular dynamics simulation and the normal distribution model shows that not only the cation sublattice but also the anion sublattice can be described by the normal distribution model. The cations keep their original fcc frame, but in the anion case there is a great deal of random distortions from the original anion sublattice.     

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.     

Unveiling the Pathways of Dioxygen Through the C2 Component of the Environmentally Relevant Monooxygenase p‐Hydroxyphenylacetate Hydroxylase from Acinetobacter baumannii: A Molecular Dynamics Investigation

In this work, models of the homotetrameric C2 component of the monooxygenase p‐hydroxyphenylacetate hydroxylase from Acinetobacter baumannii, in complex with dioxygen (O₂) and, or not, the substrate p‐hydroxyphenylacetate (HPA) were built. Both models proved to be amenable to random‐acceleration molecular dynamics (RAMD) simulations, whereby a tiny randomly oriented external force, acting on O₂ at the active site in front of flavin mononucleotide (FMNH^?), accelerated displacement of O₂ toward the bulk solvent. This allowed us to carry out a sufficiently large number of RAMD simulations to be of statistical significance. The two systems behaved very similarly under RAMD, except for O₂ leaving the active site more easily in the absence of HPA, but then finding similar obstacles in getting to the gate as when the active site was sheltered by HPA. This challenges previous conclusions that HPA can only reach the active center after that the C4aOOH derivative of FMNH^? is formed, requiring uptake of O₂ at the active site before HPA. According to these RAMD simulations, O₂ could well get to FMNH^? also in the presence of the substrate at the active site.     

The Influence of Charge Calculation on Molecular Dynamics Simulation of Adenine in Water

Abstract

Molecular dynamics simulations of an aqueous solution of adenine have been performed using different methods of charge calculation to evaluate the influence of the values of the atomic charges on the dynamical results and to incorporate new information about the interaction between adenine and water. Four sets of partial charges where computed using ab-initio methods. In all cases the hydration properties of adenine were similar. These results support the view that the simulations by molecular dynamics, at least for the regime of infinite dilution, are not sensitive with respect to the different sets of partial charges used. A net hydrophobic behavior of the adenine molecule, on the water was observed.     

Molecular Dynamics Simulations of Carbon Nanotube Rolling and Sliding on Graphite

Abstract

Molecular dynamics simulations were carried out to investigate the origin of friction for carbon nanotubes on graphite substrates. In an initial simulation, a (10,10) nanotube was placed in an ‘in-registry’ starting position where the hexagonal lattice of the substrate matched that of the nanotube. In a second simulation, the substrate was oriented 90 degrees to the nanotube. A uniform force was applied to the nanotubes for 500 fs to set them into motion. The simulation was then run until the nanotubes stopped moving relative to the substrate. Only sliding was observed in the out-of-registry simulation, while periodic sliding and rolling was observed in the in-registry simulation. The latter is a result of the relatively larger surface corrugation for the in-registry case and occurs to avoid direct atomic collisions between nanotube and substrate atoms as the nanotube is moved along the substrate. Analysis of the kinetic energy suggests that the transition between sliding and rolling contributes to enhanced energy dissipation and higher net friction. These results are consistent with preliminary experimental observations by Superfine and coworkers.     

Binding of Net-Fla,a Netropsin-Flavin Hybrid Molecule,to DNA: Molecular Mechanics and Dynamics Studies in vacuo and in Water Solution

Abstract

We have studied the binding of the hybrid netropsin-flavin (Net-Fla) molecule onto four sequences containing four A.T base pairs. Molecular mechanics minimizations in vacuo show numerous minimal conformations separated by one base pair. 400 ps molecular dynamics simulations in vacuo have been performed using the lowest minima as the starting conformations. During these simulations, the flavin moiety of the drug makes two hydrogen bonds with an amino group of a neighboring guanine. A 200 ps molecular dynamics simulation in explicit water solution suggests that the binding of Net-Fla upon the DNA substrate is enhanced by water bridges. A water molecule bridging the amidinium of Net-Fla to the N3 atom of an adenine seems to be stuck in the dmg-DNA complex during the whole simulation. The fluctuations of the DNA helical parameters and of the torsion angles of the sugar-phosphate backbone are very similar in the simulations in vacuo and in water. The time auto-correlation functions for the DNA helical parameters decrease rapidly in the picosecond range in vacuo. The same functions computed from the water solution molecular dynamics simulations seem to have two modes: the rapid mode is similar to the behavior in vacuo, and is followed by a slower mode in the 10 ps range.     

On the Permeation by Dioxygen of Urate Oxidase from Aspergillus flavus in Complex with Xanthine Anion: Dioxygen Pathways and a Portrait of the Enzyme Cavities from Molecular Dynamics Simulations in Water Solution

This work describes molecular dynamics (MD) simulations in aqueous media for the complex of the homotetrameric urate oxidase (UOX) from Aspergillus flavus with xanthine anion ( 5 ) in the presence of dioxygen (O₂). After 196.6 ns of trajectory from unrestrained MD, a O₂ molecule was observed leaving the bulk solvent to penetrate the enzyme between two subunits, A/C. From here, the same O₂ molecule was observed migrating, across subunit C, to the hydrophobic cavity that shares residue V227 with the active site. The latter was finally attained, after 378.3 ns of trajectory, with O₂ at a bonding distance from 5 . The reverse same O₂ pathway, from 5 to the bulk solvent, was observed as preferred pathway under random acceleration MD (RAMD), where an external, randomly oriented force was acting on O₂. Both MD and RAMD simulations revealed several cavities populated by O₂ during its migration from the bulk solvent to the active site or backwards. Paying attention to the last hydrophobic cavity that apparently serves as O₂ reservoir for the active site, it was noticed that its volume undergoes ample fluctuations during the MD simulation, as expected from the thermal motion of a flexible protein, independently from the particular subunit and no matter whether the cavity is filled or not by O₂.     

Gibbs-Ensemble Molecular Dynamics: Liquid-Gas Equilibria for Lennard-Jones Spheres and n-Hexane

Abstract

We present a novel method to simulate phase equilibria in atomic and molecular systems. The method is a Molecular Dynamics version of the Gibbs-Ensemble Monte Carlo technique, which has been developed some years ago for the direct simulation of phase equilibria in fluid systems. The idea is to have two separate simulation boxes, which can exchange particles (or molecules) in a thermodynamically consistent fashion. Here we pres the derivation of the generalized equations of motion and discuss the relation of the resulting trajectory averages to the relevant ensemble. We test this Gibbs-Ensemble Molecular Dynamics algorithm by applying it to an atomic and a molecular system, i.e. to the liquid-gas coexistence in a Lennard-Jones fluid and in n-hexane. In both cases our results are in good accord with previous mean field and Gibbs-Ensemble Monte Carlo results as well as with the experimental data in the case of hexane. We also show that our Gibbs-Ensemble Molecular Dynamics algorithm like other Molecular Dynamics techniques can be used to study the dynamics of the system. Self-diffusion coefficients calculated with this method are in agreement with the result of conventional constant temperature Molecular Dynamics.     

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.     

On Dioxygen Permeation of MaL Laccase from the Thermophilic Fungus Melanocarpus albomyces: An all‐Atom Molecular Dynamics Investigation

In this article, biased molecular dynamics (MD) simulations of O₂ egress from the active center of MaL laccase toward the bulk solvent were described. Parameterization of the set of four Cu(II) ions, in the framework of CHARMM‐36 FF, was carried out on a recent dummy‐atom model that takes into account the Jahn–Teller effect. By carrying out a number of statistically relevant MD simulations, under a tiny randomly oriented external force applied to the molecule O₂, three preferred gates for O₂ egress from the enzyme and a few intermediate binding pockets (BPs) for O₂ were visible; all the gates and pockets were located on two of the three domains. This wide distribution of preferred gates notwithstanding the molecule O₂ was seen to follow specific pathways, exploiting consistently the interstices created by the enzyme thermal fluctuations. These are features that can be imagined to have evolved to make MaL laccase extremely efficient as catalysts in various reactions that require O₂.     

Diffusion of N2, O2, H2S and SO2 in MFI and 4A zeolites by molecular dynamics simulations

This paper presents diffusion data of N₂, O₂, H₂S and SO₂ in MFI and 4A zeolites obtained by molecular dynamics simulations, especially its dependence on temperature and loading. At high loadings and temperatures, the order of self-diffusivity of guests in two zeolites is O₂ > N₂ > H₂S>SO₂. The diffusion behaviour is different in different zeolites at lower loadings, reflecting different influences from straight channels (MFI) and α-cages (4A). Furthermore, with increasing loading, the self-diffusivity of guest molecules decreases in MFI but generally increases in 4A. The centre of mass (COM) probability densities and diffusion trajectories of guests give insight into molecular-level diffusion process. The simulation results reveal that with increasing loading, the diffusion mechanism would change from the inter-pore to intra-pore diffusion in MFI. However, in 4A, the intra-pore diffusion is predominant at low and high loadings, but inter-pore diffusion is more important at moderate loadings.