Analysis of protein kinase autophosphorylation using expressed protein ligation and computational modeling

Protein kinases represent a family of enzymes that are critical in cell signaling. One mechanism by which protein kinases are regulated is via autophosphorylation. In the studies described here, we have examined the mechanism of autophosphosphorylation at serine 338 in the regulation of protein kinase A (PKA). Expressed protein ligation allowed for the covalent linkage of an ATP moiety to a Ser mimic at this phosphorylation site. Using a combination of size exclusion chromatography, fluorescence nucleotide binding, kinase measurements, and limited proteolysis assays on this semisynthetic ATP-linked protein, we have obtained unique evidence for an intramolecular autophosphorylation mechanism in PKA regulation. Computational analysis provided a plausible model for a PKA conformation consistent with intramolecular phosphoryl transfer. This approach could be applied to other autoprocessing enzymes by exploiting appropriate transition state analogue motifs in the context of protein semisynthesis.     

Geometrically associative yet electronically dissociative character in the transition state of enzymatic reversible phosphorylation

Reversible phosphorylation of proteins is a post‐translational modification that regulates diverse biological processes. The molecular mechanism underlying phosphoryl transfer catalyzed by enzymes remains a subject of active debate. In particular, the nature of transition state (TS), whether it has an associative or dissociative character, has been one of the most controversial issues. Structural evidence supports associative TS, whereas physical organic studies point to a dissociative character. Here we perform hybrid quantum mechanics/molecular mechanics simulations for the reversible phosphorylation of phosphoserine phosphatase (PSP) to study the nature of the TS. Both phosphorylation and dephosphorylation reactions are investigated based on the two‐dimensional energy surfaces along phosphoryl and proton transfer coordinates. The structures of the active site at TS in both reactions reveal compact geometries, consistent with crystal structures of PSP with phosphate analogues. However, the electron density of the phosphoryl group in both TS structures slightly decreases compared with that in the reactant states. These findings suggest that the TS of PSP has a geometrically associative yet electronically dissociative character and strongly depends on proton transfer being coupled with phosphoryl transfer. Structure and literature database, which searches on phosphotransferases, suggest that such a hybrid TS is consistent with many structures and physical organic studies and likely holds for most enzymes catalyzing phosphoryl transfer. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Density functional study of the enzymatic reaction catalyzed by a cyclin-dependent kinase

Density functional theory (DFT) calculations were carried out to study the molecular mechanism of the phosphoryl transfer reaction catalyzed by cyclin-dependent kinases (CDKs). The DFT study presented here shows that CDKs catalyze the phosphoryl transfer reaction from ATP to the serine substrate through a single step mechanism with a SN2-like transition state.     

Active Site Model for Phosphoryl Transfer Enzymes Via Hexacoordination

Previous work has demonstrated the ease with which phosphorus can increase its coordination geometry. The present study has more closely modeled active sites of phosphoryl transfer enzymes by the inclusion of anionic phosphorus, which allows for oxygen atom donor coordination at phosphorus in the presence of a hydrogen bonding network. The resulting increase in phosphorus-oxygen donor coordination compared to analogous systems containing neutral phosphorus compounds serves as a model applicable to proposed mechanisms at active sites of phosphoryl transfer enzymes.     

Theoretical comparison of p-nitrophenyl phosphate and sulfate hydrolysis in aqueous solution: implications for enzyme-catalyzed sulfuryl transfer

