Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.     

Analysis of the isolated SecA DEAD motor suggests a mechanism for chemical-mechanical coupling

The preprotein cross-linking domain and C-terminal domains of Escherichia coli SecA were removed to create a minimal DEAD motor, SecA-DM. SecA-DM hydrolyzes ATP and has the same affinity for ADP as full-length SecA. The crystal structure of SecA-DM in complex with ADP was solved and shows the DEAD motor in a closed conformation. Comparison with the structure of the E. coli DEAD motor in an open conformation (Protein Data Bank ID 2FSI) indicates main-chain conformational changes in two critical sequences corresponding to Motif III and Motif V of the DEAD helicase family. The structures that the Motif III and Motif V sequences adopt in the DEAD motor open conformation are incompatible with the closed conformation. Therefore, when the DEAD motor makes the transition from open to closed, Motif III and Motif V are forced to change their conformations, which likely functions to regulate passage through the transition state for ATP hydrolysis. The transition state for ATP hydrolysis for the SecA DEAD motor was modeled based on the conformation of the Vasa helicase in complex with adenylyl imidodiphosphate and RNA (Protein Data Bank ID 2DB3). A mechanism for chemical-mechanical coupling emerges, where passage through the transition state for ATP hydrolysis is hindered by the conformational changes required in Motif III and Motif V, and may be promoted by binding interactions with the preprotein substrate and/or other translocase domains and subunits.     

Insights into the chemomechanical coupling of the myosin motor from simulation of its ATP hydrolysis mechanism

The molecular motor myosin converts chemical energy from ATP hydrolysis into mechanical work, thus driving a variety of essential motility processes. Although myosin function has been studied extensively, the catalytic mechanism of ATP hydrolysis and its chemomechanical coupling to the motor cycle are not completely understood. Here, the catalysis mechanism in myosin II is examined using quantum mechanical/molecular mechanical reaction path calculations. The resulting reaction pathways, found in the catalytically competent closed/closed conformation of the Switch-1/Switch-2 loops of myosin, are all associative with a pentavalent bipyramidal oxyphosphorane transition state but can vary in the activation mechanism of the attacking water molecule and in the way the hydrogens are transferred between the heavy atoms. The coordination bond between the Mg2+ metal cofactor and Ser237 in the Switch-1 loop is broken in the product state, thereby facilitating the opening of the Switch-1 loop after hydrolysis is completed, which is required for subsequent strong rebinding to actin. This reveals a key element of the chemomechanical coupling that underlies the motor cycle, namely, the modulation of actin unbinding or binding in response to the ATP or ADP x P(i) state of nucleotide-bound myosin.     

Crystal structure of tryptophanyl-tRNA synthetase complexed with adenosine-5' tetraphosphate: evidence for distributed use of catalytic binding energy in amino acid activation by class I aminoacyl-tRNA synthetases

Tryptophanyl-tRNA synthetase (TrpRS) is a functionally dimeric ligase, which specifically couples hydrolysis of ATP to AMP and pyrophosphate to the formation of an ester bond between tryptophan and the cognate tRNA. TrpRS from Bacillus stearothermophilus binds the ATP analogue, adenosine-5' tetraphosphate (AQP) competitively with ATP during pyrophosphate exchange. Estimates of binding affinity from this competitive inhibition and from isothermal titration calorimetry show that AQP binds 200 times more tightly than ATP both under conditions of induced-fit, where binding is coupled to an unfavorable conformational change, and under exchange conditions, where there is no conformational change. These binding data provide an indirect experimental measurement of +3.0 kcal/mol for the conformational free energy change associated with induced-fit assembly of the active site. Thermodynamic parameters derived from the calorimetry reveal very modest enthalpic changes, consistent with binding driven largely by a favorable entropy change. The 2.5 A structure of the TrpRS:AQP complex, determined de novo by X-ray crystallography, resembles that of the previously described, pre-transition state TrpRS:ATP complexes. The anticodon-binding domain untwists relative to the Rossmann-fold domain by 20% of the way toward the orientation observed for the Products complex. An unexpected tetraphosphate conformation allows the gamma and deltad phosphate groups to occupy positions equivalent to those occupied by the beta and gamma phosphates of ATP. The beta-phosphate effects a 1.11 A extension that relocates the alpha-phosphate toward the tryptophan carboxylate while the PPi mimic moves deeper into the KMSKS loop. This configuration improves interactions between enzyme and nucleotide significantly and uniformly in the adenosine and PPi binding subsites. A new hydrogen bond forms between S194 from the class I KMSKS signature sequence and the PPi mimic. These complementary thermodynamic and structural data are all consistent with the conclusion that the tetraphosphate mimics a transition-state in which the KMSKS loop develops increasingly tight bonds to the PPi leaving group, weakening linkage to the Palpha as it is relocated by an energetically favorable domain movement. Consistent with extensive mutational data on Tyrosyl-tRNA synthetase, this aspect of the mechanism develops high transition-state affinity for the adenosine and pyrophosphate moieties, which move significantly, relative to one another, during the catalytic step.     

