1
|
Molecular and Biological Investigation of Isolated Marine Fungal Metabolites as Anticancer Agents: A Multi-Target Approach. Metabolites 2023; 13:metabo13020162. [PMID: 36837781 PMCID: PMC9964656 DOI: 10.3390/metabo13020162] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2022] [Revised: 01/13/2023] [Accepted: 01/19/2023] [Indexed: 01/25/2023] Open
Abstract
Cancer is the leading cause of death globally, with an increasing number of cases being annually reported. Nature-derived metabolites have been widely studied for their potential programmed necrosis, cytotoxicity, and anti-proliferation leading to enrichment for the modern medicine, particularly within the last couple of decades. At a more rapid pace, the concept of multi-target agents has evolved from being an innovative approach into a regular drug development procedure for hampering the multi-fashioned pathophysiology and high-resistance nature of cancer cells. With the advent of the Red Sea Penicillium chrysogenum strain S003-isolated indole-based alkaloids, we thoroughly investigated the molecular aspects for three major metabolites: meleagrin (MEL), roquefortine C (ROC), and isoroquefortine C (ISO) against three cancer-associated biological targets Cdc-25A, PTP-1B, and c-Met kinase. The study presented, for the first time, the detailed molecular insights and near-physiological affinity for these marine indole alkaloids against the assign targets through molecular docking-coupled all-atom dynamic simulation analysis. Findings highlighted the superiority of MEL's binding affinity/stability being quite in concordance with the in vitro anticancer activity profile conducted via sulforhodamine B bioassay on different cancerous cell lines reaching down to low micromolar or even nanomolar potencies. The advent of lengthy structural topologies via the metabolites' extended tetracyclic cores and aromatic imidazole arm permitted multi-pocket accommodation addressing the selectivity concerns. Additionally, the presence decorating polar functionalities on the core hydrophobic tetracyclic ring contributed compound's pharmacodynamic preferentiality. Introducing ionizable functionality with more lipophilic characters was highlighted to improve binding affinities which was also in concordance with the conducted drug-likeness/pharmacokinetic profiling for obtaining a balanced pharmacokinetic/dynamic profile. Our study adds to the knowledge regarding drug development and optimization of marine-isolated indole-based alkaloids for future iterative synthesis and pre-clinical investigations as multi-target anticancer agents.
Collapse
|
2
|
Verma AK, Ahmed SF, Hossain MS, Bhojiya AA, Upadhyay SK, Srivastava AK, Singh N, Harina H, Rahaman MM, Bahadur NM. Unlocking SGK1 inhibitor potential of bis-[1-N,7-N, pyrazolo tetraethoxyphthalimido{-4-(3,5-Dimethyl-4-(spiro-3-methylpyazolo)-1,7-dihydro-1H-dipyrazolo[3,4-b;4',3'-e]pyridin-8-yl)}]p-disubstituted phenyl compounds: a computational study. J Biomol Struct Dyn 2022; 40:13412-13431. [PMID: 34696688 DOI: 10.1080/07391102.2021.1988711] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
SGK1 (Serum and Glucocorticoid Regulated Kinase 1), a serine/threonine kinase that is activated by various stimuli, including serum and glucocorticoids. It controls inflammation, apoptosis, hormone release, neuro-excitability and cell proliferation, all of which play an important role in cancer progression and metastasis. SGK1 was recently proposed as a potential drug target for cancer, diabetes, and neurodegenerative diseases. In this study, molecular docking, physiochemical, toxicological properties and molecular dynamic simulation of the Bis-[1-N,7-N, Pyrazolo tetraethoxyphthalimido{-4-(3,5-Dimethyl-4-(spiro-3-methylpyazolo)-1,7-dihydro-1H-dipyrazolo[3,4-b;4',3'-e]pyridin-8-yl)}]p-disubstituted phenyl compoundsand reference EMD638683 against new SGK1 target protein. Compared to the reference inhibitor EMD638683, we choose the best compounds (series 2-6) based on the binding energy (in the range from -11.0 to -10.6 kcal/mol). With the exception of compounds 2 and 6, none of the compounds posed a risk for AMES toxicity or carcinogenicity due to their toxicological properties. 100 ns MD simulation accompanied by MM/PBSA energy calculations and PCA. According to MD simulation results, the binding of compounds 3, 4 and 5 stabilizes the SGK1 structure and causes febrile conformational changes compared to EMD638683. As a result of this research, the final selected compounds 3, 4 and 5 can be used as scaffolds to develop promising SGK1 inhibitors for the treatment of related diseases such as cancer.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Abhishek Kumar Verma
- Department of Life Sciences, Faculty of Science and Technology, Mewar University, Gangrar, Chittorgarh, Rajasthan, India
| | - Sk Faisal Ahmed
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Sonapur, Noakhali, Bangladesh
| | - Md Shahadat Hossain
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Sonapur, Noakhali, Bangladesh
| | - Ali Asger Bhojiya
- Faculty of Agriculture and Veterinary Sciences, Mewar University, Gangrar, Chittorgarh, Rajasthan, India
| | - Sudhir K Upadhyay
- Department of Environmental Science, V.B.S. Purvanchal University, Jaunpur, Uttar Pradesh, India
| | | | - Nripendra Singh
- Department of Pharmacy, V.B.S, Purvanchal University, Jaunpur, Uttar Pradesh, India
| | - Harina Harina
- Department of Life Sciences, Faculty of Science and Technology, Mewar University, Gangrar, Chittorgarh, Rajasthan, India
| | | | - Newaz Mohammed Bahadur
- Department of Applied Chemistry and Chemical Engineering, Noakhali Science and Technology University, Sonapur, Noakhali, Bangladesh
| |
Collapse
|
3
|
Ajjarapu SM, Tiwari A, Taj G, Singh DB, Singh S, Kumar S. Simulation studies, 3D QSAR and molecular docking on a point mutation of protein kinase B with flavonoids targeting ovarian Cancer. BMC Pharmacol Toxicol 2021; 22:68. [PMID: 34727985 PMCID: PMC8564994 DOI: 10.1186/s40360-021-00512-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Accepted: 07/09/2021] [Indexed: 01/17/2023] Open
Abstract
BACKGROUND Ovarian cancer is the world's dreaded disease and its prevalence is expanding globally. The study of integrated molecular networks is crucial for the basic mechanism of cancer cells and their progression. During the present investigation, we have examined different flavonoids that target protein kinases B (AKT1) protein which exerts their anticancer efficiency intriguing the role in cross-talk cell signalling, by metabolic processes through in-silico approaches. METHOD Molecular dynamics simulation (MDS) was performed to analyze and evaluate the stability of the complexes under physiological conditions and the results were congruent with molecular docking. This investigation revealed the effect of a point mutation (W80R), considered based on their frequency of occurrence, with AKT1 protein. RESULTS The ligand with high docking scores and favourable behaviour on dynamic simulations are proposed as potential W80R inhibitors. A virtual screening analysis was performed with 12,000 flavonoids satisfying Lipinski's rule of five according to which drug-likeness is predicted based on its pharmacological and biological properties to be active and taken orally. The pharmacokinetic ADME (adsorption, digestion, metabolism, and excretion) studies featured drug-likeness. Subsequently, a statistically significant 3D-QSAR model of high correlation coefficient (R2) with 0.822 and cross-validation coefficient (Q2) with 0.6132 at 4 component PLS (partial least square) were used to verify the accuracy of the models. Taxifolin holds good interactions with the binding domain of W80R, highest Glide score of - 9.63 kcal/mol with OH of GLU234 and H bond ASP274 and LEU156 amino acid residues and one pi-cation interaction and one hydrophobic bond with LYS276. CONCLUSION Natural compounds have always been a richest source of active compounds with a wide variety of structures, therefore, these compounds showed a special inspiration for medical chemists. The present study has aimed molecular docking and molecular dynamics simulation studies on taxifolin targeting W80R mutant protein of protein kinase B/serine- threonine kinase/AKT1 (EC:2.7.11.1) protein of ovarian cancer for designing therapeutic intervention. The expected result supported the molecular cause in a mutant form which resulted in a gain of ovarian cancer. Here we discussed validations computationally and yet experimental evaluation or in vivo studies are endorsed for further study. Several of these compounds should become the next marvels for early detection of ovarian cancer.
Collapse
Affiliation(s)
- Suchitra Maheswari Ajjarapu
- Bioinformatics Sub-DIC, Department of Molecular Biology & Genetic Engineering, College of Basic Science and Humanities, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Udham Singh Nagar, 263145, Uttarakhand, India
- Department of Biotechnology, Andhra University, Vishakhapatnam, 530003, Andhra Pradesh, India
| | - Apoorv Tiwari
- Bioinformatics Sub-DIC, Department of Molecular Biology & Genetic Engineering, College of Basic Science and Humanities, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Udham Singh Nagar, 263145, Uttarakhand, India
- Department of Computational Biology and Bioinformatics, Jacob Institute of Biotechnology and Bio-Engineering, Sam Higginbottom University of Agriculture, Technology and Sciences, Prayagraj, Uttar Pradesh, 211007, India
| | - Gohar Taj
- Bioinformatics Sub-DIC, Department of Molecular Biology & Genetic Engineering, College of Basic Science and Humanities, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Udham Singh Nagar, 263145, Uttarakhand, India
| | - Dev Bukhsh Singh
- Department of Biotechnology, Siddharth University, Kapilvastu, Siddharth Nagar, 272202, Uttar Pradesh, India
| | - Sakshi Singh
- Department of Molecular and Human Genetics, Banaras Hindu University, Varanasi, 221005, India
| | - Sundip Kumar
- Bioinformatics Sub-DIC, Department of Molecular Biology & Genetic Engineering, College of Basic Science and Humanities, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Udham Singh Nagar, 263145, Uttarakhand, India.
| |
Collapse
|
4
|
Oliveira ASF, Ciccotti G, Haider S, Mulholland AJ. Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:144. [PMID: 34720710 PMCID: PMC8549953 DOI: 10.1140/epjb/s10051-021-00157-0] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 07/05/2021] [Indexed: 05/03/2023]
Abstract
ABSTRACT A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in β -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol, BS8 1TQ UK
| | - Giovanni Ciccotti
- Institute for Applied Computing “Mauro Picone” (IAC), CNR, Via dei Taurini 19, 00185 Rome, Italy
- School of Physics, University College of Dublin, UCD-Belfield, Dublin 4, Ireland
- Università di Roma La Sapienza, Ple. A. Moro 5, 00185 Rome, Italy
| | - Shozeb Haider
- School of Pharmacy, University College London, London, WC1N 1AX UK
| | - Adrian J. Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
| |
Collapse
|
5
|
Uko NE, Güner OF, Matesic DF, Bowen JP. Akt Pathway Inhibitors. Curr Top Med Chem 2020; 20:883-900. [DOI: 10.2174/1568026620666200224101808] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2019] [Revised: 12/24/2019] [Accepted: 12/24/2019] [Indexed: 12/14/2022]
Abstract
Cancer is a devastating disease that has plagued humans from ancient times to this day. After
decades of slow research progress, promising drug development, and the identification of new targets,
the war on cancer was launched, in 1972. The P13K/Akt pathway is a growth-regulating cellular signaling
pathway, which in many human cancers is over-activated. Studies have demonstrated that a decrease
in Akt activity by Akt inhibitors is associated with a reduction in tumor cell proliferation. There have
been several promising drug candidates that have been studied, including but not limited to ipatasertib
(RG7440), 1; afuresertib (GSK2110183), 2; uprosertib (GSK2141795), 3; capivasertib (AZD5363), 4;
which reportedly bind to the ATP active site and inhibit Akt activity, thus exerting cytotoxic and antiproliferative
activities against human cancer cells. For most of the compounds discussed in this review,
data from preclinical studies in various cancers suggest a mechanistic basis involving hyperactivated
Akt signaling. Allosteric inhibitors are also known to alter the activity of kinases. Perifosine (KRX-
0401), 5, an alkylphospholipid, is known as the first allosteric Akt inhibitor to enter clinical development
and is mechanistically characterized as a PH-domain dependent inhibitor, non-competitive with
ATP. This results in a reduction in Akt enzymatic and cellular activities. Other small molecule (MK-
2206, 6, PHT-427, Akti-1/2) inhibitors with a similar mechanism of action, alter Akt activity through the
suppression of cell growth mediated by the inhibition of Akt membrane localization and subsequent activation.
The natural product solenopsin has been identified as an inhibitor of Akt. A few promising solenopsin
derivatives have emerged through pharmacophore modeling, energy-based calculations, and
property predictions.
Collapse
Affiliation(s)
- Nne E. Uko
- Department of Pharmaceutical Sciences, College of Pharmacy, Mercer University, Atlanta, GA 30341, United States
| | - Osman F. Güner
- Department of Chemistry and Physics, Santa Rosa Junior College, Santa Rosa, CA, United States
| | - Diane F. Matesic
- Department of Pharmaceutical Sciences, College of Pharmacy, Mercer University, Atlanta, GA 30341, United States
| | - J. Phillip Bowen
- Department of Pharmaceutical Sciences, College of Pharmacy, Mercer University, Atlanta, GA 30341, United States
| |
Collapse
|
6
|
Akhoon BA, Gandhi NS, Pandey R. Computational insights into the active structure of SGK1 and its implication for ligand design. Biochimie 2019; 165:57-66. [DOI: 10.1016/j.biochi.2019.07.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2019] [Accepted: 07/08/2019] [Indexed: 11/27/2022]
|
7
|
Botello-Smith WM, Luo Y. Robust Determination of Protein Allosteric Signaling Pathways. J Chem Theory Comput 2019; 15:2116-2126. [DOI: 10.1021/acs.jctc.8b01197] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Wesley M. Botello-Smith
- Department of Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| | - Yun Luo
- Department of Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| |
Collapse
|
8
|
Liu L, Fan S, Li W, Tao W, Shi T, Zhao YL. Theoretical Investigation of the Structural Characteristics in the Active State of Akt1 Kinase. J Chem Inf Model 2019; 59:316-325. [PMID: 30571108 DOI: 10.1021/acs.jcim.8b00506] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Akt (known as protein kinase B or PKB) is a serine/threonine kinase that regulates multiple biological processes, including cell growth, survival, and differentiation. Akt plays a critical role in the intracellular signaling network through conformational changes responsive to diverse signal inputs, and dysregulation of Akt activity could give rise to a number of diseases. However, understanding of Akt's dynamic structures and conformational transitions between active and inactive states remains unclear. In this work, classical MD simulations and QM/MM calculations were carried out to unveil the structural characteristics of Akt1, especially in its active state. The doubly protonated H194 was investigated, and both ATP-Akt1 and ADP-Akt1 complexes were constructed to demonstrate the significance of ATP in maintaining the ATP-K179-E198 salt bridge and the corresponding allosteric pathway. Besides, conformational transitions from the inactive state to the active state showed different permeation patterns of water molecules in the ATP pocket. The coordination modes of Mg2+ in the dominant representative conformations ( I and I') are presented. Unlike the water-free conformation I', three water molecules appear around Mg2+ in the water-occupied conformation I, which can finally exert an influence on the catalytic mechanism of Akt1. Furthermore, QM/MM calculations were performed to study the phosphoryl-transfer reaction of Akt1. The transfer of ATP γ-phosphate was achieved through a reversible conformational change from the reactant to a critical prereaction state, with a water molecule moving into the reaction center to coordinate with Mg2+, after which the γ-phosphate group was transferred from ATP to the substrate. Taken together, our results elucidate the structural characteristics of the Akt1 active state and shed new light on the catalytic mechanism of Akt kinases.
Collapse
Affiliation(s)
- Lanxuan Liu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, School of Life Sciences and Biotechnology , Shanghai Jiao Tong University , 800 Dongchuan Road , Shanghai 200240 , China
| | - Shuobing Fan
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, School of Life Sciences and Biotechnology , Shanghai Jiao Tong University , 800 Dongchuan Road , Shanghai 200240 , China
| | - Wenjuan Li
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, School of Life Sciences and Biotechnology , Shanghai Jiao Tong University , 800 Dongchuan Road , Shanghai 200240 , China
| | - Wentao Tao
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, School of Life Sciences and Biotechnology , Shanghai Jiao Tong University , 800 Dongchuan Road , Shanghai 200240 , China
| | - Ting Shi
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, School of Life Sciences and Biotechnology , Shanghai Jiao Tong University , 800 Dongchuan Road , Shanghai 200240 , China
| | - Yi-Lei Zhao
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, School of Life Sciences and Biotechnology , Shanghai Jiao Tong University , 800 Dongchuan Road , Shanghai 200240 , China
| |
Collapse
|
9
|
Tekpinar M, Yildirim A. Only a Subset of Normal Modes is Sufficient to Identify Linear Correlations in Proteins. J Chem Inf Model 2018; 58:1947-1961. [DOI: 10.1021/acs.jcim.8b00486] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
| | - Ahmet Yildirim
- Department of Physics, Siirt University, 56100 Siirt, Turkey
| |
Collapse
|
10
|
NDRG2 suppresses proliferation, migration, invasion and epithelial-mesenchymal transition of esophageal cancer cells through regulating the AKT/XIAP signaling pathway. Int J Biochem Cell Biol 2018. [PMID: 29530788 DOI: 10.1016/j.biocel.2018.03.003] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
N-Myc downstream-regulated gene 2 (NDRG2) has recently revealed as a candidate tumor suppressor gene. To inhibit tumor growth and decrease morbidity of esophageal cancer (EC), this study aims to test the hypothesis that the upregulation of NDRG2 may suppress proliferation, invasion, migration and epithelial-mesenchymal transition (EMT) of EC cells by regulating the AKT/XIAP signaling pathway. Immunohistochemistry was conducted for the identification of NDRG2, protein kinase B (p-AKT), X-linked inhibitor of apoptosis protein (XIAP) in EC tissues. To identify the regulatory mechanism of NDRG2 on the AKT/XIAP signaling pathway and EMT in EC, over-expressed lentiviral vector and shRNA were applied for up-regulating and interfering NDRG2 expression, and a series of determinations on the biological behavior of EC cells were performed to validate this regulation action. The results of immunohistochemistry showed NDRG2 was lowly expressed in EC tissues while p-AKT and XIAP are highly expressed. Over-expression of NDRG2 suppresses the proteins related to AKT/XIAP signaling pathway and EMT. Besides, a series of determinations shows the proliferation, migration and invasion of TE-13 cells were suppressed by over-expressed NDRG2, while the cell cycle progression was blocked and cell apoptosis was promoted. And in vivo experiment also demonstrated NDRG2 could inhibit tumor growth. Our findings demonstrate over-expression of NDRG2 works as tumor suppressive role in EC through its effects on inhibition of cell migration, invasion, and EMT by inhibiting the AKT/XIAP signaling pathway.
Collapse
|
11
|
Wieder M, Garon A, Perricone U, Boresch S, Seidel T, Almerico AM, Langer T. Common Hits Approach: Combining Pharmacophore Modeling and Molecular Dynamics Simulations. J Chem Inf Model 2017; 57:365-385. [PMID: 28072524 DOI: 10.1021/acs.jcim.6b00674] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We present a new approach that incorporates flexibility based on extensive MD simulations of protein-ligand complexes into structure-based pharmacophore modeling and virtual screening. The approach uses the multiple coordinate sets saved during the MD simulations and generates for each frame a pharmacophore model. Pharmacophore models with the same pharmacophore features are pooled. In this way the high number of pharmacophore models that results from the MD simulation is reduced to only a few hundred representative pharmacophore models. Virtual screening runs are performed with every representative pharmacophore model; the screening results are combined and rescored to generate a single hit-list. The score for a particular molecule is calculated based on the number of representative pharmacophore models which classified it as active. Hence, the method is called common hits approach (CHA). The steps between the MD simulation and the final hit-list are performed automatically and without user interaction. We test the performance of CHA for virtual screening using screening databases with active and inactive compounds for 40 protein-ligand systems. The results of the CHA are compared to the (i) median screening performance of all representative pharmacophore models of protein-ligand systems, as well as to the virtual screening performance of (ii) a random classifier, (iii) the pharmacophore model derived from the experimental structure in the PDB, and (iv) the representative pharmacophore model appearing most frequently during the MD simulation. For the 34 (out of 40) protein-ligand complexes, for which at least one of the approaches was able to perform better than a random classifier, the highest enrichment was achieved using CHA in 68% of the cases, compared to 12% for the PDB pharmacophore model and 20% for the representative pharmacophore model appearing most frequently. The availabilithy of diverse sets of different pharmacophore models is utilized to analyze some additional questions of interest in 3D pharmacophore-based virtual screening.
Collapse
Affiliation(s)
- Marcus Wieder
- Faculty of Life Sciences, Department of Pharmaceutical Chemistry, University of Vienna , Althanstraße 14, 1090 Vienna, Austria.,Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna , Währingerstraße 17, 1090 Vienna, Austria
| | - Arthur Garon
- Faculty of Life Sciences, Department of Pharmaceutical Chemistry, University of Vienna , Althanstraße 14, 1090 Vienna, Austria
| | - Ugo Perricone
- Faculty of Life Sciences, Department of Pharmaceutical Chemistry, University of Vienna , Althanstraße 14, 1090 Vienna, Austria.,Department of Biological, Chemical and Pharmaceutical Sciences and Technologies (STEBICEF), University of Palermo , Via Archirafi 32, Palermo, Italy
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna , Währingerstraße 17, 1090 Vienna, Austria
| | - Thomas Seidel
- Faculty of Life Sciences, Department of Pharmaceutical Chemistry, University of Vienna , Althanstraße 14, 1090 Vienna, Austria
| | - Anna Maria Almerico
- Department of Biological, Chemical and Pharmaceutical Sciences and Technologies (STEBICEF), University of Palermo , Via Archirafi 32, Palermo, Italy
| | - Thierry Langer
- Faculty of Life Sciences, Department of Pharmaceutical Chemistry, University of Vienna , Althanstraße 14, 1090 Vienna, Austria
| |
Collapse
|
12
|
Mou L, Cui T, Liu W, Zhang H, Cai Z, Lu S, Gao G. Microsecond molecular dynamics simulations provide insight into the ATP-competitive inhibitor-induced allosteric protection of Akt kinase phosphorylation. Chem Biol Drug Des 2016; 89:723-731. [PMID: 27797456 DOI: 10.1111/cbdd.12895] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2016] [Revised: 10/02/2016] [Accepted: 10/25/2016] [Indexed: 11/28/2022]
Abstract
Akt is a serine/threonine protein kinase, a critical mediator of growth factor-induced survival in key cellular pathways. Allosteric signaling between protein intramolecular domains requires long-range communication mediated by hotspot residues, often triggered by ligand binding. Here, based on extensive 3 μs explicit solvent molecular dynamics (MD) simulations of Akt1 kinase domain in the unbound (apo) and ATP-competitive inhibitor, GDC-0068-bound states, we propose a molecular mechanism for allosteric regulation of Akt1 kinase phosphorylation by GDC-0068 binding to the ATP-binding site. MD simulations revealed that the apo Akt1 is flexible with two disengaged N- and C-lobes, equilibrated between the open and closed conformations. GDC-0068 occupancy of the ATP-binding site shifts the conformational equilibrium of Akt1 from the open conformation toward the closed conformation and stabilizes the closed state. This effect enables allosteric signal propagation from the GDC-0068 to the phosphorylated T308 (pT308) in the activation loop and restrains phosphatase access to pT308, thereby protecting the pT308 in the GDC-0068-bound Akt1. Importantly, functional hotspots involved in the allosteric communication from the GDC-0068 to the pT308 are identified. Our analysis of GDC-0068-induced allosteric protection of Akt kinase phosphorylation yields important new insights into the molecular mechanism of allosteric regulation of Akt kinase activity.
Collapse
Affiliation(s)
- Linkai Mou
- Department of Urology, Affiliated Hospital of Weifang Medical University, Weifang, Shandong, China
| | - Tongwei Cui
- Department of Urology, Affiliated Hospital of Weifang Medical University, Weifang, Shandong, China
| | - Weiguang Liu
- Department of Urology, Affiliated Hospital of Weifang Medical University, Weifang, Shandong, China
| | - Hong Zhang
- Department of Urology, Affiliated Hospital of Weifang Medical University, Weifang, Shandong, China
| | - Zhanxiu Cai
- Department of Electrical and Computer Engineering, College of Bioscience and Technology, Weifang Medical University, Weifang, Shandong, China
| | - Shaoyong Lu
- Department of Pathophysiology, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
| | - Guojun Gao
- Department of Urology, Affiliated Hospital of Weifang Medical University, Weifang, Shandong, China
| |
Collapse
|
13
|
Li T, Wang G. Computer-aided targeting of the PI3K/Akt/mTOR pathway: toxicity reduction and therapeutic opportunities. Int J Mol Sci 2014; 15:18856-91. [PMID: 25334061 PMCID: PMC4227251 DOI: 10.3390/ijms151018856] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2014] [Revised: 09/21/2014] [Accepted: 10/08/2014] [Indexed: 12/14/2022] Open
Abstract
The PI3K/Akt/mTOR pathway plays an essential role in a wide range of biological functions, including metabolism, macromolecular synthesis, cell growth, proliferation and survival. Its versatility, however, makes it a conspicuous target of many pathogens; and the consequential deregulations of this pathway often lead to complications, such as tumorigenesis, type 2 diabetes and cardiovascular diseases. Molecular targeted therapy, aimed at modulating the deregulated pathway, holds great promise for controlling these diseases, though side effects may be inevitable, given the ubiquity of the pathway in cell functions. Here, we review a variety of factors found to modulate the PI3K/Akt/mTOR pathway, including gene mutations, certain metabolites, inflammatory factors, chemical toxicants, drugs found to rectify the pathway, as well as viruses that hijack the pathway for their own synthetic purposes. Furthermore, this evidence of PI3K/Akt/mTOR pathway alteration and related pathogenesis has inspired the exploration of computer-aided targeting of this pathway to optimize therapeutic strategies. Herein, we discuss several possible options, using computer-aided targeting, to reduce the toxicity of molecularly-targeted therapy, including mathematical modeling, to reveal system-level control mechanisms and to confer a low-dosage combination therapy, the potential of PP2A as a therapeutic target, the formulation of parameters to identify patients who would most benefit from specific targeted therapies and molecular dynamics simulations and docking studies to discover drugs that are isoform specific or mutation selective so as to avoid undesired broad inhibitions. We hope this review will stimulate novel ideas for pharmaceutical discovery and deepen our understanding of curability and toxicity by targeting the PI3K/Akt/mTOR pathway.
Collapse
Affiliation(s)
- Tan Li
- Department of Biology, South University of Science and Technology of China, 1088 Xueyuan Rd., Shenzhen 518055, China.
| | - Guanyu Wang
- Department of Biology, South University of Science and Technology of China, 1088 Xueyuan Rd., Shenzhen 518055, China.
| |
Collapse
|
14
|
Lavoie H, Li JJ, Thevakumaran N, Therrien M, Sicheri F. Dimerization-induced allostery in protein kinase regulation. Trends Biochem Sci 2014; 39:475-86. [PMID: 25220378 DOI: 10.1016/j.tibs.2014.08.004] [Citation(s) in RCA: 74] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2014] [Revised: 08/12/2014] [Accepted: 08/14/2014] [Indexed: 12/12/2022]
Abstract
The ability of protein kinases to switch between inactive and active states is critical to control the outputs of cellular signaling pathways. In several protein kinases, the conformation of helix αC is a key hub on which regulatory inputs converge to induce catalytic switching. An emerging mechanism involved in regulating helix αC orientation is the allosteric coupling with kinase domain surfaces involved in homo- or heterodimerization. In this review, we discuss dimerization-mediated regulation of the rapidly accelerated fibrosarcoma (RAF) and eIF2α kinase families and draw parallels with the analogous behavior of the epidermal growth factor receptor (EGFR) and serine/threonine-protein kinase endoribonuclease 1 (IRE1)/ribonuclease L (RNAse L) kinase families. Given that resistance to RAF-targeted therapeutics often stems from dimerization-dependent mechanisms, we suggest that a better understanding of dimerization-induced allostery may assist in developing alternate therapeutic strategies.
Collapse
Affiliation(s)
- Hugo Lavoie
- Institute for Research in Immunology and Cancer, Université de Montréal, Montréal, Québec, H3C 3J7, Canada
| | - John J Li
- The Lunenfeld-Tanenbaum Research Institute, Mount Sinai Hospital, Toronto, Ontario, M5G 1X5, Canada; Department of Molecular Genetics, University of Toronto, Toronto, Ontario, M5S 1A8, Canada
| | - Neroshan Thevakumaran
- The Lunenfeld-Tanenbaum Research Institute, Mount Sinai Hospital, Toronto, Ontario, M5G 1X5, Canada; Department of Biochemistry, University of Toronto, Toronto, Ontario, M5S 1A8, Canada
| | - Marc Therrien
- Institute for Research in Immunology and Cancer, Université de Montréal, Montréal, Québec, H3C 3J7, Canada; Département de Pathologie et Biologie Cellulaire, Université de Montréal, Montréal, Québec, H3C 3J7, Canada.
| | - Frank Sicheri
- The Lunenfeld-Tanenbaum Research Institute, Mount Sinai Hospital, Toronto, Ontario, M5G 1X5, Canada; Department of Molecular Genetics, University of Toronto, Toronto, Ontario, M5S 1A8, Canada; Department of Biochemistry, University of Toronto, Toronto, Ontario, M5S 1A8, Canada.
| |
Collapse
|
15
|
Chauvot de Beauchêne I, Allain A, Panel N, Laine E, Trouvé A, Dubreuil P, Tchertanov L. Hotspot mutations in KIT receptor differentially modulate its allosterically coupled conformational dynamics: impact on activation and drug sensitivity. PLoS Comput Biol 2014; 10:e1003749. [PMID: 25079768 PMCID: PMC4117417 DOI: 10.1371/journal.pcbi.1003749] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2013] [Accepted: 06/12/2014] [Indexed: 12/03/2022] Open
Abstract
Receptor tyrosine kinase KIT controls many signal transduction pathways and represents a typical allosterically regulated protein. The mutation-induced deregulation of KIT activity impairs cellular physiological functions and causes serious human diseases. The impact of hotspots mutations (D816H/Y/N/V and V560G/D) localized in crucial regulatory segments, the juxtamembrane region (JMR) and the activation (A-) loop, on KIT internal dynamics was systematically studied by molecular dynamics simulations. The mutational outcomes predicted in silico were correlated with in vitro and in vivo activation rates and drug sensitivities of KIT mutants. The allosteric regulation of KIT in the native and mutated forms is described in terms of communication between the two remote segments, JMR and A-loop. A strong correlation between the communication profile and the structural and dynamical features of KIT in the native and mutated forms was established. Our results provide new insight on the determinants of receptor KIT constitutive activation by mutations and resistance of KIT mutants to inhibitors. Depiction of an intra-molecular component of the communication network constitutes a first step towards an integrated description of vast communication pathways established by KIT in physiopathological contexts. Receptor tyrosine kinase KIT plays a crucial role in the regulation of cell signaling. This allosterically controlled activity may be affected by gain-of-function mutations that promote the development of several cancers. Identification of the molecular basis of KIT constitutive activation and allosteric regulation has inspired computational study of KIT hotspot mutations. In the present contribution, we investigated the mutation-induced effects on KIT conformational dynamics and intra-protein communication conditionally on the mutation location and the nature of the substituting amino acid. Our data elucidate that all studied mutations stabilize an inactive non-autoinhibited state of KIT over the inactive auto-inhibited state prevalent for the native protein. This shift in the protein conformational landscape promotes KIT constitutive activation. Our in silico analysis established correlations between the structural and dynamical effects induced by oncogenic mutations and the mutants auto-activation rates and drug sensitivities measured in vitro and in vivo. Particularly, the A-loop mutations stabilize the drug-resistant forms, while the JMR mutations may facilitate inhibitors binding to the active site. Cross-correlations established between local and long-range structural and dynamical effects demonstrate the allosteric character of the gain-of-function mutations mode of action.
Collapse
Affiliation(s)
- Isaure Chauvot de Beauchêne
- Bioinformatics, Molecular Dynamics & Modeling (BiMoDyM), Laboratoire de Biologie et Pharmacologie Appliqués (LBPA-CNRS), Ecole Normale Supérieure de Cachan, Cachan, France
| | - Ariane Allain
- Bioinformatics, Molecular Dynamics & Modeling (BiMoDyM), Laboratoire de Biologie et Pharmacologie Appliqués (LBPA-CNRS), Ecole Normale Supérieure de Cachan, Cachan, France
| | - Nicolas Panel
- Bioinformatics, Molecular Dynamics & Modeling (BiMoDyM), Laboratoire de Biologie et Pharmacologie Appliqués (LBPA-CNRS), Ecole Normale Supérieure de Cachan, Cachan, France
| | - Elodie Laine
- Bioinformatics, Molecular Dynamics & Modeling (BiMoDyM), Laboratoire de Biologie et Pharmacologie Appliqués (LBPA-CNRS), Ecole Normale Supérieure de Cachan, Cachan, France
| | - Alain Trouvé
- Centre de Mathématiques et de Leurs Applications (CMLA-CNRS), Ecole Normale Supérieure de Cachan, Cachan, France
| | - Patrice Dubreuil
- Inserm, U1068, Signaling, Hematopoiesis and Mechanism of Oncogenesis (CRCM); Institut Paoli-Calmettes; Aix-Marseille University; CNRS, UMR7258, Marseille, France
| | - Luba Tchertanov
- Bioinformatics, Molecular Dynamics & Modeling (BiMoDyM), Laboratoire de Biologie et Pharmacologie Appliqués (LBPA-CNRS), Ecole Normale Supérieure de Cachan, Cachan, France
- Centre de Mathématiques et de Leurs Applications (CMLA-CNRS), Ecole Normale Supérieure de Cachan, Cachan, France
- * E-mail:
| |
Collapse
|
16
|
Structural dynamic analysis of apo and ATP-bound IRAK4 kinase. Sci Rep 2014; 4:5748. [PMID: 25034608 PMCID: PMC4103033 DOI: 10.1038/srep05748] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2014] [Accepted: 07/02/2014] [Indexed: 12/15/2022] Open
Abstract
Interleukin-1 receptor-associated kinases (IRAKs) are Ser/Thr protein kinases that play an important role as signaling mediators in the signal transduction facilitated by the Toll-like receptor (TLR) and interleukin-1 receptor families. Among IRAK family members, IRAK4 is one of the drug targets for diseases related to the TLR and IL-1R signaling pathways. Experimental evidence suggests that the IRAK4 kinase domain is phosphorylated in its activation loop at T342, T345, and S346 in the fully activated state. However, the molecular interactions of subdomains within the active and inactive IRAK4 kinase domain are poorly understood. Hence, we employed a long-range molecular dynamics (MD) simulation to compare apo IRAK4 kinase domains (phosphorylated and unphosphorylated) and ATP-bound phosphorylated IRAK4 kinase domains. The MD results strongly suggested that lobe uncoupling occurs in apo unphosphorylated IRAK4 kinase via the disruption of the R334/T345 and R310/T345 interaction. In addition, apo unphosphorylated trajectory result in high mobility, particularly in the N lobe, activation segment, helix αG, and its adjoining loops. The Asp-Phe-Gly (DFG) and His-Arg-Asp (HRD) conserved kinase motif analysis showed the importance of these motifs in IRAK4 kinase activation. This study provides important information on the structural dynamics of IRAK4 kinase, which will aid in inhibitor development.
Collapse
|
17
|
Mou L, Li M, Lu SY, Li S, Shen Q, Zhang J, Li C, Lu X. Unraveling the Role of Arg4 and Arg6 in the Auto-Inhibition Mechanism of GSK3βFrom Molecular Dynamics Simulation. Chem Biol Drug Des 2014; 83:721-30. [PMID: 24444018 DOI: 10.1111/cbdd.12286] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2013] [Revised: 12/13/2013] [Accepted: 01/06/2014] [Indexed: 12/18/2022]
Affiliation(s)
- Linkai Mou
- Department of Urology; The Second Affiliated Hospital of Dalian Medical University; Dalian China
| | - Molin Li
- Department of Pathophysiology; Dalian Medical University; Dalian China
| | - Shao-Yong Lu
- Department of Pathophysiology and Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education; School of Medicine; Shanghai Jiao-Tong University; Shanghai China
| | - Shuai Li
- Department of Pathophysiology and Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education; School of Medicine; Shanghai Jiao-Tong University; Shanghai China
- Department of Obstetrics and Gynecology; Institute of Obstetrics and Gynecology; Renji Hospital; School of Medicine; Shanghai Jiao Tong University; Shanghai 200001 China
| | - Qiancheng Shen
- Department of Pathophysiology and Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education; School of Medicine; Shanghai Jiao-Tong University; Shanghai China
- Department of Obstetrics and Gynecology; Institute of Obstetrics and Gynecology; Renji Hospital; School of Medicine; Shanghai Jiao Tong University; Shanghai 200001 China
| | - Jian Zhang
- Department of Pathophysiology and Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education; School of Medicine; Shanghai Jiao-Tong University; Shanghai China
| | - Chuangang Li
- Department of Urology; The Second Affiliated Hospital of Dalian Medical University; Dalian China
| | - Xuefeng Lu
- Department of Obstetrics and Gynecology; Institute of Obstetrics and Gynecology; Renji Hospital; School of Medicine; Shanghai Jiao Tong University; Shanghai 200001 China
| |
Collapse
|
18
|
Chapnik N, Genzer Y, Ben-Shimon A, Niv MY, Froy O. AMPK-derived peptides reduce blood glucose levels but lead to fat retention in the liver of obese mice. J Endocrinol 2014; 221:89-99. [PMID: 24478381 DOI: 10.1530/joe-13-0625] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
AMP-activated protein kinase (AMPK) is a regulator of energy balance at both the cellular and the whole-body levels. Direct activation of AMPK has been highlighted as a potential novel, and possibly safer, alternative to treat type II diabetes and obesity. In this study, we aimed to design and characterize novel peptides that mimic the αG region of the α2 AMPK catalytic domain to modulate its activity by inhibiting interactions between AMPK domains or other interacting proteins. The derived peptides were tested in vivo and in tissue culture. The computationally predicted structure of the free peptide with the addition of the myristoyl (Myr) or acetyl (Ac) moiety closely resembled the protein structure that it was designed to mimic. Myr-peptide and Ac-peptide activated AMPK in muscle cells and led to reduced adipose tissue weight, body weight, blood glucose levels, insulin levels, and insulin resistance index, as expected from AMPK activation. In addition, triglyceride, cholesterol, leptin, and adiponectin levels were also lower, suggesting increased adipose tissue breakdown, a result of AMPK activation. On the other hand, liver weight and liver lipid content increased due to fat retention. We could not find an elevated pAMPK:AMPK ratio in the liver in vivo or in hepatocytes ex vivo, suggesting that the peptide does not lead to AMPK activation in hepatocytes. The finding that an AMPK-derived peptide leads to the activation of AMPK in muscle cells and in adipose tissue and leads to reduced glucose levels in obese mice, but to fat accumulation in the liver, demonstrates the differential effect of AMPK modulation in various tissues.
Collapse
Affiliation(s)
- Nava Chapnik
- Robert H. Smith Faculty of Agriculture, Food and Environment, Institute of Biochemistry, Food Science and Nutrition, The Hebrew University of Jerusalem, Rehovot 76100, Israel Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem, Israel
| | | | | | | | | |
Collapse
|
19
|
Yang L, He JT, Guan H, Sun YD. AKT1 Inhibitory DNAzymes Inhibit Cell Proliferation and Migration of Thyroid Cancer Cells. Asian Pac J Cancer Prev 2013; 14:2571-5. [DOI: 10.7314/apjcp.2013.14.4.2571] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
|
20
|
Papaleo E, Renzetti G. Coupled motions during dynamics reveal a tunnel toward the active site regulated by the N-terminal α-helix in an acylaminoacyl peptidase. J Mol Graph Model 2012; 38:226-34. [PMID: 23085164 DOI: 10.1016/j.jmgm.2012.06.014] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2012] [Revised: 06/15/2012] [Accepted: 06/26/2012] [Indexed: 10/28/2022]
Abstract
Acylaminoacyl peptidase (AAP) subfamily belongs to the prolyl oligopeptidase (POP) family of serine-proteases. There is a great interest in the definition of molecular mechanisms related to the activity and substrate recognition of these complex multi-domain enzymes. The active site relies at the interface between the C-terminal catalytic domain and the β-propeller domain, whose N-terminal region acts as a bridge to the hydrolase domain. In AAP, the N-terminal extension is characterized by a structurally conserved α1-helix, which is known to affect thermal stability and thermal dependence of the catalytic activity. In the present contribution, results from hundreds nanosecond all-atom molecular dynamics simulations, along with analyses of the networks of cross-correlated motions of a member of the AAP subfamily are discussed. The MD investigation identifies a tunnel that from the surrounding of the N-terminal α1-helix bring to the catalytic site. This cavity seems to be regulated by conformational changes of the α1-helix itself during the dynamics. The evidence here provided can be a useful guide for a better understanding of the mechanistic aspects related to AAP activity, but also for drug design purposes.
Collapse
Affiliation(s)
- Elena Papaleo
- Department of Biotechnology and Biosciences, University of Milano-Bicocca, P.zza della Scienza 2, 20126 Milan, Italy.
| | | |
Collapse
|
21
|
Cheng S, Prot JM, Leclerc E, Bois FY. Zonation related function and ubiquitination regulation in human hepatocellular carcinoma cells in dynamic vs. static culture conditions. BMC Genomics 2012; 13:54. [PMID: 22296956 PMCID: PMC3295679 DOI: 10.1186/1471-2164-13-54] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2011] [Accepted: 02/01/2012] [Indexed: 01/19/2023] Open
Abstract
Background Understanding hepatic zonation is important both for liver physiology and pathology. There is currently no effective systemic chemotherapy for human hepatocellular carcinoma (HCC) and its pathogenesis is of special interest. Genomic and proteomic data of HCC cells in different culture models, coupled to pathway-based analysis, can help identify HCC-related gene and pathway dysfunctions. Results We identified zonation-related expression profiles contributing to selective phenotypes of HCC, by integrating relevant experimental observations through gene set enrichment analysis (GSEA). Analysis was based on gene and protein expression data measured on a human HCC cell line (HepG2/C3A) in two culture conditions: dynamic microfluidic biochips and static Petri dishes. Metabolic activity (HCC-related cytochromes P450) and genetic information processing were dominant in the dynamic cultures, in contrast to kinase signaling and cancer-specific profiles in static cultures. That, together with analysis of the published literature, leads us to propose that biochips culture conditions induce a periportal-like hepatocyte phenotype while standard plates cultures are more representative of a perivenous-like phenotype. Both proteomic data and GSEA results further reveal distinct ubiquitin-mediated protein regulation in the two culture conditions. Conclusions Pathways analysis, using gene and protein expression data from two cell culture models, confirmed specific human HCC phenotypes with regard to CYPs and kinases, and revealed a zonation-related pattern of expression. Ubiquitin-mediated regulation mechanism gives plausible explanations of our findings. Altogether, our results suggest that strategies aimed at inhibiting activated kinases and signaling pathways may lead to enhanced metabolism-mediated drug resistance of treated tumors. If that were the case, mitigating inhibition or targeting inactive forms of kinases would be an alternative.
Collapse
Affiliation(s)
- Shu Cheng
- Université de Technologie de Compiègne, BP 20529, 60205 Compiègne Cedex, France
| | | | | | | |
Collapse
|
22
|
Hu G, Michielssens S, Moors SLC, Ceulemans A. The harmonic analysis of cylindrically symmetric proteins: a comparison of Dronpa and a DNA sliding clamp. J Mol Graph Model 2011; 34:28-37. [PMID: 22306411 DOI: 10.1016/j.jmgm.2011.12.005] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2011] [Revised: 12/21/2011] [Accepted: 12/22/2011] [Indexed: 11/24/2022]
Abstract
The harmonic analysis of two types of proteins with cylindrical symmetry is performed by the Standard Force Field Normal Mode Analysis and by the elastic network model. For both proteins the global elastic modes are assigned to their characteristic topologies. Dronpa is a rigid β-barrel structure, presenting the twisting, bending and breathing motion of a cylindrical rod. The β sliding clamp of Escherichia coli is a hexagonal β-wheel, consisting of rigid segments. In its spectrum four classes of vibrations are identified which are characteristic of an elastic torus. Correlation diagrams and RMSF analysis are compared. The results provide not only a comprehensive validation of the use of both methods to describe the elastic behavior according to the low-frequency normal modes, but also depict the correlated motions of β-barrel and β-wheel proteins. The harmonic flexibility of the Dronpa protein is compared to the principal components of molecular dynamics (MD) simulation. A functionally important localized cleft opening mode is found, which is not detected by harmonic analysis.
Collapse
Affiliation(s)
- Guang Hu
- Department of Chemistry and INPAC Institute for Nanoscale Physics and Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| | | | | | | |
Collapse
|
23
|
Seco J, Ferrer-Costa C, Campanera JM, Soliva R, Barril X. Allosteric regulation of PKCθ: understanding multistep phosphorylation and priming by ligands in AGC kinases. Proteins 2011; 80:269-80. [PMID: 22072623 DOI: 10.1002/prot.23205] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2011] [Revised: 09/12/2011] [Accepted: 09/18/2011] [Indexed: 01/13/2023]
Abstract
Protein kinases play critical roles in cellular activation and differentiation, and are involved in numerous pathophysiological processes. As a critical component of the regulatory circuitry of the cell, the kinase domain has the ability to integrate multiple signals, yielding a predetermined output. In PKC and other protein kinases of the AGC family, several phosphorylation sites control the activity, but these are in turn influenced by the presence of ligands in the binding pocket, which promotes phosphorylation. Here, we take PKC-theta as a prototypical member of the family and use molecular dynamics simulations to investigate the cross-talk that exists between regulatory and functional sites. We first show how the apo-unphosphorylated form of the kinase is populating a conformational space in which access to the ATP binding site and to the activation loop (AL) are simultaneously hindered. This could explain why the inactive state is not only catalytically incompetent but also resistant to activation. AL phosphorylation induces ATP binding site opening, which can then readily accept the cofactor. But the signal transmission mechanism works both ways, and if ligand binding to the unphosphorylated form occurs first, the AL is de-protected and becomes exposed to phosphorylation, thus providing an explanation for the paradoxical activation of PKCs by their inhibitors.
Collapse
Affiliation(s)
- Jesus Seco
- Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona, Barcelona, Spain
| | | | | | | | | |
Collapse
|
24
|
Dixit A, Verkhivker GM. Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases. PLoS Comput Biol 2011; 7:e1002179. [PMID: 21998569 PMCID: PMC3188506 DOI: 10.1371/journal.pcbi.1002179] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2011] [Accepted: 07/16/2011] [Indexed: 12/15/2022] Open
Abstract
The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level. Despite recent progress in computational and experimental studies of dynamic regulation in protein kinases, a mechanistic understanding of long-range communication and mechanisms of mutation-induced signaling controlling kinase activity remains largely qualitative. In this study, we have performed a systematic modeling and analysis of allosteric activation in ABL and EGFR kinases at the increasing level of complexity - from catalytic domain to multi-domain regulatory complexes. The results of this study have revealed organizing structural and mechanistic principles of allosteric signaling in protein kinases. Although activation mechanisms in ABL and EGFR kinases have evolved through acquisition of structurally different regulatory complexes, we have found that long-range interdomain communication between common functional segments (αF-helix and αC-helix) may be important for allosteric activation. The results of study have revealed molecular signatures of activating cancer mutations and have shed the light on general mechanistic aspects of mutation-induced signaling in protein kinases. An advanced understanding and further characterization of molecular signatures of kinase mutations may aid in a better rationalization of mutational effects on clinical outcomes and facilitate molecular-based therapeutic strategies to combat kinase mutation-dependent tumorigenesis.
Collapse
Affiliation(s)
- Anshuman Dixit
- Department of Pharmaceutical Chemistry, School of Pharmacy, The University of Kansas, Lawrence, Kansas, United States of America
| | - Gennady M. Verkhivker
- Department of Pharmaceutical Chemistry, School of Pharmacy, The University of Kansas, Lawrence, Kansas, United States of America
- Department of Pharmacology, University of California San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
25
|
Hu G, Michielssens S, Moors SLC, Ceulemans A. Normal Mode Analysis of Trp RNA Binding Attenuation Protein: Structure and Collective Motions. J Chem Inf Model 2011; 51:2361-71. [DOI: 10.1021/ci200268y] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Guang Hu
- Department of Chemistry and INPAC Institute for Nanoscale Physics and Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| | - Servaas Michielssens
- Department of Chemistry and INPAC Institute for Nanoscale Physics and Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| | - Samuel L. C. Moors
- Department of Chemistry and INPAC Institute for Nanoscale Physics and Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| | - Arnout Ceulemans
- Department of Chemistry and INPAC Institute for Nanoscale Physics and Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| |
Collapse
|
26
|
Tal-Gan Y, Hurevich M, Klein S, Ben-Shimon A, Rosenthal D, Hazan C, Shalev DE, Niv MY, Levitzki A, Gilon C. Backbone cyclic peptide inhibitors of protein kinase B (PKB/Akt). J Med Chem 2011; 54:5154-64. [PMID: 21650457 DOI: 10.1021/jm2003969] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Elevated levels of activated protein kinase B (PKB/Akt) have been detected in many types of cancer. Substrate-based peptide inhibitors have the advantage of selectivity due to their extensive interactions with the kinase-specific substrate binding site but often lack necessary pharmacological properties. Chemical modifications of potent peptide inhibitors, such as cyclization, may overcome these drawbacks while maintaining potency. We present an extensive structure-activity relationship (SAR) study of a potent peptide-based PKB/Akt inhibitor. Two backbone cyclic (BC) peptide libraries with varying modes of cyclization, bridge chemistry, and ring size were synthesized and evaluated for in vitro PKB/Akt inhibition. Backbone-to-backbone urea BC peptides were more potent than N-terminus-to-backbone amide BC peptides. Several analogues were up to 10-fold more active than the parent linear peptide. Some activity trends could be rationalized using computational surface mapping of the PKB/Akt kinase catalytic domain. The novel molecules have enhanced pharmacological properties which make them promising lead candidates.
Collapse
Affiliation(s)
- Yftah Tal-Gan
- Institute of Chemistry, The Hebrew University of Jerusalem, 91904 Jerusalem, Israel
| | | | | | | | | | | | | | | | | | | |
Collapse
|
27
|
Laine E, Chauvot de Beauchêne I, Perahia D, Auclair C, Tchertanov L. Mutation D816V alters the internal structure and dynamics of c-KIT receptor cytoplasmic region: implications for dimerization and activation mechanisms. PLoS Comput Biol 2011; 7:e1002068. [PMID: 21698178 PMCID: PMC3116893 DOI: 10.1371/journal.pcbi.1002068] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2010] [Accepted: 04/11/2011] [Indexed: 12/02/2022] Open
Abstract
The type III receptor tyrosine kinase (RTK) KIT plays a crucial role in the transmission of cellular signals through phosphorylation events that are associated with a switching of the protein conformation between inactive and active states. D816V KIT mutation is associated with various pathologies including mastocytosis and cancers. D816V-mutated KIT is constitutively active, and resistant to treatment with the anti-cancer drug Imatinib. To elucidate the activating molecular mechanism of this mutation, we applied a multi-approach procedure combining molecular dynamics (MD) simulations, normal modes analysis (NMA) and binding site prediction. Multiple 50-ns MD simulations of wild-type KIT and its mutant D816V were recorded using the inactive auto-inhibited structure of the protein, characteristic of type III RTKs. Computed free energy differences enabled us to quantify the impact of D816V on protein stability in the inactive state. We evidenced a local structural alteration of the activation loop (A-loop) upon mutation, and a long-range structural re-organization of the juxta-membrane region (JMR) followed by a weakening of the interaction network with the kinase domain. A thorough normal mode analysis of several MD conformations led to a plausible molecular rationale to propose that JMR is able to depart its auto-inhibitory position more easily in the mutant than in wild-type KIT and is thus able to promote kinase mutant dimerization without the need for extra-cellular ligand binding. Pocket detection at the surface of NMA-displaced conformations finally revealed that detachment of JMR from the kinase domain in the mutant was sufficient to open an access to the catalytic and substrate binding sites. Protein kinases are involved in a huge amount of cellular processes through phosphorylation, a crucial mechanism in cell signaling, and their misregulation often results in disease. The deactivation of protein tyrosine kinases (PTKs) or their oncogenic activation arises from mutations which affect the protein primary structure and the configuration of the enzymatic site apparently by stabilizing the activation loop (A-loop) extended conformation. Particularly, mutation D816V of receptor tyrosine kinase (RTK) KIT, found in patients with pediatric mastocytosis, acute leukemia or germ cell tumors, can be considered as the archetype of mutation inducing a displacement of the population equilibrium toward the active conformation. We present a comprehensive computational study of the activating mechanism(s) of this mutation. Our multi-approach in silico procedure evidenced a local alteration of the A-loop structure, and a long-range structural re-organization of the juxta-membrane region (JMR) followed by a weakening of the interaction network with the kinase domain. Our results provided a plausible conception of how the observed departure of JMR from kinase domain in the mutant promotes kinase mutant dimerization without requiring extra-cellular ligand binding. The pocket profiles we obtained suggested putative allosteric binding sites that could be targeted by ligands/modulators that trap the mutated enzyme.
Collapse
Affiliation(s)
- Elodie Laine
- LBPA, CNRS - ENS de Cachan, Cachan, France
- * E-mail: (EL); (LT)
| | | | | | | | - Luba Tchertanov
- LBPA, CNRS - ENS de Cachan, Cachan, France
- * E-mail: (EL); (LT)
| |
Collapse
|