Activation and inhibition of cyclin-dependent kinase-2 by phosphorylation; a molecular dynamics study reveals the functional importance of the glycine-rich loop

Exploring the mechanism of action of spirooxindoles as a class of CDK2 inhibitors: a structure-based computational approach

Cyclin-dependent kinase 2 (CDK2) regulates cell cycle checkpoints in the synthesis and mitosis phases and plays a pivotal role in cancerous cell proliferation. The activation of CDK2, influenced by various protein signaling pathways, initiates the phosphorylation process. Due to its crucial role in carcinogenesis, CDK2 is a druggable hotspot target to suppress cancer cell proliferation. In this context, several studies have identified spirooxindoles as an effective class of CDK2 inhibitors. In the present study, three spirooxindoles (SOI1, SOI2, and SOI3) were studied to understand their inhibitory mechanism against CDK2 through a structure-based approach. Molecular docking and molecular dynamics (MD) simulations were performed to explore their interactions with CDK2 at the molecular level. The calculated binding free energy for the spirooxindole-based CDK2 inhibitors aligned well with experimental results regarding CDK2 inhibition. Energy decomposition (ED) analysis identified key binding residues, including I10, G11, T14, R36, F82, K89, L134, P155, T158, Y159, and T160, in the CDK2 active site and T-loop phosphorylation. Molecular mechanics (MM) energy was identified as the primary contributor to stabilizing inhibitor binding in the CDK2 protein structure. Furthermore, the analysis of binding affinity revealed that the inhibitor SOI1 binds more strongly to CDK2 compared to the other inhibitors under investigation. It demonstrated a robust interaction with the crucial residue T160 in the T-loop phosphorylation site, responsible for kinase activation. These insights into the inhibitory mechanism are anticipated to contribute to the development of potential CDK2 inhibitors using the spirooxindole scaffold.

EGR1 and EGR2 positively regulate plant ABA signaling by modulating the phosphorylation of SnRK2.2

During abscisic acid (ABA) signaling, reversible phosphorylation controls the activity and accumulation of class III SNF1-RELATED PROTEIN KINASE 2s (SnRK2s). While protein phosphatases that negatively regulate SnRK2s have been identified, those that positively regulate ABA signaling through SnRK2s are less understood. In this study, Arabidopsis thaliana mutants of Clade E Growth-Regulating 1 and 2 (EGR1/2), which belong to the protein phosphatase 2C family, exhibited reduced ABA sensitivity in terms of seed germination, cotyledon greening, and ABI5 accumulation. Conversely, overexpression increased these ABA-induced responses. Transcriptomic data revealed that most ABA-regulated genes in egr1 egr2 plants were expressed at reduced levels compared with those in Col-0 after ABA treatment. Abscisic acid up-regulated EGR1/2, which interact directly with SnRK2.2 through its C-terminal domain I. Genetic analysis demonstrated that EGR1/2 function through SnRK2.2 during ABA response. Furthermore, SnRK2.2 de-phosphorylation by EGR1/2 was identified at serine 31 within the ATP-binding pocket. A phospho-mimic mutation confirmed that phosphorylation at serine 31 inhibited SnRK2.2 activity and reduced ABA responsiveness in plants. Our findings highlight the positive role of EGR1/2 in regulating ABA signaling, they reveal a new mechanism for modulating SnRK2.2 activity, and provide novel insight into how plants fine-tune their responses to ABA.

Rapid adaptation to CDK2 inhibition exposes intrinsic cell-cycle plasticity

CDK2 is a core cell-cycle kinase that phosphorylates many substrates to drive progression through the cell cycle. CDK2 is hyperactivated in multiple cancers and is therefore an attractive therapeutic target. Here, we use several CDK2 inhibitors in clinical development to interrogate CDK2 substrate phosphorylation, cell-cycle progression, and drug adaptation in preclinical models. Whereas CDK1 is known to compensate for loss of CDK2 in Cdk2-/- mice, this is not true of acute inhibition of CDK2. Upon CDK2 inhibition, cells exhibit a rapid loss of substrate phosphorylation that rebounds within several hours. CDK4/6 activity backstops inhibition of CDK2 and sustains the proliferative program by maintaining Rb1 hyperphosphorylation, active E2F transcription, and cyclin A2 expression, enabling re-activation of CDK2 in the presence of drug. Our results augment our understanding of CDK plasticity and indicate that co-inhibition of CDK2 and CDK4/6 may be required to suppress adaptation to CDK2 inhibitors currently under clinical assessment.

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.

Prediction of allosteric druggable pockets of cyclin-dependent kinases

