Conformational Equilibrium of CDK/Cyclin Complexes by Molecular Dynamics with Excited Normal Modes

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.

Slower CDK4 and faster CDK2 activation in the cell cycle

Dysregulation of cyclin-dependent kinases (CDKs) impacts cell proliferation, driving cancer. Here, we ask why the cyclin-D/CDK4 complex governs cell cycle progression through the longer G1 phase, whereas cyclin-E/CDK2 regulates the shorter G1/S phase transition. We consider available experimental cellular and structural data including cyclin-E's high-level burst, sustained duration of elevated cyclin-D expression, and explicit solvent molecular dynamics simulations of the inactive monomeric and complexed states, to establish the conformational tendencies along the landscape of the distinct activation scenarios of cyclin-D/CDK4 and cyclin-E/CDK2 in the G1 phase and G1/S transition of the cell cycle, respectively. These lead us to propose slower activation of cyclin-D/CDK4 and rapid activation of cyclin-E/CDK2. We provide the mechanisms through which this occurs, offering innovative CDK4 drug design considerations. Our insightful mechanistic work addresses a compelling cell cycle regulation question and illuminates the distinct activation speeds between the G1 and the G1/S phases, which are crucial for function.

Cell cycle progression mechanisms: slower cyclin-D/CDK4 activation and faster cyclin-E/CDK2

Dysregulation of cyclin-dependent kinases (CDKs) impacts cell proliferation, driving cancer. Here, we ask why the cyclin-D/CDK4 complex governs cell cycle progression through the longer G1 phase, whereas cyclin-E/CDK2 regulates the short G1/S phase transition. We consider the experimentally established high-level bursting of cyclin-E, and sustained duration of elevated cyclin-D expression in the cell, available experimental cellular and structural data, and comprehensive explicit solvent molecular dynamics simulations to provide the mechanistic foundation of the distinct activation scenarios of cyclin-D/CDK4 and cyclin-E/CDK2 in the G1 phase and G1/S transition of the cell cycle, respectively. These lead us to propose slower activation of cyclin-D/CDK4 and rapid activation of cyclin-E/CDK2. Importantly, we determine the mechanisms through which this occurs, offering innovative CDK4 drug design considerations. Our insightful mechanistic work addresses the compelling cell cycle regulation question and illuminates the distinct activation speeds in the G1 versus G1/S phases, which are crucial for cell function.

MDexciteR: Enhanced Sampling Molecular Dynamics by Excited Normal Modes or Principal Components Obtained from Experiments

