Proton transfer 200 years after von Grotthuss: insights from ab initio simulations.

In the last decade, ab initio simulations and especially Car-Parrinello molecular dynamics have significantly contributed to the improvement of our understanding of both the physical and chemical properties of water, ice, and hydrogen-bonded systems in general. At the heart of this family of in silico techniques lies the crucial idea of computing the many-body interactions by solving the electronic structure problem "on the fly" as the simulation proceeds, which circumvents the need for pre-parameterized potential models. In particular, the field of proton transfer in hydrogen-bonded networks greatly benefits from these technical advances. Here, several systems of seemingly quite different nature and of increasing complexity, such as Grotthuss diffusion in water, excited-state proton-transfer in solution, phase transitions in ice, and protonated water networks in the membrane protein bacteriorhodopsin, are discussed in the realms of a unifying viewpoint.     

Structures and vibrational spectra of the sulfur-rich oxides SnO (n = 4-9): the importance of pi*-pi* interactions

The structures of a large number of isomers of the sulfur oxides S(n)O with n = 4-9 have been calculated at the G3X(MP2) level of theory. In most cases, homocyclic molecules with exocyclic oxygen atoms in an axial position are the global minimum structures. Perfect agreement is obtained with experimentally determined structures of S(7)O and S(8)O. The most stable S(4)O isomer as well as some less stable isomers of S(5)O and S(6)O are characterized by a strong pi*-pi* interaction between S==O and S==S groups, which results in relatively long S--S bonds with internuclear distances of 244-262 pm. Heterocyclic isomers are less stable than the global minimum structures, and this energy difference approximately increases with the ring size: 17 (S(4)O), 40 (S(5)O), 32 (S(6)O), 28 (S(7)O), 45 (S(8)O), and 54 kJ mol(-1) (S(9)O). Owing to a favorable pi*-pi* interaction, preference for an axial (or endo) conformation is calculated for the global energy minima of S(7)O, S(8)O, and S(9)O. Vapor-phase decomposition of S(n)O molecules to SO(2) and S(8) is strongly exothermic, whereas the formation of S(2)O and S(8) is exothermic if n<7, but slightly endothermic for S(7)O, S(8)O, and S(9)O. The calculated vibrational spectra of the most stable isomers of S(6)O, S(7)O, and S(8)O are in excellent agreement with the observed data.     

Berry pseudorotation mechanism for the interpretation of the 19F NMR spectrum in PF5 by ab initio molecular dynamics simulations.

For the first time, theoretical evidence that confirms the importance of the Berry pseudorotation process in the interpretation of the 19F NMR spectrum of phosphorus pentafluoride (PF5) is presented. Ab initio molecular dynamics simulations have been performed to generate a large number of configurations used for NMR parameter computations at the density functional theory level. Two different temperatures were set to highlight the effect of pseudorotation process on the NMR spectrum. Average 19F chemical shifts and spin-spin coupling constants calculated for the five fluorine atoms converge towards the NMR equivalence of the five atoms when the Berry pseudorotation mechanism is accounted for.     

The mechanism of formamide hydrolysis in water from ab initio calculations and simulations

The neutral hydrolysis of formamide in water is a suitable reference to quantify the efficiency of proteolytic enzymes. However, experimental data for this reaction has only very recently been obtained and the kinetic constant determined experimentally is significantly higher than that predicted by previous theoretical estimations. In this work, we have investigated in detail the possible mechanisms of this reaction. Several solvent models have been considered that represent a considerable improvement on those used in previous studies. Density functional and ab initio calculations have been carried out on a system which explicitly includes the first solvation shell of the formamide molecule. Its interaction with the bulk has been treated with the aid of a dielectric continuum model. Molecular dynamics simulations at the combined density functional/molecular mechanics level have been carried out in parallel to better understand the structure of the reaction intermediates in aqueous solution. Overall, the most favored mechanism predicted by our study involves two reaction steps. In the first step, the carbonyl group of the formamide molecule is hydrated to form a diol intermediate. The corresponding transition structure involves two water molecules. From this intermediate, a water-assisted proton transfer occurs from one of the hydroxy groups to the amino group. This reaction step may lead either to the formation of a new reaction intermediate with a marked zwitterionic character or to dissociation of the system into ammonia and formic acid. The zwitterionic intermediate dissociates quite easily but its lifetime is not negligible and it could play a role in the hydrolysis of substituted amides or peptides. The predicted pseudo-first-order kinetic constant for the rate-limiting step (the first step) of the hydrolysis reaction at 25 degrees C (3.9x10(-10) s(-1)) is in excellent agreement with experimental data (1.1x10(-10) s(-1)).     

