Exploring Chloride Selectivity and Halogenase Regioselectivity of the SalL Enzyme through Quantum Mechanical/Molecular Mechanical Modeling

Molecular Insights into Converting Hydroxide Adenosyltransferase into Halogenase

Halogenation plays a unique role in the design of agrochemicals. Enzymatic halogenation reactions have attracted great attention due to their excellent specificity and mild reaction conditions. S-adenosyl-l-methionine (SAM)-dependent halogenases mediate the nucleophilic attack of halide ions (X-) to SAM to produce 5'-XDA. However, only 11 SAM-dependent fluorinases and 3 chlorinases have been reported, highlighting the desire for additional halogenases. SAM-dependent hydroxide adenosyltransferase (HATase) has a similar reaction mechanism as halogenases but uses water as a substrate instead of halide ions. Here, we explored a HATase from the thermophile Thermotoga maritima MSB8 and transformed it into a halogenase. We identified a key dyad W8L/V71T for the halogenation reaction. We also obtained the best performing mutants for each halogenation reaction: M1, M2 and M4 for Cl-, Br- and I-, respectively. The M4 mutant retained the thermostability of HATase in the iodination reaction at 80 °C, which surpasses the natural halogenase SalL. QM/MM revealed that these mutants bind halide ions with more suitable angles for nucleophilic attack of C5' of SAM, thus conferring halogenation capabilities. Our work achieved the halide ion specificity of halogenases and generated thermostable halogenases for the first time, which provides new opportunities to expand the halogenase repertoire from hydroxylase.

Computational Analysis of SAM Analogs as Methyltransferase Inhibitors of nsp16/nsp10 Complex from SARS-CoV-2

Methyltransferases (MTases) enzymes, responsible for RNA capping into severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), are emerging important targets for the design of new anti-SARS-CoV-2 agents. Here, analogs of S-adenosylmethionine (SAM), obtained from the bioisosteric substitution of the sulfonium and amino acid groups, were evaluated by rigorous computational modeling techniques such as molecular dynamics (MD) simulations followed by relative binding free analysis against nsp16/nsp10 complex from SARS-CoV-2. The most potent inhibitor (2a) shows the lowest binding free energy (-58.75 Kcal/mol) and more potency than Sinefungin (SFG) (-39.8 Kcal/mol), a pan-MTase inhibitor, which agrees with experimental observations. Besides, our results suggest that the total binding free energy of each evaluated SAM analog is driven by van der Waals interactions which can explain their poor cell permeability, as observed in experimental essays. Overall, we provide a structural and energetic analysis for the inhibition of the nsp16/nsp10 complex involving the evaluated SAM analogs as potential inhibitors.

A Nonconventional Archaeal Fluorinase Identified by In Silico Mining for Enhanced Fluorine Biocatalysis

Fluorinases, the only enzymes known to catalyze the transfer of fluorine to an organic molecule, are essential catalysts for the biological synthesis of valuable organofluorines. However, the few fluorinases identified so far have low turnover rates that hamper biotechnological applications. Here, we isolated and characterized putative fluorinases retrieved from systematic in silico mining and identified a nonconventional archaeal enzyme from Methanosaeta sp. that mediates the fastest S_N2 fluorination rate reported to date. Furthermore, we demonstrate enhanced production of fluoronucleotides in vivo in a bacterial host engineered with this archaeal fluorinase, paving the way toward synthetic metabolism for efficient biohalogenation.

A Computational Study of the Promiscuity of the SAM-Dependent Methyltransferase AtHTMT1

A two-pronged computational approach was taken to study the promiscuity of the SAM⁺-dependent methyl transferase AtHTMT1 from thale cress with several nucleophiles (Cl^-, Br^-, I^-, NCO^-, NCS^-). First, enzyme-free methyl transfer reactions were studied with M05/6-311+G(2d,p) DFT calculations and electrostatic continuum models (PCM/SMD) for various chemical environments. Second, QM/MM MD simulations with semiempirical Hamiltonians (PM7, PM6-D3, AM1, PM6-D3H4) and the AMBER 14SB force field were used to study the enzyme catalyzed reaction in silico. The combination of the DFT and MD results shows that reactant desolvation generally accelerates the reaction, but it cannot explain the selectivity of the enzyme. The critical position of H₂O molecules at the reactive site favors the reaction of NCS^- over Cl^- and Br^- in agreement with experiments, but not observed in the quantum calculations for the cytosol. The addition of selected H₂O molecules to the N terminus of NCS^- greatly increases its reactivity, while H₂O molecules attached to Cl^- slow the reaction. The partial solvation of the nucleophiles in the reactive pouch holds the key to understanding the reactivity of AtHTMT1.

Exploring the Catalytic Mechanism of the RNA Cap Modification by nsp16-nsp10 Complex of SARS-CoV-2 through a QM/MM Approach

The inhibition of key enzymes that may contain the viral replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have assumed central importance in drug discovery projects. Nonstructural proteins (nsps) are essential for RNA capping and coronavirus replication since it protects the virus from host innate immune restriction. In particular, nonstructural protein 16 (nsp16) in complex with nsp10 is a Cap-0 binding enzyme. The heterodimer formed by nsp16-nsp10 methylates the 5′-end of virally encoded mRNAs to mimic cellular mRNAs and thus it is one of the enzymes that is a potential target for antiviral therapy. In this study, we have evaluated the mechanism of the 2′-O methylation of the viral mRNA cap using hybrid quantum mechanics/molecular mechanics (QM/MM) approach. It was found that the calculated free energy barriers obtained at M062X/6-31+G(d,p) is in agreement with experimental observations. Overall, we provide a detailed molecular analysis of the catalytic mechanism involving the 2′-O methylation of the viral mRNA cap and, as expected, the results demonstrate that the TS stabilization is critical for the catalysis.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Drug repurposing and computational modeling for discovery of inhibitors of the main protease (Mpro) of SARS-CoV-2