Cyclin-dependent kinase (Cdk) proteins play crucial roles in the cell cycle progression and are thus attractive drug targets for therapy against such aberrant cell cycle processes as cancer. Since most of the available Cdk inhibitors target the highly conserved catalytic ATP pocket and their lack of specificity often lead to side effects, it is imperative to identify and characterize less conserved non-catalytic pockets capable of interfering with the kinase activity allosterically. However, a systematic analysis of these allosteric druggable pockets is still in its infancy. Here, we summarize the existing Cdk pockets and their selectivity. Then, we outline a network-based pocket prediction approach (NetPocket) and illustrate its utility for systematically identifying the allosteric druggable pockets with case studies. Finally, we discuss potential future directions and their challenges.

Influenza A virus uses PSMA2 for downregulation of NRF2-mediated oxidative stress response

Influenza A virus (IAV), an obligatory intracellular parasite, uses host cellular molecules to complete its replication cycle and suppress immune responses. Proteasome subunit alpha type 2 (PSMA2) is a cellular protein highly expressed in IAV-infected human lung epithelial A549 cells. PSMA2 is part of the 20S proteasome complex that degrades or recycles defective proteins and involves proteolytic modification of many cellular regulatory proteins. However, the role of PSMA2 in IAV replication is not well understood. In this study, PSMA2 knockdown (KD) in A549 cells caused a significant reduction in extracellular progeny IAV, but intracellular viral protein translation and viral RNA transcription were not affected. This indicates that PSMA2 is a critical host factor for IAV maturation. To better understand the interplay between PSMA2 KD and IAV infection at the proteomic level, we used the SomaScan 1.3K version, which measures 1,307 proteins to analyze alterations induced by these treatments. We found seven cellular signaling pathways, including phospholipase C signaling, Pak signaling, and nuclear factor erythroid 2p45-related factor 2 (NRF2)-mediated oxidative stress response signaling, that were inhibited by IAV infection but significantly activated by PSMA2 KD. Further analysis of NRF2-mediated oxidative stress response signaling indicated IAV inhibits accumulation of reactive oxygen species (ROS), but ROS levels significantly increased during IAV infection in PSMA2 KD cells. However, IAV infection caused significantly higher NFR2 nuclear translocation that was inhibited in PSMA2 KD cells. This indicates that PSMA2 is required for NRF2-mediated ROS neutralization and that IAV uses PSMA2 to escape viral clearance via the NRF2-mediated cellular oxidative response.

IMPORTANCE Influenza A virus (IAV) remains one of the most significant infectious agents, responsible for 3 million to 5 million illnesses each year and more than 50 million deaths during the 20th century. The cellular processes that promote and inhibit IAV infection and pathogenesis remain only partially understood. PSMA2 is a critical component of the 20S proteasome and ubiquitin-proteasome system, which is important in the replication of numerous viruses. This study examined host protein responses to IAV infection alone, PSMA2 knockdown alone, and IAV infection in the presence of PSMA2 knockdown and determined that interfering with PSMA2 function affected IAV maturation. These results help us better understand the importance of PSMA2 in IAV replication and may pave the way for designing additional IAV antivirals targeting PSMA2 or the host proteasome for the treatment of seasonal flu.

Targeted Protein Profiling of In Vivo NIPP-Treated Tissues Using DigiWest Technology

Non-invasive physical plasma (NIPP) is a novel therapeutic tool, currently being evaluated for the treatment of cancer and precancerous lesions in gynecology and other disciplines. Additionally, patients with cervical intraepithelial neoplasia (CIN) may benefit from NIPP treatment due to its non-invasive, side-effect-free, and tissue-sparing character. However, the molecular impact of in vivo NIPP treatment needs to be further investigated. For this purpose, usually only very small tissue biopsies are available after NIPP treatment. Here, we adapted DigiWest technology, a high-throughput bead-based Western blot, for the analysis of formalin-fixed paraffin-embedded (FFPE) cervical punch biopsies with a minimal sample amount. We investigated the molecular effects of NIPP treatment directly after (0 h) and 24 h after in vivo application. Results were compared to in vitro NIPP-treated human malignant cervical cells. NIPP effects were primarily based on an inhibitory impact on the cell cycle and cell growth factors. DigiWest technology was suitable for detailed protein profiling of small, primary FFPE biopsies.

Pharmacological relevance of CDK inhibitors in Alzheimer's disease

