Secondary Amine Catalysis in Enzyme Design: Broadening Protein Template Diversity through Genetic Code Expansion

Secondary amines, due to their reactivity, can transform protein templates into catalytically active entities, accelerating the development of artificial enzymes. However, existing methods, predominantly reliant on modified ligands or N-terminal prolines, impose significant limitations on template selection. In this study, genetic code expansion was used to break this boundary, enabling secondary amines to be incorporated into alternative proteins and positions of choice. Pyrrolysine analogues carrying different secondary amines could be incorporated into superfolder green fluorescent protein (sfGFP), multidrug-binding LmrR and nucleotide-binding dihydrofolate reductase (DHFR). Notably, the analogue containing a D-proline moiety demonstrated both proteolytic stability and catalytic activity, conferring LmrR and DHFR with the desired transfer hydrogenation activity. While the LmrR variants were confined to the biomimetic 1-benzyl-1,4-dihydronicotinamide (BNAH) as the hydride source, the optimal DHFR variant favorably used the pro-R hydride from NADPH for stereoselective reactions (e.r. up to 92 : 8), highlighting that a switch of protein template could broaden the nucleophile option for catalysis. Owing to the cofactor compatibility, the DHFR-based secondary amine catalysis could be integrated into an enzymatic recycling scheme. This established method shows substantial potential in enzyme design, applicable from studies on enzyme evolution to the development of new biocatalysts.

A Single Amino Acid Able to Promote High-Temperature Ring-Opening Polymerization by Dual Activation

Amino acids are indispensable compounds in the body, performing several biological processes that enable proper functioning. In this work, it is demonstrated that a single amino acid, taurine, is also able to promote the ring-opening polymerization (ROP) of several cyclic monomers under industrially relevant conditions. It is shown that the unique zwitterionic structure of taurine, where the negatively charged sulfonic acid group and the protonated amine group are separated by two methylene groups, not only provides high thermal stability but also leads to a dual activation mechanism, which is corroborated by quantum mechanical calculations. This unique mechanism allows for the synthesis of polylactide of up to 50 kDa in bulk at 180 °C with good end-group fidelity using a highly abundant catalyst. Furthermore, cytotoxicity tests confirm that PLLA synthesized with taurine is non-toxic. Moreover, it is demonstrated that the presence of taurine does not have any detrimental effect on the thermal stability of polylactide, and therefore polymers can be used directly without any post-polymerization purification. It is believed that the demonstration that a simple structure composed of a single amino acid can promote polymerization can bring a paradigm shift in the preparation of polymers.

Electrostatics as a Guiding Principle in Understanding and Designing Enzymes

Enzyme design faces challenges related to the implementation of the basic principles that govern the catalytic activity in natural enzymes. In this work, we revisit basic electrostatic concepts that have been shown to explain the origin of enzymatic efficiency like preorganization and reorganization. Using magnitudes such as the electrostatic potential and the electric field generated by the protein, we explain how these concepts work in different enzymes and how they can be used to rationalize the consequences of point mutations. We also discuss examples of protein design in which electrostatic effects have been implemented. For the near future, molecular simulations, coupled with the use of machine learning methods, can be used to implement electrostatics as a guiding principle for enzyme designs.

Peptidyl nitroalkene inhibitors of main protease rationalized by computational and crystallographic investigations as antivirals against SARS-CoV-2

The coronavirus disease 2019 (COVID-19) pandemic continues to represent a global public health issue. The viral main protease (Mpro) represents one of the most attractive targets for the development of antiviral drugs. Herein we report peptidyl nitroalkenes exhibiting enzyme inhibitory activity against Mpro (Ki: 1-10 μM) good anti-SARS-CoV-2 infection activity in the low micromolar range (EC50: 1-12 μM) without significant toxicity. Additional kinetic studies of compounds FGA145, FGA146 and FGA147 show that all three compounds inhibit cathepsin L, denoting a possible multitarget effect of these compounds in the antiviral activity. Structural analysis shows the binding mode of FGA146 and FGA147 to the active site of the protein. Furthermore, our results illustrate that peptidyl nitroalkenes are effective covalent reversible inhibitors of the Mpro and cathepsin L, and that inhibitors FGA145, FGA146 and FGA147 prevent infection against SARS-CoV-2.

Impact of the Warhead of Dipeptidyl Keto Michael Acceptors on the Inhibition Mechanism of Cysteine Protease Cathepsin L

