Quantum Chemical Microsolvation by Automated Water Placement

We developed a quantitative approach to quantum chemical microsolvation. Key in our methodology is the automatic placement of individual solvent molecules based on the free energy solvation thermodynamics derived from molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). This protocol enabled us to rigorously define the number, position, and orientation of individual solvent molecules and to determine their interaction with the solute based on physical quantities. The generated solute–solvent clusters served as an input for subsequent quantum chemical investigations. We showcased the applicability, scope, and limitations of this computational approach for a number of small molecules, including urea, 2-aminobenzothiazole, (+)-syn-benzotriborneol, benzoic acid, and helicene. Our results show excellent agreement with the available ab initio molecular dynamics data and experimental results.     

Probing the transition from hydrophilic to hydrophobic solvation with atomic scale resolution

Picosecond and femtosecond X-ray absorption spectroscopy is used to probe the changes of the solvent shell structure upon electron abstraction of aqueous iodide using an ultrashort laser pulse. The transient L(1,3) edge EXAFS at 50 ps time delay points to the formation of an expanded water cavity around the iodine atom, in good agreement with classical and quantum mechanical/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. These also show that while the hydrogen atoms pointed toward iodide, they predominantly point toward the bulk solvent in the case of iodine, suggesting a hydrophobic behavior. This is further confirmed by quantum chemical (QC) calculations of I(-)/I(0)(H(2)O)(n=1-4) clusters. The L(1) edge sub-picosecond spectra point to the existence of a transient species that is not present at 50 ps. The QC calculations and the QM/MM MD simulations identify this transient species as an I(0)(OH(2)) complex inside the cavity. The simulations show that upon electron abstraction most of the water molecules move away from iodine, while one comes closer to form the complex that lives for 3-4 ps. This time is governed by the reorganization of the main solvation shell, basically the time it takes for the water molecules to reform an H-bond network. Only then is the interaction with the solvation shell strong enough to pull the water molecule of the complex toward the bulk solvent. Overall, much of the behavior at early times is determined by the reorientational dynamics of water molecules and the formation of a complete network of hydrogen bonded molecules in the first solvation shell.     

Molecular probes of solvation phenomena

The properties of the molecules present in any chemical or biological system are dependent on interactions with the environment, and a quantitative understanding of solvation phenomena remains a major challenge. Molecular recognition probes provide a new approach to quantitatively measure the properties of solvents. Traditionally, solvent polarity scales have been based on spectroscopic probes that provide insight into the nature of solvent-solute interactions. This review compares the solvent polarity parameters obtained from the wavelengths of UV/Visible absorption maxima with solute H-bond parameters obtained from the free energies of solution equilibria. The similarity of the solvent and solute H-bond scales leads to a general H-bond scale that uses the same parameters to describe both solvent and solute. The general H-bond scale provides a framework for understanding the relationship between local intermolecular interactions and the properties of the bulk medium. Intermolecular interactions are sensitive to solvation equilibria, so molecular recognition probes provide fundamentally different information from spectroscopic probes that are sensitive to the populations of different solvation states of the solute. Studies of mixed solvents demonstrate the potential of molecular recognition probes for providing new insights into solvation phenomena.     

A quantum mechanical strategy to investigate the structure of liquids: the cases of acetonitrile, formamide, and their mixture

A computational strategy based on quantum mechanical (QM) calculations and continuum solvation models is used to investigate the structure of liquids (either neat liquids or mixtures). The strategy is based on the comparison of calculated and experimental spectroscopic properties (IR-Raman vibrational frequencies and Raman intensities). In particular, neat formamide, neat acetonitrile, and their equimolar mixture are studied comparing isolated and solvated clusters of different nature and size. In all cases, the study seems to indicate that liquids, even when strongly associated, can be effectively modeled in terms of a shell-like system in which clusters of strongly interacting molecules (the microenvironments) are solvated by a polarizable macroenvironment represented by the rest of the molecules. Only taking into proper account both these effects can a correct picture of the liquid structure be achieved.     

Self-consistent combination of the three-dimensional RISM theory of molecular solvation with analytical gradients and the Amsterdam density functional package

The three-dimensional reference interaction site model with the closure relation by Kovalenko and Hirata (3D-RISM-KH) in combination with the density functional theory (DFT) method has been implemented in the Amsterdam density functional (ADF) software package. The analytical first derivatives of the free energy with respect to displacements of the solute nuclear coordinates have also been developed. This enables study of chemical reactions, including reaction coordinates and transition state search, with the molecular solvation described from the first principles. The method yields all of the features available by using other solvation approaches, for instance infrared spectra of solvated molecules. To evaluate the accuracy of the present method, test calculations have been carried out for a number of small molecules, including four glycine conformers, a set of small organic compounds, and carbon nanotubes of various lengths in aqueous solution. Our predictions for the solvation free energy agree well with other approaches as well as experiment. This new development makes it possible to calculate at modest computational cost the electronic properties and molecular solvation structure of a solute molecule in a given molecular liquid or mixture from the first principles.     

Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