Evidence suggests that cell cycle activation plays a role in the pathophysiology of neurodegenerative diseases. Alzheimer's disease is a progressive, terminal neurodegenerative disease that affects memory and other important mental functions. Intracellular deposition of Tau protein, a hyperphosphorylated form of a microtubule-associated protein, and extracellular aggregation of Amyloid β protein, which manifests as neurofibrillary tangles (NFT) and senile plaques, respectively, characterize this condition. In recent years, however, several studies have concluded that cell cycle re-entry is one of the key causes of neuronal death in the pathogenesis of Alzheimer's disease. The eukaryotic cell cycle is well-coordinated machinery that performs critical functions in cell replenishment, such as DNA replication, cell creation, repair, and the birth of new daughter cells from the mother cell. The complex interplay between the levels of various cyclins and cyclin-dependent kinases (CDKs) at different checkpoints is needed for cell cycle synchronization. CDKIs (cyclin-dependent kinase inhibitors) prevent cyclin degradation and CDK inactivation. Different external and internal factors regulate them differently, and they have different tissue expression and developmental functions. The checkpoints ensure that the previous step is completed correctly before starting the new cell cycle phase, and they protect against the transfer of defects to the daughter cells. Due to the development of more selective and potent ATP-competitive CDK inhibitors, CDK inhibitors appear to be on the verge of having a clinical impact. This avenue is likely to yield new and effective medicines for the treatment of cancer and other neurodegenerative diseases. These new methods for recognizing CDK inhibitors may be used to create non-ATP-competitive agents that target CDK4, CDK5, and other CDKs that have been recognized as important therapeutic targets in Alzheimer's disease treatment.

The role of cyclin-dependent kinase 5 in neuropathic pain

The chronification of pain can be attributed to changes in membrane receptors and channels underlying neuronal plasticity and signal transduction largely within nociceptive neurons that initiate and maintain pathological pain states. These proteins are subject to dynamic modification by posttranslational modifications, creating a code that controls protein function in time and space. Phosphorylation is an important posttranslational modification that affects ∼30% of proteins in vivo. Increased phosphorylation of various nociceptive ion channels and of their modulators underlies sensitization of different pain states. Cyclin-dependent kinases are proline-directed serine/threonine kinases that impact various biological and cellular systems. Cyclin-dependent kinase 5 (Cdk5), one member of this kinase family, and its activators p35 and p39 are expressed in spinal nerves, dorsal root ganglia, and the dorsal horn of the spinal cord. In neuropathic pain conditions, expression and/or activity of Cdk5 is increased, implicating Cdk5 in nociception. Experimental evidence suggests that Cdk5 is regulated through its own phosphorylation, through increasing p35's interaction with Cdk5, and through cleavage of p35 into p25. This narrative review discusses the molecular mechanisms of Cdk5-mediated regulation of target proteins involved in neuropathic pain. We focus on Cdk5 substrates that have been linked to nociceptive pathways, including channels (eg, transient receptor potential cation channel and voltage-gated calcium channel), proteins involved in neurotransmitter release (eg, synaptophysin and collapsin response mediator protein 2), and receptors (eg, glutamate, purinergic, and opioid). By altering the phosphoregulatory "set point" of proteins involved in pain signaling, Cdk5 thus appears to be an attractive target for treating neuropathic pain conditions.

KinaseMD: kinase mutations and drug response database

Mutations in kinases are abundant and critical to study signaling pathways and regulatory roles in human disease, especially in cancer. Somatic mutations in kinase genes can affect drug treatment, both sensitivity and resistance, to clinically used kinase inhibitors. Here, we present a newly constructed database, KinaseMD (kinase mutations and drug response), to structurally and functionally annotate kinase mutations. KinaseMD integrates 679 374 somatic mutations, 251 522 network-rewiring events, and 390 460 drug response records curated from various sources for 547 kinases. We uniquely annotate the mutations and kinase inhibitor response in four types of protein substructures (gatekeeper, A-loop, G-loop and αC-helix) that are linked to kinase inhibitor resistance in literature. In addition, we annotate functional mutations that may rewire kinase regulatory network and report four phosphorylation signals (gain, loss, up-regulation and down-regulation). Overall, KinaseMD provides the most updated information on mutations, unique annotations of drug response especially drug resistance and functional sites of kinases. KinaseMD is accessible at https://bioinfo.uth.edu/kmd/, having functions for searching, browsing and downloading data. To our knowledge, there has been no systematic annotation of these structural mutations linking to kinase inhibitor response. In summary, KinaseMD is a centralized database for kinase mutations and drug response.

Identification of key interactions between SARS-CoV-2 main protease and inhibitor drug candidates

The number of cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (COVID-19) has reached over 114,000. SARS-CoV-2 caused a pandemic in Wuhan, China, in December 2019 and is rapidly spreading globally. It has been reported that peptide-like anti-HIV-1 drugs are effective against SARS-CoV Main protease (Mpro). Due to the close phylogenetic relationship between SARS-CoV and SARS-CoV-2, their main proteases share many structural and functional features. Thus, these drugs are also regarded as potential drug candidates targeting SARS-CoV-2 Mpro. However, the mechanism of action of SARS-CoV-2 Mpro at the atomic-level is unknown. In the present study, we revealed key interactions between SARS-CoV-2 Mpro and three drug candidates by performing pharmacophore modeling and 1 μs molecular dynamics (MD) simulations. His41, Gly143, and Glu166 formed interactions with the functional groups that were common among peptide-like inhibitors in all MD simulations. These interactions are important targets for potential drugs against SARS-CoV-2 Mpro.

