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

Self-Uptake Mechanism of Polymyxin-Based Lipopeptide against Gram-Negative Bacterial Membrane: Role of the First Adsorbed Lipopeptide

Research in the continuously increasing threat of polymyxin-resistant multidrug-resistant Pseudomonas aeruginosa, which causes severe infection in immunocompromised patients, has resulted in the development of several polymyxin-derived cyclic lipopeptides containing l-α-γ- diamino butyric acid-like FADDI-019 (F19). In this work, F19's insertion into a minimal model of the asymmetric outer membrane of the bacterium, which contained only penta-acylated lipid A (LipA) and lacked keto-d-octulosonic acid and O-antigens, in the top leaflet and phospholipids in the bottom leaflet, was studied. F19 exhibited all of the hallmarks of the self-uptake mechanism into the asymmetric bilayer. While a single monomer of the lipopeptide did not get partitioned into the inside of the bilayer, it competitively displaced Ca²⁺ from the membrane surface, observed as a decrease in Ca²⁺ coordination number with phosphate groups (1.89 vs 1.718), resulting in membrane destabilization. This resulted in an increment of the average defect size and the probability of interplay between lipid tails and hydrophobic residues of another F19. When more than one monomer was present in the system, the first monomer remained docked on the surface, while other monomers intercalated into the bilayer interior with their hydrophobic moieties "sleeved" by lipid acyl chains. The free energy barrier for partial insertion of the lipopeptide into a bilayer in the presence of surface-docked second F19 was recorded at ∼1.3 kcal/mol using two-dimensional (2D) well-tempered metadynamics, making it a low barrier process at 300 K. This study is an attempt to demonstrate the self-uptake mechanism of F19 during intercalation process into the bilayer interior, which may help in the design of better alternates for polymyxins to work against polymyxin resistance.

Reversible assembly of nanoparticles: theory, strategies and computational simulations

The significant advances in synthesis and functionalization have enabled the preparation of high-quality nanoparticles that have found a plethora of successful applications. The unique physicochemical properties of nanoparticles can be manipulated through the control of size, shape, composition, and surface chemistry, but their technological application possibilities can be further expanded by exploiting the properties that emerge from their assembly. The ability to control the assembly of nanoparticles not only is required for many real technological applications, but allows the combination of the intrinsic properties of nanoparticles and opens the way to the exploitation of their complex interplay, giving access to collective properties. Significant advances and knowledge gained over the past few decades on nanoparticle assembly have made it possible to implement a growing number of strategies for reversible assembly of nanoparticles. In addition to being of interest for basic studies, such advances further broaden the range of applications and the possibility of developing innovative devices using nanoparticles. This review focuses on the reversible assembly of nanoparticles and includes the theoretical aspects related to the concept of reversibility, an up-to-date assessment of the experimental approaches applied to this field and the advanced computational schemes that offer key insights into the assembly mechanisms. We aim to provide readers with a comprehensive guide to address the challenges in assembling reversible nanoparticles and promote their applications.

Vibrational Circular Dichroism from DFT Molecular Dynamics: The AWV Method

The paper illustrates the Activity Weighted Velocities (AWV) methodology to compute Vibrational Circular Dichroism (VCD) anharmonic spectra from Density Functional Theory (DFT) molecular dynamics. AWV calculates the spectra by the Fourier Transform of the time correlation functions of velocities, weighted by specific observables: the Atomic Polar Tensors (APTs) and the Atomic Axial Tensors (AATs). Indeed, AWV shows to correctly reproduce the experimental spectra for systems in the gas and liquid phases, both in the case of weakly and strongly interacting systems. The comparison with the experimental spectra is striking especially in the fingerprint region, as demonstrated by the three benchmark systems discussed: (1S)-Fenchone in the gas phase, (S)-(-)-Propylene oxide in the liquid phase, and (R)-(-)-2-butanol in the liquid phase. The time evolution of APTs and AATs can be adequately described by a linear combination of the tensors of a small set of appropriate reference structures, strongly reducing the computational cost without compromising accuracy. Additionally, AWV allows the partition of the spectral signal in its molecular components without any expensive postprocessing and any localization of the charge density or the wave function.

O-Acetyl Migration within the Sialic Acid Side Chain: A Mechanistic Study Using the Ab Initio Nanoreactor

Many disease-causing viruses target sialic acids on the surface of host cells. Some viruses bind preferentially to sialic acids with O-acetyl modification at the hydroxyl group of C7, C8, or C9 on the glycerol-like side chain. Studies of proteins binding to sialosides containing O-acetylated sialic acids are crucial in understanding the related diseases but experimentally difficult due to the lability of the ester group. We recently showed that O-acetyl migration among hydroxyl groups of C7, C8, and C9 in sialic acids occurs in all directions in a pH-dependent manner. In the current study, we elucidate a full mechanistic pathway for the migration of O-acetyl among C7, C8, and C9. We used an ab initio nanoreactor to explore potential reaction pathways and density functional theory, pK_a calculations, and umbrella sampling to investigate elementary steps of interest. We found that when a base is present, migration is easy in any direction and involves three key steps: deprotonation of the hydroxyl group, cyclization between the two carbons, and the migration of the O-acetyl group. This dynamic equilibrium may play a defensive role against pathogens that evolve to gain entry to the cell by binding selectively to one acetylation state.

Free Energy Landscape of the Adenylation Reaction of the Aminoacylation Process at the Active Site of Aspartyl tRNA Synthetase

