Computational design and experimental study of tighter binding peptides to an inactivated mutant of HIV‐1 protease

Drug resistance in HIV‐1 protease, a barrier to effective treatment, is generally caused by mutations in the enzyme that disrupt inhibitor binding but still allow for substrate processing. Structural studies with mutant, inactive enzyme, have provided detailed information regarding how the substrates bind to the protease yet avoid resistance mutations; insights obtained inform the development of next generation therapeutics. Although structures have been obtained of complexes between substrate peptide and inactivated (D25N) protease, thermodynamic studies of peptide binding have been challenging due to low affinity. Peptides that bind tighter to the inactivated protease than the natural substrates would be valuable for thermodynamic studies as well as to explore whether the structural envelope observed for substrate peptides is a function of weak binding. Here, two computational methods—namely, charge optimization and protein design—were applied to identify peptide sequences predicted to have higher binding affinity to the inactivated protease, starting from an RT–RH derived substrate peptide. Of the candidate designed peptides, three were tested for binding with isothermal titration calorimetry, with one, containing a single threonine to valine substitution, measured to have more than a 10‐fold improvement over the tightest binding natural substrate. Crystal structures were also obtained for the same three designed peptide complexes; they show good agreement with computational prediction. Thermodynamic studies show that binding is entropically driven, more so for designed affinity enhanced variants than for the starting substrate. Structural studies show strong similarities between natural and tighter‐binding designed peptide complexes, which may have implications in understanding the molecular mechanisms of drug resistance in HIV‐1 protease. Proteins 2008. © 2007 Wiley‐Liss, Inc.

[1]  K. N. Trueblood,et al.  On the rigid-body motion of molecules in crystals , 1968 .

[2]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[3]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[4]  Conrad C. Huang,et al.  The MIDAS display system , 1988 .

[5]  B. Honig,et al.  Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis , 1988, Proteins.

[6]  M Karplus,et al.  Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison , 1988, Proteins.

[7]  K. Sharp,et al.  Calculating the electrostatic potential of molecules in solution: Method and error assessment , 1988 .

[8]  S D Kemp,et al.  Multiple mutations in HIV-1 reverse transcriptase confer high-level resistance to zidovudine (AZT). , 1989, Science.

[9]  Charles E. Bugg,et al.  Crystallographic and Modeling Methods in Molecular Design , 1990, Springer New York.

[10]  M. Kjeldgaard,et al.  O: A Macromolecule Modeling Environment , 1990 .

[11]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[12]  Barry Honig,et al.  Calculating total electrostatic energies with the nonlinear Poisson-Boltzmann equation , 1990 .

[13]  J Kuriyan,et al.  Rigid protein motion as a model for crystallographic temperature factors. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[14]  K. Sharp,et al.  Protein folding and association: Insights from the interfacial and thermodynamic properties of hydrocarbons , 1991, Proteins.

[15]  T. Meek,et al.  Human immunodeficiency virus-1 protease. 2. Use of pH rate studies and solvent kinetic isotope effects to elucidate details of chemical mechanism. , 1991, Biochemistry.

[16]  R. L. Baldwin,et al.  The mechanism of alpha-helix formation by peptides. , 1992, Annual review of biophysics and biomolecular structure.

[17]  Johan Desmet,et al.  The dead-end elimination theorem and its use in protein side-chain positioning , 1992, Nature.

[18]  C. Debouck,et al.  The HIV-1 protease as a therapeutic target for AIDS. , 1992, AIDS research and human retroviruses.

[19]  J. Thornton,et al.  PROCHECK: a program to check the stereochemical quality of protein structures , 1993 .

[20]  D. P. Brunner,et al.  Large scale purification and refolding of HIV-1 protease fromEscherichia coli inclusion bodies , 1993, Journal of protein chemistry.

[21]  Collaborative Computational,et al.  The CCP4 suite: programs for protein crystallography. , 1994, Acta crystallographica. Section D, Biological crystallography.

[22]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[23]  Roland L. Dunbrack,et al.  Conformational analysis of the backbone-dependent rotamer preferences of protein sidechains , 1994, Nature Structural Biology.

[24]  W N Hunter,et al.  Structure of trypanothione reductase from Crithidia fasciculata at 2.6 A resolution; enzyme-NADP interactions at 2.8 A resolution. , 1994, Acta crystallographica. Section D, Biological crystallography.

