Hydrogen bond lifetimes and energetics for solute/solvent complexes studied with 2D-IR vibrational echo spectroscopy

Weak pi hydrogen-bonded solute/solvent complexes are studied with ultrafast two-dimensional infrared (2D-IR) vibrational echo chemical exchange spectroscopy, temperature-dependent IR absorption spectroscopy, and density functional theory calculations. Eight solute/solvent complexes composed of a number of phenol derivatives and various benzene derivatives are investigated. The complexes are formed between the phenol derivative (solute) in a mixed solvent of the benzene derivative and CCl4. The time dependence of the 2D-IR vibrational echo spectra of the phenol hydroxyl stretch is used to directly determine the dissociation and formation rates of the hydrogen-bonded complexes. The dissociation rates of the weak hydrogen bonds are found to be strongly correlated with their formation enthalpies. The correlation can be described with an equation similar to the Arrhenius equation. The results are discussed in terms of transition state theory.     

Ultrafast two-dimensional infrared vibrational echo chemical exchange experiments and theory

Ultrafast two-dimensional (2D) infrared vibrational echo experiments and theory are used to examine chemical exchange between solute-solvent complexes and the free solute for the solute phenol and three solvent complex partners, p-xylene, benzene, and bromobenzene, in mixed solvents of the partner and CCl4. The experiments measure the time evolution of the 2D spectra of the hydroxyl (OD) stretching mode of the phenol. The time-dependent 2D spectra are analyzed using time-dependent diagrammatic perturbation theory with a model that includes the chemical exchange (formation and dissociation of the complexes), spectral diffusion of both the complex and the free phenol, orientational relaxation of the complexes and free phenol, and the vibrational lifetimes. The detailed calculations are able to reproduce the experimental results and demonstrate that a method employed previously that used a kinetic model for the volumes of the peaks is adequate to extract the exchange kinetics. The current analysis also yields the spectral diffusion (time evolution of the dynamic line widths) and shows that the spectral diffusion is significantly different for phenol complexes and free phenol.     

Spatial hydration structures and dynamics of phenol in sub- and supercritical water

The hydration structures and dynamics of phenol in aqueous solution at infinite dilution are investigated using molecular-dynamics simulation technique. The simulations are performed at several temperatures along the coexistence curve of water up to the critical point, and above the critical point with density fixed at 0.3 g/cm3. The hydration structures of phenol are characterized using the radial, cylindrical, and spatial distribution functions. In particular, full spatial maps of local atomic (solvent) density around a solute molecule are presented. It is demonstrated that in addition to normal H bonds with hydroxyl group of phenol, water forms pi-type complexes with the center of the benzene ring, in which H2O molecules act as H-bond donor. At ambient conditions phenol is solvated by 38 water molecules, which make up a large hydrophobic cavity, and forms on average 2.39 H bonds (1.55 of which are due to the hydroxyl group-water interactions and 0.84 are due to the pi complex) with its hydration shell. As temperature increases, the hydration structure of phenol undergoes significant changes. The disappearance of the pi-type H bonding is observed near the critical point. Self-diffusion coefficients of water and phenol are also calculated. Dramatic increase in the diffusivity of phenol in aqueous solution is observed near the critical point of simple point-charge-extended water and is related to the changes in water structure at these conditions.     

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.     

Solvent effects on hydrogen bond formation

Enthalpies of solution have been used to calculate transfer enthalpies for phenol, pyridine, and DMSO between the solvent cyclohexane and the solvents CCl₄, benzene, and CHCl₃. By use of model compounds, enthalpies due to interactions with phenol, pyridine, and DMSO have been determined. These enthalpies are used to calculate the effect of solvation relative to cyclohexane on hydrogen bonded complexes in CCl₄ and benzene solvents. Correlations with enthalpies due to interactions and frequency shifts for the hydroxyl stretch in these solvents have also been made.     

Structure and dynamics of beta-cyclodextrin in aqueous solution at the density-functional tight binding level

Structure and dynamics of beta-cyclodextrin (beta-CyD), a prototype host for inclusion compounds of biological interest, is investigated by means of density-functional based tight-binding molecular dynamics (MD) simulations. The computational protocol is benchmarked against available experimental data and first-principles calculations. Solvent-solute interactions, including the diffusion into and dwell time of the solvent in the cavity of beta-CyD, are studied with a hybrid QM/MM method. Comparison of MD simulations of beta-CyD in the gas phase and in water shows that the solvent reduces the flexibility of the structure framework, while the terminal hydroxyl groups become more flexible and are embedded in a network of hydrogen bonds. Our 160 ps MD simulations, provide enough sampling to discuss the dynamics of the water inside the cavity. The dwell time of the encapsulated water molecule has a wide distribution with a peak at 70 fs. Surprisingly, despite only the 17% difference between the "top" and "bottom" opening area of the beta-CyD cone, 64% of the water molecules enter the cavity through the slightly bigger "bottom" aperture.     

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.     