The process of protein biosynthesis is initiated by the aminoacylation process where a transfer ribonucleic acid (tRNA) is charged by the attachment of its cognate amino acid at the active site of the corresponding aminoacyl tRNA synthetase enzyme. The first step of the aminoacylation process, known as the adenylation reaction, involves activation of the cognate amino acid where it reacts with a molecule of adenosine triphosphate (ATP) at the active site of the enzyme to form the aminoacyl adenylate and inorganic pyrophosphate. In the current work, we have investigated the adenylation reaction between aspartic acid and ATP at the active site of the fully solvated aspartyl tRNA synthetase (AspRS) from Escherichia coli in aqueous medium at room temperature through hybrid quantum mechanical/molecular mechanical (QM/MM) simulations combined with enhanced sampling methods of well-tempered and well-sliced metadynamics. The objective of the present work is to study the associated free energy landscape and reaction barrier and also to explore the effects of active site mutation on the free energy surface of the reaction. The current calculations include finite temperature effects on free energy profiles. In particular, apart from contributions of interaction energies, the current calculations also include contributions of conformational, vibrational, and translational entropy of active site residues, substrates, and also the rest of the solvated protein and surrounding water into the free energy calculations. The present QM/MM metadynamics simulations predict a free energy barrier of 23.35 and 23.5 kcal/mol for two different metadynamics methods used to perform the reaction at the active site of the wild type enzyme. The free energy barrier increases to 30.6 kcal/mol when Arg217, which is an important conserved residue of the wild type enzyme at its active site, is mutated by alanine. These free energy results including the effect of mutation compare reasonably well with those of kinetic experiments that are available in the literature. The current work also provides molecular details of structural changes of the reactants and surroundings as the system dynamically evolves along the reaction pathway from reactant to the product state through QM/MM metadynamics simulations.

The chemistry of AlF and CaF production in buffer gas sources

In this work, we explore the role of chemical reactions on the properties of buffer gas cooled molecular beams. In particular, we focus on scenarios relevant to the formation of AlF and CaF via chemical reactions between the Ca and Al atoms ablated from a solid target in an atmosphere of a fluorine-containing gas, in this case, SF6 and NF3. Reactions are studied following an ab initio molecular dynamics approach, and the results are rationalized following a tree-shaped reaction model based on Bayesian inference. We find that NF3 reacts more efficiently with hot metal atoms to form monofluoride molecules than SF6. In addition, when using NF3, the reaction products have lower kinetic energy, requiring fewer collisions to thermalize with the cryogenic helium. Furthermore, we find that the reaction probability for AlF formation is much higher than for CaF across a broad range of temperatures.

Exploration of high dimensional free energy landscapes by a combination of temperature-accelerated sliced sampling and parallel biasing

Temperature-accelerated sliced sampling (TASS) is an enhanced sampling method for achieving accelerated and controlled exploration of high-dimensional free energy landscapes in molecular dynamics simulations. With the aid of umbrella bias potentials, the TASS method realizes a controlled exploration and divide-and-conquer strategy for computing high-dimensional free energy surfaces. In TASS, diffusion of the system in the collective variable (CV) space is enhanced with the help of metadynamics bias and elevated-temperature of the auxiliary degrees of freedom (DOF) that are coupled to the CVs. Usually, a low-dimensional metadynamics bias is applied in TASS. In order to further improve the performance of TASS, we propose here to use a highdimensional metadynamics bias, in the same form as in a parallel bias metadynamics scheme. Here, a modified reweighting scheme, in combination with artificial neural network is used for computing unbiased probability distribution of CVs and projections of high-dimensional free energy surfaces. We first validate the accuracy and efficiency of our method in computing the four-dimensional free energy landscape for alanine tripeptide in vacuo. Subsequently, we employ the approach to calculate the eight-dimensional free energy landscape of alanine pentapeptide in vacuo. Finally, the method is applied to a more realistic problem wherein we compute the broad four-dimensional free energy surface corresponding to the deacylation of a drug molecule which is covalently complexed with a β-lactamase enzyme. We demonstrate that using parallel bias in TASS improves the efficiency of exploration of high-dimensional free energy landscapes.

Deciphering Ethanol-Driven Swelling, Rupturing, Aggregation, and Fusion of Lipid Vesicles Using Coarse-Grained Molecular Dynamics Simulations

Traditionally, liquid ethanol is known to enhance the permeability of lipid membranes and causes vesicle aggregation and fusion. However, how the amphiphilic ethanol molecules perturb the lipid vesicles to facilitate their aggregation or fusion has not been addressed at any level of molecular simulations. Herein, not only have we developed a coarse-grained (CG) model for liquid ethanol, its aqueous mixture, and hydrated lipid membranes for molecular dynamics (MD) simulations, but also utilized it to delineate the aggregation and fusion of lipid vesicles using CG-MD simulations with multimillion particles. We have systematically parametrized the force-field for pure ethanol and its interactions with hydrated POPC and POPE model lipid membranes. In this process, we have successfully reproduced the bulk ethanol structure and concentration-dependent density of aqueous ethanol. To quantify the interaction of ethanol with lipid membranes, we have reproduced the transfer free energy of the ethanol molecule across the hydrated bilayers, and the concentration-dependent distribution of ethanol molecules across the lipid bilayers. After having acceptable force-field parameters for ethanol-membrane interactions, we have checked the effect of ethanol toward the vesicles comprising POPC lipids. We observe a rapid increase in the size of the POPC lipid vesicles with increasing amounts of ethanol up to 30 mol %. We unambiguously observe swelling and decrease in the thickness of the POPC vesicles with increasing amounts of ethanol up to 30 mol %, beyond which the vesicles begin to lose their integrity and rupture at higher mol % of ethanol. The fusion study of two vesicles demonstrates that fused vesicles can be obtained from 20 to 30 mol % of ethanol provided that they are brought closer than a critical distance at a particular mol %. The multivesicle simulations show that along with the increase in the sizes of vesicles the propensity of vesicle aggregation increases as the mol % of ethanol increases.

