The mechanism of inhibition of the cyclin-dependent kinase-2 as revealed by the molecular dynamics study on the complex CDK2 with the peptide substrate HHASPRK

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Identification of novel CDK2 inhibitors by a multistage virtual screening method based on SVM, pharmacophore and docking model

Cyclin-dependent kinase 2 (CDK2) is the family of Ser/Thr protein kinases that has emerged as a highly selective with low toxic cancer therapy target. A multistage virtual screening method combined by SVM, protein-ligand interaction fingerprints (PLIF) pharmacophore and docking was utilised for screening the CDK2 inhibitors. The evaluation of the validation set indicated that this method can be used to screen large chemical databases because it has a high hit-rate and enrichment factor (80.1% and 332.83 respectively). Six compounds were screened out from NCI, Enamine and Pubchem database. After molecular dynamics and binding free energy calculation, two compounds had great potential as novel CDK2 inhibitors and they also showed selective inhibition against CDK2 in the kinase activity assay.

Computational Strategy Revealing the Structural Determinant of Ligand Selectivity towards Highly Similar Protein Targets

BACKGROUND

Poor selectivity of drug candidates may lead to toxicity and side effects accounting for as high as 60% failure rate, thus, the selectivity is consistently significant and challenging for drug discovery.

OBJECTIVE

To find highly specific small molecules towards very similar protein targets, multiple strategies are always employed, including (1) To make use of the diverse shape of binding pocket to avoid steric bump; (2) To increase binding affinities for favorite residues; (3) To achieve selectivity through allosteric regulation of target; (4) To stabalize the inactive conformation of protein target and (5) To occupy dual binding pockets of single target.

CONCLUSION

In this review, we summarize computational strategies along with examples of their successful applications in designing selective ligands, with the aim to provide insights into everdiversifying drug development practice and inspire medicinal chemists to utilize computational strategies to avoid potential side effects due to low selectivity of ligands.

Ensemble-based modeling and rigidity decomposition of allosteric interaction networks and communication pathways in cyclin-dependent kinases: Differentiating kinase clients of the Hsp90-Cdc37 chaperone

The overarching goal of delineating molecular principles underlying differentiation of protein kinase clients and chaperone-based modulation of kinase activity is fundamental to understanding activity of many oncogenic kinases that require chaperoning of Hsp70 and Hsp90 systems to attain a functionally competent active form. Despite structural similarities and common activation mechanisms shared by cyclin-dependent kinase (CDK) proteins, members of this family can exhibit vastly different chaperone preferences. The molecular determinants underlying chaperone dependencies of protein kinases are not fully understood as structurally similar kinases may often elicit distinct regulatory responses to the chaperone. The regulatory divergences observed for members of CDK family are of particular interest as functional diversification among these kinases may be related to variations in chaperone dependencies and can be exploited in drug discovery of personalized therapeutic agents. In this work, we report the results of a computational investigation of several members of CDK family (CDK5, CDK6, CDK9) that represented a broad repertoire of chaperone dependencies—from nonclient CDK5, to weak client CDK6, and strong client CDK9. By using molecular simulations of multiple crystal structures we characterized conformational ensembles and collective dynamics of CDK proteins. We found that the elevated dynamics of CDK9 can trigger imbalances in cooperative collective motions and reduce stability of the active fold, thus creating a cascade of favorable conditions for chaperone intervention. The ensemble-based modeling of residue interaction networks and community analysis determined how differences in modularity of allosteric networks and topography of communication pathways can be linked with the client status of CDK proteins. This analysis unveiled depleted modularity of the allosteric network in CDK9 that alters distribution of communication pathways and leads to impaired signaling in the client kinase. According to our results, these network features may uniquely define chaperone dependencies of CDK clients. The perturbation response scanning and rigidity decomposition approaches identified regulatory hotspots that mediate differences in stability and cooperativity of allosteric interaction networks in the CDK structures. By combining these synergistic approaches, our study revealed dynamic and network signatures that can differentiate kinase clients and rationalize subtle divergences in the activation mechanisms of CDK family members. The therapeutic implications of these results are illustrated by identifying structural hotspots of pathogenic mutations that preferentially target regions of the increased flexibility to enable modulation of activation changes. Our study offers a network-based perspective on dynamic kinase mechanisms and drug design by unravelling relationships between protein kinase dynamics, allosteric communications and chaperone dependencies.