Cathepsin L (CatL) is a lysosomal cysteine protease whose activity has been related to several human pathologies. However, although preclinical trials using CatL inhibitors were promising, clinical trials have been unsuccessful up to now. We are presenting a study of two designed dipeptidyl keto Michael acceptor potential inhibitors of CatL with either a keto vinyl ester or a keto vinyl sulfone (KVS) warhead. The compounds were synthesized and experimentally assayed in vitro, and their inhibition molecular mechanism was explored based on molecular dynamics simulations at the density functional theory/molecular mechanics level. The results confirm that both compounds inhibit CatL in the nanomolar range and show a time-dependent inhibition. Interestingly, despite both presenting almost equivalent equilibrium constants for the reversible formation of the noncovalent enzyme/inhibitor complex, differences are observed in the chemical step corresponding to the enzyme-inhibitor covalent bond formation, results that are mirrored by the computer simulations. Theoretically determined kinetic and thermodynamic results, which are in very good agreement with the experiments, afford a detailed explanation of the relevance of the different structural features of both compounds having a significant impact on enzyme inhibition. The unprecedented binding interactions of both inhibitors in the P1' site of CatL represent valuable information for the design of inhibitors. In particular, the peptidyl KVS can be used as a starting lead compound in the development of drugs with medical applications for the treatment of cancerous pathologies since sulfone warheads have previously shown promising cell stability compared to other functions such as carboxylic esters. Future improvements can be guided by the atomistic description of the enzyme-inhibitor interactions established along the inhibition reaction derived from computer simulations.

Mechanistic studies of a lipase unveil effect of pH on hydrolysis products of small PET modules

Biocatalysis is a key technology enabling plastic recycling. However, despite advances done in the development of plastic-degrading enzymes, the molecular mechanisms that govern their catalytic performance are poorly understood, hampering the engineering of more efficient enzyme-based technologies. In this work, we study the hydrolysis of PET-derived diesters and PET trimers catalyzed by the highly promiscuous lipase B from Candida antarctica (CALB) through QM/MM molecular dynamics simulations supported by experimental Michaelis-Menten kinetics. The computational studies reveal the role of the pH on the CALB regioselectivity toward the hydrolysis of bis-(hydroxyethyl) terephthalate (BHET). We exploit this insight to perform a pH-controlled biotransformation that selectively hydrolyzes BHET to either its corresponding diacid or monoesters using both soluble and immobilized CALB. The discoveries presented here can be exploited for the valorization of BHET resulting from the organocatalytic depolymerization of PET.

Impact of the Recognition Part of Dipeptidyl Nitroalkene Compounds on the Inhibition Mechanism of Cysteine Proteases Cruzain and Cathepsin L

Cysteine proteases (CPs) are an important class of enzymes, many of which are responsible for several human diseases. For instance, cruzain of protozoan parasite Trypanosoma cruzi is responsible for the Chagas disease, while the role of human cathepsin L is associated with some cancers or is a potential target for the treatment of COVID-19. However, despite paramount work carried out during the past years, the compounds that have been proposed so far show limited inhibitory action against these enzymes. We present a study of proposed covalent inhibitors of these two CPs, cruzain and cathepsin L, based on the design, synthesis, kinetic measurements, and QM/MM computational simulations on dipeptidyl nitroalkene compounds. The experimentally determined inhibition data, together with the analysis and the predicted inhibition constants derived from the free energy landscape of the full inhibition process, allowed describing the impact of the recognition part of these compounds and, in particular, the modifications on the P2 site. The designed compounds and, in particular, the one with a bulky group (Trp) at the P2 site show promising in vitro inhibition activities against cruzain and cathepsin L for use as a starting lead compound in the development of drugs with medical applications for the treatment of human diseases and future designs.

Molecular Basis of RNA-Driven ATP Hydrolysis in DExH-Box Helicases

The spliceosome machinery catalyzes precursor messenger (pre-m)RNA splicing. In each cycle, the spliceosome experiences massive compositional and conformational remodeling fueled by the concerted action of specific RNA-dependent ATPases/helicases. Intriguingly, these enzymes are allosterically activated to perform ATP hydrolysis and trigger helicase activity only upon pre-mRNA binding. Yet, the molecular mechanism underlying the RNA-driven regulation of their ATPase function remains elusive. Here, we focus on the Prp2 ATPase/helicase which contributes to reshaping the spliceosome into its catalytic competent state. By performing classical and quantum-classical molecular dynamics simulations, we unprecedentedly unlock the molecular terms governing the Prp2 ATPase/helicase function. Namely, we dissect the molecular mechanism of ATP hydrolysis, and we disclose that RNA binding allosterically triggers the formation of a set of interactions linking the RNA binding tunnel to the catalytic site. This activates the Prp2's ATPase function by optimally placing the nucleophilic water and the general base of the enzymatic process to perform ATP hydrolysis. The key structural motifs, mechanically coupling RNA gripping and the ATPase/helicase functions, are conserved across all DExH-box helicases. This mechanism could thus be broadly applicable to all DExH-box helicase family.

Engineering a Highly Regioselective Fungal Peroxygenase for the Synthesis of Hydroxy Fatty Acids

The hydroxylation of fatty acids is an appealing reaction in synthetic chemistry, although the lack of selective catalysts hampers its industrial implementation. In this study, we have engineered a highly regioselective fungal peroxygenase for the ω-1 hydroxylation of fatty acids with quenched stepwise over-oxidation. One single mutation near the Phe catalytic tripod narrowed the heme cavity, promoting a dramatic shift toward subterminal hydroxylation with a drop in the over-oxidation activity. While crystallographic soaking experiments and molecular dynamic simulations shed light on this unique oxidation pattern, the selective biocatalyst was produced by Pichia pastoris at 0.4 g L^-1 in a fed-batch bioreactor and used in the preparative synthesis of 1.4 g of (ω-1)-hydroxytetradecanoic acid with 95 % regioselectivity and 83 % ee for the S enantiomer.