[25]  Charles J. Eyermann,et al.  NMR and X-ray Evidence That the HIV Protease Catalytic Aspartyl Groups Are Protonated in the Complex Formed by the Protease and a Non-Peptide Cyclic Urea-Based Inhibitor , 1994 .

[26]  J. Navaza,et al.  AMoRe: an automated package for molecular replacement , 1994 .

[27]  J. Coffin,et al.  HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy , 1995, Science.

[28]  J. Condra,et al.  In vivo emergence of HIV-1 variants resistant to multiple protease inhibitors , 1995, Nature.

[29]  M. Karplus,et al.  Simulation of activation free energies in molecular systems , 1996 .

[30]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[31]  D. Ho,et al.  Ordered accumulation of mutations in HIV protease confers resistance to ritonavir , 1996, Nature Medicine.

[32]  S. L. Mayo,et al.  Protein design automation , 1996, Protein science : a publication of the Protein Society.

[33]  Andrew Carr,et al.  HIV protease inhibitors , 1996, AIDS.

[34]  I. Luque,et al.  Structure-based thermodynamic analysis of HIV-1 protease inhibitors. , 1997, Biochemistry.

[35]  S L Mayo,et al.  De novo protein design: towards fully automated sequence selection. , 1997, Journal of molecular biology.

[36]  Bruce Tidor,et al.  Optimization of electrostatic binding free energy , 1997 .

[37]  G. Murshudov,et al.  Refinement of macromolecular structures by the maximum-likelihood method. , 1997, Acta crystallographica. Section D, Biological crystallography.

[38]  Crystallographic analysis of human immunodeficiency virus 1 protease with an analog of the conserved CA-p2 substrate -- interactions with frequently occurring glutamic acid residue at P2' position of substrates. , 1997 .

[39]  E A Merritt,et al.  Raster3D: photorealistic molecular graphics. , 1997, Methods in enzymology.

[40]  D. Lamarre,et al.  Impaired fitness of human immunodeficiency virus type 1 variants with high-level resistance to protease inhibitors , 1997, Journal of virology.

[41]  P. S. Kim,et al.  High-resolution protein design with backbone freedom. , 1998, Science.

[42]  A R Leach,et al.  Exploring the conformational space of protein side chains using dead‐end elimination and the A* algorithm , 1998, Proteins.

[43]  I. Luque,et al.  Molecular basis of resistance to HIV-1 protease inhibition: a plausible hypothesis. , 1998, Biochemistry.

[44]  Bruce Tidor,et al.  Optimizing electrostatic affinity in ligand-receptor binding: Theory, computation, and ligand properties , 1998 .

[45]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[46]  J. Martinez-Picado,et al.  Replicative Fitness of Protease Inhibitor-Resistant Mutants of Human Immunodeficiency Virus Type 1 , 1999, Journal of Virology.

[47]  R. Vanderbei LOQO user's manual — version 3.10 , 1999 .

[48]  Joanna Trylska,et al.  Thermodynamic linkage between the binding of protons and inhibitors to HIV‐1 protease , 2008, Protein science : a publication of the Protein Society.

[49]  B Tidor,et al.  Charge optimization leads to favorable electrostatic binding free energy. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[50]  R. Vanderbei LOQO:an interior point code for quadratic programming , 1999 .

[51]  A. Velázquez‐Campoy,et al.  HIV-1 protease inhibitors: enthalpic versus entropic optimization of the binding affinity. , 2000, Biochemistry.

[52]  Stephen L. Mayo,et al.  Conformational splitting: A more powerful criterion for dead-end elimination , 2000, J. Comput. Chem..

[53]  Bruce Tidor,et al.  Electrostatic specificity in molecular ligand design , 2000 .

[54]  C. Schiffer,et al.  How does a symmetric dimer recognize an asymmetric substrate? A substrate complex of HIV-1 protease. , 2000, Journal of molecular biology.

[55]  J. Louis,et al.  Structural implications of drug‐resistant mutants of HIV‐1 protease: High‐resolution crystal structures of the mutant protease/substrate analogue complexes , 2001, Proteins.

[56]  Bruce Tidor,et al.  Optimization of binding electrostatics: Charge complementarity in the barnase‐barstar protein complex , 2001, Protein science : a publication of the Protein Society.