Landscape of drug-resistance mutations in kinase regulatory hotspots

More than 48 kinase inhibitors (KIs) have been approved by Food and Drug Administration. However, drug-resistance (DR) eventually occurs, and secondary mutations have been found in the previously targeted primary-mutated cancer cells. Cancer and drug research communities recognize the importance of the kinase domain (KD) mutations for kinasopathies. So far, a systematic investigation of kinase mutations on DR hotspots has not been done yet. In this study, we systematically investigated four types of representative mutation hotspots (gatekeeper, G-loop, αC-helix and A-loop) associated with DR in 538 human protein kinases using large-scale cancer data sets (TCGA, ICGC, COSMIC and GDSC). Our results revealed 358 kinases harboring 3318 mutations that covered 702 drug resistance hotspot residues. Among them, 197 kinases had multiple genetic variants on each residue. We further computationally assessed and validated the epidermal growth factor receptor mutations on protein structure and drug-binding efficacy. This is the first study to provide a landscape view of DR-associated mutation hotspots in kinase's secondary structures, and its knowledge will help the development of effective next-generation KIs for better precision medicine.

Molecular Dynamics Simulation reveals the mechanism by which the Influenza Cap-dependent Endonuclease acquires resistance against Baloxavir marboxil

Baloxavir marboxil (BXM), an antiviral drug for influenza virus, inhibits RNA replication by binding to RNA replication cap-dependent endonuclease (CEN) of influenza A and B viruses. Although this drug was only approved by the FDA in October 2018, drug resistant viruses have already been detected from clinical trials owing to an I38 mutation of CEN. To investigate the reduction of drug sensitivity by the I38 mutant variants, we performed a molecular dynamics (MD) simulation on the CEN-BXM complex structure to analyze variations in the mode of interaction. Our simulation results suggest that the side chain methyl group of I38 in CEN engages in a CH-pi interaction with the aromatic ring of BXM. This interaction is abolished in various I38 mutant variants. Moreover, MD simulation on various mutation models and binding free energy prediction by MM/GBSA method suggest that the I38 mutation precludes any interaction with the aromatic ring of BXA and thereby reduces BXA sensitivity.

Insights into evolutionary interaction patterns of the 'Phosphorylation Activation Segment' in kinase

We are interested in studying the phosphorylation of the kinase activation loop, distinguishing the passage from the unphosphorylated to the phosphorylated form without allostery. We performed an interaction study to trace the change of interactions between the activation segment and the kinase catalytic core, before and after phosphorylation. Results show that the structural changes are mainly due to the attraction between the phosphate group and guanidine groups of the arginine side chains of RD-pocket, which are constituted mainly of guanidine groups of the catalytic loop, the β9, and the αC helix. This attraction causes propagation of structural variation of the activation segment, principally towards the N-terminal. The structural variations are not made on all the amino acids of the activation segment; they are conditioned by the existence of two beta sheets stabilizing the loop during phosphorylation. The first,β6-β9 sheet is usually present in most of the kinases; the second, β10-β11 is formed due to the interaction between the main chain amino acids of the activation loop and the αEF/αF loop.

Goniolanceolatins A-H, Cytotoxic Bis-styryllactones from Goniothalamus lanceolatus

Eight new bis-styryllactones, goniolanceolatins A-H (1-8), possessing a rare α,β-unsaturated δ-lactone moiety with a (6S)-configuration, were isolated from the CH₂Cl₂ extract of the stembark and roots of Goniothalamus lanceolatus Miq., a plant endemic to Malaysia. Absolute structures were established through extensive 1D- and 2D-NMR data analysis, in combination with electronic dichroism (ECD) data. All of the isolates were evaluated for their cytotoxicity against human lung and colorectal cancer cell lines. Compounds 2 and 4 showed cytotoxicity, with IC₅₀ values ranging from 2.3 to 4.2 μM, and were inactive toward human noncancerous lung and colorectal cells. Compounds 1, 3, 6, 7, and 8 showed moderate to weak cytotoxicity. Docking studies of compounds 2 and 4 showed that they bind with EGFR tyrosine kinase and cyclin-dependent kinase 2 through hydrogen bonding interactions with the important amino acids, including Lys721, Met769, Asn818, Arg157, Ile10, and Glu12.

SnRK1 Kinase and the NAC Transcription Factor SOG1 Are Components of a Novel Signaling Pathway Mediating the Low Energy Response Triggered by ATP Depletion

