Hybrid MC/QC simulations of water-assisted proton transfer in nucleosides. Guanosine and its analog acyclovir

To provide an in-depth insight into the molecular basis of spontaneous tautomerism in DNA and RNA base pairs, a hybrid Monte Carlo (MC)–quantum chemical (QC) methodology is implemented to map two-dimensional potential energy surfaces along the reaction coordinates of solvent-assisted proton transfer processes in guanosine and its analog acyclovir in aqueous solution. The solvent effects were simulated by explicit inclusion of water molecules that model the relevant part of the first hydration shell around the solute. The position of these water molecules was estimated by carrying out a classical Metropolis Monte Carlo simulation of dilute water solutions of the guanosine (Gs) and acyclovir (ACV) and subsequently analyzing solute–solvent intermolecular interactions in the statistically-independent MC-generated configurations. The solvent-assisted proton transfer processes were further investigated using two different ab initio MP2 quantum chemical approaches. In the first one, potential energy surfaces of the ‘bare’ finite solute–solvent clusters containing Gs/ACV and four water molecules (MP2/6-31+G(d,p) level) were explored, while within the second approach, these clusters were embedded in ‘bulk’ solvent treated as polarizable continuum (C-PCM/MP2/6-31+G(d,p) level of theory). It was found that in the gas phase and in water solution, the most stable tautomer for guanosine and acyclovir is the 1H-2-amino-6-oxo form followed by the 2-amino-6-(sZ)-hydroxy form. The energy barriers of the water-assisted proton transfer reaction in guanosine and in acyclovir are found to be very similar – 11.74 kcal mol^?1 for guanosine and 11.16 kcal mol^?1 for acyclovir, and the respective rate constants (k = 1.5?×?10¹ s^?1, guanosine and k = 4.09?×?10¹ s^?1, acyclovir), are sufficiently large to generate the 2-amino-6-(sZ)-hydroxy tautomer. The analysis of the reaction profiles in both compounds shows that the proton transfer processes occur through the asynchronous concerted mechanism.     

Test of the Inverse Monte Carlo Method for the Calculation of Interatomic Potential Energies in Atomic Liquids

Abstract

The principle purpose of this paper is to demonstrate the use of the Inverse Monte Carlo technique for calculating pair interaction energies in monoatomic liquids from a given equilibrium property. This method is based on the mathematical relation between transition probability and pair potential given by the fundamental equation of the “importance sampling” Monte Carlo method. In order to have well defined conditions for the test of the Inverse Monte Carlo method a Metropolis Monte Carlo simulation of a Lennard Jones liquid is carried out to give the equilibrium pair correlation function determined by the assumed potential. Because an equilibrium configuration is prerequisite for an Inverse Monte Carlo simulation a model system is generated reproducing the pair correlation function, which has been calculated by the Metropolis Monte Carlo simulation and therefore representing the system in thermal equilibrium. This configuration is used to simulate virtual atom displacements. The resulting changes in atom distribution for each single simulation step are inserted in a set of non-linear equations defining the transition probability for the virtual change of configuration. The solution of the set of equations for pair interaction energies yields the Lennard Jones potential by which the equilibrium configuration has been determined.     

A Monte Carlo scheme based on mid-density in a hysteresis loop to determine equilibrium phase transition

A new and simple method to determine equilibrium phase transition in adsorption systems exhibiting a hysteresis loop is presented as an alternative to methods such as multiple histogram reweighting, gauge cell method and thermodynamic integration. This method is based on the NVT-grand canonical Monte Carlo mid-density scheme to determine the coexistence chemical potential and coexistence densities of an adsorption system. We illustrate this new scheme with argon and methane adsorption in a number of model solids having slit and cylindrical pores. This method does not have a strong basis on thermodynamic ground, but it does provide a simple heuristic approach that is simpler to understand physically.     

A Two-Ellipsoid Model Based Upon A Gaussian Overlap Potential

Abstract

The two-ellipsoid model (TEM) is proposed as a versatile single-site model which can be used in the study of liquid crystal phases. This TEM uses two ellipsoids to describe a molecule, one ellipsoid for the geometry and the other for the interaction strengths of the molecule. The present TEM can mimic asymmetric interactions of a liquid crystal molecule by separating the center of the interaction ellipsoid from that of the geometry ellipsoid. The potential energy surfaces of the present TEMs compare favorably with those of the corresponding Gay-Berne and the site–site models.

Monte Carlo simulations with 320 particles are performed for a symmetric interaction TEM and an asymmetric interaction TEM. The asymmetric interaction TEM displays a slightly higher transition temperature than the symmetric interaction TEM indicating that asymmetric interactions can be a driving force in a phase transition. Radial and cylindrical distribution functions of two models in the isotropic phase are similar, but those in the nematic phase are quite different.     

Monte Carlo Simulations of the Chemical Potential and Free Energy for Trimer and Hexamer Rings

Abstract

The chemical potential of a trimer and hexamer model ring system was determined by computer simulation over a range of temperatures and densities. Such ring molecules are important as model aromatic and naphthenic hydrocarbons. Thermodynamic integration of the pressure along a reversible path, Widom's ghost particle insertion method and Kirkwood's charging parameter method were used over a molecular density range of 0.05 to 0.30. Data were obtained by Monte Carlo simulation of a 96 molecule system that was modelled with a Lennard-Jones 6-12 truncated potential. The original insertion method, which does not take into account the orientation of the molecule when it is inserted, gives results for the chemical potential which deviate from that obtained using the thermodynamic pressure integration. At high density or temperature the deviation is significant. We have modified the Widom insertion technique to account for this short range orientation and find good agreement between this technique and the thermodynamic integration method for the chemical potential. We also calculated the free energy difference between our model ring molecules and ring molecules made up of hard spheres.     