Palbociclib can overcome mutations in cyclin dependent kinase 6 that break hydrogen bonds between the drug and the protein

Inhibition of cyclin dependent kinases (CDKs) 4 and 6 prevent cells from entering the synthesis phase of the cell cycle. CDK4 and 6 are therefore important drug targets in various cancers. The selective CDK4/6 inhibitor palbociclib is approved for the treatment of breast cancer and has shown activity in a cellular model of mixed lineage leukaemia (MLL)‐rearranged acute myeloid leukaemia (AML). We studied the interactions of palbociclib and CDK6 using molecular dynamics simulations. Analysis of the simulations suggested several interactions that stabilized the drug in its binding site and that were not observed in the crystal structure of the protein‐drug complex. These included a hydrogen bond to His 100 that was hitherto not reported and several hydrophobic contacts. Evolutionary‐based bioinformatic analysis was used to suggest two mutants, D163G and H100L that would potentially yield drug resistance, as they lead to loss of important protein–drug interactions without hindering the viability of the protein. One of the mutants involved a change in the glycine of the well‐conserved DFG motif of the kinase. Interestingly, CDK6‐dependent human AML cells stably expressing either mutant retained sensitivity to palbociclib, indicating that the protein‐drug interactions are not affected by these. Furthermore, the cells were proliferative in the absence of palbociclib, indicating that the Asp to Gly mutation in the DFG motif did not interfere with the catalytic activity of the protein.

PDB Code(s): 2EUF

Network-based modelling and percolation analysis of conformational dynamics and activation in the CDK2 and CDK4 proteins: dynamic and energetic polarization of the kinase lobes may determine divergence of the regulatory mechanisms

Network modeling and percolation analysis of conformational dynamics and energetics of regulatory mechanisms in cyclin-dependent kinases.

Mechanism of p27 Unfolding for CDK2 Reactivation

Cell-cycle regulatory protein, CDK2 is active when bound to its complementary partner protein, CyclinA or E. Recent discovery of the Kip/Cip family of proteins has indicated that the activity of CDK2 is also regulated by a member protein, p27. Although, the mechanism of CDK2 inhibition by p27 binding is known from crystal structure, little is known about the mechanism of CDK2 reactivation. Here, we execute classical and accelerated molecular dynamics simulations of unphosphorylated- and phosphorylated-p27 bound CDK2/CyclinA to unravel the CDK2 reactivation mechanism at molecular-to-atomic detail. Results suggest that the phosphorylation of p27 Y88 residue (pY88-p27) first disrupts the p27/CDK2 hybrid β-sheet and subsequently ejects the p27 3₁₀ helix from CDK2 catalytic cleft. The unbinding of p27 from CDK2/CyclinA complex, thus, follows a two-step unfolding mechanism, where the 3₁₀ helix ejection constitutes the rate-limiting step. Interestingly, the unfolding of p27 leaves CDK2/CyclinA in an active state, where the prerequisite CDK2-CyclinA interfacial contacts were regained and ATP achieved its native position for smooth transfer of phosphate. Our findings match very well with NMR chemical shift data that indicated the flip-out of p27 3₁₀ helix from CDK2 pocket and kinetic experiments that exhibited significant kinase activity of CDK2 when saturated with pY88-p27.

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.

Atrial natriuretic factor receptor guanylate cyclase, ANF-RGC, transduces two independent signals, ANF and Ca(2+)