Plant growth is strictly controlled by cell division, elongation, and differentiation for which adequate supplies of intracellular ATP are required. However, it is unclear how changes in the amount of intracellular ATP affect cell division and growth. To reveal the specific pathway dependent on ATP concentration, we performed analyses on the Arabidopsis mitochondria mutation sd3. The mutant is tiny, a result of a low amount of ATP caused by the disruption of Tim21, a subunit of the TIM23 protein complex localized in the inner membrane of the mitochondria. Loss of function of suppressor of gamma response 1 (SOG1) also restored the dwarf phenotype of wild type treated with antimycin A, a blocker of ATP synthesis in mitochondria. The sd3 phenotype is partially restored by the introduction of sog1, suppressor of gamma response 1, and kin10/kin11, subunits of Snf1-related kinase 1 (SnRK1). Additionally, SOG1 interacted with SnRK1, and was modified by phosphorylation in planta only after treatment with antimycin A. Transcripts of several negative regulators of the endocycle were up-regulated in the sd3 mutant, and this high expression was not observed in sd3sog1 and sd3kin11. We suggest that there is a novel regulatory mechanism for the control of plant cell cycle involving SnRK1 and SOG1, which is induced by low amounts of intracellular ATP, and controls plant development.

Targeting p53 as a promising therapeutic option for cancer by re-activating the wt or mutant p53’s tumor suppression

p53 protein, a product of the TP53 tumor suppressor gene, controls the cellular genome’s integrity and is an important regulator of cell cycling, proliferation, apoptosis and metabolism. Mutations of TP53 or inactivation of its gene product are among the first events initiating malignant transformation. The consequent loss of control over the cell cycle, resulting in accelerated cell proliferation and facilitating metabolic reprogramming, gives the initiated (premalignant) cells numerous advantages over healthy cells. Interestingly, p53 status is not only an important marker in cancer diagnosis; it has also become a promising target of personalized therapy. Depending on the TP53 status different therapeutic options have been developed. (Re)-activation of p53 functionality in cancer cells offers promising new alternatives to existing oncological therapies.

Exploring the selectivity of inhibitor complexes with Bcl-2 and Bcl-XL: A molecular dynamics simulation approach

B-cell lymphoma 2 (Bcl-2) family proteins are potential drug targets in cancer and have a relatively flat and flexible binding site. ABT-199 is one of the most promising selective Bcl-2 inhibitors, and A-1155463 selectively inhibits Bcl-XL. Although the amino acid sequences of the binding sites of these two inhibitors are similar, the inhibitors selectively bind the target protein. In order to determine the origin of the selectivity of these inhibitors, we conducted molecular dynamics simulations using protein-inhibitor modeling. We confirmed that ASP103 of Bcl-2 is a key residue and that hydrogen bonding between ASP103 and ABT-199 confers the Bcl-2 selectivity of this inhibitor. For Bcl-XL selectivity, the secondary structure of α-helix 3 is a key factor. PHE105, SER106, and LEU108 in the loose α-helix 3 interact with A-1155463 to confer Bcl-XL selectivity. These findings provide important insights into the molecular mechanisms of selective inhibitors of Bcl-2 family proteins.

Reduced state transition barrier of CDK6 from open to closed state induced by Thr177 phosphorylation and its implication in binding modes of inhibitors

BACKGROUND

CDK6 is considered as a highly validated anticancer drug target due to its essential role in regulating cell cycle progression at G1 restriction point. Activation of CDK6 requires the phosphorylation of Thr177 on A-loop, but the structural insights of the activation mechanism remain unclear.

METHODS

Herein, all-atoms molecular dynamics (MD) simulations were used to study the effects of Thr177 phosphorylation on the dynamic structure of CDK6-Vcyclin complex.

RESULTS

MD results indicated that the free energy barrier of the transition from open to closed state decreased ~47.2% after Thr177 phosphorylation. Key steps along the state transition process were obtained from a cluster analysis. Binding preference of ten different inhibitors to open or closed state were also investigated through molecular docking along with MD simulations methods.

CONCLUSIONS

Our results indicated that Thr177 phosphorylation increased the flexibility around the ATP-binding pocket. The transition of the ATP-binding pocket between open and closed states should be considered for understanding the binding of CDK6 inhibitors.

GENERAL SIGNIFICANCE

This work could deepen the understanding of CDKs activation mechanism, and provide useful information for the discovery of new CDKs inhibitors with high affinity and specificity.

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

Variability of the Cyclin-Dependent Kinase 2 Flexibility Without Significant Change in the Initial Conformation of the Protein or Its Environment; a Computational Study

BACKGROUND

