Protein modeling of human prorenin using the molecular dynamics method

To study the activation-inactivation mechanism of the renin zymogen, prorenin, a tertiary structural model of human prorenin was constructed using computer graphics and molecular dynamics calculations, based on the pepsinogen structure. This prorenin model shows that the folded prosegment polypeptide can fit into the substrate binding cleft of the renin moiety. The three positively charged residues, Arg 10, Arg 15, and Arg 20, in the prosegment make salt bridges with Asp 225, Glu 331, and Asp 60, respectively, in renin. Arg 43, which is in the processing site, forms salt bridges with the catalytic residues of Asp 81 and Asp 269. These ionic interactions between the prosegment and the renin may contribute to keeping the prorenin structure as an inactive form.     

Computational analysis of plasmepsin IV bound to an allophenylnorstatine inhibitor

The plasmepsin proteases from the malaria parasite Plasmodium falciparum are attracting attention as putative drug targets. A recently published crystal structure of Plasmodium malariae plasmepsin IV bound to an allophenylnorstatine inhibitor [Clemente, J.C. et al. (2006) Acta Crystallogr. D 62, 246-252] provides the first structural insights regarding interactions of this family of inhibitors with plasmepsins. The compounds in this class are potent inhibitors of HIV-1 protease, but also show nM binding affinities towards plasmepsin IV. Here, we utilize automated docking, molecular dynamics and binding free energy calculations with the linear interaction energy LIE method to investigate the binding of allophenylnorstatine inhibitors to plasmepsin IV from two different species. The calculations yield excellent agreement with experimental binding data and provide new information regarding protonation states of active site residues as well as conformational properties of the inhibitor complexes.     

Human prorenin structure sheds light on a novel mechanism of its autoinhibition and on its non-proteolytic activation by the (pro)renin receptor

Antibodies and prorenin mutants have long been used to structurally characterize prorenin, the inactive proenzyme form of renin. They were designed on the basis of homology models built using other aspartyl protease proenzyme structures since no structure was available for prorenin. Here, we present the first X-ray structure of a prorenin. The current structure of prorenin reveals that, in this zymogene, the active site of renin is blocked by the N-terminal residues of the mature version of the renin molecule, which are, in turn, covered by an Ω-shaped prosegment. This prevents access of substrates to the active site. The departure of the prosegment on activation induces an important global conformational change in the mature renin molecule with respect to prorenin: similar to other related enzymes such as pepsin or gastricsin, the segment that constitutes the N-terminal β-strand in renin is displaced from the renin active site by about 180° straight into the position that corresponds to the N-terminal β-strand of the prorenin prosegment. This way, the renin active site will become completely exposed and capable of carrying out its catalytic functions. A unique inactivation mechanism is also revealed, which does not make use of a lysine against the catalytic aspartates, probably in order to facilitate pH-independent activation [e.g., by the (pro)renin receptor].     

The protonation state of the catalytic aspartates in plasmepsin II

Assigning the correct protonation state to the catalytic residues is essential for a realistic modelling of an enzyme's active site. Plasmepsins are pharmaceutically relevant aspartic proteases involved in haemoglobin degradation by Plasmodium spp. In aspartic proteases, one of the two catalytic aspartates is protonated, while the other is negatively charged. Here, multiple explicit-water molecular dynamics simulations of plasmepsin II, uncomplexed and with a hydroxypropylamine peptidomimetic inhibitor, indicate that protonation of Asp214 favours a stable active site structure. Moreover, the protonation state of the catalytic aspartate has a strong influence on a linear chain of hydrogen bonds with the adjacent side chains.     

Mechanism of intramolecular activation of pepsinogen. Evidence for an intermediate delta and the involvement of the active site of pepsin in the intramolecular activation of pepsinogen.

Intramolecular pepsinogen activation is inhibited either by pepstatin, a potent pepsin inhibitor, or by purified globin from hemoglobin, a good pepsin substrate. Also, pepsinogen at pH 2 can be bound to a pepstatin-Sepharose column and recovered as native zymogen upon elution in pH 8 buffer. Kinetic studies of the globin inhibition of pepsinogen activation show that globin binds to a pepsinogen intermediate. This interaction gives rise to competitive inhibition of intramolecular pepsinogen activation. The evidence presented in this paper suggests that pepsinogen is converted rapidly upon acidification to the pepsinogen intermediate delta. In the absence of an inhibitor, the intermediate undergoes conformational change to bind the activation peptide portion of this same pepsinogen molecule in the active center to form an intramolecular enzyme-substrate complex (intermediate theta). This is followed by the intramolecular hydrolysis of the peptide bond between residues 44 and 45 of the pepsinogen molecule and the dissociation of the activation peptide from the pepsin. Intermediate delta apparently does not activate another pepsinogen molecule via an intermolecular process. Neither does intermediate delta hydrolyze globin substrate.     

