Kinetic and catalytic mechanisms of protein kinases

Catalytic mechanisms and regulation of protein kinases

Protein kinases transfer a phosphoryl group from ATP onto target proteins and play a critical role in signal transduction and other cellular processes. Here, we review the kinase kinetic and chemical mechanisms and their application in understanding kinase structure and function. Aberrant kinase activity has been implicated in many human diseases, in particular cancer. We highlight applications of technologies and concepts derived from kinase mechanistic studies that have helped illuminate how kinases are regulated and contribute to pathophysiology.

A structural atlas of kinases inhibited by clinically approved drugs

The aberrant activation of protein kinases is associated with many human diseases, most notably cancer. Due to this link between kinase deregulation and disease progression, kinases are one of the most targeted protein families for small-molecule inhibition. Within the last 15 years, the U.S. Food and Drug Administration has approved over 20 small-molecule inhibitors of protein kinases for use in the clinic. These inhibitors target the kinase active site and represent the successful hurdling by medicinal chemists of the formidable challenge posed by the high similarity among the active sites of the approximately 500 human kinases. We review the conserved structural features of kinases that are important for inhibitor binding as well as for catalysis. Many clinically approved drugs elicit selectivity by exploiting subtle variation within the kinase active site. We highlight some of the crystallographic studies on the kinase-inhibitor complexes that have provided valuable guidance for the development of these drugs as well as for future drug design efforts.

Discovery of nitroaryl urea derivatives with antiproliferative properties

A series of urea derivatives bearing nitroaryl moiety has been synthesized and assayed for their potential antiproliferative activities. Some of the tested compounds displayed activity in RK33 laryngeal cancer cells and TE671 rhabdomyosarcoma cells while being generally less toxic to healthy HSF human fibroblasts cells. One compound was demonstrated to be a moderate CDK2 inhibitor with IC50 = 14.3 µM. Its structure was solved by an X-ray crystallography and molecular modelling was performed to determine structure-activity relationship. Obtained compounds constitute novel structures and generally demonstrated greater cytotoxicity in comparison to cisplatin. This study offers new structural motifs with potential for further development.

Quantitative Circadian Phosphoproteomic Analysis of Arabidopsis Reveals Extensive Clock Control of Key Components in Physiological, Metabolic, and Signaling Pathways

The circadian clock provides adaptive advantages to an organism, resulting in increased fitness and survival. The phosphorylation events that regulate circadian-dependent signaling and the processes which post-translationally respond to clock-gated signals are largely unknown. To better elucidate post-translational events tied to the circadian system we carried out a survey of circadian-regulated protein phosphorylation events in Arabidopsis seedlings. A large-scale mass spectrometry-based quantitative phosphoproteomics approach employing TiO2-based phosphopeptide enrichment techniques identified and quantified 1586 phosphopeptides on 1080 protein groups. A total of 102 phosphopeptides displayed significant changes in abundance, enabling the identification of specific patterns of response to circadian rhythms. Our approach was sensitive enough to quantitate oscillations in the phosphorylation of low abundance clock proteins (early flowering4; ELF4 and pseudoresponse regulator3; PRR3) as well as other transcription factors and kinases. During constant light, extensive cyclic changes in phosphorylation status occurred in critical regulators, implicating direct or indirect regulation by the circadian system. These included proteins influencing transcriptional regulation, translation, metabolism, stress and phytohormones-mediated responses. We validated our analysis using the elf4-211 allele, in which an S45L transition removes the phosphorylation herein identified. We show that removal of this phosphorylatable site diminishes interaction with early flowering3 (ELF3), a key partner in a tripartite evening complex required for circadian cycling. elf4-211 lengthens period, which increases with increasing temperature, relative to the wild type, resulting in a more stable temperature compensation of circadian period over a wider temperature range.

Potential-based dynamical reweighting for Markov state models of protein dynamics

As simulators attempt to replicate the dynamics of large cellular components in silico, problems related to sampling slow, glassy degrees of freedom in molecular systems will be amplified manyfold. It is tempting to augment simulation techniques with external biases to overcome such barriers with ease; biased simulations, however, offer little utility unless equilibrium properties of interest (both kinetic and thermodynamic) can be recovered from the data generated. In this Article, we present a general scheme that harnesses the power of Markov state models (MSMs) to extract equilibrium kinetic properties from molecular dynamics trajectories collected on biased potential energy surfaces. We first validate our reweighting protocol on a simple two-well potential, and we proceed to test our method on potential-biased simulations of the Trp-cage miniprotein. In both cases, we find that equilibrium populations, time scales, and dynamical processes are reliably reproduced as compared to gold standard, unbiased data sets. We go on to discuss the limitations of our dynamical reweighting approach, and we suggest auspicious target systems for further application.

Structural basis of nSH2 regulation and lipid binding in PI3Kα

We report two crystal structures of the wild-type phosphatidylinositol 3-kinase α (PI3Kα) heterodimer refined to 2.9 Å and 3.4 Å resolution: the first as the free enzyme, the second in complex with the lipid substrate, diC4-PIP₂, respectively. The first structure shows key interactions of the N-terminal SH2 domain (nSH2) and iSH2 with the activation loop that suggest a mechanism by which the enzyme is inhibited in its basal state. In the second structure, the lipid substrate binds in a positively charged pocket adjacent to the ATP-binding site, bordered by the P-loop, the activation loop and the iSH2 domain. An additional lipid-binding site was identified at the interface of the ABD, iSH2 and kinase domains. The ability of PI3Kα to bind an additional PIP₂ molecule was confirmed in vitro by fluorescence quenching experiments. The crystal structures reveal key differences in the way the nSH2 domain interacts with wild-type p110α and with the oncogenic mutant p110αH1047R. Increased buried surface area and two unique salt-bridges observed only in the wild-type structure suggest tighter inhibition in the wild-type PI3Kα than in the oncogenic mutant. These differences may be partially responsible for the increased basal lipid kinase activity and increased membrane binding of the oncogenic mutant.

Synthesis, characterization, biological evaluation and docking studies of macrocyclic binuclear manganese(II) complexes containing 3,5-dinitrobenzoyl pendant arms

