Transition State Binding Research Articles

How Robust Is the Ligand Binding Transition State?

For many drug targets, it has been shown that the kinetics of drug binding (e.g., on rate and off rate) is more predictive of drug efficacy than thermodynamic quantities alone. This motivates the development of predictive computational models that can be used to optimize compounds on the basis of their kinetics. The structural details underpinning these computational models are found not only in the bound state but also in the short-lived ligand binding transition states. Although transition states cannot be directly observed experimentally due to their extremely short lifetimes, recent successes have demonstrated that modeling the ligand binding transition state is possible with the help of enhanced sampling molecular dynamics methods. Previously, we generated unbinding paths for an inhibitor of soluble epoxide hydrolase (sEH) with a residence time of 11 min. Here, we computationally modeled unbinding events with the weighted ensemble method REVO (resampling of ensembles by variation optimization) for five additional inhibitors of sEH with residence times ranging from 14.25 to 31.75 min, with average prediction accuracy within an order of magnitude. The unbinding ensembles are analyzed in detail, focusing on features of the ligand binding transition state ensembles (TSEs). We find that ligands with similar bound poses can show significant differences in their ligand binding TSEs, in terms of their spatial distribution and protein-ligand interactions. However, we also find similarities across the TSEs when examining more general features such as ligand degrees of freedom. Together these findings show significant challenges for rational, kinetics-based drug design.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity.

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs—common to all P-loop NTPases—could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Creating Maps of the Ligand Binding Landscape for Kinetics-Based Drug Discovery.

Simulations of ligand-protein interactions can be very useful for drug design and to gain biological insight. Full pathways of ligand-protein binding can be used to get information about ligand binding transition states, which form the rate-limiting step of the binding and release processes. However, these simulations are typically limited by the presence of large energy barriers that separate stable poses of interest. Here we describe a simulation protocol for exploring and analyzing landscapes of ligand-protein interactions that makes use of molecular docking, enhanced molecular simulation with the weighted ensemble algorithm, and network analysis. It can be accomplished using a modest cluster of graphics processing units and freely accessible software. This protocol focuses on the construction and analysis of a network model of ligand binding poses and provides links to resources that describe the other steps in more detail. The end result of this protocol is a map of the ligand-protein binding landscape that identifies transition states of the ligand binding pathway, as well as alternative bound poses that could be stabilized with modifications to the ligand.

Understanding the Binding Transition State After the Conformational Selection Step: The Second Half of the Molecular Recognition Process Between NS1 of the 1918 Influenza Virus and Host p85β.

Biomolecular recognition often involves conformational changes as a prerequisite for binding (i.e., conformational selection) or concurrently with binding (i.e., induced-fit). Recent advances in structural and kinetic approaches have enabled the detailed characterization of protein motions at atomic resolution. However, to fully understand the role of the conformational dynamics in molecular recognition, studies on the binding transition state are needed. Here, we investigate the binding transition state between nonstructural protein 1 (NS1) of the pandemic 1918 influenza A virus and the human p85β subunit of PI3K. 1918 NS1 binds to p85β via conformational selection. We present the free-energy mapping of the transition and bound states of the 1918 NS1:p85β interaction using linear free energy relationship and ϕ-value analyses. We find that the binding transition state of 1918 NS1 and p85β is structurally similar to the bound state with well-defined binding orientation and hydrophobic interactions. Our finding provides a detailed view of how protein motion contributes to the development of intermolecular interactions along the binding reaction coordinate.

Membrane-Mediated Ligand Unbinding of the PK-11195 Ligand from TSPO

The translocator protein (TSPO), previously known as the peripheral benzodiazepine receptor, is of longstanding medical interest as both a biomarker for neuroinjury and a potential drug target for neuroinflammation and other disorders. Recently, it was shown that ligand residence time is a key factor determining steroidogenic efficacy of TSPO-binding compounds. This spurs interest in simulations of (un)binding pathways of TSPO ligands, which could reveal the molecular interactions governing ligand residence time. In this study, we use a weighted ensemble algorithm to determine the unbinding pathway for different poses of PK-11195, a TSPO ligand used in neuroimaging. In contrast with previous studies, our results show that PK-11195 does not dissociate directly into the solvent but instead dissociates via the lipid membrane by going between the transmembrane helices. We analyze this path ensemble in detail, constructing descriptors that can facilitate a general understanding of membrane-mediated ligand binding. We construct a set of Markov state models augmented with additional straightforward simulations to determine pose-specific ligand residence times. Together, we combine over 40 μs of trajectory data to form a coherent picture of the ligand binding landscape. We find that multiple starting poses yield residence times that roughly agree with the experimental quantity. The ligand binding transition states predicted by these Markov state models occur when PK-11195 is already in the membrane and involves only minimal ligand-protein interactions. This has implications for the design of new long-residence-time TSPO ligands.