Understanding High-Temperature Chemical Reactions on Metal Surfaces: A Case Study on Equilibrium Concentration and Diffusivity of C x H y on a Cu(111) Surface

Chemical reactions on metal surfaces are important in various processes such as heterogeneous catalysis and nanostructure growth. At moderate or lower temperatures, these reactions generally follow the minimum energy path, and temperature effects can be reasonably described by a harmonic oscillator model. At a high temperature approaching the melting point of the substrate, general behaviors of surface reactions remain elusive. In this study, by taking hydrocarbon species adsorbed on Cu(111) as a model system and performing extensive molecular dynamics simulations powered by machine learning potentials, we identify several important high-temperature effects, including local chemical environment, surface atom mobility, and substrate thermal expansion. They affect different aspects of a high-temperature surface reaction in different ways. These results deepen our understanding of high-temperature reactions.

Site dependent catalytic water dissociation on an anisotropic buckled black phosphorus surface

Black phosphorus (BP) is unique among 2D materials due to its anisotropic puckered structure. It has been used as a multifunctional catalyst for various purposes. In this study, we performed first principles molecular dynamics simulations to understand the water-splitting reaction on a bi-layer BP surface. We focused on the site-specific aqueous reactivity of the buckled surface. A difference in the axis-dependent reactivity is observed owing to edge defects and exposed sites. Thus, we believe that BP edges, which significantly affect the interfacial water or organic solvent molecules, must exhibit very different edge-dependent reactivity. Experiments suggested the increasing catalytic efficiency of undisturbed BP in the order bulk, few-layered BP, and BP quantum dots. We choose three active sites to investigate the mechanistic details of the OER: the zigzag (ZZ), armchair (AC), and bulk sites. This study will provide insight into the enhanced catalytic activity when more edges are exposed as the active surface. We hope to clarify the reactive pathway in an aqueous solution supported by bi-layer BP by exploring the two different mechanisms for forming the OOH* complex. We explore and report two mechanisms: a simple push-pull reaction for oxygen-oxygen bond formation, the nucleophilic attack by formed OH^- and an attack by a water molecule. The free energy barriers procured for mechanism 1 taking place at the zigzag, armchair, and bulk sites are 7.59 ± 0.33, 9.04 ± 0.01, and 12.80 ± 0.09 kcal mol^-1, respectively. For mechanism 2 the free energy barriers are 7.62 ± 0.11, 9.15 ± 0.16, and 11.63 ± 0.11 kcal mol^-1 for the ZZ, AC, and bulk sites. The interlink between both the mechanisms is established concerning the reported free energy barriers for OOH* formation. The ZZ site is found to lower the activation barrier for the rate-determining step, followed by the AC and bulk.

High-pressure phases of a Mn-N system

Highly compressed extended states of light elemental solids have emerged recently as a novel group of energetic materials. The application of these materials is seriously limited by the energy-safety contradiction, because the material with high energy density is highly metastable and can hardly be recovered under ambient conditions. Recently, it has been found that high-energy density transition metal polynitrides could be synthesized at ∼100 GPa and recovered at ∼20 GPa. Inspired by these findings, we have studied a high-pressure Mn-N system from the aspects of structure, stability, phase transition, energy density and electronic structure theoretically for the first time. The results reveal that MnN₄_P1̄ consisting of [N₄]_∞^2- is thermodynamically stable at 36.9-100 GPa, dynamically stable at 0 GPa and has a noticeably high volumetric energy density of 15.71 kJ cm^-3. Upon decompression, this structure will transform to MnN₄_C2/m with the transition barrier declining sharply at 5-10 GPa due to the switching of transition pathways. Hence, we propose MnN₄_P1̄ as a potential energetic material that is synthesizable above 40 GPa and recoverable until 10 GPa.

Catalytic properties of the ferryl ion in the solid state: a computational review

This review summarises the last findings in the emerging field of heterogeneous catalytic oxidation of light alkanes by ferryl species supported on solid-state systems such as the conversion of methane into methanol by FeO-MOF74.

Discovering surface reaction pathways using accelerated molecular dynamics and network analysis tools

We present an automated method that maps surface reaction pathways with no experimental data and with minimal human interventions. In this method, bias potentials promoting surface reactions are applied to enable statistical samplings of the surface reaction within the timescale of ab initio molecular dynamics (MD) simulations, and elementary reactions are extracted automatically using an extension of the method constructed for gas- or liquid-phase reactions. By converting the extracted elementary reaction data to directed graph data, MD trajectories can be efficiently mapped onto reaction pathways using a network analysis tool. To demonstrate the power of the method, it was applied to the steam reforming of methane on the Rh(111) surface and to propane reforming on the Pt(111) and Pt₃Sn(111) surfaces. We discover new energetically favorable pathways for both reactions and reproduce the experimentally-observed materials-dependence of the surface reaction activity and the selectivity for the propane reforming reactions.

We present an automated method that maps surface reaction pathways with no experimental data and with minimal human interventions.

Mechanistic Insight on BioIL-Induced Structural Alterations in DMPC Lipid Bilayer

