Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution

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.

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.

Hydrophobicity-Based Force Field In Enzymes

The biocatalysis process takes place with the participation of enzymes, which, depending on the reaction carried out, require, apart from the appropriate arrangement of catalytic residues, an appropriate external force field. It is generated by the protein body. The relatively small size of the part directly involved in the process itself is supported by the presence of an often complex structure of the protein body, the purpose of which is to provide an appropriate local force field, eliminating the influence of water. Very often, the large size of the enzyme is an expression of the complex form of this field. In this paper, a comparative analysis of arbitrarily selected enzymes, representatives of different enzyme classes, was carried out, focusing on the measurement of the diversity of the force field provided by a given protein. This analysis was based on the fuzzy oil drop model (FOD) and its modified version (FOD-M), which takes into account the participation of nonaqueous external factors in shaping the structure and thus the force field within the protein. The degree and type of ordering of the hydrophobicity distribution in the protein molecule is the result of the influence of the environment but also the supplier of the local environment for a given process, including the catalysis process in particular. Determining the share of a nonaqueous environment is important due to the ubiquity of polar water, whose participation in processes with high specificity requires control. It can be assumed that some enzymes in their composition have a permanently built-in part, the role of which is reduced to that of a permanent chaperone. It provides a specific external force field needed for the process. The proposed model, generalized to other types of proteins, may also provide a form of recording the environment model for the simulation of the in silico protein folding process, taking into account the impact of its differentiation.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Dynamics of receptor activation by agonists

How do agonists turn on receptors? The model system we have used to address this question is the adult-type skeletal muscle nicotinic acetylcholine receptor. This ligand-gated ion channel has two orthosteric sites (for neurotransmitters) in the extracellular domain linked to an allosteric site (a gate) in the transmembrane domain. The goal of this perspective is to summarize how measurements of agonist binding energy reveal the dynamics of the neurotransmitter sites and the fundamental link between binding and gating.

Decreased Transition-State Analogue Affinity in Isotopically Heavy MTAN with Increased Catalysis

5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase from Helicobacter pylori (HpMTAN) demonstrated faster chemistry when expressed as an isotopically heavy protein, with 2H, 13C, and 15N replacing the bulk of normal isotopes. The inverse heavy enzyme isotope effect has been attributed to improved enzyme-reactant interactions causing more frequent transition-state formation ( Proc. Natl. Acad. Sci. U.S.A. 2021, 118, e2109118118). Transition-state analogues stabilize the transient dynamic geometry of the transition state and inform on transition-state dynamics. Here, a slow-onset, tight-binding transition-state analogue of HpMTAN is characterized with heavy and light enzymes. Dissociation constants for the initial encounter complex (Ki) and for the tightly bound complex after slow-onset inhibition (Ki*) with hexylthio-DADMe-Immucillin-A (HTDIA) gave Ki values for light and heavy HpMTAN = 52 ± 10 and 85 ± 13 pM and Ki* values = 5.9 ± 0.3 and 10.0 ± 1.2 pM, respectively. HTDIA dissociates from heavy HpMTAN at 0.063 ± 0.002 min-1, faster than that from light HpMTAN at 0.032 ± 0.004 min-1. These values are consistent with transition-state formation by an improved catalytic site dynamic search and inconsistent with catalytic efficiency proportional to tight binding of the transition state.

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].

Dual Efficacy of a Catalytic Anti-Oligomeric Aβ42 scFv Antibody in Clearing Aβ42 Aggregates and Reducing Aβ Burden in the Brains of Alzheimer's Disease Mice

One of the primary pathological mechanisms underlying Alzheimer's disease (AD) is the deposition of amyloid β-protein (Aβ42) aggregates in the brain. In this study, a catalytic anti-oligomeric Aβ42 scFv antibody, HS72, was identified by screening a human antibody library, its ability to degrade Aβ42 aggregates was defined, and its role in the reduction of Aβ burden in the AD mouse brain was evaluated. HS72 specifically targeted Aβ42 aggregates with an approximately 14-68 kDa range. Based on molecular docking simulations, HS72 likely catalyzed the hydrolytic cleavage of the His13-His14 bond of Aβ42 chains in an Aβ42 aggregate unit, releasing N/C-terminal fragments and Aβ42 monomers. Degradation of Aβ42 aggregates by HS72 triggered a considerable disassembly or breakdown of the Aβ42 aggregates and greatly reduced their neurotoxicity. Aβ deposit/plaque load in the hippocampus of AD mice was reduced by approximately 27% after 7 days (once daily) of intravenous HS72 administration, while brain neural cells were greatly restored and their morphology was drastically improved. The above efficacies of HS72 were all greater than those of HT7, a simple anti-oligomeric Aβ42 scFv antibody. Although a catalytic anti-oligomeric Aβ42 antibody may have a slightly lower affinity for Aβ42 aggregates than a simple anti-oligomeric Aβ42 antibody, the former may display a stronger overall efficacy (dual efficacy of induction and catalysis) than the latter (induction alone) in clearing Aβ42 aggregates and improving histopathological changes in AD brain. Our findings on the catalytic antibody HS72 indicate the possibility of functional evolution of anti-oligomeric Aβ42 antibodies and provide novel insights into the immunotherapy of AD.

