Protein Flexibility and Stiffness Enable Efficient Enzymatic Catalysis

A comparative study of two α-L-rhamnosidases with high sequence identity

The GH78 α-L-rhamnosidase from Aspergillus tubingensis (AT-Rha) was proved to be a new clade of Aspergillus α-L-rhamnosidases in the previous study. A putative α-L-rhamnosidase from A. kawachii IFO 4308 (AK-Rha) has 92 % identity in amino acid sequence with AT-Rha. In this study, AK-Rha was expressed in P. pastoris and characterized. Similar to AT-rRha, the recombinant AK-Rha (AK-rRha) showed a narrow substrate specificity to naringin. Interestingly, the enzyme activity of AK-rRha was 0.816 U/mg toward naringin, significantly lower than 125.142 U/mg of AT-rRha. Their large differences in catalytic efficiency was mainly due to their differences in kcat values between AK-rRha (0.67 s-1) and AT-rRha (4.89 × 104 s-1). The molecular dynamics simulation exhibited that the overall conformation of AK-Rha was rigid and that of AT-Rha was flexible; the Loop Y-L located above the catalytic domain formed different steric hindrances to naringin, and interacted with the flavonoid matrices at different strengths. The polar solvation energy analysis implied that the glycosidic bond was more easily hydrolysed in AT-Rha. The comparative study verified that the main feature of AK-Rha and AT-Rha represented Aspergillus α-L-rhamnosidase was the narrow substrate specificity toward naringin, and provided an insight of the relationships between their catalytic abilities and structures.

Molecular Mechanism of Double-Displacement Retaining β-Kdo Glycosyltransferase WbbB

Glycosyltransferases (GTs) are pivotal enzymes involved in glycosidic bond synthesis, which can lead to either retention or inversion of the glycosyl moiety's anomeric configuration. However, the catalytic mechanism for retaining GTs remains a subject of controversy. In this study, we employ MD and QM/MM metadynamics to investigate the double-displacement catalytic mechanism of the retaining β-Kdo transferase WbbB. Our findings demonstrate that the nucleophile Asp232 initiates the reaction by attacking the sugar ring containing a carboxylate at the anomeric position, forming a covalent adduct. Subsequently, the adduct undergoes a rotational rearrangement, ensuring proper orientation of the anomeric carbon for the acceptor substrate. In the second step, Glu158 acts as the catalytic base to abstract the proton of the acceptor substrate to complete the transglycosylation reaction. Notably, His265 does not function as the anticipated catalytic acid; instead, it stabilizes the phosphate group through H-bonding interactions. Our simulations support the double-displacement mechanism implicated from the crystallographic studies of WbbB. This mechanism deviates from the common SNi-type and retaining glycoside hydrolase mechanisms, thereby expanding our understanding of GT catalytic mechanisms.

Biocatalytic approach for the synthesis of chiral alcohols for the development of pharmaceutical intermediates and other industrial applications: A review

Biocatalysis has emerged as a strong tool for the synthesis of active pharmaceutical ingredients (APIs). In the early twentieth century, whole cell biocatalysis was used to develop the first industrial biocatalytic processes, and the precise work of enzymes was unknown. Biocatalysis has evolved over the years into an essential tool for modern, cost-effective, and sustainable pharmaceutical manufacturing. Meanwhile, advances in directed evolution enable the rapid production of process-stable enzymes with broad substrate scope and high selectivity. Large-scale synthetic pathways incorporating biocatalytic critical steps towards >130 APIs of authorized pharmaceuticals and drug prospects are compared in terms of steps, reaction conditions, and scale with the corresponding chemical procedures. This review is designed on the functional group developed during the reaction forming alcohol functional groups. Some important biocatalyst sources, techniques, and challenges are described. A few APIs and their utilization in pharmaceutical drugs are explained here in this review. Biocatalysis has provided shorter, more efficient, and more sustainable alternative pathways toward existing small molecule APIs. Furthermore, non-pharmaceutical applications of biocatalysts are also mentioned and discussed. Finally, this review includes the future outlook and challenges of biocatalysis. In conclusion, Further research and development of promising enzymes are required before they can be used in industry.

Four Parallel Pathways in T4 Ligase-Catalyzed Repair of Nicked DNA with Diverse Bending Angles

The structural diversity of biological macromolecules in different environments contributes complexity to enzymological processes vital for cellular functions. Fluorescence resonance energy transfer and electron microscopy are used to investigate the enzymatic reaction of T4 DNA ligase catalyzing the ligation of nicked DNA. The data show that both the ligase-AMP complex and the ligase-AMP-DNA complex can have four conformations. This finding suggests the parallel occurrence of four ligation reaction pathways, each characterized by specific conformations of the ligase-AMP complex that persist in the ligase-AMP-DNA complex. Notably, these complexes have DNA bending angles of ≈0°, 20°, 60°, or 100°. The mechanism of parallel reactions challenges the conventional notion of simple sequential reaction steps occurring among multiple conformations. The results provide insights into the dynamic conformational changes and the versatile attributes of T4 DNA ligase and suggest that the parallel multiple reaction pathways may correspond to diverse T4 DNA ligase functions. This mechanism may potentially have evolved as an adaptive strategy across evolutionary history to navigate complex environments.

A curated rotamer library for common post-translational modifications of proteins

MOTIVATION

Sidechain rotamer libraries of the common amino acids of a protein are useful for folded protein structure determination and for generating ensembles of intrinsically disordered proteins (IDPs). However, much of protein function is modulated beyond the translated sequence through the introduction of post-translational modifications (PTMs).

RESULTS

In this work, we have provided a curated set of side chain rotamers for the most common PTMs derived from the RCSB PDB database, including phosphorylated, methylated, and acetylated sidechains. Our rotamer libraries improve upon existing methods such as SIDEpro, Rosetta, and AlphaFold3 in predicting the experimental structures for PTMs in folded proteins. In addition, we showcase our PTM libraries in full use by generating ensembles with the Monte Carlo Side Chain Entropy (MCSCE) for folded proteins, and combining MCSCE with the Local Disordered Region Sampling algorithms within IDPConformerGenerator for proteins with intrinsically disordered regions.

AVAILABILITY AND IMPLEMENTATION

The codes for dihedral angle computations and library creation are available at https://github.com/THGLab/ptm_sc.git.

Exploring Enzyme-Mimicking Metal-Organic Frameworks for CO2 Conversion through Vibrational Spectra-Based Machine Learning

In pursuing the benefits of natural enzyme catalysts while overcoming their limitations, we find metal-organic frameworks (MOFs), renowned for their highly tunable functionalities, stand out in biomimetic applications. We used unsupervised machine learning on density functional theory-computed vibrational infrared and Raman spectral features to screen 300 Zn-MOFs for CO2 conversion, similar to carbonic anhydrase (CA). Our findings confirmed that MOFs with spectroscopic attributes closely resembling those of CA hold the potential for replicating CA's electronic and catalytic properties. Unlike previous studies that relied on heuristic or trial-and-error methods and focused on geometric configurations, our research uses vibrational spectral features to explore structure-property relationships, making them more accessible through spectroscopy. Moreover, we highlight vibrational spectral features as efficient carriers for highly dimensional chemical information, enabling the simultaneous optimization of multiple performance parameters. These findings pave the way for pioneering designs of enzyme-mimetic MOFs and concurrently expand the application scope of spectroscopic tools in biomimetic catalysis.