Solvation dynamics in acetonitrile: a study incorporating solute electronic response and nuclear relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.     

Stability of ion triplets in ionic liquid/lithium salt solutions: Insights from implicit and explicit solvent models and molecular dynamics simulations

Binding energies of ion triplets formed in ionic liquids by Li⁺ with two anions have been studied using quantum‐chemical calculations with implicit and explicit solvent supplemented by molecular dynamics (MD) simulations. Explicit solvent approach confirms variation of solute‐ionic liquid interactions at distances up to 2 nm, resulting from structure of solvation shells induced by electric field of the solute. Binding energies computed in explicit solvent and from the polarizable continuum model approach differ largely, even in sign, but relative values generally agree between these two models. Stabilities of ion triplets obtained in quantum‐chemical calculations for some systems disagree with MD results; the discrepancy is attributed to the difference between static optimized geometries used in quantum chemical modeling and dynamic structures of triplets in MD simulations. © 2015 Wiley Periodicals, Inc.     

Solvation of tetraalkylammonium chlorides in acetonitrile-water mixtures: mass spectrometry and molecular dynamics simulations.

The solvation of tetramethylammonium chloride (Me4NCl) and tetra-n-butylammonium chloride (Bu4NCl) in water-acetonitrile mixtures was investigated by mass spectrometry of clusters isolated from the solution. As far as the positive ions are concerned, clusters composed of alkylammonium ions and acetonitrile molecules only were observed, even for mixtures with high water content. In contrast, for the negative ions, clusters composed of chloride with both water and/or acetonitrile molecules were observed. For the smaller system (Me4NCl) we performed quantum chemical calculations and molecular dynamics simulations. It was found that even though water is present in the solvation shell of Me4N+, only acetonitrile has a strong electrostatic interaction with the cation. Water molecules around Me4N+ form hydrogen bonds with other water molecules, and they interact with Me4N+ mainly via dispersive interactions. These results indicate that Me4N+ behaves like a hydrophobic solute. On the other hand, the interaction of Cl- with water and acetonitrile is of comparable strength and, in both cases, the electrostatic interaction dominates. Herein we demonstrate experimentally and theoretically that positive and negative ions give rise to characteristic solvation structures in mixed solvents: even a relatively small organic cation, such as Me4N+, exhibits a hydrophobic-like solvation shell.     

A noniterative mean-field QM/MM-type approach with a linear response approximation toward an efficient free-energy evaluation

Mean-field treatment of solvent provides an efficient technique to investigate chemical processes in solution in quantum mechanics/molecular mechanics (QM/MM) framework. In the algorithm, an iterative calculation is required to obtain the self-consistency between QM and MM regions, which is a time-consuming step. In the present study, we have proposed a noniterative approach by introducing a linear response approximation (LRA) into the solvation term in the one-electron part of Fock matrix in a hybrid approach between molecular-orbital calculations and a three-dimensional (3D) integral equation theory for molecular liquids (multicenter molecular Ornstein–Zernike self-consistent field [MC-MOZ-SCF]; Kido et al., J. Chem. Phys. 2015, 143, 014103). To save the computational time, we have also developed a fast method to generate electrostatic potential map near solute and the solvation term in Fock matrix, using Fourier transformation (FT) and real spherical harmonics expansion (RSHE). To numerically validate the LRA and FT-RSHE method, we applied the present approach to water, carbonic acid, and their ionic species in aqueous solution. Molecular properties of the solutes were evaluated by the present approach with four different types of initial wave functions and compared with those by the original (MC-MOZ-SCF). We found that an initial wave function considering solvation effects is needed to appropriately reproduce the properties by MC-MOZ-SCF. Furthermore, a benchmark test for 32 solute molecules was performed to evaluate the accuracy of the present approach for solvation free energy (SFE) and measure the speedup ratio for MC-MOZ-SCF. The error of SFE for MC-MOZ-SCF does not correlate with the SFE but increases in proportion to the electronic reorganization energy. Similar to water and carbonic acid, an initial wave function with solvation effects is also important to make the error small. From the averaged speed up ratio, the present approach is 13.5 times faster than MC-MOZ-SCF. © 2019 Wiley Periodicals, Inc.     

