Explaining the inhibition of cyclin-dependent kinase 5 by peptides derived from p25 with molecular dynamics simulations and MM-PBSA

In silico design of potential Mcl-1 peptide-based inhibitors

CONTEXT

BIM (Bcl-2 interacting mediator of apoptosis)-derived peptides that specifically target over-expressed Mcl-1 (myeloid cell leukemia-1) protein and induce apoptosis are potentially anti-cancer agents. Since the helicity of BIM-derived peptides has a crucial role in their functionality, a range of strategies have been used to increase the helicity including the introduction of unnatural residues and stapling methods that have some drawbacks such as the accumulation in the liver. To avoid these drawbacks, this study aimed to design a more helical peptide by utilizing bioinformatics algorithms and molecular dynamics simulations without exploiting unnatural residues and stapling methods. MM-PBSA results showed that the mutations of A4fE and A2eE in analogue 5 demonstrate a preference towards binding with Mcl-1. As evidenced by Circular dichroism results, the helicity increases from 18 to 34%, these findings could enhance the potential of analogue 5 as an anti-cancer agent targeting Mcl-1. The applied strategies in this research could shed light on the in silico peptide design. Moreover, analogue 5 as a drug candidate can be evaluated in vitro and in vivo studies.

METHODS

The sequence of the lead peptide was determined using the ApInAPDB database and PRALINE program. Contact finder and PDBsum web server softwares were used to determine the contact involved amino acids in complex with Mcl-1. All identified salt bridge contributing residues were unaltered to preserve the binding affinity. After proposing novel analogues, their secondary structures were predicted by Cham finder web server software and GOR, Neural Network, and Chou-Fasman algorithms. Finally, molecular dynamics simulations run for 100 ns were done using the GROMACS, version 5.0.7, with the CHARMM36 force field. MM-PBSA was used to assess binding affinity specificity in targeting Mcl-1 and Bcl-xL (B-cell lymphoma extra-large).

Effect of lipid saturation on amyloid-beta peptide partitioning and aggregation in neuronal membranes: molecular dynamics simulations

Aggregation of amyloid-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β (Aβ) peptides, cleaved from the amyloid precursor protein, is known as a precursor of the Alzheimer’s disease (AD). It is also known that Alzheimer’s disease is characterized by a substantial decrease of the amount of polyunsaturated lipids in the neuronal membranes of the frontal gray matter. To get insight into possible interconnection of these phenomena, we have carried out molecular dynamics simulations of two fragments of A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β peptide, A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{1-28}$$\end{document}1-28 and A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{26-40}$$\end{document}26-40, in four different lipid bilayers: two monocomponent ones (14:0-14:0 PC, 18:0-22:6 PC), and two bilayers containing mixtures of 18:0-18:0 PE, 22:6-22:6 PE, 16:0-16:0 PC and 18:1-18:1 PC lipids of composition mimicking neuronal membranes in a “healthy” and “AD” brain. The simulations showed that the presence of lipids with highly unsaturated 22:6cis fatty acids chains strongly affects the interaction of amyloid-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β peptides with lipid membranes. The polyunsaturated lipids cause stronger adsorption of A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β-peptides by the membrane and lead to weaker binding between peptides when the latter form aggregates. This difference in the behaviour observed in monocomponent bilayers is propagated in a similar fashion to the mixed membranes mimicking composition of neuronal membranes in “healthy” and “AD” brains, with “healthy” membrane having higher fraction of polyunsaturated lipids. Our simulations give strong indication that it can be physical–chemical background of the interconnection between amyloid fibrillization causing Alzheimer’s disease, and content of polyunsaturated lipids in the neuronal membranes.

The interaction of Munc 18 (p67) with the p10 domain of p35 protects in vivo Cdk5/p35 activity from inhibition by TFP5, a peptide derived from p35

In a series of studies, we have identified TFP5, a truncated fragment of p35, the Cdk5 kinase regulatory protein, which inhibits Cdk5/p35 and the hyperactive Cdk5/p25 activities in test tube experiments. In cortical neurons, however, and in vivo in Alzheimer's disease (AD) model mice, the peptide specifically inhibits the Cdk5/p25 complex and not the endogenous Cdk5/p35. To account for the selective inhibition of Cdk5/p25 activity, we propose that the "p10" N-terminal domain of p35, absent in p25, spares Cdk5/p35 because p10 binds to macromolecules (e.g., tubulin and actin) as a membrane-bound multimeric complex that favors p35 binding to Cdk5 and catalysis. To test this hypothesis, we focused on Munc 18, a key synapse-associated neuronal protein, one of many proteins copurifying with Cdk5/p35 in membrane-bound multimeric complexes. Here we show that, in vitro, the addition of p67 protects Cdk5/p35 and has no effect on Cdk5/p25 activity in the presence of TFP5. In cortical neurons transfected with p67siRNA, we also show that TFP5 inhibits Cdk5/p35 activity, whereas in the presence of p67 the activity is protected. It does so without affecting any other kinases of the Cdk family of cyclin kinases. This difference may be of significant therapeutic value because the accumulation of the deregulated, hyperactive Cdk5/p25 complex in human brains has been implicated in pathology of AD and other neurodegenerative disorders.

Cervical cancer knowledge and screening behaviors among female university graduates of year 2012 attending national graduate orientation program, Bhutan

BACKGROUND