Physical�Cchemical determinants of coil conformations in globular proteins

We present a method with the potential to generate a library of coil segments from first principles. Proteins are built from α‐helices and/or β‐strands interconnected by these coil segments. Here, we investigate the conformational determinants of short coil segments, with particular emphasis on chain turns. Toward this goal, we extracted a comprehensive set of two‐, three‐, and four‐residue turns from X‐ray–elucidated proteins and classified them by conformation. A remarkably small number of unique conformers account for most of this experimentally determined set, whereas remaining members span a large number of rare conformers, many occurring only once in the entire protein database. Factors determining conformation were identified via Metropolis Monte Carlo simulations devised to test the effectiveness of various energy terms. Simulated structures were validated by comparison to experimental counterparts. After filtering rare conformers, we found that 98% of the remaining experimentally determined turn population could be reproduced by applying a hydrogen bond energy term to an exhaustively generated ensemble of clash‐free conformers in which no backbone polar group lacks a hydrogen‐bond partner. Further, at least 90% of longer coil segments, ranging from 5‐ to 20 residues, were found to be structural composites of these shorter primitives. These results are pertinent to protein structure prediction, where approaches can be divided into either empirical or ab initio methods. Empirical methods use database‐derived information; ab initio methods rely on physical–chemical principles exclusively. Replacing the database‐derived coil library with one generated from first principles would transform any empirically based method into its corresponding ab initio homologue.     

Potentials for Molecular Dynamics Simulation of Silicate Glasses

A new potential model has been developed for the simulation of amorphous silica based on the ab initio potential model of Pyper. This model promises to be of value in the simulation of silica at high pressures.     

Bias in Markov models of disease

We examine bias in Markov models of diseases, including both chronic and infectious diseases. We consider two common types of Markov disease models: ones where disease progression changes by severity of disease, and ones where progression of disease changes in time or by age. We find sufficient conditions for bias to exist in models with aggregated transition probabilities when compared to models with state/time dependent transition probabilities. We also find that when aggregating data to compute transition probabilities, bias increases with the degree of data aggregation. We illustrate by examining bias in Markov models of Hepatitis C, Alzheimer’s disease, and lung cancer using medical data and find that the bias is significant depending on the method used to aggregate the data. A key implication is that by not incorporating state/time dependent transition probabilities, studies that use Markov models of diseases may be significantly overestimating or underestimating disease progression. This could potentially result in incorrect recommendations from cost-effectiveness studies and incorrect disease burden forecasts.     

Study of Liquid Dimethyl Sulfoxide by Computer Simulation

The paper reports Monte Carlo and molecular dynamicsresults for pure liquid dimethyl sulfoxide (DMSO) at298 K and 1 atm. The classical 6–12 Lennard–Jones plusCoulomb pairwise potential was used to calculateintermolecular interaction energy. Potentialparameters for the liquid were optimized in this work.Some thermodynamic and dynamical properties obtained,such as heat of vaporization, density and diffusioncoefficient, are in good agreement with theexperimental values. The present model is comparedwith other models for DMSO reported previously. It isshown to be an improvement over earlier potentials.The structure factors and the radial distributionfunctions (rdf), are compared with experimentalresults for the liquid. The analysis shows that thestructure of DMSO is not completely understood yet anddeserves deeper investigation. The geometry of thedimer that corresponds to the rdf plots obtained, isreported. The results suggest that the dipole momentof this dimer plays an important role in the structureof the liquid.     

A New Version of the Insertion Particle Method for Determining the Chemical Potential by Monte Carlo Simulation

A new version of the test particle method for determining the chemical potential by Monte Carlo simulations is proposed. The method, applicable to any fluid at any density, combines the Widom's test particle insertion method with the ideas of the scaled particle theory, gradual insertion method and multistage sampling. Its applicability is exemplified by evaluating the chemical potential of the hard sphere fluid at a very high density in semi-grand-canonical and grand-canonical ensembles. A theory estimating the efficiency (i.e. statistical errors) of the method is proposed and the results are compared with the Widom's and gradual insertion methods, and the analytic results.     

Intermolecular Interaction of Fluoro Propanes

Abstract

Intermolecular interaction is investigated for an isomeric pair of fluoro propane, CH₃CF₂CF₃ (HFC–245cb, CB) and CH₂FCF₂CHF₂ (HFC–245ca, CA). CB has a larger dipole moment than CA. This may suggest that CB has a larger intermolecular attractive interaction than CA; the reverse is, however, found from the experimental data: normal boiling point, critical temperature, and heat of vaporization. Systematic ab initio calculations have been done for both CB dimer and CA dimer, and confirmed that the former has a smaller attractive interaction than the latter.

On the basis of these calculations, analytic functions have been constructed as the pair potential models for the two isomers. Each of these models has 11 Lennard-Jones and Coulomb interaction sites in the molecule. The present models can explain why CB dimer has a smaller attractive interaction than CA dimer, and they will easily be extended to a series of fluoro propanes, and make it possible to perform the systematic molecular simulation studies.