Protein flexibility, which has been referred as a dynamic behavior has various roles in proteins' functions. Furthermore, for some developed tools in bioinformatics, such as protein-protein docking software, considering the protein flexibility, causes a higher degree of accuracy. Through undertaking the present work, we have accomplished the quantification plus analysis of the variations in the human Cyclin Dependent Kinase 2 (hCDK2) protein flexibility without affecting a significant change in its initial environment or the protein per se.

OBJECTIVES

The main goal of the present research was to calculate variations in the flexibility for each residue of the hCDK2, analysis of their flexibility variations through clustering, and to investigate the functional aspects of the residues with high flexibility variations.

MATERIALS AND METHODS

Using Gromacs package (version 4.5.4), three independent molecular dynamics (MD) simulations of the hCDK2 protein (PDB ID: 1HCL) was accomplished with no significant changes in their initial environments, structures, or conformations, followed by Root Mean Square Fluctuations (RMSF) calculation of these MD trajectories. The amount of variations in these three curves of RMSF was calculated using two formulas.

RESULTS

More than 50% of the variation in the flexibility (the distance between the maximum and the minimum amount of the RMSF) was found at the region of Val-154. As well, there are other major flexibility fluctuations in other residues. These residues were mostly positioned in the vicinity of the functional residues. The subsequent works were done, as followed by clustering all hCDK2 residues into four groups considering the amount of their variability with respect to flexibility and their position in the RMSF curves.

CONCLUSIONS

This work has introduced a new class of flexibility aspect of the proteins' residues. It could also help designing and engineering proteins, with introducing a new dynamic aspect of hCDK2, and accordingly, for the other similar globular proteins. In addition, it could provide a better computational calculation of the protein flexibility, which is, especially important in the comparative studies of the proteins' flexibility.

Cyclin-dependent kinase 2 protects podocytes from apoptosis

Loss of podocytes is an early feature of diabetic nephropathy (DN) and predicts its progression. We found that treatment of podocytes with sera from normoalbuminuric type 1 diabetes patients with high lipopolysaccharide (LPS) activity, known to predict progression of DN, downregulated CDK2 (cyclin-dependent kinase 2). LPS-treatment of mice also reduced CDK2 expression. LPS-induced downregulation of CDK2 was prevented in vitro and in vivo by inhibiting the Toll-like receptor (TLR) pathway using immunomodulatory agent GIT27. We also observed that CDK2 is downregulated in the glomeruli of obese Zucker rats before the onset of proteinuria. Knockdown of CDK2, or inhibiting its activity with roscovitine in podocytes increased apoptosis. CDK2 knockdown also reduced expression of PDK1, an activator of the cell survival kinase Akt, and reduced Akt phosphorylation. This suggests that CDK2 regulates the activity of the cell survival pathway via PDK1. Furthermore, PDK1 knockdown reduced the expression of CDK2 suggesting a regulatory loop between CDK2 and PDK1. Collectively, our data show that CDK2 protects podocytes from apoptosis and that reduced expression of CDK2 associates with the development of DN. Preventing downregulation of CDK2 by blocking the TLR pathway with GIT27 may provide a means to prevent podocyte apoptosis and progression of DN.

Contribution of EGFR and ErbB-3 Heterodimerization to the EGFR Mutation-Induced Gefitinib- and Erlotinib-Resistance in Non-Small-Cell Lung Carcinoma Treatments

EGFR mutation-induced drug resistance has become a major threat to the treatment of non-small-cell lung carcinoma. Essentially, the resistance mechanism involves modifications of the intracellular signaling pathways. In our work, we separately investigated the EGFR and ErbB-3 heterodimerization, regarded as the origin of intracellular signaling pathways. On one hand, we combined the molecular interaction in EGFR heterodimerization with that between the EGFR tyrosine kinase and its inhibitor. For 168 clinical subjects, we characterized their corresponding EGFR mutations using molecular interactions, with three potential dimerization partners (ErbB-2, IGF-1R and c-Met) of EGFR and two of its small molecule inhibitors (gefitinib and erlotinib). Based on molecular dynamics simulations and structural analysis, we modeled these mutant-partner or mutant-inhibitor interactions using binding free energy and its components. As a consequence, the mutant-partner interactions are amplified for mutants L858R and L858R_T790M, compared to the wild type EGFR. Mutant delL747_P753insS represents the largest difference between the mutant-IGF-1R interaction and the mutant-inhibitor interaction, which explains the shorter progression-free survival of an inhibitor to this mutant type. Besides, feature sets including different energy components were constructed, and efficient regression trees were applied to map these features to the progression-free survival of an inhibitor. On the other hand, we comparably examined the interactions between ErbB-3 and its partners (EGFR mutants, IGF-1R, ErbB-2 and c-Met). Compared to others, c-Met shows a remarkably-strong binding with ErbB-3, implying its significant role in regulating ErbB-3 signaling. Moreover, EGFR mutants corresponding to poor clinical outcomes, such as L858R_T790M, possess lower binding affinities with ErbB-3 than c-Met does. This may promote the communication between ErbB-3 and c-Met in these cancer cells. The analysis verified the important contribution of IGF-1R or c-Met in the drug resistance mechanism developed in lung cancer treatments, which may bring many benefits to specialized therapy design and innovative drug discovery.

