Characterizing the Diverse Mutational Pathways Associated with R5-Tropic Maraviroc Resistance: HIV-1 That Uses the Drug-Bound CCR5 Coreceptor

ABSTRACT Entry inhibitors represent a potent class of antiretroviral drugs that target a host cell protein, CCR5, an HIV-1 entry coreceptor, and not viral protein. Lack of sensitivity can occur due to preexisting virus that uses the CXCR4 coreceptor, while true resistance occurs through viral adaptation to use a drug-bound CCR5 coreceptor. To understand this R5 resistance pathway, we analyzed >500 envelope protein sequences and phenotypes from viruses of 20 patients from the clinical trials MOTIVATE 1 and 2, in which treatment-experienced patients received maraviroc plus optimized background therapy. The resistant viral population was phylogenetically distinct and associated with a genetic bottleneck in each patient, consistent with de novo emergence of resistance. Recombination analysis showed that the C2-V3-C3 region tends to genotypically correspond to the recombinant's phenotype, indicating its primary importance in conferring resistance. Between patients, there was a notable lack of commonality in the specific sites conferring resistance, confirming the unusual nature of R5-tropic resistance. We used coevolutionary and positive-selection analyses to characterize the genotypic determinants of resistance and found that (i) there are complicated covariation networks, indicating frequent coevolutionary/compensatory changes in the context of protein structure; (ii) covarying sites under positive selection are enriched in resistant viruses; (iii) CD4 binding sites form part of a unique covariation network independent of the V3 loop; and (iv) the covariation network formed between the V3 loop and other regions of gp120 and gp41 intersects sites involved in glycosylation and protein secretion. These results demonstrate that while envelope sequence mutations are the key to conferring maraviroc resistance, the specific changes involved are context dependent and thus inherently unpredictable. IMPORTANCE The entry inhibitor drug maraviroc makes the cell coreceptor CCR5 unavailable for use by HIV-1 and is now used in combination antiretroviral therapy. Treatment failure with drug-resistant virus is particularly interesting because it tends to be rare, with lack of sensitivity usually associated with the presence of CXCR4-using virus (CXCR4 is the main alternative coreceptor HIV-1 uses, in addition to CD4). We analyzed envelope sequences from HIV-1, obtained from 20 patients who enrolled in maraviroc clinical trials and experienced treatment failure, without detection of CXCR4-using virus. Evolutionary analysis was employed to identify molecular changes that confer maraviroc resistance. We found that in these individuals, resistant viruses form a distinct population that evolved once and was successful as a result of drug pressure. Further evolutionary analysis placed the complex network of interdependent mutational changes into functional groups that help explain the impediments to the emergence of maraviroc-associated R5 drug resistance.

[1]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[2]  M. Churchill,et al.  HIV-1 Escape from the CCR5 Antagonist Maraviroc Associated with an Altered and Less-Efficient Mechanism of gp120-CCR5 Engagement That Attenuates Macrophage Tropism , 2011, Journal of Virology.

[3]  J. Margolick,et al.  Improved Coreceptor Usage Prediction and GenotypicMonitoring of R5-to-X4 Transition by Motif Analysis of HumanImmunodeficiency Virus Type 1 env V3 LoopSequences , 2003, Journal of Virology.

[4]  John P. Moore,et al.  Resistance to CCR5 inhibitors caused by sequence changes in the fusion peptide of HIV-1 gp41 , 2009, Proceedings of the National Academy of Sciences.

[5]  G. von Heijne The signal peptide. , 1990, The Journal of membrane biology.

[6]  Y. Levy,et al.  Effect of glycosylation on protein folding: A close look at thermodynamic stabilization , 2008, Proceedings of the National Academy of Sciences.

[7]  Simon A. A. Travers,et al.  The Influence of N-Linked Glycans on the Molecular Dynamics of the HIV-1 gp120 V3 Loop , 2013, PloS one.

[8]  S. Inouye,et al.  Role of positive charge on the amino-terminal region of the signal peptide in protein secretion across the membrane. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[9]  M. Alizon,et al.  Cooperation of the V1/V2 and V3 Domains of Human Immunodeficiency Virus Type 1 gp120 for Interaction with the CXCR4 Receptor , 2001, Journal of Virology.

[10]  Madolyn L. MacDonald,et al.  Rapid Evolutionary Dynamics of Structural Disorder as a Potential Driving Force for Biological Divergence in Flaviviruses , 2013, Genome biology and evolution.

[11]  G. Ciaramella,et al.  Reduced Maximal Inhibition in Phenotypic Susceptibility Assays Indicates that Viral Strains Resistant to the CCR5 Antagonist Maraviroc Utilize Inhibitor-Bound Receptor for Entry , 2006, Journal of Virology.

[12]  G. Vonheijne The signal peptide. , 1990 .