Catalytic mechanism of the tryptophan activation reaction revealed by crystal structures of human tryptophanyl-tRNA synthetase in different enzymatic states

Human tryptophanyl-tRNA synthetase (hTrpRS) differs from its bacterial counterpart at several key positions of the catalytic active site and has an extra N-terminal domain, implying possibly a different catalytic mechanism. We report here the crystal structures of hTrpRS in complexes with Trp, tryptophanamide and ATP and tryptophanyl-AMP, respectively, which represent three different enzymatic states of the Trp activation reaction. Analyses of these structures reveal the molecular basis of the mechanisms of the substrate recognition and the activation reaction. The dimeric hTrpRS is structurally and functionally asymmetric with half-of-the-sites reactivity. Recognition of Trp is by an induced-fit mechanism involving conformational change of the AIDQ motif that creates a perfect pocket for the binding and activation of Trp and causes coupled movements of the N-terminal and C-terminal domains. The KMSAS loop appears to have an inherent flexibility and the binding of ATP stabilizes it in a closed conformation that secures the position of ATP for catalysis. Our structural data indicate that the catalytic mechanism of the Trp activation reaction by hTrpRS involves more moderate conformational changes of the structural elements at the active site to recognize and bind the substrates, which is more complex and fine-tuned than that of bacterial TrpRS.     

Electrostatic interactions determine entrance/release order of substrates in the catalytic cycle of adenylate kinase

Adenylate kinase is a monomeric phosphotransferase with important biological function in regulating concentration of adenosine triphosphate (ATP) in cells, by transferring the terminal phosphate group from ATP to adenosine monophosphate (AMP) and forming two adenosine diphosphate (ADP) molecules. During this reaction, the kinase may undergo a large conformational transition, forming different states with its substrates. Although many structures of the protein are available, atomic details of the whole process remain unclear. In this article, we use both conventional molecular dynamics (MD) simulation and an enhanced sampling technique called parallel cascade selection MD simulation to explore different conformational states of the Escherichia coli adenylate kinase. Based on the simulation results, we propose a possible entrance/release order of substrates during the catalytic cycle. The substrate-free protein prefers an open conformation, but changes to a closed state once ATP·Mg enters into its binding pocket first and then AMP does. After the reaction of ATP transferring the terminal phosphate group to AMP, ADP·Mg and ADP are released sequentially, and finally the whole catalyze cycle is completed. Detailed contact and distance analysis reveals that the entrance/release order of substrates may be largely controlled by electrostatic interactions between the protein and the substrates.     

Fine structure of conformational ensembles in adenylate kinase