A role for conformational changes in enzyme catalysis

The role played by conformational changes in enzyme catalysis is controversial. In addition to examining specific enzymes, studying formal models can help identify the conditions under which conformational changes promote catalysis. Here, we present a model demonstrating how conformational changes can break a generic trade-off due to the conflicting requirements of successive steps in catalytic cycles, namely high specificity for the transition state to accelerate the chemical transformation and low affinity for the products to favor their release. The mechanism by which the trade-off is broken is a transition between conformations with different affinities for the substrate. The role of the effector that induces the transition is played by a substrate "handle," a part of the substrate that is not chemically transformed but whose interaction with the enzyme is nevertheless essential to rapidly complete the catalytic cycle. A key element of the model is the formalization of the constraints causing the trade-off that the presence of multiple states breaks, which we attribute to the strong chemical similarity between successive reaction states-substrates, transition states, and products. For the sake of clarity, we present our model for irreversible one-step unimolecular reactions. In this context, we demonstrate how the different forms that chemical similarities between reaction states can take impose limits on the overall catalytic turnover. We first analyze catalysts without internal degrees of freedom and then show how two-state catalysts can overcome their limitations. Our results recapitulate previous proposals concerning the role of conformational changes and substrate handles in a formalism that makes explicit the constraints that elicit these features. In addition, our approach establishes links with studies in the field of heterogeneous catalysis, where the same trade-offs are observed and where overcoming them is a well-recognized challenge.

A Curated Rotamer Library for Common Post-Translational Modifications of Proteins

Sidechain rotamer libraries of the common amino acids of a protein are useful for folded protein structure determination and for generating ensembles of intrinsically disordered proteins (IDPs). However much of protein function is modulated beyond the translated sequence through thFiguree introduction of post-translational modifications (PTMs). In this work we have provided a curated set of side chain rotamers for the most common PTMs derived from the RCSB PDB database, including phosphorylated, methylated, and acetylated sidechains. Our rotamer libraries improve upon existing methods such as SIDEpro and Rosetta in predicting the experimental structures for PTMs in folded proteins. In addition, we showcase our PTM libraries in full use by generating ensembles with the Monte Carlo Side Chain Entropy (MCSCE) for folded proteins, and combining MCSCE with the Local Disordered Region Sampling algorithms within IDPConformerGenerator for proteins with intrinsically disordered regions.

Glycerol 3-Phosphate Dehydrogenase: Role of the Protein Conformational Change in Activation of a Readily Reversible Enzyme-Catalyzed Hydride Transfer Reaction

Kinetic parameters are reported for glycerol 3-phosphate dehydrogenase (GPDH)-catalyzed hydride transfer from the whole substrate glycerol 3-phosphate (G3P) or truncated substrate ethylene glycol (EtG) to NAD, and for activation of the hydride transfer reaction of EtG by phosphite dianion. These kinetic parameters were combined with parameters for enzyme-catalyzed hydride transfer in the microscopic reverse direction to give the reaction equilibrium constants Keq. Hydride transfer from G3P is favored in comparison to EtG because the carbonyl product of the former reaction is stabilized by hyperconjugative electron donation from the -CH2R keto substituent. The kinetic data show that the phosphite dianion provides the same 7.6 ± 0.1 kcal/mol stabilization of the transition states for enzyme-catalyzed reactions in the forward [reduction of NAD by EtG] and reverse [oxidation of NADH by glycolaldehyde] directions. The experimental evidence that supports a role for phosphite dianion in stabilizing the active closed form of the GPDH (EC) relative to the ca. 6 kcal/mol more unstable open form (EO) is summarized.

Activation and friction in enzymatic loop opening and closing dynamics

Protein loop dynamics have recently been recognized as central to enzymatic activity, specificity and stability. However, the factors controlling loop opening and closing kinetics have remained elusive. Here, we combine molecular dynamics simulations with string-method determination of complex reaction coordinates to elucidate the molecular mechanism and rate-limiting step for WPD-loop dynamics in the PTP1B enzyme. While protein conformational dynamics is often represented as diffusive motion hindered by solvent viscosity and internal friction, we demonstrate that loop opening and closing is activated. It is governed by torsional rearrangement around a single loop peptide group and by significant friction caused by backbone adjustments, which can dynamically trap the loop. Considering both torsional barrier and time-dependent friction, our calculated rate constants exhibit very good agreement with experimental measurements, reproducing the change in loop opening kinetics between proteins. Furthermore, we demonstrate the applicability of our results to other enzymatic loops, including the M20 DHFR loop, thereby offering prospects for loop engineering potentially leading to enhanced designs.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model. Here, we investigate protein folding, unfolding, and misfolding of fast-folding proteins by examining different force fields with their recommended water models, i.e., ff14SB with the TIP3P model and ff19SB with the OPC model. To this end, we generated long conventional MD simulations highlighting the perks and pitfalls of these setups. Using Markov state models, we defined kinetically independent conformational substates and emphasized their distinct characteristics, as well as their corresponding state probabilities. Surprisingly, we found substantial differences in thermodynamics and kinetics of protein folding, depending on the combination of the protein force field and water model, originating primarily from the different water models. These results emphasize the importance of carefully choosing the force field and the respective water model as they determine the accuracy of the observed dynamics of folding events. Thus, the findings support the hypothesis that the water model is at least equally important as the force field and hence needs to be considered in future studies investigating protein dynamics and folding in all areas of biophysics.

Stepwise introduction of stabilizing mutations reveals nonlinear additive effects in de novo TIM barrels

Over the past decades, the TIM-barrel fold has served as a model system for the exploration of how changes in protein sequences affect their structural, stability, and functional characteristics, and moreover, how this information can be leveraged to design proteins from the ground up. After numerous attempts to design de novo proteins with this specific fold, sTIM11 was the first validated de novo design of an idealized four-fold symmetric TIM barrel. Subsequent efforts to enhance the stability of this initial design resulted in the development of DeNovoTIMs, a family of de novo TIM barrels with various stabilizing mutations. In this study, we present an investigation into the biophysical and thermodynamic effects upon introducing a varying number of stabilizing mutations per quarter along the sequence of a four-fold symmetric TIM barrel. We compared the base design DeNovoTIM0 without any stabilizing mutations with variants containing mutations in one, two, three, and all four quarters-designated TIM1q, TIM2q, TIM3q, and DeNovoTIM6, respectively. This analysis revealed a stepwise and nonlinear change in the thermodynamic properties that correlated with the number of mutated quarters, suggesting positive nonadditive effects. To shed light on the significance of the location of stabilized quarters, we engineered two variants of TIM2q which contain the same number of mutations but positioned in different quarter locations. Characterization of these TIM2q variants revealed that the mutations exhibit varying effects on the overall protein stability, contingent upon the specific region in which they are introduced. These findings emphasize that the amount and location of stabilized interfaces among the four quarters play a crucial role in shaping the conformational stability of these four-fold symmetric TIM barrels. Analysis of de novo proteins, as described in this study, enhances our understanding of how sequence variations can finely modulate stability in both naturally occurring and computationally designed proteins.

