Characterizing Active Site Conformational Heterogeneity along the Trajectory of an Enzymatic Phosphoryl Transfer Reaction

Magnesium induced structural reorganization in the active site of adenylate kinase

Phosphoryl transfer is a fundamental reaction in cellular signaling and metabolism that requires Mg2+ as an essential cofactor. While the primary function of Mg2+ is electrostatic activation of substrates, such as ATP, the full spectrum of catalytic mechanisms exerted by Mg2+ is not known. In this study, we integrate structural biology methods, molecular dynamic (MD) simulations, phylogeny, and enzymology assays to provide molecular insights into Mg2+-dependent structural reorganization in the active site of the metabolic enzyme adenylate kinase. Our results demonstrate that Mg2+ induces a conformational rearrangement of the substrates (ATP and ADP), resulting in a 30° adjustment of the angle essential for reversible phosphoryl transfer, thereby optimizing it for catalysis. MD simulations revealed transitions between conformational substates that link the fluctuation of the angle to large-scale enzyme dynamics. The findings contribute detailed insight into Mg2+ activation of enzymes and may be relevant for reversible and irreversible phosphoryl transfer reactions.

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.

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.

Mechanistic Basis for a Connection between the Catalytic Step and Slow Opening Dynamics of Adenylate Kinase

Escherichia coli adenylate kinase (AdK) is a small, monomeric enzyme that synchronizes the catalytic step with the enzyme's conformational dynamics to optimize a phosphoryl transfer reaction and the subsequent release of the product. Guided by experimental measurements of low catalytic activity in seven single-point mutation AdK variants (K13Q, R36A, R88A, R123A, R156K, R167A, and D158A), we utilized classical mechanical simulations to probe mutant dynamics linked to product release, and quantum mechanical and molecular mechanical calculations to compute a free energy barrier for the catalytic event. The goal was to establish a mechanistic connection between the two activities. Our calculations of the free energy barriers in AdK variants were in line with those from experiments, and conformational dynamics consistently demonstrated an enhanced tendency toward enzyme opening. This indicates that the catalytic residues in the wild-type AdK serve a dual role in this enzyme's function─one to lower the energy barrier for the phosphoryl transfer reaction and another to delay enzyme opening, maintaining it in a catalytically active, closed conformation for long enough to enable the subsequent chemical step. Our study also discovers that while each catalytic residue individually contributes to facilitating the catalysis, R36, R123, R156, R167, and D158 are organized in a tightly coordinated interaction network and collectively modulate AdK's conformational transitions. Unlike the existing notion of product release being rate-limiting, our results suggest a mechanistic interconnection between the chemical step and the enzyme's conformational dynamics acting as the bottleneck of the catalytic process. Our results also suggest that the enzyme's active site has evolved to optimize the chemical reaction step while slowing down the overall opening dynamics of the enzyme.

Protein dynamics: The future is bright and complicated!

Biological life depends on motion, and this manifests itself in proteins that display motion over a formidable range of time scales spanning from femtoseconds vibrations of atoms at enzymatic transition states, all the way to slow domain motions occurring on micro to milliseconds. An outstanding challenge in contemporary biophysics and structural biology is a quantitative understanding of the linkages among protein structure, dynamics, and function. These linkages are becoming increasingly explorable due to conceptual and methodological advances. In this Perspective article, we will point toward future directions of the field of protein dynamics with an emphasis on enzymes. Research questions in the field are becoming increasingly complex such as the mechanistic understanding of high-order interaction networks in allosteric signal propagation through a protein matrix, or the connection between local and collective motions. In analogy to the solution to the "protein folding problem," we argue that the way forward to understanding these and other important questions lies in the successful integration of experiment and computation, while utilizing the present rapid expansion of sequence and structure space. Looking forward, the future is bright, and we are in a period where we are on the doorstep to, at least in part, comprehend the importance of dynamics for biological function.

Minimizing Pervasive Artifacts in 4D Covariance Maps for Protein Side Chain NMR Assignments