[57]  Kam Y. J. Zhang,et al.  Conversion of monomeric protein L to an obligate dimer by computational protein design , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[58]  Bruce Tidor,et al.  Electrostatic Complementarity at Ligand Binding Sites: Application to Chorismate Mutase , 2001 .

[59]  Bruce Tidor,et al.  Barstar is electrostatically optimized for tight binding to barnase , 2001, Nature Structural Biology.

[60]  Enrico O. Purisima,et al.  Optimizing Ligand Charges for Maximum Binding Affinity. A Solvated Interaction Energy Approach , 2001 .

[61]  Jae Young Lee,et al.  Crystal structure and functional analysis of the SurE protein identify a novel phosphatase family , 2001, Nature Structural Biology.

[62]  C. Schiffer,et al.  Substrate shape determines specificity of recognition for HIV-1 protease: analysis of crystal structures of six substrate complexes. , 2002, Structure.

[63]  Jeffery G Saven,et al.  Combinatorial protein design. , 2002, Current opinion in structural biology.

[64]  Julia M. Shifman,et al.  Modulating calmodulin binding specificity through computational protein design. , 2002, Journal of molecular biology.

[65]  Richard J Morris,et al.  ARP/wARP's model-building algorithms. I. The main chain. , 2002, Acta crystallographica. Section D, Biological crystallography.

[66]  L. Looger,et al.  Computational design of receptor and sensor proteins with novel functions , 2003, Nature.

[67]  D. Benjamin Gordon,et al.  Exact rotamer optimization for protein design , 2003, J. Comput. Chem..

[68]  C. Schiffer,et al.  Viability of a Drug-Resistant Human Immunodeficiency Virus Type 1 Protease Variant: Structural Insights for Better Antiviral Therapy , 2003, Journal of Virology.

[69]  D. Baker,et al.  A large scale test of computational protein design: folding and stability of nine completely redesigned globular proteins. , 2003, Journal of molecular biology.

[70]  D. Baker,et al.  Design of a Novel Globular Protein Fold with Atomic-Level Accuracy , 2003, Science.

[71]  P. Harbury,et al.  Automated design of specificity in molecular recognition , 2003, Nature Structural Biology.

[72]  Celia A. Schiffer,et al.  Structural Basis for Coevolution of a Human Immunodeficiency Virus Type 1 Nucleocapsid-p1 Cleavage Site with a V82A Drug-Resistant Mutation in Viral Protease , 2004, Journal of Virology.

[73]  Celia A. Schiffer,et al.  Structural and Thermodynamic Basis for the Binding of TMC114, a Next-Generation Human Immunodeficiency Virus Type 1 Protease Inhibitor , 2004, Journal of Virology.

[74]  B. Tidor,et al.  Escherichia coli glutaminyl-tRNA synthetase is electrostatically optimized for binding of its cognate substrates. , 2004, Journal of molecular biology.

[75]  Antiviral drugs: Inside the envelope , 2004, Nature Reviews Drug Discovery.

[76]  James Andrew McCammon,et al.  Charge optimization of the interface between protein kinases and their ligands , 2004, J. Comput. Chem..

[77]  Celia A Schiffer,et al.  Combating susceptibility to drug resistance: lessons from HIV-1 protease. , 2004, Chemistry & biology.

[78]  Peter Kirkpatrick Antiviral drugs: Inside the envelope , 2004, Nature Reviews Drug Discovery.

[79]  C. Schiffer,et al.  Combating Susceptibility to Drug ResistanceLessons from HIV-1 Protease , 2004 .

[80]  Roger A. Jones,et al.  Structures of HIV-1 RT–DNA complexes before and after incorporation of the anti-AIDS drug tenofovir , 2004, Nature Structural &Molecular Biology.

[81]  Bruce Tidor,et al.  Design of improved protein inhibitors of HIV‐1 cell entry: Optimization of electrostatic interactions at the binding interface , 2005, Proteins.

[82]  Irene T Weber,et al.  Molecular basis for substrate recognition and drug resistance from 1.1 to 1.6 Å resolution crystal structures of HIV‐1 protease mutants with substrate analogs , 2005, The FEBS journal.

[83]  Michael K Gilson Sensitivity Analysis and Charge-Optimization for Flexible Ligands:  Applicability to Lead Optimization. , 2006, Journal of chemical theory and computation.

[84]  D. Baker,et al.  Computational design of a new hydrogen bond network and at least a 300-fold specificity switch at a protein-protein interface. , 2006, Journal of molecular biology.