The emerging application risks of traditional ionic liquids (ILs) toward the ecosystem have changed the perception regarding their greenness. This resulted in the exploration of their more biocompatible alternatives known as biocompatible ILs (BioILs). Here, we have investigated the impact of two such biocompatible cholinium amino acid-based ILs on the structural behavior of model homogeneous DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) lipid bilayer using all-atom molecular dynamics simulation technique. Two classic cholinium-amino acid-based ILs, cholinium glycinate ([Ch][Gly]) and cholinium phenylalaninate ([Ch][Phe]), which differ only by the side chain lengths and hydrophobicity of the anions, have been utilized in the present work. Simultaneous analysis of the bilayer structural properties reveals that the existence of [Ch][Gly] BioIL above a particular concentration induces phase transition from fluid phase to gel phase in the DMPC lipid bilayer. Such a freezing of lipid bilayer upon the exposure to concentrated aqueous solution of [Ch][Gly] BioIL indicates the harmfulness of this BioIL toward the cell membranes majorly containing DMPC lipids, as the cell freezing can negatively affect its stability and functionality. Despite having a more hydrophobic amino acid side chain of [Phe]^- anion in [Ch][Phe], in the case of bilayer-[Ch][Phe] systems we observe the minimal impact of [Ch][Phe] BioIL on the DMPC bilayer properties up to 10 mol % concentration. In the presence of these BioIL, we observe the thickening of the bilayer and accumulation of the cations and anions of the BioILs at the interface of DMPC lipid heads and tails. The transfer free-energy profile of a [Phe]^- anion from aqueous phase to membrane center also indicates the anion partitioning at lipid head-tail interface and its inability to penetrate in the lipid membrane tail region. In contrast, the free-energy profile for a [Gly]^- anion offers a very high energy barrier to the insertion of [Gly]^- into the membrane interior, leading to accumulation of [Gly]^- anions at the lipid head-water region.

Free Energy Landscape and Proton Transfer Pathways of the Transimination Reaction at the Active site of the Serine Hydroxymethyltransferase Enzyme in Aqueous Medium

Serine hydroxymethyltransferase (SHMT) is a ubiquitous enzyme belonging to the fold type I or aspartate aminotransferase (AspAT) family of the pyridoxal 5'-phosphate (PLP)-dependent enzymes. Like other PLP-dependent enzymes, SHMT also undergoes the so-called transimination reaction before exhibiting its enzymatic activity. The transimination process constitutes an important pre-step for all PLP-dependent enzymes, where an internal aldimine of the PLP-enzyme complex gets converted to an external aldimine of the substrate-PLP complex at the active site of the enzyme. In case of the transimination reaction involving SHMT, the PLP molecule bound to the active site lysine residue of SHMT (internal aldimine) gets detached from the enzyme by a serine substrate to produce an external aldimine complex, where the PLP is now bound to the serine substrate. In the current study, the free energy surfaces and reaction pathways of different steps of the transimination reaction at the active site of SHMT are investigated by employing hybrid quantum mechanical/molecular mechanical (QM/MM) simulations combined with metadynamics methods of rare event sampling. It is found that the process of transimination involving serine and PLP at the active site of the SHMT enzyme takes place through different elementary steps such as the formation of the first geminal diamine intermediate (GDI1), transfer of a proton from the substrate serine to the phenolic oxygen of PLP, followed by another proton transfer from PLP to the amine nitrogen of lysine with the formation of the second geminal diamine intermediate (GDI2), and finally, detachment of the active site lysine residue from PLP to produce the external aldimine.

Collective variable-based enhanced sampling and machine learning

ABSTRACT

Collective variable-based enhanced sampling methods have been widely used to study thermodynamic properties of complex systems. Efficiency and accuracy of these enhanced sampling methods are affected by two factors: constructing appropriate collective variables for enhanced sampling and generating accurate free energy surfaces. Recently, many machine learning techniques have been developed to improve the quality of collective variables and the accuracy of free energy surfaces. Although machine learning has achieved great successes in improving enhanced sampling methods, there are still many challenges and open questions. In this perspective, we shall review recent developments on integrating machine learning techniques and collective variable-based enhanced sampling approaches. We also discuss challenges and future research directions including generating kinetic information, exploring high-dimensional free energy surfaces, and efficiently sampling all-atom configurations.

GRAPHIC ABSTRACT

Nanoconfinement and mass transport in metal-organic frameworks

The ubiquity of metal-organic frameworks in recent scientific literature underscores their highly versatile nature. MOFs have been developed for use in a wide array of applications, including: sensors, catalysis, separations, drug delivery, and electrochemical processes. Often overlooked in the discussion of MOF-based materials is the mass transport of guest molecules within the pores and channels. Given the wide distribution of pore sizes, linker functionalization, and crystal sizes, molecular diffusion within MOFs can be highly dependent on the MOF-guest system. In this review, we discuss the major factors that govern the mass transport of molecules through MOFs at both the intracrystalline and intercrystalline scale; provide an overview of the experimental and computational methods used to measure guest diffusivity within MOFs; and highlight the relevance of mass transfer in the applications of MOFs in electrochemical systems, separations, and heterogeneous catalysis.

Dynamics and Site Isolation: Keys to High Propane Dehydrogenation Performance of Silica-Supported PtGa Nanoparticles

