Upgrading of the general AMBER force field 2 for fluorinated alcohol biosolvents: A validation for water solutions and melittin solvation.

The standard GAFF2 force field parameterization has been refined for the fluorinated alcohols 2,2,2-trifluoroethanol (TFE), 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP), and 1,1,1,3,3,3-hexafluoropropan-2-one (HFA), which are commonly used to study proteins and peptides in biomimetic media. The structural and dynamic properties of both proteins and peptides are significantly influenced by the biomimetic environment created by the presence of these cosolvents in aqueous solutions. Quantum mechanical calculations on stable conformers were used to parameterize the atomic charges. Different systems, such as pure liquids, aqueous solutions, and systems formed by melittin protein and cosolvent/water solutions, have been used to validate the new models. The calculated macroscopic and structural properties are in agreement with experimental findings, supporting the validity of the newly proposed models.

Concerted versus stepwise proton transfer reactions in the [2, 2'-bipyridyl]-3-3'-diol molecule: A static and dynamic ab-initio investigation.

The double proton transfer (PT) reaction has been investigated in the [2,2'-bipyridyl]-3-3'-diol, a complex molecule where the proton movements is coupled to significant rearrangement of the electronic structure. Moreover, the reaction could be concerted, that is the two protons are exchanged simultaneously, or stepwise, where the two protons are transferred sequentially. To this end, a static exploration of the potential energy surface (PES) was carried together with the analysis of the free-energy surface (FES), both surfaces being evaluated at density functional theory level and different exchange-correlation functionals. While the concerted mechanism has been clearly discharged, the characteristics of the stepwise PT significantly depends on the chosen functionals, some suggesting a clear stepwise mechanism characterized by a stable reaction intermediates and two transitions states, whereas other approaches propend for a asynchronous PT, with a single TS. These features appear on both PES and FES, albeit some differences appears due to their different nature.

Absorption Spectra of Flexible Fluorescent Probes by a Combined Computational Approach: Molecular Dynamics Simulations and Time-Dependent Density Functional Theory.

A detailed understanding and interpretation of absorption spectra of molecular systems, especially in condensed phases, requires computational models that allow their structural and electronic features to be connected to the observed macroscopic spectra. This work is focused on modeling the electronic absorption spectrum of a fluorescent probe, namely, the 9-(4-((bis(2-((2-(ethylthio)ethyl)thio)ethyl)amino)methyl)phenyl)-6-(pyrrolidin-1-yl)-3H-xanthen-3-one molecule, depicted by a combined classical-quantum chemical approach. Particularly, first classical molecular dynamics (MD) has been used to explore the configurational space, and next, the absorption spectrum has been reconstructed by averaging the results of time-dependent density functional theory (TD-DFT) calculations performed on equispaced molecular conformations extracted from MD to properly sample the configurational space explored at finite temperature. To verify the effect of molecular conformation on the spectral profile, the generated electronic absorption spectra were compared with those obtained considering a single structure corresponding to the optimized one, an approach also referred to as static. This comparison allows one to highlight a sizable though small shift between the maxima of the corresponding reconstructed absorption spectra, highlighting the importance of conformational sampling in the case of this rather flexible molecule. Four different exchange and correlation functionals (PBE, BLYP, PBE0, B3LYP) were considered to compute vertical transition via TD-DFT calculations. From the results obtained in gas and in condensed, here solution, phases, it appears that the magnitude of the shift is actually more affected by the phase in which the system is found than by the functional used. This fact underlines the central importance of conformational mobility, that is flexibility, of this molecule. From a more quantitative point of view, a comparison with available experimental data shows that hybrid functionals, such as PBE0 and B3LYP, enable one to faithfully reproduce the observed absorption maxima.

Exploring the effect of Mg2+ substitution on amorphous calcium phosphate nanoparticles.

HYPOTHESIS: The study of Amorphous Calcium Phosphate (ACP) has become a hot topic due to its relevance in living organisms and as a material for biomedical applications. The preparation and characterization of Mg-substituted ACP nanoparticles (AMCP) with tunable Ca/Mg ratio is reported in the present study to address the effect of Mg2+ on their structure and stability. EXPERIMENTS: AMCPs particles were synthesized by precipitation of the precursors from aqueous solutions. The particles were analyzed in terms of morphology, crystallinity, and thermal stability, to get a complete overview of their physico-chemical characteristics. Computational methods were also employed to simulate the structure of ACP clusters at different levels of Mg2+ substitution. FINDINGS: Our results demonstrate that AMCP particles with tunable composition and crystallinity can be obtained. The analysis of the heat-induced crystallization of AMCP shows that particles' stability depends on the degree of Mg2+ substitution in the cluster, as confirmed by computational analyses. The presented results shed light on the effect of Mg2+ on ACP features at different structural levels and may be useful guidelines for the preparation and design of AMCP particles with a specific Ca/Mg ratio.

Toward an Understanding of the Pressure Effect on the Intramolecular Vibrational Frequencies of Sulfur Hexafluoride.

The structural and vibrational properties of the molecular units of sulfur hexafluoride crystal as a function of pressure have been studied by the Extreme Pressure Polarizable Continuum Model (XP-PCM) method. Within the XP-PCM model, single molecule calculations allow a consistent interpretation of the experimental measurements when considering the effect of pressure on both the molecular structure and the vibrational normal modes. This peculiar aspect of XP-PCM provides a detailed description of the electronic origin of normal modes variations with pressure, via the curvature of the potential energy surface and via the anharmonicity of the normal modes. When applied to the vibrational properties of the sulfur hexafluoride crystal, the XP-PCM method reveals a hitherto unknown interpretation of the effects of the pressure on the vibrational normal modes of the molecular units of this crystal.

"Cyclopropylidene Effect" in the 1,3-Dipolar Cycloaddition of Nitrones to Alkylidene Cyclopropanes: A Computational Rationalization.

The regioselectivity in the 1,3-dipolar cycloaddition (1,3-DC) between five-membered cyclic nitrone and methylenecyclopropane (MCP) has been studied through density functional theory (DFT) calculations. The computational study of 1,3-DC with different 1-alkyl- (or 1,1-dialkyl)-substituted alkenes and the comparison with MCP have evidenced that the electrostatic interaction has a central role in the regioselectivity of the reactions. It has been observed that the electronic effect of the substituent (donor or attractor groups) determines the polarization of the alkene double bond and the reaction mechanism, consequently determining the interaction with nitrones and favoring an orientation between this moiety and the dipolarophile.

Evidence of a Low-High Density Turning Point in Liquid Water at Ordinary Temperature under Pressure: A Molecular Dynamics Study.

Water has a fundamental role in important processes spanning a wide range of pressure and temperature conditions. Knowledge of structural, dynamic and thermodynamic properties of water at nonstandard conditions is a primary concern since interest in astronomical, geological, and technological processes is continuously growing. Molecular dynamics simulations allow us to study thermodynamic conditions that require sophisticated techniques and instruments, while at the same time offering the interpretation of properties at the atomic level. It is established that the behavior of water is strongly affected by the temperature and pressure conditions, determining the existence of low and high density regimes. For the first time, a thermodynamic property, isothermal compressibility, has been adopted to detect the low-high density turning point at ambient temperature in liquid water due to pressure. Molecular dynamics simulations have been performed with five three-site models, allowing us to characterize the complexity of water nature at these conditions at the atomic level.

Computational Investigation of the Selective Cleavage of Diastereotopic Cyclopropane Bonds in 5-Spirocyclopropane Isoxazolidines Rearrangement.

The complete path of the Brandi-Guarna rearrangement of 5-spirocyclopropane isoxazolidines has been investigated by means of density functional theory calculations to rationalize the competing formation of tetrahydropyridones and enaminones by the determination of the minimum energy reaction paths. Our calculations confirm that the rearrangement is triggered by the homolysis of the isoxazolidine N-O bond followed by cleavage of one of the two C-CH2 cyclopropane bonds as previously proposed by the Fabian group [ Eur. J. Org. Chem. 2001, 2001, 4223]. In addition, the results of this work suggest that in the presence of a stereogenic center at isoxazolidine C-4', the formation of a piperidinone or an enaminone as the final product depends on which of the two diastereotopic C-CH2 bonds of cyclopropane is cleaved in the second step of the process. The result can be of great interest for the understanding of other processes involving the opening of a cyclopropane ring.

Imidazole in Aqueous Solution: Hydrogen Bond Interactions and Structural Reorganization with Concentration.

The structural and dynamic properties of imidazole in aqueous solution have been studied by means of classical and ab initio molecular dynamics simulations. We developed a new force field for the imidazole molecule with improved modeling of the electrostatic interactions, specifically tailored to address the well-known drawbacks of existing force fields based on the atomic fractional charge approach. To this end, we reparametrized the charge distribution on the heterocyclic ring, introducing an extra site accounting for the lone pair on the deprotonated nitrogen. The accuracy of the model in describing the hydrogen bond pattern in the aqueous solvent has been confirmed by comparing the classical results on imidazole-water interactions to accurate Car-Parrinello molecular dynamics simulations. It reproduces satisfactorily the experimental water/octanol partition coefficient of imidazole, as well as the structure of the imidazole molecular crystal. The force field has been finally applied to simulate aqueous solutions at various imidazole concentrations to obtain information on both imidazole-water and imidazole-imidazole interactions, providing a description of the different molecular arrangements in solution.

Correspondence between light-absorption spectrum and nonequilibrium work distribution as a mean to access free energy differences between electronic states.

The problem of recovering the free energy difference between two electronic states has been investigated by Frezzato [Chem. Phys. Lett. 533, 106 (2012)], exploring the equivalence between light-absorption spectra and work distribution, hence opening to the application of a spectroscopic version of the Jarzynski equality (JE) [Phys. Rev. Lett. 78, 2690 (1997)]. Here, assuming the validity of the time-dependent perturbation theory, we demonstrate that such equivalence does not lead to the known form of the JE. This is ascribed to the fact that light-absorption processes cannot be described as stochastic processes. To emphasize such an aspect, we devise a stochastic model for the UV-vis (ultraviolet and visible) absorption, suitable for determining the free energy difference between two generic quantum manifolds in a JE-like fashion. However, the model would require explicit knowledge of the transition dipole moments, which are in general not available. Nonetheless, we derive a spectroscopic version of the JE that allows us to recover the free energy difference between the ground and an excited electronic state when the latter state is the only one observed in the spectrum.

Binding Free Energies of Host-Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Illustrative Calculations and Numerical Validation.

In the companion article (Giovannelli et al., 10.1021/acs.jctc.7b00594), we presented an alchemical approach, based on nonequilibrium molecular dynamics simulations, to compute absolute binding free energies of a generic host-guest system. Two alternative computational routes, called binded-domain and single-point alchemical-path schemes, have been proposed. This study is addressed to furnish numerical validation and illustrative examples of the above-mentioned alchemical schemes. Validation is provided by comparing binding free-energy data relative to two poses of a Zn(II)·anion complex with those recovered from an alternative approach, based on steered molecular dynamics simulations. We illustrate important technical and theoretical aspects for a good practice in applying both alchemical schemes, not only through the calculations on the Zn(II)·anion complex, but also estimating absolute binding free energies of 1:1 complexes of ß-cyclodextrin with aromatic compounds (benzene and naphthalene). Comparison with experimental data and previous molecular dynamics simulation studies further confirms the validity of the present nonequilibrium-alchemical methodology.

Binding Free Energies of Host-Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Theoretical Framework.

The fast-switching decoupling method is a powerful nonequilibrium technique to compute absolute binding free energies of ligand-receptor complexes (Sandberg et al., J. Chem. Theory Comput. 2014, 11, 423-435). Inspired by the theory of noncovalent binding association of Gilson and co-workers (Biophys. J. 1997, 72, 1047-1069), we develop two approaches, termed binded-domain and single-point alchemical-path schemes (BiD-AP and SiP-AP), based on the possibility of performing alchemical trajectories during which the ligand is constrained to fixed positions relative to the receptor. The BiD-AP scheme exploits a recent generalization of nonequilibrium work theorems to estimate the free energy difference between the coupled and uncoupled states of the ligand-receptor complex. With respect to the fast-switching decoupling method without constraints, BiD-AP prevents the ligand from leaving the binding site, but still requires an estimate of the positional binding-site volume, which may not be a simple task. On the other side, the SiP-AP scheme allows avoidance of the calculation of the binding-site volume by introducing an additional equilibrium simulation of ligand and receptor in the bound state. In the companion article (DOI: 10.1021/acs.jctc.7b00595), we show that the extra computational effort required by SiP-AP leads to a significant improvement of accuracy in the free energy estimates.

Insights on the Realgar Crystal Under Pressure from XP-PCM and Periodic Model Calculations.

The spectroscopic properties of As4S4 with pressure have been computed by the quantum mechanical XP-PCM method and by density functional theory periodic calculations. The comparison has allowed the interpretation of the available experimental data. By comparison of the two methods and with experiments, we show that the XP-PCM method is able to reproduce the same behavior of the periodic calculations with much lower computational cost allowing to be adopted as a first choice computational tool for a qualitative interpretation of molecular crystals properties under pressure.

XP-PCM Calculations of High Pressure Structural and Vibrational Properties of P4S3.

The structure and the vibrational properties of the P4S3 crystal at high pressures are discussed by application of the XP-PCM method. The vibrational assignment has been clarified. The structure and the electron distribution changes as a function of pressure are analyzed. The pressure effect on the vibrational frequencies is satisfactorily reproduced and discussed in terms of confinement and structure relaxation contributions.

Elastic Barrier Dynamical Freezing in Free Energy Calculations: A Way To Speed Up Nonequilibrium Molecular Dynamics Simulations by Orders of Magnitude.

An important issue concerning computer simulations addressed to free energy estimates via nonequilibrium work theorems, such as the Jarzynski equality [Phys. Rev. Lett. 1997, 78, 2690], is the computational effort required to achieve results with acceptable accuracy. In this respect, the dynamical freezing approach [Phys. Rev. E 2009, 80, 041124] has been shown to improve the efficiency of this kind of simulations, by blocking the dynamics of particles located outside an established mobility region. In this report, we show that dynamical freezing produces a systematic spurious decrease of the particle density inside the mobility region. As a consequence, the requirements to apply nonequilibrium work theorems are only approximately met. Starting from these considerations, we have developed a simulation scheme, called "elastic barrier dynamical freezing", according to which a stiff potential-energy barrier is enforced at the boundaries of the mobility region, preventing the particles from leaving this region of space during the nonequilibrium trajectories. The method, tested on the calculation of the distance-dependent free energy of a dimer immersed into a Lennard-Jones fluid, provides an accuracy comparable to the conventional steered molecular dynamics, with a computational speedup exceeding a few orders of magnitude.

Simulations in generalized ensembles through noninstantaneous switches.

Generalized-ensemble simulations, such as replica exchange and serial generalized-ensemble methods, are powerful simulation tools to enhance sampling of free energy landscapes in systems with high energy barriers. In these methods, sampling is enhanced through instantaneous transitions of replicas, i.e., copies of the system, between different ensembles characterized by some control parameter associated with thermodynamical variables (e.g., temperature or pressure) or collective mechanical variables (e.g., interatomic distances or torsional angles). An interesting evolution of these methodologies has been proposed by replacing the conventional instantaneous (trial) switches of replicas with noninstantaneous switches, realized by varying the control parameter in a finite time and accepting the final replica configuration with a Metropolis-like criterion based on the Crooks nonequilibrium work (CNW) theorem. Here we revise these techniques focusing on their correlation with the CNW theorem in the framework of Markovian processes. An outcome of this report is the derivation of the acceptance probability for noninstantaneous switches in serial generalized-ensemble simulations, where we show that explicit knowledge of the time dependence of the weight factors entering such simulations is not necessary. A generalized relationship of the CNW theorem is also provided in terms of the underlying equilibrium probability distribution at a fixed control parameter. Illustrative calculations on a toy model are performed with serial generalized-ensemble simulations, especially focusing on the different behavior of instantaneous and noninstantaneous replica transition schemes.

Computing Free Energy Differences of Configurational Basins.

A simulation-based approach is proposed to estimate free energy differences between configurational states A and B, defined in terms of collective coordinates of the molecular system. The computational protocol is organized into three stages that can be carried on simultaneously. Two of them consist of independent simulations aimed at sampling, in turn, A and B states. In order to limit the evolution of the system around A and B, biased sampling simulations such as umbrella sampling can be employed. These simulations allow us to estimate local configuration integrals associated with A and B, which can be viewed as vibrational contributions to the free energy. Free energy evaluation is completed by the linking-path stage, in which the potential of mean force difference is estimated between two arbitrary points of the configurational surface, located the first around A and the second around B. The linking path in the space of the collective coordinates is arbitrary and can be computed with any method, starting from adaptive biasing potential/force approaches to nonequilibrium techniques. As an illustrative example, we present the calculation of free energy differences between conformational states of the alanine dipeptide in the space of backbone dihedral angles. The basic advantage of this method, that we term "path-linked domains" scheme, is to prevent accurate calculation of the whole free energy hypersurface in the space of the collective coordinates, thus limiting the statistical sampling to a minimum. Path-linked domains schemes can be applied to a variety of biochemical processes, such as protein-ligand complexation or folding-unfolding interconversion.

Identification of Di(oxymethylene)glycol in the Raman Spectrum of Formaldehyde Aqueous Solutions by ab Initio Molecular Dynamics Simulations and Quantum Chemistry Calculations.

Di(oxymethylene)glycol forms in formaldehyde aqueous solutions by polymerization of methanediol. The structure and hydrogen bond interactions of di(oxymethylene)glycol with water were characterized by performing Car-Parrinello molecular dynamics simulations. The anharmonic vibrational frequencies of di(oxymethylene)glycol in solution were determined with ab initio calculations considering explicitly the hydrogen-bonded water molecules, while other interactions with solvent were described within a polarizable continuum model approach. The calculations allow for a detailed interpretation of the experimental Raman spectrum of formaldehyde aqueous solutions, leading to the assignment of the band at 920 cm(-1) to the symmetric CO stretching mode of di(oxymethylene)glycol.

Annealed importance sampling with constant cooling rate.

Annealed importance sampling is a simulation method devised by Neal [Stat. Comput. 11, 125 (2001)] to assign weights to configurations generated by simulated annealing trajectories. In particular, the equilibrium average of a generic physical quantity can be computed by a weighted average exploiting weights and estimates of this quantity associated to the final configurations of the annealed trajectories. Here, we review annealed importance sampling from the perspective of nonequilibrium path-ensemble averages [G. E. Crooks, Phys. Rev. E 61, 2361 (2000)]. The equivalence of Neal's and Crooks' treatments highlights the generality of the method, which goes beyond the mere thermal-based protocols. Furthermore, we show that a temperature schedule based on a constant cooling rate outperforms stepwise cooling schedules and that, for a given elapsed computer time, performances of annealed importance sampling are, in general, improved by increasing the number of intermediate temperatures.

Structural and spectroscopic properties of methanediol in aqueous solutions from quantum chemistry calculations and ab initio molecular dynamics simulations.

The structural, electronic, and spectroscopic properties of methanediol in aqueous solutions have been studied by a combined approach based on Car-Parrinello molecular dynamics simulations and ab initio calculations. The hydrogen bond interactions between the solute and water have been characterized, showing the important role of the solvent in the stabilization of the methanediol conformers in solution. First insights on the experimental vibrational spectra have been obtained by the analysis of the simulation results, with particular regard to the most prominent band at 1050 cm(-1) that has been attributed to both the symmetric and antisymmetric CO stretching modes. The assignment has been completed adopting both electric and mechanical anharmonic calculations considering the interactions with the solvent using a polarizable continuum model.