Computational study of the mechanism of a polyurethane esterase A (PueA) from Pseudomonas chlororaphis.

The effective management of plastic waste has become a global imperative, given our reliance on a linear model in which plastics are manufactured, used once, and then discarded. This has led to the pervasive accumulation of plastic debris in landfills and environmental contamination. Recognizing this issue, numerous initiatives are underway to address the environmental repercussions associated with plastic disposal. In this study, we investigate the possible molecular mechanism of polyurethane esterase A (PueA), which has been previously identified as responsible for the degradation of a polyester polyurethane (PU) sample in Pseudomonas chlororaphis, as an effort to develop enzymatic biodegradation solutions. After generating the unsolved 3D structure of the protein by AlphaFold2 from its known genome, the enzymatic hydrolysis of the same model PU compound previously used in experiments has been explored employing QM/MM molecular dynamics simulations. This required a preliminary analysis of the 3D structure of the apo-enzyme, identifying the putative active site, and the search for the optimal protein-substrate binding site. Finally, the resulting free energy landscape indicates that wild-type PueA can degrade PU chains, although with low-level activity. The reaction takes place by a characteristic four-step path of the serine hydrolases, involving an acylation followed by a diacylation step. Energetics and structural analysis of the evolution of the active site along the reaction suggests that PueA can be considered a promising protein scaffold for further development to achieve efficient biodegradation of PU.

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.

Unlocking a (Pseudo)-Mechanically Interlocked Molecule with a Coronene "Shoehorn".

Mechanically interlocked molecules (MIMs) have gained increasing interest during the last decades, not only because of their aesthetic appeal, but also because their unique properties have allowed them to find applications in nanotechnology, catalysis, chemosensing and biomedicine. Herein we describe how a pyrene molecule with four octynyl substituents can be easily encapsulated within the cavity of a tetragold(I) rectangle-like metallobox, by template formation of the metallo-assembly in the presence of the guest. The resulting assembly behaves as a mechanically interlocked molecule (MIM), in which the four long limbs of the guest protrude from the entrances of the metallobox, thus locking the guest inside the cavity of the metallobox. The new assembly resembles a metallo-suit[4]ane, given the number of protruding long limbs and the presence of the metal atoms in the host molecule. However, unlike normal MIMs, this molecule can release the tetra-substituted pyrene guest by the addition of coronene, which can smoothly replace the guest in the cavity of the metallobox. Combined experimental and computational studies allowed the role of the coronene molecule in facilitating the release of the tetrasubstituted pyrene guest to be explained, through a process that we named "shoehorning", as the coronene compresses the flexible limbs of the guest so that it can reduce its size to slide in and out the metallobox.

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.

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.

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.

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.

Parallel reaction pathways and noncovalent intermediates in thymidylate synthase revealed by experimental and computational tools.

Thymidylate synthase was one of the most studied enzymes due to its critical role in molecular pathogenesis of cancer. Nevertheless, many atomistic details of its chemical mechanism remain unknown or debated, thereby imposing limits on design of novel mechanism-based anticancer therapeutics. Here, we report unprecedented isolation and characterization of a previously proposed intact noncovalent bisubstrate intermediate formed in the reaction catalyzed by thymidylate synthase. Free-energy surfaces of the bisubstrate intermediates interconversions computed with quantum mechanics/molecular mechanics (QM/MM) methods and experimental assessment of the corresponding kinetics indicate that the species is the most abundant productive intermediate along the reaction coordinate, whereas accumulation of the covalent bisubstrate species largely occurs in a parallel nonproductive pathway. Our findings not only substantiate relevance of the previously proposed noncovalent intermediate but also support potential implications of the overstabilized covalent intermediate in drug design targeting DNA biosynthesis.

Influence of Dielectric Environment upon Isotope Effects on Glycoside Heterolysis: Computational Evaluation and Atomic Hessian Analysis.

Isotope effects depend upon the polarity of the bulk medium in which a chemical process occurs. Implicit solvent calculations with molecule-shaped cavities show that the equilibrium isotope effect (EIE) for heterolysis of the glycosidic bonds in 5'-methylthioadenosine and in 2-(p-nitrophenoxy)tetrahydropyran, both in water, are very sensitive in the range 2 ≤ ε ≤ 10 to the relative permittivity of the continuum surrounding the oxacarbenium ion. However, different implementations of nominally the same PCM method can lead to opposite trends being predicted for the same molecule. Computational modeling of the influence of the inhomogeneous effective dielectric surrounding a substrate within the protein environment of an enzymic reaction requires an explicit treatment. The EIE (KH/KD) for transfer of cyclopentyl, cyclohexyl, tetrahydrofuranyl and tetrahydropyranyl cations from water to cyclohexane is predicted by B3LYP/6-31+G(d) calculations with implicit solvation and confirmed by B3LYP/6-31+G(d)/OPLS-AA calculations with averaging over many explicit solvation configurations. Atomic Hessian analysis, whereby the full Hessian is reduced to the elements belonging to a single atom at the site of isotopic substitution, reveals a remarkable result for both implicit and explicit solvation: the influence of the solvent environment on these EIEs is essentially captured completely by only a 3 × 3 block of the Hessian, although these values must correctly reflect the influence of the whole environment. QM/MM simulation with ensemble averaging has an important role to play in assisting the meaningful interpretation of observed isotope effects for chemical reactions both in solution and catalyzed by enzymes.