A series of bis(phenoxo) bridged binuclear manganese(II) complexes of the type [Mn2L(1-3)](ClO4)2 (1-3) containing 3,5-dinitrobenzoyl pendant-arms have been synthesized by cyclocondensation of 2,6-diformyl-4-R-phenols (where R=CH3, C(CH3)3 or Br) with 2,2'-3,5-dinitrobenzoyliminodi(ethylamine) trihydrochloride in the presence of manganese(II) perchlorate. The IR spectra of complexes indicate the presence of uncoordinated perchlorate anions. The UV-Vis spectra of complexes suggest the distorted octahedral geometry around manganese(II) nuclei. The EPR spectra of Mn(II) complexes show a broad signal with g value 2.03-2.04, which is characteristic for octahedral high spin Mn(2+) complex. The observed room temperature magnetic moment values of the Mn(II) complexes (5.60-5.62B.M.) are less than the normal value (5.92B.M.), indicating weak antiferromagnetic coupling interaction between the two metal ions. Electrochemical studies of the complexes show two distinct quasi-reversible one electron transfer processes in the cathodic (E(1)pc=-0.73 to -0.76V, E(2)pc=-1.30 to -1.36V), and anodic (E(1)pa=1.02-1.11V, E(2)pa=1.32-1.79V) potential regions. Antibacterial efficacy of complexes have been screened against four Gram (-ve) and two Gram (+ve) bacterial strains. The DNA interaction studies suggest that these complexes bind with CT-DNA by intercalation, giving the binding affinity in the order 1>2>3. All the complexes display significant cleavage activity against circular plasmid pBR322 DNA. Docking simulation was performed to insert complexes into the crystal structure of EGFR tyrosine kinase and B-DNA at active site to determine the probable binding mode.

FDA-approved small-molecule kinase inhibitors

Kinases have emerged as one of the most intensively pursued targets in current pharmacological research, especially for cancer, due to their critical roles in cellular signaling. To date, the US FDA has approved 28 small-molecule kinase inhibitors, half of which were approved in the past 3 years. While the clinical data of these approved molecules are widely presented and structure-activity relationship (SAR) has been reported for individual molecules, an updated review that analyzes all approved molecules and summarizes current achievements and trends in the field has yet to be found. Here we present all approved small-molecule kinase inhibitors with an emphasis on binding mechanism and structural features, summarize current challenges, and discuss future directions in this field.

Regulation and function of tetrapyrrole biosynthesis in plants and algae

Tetrapyrroles are macrocyclic molecules with various structural variants and multiple functions in Prokaryotes and Eukaryotes. Present knowledge about the metabolism of tetrapyrroles reflects the complex evolution of the pathway in different kingdoms of organisms, the complexity of structural and enzymatic variations of enzymatic steps, as well as a wide range of regulatory mechanisms, which ensure adequate synthesis of tetrapyrrole end-products at any time of development and environmental condition. This review intends to highlight new findings of research on tetrapyrrole biosynthesis in plants and algae. In the course of the heme and chlorophyll synthesis in these photosynthetic organisms, glutamate, one of the central and abundant metabolites, is converted into highly photoreactive tetrapyrrole intermediates. Thereby, several mechanisms of posttranslational control are thought to be essential for a tight regulation of each enzymatic step. Finally, we wish to discuss the potential role of tetrapyrroles in retrograde signaling and point out perspectives of the formation of macromolecular protein complexes in tetrapyrrole biosynthesis as an efficient mechanism to ensure a fine-tuned metabolic flow in the pathway. This article is part of a Special Issue entitled: Chloroplast Biogenesis.

Functional Role of Histidine in the Conserved His-x-Asp Motif in the Catalytic Core of Protein Kinases

The His-x-Asp (HxD) motif is one of the most conserved structural components of the catalytic core of protein kinases; however, the functional role of the conserved histidine is unclear. Here we report that replacement of the HxD-histidine with Arginine or Phenylalanine in Aurora A abolishes both the catalytic activity and auto-phosphorylation, whereas the Histidine-to-tyrosine impairs the catalytic activity without affecting its auto-phosphorylation. Comparisons of the crystal structures of wild-type (WT) and mutant Aurora A demonstrate that the impairment of the kinase activity is accounted for by (1) disruption of the regulatory spine in the His-to-Arg mutant, and (2) change in the geometry of backbones of the Asp-Phe-Gly (DFG) motif and the DFG-1 residue in the His-to-Tyr mutant. In addition, bioinformatics analyses show that the HxD-histidine is a mutational hotspot in tumor tissues. Moreover, the H174R mutation of the HxD-histidine, in the tumor suppressor LKB1 abrogates the inhibition of anchorage-independent growth of A549 cells by WT LKB1. Based on these data, we propose that the HxD-histidine is involved in a conserved inflexible organization of the catalytic core that is required for the kinase activity. Mutation of the HxD-histidine may also be involved in the pathogenesis of some diseases including cancer.

Phosphoryl Transfer Reaction Snapshots in Crystals: INSIGHTS INTO THE MECHANISM OF PROTEIN KINASE A CATALYTIC SUBUNIT

To study the catalytic mechanism of phosphorylation catalyzed by cAMP-dependent protein kinase (PKA) a structure of the enzyme-substrate complex representing the Michaelis complex is of specific interest as it can shed light on the structure of the transition state. However, all previous crystal structures of the Michaelis complex mimics of the PKA catalytic subunit (PKAc) were obtained with either peptide inhibitors or ATP analogs. Here we utilized Ca(2+) ions and sulfur in place of the nucleophilic oxygen in a 20-residue pseudo-substrate peptide (CP20) and ATP to produce a close mimic of the Michaelis complex. In the ternary reactant complex, the thiol group of Cys-21 of the peptide is facing Asp-166 and the sulfur atom is positioned for an in-line phosphoryl transfer. Replacement of Ca(2+) cations with Mg(2+) ions resulted in a complex with trapped products of ATP hydrolysis: phosphate ion and ADP. The present structural results in combination with the previously reported structures of the transition state mimic and phosphorylated product complexes complete the snapshots of the phosphoryl transfer reaction by PKAc, providing us with the most thorough picture of the catalytic mechanism to date.

Determining the Functions of HIV-1 Tat and a Second Magnesium Ion in the CDK9/Cyclin T1 Complex: A Molecular Dynamics Simulation Study