Both phosphoryl and sulfuryl transfers are ubiquitous in biology, being involved in a wide range of processes, ranging from cell division to apoptosis. Additionally, it is becoming increasingly clear that enzymes that can catalyze phosphoryl transfer can often cross-catalyze sulfuryl transfer (and vice versa). However, while there have been extensive experimental and theoretical studies performed on phosphoryl transfer, the body of available research on sulfuryl transfer is comparatively much smaller. The present work presents a direct theoretical comparison of p-nitrophenyl phosphate and sulfate monoester hydrolysis, both of which are considered prototype systems for probing phosphoryl and sulfuryl transfer, respectively. Specifically, free energy surfaces have been generated using density functional theory, by initial geometry optimization in PCM using the 6-31+G* basis set and the B3LYP density functional, followed by single-point calculations using the larger 6-311+G** basis set and the COSMO continuum model. The resulting surfaces have been then used to identify the relevant transition states, either by further unconstrained geometry optimization or from the surface itself where possible. Additionally, configurational entropies were evaluated using a combination of the quasiharmonic approximation and the restraint release approach and added to the calculated activation barriers as a correction. Finally, the overall activation entropy was estimated by approximating the solvent contribution to the total activation entropy using the Langevin dipoles solvation model. We have reproduced both the experimentally observed activation barriers and the observed trend in the activation entropies with reasonable accuracy, as well as providing a comparison of calculated and observed (15)N and (18)O kinetic isotope effects. We demonstrate that, counterintuitively, the hydrolysis of the p-nitrophenyl sulfate proceeds through a more expansive pathway than its phosphate analogue. Additionally, we show that the solvation effects upon moving from the ground state to the transition state are quite different for both reactions, suggesting that the enzymes that catalyze these reactions would need active sites with quite different electrostatic preorganization for the efficient catalysis of either reaction (despite which many enzymes can catalyze both phosphoryl and sulfuryl transfer). We believe that such a comparative study is an important foundation for understanding the molecular basis for phosphate-sulfate cross-promiscuity within members of the alkaline phosphatase superfamily.     

Visualizing an Ultra‐Weak Protein–Protein Interaction in Phosphorylation Signaling

Proteins interact with each other to fulfill their functions. The importance of weak protein–protein interactions has been increasingly recognized. However, owing to technical difficulties, ultra‐weak interactions remain to be characterized. Phosphorylation can take place via a K_D≈25 mM interaction between two bacterial enzymes. Using paramagnetic NMR spectroscopy and with the introduction of a novel Gd^III‐based probe, we determined the structure of the resulting complex to atomic resolution. The structure accounts for the mechanism of phosphoryl transfer between the two enzymes and demonstrates the physical basis for their ultra‐weak interaction. Further, molecular dynamics (MD) simulations suggest that the complex has a lifetime in the micro‐ to millisecond regimen. Hence such interaction is termed a fleeting interaction. From mathematical modeling, we propose that an ultra‐weak fleeting interaction enables rapid flux of phosphoryl signal, providing a high effective protein concentration.     

Insight into the phosphoryl transfer of the Escherichia coli glucose phosphotransferase system from QM/MM simulations

Phosphoryl transfer is a key reaction in many aspects of metabolism, gene regulation, and signal transduction. One prominent example is the phosphoenolpyruvate:sugar phosphotransferase system (PTS), which represents an integral part of the bacterial sugar metabolism. The transfer between the enzymes IIA (Glc) and IIB (Glc) in the glucose-specific branch of the PTS is of particular interest due to the unusual combination of donor and acceptor residues involved in phosphoryl transfer: The phosphoryl group is initially attached to the Nepsilon2 atom of His 90 in IIA (Glc) and then transferred to the S gamma atom of Cys 35 in IIB (Glc). To gain insight into the details of the transfer mechanism, we have performed a QM/MM simulation which treats the entire active site quantum-mechanically. The transfer has a high dissociative character, and the Nepsilon2-P bond gets immediately destabilized after complex formation by numerous interactions formed between residues of IIB (Glc) and the phosphoryl group. The final formation of a tight S gamma-P bond is accompanied by a reorientation of the side chain of the phosphoryl donor. This reorientation results in the loss of interaction between the imidazole ring and the phosphate group thus hindering the reverse transfer. A comparison to the transfer in protein tyrosine phosphatases, which also use a cysteine as acceptor of the phosphoryl group, reveals significant similarities in the conformation of the active site, the energy profile of the reaction, and in the pattern of interactions that stabilize the phosphoryl group during the transfer.     