[13]  Thomas Lengauer,et al.  Recco: recombination analysis using cost optimization , 2006, Bioinform..

[14]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[15]  O. Pybus,et al.  Bayesian coalescent inference of past population dynamics from molecular sequences. , 2005, Molecular biology and evolution.

[16]  B. Murrell,et al.  RDP4: Detection and analysis of recombination patterns in virus genomes , 2015, Virus evolution.

[17]  Mike Youle,et al.  Emergence of CXCR4-Using Human Immunodeficiency Virus Type 1 (HIV-1) Variants in a Minority of HIV-1-Infected Patients following Treatment with the CCR5 Antagonist Maraviroc Is from a Pretreatment CXCR4-Using Virus Reservoir , 2006, Journal of Virology.

[18]  D. Y. Thomas,et al.  Control of expression, glycosylation, and secretion of HIV-1 gp120 by homologous and heterologous signal sequences. , 1994, Virology.

[19]  M. Suchard,et al.  Bayesian Phylogenetics with BEAUti and the BEAST 1.7 , 2012, Molecular biology and evolution.

[20]  K. Crandall,et al.  The Effect of Recombination on the Accuracy of Phylogeny Estimation , 2002, Journal of Molecular Evolution.

[21]  J. Montaner,et al.  Enfuvirtide, an HIV-1 fusion inhibitor, for drug-resistant HIV infection in North and South America. , 2003, The New England journal of medicine.

[22]  W. Atchley,et al.  Solving the protein sequence metric problem. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Yu Xue,et al.  DOG 1.0: illustrator of protein domain structures , 2009, Cell Research.

[24]  Bert van den Berg,et al.  X-ray structure of a protein-conducting channel , 2004, Nature.

[25]  A. Heredia,et al.  An Alteration of Human Immunodeficiency Virus gp41 Leads to Reduced CCR5 Dependence and CD4 Independence , 2008, Journal of Virology.

[26]  Elias Lolis,et al.  Structure, function, and inhibition of chemokines. , 2002, Annual review of pharmacology and toxicology.

[27]  F. Bushman,et al.  A Maraviroc-Resistant HIV-1 with Narrow Cross-Resistance to Other CCR5 Antagonists Depends on both N-Terminal and Extracellular Loop Domains of Drug-Bound CCR5 , 2010, Journal of Virology.

[28]  Trey Ideker,et al.  Cytoscape 2.8: new features for data integration and network visualization , 2010, Bioinform..

[29]  M. Aebi,et al.  Mechanisms and principles of N-linked protein glycosylation. , 2011, Current opinion in structural biology.

[30]  Sergei L. Kosakovsky Pond,et al.  Evolutionary Interactions between N-Linked Glycosylation Sites in the HIV-1 Envelope , 2006, PLoS Comput. Biol..

[31]  John P. Moore,et al.  Two HIV-1 Variants Resistant to Small Molecule CCR5 Inhibitors Differ in How They Use CCR5 for Entry , 2009, PLoS pathogens.

[32]  M. Fares,et al.  IDENTIFYING COEVOLUTIONARY PATTERNS IN HUMAN LEUKOCYTE ANTIGEN (HLA) MOLECULES , 2009, Evolution; international journal of organic evolution.

[33]  O. Gascuel,et al.  A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. , 2003, Systematic biology.

[34]  A. J. Parodi,et al.  Protein glucosylation and its role in protein folding. , 2000, Annual review of biochemistry.

[35]  R. Doms,et al.  HIV type 1 from a patient with baseline resistance to CCR5 antagonists uses drug-bound receptor for entry. , 2010, AIDS research and human retroviruses.

[36]  Mario A. Fares,et al.  CAPS: coevolution analysis using protein sequences , 2006, Bioinform..

[37]  Sergei L. Kosakovsky Pond,et al.  HyPhy: hypothesis testing using phylogenies , 2005, Bioinform..

[38]  E. Poveda,et al.  Primary resistance to maraviroc in a large set of R5-V3 viral sequences from HIV-1-infected patients. , 2010, The Journal of antimicrobial chemotherapy.

[39]  N. Goldman,et al.  Codon-substitution models for heterogeneous selection pressure at amino acid sites. , 2000, Genetics.

[40]  Simon A. A. Travers,et al.  A study of the coevolutionary patterns operating within the env gene of the HIV-1 group M subtypes. , 2007, Molecular biology and evolution.

[41]  E. Arts,et al.  HIV-1 Resistance to Maraviroc Conferred by a CD4 Binding Site Mutation in the Envelope Glycoprotein gp120 , 2012, Journal of Virology.

[42]  David A. Price,et al.  Maraviroc (UK-427,857), a Potent, Orally Bioavailable, and Selective Small-Molecule Inhibitor of Chemokine Receptor CCR5 with Broad-Spectrum Anti-Human Immunodeficiency Virus Type 1 Activity , 2005, Antimicrobial Agents and Chemotherapy.

