Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics

Reaction path ensemble of the B-Z-DNA transition: a comprehensive atomistic study

Since its discovery in 1979, left-handed Z-DNA has evolved from an in vitro curiosity to a challenging DNA structure with crucial roles in gene expression, regulation and recombination. A fundamental question that has puzzled researchers for decades is how the transition from B-DNA, the prevalent right-handed form of DNA, to Z-DNA is accomplished. Due to the complexity of the B–Z-DNA transition, experimental and computational studies have resulted in several different, apparently contradictory models. Here, we use molecular dynamics simulations coupled with state-of-the-art enhanced sampling techniques operating through non-conventional reaction coordinates, to investigate the B–Z-DNA transition at the atomic level. Our results show a complex free energy landscape, where several phenomena such as over-stretching, unpeeling, base pair extrusion and base pair flipping are observed resulting in interconversions between different DNA conformations such as B-DNA, Z-DNA and S-DNA. In particular, different minimum free energy paths allow for the coexistence of different mechanisms (such as zipper and stretch–collapse mechanisms) that previously had been proposed as independent, disconnected models. We find that the B–Z-DNA transition—in absence of other molecular partners—can encompass more than one mechanism of comparable free energy, and is therefore better described in terms of a reaction path ensemble.

Proton transport in functionalised additives for PEM fuel cells: contributions from atomistic simulations

The conventional polymer electrolyte membrane (PEM) materials for fuel cell applications strongly rely on temperature and pressure conditions for optimal performance. In order to expand the range of operating conditions of these conventional PEM materials, mesoporous functionalised SiO(2) additives are developed. It has been demonstrated that these additives themselves achieve proton conductivities approaching those of conventional materials. However, the proton conduction mechanisms and especially factors influencing charge carrier mobility under different hydration conditions are not well known and difficult to separate from concentration effects in experiments. This tutorial review highlights contributions of atomistic computer simulations to the basic understanding and eventual design of these materials. Some basic introduction to the theoretical and computational framework is provided to introduce the reader to the field, the techniques are in principle applicable to a wide range of other situations as well. Simulation results are directly compared to experimental data as far as possible.

Mechanistic insights into the hydrolysis of a nucleoside triphosphate model in neutral and acidic solution

Nucleoside triphosphate hydrolysis is an essential component of all living systems. Despite extensive research, the exact modus and mechanism of this ubiquitous reaction still remain elusive. In this work, we examined the detailed hydrolysis mechanisms of a model nucleoside triphosphate in acidic and neutral solution by means of ab initio simulations. The timescale of the reaction was accessed through use of an accelerated sampling method, metadynamics. Both hydrolyses were found to proceed via different mechanisms; the acidic system reacted by means of concerted general acid catalysis (found to be a so-called D(N)A(N)A(H)D(xh) mechanism), whereas the neutral system reacted by way of a different mechanism (namely, D(N)*A(N)D(xh)A(H)). A neighboring water molecule took on the role of a general base in both systems, which has not been seen before but is a highly plausible reaction path, meaning that substrate-assisted catalysis was not observed in the bulk water environment.

The reaction mechanisms of heme catalases: an atomistic view by ab initio molecular dynamics

Catalases are ubiquitous enzymes that prevent cell oxidative damage by degrading hydrogen peroxide to water and oxygen (2H(2)O(2) → 2H(2)O+O(2)) with high efficiency. The enzyme is first oxidized to a high-valent iron intermediate, known as Compound I (Cpd I, Por(·+)-Fe(IV)=O) which, at difference from other hydroperoxidases, is reduced back to the resting state by further reacting with H(2)O(2). The normal catalase activity is reduced if Cpd I is consumed in a competing side reaction, forming a species named Cpd I*. In recent years, Density Functional Theory (DFT) methods have unraveled the electronic configuration of these high-valent iron species, helping to assign the intermediates trapped in the crystal structures of oxidized catalases. It has been demonstrated that the a priori assumption that the H(+)/H(-) type of mechanism for Cpd I reduction leads to the generation of singlet oxygen is not justified. Moreover, it has been shown by ab initio metadynamics simulations that two pathways are operative for Cpd I reduction: a His-mediated mechanism (described as H·/H(+) + e(-)) in which the distal His acts as an acid-base catalyst and a direct mechanism (described as H·/H·) in which the distal His does not play a direct role. Independently of the mechanism, the reaction proceeds by two one-electron transfers rather than one two-electron transfer, as previously assumed. Electron transfer to Cpd I, regardless of whether the electron is exogenous or endogenous, facilitates protonation of the oxoferryl group, to the point that formation of Cpd I* may be controlled by the easiness of protonation of reduced Cpd I.