The role of the C8 proton of ATP in the catalysis of shikimate kinase and adenylate kinase

ABSTRACT: BACKGROUND: It has been demonstrated that the adenyl moiety of ATP plays a direct role in the regulation of ATP binding and/or phosphoryl transfer within a range of kinase and synthetase enzymes. The role of the C8-H of ATP in the binding and/or phosphoryl transfer on the enzyme activity of a number of kinase and synthetase enzymes has been elucidated. The intrinsic catalysis rate mediated by each kinase enzyme is complex, yielding apparent KM values ranging from less than 0.4 muM to more than 1 mM for ATP in the various kinases. Using a combination of ATP deuterated at the C8 position (C8D-ATP) as a molecular probe with site directed mutagenesis (SDM) of conserved amino acid residues in shikimate kinase and adenylate kinase active sites, we have elucidated a mechanism by which the ATP C8-H is induced to be labile in the broader kinase family. We have demonstrated the direct role of the C8-H in the rate of ATP consumption, and the direct role played by conserved Thr residues interacting with the C8-H. The mechanism by which the vast range in KM might be achieved is also suggested by these findings. RESULTS: We have demonstrated the mechanism by which the enzyme activities of Group 2 kinases, shikimate kinase (SK) and adenylate kinase 1 (AK1), are controlled by the C8-H of ATP. Mutations of the conserved threonine residues associated with the labile C8-H cause the enzymes to lose their saturation kinetics over the concentration range tested. The relationship between the role C8-H of ATP in the reaction mechanism and the ATP concentration as they influence the saturation kinetics of the enzyme activity is also shown. The SDM clearly identified the amino acid residues involved in both the catalysis and regulation of phosphoryl transfer in SK and AK1 as mediated by C8H-ATP. CONCLUSIONS: The data outlined serves to demonstrate the "push" mechanism associated with the control of the saturation kinetics of Group 2 kinases mediated by ATP C8-H. It is therefore conceivable that kinase enzymes achieve the observed 2,500-fold variation in KM through a combination of the various conserved "push" and "pull" mechanisms associated with the release of C8-H, the proton transfer cascades unique to the class of kinase in question and the resultant/concomitant creation of a pentavalent species from the gamma-phosphate group of ATP. Also demonstrated is the interplay between the role of the C8-H of ATP and the ATP concentration in the observed enzyme activity. The lability of the C8-H mediated by active site residues co-ordinated to the purine ring of ATP therefore plays a significant role in explaining the broad KM range associated with kinase steady state enzyme activities.     

Metal Fluorides as Analogues for Studies on Phosphoryl Transfer Enzymes

The 1994 structure of a transition‐state analogue with AlF₄⁻ and GDP complexed to G1α, a small G protein, heralded a new field of research into the structure and mechanism of enzymes that manipulate the transfer of phosphoryl (PO₃⁻) groups. The number of enzyme structures in the PDB containing metal fluorides (MF_x) as ligands that imitate either a phosphoryl or a phosphate group was 357 at the end of 2016. They fall into three distinct geometrical classes: 1) Tetrahedral complexes based on BeF₃⁻ that mimic ground‐state phosphates; 2) octahedral complexes, primarily based on AlF₄⁻, which mimic “in‐line” anionic transition states for phosphoryl transfer; and 3) trigonal bipyramidal complexes, represented by MgF₃⁻ and putative AlF₃⁰ moieties, which mimic the geometry of the transition state. The interpretation of these structures provides a deeper mechanistic understanding into the behavior and manipulation of phosphate monoesters in molecular biology. This Review provides a comprehensive overview of these structures, their uses, and their computational development.     

Biologically relevant phosphoranes: structural characterization of glucofuranose- and xylofuranose-based phosphoranes