Molecular Insights into the Substrate-Assisted Mechanism of Viral DNA 3'-End Processing in Intasome of Prototype Foamy Virus Integrase from Molecular Dynamic and QM/MM Studies.

Integrases participate in two important steps for virus replication such as 3'-end processing of viral DNA ( vDNA) and nuclear entry of host DNA ( hDNA). In this work, insight into the structural changes in intasome of prototype foamy virus integrase (PFV-IN) complexed with vDNA from classical molecular dynamic (MD) simulations are done. Analysis of the results reveal the existence of alternative conformations of the enzyme active site indicating that the 3'-end processing reaction can occur according to three different pathways and taking place with the possible participation of aspartate 185 of a neighboring phosphate group or involving an internal phosphate group of the substrate. In this work, one of them, the so-called substrate-assisted mechanism was explored by QM/MM methods. The free energy barriers of 34.4 kcal mol-1 for the first and 35.3 kcal mol-1 for the second step of reaction computed with free energy perturbation (FEP) methods at the M06-2X/AMBER level show that 3'-end processing has to proceed via a different mechanism than studied herein. Nevertheless, the obtained results are in good agreement with the experimental observations that the substitution of the key atom for this mechanism, oxygen by sulfur, did not influence the catalysis. Additionally, the obtained mechanism reveals significant similarities to the previously studied substrate-assisted mechanism in twister ribozyme. The possible role of Mg2+ in the active site is discussed.

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.

Reaction Mechanism of Organocatalytic Michael Addition of Nitromethane to Cinnamaldehyde: A Case Study on Catalyst Regeneration and Solvent Effects.

The Michael addition of nitromethane to cinnamaldehyde has been computationally studied in the absence of a catalyst and the presence of a biotinylated secondary amine by a combined computational and experimental approach. The calculations were performed at the density functional theory (DFT) level with the M06-2X hybrid functional, and a polarizable continuum model has been employed to mimic the effect of two different solvents: dichloromethane (DCM) and water. Contrary to common assumption, the product-derived iminium intermediate was absent in both of the solvents tested. Instead, hydrating the C1-C2 double bond in the enamine intermediate directly yields the tetrahedral intermediate, which is key for forming the product and regenerating the catalyst. Enamine hydration is concerted and found to be rate-limiting in DCM but segregated into two non-rate-limiting steps when the solvent is replaced with water. However, further analysis revealed that the use of water as solvent also raises the energy barriers for other chemical steps, particularly the critical step of C-C bond formation between the iminium intermediate and nucleophile; this consequently lowers both the reaction yield and enantioselectivity of this LUMO-lowering reaction, as experimentally detected. These findings provide a logical explanation to why water often enhances organocatalysis when used as an additive but hampers the reaction progress when employed as a solvent.

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.

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.

Binding Analysis of Some Classical Acetylcholinesterase Inhibitors: Insights for a Rational Design Using Free Energy Perturbation Method Calculations with QM/MM MD Simulations.

In the present study, the binding free energy of some classical inhibitors (DMT, DNP, GNT, HUP, THA) with acetylcholinesterase (AChE) is calculated by means of the free energy perturbation (FEP) method based on hybrid quantum mechanics and molecular mechanics (QM/MM) potentials. The results highlight the key role of the van der Waals interaction for the inhibition process, since the contribution of this term to the binding free energy is almost as decisive as the electrostatic one. The analysis of the geometrical parameters and the interaction energy per residue along the QM/MM molecular dynamics (MD) simulations highlights the most relevant interactions in the different AChE-ligand systems, showing that the charged residues with a more prominent contribution to the interaction energy are Asp72 and Glu199, although the relative importance depends on the molecular size of the ligand. A correlation between the binding free energy and the number of cation-π interactions present in the systems has been established, DMT being the most potent inhibitor, capable of forming four cation-π interactions. A layer of water molecules surrounding the inhibitors has been observed, which act as bridges along a network formed by the ligands and the residues of the gorge and also between different residues. Although several hydrogen bonds between ligands and AChE do appear, no significant values of BIEs have been recorded. This behavior can be accounted for by the special features of AChE, such as the presence of several subsites of different natures in the gorge or the existence of several water molecules that act as bridges in the electrostatic interactions.