The nature of [PdCl(2)(C(2)H(4))(H(2)O)] as an active species in the Wacker process: new insights from ab initio molecular dynamics simulations

First-principles molecular dynamics coupled with metadynamics have been used to gain a deeper insight into the reaction mechanism of the Wacker process by determining the nature of the active species. An explicit and dynamic representation of the aqueous solvent, which was essential for modeling this reaction, was efficiently included into the simulations. Prompted by our earlier results, which showed that the configuration of the catalytically active species [PdCl(2)(H(2)O)(C(2)H(4))] was crucial in the subsequent steps of the Wacker process, herein we focused on the preceding equilibria that led to the formation of both the cis and trans isomers. Starting from the initial catalyst, [PdCl(4)](2-), the free-energy barriers for the forward and backward reactions were calculated. These results confirmed the relevance of the trans intermediate in the reaction mechanism, whilst conversely, they showed that the cis configuration played no role in it. This sole participation of the trans intermediate has some very important implications; besides the mechanistic interpretation of the initial steps in the Wacker reaction mechanism, the analysis of these equilibria provides additional information about the chemical nature of these ligand-substitution processes.

A first-principles theoretical approach to heterogeneous nanocatalysis

A theoretical approach to heterogeneous catalysis by sub-nanometre supported metal clusters and alloys is presented and discussed. Its goal is to perform a computational sampling of the reaction paths in nanocatalysis via a global search in the phase space of structures and stoichiometry combined with filtering which takes into account the given experimental conditions (catalytically relevant temperature and reactant pressure), and corresponds to an incremental exploration of the disconnectivity diagram of the system. The approach is implemented and applied to the study of propylene partial oxidation by Ag(3) supported on MgO(100). First-principles density-functional theory calculations coupled with a Reactive Global Optimization algorithm are performed, finding that: (1) the presence of an oxide support drastically changes the potential energy landscape of the system with respect to the gas phase, favoring configurations which interact positively with the electrostatic field generated by the surface; (2) the reaction energy barriers for the various mechanisms are crucial in the competition between thermodynamically and kinetically favored reaction products; (3) a topological database of structures and saddle points is produced which has general validity and can serve for future studies or for deriving general trends; (4) the MgO(100) surface captures some major features of the effect of an oxide support and appears to be a good model of a simple oxide substrate; (5) strong cooperative effects are found in the co-adsorption of O(2) and other ligands on small metal clusters. The proposed approach appears as a viable route to advance the role of predictive computational science in the field of heterogeneous nanocatalysis.

Simulating microbial systems: addressing model uncertainty/incompleteness via multiscale and entropy methods

Most systems of interest in the natural and engineering sciences are multiscale in character. Typically available models are incomplete or uncertain. Thus, a probabilistic approach is required. We present a deductive multiscale approach to address such problems, focusing on virus and cell systems to demonstrate the ideas. There is usually an underlying physical model, all factors in which (e.g., particle masses, charges, and force constants) are known. For example, the underlying model can be cast in terms of a collection of N-atoms evolving via Newton's equations. When the number of atoms is 10(6) or more, these physical models cannot be simulated directly. However, one may only be interested in a coarse-grained description, e.g., in terms of molecular populations or overall system size, shape, position, and orientation. The premise of this chapter is that the coarse-grained equations should be derived from the underlying model so that a deductive calibration-free methodology is achieved. We consider a reduction in resolution from a description for the state of N-atoms to one in terms of coarse-grained variables. This implies a degree of uncertainty in the underlying microstates. We present a methodology for modeling microbial systems that integrates equations for coarse-grained variables with a probabilistic description of the underlying fine-scale ones. The implementation of our strategy as a general computational platform (SimEntropics™) for microbial modeling and prospects for developments and applications are discussed.

