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

Metadynamics combined with auxiliary density functional and density functional tight-binding methods: alanine dipeptide as a case study

Application of ab initio molecular dynamics to study free energy surfaces (FES) is still not commonly performed because of the extensive sampling required. Indeed, it generally necessitates computationally costly simulations of more than several hundreds of picoseconds. To achieve such studies, efficient density functional theory (DFT) formalisms, based on various levels of approximate computational schemes, have been developed, and provide a good alternative to commonly used DFT implementations. We report benchmark results on the conformational change FES of alanine dipeptide obtained with auxiliary density functional theory (ADFT) and second- and third-order density functional tight-binding (DFTB) methods coupled to metadynamics simulations. The influence of an explicit water solvent is also studied with DFTB, which was made possible by its lower computational cost compared to ADFT. Simulations lengths of 2.1 and 15 ns were achieved with ADFT and DFTB, respectively, in a reasonably short computational time. ADFT leads to a free energy difference (ΔF _eq-ax) of ∼ -3 kcal mol^-1 between the two low energy conformers, C_7eq and C_7ax, which is lower by only 1.5 kcal mol^-1 than the ΔF _eq-ax computed with DFTB. The two minima in ADFT FES are separated by an energy barrier of 9 kcal mol^-1, which is higher than the DFTB barriers by 2-4 kcal mol^-1. Despite these small quantitative differences, the DFTB method reveals FES shapes, confor-mation geometries and energies of the stationary points in good agreement with these found with ADFT. This validates the promising applicability of DFTB to FES of reactions occurring in larger-size systems placed in complex environments.

Full atomistic reaction mechanism with kinetics for CO reduction on Cu(100) from ab initio molecular dynamics free-energy calculations at 298 K

A critical step toward the rational design of new catalysts that achieve selective and efficient reduction of CO₂ to specific hydrocarbons and oxygenates is to determine the detailed reaction mechanism including kinetics and product selectivity as a function of pH and applied potential for known systems. To accomplish this, we apply ab initio molecular metadynamics simulations (AIMμD) for the water/Cu(100) system with five layers of the explicit solvent under a potential of -0.59 V [reversible hydrogen electrode (RHE)] at pH 7 and compare with experiment. From these free-energy calculations, we determined the kinetics and pathways for major products (ethylene and methane) and minor products (ethanol, glyoxal, glycolaldehyde, ethylene glycol, acetaldehyde, ethane, and methanol). For an applied potential (U) greater than -0.6 V (RHE) ethylene, the major product, is produced via the Eley-Rideal (ER) mechanism using H₂O + e^- The rate-determining step (RDS) is C-C coupling of two CO, with ΔG^‡ = 0.69 eV. For an applied potential less than -0.60 V (RHE), the rate of ethylene formation decreases, mainly due to the loss of CO surface sites, which are replaced by H*. The reappearance of C₂H₄ along with CH₄ at U less than -0.85 V arises from *CHO formation produced via an ER process of H* with nonadsorbed CO (a unique result). This *CHO is the common intermediate for the formation of both CH₄ and C₂H₄ These results suggest that, to obtain hydrocarbon products selectively and efficiency at pH 7, we need to increase the CO concentration by changing the solvent or alloying the surface.

Reliable and efficient reaction path and transition state finding for surface reactions with the growing string method

The computational challenge of fast and reliable transition state and reaction path optimization requires new methodological strategies to maintain low cost, high accuracy, and systematic searching capabilities. The growing string method using internal coordinates has proven to be highly effective for the study of molecular, gas phase reactions, but difficulties in choosing a suitable coordinate system for periodic systems has prevented its use for surface chemistry. New developments are therefore needed, and presented herein, to handle surface reactions which include atoms with large coordination numbers that cannot be treated using standard internal coordinates. The double-ended and single-ended growing string methods are implemented using a hybrid coordinate system, then benchmarked for a test set of 43 elementary reactions occurring on surfaces. These results show that the growing string method is at least 45% faster than the widely used climbing image-nudged elastic band method, which also fails to converge in several of the test cases. Additionally, the surface growing string method has a unique single-ended search method which can move outward from an initial structure to find the intermediates, transition states, and reaction paths simultaneously. This powerful explorative feature of single ended-growing string method is demonstrated to uncover, for the first time, the mechanism for atomic layer deposition of TiN on Cu(111) surface. This reaction is found to proceed through multiple hydrogen-transfer and ligand-exchange events, while formation of H-bonds stabilizes intermediates of the reaction. Purging gaseous products out of the reaction environment is the driving force for these reactions. © 2017 Wiley Periodicals, Inc.