Enzymes of the crotonase superfamily: Diverse assembly and diverse function

The crotonase fold is generated by a framework of four repeats of a ββα-unit, extended by two helical regions. The active site of crotonase superfamily (CS) enzymes is located at the N-terminal end of the helix of the third repeat, typically being covered by a C-terminal helix. A major subset of CS-enzymes catalyzes acyl-CoA-dependent reactions, allowing for a diverse range of acyl-tail modifications. Most of these enzymes occur as trimers or hexamers (dimers of trimers), but monomeric forms are also observed. A common feature of the active sites of CS-enzymes is an oxyanion hole, formed by two peptide-NH hydrogen bond donors, which stabilises the negatively charged thioester oxygen atom of the reaction intermediate. Structural properties and possible use of these enzymes for biotechnological applications are discussed.

Human 2'-Deoxynucleoside 5'-Phosphate N-Hydrolase 1: Mechanism of 2'-Deoxyuridine 5'-Monophosphate Hydrolysis

The enzyme 2'-deoxynucleoside 5'-phosphate N-hydrolase 1 (DNPH1) catalyzes the N-ribosidic bond cleavage of 5-hydroxymethyl-2'-deoxyuridine 5'-monophosphate to generate 2-deoxyribose 5-phosphate and 5-hydroxymethyluracil. DNPH1 accepts other 2'-deoxynucleoside 5'-monophosphates as slow-reacting substrates. DNPH1 inhibition is a promising strategy to overcome resistance to and potentiate anticancer poly(ADP-ribose) polymerase inhibitors. We solved the crystal structure of unliganded human DNPH1 and took advantage of the slow reactivity of 2'-deoxyuridine 5'-monophosphate (dUMP) as a substrate to obtain a crystal structure of the DNPH1:dUMP Michaelis complex. In both structures, the carboxylate group of the catalytic Glu residue, proposed to act as a nucleophile in covalent catalysis, forms an apparent low-barrier hydrogen bond with the hydroxyl group of a conserved Tyr residue. The crystal structures are supported by functional data, with liquid chromatography-mass spectrometry analysis showing that DNPH1 incubation with dUMP leads to slow yet complete hydrolysis of the substrate. A direct UV-vis absorbance-based assay allowed characterization of DNPH1 kinetics at low dUMP concentrations. A bell-shaped pH-rate profile indicated that acid-base catalysis is operational and that for maximum kcat/KM, two groups with an average pKa of 6.4 must be deprotonated, while two groups with an average pKa of 8.2 must be protonated. A modestly inverse solvent viscosity effect rules out diffusional processes involved in dUMP binding to and possibly uracil release from the enzyme as rate limiting to kcat/KM. Solvent deuterium isotope effects on kcat/KM and kcat were inverse and unity, respectively. A reaction mechanism for dUMP hydrolysis is proposed.

History of advances in enzyme kinetic methods: From minutes to milliseconds

The last review on transient-state kinetic methods in The Enzymes was published three decades ago (Johnson, K.A., 1992. The Enzymes, XX, 1-61). In that review the foundations were laid out for the logic behind the design and interpretation of experiments. In the intervening years the instrumentation has improved mainly in providing better sample economy and shorter dead times. More significantly, in 1992 we were just introducing methods for fitting data based on numerical integration of rate equations, but the software was slow and difficult to use. Today, advances in numerical integration methods for data fitting have led to fast and dynamic software, making it easy to fit data without simplifying approximations. This approach overcomes multiple disadvantages of traditional data fitting based on equations derived by analytical integration of rate equations, requiring simplifying approximations. Mechanism-based fitting using computer simulation resolves mechanisms by accounting for the concentration dependence of the rates and amplitudes of the reaction to find a set of intrinsic rate constants that reproduce the experimental data, including details about how the experiment was performed in modeling the data. Rather than discuss how to design and interpret rapid-quench and stopped-flow experiments individually, we now focus on how to fit them simultaneously so that the quench-flow data anchor the interpretation of fluorescence signals. Computer simulation streamlines the fitting of multiple experiments globally to yield a single unifying model to account for all available data.

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.