Nonoxidative dehydrogenation of light alkanes has seen a renewed interest in recent years. While PtGa systems appear among the most efficient catalyst for this reaction and are now implemented in production plants, the origin of the high catalytic performance in terms of activity, selectivity, and stability in PtGa-based catalysts is largely unknown. Here we use molecular modeling at the DFT level on three different models: (i) periodic surfaces, (ii) clusters using static calculations, and (iii) realistic size silica-supported nanoparticles (1 nm) using molecular dynamics and metadynamics. The combination of the models with experimental data (XAS, TEM) allowed the refinement of the structure of silica-supported PtGa nanoparticles synthesized via surface organometallic chemistry and provided a structure-activity relationship at the molecular level. Using this approach, the key interaction between Pt and Ga was evidenced and analyzed: the presence of Ga increases (i) the interaction between the oxide surface and the nanoparticles, which reduces sintering, (ii) the Pt site isolation, and (iii) the mobility of surface atoms which promotes the high activity, selectivity, and stability of this catalyst. Considering the complete system for modeling that includes the silica support as well as the dynamics of the PtGa nanoparticle is essential to understand the catalytic performances.

Molecular Dynamics and Machine Learning in Catalysts

Given the importance of catalysts in the chemical industry, they have been extensively investigated by experimental and numerical methods. With the development of computational algorithms and computer hardware, large-scale simulations have enabled influential studies with more atomic details reflecting microscopic mechanisms. This review provides a comprehensive summary of recent developments in molecular dynamics, including ab initio molecular dynamics and reaction force-field molecular dynamics. Recent research on both approaches to catalyst calculations is reviewed, including growth, dehydrogenation, hydrogenation, oxidation reactions, bias, and recombination of carbon materials that can guide catalyst calculations. Machine learning has attracted increasing interest in recent years, and its combination with the field of catalysts has inspired promising development approaches. Its applications in machine learning potential, catalyst design, performance prediction, structure optimization, and classification have been summarized in detail. This review hopes to shed light and perspective on ML approaches in catalysts.

Chemically Accurate Vibrational Free Energies of Adsorption from Density Functional Theory Molecular Dynamics: Alkanes in Zeolites

We present a methodology to compute, at reduced computational cost, Gibbs free energies, enthalpies, and entropies of adsorption from molecular dynamics. We calculate vibrational partition functions from vibrational energies, which we obtain from the vibrational density of states by projection on the normal modes. The use of a set of well-chosen reference structures along the trajectories accounts for the anharmonicities of the modes. For the adsorption of methane, ethane, and propane in the H-CHA zeolite, we limit our treatment to a set of vibrational modes localized at the adsorption site (zeolitic OH group) and the alkane molecule interacting with it. Only two short trajectories (1–20 ps) are required to reach convergence (<1 kJ/mol) for the thermodynamic functions. The mean absolute deviations from the experimentally measured values are 2.6, 2.8, and 4.7 kJ/mol for the Gibbs free energy, the enthalpy, and the entropy term (−TΔS), respectively. In particular, the entropy terms show a major improvement compared to the harmonic approximation and almost reach the accuracy of the previous use of anharmonic frequencies obtained with curvilinear distortions of individual modes. The thermodynamic functions so obtained follow the trend of the experimental values for methane, ethane, and propane, and the Gibbs free energy of adsorption at experimental conditions is correctly predicted to change from positive for methane (5.9 kJ/mol) to negative for ethane (−4.8 kJ/mol) and propane (−7.1 kJ/mol).

Mechanistic Insight into the O2 Evolution Catalyzed by Copper Complexes with Tetra- and Pentadentate Ligands

The mononuclear complexes ([(bztpen)Cu] (BF₄)₂ (bztpen = N-benzyl-N,N',N'-tris (pyridin-2-yl methyl ethylenediamine))) and ([(dbzbpen)Cu(OH₂)] (BF₄)₂ (dbzbpen = N,N'-dibenzyl-N,N'-bis(pyridin-2-ylmethyl) ethylenediamine)) have been reported as water oxidation catalysts in basic medium (pH = 11.5). We explore the O₂ evolution process catalyzed by these copper catalysts with various ligands (L) by applying the first-principles molecular dynamics simulations. First, the oxidation of catalysts to the metal-oxo intermediates [LCu(O)]²⁺ occurs through the proton-coupled electron transfer (PCET) process. These intermediates are involved in the oxygen-oxygen bond formation through the water-nucleophilic addition process. Here, we have considered two types of oxygen-oxygen bond formation. The first one is the transfer of the hydroxide of the water molecule to the Cu═O moiety; the proton transfer to the solvent leads to the formation of the peroxide complex ([LCu(OOH)]⁺). The other is the formation of the hydrogen peroxide complex ([LCu(HOOH)]²⁺) by the transfer of proton and hydroxide of the water molecule to the metal-oxo intermediate. The formation of the peroxide complex requires less activation free energy than hydrogen peroxide formation for both catalysts. We found two transition states in the well-tempered metadynamics simulations: one for proton transfer and another for hydroxide transfer. In both cases, the proton transfer requires higher free energy. Following the formation of the oxygen-oxygen bond, we study the release of the dioxygen molecule. The formed peroxide and hydrogen peroxide complexes are converted into the superoxide complex ([LCu(OO)]²⁺) through the transfer of proton, electron, and PCET processes. The superoxide complex releases an oxygen molecule upon the addition of a water molecule. The free energy of activation for the release of the dioxygen molecule is lesser than that of the oxygen-oxygen bond formation. When we observe the entire water oxidation process, the oxygen-oxygen bond formation is the rate-determining step. We calculated the rates of reaction by using the Eyring equation and found them to be close to the experimental values.

Overcoming Free-Energy Barriers with a Seamless Combination of a Biasing Force and a Collective Variable-Independent Boost Potential

