1
|
Handa Y, Okuwaki K, Kawashima Y, Hatada R, Mochizuki Y, Komeiji Y, Tanaka S, Furuishi T, Yonemochi E, Honma T, Fukuzawa K. Prediction of Binding Pose and Affinity of Nelfinavir, a SARS-CoV-2 Main Protease Repositioned Drug, by Combining Docking, Molecular Dynamics, and Fragment Molecular Orbital Calculations. J Phys Chem B 2024; 128:2249-2265. [PMID: 38437183 PMCID: PMC10946393 DOI: 10.1021/acs.jpcb.3c05564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 02/09/2024] [Accepted: 02/14/2024] [Indexed: 03/06/2024]
Abstract
A novel in silico drug design procedure is described targeting the Main protease (Mpro) of the SARS-CoV-2 virus. The procedure combines molecular docking, molecular dynamics (MD), and fragment molecular orbital (FMO) calculations. The binding structure and properties of Mpro were predicted for Nelfinavir (NFV), which had been identified as a candidate compound through drug repositioning, targeting Mpro. Several poses of the Mpro and NFV complexes were generated by docking, from which four docking poses were selected by scoring with FMO energy. Then, each pose was subjected to MD simulation, 100 snapshot structures were sampled from each of the generated MD trajectories, and the structures were evaluated by FMO calculations to rank the pose based on binding energy. Several residues were found to be important in ligand recognition, including Glu47, Asp48, Glu166, Asp187, and Gln189, all of which interacted strongly with NFV. Asn142 is presumably regarded to form hydrogen bonds or CH/π interaction with NFV; however, in the present calculation, their interactions were transient. Moreover, the tert-butyl group of NFV had no interaction with Mpro. Identifying such strong and weak interactions provides candidates for maintaining and substituting ligand functional groups and important suggestions for drug discovery using drug repositioning. Besides the interaction between NFV and the amino acid residues of Mpro, the desolvation effect of the binding pocket also affected the ranking order. A similar procedure of drug design was applied to Lopinavir, and the calculated interaction energy and experimental inhibitory activity value trends were consistent. Our approach provides a new guideline for structure-based drug design starting from a candidate compound whose complex crystal structure has not been obtained.
Collapse
Affiliation(s)
- Yuma Handa
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Koji Okuwaki
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Yusuke Kawashima
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Ryo Hatada
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Yuji Mochizuki
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Institute
of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Yuto Komeiji
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Health
and Medical Research Institute, AIST, Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
- RIKEN
Center
for Biosystems Dynamics Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Shigenori Tanaka
- Graduate
School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Takayuki Furuishi
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Etsuo Yonemochi
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Teruki Honma
- RIKEN
Center
for Biosystems Dynamics Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Kaori Fukuzawa
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
- Department
of Biomolecular Engineering, Graduate School of Engineering, Tohoku University, 6-6-11 Aoba, Aramaki, Aoba-ku, Sendai 980-8579, Japan
| |
Collapse
|
2
|
Ren P, Li S, Wang S, Zhang X, Bai F. Computer-Aided Prediction of the Interactions of Viral Proteases with Antiviral Drugs: Antiviral Potential of Broad-Spectrum Drugs. Molecules 2023; 29:225. [PMID: 38202808 PMCID: PMC10780089 DOI: 10.3390/molecules29010225] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 12/27/2023] [Accepted: 12/29/2023] [Indexed: 01/12/2024] Open
Abstract
Human society is facing the threat of various viruses. Proteases are promising targets for the treatment of viral infections. In this study, we collected and profiled 170 protease sequences from 125 viruses that infect humans. Approximately 73 of them are viral 3-chymotrypsin-like proteases (3CLpro), and 11 are pepsin-like aspartic proteases (PAPs). Their sequences, structures, and substrate characteristics were carefully analyzed to identify their conserved nature for proposing a pan-3CLpro or pan-PAPs inhibitor design strategy. To achieve this, we used computational prediction and modeling methods to predict the binding complex structures for those 73 3CLpro with 4 protease inhibitors of SARS-CoV-2 and 11 protease inhibitors of HCV. Similarly, the complex structures for the 11 viral PAPs with 9 protease inhibitors of HIV were also obtained. The binding affinities between these compounds and proteins were also evaluated to assess their pan-protease inhibition via MM-GBSA. Based on the drugs targeting viral 3CLpro and PAPs, repositioning of the active compounds identified several potential uses for these drug molecules. As a result, Compounds 1-2, modified based on the structures of Ray1216 and Asunaprevir, indicate potential inhibition of DENV protease according to our computational simulation results. These studies offer ideas and insights for future research in the design of broad-spectrum antiviral drugs.
Collapse
Affiliation(s)
- Pengxuan Ren
- School of Life Science and Technology, Shanghai Institute for Advanced Immunochemical Studies, ShanghaiTech University, Shanghai 201210, China; (P.R.); (S.L.); (S.W.)
| | - Shiwei Li
- School of Life Science and Technology, Shanghai Institute for Advanced Immunochemical Studies, ShanghaiTech University, Shanghai 201210, China; (P.R.); (S.L.); (S.W.)
| | - Shihang Wang
- School of Life Science and Technology, Shanghai Institute for Advanced Immunochemical Studies, ShanghaiTech University, Shanghai 201210, China; (P.R.); (S.L.); (S.W.)
| | - Xianglei Zhang
- School of Life Science and Technology, Shanghai Institute for Advanced Immunochemical Studies, ShanghaiTech University, Shanghai 201210, China; (P.R.); (S.L.); (S.W.)
| | - Fang Bai
- School of Life Science and Technology, Shanghai Institute for Advanced Immunochemical Studies, ShanghaiTech University, Shanghai 201210, China; (P.R.); (S.L.); (S.W.)
- School of Information Science and Technology, ShanghaiTech University, Shanghai 201210, China
- Shanghai Clinical Research and Trial Center, Shanghai 201210, China
| |
Collapse
|
3
|
Dakshinamoorthy A, Asmita A, Senapati S. Comprehending the Structure, Dynamics, and Mechanism of Action of Drug-Resistant HIV Protease. ACS OMEGA 2023; 8:9748-9763. [PMID: 36969469 PMCID: PMC10034783 DOI: 10.1021/acsomega.2c08279] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Accepted: 02/23/2023] [Indexed: 06/18/2023]
Abstract
Since the emergence of the Human Immunodeficiency Virus (HIV) in the 1980s, strategies to combat HIV-AIDS are continuously evolving. Among the many tested targets to tackle this virus, its protease enzyme (PR) was proven to be an attractive option that brought about numerous research publications and ten FDA-approved drugs to inhibit the PR activity. However, the drug-induced mutations in the enzyme made these small molecule inhibitors ineffective with prolonged usage. The research on HIV PR, therefore, remains a thrust area even today. Through this review, we reiterate the importance of understanding the various structural and functional components of HIV PR in redesigning the structure-based small molecule inhibitors. We also discuss at length the currently available FDA-approved drugs and how these drug molecules induced mutations in the enzyme structure. We then recapitulate the reported mechanisms on how these drug-resistant variants remain sufficiently active to cleave the natural substrates. We end with the future scope covering the recently proposed strategies that show promise to deal with the mutations.
Collapse
|
4
|
Marie V, Gordon M. Understanding the co-evolutionary molecular mechanisms of resistance in the HIV-1 Gag and protease. J Biomol Struct Dyn 2022; 40:10852-10861. [PMID: 34253143 DOI: 10.1080/07391102.2021.1950569] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Due to high human immunodeficiency virus type 1 (HIV-1) subtype C infections coupled with increasing antiretroviral treatment failure, the elucidation of complex drug resistance mutational patterns arising through protein co-evolution is required. Despite the inclusion of potent protease inhibitors Lopinavir (LPV) and Darunavir (DRV) in second- and third-line therapies, many patients still fail treatment due to the accumulation of mutations in protease (PR) and recently, Gag. To understand the co-evolutionary molecular mechanisms of resistance in the HIV-1 PR and Gag, we performed 100 ns molecular dynamic simulations on multidrug resistant PR's when bound to LPV, DRV or a mutated A431V NC|p1 Gag cleavage site (CS). Here we showed that distinct changes in PR's active site, flap and elbow regions due to several PR resistance mutations (L10F, M46I, I54V, L76V, V82A) were found to alter LPV and DRV drug binding. However, binding was significantly exacerbated when the mutant PRs were bound to the NC|p1 Gag CS. Although A431V was shown to coordinate several residues in PR, the L76V PR mutation was found to have a significant role in substrate recognition. Consequently, a greater binding affinity was observed when the mutated substrate was bound to an L76V-inclusive PR mutant (Gbind: -62.46 ± 5.75 kcal/mol) than without (Gbind: -50.34 ± 6.28 kcal/mol). These data showed that the co-selection of resistance mutations in the enzyme and substrate can simultaneously constrict regions in PR's active site whilst flexing the flaps to allow flexible movement of the substrate and multiple, complex mechanisms of resistance to occur. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Veronna Marie
- KwaZulu-Natal Research Innovation & Sequencing Platform (KRISP), Nelson R Mandela School of Medicine, University of KwaZulu-Natal, DurbanSouth Africa
| | - Michelle Gordon
- KwaZulu-Natal Research Innovation & Sequencing Platform (KRISP), Nelson R Mandela School of Medicine, University of KwaZulu-Natal, DurbanSouth Africa
| |
Collapse
|
5
|
Wong-Sam A, Wang YF, Kneller DW, Kovalevsky AY, Ghosh AK, Harrison RW, Weber IT. HIV-1 protease with 10 lopinavir and darunavir resistance mutations exhibits altered inhibition, structural rearrangements and extreme dynamics. J Mol Graph Model 2022; 117:108315. [PMID: 36108568 PMCID: PMC10091457 DOI: 10.1016/j.jmgm.2022.108315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2022] [Revised: 08/17/2022] [Accepted: 08/21/2022] [Indexed: 01/14/2023]
Abstract
Antiretroviral drug resistance is a therapeutic obstacle for people with HIV. HIV protease inhibitors darunavir and lopinavir are recommended for resistant infections. We characterized a protease mutant (PR10x) derived from a highly resistant clinical isolate including 10 mutations associated with resistance to lopinavir and darunavir. Compared to the wild-type protease, PR10x exhibits ∼3-fold decrease in catalytic efficiency and Ki values of 2-3 orders of magnitude worse for darunavir, lopinavir, and potent investigational inhibitor GRL-519. Crystal structures of the mutant were solved in a ligand-free form and in complex with GRL-519. The structures show altered interactions in the active site, flap-core interface, hydrophobic core, hinge region, and 80s loop compared to the corresponding wild-type protease structures. The ligand-free crystal structure exhibits a highly curled flap conformation which may amplify drug resistance. Molecular dynamics simulations performed for 1 μs on ligand-free dimers showed extremely large fluctuations in the flaps for PR10x compared to equivalent simulations on PR with a single L76V mutation or wild-type protease. This analysis offers insight about the synergistic effects of mutations in highly resistant variants.
Collapse
Affiliation(s)
- Andres Wong-Sam
- Department of Biology, Georgia State University, Atlanta, GA, 30303, USA
| | - Yuan-Fang Wang
- Department of Biology, Georgia State University, Atlanta, GA, 30303, USA
| | - Daniel W Kneller
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Andrey Y Kovalevsky
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Arun K Ghosh
- Department of Chemistry and Department of Medicinal Chemistry, Purdue University, West Lafayette, IN, 47907, USA
| | - Robert W Harrison
- Department of Biology, Georgia State University, Atlanta, GA, 30303, USA; Department of Computer Science, Georgia State University, Atlanta, GA, 30303, USA
| | - Irene T Weber
- Department of Biology, Georgia State University, Atlanta, GA, 30303, USA; Department of Chemistry, Georgia State University, Atlanta, GA, 30303, USA.
| |
Collapse
|
6
|
Truong TT, Hayn M, Frich CK, Olari L, Ladefoged LK, Jarlstad Olesen MT, Jakobsen JH, Lunabjerg‐Vestergaard CK, Schiøtt B, Münch J, Zelikin AN. Potentiation of Drug Toxicity Through Virus Latency Reversal Promotes Preferential Elimination of HIV Infected Cells. ADVANCED THERAPEUTICS 2022. [DOI: 10.1002/adtp.202200113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Thanh Tung Truong
- Department of Chemistry Aarhus University Langelandsgade 140 Aarhus C 8000 Denmark
| | - Manuel Hayn
- Institute of Molecular Virology Ulm University Medical Center 89081 Ulm Germany
| | - Camilla Kaas Frich
- Department of Chemistry Aarhus University Langelandsgade 140 Aarhus C 8000 Denmark
| | - Lia‐Raluca Olari
- Institute of Molecular Virology Ulm University Medical Center 89081 Ulm Germany
| | | | | | - Josefine H. Jakobsen
- Department of Chemistry Aarhus University Langelandsgade 140 Aarhus C 8000 Denmark
| | | | - Birgit Schiøtt
- Department of Chemistry Aarhus University Langelandsgade 140 Aarhus C 8000 Denmark
- iNano Interdisciplinary Nanoscience Centre Aarhus University Aarhus 8000 Denmark
| | - Jan Münch
- Institute of Molecular Virology Ulm University Medical Center 89081 Ulm Germany
- iNano Interdisciplinary Nanoscience Centre Aarhus University Aarhus 8000 Denmark
| | - Alexander N. Zelikin
- Department of Chemistry Aarhus University Langelandsgade 140 Aarhus C 8000 Denmark
- iNano Interdisciplinary Nanoscience Centre Aarhus University Aarhus 8000 Denmark
| |
Collapse
|
7
|
Ahmad B, Batool M, Ain QU, Kim MS, Choi S. Exploring the Binding Mechanism of PF-07321332 SARS-CoV-2 Protease Inhibitor through Molecular Dynamics and Binding Free Energy Simulations. Int J Mol Sci 2021; 22:9124. [PMID: 34502033 PMCID: PMC8430524 DOI: 10.3390/ijms22179124] [Citation(s) in RCA: 44] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2021] [Revised: 08/20/2021] [Accepted: 08/21/2021] [Indexed: 12/23/2022] Open
Abstract
The novel coronavirus disease, caused by severe acute respiratory coronavirus 2 (SARS-CoV-2), rapidly spreading around the world, poses a major threat to the global public health. Herein, we demonstrated the binding mechanism of PF-07321332, α-ketoamide, lopinavir, and ritonavir to the coronavirus 3-chymotrypsin-like-protease (3CLpro) by means of docking and molecular dynamic (MD) simulations. The analysis of MD trajectories of 3CLpro with PF-07321332, α-ketoamide, lopinavir, and ritonavir revealed that 3CLpro-PF-07321332 and 3CLpro-α-ketoamide complexes remained stable compared with 3CLpro-ritonavir and 3CLpro-lopinavir. Investigating the dynamic behavior of ligand-protein interaction, ligands PF-07321332 and α-ketoamide showed stronger bonding via making interactions with catalytic dyad residues His41-Cys145 of 3CLpro. Lopinavir and ritonavir were unable to disrupt the catalytic dyad, as illustrated by increased bond length during the MD simulation. To decipher the ligand binding mode and affinity, ligand interactions with SARS-CoV-2 proteases and binding energy were calculated. The binding energy of the bespoke antiviral PF-07321332 clinical candidate was two times higher than that of α-ketoamide and three times than that of lopinavir and ritonavir. Our study elucidated in detail the binding mechanism of the potent PF-07321332 to 3CLpro along with the low potency of lopinavir and ritonavir due to weak binding affinity demonstrated by the binding energy data. This study will be helpful for the development and optimization of more specific compounds to combat coronavirus disease.
Collapse
Affiliation(s)
- Bilal Ahmad
- Department of Molecular Science and Technology, Ajou University, Suwon 16499, Korea; (B.A.); (M.B.); (Q.u.A.); (M.S.K.)
| | - Maria Batool
- Department of Molecular Science and Technology, Ajou University, Suwon 16499, Korea; (B.A.); (M.B.); (Q.u.A.); (M.S.K.)
- S&K Therapeutics, Campus Plaza 418, Ajou University, Suwon 16502, Korea
| | - Qurat ul Ain
- Department of Molecular Science and Technology, Ajou University, Suwon 16499, Korea; (B.A.); (M.B.); (Q.u.A.); (M.S.K.)
| | - Moon Suk Kim
- Department of Molecular Science and Technology, Ajou University, Suwon 16499, Korea; (B.A.); (M.B.); (Q.u.A.); (M.S.K.)
| | - Sangdun Choi
- Department of Molecular Science and Technology, Ajou University, Suwon 16499, Korea; (B.A.); (M.B.); (Q.u.A.); (M.S.K.)
- S&K Therapeutics, Campus Plaza 418, Ajou University, Suwon 16502, Korea
| |
Collapse
|
8
|
Ota R, So K, Tsuda M, Higuchi Y, Yamashita F. Prediction of HIV drug resistance based on the 3D protein structure: Proposal of molecular field mapping. PLoS One 2021; 16:e0255693. [PMID: 34347839 PMCID: PMC8336827 DOI: 10.1371/journal.pone.0255693] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2021] [Accepted: 07/21/2021] [Indexed: 11/19/2022] Open
Abstract
A method for predicting HIV drug resistance by using genotypes would greatly assist in selecting appropriate combinations of antiviral drugs. Models reported previously have had two major problems: lack of information on the 3D protein structure and processing of incomplete sequencing data in the modeling procedure. We propose obtaining the 3D structural information of viral proteins by using homology modeling and molecular field mapping, instead of just their primary amino acid sequences. The molecular field potential parameters reflect the physicochemical characteristics associated with the 3D structure of the proteins. We also introduce the Bayesian conditional mutual information theory to estimate the probabilities of occurrence of all possible protein candidates from an incomplete sequencing sample. This approach allows for the effective use of uncertain information for the modeling process. We applied these data analysis techniques to the HIV-1 protease inhibitor dataset and developed drug resistance prediction models with reasonable performance.
Collapse
Affiliation(s)
- Ryosaku Ota
- Department of Drug Delivery Research, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Kanako So
- Department of Applied Pharmaceutics and Pharmacokinetics, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Masahiro Tsuda
- Department of Applied Pharmaceutics and Pharmacokinetics, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Yuriko Higuchi
- Department of Drug Delivery Research, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Fumiyoshi Yamashita
- Department of Drug Delivery Research, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
- Department of Applied Pharmaceutics and Pharmacokinetics, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
- * E-mail:
| |
Collapse
|
9
|
Eche S, Kumar A, Sonela N, Gordon ML. Acquired HIV-1 Protease Conformational Flexibility Associated with Lopinavir Failure May Shape the Outcome of Darunavir Therapy after Antiretroviral Therapy Switch. Biomolecules 2021; 11:489. [PMID: 33805099 PMCID: PMC8064090 DOI: 10.3390/biom11040489] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Revised: 03/08/2021] [Accepted: 03/12/2021] [Indexed: 11/16/2022] Open
Abstract
Understanding the underlying molecular interaction during a therapy switch from lopinavir (LPV) to darunavir (DRV) is essential to achieve long-term virological suppression. We investigated the kinetic and structural characteristics of multidrug-resistant South African HIV-1 subtype C protease (HIV-1 PR) during therapy switch from LPV to DRV using enzyme activity and inhibition assay, fluorescence spectroscopy, and molecular dynamic simulation. The HIV-1 protease variants were from clinical isolates with a combination of drug resistance mutations; MUT-1 (M46I, I54V, V82A, and L10F), MUT-2 (M46I, I54V, L76V, V82A, L10F, and L33F), and MUT-3 (M46I, I54V, L76V, V82A, L90M, and F53L). Enzyme kinetics analysis shows an association between increased relative resistance to LPV and DRV with the progressive decrease in the mutant HIV-1 PR variants' catalytic efficiency. A direct relationship between high-level resistance to LPV and intermediate resistance to DRV with intrinsic changes in the three-dimensional structure of the mutant HIV-1 PR as a function of the multidrug-resistance mutation was observed. In silico analysis attributed these structural adjustments to the multidrug-resistance mutations affecting the LPV and DRV binding landscape. Though DRV showed superiority to LPV, as a lower concentration was needed to inhibit the HIV-1 PR variants, the inherent structural changes resulting from mutations selected during LPV therapy may dynamically shape the DRV treatment outcome after the therapy switch.
Collapse
Affiliation(s)
- Simeon Eche
- Discipline of Virology, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban 4001, South Africa;
| | - Ajit Kumar
- Discipline of Microbiology, School of Life Sciences, University of KwaZulu-Natal (Westville Campus), Durban 4000, South Africa;
| | - Nelson Sonela
- School of Medicine, Physical and Natural Sciences, University of Rome Tor Vegata, 1-00133 Rome, Italy;
- Chantal Biya International Reference Center for Research on the Management and Prevention of HIV/AIDS (CIRCB), Yaoundé P.O. Box 3077, Cameroon
| | - Michelle L. Gordon
- Discipline of Virology, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban 4001, South Africa;
| |
Collapse
|
10
|
Affiliation(s)
- Mark Aldren M. Feliciano
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Brian Gold
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
11
|
Darunavir-Resistant HIV-1 Protease Constructs Uphold a Conformational Selection Hypothesis for Drug Resistance. Viruses 2020; 12:v12111275. [PMID: 33171603 PMCID: PMC7695139 DOI: 10.3390/v12111275] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2020] [Revised: 10/23/2020] [Accepted: 10/30/2020] [Indexed: 01/24/2023] Open
Abstract
Multidrug resistance continues to be a barrier to the effectiveness of highly active antiretroviral therapy in the treatment of human immunodeficiency virus 1 (HIV-1) infection. Darunavir (DRV) is a highly potent protease inhibitor (PI) that is oftentimes effective when drug resistance has emerged against first-generation inhibitors. Resistance to darunavir does evolve and requires 10–20 amino acid substitutions. The conformational landscapes of six highly characterized HIV-1 protease (PR) constructs that harbor up to 19 DRV-associated mutations were characterized by distance measurements with pulsed electron double resonance (PELDOR) paramagnetic resonance spectroscopy, namely double electron–electron resonance (DEER). The results show that the accumulated substitutions alter the conformational landscape compared to PI-naïve protease where the semi-open conformation is destabilized as the dominant population with open-like states becoming prevalent in many cases. A linear correlation is found between values of the DRV inhibition parameter Ki and the open-like to closed-state population ratio determined from DEER. The nearly 50% decrease in occupancy of the semi-open conformation is associated with reduced enzymatic activity, characterized previously in the literature.
Collapse
|
12
|
Kneller DW, Agniswamy J, Harrison RW, Weber IT. Highly drug-resistant HIV-1 protease reveals decreased intra-subunit interactions due to clusters of mutations. FEBS J 2020; 287:3235-3254. [PMID: 31920003 PMCID: PMC7343616 DOI: 10.1111/febs.15207] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2019] [Revised: 11/16/2019] [Accepted: 01/08/2020] [Indexed: 01/07/2023]
Abstract
Drug-resistance is a serious problem for treatment of the HIV/AIDS pandemic. Potent clinical inhibitors of HIV-1 protease show several orders of magnitude worse inhibition of highly drug-resistant variants. Hence, the structure and enzyme activities were analyzed for HIV protease mutant HIV-1 protease (EC 3.4.23.16) (PR) with 22 mutations (PRS5B) from a clinical isolate that was selected by machine learning to represent high-level drug-resistance. PRS5B has 22 mutations including only one (I84V) in the inhibitor binding site; however, clinical inhibitors had poor inhibition of PRS5B activity with kinetic inhibition value (Ki ) values of 4-1000 nm or 18- to 8000-fold worse than for wild-type PR. High-resolution crystal structures of PRS5B complexes with the best inhibitors, amprenavir (APV) and darunavir (DRV) (Ki ~ 4 nm), revealed only minor changes in protease-inhibitor interactions. Instead, two distinct clusters of mutations in distal regions induce coordinated conformational changes that decrease favorable internal interactions across the entire protein subunit. The largest structural rearrangements are described and compared to other characterized resistant mutants. In the protease hinge region, the N83D mutation eliminates a hydrogen bond connecting the hinge and core of the protease and increases disorder compared to highly resistant mutants PR with 17 mutations and PR with 20 mutations with similar hinge mutations. In a distal β-sheet, mutations G73T and A71V coordinate with accessory mutations to bring about shifts that propagate throughout the subunit. Molecular dynamics simulations of ligand-free dimers show differences consistent with loss of interactions in mutant compared to wild-type PR. Clusters of mutations exhibit both coordinated and antagonistic effects, suggesting PRS5B may represent an intermediate stage in the evolution of more highly resistant variants. DATABASES: Structural data are available in Protein Data Bank under the accession codes 6P9A and 6P9B for PRS5B/DRV and PRS5B/APV, respectively.
Collapse
Affiliation(s)
- Daniel W. Kneller
- Department of Biology, Georgia State University, Atlanta, Georgia 30303, United States of America
| | - Johnson Agniswamy
- Department of Biology, Georgia State University, Atlanta, Georgia 30303, United States of America
| | - Robert W. Harrison
- Department of Computer Science, Georgia State University, Atlanta, Georgia 30303, United States of America
| | - Irene T. Weber
- Department of Biology, Georgia State University, Atlanta, Georgia 30303, United States of America,Department of Chemistry, Georgia State University, Atlanta, Georgia 30303, United States of America,Author of correspondence:
| |
Collapse
|
13
|
Marie V, Gordon M. Gag-protease coevolution shapes the outcome of lopinavir-inclusive treatment regimens in chronically infected HIV-1 subtype C patients. Bioinformatics 2020; 35:3219-3223. [PMID: 30753326 DOI: 10.1093/bioinformatics/btz076] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2018] [Revised: 01/03/2019] [Accepted: 02/11/2019] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Commonly, protease inhibitor failure is characterized by the development of multiple protease resistance mutations (PRMs). While the impact of PRMs on therapy failure are understood, the introduction of Gag mutations with protease remains largely unclear. RESULTS Here, we utilized phylogenetic analyses and Bayesian network learning as tools to understand Gag-protease coevolution and elucidate the pathways leading to Lopinavir failure in HIV-1 subtype C infected patients. Our analyses indicate that while PRMs coevolve in response to drug selection pressure within protease, the Gag mutations added to the existing network while specifically interacting with known Lopinavir failure PRMs. Additionally, the selection of mutations at specific positions in Gag-protease suggests that these coevolving mutational changes occurs to maintain structural integrity during Gag cleavage. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- V Marie
- KwaZulu-Natal Research Innovation and Sequencing Platform, University of KwaZulu-Natal, Durban, South Africa
| | - M Gordon
- KwaZulu-Natal Research Innovation and Sequencing Platform, University of KwaZulu-Natal, Durban, South Africa
| |
Collapse
|
14
|
Bastys T, Gapsys V, Walter H, Heger E, Doncheva NT, Kaiser R, de Groot BL, Kalinina OV. Non-active site mutants of HIV-1 protease influence resistance and sensitisation towards protease inhibitors. Retrovirology 2020; 17:13. [PMID: 32430025 PMCID: PMC7236880 DOI: 10.1186/s12977-020-00520-6] [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: 03/06/2020] [Accepted: 05/04/2020] [Indexed: 02/07/2023] Open
Abstract
Background HIV-1 can develop resistance to antiretroviral drugs, mainly through mutations within the target regions of the drugs. In HIV-1 protease, a majority of resistance-associated mutations that develop in response to therapy with protease inhibitors are found in the protease’s active site that serves also as a binding pocket for the protease inhibitors, thus directly impacting the protease-inhibitor interactions. Some resistance-associated mutations, however, are found in more distant regions, and the exact mechanisms how these mutations affect protease-inhibitor interactions are unclear. Furthermore, some of these mutations, e.g. N88S and L76V, do not only induce resistance to the currently administered drugs, but contrarily induce sensitivity towards other drugs. In this study, mutations N88S and L76V, along with three other resistance-associated mutations, M46I, I50L, and I84V, are analysed by means of molecular dynamics simulations to investigate their role in complexes of the protease with different inhibitors and in different background sequence contexts. Results Using these simulations for alchemical calculations to estimate the effects of mutations M46I, I50L, I84V, N88S, and L76V on binding free energies shows they are in general in line with the mutations’ effect on \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$IC_{50}$$\end{document}IC50 values. For the primary mutation L76V, however, the presence of a background mutation M46I in our analysis influences whether the unfavourable effect of L76V on inhibitor binding is sufficient to outweigh the accompanying reduction in catalytic activity of the protease. Finally, we show that L76V and N88S changes the hydrogen bond stability of these residues with residues D30/K45 and D30/T31/T74, respectively. Conclusions We demonstrate that estimating the effect of both binding pocket and distant mutations on inhibitor binding free energy using alchemical calculations can reproduce their effect on the experimentally measured \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$IC_{50}$$\end{document}IC50 values. We show that distant site mutations L76V and N88S affect the hydrogen bond network in the protease’s active site, which offers an explanation for the indirect effect of these mutations on inhibitor binding. This work thus provides valuable insights on interplay between primary and background mutations and mechanisms how they affect inhibitor binding.
Collapse
Affiliation(s)
- Tomas Bastys
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics, 66123, Saarbrücken, Germany.,Saarbrücken Graduate School of Computer Science, University of Saarland, 66123, Saarbrücken, Germany
| | - Vytautas Gapsys
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, 37077, Göttingen, Germany
| | - Hauke Walter
- Medizinisches Labor Stendal, 39576, Stendal, Germany
| | - Eva Heger
- Institute of Virology, University of Cologne, 50935, Cologne, Germany
| | - Nadezhda T Doncheva
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics, 66123, Saarbrücken, Germany.,Faculty of Health and Medical Sciences, University of Copenhagen, 2200, Copenhagen, Denmark
| | - Rolf Kaiser
- Institute of Virology, University of Cologne, 50935, Cologne, Germany
| | - Bert L de Groot
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, 37077, Göttingen, Germany
| | - Olga V Kalinina
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics, 66123, Saarbrücken, Germany. .,Helmholtz Institute for Pharmaceutical Research Saarland (HIPS), Helmholtz Centre for Infection Research (HZI), 66123, Saarbrücken, Germany. .,Faculty of Medicine, Saarland University, 66421, Homburg, Germany.
| |
Collapse
|
15
|
Stefanik M, Valdes JJ, Ezebuo FC, Haviernik J, Uzochukwu IC, Fojtikova M, Salat J, Eyer L, Ruzek D. FDA-Approved Drugs Efavirenz, Tipranavir, and Dasabuvir Inhibit Replication of Multiple Flaviviruses in Vero Cells. Microorganisms 2020; 8:microorganisms8040599. [PMID: 32326119 PMCID: PMC7232190 DOI: 10.3390/microorganisms8040599] [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: 02/10/2020] [Revised: 03/16/2020] [Accepted: 04/15/2020] [Indexed: 02/01/2023] Open
Abstract
Vector-borne flaviviruses (VBFs) affect human health worldwide, but no approved drugs are available specifically to treat VBF-associated infections. Here, we performed in silico screening of a library of U.S. Food and Drug Administration-approved antiviral drugs for their interaction with Zika virus proteins. Twelve hit drugs were identified by the docking experiments and tested in cell-based antiviral assay systems. Efavirenz, tipranavir, and dasabuvir at micromolar concentrations were identified to inhibit all VBFs tested; i.e., two representatives of mosquito-borne flaviviruses (Zika and West Nile viruses) and one representative of flaviviruses transmitted by ticks (tick-borne encephalitis virus). The results warrant further research into these drugs, either individually or in combination, as possible pan-flavivirus inhibitors.
Collapse
Affiliation(s)
- Michal Stefanik
- Department of Virology, Veterinary Research Institute, Hudcova 70, CZ-62100 Brno, Czech Republic; (M.S.); (J.J.V.); (J.H.); (M.F.); (J.S.); (L.E.)
- Department of Chemistry and Biochemistry, Mendel University in Brno, Zemedelska 1, CZ-613 00 Brno, Czech Republic
| | - James J. Valdes
- Department of Virology, Veterinary Research Institute, Hudcova 70, CZ-62100 Brno, Czech Republic; (M.S.); (J.J.V.); (J.H.); (M.F.); (J.S.); (L.E.)
- Institute of Parasitology, Biology Centre of the Czech Academy of Sciences, Branisovska 31, CZ-37005 Ceske Budejovice, Czech Republic
| | - Fortunatus C. Ezebuo
- Department of Pharmaceutical and Medicinal Chemistry, Faculty of Pharmaceutical Sciences, Nnamdi Azikiwe University, PMB 5025 Awka 420281, Nigeria; (F.C.E.); (I.C.U.)
| | - Jan Haviernik
- Department of Virology, Veterinary Research Institute, Hudcova 70, CZ-62100 Brno, Czech Republic; (M.S.); (J.J.V.); (J.H.); (M.F.); (J.S.); (L.E.)
- Faculty of Science, Masaryk University, Kamenice 5, CZ-625 00 Brno, Czech Republic
| | - Ikemefuna C. Uzochukwu
- Department of Pharmaceutical and Medicinal Chemistry, Faculty of Pharmaceutical Sciences, Nnamdi Azikiwe University, PMB 5025 Awka 420281, Nigeria; (F.C.E.); (I.C.U.)
| | - Martina Fojtikova
- Department of Virology, Veterinary Research Institute, Hudcova 70, CZ-62100 Brno, Czech Republic; (M.S.); (J.J.V.); (J.H.); (M.F.); (J.S.); (L.E.)
| | - Jiri Salat
- Department of Virology, Veterinary Research Institute, Hudcova 70, CZ-62100 Brno, Czech Republic; (M.S.); (J.J.V.); (J.H.); (M.F.); (J.S.); (L.E.)
| | - Ludek Eyer
- Department of Virology, Veterinary Research Institute, Hudcova 70, CZ-62100 Brno, Czech Republic; (M.S.); (J.J.V.); (J.H.); (M.F.); (J.S.); (L.E.)
- Institute of Parasitology, Biology Centre of the Czech Academy of Sciences, Branisovska 31, CZ-37005 Ceske Budejovice, Czech Republic
| | - Daniel Ruzek
- Department of Virology, Veterinary Research Institute, Hudcova 70, CZ-62100 Brno, Czech Republic; (M.S.); (J.J.V.); (J.H.); (M.F.); (J.S.); (L.E.)
- Institute of Parasitology, Biology Centre of the Czech Academy of Sciences, Branisovska 31, CZ-37005 Ceske Budejovice, Czech Republic
- Correspondence:
| |
Collapse
|
16
|
Structural studies of antiviral inhibitor with HIV-1 protease bearing drug resistant substitutions of V32I, I47V and V82I. Biochem Biophys Res Commun 2019; 514:974-978. [PMID: 31092330 DOI: 10.1016/j.bbrc.2019.05.064] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2019] [Accepted: 05/07/2019] [Indexed: 12/23/2022]
Abstract
HIV-1 protease inhibitors are effective in HIV/AIDS therapy, although drug resistance is a severe problem. This study examines the effects of four investigational inhibitors against HIV-1 protease with drug resistant mutations of V32I, I47V and V82I (PRTri) that model the inhibitor-binding site of HIV-2 protease. These inhibitors contain diverse chemical modifications on the darunavir scaffold and form new interactions with wild type protease, however, the measured inhibition constants for PRTri mutant range from 17 to 40 nM or significantly worse than picomolar values reported for wild type enzyme. The X-ray crystal structure of PRTri mutant in complex with inhibitor 1 at 1.5 Å resolution shows minor changes in interactions with inhibitor compared with the corresponding wild type PR complex. Instead, the basic amine at P2 of inhibitor together with mutation V82I induces two alternate conformations for the side chain of Arg8 with new interactions with inhibitor and Leu10. Hence, inhibition is influenced by small coordinated changes in hydrophobic interactions.
Collapse
|