Nuclear magnetic resonance (NMR) is a mainstay of biophysical studies that provides atomic level readouts to formulate molecular mechanisms. Side chains are particularly important to derive mechanisms involving proteins as they carry functional groups, but NMR studies of side chains are often limited by challenges in assigning their signals. Here, we designed a novel computational method that combines spectral derivatives and matrix square-rooting to produce reliable 4D covariance maps from routinely acquired 3D spectra and facilitates side chain resonance assignments. Thus, we generate two 4D maps from 3D-HcccoNH and 3D-HCcH-TOCSY spectra that each help overcome signal overlap or sensitivity losses. These 4D maps feature HC-HSQCs of individual side chains that can be paired to assigned backbone amide resonances of individual aliphatic signals, and both are obtained from a single modified covariance calculation. Further, we present 4D maps produced using conventional triple resonance experiments to easily assign asparagine side chain amide resonances. The 4D covariance maps encapsulate the lengthy manual pattern recognition used in traditional assignment methods and distill the information as correlations that can be easily visualized. We showcase the utility of the 4D covariance maps with a 10 kDa peptidyl carrier protein and a 52 kDa cyclization domain from a nonribosomal peptide synthetase.

Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions

Enzymes employ a wide range of protein motions to achieve efficient catalysis of chemical reactions. While the role of collective protein motions in substrate binding, product release, and regulation of enzymatic activity is generally understood, their roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations, enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the roles of catalytic residues in catalysis and the conformational change in the enzyme. This study reveals that the motions in the active site, which occur on a time scale of picoseconds to nanoseconds, link the catalytic reaction to the slow conformational dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In particular, substantial conformational rearrangement occurs in the active site following the catalytic reaction. This rearrangement not only affects the reaction barrier but also promotes a more open conformation of the enzyme after the reaction, which then results in an accelerated opening of the enzyme compared to that of the reactant state. The results illustrate a linkage between enzymatic catalysis and collective protein motions, whereby the disparate time scales between the two processes are bridged by a cascade of intermediate-scale motion of catalytic residues modulating the free energy landscapes of the catalytic and conformational change processes.

Serine protease dynamics revealed by NMR analysis of the thrombin-thrombomodulin complex

Serine proteases catalyze a multi-step covalent catalytic mechanism of peptide bond cleavage. It has long been assumed that serine proteases including thrombin carry-out catalysis without significant conformational rearrangement of their stable two-β-barrel structure. We present nuclear magnetic resonance (NMR) and hydrogen deuterium exchange mass spectrometry (HDX-MS) experiments on the thrombin-thrombomodulin (TM) complex. Thrombin promotes procoagulative fibrinogen cleavage when fibrinogen engages both the anion binding exosite 1 (ABE1) and the active site. It is thought that TM promotes cleavage of protein C by engaging ABE1 in a similar manner as fibrinogen. Thus, the thrombin-TM complex may represent the catalytically active, ABE1-engaged thrombin. Compared to apo- and active site inhibited-thrombin, we show that thrombin-TM has reduced μs-ms dynamics in the substrate binding (S1) pocket consistent with its known acceleration of protein C binding. Thrombin-TM has increased μs-ms dynamics in a β-strand connecting the TM binding site to the catalytic aspartate. Finally, thrombin-TM had doublet peaks indicative of dynamics that are slow on the NMR timescale in residues along the interface between the two β-barrels. Such dynamics may be responsible for facilitating the N-terminal product release and water molecule entry that are required for hydrolysis of the acyl-enzyme intermediate.

Reactions of aryl dimethylphosphinothioate esters with anionic oxygen nucleophiles: transition state structure in 70% water-30% ethanol