Atrial natriuretic factor receptor guanylate cyclase (ANF-RGC), was the first discovered member of the mammalian membrane guanylate cyclase family. The hallmark feature of the family is that a single protein contains both the site for recognition of the regulatory signal and the ability to transduce it into the production of the second messenger, cyclic GMP. For over two decades, the family has been classified into two subfamilies, the hormone receptor subfamily with ANF-RGC being its paramount member, and the Ca²⁺ modulated subfamily, which includes the rod outer segment guanylate cyclases, ROS-GC1 and 2, and the olfactory neuroepithelial guanylate cyclase. ANF-RGC is the receptor and the signal transducer of the most hypotensive hormones, ANF– and B-type natriuretic peptide (BNP). After binding these hormones at the extracellular domain it, at its intracellular domain, signals activation of the C-terminal catalytic module and accelerates the production of cyclic GMP. Cyclic GMP then serves the second messenger role in biological responses of ANF and BNP such as natriuresis, diuresis, vasorelaxation, and anti-proliferation. Very recently another modus operandus for ANF-RGC was revealed. Its crux is that ANF-RGC activity is also regulated by Ca²⁺. The Ca²⁺ sensor neurocalcin d mediates this signaling mechanism. Strikingly, the Ca²⁺ and ANF signaling mechanisms employ separate structural motifs of ANF-RGC in modulating its core catalytic domain in accelerating the production of cyclic GMP. In this review the biochemistry and physiology of these mechanisms with emphasis on cardiovascular regulation will be discussed.

Tyrosine phosphorylation allows integration of multiple signaling inputs by IKKβ

Signaling regulated by NFκB and related transcription factors is centrally important to many inflammatory and autoimmune diseases, cancer, and stress responses. The kinase that directly regulates the canonical NFκB transcriptional pathway, Inhibitor of κB kinase β (IKKβ), undergoes activation by Ser phosphorylation mediated by NIK or TAK1 in response to inflammatory signals. Using titanium dioxide-based phosphopeptide enrichment (TiO2)-liquid chromatography (LC)-high mass accuracy tandem mass spectrometry (MS/MS), we analyzed IKKβ phosphorylation in human HEK293 cells expressing IKKβ and FGFR2, a Receptor tyrosine kinase (RTK) essential for embryonic differentiation and dysregulated in several cancers. We attained unusually high coverage of IKKβ, identifying an abundant site of Tyr phosphorylation at Tyr169 within the Activation Loop. The phosphomimic at this site confers a level of kinase activation and NFκB nuclear localization exceeding the iconic mutant S177E/S181E, demonstrating that RTK-mediated Tyr phosphorylation of IKKβ has the potential to directly regulate NFκB transcriptional activation.

Molecular basis of differential selectivity of cyclobutyl-substituted imidazole inhibitors against CDKs: insights for rational drug design

Cyclin-dependent kinases (CDKs) belong to the CMGC subfamily of protein kinases and play crucial roles in eukaryotic cell division cycle. At least seven different CDKs have been reported to be implicated in the cell cycle regulation in vertebrates. These CDKs are highly homologous and contain a conserved catalytic core. This makes the design of inhibitors specific for a particular CDK difficult. There is, however, growing need for CDK5 specific inhibitors to treat various neurodegenerative diseases. Recently, cis-substituted cyclobutyl-4-aminoimidazole inhibitors have been identified as potent CDK5 inhibitors that gave up to 30-fold selectivity over CDK2. Available IC50 values also indicate a higher potency of this class of inhibitors over commercially available drugs, such as roscovitine. To understand the molecular basis of higher potency and selectivity of these inhibitors, here, we present molecular dynamics simulation results of CDK5/p25 and CDK2/CyclinE complexed with a series of cyclobutyl-substituted imidazole inhibitors and roscovitine. The atomic details of the stereospecificity and selectivity of these inhibitors are obtained from energetics and binding characteristics to the CDK binding pocket. The study not only complements the experimental findings, but also provides a wealth of detailed information that could help the structure-based drug designing processes.

How calcium inhibits the magnesium-dependent kinase gsk3β: a molecular simulation study

