How mitogen-activated protein kinases recognize and phosphorylate their targets: A QM/MM study

Activation Loop Plasticity and Active Site Coupling in the MAP Kinase, ERK2

Previous studies of the protein kinase, ERK2, using NMR and hydrogen-exchange measurements have shown changes in dynamics accompanying its activation by phosphorylation. However, knowledge about the conformational motions involved is incomplete. Here, we examined ERK2 using long conventional molecular dynamics (MD) simulations starting from crystal structures of phosphorylated (2P) and unphosphorylated (0P) forms. Individual trajectories were run for (5 to 25) μs, totaling 727 μs. The results show unexpected flexibility of the A-loop, with multiple long-lived (>5 μs) conformational states in both 2P- and 0P-ERK2. Differential contact network and principal component analyses reveal coupling between the A-loop fold and active site dynamics, with evidence for conformational selection in the kinase core of 2P-ERK2 but not 0P-ERK2. Simulations of 2P-ERK2 show A-loop states corresponding to restrained dynamics within the N-lobe, including regions around catalytic residues. One A-loop conformer forms lasting interactions with the L16 segment, leading to reduced RMSF and greater compaction in the active site. By contrast, simulations of 0P-ERK2 reveal excursions of A-loop residues away from the C-lobe, leading to greater active site mobility. Thus, the A-loop in ERK2 switches between distinct conformations that reflect coupling with the active site, possibly via the L16 segment. Crystal packing interactions suggest that lattice contacts with the A-loop may restrain its structural variation in X-ray structures of ERK2. The novel conformational states identified by MD expand our understanding of ERK2 regulation, by linking the activated state of the kinase to reduced dynamics and greater compaction surrounding the catalytic site.

Activation loop plasticity and active site coupling in the MAP kinase, ERK2

Changes in the dynamics of the protein kinase, ERK2, have been shown to accompany its activation by dual phosphorylation. However, our knowledge about the conformational changes represented by these motions is incomplete. Previous NMR relaxation dispersion studies showed that active, dual-phosphorylated ERK2 undergoes global exchange between at least two energetically similar conformations. These findings, combined with measurements by hydrogen exchange mass spectrometry (HX-MS), suggested that the global conformational exchange involves motions of the activation loop (A-loop) that are coupled to regions surrounding the kinase active site. In order to better understand the contribution of dynamics to the activation of ERK2, we applied long conventional molecular dynamics (MD) simulations starting from crystal structures of active, phosphorylated (2P), and inactive, unphosphorylated (0P) ERK2. Individual trajectories were run for (5 to 25) µ s and totaled 727 µ s. The results showed that the A-loop is unexpectedly flexible in both 2P- and 0P-ERK2, and able to adopt multiple long-lived (>5 µ s) conformational states. Simulations starting from the X-ray structure of 2P-ERK2 (2ERK) revealed A-loop states corresponding to restrained dynamics within the N-lobe, including regions surrounding catalytic residues. One A-loop conformer forms lasting interactions with the C-terminal L16 segment and shows reduced RMSF and greater compaction in the active site. By contrast, simulations starting from the most common X-ray conformation of 0P-ERK2 (5UMO) reveal frequent excursions of A-loop residues away from a C-lobe docking site pocket and towards a new state that shows greater dynamics in the N-lobe and disorganization around the active site. Thus, the A-loop in ERK2 appears to switch between distinct conformational states that reflect allosteric coupling with the active site, likely occurring via the L16 segment. Analyses of crystal packing interactions across many structural datasets suggest that the A-loop observed in X-ray structures of ERK2 may be driven by lattice contacts and less representative of the solution structure. The novel conformational states identified by MD expand our understanding of ERK2 regulation, by linking the activated state of the kinase to reduced dynamics and greater compaction surrounding the catalytic site.

Specificity and Reactivity of Mycobacterium tuberculosis Serine/Threonine Kinases PknG and PknB