Kinetics and thermodynamics of catalysis and thermal inactivation of a novel α-amylase (Tp-AmyS) from Thermotoga petrophila

This research is focussed on kinetic, thermodynamic and thermal inactivation of a novel thermostable recombinant α-amylase (Tp-AmyS) from Thermotoga petrophila. The amylase gene was cloned in pHIS-parallel1 expression vector and overexpressed in Escherichia coli. The steady-state kinetic parameters (Vmax, Km, kcat and kcat/Km) for the hydrolysis of amylose (1.39 mg/min, 0.57 mg, 148.6 s−1, 260.7), amylopectin (2.3 mg/min, 1.09 mg, 247.1 s−1, 226.7), soluble starch (2.67 mg/min, 2.98 mg, 284.2 s−1, 95.4) and raw starch (2.1 mg/min, 3.6 mg, 224.7 s−1, 61.9) were determined. The activation energy (Ea), free energy (ΔG), enthalpy (ΔH) and entropy of activation (ΔS) at 98 °C were 42.9 kJ mol−1, 74 kJ mol−1, 39.9 kJ mol−1 and −92.3 J mol−1 K−1, respectively, for soluble starch hydrolysis. While ΔG of substrate binding (ΔGE-S) and ΔG of transition state binding (ΔGE-T) were 3.38 and −14.1 kJ mol−1, respectively. Whereas, EaD, Gibbs free energy (ΔG*), increase in the enthalpy (ΔH*) and activation entropy (ΔS*) for activation of the unfolding of transition state were 108, 107, 105 kJ mol−1 and −4.1 J mol−1 K−1. The thermodynamics of irreversible thermal inactivation of Tp-AmyS revealed that at high temperature the process involves the aggregation of the protein.

Adaptive landscape flattening allows the design of both enzyme: Substrate binding and catalytic power.

Designed enzymes are of fundamental and technological interest. Experimental directed evolution still has significant limitations, and computational approaches are a complementary route. A designed enzyme should satisfy multiple criteria: stability, substrate binding, transition state binding. Such multi-objective design is computationally challenging. Two recent studies used adaptive importance sampling Monte Carlo to redesign proteins for ligand binding. By first flattening the energy landscape of the apo protein, they obtained positive design for the bound state and negative design for the unbound. We have now extended the method to design an enzyme for specific transition state binding, i.e., for its catalytic power. We considered methionyl-tRNA synthetase (MetRS), which attaches methionine (Met) to its cognate tRNA, establishing codon identity. Previously, MetRS and other synthetases have been redesigned by experimental directed evolution to accept noncanonical amino acids as substrates, leading to genetic code expansion. Here, we have redesigned MetRS computationally to bind several ligands: the Met analog azidonorleucine, methionyl-adenylate (MetAMP), and the activated ligands that form the transition state for MetAMP production. Enzyme mutants known to have azidonorleucine activity were recovered by the design calculations, and 17 mutants predicted to bind MetAMP were characterized experimentally and all found to be active. Mutants predicted to have low activation free energies for MetAMP production were found to be active and the predicted reaction rates agreed well with the experimental values. We suggest the present method should become the paradigm for computational enzyme design.

Protein Flexibility and Stiffness Enable Efficient Enzymatic Catalysis.

The enormous rate accelerations observed for many enzyme catalysts are due to strong stabilizing interactions between the protein and reaction transition state. The defining property of these catalysts is their specificity for binding the transition state with a much higher affinity than substrate. Experimental results are presented which show that the phosphodianion-binding energy of phosphate monoester substrates is used to drive conversion of their protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes. These results are generalized to other enzyme-catalyzed reactions. The existence of many enzymes in flexible, entropically rich, and inactive ground states provides a mechanism for utilization of ligand-binding energy to mold these catalysts into stiff and active forms. This reduces the substrate-binding energy expressed at the Michaelis complex, while enabling the full and specific expression of large transition-state binding energies. Evidence is presented that the complexity of enzyme conformational changes increases with increases in the enzymatic rate acceleration. The requirement that a large fraction of the total substrate-binding energy be utilized to drive conformational changes of floppy enzymes is proposed to favor the selection and evolution of protein folds with multiple flexible unstructured loops, such as the TIM-barrel fold. The effect of protein motions on the kinetic parameters for enzymes that undergo ligand-driven conformational changes is considered. The results of computational studies to model the complex ligand-driven conformational change in catalysis by triosephosphate isomerase are presented.