Glycogen synthase kinase 3β (GSK3β) is a ubiquitous serine/threonine kinase that plays a pivotal role in many biological processes. GSK3β catalyzes the transfer of γ-phosphate of ATP to the unique substrate Ser/Thr residues with the assistance of two natural activating cofactors Mg(2+). Interestingly, the biological observation reveals that a non-native Ca(2+) ion can inhibit the GSK3β catalytic activity. Here, the inhibitory mechanism of GSK3β by the displacement of native Mg(2+) at site 1 by Ca(2+) was investigated by means of 80 ns comparative molecular dynamics (MD) simulations of the GSK3β···Mg(2+)-2/ATP/Mg(2+) -1 and GSK3β···Mg(2+)-2/ATP/Ca(2+)-1 systems. MD simulation results revealed that using the AMBER point charge model force field for Mg(2+) was more appropriate in the reproduction of the active site architectural characteristics of GSK3β than using the magnesium-cationic dummy atom model force field. Compared with the native Mg(2+) bound system, the misalignment of the critical triphosphate moiety of ATP, the erroneous coordination environments around the Mg(2+) ion at site 2, and the rupture of the key hydrogen bond between the invariant Lys85 and the ATP O(β2) atom in the Ca(2+) substituted system were observed in the MD simulation due to the Ca(2+) ion in active site in order to achieve its preferred sevenfold coordination geometry, which adequately abolish the enzymatic activity. The obtained results are valuable in understanding the possible mechanism by why Ca(2+) inhibits the GSK3β activity and also provide insights into the mechanism of Ca(2+) inhibition in other structurally related protein kinases.

Insights into the role of magnesium triad in myo-inositol monophosphatase: metal mechanism, substrate binding, and lithium therapy

myo-Inositol monophosphatase (IMPase) plays a pivotal role in the intracellular phosphatidylinositol cell signaling pathway. It has attracted considerable attention as a putative therapeutic target for lithium therapy in the treatment of bipolar disorder. A trio of activated cofactor Mg²⁺ ions is required for inositol monophosphate hydrolysis by IMPase. In the present study, computational studies, including two-layered ONIOM-based quantum mechanics/mechanical mechanics (QM/MM) calculations, molecular modeling, and molecular dynamics (MD) simulations, were performed to ascertain the role of the Mg²⁺ triad in the IMPase active site. The QM/MM calculations show that the structural identity of the nucleophilic water molecule W1 shared by Mg²⁺-1 and Mg²⁺-3, activated by Thr95/Asp47 dyad, is a hydroxide ion. Moreover, Mg²⁺-3 needs to be conjugated with Mg²⁺-1 in the binding site to create the activated nucleophilic hydroxide ion in accordance with the three-metal ion catalytic mechanism. The MD simulation of the IMPase-substrate-Mg²⁺ complex shows that the three Mg²⁺ ions promote substrate binding and help fix the phosphate moiety of the substrate for nucleophilic attack by the hydroxide ion. When Mg²⁺-2 is displaced with Li⁺, the MD simulations of the postreaction complex indicate that the conformation of the catalytic loop (residues 33 to 44) is disrupted and water molecule W2 does not coordinate with Li⁺. This disruption traps the inorganic phosphate and inositolate in the active site, which lead to IMPase inhibition. By contrast, in the native Mg²⁺ system, the W2 ligated by Mg²⁺-2 and Asp200 aids in protonation of the leaving inositolate moiety.

Molecular modeling and molecular dynamics simulation studies of the GSK3β/ATP/substrate complex: understanding the unique P+4 primed phosphorylation specificity for GSK3β substrates

Substrate specificity of protein kinases is of fundamental importance for the integrity and fidelity of signaling pathways. Glycogen synthase kinase 3β (GSK3β) has a unique substrate specificity that prefers phosphorylation of its substrates at the P+4 serine before it can further phosphorylate the substrate at the P0 serine in the canonical motif SXXXS(p), where S(p) is the primed phosphorylation site. The detailed phosphorylation mechanism, however, is not clearly understood. In this study, a three-dimensional (3D) model of the ternary complex of GSK3β, ATP, and the phosphorylated glycogen synthase (pGS), termed GSK3β/ATP/pGS, is constructed using a hierarchical approach and by integrating molecular modeling and molecular dynamics (MD) simulations. Based on the 3D model, the substrate primed phosphorylation mechanism is investigated via two 12 ns comparative MD simulations of the GSK3β/ATP/pGS and GSK3β/ATP/GS systems, which differ in the phosphate group bound to the P+4 serine of GS. In agreement with structural analysis, computed binding free energies reveal that the binding of pGS to GSK3β is favored in the prephosphorylated state compared with the GS native state. More importantly, comparison with the system simulated without primed phosphorylation in the GSK3β/ATP/GS complex shows that for an optimal phosphorylation reaction to occur, the pGS priming phosphate in the GSK3β/ATP/pGS system optimizes the proper orientation of the GSK3β N- and C-terminal domains and clamps the P0 serine of pGS in the appropriate configuration for interaction with the ATP γ-phosphate within the catalytic groove.

