1
|
Alam A, Khan MS, Mathur Y, Sulaimani MN, Farooqui N, Ahmad SF, Nadeem A, Yadav DK, Mohammad T. Structure-based identification of potential inhibitors of ribosomal protein S6 kinase 1, targeting cancer therapy: a combined docking and molecular dynamics simulations approach. J Biomol Struct Dyn 2024; 42:5758-5769. [PMID: 37365756 DOI: 10.1080/07391102.2023.2228912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 06/17/2023] [Indexed: 06/28/2023]
Abstract
Ribosomal protein S6 kinase 1 (S6K1), commonly known as P70-S6 kinase 1 (p70S6), is a key protein kinase involved in cellular signaling pathways that regulate cell growth, proliferation, and metabolism. Its significant role is reported in the PIK3/mTOR signaling pathway and is associated with various complex diseases, including diabetes, obesity, and different types of cancer. Due to its involvement in various physiological and pathological conditions, S6K1 is considered as an attractive target for drug design and discovery. One way to target S6K1 is by developing small molecule inhibitors that specifically bind to its ATP-binding site, preventing its activation and thus inhibiting downstream signaling pathways necessary for cell growth and survival. In this study, we have conducted a multitier virtual screening of a pool of natural compounds to identify potential S6K1 inhibitors. We performed molecular docking on IMPPAT 2.0 library and selected top hits based on their binding affinity, ligand efficiency, and specificity towards S6K1. The selected hits were further assessed based on different filters of drug-likeliness where two compounds (Hecogenin and Glabrene) were identified as potential leads for S6K1 inhibition. Both compounds showed appreciable affinity, ligand efficiency and specificity towards S6K1 binding pocket, drug-like properties, and stable protein-ligand complexes in molecular dynamics (MD) simulations. Finally, our study has suggested that Hecogenin and Glabrene can be potential S6K1 inhibitors which are presumably implicated in the therapeutic management of associated diseases such as diabetes, obesity, and varying types of cancer.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Afsar Alam
- Department of Computer Science, Jamia Millia Islamia, New Delhi, India
| | - Mohammad Shahzeb Khan
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi, India
| | - Yash Mathur
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi, India
| | - Md Nayab Sulaimani
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi, India
| | - Naqiya Farooqui
- Department of Biotechnology, Jamia Millia Islamia, New Delhi, India
| | - Sheikh F Ahmad
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Ahmed Nadeem
- Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Dharmendra Kumar Yadav
- Gachon Institute of Pharmaceutical Science and Department of Pharmacy, College of Pharmacy, Gachon University, Incheon, Republic of Korea
| | - Taj Mohammad
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, New Delhi, India
| |
Collapse
|
2
|
Shahab M, de Farias Morais GC, Akash S, Fulco UL, Oliveira JIN, Zheng G, Akter S. A robust computational quest: Discovering potential hits to improve the treatment of pyrazinamide-resistant Mycobacterium tuberculosis. J Cell Mol Med 2024; 28:e18279. [PMID: 38634203 PMCID: PMC11024510 DOI: 10.1111/jcmm.18279] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 03/12/2024] [Accepted: 03/18/2024] [Indexed: 04/19/2024] Open
Abstract
The rise of pyrazinamide (PZA)-resistant strains of Mycobacterium tuberculosis (MTB) poses a major challenge to conventional tuberculosis (TB) treatments. PZA, a cornerstone of TB therapy, must be activated by the mycobacterial enzyme pyrazinamidase (PZase) to convert its active form, pyrazinoic acid, which targets the ribosomal protein S1. Resistance, often associated with mutations in the RpsA protein, complicates treatment and highlights a critical gap in the understanding of structural dynamics and mechanisms of resistance, particularly in the context of the G97D mutation. This study utilizes a novel integration of computational techniques, including multiscale biomolecular and molecular dynamics simulations, physicochemical and medicinal chemistry predictions, quantum computations and virtual screening from the ZINC and Chembridge databases, to elucidate the resistance mechanism and identify lead compounds that have the potential to improve treatment outcomes for PZA-resistant MTB, namely ZINC15913786, ZINC20735155, Chem10269711, Chem10279789 and Chem10295790. These computational methods offer a cost-effective, rapid alternative to traditional drug trials by bypassing the need for organic subjects while providing highly accurate insight into the binding sites and efficacy of new drug candidates. The need for rapid and appropriate drug development emphasizes the need for robust computational analysis to justify further validation through in vitro and in vivo experiments.
Collapse
Affiliation(s)
- Muhammad Shahab
- State key laboratories of Chemical Resources Engineering Beijing, University of Chemical TechnologyBeijingChina
| | | | - Shopnil Akash
- Department of PharmacyDaffodil International UniversityDhakaBangladesh
| | - Umberto Laino Fulco
- Department of Biophysics and Pharmacology, Bioscience CenterFederal University of Rio Grande do NorteNatalRio Grande do NorteBrazil
| | - Jonas Ivan Nobre Oliveira
- Department of Biophysics and Pharmacology, Bioscience CenterFederal University of Rio Grande do NorteNatalRio Grande do NorteBrazil
| | - Guojun Zheng
- State key laboratories of Chemical Resources Engineering Beijing, University of Chemical TechnologyBeijingChina
| | - Shahina Akter
- Bangladesh Council of Scientific and Industrial ResearchDhakaBangladesh
| |
Collapse
|
3
|
Majeed A, Tahir Ul Qamar M, Maryam A, Mirza MU, Alhussain L, Al Otaibi SO, Almatroudi A, Allemailem KS, Alrumaihi F, Aloliqi AA, Alshehri FF. Structural insights into the mechanism of resistance to bicalutamide by the clinical mutations in androgen receptor in chemo-treatment resistant prostate cancer. J Biomol Struct Dyn 2024; 42:1181-1190. [PMID: 37144757 DOI: 10.1080/07391102.2023.2208203] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Accepted: 03/28/2023] [Indexed: 05/06/2023]
Abstract
Despite advanced diagnosis and detection technologies, prostate cancer (PCa) is the most prevalent neoplasms in males. Dysregulation of the androgen receptor (AR) is centrally involved in the tumorigenesis of PCa cells. Acquisition of drug resistance due to modifications in AR leads to therapeutic failure and relapse in PCa. An overhaul of comprehensive catalogues of cancer-causing mutations and their juxta positioning on 3D protein can help in guiding the exploration of small drug molecules. Among several well-studied PCa-specific mutations, T877A, T877S and H874Y are the most common substitutions in the ligand-binding domain (LBD) of the AR. In this study, we combined structure as well as dynamics-based in silico approaches to infer the mechanistic effect of amino acid substitutions on the structural stability of LBD. Molecular dynamics simulations allowed us to unveil a possible drug resistance mechanism that acts through structural alteration and changes in the molecular motions of LBD. Our findings suggest that the resistance to bicalutamide is partially due to increased flexibility in the H12 helix, which disturbs the compactness, thereby reducing the affinity for bicalutamide. In conclusion, the current study helps in understanding the structural changes caused by mutations and could assist in the drug development process.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Abdul Majeed
- Integrative Omics and Molecular Modeling Laboratory, Department of Bioinformatics and Biotechnology, Government College University, Faisalabad, Pakistan
| | - Muhammad Tahir Ul Qamar
- Integrative Omics and Molecular Modeling Laboratory, Department of Bioinformatics and Biotechnology, Government College University, Faisalabad, Pakistan
| | - Arooma Maryam
- Department of Biochemistry and Molecular Biotechnology, UMass Chan Medical School, Worcester, MA, USA
| | - Muhammad Usman Mirza
- Department of Chemistry and Biochemistry, University of Windsor, Windsor, ON, Canada
| | - Laila Alhussain
- Department of Biology, College of Science, Qassim University, Buraydah, Saudi Arabia
| | - Seham Obaid Al Otaibi
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Ahmad Almatroudi
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Khaled S Allemailem
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Faris Alrumaihi
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Abdulaziz A Aloliqi
- Department of Medical Biotechnology, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| | - Faez Falah Alshehri
- Department of Medical Laboratories, College of Applied Medical Sciences, Shaqra University, Aldawadmi, Saudi Arabia
| |
Collapse
|
4
|
Pitaloka DAE, Arfan A, Ramadhan DSF, Chaidir L. Insights from the molecular mechanism of pyrazinamide to mutated pyrazinamidase linked to the pncA gene in clinical isolates of Mycobacterium tuberculosis. J Biomol Struct Dyn 2024; 42:759-765. [PMID: 37096659 DOI: 10.1080/07391102.2023.2195002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 03/18/2023] [Indexed: 04/26/2023]
Abstract
This study aims to conduct a comprehensive molecular dynamics strategy to evaluate whether mutations found in pyrazinamide monoresistant (PZAMR) strains of Mycobacterium tuberculosis (MTB) can potentially reduce the effectiveness of pyrazinamide (PZA) for tuberculosis (TB) treatment. Five single point mutations of pyrazinamidase (PZAse), an enzyme which is responsible for the activation of prodrug PZA into pyrazinoic acid, found in MTB clinical isolates, namely His82Arg, Thr87Met, Ser66Pro, Ala171Val, and Pro62Leu, were analyzed by the dynamics simulations both in the apo state (unbound state) and in the PZA bound state. The results showed that the mutation of His82 to Arg, Thr87 to Met, and Ser66 to Pro in PZAse affects the coordination state of the Fe2+ ion, which is a cofactor required for enzyme activity. These mutations change the flexibility, stability, and fluctuation of His51, His57, and ASP49 amino acid residues around the Fe2+ ion, culminating in an unstable complex and dissociation of PZA from the PZAse binding site. However, mutations of Ala171 to Val and Pro62 to Leu were found to have no effect on the complex's stability. Based on the results, PZAse mutations of His82Arg, Thr87Met, and Ser66Pro culminated in weak binding affinity for PZA and caused significant structural deformations that led to PZA resistance. Future structural and functional studies, as well as investigations into other aspects of drug resistance in PZAse, will require experimental clarification.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Dian Ayu Eka Pitaloka
- Department of Pharmacology and Clinical Pharmacy, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
- Center for Translational Biomarker Research, Universitas Padjadjaran, Bandung, Indonesia
| | - Arfan Arfan
- Department of Medicinal Chemistry, Faculty of Pharmacy, Universitas Halu Oleo, Kendari, Indonesia
| | - Dwi Syah Fitra Ramadhan
- Department of Pharmaceutical Chemistry, Sekolah Tinggi Ilmu Kesehatan Mandala Waluya, Kendari, Indonesia
| | - Lidya Chaidir
- Center for Translational Biomarker Research, Universitas Padjadjaran, Bandung, Indonesia
- Department of Biomedical Sciences, Faculty of Medicine, Universitas Padjadjaran, Sumedang, Indonesia
| |
Collapse
|
5
|
Modlin SJ, Mansjö M, Werngren J, Ejike CM, Hoffner SE, Valafar F. Pyrazinamide-resistant Tuberculosis Obscured From Common Targeted Molecular Diagnostics. Drug Resist Updat 2023; 68:100959. [PMID: 37043916 DOI: 10.1016/j.drup.2023.100959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Revised: 03/10/2023] [Accepted: 03/21/2023] [Indexed: 04/08/2023]
Abstract
Here, we describe a clinical case of pyrazinamide-resistant (PZA-R) tuberculosis (TB) reported as PZA-susceptible (PZA-S) by common molecular diagnostics. Phenotypic susceptibility testing (pDST) indicated PZA-R TB. Targeted Sanger sequencing reported wild-type PncA, indicating PZA-S TB. Whole Genome Sequencing (WGS) by PacBio and IonTorrent both detected deletion of a large portion of pncA, indicating PZA-R. Importantly, both WGS methods showed deletion of part of the primer region targeted by Sanger sequencing. Repeating Sanger sequencing from a culture in presence of PZA returned no result, revealing that 1) two minority susceptible subpopulations had vanished, 2) the PZA-R majority subpopulation harboring the pncA deletion could not be amplified by Sanger primers, and was thus obscured by amplification process. This case demonstrates how a small susceptible subpopulation can entirely obscure majority resistant populations from targeted molecular diagnostics and falsely imply homogenous susceptibility, leading to incorrect diagnosis. To our knowledge, this is the first report of a minority susceptible subpopulation masking a majority resistant population, causing targeted molecular diagnostics to call false susceptibility. The consequence of such genomic events is not limited to PZA. This phenomenon can impact molecular diagnostics' sensitivity whenever the resistance-conferring mutation is not fully within primer-targeted regions. This can be caused by structural changes of genomic context with phenotypic consequence as we report here, or by uncommon mechanisms of resistance. Such false susceptibility calls promote suboptimal treatment and spread of strains that challenge targeted molecular diagnostics. This motivates development of molecular diagnostics unreliant on primer conservation, and impels frequent WGS surveillance for variants that evade prevailing molecular diagnostics.
Collapse
|
6
|
Alshabrmi FM, Alatawi EA. Deciphering the mechanism of resistance by novel double mutations in pncA in Mycobacterium tuberculosis using protein structural graphs (PSG) and structural bioinformatic approaches. Comput Biol Med 2023; 154:106599. [PMID: 36731361 DOI: 10.1016/j.compbiomed.2023.106599] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 12/17/2022] [Accepted: 01/22/2023] [Indexed: 01/29/2023]
Abstract
The evolution of MDR and XDR-TB is a growing concern and public health safety threat around the world. Gene mutations are the prime cause of drug resistance in tuberculosis, however the reports of double mutations further aggravated the situation. Despite the large-scale genomic sequencing and identification of novel mutations, structure investigation of the protein is still required to structurally and functionally characterize these novel mutations to design novel drugs for improved clinical outcome. Hence, we used structural bioinformatics approaches i.e. molecular modeling, residues communication and molecular simulation to understand the impact of novel double S59Y-L85P, D86G-V180F and S104G-V130 M mutation on the structure, function of pncA encoded Pyrazinamidase (PZase) and resistance of Pyrazinamide (PZA). Our results revealed that these mutations alter the binding paradigm and destabilize the protein to release the drug. Protein commination network (PCN) revealed variations in the hub residues and sub-networks which consequently alter the internal communication and signaling. The region 1-75 demonstrated higher flexibility in the mutant structures and minimal by the wild type which destabilize of the internally arranged beta-sheets which consequently reduce the binding of PZA and potentially Fe ion in the mutants. Hydrogen bonding analysis further validated the findings. The total binding free energy (ΔG) for each complex i.e. wild type -7.46 kcal/mol, S59Y-L85P -5.21 kcal/mol, S104G-V130 M -5.33 kcal/mol while for the D86G-V180F mutant the TBE was calculated to be -6.26 kcal/mol. This further confirms that these mutations reduce the binding energy of PZA for PZase and causes resistance in the effective therapy for TB. The trajectories motion was also observed to be affected by these mutations. In conclusion, these mutations use destabilizing approach to reduce the binding of PZA and causes resistance. These features can be used to design novel structure-based drugs against Tuberculosis.
Collapse
Affiliation(s)
- Fahad M Alshabrmi
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, 51452, Saudi Arabia.
| | - Eid A Alatawi
- Department of Medical Laboratory Technology, Faculty of Applied Medical Sciences, University of Tabuk, Tabuk, 71491, Saudi Arabia.
| |
Collapse
|
7
|
Ullah I, Hassan M, Khan KM, Sajid M, Umar M, Hassan S, Ullah A, El-Serehy HA, Charifi W, Yasmin H. Thiourea derivatives inhibit key diabetes-associated enzymes and advanced glycation end-product formation as a treatment for diabetes mellitus. IUBMB Life 2023; 75:161-180. [PMID: 36565478 DOI: 10.1002/iub.2699] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2022] [Accepted: 11/14/2022] [Indexed: 12/25/2022]
Abstract
This study was designed to screen novel thiourea derivatives against different enzymes, such as α-amylase, α-glucosidase, protein tyrosine phosphatase 1 B, and advanced glycated end product (AGEs). A cytotoxicity analysis was performed using rat L6 myotubes and molecular docking analysis was performed to map the binding interactions between the active compounds and α-amylase and α-glucosidase. The data revealed the potency of five compounds, including E (1-(2,4-difluorophenyl)-3-(3,4-dimethyl phenyl) thiourea), AG (1-(2-methoxy-5-(trifluoromethyl) phenyl)-3-(3-methoxy phenyl) thiourea), AF (1-(2,4-dichlorophenyl)-3-(4-ethylphenyl) thiourea), AD (1-(2,4-dichlorophenyl)-3-(4-ethylphenyl) thiourea), and AH (1-(2,4-difluorophenyl)-3-(2-iodophenyl) thiourea), showed activity against α-amylase. The corresponding percentage inhibitions were found to be 85 ± 1.9, 82 ± 0.7, 75 ± 1.2, 72 ± 0.4, and 65 ± 1.1%, respectively. These compounds were then screened using in vitro assays. Among them, AH showed the highest activity against α-glucosidase, AGEs, and PTP1B, with percentage inhibitions of 86 ± 0.4% (IC50 = 47.9 μM), 85 ± 0.7% (IC50 = 49.51 μM), and 85 ± 0.5% (IC50 = 79.74 μM), respectively. Compound AH showed an increased glucose uptake at a concentration of 100 μM. Finally, an in vivo study was conducted using a streptozotocin-induced diabetic mouse model and PTP1B expression was assessed using real-time PCR. Additionally, we examined the hypoglycemic effect of compound AH in diabetic rats compared to the standard drug glibenclamide.
Collapse
Affiliation(s)
- Imran Ullah
- Department of Biochemistry, Hazara University Mansehra, Mansehra, Khyber Pakhtunkhwa, Pakistan
| | - Mukhtiar Hassan
- Department of Biochemistry, Hazara University Mansehra, Mansehra, Khyber Pakhtunkhwa, Pakistan
| | - Khalid M Khan
- H.E.J Research Institute of Chemistry, University of Karachi, Karachi, Pakistan
| | - Muhammad Sajid
- Department of Biochemistry, Hazara University Mansehra, Mansehra, Khyber Pakhtunkhwa, Pakistan
| | - Muhammad Umar
- Department of Orthopaedics, Shenzhen University General Hospital, Shenzhen, Guangdong, People's Republic of China
| | - Said Hassan
- Institute of Biotechnology and Microbiology, Bacha Khan University Charsadda, Charsadda, Khyber Pakhtunkhwa, Pakistan
| | - Amin Ullah
- Department of Health and Biological Sciecnes, Abasyn University Peshawar, Peshawar, Khyber Pakhtunkhwa, Pakistan
| | - Hamed A El-Serehy
- Department of Zoology, College of Science, King Saud University, Riyadh, Saudi Arabia
| | - Wafa Charifi
- Cochin Institute, University of Paris, INSERM, U1016, Paris, France
| | - Humaira Yasmin
- Department of Biosciences, COMSATS University Islamabad, Islamabad, Pakistan
| |
Collapse
|
8
|
Validation of catalytic site residues of Ubiquitin Specific Protease 2 (USP2) by molecular dynamic simulation and novel kinetics assay for rational drug design. Mol Divers 2022:10.1007/s11030-022-10499-1. [PMID: 35932436 DOI: 10.1007/s11030-022-10499-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Accepted: 06/29/2022] [Indexed: 10/15/2022]
Abstract
Post-translational modifications of proteins such as protein ubiquitination are crucial for regulating conformation, stability and localization of the modified protein. Ubiquitin-specific protease 2 (USP2), a multifunctional cysteine protease is reported to be a key regulator of ubiquitylation events in numerous oncogenic proteins e.g., fatty acid synthetase, Mdm2, EGFR, cyclin A1, and cyclin-D1, etc. Thus targeting USP2 is a promising strategy for cancer therapy. USP2 is characterized by a catalytic triad comprising of cysteine, histidine and aspartic acid residues. Five residues including three from the catalytic triad and two from outside of the catalytic triad have been reported as a catalytic site of USP2 that catalyze hydrolysis and stabilizes the oxyanion formed in the intermediate step of catalysis. Here, we report two more novel residues (L269 and Y558) on USP2 involved in the catalysis of Ubiquitin using computational alanine scanning (CAS) followed by molecular dynamic simulation studies. The results obtained from CAS were further validated by a highly reliable, time- and cost-effective SDS-PAGE-based kinetics assay using UBA52 which is a natural substrate of USP2. Our results showed that mutating L269 and Y558 significantly compromised the catalytic efficiency of USP2 in hydrolyzing UBA52 which can further be extended to rational drug design of USP2 selective inhibitors and to explore the catalytic sites of other USPs. Two novel residues take part in catalytic activity of USP2 which were depicted by MD Simulations and were further validated by novel SDS-PAGE-based reliable time- and cost-effective kinetics assay.
Collapse
|
9
|
Mou K, Abdalla M, Wei DQ, Khan MT, Lodhi MS, Darwish DB, Sharaf M, Tu X. Emerging mutations in envelope protein of SARS-CoV-2 and their effect on thermodynamic properties. INFORMATICS IN MEDICINE UNLOCKED 2021; 25:100675. [PMID: 34337139 PMCID: PMC8314890 DOI: 10.1016/j.imu.2021.100675] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Revised: 07/11/2021] [Accepted: 07/19/2021] [Indexed: 12/23/2022] Open
Abstract
Structural proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are potential drug targets due to their role in the virus life cycle. The envelope (E) protein is one of the structural proteins; plays a critical role in virulency. However, the emergence of mutations oftenly leads to drug resistance and may also play a vital role in virus stabilization and evolution. In this study, we aimed to identify mutations in E proteins that affect the protein stability. About 0.3 million complete whole genome sequences were analyzed to screen mutations in E protein. All these mutations were subjected to stability prediction using the DynaMut server. The most common mutations that were detected at the C-terminal domain, Ser68Phe, Pro71Ser, and Leu73Phe, were examined through molecular dynamics (MD) simulations for a 100ns period. The sequence analysis shows the existence of 259 mutations in E protein. Interestingly, 16 of them were detected in the DFLV amino acid (aa) motif (aa72-aa75) that binds the host PALS1 protein. The results of root mean square deviation, fluctuations, radius of gyration, and free energy landscape show that Ser68Phe, Pro71Ser, and Leu73Phe are exhibiting a more stabilizing effect. However, a more comprehensive experimental study may be required to see the effect on virus pathogenicity. Potential antiviral drugs, and vaccines may be developed used after screening the genomic variations for better management of SARS-CoV-2 infections.
Collapse
Affiliation(s)
- Kejie Mou
- Department of Neurosurgery, Bishan Hospital of Chongqing, Chongqing, China
| | - Mohnad Abdalla
- Key Laboratory of Chemical Biology (Ministry of Education), Department of Pharmaceutics, School of Pharmaceutical Sciences, Cheeloo College of Medicine, Shandong University, 44 Cultural West Road, Shandong Province, 250012, PR China
| | - Dong Qing Wei
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200030, PR China
- Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nashan District, Shenzhen, Guangdong, 518055, PR China
| | - Muhammad Tahir Khan
- Institute of Molecular Biology and Biotechnology (IMBB), The University of Lahore, KM Defence Road, Lahore, Pakistan, 58810
| | - Madeeha Shahzad Lodhi
- Institute of Molecular Biology and Biotechnology (IMBB), The University of Lahore, KM Defence Road, Lahore, Pakistan, 58810
| | - Doaa B Darwish
- Department of Biology, Faculty of Science, University of Tabuk, 71491, Saudi Arabia
| | - Mohamed Sharaf
- Department of Biochemistry and Molecular Biology, College of Marine Life Sciences, Ocean University of China, Qingdao, 266003, PR China
- Department of Biochemistry, Faculty of Agriculture, AL-Azhar University, Nasr City, Cairo, 11751, Egypt
| | - Xudong Tu
- Chongqing Medical and Pharmaceutical College, Chongqing, China
| |
Collapse
|
10
|
Nangraj AS, Khan A, Umbreen S, Sahar S, Arshad M, Younas S, Ahmad S, Ali S, Ali SS, Ali L, Wei DQ. Insights Into Mutations Induced Conformational Changes and Rearrangement of Fe 2+ Ion in pncA Gene of Mycobacterium tuberculosis to Decipher the Mechanism of Resistance to Pyrazinamide. Front Mol Biosci 2021; 8:633365. [PMID: 34095218 PMCID: PMC8174790 DOI: 10.3389/fmolb.2021.633365] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 04/07/2021] [Indexed: 11/15/2022] Open
Abstract
Pyrazinamide (PZA) is the first-line drug commonly used in treating Mycobacterium tuberculosis (Mtb) infections and reduces treatment time by 33%. This prodrug is activated and converted to an active form, Pyrazinoic acid (POA), by Pyrazinamidase (PZase) enzyme. Mtb resistance to PZA is the outcome of mutations frequently reported in pncA, rpsA, and panD genes. Among the mentioned genes, pncA mutations contribute to 72-99% of the total resistance to PZA. Thus, considering the vital importance of this gene in PZA resistance, its frequent mutations (D49N, Y64S, W68G, and F94A) were investigated through in-depth computational techniques to put conclusions that might be useful for new scaffolds design or structure optimization to improve the efficacy of the available drugs. Mutants and wild type PZase were used in extensive and long-run molecular dynamics simulations in triplicate to disclose the resistance mechanism induced by the above-mentioned point mutations. Our analysis suggests that these mutations alter the internal dynamics of PZase and hinder the correct orientation of PZA to the enzyme. Consequently, the PZA has a low binding energy score with the mutants compared with the wild type PZase. These mutations were also reported to affect the binding of Fe2+ ion and its coordinated residues. Conformational dynamics also revealed that β-strand two is flipped, which is significant in Fe2+ binding. MM-GBSA analysis confirmed that these mutations significantly decreased the binding of PZA. In conclusion, these mutations cause conformation alterations and deformities that lead to PZA resistance.
Collapse
Affiliation(s)
- Asma Sindhoo Nangraj
- Department of Bioinformatics and Biological Statistics, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| | - Abbas Khan
- Department of Bioinformatics and Biological Statistics, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| | | | - Sana Sahar
- The Islamia University of Bahawalpur, Bahawalpur, Pakistan
| | - Maryam Arshad
- Government College University Faisalabad, Sahiwal, Pakistan
| | | | - Sajjad Ahmad
- Department of Health and Biological Sciences, Abasyn University, Peshawar, Pakistan
| | - Shahid Ali
- Center for Biotechnology and Microbiology, University of Swat, Swat, Pakistan
| | - Syed Shujait Ali
- Center for Biotechnology and Microbiology, University of Swat, Swat, Pakistan
| | - Liaqat Ali
- Department of Biological Sciences, National University of Medical Sciences, Islamabad, Pakistan
| | - Dong-Qing Wei
- Department of Bioinformatics and Biological Statistics, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
- Peng Cheng Laboratory, Shenzhen, China
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic and Developmental Sciences, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| |
Collapse
|
11
|
Chen J, Wang W, Pang L, Zhu W. Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics. Phys Chem Chem Phys 2021; 22:21238-21250. [PMID: 32930679 DOI: 10.1039/d0cp03766d] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China.
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China.
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan 250357, China.
| | - Weiliang Zhu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China.
| |
Collapse
|
12
|
Khan MT, Chinnasamy S, Cui Z, Irfan M, Wei DQ. Mechanistic analysis of A46V, H57Y, and D129N in pyrazinamidase associated with pyrazinamide resistance. Saudi J Biol Sci 2020; 27:3150-3156. [PMID: 33100877 PMCID: PMC7569123 DOI: 10.1016/j.sjbs.2020.07.015] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 07/10/2020] [Accepted: 07/12/2020] [Indexed: 12/20/2022] Open
Abstract
Pyrazinamide (PZA) is a component of first-line drugs, active against latent Mycobacterium tuberculosis (MTB) isolates. The prodrug is activated into the active form, pyrazinoic acid (POA) via pncA gene-encoded pyrazinamidase (PZase). Mutations in pncA have been reported, most commonly responsible for PZA-resistance in more than 70% of the resistant cases. In our previous study, we detected many mutations in PZase among PZA-resistance MTB isolates including A46V, H71Y, and D129N. The current study was aimed to investigate the molecular mechanism of PZA-resistance behind mutants (MTs) A46V, H71Y, and D129N in comparison with the wild type (WT) through molecular dynamic (MD) simulation. MTB positive samples were subjected to PZA drug susceptibility testing (DST) against critical concentration (100ug/ml). The resistant samples were subjected to pncA sequencing. Thirty-six various mutations have been observed in the coding region of pncA of PZA-resistant isolates (GenBank accession No. MH461111) including A46V, H71Y, and D129N. The post-simulation analysis revealed a significant variation in MTs structural dynamics as compared to the WT. Root means square deviations (RMSD) and Root means square fluctuation (RMSF) has been found in variation between WT and MTs. Folding effect and pocket volume were altered in MTs when compared with WT. Geometric matching supports the effect of mutation A46V, H71Y, and D129N on PZase structure that may have an insight effect on PZase dynamics, making them vulnerable to convert pro-PZA into active form, POA. In conclusion, the current analyses will provide useful information behind PZA-resistance for better management of drug-resistant TB.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Pakistan
| | - Sathishkumar Chinnasamy
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Zhilei Cui
- Department of Respiratory Medicine, XinHua Hospital Affiliated to Shanghai Jiao Tong University School of Medicine, China
| | - Muhammad Irfan
- Department of Microbiology and Cell Science, Genetics Institute and Institute of Food and Agricultural Sciences, University of Florida, Gainesville, FL, USA
| | - Dong-Qing Wei
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240, China
- Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District, Shenzhen, Guangdong 518055, China
| |
Collapse
|
13
|
Ali A, Khan MT, Khan A, Ali S, Chinnasamy S, Akhtar K, Shafiq A, Wei DQ. Pyrazinamide resistance of novel mutations in pncA and their dynamic behavior. RSC Adv 2020; 10:35565-35573. [PMID: 35515677 PMCID: PMC9056903 DOI: 10.1039/d0ra06072k] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2020] [Accepted: 08/25/2020] [Indexed: 11/21/2022] Open
Abstract
Pyrazinamide (PZA) is one of the essential anti-mycobacterium drugs, active against non-replicating Mycobacterium tuberculosis (MTB) isolates. PZA is converted into its active state, called pyrazinoic acid (POA), by action of pncA encoding pyrazinamidase (PZase). In the majority of PZA-resistance isolates, pncA harbored mutations in the coding region. In our recent report, we detected a number of novel variants in PZA-resistance (PZAR) MTB isolates, whose resistance mechanisms were yet to be determined. Here we performed several analyses to unveil the PZAR mechanism of R123P, T76P, G150A, and H71R mutants (MTs) through molecular dynamics (MD) simulations. In brief, culture positive MTB isolates were subjected to PZA susceptibility tests using the WHO recommended concentration of PZA (100 μg ml−1). The PZAR samples were screened for mutations in pncA along sensitive isolates through polymerase chain reactions and sequencing. A large number of variants (GeneBank accession no. MH461111), including R123P, T76P, G150A, and H71R, have been spotted in more than 70% of isolates. However, the mechanism of PZAR for mutants (MTs) R123P, T76P, G150A, and H71R was unknown. For the MTs and native PZase structures (WT), thermodynamic properties were compared using molecular dynamics simulations for 100 ns. The MTs structural activity was compared to the WT. Folding effect and pocket volume variations have been detected when comparing between WT and MTs. Geometric matching further confirmed the effect of R123P, T76P, G150A, and H71R mutations on PZase dynamics, making them vulnerable for activating the pro-drug into POA. This study offers a better understanding for management of PZAR TB. The results may be used as alternative diagnostic tools to infer PZA resistance at a structural dynamics level. We performed several analyses to unveil the pyrazinamide-resistance mechanism of R123P, T76P, G150A, and H71R mutants through molecular dynamics simulations.![]()
Collapse
Affiliation(s)
- Arif Ali
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University 800 Dongchuan Road Shanghai, Minhang District Shanghai 200240 China +86-21-3420-4573
| | - Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology Pakistan
| | - Abbas Khan
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University 800 Dongchuan Road Shanghai, Minhang District Shanghai 200240 China +86-21-3420-4573
| | - Sajid Ali
- Quaid-i-Azam University Islamabad, Provincial Tuberculosis Reference Laboratory Hayatabad Medical Complex Peshawar Pakistan
| | - Sathishkumar Chinnasamy
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University 800 Dongchuan Road Shanghai, Minhang District Shanghai 200240 China +86-21-3420-4573
| | - Khalid Akhtar
- National University of Science and Technology Pakistan
| | - Athar Shafiq
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University 800 Dongchuan Road Shanghai, Minhang District Shanghai 200240 China +86-21-3420-4573
| | - Dong-Qing Wei
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University 800 Dongchuan Road Shanghai, Minhang District Shanghai 200240 China +86-21-3420-4573.,Peng Cheng Laboratory Vanke Cloud City Phase I Building 8, Xili Street, Nanshan District Shenzhen Guangdong 518055 China
| |
Collapse
|
14
|
Supo-Escalante RR, Médico A, Gushiken E, Olivos-Ramírez GE, Quispe Y, Torres F, Zamudio M, Antiparra R, Amzel LM, Gilman RH, Sheen P, Zimic M. Prediction of Mycobacterium tuberculosis pyrazinamidase function based on structural stability, physicochemical and geometrical descriptors. PLoS One 2020; 15:e0235643. [PMID: 32735615 PMCID: PMC7394417 DOI: 10.1371/journal.pone.0235643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Accepted: 06/19/2020] [Indexed: 12/02/2022] Open
Abstract
BACKGROUND Pyrazinamide is an important drug against the latent stage of tuberculosis and is used in both first- and second-line treatment regimens. Pyrazinamide-susceptibility test usually takes a week to have a diagnosis to guide initial therapy, implying a delay in receiving appropriate therapy. The continued increase in multi-drug resistant tuberculosis and the prevalence of pyrazinamide resistance in several countries makes the development of assays for prompt identification of resistance necessary. The main cause of pyrazinamide resistance is the impairment of pyrazinamidase function attributed to mutations in the promoter and/or pncA coding gene. However, not all pncA mutations necessarily affect the pyrazinamidase function. OBJECTIVE To develop a methodology to predict pyrazinamidase function from detected mutations in the pncA gene. METHODS We measured the catalytic constant (kcat), KM, enzymatic efficiency, and enzymatic activity of 35 recombinant mutated pyrazinamidase and the wild type (Protein Data Bank ID = 3pl1). From all the 3D modeled structures, we extracted several predictors based on three categories: structural stability (estimated by normal mode analysis and molecular dynamics), physicochemical, and geometrical characteristics. We used a stepwise Akaike's information criterion forward multiple log-linear regression to model each kinetic parameter with each category of predictors. We also developed weighted models combining the three categories of predictive models for each kinetic parameter. We tested the robustness of the predictive ability of each model by 6-fold cross-validation against random models. RESULTS The stability, physicochemical, and geometrical descriptors explained most of the variability (R2) of the kinetic parameters. Our models are best suited to predict kcat, efficiency, and activity based on the root-mean-square error of prediction of the 6-fold cross-validation. CONCLUSIONS This study shows a quick approach to predict the pyrazinamidase function only from the pncA sequence when point mutations are present. This can be an important tool to detect pyrazinamide resistance.
Collapse
Affiliation(s)
- Rydberg Roman Supo-Escalante
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Aldhair Médico
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Eduardo Gushiken
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Gustavo E. Olivos-Ramírez
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Yaneth Quispe
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Fiorella Torres
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Melissa Zamudio
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Ricardo Antiparra
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - L. Mario Amzel
- Department of Biophysics and Biophysical Chemistry, Johns Hopkins University, Baltimore, MD, United States of America
| | - Robert H. Gilman
- International Health Department, Johns Hopkins School of Public Health, Baltimore, MD, United States of America
| | - Patricia Sheen
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Mirko Zimic
- Laboratorio de Bioinformática, Biología Molecular y Desarrollos Tecnológicos, Facultad de Ciencias y Filosofía, Universidad Peruana Cayetano Heredia, Lima, Peru
| |
Collapse
|
15
|
Sheik Amamuddy O, Musyoka TM, Boateng RA, Zabo S, Tastan Bishop Ö. Determining the unbinding events and conserved motions associated with the pyrazinamide release due to resistance mutations of Mycobacterium tuberculosis pyrazinamidase. Comput Struct Biotechnol J 2020; 18:1103-1120. [PMID: 32489525 PMCID: PMC7251373 DOI: 10.1016/j.csbj.2020.05.009] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2020] [Revised: 04/23/2020] [Accepted: 05/06/2020] [Indexed: 01/04/2023] Open
Abstract
Pyrazinamide (PZA) is the only first-line antitubercular drug active against latent Mycobacterium tuberculosis (Mtb). It is activated to pyrazinoic acid by the pncA-encoded pyrazinamidase enzyme (PZase). Despite the emergence of PZA drug resistance, the underlying mechanisms of resistance remain unclear. This study investigated part of these mechanisms by modelling a PZA-bound wild type and 82 mutant PZase structures before applying molecular dynamics (MD) with an accurate Fe2+ cofactor coordination geometry. After observing nanosecond-scale PZA unbinding from several PZase mutants, an algorithm was developed to systematically detect ligand release via centre of mass distances (COM) and ligand average speed calculations, before applying the statistically guided network analysis (SGNA) method to investigate conserved protein motions associated with ligand unbinding. Ligand and cofactor perspectives were also investigated. A conserved pair of lid-destabilising motions was found. These consisted of (1) antiparallel lid and side flap motions; (2) the contractions of a flanking region within the same flap and residue 74 towards the core. Mutations affecting the hinge residues (H51 and H71), nearby residues or L19 were found to destabilise the lid. Additionally, other metal binding site (MBS) mutations delocalised the Fe2+ cofactor, also facilitating lid opening. In the early stages of unbinding, a wider variety of PZA poses were observed, suggesting multiple exit pathways. These findings provide insights into the late events preceding PZA unbinding, which we found to occur in some resistant PZase mutants. Further, the algorithm developed here to identify unbinding events coupled with SGNA can be applicable to other similar problems.
Collapse
Key Words
- 3D, Three-dimensional
- ACPYPE, AnteChamber Python Parser interface
- Amber force field parameters
- CHPC, Center for High Performance Computing
- COM, Center of mass
- Drug resistance
- Drug unbinding
- FDA, Food and Drug Administration
- HTMD, High throughput molecular dynamics
- INH, Isoniazid
- MBS, Metal binding site
- MCBP, Metal Center Parameter Builder
- MD, Molecular dynamics
- MDR-TB, Multidrug-resistant tuberculosis
- Missense mutations
- Molecular dynamics simulations
- PBC, Periodic boundary conditions
- PDB, Protein Data bank
- POA, Pyrazinoic acid
- PZA, Pyrazinamide
- PZase, Pyrazinamidase
- QM, Quantum Mechanics
- RIF, Rifampicin
- SGNA, Statistically guided network analysis
- Statistically guided network analysis
- TB, Tuberculosis
- VAPOR, Variant Analysis Portal
- WHO, World Health Organization
- WT, Wild type
Collapse
Affiliation(s)
| | | | - Rita Afriyie Boateng
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| | - Sophakama Zabo
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| |
Collapse
|
16
|
Junaid M, Li CD, Li J, Khan A, Ali SS, Jamal SB, Saud S, Ali A, Wei DQ. Structural insights of catalytic mechanism in mutant pyrazinamidase of Mycobacterium tuberculosis. J Biomol Struct Dyn 2020; 39:3172-3185. [PMID: 32340563 DOI: 10.1080/07391102.2020.1761879] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Pyrazinamidase (PZase) is a member of Fe-dependent amidohydrolases that activates pyrazinamide (PZA) into active pyrazinoic acid (POA). PZA, a nicotinamide analogue, is an essential first-line drug used in Mycobacterium tuberculosis (Mtb) treatment. The active form of PZA, POA, is toxic and potently inhibits the growth of latent Mtb, which makes it possible to shorten the conventional 9-month tuberculosis treatment to 6 months. In this study, an extensive molecular dynamics simulation was carried out to the study the resistance mechanism offered by the three mutations Q10P and D12A and G97D. Our results showed that two regions Gln10-His43, Phe50-Gly75 are profoundly affected by these mutations. Among the three mutations, Q10P and D12A mutations strongly disturb the communication among the catalytic triad (Asp8, Lys98 and Cys138). The oxyanion hole is formed between the backbone nitrogen atoms of A134 and C138 which stabilizes the hydroxyl anion of nicotinamide. The D12A mutation greatly disturbs the oxyanion hole formation followed by the Q10P and G97D. Our results also showed that these mutations destabilize the interaction between Fe2+ ion and Asp49, His51, His57 and His71. The binding pocket analysis showed that these mutations increase the cavity volume, which results in loose binding of PZA. MMGBSA analyzes have shown that these mutations reduce the binding affinity to the PZA drug. Our results may provide useful information for the design of new and effective PZase inhibitors based on structural information of WT and mutant PZases.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Muhammad Junaid
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, China.,Peng Cheng Laboratory, Shenzhen, Guangdong, P.R China.,Ministry of Education, Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, P.R China
| | - Cheng-Dong Li
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, China.,Peng Cheng Laboratory, Shenzhen, Guangdong, P.R China.,Ministry of Education, Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, P.R China
| | - Jiayi Li
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, China
| | - Abbas Khan
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, China.,Peng Cheng Laboratory, Shenzhen, Guangdong, P.R China.,Ministry of Education, Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, P.R China
| | - Syed Shujait Ali
- Center for Biotechnology and Microbiology, University of Swat, Swat, Pakistan
| | - Syed Baber Jamal
- Department of biological sciences, National University of Medical Sciences, Punjab, Rawalpendi, Pakistan
| | - Shah Saud
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, China
| | - Arif Ali
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, China.,Peng Cheng Laboratory, Shenzhen, Guangdong, P.R China.,Ministry of Education, Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, P.R China
| | - Dong-Qing Wei
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, China.,Peng Cheng Laboratory, Shenzhen, Guangdong, P.R China.,Ministry of Education, Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, P.R China
| |
Collapse
|
17
|
Khan MT, Ali S, Zeb MT, Kaushik AC, Malik SI, Wei DQ. Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance. Front Mol Biosci 2020; 7:52. [PMID: 32328498 PMCID: PMC7160322 DOI: 10.3389/fmolb.2020.00052] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 03/16/2020] [Indexed: 12/16/2022] Open
Abstract
A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Sajid Ali
- Department of Microbiology, Quaid-i-Azam University Islamabad, Islamabad, Pakistan
| | | | - Aman Chandra Kaushik
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
| | - Shaukat Iqbal Malik
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Dong-Qing Wei
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
- Peng Cheng Laboratory, Shenzhen, China
| |
Collapse
|
18
|
Mehmood A, Khan MT, Kaushik AC, Khan AS, Irfan M, Wei DQ. Structural Dynamics Behind Clinical Mutants of PncA-Asp12Ala, Pro54Leu, and His57Pro of Mycobacterium tuberculosis Associated With Pyrazinamide Resistance. Front Bioeng Biotechnol 2019; 7:404. [PMID: 31921809 PMCID: PMC6914729 DOI: 10.3389/fbioe.2019.00404] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2019] [Accepted: 11/26/2019] [Indexed: 11/15/2022] Open
Abstract
Pyrazinamide (PZA) is one of the main FDA approved drugs to be used as the first line of defense against Mycobacterium Tuberculosis (MTB). It is activated into pyrazinoic acid (POA) via MTB's pncA gene-encoded pyrazinamidase (PZase). Mutations are most commonly responsible for PZA-resistance in nearly 70% of the resistant samples. In the present work, MTB positive samples were chosen for PZA drug susceptibility testing (DST) against critical concentration (100 ug/ml) of PZA. The resistant samples were subjected to pncA sequencing. As a result, 36 various mutations have been observed in the PZA resistant samples, uploaded to the NCBI (GeneBank accession no. MH461111). Here we report the mechanism of PZA resistance behind the three mutants (MTs), Asp12Ala, Pro54Leu, and His57Pro in comparison with the wild type (WT) through molecular dynamics simulation to unveil how these mutations affect the overall conformational stability. The post-simulation analyses revealed notable deviations as compared to the WT structure. Molecular docking studies of PZA with MTs and WT, pocket volume inspection and overall shape complementarity analysis confirmed the deleterious nature of these mutations and gave an insight into the mechanism behind PZA-resistance. These analyses provide vital information regarding MTB drug resistance and could be extremely useful in therapy management and overcoming its global burden.
Collapse
Affiliation(s)
- Aamir Mehmood
- The State Key Laboratory of Microbial Metabolism, College of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| | - Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | | | - Anwar Sheed Khan
- Department of Microbiology, Kohat University of Science and Technology, Kohat, Pakistan
| | - Muhammad Irfan
- Department of Microbiology and Cell Science, Genetics Institute and Institute of Food and Agricultural Sciences, University of Florida, Gainesville, FL, United States
| | - Dong-Qing Wei
- The State Key Laboratory of Microbial Metabolism, College of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| |
Collapse
|
19
|
Khan MT, Kaushik AC, Bhatti AI, Zhang YJ, Zhang S, Wei AJ, Malik SI, Wei DQ. Marine Natural Products and Drug Resistance in Latent Tuberculosis. Mar Drugs 2019; 17:md17100549. [PMID: 31561525 PMCID: PMC6836121 DOI: 10.3390/md17100549] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2019] [Revised: 09/05/2019] [Accepted: 09/06/2019] [Indexed: 12/16/2022] Open
Abstract
Pyrazinamide (PZA) is the only drug for the elimination of latent Mycobacterium tuberculosis (MTB) isolates. However, due to the increased number of PZA-resistance, the chances of the success of global TB elimination seems to be more prolonged. Recently, marine natural products (MNPs) as an anti-TB agent have received much attention, where some compounds extracted from marine sponge, Haliclona sp. exhibited strong activity under aerobic and hypoxic conditions. In this study, we screened articles from 1994 to 2019 related to marine natural products (MNPs) active against latent MTB isolates. The literature was also mined for the major regulators to map them in the form of a pathway under the dormant stage. Five compounds were found to be more suitable that may be applied as an alternative to PZA for the better management of resistance under latent stage. However, the mechanism of actions behind these compounds is largely unknown. Here, we also applied synthetic biology to analyze the major regulatory pathway under latent TB that might be used for the screening of selective inhibitors among marine natural products (MNPs). We identified key regulators of MTB under latent TB through extensive literature mining and mapped them in the form of regulatory pathway, where SigH is negatively regulated by RshA. PknB, RshA, SigH, and RNA polymerase (RNA-pol) are the major regulators involved in MTB survival under latent stage. Further studies are needed to screen MNPs active against the main regulators of dormant MTB isolates. To reduce the PZA resistance burden, understanding the regulatory pathways may help in selective targets of MNPs from marine natural sources.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad 44000, Pakistan; (M.T.K.); (S.I.M.)
| | - Aman Chandra Kaushik
- The State Key Laboratory of Microbial Metabolism, College of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China;
| | - Aamer Iqbal Bhatti
- Department of Electrical Engineering, Capital University of Science and Technology, Islamabad 44000, Pakistan;
| | - Yu-Juan Zhang
- College of Life Sciences, Chongqing Normal University, Chongqing 401331, China;
| | - Shulin Zhang
- Department of Immunology and Microbiology, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China; (S.Z.)
| | - Amie Jinghua Wei
- Department of Immunology and Microbiology, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China; (S.Z.)
| | - Shaukat Iqbal Malik
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad 44000, Pakistan; (M.T.K.); (S.I.M.)
| | - Dong Qing Wei
- The State Key Laboratory of Microbial Metabolism, College of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China;
- Correspondence: ; Tel.: +86-21-3420-4573
| |
Collapse
|
20
|
Khan MT, Khan A, Rehman AU, Wang Y, Akhtar K, Malik SI, Wei DQ. Structural and free energy landscape of novel mutations in ribosomal protein S1 (rpsA) associated with pyrazinamide resistance. Sci Rep 2019; 9:7482. [PMID: 31097767 PMCID: PMC6522564 DOI: 10.1038/s41598-019-44013-9] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2018] [Accepted: 04/29/2019] [Indexed: 02/04/2023] Open
Abstract
Resistance to key first-line drugs is a major hurdle to achieve the global end tuberculosis (TB) targets. A prodrug, pyrazinamide (PZA) is the only drug, effective in latent TB, recommended in drug resistance and susceptible Mycobacterium tuberculosis (MTB) isolates. The prodrug conversion into active form, pyrazinoic acid (POA), required the activity of pncA gene encoded pyrazinamidase (PZase). Although pncA mutations have been commonly associated with PZA resistance but a small number of resistance cases have been associated with mutationss in RpsA protein. Here in this study a total of 69 PZA resistance isolates have been sequenced for pncA mutations. However, samples that were found PZA resistant but pncA wild type (pncAWT), have been sequenced for rpsA and panD genes mutation. We repeated a drug susceptibility testing according to the WHO guidelines on 18 pncAWT MTB isolates. The rpsA and panD genes were sequenced. Out of total 69 PZA resistant isolates, 51 harbored 36 mutations in pncA gene (GeneBank Accession No. MH46111) while, fifteen different mutations including seven novel, were detected in the fourth S1 domain of RpsA known as C-terminal (MtRpsACTD) end. We did not detect any mutations in panD gene. Among the rpsA mutations, we investigated the molecular mechanism of resistance behind mutations, D342N, D343N, A344P, and I351F, present in the MtRpsACTD through molecular dynamic simulations (MD). WT showed a good drug binding affinity as compared to mutants (MTs), D342N, D343N, A344P, and I351F. Binding pocket volume, stability, and fluctuations have been altered whereas the total energy, protein folding, and geometric shape analysis further explored a significant variation between WT and MTs. In conclusion, mutations in MtRpsACTD might be involved to alter the RpsA activity, resulting in drug resistance. Such molecular mechanism behind resistance may provide a better insight into the resistance mechanism to achieve the global TB control targets.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Abbas Khan
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China
| | - Ashfaq Ur Rehman
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China
| | - Yanjie Wang
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China
| | - Khalid Akhtar
- National University of Science and Technology, Islamabad, Pakistan
| | - Shaukat Iqbal Malik
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan.
| | - Dong-Qing Wei
- College of Life Sciences and Biotechnology, The State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China.
| |
Collapse
|