Structural analysis of the inhibition of Cdk4 and Cdk6 by p16(INK4a) through molecular dynamics simulations

Clinical and molecular overview of immunotherapeutic approaches for malignant skin melanoma: Past, present and future

Traditional therapeutic approaches for malignant melanoma, have proved to be limited and/or ineffective, especially with respect to their role in improving patient survival and tumor recurrence. In this regard, immunotherapy has been demonstrated to be a promising therapeutic alternative, boosting antitumor responses through the modulation of cell signaling pathways involved in the effector mechanisms of the immune system, particularly, the so-called "immunological checkpoints". Clinical studies on the efficacy and safety of immunotherapeutic regimens, alone or in combination with other antitumor approaches, have increased dramatically in recent decades, with very encouraging results. Hence, this review will discuss the current immunotherapeutic regimens used to treat malignant melanoma, as well as the molecular and cellular mechanisms involved. In addition, current clinical studies that have investigated the use, efficacy, and adverse events of immunotherapy in melanoma will also be discussed.

A Review of Key Biological and Molecular Events Underpinning Transformation of Melanocytes to Primary and Metastatic Melanoma

Melanoma is a major public health concern that is responsible for significant morbidity and mortality, particularly in countries such as New Zealand and Australia where it is the commonest cause of cancer death in young adults. Until recently, there were no effective drug therapies for patients with advanced melanoma however significant advances in our understanding of the biological and molecular basis of melanoma in recent decades have led to the development of revolutionary treatments, including targeted molecular therapy and immunotherapy. This review summarizes our current understanding of the key events in the pathway of melanomagenesis and discusses the role of genomic analysis as a potential tool for improved diagnostic evaluation, prognostication and treatment strategies. Ultimately, it is hoped that a continued deeper understanding of the mechanisms of melanomagenesis will lead to the development of even more effective treatments that continue to provide better outcomes for patients with melanoma.

Cellular senescence and senescence-associated secretory phenotype: comparison of idiopathic pulmonary fibrosis, connective tissue disease-associated interstitial lung disease, and chronic obstructive pulmonary disease

Background

The senescence-associated secretory phenotype (SASP) develops due to cellular senescence during conditions such as chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF). However, studies comparing the degree of cellular senescence and SASP between COPD and IPF are limited. Furthermore, to the best of our knowledge, no study has examined cellular senescence and/or SASP in connective tissue disease-associated interstitial lung disease (CTD-ILD).

Methods

To compare the degree of cellular senescence among COPD, IPF, and CTD-ILD, tissue samples from surgical lung biopsies or noncancerous tissue from lobectomy specimens of patients with lung cancer were subjected to immunostaining for p16 and p21. Double-staining for p16 and phosphorylated NF-κB was performed to verify the relationship between cellular senescence and SASP.

Results

There was a greater degree of enhancement of p16 and p21 expression in patients with IPF than in those with COPD and controls. Immunostaining for p16 revealed an enhanced expression of this marker in patients with COPD compared with that in controls. No significant differences were observed in the phosphorylated NF-κB expression rate of p16-positive and p16-negative cells among patients with IPF, CTD-ILD, and COPD.

Conclusions

Epithelial cells in patients with IPF express higher levels of both cellular senescence and SASP than those in patients with COPD or controls.

New insights into the structural dynamics of the kinase JNK3

In this work, we study the dynamics and the energetics of the all-atom structure of a neuronal-specific serine/threonine kinase c-Jun N-terminal kinase 3 (JNK3) in three states: unphosphorylated, phosphorylated, and ATP-bound phosphorylated. A series of 2 µs atomistic simulations followed by a conformational landscape mapping and a principal component analysis supports the mechanistic understanding of the JNK3 inactivation/activation process and also indicates key structural intermediates. Our analysis reveals that the unphosphorylated JNK3 undergoes the ‘open-to-closed’ movement via a two-step mechanism. Furthermore, the phosphorylation and ATP-binding allow the JNK3 kinase to attain a fully active conformation. JNK3 is a widely studied target for small-drugs used to treat a variety of neurological disorders. We believe that the mechanistic understanding of the large-conformational changes upon the activation of JNK3 will aid the development of novel targeted therapeutics.

Gene regulatory network underlying the immortalization of epithelial cells

BACKGROUND