The current paradigm of cyclin-dependent kinase (CDK) regulation based on the well-established CDK2 has been recently expanded. The determination of CDK9 crystal structures suggests the requirement of an additional regulatory protein, such as human immunodeficiency virus type 1 (HIV-1) Tat, to exert its physiological functions. In most kinases, the exact number and roles of the cofactor metal ions remain unappreciated, and the repertoire has thus gained increasing attention recently. Here, molecular dynamics (MD) simulations were implemented on CDK9 to explore the functional roles of HIV-1 Tat and the second Mg²⁺ ion at site 1 (Mg₁²⁺). The simulations unveiled that binding of HIV-1 Tat to CDK9 not only stabilized hydrogen bonds (H-bonds) between ATP and hinge residues Asp104 and Cys106, as well as between ATP and invariant Lys48, but also facilitated the salt bridge network pertaining to the phosphorylated Thr186 at the activation loop. By contrast, these H-bonds cannot be formed in CDK9 owing to the absence of HIV-1 Tat. MD simulations further revealed that the Mg₁²⁺ ion, coupled with the Mg₂²⁺ ion, anchored to the triphosphate moiety of ATP in its catalytic competent conformation. This observation indicates the requirement of the Mg₁²⁺ ion for CDK9 to realize its function. Overall, the introduction of HIV-1 Tat and Mg₁²⁺ ion resulted in the active site architectural characteristics of phosphorylated CDK9. These data highlighted the functional roles of HIV-1 Tat and Mg₁²⁺ ion in the regulation of CDK9 activity, which contributes an important complementary understanding of CDK molecular underpinnings.

Dynamics of protein kinases: insights from nuclear magnetic resonance

Protein kinases are ubiquitous enzymes with critical roles in cellular processes and pathology. As a result, researchers have studied their activity and regulatory mechanisms extensively. Thousands of X-ray structures give snapshots of the architectures of protein kinases in various states of activation and ligand binding. However, the extent of and manner by which protein motions and conformational dynamics underlie the function and regulation of these important enzymes is not well understood. Nuclear magnetic resonance (NMR) methods provide complementary information about protein conformation and dynamics in solution. However, until recently, the large size of these enzymes prevented researchers from using these methods with kinases. Developments in transverse relaxation-optimized spectroscopy (TROSY)-based techniques and more efficient isotope labeling strategies are now allowing researchers to carry out NMR studies on full-length protein kinases. In this Account, we describe recent insights into the role of dynamics in protein kinase regulation and catalysis that have been gained from NMR measurements of chemical shift changes and line broadening, residual dipolar couplings, and relaxation. These findings show strong associations between protein motion and events that control kinase activity. Dynamic and conformational changes occurring at ligand binding sites and other regulatory domains of these proteins propagate to conserved kinase core regions that mediate catalytic function. NMR measurements of slow time scale (microsecond to millisecond) motions also reveal that kinases carry out global exchange processes that synchronize multiple residues and allosteric interconversion between conformational states. Activating covalent modifications or ligand binding to form the Michaelis complex can induce these global processes. Inhibitors can also exploit the exchange properties of kinases by using conformational selection to form dynamically quenched states. These investigations have revealed that kinases are highly dynamic enzymes, whose regulation by interdomain interactions, ligand binding, and covalent modifications involve changes in motion and conformational equilibrium in a manner that can be correlated with function. Thus, NMR provides a unique window into the role of protein dynamics in kinase regulation and catalysis with important implications for drug design.

The effect of a widespread cancer-causing mutation on the inactive to active dynamics of the B-Raf kinase

Protein kinases play a key role in regulating cellular processes. Kinase dysfunction can lead to disease, making them an attractive target for drug design. The B-Raf kinase is a key target for the treatment of melanoma since a single mutation (V600E) is found in more than 50% of all malignant melanomas. Despite the importance of B-Raf in melanoma treatment, the molecular mechanism by which the mutation increases kinase activity remains elusive. Since kinases are tightly regulated by a conformational transition between an active and inactive state, which is difficult to capture experimentally, large-scale enhanced-sampling simulations are performed to examine the mechanism by which the V600E mutation enhances the activity of the B-Raf monomer. The results reveal that the mutation has a twofold effect. First, the mutation increases the barrier of the active to inactive transition trapping B-Raf in the active state. The mutation also increases the flexibility of the activation loop which might speed-up the rate-limiting step of phosphorylation. Both effects can be explained by the formation of salt-bridges with the Glu600 residue.

Dysfunctional conformational dynamics of protein kinase A induced by a lethal mutant of phospholamban hinder phosphorylation

The dynamic interplay between kinases and substrates is crucial for the formation of catalytically committed complexes that enable phosphoryl transfer. However, a clear understanding on how substrates modulate kinase structural dynamics to control catalytic efficiency is still missing. Here, we used solution NMR spectroscopy to study the conformational dynamics of two complexes of the catalytic subunit of the cAMP-dependent protein kinase A with WT and R14 deletion phospholamban, a lethal human mutant linked to familial dilated cardiomyopathy. Phospholamban is a central regulator of heart muscle contractility, and its phosphorylation by protein kinase A constitutes a primary response to β-adrenergic stimulation. We found that the single deletion of arginine in phospholamban's recognition sequence for the kinase reduces its binding affinity and dramatically reduces phosphorylation kinetics. Structurally, the mutant prevents the enzyme from adopting conformations and motions committed for catalysis, with concomitant reduction in catalytic efficiency. Overall, these results underscore the importance of a well-tuned structural and dynamic interplay between the kinase and its substrates to achieve physiological phosphorylation levels for proper Ca(2+) signaling and normal cardiac function.

Investigation of structural mimetics of natural phosphate ion binding motifs

Phosphates are ubiquitous in biology and nearly half of all proteins interact with their partners by means of recognition of phosphate residues. Therefore, a better understanding of the phosphate ion binding by peptidic structures is highly desirable. Two new receptors have been designed and synthesized and their anion binding properties in an acetonitrile solution have been determined. The structure of hosts mimics a part of the kinase active site that is responsible for the recognition of the phosphate residue. New hosts contain additional free amino groups with the aim to facilitate coordination of protonated anions, such as dihydrogen phosphate. According to spectrophotometric measurements, stepwise 1:1 and 1:2 binding modes have been observed for both receptors in the presence of acetate, hydrogen sulfate and dihydrogen phosphate. Compared with the acyclic receptor, the macrocyclic receptor has demonstrated a remarkably enhanced selectivity for dihydrogen phosphate over other anions. Fluorometric measurements have revealed different responses of the acyclic and macrocyclic receptors towards anions. However, in both cases, a 5–8 nm hypsochromic shift of fluorescence maximum has been observed upon interaction of acetate and dihydrogen phosphate with receptors.