Amid collective-variable (CV)-based importance-sampling algorithms, a hybrid of the extended adaptive biasing force and the well-tempered metadynamics algorithms (WTM-eABF) has proven particularly cost-effective for exploring the rugged free-energy landscapes that underlie biological processes. However, as an inherently CV-based algorithm, this hybrid scheme does not explicitly accelerate sampling in the space orthogonal to the chosen CVs, thereby limiting its efficiency and accuracy, most notably in those cases where the slow degrees of freedom of the process at hand are not accounted for in the model transition coordinate. Here, inspired by Gaussian-accelerated molecular dynamics (GaMD), we introduce the same CV-independent harmonic boost potential into WTM-eABF, yielding a hybrid algorithm coined GaWTM-eABF. This algorithm leans on WTM-eABF to explore the transition coordinate with a GaMD-mollified potential and recovers the unbiased free-energy landscape through thermodynamic integration followed by proper reweighting. As illustrated in our numerical tests, GaWTM-eABF effectively overcomes the free-energy barriers in orthogonal space and correctly recovers the unbiased potential of mean force (PMF). Furthermore, applying both GaWTM-eABF and WTM-eABF to two biologically relevant processes, namely, the reversible folding of (i) deca-alanine and (ii) chignolin, our results indicate that GaWTM-eABF reduces the uncertainty in the PMF calculation and converges appreciably faster than WTM-eABF. Obviating the need of multiple-copy strategies, GaWTM-eABF is a robust, computationally efficient algorithm to surmount the free-energy barriers in orthogonal space and maps with utmost fidelity the free-energy landscape along selections of CVs. Moreover, our strategy that combines WTM-eABF with GaMD can be easily extended to other biasing-force algorithms.

Mechanism of the Conformational Change of the Protein Methyltransferase SMYD3: A Molecular Dynamics Simulation Study

SMYD3 is a SET-domain-containing methyltransferase that catalyzes the transfer of methyl groups onto lysine residues of substrate proteins. Methylation of MAP3K2 by SMYD3 has been implicated in Ras-driven tumorigenesis, which makes SMYD3 a potential target for cancer therapy. Of all SMYD family proteins, SMYD3 adopt a closed conformation in a crystal structure. Several studies have suggested that the conformational changes between the open and closed forms may regulate the catalytic activity of SMYD3. In this work, we carried out extensive molecular dynamics simulations on a series of complexes with a total of 21 μs sampling to investigate the conformational changes of SMYD3 and unveil the molecular mechanisms. Based on the C-terminal domain movements, the simulated models could be depicted in three different conformational states: the closed, intermediate and open states. Only in the case that both the methyl donor binding pocket and the target lysine-binding channel had bound species did the simulations show SMYD3 maintaining its conformation in the closed state, indicative of a synergetic effect of the cofactors and target lysine on regulating the conformational change of SMYD3. In addition, we performed analyses in terms of structure and energy to shed light on how the two regions might regulate the C-terminal domain movement. This mechanistic study provided insights into the relationship between the conformational change and the methyltransferase activity of SMYD3. The more complete understanding of the conformational dynamics developed here together with further work may lay a foundation for the rational drug design of SMYD3 inhibitors.

Isohemigossypolone: Antiophidic properties of a naphthoquinone isolated from Pachira aquatica Aubl

We investigated the antiophidic properties of isohemigossypolone (ISO), a naphthoquinone isolated from the outer bark of the Pachira aquatic Aubl. The inhibition of phospholipase A₂, coagulant, fibrinogenolytic, hemorrhagic and myotoxic activities induced by Bothrops pauloensis venom (Pb) was investigated. For this, we use samples resulting from the incubation of Pb with ISO in different concentrations (1:1, 1:5 and 1:10 w/w), we also evaluated the condition of treatment using ISO after 15 min of venom inoculation. The activities of phospholipase A₂, coagulant, fibrinogenolytic, hemorrhagic and myotoxic induced by the B. pauloensis venom were significantly inhibited when the ISO was pre-incubated with the crude venom. For in vivo neutralization tests, the results were observed even when the ISO was applied after 15 min of inoculation of the venom or metalloprotease (BthMP). Also, to identify the inhibition mechanism, we performed in silico assays, across simulations of molecular coupling and molecular dynamics, it was possible to identify the modes of interaction between ISO and bothropic toxins BmooMPα-I, Jararacussin-I and BNSP-7. The present study shows that naphthoquinone isohemigossypolone isolated from the P. aquatica plant inhibited part of the local and systemic damage caused by venom proteins, demonstrating the pharmacological potential of this compound in neutralizing the harmful effects caused by snakebites.

Comparative first principles-based molecular dynamics study of catalytic mechanism and reaction energetics of water oxidation reaction on 2D-surface

The study of the water-splitting process, which can proceed in 2e^- as well as 4e^- pathway, reveals that the process is entirely an uphill process, and the third step, that is, the oxooxo bond formation is the rate-determining step. The kinetic barrier of the oxygen evolution reaction (OER) on the 2D material catalysts in the presence of explicit solvents is scarcely studied. Here, we investigate the dynamics of the OER on the undoped graphene and the activation energy barrier of each step using first principles molecular dynamics simulations. Here we provide a detailed analysis of the kinetics of all the 4e^- transfer steps of OER on the graphene surface. We also compare the accuracy of one of the density functional theory (DFT) functionals and density functional based tight binding (DFTB) method in explaining the OER steps. The comparative study reveals that DFTB can be used for performing metadynamics simulations quipped with much less computational cost than DFT functionals. By both Perdew-Burke-Ernzerhof and DFTB methods, the third step is revealed to be the rate-determining step with an energy barrier of 21.19 ± 0.51 and 20.23 ± 0.20 kcal mol^-1 , respectively. DFTB gives an impression of being successful in predicting the energy barriers of OER in 4e- transfer pathway and comparable to the DFT method, and we would like to extend the use of DFTB for further studies with a sizable and complex system.

