Competing Reaction Mechanisms of Peptide Bond Formation in Water Revealed by Deep Potential Molecular Dynamics and Path Sampling.

The formation of an amide bond is an essential step in the synthesis of materials and drugs, and in the assembly of amino acids to form peptides. The mechanism of this reaction has been studied extensively, in particular to understand how it can be catalyzed, but a representation capable of explaining all the experimental data is still lacking. Numerical simulation should provide the necessary molecular description, but the solvent involvement poses a number of challenges. Here, we combine the efficiency and accuracy of neural network potential-based reactive molecular dynamics with the extensive and unbiased exploration of reaction pathways provided by transition path sampling. Using microsecond-scale simulations at the density functional theory level, we show that this method reveals the presence of two competing distinct mechanisms for peptide bond formation between alanine esters in aqueous solution. We describe how both reaction pathways, via a general base catalysis mechanism and via direct cleavage of the tetrahedral intermediate respectively, change with pH. This result contrasts with the conventional mechanism involving a single pathway in which only the barrier heights are affected by pH. We show that this new proposal involving two competing mechanisms is consistent with the experimental data, and we discuss the implications for peptide bond formation under prebiotic conditions and in the ribosome. Our work shows that integrating deep potential molecular dynamics with path sampling provides a powerful approach for exploring complex chemical mechanisms.

Understanding the Interactions of Ubiquitin-Specific Protease 7 with Its Substrates through Molecular Dynamics Simulations: Insights into the Role of Its C-Terminal Domains in Substrate Recognition.

Ubiquitin-specific protease 7 (USP7) is a deubiquitinase enzyme that plays a critical role in regulating various cellular processes by cleaving ubiquitin molecules from target proteins. The C-terminal loop (CTL) motif is a specific region at the C-terminal end of the USP7 enzyme. Recent experiments suggest that the CTL motif plays a role in USP7's catalytic activity by contributing to the enzyme's structural stability, substrate recognition, and catalytic efficiency. The objective of this work is to elucidate these roles through the utilization of computational methods for molecular simulations. For this, we conducted extensive molecular dynamics (MD) simulations to investigate the conformational dynamics and protein-protein interactions within the USP7 enzyme-substrate complex with the substrate consisting of the ubiquitin tagged with the fluorescent label rhodamine 110-gly (Ub-Rho). Our results demonstrate that the CTL motif plays a crucial role in stabilizing the Ubl domains' conformation and augmenting the stability of active conformations within the enzyme-substrate complex. Conversely, the absence of the CTL motif results in increased flexibility and variability in Ubl domains' motion, leading to a reduced percentage of active conformations. Furthermore, our analysis of protein-protein interactions highlights the significance of the CTL motif in anchoring the Ubl45 domains to the catalytic domain (CD), thereby facilitating stable interactions with the substrate. Overall, our findings provide valuable insights into the conformational dynamics and protein-protein interactions inherent in the USP7 enzyme-substrate complex. These insights shed light on some mechanistic details of USP7 concerning the substrate's recognition before its catalytic action.

Pomotrelvir and Nirmatrelvir Binding and Reactivity with SARS-CoV-2 Main Protease: Implications for Resistance Mechanisms from Computations.

We investigate the inhibition mechanism between pomotrelvir and the SARS-CoV-2 main protease using molecular mechanics and quantum mechanics/molecular mechanics simulations. Alchemical transformations where each Pi group of pomotrelvir was transformed into its counterpart in nirmatrelvir were performed to unravel the individual contribution of each group to the binding and reaction processes. We have shown that while a γ-lactam ring is preferred at position P1, a δ-lactam ring is a good alternative for the design of inhibitors for variants presenting mutations at position 166. For the P2 position, tertiary amines are preferred with respect to secondary amines. Flexible side chains at the P2 position can disrupt the preorganization of the active site, favouring the exploration of non-reactive conformations. The substitution of the P2 group of pomotrelvir by that of nirmatrelvir resulted in a compound, named as C2, that presents a better binding free energy and a higher population of reactive conformations in the Michaelis complex. Analysis of the chemical reaction to form the covalent complex has shown a similar reaction mechanism and activation free energies for pomotrelvir, nirmatrelvir and C2. We hope that these findings could be useful to design better inhibitors to fight present and future variants of the SARS-CoV-2 virus.

Insights into the Enhancement of the Poly(ethylene terephthalate) Degradation by FAST-PETase from Computational Modeling.