Src kinase activation: A switched electrostatic network

Src tyrosine kinases are essential in numerous cell signaling pathways, and improper functioning of these enzymes has been implicated in many diseases. The activity of Src kinases is regulated by conformational activation, which involves several structural changes within the catalytic domain (CD): the orientation of two lobes of CD; rearrangement of the activation loop (A-loop); and movement of an alpha-helix (alphaC), which is located at the interface between the two lobes, into or away from the catalytic cleft. Conformational activation was investigated using biased molecular dynamics to explore the transition pathway between the active and the down-regulated conformation of CD for the Src-kinase family member Lyn kinase, and to gain insight into the interdependence of these changes. Lobe opening is observed to be a facile motion, whereas movement of the A-loop motion is more complex requiring secondary structure changes as well as communication with alphaC. A key result is that the conformational transition involves a switch in an electrostatic network of six polar residues between the active and the down-regulated conformations. The exchange between interactions links the three main motions of the CD. Kinetic experiments that would demonstrate the contribution of the switched electrostatic network to the enzyme mechanism are proposed. Possible implications for regulation conferred by interdomain interactions are also discussed.     

Flap flexibility amongst plasmepsins I,II, III,IV, and V: Sequence,structural, and molecular dynamics analyses

Herein, for the first time, we comparatively report the opening and closing of apo plasmepsin I – V. Plasmepsins belong the aspartic protease family of enzymes, and are expressed during the various stages of the P. falciparum lifecycle, the species responsible for the most lethal and virulent malaria to infect humans. Plasmepsin I, II, IV and HAP degrade hemoglobin from infected red blood cells, whereas plasmepsin V transport proteins crucial to the survival of the malaria parasite across the endoplasmic reticulum. Flap‐structures covering the active site of aspartic proteases (such as HIV protease) are crucial to the conformational flexibility and dynamics of the protein, and ultimately control the binding landscape. The flap‐structure in plasmepsins is made up of a flip tip in the N‐terminal lying perpendicular to the active site, adjacent to the flexible loop region in the C‐terminal. Using molecular dynamics, we propose three parameters to better describe the opening and closing of the flap‐structure in apo plasmepsins. Namely, the distance, d₁, between the flap tip and the flexible region; the dihedral angle, ?, to account for the twisting motion; and the TriCα angle, θ₁. Simulations have shown that as the flap‐structure twists, the flap and flexible region move apart opening the active site, or move toward each other closing the active site. The data from our study indicate that of all the plasmepsins investigated in the present study, Plm IV and V display the highest conformational flexibility and are more dynamic structures versus Plm I, II, and HAP. Proteins 2015; 83:1693–1705. © 2015 Wiley Periodicals, Inc.     

p38γ activation triggers dynamical changes in allosteric docking sites

Mitogen-activated protein kinases (MAPKs) are serine-threonine kinases that participate in signal transduction pathways. p38 MAPKs have four isoforms (p38α, p38β, p38γ, and p38δ) which are involved in multiple cellular functions such as proliferation, differentiation, survival, and migration. MAPK kinases phosphorylate p38s in the dual-phosphorylation motif, Thr-Gly-Tyr, located in their activation loop, which induces a conformational change that increases ATP binding affinity and catalytic activity. Several works have proposed that MAPK dynamics is a key factor in determining their function. However, we still do not understand the dynamical changes that lead to MAPK activation. In this work we have used molecular dynamics techniques to study the dynamical changes associated with p38γ activation, the only fully active MAPK crystallized so far. We performed MD simulations of p38γ in three different states, fully active with ATP, active without ATP, and inactive. We found that the dynamical fluctuations of the docking sites, important for protein-protein interactions, are regulated allosterically by changes in the active site. Interestingly, in the phosphorylated and ATP-bound states the whole protein dynamics lead to concerted motions of whole protein domains in contrast to the inactive state. The binding/unbinding of ATP participates in the reorientation of the two domains and in the regulation of protein plasticity. Our study shows that beyond the conformational changes associated with MAPK activation their correlated dynamics are highly regulated by phosphorylation and ATP binding. This means that MAPK plasticity may have a role in their catalytic activity, specificity, and protein-protein interactions and, therefore, in the outcome of the signaling network.     