Adenylate kinase (ADK) catalyzes the reversible Mg²⁺‐dependent phosphoryl transfer reaction Mg²⁺+2ADP ?Mg²⁺+ATP + AMP in essential cellular systems. This reaction is a major player in cellular energy homeostasis and the isoform network of ADK plays an important role in AMP metabolic signaling circuits. ADK has 3 domains, the LID, NMP, and CORE domains, that undergo large conformational rearrangements during ADK's catalytic cycle. In spite of extensive experimental and computational studies, details of the conformational pathway from open to closed forms remain uncertain. In this paper we explore this pathway using coarse‐grained molecular dynamics (MD) trajectories of ADK calculated by GROMACS using a SMOG model and classify the conformations within the resultant trajectories by K‐means clustering. ADK conformations segregate naturally into open; intermediate; and closed forms with long‐term residence in the intermediate state. Structural clustering divides the intermediate conformation into 3 sub‐states that are distinguished from one another on the basis of differences in both structure and dynamics. These distinctions are defined on the basis of a number of different metrics including radius of gyration, dihedral angle fluctuation, and fluctuations of interatomic pair distances. Furthermore, differences in the sub‐states appear to correspond to the distinct ways each sub‐state contributes to the molecular mechanism of catalysis: One sub‐state acts as a gate‐way to the open conformation; one sub‐state a gate‐way to the closed conformation. A third intermediate sub‐state appears to represent a metastable off‐pathway structure that is nevertheless frequently visited during the passage from open to closed state.     

Interconversion of ATP binding and conformational free energies by tryptophanyl-tRNA synthetase: structures of ATP bound to open and closed,pre-transition-state conformations

Binding ATP to tryptophanyl-tRNA synthetase (TrpRS) in a catalytically competent configuration for amino acid activation destabilizes the enzyme structure prior to forming the transition state. This conclusion follows from monitoring the titration of TrpRS with ATP by small angle solution X-ray scattering, enzyme activity, and crystal structures. ATP induces a significantly smaller radius of gyration at pH=7 with a transition midpoint at approximately 8mM. A non-reciprocal dependence of Trp and ATP dissociation constants on concentrations of the second substrate show that Trp binding enhances affinity for ATP, while the affinity for Trp falls with the square of the [ATP] over the same concentration range ( approximately 5mM) that induces the more compact conformation. Two distinct TrpRS:ATP structures have been solved, a high-affinity complex grown with 1mM ATP and a low-affinity complex grown at 10mM ATP. The former is isomorphous with unliganded TrpRS and the Trp complex from monoclinic crystals. Reacting groups of the two individually-bound substrates are separated by 6.7A. Although it lacks tryptophan, the low-affinity complex has a closed conformation similar to that observed in the presence of both ATP and Trp analogs such as indolmycin, and resembles a complex previously postulated to form in the closely-related TyrRS upon induced-fit active-site assembly, just prior to catalysis. Titration of TrpRS with ATP therefore successively produces structurally distinct high- and low-affinity ATP-bound states. The higher quality X-ray data for the closed ATP complex (2.2A) provide new structural details likely related to catalysis, including an extension of the KMSKS loop that engages the second lysine and serine residues, K195 and S196, with the alpha and gamma-phosphates; interactions of the K111 side-chain with the gamma-phosphate; and a water molecule bridging the consensus sequence residue T15 to the beta-phosphate. Induced-fit therefore strengthens active-site interactions with ATP, substantially intensifying the interaction of the KMSKS loop with the leaving PP(i) group. Formation of this conformation in the absence of a Trp analog implies that ATP is a key allosteric effector for TrpRS. The paradoxical requirement for high [ATP] implies that Gibbs binding free energy is stored in an unfavorable protein conformation and can then be recovered for useful purposes, including catalysis in the case of TrpRS.     

Structure of substrate-free human insulin-degrading enzyme (IDE) and biophysical analysis of ATP-induced conformational switch of IDE