Circular dichroism and conformational dynamics of cephams and their carba and oxa analogues

The biological activity of bicyclic beta-lactam antibiotics depends strongly on the absolute configuration of the bridgehead carbon atom. Frelek and co-workers proposed an empirical helicity rule relating the configuration of the bridgehead carbon atom to the sign of the 220 nm band in the electronic circular dichroism (CD) spectrum of beta-lactams. Here we use synthetic organic chemistry, CD spectroscopy, and time-dependent density functional theory (TDDFT) to investigate the validity of this structure-property relationship for eight model compounds. For conformationally flexible beta-lactams, substantial thermal effects are found which must be included in calculations. To this end, we combine TDDFT calculations of CD with full quantum-mechanical Born-Oppenheimer molecular dynamics (MD) simulations for the first time. The CD spectra are sampled with ground-state density functional trajectories of up to 60 ps. The MD simulations show a surprisingly high sensitivity of the CD to the molecular conformation. On the other hand, the relation between CD and thermally averaged structural parameters is much less complex. While the helicity rule does not seem to hold for individual conformers, it is confirmed by the calculations for seven out of eight systems studied if thermally averaged CD spectra and structures are considered. Since thermal effects on CD can be larger than typical inherent inaccuracies of TDDFT, our results emphasize the need for a systematic treatment of conformational dynamics in CD calculations even for moderately flexible systems. Temperature-dependent CD measurements are very useful for this purpose. Our results also suggest that CD spectroscopy may be used as a sensitive probe of conformational dynamics if combined with electronic structure calculations.     

Microsolvation effects on the excited-state dynamics of protonated tryptophan

To better understand the complex photophysics of the amino acid tryptophan, which is widely used as a probe of protein structure and dynamics, we have measured electronic spectra of protonated, gas-phase tryptophan solvated with a controlled number of water molecules and cooled to approximately 10 K. We observe that, even at this temperature, the bare molecule exhibits a broad electronic spectrum, implying ultrafast, nonradiative decay of the excited state. Surprisingly, the addition of two water molecules sufficiently lengthens the excited-state lifetime that we obtain a fully vibrationally resolved electronic spectrum. Quantum chemical calculations at the RI-CC2/aug-cc-pVDZ level, together with TDDFT/pw based first-principles MD simulations of the excited-state dynamics, clearly demonstrate how interactions with water destabilize the photodissociative states and increase the excited-state lifetime.     

Elucidating the vibrational spectra of hydrogen-bonded aggregates in solution: electronic structure calculations with implicit solvent and first-principles molecular dynamics simulations with explicit solvent for 1-hexanol in n-hexane

Fourier transform infrared spectroscopy is a popular method for the experimental investigation of hydrogen-bonded aggregates, but linking spectral information to microscopic information on aggregate size distribution and aggregate architecture is an arduous task. Static electronic structure calculations with an implicit solvent model, Car-Parrinello molecular dynamics (CPMD) using the Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and classical molecular dynamics simulations for the all-atom version of the optimized parameters for liquid simulations (OPLS-AA) force field were carried out for an ensemble of 1-hexanol aggregates solvated in n-hexane. The initial configurations for these calculations were size-selected from a distribution of aggregates obtained from a large-scale Monte Carlo simulation. The vibrational spectra computed from the static electronic structure calculations for monomers and dimers and from the CPMD simulations for aggregates up to pentamers demonstrate the extent of the contribution of dangling or nondonating hydroxyl groups found in linear and branched aggregates to the "monomeric" peak. Furthermore, the computed spectra show that there is no simple relationship between peak shift and aggregate size nor architecture, but the effect of hydrogen-bond cooperativity is shown to differentiate polymer-like (cooperative) and dimer-like (noncooperative) hydrogen bonds in the vibrational spectrum. In contrast to the static electronic structure calculations and the CPMD simulations, the classical molecular dynamics simulations greatly underestimate the vibrational peak shift due to hydrogen bonding.     

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.     

Non-aqueous solvation of n-octanol and ethanol: spectroscopic and computational studies