Fluctuation enhancement of ion diffusivity in liquids

The diffusivity of ions in liquid solutions is known either to decrease with an increase in the ion size or to have a single maximum depending on the ion size. This article presents evidence for the appearance of multiple maxima and thus multiple ion sizes with enhanced diffusivity.

A β-Mannanase with a Lysozyme-like Fold and a Novel Molecular Catalytic Mechanism

The enzymatic cleavage of β-1,4-mannans is achieved by endo-β-1,4-mannanases, enzymes involved in germination of seeds and microbial hemicellulose degradation, and which have increasing industrial and consumer product applications. β-Mannanases occur in a range of families of the CAZy sequence-based glycoside hydrolase (GH) classification scheme including families 5, 26, and 113. In this work we reveal that β-mannanases of the newly described GH family 134 differ from other mannanase families in both their mechanism and tertiary structure. A representative GH family 134 endo-β-1,4-mannanase from a Streptomyces sp. displays a fold closely related to that of hen egg white lysozyme but acts with inversion of stereochemistry. A Michaelis complex with mannopentaose, and a product complex with mannotriose, reveal ligands with pyranose rings distorted in an unusual inverted chair conformation. Ab initio quantum mechanics/molecular mechanics metadynamics quantified the energetically accessible ring conformations and provided evidence in support of a ¹C₄ → ³H₄^‡ → ³S₁ conformational itinerary along the reaction coordinate. This work, in concert with that on GH family 124 cellulases, reveals how the lysozyme fold can be co-opted to catalyze the hydrolysis of different polysaccharides in a mechanistically distinct manner.

Smoothed Biasing Forces Yield Unbiased Free Energies with the Extended-System Adaptive Biasing Force Method

We report a theoretical description and numerical tests of the extended-system adaptive biasing force method (eABF), together with an unbiased estimator of the free energy surface from eABF dynamics. Whereas the original ABF approach uses its running estimate of the free energy gradient as the adaptive biasing force, eABF is built on the idea that the exact free energy gradient is not necessary for efficient exploration, and that it is still possible to recover the exact free energy separately with an appropriate estimator. eABF does not directly bias the collective coordinates of interest, but rather fictitious variables that are harmonically coupled to them; therefore is does not require second derivative estimates, making it easily applicable to a wider range of problems than ABF. Furthermore, the extended variables present a smoother, coarse-grain-like sampling problem on a mollified free energy surface, leading to faster exploration and convergence. We also introduce CZAR, a simple, unbiased free energy estimator from eABF trajectories. eABF/CZAR converges to the physical free energy surface faster than standard ABF for a wide range of parameters.

Conformational changes during permeation of Na+ through a modified cyclic peptide nanotube promote energy landscape roughness

Using a metadynamics approach, we investigate the potential of mean force for Na⁺ permeation inside a cyclic peptide nanotube (CPN) with modified interior as a function of ion position, coordination number, and lumen chemistry. We show that functionalizing the lumen of a CPN with a methyl-benzoic acid group introduces non-periodic variations in the internal energy of the nanotube, which dictate the overall free energy roughness during the permeation of Na⁺. These non-periodic variations arise from the structural dynamics of the functional group, where changes in the dihedral angles induced by the proximity of the ion give rise to conformational changes that increase landscape roughness and thereby decrease transport rate. Our computational framework emphasizes the advantages of using the coordination number as a collective variable to investigate the available conformations during ion permeation through CPNs, and reveals new structure-function relations for chemically tunable CPNs, paving the way for rational design of nano-porous systems with tunable selectivity and flux.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