Dynamic and Electrostatic Effects on the Reaction Catalyzed by HIV-1 Protease.

HIV-1 Protease (HIV-1 PR) is one of the three enzymes essential for the replication process of HIV-1 virus, which explains why it has been the main target for design of drugs against acquired immunodeficiency syndrome (AIDS). This work is focused on exploring the proteolysis reaction catalyzed by HIV-1 PR, with special attention to the dynamic and electrostatic effects governing its catalytic power. Free energy surfaces for all possible mechanisms have been computed in terms of potentials of mean force (PMFs) within hybrid QM/MM potentials, with the QM subset of atoms described at semiempirical (AM1) and DFT (M06-2X) level. The results suggest that the most favorable reaction mechanism involves formation of a gem-diol intermediate, whose decomposition into the product complex would correspond to the rate-limiting step. The agreement between the activation free energy of this step with experimental data, as well as kinetic isotope effects (KIEs), supports this prediction. The role of the protein dynamic was studied by protein isotope labeling in the framework of the Variational Transition State Theory. The predicted enzyme KIEs, also very close to the values measured experimentally, reveal a measurable but small dynamic effect. Our calculations show how the contribution of dynamic effects to the effective activation free energy appears to be below 1 kcal·mol-1. On the contrary, the electric field created by the protein in the active site of the enzyme emerges as being critical for the electronic reorganization required during the reaction. These electrostatic properties of the active site could be used as a mold for future drug design.

Peptide Bond Formation Mechanism Catalyzed by Ribosome.

In this paper we present a study of the peptide bond formation reaction catalyzed by ribosome. Different mechanistic proposals have been explored by means of Free Energy Perturbation methods within hybrid QM/MM potentials, where the chemical system has been described by the M06-2X functional and the environment by means of the AMBER force field. According to our results, the most favorable mechanism in the ribosome would proceed through an eight-membered ring transition state, involving a proton shuttle mechanism through the hydroxyl group of the sugar and a water molecule. This transition state is similar to that described for the reaction in solution (J. Am. Chem. Soc. 2013, 135, 8708-8719), but the reaction mechanisms are noticeably different. Our simulations reproduce the experimentally determined catalytic effect of ribosome that can be explained by the different behavior of the two environments. While the solvent reorganizes during the chemical process involving an entropic penalty, the ribosome is preorganized in the formation of the Michaelis complex and does not suffer important changes along the reaction, dampening the charge redistribution of the chemical system.

The influence of active site conformations on the hydride transfer step of the thymidylate synthase reaction mechanism.

The hydride transfer from C6 of tetrahydrofolate to the reaction's exocyclic methylene-dUMP intermediate is the rate limiting step in thymidylate synthase (TSase) catalysis. This step has been studied by means of QM/MM molecular dynamics simulations to generate the corresponding free energy surfaces. The use of two different initial X-ray structures has allowed exploring different conformational spaces and the existence of chemical paths with not only different reactivities but also different reaction mechanisms. The results confirm that this chemical conversion takes place preferentially via a concerted mechanism where the hydride transfer is conjugated to thiol-elimination from the product. The findings also confirm the labile character of the substrate-enzyme covalent bond established between the C6 of the nucleotide substrate and a conserved cysteine residue. The calculations also reproduce and rationalize a normal H/T 2° kinetic isotope effect measured for that step. From a computational point of view, the results demonstrate that the use of an incomplete number of coordinates to describe the real reaction coordinate can render biased results.

Is Promiscuous CALB a Good Scaffold for Designing New Epoxidases?

Candida Antarctica lipase B (CALB) is a well-known enzyme, especially because of its promiscuous activity. Due to its properties, CALB was widely used as a benchmark for designing new catalysts for important organic reactions. The active site of CALB is very similar to that of soluble epoxide hydrolase (sEH) formed by a nucleophile-histidine-acid catalytic triad and an oxyanion hole typical for molecular structures derived from processes of α/ß hydrolases. In this work we are exploring these similarities and proposing a Ser105Asp variant of CALB as a new catalyst for epoxide hydrolysis. In particular, the hydrolysis of the trans-diphenylpropene oxide (t-DPPO) is studied by means of quantum cluster models mimicking the active site of both enzymes. Our results, based on semi-empirical and DFT calculations, suggest that mutant Ser105Asp CALB is a good protein scaffold to be used for the bio-synthesis of chiral compounds.

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.