Molecular dynamics with excited normal modes (MDeNM) is an enhanced sampling method for exploring conformational changes in proteins with minimal biases. The excitation corresponds to injecting kinetic energy along normal modes describing intrinsic collective motions. Herein, we developed a new automated open-source implementation, MDexciteR (https://github.com/mcosta27/MDexciteR), enabling the integration of MDeNM with two commonly used simulation programs with GPU support. Second, we generalized the method to include the excitation of principal components calculated from experimental ensembles. Finally, we evaluated whether the use of coarse-grained normal modes calculated with elastic network representations preserved the performance and accuracy of the method. The advantages and limitations of these new approaches are discussed based on results obtained for three different protein test cases: two globular and a protein/membrane system.

ABCG2/BCRP transport mechanism revealed through kinetically excited targeted molecular dynamics simulations

ABCG2/BCRP is an ABC transporter that plays an important role in tissue protection by exporting endogenous substrates and xenobiotics. ABCG2 is of major interest due to its involvement in multidrug resistance (MDR), and understanding its complex efflux mechanism is essential to preventing MDR and drug-drug interactions (DDI). ABCG2 export is characterized by two major conformational transitions between inward- and outward-facing states, the structures of which have been resolved. Yet, the entire transport cycle has not been characterized to date. Our study bridges the gap between the two extreme conformations by studying connecting pathways. We developed an innovative approach to enhance molecular dynamics simulations, ‘kinetically excited targeted molecular dynamics’, and successfully simulated the transitions between inward- and outward-facing states in both directions and the transport of the endogenous substrate estrone 3-sulfate. We discovered an additional pocket between the two substrate-binding cavities and found that the presence of the substrate in the first cavity is essential to couple the movements between the nucleotide-binding and transmembrane domains. Our study shed new light on the complex efflux mechanism, and we provided transition pathways that can help to identify novel substrates and inhibitors of ABCG2 and probe new drug candidates for MDR and DDI.

A comprehensive review on acridone based derivatives as future anti-cancer agents and their structure activity relationships

The development of drug resistance and severe side-effects has reduced the clinical efficacy of the existing anti-cancer drugs available in the market. Thus, there is always a constant need to develop newer anti-cancer drugs with minimal adverse effects. Researchers all over the world have been focusing on various alternative strategies to discover novel, potent, and target specific molecules for cancer therapy. In this direction, several heterocyclic compounds are being explored but amongst them one promising heterocycle is acridone which has attracted the attention of medicinal chemists and gained huge biological importance as acridones are found to act on different therapeutically proven molecular targets, overcome ABC transporters mediated drug resistance and DNA intercalation in cancer cells. Some of these acridone derivatives have reached clinical studies as these heterocycles have shown huge potential in cancer therapeutics and imaging. Here, the authors have attempted to compile and make some recommendations of acridone based derivatives concerning their cancer biological targets and in vitro-cytotoxicity based on drug design and novelty to increase their therapeutic potential. This review also provides some important insights on the design, receptor targeting and future directions for the development of acridones as possible clinically effective anti-cancer agents.

Biotin-painted proteins have thermodynamic stability switched by kinetic folding routes

Biotin-labeled proteins are widely used as tools to study protein-protein interactions and proximity in living cells. Proteomic methods broadly employ proximity-labeling technologies based on protein biotinylation in order to investigate the transient encounters of biomolecules in subcellular compartments. Biotinylation is a post-translation modification in which the biotin molecule is attached to lysine or tyrosine residues. So far, biotin-based technologies proved to be effective instruments as affinity and proximity tags. However, the influence of biotinylation on aspects such as folding, binding, mobility, thermodynamic stability, and kinetics needs to be investigated. Here, we selected two proteins [biotin carboxyl carrier protein (BCCP) and FKBP3] to test the influence of biotinylation on thermodynamic and kinetic properties. Apo (without biotin) and holo (biotinylated) protein structures were used separately to generate all-atom structure-based model simulations in a wide range of temperatures. Holo BCCP contains one biotinylation site, and FKBP3 was modeled with up to 23 biotinylated lysines. The two proteins had their estimated thermodynamic stability changed by altering their energy landscape. In all cases, after comparison between the apo and holo simulations, differences were observed on the free-energy profiles and folding routes. Energetic barriers were altered with the density of states clearly showing changes in the transition state. This study suggests that analysis of large-scale datasets of biotinylation-based proximity experiments might consider possible alterations in thermostability and folding mechanisms imposed by the attached biotins.

Cell cycle involvement in cancer therapy; WEE1 kinase, a potential target as therapeutic strategy

Mitosis is the process of cell division and is regulated by checkpoints in the cell cycle. G1-S, S, and G2-M are the three main checkpoints that prevent initiation of the next phase of the cell cycle phase until previous phase has completed. DNA damage leads to activation of the G2-M checkpoint, which can trigger a downstream DNA damage response (DDR) pathway to induce cell cycle arrest while the damage is repaired. If the DNA damage cannot be repaired, the replication stress response (RSR) pathway finally leads to cell death by apoptosis, in this case called mitotic catastrophe. Many cancer treatments (chemotherapy and radiotherapy) cause DNA damages based on SSBs (single strand breaks) or DSBs (double strand breaks), which cause cell death through mitotic catastrophe. However, damaged cells can activate WEE1 kinase (as a part of the DDR and RSR pathways), which prevents apoptosis and cell death by inducing cell cycle arrest at G2 phase. Therefore, inhibition of WEE1 kinase could sensitize cancer cells to chemotherapeutic drugs. This review focuses on the role of WEE1 kinase (as a biological macromolecule which has a molecular mass of 96 kDa) in the cell cycle, and its interactions with other regulatory pathways. In addition, we discuss the potential of WEE1 inhibition as a new therapeutic approach in the treatment of various cancers, such as melanoma, breast cancer, pancreatic cancer, cervical cancer, etc.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Cyclin-dependent kinase 4/6 inhibitors for cancer therapy: a patent review (2015 - 2019)

INTRODUCTION

Cyclin-dependent kinases 4 and 6 (CDK4/6) along with their upstream/downstream components are pivotal regulators for the cell cycle progression. The dysfunction of CDK4/6 is the common feature and promoting factor in various cancer types. In-depth research on CDK4/6 inhibitors has afforded therapeutic agents, while new challenges and ideas are emerging concomitantly.

AREAS COVERED

This review focuses on patent publications related to CDK4/6 inhibitors which could be utilized for anti-cancer purposes during the period 2015-2019.

EXPERT OPINION

The increasingly comprehensive and thorough understanding of CDK4/6 inhibitors facilitates them to break through the current limitations. Hence the utilization of CDK4/6 inhibitors for cancer therapy in the near future is likely to be performed in diverse forms and for distinct purposes. Selectivity over kinases is still crucial to new agent development but shall be prudently dealt with. The gradually revealing of resistance and adverse events proposed another issue that calls for new tackling strategies.

Effects of pH and aggregation in the human prion conversion into scrapie form: a study using molecular dynamics with excited normal modes

There are two different prion conformations: (1) the cellular natural (PrP^C) and (2) the scrapie (PrP^Sc), an infectious form that tends to aggregate under specific conditions. PrP^C and PrP^Sc are widely different regarding secondary and tertiary structures. PrP^Sc contains more and longer β-strands compared to PrP^C. The lack of solved PrP^Sc structures precludes a proper understanding of the mechanisms related to the transition between cellular and scrapie forms, as well as the aggregation process. In order to investigate the conformational transition between PrP^C and PrP^Sc, we applied MDeNM (molecular dynamics with excited normal modes), an enhanced sampling simulation technique that has been recently developed to probe large structural changes. These simulations yielded new structural rearrangements of the cellular prion that would have been difficult to obtain with standard MD simulations. We observed an increase in β-sheet formation under low pH (≤ 4) and upon oligomerization, whose relevance was discussed on the basis of the energy landscape theory for protein folding. The characterization of intermediate structures corresponding to transition states allowed us to propose a conversion model from the cellular to the scrapie prion, which possibly ignites the fibril formation. This model can assist the design of new drugs to prevent neurological disorders related to the prion aggregation mechanism.

Computational engineering of cellulase Cel9A-68 functional motions through mutations in its linker region

Cellulase collective motions design through linker mutations leads to the enhancement of protein flexibility and function.

PLA2-like proteins myotoxic mechanism: a dynamic model description

Phospholipase A₂-like (PLA₂-like) proteins contribute to the development of muscle necrosis in Viperidae snake bites and are not efficiently neutralized by current antivenom treatments. The toxic mechanisms of PLA₂-like proteins are devoid of catalytic activity and not yet fully understood, although structural and functional experiments suggest a dimeric assembly and that the C-terminal residues are essential to myotoxicity. Herein, we characterized the functional mechanism of bothropic PLA₂-like structures related to global and local measurements using the available models in the Protein Data Bank and normal mode molecular dynamics (NM-MD). Those measurements include: (i) new geometric descriptions between their monomers, based on Euler angles; (ii) characterizations of canonical and non-canonical conformations of the C-terminal residues; (iii) accessibility of the hydrophobic channel; (iv) inspection of ligands; and (v) distance of clustered residues to toxin interface of interaction. Thus, we described the allosteric activation of PLA₂-like proteins and hypothesized that the natural movement between monomers, calculated from NM-MD, is related to their membrane disruption mechanism, which is important for future studies of the inhibition process. These methods and strategies can be applied to other proteins to help understand their mechanisms of action.

Frustration, function and folding

Natural protein molecules are exceptional polymers. Encoded in apparently random strings of amino-acids, these objects perform clear physical tasks that are rare to find by simple chance. Accurate folding, specific binding, powerful catalysis, are examples of basic chemical activities that the great majority of polypeptides do not display, and are thought to be the outcome of the natural history of proteins. Function, a concept genuine to Biology, is at the core of evolution and often conflicts with the physical constraints. Locating the frustration between discrepant goals in a recurrent system leads to fundamental insights about the chances and necessities that shape the encoding of biological information.

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.

Recent advances on CDK inhibitors: An insight by means of in silico methods

The cyclin dependent kinases (CDKs) are a small family of serine/threonine protein kinases that can act as a potential therapeutic target in several proliferative diseases, including cancer. This short review is a survey on the more recent research progresses in the field achieved by using in silico methods. All the "armamentarium" available to the medicinal chemists (docking protocols and molecular dynamics, fragment-based, de novo design, virtual screening, and QSAR) has been employed to the discovery of new, potent, and selective inhibitors of cyclin dependent kinases. The results cited herein can be useful to understand the nature of the inhibitor-target interactions, and furnish an insight on the structural/molecular requirements necessary to achieve the required selectivity against cyclin dependent kinases over other types of kinases.

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.

Bending-Twisting Motions and Main Interactions in Nucleoplasmin Nuclear Import

Alpha solenoid proteins play a key role in regulating the classical nuclear import pathway, recognizing a target protein and transporting it into the nucleus. Importin-α (Impα) is the solenoid responsible for cargo protein recognition, and it has been extensively studied by X-ray crystallography to understand the binding specificity. To comprehend the main motions of Impα and to extend the information about the critical interactions during carrier-cargo recognition, we surveyed different conformational states based on molecular dynamics (MD) and normal mode (NM) analyses. Our model of study was a crystallographic structure of Impα complexed with the classical nuclear localization sequence (cNLS) from nucleoplasmin (Npl), which was submitted to multiple 100 ns of MD simulations. Representative conformations were selected for calculating the 87 lowest frequencies NMs of vibration, and a displacement approach was applied along each NM. Based on geometric criteria, using the radius of curvature and inter-repeat angles as the reference metrics, the main motions of Impα were described. Moreover, we determined the salt bridges, hydrogen bonds and hydrophobic interactions in the Impα-NplNLS interface. Our results show the bending and twisting motions participating in the recognition of nuclear proteins, allowing the accommodation and adjustment of a classical bipartite NLS sequence. The essential contacts for the nuclear import were also described and were mostly in agreement with previous studies, suggesting that the residues in the cNLS linker region establish important contacts with Impα adjusting the cNLS backbone. The MD simulations combined with NM analysis can be applied to the Impα-NLS system to help understand interactions between Impα and cNLSs and the analysis of non-classic NLSs.

TP-2Rho Is a Sensitive Solvatochromic Red-Shifted Probe for Monitoring the Interactions between CDK4 and Cyclin D

Understanding the intricate steps of protein kinase regulation requires characterization of protein-protein interactions between the catalytic subunit, its regulatory partners and the substrate. Fluorescent probes are useful tools with which to study such interactions and to gain insight into their affinities and specificities. Solvatochromic probes, which display changes in their fluorescence emission in response to changes in the polarity of the medium, are particularly attractive. Here we describe conjugation of a switchable fluorescent dye, TP-2Rho, to peptide and protein derivatives of cyclin-dependent kinase 4 (CDK4) and its application to characterization of the interactions between the catalytic subunit of this kinase, its regulatory partner cyclin D1 and a peptide substrate. We demonstrate the sensitivity of TP-2Rho in relation to of those other dyes used for monitoring peptide-protein and protein-protein interactions. Moreover, we show that TP-Rho-labelled peptides can be introduced into living cells to probe endogenous CDK4/cyclin D.