Polyethylene terephthalate (PET) is the most abundant polyester plastic, widely used in textiles and packaging, but, unfortunately, it is also one of the most discarded plastics after one use. In the last years, the enzymatic biodegradation of PET has sparked great interest owing to the discovery and subsequent mutation of PETase-like enzymes, able to depolymerize PET. FAST-PETase is one of the best enzymes hitherto proposed to efficiently degrade PET, although the origin of its efficiency is not completely clear. To understand the molecular origin of its enhanced catalytic activity, we have carried out a thorough computational study of PET degradation by the FAST-PETase action by employing classical and hybrid (QM/MM) molecular dynamics (MD) simulations. Our findings show that the rate-limiting reaction step for FAST-PETase corresponds to the acylation stage with an estimated free energy barrier of 12.1 kcal mol-1, which is significantly smaller than that calculated for PETase (16.5 kcal mol-1) and, therefore, supports the enhanced catalytic activity of FAST-PETase. The origin of this enhancement is mainly attributed to the N233K mutation, which, although sited relatively far from the active site, induces a chain folding where the Asp206 of the catalytic triad is located, impeding that this residue sets effective H-bonds with its neighboring residues. This effect makes Asp206 hold a more basic character compared to the wild-type PETase and boosts the interaction with the protonated His237 of the catalytic triad in the transition state of acylation, with the consequent decrease of the catalytic barrier and acceleration of the PET degradation reaction.

Conformational Changes and ATP Hydrolysis in Zika Helicase: The Molecular Basis of a Biomolecular Motor Unveiled by Multiscale Simulations.

We computationally study the Zika NS3 helicase, a biological motor, using ATP hydrolysis energy for nucleic acid remodeling. Through molecular mechanics and hybrid quantum mechanics/molecular mechanics simulations, we explore the conformational landscape of motif V, a conserved loop connecting the active sites for ATP hydrolysis and nucleic acid binding. ATP hydrolysis, initiated by a meta-phosphate group formation, involves the nucleophilic attack of a water molecule activated by Glu286 proton abstraction. Motif V hydrogen bonds to this water via the Gly415 backbone NH group, assisting hydrolysis. Posthydrolysis, free energy is released when the inorganic phosphate moves away from the coordination shell of the magnesium ion, inducing a significant shift in the conformational landscape of motif V to establish a hydrogen bond between the Gly415 NH group and Glu285. According to our simulations, the Zika NS3 helicase acts as a ratchet biological motor with motif V transitions steered by Gly415's γ-phosphate sensing in the ATPase site.

Elucidation of the Active Form and Reaction Mechanism in Human Asparaginase Type III Using Multiscale Simulations.

l-asparaginases catalyze the asparagine hydrolysis to aspartate. These enzymes play an important role in the treatment of acute lymphoblastic leukemia because these cells are unable to produce their own asparagine. Due to the immunogenic response and various side effects of enzymes of bacterial origin, many attempts have been made to replace these enzymes with mammalian enzymes such as human asparaginase type III (hASNaseIII). This study investigates the reaction mechanism of hASNaseIII through molecular dynamics simulations, quantum mechanics/molecular mechanics methods, and free energy calculations. Our simulations reveal that the dimeric form of the enzyme plays a vital role in stabilizing the substrate in the active site, despite the active site residues coming from a single protomer. Protomer-protomer interactions are essential to keep the enzyme in an active conformation. Our study of the reaction mechanism indicates that the self-cleavage process that generates an N-terminal residue (Thr168) is required to activate the enzyme. This residue acts as the nucleophile, attacking the electrophilic carbon of the substrate after a proton transfer from its hydroxyl group to the N-terminal amino group. The reaction mechanism proceeds with the formation of an acyl-enzyme complex and its hydrolysis, which turns out to be the rate-determining step. Our proposal of the enzymatic mechanism sheds light on the role of different active site residues and rationalizes the studies on mutations. The insights provided here about hASNaseIII activity could contribute to the comprehension of the disparities among different ASNases and might even guide the design of new variants with improved properties for acute lymphoblastic leukemia treatment.

Spatial and Temporal Resolution of the Oxygen-Independent Photoinduced DNA Interstrand Cross-Linking by a Nitroimidazole Derivative.