Refined structure of porcine pepsinogen at 1.8 A resolution

The molecular structure of porcine pepsinogen at 1.8 A resolution has been determined by a combination of molecular replacement and multiple isomorphous phasing techniques. The resulting structure was refined by restrained-parameter least-squares methods. The final R factor [formula: see text] is 0.164 for 32,264 reflections with I greater than or equal to sigma (I) in the resolution range of 8.0 to 1.8 A. The model consists of 2785 protein atoms in 370 residues, a phosphoryl group on Ser68 and 238 ordered water molecules. The resulting molecular stereochemistry is consistent with a well-refined crystal structure with co-ordinate accuracy in the range of 0.10 to 0.15 A for the well-ordered regions of the molecule (B less than 15 A2). For the enzyme portion of the zymogen, the root-mean-square difference in C alpha atom co-ordinates with the refined porcine pepsin structure is 0.90 A (284 common atoms) and with the C alpha atoms of penicillopepsin it is 1.63 A (275 common atoms). The additional 44 N-terminal amino acids of the prosegment (Leu1p to Leu44p, using the letter p after the residue number to distinguish the residues of the prosegment) adopt a relatively compact structure consisting of a long beta-strand followed by two approximately orthogonal alpha-helices and a short 3(10)-helix. Intimate contacts, both electrostatic and hydrophobic interactions, are made with residues in the pepsin active site. The N-terminal beta-strand, Leu1p to Leu6p, forms part of the six-stranded beta-sheet common to the aspartic proteinases. In the zymogen the first 13 residues of pepsin, Ile1 to Glu13, adopt a completely different conformation from that of the mature enzyme. The C alpha atom of Ile1 must move approximately 44 A in going from its position in the inactive zymogen to its observed position in active pepsin. Electrostatic interactions of Lys36pN and hydrogen-bonding interactions of Tyr37pOH, and Tyr90H with the two catalytic aspartate groups, Asp32 and Asp215, prevent substrate access to the active site of the zymogen. We have made a detailed comparison of the mammalian pepsinogen fold with the fungal aspartic proteinase fold of penicillopepsin, used for the molecular replacement solution. A structurally derived alignment of the two sequences is presented.     

Predicting functional residues in Plasmodium falciparum plasmepsins by combining sequence and structural analysis with molecular dynamics simulations

Plasmepsins are aspartic proteases involved in the initial steps of the hemoglobin degradation pathway, a critical stage in the Plasmodium falciparum life cycle during human infection. Thus, they are attractive targets for novel therapeutic compounds to treat malaria, which remains one of the world's biggest health problems. The three-dimensional structures available for P. falciparum plasmepsins II and IV make structure-based drug design of antimalarial compounds that focus on inhibiting plasmepsins possible. However, the structural flexibility of the plasmepsin active site cavity combined with insufficient knowledge of the functional residues and of those determining the specificity of parasitic enzymes is a drawback when designing specific inhibitors. In this study, we have combined a sequence and structural analysis with molecular dynamics simulations to predict the functional residues in P. falciparum plasmepsins. The careful analysis of X-ray structures and 3D models carried out here suggests that residues Y17, V105, T108, L191, L242, Q275, and T298 are important for plasmepsin function. These seven amino acids are conserved across the malarial strains but not in human aspartic proteases. Residues V105 and T108 are localized in a flap of an interior pocket and they only establish contacts with a specific non-peptide achiral inhibitor. We also observed a rapid conformational change in the L3 region of plasmepsins that closes the active site of the enzyme, which explains earlier experimental findings. These results shed light on the role of V105 and T108 residues in plasmepsin specificities, and they should be useful in structure-based design of novel, selective inhibitors that may serve as antimalarial drugs.     

The role of the flap residue, threonine 77, in the activation and catalytic activity of pepsin A