[43]  Chi‐Huey Wong,et al.  The core trisaccharide of an N-linked glycoprotein intrinsically accelerates folding and enhances stability , 2009, Proceedings of the National Academy of Sciences.

[44]  D. Weissman,et al.  Quantification of CD4, CCR5, and CXCR4 levels on lymphocyte subsets, dendritic cells, and differentially conditioned monocyte-derived macrophages. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Steven M. Wolinsky,et al.  Changes in the V 3 region of gp 120 contribute to unusually broad coreceptor usage of an HIV-1 isolate from a CCR 5 Delta 32 heterozygote , 2007 .

[46]  John R Mascola,et al.  Changes in the V3 region of gp120 contribute to unusually broad coreceptor usage of an HIV-1 isolate from a CCR5 Delta32 heterozygote. , 2007, Virology.

[47]  D. Kuritzkes,et al.  HIV-1 Clinical Isolates Resistant to CCR5 Antagonists Exhibit Delayed Entry Kinetics That Are Corrected in the Presence of Drug , 2011, Journal of Virology.

[48]  G. Sapiro,et al.  Molecular architecture of native HIV-1 gp120 trimers , 2008, Nature.

[49]  Chi-Huey Wong,et al.  Protein Native-State Stabilization by Placing Aromatic Side Chains in N-Glycosylated Reverse Turns , 2011, Science.

[50]  N. Haigwood,et al.  N-glycosylation of HIV-gp120 may constrain recognition by T lymphocytes. , 1991, Journal of immunology.

[51]  R. Swanstrom,et al.  Improved success of phenotype prediction of the human immunodeficiency virus type 1 from envelope variable loop 3 sequence using neural networks. , 2001, Virology.

[52]  Joost Schymkowitz,et al.  The stability effects of protein mutations appear to be universally distributed. , 2007, Journal of molecular biology.

[53]  J. Moult,et al.  Loss of protein structure stability as a major causative factor in monogenic disease. , 2005, Journal of molecular biology.

[54]  Feng Gao,et al.  Recurrent Signature Patterns in HIV-1 B Clade Envelope Glycoproteins Associated with either Early or Chronic Infections , 2011, PLoS pathogens.

[55]  Ramón Doallo,et al.  ProtTest 3: fast selection of best-fit models of protein evolution , 2011, Bioinform..

[56]  Martin A. Nowak,et al.  Antibody neutralization and escape by HIV-1 , 2003, Nature.

[57]  Simon A. A. Travers,et al.  A Novel Method for Detecting Intramolecular Coevolution: Adding a Further Dimension to Selective Constraints Analyses , 2006, Genetics.

[58]  I. Braakman,et al.  Folding of the human immunodeficiency virus type 1 envelope glycoprotein in the endoplasmic reticulum. , 2001, Biochimie.

[59]  Jaap Goudsmit,et al.  N-Linked Glycosylation of the HIV Type-1 gp120 Envelope Glycoprotein as a Major Determinant of CCR5 and CXCR4 Coreceptor Utilization* , 2001, The Journal of Biological Chemistry.

[60]  Anders Karlsson,et al.  Maraviroc for previously treated patients with R5 HIV-1 infection. , 2008, The New England journal of medicine.

[61]  A. Perelson,et al.  A Signature in HIV-1 Envelope Leader Peptide Associated with Transition from Acute to Chronic Infection Impacts Envelope Processing and Infectivity , 2011, PloS one.

[62]  R. Dwek,et al.  Exploiting the defensive sugars of HIV-1 for drug and vaccine design , 2007, Nature.

[63]  Alexandros Stamatakis,et al.  RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models , 2006, Bioinform..

[64]  D. Weissman,et al.  Quantification of CD 4 , CCR 5 , and CXCR 4 levels on lymphocyte subsets , dendritic cells , and differentially conditioned monocyte-derived macrophages , 1999 .

[65]  M. Sitbon,et al.  Functional tolerance of the human immunodeficiency virus type 1 envelope signal peptide to mutations in the amino-terminal and hydrophobic regions , 1992, Journal of virology.

[66]  Geoffrey J. Barton,et al.  Jalview Version 2—a multiple sequence alignment editor and analysis workbench , 2009, Bioinform..

[67]  David L. Robertson,et al.  The Evolutionary Analysis of Emerging Low Frequency HIV-1 CXCR4 Using Variants through Time—An Ultra-Deep Approach , 2010, PLoS Comput. Biol..

[68]  S. Zolla-Pazner,et al.  Structural basis for coreceptor selectivity by the HIV type 1 V3 loop. , 2007, AIDS research and human retroviruses.

[69]  Christopher J. Oldfield,et al.  Do viral proteins possess unique biophysical features? , 2009, Trends in biochemical sciences.