Cooperative Electrostatic Interactions Drive Functional Evolution in the Alkaline Phosphatase Superfamily

Mechanistic and evolutionary insights into alkaline phosphatase superfamily through structure-function studies on Sphingomonas alkaline phosphatase

Alkaline phosphatases (APs), represented by E. coli AP (ECAP), employ an arginine residue to stabilize the phosphoryl group in the active site; whereas, AP from Sphingomonas (SPAP) shows a unique combination of substrate-binding residues; Thr89, Asn110, Lys171, and Arg173. Although such combination has been observed only in SPAP, these residues are present separately in different members of the AP superfamily. Here, we establish the presence of two distinct classes of APs; ECAP-type and SPAP-type. Bioinformatic analyses show that SPAP-type of APs are widely distributed in the bacterial kingdom. The role of active site residues in the catalytic mechanism has been delineated through a set of crystal structures reported here. These structures, representing different stages of the reaction pathway provide wealth of information for the catalytic mechanism. Despite critical differences in the substrate binding residues, SPAP follows a mechanism similar to that of ECAP-type of APs. Structure-based phylogenetic analysis suggests that SPAP and ECAP may have diverged very early during the evolution from a common ancestor. Moreover, it is proposed that the SPAP-type of APs are fundamental members of the AP superfamily and are more closely related to other members of the superfamily as compared to the ECAP-type of APs.

Insights into the importance of WPD-loop sequence for activity and structure in protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) possess a conserved mobile catalytic loop, the WPD-loop, which brings an aspartic acid into the active site where it acts as an acid/base catalyst. Prior experimental and computational studies, focused on the human enzyme PTP1B and the PTP from Yersinia pestis, YopH, suggested that loop conformational dynamics are important in regulating both catalysis and evolvability. We have generated a chimeric protein in which the WPD-loop of YopH is transposed into PTP1B, and eight chimeras that systematically restored the loop sequence back to native PTP1B. Of these, four chimeras were soluble and were subjected to detailed biochemical and structural characterization, and a computational analysis of their WPD-loop dynamics. The chimeras maintain backbone structural integrity, with somewhat slower rates than either wild-type parent, and show differences in the pH dependency of catalysis, and changes in the effect of Mg²⁺. The chimeric proteins' WPD-loops differ significantly in their relative stability and rigidity. The time required for interconversion, coupled with electrostatic effects revealed by simulations, likely accounts for the activity differences between chimeras, and relative to the native enzymes. Our results further the understanding of connections between enzyme activity and the dynamics of catalytically important groups, particularly the effects of non-catalytic residues on key conformational equilibria.

Structural Characterization of Functionally Important Chloride Binding Sites in the Marine Vibrio Alkaline Phosphatase

Enzyme stability and function can be affected by various environmental factors, such as temperature, pH, and ionic strength. Enzymes that are located outside the relatively unchanging environment of the cytosol, such as those residing in the periplasmic space of bacteria or extracellularly secreted, are challenged by more fluctuations in the aqueous medium. Bacterial alkaline phosphatases (APs) are generally affected by ionic strength of the medium, but this varies substantially between species. An AP from the marine bacterium Vibrio splendidus (VAP) shows complex pH-dependent activation and stabilization in the 0-1.0 M range of halogen salts and has been hypothesized to specifically bind chloride anions. Here, using X-ray crystallography and anomalous scattering, we have located two chloride binding sites in the structure of VAP, one in the active site and another one at a peripheral site. Further characterization of the binding sites using site-directed mutagenesis and small-angle X-ray scattering showed that upon binding of chloride to the peripheral site, structural dynamics decreased locally, resulting in thermal stabilization of the VAP active conformation. Binding of the chloride ion in the active site did not displace the bound inorganic phosphate product, but it may promote product release by facilitating rotational stabilization of the substrate-binding Arg129. Overall, these results reveal the complex nature and dynamics of chloride binding to enzymes through long-range modulation of electronic potential in the vicinity of the active site, resulting in increased catalytic efficiency and stability.

Structural insights into choline-O-sulfatase reveal the molecular determinants for ligand binding

The first structures of a choline-O-sulfatase bound to different ligands are reported.