Mechanism of kinase inactivation and nonbinding of FRATide to GSK3β due to K85M mutation: molecular dynamics simulation and normal mode analysis

As a serine/threonine protein kinase, glycogen synthase kinase 3β (GSK3β) is an essential component of several cellular processes, including insulin, growth factor, and Wnt signaling. The conserved K85 is important to GSK3β activity and FRATide binding. To elucidate the mechanisms concerning kinase inactivation and nonbinding of FRATide to GSK3β, molecular dynamics (MD) simulation, molecular mechanics generalized Born/surface area (MM_GBSA) calculation, and normal mode analysis (NMA) were performed on both the wild-type (WT) and the K85M mutation of the GSK3β-FRATide complex. The results revealed that the periodic open-closed conformational change of the G loop, together with the compact conformation of the RD pocket, was disturbed in the K85M mutant, in contrast to those in the WT. This in turn caused inhibition of GSK3β. Specifically, the correct folding pattern of GSK3β was disrupted in the K85M mutant, resulting in the loss of two key hydrogen bonds between K214 of FRATide and E290 and K292 of GSK3β, respectively. Furthermore, MM_GBSA calculations indicated that the K85M mutation could lead to a less energy-favorable GSK3β-FRATide complex. In addition, NMA demonstrated that the "rocking" of the N- and C-terminal domains of GSK3β, which coordinates the mutual movement of both lobes, inducing the opening and closing of the active site of GSK3β, which may assist the entry of ATP into the ATP binding site and the release of the ADP product. Strikingly, this phenomenon was not clearly observed in the K85M mutation. This study provides a structural basis for the effect of the K85M mutation on the GSK3β-FRATide complex.

Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact

Conventional molecular dynamics simulations on 50 ns to 1 μs time scales were used to study the effects of explicit solvent models on the conformational behavior and solvation of two oligopeptide solutes: α-helical EK-peptide (14 amino acids) and a β-hairpin chignolin (10 amino acids). The widely used AMBER force fields (ff99, ff99SB, and ff03) were combined with four of the most commonly used explicit solvent models (TIP3P, TIP4P, TIP5P, and SPC/E). Significant differences in the specific solvation of chignolin among the studied water models were identified. Chignolin was highly solvated in TIP5P, whereas reduced specific solvation was found in the TIP4P, SPC/E, and TIP3P models for kinetic, thermodynamic, and both kinetic and thermodynamic reasons, respectively. The differences in specific solvation did not influence the dynamics of structured parts of the folded peptide. However, substantial differences between TIP5P and the other models were observed in the dynamics of unfolded chignolin, stability of salt bridges, and specific solvation of the backbone carbonyls of EK-peptide. Thus, we conclude that the choice of water model may affect the dynamics of flexible parts of proteins that are solvent-exposed. On the other hand, all water models should perform similarly for well-structured folded protein regions. The merits of the TIP3P model include its high and overestimated mobility, which accelerates simulation processes and thus effectively increases sampling.

The effect of Li+ on GSK-3 inhibition: molecular dynamics simulation

Glycogen synthase kinase-3 (GSK-3) is a kind of serine-threonine protein kinase. It places important roles in several signaling pathways and it is a key therapeutic target for a number of diseases, such as diabetes, cancer, Alzheimer's disease and chronic inflammation. Mg(2+) ions which interact with ATP are conserved in GSK. They are important in phosphoryl transfer. Li(+) is an inhibitor for GSK-3. It is used to treat bipolar mood disorder. This paper illustrates the effect of Li(+) on GSK-3. When Mg(I)(2+) is replaced by Li(+), the atom fluctuation of GSK-3 will rise, and the in-line phosphoryl transfer mechanism is probably demolished and the binding of pre-phosphorylated substrates may be disturbed. All the results we obtained clearly suggest that inhibition to GSK-3 is caused by the Mg(I)(2+) replacement with Li(+).