The Arylamine N-Acetyltransferases as Therapeutic Targets in Metabolic Diseases Associated with Mitochondrial Dysfunction

In humans, there are two arylamine N-acetyltransferase genes that encode functional enzymes (NAT1 and NAT2) as well as one pseudogene, all of which are located together on chromosome 8. Although they were first identified by their role in the acetylation of drugs and other xenobiotics, recent studies have shown strong associations for both enzymes in a variety of diseases, including cancer, cardiovascular disease, and diabetes. There is growing evidence that this association may be causal. Consistently, NAT1 and NAT2 are shown to be required for healthy mitochondria. This review discusses the current literature on the role of both NAT1 and NAT2 in mitochondrial bioenergetics. It will attempt to relate our understanding of the evolution of the two genes with biologic function and then present evidence that several major metabolic diseases are influenced by NAT1 and NAT2. Finally, it will discuss current and future approaches to inhibit or enhance NAT1 and NAT2 activity/expression using small-molecule drugs. SIGNIFICANCE STATEMENT: The arylamine N-acetyltransferases (NATs) NAT1 and NAT2 share common features in their associations with mitochondrial bioenergetics. This review discusses mitochondrial function as it relates to health and disease, and the importance of NAT in mitochondrial function and dysfunction. It also compares NAT1 and NAT2 to highlight their functional similarities and differences. Both NAT1 and NAT2 are potential drug targets for diseases where mitochondrial dysfunction is a hallmark of onset and progression.

Optimizing Trilobatin Production via Screening and Modification of Glycosyltransferases

Trilobatin (TBL) is a key sweet compound from the traditional Chinese sweet tea plant (Rubus suavissimus S. Lee). Because of its intense sweetness, superior taste profile, and minimal caloric value, it serves as an exemplary natural dihydrochalcone sweetener. It also has various health benefits, including anti-inflammatory and glucose-lowering effects. It is primarily produced through botanical extraction, which impedes its scalability and cost-effectiveness. In a novel biotechnological approach, phloretin is used as a precursor that is transformed into TBL by the glycosyltransferase enzyme ph-4'-OGT. However, this enzyme's low catalytic efficiency and by-product formation limit the large-scale synthesis of TBL. In our study, the enzyme Mdph-4'-OGT was used to screen 17 sequences across species for TBL synthesis, of which seven exhibited catalytic activity. Notably, PT577 exhibited an unparalleled 97.3% conversion yield within 3 h. We then optimized the reaction conditions of PT577, attaining a peak TBL bioproduction of 163.3 mg/L. By employing virtual screening, we identified 25 mutation sites for PT577, thereby creating mutant strains that reduced by-products by up to 50%. This research enhances the enzymatic precision for TBL biosynthesis and offers a robust foundation for its industrial-scale production, with broader implications for the engineering and in silico analysis of glycosyltransferases.

Targeted Mutation of a Non-catalytic Gating Residue Increases the Rate of Pseudomonas aeruginosa d-Arginine Dehydrogenase Catalytic Turnover

Commercial food and l-amino acid industries rely on bioengineered d-amino acid oxidizing enzymes to detect and remove d-amino acid contaminants. However, the bioengineering of enzymes to generate faster biological catalysts has proven difficult as a result of the failure to target specific kinetic steps that limit enzyme turnover, kcat, and the poor understanding of loop dynamics critical for catalysis. Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) oxidizes most d-amino acids and is a good candidate for application in the l-amino acid and food industries. The side chain of the loop L2 E246 residue located at the entrance of the PaDADH active site pocket potentially favors the closed active site conformation and secures the substrate upon binding. This study used site-directed mutagenesis, steady-state, and rapid reaction kinetics to generate the glutamine, glycine, and leucine variants and investigate whether increasing the rate of product release could translate to an increased enzyme turnover rate. Upon E246 mutation to glycine, there was an increased rate of d-arginine turnover kcat from 122 to 500 s-1. Likewise, the kcat values increased 2-fold for the glutamine or leucine variants. Thus, we have engineered a faster biocatalyst for industrial applications by selectively increasing the rate of the PaDADH product release.

The Relationship Between the Cross-Linker on Chitosan-Coated Magnetic Nanoparticles and the Properties of Immobilized Papain

The immobilized enzymes' properties can be affected by cross-linkers on the surface of supports. To study how cross-linkers alter enzymes function, chitosan-coated magnetic nanoparticles (CMNPs) with immobilized papain were prepared using glutaraldehyde and or genipin, and then, the properties of the nanoparticles and the immobilized enzymes were assessed. The Scanning Electron Microscope (SEM), Fourier Transform Infrared (FTIR), and X-Ray Diffraction (XRD) results showed that the CMNPs were prepared and papain molecules were immobilized on CMNPs by glutaraldehyde (CMNP-Glu-Papain) or by genipin (CMNP-Gen-Papain). Also, the results associated with enzymes activity indicated that the immobilization by glutaraldehyde and genipin increased the pH optimum of papain from 7 to 7.5 and 9, respectively. The kinetic results indicated that the immobilization by genipin slightly affects the enzyme affinity to the substrate. The stability results showed that CMNP-Gen-Papain has more thermal stability than CMNP-Glu-Papain and papain immobilization on CMNPs by genipin leads to stabilization of the enzyme in the presence of polar solvents, probably due to the more hydroxyl groups on CMNPs activated by genipin. In conclusion, this study suggests that there is a relationship between the types of cross-linker on the surface of supports, and the mechanism of action, kinetic parameters, and the stability of immobilized papain.

Triosephosphate Isomerase: The Crippling Effect of the P168A/I172A Substitution at the Heart of an Enzyme Active Site

The P168 and I172 side chains sit at the heart of the active site of triosephosphate isomerase (TIM) and play important roles in the catalysis of the isomerization reaction. The phosphodianion of substrate glyceraldehyde 3-phosphate (GAP) drives a conformational change at the TIM that creates a steric interaction with the P168 side chain that is relieved by the movement of P168 that carries the basic E167 side chain into a clamp that consists of the hydrophobic I172 and L232 side chains. The P168A/I172A substitution at TIM from Trypanosoma brucei brucei (TbbTIM) causes a large 120,000-fold decrease in kcat for isomerization of GAP that eliminates most of the difference in the reactivity of TIM compared to the small amine base quinuclidinone for deprotonation of catalyst-bound GAP. The I172A substitution causes a > 2-unit decrease in the pKa of the E167 carboxylic acid in a complex to the intermediate analog PGA, but the P168A substitution at the I172A variant has no further effect on this pKa. The P168A/I172A substitutions cause a 5-fold decrease in Km for the isomerization of GAP from a 0.9 kcal/mol stabilization of the substrate Michaelis complexes. The results show that the P168 and I172 side chains play a dual role in destabilizing the ground-state Michaelis complex to GAP and in promoting stabilization of the transition state for substrate isomerization. This is consistent with an important role for these side chains in an induced fit reaction mechanism [Richard, J. P. (2022) Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution. Biochemistry 61, 1533-1542].

How Flexibility Can Enhance Catalysis