Conformational Heterogeneity of Bax Helix 9 Dimer for Apoptotic Pore Formation

Helix α9 of Bax protein can dimerize in the mitochondrial outer membrane (MOM) and lead to apoptotic pores. However, it remains unclear how different conformations of the dimer contribute to the pore formation on the molecular level. Thus we have investigated various conformational states of the α9 dimer in a MOM model — using computer simulations supplemented with site-specific mutagenesis and crosslinking of the α9 helices. Our data not only confirmed the critical membrane environment for the α9 stability and dimerization, but also revealed the distinct lipid-binding preference of the dimer in different conformational states. In our proposed pathway, a crucial iso-parallel dimer that mediates the conformational transition was discovered computationally and validated experimentally. The corroborating evidence from simulations and experiments suggests that, helix α9 assists Bax activation via the dimer heterogeneity and interactions with specific MOM lipids, which eventually facilitate proteolipidic pore formation in apoptosis regulation.

Quantum free energy landscapes from ab initio path integral metadynamics: Double proton transfer in the formic acid dimer is concerted but not correlated

With the goal of computing quantum free energy landscapes of reactive (bio)chemical systems in multi-dimensional space, we combine the metadynamics technique for sampling potential energy surfaces with the ab initio path integral approach to treating nuclear quantum motion. This unified method is applied to the double proton transfer process in the formic acid dimer (FAD), in order to study the nuclear quantum effects at finite temperatures without imposing a one-dimensional reaction coordinate or reducing the dimensionality. Importantly, the ab initio path integral metadynamics technique allows one to treat the hydrogen bonds and concomitant proton transfers in FAD strictly independently and thus provides direct access to the much discussed issue of whether the double proton transfer proceeds via a stepwise or concerted mechanism. The quantum free energy landscape we compute for this H-bonded molecular complex reveals that the two protons move in a concerted fashion from initial to product state, yet world-line analysis of the quantum correlations demonstrates that the protons are as quantum-uncorrelated at the transition state as they are when close to the equilibrium structure.

Enzymatic Cleavage of Glycosidic Bonds: Strategies on How to Set Up and Control a QM/MM Metadynamics Simulation

Carbohydrates play crucial roles in many biological processes, from cell-cell adhesion to chemical signaling. Their complexity and diversity, related to α/β anomeric configuration, ring substituents, and conformational variations, require a diverse set of enzymes for their processing. Among them, glycoside hydrolases (GHs) are responsible for the hydrolysis of one of the strongest bonds in nature: the glycosidic bond. These highly specialized biological catalysts select particular conformations their carbohydrate substrates to enhance catalysis. The evolution of this conformation during the reaction of glycosidic bond cleavage, known as the conformational catalytic itinerary, is of fundamental interest in glycobiology, with impact on inhibitor and drug design. Here we review some of the aspects and the main strategies one needs to take into account when simulating a reaction in a GH enzyme using QM/MM metadynamics. Several specific aspects are highlighted, from the importance of the distortion of the substrate at the Michaelis complex to the variable control during the metadynamics simulation or the analysis of the reaction mechanism and conformational itinerary. The increasing speed of computer power and methodological advances have added a vital tool to the study of GH mechanisms, as shown here and recent reviews. It is hoped that this chapter will serve as a first guide for those attempting to perform a metadynamics simulation of these relevant and fascinating enzymes.

The effect of molecular dynamics sampling on the calculated observable gas-phase structures

