1
|
Sankaran S, Krishnan SR, Sayed Y, Gromiha MM. Mechanism of drug resistance in HIV-1 protease subtype C in the presence of Atazanavir. Curr Res Struct Biol 2024; 7:100132. [PMID: 38435053 PMCID: PMC10907180 DOI: 10.1016/j.crstbi.2024.100132] [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/28/2023] [Revised: 02/09/2024] [Accepted: 02/13/2024] [Indexed: 03/05/2024] Open
Abstract
AIDS is one of the deadliest diseases in the history of humankind caused by HIV. Despite the technological development, curtailing the viral infection inside human host still remains a challenge. Therapies such as HAART uses a combination of drugs to inhibit the viral activity. One of the important targets includes HIV protease and inhibiting its activity will minimize the production of mature structural proteins. However, the genetic diversity and the occurrence of drug resistant mutations adds complexity to effective drug design. In this study, we aimed at understanding the drug binding mechanism of one such subtype, namely subtype C and its insertion variant L38HL. We performed multiple molecular dynamics simulations along with binding free energy analysis of wild-type and L38HL bound to Atazanavir (ATV). From the analysis, we revealed that the insertion alters the hydrogen bond and hydrophobic interaction networks. The alterations in the interaction networks increase flexibility at the hinge-fulcrum interface. Further, the effects of these changes affect flap tip curling. Moreover, the changes in the hinge-fulcrum-cantilever interface alters the concerted motion of the functional regions leading to change in the direction of flap movement thus causing a subtle change in the active site volume. Additionally, formation of intramolecular hydrogen bonds in the ATV docked to L38HL restricted the movement of R1 and R2 groups thereby altering the interactions. Overall, the changes in the flexibility of flap together with the changes in the active site volume and compactness of the ligand provide insights for increased binding affinity of ATV with L38HL.
Collapse
Affiliation(s)
- S.V. Sankaran
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, 600036, India
| | - Sowmya R. Krishnan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, 600036, India
| | - Yasien Sayed
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg, South Africa
| | - M. Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, 600036, India
| |
Collapse
|
2
|
Serna-Arbeláez MS, García-Cárcamo V, Rincón-Tabares DS, Guerra D, Loaiza-Cano V, Martinez-Gutierrez M, Pereañez JA, Pastrana-Restrepo M, Galeano E, Zapata W. In Vitro and In Silico Antiviral Activity of Di-Halogenated Compounds Derived from L-Tyrosine against Human Immunodeficiency Virus 1 (HIV-1). Curr Issues Mol Biol 2023; 45:8173-8200. [PMID: 37886959 PMCID: PMC10605077 DOI: 10.3390/cimb45100516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Revised: 09/28/2023] [Accepted: 09/29/2023] [Indexed: 10/28/2023] Open
Abstract
HIV-1 infection is considered one of the major public health problems worldwide. Due to the limited access to antiretroviral therapy, the associated side effects, and the resistance that the virus can generate, it has become necessary to continue the development of new antiviral agents. The study aimed to identify potential antiviral agents for HIV-1 by evaluating the in vitro and in silico activity of 16 synthetic di-halogenated compounds derived from L-Tyrosine. The compounds were tested for cytotoxicity, which was determined using MTT, and a combined antiviral screening strategy (pre- and post-infection treatment) was performed against R5 and X4 strains of HIV-1. The most promising compounds were evaluated against a pseudotyped virus (HIV-GFP-VSV-G), and the effectiveness of these compounds was measured through GFP flow cytometry. Also, the antiviral effect of these compounds was evaluated in PBMCs using flow cytometry and ELISA for p24. The TODB-2M, TODC-2M, TODC-3M, and YDC-3M compounds showed low toxicity and significant inhibitory activity against HIV-1. In silico docking and molecular dynamics assays suggest that the compounds' antiviral activity may be due to interaction with reverse transcriptase, viral protease, or envelope gp120.
Collapse
Affiliation(s)
- Maria S. Serna-Arbeláez
- Grupo Infettare, Facultad de Medicina, Universidad Cooperativa de Colombia, Medellín 050001, Colombia; (M.S.S.-A.); (V.G.-C.)
- Grupo Inmunovirología, Facultad de Medicina, Universidad de Antioquia UdeA, Medellín 050001, Colombia;
| | - Valentina García-Cárcamo
- Grupo Infettare, Facultad de Medicina, Universidad Cooperativa de Colombia, Medellín 050001, Colombia; (M.S.S.-A.); (V.G.-C.)
- Grupo Inmunovirología, Facultad de Medicina, Universidad de Antioquia UdeA, Medellín 050001, Colombia;
| | - Daniel S. Rincón-Tabares
- Grupo Inmunovirología, Facultad de Medicina, Universidad de Antioquia UdeA, Medellín 050001, Colombia;
| | - Diego Guerra
- Instituto de Parasitología y Biomedicina “López-Neyra”, Consejo Superior de Investigaciones Científicas, PTS Granada, 18016 Granada, Spain;
- Programa de Estudio y Control de Enfermedades Tropicales PECET, Faculty of Medicine, University of Antioquia, Medellín 050010, Colombia
| | - Vanessa Loaiza-Cano
- Grupo de Investigación en Ciencias Animales-GRICA, Facultad de Medicina Veterinaria y Zootecnia, Universidad Cooperativa de Colombia, Bucaramanga 680005, Colombia; (V.L.-C.); (M.M.-G.)
| | - Marlen Martinez-Gutierrez
- Grupo de Investigación en Ciencias Animales-GRICA, Facultad de Medicina Veterinaria y Zootecnia, Universidad Cooperativa de Colombia, Bucaramanga 680005, Colombia; (V.L.-C.); (M.M.-G.)
| | - Jaime A. Pereañez
- Grupo Toxinología, Alternativas Terapéuticas y Alimentarias, Facultad de Ciencias Farmacéuticas y Alimentarias, Universidad de Antioquia UdeA, Medellín 050001, Colombia;
| | - Manuel Pastrana-Restrepo
- Productos Naturales Marinos, Departamento de Farmacia, Facultad de Ciencias Farmacéuticas y Alimentarias, Universidad de Antioquia UdeA, Medellín 050001, Colombia; (M.P.-R.); (E.G.)
| | - Elkin Galeano
- Productos Naturales Marinos, Departamento de Farmacia, Facultad de Ciencias Farmacéuticas y Alimentarias, Universidad de Antioquia UdeA, Medellín 050001, Colombia; (M.P.-R.); (E.G.)
| | - Wildeman Zapata
- Grupo Infettare, Facultad de Medicina, Universidad Cooperativa de Colombia, Medellín 050001, Colombia; (M.S.S.-A.); (V.G.-C.)
- Grupo Inmunovirología, Facultad de Medicina, Universidad de Antioquia UdeA, Medellín 050001, Colombia;
| |
Collapse
|
3
|
Venkatachalam S, Murlidharan N, Krishnan SR, Ramakrishnan C, Setshedi M, Pandian R, Barh D, Tiwari S, Azevedo V, Sayed Y, Gromiha MM. Understanding Drug Resistance of Wild-Type and L38HL Insertion Mutant of HIV-1 C Protease to Saquinavir. Genes (Basel) 2023; 14:533. [PMID: 36833460 PMCID: PMC9957153 DOI: 10.3390/genes14020533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2023] [Revised: 02/16/2023] [Accepted: 02/17/2023] [Indexed: 02/25/2023] Open
Abstract
Acquired immunodeficiency syndrome (AIDS) is one of the most challenging infectious diseases to treat on a global scale. Understanding the mechanisms underlying the development of drug resistance is necessary for novel therapeutics. HIV subtype C is known to harbor mutations at critical positions of HIV aspartic protease compared to HIV subtype B, which affects the binding affinity. Recently, a novel double-insertion mutation at codon 38 (L38HL) was characterized in HIV subtype C protease, whose effects on the interaction with protease inhibitors are hitherto unknown. In this study, the potential of L38HL double-insertion in HIV subtype C protease to induce a drug resistance phenotype towards the protease inhibitor, Saquinavir (SQV), was probed using various computational techniques, such as molecular dynamics simulations, binding free energy calculations, local conformational changes and principal component analysis. The results indicate that the L38HL mutation exhibits an increase in flexibility at the hinge and flap regions with a decrease in the binding affinity of SQV in comparison with wild-type HIV protease C. Further, we observed a wide opening at the binding site in the L38HL variant due to an alteration in flap dynamics, leading to a decrease in interactions with the binding site of the mutant protease. It is supported by an altered direction of motion of flap residues in the L38HL variant compared with the wild-type. These results provide deep insights into understanding the potential drug resistance phenotype in infected individuals.
Collapse
Affiliation(s)
- Sankaran Venkatachalam
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Nisha Murlidharan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sowmya R. Krishnan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - C. Ramakrishnan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Mpho Setshedi
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg 2000, South Africa
| | - Ramesh Pandian
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg 2000, South Africa
| | - Debmalya Barh
- Department of Genetics, Ecology, and Evolution, Institute of Biological Sciences, Federal University of Minas Gerais (UFMG), Belo Horizonte 31270-901, Brazil
- Institute of Integrative Omics and Applied Biotechnology (IIOAB), Nonakuri, Purba Medinipur 721172, West Bengal, India
| | - Sandeep Tiwari
- Department of Genetics, Ecology, and Evolution, Institute of Biological Sciences, Federal University of Minas Gerais (UFMG), Belo Horizonte 31270-901, Brazil
- Institute of Biology, Federal University of Bahia, Salvador, BA 40110-909, Brazil
- Institute of Health Sciences, Federal University of Bahia, Salvador, BA 40110-909, Brazil
| | - Vasco Azevedo
- Department of Genetics, Ecology, and Evolution, Institute of Biological Sciences, Federal University of Minas Gerais (UFMG), Belo Horizonte 31270-901, Brazil
| | - Yasien Sayed
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg 2000, South Africa
| | - M. Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
4
|
Bao HY, Wang W, Sun HB, Chen JZ. Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2023; 34:65-89. [PMID: 36762439 DOI: 10.1080/1062936x.2023.2165542] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 12/31/2022] [Indexed: 06/18/2023]
Abstract
Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.
Collapse
Affiliation(s)
- H Y Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - J Z Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
5
|
Yu YX, Wang W, Sun HB, Zhang LL, Wang LF, Yin YY. Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2022; 33:805-831. [PMID: 36322686 DOI: 10.1080/1062936x.2022.2140708] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Accepted: 10/20/2022] [Indexed: 06/16/2023]
Abstract
Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L F Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Y Y Yin
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
6
|
Ritonavir and xk263 Binding-Unbinding with HIV-1 Protease: Pathways, Energy and Comparison. LIFE (BASEL, SWITZERLAND) 2022; 12:life12010116. [PMID: 35054509 PMCID: PMC8779838 DOI: 10.3390/life12010116] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 01/04/2022] [Accepted: 01/10/2022] [Indexed: 01/22/2023]
Abstract
Understanding non-covalent biomolecular recognition, which includes drug-protein bound states and their binding/unbinding processes, is of fundamental importance in chemistry, biology, and medicine. Fully revealing the factors that govern the binding/unbinding processes can further assist in designing drugs with desired binding kinetics. HIV protease (HIVp) plays an integral role in the HIV life cycle, so it is a prime target for drug therapy. HIVp has flexible flaps, and the binding pocket can be accessible by a ligand via various pathways. Comparing ligand association and dissociation pathways can help elucidate the ligand-protein interactions such as key residues directly involved in the interaction or specific protein conformations that determine the binding of a ligand under certain pathway(s). Here, we investigated the ligand unbinding process for a slow binder, ritonavir, and a fast binder, xk263, by using unbiased all-atom accelerated molecular dynamics (aMD) simulation with a re-seeding approach and an explicit solvent model. Using ritonavir-HIVp and xk263-HIVp ligand-protein systems as cases, we sampled multiple unbinding pathways for each ligand and observed that the two ligands preferred the same unbinding route. However, ritonavir required a greater HIVp motion to dissociate as compared with xk263, which can leave the binding pocket with little conformational change of HIVp. We also observed that ritonavir unbinding pathways involved residues which are associated with drug resistance and are distal from catalytic site. Analyzing HIVp conformations sampled during both ligand-protein binding and unbinding processes revealed significantly more overlapping HIVp conformations for ritonavir-HIVp rather than xk263-HIVp. However, many HIVp conformations are unique in xk263-HIVp unbinding processes. The findings are consistent with previous findings that xk263 prefers an induced-fit model for binding and unbinding, whereas ritonavir favors a conformation selection model. This study deepens our understanding of the dynamic process of ligand unbinding and provides insights into ligand-protein recognition mechanisms and drug discovery.
Collapse
|