Src protein-tyrosine kinase structure, mechanism, and small molecule inhibitors

The physiological Src proto-oncogene is a protein-tyrosine kinase that plays key roles in cell growth, division, migration, and survival signaling pathways. From the N- to C-terminus, Src contains a unique domain, an SH3 domain, an SH2 domain, a protein-tyrosine kinase domain, and a regulatory tail. The chief phosphorylation sites of human Src include an activating pTyr419 that results from phosphorylation in the kinase domain by an adjacent Src molecule and an inhibitory pTyr530 in the regulatory tail that results from phosphorylation by C-terminal Src kinase (Csk) or Chk (Csk homologous kinase). The oncogenic Rous sarcoma viral protein lacks the equivalent of Tyr530 and is constitutively activated. Inactive Src is stabilized by SH2 and SH3 domains on the rear of the kinase domain where they form an immobilizing and inhibitory clamp. Protein kinases including Src contain hydrophobic regulatory and catalytic spines and collateral shell residues that are required to assemble the active enzyme. In the inactive enzyme, the regulatory spine contains a kink or a discontinuity with a structure that is incompatible with catalysis. The conversion of inactive to active Src is accompanied by electrostatic exchanges involving the breaking and making of distinct sets of kinase domain salt bridges and hydrogen bonds. Src-catalyzed protein phosphorylation requires the participation of two Mg(2+) ions. Although nearly all protein kinases possess a common K/E/D/D signature, each enzyme exhibits its unique variations of the protein-kinase reaction template. Bosutinib, dasatinib, and ponatinib are Src/multikinase inhibitors that are approved by the FDA for the treatment of chronic myelogenous leukemia and vandetanib is approved for the treatment of medullary thyroid cancer. The Src and BCR-Abl inhibitors saracatinib and AZD0424, along with the previous four drugs, are in clinical trials for a variety of solid tumors including breast and lung cancers. Both ATP and targeted therapeutic Src protein kinase inhibitors such as dasatinib and ponatinib make hydrophobic contacts with catalytic spine residues and form hydrogen bonds with hinge residues connecting the small and large kinase lobes.

Open-close structural change upon ligand binding and two magnesium ions required for the catalysis of N-acetylhexosamine 1-kinase

Infant gut-associated bifidobacteria possess a metabolic pathway to utilize lacto-N-biose (Gal-β1,3-GlcNAc) and galacto-N-biose (Gal-β1,3-GalNAc) from human milk and glycoconjugates specifically. In this pathway, N-acetylhexosamine 1-kinase (NahK) catalyzes the phosphorylation of GlcNAc or GalNAc at the anomeric C1 position with ATP. Crystal structures of NahK have only been determined in the closed state. In this study, we determined open state structures of NahK in three different forms (apo, ADP complex, and ATP complex). A comparison of the open and closed state structures revealed an induced fit structural change defined by two rigid domains. ATP binds to the small N-terminal domain, and binding of the N-acetylhexosamine substrate to the large C-terminal domain induces a closing conformational change with a rotation angle of 16°. In the nucleotide binding site, two magnesium ions bridging the α-γ and β-γ phosphates were identified. A mutational analysis indicated that a residue coordinating both of the two magnesium ions (Asp228) is essential for catalysis. The involvement of two magnesium ions in the catalytic machinery is structurally similar to the catalytic structures of protein kinases and aminoglycoside phosphotransferases, but distinct from the structures of other anomeric kinases or sugar 6-kinases. These findings help to elucidate the possible evolutionary adaptation of substrate specificities and induced fit mechanism.

A dynamically coupled allosteric network underlies binding cooperativity in Src kinase

Protein tyrosine kinases are attractive drug targets because many human diseases are associated with the deregulation of kinase activity. However, how the catalytic kinase domain integrates different signals and switches from an active to an inactive conformation remains incompletely understood. Here we identify an allosteric network of dynamically coupled amino acids in Src kinase that connects regulatory sites to the ATP- and substrate-binding sites. Surprisingly, reactants (ATP and peptide substrates) bind with negative cooperativity to Src kinase while products (ADP and phosphopeptide) bind with positive cooperativity. We confirm the molecular details of the signal relay through the allosteric network by biochemical studies. Experiments on two additional protein tyrosine kinases indicate that the allosteric network may be largely conserved among these enzymes. Our work provides new insights into the regulation of protein tyrosine kinases and establishes a potential conduit by which resistance mutations to ATP-competitive kinase inhibitors can affect their activity.

Protein tyrosine kinases are subject to multiple regulatory mechanisms. Foda et al. show that reactants and products of the tyrosine kinase Src bind its catalytic domain with opposite cooperativity, and identify an allosteric network of dynamically coupled amino acids that underlie this behaviour.

Targeting conformational plasticity of protein kinases

The quest for ever more selective kinase inhibitors as potential future drugs has yielded a large repertoire of chemical probes that are selective for specific kinase conformations. These probes have been useful tools to obtain structural snapshots of kinase conformational plasticity. Similarly, kinetic and thermodynamic inhibitor binding experiments provide glimpses at the time scales and energetics of conformational interconversions. These experimental insights are complemented by computational predictions of conformational energy landscapes and simulations of conformational transitions and of the process of inhibitors binding to the protein kinase domain. A picture emerges in which highly selective inhibitors capitalize on the dynamic nature of kinases.

The energy landscape of adenylate kinase during catalysis

Kinases perform phosphoryl-transfer reactions in milliseconds; without enzymes, these reactions would take about 8,000 years under physiological conditions. Despite extensive studies, a comprehensive understanding of kinase energy landscapes, including both chemical and conformational steps, is lacking. Here we scrutinize the microscopic steps in the catalytic cycle of adenylate kinase, through a combination of NMR measurements during catalysis, pre-steady-state kinetics, molecular-dynamics simulations and crystallography of active complexes. We find that the Mg(2+) cofactor activates two distinct molecular events: phosphoryl transfer (>10(5)-fold) and lid opening (10(3)-fold). In contrast, mutation of an essential active site arginine decelerates phosphoryl transfer 10(3)-fold without substantially affecting lid opening. Our results highlight the importance of the entire energy landscape in catalysis and suggest that adenylate kinases have evolved to activate key processes simultaneously by precise placement of a single, charged and very abundant cofactor in a preorganized active site.

A QM/MM study of Kemptide phosphorylation catalyzed by protein kinase A. The role of Asp166 as a general acid/base catalyst