Insulin-degrading enzyme (IDE) is a zinc metalloprotease that hydrolyzes amyloid-beta (Abeta) and insulin, which are peptides associated with Alzheimer disease (AD) and diabetes, respectively. Our previous structural analysis of substrate-bound human 113-kDa IDE reveals that the N- and C-terminal domains of IDE, IDE-N and IDE-C, make substantial contact to form an enclosed catalytic chamber to entrap its substrates. Furthermore, IDE undergoes a switch between the closed and open conformations for catalysis. Here we report a substrate-free IDE structure in its closed conformation, revealing the molecular details of the active conformation of the catalytic site of IDE and new insights as to how the closed conformation of IDE may be kept in its resting, inactive conformation. We also show that Abeta is degraded more efficiently by IDE carrying destabilizing mutations at the interface of IDE-N and IDE-C (D426C and K899C), resulting in an increase in Vmax with only minimal changes to Km. Because ATP is known to activate the ability of IDE to degrade short peptides, we investigated the interaction between ATP and activating mutations. We found that these mutations rendered IDE less sensitive to ATP activation, suggesting that ATP might facilitate the transition from the closed state to the open conformation. Consistent with this notion, we found that ATP induced an increase in hydrodynamic radius, a shift in electrophoretic mobility, and changes in secondary structure. Together, our results highlight the importance of the closed conformation for regulating the activity of IDE and provide new molecular details that will facilitate the development of activators and inhibitors of IDE.     

The role of magnesium for geometry and charge in GTP hydrolysis, revealed by quantum mechanics/molecular mechanics simulations

The coordination of the magnesium ion in proteins by triphosphates plays an important role in catalytic hydrolysis of GTP or ATP, either in signal transduction or energy conversion. For example, in Ras the magnesium ion contributes to the catalysis of GTP hydrolysis. The cleavage of GTP to GDP and P(i) in Ras switches off cellular signaling. We analyzed GTP hydrolysis in water, Ras, and Ras·Ras-GTPase-activating protein using quantum mechanics/molecular mechanics simulations. By comparison of the theoretical IR-difference spectra for magnesium ion coordinated triphosphate to experimental ones, the simulations are validated. We elucidated thereby how the magnesium ion contributes to catalysis. It provides a temporary storage for the electrons taken from the triphosphate and it returns them after bond cleavage and P(i) release back to the diphosphate. Furthermore, the Ras·Mg(2+) complex forces the triphosphate into a stretched conformation in which the β- and γ-phosphates are coordinated in a bidentate manner. In this conformation, the triphosphate elongates the bond, which has to be cleaved during hydrolysis. Furthermore, the γ-phosphate adopts a more planar structure, driving the conformation of the molecule closer to the hydrolysis transition state. GTPase-activating protein enhances these changes in GTP conformation and charge distribution via the intruding arginine finger.     

Cross-linking of two beta subunits in the closed conformation in F1-ATPase

In the crystal structure of mitochondrial F1-ATPase, two beta subunits with a bound Mg-nucleotide are in "closed" conformations, whereas the third beta subunit without bound nucleotide is in an "open" conformation. In this "CCO" (beta-closed beta-closed beta-open) conformational state, Ile-390s of the two closed beta subunits, even though they are separated by an intervening alpha subunit, have a direct contact. We replaced the equivalent Ile of the alpha3beta3gamma subcomplex of thermophilic F1-ATPase with Cys and observed the formation of the beta-beta cross-link through a disulfide bond. The analysis of conditions required for the cross-link formation indicates that: (i) F1-ATPase takes the CCO conformation when two catalytic sites are filled with Mg-nucleotide, (ii) intermediate(s) with the CCO conformation are generated during catalytic cycle, (iii) the Mg-ADP inhibited form is in the CCO conformation, and (iv) F1-ATPase dwells in conformational state(s) other than CCO when only one (or none) of catalytic sites is filled by Mg-nucleotide or when catalytic sites are filled by Mg2+-free nucleotide. The alpha3beta3gamma subcomplex containing the beta-beta cross-link retained the activity of uni-site catalysis but lost that of multiple catalytic turnover, suggesting that open-closed transition of beta subunits is required for the rotation of gamma subunit but not for hydrolysis of a single ATP.     

ATP hydrolysis in the betaTP and betaDP catalytic sites of F1-ATPase