Carbohydrate-based phosphoranes were synthesized by reacting 2,2'-ethylidenebis(4,6-di-tert-butylphenyl)fluorophosphite with 1,2-O-isopropylidene-alpha-D-glucofuranose, beta-chloralose, and 1,2-isopropylidene-alpha-D-xylofuranose to form the monocyclic biophosphoranes 1-3, respectively, in the presence of N-chlorodiisopropylamine. Synthesis of the monocyclic biophosphorane 4 was achieved by reacting tris(2,6-di-isopropylphenyl)phosphite with 1,2-O-isopropylidene-alpha-D-glucofuranose in the presence of N-chlorodiisopropylamine. X-ray analysis of 1-4 revealed trigonal bipyramidal structures with the carbohydrate components occupying axial-equatorial sites. An eight-membered ring in 1-3 occupied diequatorial sites of the trigonal bipyramid. Solution and solid state 31P and solution 19F, 1H, and 13C NMR measurements including variable temperature and correlation spectroscopy studies established retention of the solid state structure in solution. A dynamic equilibrium exists among two isomeric forms. These biophosphoranes serve as models for active sites of phosphoryl transfer enzymes. The rapid exchange process reorients the carbohydrate component of the trigonal bipyramidal phosphorane. At an active site, this type of pseudorotational behavior provides a mechanism that could bring another active site residue into play and account for a means by which some phosphoryl transfer enzymes express promiscuous behavior.     

Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis

The molecular mechanism of hairpin ribozyme catalysis is studied with molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential with a recently developed semiempirical AM1/d-PhoT model for phosphoryl transfer reactions. Simulations are used to derive one- and two-dimensional potentials of mean force to examine specific reaction paths and assess the feasibility of proposed general acid and base mechanisms. Density-functional calculations of truncated active site models provide complementary insight to the simulation results. Key factors utilized by the hairpin ribozyme to enhance the rate of transphosphorylation are presented, and the roles of A38 and G8 as general acid and base catalysts are discussed. The computational results are consistent with available experimental data, provide support for a general acid/base mechanism played by functional groups on the nucleobases, and offer important insight into the ability of RNA to act as a catalyst without explicit participation by divalent metal ions.     

The role of the C8 proton of ATP in the regulation of phosphoryl transfer within kinases and synthetases

Background

The kinome comprises functionally diverse enzymes, with the current classification indicating very little about the extent of conserved regulatory mechanisms associated with phosphoryl transfer. The apparent K _m of the kinases ranges from less than 0.4 μM to in excess of 1000 μM for ATP. It is not known how this diverse range of enzymes mechanistically achieves the regulation of catalysis via an affinity range for ATP varying by three-orders of magnitude.

Results

We have demonstrated a previously undiscovered mechanism in kinase and synthetase enzymes where the overall rate of reaction is regulated via the C8-H of ATP. Using ATP deuterated at the C8 position (C8D-ATP) as a molecular probe it was shown that the C8-H plays a direct role in the regulation of the overall rate of reaction in a range of kinase and synthetase enzymes. Using comparative studies on the effect of the concentration of ATP and C8D-ATP on the activity of the enzymes we demonstrated that not only did C8D-ATP give a kinetic isotope effect (KIE) but the KIE's obtained are clearly not secondary KIE effects as the magnitude of the KIE in all cases was at least 2 fold and in most cases in excess of 7 fold.

Conclusions

Kinase and synthetase enzymes utilise C8D-ATP in preference to non-deuterated ATP. The KIE obtained at low ATP concentrations is clearly a primary KIE demonstrating strong evidence that the bond to the isotopically substituted hydrogen is being broken. The effect of the ATP concentration profile on the KIE was used to develop a model whereby the C8H of ATP plays a role in the overall regulation of phosphoryl transfer. This role of the C8H of ATP in the regulation of substrate binding appears to have been conserved in all kinase and synthetase enzymes as one of the mechanisms associated with binding of ATP. The induction of the C8H to be labile by active site residues coordinated to the ATP purine ring may play a significant role in explaining the broad range of K _m associated with kinase enzymes.     