A cancer-derived mutation in the PSTAIRE helix of cyclin-dependent kinase 2 alters the stability of cyclin binding

Cyclin-dependent kinase 2 (cdk2) is a central regulator of the mammalian cell cycle. Here we describe the properties of a mutant form of cdk2 identified during large-scale sequencing of protein kinases from cancerous tissue. The mutation substituted a leucine for a proline in the PSTAIRE helix, the central motif in the interaction of the cdk with its regulatory cyclin subunit. We demonstrate that whilst the mutant cdk2 is considerably impaired in stable cyclin association, it is still able to generate an active kinase that can functionally complement defective cdks in vivo. Molecular dynamic simulations and biophysical measurements indicate that the observed biochemical properties likely stem from increased flexibility within the cyclin-binding helix.

Glycine-rich loop of mitochondrial processing peptidase alpha-subunit is responsible for substrate recognition by a mechanism analogous to mitochondrial receptor Tom20

Tryptophan fluorescence measurements were used to characterize the local dynamics of the highly conserved glycine-rich loop (GRL) of the mitochondrial processing peptidase (MPP) alpha-subunit in the presence of the substrate precursor. Reporter tryptophan residue was introduced into the GRL of the yeast alpha-MPP (Y299W) or at a proximal site (Y303W). Time-resolved and steady-state fluorescence spectroscopy demonstrated that for Trp299, the primary contact with the yeast malate dehydrogenase precursor evokes a change of the local GRL mobility. Moreover, time-resolved measurements showed that a functionless alpha-MPP with a single-residue deletion in the loop (Y303W/DeltaG292) is defective particularly in the primary contact with substrate. Thus, the GRL was proved to be part of a contact site of the enzyme specifically recognizing the substrate. Regarding the surface exposure and presence of the hydrophobic patches within the GRL, we proposed a functional analogy between the presequence recognition by the hydrophobic binding groove of the Tom20 mitochondrial import receptor and the GRL of the alpha-MPP. A molecular dynamics (MD) simulation of the MPP-substrate peptide complex model was employed to test this hypothesis. The initial positioning and conformation of the substrate peptide in the model fitting were chosen based on the analogy of its interaction with the Tom20 binding groove. MD simulation confirmed the stability of the proposed interaction and showed also a decrease in GRL flexibility in the presence of substrate, in agreement with fluorescence measurements. Moreover, conserved substrate hydrophobic residues in positions +1 and -4 to the cleavage site remain in close contact with the side chains of the GRL during the entire production part of MD simulation as stabilizing points of the hydrophobic interaction. We conclude that the GRL of the MPP alpha-subunit is the crucial evolutional outcome of the presequence recognition by MPP and represents a functional parallel with Tom20 import receptor.

Structural basis for the complete loss of GSK3beta catalytic activity due to R96 mutation investigated by molecular dynamics study

Many Ser/Thr protein kinases, to be fully activated, are obligated to introduce a phospho-Ser/Thr in their activation loop. Presently, the similarity of activation loop between two crystal complexes, i.e. glycogen synthase kinase 3beta (GSK3beta)-AMPNP and GSK3beta-sulfate ion complex, indicates that the activation segment of GSK3beta is preformed requiring neither a phosphorylation event nor conformational changes. GSK3beta, when participated in glycogen synthesis and Wnt signaling pathways, possesses a unique feature with the preference of such substrate with a priming phosphate. Experimental mutagenesis proved that the residue arginine at amino acid 96 mutations to lysine (R96K) or alanine (R96A) selectively abolish activity on the substrates involved in glycogen synthesis signaling pathway. Based on two solved crystal structures, wild type (WT) and two mutants (R96K and R96A) GSK3beta-ATP-phospho-Serine (pSer) complexes were modeled. Molecular dynamics simulations and energy analysis were employed to investigate the effect of pSer involvement on the GSK3beta structure in WT, and the mechanisms of GSK3beta deactivation due to R96K and R96A mutations. The results indicate that the introduction of pSer to WT GSK3beta generates a slight lobe closure on GSK3beta without any remarkable changes, which may illuminate the experimental conclusion, whereas the conformations of GSK3beta and ATP undergo significant changes in two mutants. As to GSK3beta, the affected positions distribute over activation loop, alpha-helix, and glycine-rich loop. Based on coupling among the mentioned positions, the allosteric mechanisms for distorted ATP were proposed. Energy decomposition on the residues of activation loop identified the important residues Arg96 and Arg180 in anchoring the phosphate group.