Mycobacterium tuberculosis (Mtb), the causative agent of Tuberculosis, has 11 eukaryotic-like serine/threonine protein kinases, which play essential roles in cell growth, signal transduction, and pathogenesis. Protein kinase G (PknG) regulates the carbon and nitrogen metabolism by phosphorylation of the glycogen accumulation regulator (GarA) protein at Thr21. Protein kinase B (PknB) is involved in cell wall synthesis and cell shape, as well as phosphorylates GarA but at Thr22. While PknG seems to be constitutively activated and recognition of GarA requires phosphorylation in its unstructured tail, PknB activation is triggered by phosphorylation of its activation loop, which allows binding of the forkhead-associated domain of GarA. In the present work, we used molecular dynamics and quantum-mechanics/molecular mechanics simulations of the catalytically competent complex and kinase activity assays to understand PknG/PknB specificity and reactivity toward GarA. Two hydrophobic residues in GarA, Val24 and Phe25, seem essential for PknG binding and allow specificity for Thr21 phosphorylation. On the other hand, phosphorylated residues in PknB bind Arg26 in GarA and regulate its specificity for Thr22. We also provide a detailed analysis of the free energy profile for the phospho-transfer reaction and show why PknG has a constitutively active conformation not requiring priming phosphorylation in contrast to PknB. Our results provide new insights into these two key enzymes relevant for Mtb and the mechanisms of serine/threonine phosphorylation in bacteria.

MAP Kinase-Mediated Activation of RSK1 and MK2 Substrate Kinases

Mitogen-activated protein kinases (MAPKs) control essential eukaryotic signaling pathways. While much has been learned about MAPK activation, much less is known about substrate recruitment and specificity. MAPK substrates may be other kinases that are crucial to promote a further diversification of the signaling outcomes. Here, we used a variety of molecular and cellular tools to investigate the recruitment of two substrate kinases, RSK1 and MK2, to three MAPKs (ERK2, p38α, and ERK5). Unexpectedly, we identified that kinase heterodimers form structurally and functionally distinct complexes depending on the activation state of the MAPK. These may be incompatible with downstream signaling, but naturally they may also form structures that are compatible with the phosphorylation of the downstream kinase at the activation loop, or alternatively at other allosteric sites. Furthermore, we show that small-molecule inhibitors may affect the quaternary arrangement of kinase heterodimers and thus influence downstream signaling in a specific manner.

Conformational and Reaction Dynamic Coupling in Histidine Kinases: Insights from Hybrid QM/MM Simulations

Histidine kinases (HK) of bacterial two-component systems represent a hallmark of allosterism in proteins, being able to detect a signal through the sensor domain and transmit this information through the protein matrix to the kinase domain which, once active, autophosphorylates a specific histidine residue. Inactive-to-active transition results in a large conformational change that moves the kinase on top of the histidine. In the present work, we use several molecular simulation techniques (Molecular Dynamics, Hybrid QM/MM, and constant pH molecular dynamics) to study the activation and autophosphorylation reactions in L. plantarum WalK, a cis-acting HK. In agreement with previous results, we show that the chemical step requires tight coupling with the conformational step in order to maintain the histidine phosphoacceptor in the correct tautomeric state, with a reactive δ-nitrogen. During the conformational transition, the kinase domain is never released and walks along the HK helix axis, breaking and forming several conserved residue-based contacts. The phosphate transfer reaction is concerted in the transition state region and is catalyzed through the stabilization of the negative developing charge of transferring phosphate along the reaction.

Kinase Activation by Small Conformational Changes

Protein kinases (PKs) are allosteric enzymes that play an essential role in signal transduction by regulating a variety of key cellular processes. Most PKs suffer conformational rearrangements upon phosphorylation that strongly enhance the catalytic activity. Generally, it involves the movement of the phosphorylated loop toward the active site and the rotation of the whole C-terminal lobe. However, not all kinases undergo such a large configurational change: The MAPK extracellular signal-regulated protein kinases ERK1 and ERK2 achieve a 50 000 fold increase in kinase activity with only a small motion of the C-terminal region. In the present work, we used a combination of molecular simulation tools to characterize the conformational landscape of ERK2 in the active (phosphorylated) and inactive (unphosphorylated) states in solution in agreement with NMR experiments. We show that the chemical reaction barrier is strongly dependent on ATP conformation and that the "active" low-barrier configuration is subtly regulated by phosphorylation, which stabilizes a key salt bridge between the conserved Lys52 and Glu69 belonging to helix-C and promotes binding of a second Mg ion. Our study highlights that the on-off switch embedded in the kinase fold can be regulated by small, medium, and large conformational changes.

Arrestin-3 scaffolding of the JNK3 cascade suggests a mechanism for signal amplification