Computational Study of the Inhibition of RgpB Gingipain, a Promising Target for the Treatment of Alzheimer's Disease

Alzheimer's disease represents one of the most ambitious challenges for biomedical sciences due to the growing number of cases worldwide in the elderly population and the lack of efficient treatments. One of the recent attempts to develop a treatment points to the cysteine protease RgpB as a promising drug target. In this attempt, several small-molecule covalent inhibitors of this enzyme have been proposed. Here, we report a computational study at the atomic level of the inhibition mechanism of the most promising reported compounds. Molecular dynamics simulations were performed on six of them, and their binding energies in the active site of the protein were computed. Contact maps and interaction energies were decomposed by residues to disclose those key interactions with the enzyme. Finally, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations were performed to evaluate the reaction mechanism by which these drug candidates lead to covalently bound complexes, inhibiting the RgpB protease. The results provide a guide for future re-design of prospective and efficient inhibitors for the treatment of Alzheimer's disease.

Mechanistic Modeling of Lys745 Sulfonylation in EGFR C797S Reveals Chemical Determinants for Inhibitor Activity and Discriminates Reversible from Irreversible Agents

Targeted covalent inhibitors hold promise for drug discovery, particularly for kinases. Targeting the catalytic lysine of epidermal growth factor receptor (EGFR) has attracted attention as a new strategy to overcome resistance due to the emergence of C797S mutation. Sulfonyl fluoride derivatives able to inhibit EGFR^{L858R/T790M/C797S} by sulfonylation of Lys745 have been reported. However, atomistic details of this process are still poorly understood. Here, we describe the mechanism of inhibition of an innovative class of compounds that covalently engage the catalytic lysine of EGFR, through a sulfur(VI) fluoride exchange (SuFEx) process, with the help of hybrid quantum mechanics/molecular mechanics (QM/MM) and path collective variables (PCVs) approaches. Our simulations identify the chemical determinants accounting for the irreversible activity of agents targeting Lys745 and provide hints for the further optimization of sulfonyl fluoride agents.

Computational design of an amidase by combining the best electrostatic features of two promiscuous hydrolases

While there has been emerging interest in designing new enzymes to solve practical challenges, computer-based options to redesign catalytically active proteins are rather limited. Here, a rational QM/MM molecular dynamics strategy based on combining the best electrostatic properties of enzymes with activity in a common reaction is presented. The computational protocol has been applied to the re-design of the protein scaffold of an existing promiscuous esterase from Bacillus subtilis Bs2 to enhance its secondary amidase activity. After the alignment of Bs2 with a non-homologous amidase Candida antarctica lipase B (CALB) within rotation quaternions, a relevant spatial aspartate residue of the latter was transferred to the former as a means to favor the electrostatics of transition state formation, where a clear separation of charges takes place. Deep computational insights, however, revealed a significant conformational change caused by the amino acid replacement, provoking a shift in the pK a of the inserted aspartate and counteracting the anticipated catalytic effect. This prediction was experimentally confirmed with a 1.3-fold increase in activity. The good agreement between theoretical and experimental results, as well as the linear correlation between the electrostatic properties and the activation energy barriers, suggest that the presented computational-based investigation can transform in an enzyme engineering approach.

Impact of Warhead Modulations on the Covalent Inhibition of SARS-CoV-2 Mpro Explored by QM/MM Simulations

The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2, SARS-CoV-2, shows the need for effective antiviral treatments. Here, we present a simulation study of the inhibition of the SARS-CoV-2 main protease (Mpro), a cysteine hydrolase essential for the life cycle of the virus. The free energy landscape for the mechanism of the inhibition process is explored by QM/MM umbrella sampling and free energy perturbation simulations at the M06-2X/MM level of theory for two proposed peptidyl covalent inhibitors that share the same recognition motif but feature distinct cysteine-targeting warheads. Regardless of the intrinsic reactivity of the modeled inhibitors, namely a Michael acceptor and a hydroxymethyl ketone activated carbonyl, our results confirm that the inhibitory process takes place by means of a two-step mechanism, in which the formation of an ion pair C145/H41 dyad precedes the protein-inhibitor covalent bond formation. The nature of this second step is strongly dependent on the functional groups in the warhead: while the nucleophilic attack of the C145 sulfur atom on the Cα of the double bond of the Michael acceptor takes place concertedly with the proton transfer from H41 to Cβ, in the compound with an activated carbonyl, the sulfur attacks the carbonyl carbon concomitant with a proton transfer from H41 to the carbonyl oxygen via the hydroxyl group. An analysis of the free energy profiles, structures along the reaction path, and interactions between the inhibitors and the different pockets of the active site on the protein shows a measurable effect of the warhead on the kinetics and thermodynamics of the process. These results and QM/MM methods can be used as a guide to select warheads to design efficient irreversible and reversible inhibitors of SARS-CoV-2 Mpro.