Choline-O-sulfatase (COSe; EC 3.1.6.6) is a member of the alkaline phosphatase (AP) superfamily, and its natural function is to hydrolyze choline-O-sulfate into choline and sulfate. Despite its natural function, the major interest in this enzyme resides in the landmark catalytic/substrate promiscuity of sulfatases, which has led to attention in the biotechnological field due to their potential in protein engineering. In this work, an in-depth structural analysis of wild-type Sinorhizobium (Ensifer) meliloti COSe (SmeCOSe) and its C54S active-site mutant is reported. The binding mode of this AP superfamily member to both products of the reaction (sulfate and choline) and to a substrate-like compound are shown for the first time. The structures further confirm the importance of the C-terminal extension of the enzyme in becoming part of the active site and participating in enzyme activity through dynamic intra-subunit and inter-subunit hydrogen bonds (Asn146^A–Asp500^B–Asn498^B). These residues act as the ‘gatekeeper’ responsible for the open/closed conformations of the enzyme, in addition to assisting in ligand binding through the rearrangement of Leu499 (with a movement of approximately 5 Å). Trp129 and His145 clamp the quaternary ammonium moiety of choline and also connect the catalytic cleft to the C-terminus of an adjacent protomer. The structural information reported here contrasts with the proposed role of conformational dynamics in promoting the enzymatic catalytic proficiency of an enzyme.

Biomimetic Heterodimerization of Tetrapeptides to Generate Liquid Crystalline Hydrogel in A Two-Component System

Anisotropic structures made by hierarchical self-assembly and crystallization play an essential role in the living system. However, the spontaneous formation of liquid crystalline hydrogel of low molecular weight organic molecules with controlled properties remains challenging. This work describes a rational design of tetrapeptide without N-terminal modification and chemical conjugation that utilizes intermolecular interactions to drive the formation of nanofiber bundles in a two-component system, which could not be accessed by a single component. The diameter of nanofibers can be simply controlled by varying the enantiomer of electrostatic pairs. Mutation of lysine (K) to arginine (R) results in an over 30-fold increase of mechanical property. Mechanistic studies using different techniques unravel the mechanism of self-assembly and formation of anisotropic liquid crystalline domains. All-atom molecular dynamics simulations reveal that the mixture of heterochiral peptides self-assembles into a nanofiber with a larger width compared to the homochiral assemblies due to the different stacking pattern and intermolecular interactions. The intermolecular interactions show an obvious increase by substituting the K with R, facilitating a more stable assembly and further altering the assembly mechanics and bulk material properties. Moreover, we also demonstrated that the hydrogel properties can be easily controlled by incorporating a light-responsive group. This work provides a method to generate the liquid crystalline hydrogel from isotropic monomers.

Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases

Protein conformational changes can facilitate the binding of noncognate substrates and underlying promiscuous activities. However, the contribution of substrate conformational dynamics to this process is comparatively poorly understood. Here, we analyze human (hMAT2A) and Escherichia coli (eMAT) methionine adenosyltransferases that have identical active sites but different substrate specificity. In the promiscuous hMAT2A, noncognate substrates bind in a stable conformation to allow catalysis. In contrast, noncognate substrates sample stable productive binding modes less frequently in eMAT owing to altered mobility in the enzyme active site. Different cellular concentrations of substrates likely drove the evolutionary divergence of substrate specificity in these orthologues. The observation of catalytic promiscuity in hMAT2A led to the detection of a new human metabolite, methyl thioguanosine, that is produced at elevated levels in a cancer cell line. This work establishes that identical active sites can result in different substrate specificity owing to the effects of substrate and enzyme dynamics.

Loop Dynamics and Enzyme Catalysis in Protein Tyrosine Phosphatases

Protein tyrosine phosphatases (PTPs) play an important role in cellular signaling and have been implicated in human cancers, diabetes, and obesity. Despite shared catalytic mechanisms and transition states for the chemical steps of catalysis, catalytic rates within the PTP family vary over several orders of magnitude. These rate differences have been implied to arise from differing conformational dynamics of the closure of a protein loop, the WPD-loop, which carries a catalytically critical residue. The present work reports computational studies of the human protein tyrosine phosphatase 1B (PTP1B) and YopH from Yersinia pestis, for which NMR has demonstrated a link between their respective rates of WPD-loop motion and catalysis rates, which differ by an order of magnitude. We have performed detailed structural analysis, both conventional and enhanced sampling simulations of their loop dynamics, as well as empirical valence bond simulations of the chemical step of catalysis. These analyses revealed the key residues and structural features responsible for these differences, as well as the residues and pathways that facilitate allosteric communication in these enzymes. Curiously, our wild-type YopH simulations also identify a catalytically incompetent hyper-open conformation of its WPD-loop, sampled as a rare event, previously only experimentally observed in YopH-based chimeras. The effect of differences within the WPD-loop and its neighboring loops on the modulation of loop dynamics, as revealed in this work, may provide a facile means for the family of PTP enzymes to respond to environmental changes and regulate their catalytic activities.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity

Hydroxynitrile lyases (HNL's) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity.. [DOI: 10.1101/2020.04.06.027797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Hydroxynitrile lyases (HNL’s) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/ mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

Transition-State Interactions in a Promiscuous Enzyme: Sulfate and Phosphate Monoester Hydrolysis by Pseudomonas aeruginosa Arylsulfatase

Pseudomonas aeruginosa arylsulfatase (PAS) hydrolyzes sulfate and, promiscuously, phosphate monoesters. Enzyme-catalyzed sulfate transfer is crucial to a wide variety of biological processes, but detailed studies of the mechanistic contributions to its catalysis are lacking. We present linear free energy relationships (LFERs) and kinetic isotope effects (KIEs) of PAS and analyses of active site mutants that suggest a key role for leaving group (LG) stabilization. In LFERs PASWT has a much less negative Brønsted coefficient (βleaving groupobs-Enz = -0.33) than the uncatalyzed reaction (βleaving groupobs = -1.81). This situation is diminished when cationic active site groups are exchanged for alanine. The considerable degree of bond breaking during the transition state (TS) is evidenced by an 18Obridge KIE of 1.0088. LFER and KIE data for several active site mutants point to leaving group stabilization by active site K375, in cooperation with H211. 15N KIEs and the increased sensitivity to leaving group ability of the sulfatase activity in neat D2O (Δβleaving groupH-D = +0.06) suggest that the mechanism for S-Obridge bond fission shifts, with decreasing leaving group ability, from charge compensation via Lewis acid interactions toward direct proton donation. 18Ononbridge KIEs indicate that the TS for PAS-catalyzed sulfate monoester hydrolysis has a significantly more associative character compared to the uncatalyzed reaction, while PAS-catalyzed phosphate monoester hydrolysis does not show this shift. This difference in enzyme-catalyzed TSs appears to be the major factor favoring specificity toward sulfate over phosphate esters by this promiscuous hydrolase, since other features are either too similar (uncatalyzed TS) or inherently favor phosphate (charge).

Cryptic genetic variation shapes the adaptive evolutionary potential of enzymes

Genetic variation among orthologous proteins can cause cryptic phenotypic properties that only manifest in changing environments. Such variation may impact the evolvability of proteins, but the underlying molecular basis remains unclear. Here, we performed comparative directed evolution of four orthologous metallo-β-lactamases toward a new function and found that different starting genotypes evolved to distinct evolutionary outcomes. Despite a low initial fitness, one ortholog reached a significantly higher fitness plateau than its counterparts, via increasing catalytic activity. By contrast, the ortholog with the highest initial activity evolved to a less-optimal and phenotypically distinct outcome through changes in expression, oligomerization and activity. We show how cryptic molecular properties and conformational variation of active site residues in the initial genotypes cause epistasis, that could lead to distinct evolutionary outcomes. Our work highlights the importance of understanding the molecular details that connect genetic variation to protein function to improve the prediction of protein evolution.

In Silico-Directed Evolution Using CADEE

Recent years have seen an explosion of interest in both sequence- and structure-based approaches toward in silico-directed evolution. We recently developed a novel computational toolkit, CADEE, which facilitates the computer-aided directed evolution of enzymes. Our initial work (Amrein et al., IUCrJ 4:50-64, 2017) presented a pedagogical example of the application of CADEE to triosephosphate isomerase, to illustrate the CADEE workflow. In this contribution, we describe this workflow in detail, including code input/output snippets, in order to allow users to set up and execute CADEE simulations on any system of interest.

Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case

Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 "grips" the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple "open and shut" case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Epoxide hydrolysis as a model system for understanding flux through a branched reaction scheme

The epoxide hydrolase StEH1 catalyzes the hydrolysis of trans-methylstyrene oxide to 1-phenyl-propane-1,2-diol. The (S,S)-epoxide is exclusively transformed into the (1R,2S)-diol, while hydrolysis of the (R,R)-epoxide results in a mixture of product enantiomers. In order to understand the differences in the stereoconfigurations of the products, the reactions were studied kinetically during both the pre-steady-state and steady-state phases. A number of closely related StEH1 variants were analyzed in parallel, and the results were rationalized by structure-activity analysis using the available crystal structures of all tested enzyme variants. Finally, empirical valence-bond simulations were performed in order to provide additional insight into the observed kinetic behaviour and ratios of the diol product enantiomers. These combined data allow us to present a model for the flux through the catalyzed reactions. With the (R,R)-epoxide, ring opening may occur at either C atom and with similar energy barriers for hydrolysis, resulting in a mixture of diol enantiomer products. However, with the (S,S)-epoxide, although either epoxide C atom may react to form the covalent enzyme intermediate, only the pro-(R,S) alkylenzyme is amenable to subsequent hydrolysis. Previously contradictory observations from kinetics experiments as well as product ratios can therefore now be explained for this biocatalytically relevant enzyme.

Challenges and advances in the computational modeling of biological phosphate hydrolysis

Phosphate ester hydrolysis is fundamental to many life processes, and has been the topic of substantial experimental and computational research effort. However, even the simplest of phosphate esters can be hydrolyzed through multiple possible pathways that can be difficult to distinguish between, either experimentally, or computationally. Therefore, the mechanisms of both the enzymatic and non-enzymatic reactions have been historically controversial. In the present contribution, we highlight a number of technical issues involved in reliably modeling these computationally challenging reactions, as well as proposing potential solutions. We also showcase examples of our own work in this area, discussing both the non-enzymatic reaction in aqueous solution, as well as insights obtained from the computational modeling of organophosphate hydrolysis and catalytic promiscuity amongst enzymes that catalyze phosphoryl transfer.

Determinants and Prediction of Esterase Substrate Promiscuity Patterns

Esterases receive special attention because of their wide distribution in biological systems and environments and their importance for physiology and chemical synthesis. The prediction of esterases' substrate promiscuity level from sequence data and the molecular reasons why certain such enzymes are more promiscuous than others remain to be elucidated. This limits the surveillance of the sequence space for esterases potentially leading to new versatile biocatalysts and new insights into their role in cellular function. Here, we performed an extensive analysis of the substrate spectra of 145 phylogenetically and environmentally diverse microbial esterases, when tested with 96 diverse esters. We determined the primary factors shaping their substrate range by analyzing substrate range patterns in combination with structural analysis and protein-ligand simulations. We found a structural parameter that helps rank (classify) the promiscuity level of esterases from sequence data at 94% accuracy. This parameter, the active site effective volume, exemplifies the topology of the catalytic environment by measuring the active site cavity volume corrected by the relative solvent accessible surface area (SASA) of the catalytic triad. Sequences encoding esterases with active site effective volumes (cavity volume/SASA) above a threshold show greater substrate spectra, which can be further extended in combination with phylogenetic data. This measure provides also a valuable tool for interrogating substrates capable of being converted. This measure, found to be transferred to phosphatases of the haloalkanoic acid dehalogenase superfamily and possibly other enzymatic systems, represents a powerful tool for low-cost bioprospecting for esterases with broad substrate ranges, in large scale sequence data sets.

Computer simulations of the catalytic mechanism of wild-type and mutant β-phosphoglucomutase

β-Phosphoglucomutase (β-PGM) has served as an important model system for understanding biological phosphoryl transfer.

Structural plasticity mediates distinct GAP-dependent GTP hydrolysis mechanisms in Rab33 and Rab5

The classical GTP hydrolysis mechanism, as seen in Ras, employs a catalytic glutamine provided in cis by the GTPase and an arginine supplied in trans by a GTPase activating protein (GAP). The key idea emergent from a large body of research on small GTPases is that GTPases employ a variety of different hydrolysis mechanisms; evidently, these variations permit diverse rates of GTPase inactivation, crucial for temporal regulation of different biological processes. Recently, we unified these variations and argued that a steric clash between active site residues (corresponding to positions 12 and 61 of Ras) governs whether a GTPase utilizes the cis-Gln or the trans-Gln (from the GAP) for catalysis. As the cis-Gln encounters a steric clash, the Rab GTPases employ the so-called dual finger mechanism where the interacting GAP supplies a trans-Gln for catalysis. Using experimental and computational methods, we demonstrate how the cis-Gln of Rab33 overcomes the steric clash when it is stabilized by a residue in the vicinity. In effect, this demonstrates how both cis-Gln- and trans-Gln-mediated mechanisms could operate in the same GTPase in different contexts, i.e. depending on the GAP that regulates its action. Interestingly, in the case of Rab5, which possesses a higher intrinsic GTP hydrolysis rate, a similar stabilization of the cis-Gln appears to overcome the steric clash. Taken together with the mechanisms seen for Rab1, it is evident that the observed variations in Rab and their GAP partners allow structural plasticity, or in other words, the choice of different catalytic mechanisms.

Extending the Nonbonded Cationic Dummy Model to Account for Ion-Induced Dipole Interactions

Modeling metalloproteins often requires classical molecular dynamics (MD) simulations in order to capture their relevant motions, which in turn necessitates reliable descriptions of the metal centers involved. One of the most successful approaches to date is provided by the "cationic dummy model", where the positive charge of the metal ion is transferred toward dummy particles that are bonded to the central metal ion in a predefined coordination geometry. While this approach allows for ligand exchange, and captures the correct electrostatics as demonstrated for different divalent metal ions, current dummy models neglect ion-induced dipole interactions. In the present work, we resolve this weakness by taking advantage of the recently introduced 12-6-4 type Lennard-Jones potential to include ion-induced dipole interactions. We revise our previous dummy model for Mg²⁺ and demonstrate that the resulting model can simultaneously reproduce the experimental solvation free energy and metal-ligand distances without the need for artificial restraints or bonds. As ion-induced dipole interactions become particularly important for highly charged metal ions, we develop dummy models for the biologically relevant ions Al³⁺, Fe³⁺, and Cr³⁺. Finally, the effectiveness of our new models is demonstrated in MD simulations of several diverse (and highly challenging to simulate) metalloproteins.

Unraveling the Addition-Elimination Mechanism of EPSP Synthase through Computer Modeling

Enolpyruvyl transfer from phosphoenolpyruvate (PEP) to the hydroxyl group of shikimate-5-OH-3-phosphate (S3P) is catalyzed by 5-enolpyruvylshikimate 3-phosphate (EPSP) synthase in a reaction that involves breaking the C-O bond of PEP. Catalysis involves an addition-elimination mechanism with the formation of a tetrahedral intermediate (THI). Experiments have elucidated the mechanism of THI formation and breakdown. However, the catalytic action of EPSP synthase and the individual roles of catalytic residues Asp313 and Glu341 remains unclear. We have used a hybrid quantum mechanical/molecular mechanical (QM/MM) approach to explore the free energy surface in a reaction catalyzed by EPSP synthase. The Glu341 was the most favorable acid/base catalyst. Our results indicate that the protonation of PEP C3 precedes the nucleophilic attack on PEP C2 in the addition mechanism. Also, the breaking of the C-O bond of THI to form an EPSP cation intermediate must occur before proton transfer from PEP C3 to Glu341 in the elimination mechanism. Analysis of the FES supports cationic intermediate formation during the reaction catalyzed by EPSP synthase. Finally, the computational model indicates a proton transfer shift (Hammond shift) from Glu341 to C3 for an enzyme-based reaction with the shifted transition state, earlier than in the reference reaction in water.

Enzyme Architecture: Modeling the Operation of a Hydrophobic Clamp in Catalysis by Triosephosphate Isomerase

Triosephosphate isomerase (TIM) is a proficient catalyst of the reversible isomerization of dihydroxyacetone phosphate (DHAP) to d-glyceraldehyde phosphate (GAP), via general base catalysis by E165. Historically, this enzyme has been an extremely important model system for understanding the fundamentals of biological catalysis. TIM is activated through an energetically demanding conformational change, which helps position the side chains of two key hydrophobic residues (I170 and L230), over the carboxylate side chain of E165. This is critical both for creating a hydrophobic pocket for the catalytic base and for maintaining correct active site architecture. Truncation of these residues to alanine causes significant falloffs in TIM's catalytic activity, but experiments have failed to provide a full description of the action of this clamp in promoting substrate deprotonation. We perform here detailed empirical valence bond calculations of the TIM-catalyzed deprotonation of DHAP and GAP by both wild-type TIM and its I170A, L230A, and I170A/L230A mutants, obtaining exceptional quantitative agreement with experiment. Our calculations provide a linear free energy relationship, with slope 0.8, between the activation barriers and Gibbs free energies for these TIM-catalyzed reactions. We conclude that these clamping side chains minimize the Gibbs free energy for substrate deprotonation, and that the effects on reaction driving force are largely expressed at the transition state for proton transfer. Our combined analysis of previous experimental and current computational results allows us to provide an overview of the breakdown of ground-state and transition state effects in enzyme catalysis in unprecedented detail, providing a molecular description of the operation of a hydrophobic clamp in triosephosphate isomerase.

De novo active sites for resurrected Precambrian enzymes

Protein engineering studies often suggest the emergence of completely new enzyme functionalities to be highly improbable. However, enzymes likely catalysed many different reactions already in the last universal common ancestor. Mechanisms for the emergence of completely new active sites must therefore either plausibly exist or at least have existed at the primordial protein stage. Here, we use resurrected Precambrian proteins as scaffolds for protein engineering and demonstrate that a new active site can be generated through a single hydrophobic-to-ionizable amino acid replacement that generates a partially buried group with perturbed physico-chemical properties. We provide experimental and computational evidence that conformational flexibility can assist the emergence and subsequent evolution of new active sites by improving substrate and transition-state binding, through the sampling of many potentially productive conformations. Our results suggest a mechanism for the emergence of primordial enzymes and highlight the potential of ancestral reconstruction as a tool for protein engineering.

Perspective: Defining and quantifying the role of dynamics in enzyme catalysis

Enzymes control chemical reactions that are key to life processes, and allow them to take place on the time scale needed for synchronization between the relevant reaction cycles. In addition to general interest in their biological roles, these proteins present a fundamental scientific puzzle, since the origin of their tremendous catalytic power is still unclear. While many different hypotheses have been put forward to rationalize this, one of the proposals that has become particularly popular in recent years is the idea that dynamical effects contribute to catalysis. Here, we present a critical review of the dynamical idea, considering all reasonable definitions of what does and does not qualify as a dynamical effect. We demonstrate that no dynamical effect (according to these definitions) has ever been experimentally shown to contribute to catalysis. Furthermore, the existence of non-negligible dynamical contributions to catalysis is not supported by consistent theoretical studies. Our review is aimed, in part, at readers with a background in chemical physics and biophysics, and illustrates that despite a substantial body of experimental effort, there has not yet been any study that consistently established a connection between an enzyme's conformational dynamics and a significant increase in the catalytic contribution of the chemical step. We also make the point that the dynamical proposal is not a semantic issue but a well-defined scientific hypothesis with well-defined conclusions.

Active Site Hydrophobicity and the Convergent Evolution of Paraoxonase Activity in Structurally Divergent Enzymes: The Case of Serum Paraoxonase 1

Serum paraoxonase 1 (PON1) is a native lactonase capable of promiscuously hydrolyzing a broad range of substrates, including organophosphates, esters, and carbonates. Structurally, PON1 is a six-bladed β-propeller with a flexible loop (residues 70–81) covering the active site. This loop contains a functionally critical Tyr at position 71. We have performed detailed experimental and computational analyses of the role of selected Y71 variants in the active site stability and catalytic activity in order to probe the role of Y71 in PON1’s lactonase and organophosphatase activities. We demonstrate that the impact of Y71 substitutions on PON1’s lactonase activity is minimal, whereas the k_cat for the paraoxonase activity is negatively perturbed by up to 100-fold, suggesting greater mutational robustness of the native activity. Additionally, while these substitutions modulate PON1’s active site shape, volume, and loop flexibility, their largest effect is in altering the solvent accessibility of the active site by expanding the active site volume, allowing additional water molecules to enter. This effect is markedly more pronounced in the organophosphatase activity than the lactonase activity. Finally, a detailed comparison of PON1 to other organophosphatases demonstrates that either a similar “gating loop” or a highly buried solvent-excluding active site is a common feature of these enzymes. We therefore posit that modulating the active site hydrophobicity is a key element in facilitating the evolution of organophosphatase activity. This provides a concrete feature that can be utilized in the rational design of next-generation organophosphate hydrolases that are capable of selecting a specific reaction from a pool of viable substrates.

CADEE: Computer-Aided Directed Evolution of Enzymes

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

Complete sampling of an enzyme reaction pathway: a lesson from gas phase simulations

With proper sampling strategy, convergence of free energy profiles of biomolecular reactions in the gas phase can be achieved in microseconds of simulation.

Probing the mechanisms for the selectivity and promiscuity of methyl parathion hydrolase

Diverse organophosphate hydrolases have convergently evolved the ability to hydrolyse man-made organophosphates. Thus, these enzymes are attractive model systems for studying the factors shaping enzyme functional evolution. Methyl parathion hydrolase (MPH) is an enzyme from the metallo-β-lactamase superfamily, which hydrolyses a wide range of organophosphate, aryl ester and lactone substrates. In addition, MPH demonstrates metal-ion-dependent selectivity patterns. The origins of this remain unclear, but are linked to open questions about the more general role of metal ions in functional evolution and divergence within enzyme superfamilies. Here, we present detailed mechanistic studies of the paraoxonase and arylesterase activities of MPH complexed with five different transition metal ions, and demonstrate that the hydrolysis reactions proceed via similar pathways and transition states. However, while it is possible to discern a clear structural origin for the selectivity between different substrates, the selectivity between different metal ions appears to lie instead in the distinct electrostatic properties of the metal ions themselves, which causes subtle changes in transition state geometries and metal-metal distances at the transition state rather than significant structural changes in the active site. While subtle, these differences can be significant for shaping the metal-ion-dependent activity patterns observed for this enzyme.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Catalytic scaffolds for phosphoryl group transfer

A single genome encodes a large number of phosphoryl hydrolases for the purposes of phosphate recycling, primary and secondary metabolism, signal transduction and regulation, and protection from xenobiotics. Phosphate monoester hydrolysis faces a high kinetic barrier, yet there are multiple solutions to the problem both in terms of catalytic mechanisms and three-dimensional structure of the hydrolases. Recent structural and mechanistic findings highlight the trigonal-bipyramidal nature of the transition state for enzyme promoted phosphate monoester hydrolysis and the evolution and role of inserted loops/domains in governing substrate specificity and promiscuity. Important questions remain as to how electrostatics modulate water networks and critical proton-transfer events. How substrate targeting and catalysis is achieved by the independently evolved catalytic platforms is compared and contrasted in this article.

Structural snapshots of the catalytic cycle of the phosphodiesterase Autotaxin

Autotaxin (ATX) is a secreted phosphodiesterase that produces the signalling lipid lysophosphatidic acid (LPA). The bimetallic active site of ATX is structurally related to the alkaline phosphatase superfamily. Here, we present a new crystal structure of ATX in complex with orthovanadate (ATX-VO5), which binds the Oγ nucleophile of Thr209 and adopts a trigonal bipyramidal conformation, following the nucleophile attack onto the substrate. We have now a portfolio of ATX structures we discuss as intermediates of the catalytic mechanism: the new ATX-VO5 structure; a unique structure where the nucleophile Thr209 is phosphorylated (ATX-pThr). Comparing these to a complex with the LPA product (ATX-LPA) and with a complex with a phosphate ion (ATX-PO4), that represent the Michaelis complex of the reaction, we observe movements of Thr209, changes in the relative displacement of the zinc ions, and a water molecule that likely fulfils the second nucleophilic attack. We propose that ATX follows the associative two-step in-line displacement mechanism.

Leaving Group Ability Observably Affects Transition State Structure in a Single Enzyme Active Site

A reaction's transition state (TS) structure plays a critical role in determining reactivity and has important implications for the design of catalysts, drugs, and other applications. Here, we explore TS structure in the enzyme alkaline phosphatase using hybrid Quantum Mechanics/Molecular Mechanics simulations. We find that minor perturbations to the substrate have major effects on TS structure and the way the enzyme stabilizes the TS. Substrates with good leaving groups (LGs) have little cleavage of the phosphorus-LG bond at the TS, while substrates with poor LGs have substantial cleavage of that bond. The results predict nonlinear free energy relationships for a single rate-determining step, and substantial differences in kinetic isotope effects for different substrates; both trends were observed in previous experimental studies, although the original interpretations differed from the present model. Moreover, due to different degrees of phosphorus-LG bond cleavage at the TS for different substrates, the LG is stabilized by different interactions at the TS: while a poor LG is directly stabilized by an active site zinc ion, a good LG is mainly stabilized by active site water molecules. Our results demonstrate the considerable plasticity of TS structure and stabilization in enzymes. Furthermore, perturbations to reactivity that probe TS structure experimentally (i.e., substituent effects) may substantially perturb the TS they aim to probe, and thus classical experimental approaches such as free energy relations should be interpreted with care.

Promiscuity in the Enzymatic Catalysis of Phosphate and Sulfate Transfer

The enzymes that facilitate phosphate and sulfate hydrolysis are among the most proficient natural catalysts known to date. Interestingly, a large number of these enzymes are promiscuous catalysts that exhibit both phosphatase and sulfatase activities in the same active site and, on top of that, have also been demonstrated to efficiently catalyze the hydrolysis of other additional substrates with varying degrees of efficiency. Understanding the factors that underlie such multifunctionality is crucial both for understanding functional evolution in enzyme superfamilies and for the development of artificial enzymes. In this Current Topic, we have primarily focused on the structural and mechanistic basis for catalytic promiscuity among enzymes that facilitate both phosphoryl and sulfuryl transfer in the same active site, while comparing this to how catalytic promiscuity manifests in other promiscuous phosphatases. We have also drawn on the large number of experimental and computational studies of selected model systems in the literature to explore the different features driving the catalytic promiscuity of such enzymes. Finally, on the basis of this comparative analysis, we probe the plausible origins and determinants of catalytic promiscuity in enzymes that catalyze phosphoryl and sulfuryl transfer.

Insights into enzyme point mutation effect by molecular simulation: phenylethylamine oxidation catalyzed by monoamine oxidase A

The I335Y point mutation effect on the kinetics of phenylethylamine decomposition catalyzed by monoamine oxidase A was elucidated by means of molecular simulation. The established empirical valence bond methodology was used in conjunction with the free energy perturbation sampling technique and a classical force field representing the state of reactants and products. The methodology allows for the simulation of chemical reactions, in the present case the breaking of the α-C-H bond in a phenylethylamine substrate and the subsequent hydrogen transfer to the flavin cofactor, resulting in the formation of the N-H bond on flavin. The empirical parameters were calibrated against the experimental data for the simulated reaction in a wild type protein and then used for the calculation of the reaction free energy profile in the I335Y mutant. In very good agreement with the measured kinetic data, mutation increases the free energy barrier for the rate limiting step by slightly more than 1 kcal mol(-1) and consequently decreases the rate constant by about an order of magnitude. The magnitude of the computed effect slightly varies with simulation settings, but always remains in reasonable agreement with the experiment. Analysis of trajectories reveals a major change in the interaction between phenyl rings of the substrate and the neighboring Phe352 residue upon the I335Y mutation due to the increased local polarity, leading to an attenuated quadrupole interaction between the rings and destabilization of the transition state. Additionally, the increased local polarity in the mutant allows for a larger number of water molecules to be present near the active site, effectively shielding the catalytic effect of the enzyme and contributing to the increased barrier.

Catalytic mechanism of acetolactate decarboxylase from Brevibacillus brevis towards both enantiomers of α-acetolactate

QM/MM calculations detail the conversion of both enantiomers of α-acetolactate to (R)-acetoin. (S)-α-acetolactate can be directly converted, whereas (R)-α-acetolactate has to firstly rearrange to the natural (S)-enantiomer through a carboxylate migration.

Promiscuity and electrostatic flexibility in the alkaline phosphatase superfamily

Catalytic promiscuity, that is, the ability of single enzymes to facilitate the turnover of multiple, chemically distinct substrates, is a widespread phenomenon that plays an important role in the evolution of enzyme function. Additionally, such pre-existing multifunctionality can be harnessed in artificial enzyme design. The members of the alkaline phosphatase superfamily have served extensively as both experimental and computational model systems for enhancing our understanding of catalytic promiscuity. In this Opinion, we present key recent computational studies into the catalytic activity of these highly promiscuous enzymes, highlighting the valuable insight they have provided into both the molecular basis for catalytic promiscuity in general, and its implications for the evolution of phosphatase activity.

Linking coupled motions and entropic effects to the catalytic activity of 2-deoxyribose-5-phosphate aldolase (DERA)

Local mutations in the phosphate binding group of DERA alter global conformation dynamics, catalytic activities and reaction entropies.

DERA, 2-deoxyribose-5-phosphate aldolase, catalyzes the retro-aldol cleavage of 2-deoxy-ribose-5-phosphate (dR5P) into glyceraldehyde-3-phosphate (G3P) and acetaldehyde in a branch of the pentose phosphate pathway. In addition to the physiological reaction, DERA also catalyzes the reverse addition reaction and, hence, is an interesting candidate for bio-catalysis of carbo-ligation reactions, which are central to synthetic chemistry. An obstacle to overcome for this enzyme to become a truly useful biocatalyst, however, is to relax the very strict dependency of this enzyme on phosphorylated substrates. We have studied herein the role of the non-canonical phosphate-binding site of this enzyme, consisting of Ser²³⁸ and Ser²³⁹, by site-directed and site-saturation mutagenesis, coupled to kinetic analysis of mutants. In addition, we have performed molecular dynamics simulations on the wild-type and four mutant enzymes, to analyse how mutations at this phosphate-binding site may affect the protein structure and dynamics. Further examination of the S239P mutant revealed that this variant increases the enthalpy change at the transition state, relative to the wild-type enzyme, but concomitant loss in entropy causes an overall relative loss in the TS free energy change. This entropy loss, as measured by the temperature dependence of catalysed rates, was mirrored in both a drastic loss in dynamics of the enzyme, which contributes to phosphate binding, as well as an overall loss in anti-correlated motions distributed over the entire protein. Our combined data suggests that the degree of anticorrelated motions within the DERA structure is coupled to catalytic efficiency in the DERA-catalyzed retro-aldol cleavage reaction, and can be manipulated for engineering purposes.