Aryl dimethylphosphinates, 2, react with anionic oxygen nucleophiles in water via a concerted (ANDN) mechanism. With EtO- in anhydrous ethanol, the mechanism is associative (AN + DN), with rate-limiting pentacoordinate intermediate formation. This change in mechanism with solvent change has been ascribed to changes in the nucleophile and leaving group basicities accompanying solvent change. This paper reports on a kinetic analysis of the reactions of the aryl dimethylphosphinothioates, 3a-g, with oxygen nucleophiles in 70% water-30% ethanol (v/v) solvent at 25 °C, reactions known to proceed by a concerted mechanism in water, to test the rationalization stated above, since the nucleophiles and LGs of interest are more basic in aqueous ethanol than in water. The change in solvent causes an ca. 14 to 320-fold decrease in rate. Hammett and Brønsted-type correlations characterize a concerted TS with less P-LG bonding in aqueous ethanol than in water. Two opposing consequences are associated with the solvent change: (a) increased basicity of nucleophiles and LGs, which lead to a modest tightening of the TS; and (b) better stabilization of the IS relative to the TS in aqueous ethanol, which results in a slower reaction with a more product-like TS. Hammond and anti-Hammond effects on the TS arising from better stabilization of the IS over the TS dominate over the effects of increased nucleophile and LG basicity in determining the looser TS structure in aqueous ethanol. An altered TS structure is consistent with an altered reaction potential energy surface, in this case caused by a change in solvent polarity.

Determining isoleucine side-chain rotamer-sampling in proteins from 13C chemical shift

A framework is presented to derive the conformational sampling of isoleucine side chains from nuclear magnetic resonance ¹³C chemical shifts.

Chemical shifts are often the only nuclear magnetic resonance parameter that can be obtained for challenging macromolecular systems. Here we present a framework to derive the conformational sampling of isoleucine side chains from ¹³C chemical shifts and demonstrate that side-chain conformations in a low-populated folding intermediate can be determined.

Characterising side chains in large proteins by protonless 13C-detected NMR spectroscopy

Side chains cover protein surfaces and are fundamental to processes as diverse as substrate recognition, protein folding and enzyme catalysis. However, characterisation of side-chain motions has so far been restricted to small proteins and methyl-bearing side chains. Here we present a class of methods, based on ¹³C-detected NMR spectroscopy, to more generally quantify motions and interactions of side chains in medium-to-large proteins. A single, uniformly isotopically labelled sample is sufficient to characterise the side chains of six different amino acid types. Side-chain conformational dynamics on the millisecond time-scale can be quantified by incorporating chemical exchange saturation transfer (CEST) into the presented methods, whilst long-range ¹³C-¹³C scalar couplings reporting on nanosecond to millisecond motions can be quantified in proteins as large as 80 kDa. The presented class of methods promises characterisation of side-chain behaviour at a level that has so far been reserved for the protein backbone.

Analysis of side-chain motions by NMR has so far been restricted to small proteins and methyl-bearing side chains. Here, the authors present NMR methods based on ¹³C direct detection of highly deuterated protein samples that yield sharp and well-resolved signals and allow the characterisation of side-chain conformational dynamics of six different amino acid types in medium-to-large proteins.

Arginine Side-Chain Hydrogen Exchange: Quantifying Arginine Side-Chain Interactions in Solution

The rate with which labile backbone hydrogen atoms in proteins exchange with the solvent has long been used to probe protein interactions in aqueous solutions. Arginine, an essential amino acid found in many interaction interfaces, is capable of an impressive range of interactions via its guanidinium group. The hydrogen exchange rate of the guanidinium hydrogens therefore becomes an important measure to quantify side-chain interactions. Herein we present an NMR method to quantify the hydrogen exchange rates of arginine side-chain 1 Hϵ protons and thus present a method to gauge the strength of arginine side-chain interactions. The method employs 13 C-detection and the one-bond deuterium isotope shift observed for 15 Nϵ to generate two exchanging species in 1 H2 O/2 H2 O mixtures. An application to the protein T4 Lysozyme is shown, where protection factors calculated from the obtained exchange rates correlate well with the interactions observed in the crystal structure. The methodology presented provides an important step towards characterising interactions of arginine side-chains in enzymes, in phase separation, and in protein interaction interfaces in general.