Substrate binding process and mechanistic functioning of type 1 11β-hydroxysteroid dehydrogenase from enhanced sampling methods

In humans, type 1 11β-hydroxysteroid dehydrogenase (11β-HSD-1) plays a key role in the regulation of the glucocorticoids balance by converting the inactive hormone cortisone into cortisol. Numerous functional aspects of 11β-HSD-1 have been understood thanks to the availability at the Worldwide Protein Data Bank of a number of X-ray structures of the enzyme either alone or in complex with inhibitors, and to several experimental data. However at present, a complete description of the dynamic behaviour of 11β-HSD-1 upon substrate binding is missing. To this aim we firstly docked cortisone into the catalytic site of 11β-HSD-1 (both wild type and Y177A mutant), and then we used steered molecular dynamics and metadynamics to simulate its undocking. This methodology helped shedding light at molecular level on the complex relationship between the enzyme and its natural substrate. In particular, the work highlights a) the reason behind the functional dimerisation of 11β-HSD-1, b) the key role of Y177 in the cortisone binding event, c) the fine tuning of the active site degree of solvation, and d) the role of the S228-P237 loop in ligand recognition.

A Combined Metadynamics and Umbrella Sampling Method for the Calculation of Ion Permeation Free Energy Profiles

Free energy calculations are one of the most useful methods for the study of ion transport mechanisms through confined spaces such as protein ion channels. Their reliability depends on a correctly defined reaction coordinate (RC). A straight line is usually not a proper RC for such complicated processes so in this work a combined metadynamics/umbrella sampling (MTD/US) method is proposed. In the combined method, the ion transport pathway is first identified by the MTD method and then the free energy profile or potential of mean force (PMF) along the path is calculated using umbrella sampling. This combined method avoids the discontinuity problem often associated with normal umbrella sampling calculations that assume a straight line RC and it provides a more physically accurate PMF for such processes. The method is demonstrated for the proton transport process through the protein channel of aquaporin-1.

Studies of glutathione transferase P1-1 bound to a platinum(IV)-based anticancer compound reveal the molecular basis of its activation

Platinum-based cancer drugs, such as cisplatin, are highly effective chemotherapeutic agents used extensively for the treatment of solid tumors. However, their effectiveness is limited by drug resistance, which, in some cancers, has been associated with an overexpression of pi class glutathione S-transferase (GST P1-1), an important enzyme in the mercapturic acid detoxification pathway. Ethacraplatin (EA-CPT), a trans-Pt(IV) carboxylate complex containing ethacrynate ligands, was designed as a platinum cancer metallodrug that could also target cytosolic GST enzymes. We previously reported that EA-CPT was an excellent inhibitor of GST activity in live mammalian cells compared to either cisplatin or ethacrynic acid. In order to understand the nature of the drug-protein interactions between EA-CPT and GST P1-1, and to obtain mechanistic insights at a molecular level, structural and biochemical investigations were carried out, supported by molecular modeling analysis using quantum mechanical/molecular mechanical methods. The results suggest that EA-CPT preferentially docks at the dimer interface at GST P1-1 and subsequent interaction with the enzyme resulted in docking of the ethacrynate ligands at both active sites (in the H-sites), with the Pt moiety remaining bound at the dimer interface. The activation of the inhibitor by its target enzyme and covalent binding accounts for the strong and irreversible inhibition of enzymatic activity by the platinum complex.

Free energy of conformational transition paths in biomolecules: the string method and its application to myosin VI