New Irreversible α-l-Iduronidase Inhibitors and Activity-Based Probes.

Cyclophellitol aziridines are potent irreversible inhibitors of retaining glycosidases and versatile intermediates in the synthesis of activity‐based glycosidase probes (ABPs). Direct 3‐amino‐2‐(trifluoromethyl)quinazolin‐4(3H)‐one‐mediated aziridination of l‐ido‐configured cyclohexene has enabled the synthesis of new covalent inhibitors and ABPs of α‐l‐iduronidase, deficiency of which underlies the lysosomal storage disorder mucopolysaccharidosis type I (MPS I). The iduronidase ABPs react covalently and irreversibly in an activity‐based manner with human recombinant α‐l‐iduronidase (rIDUA, Aldurazyme®). The structures of IDUA when complexed with the inhibitors in a non‐covalent transition state mimicking form and a covalent enzyme‐bound form provide insights into its conformational itinerary. Inhibitors 1–3 adopt a half‐chair conformation in solution (4H3 and 3H4), as predicted by DFT calculations, which is different from the conformation of the Michaelis complex observed by crystallographic studies. Consequently, 1–3 may need to overcome an energy barrier in order to switch from the 4H3 conformation to the transition state (2, 5B) binding conformation before reacting and adopting a covalent 5S1 conformation. rIDUA can be labeled with fluorescent Cy5 ABP 2, which allows monitoring of the delivery of therapeutic recombinant enzyme to lysosomes, as is intended in enzyme replacement therapy for the treatment of MPS I patients.

Conjugation of Agrobacterium radiobacter epoxide hydrolase with ficoll: Catalytic, kinetic and thermodynamic analysis

Epoxide hydrolase-mediated biocatalysis has a great prospective in the biosynthesis of optically pure epoxides. The present work targets toward the thermo-stabilization of epoxide hydrolase by covalent conjugation with polysaccharide. An epoxide hydrolase from Agrobacterium radiobacter (ArEH) was modified by covalent coupling to seven oxidized polysaccharides with different chemical structures and characteristics. Among all conjugates, ArEH with ficoll exhibited both the highest specific activity (494 U/mg) and half-life at 60 °C (t1/2 = 183 min). The conjugated enzyme also displayed wider optimum pH and temperature ranges, higher catalytic number (kcat), and catalytic efficiency (kcat/Km) as compared to its native counterpart. The conjugation significantly decreased the enthalpy of activation (ΔH*), free energy of transition state binding (ΔG*E−T), and free energy of activation (ΔG*) for epichlorohydrin hydrolysis. Moreover, the conjugated ArEH showed enhanced thermal stability as evidenced by its longer half-life (t1/2), lower thermal deactivation constant (kd), and higher D values at 50–70 °C. Furthermore, the enzymatic activity of conjugated ArEH elevated significantly by Ca2+ and it showed increased tolerance against Co2+, Fe3+ and EDTA inhibitors. Finally, fluorescence and circular dichroism spectroscopic analysis were performed to confirm the noticeable conformational changes in ArEH structure after conjugation with ficoll, which could be responsible for the observed catalytic and stability enhancement.

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.

The transition state structure for binding between TAZ1 of CBP and the disordered Hif-1\u03b1 CAD

Intrinsically disordered proteins (IDPs) are common in eukaryotes. However, relatively few experimental studies have addressed the nature of the rate-limiting transition state for the coupled binding and folding reactions involving IDPs. By using site-directed mutagenesis in combination with kinetics measurements we have here characterized the transition state for binding between the globular TAZ1 domain of CREB binding protein and the intrinsically disordered C-terminal activation domain of Hif-1α (Hif-1α CAD). A total of 17 Hif-1α CAD point-mutations were generated and a Φ-value binding analysis was carried out. We found that native hydrophobic binding interactions are not formed at the transition state. We also investigated the effect the biologically important Hif-1α CAD Asn-803 hydroxylation has on the binding kinetics, and found that the whole destabilization effect due the hydroxylation is within the dissociation rate constant. Thus, the rate-limiting transition state is “disordered-like”, with native hydrophobic binding contacts being formed cooperatively after the rate-limiting barrier, which is clearly shown by linear free energy relationships. The same behavior was observed in a previously characterized TAZ1/IDP interaction, which may suggest common features for the rate-limiting transition state for TAZ1/IDP interactions.