Tumorigenic transformation of human epithelial cells in vitro has been described experimentally as the potential result of spontaneous immortalization. This process is characterized by a series of cell-state transitions, in which normal epithelial cells acquire first a senescent state which is later surpassed to attain a mesenchymal stem-like phenotype with a potentially tumorigenic behavior. In this paper we aim to provide a system-level mechanistic explanation to the emergence of these cell types, and to the time-ordered transition patterns that are common to neoplasias of epithelial origin. To this end, we first integrate published functional and well-curated molecular data of the components and interactions that have been found to be involved in such cell states and transitions into a network of 41 molecular components. We then reduce this initial network by removing simple mediators (i.e., linear pathways), and formalize the resulting regulatory core into logical rules that govern the dynamics of each of the network components as a function of the states of its regulators.

RESULTS

Computational dynamic analysis shows that our proposed Gene Regulatory Network model recovers exactly three attractors, each of them defined by a specific gene expression profile that corresponds to the epithelial, senescent, and mesenchymal stem-like cellular phenotypes, respectively. We show that although a mesenchymal stem-like state can be attained even under unperturbed physiological conditions, the likelihood of converging to this state is increased when pro-inflammatory conditions are simulated, providing a systems-level mechanistic explanation for the carcinogenic role of chronic inflammatory conditions observed in the clinic. We also found that the regulatory core yields an epigenetic landscape that restricts temporal patterns of progression between the steady states, such that recovered patterns resemble the time-ordered transitions observed during the spontaneous immortalization of epithelial cells, both in vivo and in vitro.

CONCLUSION

Our study strongly suggests that the in vitro tumorigenic transformation of epithelial cells, which strongly correlates with the patterns observed during the pathological progression of epithelial carcinogenesis in vivo, emerges from underlying regulatory networks involved in epithelial trans-differentiation during development.

Prospecting for novel plant-derived molecules of Rauvolfia serpentina as inhibitors of Aldose Reductase, a potent drug target for diabetes and its complications

Aldose Reductase (AR) is implicated in the development of secondary complications of diabetes, providing an interesting target for therapeutic intervention. Extracts of Rauvolfia serpentina, a medicinal plant endemic to the Himalayan mountain range, have been known to be effective in alleviating diabetes and its complications. In this study, we aim to prospect for novel plant-derived inhibitors from R. serpentina and to understand structural basis of their interactions. An extensive library of R. serpentina molecules was compiled and computationally screened for inhibitory action against AR. The stability of complexes, with docked leads, was verified using molecular dynamics simulations. Two structurally distinct plant-derived leads were identified as inhibitors: indobine and indobinine. Further, using these two leads as templates, 16 more leads were identified through ligand-based screening of their structural analogs, from a small molecules database. Thus, we obtained plant-derived indole alkaloids, and their structural analogs, as potential AR inhibitors from a manually curated dataset of R. serpentina molecules. Indole alkaloids reported herein, as a novel structural class unreported hitherto, may provide better insights for designing potential AR inhibitors with improved efficacy and fewer side effects.

Cyclin-dependent kinases 4 and 6 control tumor progression and direct glucose oxidation in the pentose cycle

Cyclin-dependent kinases CDK4 and CDK6 are essential for the control of the cell cycle through the G(1) phase. Aberrant expression of CDK4 and CDK6 is a hallmark of cancer, which would suggest that CDK4 and CDK6 are attractive targets for cancer therapy. Herein, we report that calcein AM (the calcein acetoxymethyl-ester) is a potent specific inhibitor of CDK4 and CDK6 in HCT116 human colon adenocarcinoma cells, inhibiting retinoblastoma protein (pRb) phosphorylation and inducing cell cycle arrest in the G(1) phase. The metabolic effects of calcein AM on HCT116 cells were also evaluated and the flux between the oxidative and non-oxidative branches of the pentose phosphate pathway was significantly altered. To elucidate whether these metabolic changes were due to the inhibition of CDK4 and CDK6, we also characterized the metabolic profile of a CDK4, CDK6 and CDK2 triple knockout of mouse embryonic fibroblasts. The results show that the metabolic profile associated with the depletion of CDK4, CDK6 and CDK2 coincides with the metabolic changes induced by calcein AM on HCT116 cells, thus confirming that the inhibition of CDK4 and CDK6 disrupts the balance between the oxidative and non-oxidative branches of the pentose phosphate pathway. Taken together, these results indicate that low doses of calcein can halt cell division and kill tumor cells. Thus, selective inhibition of CDK4 and CDK6 may be of greater pharmacological interest, since inhibitors of these kinases affect both cell cycle progression and the robust metabolic profile of tumors.