Two C-terminal ankyrin repeats form the minimal stable unit of the ankyrin repeat protein p18INK4c

Ankyrin repeat proteins (ARPs) appear to be abundant in organisms from all phyla, and play critical regulatory roles, mediating specific interactions with target biomolecules and thus ordering the sequence of events in diverse cellular processes. ARPs possess a non-globular scaffold consisting of repeating motifs named ankyrin (ANK) repeats, which stack on each other. The modular architecture of ARPs provides a new paradigm for understanding protein stability and folding mechanisms. In the present study, the stability of various C-terminal fragments of the ARP p18(INK4c) was investigated by all-atomic 450 ns molecular dynamics (MD) simulations in explicit water solvent. Only motifs with at least two ANK repeats made stable systems in the available timescale. All smaller fragments were unstable, readily losing their native fold and alpha-helical content. Since each non-terminal ANK repeat has two hydrophobic sides, we may hypothesize that at least one hydrophobic side must be fully covered and shielded from the water as a necessary, but not sufficient, condition to maintain ANK repeat stability. Consequently, at least two ANK repeats are required to make a stable ARP.

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15. Phosphorylation at Thr14/Tyr15 causes ATP misalignment for the phosphate-group transfer, changes in the Mg(2+) coordination sphere, and changes in the H-bond network formed by CDK2 catalytic residues (Asp127, Lys129, Asn132). The inhibitory phosphorylation causes the G-loop to shift from the ATP binding site, which leads to opening of the CDK2 substrate binding box, thus probably weakening substrate binding. All these effects explain the decrease in kinase activity observed after inhibitory phosphorylation at Thr14/Tyr15 in the G-loop. Interaction of the peptide substrate, and the phosphorylated peptide product, with CDK2 was also studied and compared. These results broaden hypotheses drawn from our previous MD studies as to why a basic residue (Arg/Lys) is preferred at the P(+2) substrate position.

Functional flexibility of human cyclin-dependent kinase-2 and its evolutionary conservation

Cyclin-dependent kinase 2 (CDK2) is the most thoroughly studied of the cyclin-dependent kinases that regulate essential cellular processes, including the cell cycle, and it has become a model for studies of regulatory mechanisms at the molecular level. This contribution identifies flexible and rigid regions of CDK2 based on temperature B-factors acquired from both X-ray data and molecular dynamics simulations. In addition, the biological relevance of the identified flexible regions and their motions is explored using information from the essential dynamics analysis related to conformational changes of CDK2 and knowledge of its biological function(s). The conserved regions of CMGC protein kinases' primary sequences are located in the most rigid regions identified in our analyses, with the sole exception of the absolutely conserved G13 in the tip of the glycine-rich loop. The conserved rigid regions are important for nucleotide binding, catalysis, and substrate recognition. In contrast, the most flexible regions correlate with those where large conformational changes occur during CDK2 regulation processes. The rigid regions flank and form a rigid skeleton for the flexible regions, which appear to provide the plasticity required for CDK2 regulation. Unlike the rigid regions (which as mentioned are highly conserved) no evidence of evolutionary conservation was found for the flexible regions.

Significance of water molecules in the inhibition of cylin-dependent kinase 2 and 5 complexes

Interest in CDK2 and CDK5 has stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. In the present paper, eight Molecular Dynamics (MD) simulations are carried out to examine the importance of structure and dynamics of water in the active site of both CDK2 and CDK5 complexes with roscovitine and indirubin analogues. Together with previous results, the current work shows a highly conserved water-involved hydrogen bonding (HB) network in both CDK2- and CDK5-indirubin combinations to complete information from the X-ray crystallography. The simulations suggest the importance of such a network for combining the inhibitor to the host protein as well as the significance of using an activated CDK as a template when designing new inhibitors. Different binding patterns of roscovitine in CDK2 and CDK5 are detected during the simulations because of the different binding conformations of the group on the C2 side chain, which might offer a clue toward finding highly selective inhibitors with regards to CDK2 and CDK5.