The enzyme F1-adenosine triphosphatase (ATPase) is a molecular motor that converts the chemical energy stored in the molecule adenosine triphosphate (ATP) into mechanical rotation of its gamma-subunit. During steady-state catalysis, the three catalytic sites of F1 operate in a cooperative fashion such that at every instant each site is in a different conformation corresponding to a different stage along the catalytic cycle. Notwithstanding a large amount of biochemical and, recently, structural data, we still lack an understanding of how ATP hydrolysis in F1 is coupled to mechanical motion and how the catalytic sites achieve cooperativity during rotatory catalysis. In this publication, we report combined quantum mechanical/molecular mechanical simulations of ATP hydrolysis in the betaTP and betaDP catalytic sites of F1-ATPase. Our simulations reveal a dramatic change in the reaction energetics from strongly endothermic in betaTP to approximately equienergetic in betaDP. The simulations identify the responsible protein residues, the arginine finger alphaR373 being the most important one. Similar to our earlier study of betaTP, we find a multicenter proton relay mechanism to be the energetically most favorable hydrolysis pathway. The results elucidate how cooperativity between catalytic sites might be achieved by this remarkable molecular motor.     

Targeted molecular dynamics simulation studies of binding and conformational changes in E. coli MurD

Enzymes involved in the biosynthesis of bacterial peptidoglycan, an essential cell wall polymer unique to prokaryotic cells, represent a highly interesting target for antibacterial drug design. Structural studies of E. coli MurD, a three-domain ATP hydrolysis driven muramyl ligase revealed two inactive open conformations of the enzyme with a distinct C-terminal domain position. It was hypothesized that the rigid body rotation of this domain brings the enzyme to its closed active conformation, a structure, which was also determined experimentally. Targeted molecular dynamics 1 ns-length simulations were performed in order to examine the substrate binding process and gain insight into structural changes in the enzyme that occur during the conformational transitions into the active conformation. The key interactions essential for the conformational transitions and substrate binding were identified. The results of such studies provide an important step toward more powerful exploitation of experimental protein structures in structure-based inhibitor design.     

Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials

The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.     

Guanylate kinase, induced fit, and the allosteric spring probe

Since the introduction of the induced-fit theory by D. E. Koshland Jr., it has been established that conformational motion invariably accompanies the execution of protein function. The catalytic activity of kinases, specifically, is associated with large conformational changes (∼1 nm amplitude). In the case of guanylate kinase, upon substrate binding, the LID and nucleotide-monophosphate-binding domains are brought together and toward the CORE with large concerted movements about the α3 (helix 3) axis. However, whether the change in conformation mostly affects the catalytic rate or mostly increases binding affinities for one or the other substrate is unclear. We investigate this question using a nanotechnology approach based on mechanical stress. Using an “allosteric spring probe”, we bias conformational states in favor of the “open” (substrate-free) conformation of the enzyme; the result is that the binding constant for the substrate guanosine monophosphate (GMP) is reduced by up to a factor of 10, whereas the binding constant for adenosine triphosphate (ATP) and the catalytic rate are essentially unaffected. The results show that the GMP-induced conformational change, which promotes catalysis, does not promote ATP binding, consistent with previous mutagenesis studies. Furthermore, they show that this conformational change is of the induced-fit type with respect to GMP binding (but not ATP binding). We elaborate on this point by proposing a quantitative criterion for the classification of conformational changes with respect to the induced-fit theory. More generally, these results show that the allosteric spring probe can be used to affect enzymatic activity in a continuously controlled manner, and also to affect specific steps of the reaction mechanism while leaving others unaffected. It is presumed that this will enable informative comparisons with the results of future molecular dynamics or statistical mechanics computations.     

The dynamics of the MgATP-driven closure of MalK, the energy-transducing subunit of the maltose ABC transporter

The nucleotide binding domains (NBDs) are the energy supplying subunits of ATP-binding cassette (ABC) proteins. They power transport by binding and hydrolyzing ATP. Tracing the pathway between different conformational states of the NBDs during ATP binding, hydrolysis, and release has, however, proven difficult. We have used molecular dynamics simulations to study the ATP-driven association of the NBDs of the maltose ABC transporter, MalK, based on the crystal structures of its open and semiopen dimers. When MgATP was introduced into the binding pockets, the semiopen dimer transitioned to a closed conformation, whereas the open dimer evolved to a semiopen state. In the absence of docked MgATP, however, the twin NBDs of both the open and semiopen starting configurations drifted further apart. Both the presence of MgATP and direct cross-interface protein-protein hydrogen bonds, primarily involving the D-loop, quite likely play a key role in initiating closure. The simulations of the MgATP-docked semiopen form indicate that completion of closure is driven mainly by cross-interface contacts between the gamma-phosphate of ATP and residues in the signature motif. Our simulations also give insight into possible interactions of MalK with the regulatory proteins MalT and enzyme IIA(glc).     