Conformational changes are observed in many enzymes, but their role in catalysis is highly controversial. Here we present a theoretical model that illustrates how rigid catalysts can be fundamentally limited and how a conformational change induced by substrate binding can overcome this limitation, ultimately enabling barrier-free catalysis. The model is deliberately minimal, but the principle it illustrates is general and consistent with unique features of proteins as well as with previous informal proposals to explain the superiority of enzymes over other classes of catalysts. Implementing the discriminative switch suggested by the model could help overcome limitations currently encountered in the design of artificial catalysts.

Utilization of Cofactor Binding Energy for Enzyme Catalysis: Formate Dehydrogenase-Catalyzed Reactions of the Whole NAD Cofactor and Cofactor Pieces

The pressure to optimize enzymatic rate accelerations has driven the evolution of the induced-fit mechanism for enzyme catalysts where the binding interactions of nonreacting phosphodianion or adenosyl substrate pieces drive enzyme conformational changes to form protein substrate cages that are activated for catalysis. We report the results of experiments to test the hypothesis that utilization of the binding energy of the adenosine 5'-diphosphate ribose (ADP-ribose) fragment of the NAD cofactor to drive a protein conformational change activates Candida boidinii formate dehydrogenase (CbFDH) for catalysis of hydride transfer from formate to NAD+. The ADP-ribose fragment provides a >14 kcal/mol stabilization of the transition state for CbFDH-catalyzed hydride transfer from formate to NAD+. This is larger than the ca. 6 kcal/mol stabilization of the ground-state Michaelis complex between CbFDH and NAD+ (KNAD = 0.032 mM). The ADP, AMP, and ribose 5'-phosphate fragments of NAD+ activate CbFDH for catalysis of hydride transfer from formate to nicotinamide riboside (NR). At a 1.0 M standard state, these activators stabilize the hydride transfer transition states by ≈5.5 (ADP), 5.5 (AMP), and 4.4 (ribose 5'-phosphate) kcal/mol. We propose that activation by these cofactor fragments is partly or entirely due to the ion-pair interaction between the guanidino side chain cation of R174 and the activator phosphate anion. This substitutes for the interaction between the α-adenosyl pyrophosphate anion of the whole NAD+ cofactor that holds CbFDH in the catalytically active closed conformation.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Drug Design in the Exascale Era: A Perspective from Massively Parallel QM/MM Simulations

The initial phases of drug discovery - in silico drug design - could benefit from first principle Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics (MD) simulations in explicit solvent, yet many applications are currently limited by the short time scales that this approach can cover. Developing scalable first principle QM/MM MD interfaces fully exploiting current exascale machines - so far an unmet and crucial goal - will help overcome this problem, opening the way to the study of the thermodynamics and kinetics of ligand binding to protein with first principle accuracy. Here, taking two relevant case studies involving the interactions of ligands with rather large enzymes, we showcase the use of our recently developed massively scalable Multiscale Modeling in Computational Chemistry (MiMiC) QM/MM framework (currently using DFT to describe the QM region) to investigate reactions and ligand binding in enzymes of pharmacological relevance. We also demonstrate for the first time strong scaling of MiMiC-QM/MM MD simulations with parallel efficiency of ∼70% up to >80,000 cores. Thus, among many others, the MiMiC interface represents a promising candidate toward exascale applications by combining machine learning with statistical mechanics based algorithms tailored for exascale supercomputers.

The Role of Asn11 in Catalysis by Triosephosphate Isomerase

Four catalytic amino acids at triosephosphate isomerase (TIM) are highly conserved: N11, K13, H95, and E167. Asparagine 11 is the last of these to be characterized in mutagenesis studies. The ND2 side chain atom of N11 is hydrogen bonded to the O-1 hydroxyl of enzyme-bound dihydroxyacetone phosphate (DHAP), and it sits in an extended chain of hydrogen-bonded side chains that includes T75' from the second subunit. The N11A variants of wild-type TIM from Trypanosoma brucei brucei (TbbTIM) and Leishmania mexicana (LmTIM) undergo dissociation from the dimer to monomer under our assay conditions. Values of Kas = 8 × 103 and 1 × 106 M-1, respectively, were determined for the conversion of monomeric N11A TbbTIM and LmTIM into their homodimers. The N11A substitution at the variant of LmTIM previously stabilized by the E65Q substitution gives the N11A/E65Q variant that is stable to dissociation under our assay conditions. The X-ray crystal structure of N11A/E65Q LmTIM shows an active site that is essentially superimposable on that for wild-type TbbTIM, which also has a glutamine at position 65. A comparison of the kinetic parameters for E65Q LmTIM and N11A/E65Q LmTIM-catalyzed reactions of (R)-glyceraldehyde 3-phosphate (GAP) and (DHAP) shows that the N11A substitution results in a (13-14)-fold decrease in kcat/Km for substrate isomerization and a similar decrease in kcat for DHAP but only a 2-fold decrease in kcat for GAP.

Probing the mechanism of flavin action in the oxidative decarboxylation catalyzed by salicylate hydroxylase

Salicylate hydroxylase (NahG) is a FAD-dependent monooxygenase in which the reduced flavin activates O₂ coupled to the oxidative decarboxylation of salicylate to catechol or uncoupled from substrate oxidation to afford H₂O₂. This chapter presents different methodologies in equilibrium studies, steady-state kinetics, and identification of reaction products, which were important to understand the S_EAr mechanism of catalysis in NahG, the role of the different FAD parts for ligand binding, the extent of uncoupled reaction, and the catalysis of salicylate's oxidative decarboxylation. These features are likely familiar to many other FAD-dependent monooxygenases and offer a potential asset for developing new tools and strategies in catalysis.

Rational Construction of Protein-Mimetic Nano-Switch Systems Based on Secondary Structure Transitions of Synthetic Polypeptides

The manipulation of the flexibility/rigidity of polymeric chains to control their function is commonly observed in natural macromolecules but largely unexplored in synthetic systems. Herein, we construct a series of protein-mimetic nano-switches consisting of a gold nanoparticle (GNP) core, a synthetic polypeptide linker, and an optically functional molecule (OFM), whose biological function can be dynamically regulated by the flexibility of the polypeptide linker. At the dormant state, the polypeptide adopts a flexible, random-coiled conformation, bringing GNP and OFM in close proximity that leads to the "turn-off" of the OFM. Once treated with alkaline phosphatase (ALP), the nano-switches are activated due to the increased separation distance between GNP and OFM driven by the coil-to-helix and flexible-to-rigid transition of the polypeptide linker. The nano-switches therefore enable selective fluorescence imaging or photodynamic therapy in response to ALP overproduced by tumor cells. The control over polymer flexibility represents an effective strategy to manipulate the optical activity of nano-switches, which mimics the delicate structure-property relationship of natural proteins.

Kinetics and mechanism for enzyme-catalyzed reactions of substrate pieces

The most important difference between enzyme and small molecule catalysts is that only enzymes utilize the large intrinsic binding energies of nonreacting portions of the substrate in stabilization of the transition state for the catalyzed reaction. A general protocol is described to determine the intrinsic phosphodianion binding energy for enzymatic catalysis of reactions of phosphate monoester substrates, and the intrinsic phosphite dianion binding energy in activation of enzymes for catalysis of phosphodianion truncated substrates, from the kinetic parameters for enzyme-catalyzed reactions of whole and truncated substrates. The enzyme-catalyzed reactions so-far documented that utilize dianion binding interactions for enzyme activation; and, their phosphodianion truncated substrates are summarized. A model for the utilization of dianion binding interactions for enzyme activation is described. The methods for the determination of the kinetic parameters for enzyme-catalyzed reactions of whole and truncated substrates, from initial velocity data, are described and illustrated by graphical plots of kinetic data. The results of studies on the effect of site-directed amino acid substitutions at orotidine 5'-monophosphate decarboxylase, triosephosphate isomerase, and glycerol-3-phosphate dehydrogenase provide strong support for the proposal that these enzymes utilize binding interactions with the substrate phosphodianion to hold the protein catalysts in reactive closed conformations.

Mutation of a distal gating residue modulates NADH binding in NADH:Quinone oxidoreductase from Pseudomonas aeruginosa PAO1

Enzymes require flexible regions to adopt multiple conformations during catalysis. The mobile regions of enzymes include gates that modulate the passage of molecules in and out of the enzyme's active site. The enzyme PA1024 from Pseudomonas aeruginosa PA01 is a recently discovered flavin-dependent NADH:quinone oxidoreductase (NQO, EC 1.6.5.9). Q80 in loop 3 (residues 75-86) of NQO is ∼15 Å away from the flavin and creates a gate that seals the active site through a hydrogen bond with Y261 upon NADH binding. In this study, we mutated Q80 to glycine, leucine, or glutamate to investigate the mechanistic significance of distal residue Q80 in NADH binding in the active site of NQO. The UV-visible absorption spectrum reveals that the mutation of Q80 minimally affects the protein microenvironment surrounding the flavin. The anaerobic reductive half-reaction of the NQO-mutants yields a ≥25-fold increase in the Kd value for NADH compared to the WT enzyme. However, we determined that the kred value was similar in the Q80G, Q80L, and wildtype enzymes and only ∼25% smaller in the Q80E enzyme. Steady-state kinetics with NQO-mutants and NQO-WT at varying concentrations of NADH and 1,4-benzoquinone establish a ≤5-fold decrease in the kcat/KNADH value. Moreover, there is no significant difference in the kcat/KBQ (∼1 × 106 M-1s-1) and kcat (∼24 s-1) values in NQO-mutants and NQO-WT. These results are consistent with the distal residue Q80 being mechanistically essential for NADH binding to NQO with minimal effect on the quinone binding to the enzyme and hydride transfer from NADH to flavin.

Underlying Role of Hydrophobic Environments in Tuning Metal Elements for Efficient Enzyme Catalysis

The catalytic functions of metalloenzymes are often strongly correlated with metal elements in the active sites. However, dioxygen-activating nonheme quercetin dioxygenases (QueD) are found with various first-row transition-metal ions when metal swapping inactivates their innate catalytic activity. To unveil the molecular basis of this seemingly promiscuous yet metal-specific enzyme, we transformed manganese-dependent QueD into a nickel-dependent enzyme by sequence- and structure-based directed evolution. Although the net effect of acquired mutations was primarily to rearrange hydrophobic residues in the active site pocket, biochemical, kinetic, X-ray crystallographic, spectroscopic, and computational studies suggest that these modifications in the secondary coordination spheres can adjust the electronic structure of the enzyme-substrate complex to counteract the effects induced by the metal substitution. These results explicitly demonstrate that such noncovalent interactions encrypt metal specificity in a finely modulated manner, revealing the underestimated chemical power of the hydrophobic sequence network in enzyme catalysis.

Enzymatische C4-Epimerisierung von UDP-Glucuronsäure: präzise gesteuerte Rotation eines transienten 4-Ketointermediats für eine invertierende Reaktion ohne Decarboxylierung

UDP‐Glucuronsäure(UDP‐GlcA)‐4‐Epimerase repräsentiert eine wichtige Fragestellung in der Enzymkatalyse: die Balance zwischen konformativer Flexibilität und genauer Positionierung. Das Enzym koordiniert die C4‐Oxidation des Substrats durch NAD+ mit der Rotation eines leicht decarboxylierbaren β‐Ketosäure‐Intermediats im aktiven Zentrum zur Ermöglichung der stereoinvertierenden Reduktion der Ketogruppe durch NADH. Wir zeigen hier die nur schwer erfassbare Rotationskoordinate des 4‐Ketointermediats. Distorsion des Zuckerrings in eine Boot‐Konformation erzeugt torsionale Mobilität in der Bindungstasche des Enzyms. Die Endpunkte der Rotation zeigen den 4‐Ketozucker in einer unverformten 4C1‐Sesselkonformation. Die äquatorial positionierte Carboxylatgruppe ist ungünstig für die 4‐Ketozucker‐Decarboxylierung. Varianten der Epimerase zeigen Decarboxylierung, wenn sie die Bindung mit der Carboxylatgruppe im entgegengesetzten Rotationsisomer des Substrats entfernen. R185A/D‐Substitutionen wandeln die Epimerase in UDP‐Xylose‐Synthasen um, welche UDP‐GlcA in stereospezifischen, konfigurationserhaltenden Reaktionen decarboxylieren.

Enzymatic C4-Epimerization of UDP-Glucuronic Acid: Precisely Steered Rotation of a Transient 4-Keto Intermediate for an Inverted Reaction without Decarboxylation

UDP-glucuronic acid (UDP-GlcA) 4-epimerase illustrates an important problem regarding enzyme catalysis: balancing conformational flexibility with precise positioning. The enzyme coordinates the C4-oxidation of the substrate by NAD⁺ and rotation of a decarboxylation-prone β-keto acid intermediate in the active site, enabling stereoinverting reduction of the keto group by NADH. We reveal the elusive rotational landscape of the 4-keto intermediate. Distortion of the sugar ring into boat conformations induces torsional mobility in the enzyme's binding pocket. The rotational endpoints show that the 4-keto sugar has an undistorted ⁴ C₁ chair conformation. The equatorially placed carboxylate group disfavors decarboxylation of the 4-keto sugar. Epimerase variants lead to decarboxylation upon removal of the binding interactions with the carboxylate group in the opposite rotational isomer of the substrate. Substitutions R185A/D convert the epimerase into UDP-xylose synthases that decarboxylate UDP-GlcA in stereospecific, configuration-retaining reactions.

Protein adsorption by nanomechanical mass spectrometry: Beyond the real-time molecular weighting

During past decades, enormous progress in understanding the mechanisms of the intermolecular interactions between the protein and surface at the single-molecule level has been achieved. These advances could only be possible by the ongoing development of highly sophisticated experimental methods such as atomic force microscopy, optical microscopy, surface plasmon resonance, ellipsometry, quartz crystal microbalance, conventional mass spectrometry, and, more recently, the nanomechanical systems. Here, we highlight the main findings of recent studies on the label-free single-molecule (protein) detection by nanomechanical systems including those focusing on the protein adsorption on various substrate surfaces. Since the nanomechanical techniques are capable of detecting and manipulating proteins even at the single-molecule level, therefore, they are expected to open a new way of studying the dynamics of protein functions. It is noteworthy that, in contrast to other experimental methods, where only given protein properties like molecular weight or protein stiffness can be determined, the nanomechanical systems enable a real-time measurement of the multiple protein properties (e.g., mass, stiffness, and/or generated surface stress), making them suitable for the study of protein adsorption mechanisms. Moreover, we also discuss the possible future trends in label-free detection and analysis of dynamics of protein complexes with these nanomechanical systems.

Biocatalytic gateway to convert glycerol into 3-hydroxypropionic acid in waste-based biorefineries: Fundamentals, limitations, and potential research strategies

Microbial conversion of bioenergy-derived waste glycerol into value-added chemicals has emerged as an important bioprocessing technology due to its eco-friendliness, feasible technoeconomics, and potential to provide sustainability in biodiesel and bioethanol production. Glycerol is an abundant liquid waste from bioenergy plants with a projected volume of 6 million tons by 2025, accounting for about 10% of biodiesel and 2.5% of bioethanol yields. 3-Hydroxypropionic acid (3-HP) is a major product of glycerol bioconversion, which is the third largest biobased platform compound with expected market size and value of 3.6 million tons/year and USD 10 billion/year, respectively. Despite these biorefinery values, 3-HP biosynthesis from glycerol is still at an immature stage of commercial exploitation. The main challenges behind this immaturity are the toxic effects of 3-HPA on cells, the distribution of carbon flux to undesirable pathways, low tolerance of cells to glycerol and 3-HP, co-factor dependence of enzymes, low enzyme activity and stability, and the problems of substrate inhibition and specificity of enzymes. To address these challenges, it is necessary to understand the fundamentals of glycerol bioconversion and 3-HP production in terms of metabolic pathways, related enzymes, cell factories, midstream process configurations, and downstream 3-HP recovery, as discussed in this review critically and comprehensively. It is equally important to know the current challenges and limitations in 3-HP production, which are discussed in detail along with recent research efforts and remaining gaps. Finally, possible research strategies are outlined considering the recent technological advances in microbial biosynthesis, aiming to attract further research efforts to achieve a sustainable and industrially exploitable 3-HP production technology. By discussing the use of advanced tools and strategies to overcome the existing challenges in 3-HP biosynthesis, this review will attract researchers from many other similar biosynthesis technologies and provide a common gateway for their further development.

Adenylate Kinase-Catalyzed Reactions of AMP in Pieces: Specificity for Catalysis at the Nucleoside Activator and Dianion Catalytic Sites

The pressure to optimize the enzymatic rate acceleration for adenylate kinase (AK)-catalyzed phosphoryl transfer has led to the evolution of an induced-fit mechanism, where the binding energy from interactions between the protein and substrate adenosyl group is utilized to drive a protein conformational change that activates the enzyme for catalysis. The adenine group of adenosine contributes 11.8 kcal mol-1 to the total ≥14.7 kcal mol-1 adenosine stabilization of the transition state for AK-catalyzed phosphoryl transfer to AMP. The relative third-order rate constants for activation of adenylate kinase, by the C-5 truncated adenosine 1-(β-d-erythrofuranosyl)adenine (EA), for catalysis of phosphoryl transfer from ATP to phosphite dianion (HP, kcat/KHPKAct = 260 M-2 s-1), fluorophosphate (47 M-2 s-1), and phosphate (9.6 M-2 s-1), show that substitution of -F for -H and of -OH for -H at HP results, respectively, in decreases in the reactivity of AK for catalysis of phosphoryl transfer due to polar and steric effects of the -F and -OH substituents. The addition of a 5'-CH2OH to the EA activator results in a 3.0 kcal mol-1 destabilization of the transition state for AK-activated phosphoryl transfer to HP due to a steric effect. This is smaller than the 8.3 kcal mol-1 steric effect of the 5'-CH2OH substituent at OMP on HP-activated OMPDC-catalyzed decarboxylation of 1-(β-d-erythrofuranosyl)orotate. The 2'-OH ribosyl substituent shows significant interactions with the transition states for AK-catalyzed phosphoryl transfer from ATP to AMP and for adenosine-activated AK-catalyzed phosphoryl transfer from ATP to HP.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

Design and optimization of enzymatic activity in a de novo β-barrel scaffold

While native scaffolds offer a large diversity of shapes and topologies for enzyme engineering, their often unpredictable behavior in response to sequence modification makes de novo generated scaffolds an exciting alternative. Here we explore the customization of the backbone and sequence of a de novo designed eight stranded β‐barrel protein to create catalysts for a retro‐aldolase model reaction. We show that active and specific catalysts can be designed in this fold and use directed evolution to further optimize activity and stereoselectivity. Our results support previous suggestions that different folds have different inherent amenability to evolution and this property could account, in part, for the distribution of natural enzymes among different folds.

Temperature-dependent hydrogen deuterium exchange shows impact of analog binding on adenosine deaminase flexibility but not embedded thermal networks

The analysis of hydrogen deuterium exchange by mass spectrometry as a function of temperature and mutation (TDHDX-MS) has emerged as a generic and efficient tool for the spatial resolution of protein networks that are proposed to function in the thermal activation of catalysis. In this work, we extend TDHDX from apo-enzyme structures to protein-ligand complexes. Using adenosine deaminase as a prototype, we compared the impacts of a substrate analog (1-deaza-adenosine or DAA) and a very tight-binding inhibitor/transition state analog (pentostatin) at single and multiple temperatures. At a single temperature, we observed different HDX-MS properties for the two ligands, as expected from their 10⁶-fold differences in strength of binding. By contrast, analogous patterns for TDHDX-MS emerge in the presence of both DAA and pentostatin, indicating similar impacts of either ligand on the enthalpic barriers for local protein unfolding. We extended TDHDX to a function-altering mutant of adenosine deaminase in the presence of pentostatin and revealed a protein thermal network that is highly similar to that previously reported for the apo-enzyme (Gao et al., 2020, JACS 142, 19936-19949). Finally, we discuss the differential impacts of pentostatin binding on overall protein flexibility vs. site-specific thermal transfer pathways in the context of models for substrate-induced changes to a distributed protein conformational landscape that act in synergy with embedded protein thermal networks to achieve efficient catalysis.

Uniform binding and negative catalysis at the origin of enzymes

Enzymes are well known for their catalytic abilities, some even reaching "catalytic perfection" in the sense that the reaction they catalyze has reached the physical bound of the diffusion rate. However, our growing understanding of enzyme superfamilies has revealed that only some share a catalytic chemistry while others share a substrate-handle binding motif, for example, for a particular phosphate group. This suggests that some families emerged through a "substrate-handle-binding-first" mechanism ("binding-first" for brevity) instead of "chemistry-first" and we are, therefore, left to wonder what the role of non-catalytic binders might have been during enzyme evolution. In the last of their eight seminal, back-to-back articles from 1976, John Albery and Jeremy Knowles addressed the question of enzyme evolution by arguing that the simplest mode of enzyme evolution is what they defined as "uniform binding" (parallel stabilization of all enzyme-bound states to the same degree). Indeed, we show that a uniform-binding proto-catalyst can accelerate a reaction, but only when catalysis is already present, that is, when the transition state is already stabilized to some degree. Thus, we sought an alternative explanation for the cases where substrate-handle-binding preceded any involvement of a catalyst. We find that evolutionary starting points that exhibit negative catalysis can redirect the reaction's course to a preferred product without need for rate acceleration or product release; that is, if they do not stabilize, or even destabilize, the transition state corresponding to an undesired product. Such a mechanism might explain the emergence of "binding-first" enzyme families like the aldolase superfamily.

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution

Many enzymes that show a large specificity in binding the enzymatic transition state with a higher affinity than the substrate utilize substrate binding energy to drive protein conformational changes to form caged substrate complexes. These protein cages provide strong stabilization of enzymatic transition states. Using part of the substrate binding energy to drive the protein conformational change avoids a similar strong stabilization of the Michaelis complex and irreversible ligand binding. A seminal step in the development of modern enzyme catalysts was the evolution of enzymes that couple substrate binding to a conformational change. These include enzymes that function in glycolysis (triosephosphate isomerase), the biosynthesis of lipids (glycerol phosphate dehydrogenase), the hexose monophosphate shunt (6-phosphogluconate dehydrogenase), and the mevalonate pathway (isopentenyl diphosphate isomerase), catalyze the final step in the biosynthesis of pyrimidine nucleotides (orotidine monophosphate decarboxylase), and regulate the cellular levels of adenine nucleotides (adenylate kinase). The evolution of enzymes that undergo ligand-driven conformational changes to form active protein-substrate cages is proposed to proceed by selection of variants, in which the selected side chain substitutions destabilize a second protein conformer that shows compensating enhanced binding interactions with the substrate. The advantages inherent to enzymes that incorporate a conformational change into the catalytic cycle provide a strong driving force for the evolution of flexible protein folds such as the TIM barrel. The appearance of these folds represented a watershed event in enzyme evolution that enabled the rapid propagation of enzyme activities within enzyme superfamilies.

Glycerol-3-Phosphate Dehydrogenase: The K120 and K204 Side Chains Define an Oxyanion Hole at the Enzyme Active Site

The cationic K120 and K204 side chains lie close to the C-2 carbonyl group of substrate dihydroxyacetone phosphate (DHAP) at the active site of glycerol-3-phosphate dehydrogenase (GPDH), and the K120 side chain is also positioned to form a hydrogen bond to the C-1 hydroxyl of DHAP. The kinetic parameters for unactivated and phosphite dianion-activated GPDH-catalyzed reduction of glycolaldehyde and acetaldehyde (AcA) show that the transition state for the former reaction is stabilized by ca 5 kcal/mole by interactions of the C-1 hydroxyl group with the protein catalyst. The K120A and K204A substitutions at wild-type GPDH result in similar decreases in k_cat, but K_m is only affected by the K120A substitution. These results are consistent with 3 kcal/mol stabilizing interactions between the K120 or K204 side chains and a negative charge at the C-2 oxygen at the transition state for hydride transfer from NADH to DHAP. This stabilization resembles that observed at oxyanion holes for other enzymes. There is no detectable rescue of the K204A variant by ethylammonium cation (EtNH₃⁺), compared with the efficient rescue of the K120A variant. This is consistent with a difference in the accessibility of the variant enzyme active sites to exogenous EtNH₃⁺. The K120A/K204A substitutions cause a (6 × 10⁶)-fold increase in the promiscuity of wild-type hlGPDH for catalysis of the reduction of AcA compared to DHAP. This may reflect conservation of the active site for an ancestral alcohol dehydrogenase, whose relative activity for catalysis of reduction of AcA increases with substitutions that reduce the activity for reduction of the specific substrate DHAP.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

An Enzyme with High Catalytic Proficiency Utilizes Distal Site Substrate Binding Energy to Stabilize the Closed State but at the Expense of Substrate Inhibition

Understanding the factors that underpin the enormous catalytic proficiencies of enzymes is fundamental to catalysis and enzyme design. Enzymes are, in part, able to achieve high catalytic proficiencies by utilizing the binding energy derived from nonreacting portions of the substrate. In particular, enzymes with substrates containing a nonreacting phosphodianion group coordinated in a distal site have been suggested to exploit this binding energy primarily to facilitate a conformational change from an open inactive form to a closed active form, rather than to either induce ground state destabilization or stabilize the transition state. However, detailed structural evidence for the model is limited. Here, we use β-phosphoglucomutase (βPGM) to investigate the relationship between binding a phosphodianion group in a distal site, the adoption of a closed enzyme form, and catalytic proficiency. βPGM catalyzes the isomerization of β-glucose 1-phosphate to glucose 6-phosphate via phosphoryl transfer reactions in the proximal site, while coordinating a phosphodianion group of the substrate(s) in a distal site. βPGM has one of the largest catalytic proficiencies measured and undergoes significant domain closure during its catalytic cycle. We find that side chain substitution at the distal site results in decreased substrate binding that destabilizes the closed active form but is not sufficient to preclude the adoption of a fully closed, near-transition state conformation. Furthermore, we reveal that binding of a phosphodianion group in the distal site stimulates domain closure even in the absence of a transferring phosphoryl group in the proximal site, explaining the previously reported β-glucose 1-phosphate inhibition. Finally, our results support a trend whereby enzymes with high catalytic proficiencies involving phosphorylated substrates exhibit a greater requirement to stabilize the closed active form.

The role of remote flavin adenine dinucleotide pieces in the oxidative decarboxylation catalyzed by salicylate hydroxylase

Salicylate hydroxylase (NahG) has a single redox site in which FAD is reduced by NADH, the O₂ is activated by the reduced flavin, and salicylate undergoes an oxidative decarboxylation by a C(4a)-hydroperoxyflavin intermediate to give catechol. We report experimental results that show the contribution of individual pieces of the FAD cofactor to the observed enzymatic activity for turnover of the whole cofactor. A comparison of the kinetic parameters and products for the NahG-catalyzed reactions of FMN and riboflavin cofactor fragments reveal that the adenosine monophosphate (AMP) and ribitol phosphate pieces of FAD act to anchor the flavin to the enzyme and to direct the partitioning of the C(4a)-hydroperoxyflavin reaction intermediate towards hydroxylation of salicylate. The addition of AMP or ribitol phosphate pieces to solutions of the truncated flavins results in a partial restoration of the enzymatic activity lost upon truncation of FAD, and the pieces direct the reaction of the C(4a)-hydroperoxyflavin intermediate towards hydroxylation of salicylate.

Single-Molecule Sampling of Dihydrofolate Reductase Shows Kinetic Pauses and an Endosteric Effect Linked to Catalysis

The ability to sample multiple reactions on the same single enzyme is important to link rare intermediates with catalysis and to unravel the role of conformational changes. Despite decades of efforts, however, the single-molecule characterization of nonfluorogenic enzymes during multiple catalytic turnovers has been elusive. Here, we show that nanopore currents allow sampling the dynamic exchange between five structural intermediates during E. coli dihydrofolate reductase (DHFR) catalysis. We found that an endosteric effect promotes the binding of the substrate to the enzyme with a specific hierarchy. The chemical step then switched the enzyme from the closed to the occluded conformation, which in turn promotes the release of the reduced cofactor NADP⁺. Unexpectedly, only a few reactive complexes lead to catalysis. Furthermore, second-long catalytic pauses were observed, possibly reflecting an off-path conformation generated during the reaction. Finally, the free energy from multiple cofactor binding events were required to release the product and switch DHFR back to the reactive conformer. This catalytic fueled concerted mechanism is likely to have evolved to improve the catalytic efficiency of DHFR under the high concentrations of NADP⁺ in E. coli and might be a general feature for complex enzymatic reactions where the binding and release of the products must be tightly controlled.

A biomimetic metal–organic framework with cuboid inner cavities for enantioselective separation

A biomimetic metal–organic framework with cuboid inner cavities and multiple recognition sites was constructed from a phenylalanine-derived ligand. It can enantioselectively separate various racemic alcohols, diols and epoxides with ee up to 99.5%.

[(Flu)-(CH2)2-(NHC)-CH2-(NHC)-(CH2)2-(Flu)]2−: an ‘all-organic’ hybrid and flexible ligand that enwraps a Ca2+ pseudo-tetrahedrally

A conformationally flexible ligand, [(Flu)-(CH2)2-(NHC)-CH2-(NHC)-(CH2)2-(Flu)]2− (L), that offers four sequential organic donor modules and enwraps a Ca2+ pseudo-tetrahedrally in a twisted ‘S’-shape is devised.

An Efficient Protein Evolution Workflow for the Improvement of Bacterial PET Hydrolyzing Enzymes

Enzymatic degradation is a promising green approach to bioremediation and recycling of the polymer poly(ethylene terephthalate) (PET). In the past few years, several PET-hydrolysing enzymes (PHEs) have been discovered, and new variants have been evolved by protein engineering. Here, we report on a straightforward workflow employing semi-rational protein engineering combined to a high-throughput screening of variant libraries for their activity on PET nanoparticles. Using this approach, starting from the double variant W159H/S238F of Ideonella sakaiensis 201-F6 PETase, the W159H/F238A-ΔIsPET variant, possessing a higher hydrolytic activity on PET, was identified. This variant was stabilized by introducing two additional known substitutions (S121E and D186H) generating the TS-ΔIsPET variant. By using 0.1 mg mL⁻¹ of TS-ΔIsPET, ~10.6 mM of degradation products were produced in 2 days from 9 mg mL⁻¹ PET microparticles (~26% depolymerization yield). Indeed, TS-ΔIsPET allowed a massive degradation of PET nanoparticles (>80% depolymerization yield) in 1.5 h using only 20 μg of enzyme mL⁻¹. The rationale underlying the effect on the catalytic parameters due to the F238A substitution was studied by enzymatic investigation and molecular dynamics/docking analysis. The present workflow is a well-suited protocol for the evolution of PHEs to help generate an efficient enzymatic toolbox for polyester degradation.

Protein-Ribofuranosyl Interactions Activate Orotidine 5'-Monophosphate Decarboxylase for Catalysis

The role of a global, substrate-driven, enzyme conformational change in enabling the extraordinarily large rate acceleration for orotidine 5'-monophosphate decarboxylase (OMPDC)-catalyzed decarboxylation of orotidine 5'-monophosphate (OMP) is examined in experiments that focus on the interactions between OMPDC and the ribosyl hydroxyl groups of OMP. The D37 and T100' side chains of OMPDC interact, respectively, with the C-3' and C-2' hydroxyl groups of enzyme-bound OMP. D37G and T100'A substitutions result in 1.4 kcal/mol increases in the activation barrier ΔG^⧧ for catalysis of decarboxylation of the phosphodianion-truncated substrate 1-(β-d-erythrofuranosyl)orotic acid (EO) but result in larger 2.1-2.9 kcal/mol increases in ΔG^⧧ for decarboxylation of OMP and for phosphite dianion-activated decarboxylation of EO. This shows that these substitutions reduce transition-state stabilization by the Q215, Y217, and R235 side chains at the dianion binding site. The D37G and T100'A substitutions result in <1.0 kcal/mol increases in ΔG^⧧ for activation of OMPDC-catalyzed decarboxylation of the phosphoribofuranosyl-truncated substrate FO by phosphite dianions. Experiments to probe the effect of D37 and T100' substitutions on the kinetic parameters for d-glycerol 3-phosphate and d-erythritol 4-phosphate activators of OMPDC-catalyzed decarboxylation of FO show that ΔG^⧧ for sugar phosphate-activated reactions is increased by ca. 2.5 kcal/mol for each -OH interaction eliminated by D37G or T100'A substitutions. We conclude that the interactions between the D37 and T100' side chains and ribosyl or ribosyl-like hydroxyl groups are utilized to activate OMPDC for catalysis of decarboxylation of OMP, EO, and FO.

Adenylate Kinase-Catalyzed Reaction of AMP in Pieces: Enzyme Activation for Phosphoryl Transfer to Phosphite Dianion

The binding of adenosine 5'-triphosphate (ATP) and adenosine 5'-monophosphate (AMP) to adenylate kinase (AdK) drives closure of lids over the substrate adenosyl groups. We test the hypothesis that this conformational change activates AdK for catalysis. The rate constants for Homo sapiens adenylate kinase 1 (HsAdK1)-catalyzed phosphoryl group transfer to AMP, k_cat/K_m = 7.0 × 10⁶ M^-1 s^-1, and phosphite dianion, (k_HPi)_obs ≤1 × 10^-4 M^-1 s^-1, show that the binding energy of the adenosyl group effects a ≥7.0 × 10¹⁰-fold rate acceleration of phosphoryl transfer from ATP. The third-order rate constant of k_cat/K_HPiK_EA = 260 M^-2 s^-1 for 1-(β-d-erythrofuranosyl)adenine (EA)-activated phosphoryl transfer to phosphite dianion was determined, and the isohypophosphate reaction product characterized by ³¹P NMR. The results demonstrate the following: (i) a ≥14.7 kcal/mol stabilization of the transition state for phosphoryl transfer by the adenosyl group of AMP and a ≥2.6 × 10⁶-fold rate acceleration from the EA-driven conformational change and (ii) the recovery of ≥8.7 kcal/mol of this transition state stabilization for EA-activated phosphoryl transfer from ATP to phosphite.

Understanding enzyme behavior in a crowded scenario through modulation in activity, conformation and dynamics

Macromolecular crowding, inside the physiological interior, modulates the energy landscape of biological macromolecules in multiple ways. Amongst these, enzymes occupy a special place and hence understanding the function of the same in the crowded interior is of utmost importance. In this study, we have investigated the manner in which the multidomain enzyme, AK3L1 (PDB ID: 1ZD8), an isoform of adenylate kinase, has its features affected in presence of commonly used crowders (PEG 8, Dextran 40, Dextran 70, and Ficoll 70). Michaelis Menten plots reveal that the crowders in general enhance the activity of the enzyme, with the K_m and V_max values showing significant variations. Ficoll 70, induced the maximum activity for AK3L1 at 100 g/L, beyond which the activity reduced. Ensemble FRET studies were performed to provide insights into the relative domain (LID and CORE) displacements in presence of the crowders. Solvation studies reveal that the protein matrix surrounding the probe CPM (7-diethylamino-3-(4-maleimido-phenyl)-4-methylcoumarin) gets restricted in presence of the crowders, with Ficoll 70 providing the maximum rigidity, the same being linked to the decrease in the activity of the enzyme. Through our multipronged approach, we have observed a distinct correlation between domain displacement, enzyme activity and associated dynamics. Thus, keeping in mind the complex nature of enzyme activity and the surrounding bath of dense soup that the biological entity remains immersed in, indeed more such approaches need to be undertaken to have a better grasp of the "enzymes in the crowd".