DNA damage is ubiquitous in nature and is at the basis of emergent treatments such as photodynamic therapy, which is based on the activation of highly oxidative reactive oxygen species by photosensitizing O2. However, hypoxia observed in solid tumors imposes the necessity to devise oxygen-independent modes of action able to induce DNA damage under a low oxygen concentration. The complexity of these DNA damage mechanisms in realistic environments grows exponentially when taking into account light absorption and subsequent excited-state population, photochemical and (photo)-redox reactions, the multiple species involved in different electronic states, noncovalent interactions, multiple reaction steps, and the large number of DNA reactive sites. This work tackles all the intricate reactivity of a photosensitizer based on a nitroimidazole derivative reacting toward DNA in solution under UV light exposition. This is performed through a combination of ground- and excited-state quantum chemistry, classical molecular dynamics, and hybrid QM/MM simulations to rationalize in detail the formation of DNA interstrand cross-links (ICLs) exerted by the noncanonical noncovalent photosensitizer. Unprecedented spatial and temporal resolution of these phenomena is achieved, revealing that the ICL is sequence-specific and that the fastest reactions take place at AT, GC, and GT steps involving either the opposite nucleobases or adjacent Watson-Crick base pairs. The N7 and O6 positions of guanine, the N7 and N3 sites of adenine, the N4 position of cytosine, and the O2 atom of thymine are deemed as the most nucleophile sites and are positively identified to participate in the ICL productions. This work provides a multiscale computational protocol to study DNA reactivity with noncovalent photosensitizers, and contributes to the understanding of therapies based on photoinduced DNA damage at molecular and electronic levels. In addition, we believe the depth understanding of these processes should assist the design of new photosensitizers considering their molecular size, electronic properties, and the observed regioselectivity toward nucleic acids.

Unraveling the Role of the Tyrosine Tetrad from the Binding Site of the Epigenetic Writer MLL3 in the Catalytic Mechanism and Methylation Multiplicity.

MLL3, also known as KMT2C, is a lysine mono-methyltransferase in charge of the writing of an epigenetic mark on lysine 4 from histone 3. The catalytic site of MLL3 is composed of four tyrosines, namely, Y44, Y69, Y128, and Y130. Tyrosine residues are highly conserved among lysine methyltransferases' catalytic sites, although their complete function is still unclear. The exploration of how modifications on these residues from the enzymatic machinery impact the enzymatic activity of MLL3 could shed light transversally into the inner functioning of enzymes with similar characteristics. Through the use of QMMM calculations, we focus on the effect of the mutation of each tyrosine from the catalytic site on the enzymatic activity and the product specificity in the current study. While we found that the mutations of Y44 and Y128 by phenylalanine inactivated the enzyme, the mutation of Y128 by alanine reactivated the enzymatic activity of MLL3. Moreover, according to our models, the Y128A mutant was even found to be capable of di- and tri-methylate lysine 4 from histone 3, what would represent a gain of function mutation, and could be responsible for the development of diseases. Finally, we were able to establish the inactivation mechanism, which involved the use of Y130 as a water occlusion structure, whose conformation, once perturbed by its mutation or Y128 mutant, allows the access of water molecules that sequester the electron pair from lysine 4 avoiding its methylation process and, thus, increasing the barrier height.

Mechanistic study of the biosynthesis of R-phenylacetylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations.

The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the Cß atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4'-aminopyrimidium (APH+) form. The calculated activation free energy for this step is 17.4 kcal mol-1 at 27 °C. From this point, the reaction continues with the abstraction of Hß atom of the HEThDP intermediate by the Oß atom of benzaldehyde to form the intermediate I. The reaction is completed with the cleavage of the bond C2α-C2 to form the product R-PAC and to regenerate the ylide intermediate under the APH+ form, allowing in this way to reinitiate to the catalytic cycle once more. The calculated activation barrier for this last step is 15.9 kcal mol-1 at 27 °C.

Mechanistic study of the biosynthesis of R-phenylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations.

The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the Cß atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4'-aminopyrimidium (APH+) form. The calculated activation free energy for this step is 17.4kcal mol-1 at 27 °C. From this point, the reaction continues with the abstraction of Hß atom of the HEThDP intermediate by the Oß atom of benzaldehyde to form the intermediate I. The reaction is completed with the cleavage of the bond C2α-C2 to form the product R-PAC and to regenerate the ylide intermediate under the APH+ form, allowing in this way to reinitiate to the catalytic cycle once more. The calculated activation barrier for this last step is 15.9kcal mol-1 at 27 °C.