The role of streptavidin and its variants in catalysis by biotinylated secondary amines

Here, we combine the use of host screening, protein crystallography and QM/MM molecular dynamics simulations to investigate how the protein structure affects iminium catalysis by biotinylated secondary amines in a model 1,4 conjugate addition reaction. Monomeric streptavidin (M-Sav) lacks a quaternary structure and the solvent-exposed reaction site resulted in poor product conversion in the model reaction with low enantio- and regioselectivities. These parameters were much improved when the tetrameric host T-Sav was used; indeed, residues at the symmetrical subunit interface were proven to be critical for catalysis through a mutagenesis study. The use of QM/MM simulations and the asymmetric dimeric variant D-Sav revealed that both Lys121 residues which are located in the hosting and neighboring subunits play a critical role in controlling the stereoselectivity and reactivity. Lastly, the D-Sav template, though providing a lower conversion than that of the symmetric tetrameric counterpart, is likely a better starting point for future protein engineering because each surrounding residue within the asymmetric scaffold can be refined for secondary amine catalysis.

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding

The main protease (Mpro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural Mpro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial Mpro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective 'peptibitors' inhibit Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

Unrevealing the Proteolytic Activity of RgpB Gingipain from Computational Simulations

Alzheimer's disease represents one of the greatest medical concerns for today's population and health services. Its multifactorial inherent nature represents a challenge for its treatment and requires the development of a broad spectrum of drugs. Recently, the cysteine protease gingipain RgpB has been related to neurodegenerative diseases, including Alzheimer's disease, and its inhibition appears to be a promising neuroprotective strategy. Given these features, a computational study that integrates molecular dynamics (MD) simulations with classical and hybrid quantum mechanics/molecular mechanics (QM/MM) potentials was carried out to unravel the atomistic details of RgpB activity. First, a preliminary study based on principal component analysis (PCA), determined the protonation state of the Cys/His catalytic dyad, as well as the crucial role of a flexible loop that favors reactive interactions of the catalytic residues and the peptide in the precatalytic state in its closed conformation. Then, different mechanisms were explored by means of QM/MM MD simulations. The most favorable mechanism consists of two stages. First is an acylation stage that takes place in two steps where, initially, the sulfur atom of the C244 residue attacks the carbonylic carbon of the peptide and the proton of the C244 residue is transferred to the amino group of the peptide in a concerted manner. Subsequently, the peptide bond is broken, and a fragment of the peptide is released. After that, the deacylation stage takes place in a single step where a water molecule attacks the carbonylic carbon of the peptide and a proton of the water is transferred to the C244 residue. The free energy barrier of the rate limiting step is in very good agreement with available experimental data. The mechanism exhibits an unusual role of H211 residue compared with other cysteine proteases but a crucial role of the peptide in triggering the catalysis. Notably, the atomic and energetic particularities found represent a significant contribution to the comprehension of the reaction mechanism and a great opportunity for the design of efficient inhibitors of gingipain RgpB.

Computational Studies Suggest Promiscuous Candida antarctica Lipase B as an Environmentally Friendly Alternative for the Production of Epoxides

Environmentally friendly processes are nowadays a trending topic to get highly desired chemical compounds and, in this sense, the use of enzyme-catalyzed routes is becoming a promising alternative to traditional synthetic methods. In the present paper, a hybrid quantum mechanics/molecular mechanics (QM/MM) computational study on the epoxidation of alkenes catalyzed by the Ser105Ala variant of the promiscuous Candida antarctica lipase B (CALB) is presented in an attempt to search for alternative paths to get useful intermediates in industries. The catalyzed reaction, described at the atomistic level with a model of the full solvated in a box of water molecules, is compared with the alternative epoxidation of alkenes by peroxy acids in chloroform. Free-energy profiles obtained at the density functional theory (DFT)/MM level show how Ser105Ala CALB is capable of epoxide short alkenes in a two-step process with free-energy barriers, in agreement with available experimental data, that are significantly lower than those of the single-step reaction in solution. The possible (R)-enantioselectivity dictated by the binding step, explored by means of alchemical QM/MM free-energy perturbation (FEP) methods, and the preference for the (S)-enantiomer derived from the free-energy landscape of the chemical steps would cancel out, thus predicting the lack of enantioselectivity experimentally observed. In general, our results provide general information on the molecular mechanism employed by a highly promiscuous enzyme, with potential applications in biotechnology.

Elucidating the Dual Mode of Action of Dipeptidyl Enoates in the Inhibition of Rhodesain Cysteine Proteases

A computational study of the two possible inhibition mechanisms of rhodesain cysteine protease by the dipeptidyl enoate Cbz-Phe-Leu-CH=CH-CO₂ C₂ H₅ has been carried out by means of molecular dynamics simulations with hybrid QM/MM potentials. The low free energy barriers confirm that the Cys25 residue can attack both C_β and C1 atoms of the inhibitor, confirming a dual mode of action in the inhibition of the rhodesain by enoates. According to the results, the inhibition process through the Cys25 attack on the C_β atom of the inhibitor is an exergonic and irreversible process, while the inhibition process when Cys25 attacks on the C1 atom of the inhibitor is and exergonic but reversible process. The interactions between the inhibitor and rhodesain suggest that P2 is the most important fragment to consider in the design of new efficient inhibitors of rhodesain. These results may be useful for the design of new inhibitors of rhodesain and other related cysteine proteases based on dipeptidyl enoates scaffolds.

Combined Theoretical and Experimental Study to Unravel the Differences in Promiscuous Amidase Activity of Two Nonhomologous Enzymes

Convergent evolution has resulted in nonhomologous enzymes that contain similar active sites that catalyze the same primary and secondary reactions. Comparing how these enzymes achieve their reaction promiscuity can yield valuable insights to develop functions from the optimization of latent activities. In this work, we have focused on the promiscuous amidase activity in the esterase from Bacillus subtilis (Bs2) and compared with the same activity in the promiscuous lipase B from Candida antarctica (CALB). The study, combining multiscale quantum mechanics/molecular mechanics (QM/MM) simulations, deep machine learning approaches, and experimental characterization of Bs2 kinetics, confirms the amidase activity of Bs2 and CALB. The computational results indicate that both enzymes offer a slightly different reaction environment reflected by electrostatic effects within the active site, thus resulting in a different reaction mechanism during the acylation step. A convolutional neural network (CNN) has been used to understand the conserved amino acids among the evolved protein family and suggest that Bs2 provides a more robust protein scaffold to perform future mutagenesis studies. Results derived from this work will help reveal the origin of enzyme promiscuity, which will find applications in enzyme (re)design, particularly in creating a highly active amidase.

Assessment of the PETase conformational changes induced by poly(ethylene terephthalate) binding

Recently, a bacterium strain of Ideonella sakaiensis was identified with the uncommon ability to degrade the poly(ethylene terephthalate) (PET). The PETase from I. sakaiensis strain 201-F6 (IsPETase) catalyzes the hydrolysis of PET converting it to mono(2-hydroxyethyl) terephthalic acid (MHET), bis(2-hydroxyethyl)-TPA (BHET), and terephthalic acid (TPA). Despite the potential of this enzyme for mitigation or elimination of environmental contaminants, one of the limitations of the use of IsPETase for PET degradation is the fact that it acts only at moderate temperature due to its low thermal stability. Besides, molecular details of the main interactions of PET in the active site of IsPETase remain unclear. Herein, molecular docking and molecular dynamics (MD) simulations were applied to analyze structural changes of IsPETase induced by PET binding. Results from the essential dynamics revealed that the β1-β2 connecting loop is very flexible. This loop is located far from the active site of IsPETase and we suggest that it can be considered for mutagenesis to increase the thermal stability of IsPETase. The free energy landscape (FEL) demonstrates that the main change in the transition between the unbound to the bound state is associated with the β7-α5 connecting loop, where the catalytic residue Asp206 is located. Overall, the present study provides insights into the molecular binding mechanism of PET into the IsPETase structure and a computational strategy for mapping flexible regions of this enzyme, which can be useful for the engineering of more efficient enzymes for recycling plastic polymers using biological systems.

QM/MM Study of the Enzymatic Biodegradation Mechanism of Polyethylene Terephthalate

The environmental problems derived from the generalized plastic consumption and disposal could find a friendly solution in enzymatic biodegradation. Recently, two hydrolases from Ideonella sakaiensis 201-F6 and the metagenome-derived leaf-branch compost cutinase (LCC), more specially the improved ICCG variant, have revealed degradation activity toward poly ethylene terephthalate (PET). In the present study, the reaction mechanism of this polymer breakage is studied at an atomic level by multiscale QM/MM molecular dynamics simulations, using semiempirical and DFT Hamiltonians to describe the QM region. The obtained free energy surfaces confirmed a characteristic four-step path for both systems, with activation energies in agreement with the experimental observations. Structural analysis of the evolution of the active site along the reaction progress and the study of electrostatic effects generated by the proteins reveal the similarity in the behavior of the active site of these two enzymes. The origin of the apparent better performance of the LCC-ICCG protein over PETase must be due to its capabilities of working at higher temperature and its intrinsic relationship with the crystallinity grade of the polymer. Our results may be useful for the development of more efficient enzymes in the biodegradation of PET.

A QM/MM study on the origin of retro-aldolase activity of a catalytic antibody

The retro-aldolase mechanism of methodol catalysed by the catalytic antibody 33F12 is described based on the exploration of the free energy landscape obtained with QM/MM methods. The amino acids involved in the reaction have been identified, as well as their specific role played in the active site and in the flexibility of the loops. Finally, the comparison with a de novo enzyme RA95.5-8F provides a deeper understanding of catalytic differences between such different protein scaffolds.

Caught in Action: X-ray Structure of Thymidylate Synthase with Noncovalent Intermediate Analog

Methylation of 2-deoxyuridine-5'-monophosphate (dUMP) at the C5 position by the obligate dimeric thymidylate synthase (TSase) in the sole de novo biosynthetic pathway to thymidine 5'-monophosphate (dTMP) proceeds by forming a covalent ternary complex with dUMP and cosubstrate 5,10-methylenetetrahydrofolate. The crystal structure of an analog of this intermediate gives important mechanistic insights but does not explain the half-of-the-sites activity of the enzyme. Recent experiments showed that the C5 proton and the catalytic Cys are eliminated in a concerted manner from the covalent ternary complex to produce a noncovalent bisubstrate intermediate. Here, we report the crystal structure of TSase with a close synthetic analog of this intermediate in which it has partially reacted with the enzyme but in only one protomer, consistent with the half-of-the-sites activity of this enzyme. Quantum mechanics/molecular mechanics simulations confirmed that the analog could undergo catalysis. The crystal structure shows a new water 2.9 Å from the critical C5 of the dUMP moiety, which in conjunction with other residues in the network, may be the elusive general base that abstracts the C5 proton of dUMP during the reaction.

Transfer hydrogenations catalyzed by streptavidin-hosted secondary amine organocatalysts

Here, the streptavidin-biotin technology was applied to enable organocatalytic transfer hydrogenation. By introducing a biotin-tethered pyrrolidine (1) to the tetrameric streptavidin (T-Sav), the resulting hybrid catalyst was able to mediate hydride transfer from dihydro-benzylnicotinamide (BNAH) to α,β-unsaturated aldehydes. Hydrogenation of cinnamaldehyde and some of its aryl-substituted analogues was found to be nearly quantitative. Kinetic measurements revealed that the T-Sav:1 assembly possesses enzyme-like behavior, whereas isotope effect analysis, performed by QM/MM simulations, illustrated that the step of hydride transfer is at least partially rate-limiting. These results have proven the concept that T-Sav can be used to host secondary amine-catalyzed transfer hydrogenations.

Mechanism of inhibition of SARS-CoV-2 Mpro by peptidyl Michael acceptor explained by QM/MM simulations and design of new derivatives with tunable chemical reactivity

The SARS-CoV-2 main protease (M^pro) is essential for replication of the virus responsible for the COVID-19 pandemic, and one of the main targets for drug design. Here, we simulate the inhibition process of SARS-CoV-2 M^pro with a known Michael acceptor (peptidyl) inhibitor, N3. The free energy landscape for the mechanism of the formation of the covalent enzyme-inhibitor product is computed with QM/MM molecular dynamics methods. The simulations show a two-step mechanism, and give structures and calculated barriers in good agreement with experiment. Using these results and information from our previous investigation on the proteolysis reaction of SARS-CoV-2 M^pro, we design two new, synthetically accessible N3-analogues as potential inhibitors, in which the recognition and warhead motifs are modified. QM/MM modelling of the mechanism of inhibition of M^pro by these novel compounds indicates that both may be promising candidates as drug leads against COVID-19, one as an irreversible inhibitor and one as a potential reversible inhibitor.

Selective oxidation of alkyl and aryl glyceryl monoethers catalysed by an engineered and immobilised glycerol dehydrogenase

Enzymes acting over glyceryl ethers are scarce in living cells, and consequently biocatalytic transformations of these molecules are rare despite their interest for industrial chemistry. In this work, we have engineered and immobilised a glycerol dehydrogenase from Bacillus stearothermophilus (BsGlyDH) to accept a battery of alkyl/aryl glyceryl monoethers and catalyse their enantioselective oxidation to yield the corresponding 3-alkoxy/aryloxy-1-hydroxyacetones. QM/MM computational studies decipher the key role of D123 in the oxidation catalytic mechanism, and reveal that this enzyme is highly enantioselective towards S-isomers (ee > 99%). Through structure-guided site-selective mutagenesis, we find that the mutation L252A sculpts the active site to accommodate a productive configuration of 3-monoalkyl glycerols. This mutation enhances the k cat 163-fold towards 3-ethoxypropan-1,2-diol, resulting in a specific activity similar to the one found for the wild-type towards glycerol. Furthermore, we immobilised the L252A variant to intensify the process, demonstrating the reusability and increasing the operational stability of the resulting heterogeneous biocatalyst. Finally, we manage to integrate this immobilised enzyme into a one-pot chemoenzymatic process to convert glycidol and ethanol into 3-ethoxy-1-hydroxyacetone and (R)-3-ethoxypropan-1,2-diol, without affecting the oxidation activity. These results thus expand the uses of engineered glycerol dehydrogenases in applied biocatalysis for the kinetic resolution of glycerol ethers and the manufacturing of substituted hydroxyacetones.

Glycoside hydrolase stabilization of transition state charge: new directions for inhibitor design

Carbasugars are structural mimics of naturally occurring carbohydrates that can interact with and inhibit enzymes involved in carbohydrate processing. In particular, carbasugars have attracted attention as inhibitors of glycoside hydrolases (GHs) and as therapeutic leads in several disease areas. However, it is unclear how the carbasugars are recognized and processed by GHs. Here, we report the synthesis of three carbasugar isotopologues and provide a detailed transition state (TS) analysis for the formation of the initial GH-carbasugar covalent intermediate, as well as for hydrolysis of this intermediate, using a combination of experimentally measured kinetic isotope effects and hybrid QM/MM calculations. We find that the α-galactosidase from Thermotoga maritima effectively stabilizes TS charge development on a remote C5-allylic center acting in concert with the reacting carbasugar, and catalysis proceeds via an exploded, or loose, S_N2 transition state with no discrete enzyme-bound cationic intermediate. We conclude that, in complement to what we know about the TS structures of enzyme-natural substrate complexes, knowledge of the TS structures of enzymes reacting with non-natural carbasugar substrates shows that GHs can stabilize a wider range of positively charged TS structures than previously thought. Furthermore, this enhanced understanding will enable the design of new carbasugar GH transition state analogues to be used as, for example, chemical biology tools and pharmaceutical lead compounds.

Positive charge stabilized on remote C5-allylic center with catalysis occurring via a loose S_N2 transition state.

Electric Field Measurements Reveal the Pivotal Role of Cofactor-Substrate Interaction in Dihydrofolate Reductase Catalysis

The contribution of ligand–ligand electrostatic interaction to transition state formation during enzyme catalysis has remained unexplored, even though electrostatic forces are known to play a major role in protein functions and have been investigated by the vibrational Stark effect (VSE). To monitor electrostatic changes along important steps during catalysis, we used a nitrile probe (T46C-CN) inserted proximal to the reaction center of three dihydrofolate reductases (DHFRs) with different biophysical properties, Escherichia coli DHFR (EcDHFR), its conformationally impaired variant (EcDHFR-S148P), and Geobacillus stearothermophilus DHFR (BsDHFR). Our combined experimental and computational approach revealed that the electric field projected by the substrate toward the probe negates those exerted by the cofactor when both are bound within the enzymes. This indicates that compared to previous models that focus exclusively on subdomain reorganization and protein–ligand contacts, ligand–ligand interactions are the key driving force to generate electrostatic environments conducive for catalysis.

Revealing the molecular mechanisms of proteolysis of SARS-CoV-2 Mpro by QM/MM computational methods

SARS-CoV-2 M^pro is one of the enzymes essential for the replication process of the virus responsible for the COVID-19 pandemic. This work is focused on exploring its proteolysis reaction by means of QM/MM methods. The resulting free energy landscape of the process provides valuable information on the species appearing along the reaction path and suggests that the mechanism of action of this enzyme, taking place in four steps, slightly differs from that of other cysteine proteases. Our predictions, which are in agreement with some recently published experimental data, can be used to guide the design of COVID-19 antiviral compounds with clinical potential.

The molecular mechanism of the proteolysis reaction catalyzed by SARS-CoV-2 M^pro, one of the enzymes essential for the replication process of the virus responsible for the COVID-19 pandemic, is described using computational QM/MM methods.

Quantum Mechanics/Molecular Mechanics Studies of the Mechanism of Cysteine Proteases Inhibition by Dipeptidyl Nitroalkenes

In this work a computational study of the mechanism of inhibition of cruzain, rhodesain, and cathepsin L cysteine proteases by the dipeptidyl nitroalkene Cbz-Phe-Ala-CH=CH-NO₂ has been carried out by means of molecular dynamics simulations with hybrid QM/MM potentials. The free-energy surfaces confirmed that the inhibition takes place by the formation of a covalent bond between the protein and the β-carbon atom of the inhibitor. According to the results, the tested inhibitor should be a much more efficient inhibitor of cruzain than of rhodesain, and little activity would be expected against cathepsin L, in total correspondence with the available experimental data. The origin of these differences may lie in the different stabilizing electrostatic interactions established between the inhibitor and the residues of the active site and S2 pocket of these enzymes. These results may be useful for the rational design of new dipeptidyl nitroalkenes with higher and more selective inhibitory activity against cysteine proteases.

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.

Loss of Hyperconjugative Effects Drives Hydride Transfer during Dihydrofolate Reductase Catalysis

Hydride transfer is widespread in nature and has an essential role in applied research. However, the mechanisms of how this transformation occurs in living organisms remain a matter of vigorous debate. Here, we examined dihydrofolate reductase (DHFR), an enzyme that catalyzes hydride from C4′ of NADPH to C6 of 7,8-dihydrofolate (H₂F). Despite many investigations of the mechanism of this reaction, the contribution of polarization of the π-bond of H₂F in driving hydride transfer remains unclear. H₂F was stereospecifically labeled with deuterium β to the reacting center, and β-deuterium kinetic isotope effects were measured. Our experimental results combined with analysis derived from QM/MM simulations reveal that hydride transfer is triggered by polarization at the C6 of H₂F. The σ C_β–H bonds contribute to the buildup of the cationic character during the chemical transformation, and hyperconjugation influences the formation of the transition state. Our findings provide key insights into the hydride transfer mechanism of the DHFR-catalyzed reaction, which is a target for antiproliferative drugs and a paradigmatic model in mechanistic enzymology.

Why are some Enzymes Dimers? Flexibility and Catalysis in Thermotoga Maritima Dihydrofolate Reductase

Dihydrofolate reductase from Thermotoga maritima (TmDFHFR) is a dimeric thermophilic enzyme that catalyzes the hydride transfer from the cofactor NADPH to dihydrofolate less efficiently than other DHFR enzymes, such as the mesophilic analogue Escherichia coli DHFR (EcDHFR). Using QM/MM potentials, we show that the reduced catalytic efficiency of TmDHFR is most likely due to differences in the amino acid sequence that stabilize the M20 loop in an open conformation, which prevents the formation of some interactions in the transition state and increases the number of water molecules in the active site. However, dimerization provides two advantages to the thermophilic enzyme: it protects its structure against denaturation by reducing thermal fluctuations and it provides a less negative activation entropy, toning down the increase of the activation free energy with temperature. Our molecular picture is confirmed by the analysis of the temperature dependence of enzyme kinetic isotope effects in different DHFR enzymes.

Convergence of theory and experiment on the role of preorganization, quantum tunneling and enzyme motions into flavoenzyme-catalyzed hydride transfer

Hydride transfer is one of the most common reactions catalyzed by enzymatic systems and it has become an object of study due to possible significant quantum tunneling effects. In the present work, we provide a combination of theoretical QM/MM simulations and experimental measurements of the rate constants and kinetic isotopic effects (KIEs) for the hydride transfer reaction catalyzed by morphinone reductase, MR. Quantum mechanical tunneling coefficients, computed in the framework of variational transition-state theory, play a significant role in this reaction, reaching values of 23.8 ± 5.5 for the lightest isotopologue; one of the largest values reported for enzymatic systems. This prediction is supported by the agreement between the theoretically predicted rate constants and the corresponding experimental values. Simulations indicate that the role of protein motions can be satisfactorily described as equilibrium fluctuations along the reaction coordinate, in line with a high degree of preorganization displayed by this enzyme.

Experimental and Computational Studies Delineate the Role of Asparagine 177 in Hydride Transfer for E. coli Thymidylate Synthase

Thymidylate synthase (TSase), an enzyme responsible for the de novo biosynthesis of 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) necessary for DNA synthesis, has been a drug target for decades. TSase is a highly conserved enzyme across species ranging from very primitive organisms to mammals. Among the many conserved active site residues, an asparagine (N177, using Escherichia coli residues numbering) appears to make direct hydrogen bonds with both the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate and its pyrimidine ring's N3. Recent studies have reassessed the TSase catalytic mechanism, focusing on the degree of negative charge accumulation at the O4 carbonyl of the substrate during two critical H-transfers - a proton abstraction and a hydride transfer. To obtain insights into the role of this conserved N177 on the hydride transfer, we examined its aspartic acid (D) and serine (S) mutants - each of which is expected to alter hydrogen bonding and charge stabilization around the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate. Steady-state kinetics, substrate binding order studies and temperature-dependency analysis of intrinsic KIEs for the hydride transfer step of the TSase catalytic cycle suggest the active site of N177D is not precisely organized for that step. A smaller disruption was observed for N177S, which could be rationalized by partial compensation by water molecules and rearrangement of other residues toward preparation of the system for the hydride transfer under study. These experimental findings are qualitatively mirrored by QM/MM computational simulations, thereby shedding light on the sequence and synchronicity of steps in the TSase-catalyzed reaction. This information could potentially inform the design of mechanism-based drugs targeting this enzyme.

Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase

The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed light on this debate, has been computed and the results have been compared with previous experimental data. The classical mechanical free energy barrier and the number of recrossing trajectories seem to be temperature-independent while the quantum vibrational corrections and the tunnelling effects are slightly temperature-dependent over the interval of 5-45 °C. The computed primary KIEs are in very good agreement with previous experimental data, being almost temperature-independent within the standard deviations. The modest dependence on the temperature is due to just the quantum vibrational correction contribution. These results, together with the analysis of the evolution of the collective variables such as the electrostatic potential or the electric field created by the protein on the key atoms involved in the reaction, confirm that while the protein is well preorganised, some changes take place along the reaction that favour the hydride transfer and the product release. Coordinates defining these movements are, in fact, part of the real reaction coordinate.

Reactivity and Selectivity of Iminium Organocatalysis Improved by a Protein Host

There has been growing interest in performing organocatalysis within a supramolecular system as a means of controlling reaction reactivity and stereoselectivity. Here, a protein is used as a host for iminium catalysis. A pyrrolidine moiety is covalently linked to biotin and introduced to the protein host streptavidin for organocatalytic activity. Whereas in traditional systems stereoselectivity is largely controlled by the substituents added to the organocatalyst, enantiomeric enrichment by the reported supramolecular system is completely controlled by the host. Also, the yield of the model reaction increases over 10-fold when streptavidin is included. A 1.1 Å crystal structure of the protein-catalyst complex and molecular simulations of a key intermediate reveal the chiral scaffold surrounding the organocatalytic reaction site. This work illustrates that proteins can be an excellent supramolecular host for driving stereoselective secondary amine organocatalysis.