A 13C-detected 15N double-quantum NMR experiment to probe arginine side-chain guanidinium 15Nη chemical shifts

Arginine side-chains are often key for enzyme catalysis, protein-ligand and protein-protein interactions. The importance of arginine stems from the ability of the terminal guanidinium group to form many key interactions, such as hydrogen bonds and salt bridges, as well as its perpetual positive charge. We present here an arginine ¹³C^ζ-detected NMR experiment in which a double-quantum coherence involving the two ¹⁵N^η nuclei is evolved during the indirect chemical shift evolution period. As the precession frequency of the double-quantum coherence is insensitive to exchange of the two ¹⁵N^η; this new approach is shown to eliminate the previously deleterious line broadenings of ¹⁵N^η resonances caused by the partially restricted rotation about the C^ζ-N^ε bond. Consequently, sharp and well-resolved ¹⁵N^η resonances can be observed. The utility of the presented method is demonstrated on the L99A mutant of the 19 kDa protein T4 lysozyme, where the measurement of small chemical shift perturbations, such as one-bond deuterium isotope shifts, of the arginine amine ¹⁵N^η nuclei becomes possible using the double-quantum experiment.

Determining rotational dynamics of the guanidino group of arginine side chains in proteins by carbon-detected NMR

A new NMR-based method is presented to determine the rotational dynamics around the N_ε–C_ζ bond of arginine to characterise the interactions mediated by arginine side chains.

Arginine residues are imperative for many active sites and protein-interaction interfaces. A new NMR-based method is presented to determine the rotational dynamics around the N_ε–C_ζ bond of arginine side chains. An application to a 19 kDa protein shows that the strengths of interactions involving arginine side chains can be characterised.

Structural studies of a hyperthermophilic thymidylate kinase enzyme reveal conformational substates along the reaction coordinate

Thymidylate kinase (TMK) is a key enzyme which plays an important role in DNA synthesis. It belongs to the family of nucleoside monophosphate kinases, several of which undergo structure-encoded conformational changes to perform their function. However, the absence of three-dimensional structures for all the different reaction intermediates of a single TMK homolog hinders a clear understanding of its functional mechanism. We herein report the different conformational states along the reaction coordinate of a hyperthermophilic TMK from Aquifex aeolicus, determined via X-ray diffraction and further validated through normal-mode studies. The analyses implicate an arginine residue in the Lid region in catalysis, which was confirmed through site-directed mutagenesis and subsequent enzyme assays on the wild-type protein and mutants. Furthermore, the enzyme was found to exhibit broad specificity toward phosphate group acceptor nucleotides. Our comprehensive analyses of the conformational landscape of TMK, together with associated biochemical experiments, provide insights into the mechanistic details of TMK-driven catalysis, for example, the order of substrate binding and the reaction mechanism for phosphate transfer. Such a study has utility in the design of potent inhibitors for these enzymes.

DATABASE

Structural data are available in the PDB under the accession numbers 2PBR, 4S2E, 5H5B, 5XAI, 4S35, 5XB2, 5H56, 5XB3, 5H5K, 5XB5, and 5XBH.

Characterizing Active Site Conformational Heterogeneity along the Trajectory of an Enzymatic Phosphoryl Transfer Reaction

States along the phosphoryl transfer reaction catalyzed by the nucleoside monophosphate kinase UmpK were captured and changes in the conformational heterogeneity of conserved active site arginine side‐chains were quantified by NMR spin‐relaxation methods. In addition to apo and ligand‐bound UmpK, a transition state analog (TSA) complex was utilized to evaluate the extent to which active site conformational entropy contributes to the transition state free energy. The catalytically essential arginine side‐chain guanidino groups were found to be remarkably rigid in the TSA complex, indicating that the enzyme has evolved to restrict the conformational freedom along its reaction path over the energy landscape, which in turn allows the phosphoryl transfer to occur selectively by avoiding side reactions.