The main protease (M^pro or 3CL^pro) is a conserved cysteine protease from the coronaviruses and started to be considered an important drug target for developing antivirals, as it produced a deadly outbreak of COVID-19. Herein, we used a combination of drug reposition and computational modeling approaches including molecular docking, molecular dynamics (MD) simulations, and the calculated binding free energy to evaluate a set of drugs in complex with the M^pro enzyme. Particularly, our results show that darunavir and triptorelin drugs have favorable binding free energy (−63.70 and −77.28 kcal mol⁻¹, respectively) in complex with the M^pro enzyme. Based on the results, the structural and energetic features that explain why some drugs can be repositioned to inhibit M^pro from SARS-CoV-2 were exposed. These features should be considered for the design of novel M^pro inhibitors.

Structural and energetic features explain why some drugs can be repositioned to inhibit M^pro from SARS-CoV-2.

QM/MM Study of the Fosfomycin Resistance Mechanism Involving FosB Enzyme

Multidrug-resistant organisms contain antibiotic-modifying enzymes that facilitate resistance to a variety of antimicrobial compounds. Particularly, the fosfomycin (FOF) drug can be structurally modified by several FOF-modifying enzymes before it reaches the biological target. Among them, FosB is an enzyme that utilizes l-cysteine or bacillithiol in the presence of a divalent metal to open the epoxide ring of FOF and, consequently, inactivate the drug. Here, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations to explore the mechanism of the reaction involving FosB and FOF. The calculated free-energy profiles show that the cost to open the epoxide ring of FOF at the C2 atom is ∼3.0 kcal/mol higher than that at the C1 atom. Besides, our QM/MM MD results revealed the critical role of conformation change of Cys9 and Asn50 to release the drug from the active site. Overall, the present study provides insights into the mechanism of FOF-resistant proteins.

Synthetic biology enabling access to designer polyketides

The full potential of polyketide discovery has yet to be reached owing to a lack of suitable technologies and knowledge required to advance engineering of polyketide biosynthesis. Recent investigations on the discovery, enhancement, and non-natural use of these biosynthetic gene clusters via computational biology, metabolic engineering, structural biology, and enzymology-guided approaches have facilitated improved access to designer polyketides. Here, we discuss recent successes in gene cluster discovery, host strain engineering, precursor-directed biosynthesis, combinatorial biosynthesis, polyketide tailoring, and high-throughput synthetic biology, as well as challenges and outlooks for rapidly generating useful target polyketides.

Evaluating the Performance of a Non-Bonded Cu2+ Model Including Jahn-Teller Effect into the Binding of Tyrosinase Inhibitors

Tyrosinase (TYR) is a metalloenzyme classified as a type-3 copper protein, which is involved in the synthesis of melanin through a catalytic process beginning with the conversion of the amino acid l-Tyrosine (l-Tyr) to l-3,4-dihydroxyphenylalanine (l-DOPA). It plays an important role in the mechanism of melanogenesis in various organisms including mammals, plants, and fungi. Herein, we used a combination of computational molecular modeling techniques including molecular dynamic (MD) simulations and the linear interaction energy (LIE) model to evaluate the binding free energy of a set of analogs of kojic acid (KA) in complex with TYR. For the MD simulations, we used a dummy model including the description of the Jahn–Teller effect for Cu²⁺ ions in the active site of this enzyme. Our results show that the LIE model predicts the TYR binding affinities of the inhibitor in close agreement to experimental results. Overall, we demonstrate that the classical model provides a suitable description of the main interactions between analogs of KA and Cu²⁺ ions in the active site of TYR.

One Pot Use of Combilipases for Full Modification of Oils and Fats: Multifunctional and Heterogeneous Substrates

Lipases are among the most utilized enzymes in biocatalysis. In many instances, the main reason for their use is their high specificity or selectivity. However, when full modification of a multifunctional and heterogeneous substrate is pursued, enzyme selectivity and specificity become a problem. This is the case of hydrolysis of oils and fats to produce free fatty acids or their alcoholysis to produce biodiesel, which can be considered cascade reactions. In these cases, to the original heterogeneity of the substrate, the presence of intermediate products, such as diglycerides or monoglycerides, can be an additional drawback. Using these heterogeneous substrates, enzyme specificity can promote that some substrates (initial substrates or intermediate products) may not be recognized as such (in the worst case scenario they may be acting as inhibitors) by the enzyme, causing yields and reaction rates to drop. To solve this situation, a mixture of lipases with different specificity, selectivity and differently affected by the reaction conditions can offer much better results than the use of a single lipase exhibiting a very high initial activity or even the best global reaction course. This mixture of lipases from different sources has been called “combilipases” and is becoming increasingly popular. They include the use of liquid lipase formulations or immobilized lipases. In some instances, the lipases have been coimmobilized. Some discussion is offered regarding the problems that this coimmobilization may give rise to, and some strategies to solve some of these problems are proposed. The use of combilipases in the future may be extended to other processes and enzymes.