Mechanistic Insights into RNA Transphosphorylation from Kinetic Isotope Effects and Linear Free Energy Relationships of Model Reactions

Phosphoryl transfer reactions are ubiquitous in biology and the understanding of the mechanisms whereby these reactions are catalyzed by protein and RNA enzymes is central to reveal design principles for new therapeutics. Two of the most powerful experimental probes of chemical mechanism involve the analysis of linear free energy relations (LFERs) and the measurement of kinetic isotope effects (KIEs). These experimental data report directly on differences in bonding between the ground state and the rate‐controlling transition state, which is the most critical point along the reaction free energy pathway. However, interpretation of LFER and KIE data in terms of transition‐state structure and bonding optimally requires the use of theoretical models. In this work, we apply density‐functional calculations to determine KIEs for a series of phosphoryl transfer reactions of direct relevance to the 2′‐O‐transphosphorylation that leads to cleavage of the phosphodiester backbone of RNA. We first examine a well‐studied series of phosphate and phosphorothioate mono‐, di‐ and triesters that are useful as mechanistic probes and for which KIEs have been measured. Close agreement is demonstrated between the calculated and measured KIEs, establishing the reliability of our quantum model calculations. Next, we examine a series of RNA transesterification model reactions with a wide range of leaving groups in order to provide a direct connection between observed Brønsted coefficients and KIEs with the structure and bonding in the transition state. These relations can be used for prediction or to aid in the interpretation of experimental data for similar non‐enzymatic and enzymatic reactions. Finally, we apply these relations to RNA phosphoryl transfer catalyzed by ribonuclease A, and demonstrate the reaction coordinate–KIE correlation is reasonably preserved. A prediction of the secondary deuterium KIE in this reaction is also provided. These results demonstrate the utility of building up knowledge of mechanism through the systematic study of model systems to provide insight into more complex biological systems such as phosphoryl transfer enzymes and ribozymes.     

Intermolecular Phosphoryl Transfer Between Serine and Histidine Residues

A novel intermolecular phosphoryl transfer from O-trimethylsilyl-N-(O, O-diisopropyl)phosphoryl serine trimethylsilyl ester to N, N‘-bis(trimethylsilyl) histidine trimethylsilyl ester wasstudied through electrospray ionization mass spectrometry (ESI-MS). It was proposed that thetransfer reaction went through penta-coordinated phosphorus intermediate.     

The Self-catalytic Esterification Reaction of O-Phosphoryl Serine Derivative

O-Phosphoryl serine derivative can perform serf-catalytic esterification reaction in the mixture of CH3OH and CHCl3 at the room temperature. The phosphoryl group participation was the key step of the esterification. This type of reactions were proposed through an intermediate of mixed phosphoric-carboxylic anhydride that might provide a clue to the function of the phosphoryl group in the phosphorylated enzymes and in the prebiotic synthesis of protein.     

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 1. Implementation

The present paper describes the extension of a recently developed smooth conductor-like screening model for solvation to a d-orbital semiempirical framework (MNDO/d-SCOSMO) with analytic gradients that can be used for geometry optimizations, transition state searches, and molecular dynamics simulations. The methodology is tested on the potential energy surfaces for separating ions and the dissociative phosphoryl transfer mechanism of methyl phosphate. The convergence behavior of the smooth COSMO method with respect to discretization level is examined and the numerical stability of the energy and gradient are compared to that from conventional COSMO calculations. The present method is further tested in applications to energy minimum and transition state geometry optimizations of neutral and charged metaphosphates, phosphates, and phosphoranes that are models for stationary points in transphosphorylation reaction pathways of enzymes and ribozymes. The results indicate that the smooth COSMO method greatly enhances the stability of quantum mechanical geometry optimization and transition state search calculations that would routinely fail with conventional solvation methods. The present MNDO/d-SCOSMO method has considerable computational advantages over hybrid quantum mechanical/molecular mechanical methods with explicit solvation, and represents a potentially useful tool in the arsenal of multi-scale quantum models used to study biochemical reactions.     

