Search

Development of multi-target drugs is becoming increasingly attractive in the repertoire of protein kinase inhibitors discovery. In this study, we carried out molecular docking, molecular dynamics simulations, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy calculations, principal component analysis (PCA), and dynamical cross-correlation matrices (DCCM) to dissect the molecular mechanism for the valmerin-19 acting as a dual inhibitor for glycogen synthase kinase 3β (GSK3β) and cyclin-dependent kinase 5 (CDK5). Detailed MM-PBSA calculations revealed that the binding free energies of the valmerin-19 to GSK3β/CDK5 were calculated to be -12.60 ± 2.28 kcal mol(-1) and -11.85 ± 2.54 kcal mol(-1), respectively, indicating that valmerin-19 has the potential to act as a dual inhibitor of GSK3β/CDK5. The analyses of PCA and DCCM results unraveled that binding of the valmerin-19 reduced the conformational dynamics of GSK3β/CDK5 and the valmerin-19 bound to GSK3β/CDK5 might occur mostly through a conformational selection mechanism. This study may be helpful for the future design of novel and potent dual GSK3β/CDK5 inhibitors.

Molecular docking, molecular dynamics (MD) simulations, and binding free energy analysis were performed to reveal differences in the binding affinities between five 2-aminothiazole inhibitors and CDK5. The hydrogen bonding and hydrophobic interactions between inhibitors and adjacent residues are analyzed and discussed. The rank of calculated binding free energies using the MM-PBSA method is consistent with experimental result. The results illustrate that hydrogen bonds with Cys83 favor inhibitor binding. The van der Waals interactions, especially the important contact with Ile10, dominate in the binding free energy and play a crucial role in distinguishing the different bioactivity of the five inhibitors.

Structure-based 3D-QSAR approaches (CoMFA and CoMSIA) were applied to understand the structural requirements of the Cyclin-dependent kinase 5/p25 inhibitors. Cyclin-dependent kinase 5 (CDK5) is believed to play an important role in the development of the central nervous system during the process of mammalian embryogenesis. Genetic algorithm based docking program (GOLD) was successfully utilized to orient the compounds inside the binding pocket of the CDK5/p25 structure. The adapted alignment method with the suitable parameters resulted in a reliable model. Furthermore, the final model was robust enough to forecast the activities of test compounds, satisfactorily. The contour maps were produced around the functional groups to understand the SAR requirements. Moreover, we also investigate the structural attributes of the inhibitors which make them selective toward CDK5/p25 over its close counterpart, i.e., CDK2. The study could be helpful to rationalize the new compounds with better inhibition and selectivity profiles against CDK5/p25.

Glycogen synthase kinase 3 (GSK-3) is an attractive target for the treatment of diabetes, and paullones have been reported to be effective inhibitors of GSK-3. However, it is still a challenging task to improve selectivity among protein kinases, especially cyclin-dependent kinases (CDKs). Here we investigated the mechanism that enables paullones to selectively inhibit GSK-3 rather than cyclin-dependent kinase 5 (CDK5) using sequence alignment, molecular dynamics simulations, free-energy calculations and free-energy decomposition analysis. The results indicate that the interaction between paullones and Val135 of GSK-3 is obviously stronger than that between paullones and Cys83 of CDK5, suggesting that paullones could be utilized as potent selective inhibitors. Meanwhile, we observed that the decrease in the interaction between paullones and the Asp86 of CDK5 favors their selectivity towards GSK-3 rather than CDK5, as demonstrated using 1-azakenpaullone as an example. Although substitution at position 9 and replacement at position 2 may influence the activity of GSK-3, they only have a minor effect on the selectivity. We expect that the information obtained here could prove useful for developing specific paullone inhibitors of GSK-3.

In the current work, CDK5/p25 complexes were pulled apart by applying external forces with steered molecular dynamics (SMD) simulations. The crucial interactions between the kinase and the activation protein were investigated and the SMD simulations showed that several activation-relevant motifs of CDK5 leave p25 in sequence during the pulling and lead to an apo-CDK2 like CDK5 structure after separation. Based on systematic examination of hydrogen bond breaking and classical MD/molecular mechanics-generalized Born/surface area) (MM-GBSA) calculations, a CDK5 activation mechanism by p25 is suggested. This is the first step towards the systemic development of CDK inhibitors and the mechanism proposed could lead to a better understanding of the protein-protein recognition characteristics between the kinase and its activator.

A cyclin-dependent kinase (CDK) 5 inhibitory peptide (CIP) from p25 was recently reported to inhibit CDK5/p25 activity in vitro but had no effect on endogenous cdc2 kinase activity. This may lead to a specific CDK5 inhibition strategy in the treatment of neurodegeneration. However, the mechanism of the inhibition remains unclear. In this work, molecular dynamics simulations and energy decomposition calculation models were set up to investigate the deregulation mechanisms of CIP on CDK5 activity. The results show that truncation of the N, and C terminals of p25 introduces important conformational changes into a hydrophobic pocket that is crucial for accommodating Ile153 on the activation loop of CDK5. In addition, such truncations lead to distortion and displacement of the activation loop and consequently affect binding of the substrate peptide. New inhibition sites for selectively inhibiting the activity of CDK5 are also suggested.