Flexible loops, often referred to as flaps, have been shown to play a role in catalytic mechanisms of different enzymes. Flaps at the active site regions have been observed in the crystal structures of aspartic proteinases and their residues implicated in the catalytic processes. This research investigated the role of the flap residue, threonine 77, in the activation of pepsinogen and the catalytic mechanism of pepsin. Three mutants, T77S, T77V and T77G, were constructed. Differences in amino acid polarity and hydrogen bonding potential were shown to have an influence on the activation and catalytic processes. T77S activated at the same rate and had similar catalytic parameters as the wild-type pepsin. The activation rates of T77V and T77G were slower and their catalytic efficiencies lower than the wild-type. The results demonstrated that the threonine 77 polar side chain played a role in a proteolysis. The contribution of the side chain to zymogen activation was associated with the proteolytic cleavage of the prosegment. It was postulated that the hydroxyl group at position 77 provided an essential hydrogen bond that contributed to proper substrate alignment and, indirectly, to a catalytically favorable geometry of the transition state.     

Understanding the Mechanism of Prosegment-catalyzed Folding by Solution NMR Spectroscopy

Multidomain protein folding is often more complex than a two-state process, which leads to the spontaneous folding of the native state. Pepsin, a zymogen-derived enzyme, without its prosegment (PS), is irreversibly denatured and folds to a thermodynamically stable, non-native conformation, termed refolded pepsin, which is separated from native pepsin by a large activation barrier. While it is known that PS binds refolded pepsin and catalyzes its conversion to the native form, little structural details are known regarding this conversion. In this study, solution NMR was used to elucidate the PS-catalyzed folding mechanism by examining the key equilibrium states, e.g. native and refolded pepsin, both in the free and PS-bound states, and pepsinogen, the zymogen form of pepsin. Refolded pepsin was found to be partially structured and lacked the correct domain-domain structure and active-site cleft formed in the native state. Analysis of chemical shift data revealed that upon PS binding refolded pepsin folds into a state more similar to that of pepsinogen than to native pepsin. Comparison of pepsin folding by wild-type and mutant PSs, including a double mutant PS, indicated that hydrophobic interactions between residues of prosegment and refolded pepsin lower the folding activation barrier. A mechanism is proposed for the binding of PS to refolded pepsin and how the formation of the native structure is mediated.     

An electrostatic network and long-range regulation of Src kinases

The regulatory mechanism of Src tyrosine kinases includes conformational activation by a change in the catalytic domain tertiary structure and in domain-domain contacts between the catalytic domain and the SH2/SH3 regulatory domains. The kinase is activated when tyrosine phosphorylation occurs on the activation loop, but without phosphorylation of the C-terminal tail. Activation also occurs by allostery when contacts between the catalytic domain (CD) and the regulatory SH3 and SH2 domains are released as a result of exogenous protein binding. The aim of this work is to examine the proposed role of an electrostatic network in the conformational transition and to elucidate the molecular mechanism for long-range, allosteric conformational activation by using a combination of experimental enzyme kinetics and nonequilibrium molecular dynamics simulations. Salt dependence of the induction phase is observed in kinetic assays and supports the role of an electrostatic network in the transition. In addition, simulations provide evidence that allosteric activation involves a concerted motion coupling highly conserved residues, and spanning several nanometers from the catalytic site to the regulatory domain interface to communicate between the CD and the regulatory domains.     

A recombinant transductor-effector system: in vitro study of G inhibitory protein (G-alpha-i1) direct activators

Most protease prosegments are co-synthesized at the N-termini of cysteine proteases and are involved in folding assistance, inhibition, and activation of their mature enzymes. By using circular dichroism, UV-difference and fluorescence spectroscopies, we studied the thermal unfolding of papain prosegment. The transition seems to be two-state and reversible, with an unfolded state prone to aggregation. Unfolding thermodynamic parameters obtained show low values both for deltaH(Tm) and deltaCp(U), indicative of a loosely packed three-dimensional conformation for the prosegment at near-neutral pH conditions. In spite of these results, fluorescence experiments demonstrate that papain prosegment is able to recognize and inhibit its cognate protease. An acid medium induces a molten globule-like state without intermediates, which in turn undergoes an irreversible thermal unfolding. Our results suggest that papain prosegment has a high degree of conformational flexibility, with the ability to form not only a molten globule-like structure in activating conditions, but also requiring an induced fit in order to be functional as inhibitor.     