Scaffold proteins tether and orient components of a signaling cascade to facilitate signaling. Although much is known about how scaffolds colocalize signaling proteins, it is unclear whether scaffolds promote signal amplification. Here, we used arrestin-3, a scaffold of the ASK1-MKK4/7-JNK3 cascade, as a model to understand signal amplification by a scaffold protein. We found that arrestin-3 exhibited >15-fold higher affinity for inactive JNK3 than for active JNK3, and this change involved a shift in the binding site following JNK3 activation. We used systems biochemistry modeling and Bayesian inference to evaluate how the activation of upstream kinases contributed to JNK3 phosphorylation. Our combined experimental and computational approach suggested that the catalytic phosphorylation rate of JNK3 at Thr-221 by MKK7 is two orders of magnitude faster than the corresponding phosphorylation of Tyr-223 by MKK4 with or without arrestin-3. Finally, we showed that the release of activated JNK3 was critical for signal amplification. Collectively, our data suggest a "conveyor belt" mechanism for signal amplification by scaffold proteins. This mechanism informs on a long-standing mystery for how few upstream kinase molecules activate numerous downstream kinases to amplify signaling.

Nature of the Pre-Chemistry Ensemble in Mitogen-Activated Protein Kinases

In spite of the availability of a significant amount of structural detail on docking interactions involving mitogen-activated protein kinases (MAPKs) and their substrates, the mechanism by which the disordered phospho-acceptor on the substrate transiently interacts with the kinase catalytic elements and is phosphorylated, often with high efficiency, remains poorly understood. Here, this dynamic interaction is analyzed in the context of available biophysical and biochemical data for ERK2, an archetypal MAPK. A hypothesis about the nature of the ternary complex involving a MAPK, its substrate, and ATP immediately prior to the chemical step (the pre-chemistry complex) is proposed. It is postulated that the solution ensemble (the pre-chemistry ensemble) representing the pre-chemistry complex comprises several conformations that are linked by dynamics on multiple timescales. These individual conformations possess different intrinsic abilities to proceed through the chemical step. The overall rate of chemistry is therefore related to the microscopic nature of the pre-chemistry ensemble, its constituent conformational microstates, and their intrinsic abilities to yield a phosphorylated product. While characterizing these microstates within the pre-chemistry ensemble in atomic or near-atomic detail is an extremely challenging proposition, recent developments in hybrid methodologies that employ computational approaches driven by experimental data appear to provide the most promising path forward toward achieving this goal.

Local destabilization, rigid body, and fuzzy docking facilitate the phosphorylation of the transcription factor Ets-1 by the mitogen-activated protein kinase ERK2

Mitogen-activated protein (MAP) kinase substrates are believed to require consensus docking motifs (D-site, F-site) to engage and facilitate efficient site-specific phosphorylation at specific serine/threonine-proline sequences by their cognate kinases. In contrast to other MAP kinase substrates, the transcription factor Ets-1 has no canonical docking motifs, yet it is efficiently phosphorylated by the MAP kinase ERK2 at a consensus threonine site (T38). Using NMR methodology, we demonstrate that this phosphorylation is enabled by a unique bipartite mode of ERK2 engagement by Ets-1 and involves two suboptimal noncanonical docking interactions instead of a single canonical docking motif. The N terminus of Ets-1 interacts with a part of the ERK2 D-recruitment site that normally accommodates the hydrophobic sidechains of a canonical D-site, retaining a significant degree of disorder in its ERK2-bound state. In contrast, the C-terminal region of Ets-1, including its Pointed (PNT) domain, engages in a largely rigid body interaction with a section of the ERK2 F-recruitment site through a binding mode that deviates significantly from that of a canonical F-site. This latter interaction is notable for the destabilization of a flexible helix that bridges the phospho-acceptor site to the rigid PNT domain. These two spatially distinct, individually weak docking interactions facilitate the highly specific recognition of ERK2 by Ets-1, and enable the optimal localization of its dynamic phospho-acceptor T38 at the kinase active site to enable efficient phosphorylation.

Functional roles of magnesium binding to extracellular signal-regulated kinase 2 explored by molecular dynamics simulations and principal component analysis

Molecular dynamics (MD) simulations coupled with principal component (PC) analysis were carried out to study functional roles of Mg²⁺ binding to extracellular signal-regulated kinase 2 (ERK2). The results suggest that Mg²⁺ binding heavily decreases eigenvalue of the first principal component and totally inhibits motion strength of ERK2, which favors stabilization of ERK2 structure. Binding free energy predictions indicate that Mg²⁺ binding produces an important effect on binding ability of adenosine triphosphate (ATP) to ERK2 and strengthens the ATP binding. The calculations of residue-based free energy decomposition show that lack of Mg²⁺ weakens interactions between the hydrophobic rings of ATP and five residues I29, V37, A50, L105, and L154. Hydrogen bond analyses also prove that Mg²⁺ binding increases occupancies of hydrogen bonds formed between ATP and residues K52, Q103, D104, and M106. We expect that this study can provide a significant theoretical hint for designs of anticancer drugs targeting ERK2.