A set of techniques developed under the umbrella of the string method is used in combination with all-atom molecular dynamics simulations to analyze the conformation change between the prepowerstroke (PPS) and rigor (R) structures of the converter domain of myosin VI. The challenges specific to the application of these techniques to such a large and complex biomolecule are addressed in detail. These challenges include (i) identifying a proper set of collective variables to apply the string method, (ii) finding a suitable initial string, (iii) obtaining converged profiles of the free energy along the transition path, (iv) validating and interpreting the free energy profiles, and (v) computing the mean first passage time of the transition. A detailed description of the PPS↔R transition in the converter domain of myosin VI is obtained, including the transition path, the free energy along the path, and the rates of interconversion. The methodology developed here is expected to be useful more generally in studies of conformational transitions in complex biomolecules.

Order parameters for macromolecules: application to multiscale simulation

Order parameters (OPs) characterizing the nanoscale features of macromolecules are presented. They are generated in a general fashion so that they do not need to be redesigned with each new application. They evolve on time scales much longer than 10(-14) s typical for individual atomic collisions/vibrations. The list of OPs can be automatically increased, and completeness can be determined via a correlation analysis. They serve as the basis of a multiscale analysis that starts with the N-atom Liouville equation and yields rigorous Smoluchowski/Langevin equations of stochastic OP dynamics. Such OPs and the multiscale analysis imply computational algorithms that we demonstrate in an application to ribonucleic acid structural dynamics for 50 ns.

First-principles molecular dynamics study on the atomistic behavior of His503 in bovine cytochrome c oxidase

We report first-principles molecular dynamics calculations based on density functional theory performed on the entrance part of the D-path pathway in bovine cytochrome c oxidase. Our models, which are extracted from the fully reduced and oxidized X-ray structures, include His503 as a protonatable site. We find that the protonated His503 with the deprotonated Asp91 [H503-N(δ1)H(+) and D91-C(γ)OO(γ)] are more energetically favorable than other protonation states, [H503-N(δ1) and D91-C(γ)OOH], with an energy difference of about -5kcal/mol in reduced case, while the [H503-N(δ1)H+ and D91-C(γ)OO(-)] state is energetically unstable, about +3kcal/mol higher in energy in the oxidized case. The local interaction of His503 with the surrounding polar residues is necessary and sufficient for determining the energetics. The redox-coupled rotation of His503 is found to change the energetics of the protonation states. We also find that this rotation is coupled with the proton transfer from His503 and Asp91, which leads to the transition between the two different protonation states. This study suggests that His503 is involved in the proton supply to the D-path as a proton acceptor and that the redox-controlled proton-transfer-coupled rotation of His503 is a key process for an effective proton supply to the D-path from water bulk. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Water exchange of a ProHance MRI contrast agent: isomer-dependent free-energy landscapes and mechanisms

The water-exchange reaction in two diastereoisomers of the clinical magnetic resonance imaging contrast agent [Gd(HP-DO3A)(H(2)O)] (also known as ProHance) has been studied using ab initio simulations. On the basis of the molecular-level details of the mechanism derived from these simulations in aqueous solution, we unravel the underlying difference in the free energies and mechanisms of water exchange in the two diastereoisomers. These findings reveal the crucial role played by hydrogen-bonding dynamics and thus suggest their appropriate control in tailoring improved gadolinium-based constrast agents.

An AIMD study of the CPD repair mechanism in water: reaction free energy surface and mechanistic implications

In a series of two papers, we report the detailed mechanism of cyclobutane pyrimidine dimer repair in aqueous solvent using ab initio molecular dynamics simulations (AIMD). Umbrella sampling is used to determine the free energy surface for dimer splitting. The two-dimensional free energy surface for splitting of the C5-C5' and C6-C6' bonds on the anion surface is reported. The splitting of the C5-C5' and C6-C6' bonds occurs on a picosecond time scale. The transition state along the splitting coordinate in the anion state coincides with a maximum in the free energy along the same coordinate on the neutral surface. The implication is that back electron transfer occurring before the anion reaches the transition state leads to reformation of the cyclobutane dimer, while back electron transfer after transit through the transition state, leads to successful repair. On the basis of our calculations for CPD splitting in water, we propose a framework for understanding how various factors, such as solvent polarity, can control repair efficiency. This framework explains why back electron transfer leads predominantly to unsuccessful repair in some situations, and successful repair in others. A key observation is that the same free energy surfaces that control dimer splitting also govern how the back electron transfer rate changes during the splitting process. Configurational changes of the dimer along the splitting coordinate are also documented.

