1
|
Ye X, Zhang P, Tao J, Wang JCK, Mafi A, Grob NM, Quartararo AJ, Baddock HT, Chan LJG, McAllister FE, Foe I, Loas A, Eaton DL, Hao Q, Nile AH, Pentelute BL. Discovery of reactive peptide inhibitors of human papillomavirus oncoprotein E6. Chem Sci 2023; 14:12484-12497. [PMID: 38020382 PMCID: PMC10646941 DOI: 10.1039/d3sc02782a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2023] [Accepted: 08/22/2023] [Indexed: 12/01/2023] Open
Abstract
Human papillomavirus (HPV) infections account for nearly all cervical cancer cases, which is the fourth most common cancer in women worldwide. High-risk variants, including HPV16, drive tumorigenesis in part by promoting the degradation of the tumor suppressor p53. This degradation is mediated by the HPV early protein 6 (E6), which recruits the E3 ubiquitin ligase E6AP and redirects its activity towards ubiquitinating p53. Targeting the protein interaction interface between HPV E6 and E6AP is a promising modality to mitigate HPV-mediated degradation of p53. In this study, we designed a covalent peptide inhibitor, termed reactide, that mimics the E6AP LXXLL binding motif by selectively targeting cysteine 58 in HPV16 E6 with quantitative conversion. This reactide provides a starting point in the development of covalent peptidomimetic inhibitors for intervention against HPV-driven cancers.
Collapse
Affiliation(s)
- Xiyun Ye
- Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
| | - Peiyuan Zhang
- Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
| | - Jason Tao
- Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
| | - John C K Wang
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Amirhossein Mafi
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Nathalie M Grob
- Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
| | - Anthony J Quartararo
- Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
| | - Hannah T Baddock
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Leanne J G Chan
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Fiona E McAllister
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Ian Foe
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Andrei Loas
- Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
| | - Dan L Eaton
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Qi Hao
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Aaron H Nile
- Calico Life Sciences LLC 1170 Veterans Boulevard South San Francisco CA 94080 USA
| | - Bradley L Pentelute
- Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
- The Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology 500 Main Street Cambridge MA 02142 USA
- Center for Environmental Health Sciences, Massachusetts Institute of Technology 77 Massachusetts Avenue Cambridge MA 02139 USA
- Broad Institute of MIT and Harvard 415 Main Street Cambridge MA 02142 USA
| |
Collapse
|
2
|
Deeks HM, Zinovjev K, Barnoud J, Mulholland AJ, van der Kamp MW, Glowacki DR. Free energy along drug-protein binding pathways interactively sampled in virtual reality. Sci Rep 2023; 13:16665. [PMID: 37794083 PMCID: PMC10551034 DOI: 10.1038/s41598-023-43523-x] [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: 06/30/2023] [Accepted: 09/25/2023] [Indexed: 10/06/2023] Open
Abstract
We describe a two-step approach for combining interactive molecular dynamics in virtual reality (iMD-VR) with free energy (FE) calculation to explore the dynamics of biological processes at the molecular level. We refer to this combined approach as iMD-VR-FE. Stage one involves using a state-of-the-art 'human-in-the-loop' iMD-VR framework to generate a diverse range of protein-ligand unbinding pathways, benefitting from the sophistication of human spatial and chemical intuition. Stage two involves using the iMD-VR-sampled pathways as initial guesses for defining a path-based reaction coordinate from which we can obtain a corresponding free energy profile using FE methods. To investigate the performance of the method, we apply iMD-VR-FE to investigate the unbinding of a benzamidine ligand from a trypsin protein. The binding free energy calculated using iMD-VR-FE is similar for each pathway, indicating internal consistency. Moreover, the resulting free energy profiles can distinguish energetic differences between pathways corresponding to various protein-ligand conformations (e.g., helping to identify pathways that are more favourable) and enable identification of metastable states along the pathways. The two-step iMD-VR-FE approach offers an intuitive way for researchers to test hypotheses for candidate pathways in biomolecular systems, quickly obtaining both qualitative and quantitative insight.
Collapse
Affiliation(s)
- Helen M Deeks
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | - Kirill Zinovjev
- Departamento de Química Física, Universidad de Valencia, 46100, Burjassot, Spain
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK
| | - Jonathan Barnoud
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
- CiTIUS | Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, Rúa de Jenaro de la Fuente, s/n, 15705, Santiago de Compostela, A Coruña, Spain
| | - Adrian J Mulholland
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | - Marc W van der Kamp
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK.
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK.
| | - David R Glowacki
- CiTIUS | Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, Rúa de Jenaro de la Fuente, s/n, 15705, Santiago de Compostela, A Coruña, Spain.
| |
Collapse
|
3
|
Sherry D, Worth R, Ismail ZS, Sayed Y. Cantilever-centric mechanism of cooperative non-active site mutations in HIV protease: Implications for flap dynamics. J Mol Graph Model 2021; 106:107931. [PMID: 34030114 DOI: 10.1016/j.jmgm.2021.107931] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2021] [Revised: 03/26/2021] [Accepted: 04/14/2021] [Indexed: 11/16/2022]
Abstract
The HIV-1 protease is an important drug target in antiretroviral therapy due to the crucial role it plays in viral maturation. A greater understanding of the dynamics of the protease as a result of drug-induced mutations has been successfully elucidated using computational models in the past. We performed induced-fit docking studies and molecular dynamics simulations on the wild-type South African HIV-1 subtype C protease and two non-active site mutation-containing protease variants; HP3 PR and HP4 PR. The HP3 PR contained the I13V, I62V, and V77I mutations while HP4 PR contained the same mutations with the addition of the L33F mutation. The simulations were initiated in a cubic cell universe containing explicit solvent, with the protease variants beginning in the fully closed conformation. The trajectory for each simulation totalled 50 ns. The results indicate that the mutations increase the dynamics of the flap, hinge, fulcrum and cantilever regions when compared to the wild-type protease while in complex with protease inhibitors. Specifically, these mutations result in the protease favouring the semi-open conformation when in complex with inhibitors. Moreover, the HP4 PR adopted curled flap tip conformers which coordinated several water molecules into the active site in a manner that may reduce inhibitor binding affinity. The mutations affected the thermodynamic landscape of inhibitor binding as there were fewer observable chemical contacts between the mutated variants and saquinavir, atazanavir and darunavir. These data help to elucidate the biophysical basis for the selection of cooperative non-active site mutations by the HI virus.
Collapse
Affiliation(s)
- Dean Sherry
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg, 2050, South Africa
| | - Roland Worth
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg, 2050, South Africa
| | - Zaahida Sheik Ismail
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg, 2050, South Africa
| | - Yasien Sayed
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Johannesburg, 2050, South Africa.
| |
Collapse
|
4
|
Deeks HM, Walters RK, Hare SR, O’Connor MB, Mulholland AJ, Glowacki DR. Interactive molecular dynamics in virtual reality for accurate flexible protein-ligand docking. PLoS One 2020; 15:e0228461. [PMID: 32160194 PMCID: PMC7065745 DOI: 10.1371/journal.pone.0228461] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2019] [Accepted: 01/15/2020] [Indexed: 12/11/2022] Open
Abstract
Simulating drug binding and unbinding is a challenge, as the rugged energy landscapes that separate bound and unbound states require extensive sampling that consumes significant computational resources. Here, we describe the use of interactive molecular dynamics in virtual reality (iMD-VR) as an accurate low-cost strategy for flexible protein-ligand docking. We outline an experimental protocol which enables expert iMD-VR users to guide ligands into and out of the binding pockets of trypsin, neuraminidase, and HIV-1 protease, and recreate their respective crystallographic protein-ligand binding poses within 5-10 minutes. Following a brief training phase, our studies shown that iMD-VR novices were able to generate unbinding and rebinding pathways on similar timescales as iMD-VR experts, with the majority able to recover binding poses within 2.15 Å RMSD of the crystallographic binding pose. These results indicate that iMD-VR affords sufficient control for users to carry out the detailed atomic manipulations required to dock flexible ligands into dynamic enzyme active sites and recover crystallographic poses, offering an interesting new approach for simulating drug docking and generating binding hypotheses.
Collapse
Affiliation(s)
- Helen M. Deeks
- Intangible Realities Laboratory, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
- Department of Computer Science, University of Bristol, Bristol, England, United Kingdom
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
| | - Rebecca K. Walters
- Intangible Realities Laboratory, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
- Department of Computer Science, University of Bristol, Bristol, England, United Kingdom
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
| | - Stephanie R. Hare
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
| | - Michael B. O’Connor
- Intangible Realities Laboratory, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
- Department of Computer Science, University of Bristol, Bristol, England, United Kingdom
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
| | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
- * E-mail: (AJM); (DRG)
| | - David R. Glowacki
- Intangible Realities Laboratory, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
- Department of Computer Science, University of Bristol, Bristol, England, United Kingdom
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, England, United Kingdom
- * E-mail: (AJM); (DRG)
| |
Collapse
|
5
|
IJzerman AP, Guo D. Drug-Target Association Kinetics in Drug Discovery. Trends Biochem Sci 2019; 44:861-871. [PMID: 31101454 DOI: 10.1016/j.tibs.2019.04.004] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2019] [Revised: 03/28/2019] [Accepted: 04/08/2019] [Indexed: 02/07/2023]
Abstract
The important role of ligand-receptor binding kinetics in drug design and discovery is increasingly recognized by the drug research community. Over the past decade, accumulating evidence has shown that optimizing the ligand's dissociation rate constant can lead to desirable duration of in vivo target occupancy and, hence, improved pharmacodynamic properties. However, the association rate constant as a pharmacological principle remains less investigated, whereas it can play an equally important role in the selection of drug candidates. This review provides a compilation and discussion of otherwise scarce and dispersed information on this topic, bringing to light the importance of drug-target association in kinetics-directed drug design and discovery.
Collapse
Affiliation(s)
- Adriaan P IJzerman
- Division of Drug Discovery and Safety, Leiden Academic Centre for Drug Research (LACDR), Leiden University, P.O. Box 9502, 2300, RA, Leiden, The Netherlands
| | - Dong Guo
- Jiangsu Key Laboratory of New Drug Research and Clinical Pharmacy, Xuzhou Medical University, 209 Tongshan Road, Xuzhou, 221004, Jiangsu, China.
| |
Collapse
|
6
|
Huang S, Zhang D, Mei H, Kevin M, Qu S, Pan X, Lu L. SMD-Based Interaction-Energy Fingerprints Can Predict Accurately the Dissociation Rate Constants of HIV-1 Protease Inhibitors. J Chem Inf Model 2018; 59:159-169. [DOI: 10.1021/acs.jcim.8b00567] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
| | | | | | | | | | - Xianchao Pan
- Department of Medicinal Chemistry, College of Pharmacy, Southwest Medical University, Luzhou, Sichuan 646000, China
| | | |
Collapse
|
7
|
Ganotra GK, Wade RC. Prediction of Drug-Target Binding Kinetics by Comparative Binding Energy Analysis. ACS Med Chem Lett 2018; 9:1134-1139. [PMID: 30429958 PMCID: PMC6231175 DOI: 10.1021/acsmedchemlett.8b00397] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2018] [Accepted: 10/04/2018] [Indexed: 12/02/2022] Open
Abstract
![]()
A growing
consensus is emerging that optimizing the drug–target
affinity alone under equilibrium conditions does not necessarily translate
into higher potency in vivo and that instead binding kinetic parameters
should be optimized to ensure better efficacy. Therefore, in silico
methods are needed to predict the kinetic parameters and the mechanistic
determinants of drug–protein binding. Here we demonstrate the
application of COMparative BINding Energy (COMBINE) analysis to derive
quantitative structure–kinetics relationships (QSKRs) for the
dissociation rate constants (koff) of
inhibitors of heat shock protein 90 (HSP90) and HIV-1 protease. We
derived protein-specific scoring functions by correlating koff rate constants with a subset of weighted
interaction energy components determined from the energy-minimized
structures of drug–protein complexes. As the QSKRs derived
for these sets of chemically diverse compounds have good predictive
ability and provide insights into important drug–protein interactions
for optimizing koff, COMBINE analysis
offers a promising approach for binding kinetics-guided lead optimization.
Collapse
Affiliation(s)
- Gaurav K. Ganotra
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Rebecca C. Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
8
|
Kumar PKR. Systematic screening of viral entry inhibitors using surface plasmon resonance. Rev Med Virol 2017; 27. [PMID: 29047180 DOI: 10.1002/rmv.1952] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Revised: 08/02/2017] [Accepted: 09/07/2017] [Indexed: 12/12/2022]
Abstract
Viral binding and entry into host cells for various viruses have been studied extensively, yielding a detailed understanding of the overall viral entry process. As cell entry is an essential and requisite process by which a virus initiates infection, it is an attractive target for therapeutic intervention. The advantages of targeting viral entry are an extracellular target site, relatively easy access for biological interventions, and lower toxicity. Several cell-based strategies and biophysical techniques have been used to screen compounds that block viral entry. These studies led to the discovery of inhibitors against HIV, HCV, influenza, Ebola, and RSV. In recent years, several compounds screened by fragment-based drug discovery have been approved as drugs or are in the final stages of clinical trials. Among fragment screening technologies, surface plasmon resonance has been widely used because it provides accurate information on binding kinetics, allows real-time monitoring of ligand-drug interactions, requires very small sample amounts to perform analyses, and requires no modifications to or labeling of ligands. This review focuses on surface plasmon resonance-based schemes for screening viral entry inhibitors.
Collapse
Affiliation(s)
- Penmetcha K R Kumar
- National Institute of Advanced Industrial Science and Technology, Tsukuba City, Ibaraki, Japan
| |
Collapse
|
9
|
Louis JM, Roche J. Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease. J Mol Biol 2016; 428:2780-92. [PMID: 27170547 PMCID: PMC4905781 DOI: 10.1016/j.jmb.2016.05.005] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2016] [Revised: 04/27/2016] [Accepted: 05/02/2016] [Indexed: 01/08/2023]
Abstract
Using high-pressure NMR spectroscopy and differential scanning calorimetry, we investigate the folding landscape of the mature HIV-1 protease homodimer. The cooperativity of unfolding was measured in the absence or presence of a symmetric active site inhibitor for the optimized wild type protease (PR), its inactive variant PRD25N, and an extremely multidrug-resistant mutant, PR20. The individual fit of the pressure denaturation profiles gives rise to first order, ∆GNMR, and second order, ∆VNMR (the derivative of ∆GNMR with pressure); apparent thermodynamic parameters for each amide proton considered. Heterogeneity in the apparent ∆VNMR values reflects departure from an ideal cooperative unfolding transition. The narrow to broad distribution of ∆VNMR spanning the extremes from inhibitor-free PR20D25N to PR-DMP323 complex, and distinctively for PRD25N-DMP323 complex, indicated large variations in folding cooperativity. Consistent with this data, the shape of thermal unfolding transitions varies from asymmetric for PR to nearly symmetric for PR20, as dimer-inhibitor ternary complexes. Lack of structural cooperativity was observed between regions located close to the active site, including the hinge and tip of the glycine-rich flaps, and the rest of the protein. These results strongly suggest that inhibitor binding drastically decreases the cooperativity of unfolding by trapping the closed flap conformation in a deep energy minimum. To evade this conformational trap, PR20 evolves exhibiting a smoother folding landscape with nearly an ideal two-state (cooperative) unfolding transition. This study highlights the malleability of retroviral protease folding pathways by illustrating how the selection of mutations under drug pressure remodels the free-energy landscape as a primary mechanism.
Collapse
Affiliation(s)
- John M Louis
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA.
| | - Julien Roche
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA; Roy J. Carver Department of Biochemistry, Biophysics and Molecular Biology, Iowa State University, Ames, IA 50011, USA.
| |
Collapse
|
10
|
Bauer RA. Covalent inhibitors in drug discovery: from accidental discoveries to avoided liabilities and designed therapies. Drug Discov Today 2015; 20:1061-73. [PMID: 26002380 DOI: 10.1016/j.drudis.2015.05.005] [Citation(s) in RCA: 346] [Impact Index Per Article: 38.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2014] [Revised: 03/25/2015] [Accepted: 05/11/2015] [Indexed: 02/07/2023]
Abstract
Drugs that covalently bond to their biological targets have a long history in drug discovery. A look at drug approvals in recent years suggests that covalent drugs will continue to make impacts on human health for years to come. Although fraught with concerns about toxicity, the high potencies and prolonged effects achievable with covalent drugs may result in less-frequent drug dosing and in wide therapeutic margins for patients. Covalent inhibition can also dissociate drug pharmacodynamics (PD) from pharmacokinetics (PK), which can result in desired drug efficacy for inhibitors that have short systemic exposure. Evidence suggests that there is a reduced risk for the development of resistance against covalent drugs, which is a major challenge in areas such as oncology and infectious disease.
Collapse
Affiliation(s)
- Renato A Bauer
- Lilly Research Laboratories, Eli Lilly and Co., Indianapolis, IN 46285, USA.
| |
Collapse
|
11
|
Mah R, Thomas JR, Shafer CM. Drug discovery considerations in the development of covalent inhibitors. Bioorg Med Chem Lett 2013; 24:33-9. [PMID: 24314671 DOI: 10.1016/j.bmcl.2013.10.003] [Citation(s) in RCA: 148] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2013] [Revised: 09/27/2013] [Accepted: 10/01/2013] [Indexed: 12/31/2022]
Abstract
In recent years, the number of drug candidates with a covalent mechanism of action progressing through clinical trials or being approved by the FDA has increased significantly. And as interest in covalent inhibitors has increased, the technical challenges for characterizing and optimizing these inhibitors have become evident. A number of new tools have been developed to aid this process, but these have not gained wide-spread use. This review will highlight a number of methods and tools useful for prosecuting covalent inhibitor drug discovery programs.
Collapse
Affiliation(s)
- Robert Mah
- Global Discovery Chemistry/Oncology & Exploratory Chemistry, Novartis Institutes for Biomedical Research, Klybeckstrasse 141, CH-4057 Basel, Switzerland
| | - Jason R Thomas
- Developmental and Molecular Pathways, Novartis Institutes for Biomedical Research, 250 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Cynthia M Shafer
- Global Discovery Chemistry/Oncology & Exploratory Chemistry, Novartis Institutes for Biomedical Research, 4560 Horton Street, Emeryville, CA 94608, USA.
| |
Collapse
|
12
|
Foulkes-Murzycki JE, Rosi C, Kurt Yilmaz N, Shafer RW, Schiffer CA. Cooperative effects of drug-resistance mutations in the flap region of HIV-1 protease. ACS Chem Biol 2013; 8:513-8. [PMID: 23252515 DOI: 10.1021/cb3006193] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Understanding the interdependence of multiple mutations in conferring drug resistance is crucial to the development of novel and robust inhibitors. As HIV-1 protease continues to adapt and evade inhibitors while still maintaining the ability to specifically recognize and efficiently cleave its substrates, the problem of drug resistance has become more complicated. Under the selective pressure of therapy, correlated mutations accumulate throughout the enzyme to compromise inhibitor binding, but characterizing their energetic interdependency is not straightforward. A particular drug resistant variant (L10I/G48V/I54V/V82A) displays extreme entropy-enthalpy compensation relative to wild-type enzyme but a similar variant (L10I/G48V/I54A/V82A) does not. Individual mutations of sites in the flaps (residues 48 and 54) of the enzyme reveal that the thermodynamic effects are not additive. Rather, the thermodynamic profile of the variants is interdependent on the cooperative effects exerted by a particular combination of mutations simultaneously present.
Collapse
Affiliation(s)
- Jennifer E. Foulkes-Murzycki
- Department of Biochemistry and
Molecular Pharmacology, University of Massachusetts Medical School, 364 Plantation Street, Worcester, Massachusetts
01605, United States
| | - Christina Rosi
- Department of Biochemistry and
Molecular Pharmacology, University of Massachusetts Medical School, 364 Plantation Street, Worcester, Massachusetts
01605, United States
| | - Nese Kurt Yilmaz
- Department of Biochemistry and
Molecular Pharmacology, University of Massachusetts Medical School, 364 Plantation Street, Worcester, Massachusetts
01605, United States
| | - Robert W. Shafer
- Division
of Infectious Diseases,
Department of Medicine, Stanford University, Stanford, California 94305, United States
| | - Celia A. Schiffer
- Department of Biochemistry and
Molecular Pharmacology, University of Massachusetts Medical School, 364 Plantation Street, Worcester, Massachusetts
01605, United States
| |
Collapse
|
13
|
How conformational changes can affect catalysis, inhibition and drug resistance of enzymes with induced-fit binding mechanism such as the HIV-1 protease. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:867-73. [PMID: 23376188 DOI: 10.1016/j.bbapap.2013.01.027] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2012] [Revised: 01/21/2013] [Accepted: 01/24/2013] [Indexed: 11/21/2022]
Abstract
A central question is how the conformational changes of proteins affect their function and the inhibition of this function by drug molecules. Many enzymes change from an open to a closed conformation upon binding of substrate or inhibitor molecules. These conformational changes have been suggested to follow an induced-fit mechanism in which the molecules first bind in the open conformation in those cases where binding in the closed conformation appears to be sterically obstructed such as for the HIV-1 protease. In this article, we present a general model for the catalysis and inhibition of enzymes with induced-fit binding mechanism. We derive general expressions that specify how the overall catalytic rate of the enzymes depends on the rates for binding, for the conformational changes, and for the chemical reaction. Based on these expressions, we analyze the effect of mutations that mainly shift the conformational equilibrium on catalysis and inhibition. If the overall catalytic rate is limited by product unbinding, we find that mutations that destabilize the closed conformation relative to the open conformation increase the catalytic rate in the presence of inhibitors by a factor exp(ΔΔGC/RT) where ΔΔGC is the mutation-induced shift of the free-energy difference between the conformations. This increase in the catalytic rate due to changes in the conformational equilibrium is independent of the inhibitor molecule and, thus, may help to understand how non-active-site mutations can contribute to the multi-drug-resistance that has been observed for the HIV-1 protease. A comparison to experimental data for the non-active-site mutation L90M of the HIV-1 protease indicates that the mutation slightly destabilizes the closed conformation of the enzyme. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.
Collapse
|
14
|
ElSawy KM, Twarock R, Lane DP, Verma CS, Caves LSD. Characterization of the Ligand Receptor Encounter Complex and Its Potential for in Silico Kinetics-Based Drug Development. J Chem Theory Comput 2011; 8:314-21. [PMID: 26592892 DOI: 10.1021/ct200560w] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The study of drug-receptor interactions has largely been framed in terms of the equilibrium thermodynamic binding affinity, an in vitro measure of the stability of the drug-receptor complex that is commonly used as a proxy measure of in vivo biological activity. In response to the growing realization of the importance of binding kinetics to in vivo drug activity we present a computational methodology for the kinetic characterization of drug-receptor interactions in terms of the encounter complex. Using trajectory data from multiple Brownian dynamics simulations of ligand diffusion, we derive the spatial density of the ligand around the receptor and show how it can be quantitatively partitioned into different basins of attraction. Numerical integration of the ligand densities within the basins can be used to estimate the residence time of the ligand within these diffusive binding sites. Simulations of two structurally similar inhibitors of Hsp90 exhibit diffusive binding sites with similar spatial structure but with different ligand residence times. In contrast, a pair of structurally dissimilar inhibitors of MDM2, a peptide and a small molecule, exhibit spatially distinct basins of attraction around the receptor, which in turn reveal differences in ligand orientational order. Thus, our kinetic approach provides microscopic details of drug-receptor dynamics that provide novel insight into the observed differences in the thermodynamic binding affinities for the two inhibitors, such as the differences in the entropic contributions to binding. The characterization of the encounter complex, in terms of the structure, topology, and dynamics of diffusive binding sites, offers a new perspective on ligand-receptor interactions and the potential for greater insight into drug action. The method, which requires no prior knowledge of the bound state, is a first step toward the incorporation of ligand kinetics into in silico drug development protocols.
Collapse
Affiliation(s)
| | | | - David P Lane
- P53 Laboratory (p53Lab, A* STAR), 8A Biomedical Grove 06-06, Immunos, Singapore 138648
| | - Chandra S Verma
- Bioinformatics Institute (A*STAR), 30 Biopolis Str., 07-01 Matrix , Singapore 138671
| | | |
Collapse
|
15
|
Kang M, Roberts C, Cheng Y, Chang CEA. Gating and Intermolecular Interactions in Ligand-Protein Association: Coarse-Grained Modeling of HIV-1 Protease. J Chem Theory Comput 2011; 7:3438-46. [PMID: 26598172 DOI: 10.1021/ct2004885] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Most biological processes are initiated or mediated by the association of ligands and proteins. This work studies multistep, ligand-protein association processes by Brownian dynamics simulations with coarse-grained models for HIV-1 protease (HIVp) and its neutral ligands. We report the average association times when the ligand concentration is 100 μM. The influence of crowding on the simulated binding time was also studied. HIVp has flexible loops that serve as a gate during the ligand binding processes. It is believed that the flaps are partially closed most of the time in its free state. To accelerate our simulations, we fixed a part of the HIVp and reparameterized our coarse-grained model, using atomistic molecular dynamics simulations, to reproduce the "gating" motions of HIVp. HIVp-ligand interactions changed the gating behavior of HIVp and helped ligands diffuse on HIVp surface to accelerate binding. The structural adjustment of the ligand toward its final stable state was the limiting step in the binding processes, which is highly system dependent. The intermolecular attraction between the ligands and crowder proteins contributes the most to the crowding effects. The results highlight broader implications in recognition pathways under more complex environment that considers molecular dynamics and conformational changes. This work brings insights into ligand-protein associations and is helpful in the design of targeted ligands.
Collapse
Affiliation(s)
- Myungshim Kang
- Department of Chemistry, University of California , Riverside, California, United States
| | - Christopher Roberts
- Department of Chemistry, University of California , Riverside, California, United States
| | - Yuhui Cheng
- Pacific Northwest National Laboratory , Richland, Washington, United States
| | - Chia-En A Chang
- Department of Chemistry, University of California , Riverside, California, United States
| |
Collapse
|
16
|
Ali A, Bandaranayake RM, Cai Y, King NM, Kolli M, Mittal S, Murzycki JF, Nalam MN, Nalivaika EA, Özen A, Prabu-Jeyabalan MM, Thayer K, Schiffer CA. Molecular Basis for Drug Resistance in HIV-1 Protease. Viruses 2010; 2:2509-2535. [PMID: 21994628 PMCID: PMC3185577 DOI: 10.3390/v2112509] [Citation(s) in RCA: 98] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2010] [Revised: 10/22/2010] [Accepted: 10/28/2010] [Indexed: 02/01/2023] Open
Abstract
HIV-1 protease is one of the major antiviral targets in the treatment of patients infected with HIV-1. The nine FDA approved HIV-1 protease inhibitors were developed with extensive use of structure-based drug design, thus the atomic details of how the inhibitors bind are well characterized. From this structural understanding the molecular basis for drug resistance in HIV-1 protease can be elucidated. Selected mutations in response to therapy and diversity between clades in HIV-1 protease have altered the shape of the active site, potentially altered the dynamics and even altered the sequence of the cleavage sites in the Gag polyprotein. All of these interdependent changes act in synergy to confer drug resistance while simultaneously maintaining the fitness of the virus. New strategies, such as incorporation of the substrate envelope constraint to design robust inhibitors that incorporate details of HIV-1 protease’s function and decrease the probability of drug resistance, are necessary to continue to effectively target this key protein in HIV-1 life cycle.
Collapse
Affiliation(s)
- Akbar Ali
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Rajintha M. Bandaranayake
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Yufeng Cai
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Nancy M. King
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Madhavi Kolli
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Seema Mittal
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Jennifer F. Murzycki
- Department of Pediatrics, University of Rochester, Rochester, NY 14627, USA; E-Mail:
| | - Madhavi N.L. Nalam
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Ellen A. Nalivaika
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Ayşegül Özen
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Moses M. Prabu-Jeyabalan
- Division of Basic Sciences, The Commonwealth Medical College, 150 N. Washington Avenue, Scranton, PA 18503, USA; E-Mail:
| | - Kelly Thayer
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
| | - Celia A. Schiffer
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA; E-Mails: (A.A.); (R.M.B.); (Y.C.); (N.M.K.); (M.K.); (S.M.), (M.N.L.N.); (E.A.N.); (A.Ö.); (K.T.)
- Author to whom correspondence should be addressed; E-Mail: ; Tel.: +1-508-856-8008; Fax: +1-508-856-6464
| |
Collapse
|
17
|
Pietrucci F, Marinelli F, Carloni P, Laio A. Substrate binding mechanism of HIV-1 protease from explicit-solvent atomistic simulations. J Am Chem Soc 2009; 131:11811-8. [PMID: 19645490 DOI: 10.1021/ja903045y] [Citation(s) in RCA: 121] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The binding mechanism of a peptide substrate (Thr-Ile-Met-Met-Gln-Arg, cleavage site p2-NC of the viral polyprotein) to wild-type HIV-1 protease has been investigated by 1.6 micros biased all-atom molecular dynamics simulations in explicit water. The configuration space has been explored biasing seven reaction coordinates by the bias-exchange metadynamics technique. The structure of the Michaelis complex is obtained starting from the substrate outside the enzyme within a backbone rmsd of 0.9 A. The calculated free energy of binding is -6 kcal/mol, and the kinetic constants for association and dissociation are 1.3 x 10(6) M(-1) s(-1) and 57 s(-1), respectively, consistent with experiments. In the main binding pathway, the flaps of the protease do not open sizably. The substrate slides inside the enzyme cavity from the tight lateral channel. This may contrast with the natural polyprotein substrate which is expected to bind by opening the flaps. Thus, mutations might influence differently the binding kinetics of peptidomimetic ligands and of the natural substrate.
Collapse
Affiliation(s)
- Fabio Pietrucci
- International School for Advanced Studies (SISSA-ISAS), via Beirut 2-4, I-34014 Trieste, Italy
| | | | | | | |
Collapse
|
18
|
Liu X, Xiu Z, Hao C. Drug-resistant molecular mechanism of CRF01_AE HIV-1 protease due to V82F mutation. J Comput Aided Mol Des 2009; 23:261-72. [PMID: 19219633 DOI: 10.1007/s10822-008-9256-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2008] [Accepted: 12/15/2008] [Indexed: 11/28/2022]
Abstract
Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the major targets of anti-AIDS drug discovery. The circulating recombinant form 01 A/E (CRF01_AE, abbreviated AE) subtype is one of the most common HIV-1 subtypes, which is infecting more humans and is expanding rapidly throughout the world. It is, therefore, necessary to develop inhibitors against subtype AE HIV-1 PR. In this work, we have performed computer simulation of subtype AE HIV-1 PR with the drugs lopinavir (LPV) and nelfinavir (NFV), and examined the mechanism of resistance of the V82F mutation of this protease against LPV both structurally and energetically. The V82F mutation at the active site results in a conformational change of 79's loop region and displacement of LPV from its proper binding site, and these changes lead to rotation of the side-chains of residues D25 and I50'. Consequently, the conformation of the binding cavity is deformed asymmetrically and some interactions between PR and LPV are destroyed. Additionally, by comparing the interactive mechanisms of LPV and NFV with HIV-1 PR we discovered that the presence of a dodecahydroisoquinoline ring at the P1' subsite, a [2-(2,6-dimethylphenoxy)acetyl]amino group at the P2' subsite, and an N2 atom at the P2 subsite could improve the binding affinity of the drug with AE HIV-1 PR. These findings are helpful for promising drug design.
Collapse
Affiliation(s)
- Xiaoqing Liu
- Department of Bioscience and Biotechnology, School of Environmental and Biological Science and Technology, Dalian University of Technology, Dalian 116024, People's Republic of China.
| | | | | |
Collapse
|
19
|
Rajagopalan R, Misialek S, Stevens SK, Myszka DG, Brandhuber BJ, Ballard JA, Andrews SW, Seiwert SD, Kossen K. Inhibition and Binding Kinetics of the Hepatitis C Virus NS3 Protease Inhibitor ITMN-191 Reveals Tight Binding and Slow Dissociative Behavior. Biochemistry 2009; 48:2559-68. [DOI: 10.1021/bi900038p] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Ravi Rajagopalan
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - Shawn Misialek
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - Sarah K. Stevens
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - David G. Myszka
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - Barbara J. Brandhuber
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - Joshua A. Ballard
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - Steven W. Andrews
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - Scott D. Seiwert
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| | - Karl Kossen
- InterMune Inc., 3280 Bayshore Boulevard, Brisbane, California 94005, Biosensor Tools, 1588 East Connecticut Drive, Salt Lake City, Utah 84132, and Array Biopharma, 3200 Walnut Street, Boulder, Colorado 80301
| |
Collapse
|
20
|
Anomalous adsorptive properties of HIV protease: Indication of two-dimensional crystallization? Colloids Surf B Biointerfaces 2008; 64:145-9. [DOI: 10.1016/j.colsurfb.2008.01.011] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2007] [Accepted: 01/15/2008] [Indexed: 11/23/2022]
|
21
|
Abstract
PURPOSE OF REVIEW Individuals with HIV are at greater risk of human papillomavirus-related cancers. This report will assess the potential and limitations of vaccines against human papillomavirus in HIV-positive individuals. RECENT FINDINGS A worldwide meta-analysis of published data established the under-representation of HPV16, and increased prevalence of multiple-type human papillomavirus infections in HIV-positive women. Associations between HIV-related immunodepression, the progression of human papillomavirus infection to cervical lesions, and an increased risk of cervical cancer in women with HIV have also been shown. An increased incidence of human papillomavirus infection in anal and vulvar/vaginal neoplasia has been reported in individuals with HIV. A prophylactic vaccine against HPV6, 11, 16 and 18 has been licensed, and one against HPV16 and 18 is under evaluation. Both have shown efficacy against persistent infection, as well as related human papillomavirus cervical lesions for up to 5 years. Preliminary results have also been reported on therapeutic vaccines, notably for the treatment of cervical intraepithelial neoplasia grades 2 and 3. SUMMARY The safety and efficacy of human papillomavirus vaccines in individuals with HIV need to be assessed to prevent cervical cancer in current and future generations. Screening for human papillomavirus-related cancer among HIV patients should be undertaken.
Collapse
Affiliation(s)
- Hugo De Vuyst
- International Agency for Research on Cancer, Lyon, France.
| | | |
Collapse
|
22
|
Saen-oon S, Aruksakunwong O, Wittayanarakul K, Sompornpisut P, Hannongbua S. Insight into analysis of interactions of saquinavir with HIV-1 protease in comparison between the wild-type and G48V and G48V/L90M mutants based on QM and QM/MM calculations. J Mol Graph Model 2007; 26:720-7. [PMID: 17543558 DOI: 10.1016/j.jmgm.2007.04.009] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2006] [Revised: 04/19/2007] [Accepted: 04/28/2007] [Indexed: 11/26/2022]
Abstract
Saquinavir (SQV) was the first HIV-1 PR inhibitor licensed for clinical use and widely used for acquired immunodeficiency syndrome (AIDS) therapy. Its effectiveness, however, has been hindered by the emergence of resistant mutations. The two most important HIV-1 PR mutants are G48V and G48V/L90M. Inhibition studies of SQV on these mutants demonstrated 13.5- and 419-fold reductions of susceptibility, respectively. In this study, an analysis of energetic binding affinity between saquinavir and the HIV-1 PR wild-type and these two mutants has been performed in detail based on density functional theory and the hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. We have found that the interaction of SQV with flap residue 48 of the mutants is significantly perturbed, as shown by the reduced stability of binding between SQV and residue 48 for the G48V and G48V/L90M mutants over the wild-type. This was associated with conformational changes of the inhibitor and the enzyme, leading to the loss of hydrogen bonding between the binding subsite P2 and the backbone carbonyl of residue 48. Moreover, the G48V/L90M mutations cause the repositioning of the residues close to residues 48 and 90, at important locations as a part of the flap and catalytic regions, respectively. The repositioning of these residues consequently perturbed the binding affinity of SQV in the pocket as indicated by the decreasing interaction energies. In addition to the loss of inhibitor/enzyme binding, it is interesting to observe that the mutation leads significantly to an increase of the stability of the enzyme.
Collapse
Affiliation(s)
- Suwipa Saen-oon
- Department of Chemistry, Faculty of Science, Chulalongkorn University, Prathumwan, Bangkok 10330, Thailand
| | | | | | | | | |
Collapse
|
23
|
Chang CEA, Trylska J, Tozzini V, McCammon JA. Binding pathways of ligands to HIV-1 protease: coarse-grained and atomistic simulations. Chem Biol Drug Des 2007; 69:5-13. [PMID: 17313452 DOI: 10.1111/j.1747-0285.2007.00464.x] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Multiscale simulations (coarse-grained Brownian dynamics simulations and all-atom molecular dynamics simulations in implicit solvent) were applied to reveal the binding processes of ligands as they enter the binding site of the HIV-1 protease. The initial structures used for the molecular dynamics simulations were generated based on the Brownian dynamics trajectories, and this is the first molecular dynamics simulation of modeling the association of a ligand with the protease. We found that a protease substrate successfully binds to the protein when the flaps are fully open. Surprisingly, a smaller cyclic urea inhibitor (XK263) can reach the binding site when the flaps are not fully open. However, if the flaps are nearly closed, the inhibitor must rearrange or binding can fail because the inhibitor cannot attain proper conformations to enter the binding site. Both the peptide substrate and XK263 can also affect the protein's internal motion, which may help the flaps to open. Simulations allow us to efficiently study the ligand binding processes and may help those who study drug discovery to find optimal association pathways and to design those ligands with the best binding kinetics.
Collapse
Affiliation(s)
- Chia-En A Chang
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California 92093, USA.
| | | | | | | |
Collapse
|
24
|
Kerman K, Mahmoud KA, Kraatz HB. An electrochemical approach for the detection of HIV-1 protease. Chem Commun (Camb) 2007:3829-31. [DOI: 10.1039/b707140j] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
25
|
Copeland RA, Pompliano DL, Meek TD. Drug-target residence time and its implications for lead optimization. Nat Rev Drug Discov 2006; 5:730-9. [PMID: 16888652 DOI: 10.1038/nrd2082] [Citation(s) in RCA: 1027] [Impact Index Per Article: 57.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Much of drug discovery today is predicated on the concept of selective targeting of particular bioactive macromolecules by low-molecular-mass drugs. The binding of drugs to their macromolecular targets is therefore seen as paramount for pharmacological activity. In vitro assessment of drug-target interactions is classically quantified in terms of binding parameters such as IC(50) or K(d). This article presents an alternative perspective on drug optimization in terms of drug-target binary complex residence time, as quantified by the dissociative half-life of the drug-target binary complex. We describe the potential advantages of long residence time in terms of duration of pharmacological effect and target selectivity.
Collapse
Affiliation(s)
- Robert A Copeland
- Department of Enzymology and Mechanistic Pharmacology, at GlaxoSmithKline, UP1345, 1250 South Collegeville Road, Collegeville, Pennsylvania 19426, USA.
| | | | | |
Collapse
|
26
|
Chang CE, Shen T, Trylska J, Tozzini V, McCammon JA. Gated binding of ligands to HIV-1 protease: Brownian dynamics simulations in a coarse-grained model. Biophys J 2006; 90:3880-5. [PMID: 16533835 PMCID: PMC1459512 DOI: 10.1529/biophysj.105.074575] [Citation(s) in RCA: 73] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The internal motions of proteins may serve as a "gate" in some systems, which controls ligand-protein association. This study applies Brownian dynamics simulations in a coarse-grained model to study the gated association rate constants of HIV-1 proteases and drugs. The computed gated association rate constants of three protease mutants, G48V/V82A/I84V/L90M, G48V, and L90M with three drugs, amprenavir, indinavir, and saquinavir, yield good agreements with experiments. The work shows that the flap dynamics leads to "slow gating". The simulations suggest that the flap flexibility and the opening frequency of the wild-type, the G48V and L90M mutants are similar, but the flaps of the variant G48V/V82A/I84V/L90M open less frequently, resulting in a lower gated rate constant. The developed methodology is fast and provides an efficient way to predict the gated association rate constants for various protease mutants and ligands.
Collapse
Affiliation(s)
- Chia-En Chang
- Department of Chemistry and Biochemistry, Center for Theoretical Biological Physics, Howard Hughes Medical Institute, University of California at San Diego, La Jolla, CA, USA.
| | | | | | | | | |
Collapse
|
27
|
Aruksakunwong O, Wittayanarakul K, Sompornpisut P, Sanghiran V, Parasuk V, Hannongbua S. Structural and dynamical properties of different protonated states of mutant HIV-1 protease complexed with the saquinavir inhibitor studied by molecular dynamics simulations. J Mol Graph Model 2006; 25:324-32. [PMID: 16504560 DOI: 10.1016/j.jmgm.2006.01.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2005] [Revised: 01/17/2006] [Accepted: 01/18/2006] [Indexed: 11/22/2022]
Abstract
To understand the basis of drug resistance, particularly of the HIV-1 PR, three molecular dynamics (MD) simulations of HIV-1 PR mutant species, G48V, complexed with saquinavir (SQV) in explicit aqueous solution with three protonation states, diprotonation on Asp25 and Asp25' (Di-pro) and monoprotonation on each Asp residue (Mono-25 and Mono-25'). For all three states, H-bonds between saquinavir and HIV-1 PR were formed only in the two regions, flap and active site. It was found that conformation of P2 subsite of SQV in the Mono-25 state differs substantially from the other two states. The rotation about 177 degrees from the optimal structure of the wild type was observed, the hydrogen bond between P2 and the flap residue (Val48) was broken and indirect hydrogen bonds with the three residues (Asp29, Gly27, and Asp30) were found instead. In terms of complexation energies, interaction energy of -37.3 kcal/mol for the Mono-25 state is significantly lower than those of -30.7 and -10.7kcal/mol for the Mono-25' and Di-pro states, respectively. It was found also that protonation at the Asp25 leads to a better arrangement in the catalytic dyad, i.e., the Asp25-Asp25' interaction energy of -8.8 kcal/mol of the Mono-25 is significantly lower than that of -2.6kcal/mol for the Mono-25' state. The above data suggest us to conclude that interaction in the catalytic area should be used as criteria to enhance capability in drug designing and drug screening instead of using the total inhibitor/enzyme interaction.
Collapse
Affiliation(s)
- Ornjira Aruksakunwong
- Department of Chemistry, Faculty of Science, Chulalongkorn University, Prathumwan, Bangkok 10330, Thailand
| | | | | | | | | | | |
Collapse
|
28
|
Cimitan S, Lindgren MT, Bertucci C, Danielson UH. Early Absorption and Distribution Analysis of Antitumor and Anti-AIDS Drugs: Lipid Membrane and Plasma Protein Interactions. J Med Chem 2005; 48:3536-46. [PMID: 15887962 DOI: 10.1021/jm049343o] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The interactions of a set of compounds of potential importance for anticancer and AIDS chemotherapy with lipid membranes and plasma proteins were studied with a surface plasmon resonance (SPR) based optical biosensor, giving valuable information on the absorption and distribution of the compounds. The technique allowed both effective screening of compounds and more detailed kinetic and mechanistic analysis of specific interactions. The interaction with two different types of lipid membranes could be reliably measured at a drug concentration as low as 20 microM, allowing analysis of poorly soluble compounds. Distribution was evaluated by investigation of the interactions with two human plasma proteins, human serum albumin (HSA) and alpha(1)-acid glycoprotein (AGP). Two apparent binding sites were clearly defined for HSA: one with rapid and one with slow association and dissociation rates. The sites appear to differ in accessibility and recognition characteristics rather than in their capacities to form strong complexes with drugs.
Collapse
Affiliation(s)
- Samanta Cimitan
- Department of Biochemistry, Uppsala University, Box 576, SE-751 23 Uppsala, Sweden
| | | | | | | |
Collapse
|
29
|
Abstract
In the year 2003 there was a 17% increase in the number of publications citing work performed using optical biosensor technology compared with the previous year. We collated the 962 total papers for 2003, identified the geographical regions where the work was performed, highlighted the instrument types on which it was carried out, and segregated the papers by biological system. In this overview, we spotlight 13 papers that should be on everyone's 'must read' list for 2003 and provide examples of how to identify and interpret high-quality biosensor data. Although we still find that the literature is replete with poorly performed experiments, over-interpreted results and a general lack of understanding of data analysis, we are optimistic that these shortcomings will be addressed as biosensor technology continues to mature.
Collapse
Affiliation(s)
- Rebecca L Rich
- Center for Biomolecular Interaction Analysis, University of Utah, Salt Lake City, UT 84132, USA
| | | |
Collapse
|
30
|
Huber W. A new strategy for improved secondary screening and lead optimization using high-resolution SPR characterization of compound-target interactions. J Mol Recognit 2005; 18:273-81. [PMID: 15997470 DOI: 10.1002/jmr.744] [Citation(s) in RCA: 77] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Biophysical label-free assays such as those based on SPR are essential tools in generating high-quality data on affinity, kinetic, mechanistic and thermodynamic aspects of interactions between target proteins and potential drug candidates. Here we show examples of the integration of SPR with bioinformatic approaches and mutation studies in the early drug discovery process. We call this combination 'structure-based biophysical analysis'. Binding sites are identified on target proteins using information that is either extracted from three-dimensional structural analysis (X-ray crystallography or NMR), or derived from a pharmacore model based on known binders. The binding site information is used for in silico screening of a large substance library (e.g. available chemical directory), providing virtual hits. The three-dimensional structure is also used for the design of mutants where the binding site has been impaired. The wild-type target and the impaired mutant are then immobilized on different spots of the sensor chip and the interactions of compounds with the wild-type and mutant are compared in order to identify selective binders for the binding site of the target protein. This method can be used as a cost-effective alternative to high-throughput screening methods in cases when detailed binding site information is available. Here, we present three examples of how this technique can be applied to provide invaluable data during different phases of the drug discovery process.
Collapse
Affiliation(s)
- Walter Huber
- F. Hoffman-La Roche Ltd, Pharmaceutical Research, Discovery Chemistry, Basel, Switzerland.
| |
Collapse
|
31
|
Shuman CF, Vrang L, Danielson UH. Improved Structure−Activity Relationship Analysis of HIV-1 Protease Inhibitors Using Interaction Kinetic Data. J Med Chem 2004; 47:5953-61. [PMID: 15537350 DOI: 10.1021/jm0499110] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Despite the availability of large amounts of data for HIV-protease inhibitors and their effectiveness with wild type and resistant enzyme, there is limited knowledge about how this and other information can be systematically applied to the development of new antiviral compounds. To identify in vitro parameters that correlate with the efficacy of HIV inhibitors in cell culture, the relationships between inhibition, interaction kinetic, and cell culture parameters for HIV-1 protease inhibitors were analyzed. Correlation, cluster, and principal component analysis of data for 37 cyclic and linear compounds revealed that the affinities (K(D)) determined from SPR-biosensor binding studies correlated better to cell culture efficacy (ED(50)) than that of the inhibition constants (K(i)), indicating that the conventional use of K(i) values for structure-activity relationship analysis of HIV-1 inhibitors should be seriously reconsidered. The association and dissociation kinetic rate constants (k(on) and k(off)) alone showed weak correlations with ED(50) values. However, ED(50) values were most related to the free enzyme concentration in the viral particle ([E]), calculated from the rate constants and the total enzyme concentration in a viral particle. A structure-activity relationship analysis of the current data set was found to be valid for all classes of compounds analyzed. In summary, use of affinity, based on interaction kinetic rate constants, rather than inhibition constants, and theoretical consideration of the physiological conditions in the virus particle provide improved structure-activity relationship analysis of HIV-1 protease inhibitors.
Collapse
Affiliation(s)
- Cynthia F Shuman
- Department of Biochemistry, Uppsala University, BMC, Box 576, SE-751 23 Uppsala, Sweden
| | | | | |
Collapse
|