In this study, we compare the performance of various ab initio molecular dynamics (MD) sampling methods for the calculation of the observable vibrationally-averaged gas-phase structures of benzene, naphthalene and anthracene molecules. Nose-Hoover (NH), canonical and quantum generalized-Langevin-equation (GLE) thermostats as well as the a posteriori quantum correction to the classical trajectories have been tested and compared to the accurate path-integral molecular dynamics (PIMD), static anharmonic vibrational calculations as well as to the experimental gas electron diffraction data. Classical sampling methods neglecting quantum effects (NH and canonical GLE thermostats) dramatically underestimate vibrational amplitudes for the bonded atom pairs, both C-H and C-C, the resulting radial distribution functions exhibit nonphysically narrow peaks. This deficiency is almost completely removed by taking the quantum effects on the nuclei into account. The quantum GLE thermostat and a posteriori correction to the canonical GLE and NH thermostatted trajectories capture most vibrational quantum effects and closely reproduce computationally expensive PIMD and experimental radial distribution functions. These methods are both computationally feasible and accurate and are therefore recommended for calculations of the observable gas-phase structures. A good performance of the quantum GLE thermostat for the gas-phase calculations is encouraging since its parameters have been originally fitted for the condensed-phase calculations. Very accurate molecular structures can be predicted by combining the equilibrium geometry obtained at a high level of electronic structure theory with vibrational amplitudes and corrections calculated using MD driven by a lower level of electronic structure theory.

First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes

Computational chemistry is a valuable aid to complement experimental studies of organometallic systems and their reactivity. It allows probing mechanistic hypotheses and investigating molecular structures, shedding light on the behavior and properties of molecular assemblies at the atomic scale. When approaching a chemical problem, the computational chemist has to decide on the theoretical approach needed to describe electron/nuclear interactions and the composition of the model used to approximate the actual system. Both factors determine the reliability of the modeling study. The community dedicated much effort to developing and improving the performance and accuracy of theoretical approaches for electronic structure calculations, on which the description of (inter)atomic interactions rely. Here, the importance of the model system used in computational studies is highlighted through examples from our recent research focused on organometallic systems and homogeneous catalytic processes. We show how the inclusion of explicit solvent allows the characterization of molecular events that would otherwise not be accessible in reduced model systems (clusters). These include the stabilization of nascent charged fragments via microscopic solvation (notably, hydrogen bonding), transfer of charge (protons) between distant fragments mediated by solvent molecules, and solvent coordination to unsaturated metal centers. Furthermore, when weak interactions are involved, we show how conformational and solvation properties of organometallic complexes are also affected by the explicit inclusion of solvent molecules. Such extended model systems may be treated under periodic boundary conditions, thus removing the cluster/continuum (or vacuum) boundary, and require a statistical mechanics simulation technique to sample the accessible configurational space. First-principles molecular dynamics, in which atomic forces are computed from electronic structure calculations (namely, density functional theory), is certainly the technique of choice to investigate chemical events in solution. This methodology is well established and thanks to advances in both algorithms and computational resources simulation times required for the modeling of chemical events are nowadays accessible, though the computational requirements use to be high. Specific applications reviewed here include mechanistic studies of the Shilov and Wacker processes, speciation in Pd chemistry, hydrogen bonding to metal centers, and the dynamics of agostic interactions.

Impact of Ser17 Phosphorylation on the Conformational Dynamics of the Oncoprotein MDM2

MDM2 is an important oncoprotein that downregulates the activity of the tumor suppressor protein p53 via binding of its N-terminal domain to the p53 transactivation domain. The first 24 residues of the MDM2 N-terminal domain form an intrinsically disordered "lid" region that interconverts on a millisecond time scale between "open" and "closed" states in unliganded MDM2. While the former conformational state is expected to facilitate p53 binding, the latter competes in a pseudo-substrate manner with p53 for its binding site. Phosphorylation of serine 17 in the MDM2 lid region is thought to modulate the equilibrium between "open" and "closed" lid states, but contradictory findings on the favored lid conformational state upon phosphorylation have been reported. Here, the nature of the conformational states of MDM2 pSer17 and Ser17Asp variants was addressed by means of enhanced sampling molecular dynamics simulations. Detailed analyses of the computed lid conformational ensembles indicate that both lid variants stabilize a "closed" state, with respect to wild type. Nevertheless, the nature of the closed-state conformational ensembles differs significantly between the pSer17 and Ser17Asp variants. Thus, care should be applied in the interpretation of biochemical experiments that use phosphomimetic variants to model the effects of phosphorylation on the structure and dynamics of this disordered protein region.