All-atom models of the membrane-spanning domain of HIV-1 gp41 from metadynamics

The 27-residue membrane-spanning domain (MSD) of the HIV-1 glycoprotein gp41 bears conserved sequence elements crucial to the biological function of the virus, in particular a conserved GXXXG motif and a midspan arginine. However, structure-based explanations for the roles of these and other MSD features remain unclear. Using molecular dynamics and metadynamics calculations of an all-atom, explicit solvent, and membrane-anchored model, we study the conformational variability of the HIV-1 gp41 MSD. We find that the MSD peptide assumes a stable tilted α-helical conformation in the membrane. However, when the side chain of the midspan Arg (694) "snorkels" to the outer leaflet of the viral membrane, the MSD assumes a metastable conformation where the highly-conserved N-terminal core (between Lys(681) and Arg(694) and containing the GXXXG motif) unfolds. In contrast, when the Arg(694) side chain snorkels to the inner leaflet, the MSD peptide assumes a metastable conformation consistent with experimental observations where the peptide kinks at Phe(697) to facilitate Arg(694) snorkeling. Both of these models suggest specific ways that gp41 may destabilize viral membrane, priming the virus for fusion with a target cell.

Aggregation-induced chemical reactions: acid dissociation in growing water clusters

Understanding chemical reactivity at ultracold conditions, thus enabling molecular syntheses via interstellar and atmospheric processes, is a key issue in cryochemistry. In particular, acid dissociation and proton transfer reactions are ubiquitous in aqueous microsolvation environments. Here, the full dissociation of a HCl molecule upon stepwise solvation by a small number of water molecules at low temperatures, as relevant to helium nanodroplet isolation (HENDI) spectroscopy, is analyzed in mechanistic detail. It is found that upon successive aggregation of HCl with H(2)O molecules, a series of cyclic heteromolecular structures, up to and including HCl(H(2)O)(3), are initially obtained before a precursor state for dissociation, HCl(H(2)O)(3)···H(2)O, is observed upon addition of a fourth water molecule. The latter partially aggregated structure can be viewed as an "activated species", which readily leads to dissociation of HCl and to the formation of a solvent-shared ion pair, H(3)O(+)(H(2)O)(3)Cl(-). Overall, the process is mostly downhill in potential energy, and, in addition, small remaining barriers are overcome by using kinetic energy released as a result of forming hydrogen bonds due to aggregation. The associated barrier is not ruled by thermal equilibrium but is generated by athermal non-equilibrium dynamics. These "aggregation-induced chemical reactions" are expected to be of broad relevance to chemistry at ultralow temperature much beyond HENDI spectroscopy.

Comparison of free energy surfaces calculations from ab initio molecular dynamic simulations at the example of two transition metal catalyzed reactions

We carried out ab initio molecular dynamic simulations in order to determine the free energy surfaces of two selected reactions including solvents, namely a rearrangement of a ruthenium oxoester in water and a carbon dioxide addition to a palladium complex in carbon dioxide. For the latter reaction we also investigated the gas phase reaction in order to take solvent effects into account. We used two techniques to reconstruct the free energy surfaces: thermodynamic integration and metadynamics. Furthermore, we gave a reasonable error estimation of the computed free energy surface. We calculated a reaction barrier of ΔF = 59.5 ± 8.5 kJ mol(-1) for the rearrangement of a ruthenium oxoester in water from thermodynamic integration. For the carbon dioxide addition to the palladium complex in carbon dioxide we found a ΔF = 44.9 ± 3.3 kJ mol(-1) from metadynamics simulations with one collective variable. The investigation of the same reactions in the gas phase resulted in ΔF = 24.9 ± 6.7 kJ mol(-1) from thermodynamic integration, in ΔF = 26.7 ± 2.3 kJ mol(-1) from metadynamics simulations with one collective variable, and in ΔF = 27.1 ± 5.9 kJ mol(-1) from metadynamics simulations with two collective variables.