Unveiling Zwitterionization of Glycine in the Microhydration Limit

Charge separation under solvation stress conditions is a fundamental process that comes in many forms in doped water clusters. Yet, the mechanism of intramolecular charge separation, where constraints due to the molecular structure might be intricately tied to restricted solvation structures, remains largely unexplored. Microhydrated amino acids are such paradigmatic molecules. Ab initio simulations are carried out at 300 K in the frameworks of metadynamics sampling and thermodynamic integration to map the thermal mechanisms of zwitterionization using Gly(H₂O) _n with n = 4 and 10. In both cases, a similar water-mediated proton transfer chain mechanism is observed; yet, detailed analyses of thermodynamics and kinetics demonstrate that the charge-separated zwitterion is the preferred species only for n = 10 mainly due to kinetic stabilization. Structural analyses disclose that bifurcated H-bonded water bridges, connecting the cationic and anionic sites in the fluctuating microhydration network at room temperature, are enhanced in the transition-state ensemble exclusively for n = 10 and become overwhelmingly abundant in the stable zwitterion. The findings offer potential insights into charge separation under solvation stress conditions beyond the present example.

Chemically Induced Mismatch of Rings and Stations in [3]Rotaxanes

The mechanical interlocking of molecular components can lead to the appearance of novel and unconventional properties and processes, with potential relevance for applications in nanoscience, sensing, catalysis, and materials science. We describe a [3]rotaxane in which the number of recognition sites available on the axle component can be changed by acid–base inputs, encompassing cases in which this number is larger, equal to, or smaller than the number of interlocked macrocycles. These species exhibit very different properties and give rise to a unique network of acid–base reactions that leads to a fine pK_a tuning of chemically equivalent acidic sites. The rotaxane where only one station is available for two rings exhibits a rich coconformational dynamics, unveiled by an integrated experimental and computational approach. In this compound, the two crown ethers compete for the sole recognition site, but can also come together to share it, driven by the need to minimize free energy without evident inter-ring interactions.

Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms

We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.

Reversible Hydration of α-Dicarbonyl Compounds from Ab Initio Metadynamics Simulations: Comparison between Pyruvic and Glyoxylic Acids in Aqueous Solutions

Glyoxylic and pyruvic oxoacids are widely available in the atmosphere as gas-phase clusters and particles or in wet aerosols. In aqueous conditions, they undergo interconversion between the unhydrated oxo and gem-diol forms, where two hydroxyl groups replace the carbonyl group. We here examine the hydration equilibrium of glyoxylic and pyruvic acids with first-principles simulations in water at ambient conditions using ab initio metadynamics to reconstruct the corresponding free-energy landscapes. The main results are as follows: (i) our simulations reveal the high conformational diversity of these species in aqueous solutions. (ii) We show that gem-diol is strongly favored in water compared to its oxo counterpart by 29 and 16 kJ/mol for glyoxylic and pyruvic acids, respectively. (iii) From our atomic-scale simulations, we present new insights into the reaction mechanisms with a special focus on hydrogen-bond arrangements and the electronic structure of the transition state.

Effects of Doped N, B, P, and S Atoms on Graphene toward Oxygen Evolution Reactions

Molecular oxygen and hydrogen can be obtained from the water-splitting process through the electrolysis technique. However, harnessing energy is very challenging in this way due to the involvement of the 4e- reaction pathway, which is associated with a substantial amount of reaction barrier. After the report of the first N-doped graphene acting as an oxygen reduction reaction catalyst, the scientific community set out on exploring more reliable doping materials, better material engineering techniques, and developing computational models to explain the interfacial reactions. In this study, we modeled the graphene surface with four different nonmetal doping atoms N, B, P, and S individually by replacing a carbon atom from one of the graphitic positions. We report the mechanism of the complete catalytic cycle for each of the doped surfaces by the doping atom. The energy barriers for individual steps were explored using the biased first-principles molecular dynamics simulations to overcome the high reaction barrier. We explain the active sites and provide a comparison between the activation energy obtained by the application of two computational methods. Observing the rate-determining step, that is, oxo-oxo bond formation, S-doped graphene is the most effective. In contrast, N-doped graphene seems to be the least useful for oxygen evolution catalysis compared to the undoped graphene surface. B-doped graphene and P-doped graphene have an equivalent impact on the catalytic cycle.

Nitromethane Decomposition via Automated Reaction Discovery and an Ab Initio Corrected Kinetic Model

We explore the systematic construction of kinetic models from in silico reaction data for the decomposition of nitromethane. Our models are constructed in a computationally affordable manner by using reactions discovered through accelerated molecular dynamics simulations using the ReaxFF reactive force field. The reaction paths are then optimized to determine reaction rate parameters. We introduce a reaction barrier correction scheme that combines accurate thermochemical data from density functional theory with ReaxFF minimal energy paths. We validate our models across different thermodynamic regimes, showing predictions of gas phase CO and NO concentrations and high-pressure induction times that are similar to experimental data. The kinetic models are analyzed to find fundamental decomposition reactions in different thermodynamic regimes.

Ab Initio Molecular Dynamics Simulations of the SN1/SN2 Mechanistic Continuum in Glycosylation Reactions