Raman spectroscopy was used to examine the interactions of the free O-H bonds in n-octanol and ethanol with the organic solvents carbon tetrachloride (CCl(4)), cyclohexane, and benzene. These spectra reveal that the solvents CCl(4) and cyclohexane have a small effect on the free O-H peak of alcohols, whereas benzene as a solvent significantly red-shifts the free O-H band. Calculated spectra were generated via MP2/6-31G* calculations and the B3LYP/6-31+G**//MP2/6-31G*-derived Boltzmann populations of each ethanol complex and are consistent with the experimental results. Additional spectra were calculated using Boltzmann populations derived from single-point energies at the polarizable continuum model (PCM) level with the B3LYP/6-31+G** level of theory to take overall solvent effects into account, and these simulated spectra are also largely consistent with the experimental results. Analysis of the computational results reveals a lengthening of the O-H bond from the O-H interaction with the delocalized electronic structure of benzene as well as a bimodal distribution of the free O-H peak of the alcohol/benzene mixtures due to two distinctly different types of alcohol/benzene complexes.     

Accidental vibrational degeneracy in vibrational excited states observed with ultrafast two-dimensional IR vibrational echo spectroscopy

The coupling between the OD stretch v=2 level and benzene-ring modes in 2-methoxyphenol-OD (hydroxyl H replaced by D) is observed with ultrafast two-dimensional (2D) IR vibrational echo spectroscopy. Because of this coupling, the 1-2 transition peak in the 2D spectrum is split into a doublet with peaks of approximately equal amplitudes. Several molecules and solvents were used to study this phenomenon. Near-IR (NIR) spectroscopy measurements and density-functional theory calculations (B3LYP6-31+G(d,p) level) were also applied. Experimental results and calculations show that the OD stretch 1-2 transition is coupled to a combination band related to the benzene-ring motions. A simple quantum-mechanical model indicates that the combination band has a frequency of 5172 and 5176.5 cm(-1) in CCl4 and hexane, respectively. The transition between this combination band and the ground state is too weak to detect by NIR. The transition between this band and the OD stretch first excited state is also so weak that most of the intensity of the doublet comes from the oscillator strength produced by coupling to the OD stretch. The model gives the coupling strengths as 6.5 and 7 cm(-1) in CCl4 and hexane, respectively.     

Extensive Numerical Tests of Leapfrog Integrator in Middle Thermostat Scheme in Molecular Simulations

Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics (MD) simulations. Despite the massive use of the conventional leapfrog (LF) integrator in modern computational tools within the framework of MD propagation, further development for better performance is still possible. The alternative version of LF in the middle thermostat scheme (LF-middle) achieves a higher order of accuracy and efficiency and maintains stable dynamics even with the integration time stepsize extended by several folds. In this work, we perform a benchmark test of the two integrators (LF and LF-middle) in extensive conventional and enhanced sampling simulations, aiming at quantifying the time-stepsize-induced variations of global properties (e.g., detailed potential energy terms) as well as of local observables (e.g., free energy changes or bond-lengths) in practical simulations of complex systems. The test set is composed of six chemically and biologically relevant systems, including the conformational change of dihedral flipping in the N-methylacetamide and an AT (Adenine-Thymine) tract, the intra-molecular proton transfer inside malonaldehyde, the binding free energy calculations of benzene and phenol targeting T4 lysozyme L99A, the hydroxyl bond variations in ethaline deep eutectic solvent, and the potential energy of the blue-light using flavin photoreceptor. It is observed that the time-step-induced error is smaller for the LF-middle scheme. The outperformance of LF-middle over the conventional LF integrator is much more significant for global properties than local observables. Overall, the current work demonstrates that the LF-middle scheme should be preferably applied to obtain accurate thermodynamics in the simulation of practical chemical and biological systems.     

Solvent effects on electronic structures and chain conformations of alpha-oligothiophenes in polar and apolar solutions

Solvent effects on electronic structures and chain conformations of alpha-oligothiophenes nTs (n = 1 to 10) are investigated in solvents of n-hexane, 1,4-dioxane, carbon tetrachloride, chloroform, and water by using density functional theory (DFT) and molecular dynamics (MD) simulations. Both implicit and explicit solvent models are employed. The polarized continuum model (PCM) calculations and MD simulations demonstrate the weak solvent effects on the electronic structures of alpha-oligothiophenes. The lowest dipole-allowed vertical excitation energies of nTs, obtained from time-dependent DFT/PCM calculations at the B3LYP/6-31G(d) level, exhibit a red shift as the solvent polarity increases, in agreement with experiments. The studied solvents have little impact on the state order of the low-lying excited states provided that the nTs are kept in C2h or C2v symmetry. The MD simulations demonstrate that the chain conformations are distorted to some extent in polar and nonpolar solvents. A qualitative picture of the distribution of solvent molecules around the solvated nTs is drawn by means of radial and spatial distribution functions. The S...H-O and pi...H-O solute-solvent interactions are insignificant in aqueous solution.     