Correlation of the enzyme activities of Bacillus stearothermophilus lactate dehydrogenase on three substrates with the results of molecular dynamics/energy minimization conformational searching

Current methods for reengineering enzyme substrate specificities rely heavily on the use of static x-ray crystallographic models. In this article we detail the use of a molecular mechanics approach for suggesting regions of Bacillus stearothermophilus L-lactate dehydrogenase (EC 1.1.1.27) involved in substrate specificity, and hence areas of interest for protein engineers. The approach combines molecular dynamics with energy minimization (MD/EM) to search the conformational space available to a 15-Å sphere of the ternary complex centered on the catalytic histidine. The search is carried out by calculating a 30-ps dynamics trajectory at 300 K and minimizing structures at 1-ps intervals. The protocol has been performed on 14 systems containing different combinations of substrate and mutant /wt LDH. In order to discover which interactions are important in defining substrate specificity, eight conformational parameters representing substrate–active site interactions were measured in each of the 420 minimized structures. These parameters were then compared to the measured catalytic activity of the protein–substrate combinations. These comparisons show that arginine 109 orientation is a major determining factor in LDH specificity. Using this methodolgy it is possible to estimate the catalytic activity of proteins of varied sequence by computer simulation before synthesis. Proteins 29:228–239, 1997. © 1997 Wiley-Liss, Inc.     

Effect of molecular confinement on internal enzyme dynamics: frequency domain fluorometry and molecular dynamics simulation studies

The tryptophanyl emission decay of the mesophilic beta-galactosidase from Aspergillus oryzae free in buffer and entrapped in agarose gel is investigated as a function of temperature and compared to that of the hyperthermophilic enzyme from Sulfolobus solfataricus. Both enzymes are tetrameric proteins with a large number of tryptophanyl residues, so the fluorescence emission can provide information on the conformational dynamics of the overall protein structure rather than that of the local environment. The tryptophanyl emission decays are best fitted by bimodal Lorentzian distributions. The long-lived component is ascribed to close, deeply buried tryptophanyl residues with reduced mobility; the short-lived one arises from tryptophanyl residues located in more flexible external regions of each subunit, some of which are involved in forming the catalytic site. The center of both lifetime distribution components at each temperature increases when going from the free in solution mesophilic enzyme to the gel-entrapped and hyperthermophilic enzyme, thus indicating that confinement of the mesophilic enzyme in the agarose gel limits the freedom of the polypeptide chain. A more complex dependence is observed for the distribution widths. Computer modeling techniques are used to recognize that the catalytic sites are similar for the mesophilic and hyperthermophilic beta-galactosidases. The effect due to gel entrapment is considered in dynamic simulations by imposing harmonic restraints to solvent-exposed atoms of the protein with the exclusion of those around the active site. The temperature dependence of the tryptophanyl fluorescence emission decay and the dynamic simulation confirm that more rigid structures, as in the case of the immobilized and/or hyperthermophilic enzyme, require higher temperatures to achieve the requisite conformational dynamics for an effective catalytic action and strongly suggest a link between conformational rigidity and enhanced thermal stability.     

Catalytic loop motion in human glutathione synthetase: A molecular modeling approach

Conformational changes of three flexible loops (G, A, and S) in human glutathione synthetase (hGS) arise to accommodate the substrates inside the active site. The crystal structure of hGS, a member of the ATP-grasp superfamily, has been reported only for the product-enzyme complex. To study the function of the hGS loops, molecular dynamics simulations are performed on three different conformational models: unbound enzyme, reactant-enzyme, and product-enzyme complex of hGS. The conformational changes among the three models are analyzed and the roles of the loops during the catalytic process are described. The modeled structures of hGS show that the central portions of the G- and A-loop have a double role in the reactant complex conformation: they bind the substrates and simultaneously interact with each other through an extensive network of hydrogen bonds. The present study proposes that these favorable loop-ligand and loop-loop interactions are required for opening and closing of the active site of hGS. Additionally, this research identifies important amino acid residues and explains their function within the catalytic loops of hGS.     

Active-site residues move independently from the rest of the protein in a 200 ns molecular dynamics simulation of cytochrome P450 CYP119