Active subsite properties, subsite residues and targeting to lysosomes or midgut lumen of cathepsins L from the beetle Tenebrio molitor

Cathepsins L are the major digestive peptidases in the beetle Tenebrio molitor. Two digestive cathepsins L (TmCAL2 and TmCAL3) from it had their 3D structures solved. The aim of this paper was to study in details TmCAL3 specificity and properties and relate them to its 3D structure. Recombinant TmCAL3 was assayed with 64 oligopeptides with different amino acid replacements in positions P2, P1, P1′ and P2'. Results showed that TmCAL3 S2 specificity differs from the human enzyme and that its specificities also explain why on autoactivation two propeptide residues remain in the enzyme. Data on free energy of binding and of activation showed that S1 and S2′ are mainly involved in substrate binding, S1′ acts in substrate binding and catalysis, whereas S2 is implied mainly in catalysis. Enzyme subsite residues were identified by docking with the same oligopeptide used for kinetics. The subsite hydrophobicities were calculated from the efficiency of hydrolysis of different amino acid replacements in the peptide and from docking data. The results were closer for S1 and S2′ than for S1' and S2, indicating that the residue subsites that were more involved in transition state binding are different from those binding the substrate seen in docking. Besides TmCAL1-3, there are nine other cathepsins L, most of them more expressed at midgut. They are supposed to be directed to lysosomes by a Drosophila-like Lerp receptor and/or motifs in their prodomains. The mannose 6-phosphate lysosomal sorting machinery is absent from T. molitor transcriptome. Cathepsin L direction to midgut contents seems to depend on overexpression.

Native Hydrophobic Binding Interactions at the Transition State for Association between the TAZ1 Domain of CBP and the Disordered TAD-STAT2 Are Not a Requirement.

A significant fraction of the eukaryotic proteome consists of proteins that are either partially or completely disordered under native-like conditions. Intrinsically disordered proteins (IDPs) are common in protein-protein interactions and are involved in numerous cellular processes. Although many proteins have been identified as disordered, much less is known about the binding mechanisms of the coupled binding and folding reactions involving IDPs. Here we have analyzed the rate-limiting transition state for binding between the TAZ1 domain of CREB binding protein and the intrinsically disordered transactivation domain of STAT2 (TAD-STAT2) by site-directed mutagenesis and kinetic experiments (Φ-value analysis) and found that the native protein-protein binding interface is not formed at the transition state for binding. Instead, native hydrophobic binding interactions form late, after the rate-limiting barrier has been crossed. The association rate constant in the absence of electrostatic enhancement was determined to be rather high. This is consistent with the Φ-value analysis, which showed that there are few or no obligatory native contacts. Also, linear free energy relationships clearly demonstrate that native interactions are cooperatively formed, a scenario that has usually been observed for proteins that fold according to the so-called nucleation-condensation mechanism. Thus, native hydrophobic binding interactions at the rate-limiting transition state for association between TAD-STAT2 and TAZ1 are not a requirement, which is generally in agreement with previous findings on other IDP systems and might be a common mechanism for IDPs.

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.

Ordered yttrium concentration effects on barium zirconate structure, proton binding sites and transition states

Acceptor doped perovskites show promise as fuel cell proton conducting ceramic membranes and continue to be rapidly developed. Proton conduction in barium zirconate with 6.25% ordered yttrium doping at the zirconium site is explored. The dopant increases octahedral distortions in the lattice similar to our earlier studied 12.5% dopant system though the angle distribution is slightly broader. Normal modes promoting vibrations about a non-transient octahedral tilt are found for both concentrations at frequencies thermally accessible at typical proton conduction temperatures. The energy of proton binding sites and transition states roughly increases with distance away from the dopant though there are many oscillations in the trend due to local proton induced distortions. Protons can both locally enhance and decrease the octahedral distortions induced by the dopant. Transition state energies roughly increase with proton-dopant distance. Rotational and intraoctahedral transfer transition state energies vary through a similar magnitude range though the intraoctahedral transition states are often higher in energy. Interoctahedral transfer transition states tend to have the highest energies. The diversity of minima and transition states suggests that proton conduction in 6.25% yttrium doped barium zirconate is more complex than conduction in 12.5% doped barium zirconate.

Evaluating the Contribution of Transition-State Destabilization to Changes in the Residence Time of Triazole-Based InhA Inhibitors.