Tumor Necrosis Factor Alpha Induces Neural Stem Cell Apoptosis Through Activating p38 MAPK Pathway

Tumor necrosis factor alpha (TNF-α) is an essential cytokine that mediates cell death and has been shown to play a potential role in inducing neural stem cell (NSC) apoptosis. We have previously shown that TNF-α antagonist etanercept can suppress the transplanted NSC apoptosis induced by TNF-α in spinal cord injury (SCI) sites; however, the precise molecular mechanism remains unclear. This study aimed to investigate the signaling pathways responsible for TNF-α-induced apoptosis in NSCs. TNF-α treatment impairs cell viability and increases apoptosis of NSCs in concentration- and time-dependent manners. This is embodied in an increase in Bax and cleaved caspase-3 production, coupled with decreased Bcl-2 levels. Additionally, TNF-α remarkably increased the expression of phosphatidylinositol p38 Mitogen-activated protein kinase (p38 MAPK) in NSCs. p38 MAPK regulates apoptosis, acting as an apoptotic signal due to TNF-α exposure. TNF-α-induced apoptosis was significantly alleviated by the p38 MAPK pathway inhibitor SB203580, as well as targeted inhibition of p38 gene in NSCs, or TNF-α antagonist etanercept. These results suggest that TNF-α induces NSCs apoptosis by activating the p38 MAPK signaling pathway and etanercept acts as an effective TNF-α antagonist to prevent p38 MAPK-dependent apoptosis induced by TNF-α in NSCs. Our research represents a potential gene targeting that can prevent unnecessary grafted cell death after transplantation into the SCI models.

Phosphoryl transfer reaction catalyzed by membrane diacylglycerol kinase: a theoretical mechanism study

Diacylglycerol kinase is an integral membrane protein which catalyzes phosphoryl transfer from ATP to diacylglycerol. As the smallest kinase known, it shares no sequence homology with conventional kinases and possesses a distinct trimer structure. Thus far, its catalytic mechanism remains elusive. Using molecular dynamics and quantum mechanics calculations, we investigated the co-factor and the substrate binding and phosphoryl transfer mechanism. Based on the analysis of density functional theory calculations, we reveal that the phosphorylation reaction of diacylglycerol kinase features the same phosphoryl transfer mechanism as other kinases, despite its unique structural properties. Our results further show that the active site is relatively open and able to accommodate ligands in multiple orientations, suggesting that the optimization of binding orientations and conformational changes would occur prior to actual phosphoryl transfer.

Improving Efficiency in SMD Simulations Through a Hybrid Differential Relaxation Algorithm