We report a computational approach to evaluate the reaction mechanisms of glycosylation using ab initio molecular dynamics (AIMD) simulations in explicit solvent. The reaction pathways are simulated via free energy calculations based on metadynamics and trajectory simulations using Born-Oppenheimer molecular dynamics. We applied this approach to investigate the mechanisms of the glycosylation of glucosyl α-trichloroacetimidate with three acceptors (EtOH, i-PrOH, and t-BuOH) in three solvents (ACN, DCM, and MTBE). The reactants and the solvents are treated explicitly using density functional theory. We show that the profile of the free energy surface, the synchronicity of the transition state structure, and the time gap between leaving group dissociation and nucleophile association can be used as three complementary indicators to describe the glycosylation mechanism within the SN1/SN2 continuum for a given reaction. This approach provides a reliable means to rationalize and predict reaction mechanisms and to estimate lifetimes of oxocarbenium intermediates and their dependence on the glycosyl donor, acceptor, and solvent environment.

Influence of surface defects on activity and selectivity: a quantitative study of structure sensitivity of Pd catalysts for acetylene hydrogenation

The structure sensitivity of Pd catalysed acetylene hydrogenation is quantitatively examined using a coverage-dependent microkinetic model. Pd(211) was found to be more active than Pd(111), but present a poorer selectivity toward ethylene.

Proton transfer from water to aromatic N-heterocyclic anions from DFT-MD simulations

The phenomenon of proton transfer from water to six N-heterocyclic anions and free energy landscapes of this process are studied using both electronic structure calculations and first principles molecular metadynamics simulations. Our investigation involves microhydrated and aqueous phase interaction of water with six aromatic heterocyclic anions relevant to chemistry and biology: imidazolide, pyrrolide, benzimidazolide, 2-cyanopyrrolide, indolide, and indazolide. The basic structures of all these heterocyclic anions differ by substituted functional groups as well as fused rings. We study the proton transfer reaction and the minimum number of required water molecules for the reaction in hydrated microclusters. We find out that at least four water molecules are necessary for hydrated clusters to facilitate the intracluster proton transfer reaction from water to anions except for pyrrolide, for which this magic number is 3. To obtain the reaction free energy and activation barrier of the proton transfer process in an aqueous solution, the metadynamics method based first principles molecular dynamics simulations were performed. The complete proton transfer was observed in aqueous solutions for all the anions. The water molecule directly involved in proton transfer becomes acidic due to the cooperative effect of neighboring water molecules. From the metadynamics simulation, we obtain the values of activation barrier for the proton transfer processes from neutral water to anions, and the highest activation barrier is obtained for benzimidazolide, whereas the lowest activation barrier is obtained for pyrrolide. The structures and free energy profiles of the process for all the anions are discussed, and a comparative outlook of the study is presented here.

A QM/MM simulation study of transamination reaction at the active site of aspartate aminotransferase: Free energy landscape and proton transfer pathways

Transaminase is a key enzyme for amino acid metabolism, which reversibly catalyzes the transamination reaction with the help of PLP (pyridoxal 5^' -phosphate) as its cofactor. Here we have investigated the mechanism and free energy landscape of the transamination reaction involving the aspartate transaminase (AspTase) enzyme and aspartate-PLP (Asp-PLP) complex using QM/MM simulation and metadynamics methods. The reaction is found to follow a stepwise mechanism where the active site residue Lys258 acts as a base to shuttle a proton from α-carbon (CA) to imine carbon (C4A) of the PLP-Asp Schiff base. In the first step, the Lys258 abstracts the CA proton of the substrate leading to the formation of a carbanionic intermediate which is followed by the reprotonation of the Asp-PLP Schiff base at C4A atom by Lys258. It is found that the free energy barrier for the proton abstraction by Lys258 and that for the reprotonation are 17.85 and 3.57 kcal/mol, respectively. The carbanionic intermediate is 7.14 kcal/mol higher in energy than the reactant. Hence, the first step acts as the rate limiting step. The present calculations also show that the Lys258 residue undergoes a conformational change after the first step of transamination reaction and becomes proximal to C4A atom of the Asp-PLP Schiff base to favor the second step. The active site residues Tyr70* and Gly38 anchor the Lys258 in proper position and orientation during the first step of the reaction and stabilize the positive charge over Lys258 generated at the intermediate step.

Exploring Free Energy Profiles of Enantioselective Organocatalytic Aldol Reactions under Full Solvent Influence

We present a computational study on the enantioselectivity of organocatalytic proline-catalyzed aldol reactions between aldehydes in dimethylformamide (DMF). To explore the free energy surface of the reaction, we apply two-dimensional metadynamics on top of ab initio molecular dynamics (AIMD) simulations with explicit solvent description on the DFT level of theory. We avoid unwanted side reactions by utilizing our newly developed hybrid AIMD (HyAIMD) simulation scheme, which adds a simple force field to the AIMD simulation to prevent unwanted bond breaking and formation. Our condensed phase simulation results are able to nicely reproduce the experimental findings, including the main stereoisomer that is formed, and give a correct qualitative prediction of the change in syn:anti product ratio with different substituents. Furthermore, we give a microscopic explanation for the selectivity. We show that both the explicit description of the solvent and the inclusion of entropic effects are vital to a good outcome-metadynamics simulations in vacuum and static nudged elastic band (NEB) calculations yield significantly worse predictions when compared to the experiment. The approach described here can be applied to a plethora of other enantioselective or organocatalytic reactions, enabling us to tune the catalyst or determine the solvent with the highest stereoselectivity.