The Jahn-Teller effect of the TiIII ion in aqueous solution: extended ab initio QM/MM molecular dynamics simulations.

Combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, including only the first and the first and second hydration shells in the QM region, were performed for TiIII in aqueous solution. The hydration structure of TiIII is discussed in terms of radial distribution functions, coordination-number distributions and several angle distributions. Dynamical properties, such as librational and vibrational motions and TiIII-O vibrations, were evaluated. A fast dynamical Jahn-Teller effect of TiIII(aq) was observed in the QM/MM simulations, in particular when the second hydration shell was included into the QM region. The results justify the computational effort required for the inclusion of the second hydration shell into the QM region and show the importance of this effort for obtaining accurate hydration-shell geometries, dynamical properties, and details of the Jahn-Teller effect.     

Ab initio static and molecular dynamics study of 4-styrylpyridine.

We report an in‐depth theoretical study of 4‐styrylpyridine in its singlet S₀ ground state. The geometries and the relative stabilities of the trans and cis isomers were investigated within density functional theory (DFT) as well as within Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and coupled cluster (CC) theories. The DFT calculations were performed using the B3LYP and PBE functionals, with basis sets of different qualities, and gave results that are very consistent with each other. The molecular structure is thus predicted to be planar at the energy minimum, which is associated with the trans conformation, and to become markedly twisted at the minimum of higher energy, which is associated with the cis conformation. The results of the calculations performed with the post‐HF methods approach those obtained with the DFT methods, provided that the level of treatment of the electronic correlation is high enough and that sufficiently flexible basis sets are used. Calculations carried out within DFT also allowed the determination of the geometry and the energy of the molecule at the biradicaloid transition state associated with the thermal cis?trans isomerization and at the transition states associated with the enantiomerization of the cis isomer and with the rotations of the pyridinyl and phenyl groups in the trans and cis isomers. Car–Parrinello molecular dynamics simulations were also performed at 50, 150, and 300 K using the PBE functional. The studies allowed us to evidence the highly flexible nature of the molecule in both conformations. In particular, the trans isomer was found to exist mainly in a nonplanar form at finite temperatures, while the rotation of the pyridinyl ring in the cis isomer was incidentally observed to take place within ≈1 ps during the simulation carried out at 150 K on this isomer.     

Exceptional thermodynamic stability of DNA duplexes modified by nonpolar base analogues is due to increased stacking interactions and favorable solvation: Correlated ab initio calculations and molecular dynamics simulations

The geometries of DNA hexamer (5'-GGAACC-3') and DNA 13-mer (5'-GCGTACACATGCG-3') have been determined by molecular dynamics (MD) simulations using an empirical force field. The central canonical base pair was replaced by a pair of nonpolar base analogues, 2,2'-bipyridyl and 3-methylisocarbostyril. The stabilization energy of the model system (model A) consisting of a central base pair (base-analogue pair) and two neighboring base pairs was determined by the RI-MP2 method using an extended aug-cc-pVDZ basis set. The geometry of the model was averaged from structures determined by MD simulations. The role of the solvent was covered by the COSMO continuum solvent model and calculations were performed for a larger model system (model B) which also contained a sugar-phosphate backbone. The total stabilization energies of the unperturbed system and the system perturbed by a base-analogue pair (model A) were comparable to the stability of both duplexes experimentally determined. This is due to large stacking interaction energy of the base-analogue self-pair which compensates for the missing hydrogen-bonding energy of the replaced adenine...thymine base pair. The selectivity of the base-analogue pair was reproduced (model B) when their desolvation energy was included with the interaction energy of both strands determined by the approximate SCC-DFTB-D method.     