p-π共轭分子的氨基面外弯曲振动模的拉曼光谱

拉曼光谱是一种用途广泛的无损分子检测技术,其能够提供化学物质的分子结构指纹信息.一种面外弯曲振动模被称作wagging振动,它的信号尤为特殊,其频率和强度都非常依赖于检测环境.以乙烯胺和苯胺为例,采用密度泛函理论计算研究了p-π共轭分子分别与水簇和银簇作用的平衡结构、成键作用和拉曼光谱.结果表明,弱相互作用,如分子与金属表面的弱吸附作用以及分子与水之间的氢键作用,均使氨基面外弯曲振动模(ωNH2)的拉曼信号发生显著的变化.考虑溶剂化效应后,氢键作用减弱,计算拉曼光谱趋于一致.通过进一步对电子结构的分析,解释了面外弯曲振动信号显著增强的原因,揭示了面外弯曲振动模与分子p-π共轭作用之间的关系.     

Solute-solvent friction kernels and solution properties of methyl oxazoline-phenyl oxazoline (MeOx-PhOx) copolymers in binary ethanol-water mixtures

Solvent mixtures often alter the solubility of polymeric substances. Statistical copolymers made from 2-methyl-2-oxazoline (MeOx) and 2-phenyl-2-oxazoline (PhOx) are known for their varying solubilities in pure ethanol, pure water and in binary mixtures of ethanol-water. Constrained Molecular Dynamics (MD) simulations have been carried out with an aim to explain the varying solubilities of the statistical MeOx-PhOx copolymers. The solute-solvent dynamic friction kernels calculated through constrained MD simulations corroborate the solubility pattern in these copolymers. The solvation characteristics have been analyzed in terms of the solute-solvent radial distribution functions (RDFs). The ethanol-soluble MeOx-PhOx copolymers exhibit characteristic solute-composition dependence in the dynamic solute-solvent friction kernels, indicating the strength of the solute-solvent correlations. The aggressive solvation by the ethanol molecules in the binary solvent mixtures has been brought out by the O(solute)-H(ethanol) RDFs which exhibit a characteristic dependence on the ethanol content in the solvent composition. The corresponding O(solute)-H(water) RDFs are devoid of any such composition dependence. For all the MeOx-PhOx copolymers, the O-site solvation is strongly dominated by the water molecules and the N-sites are solvated equally by both ethanol and water molecules.     

Effects of solute structure on local solvation and solvent interactions: results from UV/vis spectroscopy and molecular dynamics simulations

Solvation of heterocyclic amines in CO(2)-expanded methanol (MeOH) has been explored with UV/vis spectroscopy and molecular dynamics (MD) simulations. A synergistic study of experiments and simulations allows exploration of solute and solvent effects on solvation and the molecular interactions that affect absorption. MeOH-nitrogen hydrogen bonds hinder the n-pi* transition; however, CO(2) addition causes a blue shift relative to MeOH because of Lewis acid/base interactions with nitrogen. Effects of solute structure are considered, and very different absorption spectra are obtained as nitrogen positions change. MD simulations provide detailed solvent clustering behavior around the solute molecules and show that the local solvent environment and ultimately the spectra are sensitive to the solute structure. This work demonstrates the importance of atomic-level information in determining the structure-property relationships between solute structure, local salvation, and solvatochromism.     

Infrared photodissociation spectroscopy of Na(NH3)n clusters: probing the solvent coordination

The first mass-selective vibrational spectra have been recorded for Na(NH3)n clusters. Infrared spectra have been obtained for n = 3-8 in the N-H stretching region. The spectroscopic work has been supported by ab initio calculations carried out at both the DFT(B3LYP) and MP2 levels, using a 6-311++G(d,p) basis set. The calculations reveal that the lowest energy isomer for n or= 7 is indicative of molecules entering a second solvation shell, i.e., the inner solvation shell around the sodium atom can accommodate a maximum of six NH3 molecules.     

The structure of the hydrated electron. Part 2. A mixed quantum/classical molecular dynamics embedded cluster density functional theory: single-excitation configuration interaction study