The conformational dynamics of cytochrome P450 enzymes are critical to their catalytic activity. In this study, the correlated motion between residues in a 200 ns molecular dynamics trajectory of the thermophilic CYP119 was analyzed to parse out conformational relationships. Residues that are structurally related, for example residues within a helix, generally have highly correlated motion. In addition, clusters of non-adjacent residues that show correlated motion (“hot spots”) are seen in various regions, including at the base of the F and G helices that make up the most dynamic region of the enzyme. A modified k-means algorithm that clusters residues based on their correlated motion indicates that functionally related residues are in the same cluster (e.g., the catalytic threonines and the heme). Tightly coupled clusters form a solvent-exposed “shell” around the enzyme, whereas less coupling between clusters is seen in regions that are critical to ligand interactions, redox partner interactions, and catalysis. Most notably, we find that residues in the active site move independently from the rest of the enzyme, effectively insulating the catalytic machinery from other regions of the protein.     

Hierarchical Modeling of Activation Mechanisms in the ABL and EGFR Kinase Domains: Thermodynamic and Mechanistic Catalysts of Kinase Activation by Cancer Mutations

Structural and functional studies of the ABL and EGFR kinase domains have recently suggested a common mechanism of activation by cancer-causing mutations. However, dynamics and mechanistic aspects of kinase activation by cancer mutations that stimulate conformational transitions and thermodynamic stabilization of the constitutively active kinase form remain elusive. We present a large-scale computational investigation of activation mechanisms in the ABL and EGFR kinase domains by a panel of clinically important cancer mutants ABL-T315I, ABL-L387M, EGFR-T790M, and EGFR-L858R. We have also simulated the activating effect of the gatekeeper mutation on conformational dynamics and allosteric interactions in functional states of the ABL-SH2-SH3 regulatory complexes. A comprehensive analysis was conducted using a hierarchy of computational approaches that included homology modeling, molecular dynamics simulations, protein stability analysis, targeted molecular dynamics, and molecular docking. Collectively, the results of this study have revealed thermodynamic and mechanistic catalysts of kinase activation by major cancer-causing mutations in the ABL and EGFR kinase domains. By using multiple crystallographic states of ABL and EGFR, computer simulations have allowed one to map dynamics of conformational fluctuations and transitions in the normal (wild-type) and oncogenic kinase forms. A proposed multi-stage mechanistic model of activation involves a series of cooperative transitions between different conformational states, including assembly of the hydrophobic spine, the formation of the Src-like intermediate structure, and a cooperative breakage and formation of characteristic salt bridges, which signify transition to the active kinase form. We suggest that molecular mechanisms of activation by cancer mutations could mimic the activation process of the normal kinase, yet exploiting conserved structural catalysts to accelerate a conformational transition and the enhanced stabilization of the active kinase form. The results of this study reconcile current experimental data with insights from theoretical approaches, pointing to general mechanistic aspects of activating transitions in protein kinases.     

Computational studies of tryptophanyl-tRNA synthetase: activation of ATP by induced-fit

Catalysis of amino acid activation by Bacillus stearothermophilus tryptophanyl-tRNA synthetase involves three allosteric states: (1) Open; (2) closed pre-transition state (PreTS); and (3) closed products (Product). The interconversions of these states entail significant domain motions driven by ligand binding. We explore the application of molecular dynamics simulations to investigate ligand-linked conformational stability changes associated with this catalytic cycle. Multiple molecular dynamics trajectories (5 ns) for 11 distinct liganded and unliganded monomer configurations show that the PreTS conformation is unstable in the absence of ATP, reverting within approximately 600 ps nearly to the Open conformation. In contrast, Open and Product state trajectories were stable, even without ligands, confirming the previous suggestion that catalysis entails destabilization of the protein conformation, driven by ATP-binding energies developed as the PreTS state assembles during induced-fit. The simulations suggest novel mechanistic details associated with both induced-fit (Open-PreTS) and catalysis (PreTS-Product). Notably, Mg2+ -ATP interactions are coupled to interactions between ATP and active-site lysine side-chains via mechanisms that cannot be captured under the molecular mechanics approximations, and which therefore require restraining potentials for stable simulation. Simulations of Mg2+. ATP-bound PreTS complexes with restraining potentials and with a virtual K111A mutant confirm that these coupling interactions are necessary to sustain the PreTS conformation and, in turn, provide a new model for how the PreTS conformation activates ATP for catalysis. These results emphasize the central role of the PreTS state as a high-energy intermediate structure along the catalytic pathway and suggest that Mg2+ and the KMSKS loop function cooperatively during catalysis.