A critical goal of lead compound selection and optimization is to maximize target engagement while minimizing off-target binding. Since target engagement is a function of both the thermodynamics and kinetics of drug-target interactions, it follows that the structures of both the ground states and transition states on the binding reaction coordinate are needed to rationally modulate the lifetime of the drug-target complex. Previously, we predicted the structure of the rate-limiting transition state that controlled the time-dependent inhibition of the enoyl-ACP reductase InhA. This led to the discovery of a triazole-containing diphenyl ether with an increased residence time on InhA due to transition-state destabilization rather than ground-state stabilization. In the present work, we evaluate the inhibition of InhA by 14 triazole-based diphenyl ethers and use a combination of enzyme kinetics and X-ray crystallography to generate a structure-kinetic relationship for time-dependent binding. We show that the triazole motif slows the rate of formation for the final drug-target complex by up to 3 orders of magnitude. In addition, we identify a novel inhibitor with a residence time on InhA of 220 min, which is 3.5-fold longer than that of the INH-NAD adduct formed by the tuberculosis drug, isoniazid. This study provides a clear example in which the lifetime of the drug-target complex is controlled by interactions in the transition state for inhibitor binding rather than the ground state of the enzyme-inhibitor complex, and demonstrates the important role that on-rates can play in drug-target residence time.

Kinetic Insights into the Binding between the nSH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The interaction between CrkII and cAbl is implicated in diverse cellular processes. This interaction starts with the binding of the N-terminal Src homology 3 (nSH3) domain of CrkII to the proline-rich motifs of cAbl (PRMscAbl). Despite its critical importance, the detailed binding mechanism between the nSH3 domain and PRMs remains elusive. In this study, we used nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiment to study the binding kinetics between the nSH3 domain of CrkII and PRMscAbl. Our results highlight that the nSH3 domain binds to three PRMscAbl with very high on- and off-rate constants, indicating the transient nature of the binding. To further characterize the binding transition state, we conducted the Eyring and linear free energy relationship analyses using temperature-dependent kinetic data. These data indicate that the binding transition state of the nSH3 domain and PRM is accompanied by small activation enthalpy, owing to partial desolvation of the transition state. These results also highlight the similarity between the transition and free states, in terms of structure and energetics. Although the binding of the nSH3 domain and PRM displays the features consistent with a diffusion-limited process within our experimental conditions, further tests are necessary to determine if the binding is a true diffusion-limited process.

Methanol Partial Oxidation on Ag(1 1 1) from First Principles

AbstractIn this work, we examine the thermochemistry and kinetics of the partial oxidation of methanol to formaldehyde on silver surfaces. Periodic density functional theory calculations employing the BEEF‐vdW functional are used to identify the most stable phases of the silver surface under relevant reaction conditions and the reaction energetics are obtained on these surfaces. The calculated binding energies and transition state energies are used as input in a mean‐field microkinetic model providing the reaction kinetics on silver surfaces under different reaction conditions. Our results show that, under conditions pertaining to methanol partial oxidation, oxygen is present at low concentrations and it plays a critical role in the catalytic reaction. Surface oxygen promotes the reaction by activating the OH bond in methanol, thus forming a methoxy intermediate, which can react further to form formaldehyde. The dissociation of molecular oxygen is identified as the most critical step.

Substrate and Transition State Binding in Alkaline Phosphatase Analyzed by Computation of Oxygen Isotope Effects.

Enzymes are powerful catalysts, and a thorough understanding of the sources of their catalytic power will facilitate many medical and industrial applications. Here we have studied the catalytic mechanism of alkaline phosphatase (AP), which is one of the most catalytically proficient enzymes known. We have used quantum mechanics calculations and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to model a variety of isotope effects relevant to the reaction of AP. We have calculated equilibrium isotope effects (EIEs), binding isotope effects (BIEs), and kinetic isotope effects (KIEs) for a range of phosphate mono- and diester substrates. The results agree well with experimental values, but the model for the reaction's transition state (TS) differs from the original interpretation of those experiments. Our model indicates that isotope effects on binding make important contributions to measured KIEs on V/K, which complicated interpretation of the measured values. Our results provide a detailed interpretation of the measured isotope effects and make predictions that can test the proposed model. The model indicates that the substrate is deformed in the ground state (GS) of the reaction and partially resembles the TS. The highly preorganized active site preferentially binds conformations that resemble the TS and not the GS, which induces the substrate to adapt to the enzyme, rather than the other way around-as with classic "induced fit" models. The preferential stabilization of the TS over the GS is what lowers the barrier to the chemical step.