Deacylation Mechanism and Kinetics of Acyl-Enzyme Complex of Class C β-Lactamase and Cephalothin

Understanding the molecular details of antibiotic resistance by the bacterial enzymes β-lactamases is vital for the development of novel antibiotics and inhibitors. In this spirit, the detailed mechanism of deacylation of the acyl-enzyme complex formed by cephalothin and class C β-lactamase is investigated here using hybrid quantum-mechanical/molecular-mechanical molecular dynamics methods. The roles of various active-site residues and substrate in the deacylation reaction are elucidated. We identify the base that activates the hydrolyzing water molecule and the residue that protonates the catalytic serine (Ser64). Conformational changes in the active sites and proton transfers that potentiate the efficiency of the deacylation reaction are presented. We have also characterized the oxyanion holes and other H-bonding interactions that stabilize the reaction intermediates. Together with the kinetic and mechanistic details of the acylation reaction, we analyze the complete mechanism and the overall kinetics of the drug hydrolysis. Finally, the apparent rate-determining step in the drug hydrolysis is scrutinized.

Structural insights into the mechanism of activation of the TRPV1 channel by a membrane-bound tarantula toxin

Venom toxins are invaluable tools for exploring the structure and mechanisms of ion channels. Here, we solve the structure of double-knot toxin (DkTx), a tarantula toxin that activates the heat-activated TRPV1 channel. We also provide improved structures of TRPV1 with and without the toxin bound, and investigate the interactions of DkTx with the channel and membranes. We find that DkTx binds to the outer edge of the external pore of TRPV1 in a counterclockwise configuration, using a limited protein-protein interface and inserting hydrophobic residues into the bilayer. We also show that DkTx partitions naturally into membranes, with the two lobes exhibiting opposing energetics for membrane partitioning and channel activation. Finally, we find that the toxin disrupts a cluster of hydrophobic residues behind the selectivity filter that are critical for channel activation. Collectively, our findings reveal a novel mode of toxin-channel recognition that has important implications for the mechanism of thermosensation.

Prebiotic NH3 Formation: Insights from Simulations

Simulations of prebiotic NH₃ synthesis from NO₃⁻ and NO₂⁻ on pyrite surfaces under hydrothermal conditions are reported. Ab initio metadynamics calculations have successfully explored the full reaction path which explains earlier experimental observations. We have found that the reaction mechanism can be constructed from stepwise single atom transfers which are compatible with the expected reaction time scales. The roles of the hot-pressurized water and of the pyrite surfaces have been addressed. The mechanistic picture that emerged from the simulations strengthens the theory of chemoautotrophic origin of life by providing plausible reaction pathways for the formation of ammonia within the iron-sulfur-world scenario.

Automated Discovery and Refinement of Reactive Molecular Dynamics Pathways

We describe a flexible and broadly applicable energy refinement method, "nebterpolation," for identifying and characterizing the reaction events in a molecular dynamics (MD) simulation. The new method is applicable to ab initio simulations with hundreds of atoms containing complex and multimolecular reaction events. A key aspect of nebterpolation is smoothing of the reactive MD trajectory in internal coordinates to initiate the search for the reaction path on the potential energy surface. We apply nebterpolation to analyze the reaction events in an ab initio nanoreactor simulation that discovers new molecules and mechanisms, including a C-C coupling pathway for glycolaldehyde synthesis. We find that the new method, which incorporates information from the MD trajectory that connects reactants with products, produces a dramatically distinct set of minimum energy paths compared to existing approaches that start from information for the reaction end points alone. The energy refinement method described here represents a key component of an emerging simulation paradigm where molecular dynamics simulations are applied to discover the possible reaction mechanisms.

Predictive Sampling of Rare Conformational Events in Aqueous Solution: Designing a Generalized Orthogonal Space Tempering Method

