Prediction of octanol-water partition coefficients for the SAMPL6-[Formula: see text] molecules using molecular dynamics simulations with OPLS-AA, AMBER and CHARMM force fields.

All-atom molecular dynamics simulations with stratified alchemical free energy calculations were used to predict the octanol-water partition coefficient [Formula: see text] of eleven small molecules as part of the SAMPL6-[Formula: see text] blind prediction challenge using four different force field parametrizations: standard OPLS-AA with transferable charges, OPLS-AA with non-transferable CM1A charges, AMBER/GAFF, and CHARMM/CGenFF. Octanol parameters for OPLS-AA, GAFF and CHARMM were validated by comparing the density as a function of temperature, the chemical potential, and the hydration free energy to experimental values. The partition coefficients were calculated from the solvation free energy for the compounds in water and pure ("dry") octanol or "wet" octanol with 27 mol% water dissolved. Absolute solvation free energies were computed by thermodynamic integration (TI) and the multistate Bennett acceptance ratio with uncorrelated samples from data generated by an established protocol using 5-ns windowed alchemical free energy perturbation (FEP) calculations with the Gromacs molecular dynamics package. Equilibration of sets of FEP simulations was quantified by a new measure of convergence based on the analysis of forward and time-reversed trajectories. The accuracy of the [Formula: see text] predictions was assessed by descriptive statistical measures such as the root mean square error (RMSE) of the data set compared to the experimental values. Discarding the first 1 ns of each 5-ns window as an equilibration phase had a large effect on the GAFF data, where it improved the RMSE by up to 0.8 log units, while the effect for other data sets was smaller or marginally worsened the agreement. Overall, CGenFF gave the best prediction with RMSE 1.2 log units, although for only eight molecules because the current CGenFF workflow for Gromacs does not generate files for certain halogen-containing compounds. Over all eleven compounds, GAFF gave an RMSE of 1.5. The effect of using a mixed water/octanol solvent slightly decreased the accuracy for CGenFF and GAFF and slightly increased it for OPLS-AA. The GAFF and OPLS-AA results displayed a systematic error where molecules were too hydrophobic whereas CGenFF appeared to be more balanced, at least on this small data set.

Prediction of cyclohexane-water distribution coefficients for the SAMPL5 data set using molecular dynamics simulations with the OPLS-AA force field.

All-atom molecular dynamics simulations were used to predict water-cyclohexane distribution coefficients [Formula: see text] of a range of small molecules as part of the SAMPL5 blind prediction challenge. Molecules were parameterized with the transferable all-atom OPLS-AA force field, which required the derivation of new parameters for sulfamides and heterocycles and validation of cyclohexane parameters as a solvent. The distribution coefficient was calculated from the solvation free energies of the compound in water and cyclohexane. Absolute solvation free energies were computed by an established protocol using windowed alchemical free energy perturbation with thermodynamic integration. This protocol resulted in an overall root mean square error in [Formula: see text] of almost 4 log units and an overall signed error of -3 compared to experimental data. There was no substantial overall difference in accuracy between simulating in NVT and NPT ensembles. The signed error suggests a systematic error but the experimental [Formula: see text] data on their own are insufficient to uncover the source of this error. Preliminary work suggests that the major source of error lies in the hydration free energy calculations.

13C longitudinal relaxation time measurements and DFT-GIAO NMR computations for two ammonium ions of a tetraazamacrocyclic scorpiand system.

ABSTRACT: Spin-lattice relaxation times, T1s, for 13C nuclei in two cations H n1n+ (n = 1, 5) of N-(2-aminoethyl)-cyclam (1, scorpiand) were determined by means of 13C{1H} NMR experiments in aqueous solution at pH 11.5 and 0.2. The theoretical study [modeling with OPLS-AA, B3LYP/6-31G(d) geometry optimizations, dispersion-corrected energies (DFT-D3), and DFT-GIAO predictions of the NMR chemical shifts (including an IEF-PCM simulation of hydration)] was also done for several conformers of the tautomer iso-H414+ not investigated before. The binding directions in protonated polyamino receptors necessary for efficient complexation of the nitrate anion(s) were briefly outlined, as well. All these results were discussed in terms of 'abnormal' 13C chemical shift changes found previously for the side-chain carbons of amine 1 in strongly acidic solution (HNO3). In conclusion, an earlier proposal of its association with NO3- at pH <1 was rejected. Instead, the participation of small amounts of a micro-species iso-H414+Dhydr under such conditions can be proposed. GRAPHICAL ABSTRACT: A small contribution of iso-H414+Dhydr (see figure) to an ionic mixture of pentamine 1 was proposed to explain the 'abnormal' 13C NMR shifts observed for atoms C11 and C12 in its side-chain arm, at pH <1.

Molecular dynamics simulation of salt diffusion in constituting phosphazene-based polymer electrolyte.

A growing demand to visualize polymer models in liquid poses a computational challenge in molecular dynamics (MD) simulation, as this requires emerging models under suitable force fields (FFs) to capture the underlying molecular behaviour accurately. In our present study, we have employed TIP3P potential on water and all atomistic optimized potentials for liquid simulations FFs to study the liquid electrolyte behavior of phosphazene-based polymer by considering its potential use in lithium-ion polymer batteries. We have explored the polymer's local structure, chain packing, wettability, and hydrophobic tendencies against the silicon surface using a combination of a pseudocontinuum model in MD simulation, and surface-sensitive sum frequency generation (SFG) vibrational spectroscopy. The finding yields invaluable insights into the molecular architecture of phosphazene. This approach identifies the importance of hydrophobic interactions with air and hydrophilic units with water molecules in understanding the behavior and properties of phosphazene-based polymers at interfaces, contributing to its advancements in materials science. The MD study uniquely captures traces of the polymer-ion linkage, which is observed to become more pronounced with the increase in polymer weight fraction. The theoretical observation of this linkage's influence on lithium-ion diffusion motion offers valuable insights into the fundamental physics governing the behavior of atoms and molecules within phosphazene-based polymer electrolytes in aqueous environments. Further these predictions are corroborated in the molecular-level depiction at the air-aqueous interface, as evidenced from the OH-oscillator strength variation measured by the SFG spectroscopy.The fundamental findings from this study open new avenues for utilizing MD simulation as a versatile methodology to gain profound insights into intermolecular interactions of polymer. It could be useful in the application of biomedical and energy-related research, such as polymer lithium-ion batteries, fuel cells, and organic solar cells.

Cis and trans unsaturated phosphatidylcholine bilayers: A molecular dynamics simulation study.

Trans unsaturated lipids are uncommon in nature. In the human diet, they occur as natural products of ruminal bacteria or from industrial food processing like hydrogenation of vegetable oils. Consumption of trans unsaturated lipids has been shown to have a negative influence on human health; in particular, the risk of cardiovascular disease is higher when the amount of trans unsaturated lipids in the diet is elevated. In this study, we first performed quantum mechanical calculations to specifically and accurately parameterize cis and trans mono-unsaturated lipids and subsequently validated the newly derived parameter set. Then, we carried out molecular dynamics (MD) simulations of lipid bilayers composed of cis or trans unsaturated lipids with and without cholesterol. Our results show that trans mono-unsaturated chains are more flexible than cis mono-unsaturated chains due to lower barriers for rotation around the single bonds next to the trans double bond than those next to the cis double bond. In effect, interactions between cholesterol and trans unsaturated chains are stronger than cis unsaturated chains, which results in a higher ordering effect of cholesterol in trans unsaturated bilayers.