Azulene-to-naphthalene rearrangement: the Car-Parrinello metadynamics method explores various reaction mechanisms.

We studied the thermal intramolecular and radical rearrangement of azulene to naphthalene by employing a novel metadynamics method based on Car-Parrinello molecular dynamics. We demonstrate that relatively short simulations can provide us with several possible reaction mechanisms for the rearrangement. We show that different choices of the collective coordinates can steer the reaction along different pathways, thus offering the possibility of choosing the most probable mechanism. We consider herein three intramolecular mechanisms and two radical pathways. We found the norcaradiene pathway to be the preferable intramolecular mechanism, whereas the spiran mechanism is the favored radical route. We obtained high activation energies for all the intramolecular pathways (81.5-98.6 kcal mol(-1)), whereas the radical routes have activation energies of 24-39 kcal mol(-1). The calculations have also resulted in elementary steps and intermediates not yet considered. A few attractive features of the metadynamics method in studying chemical reactions are pointed out.     

Cover Feature: Ab Initio Molecular Dynamics Study of H2 Dissociation Mechanisms on Cu13 and Defective Graphene-supported Cu13 Clusters: Active Sites,Energy Barriers and Adsorption States (ChemPhysChem 19/2023)

Ab initio molecular dynamics calculations were performed to study H₂ dissociation mechanisms on Cu₁₃ and defective graphene-supported Cu₁₃ clusters. The study reveals that seven types of corresponding dissociation processes are found on the two clusters. The average dissociation energy barriers are 0.51 eV on the Cu₁₃ cluster and 0.12 eV on the defective graphene-supported Cu₁₃ cluster, which are lowered by ~19 % and ~81 % compared with the pristine Cu(111) surface, respectively. Furthermore, compared with the pure Cu₁₃ cluster, the average dissociation energy barrier on the defective graphene-supported Cu₁₃ cluster is substantially reduced by about 76 %. The preferred dissociation mechanisms on the two clusters are H₂ located at a top-bridge site with the barrier heights of 0.30 eV on the Cu₁₃ cluster and −0.31 eV on the defective graphene-supported Cu₁₃ cluster_. Analysis of the H−Cu bond interactions in the transition states shows that the antibonding-orbital center shifts upward on the defective graphene-supported Cu₁₃ cluster compared with the one on the Cu₁₃ cluster, which explains the reduction of the dissociation energy barrier. The average adsorption energy of dissociated H atoms is also greatly enhanced on the defective graphene-supported Cu₁₃ cluster, about twice that on the Cu₁₃ cluster.     

Interaction of carboranes with biomolecules: formation of dihydrogen bonds.

Noncovalent interactions of the polyhedral carborane 1-carba-closo-dodecaborane (CB(11)H(12))(-) with building blocks of biomolecules, modelled by glycine (GLY), serine (SER), phenylalanine (PHE), glutamic acid (GLU), lysine (LYS) and arginine (ARG), were investigated in vacuo by molecular dynamics simulations with the UFF empirical potential. Selected structures were further studied by accurate ab initio quantum chemical procedures. Interactions with a peptide bond (GLY-SER dipeptide) and a nucleic acid building block (guanine) were also considered. The RESP and NPA charges of carboranes and small model systems are compared and their use is discussed. The dominant interaction between carboranes and biomolecules is the formation of unconventional proton-hydride hydrogen bonds (dihydrogen bonds) characterized by a short distance between hydrogen atoms (as close as 1.8 A) and an average strength in the range of 4.2-5.8 kcal mol(-1). The total stabilization energy of complexes investigated is rather large, and the largest value (approximately 15 kcal mol(-1)) was found for the carborane complexes with ARG and the GLY-SER dipeptide. These interactions are ubiquitous under geometrical constraints influencing the strength of the interaction. The carborane forms dihydrogen bonds with biomolecules preferably with the hydrogen atoms of its lower hemisphere (i.e. the part of the cage opposite to the carbon atom). These two geometrical factors can be used to explain the specificity of inhibition of HIV protease by carboranes.     