In aqueous solution, solute conformational transitions are governed by intimate interplays of the fluctuations of solute-solute, solute-water, and water-water interactions. To promote molecular fluctuations to enhance sampling of essential conformational changes, a common strategy is to construct an expanded Hamiltonian through a series of Hamiltonian perturbations and thereby broaden the distribution of certain interactions of focus. Due to a lack of active sampling of configuration response to Hamiltonian transitions, it is challenging for common expanded Hamiltonian methods to robustly explore solvent mediated rare conformational events. The orthogonal space sampling (OSS) scheme, as exemplified by the orthogonal space random walk and orthogonal space tempering methods, provides a general framework for synchronous acceleration of slow configuration responses. To more effectively sample conformational transitions in aqueous solution, in this work, we devised a generalized orthogonal space tempering (gOST) algorithm. Specifically, in the Hamiltonian perturbation part, a solvent-accessible-surface-area-dependent term is introduced to implicitly perturb near-solute water-water fluctuations; more importantly in the orthogonal space response part, the generalized force order parameter is generalized as a two-dimension order parameter set, in which essential solute-solvent and solute-solute components are separately treated. The gOST algorithm is evaluated through a molecular dynamics simulation study on the explicitly solvated deca-alanine (Ala10) peptide. On the basis of a fully automated sampling protocol, the gOST simulation enabled repetitive folding and unfolding of the solvated peptide within a single continuous trajectory and allowed for detailed constructions of Ala10 folding/unfolding free energy surfaces. The gOST result reveals that solvent cooperative fluctuations play a pivotal role in Ala10 folding/unfolding transitions. In addition, our assessment analysis suggests that because essential conformational events are mainly driven by the compensating fluctuations of essential solute-solvent and solute-solute interactions, commonly employed "predictive" sampling methods are unlikely to be effective on this seemingly "simple" system. The gOST development presented in this paper illustrates how to employ the OSS scheme for physics-based sampling method designs.

Simple, yet powerful methodologies for conformational sampling of proteins

Several biological functions, such as molecular recognition, enzyme catalysis, signal transduction, allosteric regulation, and protein folding, are strongly related to conformational transitions of proteins. These conformational transitions are generally induced as slow dynamics upon collective motions, including biologically relevant large-amplitude fluctuations of proteins. Although molecular dynamics (MD) simulation has become a powerful tool for extracting conformational transitions of proteins, it might still be difficult to reach time scales of the biological functions because the accessible time scales of MD simulations are far from biological time scales, even if straightforward conventional MD (CMD) simulations using massively parallel computers are employed. Thus, it is desirable to develop efficient methods to achieve canonical ensembles with low computational costs. From this perspective, we review several enhanced conformational sampling techniques of biomolecules developed by us. In our methods, multiple independent short-time MD simulations are employed instead of single straightforward long-time CMD simulations. Our basic strategy is as follows: (i) selection of initial seeds (initial structures) for the conformational sampling in restarting MD simulations. Here, the seeds should be selected as candidates with high potential to transit. (ii) Resampling from the selected seeds by initializing velocities in restarting short-time MD simulations. A cycle of these simple protocols might drastically promote the conformational transitions of biomolecules. (iii) Once reactive trajectories extracted from the cycles of short-time MD simulations are obtained, a free energy profile is evaluated by means of umbrella sampling (US) techniques with the weighted histogram analysis method (WHAM) as a post-processing technique. For the selection of the initial seeds, we proposed four different choices: (1) Parallel CaScade molecular dynamics (PaCS-MD), (2) Fluctuation Flooding Method (FFM), (3) Outlier FLOODing (OFLOOD) method, and (4) TaBoo SeArch (TBSA) method. We demonstrate applications of our methods to several biological systems, such as domain motions of proteins with large-amplitude fluctuations, conformational transitions upon ligand binding, and protein folding/refolding to native structures of proteins. Finally, we show the conformational sampling efficiencies of our methods compared with those by CMD simulations and other previously developed enhanced conformational sampling methods.

Isolation of pristine MXene from Nb4AlC3 MAX phase: a first-principles study