Cervical cancer is the leading female cancer in Bhutan. This study describes the level of cervical cancer knowledge and screening behaviors among female university graduates attending the National Graduate Orientation Program (NGOP), 2012.

METHODS

A cross-sectional study of female graduates attending NGOP was conducted using self-administered anonymous questionnaire developed through literature reviews and expert discussions to elicit information on demographic characteristics, knowledge, screening behaviors and determinants of cervical cancer. The association of demographic and other important study characteristics with uptake of Pap test was investigated using cross tabulation and Fischer Exact test. Frequencies and percentages were calculated for all the questions.

RESULTS

The average age of the participants was 23.43 ± SD 2.73. About 92% (n = 513) of the respondents were aged 25 years or less and 7.9% (n = 44) of the respondents were aged 26 or more. The study revealed low cervical cancer knowledge and poor screening behavior among the graduates. The mean knowledge score was 3.571 (SD1.75, Range 0-8). About 6% (n=34) of the respondents reported undergoing Pap test at least once and 94% reported as never having done Pap test. The most commonly cited reasons for not doing Pap test included "never thought I needed one" (57%, n = 320), "embarrassment of being examined by male health professional" and "fear of finding out cancer". The study revealed evidence of significant association between increasing age, those who are married, knowledge score and those recommended for screening by health professionals with the uptake of Pap test.

CONCLUSION

Our study revealed poor knowledge and screening behaviors among female university graduates in Bhutan. This may be suggestive of even poorer awareness and screening practices among young unmarried women who are less educated or with no education. Although our study group is not appropriate for measuring practice of cervical cancer screening in the country, the findings are expected to highlight the shortcomings and trigger development of comprehensive cervical cancer control programs in Bhutan.

Molecular dynamic simulations reveal the mechanism of binding between xanthine inhibitors and DPP-4

We apply molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation to investigate and reveal the binding mechanism between five xanthine inhibitors and DPP-4. The electrostatic and van der Waals interactions of the five inhibitors with DPP-4 are analyzed and discussed. The computed binding free energies using MM-PBSA method are in qualitatively agreement with experimental inhibitory potency of five inhibitors. The hydrogen bonds of inhibitors with Ser630 and Asp663 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the key contacts with His740, Asn710, Trp629, and Tyr666 have larger contributions to the binding free energy and play important roles in distinguishing the variant bioactivity of five inhibitors.

Molecular dynamic simulations give insight into the mechanism of binding between 2-aminothiazole inhibitors and CDK5

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.

Cyclin-dependent kinases 5 template: useful for virtual screening

The present study reports the development of a template for the active binding site of Cdk5 for structure-based drug design. The developed template of Cdk5 was validated by redocking with ligands I (PBD code 1UNG), II (PBD code 1UNL) and III (PBD code 1UNH). The results demonstrate a good match of the docked and the crystallographic binding orientations with RMSD less than 2.0Å. The validation results show that the constructed Cdk5 template is a good model system for predicting ligand binding orientations and binding affinities. Furthermore, the developed template was applied to predict binding mode and binding affinity of thirty-six known Cdk5 inhibitors. The results showed that the binding energy of almost Cdk5 inhibitors related to their biological evaluation.

First BRET-based screening assay performed in budding yeast leads to the discovery of CDK5/p25 interaction inhibitors

The protein kinase CDK5 (cyclin-dependent kinase 5) is activated through its association with a cyclin-like protein p35 or p39. In pathological conditions (such as Alzheimer's disease and various other neuropathies), truncation of p35 leads to the appearance of the p25 protein. The interaction of p25 with CDK5 up-regulates the kinase activity and modifies the substrate specificity. ATP-mimetic inhibitors of CDK5 have already been developed. However, the lack of selectivity of such inhibitors is often a matter of concern. An alternative approach can be used to identify highly specific inhibitors that disrupt protein interactions involving protein kinases. We have developed a bioluminescence resonance energy transfer (BRET)-based screening assay in yeast to discover protein-protein interaction inhibitors (P2I2). Here, we present the first use of BRET in yeast for the screening of small molecule libraries. This screening campaign led to the discovery of one molecule that prevents the interaction between CDK5 and p25, thus inhibiting the protein kinase activity. This molecule may give rise to high-specificity drug candidates.

Structural and dynamic determinants of ligand binding and regulation of cyclin-dependent kinase 5 by pathological activator p25 and inhibitory peptide CIP

The crystal structure of the cdk5/p25 complex has provided information on possible molecular mechanisms of the ligand binding, specificity, and regulation of the kinase. Comparative molecular dynamics simulations are reported here for physiological conditions. This study provides new insight on the mechanisms that modulate such processes, which may be exploited to control pathological activation by p25. The structural changes observed in the kinase are stabilized by a network of interactions involving highly conserved residues within the cyclin-dependent kinase (cdk) family. Collective motions of the proteins (cdk5, p25, and CIP) and their complexes are identified by principal component analysis, revealing two conformational states of the activation loop upon p25 complexation, which are absent in the uncomplexed kinase and not apparent from the crystal. Simulations of the uncomplexed inhibitor CIP show structural rearrangements and increased flexibility of the interfacial loop containing the critical residue E240, which becomes fully hydrated and available for interactions with one of several positively charged residues in the kinase. These changes provide a rationale for the observed high affinity and enhanced inhibitory action of CIP when compared to either p25 or the physiological activators of cdk5.