Local structure and thermodynamics of a core-softened potential fluid: theory and simulation.

Phase behavior and structural properties of homogeneous and inhomogeneous core-softened (CS) fluid consisting of particles interacting via the potential, which combines the hard-core repulsion and double attractive well interaction, are investigated. The vapour-liquid coexistence curves and critical points for various interaction ranges of the potential are determined by discrete molecular dynamics simulations to provide guidance for the choice of the bulk density and potential parameters for the study of homogeneous and inhomogeneous structures. Spatial correlations in the homogeneous CS system are studied by the Ornstein-Zernike integral equation in combination with the modified hypernetted chain (MHNC) approximation. The local structure of CS fluid subjected to diverse external fields maintaining the equilibrium with the bulk CS fluid are studied on the basis of a recently proposed third order+second order perturbation density functional approximation (DFA). The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. Reasonable agreement between the results of both methods proves that the DFA theory applied in this work is a convenient theoretical tool for the investigation of the CS fluid, which is practically applicable for modeling numerous real systems.     

Ab initio molecular dynamics study of the keto-enol tautomerism of acetone in solution.

We have studied the keto-enol interconversion of acetone to understand the mechanism of tautomerism relevant to numerous organic and biochemical processes. Applying the ab initio metadynamics method, we simulated the keto-enol isomerism both in the gas phase and in the presence of water. For the gas-phase intramolecular mechanism we show that no other hydrogen-transfer reactions can compete with the simple keto-enol tautomerism. We obtain an intermolecular mechanism and remarkable participation of water when acetone is solvated by neutral water. The simulations reveal that C deprotonation is the kinetic bottleneck of the keto-enol transformation, in agreement with experimental observations. The most interesting finding is the formation of short H-bonded chains of water molecules that provide the route for proton transfer from the carbon to the oxygen atom of acetone. The mechanistic picture that emerged from the present study involves proton migration and emphasizes the importance of active solvent participation in tautomeric interconversion.     

Spectroscopic Signature of Dynamical Instability of the Aqueous Complex in the Brown-Ring Nitrate Test**

The chemistry of the brown-ring test has been investigated for nearly a century. Though recent studies have focused on solid state structure determination and measurement of spectra, mechanistic details and kinetics, the aspects of solution structure and dynamics remain unknown. We have studied structural fluctuations of the brown-ring complex in aqueous solution with ab-initio molecular dynamics simulations, from which we identified that the classically established pseudo-octahedral [Fe(H₂O)₅(NO)]²⁺ complex is present along with a square-pyramidal [Fe(H₂O)₄(NO)]²⁺ complex. Based on the inability in multi-reference calculations to reproduce the experimental UV-vis spectra in aqueous solution by inclusion of thermal fluctuations of the [Fe(H₂O)₅(NO)]²⁺ complex alone, we propose the existence of an equilibrium between pseudo-octahedral and square-pyramidal complexes. Despite challenges in constructing models reproducing the solid-state UV-vis spectrum, the advanced spectrum simulation tool motivates us to challenge the established picture of a sole pseudo-octahedral complex in solution.     

HNO‐Binding in Heme Proteins: Effects of Iron Oxidation State,Axial Ligand,and Protein Environment

HNO plays significant roles in many biological processes. Numerous heme proteins bind HNO, an important step for its biological functions. A systematic computational study was performed to provide the first detailed trends and origins of the effects of iron oxidation state, axial ligand, and protein environment on HNO binding. The results show that HNO binds much weaker with ferric porphyrins than corresponding ferrous systems, offering strong thermodynamic driving force for experimentally observed reductive nitrosylation. The axial ligand was found to influence HNO binding through its trans effect and charge donation effect. The protein environment significantly affects the HNO hydrogen bonding structures and properties. The predicted NMR and vibrational data are in excellent agreement with experiment. This broad range of results shall facilitate studies of HNO binding in many heme proteins, models, and related metalloproteins.     