Complete chemical transformation of MAX (Nb₄AlC₃) into pristine MXene (Nb₄C₃) in the presence of LiF.

A delicate case of unidirectional proton transfer from water to an aromatic heterocyclic anion

One of the hydroxyl modes of the surrounding water molecules with the lowest stretching frequency facilitates the proton transfer from water to an anion.

Deprotonation of formic acid in collisions with a liquid water surface studied by molecular dynamics and metadynamics simulations

Formic acid has a lower barrier to deprotonation at the air–water interface than in bulk liquid water.

Coordination properties of a metal chelator clioquinol to Zn(2+) studied by static DFT and ab initio molecular dynamics

Several lines of evidence supporting the role of metal ions in amyloid aggregation, one of the hallmarks of Alzheimer's disease (AD), have turned metal ion chelation into a promising therapeutic treatment. The design of efficient chelating ligands requires proper knowledge of the electronic and molecular structure of the complexes formed, including their hydration properties. Among various potential chelators, clioquinol (5-chloro-7-iodo-8-hydroxyquinoline, CQH) has been evaluated with relative success in in vitro experiments and even in phase 2 clinical trials. Clioquinol interacts with Zn(ii) to lead to a binary metal/ligand 1 : 2 stoichiometric complex in which the phenolic group of CQH is deprotonated, resulting in Zn(CQ)2 neutral complexes, to which additional water molecules may coordinate. In the present work, the coordinative properties of clioquinol in aqueous solution have been analyzed by means of static, minimal cluster based DFT calculations and explicit solvent ab initio molecular dynamics simulations. Results from static calculations accounting for solvent effects by means of polarized continuum models suggest that the preferred metal coordination environment is tetrahedral Zn(CQ)2, whereas ab initio molecular dynamics simulations point to quasi degenerate penta Zn(CQ)2(H2O) and hexa Zn(CQ)2(H2O)2 coordinated complexes. The possible reasons for these discrepant results are discussed.

Cationic Methylene-Pyrene Isomers and Isomerization Pathways: Finite Temperature Theoretical Studies

This paper provides spectral characterizations of the two isomers of the 1-methylenepyrene cation, namely, the 1-pyrenemethylium and a pyrene-like isomer owing a tropylium cycle. Both are possible photodissociation products of the 1-methylpyrene cation and were proposed as potential contributors to the diffuse interstellar bands. In that respect, vibrational and electronic spectra are computed for the optimized structures at the density functional theory (DFT) and time-dependent (TD-)DFT levels. Finite temperature effects on these spectra are estimated from molecular dynamics simulations within the density functional-based tight-binding (DFTB) and TD-DFTB frameworks, these methods being first benchmarked against DFT and TD-DFT calculations. The computed spectra allow discrimination of the two isomers. When the temperature increases, bands are observed to redshift and merge. The isomerization mechanism is investigated with the metadynamics technique, a biased dynamics scheme allowing to probe reaction mechanisms with high energy barriers by investigating the free energy surface at various temperatures. Four pathways with similar barrier heights (3.5-4 eV) are found, showing that the interconversion process would only occur in interstellar clouds under photoactivation. The present study opens the way to simulations on larger methyl- and methylenePAHs of astrophysical interest and their experimental investigation.

Free-Energy Barriers and Reaction Mechanisms for the Electrochemical Reduction of CO on the Cu(100) Surface, Including Multiple Layers of Explicit Solvent at pH 0

The great interest in the photochemical reduction from CO2 to fuels and chemicals has focused attention on Cu because of its unique ability to catalyze formation of carbon-containing fuels and chemicals. A particular goal is to learn how to modify the Cu catalysts to enhance the production selectivity while reducing the energy requirements (overpotential). To enable such developments, we report here the free-energy reaction barriers and mechanistic pathways on the Cu(100) surface, which produces only CH4 (not C2H4 or CH3OH) in acid (pH 0). We predict a threshold potential for CH4 formation of -0.52 V, which compares well to experiments at low pH, -0.45 to -0.50 V. These quantum molecular dynamics simulations included ∼5 layers of explicit water at the water/electrode interface using enhanced sampling methodology to obtain the free energies. We find that that chemisorbed hydroxyl-methylene (CH-OH) is the key intermediate determining the selectivity for methane over methanol.

