Ligand-driven Conformational Changes Research Articles

Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release.

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

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

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution.

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

Hydride Transfer Catalyzed by Glycerol Phosphate Dehydrogenase: Recruitment of an Acidic Amino Acid Side Chain to Rescue a Damaged Enzyme.

K120 of glycerol 3-phosphate dehydrogenase (GPDH) lies close to the carbonyl group of the bound dihydroxyacetone phosphate (DHAP) dianion. pH rate (pH 4.6–9.0) profiles are reported for kcat and (kcat/Km)dianion for wild type and K120A GPDH-catalyzed reduction of DHAP by NADH, and for (kcat/KdKam) for activation of the variant-catalyzed reduction by CH3CH2NH3+, where Kam and Kd are apparent dissociation constants for CH3CH2NH3+ and DHAP, respectively. These profiles provide evidence that the K120 side chain cation, which is stabilized by an ion-pairing interaction with the D260 side chain, remains protonated between pH 4.6 and 9.0. The profiles for wild type and K120A variant GPDH show downward breaks at a similar pH value (7.6) that are attributed to protonation of the K204 side chain, which also lies close to the substrate carbonyl oxygen. The pH profiles for (kcat/Km)dianion and (kcat/KdKam) for the K120A variant show that the monoprotonated form of the variant is activated for catalysis by CH3CH2NH3+ but has no detectable activity, compared to the diprotonated variant, for unactivated reduction of DHAP. The pH profile for kcat shows that the monoprotonated K120A variant is active toward reduction of enzyme-bound DHAP, because of activation by a ligand-driven conformational change. Upward breaks in the pH profiles for kcat and (kcat/Km)dianion for K120A GPDH are attributed to protonation of D260. These breaks are consistent with the functional replacement of K120 by D260, and a plasticity in the catalytic roles of the active site side chains.

Orotidine 5'-Monophosphate Decarboxylase: The Operation of Active Site Chains Within and Across Protein Subunits.

The D37 and T100′ side chains of orotidine 5′-monophosphate decarboxylase (OMPDC) interact with the C-3′ and C-2′ ribosyl hydroxyl groups, respectively, of the bound substrate. We compare the intra-subunit interactions of D37 with the inter-subunit interactions of T100′ by determining the effects of the D37G, D37A, T100′G, and T100′A substitutions on the following: (a) kcat and kcat/Km values for the OMPDC-catalyzed decarboxylations of OMP and 5-fluoroorotidine 5′-monophosphate (FOMP) and (b) the stability of dimeric OMPDC relative to the monomer. The D37G and T100′A substitutions resulted in 2 kcal mol–1 increases in ΔG† for kcat/Km for the decarboxylation of OMP, while the D37A and T100′G substitutions resulted in larger 4 and 5 kcal mol–1 increases, respectively, in ΔG†. The D37G and T100′A substitutions both resulted in smaller 2 kcal mol–1 decreases in ΔG† for the decarboxylation of FOMP compared to that of OMP. These results show that the D37G and T100′A substitutions affect the barrier to the chemical decarboxylation step while the D37A and T100′G substitutions also affect the barrier to a slow, ligand-driven enzyme conformational change. Substrate binding induces the movement of an α-helix (G′98–S′106) toward the substrate C-2′ ribosyl hydroxy bound at the main subunit. The T100′G substitution destabilizes the enzyme dimer by 3.5 kcal mol–1 compared to the monomer, which is consistent with the known destabilization of α-helices by the internal Gly side chains [Serrano, L., et al. (1992) Nature, 356, 453–455]. We propose that the T100′G substitution weakens the α-helical contacts at the dimer interface, which results in a decrease in the dimer stability and an increase in the barrier to the ligand-driven conformational change.

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.

Role of Ligand-Driven Conformational Changes in Enzyme Catalysis: Modeling the Reactivity of the Catalytic Cage of Triosephosphate Isomerase.

We have previously performed empirical valence bond calculations of the kinetic activation barriers, ΔG‡calc, for the deprotonation of complexes between TIM and the whole substrate glyceraldehyde-3-phosphate (GAP, Kulkarni et al.J. Am. Chem. Soc.2017, 139, 10514–1052528683550). We now extend this work to also study the deprotonation of the substrate pieces glycolaldehyde (GA) and GA·HPi [HPi = phosphite dianion]. Our combined calculations provide activation barriers, ΔG‡calc, for the TIM-catalyzed deprotonation of GAP (12.9 ± 0.8 kcal·mol–1), of the substrate piece GA (15.0 ± 2.4 kcal·mol–1), and of the pieces GA·HPi (15.5 ± 3.5 kcal·mol–1). The effect of bound dianion on ΔG‡calc is small (≤2.6 kcal·mol–1), in comparison to the much larger 12.0 and 5.8 kcal·mol–1 intrinsic phosphodianion and phosphite dianion binding energy utilized to stabilize the transition states for TIM-catalyzed deprotonation of GAP and GA·HPi, respectively. This shows that the dianion binding energy is essentially fully expressed at our protein model for the Michaelis complex, where it is utilized to drive an activating change in enzyme conformation. The results represent an example of the synergistic use of results from experiments and calculations to advance our understanding of enzymatic reaction mechanisms.