Targeting the JAK/STAT Pathway: A Combined Ligand- and Target-Based Approach.

Janus kinases (JAKs) are a family of proinflammatory enzymes able to mediate the immune responses and the inflammatory cascade by modulating multiple cytokine expressions as well as various growth factors. In the present study, the inhibition of the JAK-signal transducer and activator of transcription (STAT) signaling pathway is explored as a potential strategy for treating autoimmune and inflammatory disorders. A computationally driven approach aimed at identifying novel JAK inhibitors based on molecular topology, docking, and molecular dynamics simulations was carried out. For the best candidates selected, the inhibitory activity against JAK2 was evaluated in vitro. Two hit compounds with a novel chemical scaffold, 4 (IC50 = 0.81 µM) and 7 (IC50 = 0.64 µM), showed promising results when compared with the reference drug Tofacitinib (IC50 = 0.031 µM).

Inhibition Mechanism of SARS-CoV-2 Main Protease with Ketone-Based Inhibitors Unveiled by Multiscale Simulations: Insights for Improved Designs*.

We present the results of classical and QM/MM simulations for the inhibition of SARS-CoV-2 3CL protease by a hydroxymethylketone inhibitor, PF-00835231. In the noncovalent complex the carbonyl oxygen atom of the warhead is placed in the oxyanion hole formed by residues 143 to 145, while P1-P3 groups are accommodated in the active site with interactions similar to those observed for the peptide substrate. According to alchemical free energy calculations, the P1' hydroxymethyl group also contributes to the binding free energy. Covalent inhibition of the enzyme is triggered by the proton transfer from Cys145 to His41. This step is followed by the nucleophilic attack of the Sγ atom on the carbonyl carbon atom of the inhibitor and a proton transfer from His41 to the carbonyl oxygen atom mediated by the P1' hydroxyl group. Computational simulations show that the addition of a chloromethyl substituent to the P1' group may lower the activation free energy for covalent inhibition.

Quantifying the limits of transition state theory in enzymatic catalysis.

While being one of the most popular reaction rate theories, the applicability of transition state theory to the study of enzymatic reactions has been often challenged. The complex dynamic nature of the protein environment raised the question about the validity of the nonrecrossing hypothesis, a cornerstone in this theory. We present a computational strategy to quantify the error associated to transition state theory from the number of recrossings observed at the equicommittor, which is the best possible dividing surface. Application of a direct multidimensional transition state optimization to the hydride transfer step in human dihydrofolate reductase shows that both the participation of the protein degrees of freedom in the reaction coordinate and the error associated to the nonrecrossing hypothesis are small. Thus, the use of transition state theory, even with simplified reaction coordinates, provides a good theoretical framework for the study of enzymatic catalysis.

Heavy Enzymes and the Rational Redesign of Protein Catalysts.

An unsolved mystery in biology concerns the link between enzyme catalysis and protein motions. Comparison between isotopically labelled "heavy" dihydrofolate reductases and their natural-abundance counterparts has suggested that the coupling of protein motions to the chemistry of the catalysed reaction is minimised in the case of hydride transfer. In alcohol dehydrogenases, unnatural, bulky substrates that induce additional electrostatic rearrangements of the active site enhance coupled motions. This finding could provide a new route to engineering enzymes with altered substrate specificity, because amino acid residues responsible for dynamic coupling with a given substrate present as hotspots for mutagenesis. Detailed understanding of the biophysics of enzyme catalysis based on insights gained from analysis of "heavy" enzymes might eventually allow routine engineering of enzymes to catalyse reactions of choice.

A first peek into sub-picosecond dynamics of spin energy levels in magnetic biomolecules.

We estimate the time- and temperature-evolution of spin energy levels in a metallopeptide by combining molecular dynamics with crystal field analysis. Fluctuations of tens of cm-1 for spin energy levels at fs times gradually average out at longer times. We confirm that local vibrations are key in spin dynamics.

Insights on the Origin of Catalysis on Glycine N-Methyltransferase from Computational Modeling.