Biased Molecular Simulations for Free-Energy Mapping: A Comparison on the KcsA Channel as a Test Case

The calculation of free-energy landscapes in proteins is a challenge for modern numerical simulations. As to the case of potassium ion channels is concerned, it is particularly interesting because of the nanometric dimensions of the selectivity filter, where the complex electrostatics is highly relevant. The present study aims at comparing three different techniques used to bias molecular dynamics simulations, namely Umbrella Sampling, Steered Molecular Dynamics, and Metadynamics, never applied all together in the past to the same channel protein. Our test case is represented by potassium ions permeating the selectivity filter of the KcsA channel.

Free Energy Barriers for the N-Terminal Asparagine to Succinimide Conversion: Quantum Molecular Dynamics Simulations for the Fully Solvated Model

Deamidation of asparagine residues represents one of the main routes for the post-translational modification of protein sequences. We computed the estimates of the free energy barriers for three stages of the deamidation process, deprotonation, cyclization, and deamination, of the conversion of asparagine to the succinimide intermediate within the fully solvated model with explicit water molecules. The Born-Oppenheimer molecular dynamics in the Gaussian and Plane Wave (GPW) approximation as implemented in the CP2K quantum chemistry package was utilized to sample the configurational space. By applying the metadynamics technique, the estimates of the free energy barriers were obtained for three separated stages of the reaction. In agreement with the experimental kinetic measurements, the estimated activation barriers do not exceed 21 kcal/mol. We demonstrate that the use of fully solvated models is the critical issue in theoretical studies of these reactions. We also conclude that more extensive sampling is necessary to obtain full free energy profiles and accurate barriers for the reaction stages.

Ab Initio MD Simulations of the Brønsted Acidity of Glutathione in Aqueous Solutions: Predicting pKa Shifts of the Cysteine Residue

The tripeptide glutathione (GSH) is one of the most abundant peptides and the major repository for nonprotein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thiol-disulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influence the redox couple, and hence the pKa value of the cysteine residue of GSH is critical to its functioning. Here we report ab initio Car-Parrinello molecular dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. We show that the free-energy landscape for the protonation-deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides accurate estimates of the pKa and correctly predicts the shift in the dissociation constant values as compared with the isolated cysteine amino acid.

Computational metallomics of the anticancer drug cisplatin

Cisplatin, cis-diamminedichlorido-platinum(II), is an important therapeutic tool in the struggle against different tumors, yet it is plagued with the emergence of resistance mechanisms after repeated administrations. This hampers greatly its efficacy. Overcoming resistance problems requires first and foremost an integrated and systematic understanding of the structural determinants and molecular recognition processes involving the drug and its cellular targets. Here we review a strategy that we have followed for the last few years, based on the combination of modern tools from computational chemistry with experimental biophysical methods. Using hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) simulations, validated by spectroscopic experiments (including NMR, and CD), we have worked out for the first time at atomic level the structural determinants in solution of platinated cellular substrates. These include the copper homeostasis proteins Ctr1, Atox1, and ATP7A. All of these proteins have been suggested to influence the pre-target resistance mechanisms. Furthermore, coupling hybrid QM/MM simulations with classical Molecular Dynamics (MD) and free energy calculations, based on force field parameters refined by the so-called "Force Matching" procedure, we have characterized the structural modifications and the free energy landscape associated with the recognition between platinated DNA and the protein HMGB1, belonging to the chromosomal high-mobility group proteins HMGB that inhibit the repair of platinated DNA. This may alleviate issues relative to on-target resistance process. The elucidation of the mechanisms by which tumors are sensitive or refractory to cisplatin may lead to the discovery of prognostic biomarkers. The approach reviewed here could be straightforwardly extended to other metal-based drugs.