In this work a theoretical study of the γ-phosphoryl group transfer from ATP to Ser17 of the synthetic substrate Kemptide (LRRASLG) in protein kinase A (PKA) has been carried out with a solvated model of the PKA-Mg2ATP-Kemptide system based on the X-ray crystallographic structure. We have used high levels (B3LYP/MM and MP2/MM) of theory to determine the overall reaction paths of the so-called concerted loose mechanism trying to clarify some aspects of that mechanism still under debate. Our calculations demonstrate for the first time in a complete model of the ternary system the viability of the final step of the catalytic mechanism in which the protonation of the phosphokemptide product by Asp166 takes place. Asp166 is a base catalyst that abstracts the HγSer17 of Kemptide thus facilitating the phosphoryl transfer, but it also acts as an acid catalyst by donating the proton just accepted from Ser17 to the O2γATP atom of the phosphoryl group.

dbPPT: a comprehensive database of protein phosphorylation in plants

As one of the most important protein post-translational modifications, the reversible phosphorylation is critical for plants in regulating a variety of biological processes such as cellular metabolism, signal transduction and responses to environmental stress. Numerous efforts especially large-scale phosphoproteome profiling studies have been contributed to dissect the phosphorylation signaling in various plants, while a large number of phosphorylation events were identified. To provide an integrated data resource for further investigations, here we present a comprehensive database of dbPPT (database of Phosphorylation site in PlanTs, at http://dbppt.biocuckoo.org), which contains experimentally identified phosphorylation sites in proteins from plants. The phosphorylation sites in dbPPT were manually curated from the literatures, whereas datasets in other public databases were also integrated. In total, there were 82 175 phosphorylation sites in 31 012 proteins from 20 plant organisms in dbPPT, presenting a larger quantity of phosphorylation sites and a higher coverage of plant species in comparison with other databases. The proportions of residue types including serine, threonine and tyrosine were 77.99, 17.81 and 4.20%, respectively. All the phosphoproteins and phosphorylation sites in the database were critically annotated. Since the phosphorylation signaling in plants attracted great attention recently, such a comprehensive resource of plant protein phosphorylation can be useful for the research community.

Database URL:http://dbppt.biocuckoo.org

Synchronous opening and closing motions are essential for cAMP-dependent protein kinase A signaling

Conformational fluctuations play a central role in enzymatic catalysis. However, it is not clear how the rates and the coordination of the motions affect the different catalytic steps. Here, we used NMR spectroscopy to analyze the conformational fluctuations of the catalytic subunit of the cAMP-dependent protein kinase (PKA-C), a ubiquitous enzyme involved in a myriad of cell signaling events. We found that the wild-type enzyme undergoes synchronous motions involving several structural elements located in the small lobe of the kinase, which is responsible for nucleotide binding and release. In contrast, a mutation (Y204A) located far from the active site desynchronizes the opening and closing of the active cleft without changing the enzyme's structure, rendering it catalytically inefficient. Since the opening and closing motions govern the rate-determining product release, we conclude that optimal and coherent conformational fluctuations are necessary for efficient turnover of protein kinases.

Solution NMR of MPS-1 reveals a random coil cytosolic domain structure

Caenorhabditis elegans MPS1 is a single transmembrane helical auxiliary subunit that co-localizes with the voltage-gated potassium channel KVS1 in the nematode nervous system. MPS-1 shares high homology with KCNE (potassium voltage-gated channel subfamily E member) auxiliary subunits, and its cytosolic domain was reported to have a serine/threonine kinase activity that modulates KVS1 channel function via phosphorylation. In this study, NMR spectroscopy indicated that the full length and truncated MPS-1 cytosolic domain (134–256) in the presence or absence of n-dodecylphosphocholine detergent micelles adopted a highly flexible random coil secondary structure. In contrast, protein kinases usually adopt a stable folded conformation in order to implement substrate recognition and phosphoryl transfer. The highly flexible random coil secondary structure suggests that MPS-1 in the free state is unstructured but may require a substrate or binding partner to adopt stable structure required for serine/threonine kinase activity.

The SH2 domain regulates c-Abl kinase activation by a cyclin-like mechanism and remodulation of the hinge motion

Regulation of the c-Abl (ABL1) tyrosine kinase is important because of its role in cellular signaling, and its relevance in the leukemiogenic counterpart (BCR-ABL). Both auto-inhibition and full activation of c-Abl are regulated by the interaction of the catalytic domain with the Src Homology 2 (SH2) domain. The mechanism by which this interaction enhances catalysis is not known. We combined computational simulations with mutagenesis and functional analysis to find that the SH2 domain conveys both local and global effects on the dynamics of the catalytic domain. Locally, it regulates the flexibility of the αC helix in a fashion reminiscent of cyclins in cyclin-dependent kinases, reorienting catalytically important motifs. At a more global level, SH2 binding redirects the hinge motion of the N and C lobes and changes the conformational equilibrium of the activation loop. The complex network of subtle structural shifts that link the SH2 domain with the activation loop and the active site may be partially conserved with other SH2-domain containing kinases and therefore offer additional parameters for the design of conformation-specific inhibitors.

The Abl kinase is a key player in many crucial cellular processes. It is also an important anti-cancer drug target, because a mutation leading to the fusion protein Bcr-Abl is the main cause for chronic myeloid leukemia (CML). Abl inhibitors are currently the only pharmaceutical treatment for CML. There are two main difficulties associated with the development of kinase inhibitors: the high similarity between active sites of different kinases, which makes selectivity a challenge, and mutations leading to resistance, which make it mandatory to search for alternative drugs. One important factor controlling Abl is the interplay between the catalytic domain and an SH2 domain. We used computer simulations to understand how the interactions between the domains modify the dynamic of the kinase and detected both local and global effects. Based on our computer model, we suggested mutations that should alter the domain-domain interplay. Consequently, we tested the mutants experimentally and found that they support our hypothesis. We propose that our findings can be of help for the development of new classes of Abl inhibitors, which would modify the domain-domain interplay instead of interfering directly with the active site.

The structural basis for mTOR function

The phosphoinositide 3-kinase (PI3K) related protein kinases (PIKKs) are a family of protein kinases with a diverse range of vital cellular functions. Recent high-resolution crystal structures of the protein kinase mTOR suggest general architectural principles that are likely to be common to all of the PIKKs. Furthermore, the structures make clear the close relationship of the PIKKs to the PI3Ks. However, the structures also make clear the unique features of mTOR that enable its substrate specificity. The active site is deeply recessed and flanked by structural elements unique to the PIKKs, namely, the FRB domain, the LST8 binding element, and a C-terminal stretch of helices known as the FATC domain. The FRB has a conserved element in it that is part of a bipartite substrate recognition mechanism that is probably characteristic of all of the PIKKs. The FRB also binds the mTOR inhibitor rapamycin that has been referred to as an allosteric inhibitor, implying that this inhibitor is actually a competitive inhibitor of the protein substrate. This bipartite substrate-binding site also helps clarify how rapamycin can result in substrate-specific inhibition.

Study of the affinity between the protein kinase PKA and peptide substrates derived from kemptide using molecular dynamics simulations and MM/GBSA

We have carried out a protocol in computational biochemistry including molecular dynamics (MD) simulations and MM/GBSA free energy calculations on the complex between the protein kinase A (PKA) and the specific peptide substrate Kemptide (LRRASLG). We made the same calculations on other PKA complexes that contain Kemptide derivatives (with mutations of the arginines, and with deletions of N and C-terminal amino acids). We predicted shifts in the free energy changes from the free PKA to PKA-substrate complex (ΔΔG_E→ES) when Kemptide structure is modified (we consider that the calculated shifts correlate with the experimental shifts of the free energy changes from the free PKA to the transition states (ΔΔG_E→TS) determined by the catalytic efficiency (k_cat/K_M) changes). Our results demonstrate that it is possible to predict the kinetic properties of protein kinases using simple computational biochemistry methods. As an additional benefit, these methods give detailed molecular information that permit the analysis of the atomic forces that contribute to the affinity between protein kinases and their substrates.

N-terminus of the protein kinase CLK1 induces SR protein hyperphosphorylation

SR proteins are essential splicing factors that are regulated through multisite phosphorylation of their RS (arginine/serine-rich) domains by two major families of protein kinases. The SRPKs (SR-specific protein kinases) efficiently phosphorylate the arginine/serine dipeptides in the RS domain using a conserved docking groove in the kinase domain. In contrast, CLKs (Cdc2-like kinases) lack a docking groove and phosphorylate both arginine/serine and serine-proline dipeptides, modifications that generate a hyperphosphorylated state important for unique SR protein-dependent splicing activities. All CLKs contain long flexible N-terminal extensions (140-300 residues) that resemble the RS domains present in their substrate SR proteins. We showed that the N-terminus in CLK1 contacts both the kinase domain and the RS domain of the SR protein SRSF1 (SR protein splicing factor 1). This interaction not only is essential for facilitating hyperphosphorylation, but also induces co-operative binding of SRSF1 to RNA. The N-terminus of CLK1 enhances the total phosphoryl contents of a panel of physiological substrates including SRSF1, SRSF2, SRSF5 and Tra2β1 (transformer 2β1) by 2-3-fold. These findings suggest that CLK1-dependent hyperphosphorylation is the result of a general mechanism in which the N-terminus acts as a bridge connecting the kinase domain and the RS domain of the SR protein.

Phosphorylation- and nucleotide-binding-induced changes to the stability and hydrogen exchange patterns of JNK1β1 provide insight into its mechanisms of activation

Many studies have characterized how changes to the stability and internal motions of a protein during activation can contribute to their catalytic function, even when structural changes cannot be observed. Here, unfolding studies and hydrogen-deuterium exchange (HX) mass spectrometry were used to investigate the changes to the stability and conformation/conformational dynamics of JNK1β1 induced by phosphorylative activation. Equivalent studies were also employed to determine the effects of nucleotide binding on both inactive and active JNK1β1 using the ATP analogue, 5'-adenylyl-imidodiphosphate (AMP-PNP). JNK1β1 phosphorylation alters HX in regions involved in catalysis and substrate binding, changes that can be ascribed to functional modifications in either structure and/or backbone flexibility. Increased HX in the hinge between the N- and C-terminal domains implied that it acquires enhanced flexibility upon phosphorylation that may be a prerequisite for interdomain closure. In combination with the finding that nucleotide binding destabilizes the kinase, the patterns of solvent protection by AMP-PNP were consistent with a novel mode of nucleotide binding to the C-terminal domain of a destabilized and open domain conformation of inactive JNK1β1. Solvent protection by AMP-PNP of both N- and C-terminal domains in active JNK1β1 revealed that the domains close around nucleotide upon phosphorylation, concomitantly stabilizing the kinase. This suggests that phosphorylation activates JNK1β1 in part by increasing hinge flexibility to facilitate interdomain closure and the creation of a functional active site. By uncovering the complex interplay that occurs between nucleotide binding and phosphorylation, we present new insight into the unique mechanisms by which JNK1β1 is regulated.

A QM/MM study of the associative mechanism for the phosphorylation reaction catalyzed by protein kinase A and its D166A mutant

Here we analyze in detail the possible catalytic role of the associative mechanism in the γ-phosphoryl transfer reaction in the catalytic subunit of the mammalian cyclic AMP-dependent protein kinase (PKA) enzyme and its D166A mutant. We have used a complete solvated model of the ATP-Mg2-Kemptide/PKA system and good levels of theory (B3LYP/MM and MP2/MM) to determine several potential energy paths from different MD snapshots, and we present a deep analysis of the interaction distances and energies between ligands, metals and enzyme residues. We have also tested the electrostatic stabilization of the transition state structures localized herein with the charge balance hypothesis. Overall, the results obtained in this work reopen the discussion about the plausibility of the associative reaction pathway and highlight the proposed role of the catalytic triad Asp166-Lys168-Thr201.

Molecular features of product release for the PKA catalytic cycle

Although ADP release is the rate limiting step in product turnover by protein kinase A, the steps and motions involved in this process are not well resolved. Here we report the apo and ADP bound structures of the myristylated catalytic subunit of PKA at 2.9 and 3.5 Å resolution, respectively. The ADP bound structure adopts a conformation that does not conform to the previously characterized open, closed, or intermediate states. In the ADP bound structure, the C-terminal tail and Gly-rich loop are more closed than in the open state adopted in the apo structure but are also much more open than the intermediate or closed conformations. Furthermore, ADP binds at the active site with only one magnesium ion, termed Mg2 from previous structures. These structures thus support a model where ADP release proceeds through release of the substrate and Mg1 followed by lifting of the Gly-rich loop and disengagement of the C-terminal tail. Coupling of these two structural elements with the release of the first metal ion fills in a key step in the catalytic cycle that has been missing and supports an ensemble of correlated conformational states that mediate the full catalytic cycle for a protein kinase.

Detecting enzyme activities with exogenous MRI contrast agents

This review focuses on exogenous magnetic resonance imaging (MRI) contrast agents that are responsive to enzyme activity. Enzymes can catalyze a change in water access, rotational tumbling time, the proximity of a (19)F-labeled ligand, the aggregation state, the proton chemical-exchange rate between the agent and water, or the chemical shift of (19)F, (31)P, (13)C or a labile (1)H of an agent, all of which can be used to detect enzyme activity. The variety of agents attests to the creativity in developing enzyme-responsive MRI contrast agents.

Structural dynamic analysis of apo and ATP-bound IRAK4 kinase

Interleukin-1 receptor-associated kinases (IRAKs) are Ser/Thr protein kinases that play an important role as signaling mediators in the signal transduction facilitated by the Toll-like receptor (TLR) and interleukin-1 receptor families. Among IRAK family members, IRAK4 is one of the drug targets for diseases related to the TLR and IL-1R signaling pathways. Experimental evidence suggests that the IRAK4 kinase domain is phosphorylated in its activation loop at T342, T345, and S346 in the fully activated state. However, the molecular interactions of subdomains within the active and inactive IRAK4 kinase domain are poorly understood. Hence, we employed a long-range molecular dynamics (MD) simulation to compare apo IRAK4 kinase domains (phosphorylated and unphosphorylated) and ATP-bound phosphorylated IRAK4 kinase domains. The MD results strongly suggested that lobe uncoupling occurs in apo unphosphorylated IRAK4 kinase via the disruption of the R334/T345 and R310/T345 interaction. In addition, apo unphosphorylated trajectory result in high mobility, particularly in the N lobe, activation segment, helix αG, and its adjoining loops. The Asp-Phe-Gly (DFG) and His-Arg-Asp (HRD) conserved kinase motif analysis showed the importance of these motifs in IRAK4 kinase activation. This study provides important information on the structural dynamics of IRAK4 kinase, which will aid in inhibitor development.

Exploration of type II binding mode: A privileged approach for kinase inhibitor focused drug discovery?

The ATP site of kinases displays remarkable conformational flexibility when accommodating chemically diverse small molecule inhibitors. The so-called activation segment, whose conformation controls catalytic activity and access to the substrate binding pocket, can undergo a large conformational change with the active state assuming a ‘DFG-in’ and an inactive state assuming a ‘DFG-out’ conformation. Compounds that preferentially bind to the DFG-out conformation are typically called ‘type II’ inhibitors in contrast to ‘type I’ inhibitors that bind to the DFG-in conformation. This review surveys the large number of type II inhibitors that have been developed and provides an analysis of their crystallographically determined binding modes. Using a small library of type II inhibitors, we demonstrate that more than 200 kinases can be targeted, suggesting that type II inhibitors may not be intrinsically more selective than type I inhibitors.

ErbB/HER protein-tyrosine kinases: Structures and small molecule inhibitors

The epidermal growth factor receptor (EGFR) family consists of four members that belong to the ErbB lineage of proteins (ErbB1-4). These receptors consist of an extracellular domain, a single hydrophobic transmembrane segment, and an intracellular portion with a juxtamembrane segment, a protein kinase domain, and a carboxyterminal tail. The ErbB proteins function as homo and heterodimers. Growth factor binding to EGFR induces a large conformational change in the extracellular domain. Two ligand-EGFR complexes unite to form a back-to-back dimer in which the ligands are on opposite sides of the aggregate. Following ligand binding, EGFR intracellular kinase domains form an asymmetric dimer. The carboxyterminal lobe of the activator kinase of the dimer interacts with the amino-terminal lobe of the receiver kinase thereby leading to its allosteric stimulation. Several malignancies are associated with the mutation or increased expression of members of the ErbB family including lung, breast, stomach, colorectal, head and neck, and pancreatic carcinomas. Gefitinib, erlotinib, and afatinib are orally effective protein-kinase targeted quinazoline derivatives that are used in the treatment of ERBB1-mutant lung cancer and lapatinib is an orally effective quinazoline derivative used in the treatment of ErbB2-overexpressing breast cancer. Moreover, monoclonal antibodies that target the extracellular domain of ErbB2 are used for the treatment of ErbB2-positive breast cancer and monoclonal antibodies that target ErbB1 and are used for the treatment of colorectal cancer. Cancers treated with these targeted drugs eventually become resistant to them, and a current goal of research is to develop drugs that are effective against drug-resistant tumors.

Day of the dead: pseudokinases and pseudophosphatases in physiology and disease

Pseudophosphatases and pseudokinases are increasingly viewed as integral elements of signaling pathways, and there is mounting evidence that they have frequently retained the ability to interact with cellular 'substrates', and can exert important roles in different diseases. However, these pseudoenzymes have traditionally received scant attention compared to classical kinases and phosphatases. In this review we explore new findings in the emerging pseudokinase and pseudophosphatase fields, and discuss their different modes of action which include exciting new roles as scaffolds, anchors, spatial modulators, traps, and ligand-driven regulators of canonical kinases and phosphatases. Thus, it is now apparent that pseudokinases and pseudophosphatases both support and drive a panoply of signaling networks. Finally, we highlight recent evidence on their involvement in human pathologies, marking them as potential novel drug targets.

Metal-free cAMP-dependent protein kinase can catalyze phosphoryl transfer

X-ray structures of several ternary product complexes of the catalytic subunit of cAMP-dependent protein kinase (PKAc) have been determined with no bound metal ions and with Na(+) or K(+) coordinated at two metal-binding sites. The metal-free PKAc and the enzyme with alkali metals were able to facilitate the phosphoryl transfer reaction. In all studied complexes, the ATP and the substrate peptide (SP20) were modified into the products ADP and the phosphorylated peptide. The products of the phosphotransfer reaction were also found when ATP-γS, a nonhydrolyzable ATP analogue, reacted with SP20 in the PKAc active site containing no metals. Single turnover enzyme kinetics measurements utilizing (32)P-labeled ATP confirmed the phosphotransferase activity of the enzyme in the absence of metal ions and in the presence of alkali metals. In addition, the structure of the apo-PKAc binary complex with SP20 suggests that the sequence of binding events may become ordered in a metal-free environment, with SP20 binding first to prime the enzyme for subsequent ATP binding. Comparison of these structures reveals conformational and hydrogen bonding changes that might be important for the mechanism of catalysis.

Computational delineation of tyrosyl-substrate recognition and catalytic landscapes by the epidermal growth factor receptor tyrosine kinase domain

The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase (RTK), which catalyzes protein phosphorylation reactions by transferring the γ-phosphoryl group from an ATP molecule to the hydroxyl group of tyrosine residues in protein substrates. EGFR is an important drug target in the treatment of cancers and a better understanding of the receptor function is critical to discern cancer mechanisms. We employ a suite of molecular simulation methods to explore the mechanism of substrate recognition and to delineate the catalytic landscape of the phosphoryl transfer reaction. Based on our results, we propose that a highly conserved region corresponding to Val852-Pro853-Ile854-Lys855-Trp856 in the EGFR tyrosine kinase domain (TKD) is essential for substrate binding. We also provide a possible explanation for the established experimental observation that protein tyrosine kinases (including EGFR) select substrates with a glutamic acid at the P - 1 position and a large hydrophobic amino acid at the P + 1 position. Furthermore, our mixed quantum mechanics/molecular mechanics (QM/MM) simulations show that the EGFR protein kinase favors the dissociative mechanism, although an alternative channel through the formation of an associative transition state is also possible. Our simulations establish some key molecular rules in the operation for substrate-recognition and for phosphoryl transfer in the EGFR TKD.

The structure and substrate specificity of human Cdk12/Cyclin K

Phosphorylation of the RNA polymerase II C-terminal domain (CTD) by cyclin-dependent kinases is important for productive transcription. Here we determine the crystal structure of Cdk12/CycK and analyse its requirements for substrate recognition. Active Cdk12/CycK is arranged in an open conformation similar to that of Cdk9/CycT but different from those of cell cycle kinases. Cdk12 contains a C-terminal extension that folds onto the N- and C-terminal lobes thereby contacting the ATP ribose. The interaction is mediated by an HE motif followed by a polybasic cluster that is conserved in transcriptional CDKs. Cdk12/CycK showed the highest activity on a CTD substrate prephosphorylated at position Ser7, whereas the common Lys7 substitution was not recognized. Flavopiridol is most potent towards Cdk12 but was still 10-fold more potent towards Cdk9. T-loop phosphorylation of Cdk12 required coexpression with a Cdk-activating kinase. These results suggest the regulation of Pol II elongation by a relay of transcriptionally active CTD kinases.

Cyclin-dependent kinase 12 (Cdk12) phosphorylates the C-terminal domain (CTD) of RNA polymerase II to regulate transcription. Here, the authors solve the crystal structure of the Cdk12 kinase domain and show that Cdk12 has its highest activity on a CTD substrate that carries a serine 7 phosphorylation.

Activation pathway of Src kinase reveals intermediate states as targets for drug design

Unregulated activation of Src kinases leads to aberrant signaling, uncontrolled growth, and differentiation of cancerous cells. Reaching a complete mechanistic understanding of large scale conformational transformations underlying the activation of kinases could greatly help in the development of therapeutic drugs for the treatment of these pathologies. In principle, the nature of conformational transition could be modeled in silico via atomistic molecular dynamics simulations, although this is very challenging due to the long activation timescales. Here, we employ a computational paradigm that couples transition pathway techniques and Markov state model-based massively distributed simulations for mapping the conformational landscape of c-src tyrosine kinase. The computations provide the thermodynamics and kinetics of kinase activation for the first time, and help identify key structural intermediates. Furthermore, the presence of a novel allosteric site in an intermediate state of c-src that could be potentially utilized for drug design is predicted.

Theoretical studies of cyclic adenosine monophosphate dependent protein kinase: native enzyme and ground-state and transition-state analogues

The mechanisms of phosphoryl transfer enzymes have garnered considerable attention. Cyclic AMP-dependent protein kinase (cAPK) catalyzes the transfer of the γ phosphoryl group of ATP to the serine hydroxyl group of a peptide chain. Metal-containing fluoro species have been used as transition-state and ground-state analogues in a variety of phosphoryl transfer enzymes and have shed light on the nature of the requirements in the active site to catalyze phosphoryl transfer. For cAPK, we present computational studies of the mechanism of phosphoryl transfer and the structure and (19)F NMR spectra of various ground- (BeF3(-)) and transition-state (MgF3(-), AlF4(-), and AlF3(0)) analogues. With native substrate, the phosphoryl transfer proceeds through a five-coordinate phosphorane transition state, i.e., there is not a five-coordinate phosphorane intermediate. Comparisons of simulated and experimental (19)F NMR spectra show cAPK prefers a monoanionic analogue (MgF3(-) or AlF4(-)) over a neutral analogue (AlF3), supporting the charge balance hypothesis.

Phosphoproteomic analysis of chromoplasts from sweet orange during fruit ripening

Like other types of plastids, chromoplasts have essential biosynthetic and metabolic activities which may be regulated via post-translational modifications, such as phosphorylation, of their resident proteins. We here report a proteome-wide mapping of in vivo phosphorylation sites in chromoplast-enriched samples prepared from sweet orange [Citrus sinensis (L.) Osbeck] at different ripening stages by titanium dioxide-based affinity chromatography for phosphoprotein enrichment with LC-MS/MS. A total of 109 plastid-localized phosphoprotein candidates were identified that correspond to 179 unique phosphorylation sites in 135 phosphopeptides. On the basis of Motif-X analysis, two distinct types of phosphorylation sites, one as proline-directed phosphorylation motif and the other as casein kinase II motif, can be generalized from these identified phosphopeptides. While most identified phosphoproteins show high homology to those already identified in plastids, approximately 22% of them are novel based on BLAST search using the public databases PhosPhAt and P(3) DB. A close comparative analysis showed that approximately 50% of the phosphoproteins identified in citrus chromoplasts find obvious counterparts in the chloroplast phosphoproteome, suggesting a rather high-level of conservation in basic metabolic activities in these two types of plastids. Not surprisingly, the phosphoproteome of citrus chromoplasts is also characterized by the lack of phosphoproteins involved in photosynthesis and by the presence of more phosphoproteins implicated in stress/redox responses. This study presents the first comprehensive phosphoproteomic analysis of chromoplasts and may help to understand how phosphorylation regulates differentiation of citrus chromoplasts during fruit ripening.