Binding of Cationic and Neutral Phenanthridine Intercalators to a DNA Oligomer Is Controlled by Dispersion Energy: Quantum Chemical Calculations and Molecular Mechanics Simulations

Correlated ab initio as well as semiempirical quantum chemical calculations and molecular dynamics simulations were used to study the intercalation of cationic ethidium, cationic 5‐ethyl‐6‐phenylphenanthridinium and uncharged 3,8‐diamino‐6‐phenylphenanthridine to DNA. The stabilization energy of the cationic intercalators is considerably larger than that of the uncharged one. The dominant energy contribution with all intercalators is represented by dispersion energy. In the case of the cationic intercalators, the electrostatic and charge‐transfer terms are also important. The ΔG of ethidium intercalation to DNA was estimated at ?4.5 kcal mol^?1 and this value agrees well with the experimental result. Of six contributions to the final free energy, the interaction energy value is crucial. The intercalation process is governed by the non‐covalent stacking (including charge‐transfer) interaction while the hydrogen bonding between the ethidium amino groups and the DNA backbone is less important. This is confirmed by the evaluation of the interaction energy as well as by the calculation of the free energy change. The intercalation affects the macroscopic properties of DNA in terms of its flexibility. This explains the easier entry of another intercalator molecule in the vicinity of an existing intercalation site.     

Complex Reaction Environments and Competing Reaction Mechanisms in Zeolite Catalysis: Insights from Advanced Molecular Dynamics

The methanol‐to‐olefin process is a showcase example of complex zeolite‐catalyzed chemistry. At real operating conditions, many factors affect the reactivity, such as framework flexibility, adsorption of various guest molecules, and competitive reaction pathways. In this study, the strength of first principle molecular dynamics techniques to capture this complexity is shown by means of two case studies. Firstly, the adsorption behavior of methanol and water in H‐SAPO‐34 at 350 °C is investigated. Hereby an important degree of framework flexibility and proton mobility was observed. Secondly, the methylation of benzene by methanol through a competitive direct and stepwise pathway in the AFI topology was studied. Both case studies clearly show that a first‐principle molecular dynamics approach enables unprecedented insights into zeolite‐catalyzed reactions at the nanometer scale to be obtained.     

Understanding the Role of Water in Aqueous Ruthenium‐Catalyzed Transfer Hydrogenation of Ketones.

We report an accurate computational study of the role of water in transfer hydrogenation of formaldehyde with a ruthenium‐based catalyst using a water‐specific model. Our results suggest that the reaction mechanism in aqueous solution is significantly different from that in the gas phase or in methanol solution. Previous theoretical studies have shown a concerted hydride and proton transfer in the gas phase (M. Yamakawa, H. Ito, R. Noyori, J. Am. Chem. Soc. 2000 , 122, 1466–1478;J.‐W. Handgraaf, J. N. H. Reek, E. J. Meijer, Organometallics 2003 , 22, 3150–3157; D. A. Alonso, P. Brandt, S. J. M. Nordin, P. G. Andersson, J. Am. Chem. Soc. 1999 , 121, 9580–9588; D. G. I. Petra, J. N. H. Reek, J.‐W. Handgraaf, E. J. Meijer, P. Dierkes, P. C. J. Kamer, J. Brussee, H. E. Schoemaker, P. W. N. M. van Leeuwen, Chem. Eur. J. 2000 , 6, 2818–2829), whereas a delayed, solvent‐mediated proton transfer has been observed in methanol solution (J.‐W. Handgraaf, E. J. Meijer, J. Am. Chem. Soc. 2007 , 129, 3099–3103). In aqueous solution, a concerted transition state is observed, as in the previous studies. However, only the hydride is transferred at that point, whereas the proton is transferred later by a water molecule instead of the catalyst.