Break CDK2/Cyclin E1 interface allosterically with small peptides

Most inhibitors of Cyclin-dependent kinase 2 (CDK2) target its ATP-binding pocket. It is difficult, however, to use this pocket to design very specific inhibitors because this catalytic pocket is highly conserved in the protein family of CDKs. Here we report some short peptides targeting a noncatalytic pocket near the interface of the CDK2/Cyclin complex. Docking and molecular dynamics simulations were used to select the peptides, and detailed dynamical network analysis revealed that these peptides weaken the complex formation via allosteric interactions. Our experiments showed that upon binding to the noncatalytic pocket, these peptides break the CDK2/Cyclin complex partially and diminish its kinase activity in vitro. The binding affinity of these peptides measured by Surface Plasmon Resonance can reach as low as 0.5 µM.

Multiple steps to activate FAK's kinase domain: adaptation to confined environments?

Protein kinases regulate cell signaling by phosphorylating their substrates in response to environment-specific stimuli. Using molecular dynamics, we studied the catalytically active and inactive conformations of the kinase domain of the focal adhesion kinase (FAK), which are distinguished by displaying a structured or unstructured activation loop, respectively. Upon removal of an ATP analog, we show that the nucleotide-binding pocket in the catalytically active conformation is structurally unstable and fluctuates between an open and closed configuration. In contrast, the pocket remains open in the catalytically inactive form upon removal of an inhibitor from the pocket. Because temporal pocket closures will slow the ATP on-rate, these simulations suggest a multistep process in which the kinase domain is more likely to bind ATP in the catalytically inactive than in the active form. Transient closures of the ATP-binding pocket might allow FAK to slow down its catalytic cycle. These short cat naps could be adaptions to crowded or confined environments by giving the substrate sufficient time to diffuse away. The simulations show further how either the phosphorylation of the activation loop or the activating mutations of the so-called SuperFAK influence the electrostatic switch that controls kinase activity.

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.

CDC25A targeting by miR-483-3p decreases CCND-CDK4/6 assembly and contributes to cell cycle arrest

Disruption of contact inhibition and serum afflux that occur after a tissue injury activate cell cycle, which then stops when confluence is reached again. Although the events involved in cell cycle entry have been widely documented, those managing cell cycle exit have remained so far ill defined. We have identified that the final stage of wound closure is preceded in keratinocytes by a strong accumulation of miR-483-3p, which acts as a mandatory signal triggering cell cycle arrest when confluence is reached. Blocking miR-483-3p accumulation strongly delays cell cycle exit, maintains cells into a proliferative state and retards their differentiation program. Using two models of cell cycle synchronization (i.e. mechanical injury and serum addition), we show that an ectopic upregulation of miR-483-3p blocks cell cycle progression in early G1 phase. This arrest results from a direct targeting of the CDC25A phosphatase by miR-483-3p, which can be impeded using an anti-miRNA against miR-483-3p or a protector that blocks the complex formation between miR-483-3p and the 3'-untranslated region (UTR) of CDC25A transcript. We show that the miRNA-induced silencing of CDC25A increases the tyrosine phosphorylation status of CDK4/6 cyclin-dependent kinases which, in turn, abolishes CDK4/6 capacity to associate with D-type cyclins. This prevents CDK4/6 kinases' activation, impairs downstream events such as cyclin E stimulation and sequesters cells in early G1. We propose this new regulatory process of cyclin-CDK association as a general mechanism coupling miRNA-mediated CDC25A invalidation to CDK post-transcriptional modifications and cell cycle control.

Localization and function of the Ska complex during mouse oocyte meiotic maturation

The Ska (spindle and kinetochore-associated) complex is composed of three proteins: Ska1, Ska2 and Ska3. It is required for stabilizing kinetochore-microtubule (KT-MT) interactions and silencing spindle checkpoint during mitosis. However, its roles in meiosis remain unclear. The present study was designed to investigate the localization and function of the Ska complex during mouse oocyte meiotic maturation. Our results showed that the localization and function of Ska complex in mouse oocyte meiosis differ in part from those in mitosis. Injection of low dose exogenous Myc-Ska mRNA showed that, instead of localizing to the kinetochores (KTs) and mediating KT-MT interactions from pro-metaphase to mid-anaphase stages as in mitosis, the members of the Ska complex were only localized on spindle microtubules from the Pro-MI to MII stages in mouse oocyte meiosis. Time-lapse live imaging analysis showed that knockdown of any member of the Ska complex by Morpholino injection into mouse oocytes resulted in spindle movement defects and enlarged polar bodies. Depletion of the whole Ska complex disrupted the stability of the anaphase spindle and influenced the extrusion of the first polar body. Taken together, these results show that the Ska complex plays an important role in meiotic spindle migration and anaphase spindle stability during mouse oocyte maturation.

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.