Molecular dynamics simulations on the inhibition of cyclin-dependent kinases 2 and 5 in the presence of activators

Interests in CDK2 and CDK5 have stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. Molecular dynamics (MD) simulations have not only become a viable approach to drug design because of advances in computer technology but are increasingly an integral part of drug discovery processes. It is common in MD simulations of inhibitor/CDK complexes to exclude the activator of the CDKs in the structural models to keep computational time tractable. In this paper, we present simulation results of CDK2 and CDK5 with roscovitine using models with and without their activators (cyclinA and p25). While p25 was found to induce slight changes in CDK5, the calculations support that cyclinA leads to significant conformational changes near the active site of CDK2. This suggests that detailed and structure-based inhibitor design targeted at these CDKs should employ activator-included models of the kinases. Comparisons between P/CDK2/cyclinA/roscovitine and CDK5/p25/roscovitine complexes reveal differences in the conformations of the glutamine around the active sites, which may be exploited to find highly selective inhibitors with respect to CDK2 and CDK5.

Study of the inhibition of cyclin-dependent kinases with roscovitine and indirubin-3'-oxime from molecular dynamics simulations

Molecular dynamics simulations were performed to elucidate the interactions of CDK2 and CDK5 complexes with three inhibitors: R-roscovitine, S-roscovitine, and indirubin-3'-oxime. The preference of the two complexes for R-roscovitine over the S enantiomer, as reported by the experiment, was also found by the simulations. More importantly, the simulations showed that the cause of the stronger affinity for the R enantiomer is the presence of an important hydrogen bond between R-roscovitine and the kinases not found with S-roscovitine. The simulations also showed two amino acid mutations in the active site of CDK5/R-roscovitine that favor binding-enhanced electrostatic contributions, making the inhibitor more effective for CDK5 than for CDK2. This suggests that the effectiveness of roscovitine-like inhibitors can be improved by enhancing their electrostatic interaction with the kinases. Finally, molecular mechanics-Possion-Boltzmann/surface area calculations of the CDK5/indirubin-3'-oxime system in both water-excluded and water-included environments gave significantly different electrostatic contributions to the binding. The simulations detected the displacement of a water molecule in the active site of the water-included CDK/indirubin-3'-oxime system. This resulted in a more conserved binding pattern than the water-excluded structure. Hence, in the design of new indirubin-like inhibitors, it is important to include the water molecule in the analysis.

Different Mechanisms of CDK5 and CDK2 Activation as Revealed by CDK5/p25 and CDK2/Cyclin A Dynamics

A detailed analysis is presented of the dynamics of human CDK5 in complexes with the protein activator p25 and the purine-like inhibitor roscovitine. These and other findings related to the activation of CDK5 are critically reviewed from a molecular perspective. In addition, the results obtained on the behavior of CDK5 are compared with data on CDK2 to assess the differences and similarities between the two kinases in terms of (i) roscovitine binding, (ii) regulatory subunit association, (iii) conformational changes in the T-loop following CDK/regulatory subunit complex formation, and (iv) specificity in CDK/regulatory subunit recognition. An energy decomposition analysis, used for these purposes, revealed why the binding of p25 alone is sufficient to stabilize the extended active T-loop conformation of CDK5, whereas the equivalent conformational change in CDK2 requires both the binding of cyclin A and phosphorylation of the Thr(160) residue. The interaction energy of the CDK5 T-loop with p25 is about 26 kcal.mol(-1) greater than that of the CDK2 T-loop with cyclin A. The binding pattern between CDK5 and p25 was compared with that of CDK2/cyclin A to find specific regions involved in CDK/regulatory subunit recognition. The analyses performed revealed that the alphaNT-helix of cyclin A interacts with the alpha6-alpha7 loop and the alpha7 helix of CDK2, but these regions do not interact in the CDK5/p25 complex. Further differences between the CDK5/p25 and CDK2/cyclin A systems studied are discussed with respect to their specific functionality.