Search

[This corrects the article DOI: 10.1371/journal.pgen.1005310.].

Viruses often encode proteins with multiple functions due to their compact genomes. Existing approaches to identify functional residues largely rely on sequence conservation analysis. Inferring functional residues from sequence conservation can produce false positives, in which the conserved residues are functionally silent, or false negatives, where functional residues are not identified since they are species-specific and therefore non-conserved. Furthermore, the tedious process of constructing and analyzing individual mutations limits the number of residues that can be examined in a single study. Here, we developed a systematic approach to identify the functional residues of a viral protein by coupling experimental fitness profiling with protein stability prediction using the influenza virus polymerase PA subunit as the target protein. We identified a significant number of functional residues that were influenza type-specific and were evolutionarily non-conserved among different influenza types. Our results indicate that type-specific functional residues are prevalent and may not otherwise be identified by sequence conservation analysis alone. More importantly, this technique can be adapted to any viral (and potentially non-viral) protein where structural information is available.

Widely used chemical genetic screens have greatly facilitated the identification of many antiviral agents. However, the regions of interaction and inhibitory mechanisms of many therapeutic candidates have yet to be elucidated. Previous chemical screens identified Daclatasvir (BMS-790052) as a potent nonstructural protein 5A (NS5A) inhibitor for Hepatitis C virus (HCV) infection with an unclear inhibitory mechanism. Here we have developed a quantitative high-resolution genetic (qHRG) approach to systematically map the drug-protein interactions between Daclatasvir and NS5A and profile genetic barriers to Daclatasvir resistance. We implemented saturation mutagenesis in combination with next-generation sequencing technology to systematically quantify the effect of every possible amino acid substitution in the drug-targeted region (domain IA of NS5A) on replication fitness and sensitivity to Daclatasvir. This enabled determination of the residues governing drug-protein interactions. The relative fitness and drug sensitivity profiles also provide a comprehensive reference of the genetic barriers for all possible single amino acid changes during viral evolution, which we utilized to predict clinical outcomes using mathematical models. We envision that this high-resolution profiling methodology will be useful for next-generation drug development to select drugs with higher fitness costs to resistance, and also for informing the rational use of drugs based on viral variant spectra from patients.

Trade-offs between throughput, read length, and error rates in high-throughput sequencing limit certain applications such as monitoring viral quasispecies. Here, we describe a molecular-based tag linkage method that allows assemblage of short sequence reads into long DNA fragments. It enables haplotype phasing with high accuracy and sensitivity to interrogate individual viral sequences in a quasispecies. This approach is demonstrated to deduce ∼ 2000 unique 1.3 kb viral sequences from HIV-1 quasispecies in vivo and after passaging ex vivo with a detection limit of ∼ 0.005% to ∼ 0.001%. Reproducibility of the method is validated quantitatively and qualitatively by a technical replicate. This approach can improve monitoring of the genetic architecture and evolution dynamics in any quasispecies population.

HIV-1 integrase (IN) is an important target for contemporary antiretroviral drug design research. Historically, efforts at inactivating the enzyme have focused upon blocking its active site. However, it has become apparent that new classes of allosteric inhibitors will be necessary to advance the antiretroviral field in light of the emergence of viral strains resistant to contemporary clinically used IN drugs. In this study we have characterized the importance of a close network of IN residues, distant from the active site, as important for the obligatory multimerization of the enzyme and viral replication as a whole. Specifically, we have determined that the configuration of six residues within a highly symmetrical region at the IN dimerization interface, composed of a four-tiered aromatic interaction flanked by two salt bridges, significantly contributes to proper HIV-1 replication. Additionally, we have utilized a quantitative luminescence assay to examine IN oligomerization and have determined that there is a very low tolerance for amino acid substitutions along this region. Even conservative residue substitutions negatively impacted IN multimerization, resulting in an inactive viral enzyme and a non-replicative virus. We have shown that there is a very low tolerance for amino acid variation at the symmetrical dimeric interface region characterized in this study, and therefore drugs designed to target the amino acid network detailed here could be expected to yield a significantly reduced number of drug-resistant escape mutations compared to contemporary clinically-evaluated antiretrovirals.

Two decades of intensive research efforts have led to the discovery of a large number of HIV-1 integrase (IN) inhibitors. Recently, the United States Food and Drug Administration (US FDA) approved MK-0518, or raltegravir ( 1), as the first IN inhibitor for HIV/AIDS treatment. Growing clinical evidence also demonstrates that the emergence of HIV-1 virus strains bearing IN amino acid substitutions that confer resistance to IN inhibitors is inevitable. The discovery of second generation inhibitors with potency against viral strains bearing drug resistant IN substitutions is necessary for ongoing effective treatment of viral infections. We generated common feature pharmacophore hypotheses using a training set of quinolone 3-carboxylic acid IN inhibitors, including the clinical candidate GS-9137 ( 2). A database search of small molecules using the quinolone 3-carboxylic acid pharmacophore model, followed by in vitro evaluation of selected hits in an assay specific to IN, resulted in the discovery of potential leads with diverse structural scaffolds useful for the design of second generation IN inhibitors. [ABSTRACT FROM AUTHOR]