Recombinant Leishmania mexicana CRK3:CYCA has protein kinase activity in the absence of phosphorylation on the T-loop residue Thr178

The activity of cyclin-dependent kinases (CDKs), which are key regulators of the eukaryotic cell cycle, is regulated through post-translational mechanisms, including binding of a cyclin and phosphorylation. Previously studies have shown that Leishmania mexicana CRK3 is an essential CDK that is a functional homologue of human CDK1. In this study, recombinant histidine tagged L. mexicana CRK3 and the cyclin CYCA were combined in vitro to produce an active histone H1 kinase that was inhibited by the CDK inhibitors, flavopiridol and indirubin-3'-monoxime. Protein kinase activity was observed in the absence of phosphorylation of the T-loop residue Thr178, but increased 5-fold upon phosphorylation by the CDK activating kinase Civ1 of Saccharomyces cerevisiae. Seven recombinant L. major CRKs (1, 2, 3, 4, 6, 7 and 8) were also expressed and purified, none of which were active as monomers. Moreover, only CRK3 was phosphorylated by Civ1. HA-tagged CYCA expressed in L. major procyclic promastigotes was co-precipitated with CRK3 and exhibited histone H1 kinase activity. These data indicate that in Leishmania CYCA interacts with CRK3 to form an active protein kinase, confirm the conservation of the regulatory mechanisms that control CDK activity in other eukaryotes, but identifies biochemical differences to human CDK1.

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.

Why pyridine containing pyrido[2,3-d]pyrimidin-7-ones selectively inhibit CDK4 than CDK2: insights from molecular dynamics simulation

Designing selective cyclin-dependent kinase 4 (CDK4) inhibitors is an area of intense research to develop potential anticancer drugs. The molecular basis governing the selective inhibition of CDK4 by lig17 (6-bromo-8-cyclopentyl-2-(5-piperazin-1-yl-pyridin-2-ylamino)-8H-pyrido[2,3-d]pyrimidin-7-one) has been investigated using molecular dynamics simulation. The positive charge on the ligand was determined to be an important contributor for CDK4 selectivity due to the electronegative nature of its active site. Similar studies on CDK2 indicated that Lys89 intrudes into the active site displacing the positive charge on lig17 away from the active center. This intrusion was observed to propel a drastic conformational change in lig17, weakening its binding interactions with the protein. The pyridine nitrogen (N(AR)) of lig17 was capable of interacting with His95 (CDK4) through hydrogen bonding. N(AR) also showed a strong tendency to mediate protein-ligand interactions through a bridged water molecule, only when bound to CDK4. The G-loop of CDK4 was observed to fluctuate extensively when complexed with lig17 and a novel "flipping-out" mechanism exhibited by Tyr17(CDK4/CDK4-17) is reported in this study. Although these proteins have similar folds, the results from principal component analysis (PCA) indicate that CDK4 and CDK2 follow an anti-correlated behavior towards the accessibility of the active site.

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.

Molecular cloning, sequence characterization and tissue transcription profile analyses of two novel genes: LCK and CDK2 from the Black-boned sheep (Ovis aries)

The complete coding sequences of two sheep genes--LCK and CDK2--were amplified using the rapid amplification of cDNA ends method based on three sheep EST sequences whose translated amino acids contain the domain PTKc_Lck_BIk and S_TKc domain, respectively. The sequence analyses of these two genes revealed that the sheep LCK gene encodes a protein of 509 amino acids which has high homology with the lymphocyte-specific protein tyrosine kinase (LCK) of eight species: bovine (99%), human (96%), dog (96%), Aotus nancymaae (95%), mouse (94%), rat (91%), horse (91%) and chicken (81%). The sheep CDK2 gene encodes a protein of 298 amino acids which has high homology with the cyclin-dependent kinase 2 (CDK2) of ten species: bovine (100%), goat (100%), rat (99%), mouse (99%), Chinese hamster (99%), dog (98%), golden hamster (98%), human (98%), horse (98%) and rhesus monkey (98%). The tissue transcription profile analyses indicated that that the Black-boned sheep LCK and CDK2 genes are generally but differentially expressed in the detected tissues including in tissues including spleen, muscle, skin, kidney, lung, liver and heart. These data serve as a foundation for further insight into these two genes.