Agonist efficiency links binding and gating in a nicotinic receptor

Receptors signal by switching between resting (C) and active (O) shapes ('gating') under the influence of agonists. The receptor's maximum response depends on the difference in agonist binding energy, O minus C. In nicotinic receptors, efficiency (η) represents the fraction of agonist binding energy applied to a local rearrangement (an induced fit) that initiates gating. In this receptor, free energy changes in gating and binding can be interchanged by the conversion factor η. Efficiencies estimated from concentration-response curves (23 agonists, 53 mutations) sort into five discrete classes (%): 0.56 (17), 0.51(32), 0.45(13), 0.41(26), and 0.31(12), implying that there are 5 C versus O binding site structural pairs. Within each class efficacy and affinity are corelated linearly, but multiple classes hide this relationship. η unites agonist binding with receptor gating and calibrates one link in a chain of coupled domain rearrangements that comprises the allosteric transition of the protein.

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.

Enzyme dynamics-a brief review

Debate has been simmering for some years regarding the importance of internal thermal motions of enzymes to catalysis. Recent developments in protein design may bring resolution of the more contentious points a little closer.

Design and interpretation of experiments to establish enzyme pathway and define the role of conformational changes in enzyme specificity

We describe the experimental methods and analysis to define the role of enzyme conformational changes in specificity based on published studies using DNA polymerases as an ideal model system. Rather than give details of how to perform transient-state and single-turnover kinetic experiments, we focus on the rationale of the experimental design and interpretation. We show how initial experiments to measure k_cat and k_cat/K_m can accurately quantify specificity but do not define its underlying mechanistic basis. We describe methods to fluorescently label enzymes to monitor conformational changes and to correlate fluorescence signals with rapid-chemical-quench flow assays to define the steps in the pathway. Measurements of the rate of product release and of the kinetics of the reverse reaction complete the kinetic and thermodynamic description of the full reaction pathway. This analysis showed that the substrate-induced change in enzyme structure from an open to a closed state was much faster than rate-limiting chemical bond formation. However, because the reverse of the conformational change was much slower than chemistry, specificity is governed solely by the product of the binding constant for the initial weak substrate binding and the rate constant for the conformational change (k_cat/K_m=K₁k₂) so that the specificity constant does not include k_cat. The enzyme conformational change leads to a closed complex in which the substrate is bound tightly and is committed to the forward reaction. In contrast, an incorrect substrate is bound weakly, and the rate of chemistry is slow, so the mismatch is released from the enzyme rapidly. Thus, the substrate-induced-fit is the major determinant of specificity. The methods outlined here should be applicable to other enzyme systems.

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.

Analyzing the Reaction of Orotidine 5'-Phosphate Decarboxylase as a Way to Examine Some Key Catalytic Proposals

This study analyzes the origin of enzyme catalysis by focusing on the reaction of orotidine 5'-phosphate decarboxylase (ODCase). This reaction involves an enormous catalytic effect of 23 kcal/mol that has been attributed to reactant state destabilization associated with the use of binding energy through the so-called Circe effect. However, our early studies and subsequent key experiments have shown that the presumed effect of the binding energy (namely, the strain exerted by a bond to a phosphate group) does not contribute to the catalysis. In this study, we perform quantitative empirical valence bond calculations that reproduce the catalytic effect of ODCase and the effect of removing the phosphate side chain. The calculations demonstrate that the effect of the phosphate is due to a change in reorganization energy and should not be described as an induced fit effect. Similarly, we show that the overall catalytic effect is due to electrostatic transition state stabilization, which again reflects the smaller reorganization energy in the enzyme than in water. We also elaborate on the problems with the induced fit proposal, including the fact that it does not serve to tell us what the actual origin of the action of the catalytic effect is. In addition to the above points, we use this paper to discuss misconceptions about the meaning of the preorganization effect, as well as other misunderstandings of what is being done in consistent calculations of enzyme catalysis.

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.