The fundamental object for studying a (bio)chemical reaction obtained from simulations is the free energy profile, which can be directly related to experimentally determined properties. Although quite accurate hybrid quantum (DFT based)-classical methods are available, achieving statistically accurate and well converged results at a moderate computational cost is still an open challenge. Here, we present and thoroughly test a hybrid differential relaxation algorithm (HyDRA), which allows faster equilibration of the classical environment during the nonequilibrium steering of a (bio)chemical reaction. We show and discuss why (in the context of Jarzynski's Relationship) this method allows obtaining accurate free energy profiles with smaller number of independent trajectories and/or faster pulling speeds, thus reducing the overall computational cost. Moreover, due to the availability and straightforward implementation of the method, we expect that it will foster theoretical studies of key enzymatic processes.

Structural and Dynamic Features of F-recruitment Site Driven Substrate Phosphorylation by ERK2

The F-recruitment site (FRS) of active ERK2 binds F-site (Phe-x-Phe-Pro) sequences found downstream of the Ser/Thr phospho-acceptor on cellular substrates. Here we apply NMR methods to analyze the interaction between active ERK2 (ppERK2), and a 13-residue F-site-bearing peptide substrate derived from its cellular target, the transcription factor Elk-1. Our results provide detailed insight into previously elusive structural and dynamic features of FRS/F-site interactions and FRS-driven substrate phosphorylation. We show that substrate F-site engagement significantly quenches slow dynamics involving the ppERK2 activation-loop and the FRS. We also demonstrate that the F-site phenylalanines make critical contacts with ppERK2, in contrast to the proline whose cis-trans isomerization has no significant effect on F-site recognition by the kinase FRS. Our results support a mechanism where phosphorylation of the disordered N-terminal phospho-acceptor is facilitated by its increased productive encounters with the ppERK2 active site due to docking of the proximal F-site at the kinase FRS.

Protein topology determines cysteine oxidation fate: the case of sulfenyl amide formation among protein families

Cysteine residues have a rich chemistry and play a critical role in the catalytic activity of a plethora of enzymes. However, cysteines are susceptible to oxidation by Reactive Oxygen and Nitrogen Species, leading to a loss of their catalytic function. Therefore, cysteine oxidation is emerging as a relevant physiological regulatory mechanism. Formation of a cyclic sulfenyl amide residue at the active site of redox-regulated proteins has been proposed as a protection mechanism against irreversible oxidation as the sulfenyl amide intermediate has been identified in several proteins. However, how and why only some specific cysteine residues in particular proteins react to form this intermediate is still unknown. In the present work using in-silico based tools, we have identified a constrained conformation that accelerates sulfenyl amide formation. By means of combined MD and QM/MM calculation we show that this conformation positions the NH backbone towards the sulfenic acid and promotes the reaction to yield the sulfenyl amide intermediate, in one step with the concomitant release of a water molecule. Moreover, in a large subset of the proteins we found a conserved beta sheet-loop-helix motif, which is present across different protein folds, that is key for sulfenyl amide production as it promotes the previous formation of sulfenic acid. For catalytic activity, in several cases, proteins need the Cysteine to be in the cysteinate form, i.e. a low pK_a Cys. We found that the conserved motif stabilizes the cysteinate by hydrogen bonding to several NH backbone moieties. As cysteinate is also more reactive toward ROS we propose that the sheet-loop-helix motif and the constraint conformation have been selected by evolution for proteins that need a reactive Cys protected from irreversible oxidation. Our results also highlight how fold conservation can be correlated to redox chemistry regulation of protein function.

Cysteine oxidation is emerging as a relevant regulatory mechanism of enzymatic function in the cell. Many proteins are protected from over oxidation by reactive oxygen species by the formation of a cyclic sulfenyl amide. Understanding how cyclic sulfenyl amide is formed and its dependence on protein structure is not only a basic question but necessary to predict which proteins may auto protect from over oxidation We describe a structural motif, which includes cysteine residues with a constrained conformation in a “forbidden” region of the Ramachandran plot plus a Beta-Cys-loop-helix motif, which has a reactive low pK_a Cysteine and also enables to form the cyclic sulfenyl amide with a low activation barrier. Our QM/MM computations show that the cyclization reaction only occurs if the “forbidden” conformation is acquired by the Cysteine residue. This structural motif was identified at least in 7 PFAM families and 145 proteins with solved structure, showing that a large number of proteins could have the ability to go through such cyclic product preventing irreversible oxidation.

FuncPatch: a web server for the fast Bayesian inference of conserved functional patches in protein 3D structures

MOTIVATION

A number of statistical phylogenetic methods have been developed to infer conserved functional sites or regions in proteins. Many methods, e.g. Rate4Site, apply the standard phylogenetic models to infer site-specific substitution rates and totally ignore the spatial correlation of substitution rates in protein tertiary structures, which may reduce their power to identify conserved functional patches in protein tertiary structures when the sequences used in the analysis are highly similar. The 3D sliding window method has been proposed to infer conserved functional patches in protein tertiary structures, but the window size, which reflects the strength of the spatial correlation, must be predefined and is not inferred from data. We recently developed GP4Rate to solve these problems under the Bayesian framework. Unfortunately, GP4Rate is computationally slow. Here, we present an intuitive web server, FuncPatch, to perform a fast approximate Bayesian inference of conserved functional patches in protein tertiary structures.

RESULTS

Both simulations and four case studies based on empirical data suggest that FuncPatch is a good approximation to GP4Rate. However, FuncPatch is orders of magnitudes faster than GP4Rate. In addition, simulations suggest that FuncPatch is potentially a useful tool complementary to Rate4Site, but the 3D sliding window method is less powerful than FuncPatch and Rate4Site. The functional patches predicted by FuncPatch in the four case studies are supported by experimental evidence, which corroborates the usefulness of FuncPatch.

AVAILABILITY AND IMPLEMENTATION

The software FuncPatch is freely available at the web site, http://info.mcmaster.ca/yifei/FuncPatch

CONTACT

golding@mcmaster.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

QM/MM study of the C-C coupling reaction mechanism of CYP121, an essential cytochrome p450 of Mycobacterium tuberculosis

Among 20 p450s of Mycobacterium tuberculosis (Mt), CYP121 has received an outstanding interest, not only due to its essentiality for bacterial viability but also because it catalyzes an unusual carbon-carbon coupling reaction. Based on the structure of the substrate bound enzyme, several reaction mechanisms were proposed involving first Tyr radical formation, second Tyr radical formation, and C-C coupling. Key and unknown features, being the nature of the species that generate the first and second radicals, and the role played by the protein scaffold each step. In the present work we have used classical and quantum based computer simulation methods to study in detail its reaction mechanism. Our results show that substrate binding promotes formation of the initial oxy complex, Compound I is the responsible for first Tyr radical formation, and that the second Tyr radical is formed subsequently, through a PCET reaction, promoted by the presence of key residue Arg386. The final C-C coupling reaction possibly occurs in bulk solution, thus yielding the product in one oxygen reduction cycle. Our results thus contribute to a better comprehension of MtCYP121 reaction mechanism, with direct implications for inhibitor design, and also contribute to our general understanding of these type of enzymes.

The role of the N-terminal tail for the oligomerization, folding and stability of human frataxin

The N-terminal stretch of human frataxin (hFXN) intermediate (residues 42–80) is not conserved throughout evolution and, under defined experimental conditions, behaves as a random-coil. Overexpression of hFXN56–210 in Escherichia coli yields a multimer, whereas the mature form of hFXN (hFXN81–210) is monomeric. Thus, cumulative experimental evidence points to the N-terminal moiety as an essential element for the assembly of a high molecular weight oligomer. The secondary structure propensity of peptide 56–81, the moiety putatively responsible for promoting protein–protein interactions, was also studied. Depending on the environment (TFE or SDS), this peptide adopts α-helical or β-strand structure. In this context, we explored the conformation and stability of hFXN56–210. The biophysical characterization by fluorescence, CD and SEC-FPLC shows that subunits are well folded, sharing similar stability to hFXN90–210. However, controlled proteolysis indicates that the N-terminal stretch is labile in the context of the multimer, whereas the FXN domain (residues 81–210) remains strongly resistant. In addition, guanidine hydrochloride at low concentration disrupts intermolecular interactions, shifting the ensemble toward the monomeric form. The conformational plasticity of the N-terminal tail might impart on hFXN the ability to act as a recognition signal as well as an oligomerization trigger. Understanding the fine-tuning of these activities and their resulting balance will bear direct relevance for ultimately comprehending hFXN function.

•

hFXN56–210 is well-folded and shares similar stability to hFXN90–210.

•

The oligomeric form of hFXN56–210 can be disassembled and reassembled in vitro.

•

Proteolysis leads to the oligomer disassembly: subunits are abridged to hFXN81–210.

•

Isolated peptide hFXN56–81 acquires structure in TFE and SDS solutions.

•

The N-terminal tail is structurally malleable and triggers oligomerization.

Molecular modelling and simulations in cancer research

The complexity of cancer and the vast amount of experimental data available have made computer-aided approaches necessary. Biomolecular modelling techniques are becoming increasingly easier to use, whereas hardware and software are becoming better and cheaper. Cross-talk between theoretical and experimental scientists dealing with cancer-research from a molecular approach, however, is still uncommon. This is in contrast to other fields, such as amyloid-related diseases, where molecular modelling studies are widely acknowledged. The aim of this review paper is therefore to expose some of the more common approaches in molecular modelling to cancer scientists in simple terms, illustrating success stories while also revealing the limitations of computational studies at the molecular level.

Molecular Dynamics Simulations Provide Atomistic Insight into Hydrogen Exchange Mass Spectrometry Experiments

It is now clear that proteins are flexible entities that in solution switch between conformations to achieve their function. Hydrogen/Deuterium Exchange Mass Spectrometry (HX/MS) is an invaluable tool to understand dynamic changes in proteins modulated by cofactor binding, post-transductional modifications, or protein-protein interactions. ERK2MAPK, a protein involved in highly conserved signal transduction pathways of paramount importance for normal cellular function, has been extensively studied by HX/MS. Experiments of the ERK2MAPK in the inactive and active states (in the presence or absence of bound ATP) have provided valuable information on the plasticity of the MAPK domain. However, interpretation of the HX/MS data is difficult, and changes are mostly explained in relation to available X-ray structures, precluding a complete atomic picture of protein dynamics. In the present work, we have used all atom Molecular Dynamics simulations (MD) to provide a theoretical framework for the interpretation of HX/MS data. Our results show that detailed analysis of protein-solvent interaction along the MD simulations allows (i) prediction of the number of protons exchanged for each peptide in the HX/MS experiments, (ii) rationalization of the experimentally observed changes in exchange rates in different protein conditions at the residue level, and (iii) that at least for ERK2MAPK, most of the functionally observed differences in protein dynamics are related to what can be considered the native state conformational ensemble. In summary, the combination of HX/MS experiments with all atom MD simulations emerges as a powerful approach to study protein native state dynamics with atomic resolution.

Modeling methods for studying post-translational and transcriptional modifying enzymes

Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.

Phosphorylated mitogen-activated protein kinase/extracellular signal-regulated kinase 1/2 may not always represent its kinase activity in a rat model of focal cerebral ischemia with or without ischemic preconditioning

The extracellular signal-regulated kinase (ERK) 1/2 protein requires a dual phosphorylation at conserved threonine and tyrosine residues to be fully activated under normal physiological conditions. Thus, ERK1/2 kinase activity is often defined by the quantity of phosphorylated kinase. However, this may not accurately represent its true activity under certain pathological conditions. We investigated whether ERK1/2 kinase activity is proportional to its phosphorylation state in a rat focal ischemia model with and without rapid ischemic preconditioning. We showed that phosphorylated-ERK1/2 protein levels were increased 2.6±0.07-fold, and ERK1/2 kinase activity was increased 10.6±1.9-fold in animals receiving ischemic preconditioning alone without test ischemia compared with sham group (P<0.05, n=6/group), suggesting that phosphorylated-ERK1/2 protein levels represent its kinase activity under these conditions. However, preconditioning plus test ischemia robustly blocked ERK1/2 kinase activity, whereas it increased phosphorylated-ERK1/2 protein levels beyond those receiving test ischemia alone, suggesting that phosphorylated-ERK1/2 protein levels were not representative of actual kinase activity in this pathological condition. In conclusion, protein phosphorylation levels of ERK1/2 do not always correspond to kinase activity, thus, measuring the true kinase activity is essential.

Insights into the phosphoryl transfer mechanism of cyclin-dependent protein kinases from ab initio QM/MM free-energy studies

Phosphorylation reactions catalyzed by kinases and phosphatases play an indispensible role in cellular signaling, and their malfunctioning is implicated in many diseases. A better understanding of the catalytic mechanism will help design novel and effective mechanism-based inhibitors of these enzymes. In this work, ab initio quantum mechanical/molecular mechanical studies are reported for the phosphoryl transfer reaction catalyzed by a cyclin-dependent kinase, CDK2. Our results suggest that an active-site Asp residue, rather than ATP as previously proposed, serves as the general base to activate the Ser nucleophile. The corresponding transition state features a dissociative, metaphosphate-like structure, stabilized by the Mg(2+) ion and several hydrogen bonds. The calculated free-energy barrier is consistent with experimental values. Implications of our results in this and other protein kinases are discussed.

Cation binding to 15-TBA quadruplex DNA is a multiple-pathway cation-dependent process

A combination of explicit solvent molecular dynamics simulation (30 simulations reaching 4 µs in total), hybrid quantum mechanics/molecular mechanics approach and isothermal titration calorimetry was used to investigate the atomistic picture of ion binding to 15-mer thrombin-binding quadruplex DNA (G-DNA) aptamer. Binding of ions to G-DNA is complex multiple pathway process, which is strongly affected by the type of the cation. The individual ion-binding events are substantially modulated by the connecting loops of the aptamer, which play several roles. They stabilize the molecule during time periods when the bound ions are not present, they modulate the route of the ion into the stem and they also stabilize the internal ions by closing the gates through which the ions enter the quadruplex. Using our extensive simulations, we for the first time observed full spontaneous exchange of internal cation between quadruplex molecule and bulk solvent at atomistic resolution. The simulation suggests that expulsion of the internally bound ion is correlated with initial binding of the incoming ion. The incoming ion then readily replaces the bound ion while minimizing any destabilization of the solute molecule during the exchange.

A conserved non-canonical motif in the pseudoactive site of the ROP5 pseudokinase domain mediates its effect on Toxoplasma virulence

The ROP5 family is a closely related set of polymorphic pseudokinases that are critical to the ability of Toxoplasma to cause disease. Polymorphisms in ROP5 also make it a major determinant of strain-specific differences in virulence. ROP5 possesses all of the major kinase motifs required for catalysis except for a substitution at the catalytic Asp. We show that this substitution in the catalytic loop of ROP5 is part of a motif conserved in other pseudokinases of both Toxoplasma and human origin, and that this motif is required for the full activity in vivo of ROP5. This suggests evolutionary selection at this site for a biochemical function, rather than simple drift away from catalysis. We present the crystal structures of a virulent isoform of ROP5 both in its ATP-bound and -unbound states and have demonstrated that despite maintaining the canonical ATP-binding motifs, ROP5 binds ATP in a distorted conformation mediated by unusual magnesium coordination sites that would not be predicted from the primary sequence. In addition, we have mapped the polymorphisms spread throughout the primary sequence of ROP5 to two major surfaces, including the activation segment of ROP5. This suggests that the pseudoactive site of this class of pseudokinases may have evolved to use the canonical ATP-binding motifs for non-catalytic signaling through allostery.

Using peptides to study protein-protein interactions

Protein-protein interactions (PPIs) govern all aspects of cell function and, as such, are a major target for research and therapeutic intervention. A major rate-limiting step in PPI research is the expression and purification of full-length proteins. The use of peptides to study PPIs significantly facilitates the structural and biophysical characterization of PPIs as well as the effort to develop drugs to control PPIs. Here we describe examples for the use of peptides to study PPI and some of the important experimental methods that are used in the field. Peptides have proved to be excellent tools to study PPIs and have been contributing both for understanding mechanisms of PPIs as well as for drug design for PPI modulation.

Prioritization of charge over geometry in transition state analogues of a dual specificity protein kinase

The direct observation of a transition state analogue (TSA) complex for tyrosine phosphorylation by a signaling kinase has been achieved using (19)F NMR analysis of MEK6 in complex with tetrafluoroaluminate (AlF(4)(-)), ADP, and p38α MAP kinase (acceptor residue: Tyr182). Solvent-induced isotope shifts and chemical shifts for the AlF(4)(-) moiety indicate that two fluorine atoms are coordinated by the two catalytic magnesium ions of the kinase active site, while the two remaining fluorides are liganded by protein residues only. An equivalent, yet distinct, AlF(4)(-) complex involving the alternative acceptor residue in p38α (Thr180) is only observed when the Tyr182 is mutated to phenylalanine. The formation of octahedral AlF(4)(-) species for both acceptor residues, rather than the trigonal bipyramidal AlF(3)(0) previously identified in the only other metal fluoride complex with a protein kinase, shows the requirement of MEK6 for a TSA that is isoelectronic with the migrating phosphoryl group. This requirement has hitherto only been demonstrated for proteins having a single catalytic magnesium ion.

p38γ activation triggers dynamical changes in allosteric docking sites

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

Noncatalytic function of ERK1/2 can promote Raf/MEK/ERK-mediated growth arrest signaling

Kinase activity is known as the key biochemical property of MAPKs. Here, we report that ERK1/2 also utilizes its noncatalytic function to mediate certain signal transductions. Sustained activation of the Raf/MEK/ERK pathway induces growth arrest, accompanied by changes in cell cycle regulators (decreased retinoblastoma phosphorylation, E2F1 down-regulation, and/or p21(CIP1) up-regulation) and cell type-specific changes in morphology and expression of c-Myc or RET in the human tumor lines LNCaP, U251, and TT. Ablation of ERK1/2 by RNA interference abrogated all these effects. However, active site-disabled ERK mutants (ERK1-K71R, ERK2-K52R, and ERK2-D147A), which competitively inhibit activation of endogenous ERK1/2, could not block Raf/MEK-induced growth arrest as well as changes in the cell cycle regulators, although they effectively blocked phosphorylation of the ERK1/2 catalytic activity readouts, p90(RSK) and ELK1, as well as the cell type-specific changes. Because this indicated a potential noncatalytic ERK1/2 function, we generated stable lines of the tumor cells in which both ERK1 and ERK2 were significantly knocked down, and we further investigated the possibility using rat-derived kinase-deficient ERK mutants (ERK2-K52R and ERK2-T183A/Y185F) that were not targeted by human small hairpin RNA. Indeed, ERK2-K52R selectively restored Raf-induced growth inhibitory signaling in ERK1/2-depleted cells, as manifested by regained cellular ability to undergo growth arrest and to control the cell cycle regulators without affecting c-Myc and morphology. However, ERK2-T183A/Y185F was less effective, indicating the requirement of TEY site phosphorylation. Our study suggests that functions of ERK1/2 other than its "canonical" kinase activity are also involved in the pathway-mediated growth arrest signaling.