Investigation of the differences in activity between hydroxycycloalkyl N1 substituted pyrazole derivatives as inhibitors of B-Raf kinase by using docking, molecular dynamics, QM/MM, and fragment-based de novo design: study of binding mode of diastereomer compounds

N1 substituted pyrazole derivatives show diverse B-Raf kinase inhibitory activities when different hydroxy-substituted cycloalkyl groups are placed at this position. Docking, molecular dynamics (MD) simulations, and hybrid calculation methods (Quantum Mechanics/Molecular Mechanics (QM/MM)) were performed on the complexes, in order to explain these differences. Docking of the inhibitors showed the same orientation that X-ray crystal structure of the analogous (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime. MD simulations of the most active diastereomer compounds containing cis- and trans-3-hydroxycyclohexyl substituents showed stable interactions with residue Ile463 at the entrance of the B-Raf active site. On the other hand, the less active diastereomer compounds containing cis- and trans-2-hydroxycyclopentyl substituents showed interactions with inner residues Asn580 and Ser465. We found that the differences in activity can be explained by considering the dynamic interactions between the inhibitors and their surrounding residues within the B-Raf binding site. We also explained the activity trend by using a testing scoring function derived from more reliable QM/MM calculations. In addition, we search for new inhibitors from a virtual screening carried out by fragment-based de novo design. We generated a set of approximately 200 virtual compounds, which interact with Ile463 and fulfill druglikeness properties according to Lipinski, Veber, and Ghose rules.

Study of differences in the VEGFR2 inhibitory activities between semaxanib and SU5205 using 3D-QSAR, docking, and molecular dynamics simulations

Semaxanib (SU5416) and 3-[4'-fluorobenzylidene]indolin-2-one (SU5205) are structurally similar drugs that are able to inhibit vascular endothelial growth factor receptor-2 (VEGFR2), but the former is 87 times more effective than the latter. Previously, SU5205 was used as a radiolabelled inhibitor (as surrogate for SU5416) and a radiotracer for positron emission tomography (PET) imaging, but the compound exhibited poor stability and only a moderate IC(50) toward VEGFR2. In the current work, the relationship between the structure and activity of these drugs as VEGFR2 inhibitors was studied using 3D-QSAR, docking and molecular dynamics (MD) simulations. First, comparative molecular field analysis (CoMFA) was performed using 48 2-indolinone derivatives and their VEGFR2 inhibitory activities. The best CoMFA model was carried out over a training set including 40 compounds, and it included steric and electrostatic fields. In addition, this model gave satisfactory cross-validation results and adequately predicted 8 compounds contained in the test set. The plots of the CoMFA fields could explain the structural differences between semaxanib and SU5205. Docking and molecular dynamics simulations showed that both molecules have the same orientation and dynamics inside the VEGFR2 active site. However, the hydrophobic pocket of VEGFR2 was more exposed to the solvent media when it was complexed with SU5205. An energetic analysis, including Embrace and MM-GBSA calculations, revealed that the potency of ligand binding is governed by van der Waals contacts.

Investigation of the structure requirement for 5-HT₆ binding affinity of arylsulfonyl derivatives: a computational study

5-HT₆ receptor has been implicated in a series of diseases including anxiety, depression, schizophrenia and cognitive dysfunctions. 5-HT₆ ligands have been reported to play a significant role in the treatment for central nervous system (CNS) diseases. Presently, a large series of 223 5-HT₆ ligands were studied using a combinational method by 3D-QSAR, molecular docking and molecular dynamics calculations for further improvement of potency. The optimal 3D models exhibit satisfying statistical results with r²_ncv, q² values of 0.85 and 0.50 for CoMFA, 0.81 and 0.53 for CoMSIA, respectively. Their predictive powers were validated by external test set, showing r²_pred of 0.71 and 0.76. The contour maps also provide a visual representation of contributions of steric, electrostatic, hydrophobic and hydrogen bond fields as well as the prospective binding models. In addition, the agreement between 3D-QSAR, molecular docking and molecular dynamics simulation proves the rationality of the developed models. These results, we hope, may be helpful in designing novel and potential 5-HT₆ ligands.

Docking and quantitative structure-activity relationship studies for 3-fluoro-4-(pyrrolo[2,1-f][1,2,4]triazin-4-yloxy)aniline, 3-fluoro-4-(1H-pyrrolo[2,3-b]pyridin-4-yloxy)aniline, and 4-(4-amino-2-fluorophenoxy)-2-pyridinylamine derivatives as c-Met kinase inhibitors

We have performed docking of 3-fluoro-4-(pyrrolo[2,1-f][1,2,4]triazin-4-yloxy)aniline (FPTA), 3-fluoro-4-(1H-pyrrolo[2,3-b]pyridin-4-yloxy)aniline (FPPA), and 4-(4-amino-2-fluorophenoxy)-2-pyridinylamine (AFPP) derivatives complexed with c-Met kinase to study the orientations and preferred active conformations of these inhibitors. The study was conducted on a selected set of 103 compounds with variations both in structure and activity. Docking helped to analyze the molecular features which contribute to a high inhibitory activity for the studied compounds. In addition, the predicted biological activities of the c-Met kinase inhibitors, measured as IC(50) values were obtained by using quantitative structure-activity relationship (QSAR) methods: Comparative molecular similarity analysis (CoMSIA) and multiple linear regression (MLR) with topological vectors. The best CoMSIA model included steric, electrostatic, hydrophobic, and hydrogen bond-donor fields; furthermore, we found a predictive model containing 2D-autocorrelation descriptors, GETAWAY descriptors (GETAWAY: Geometry, Topology and Atom-Weight AssemblY), fragment-based polar surface area (PSA), and MlogP. The statistical parameters: cross-validate correlation coefficient and the fitted correlation coefficient, validated the quality of the obtained predictive models for 76 compounds. Additionally, these models predicted adequately 25 compounds that were not included in the training set.

C-terminal domain of p16(INK4a) is adequate in inducing cell cycle arrest, growth inhibition and CDK4/6 interaction similar to the full length protein in HT-1080 fibrosarcoma cells

The tumor suppressor p16(INK4a) has earned widespread attention in cancer studies since its discovery as an inhibitor of cyclin-dependent kinases (CDKs) 4/6. Structurally, it consists of four complete ankyrin repeats, believed to be involved in CDK4 interaction. According to the previous disparities concerning the importance of domains and inactivating mutations in p16, we aimed to search for the domain possessing the functional properties of the full length protein. Upon our in silico screening analyses followed by experimental assessments, we have identified the novel minimum functional domain of p16 to be the C-terminal half including ankyrin repeats III, IV and the C-terminal flanking region accompanied by loops 2 and 3. Transfection of this truncated form into HT-1080 human fibrosarcoma cells, lacking endogenous p16, revealed that it is able to inhibit cell growth and proliferation equivalent to p16(INK4a). The functional analysis showed that this fragment like p16 can interact with CDK4/6, block the entry into S phase of the cell cycle and suppress growth as indicated by colony formation assay. Identification of p16 minimum functional domain can be of benefit to the future peptidomimetic drug design as well as gene transfer for cancer therapy.

Computational study on the interaction of N1 substituted pyrazole derivatives with B-raf kinase: an unusual water wire hydrogen-bond network and novel interactions at the entrance of the active site

Docking and molecular dynamics (MD) simulations of N1 substituted pyrazole derivatives complexed with B-Raf kinase were performed to gain insight into the structural and energetic preferences of these inhibitors. First, a comparative study of fully automated docking programs AutoDock, ICM, GLIDE, and Surflex-Dock in closely approximating the X-ray crystal structure of the inhibitor (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime was performed. Afterward, the dynamics of the above-mentioned compound and the less active analogous compounds with 1-methyl-4-piperidinyl and tetrahydro-2H-pyran-4-yl groups at position N1 of pyrazole ring inside the B-Raf active site were analyzed by MD simulations. We found that the most active compound has stable interactions with residues Ile463 and His539 at the entrance of the B-Raf active site. Those interactions were in very good agreement with more reliable quantum mechanics/molecular mechanics calculations performed on the torsional angle phi between the pyrazole ring and the substituents at position N1. In addition, we identified a water wire connecting N2 of the pyrazole ring, Cys532, and Ser536, which is composed of three water molecules for the most active compound. We found some differences in the water wire hydrogen-bond network formed by less active compounds. We suggest that the differences between these structural features are responsible for the differences in activity among the studied compounds.

Simulation of different truncated p16(INK4a) forms and in silico study of interaction with Cdk4

Protein-protein interactions studies can greatly increase the amount of structural and functional information pertaining to biologically active molecules and processes. The information obtained from such studies can lead to design and application of new modification in order to obtain a desired bioactivity. Many application packages and servers performing docking, such as HEX, DOT, AUTODOCK, and ZDOCK are now available for predicting the lowest free energy state of a protein complex. In this study, we have focused on cyclin-dependent kinase 4 (Cdk4), a key molecule in the regulation of cell cycle progression at the G₁-S phase restriction point and p16^INK4a, a tumor suppressor which inhibits Cdk4 activity. Truncated structures were created to find the more critical regions of p16 for interaction. The tertiary structures were determined by ProSAL, GENO3D Web Server. We evaluated their interactions with Cdk4 using two docking systems, HEX 4.5 and DOT 1. Calculations were performed on a high-speed computer. Minimizations and visualizations were carried out by PdbViewer 3.7. Considering shape and shape/electrostatic total energy, structures containing ANK II, III and IV motifs that lack the N-terminal region of the full length p16 molecule showed the best fit complexes among the p16 truncated forms. The free energies were compatible with that of p16 full length original form, the full length. It seems that the N-terminal of the molecule is not crucial for the interaction since the truncated structure containing only this region did not show a good total energy.

Reducing CDK4/6-p16(INK4a) interface. Computational alanine scanning of a peptide bound to CDK6 protein

The tumor suppressor gene p16INK4a is commonly found altered in numerous and different types of cancer. The encoded protein arrests cell cycle in G1 phase by binding to CDK4 and CDK6, inhibiting their kinase function. In 1995, a 20-residue peptide, extracted from p16INK4a protein sequence, was discovered that retains the cell cycle inhibition properties of the endogenous tumor suppressor. However, its structure has not been determined yet. In this article, the features of a theoretical structure of the peptide bound to CDK6 are reported. The complex was modeled from CDK6-p16INK4a X-ray crystal structure and through molecular dynamics. Final structure was assessed by comparing computed binding free energy changes, when single-alanine substitutions were brought about on the peptide, to experimental data. Better concordance was obtained when including a high level of solvation effects. Solute-solvent vdW energy and electrostatic energy between solute and first shells of water, computed through a force field and considering explicit waters, were also to be included to achieve reasonably good concordance between theoretical and experimental data.

Loop Flexibility and Solvent Dynamics as Determinants for the Selective Inhibition of Cyclin-Dependent Kinase 4: Comparative Molecular Dynamics Simulation Studies of CDK2 and CDK4

The design and discovery of selective cyclin-dependent kinase 4 (CDK4) inhibitors have been actively pursued in order to develop therapeutic cancer treatments. By means of a consecutive computational protocol involving homology modeling, docking experiments, and molecular dynamics simulations, we examine the characteristic structural and dynamic properties that distinguish CDK4 from CDK2 in its complexation with selective inhibitors. The results for all three CDK4-selective inhibitors under investigation show that the large-amplitude motion of a disordered loop of CDK4 is damped out in the presence of the inhibitors whereas their binding in the CDK2 active site has little effect on the loop flexibility. It is also found that the binding preference of CDK4- selective inhibitors for CDK4 over CDK2 stems from the reduced solvent accessibility in the active site of the former due to the formation of a stable hydrogen-bond triad by the Asp99, Arg101, and Thr102 side chains at the top of the active-site gorge. Besides the differences in loop flexibility and solvent accessibility, the dynamic stabilities of the hydrogen bonds between the inhibitors and the side chain of the lysine residue at the bottom of the active site also correlate well with the relative binding affinities of the inhibitors for the two CDKs. These results highlight the usefulness of this computational approach in evaluating the selectivity of a CDK inhibitor, and demonstrate the necessity of considering protein flexibility and solvent effects in designing new selective CDK4-selective inhibitors.