Calculations of molecules,clusters, and solids with a simplified LCAO-DFT-LDA scheme

A simplified LCAO-DFT-LDA scheme for calculations of structure and electronic structure of large molecules, clusters, and solids is presented. Forces on the atoms are calculated in a semiempirical way considering the electronic states. The small computational effort of this treatment allows one to perform molecular dynamics (MD ) simulations of molecules and clusters up to a few hundred atoms as well as corresponding simulations of condensed systems within the Born-Oppenheimer approximation. The accuracy of the method is illustrated by the results of calculations for a series of small molecules and clusters. © 1996 John Wiley & Sons, Inc.     

Molecular dynamics simulations of helix-forming, amine-functionalized m-poly(phenyleneethynylene)s

We present here the results of all-atom and united-atom molecular dynamics (MD) simulations that were used to examine the folding behavior of an amine-functionalized m-poly(phenyleneethynylene) (m-PPE) oligomer in aqueous environment. The parallelized GROMACS MD simulation code and OPLS force field were used for multiple MD simulations of m-PPE oligomers containing 24 phenyl rings in extended, coiled and helix conformations separately in water to determine the minimum energy conformation of the oligomer in aqueous solvent and what interactions are most important in determining this structure. Simulation results showed that the helix is the preferred minimum energy conformation of a single oligomer in water and that Lennard-Jones interactions are the dominant forces for the stabilization of the helix. In addition, these solvophobic interactions are strong enough to maintain the helix conformation at temperatures up to 523 K.     

Structure and dynamics of the solvated electron in alcohols from resonance Raman spectroscopy

Resonance Raman (RR) spectroscopy is used to probe the structure and excited-state dynamics of the solvated electron in the primary liquid alcohols methanol (MeOH), ethanol (EtOH), n-propanol (n-PrOH), and n-butanol (n-BuOH). The strong resonance enhancements (>or=10(4) relative to pure solvent) of the libration, CO stretch, COH bend, CH3 bend, CH2 bend, and OH stretch reveal significant Franck-Condon coupling of the intermolecular and intramolecular vibrational modes of the solvent to the electronic excitation of the solvated electron. All enhanced bands are fully accounted for by a model of the solvated electron that is comprised of several nearby solvent molecules that are only perturbed by the presence of the electron; no new molecular species are required to explain our data. The 340 cm(-1) downshift observed for the OH stretch frequency of e-(MeOH), relative to pure solvent, strongly suggests that the methanol molecules in the first solvent shell have the hydroxyl group directed linearly toward the excess electron density. The smaller downshifts observed for e-(EtOH), e-(n-PrOH), and e-(n-BuOH) are explained in terms of a OH group that is bent 28-40 degrees from linear. The Raman cross sections and absorption spectra are modeled, lending quantitative support for the inhomogeneous origin of the broad absorption spectra, the necessity of OH local motion in all enhanced Raman modes of the alcohols, and the dominant librational response of the solvent upon photoexcitation of the electron.     

Local order in aqueous NaCl solutions and pure water: X-ray scattering and molecular dynamics simulations study

The microstructures of pure water and aqueous NaCl solutions over a wide range of salt concentrations (0-4 m) under ambient conditions are characterized by X-ray scattering and molecular dynamics (MD) simulations. MD simulations are performed with the rigid SPC water model as a solvent, while the ions are treated as charged Lennard-Jones particles. Simulated data show that the first peaks in the O...O and O...H pair correlation functions clearly decrease in height with increasing salt concentration. Simultaneously, the location of the second O...O peak, the signature of the so-called tetrahedral structure of water, gradually disappears. Consequently, the degree of hydrogen bonding in liquid water decreases when compared to pure fluid. MD results also show that the hydration number around the cation decreases as the salt concentration increases, which is most likely because some water molecules in the first hydration shell are occasionally substituted by chlorine. In addition, the fraction of contact ion pairs increases and that of solvent-separated ion pairs decreases. Experimental data are analyzed to deduce the structure factors and the pair correlation functions of each system. X-ray results clearly show a perturbation of the association structure of the solvent and highlight the appearance of new interactions between ions and water. A model of intermolecular arrangement via MD results is then proposed to describe the local order in each system, as deduced from X-ray scattering data.