Comparative structural dynamics of Tyrosyl-tRNA synthetase complexed with different substrates explored by molecular dynamics

Several molecular dynamics simulations of S. aureus Tyrosyl-tRNA synthetase (TyrRS) in its free form and complexed with Tyr, ATP, tyrosyl adenylate and inhibitor respectively have been carried out to investigate the ligand-linked conformational stability changes associated with its catalytic cycle. The results show that unliganded S. aureus TyrRS samples a more relaxed conformational space than substrate-bound TyrRS. There are three high flexibility regions encompassing residues 114–118, 128–133, and 226–238 respectively. The region which includes the KMSKS motif (KFGKS in S. aureus TyrRS) shows the highest difference in fluctuations. Hydrogen bond network formed by Tyr, ATP, tyrosyl adenylate and inhibitor with S. aureus TyrRS is discussed. Our simulations suggest the induced-fit conformational changes of the KMSKS loop as follows: the KMSKS loop of substrate-free S. aureus TyrRS adopts an open conformation. The tyrosine binds in the pocket with the KMSKS loop balancing between semi-open and open forms. The ATP binding induces the KMSKS loop to the open form. After the Tyr-AMP is formed, the first two residues of KMSKS loop twists in an anticlockwise direction and drives the loop in a conformation similar to the closed one, while those of the last three GKS residues adopt a conformation between semi-open and open conformation. This conformational change may probably be involved in the initial tRNA binding. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Briefly bound to activate: transient binding of a second catalytic magnesium activates the structure and dynamics of CDK2 kinase for catalysis

We have determined high-resolution crystal structures of a CDK2/Cyclin A transition state complex bound to ADP, substrate peptide, and MgF(3)(-). Compared to previous structures of active CDK2, the catalytic subunit of the kinase adopts a more closed conformation around the active site and now allows observation of a second Mg(2+) ion in the active site. Coupled with a strong [Mg(2+)] effect on in vitro kinase activity, the structures suggest that the transient binding of the second Mg(2+) ion is necessary to achieve maximum rate enhancement of the chemical reaction, and Mg(2+) concentration could represent an important regulator of CDK2 activity in vivo. Molecular dynamics simulations illustrate how the simultaneous binding of substrate peptide, ATP, and two Mg(2+) ions is able to induce a more rigid and closed organization of the active site that functions to orient the phosphates, stabilize the buildup of negative charge, and shield the subsequently activated γ-phosphate from solvent.     

A conformational transition state accompanies tryptophan activation by B. stearothermophilus tryptophanyl-tRNA synthetase

B. stearothermophilus tryptophanyl-tRNA synthetase catalysis proceeds via high-energy protein conformations. Unliganded MD trajectories of the pretransition-state complex with Mg(2+)ATP and the (post) transition-state analog complex with adenosine tetraphosphate relax rapidly in opposite directions, the former regressing, the latter progressing along the structural reaction coordinate. The two crystal structures (rmsd 0.7 A) therefore lie on opposite sides of a conformational free-energy maximum as the chemical transition state forms. SNAPP analysis illustrates the complexity of the associated long-range conformational coupling. Switching interactions in four nonpolar core regions are locally isoenergetic throughout the transition. Different configurations, however, propagate their effects to unfavorable, longer-range interactions at the molecular surface. Designed mutation shows that switching interactions enhance the rate, perhaps by destabilizing the ground state immediately before the transition state and limiting nonproductive diffusion before and after the chemical transition state, thereby reducing the activation entropy. This paradigm may apply broadly to energy-transducing enzymes.