Adiabatic mixed quantum/classical (MQC) molecular dynamics (MD) simulations were used to generate snapshots of the hydrated electron in liquid water at 300 K. Water cluster anions that include two complete solvation shells centered on the hydrated electron were extracted from the MQC MD simulations and embedded in a roughly 18 Ax18 Ax18 A matrix of fractional point charges designed to represent the rest of the solvent. Density functional theory (DFT) with the Becke-Lee-Yang-Parr functional and single-excitation configuration interaction (CIS) methods were then applied to these embedded clusters. The salient feature of these hybrid DFT(CIS)/MQC MD calculations is significant transfer (approximately 18%) of the excess electron's charge density into the 2p orbitals of oxygen atoms in OH groups forming the solvation cavity. We used the results of these calculations to examine the structure of the singly occupied and the lower unoccupied molecular orbitals, the density of states, the absorption spectra in the visible and ultraviolet, the hyperfine coupling (hfcc) tensors, and the infrared (IR) and Raman spectra of these embedded water cluster anions. The calculated hfcc tensors were used to compute electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectra for the hydrated electron that compared favorably to the experimental spectra of trapped electrons in alkaline ice. The calculated vibrational spectra of the hydrated electron are consistent with the red-shifted bending and stretching frequencies observed in resonance Raman experiments. In addition to reproducing the visible/near IR absorption spectrum, the hybrid DFT model also accounts for the hydrated electron's 190-nm absorption band in the ultraviolet. Thus, our study suggests that to explain several important experimentally observed properties of the hydrated electron, many-electron effects must be accounted for: one-electron models that do not allow for mixing of the excess electron density with the frontier orbitals of the first-shell solvent molecules cannot explain the observed magnetic, vibrational, and electronic properties of this species. Despite the need for multielectron effects to explain these important properties, the ensemble-averaged radial wavefunctions and energetics of the highest occupied and three lowest unoccupied orbitals of the hydrated electrons in our hybrid model are close to the s- and p-like states obtained in one-electron models. Thus, one-electron models can provide a remarkably good approximation to the multielectron picture of the hydrated electron for many applications; indeed, the two approaches appear to be complementary.     

Improving the performance of the coupled reference interaction site model-hyper-netted chain (RISM-HNC)/simulation method for free energy of solvation

The coupled reference interaction site model-hyper-netted chain (RISM-HNC)/ simulation methodology determines solvation free energies as a function of the set of all radial distribution functions of solvent atoms about atomic solute sites. These functions are determined from molecular dynamics (MD) or Monte Carlo (MC) simulations rather than from solving the RISM and HNC equations iteratively. Previous applications of the method showed that it can predict relative free energies of solvation for small solutes accurately. However, the errors scale with the system size. In this study, we propose the use of the hard-sphere free energy as the reference and a linear response approximation to improve the performance, i.e., accuracy and robustness, of the method, particularly removing the size dependency of the error. The details of the new formalism are presented. To validate the proposed formalism, solvation free energies of N-methylacetamide and methylamine are computed using the new RISM-HNC-based expressions in addition to a linear response expression, which are compared to previous thermodynamic integration and thermodynamic perturbation results performed with the same force field. Additionally, free energies of solvation for cyclohexane, pyridine, benzene and derivatives, and other small organic molecules are calculated and compared to experimental values.     

Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.     

Complexes of a Zn‐metalloenzyme binding site with hydroxamate‐containing ligands. A case for detailed benchmarkings of polarizable molecular mechanics/dynamics potentials when the experimental binding structure is unknown

Zn‐metalloproteins are a major class of targets for drug design. They constitute a demanding testing ground for polarizable molecular mechanics/dynamics aimed at extending the realm of quantum chemistry (QC) to very long‐duration molecular dynamics (MD). The reliability of such procedures needs to be demonstrated upon comparing the relative stabilities of competing candidate complexes of inhibitors with the recognition site stabilized in the course of MD. This could be necessary when no information is available regarding the experimental structure of the inhibitor–protein complex. Thus, this study bears on the phosphomannose isomerase (PMI) enzyme, considered as a potential therapeutic target for the treatment of several bacterial and parasitic diseases. We consider its complexes with 5‐phospho‐d ‐arabinonohydroxamate and three analog ligands differing by the number and location of their hydroxyl groups. We evaluate the energy accuracy expectable from a polarizable molecular mechanics procedure, SIBFA. This is done by comparisons with ab initio quantum‐chemistry (QC) calculations in the following cases: (a) the complexes of the four ligands in three distinct structures extracted from the entire PMI‐ligand energy‐minimized structures, and totaling up to 264 atoms; (b) the solvation energies of several energy‐minimized complexes of each ligand with a shell of 64 water molecules; (c) the conformational energy differences of each ligand in different conformations characterized in the course of energy‐minimizations; and (d) the continuum solvation energies of the ligands in different conformations. The agreements with the QC results appear convincing. On these bases, we discuss the prospects of applying the procedure to ligand‐macromolecule recognition problems. © 2016 Wiley Periodicals, Inc.