Osmium tetroxide as a probe of RNA structure.

Structured RNAs have a central role in cellular function. The capability of structured RNAs to adopt fixed architectural structures or undergo dynamic conformational changes contributes to their diverse role in the regulation of gene expression. Although numerous biophysical and biochemical tools have been developed to study structured RNAs, there is a continuing need for the development of new methods for the investigation of RNA structures, especially methods that allow RNA structure to be studied in solution close to its native cellular conditions. Here we use osmium tetroxide (OsO4) as a chemical probe of RNA structure. In this method, we have used fluorescence-based sequencing technologies to detect OsO4 modified RNA. We characterized the requirements for OsO4 modification of RNA by investigating three known structured RNAs: the M-box, glycine riboswitch RNAs, and tRNAasp. Our results show that OsO4 predominantly modifies RNA at uracils that are conformationally exposed on the surface of the RNA. We also show that changes in OsO4 reactivity at flexible positions in the RNA correlate with ligand-driven conformational changes in the RNA structure. Osmium tetroxide modification of RNA will provide insights into the structural features of RNAs that are relevant to their underlying biological functions.

Ligand-driven conformational changes of MurD visualized by paramagnetic NMR.

Proteins, especially multi-domain proteins, often undergo drastic conformational changes upon binding to ligands or by post-translational modifications, which is a key step to regulate their function. However, the detailed mechanisms of such dynamic regulation of the functional processes are poorly understood because of the lack of an efficient tool. We here demonstrate detailed characterization of conformational changes of MurD, a 47 kDa protein enzyme consisting of three domains, by the use of solution NMR equipped with paramagnetic lanthanide probe. Quantitative analysis of pseudocontact shifts has identified a novel conformational state of MurD, named semi-closed conformation, which is found to be the key to understand how MurD regulates the binding of the ligands. The modulation of the affinity coupled with conformational changes accentuates the importance of conformational state to be evaluated in drug design.

Mechanism for Activation of Triosephosphate Isomerase by Phosphite Dianion: The Role of a Ligand-Driven Conformational Change

The L232A mutation in triosephosphate isomerase (TIM) from Trypanosoma brucei brucei results in a small 6-fold decrease in k(cat)/K(m) for the reversible enzyme-catalyzed isomerization of glyceraldehyde 3-phosphate to give dihydroxyacetone phosphate. In contrast, this mutation leads to a 17-fold increase in the second-order rate constant for the TIM-catalyzed proton transfer reaction of the truncated substrate piece [1-(13)C]glycolaldehyde ([1-(13)C]-GA) in D(2)O, a 25-fold increase in the third-order rate constant for the reaction of the substrate pieces GA and phosphite dianion (HPO(3)(2-)), and a 16-fold decrease in K(d) for binding of HPO(3)(2-) to the free enzyme. Most significantly, the mutation also results in an 11-fold decrease in the extent of activation of the enzyme toward turnover of GA by bound HPO(3)(2-). The data provide striking evidence that the L232A mutation leads to a ca. 1.7 kcal/mol stabilization of a catalytically active loop-closed form of TIM (E(c)) relative to an inactive open form (E(o)). We propose that this is due to the relief, in L232A mutant TIM, of unfavorable steric interactions between the bulky hydrophobic side chain of Leu-232 and the basic carboxylate side chain of Glu-167, the catalytic base, which destabilize E(c) relative to E(o).

Helicase Motif III in SecA is essential for coupling preprotein binding to translocation ATPase.

The SecA ATPase is a protein translocase motor and a superfamily 2 (SF2) RNA helicase. The ATPase catalytic core ('DEAD motor') contains the seven conserved SF2 motifs. Here, we demonstrate that Motif III is essential for SecA-mediated protein translocation and viability. SecA Motif III mutants can bind ligands (nucleotide, the SecYEG translocase 'channel', signal and mature preprotein domains), can catalyse basal and SecYEG-stimulated ATP hydrolysis and can be activated for catalysis. However, Motif III mutation specifically blocks the preprotein-stimulated 'translocation ATPase' at a step of the reaction pathway that lies downstream of ligand binding. A functional Motif III is required for optimal ligand-driven conformational changes and kinetic parameters that underlie optimal preprotein-modulated nucleotide cycling at the SecA DEAD motor. We propose that helicase Motif III couples preprotein binding to the SecA translocation ATPase and that catalytic activation of SF2 enzymes through Motif-III-mediated action is essential for both polypeptide and nucleic-acid substrates.