The origin of enzyme catalysis remains a question of debate despite much intense study. We report a QM/MM theoretical study of the SN2 methyl transfer reaction catalyzed by a glycine N-methyltransferase (GNMT) and three mutants to test whether recent experimental observations of rate-constant reductions and variations in inverse secondary α-3H kinetic isotope effects (KIEs) should be attributed to changes in the methyl donor-acceptor distance (DAD): Is catalysis due to a compression effect? Semiempirical (AM1) and DFT (M06-2X) methods were used to describe the QM subset of atoms, while OPLS-AA and TIP3P classical force fields were used for the protein and water molecules, respectively. The computed activation free energies and KIEs are in good agreement with experimental data, but the mutations do not meaningfully affect the DAD: Compression cannot explain the experimental variations on KIEs. On the contrary, electrostatic properties in the active site correlate with the catalytic activity of wild type and mutants. The plasticity of the enzyme moderates the effects of the mutations, explaining the rather small degree of variation in KIEs and reactivities.

Isotope Substitution of Promiscuous Alcohol Dehydrogenase Reveals the Origin of Substrate Preference in the Transition State.

The origin of substrate preference in promiscuous enzymes was investigated by enzyme isotope labelling of the alcohol dehydrogenase from Geobacillus stearothermophilus (BsADH). At physiological temperature, protein dynamic coupling to the reaction coordinate was insignificant. However, the extent of dynamic coupling was highly substrate-dependent at lower temperatures. For benzyl alcohol, an enzyme isotope effect larger than unity was observed, whereas the enzyme isotope effect was close to unity for isopropanol. Frequency motion analysis on the transition states revealed that residues surrounding the active site undergo substantial displacement during catalysis for sterically bulky alcohols. BsADH prefers smaller substrates, which cause less protein friction along the reaction coordinate and reduced frequencies of dynamic recrossing. This hypothesis allows a prediction of the trend of enzyme isotope effects for a wide variety of substrates.

Revealing the Origin of the Efficiency of the De Novo Designed Kemp Eliminase HG-3.17 by Comparison with the Former Developed HG-3.

The design of new biocatalysts is a goal in biotechnology to improve the rate, selectivity and environmental impact of industrial chemical processes. In this regard, the use of computational techniques has provided valuable assistance in the design of new enzymes with remarkable catalytic activity. In this paper, hybrid QM/MM molecular dynamics simulations have allowed insights to be gained on the origin of the limited efficiency of a computationally designed enzyme for the Kemp elimination; the HG-3. Comparison of results derived from this enzyme with those of a more evolved protein containing additional point mutations, HG-3.17, rendered important information that should be taken into account in the design of new enzymes. For this Kemp eliminase reaction, higher reactivity has been demonstrated to be related to a better electrostatic preorganisation of an environment that creates a more favourable electrostatic potential for the reaction to proceed. The limitations of HG-3 can be related to a lack of flexibility, a not well-fitted active site, and a lack of protein electrostatic preorganisation, which decrease the reorganisation around the oxyanion hole.

Adaptive Finite Temperature String Method in Collective Variables.

Here we present a modified version of the on-the-fly string method for the localization of the minimum free energy path in a space of arbitrary collective variables. In the proposed approach the shape of the biasing potential is controlled by only two force constants, defining the width of the potential along the string and orthogonal to it. The force constants and the distribution of the string nodes are optimized during the simulation, improving the convergence. The optimized parameters can be used for umbrella sampling with a path CV along the converged string as the reaction coordinate. We test the new method with three fundamentally different processes: chloride attack to chloromethane in bulk water, alanine dipeptide isomerization, and the enzymatic conversion of isochorismate to piruvate. In each case the same set of parameters resulted in a rapidly converging simulation and a precise estimation of the potential of mean force. Therefore, the default settings can be used for a wide range of processes, making the method essentially parameter free and more user-friendly.

Singlet Oxygen Attack on Guanine: Reactivity and Structural Signature within the B-DNA Helix.

Oxidatively generated DNA lesions are numerous and versatile, and have been the subject of intensive research since the discovery of 8-oxoguanine in 1984. Even for this prototypical lesion, the precise mechanism of formation remains elusive due to the inherent difficulties in characterizing high-energy intermediates. We have probed the stability of the guanine endoperoxide in B-DNA as a key intermediate and determined a unique activation free energy of around 6 kcal mol(-1) for the formation of the first C-O covalent bond upon the attack of singlet molecular oxygen ((1) O2 ) on the central guanine of a solvated 13 base-pair poly(dG-dC), described by means of quantum mechanics/molecular mechanics (QM/MM) simulations. The B-helix remains stable upon oxidation in spite of the bulky character of the guanine endoperoxide. Our modeling study has revealed the nature of the versatile (1) O2 attack in terms of free energy and shows a sensitivity to electrostatics and solvation as it involves a charge-separated intermediate.