Computational study of the phosphoryl transfer catalyzed by a cyclin-dependent kinase

A cyclin-dependent kinase, Cdk2, catalyzes the transfer of the gamma-phosphate from ATP to a threonine or serine residue of its polypeptide substrates. Here, we investigate aspects of the reaction mechanism of Cdk2 by gas-phase density functional calculations, classical molecular dynamics, and Car-Parrinello QM/MM simulations. We focus on the role of the conserved Asp127 and on the nature of the phosphoryl transfer reaction mechanism catalyzed by Cdk2. Our findings suggest that Asp127 is active in its deprotonated form by assisting the formation of the near-attack orientation of the substrate serine or threonine. Therefore, the residue does not act as a general base during the catalysis. The mechanism for the phosphoryl transfer is a single SN2-like concerted step, which shows a phosphorane-like transition state geometry. Although the resulting reaction mechanism is in agreement with a previous density functional study of the same catalytic reaction mechanism (Cavalli et al., Chem. Comm. 2003, 1308-1309), the reaction barrier is considerably lower when QM/MM calculations are performed, as in this study ( approximately 42 kcal mol(-1) QM vs. approximately 24 kcal mol(-1) QM/MM); this indicates that important roles for the catalysis are played by the protein environment and solvent waters. Because of the high amino acid sequence conservation among the whole family of cyclin-dependent kinases (CDKs), these results could be general for the CDK family.     

Phosphoryl Transfer Enzymes and Hypervalent Phosphorus Chemistry: A Keynote Lecture of the 7th International Conference on Heteroatom Chemistry

The coordination tendencies of phosphorus to form a hexacoordinated state from a pentacoordinated state which might assist in describing mechanistic action of phosphoryl transfer enzymes are delineated. The factors discussed include substrate and transition or intermediate state anionicity; hydrogen bonding; packing effects, i.e., van der Waals forces; the ease of formation of hexacoordinate phosphorus from lower coordinate states; and the pseudorotation problem common to nonrigid pentacoordinate phosphorus. In view of the work reported in this account and recent work on enzyme promiscuity and moonlighting activities, it is suggested that donor action should play a role in determining active-site interactions in phosphoryl transfer enzyme mechanisms.     

Thiophosphoryl fluoride and phosphoryl fluoride as initiators for the polymerisation of tetrahydrofuran

Thiophosphoryl fluoride and phosphoryl fluoride have been found to initiate the polymerisation of tetrahydrofuran. The living polymer formed has a high molecular weight of the order of a million and the density is found to be between 0.98 - 1.02 g/cc. A cationic mechanism for the polymerisation has been proposed.     

Benchmark calculations of proton affinities and gas-phase basicities of molecules important in the study of biological phosphoryl transfer

Benchmark calculations of proton affinities and gas-phase basicities of molecules most relevant to biological phosphoryl transfer reactions are presented and compared with available experimental results. The accuracy of proton affinity and gas-phase basicity results obtained from several multi-level model chemistries (CBS-QB3, G3B3, and G3MP2B3) and density-functional quantum models (PBE0, B1B95, and B3LYP) are assessed and compared. From these data, a set of empirical bond enthalpy, entropy, and free energy corrections are introduced that considerably improve the accuracy and predictive capability of the methods. These corrections are applied to the prediction of proton affinity and gas-phase basicity values of important biological phosphates and phosphoranes for which experimental data does not currently exist. Comparison is made with results from semiempirical quantum models that are commonly employed in hybrid quantum mechanical/molecular mechanical simulations. Data suggest that the design of improved semiempirical quantum models with increased accuracy for relative proton affinity values is necessary to obtain quantitative accuracy for phosphoryl transfer reactions in solution, enzymes, and ribozymes.