Three mutants of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (V106A, V179D, and Y181C), which occur in clinical isolates and confer resistance to nonnucleoside reverse transcriptase inhibitors (NNRTIs), were analyzed for RNA-and DNA-dependent DNA polymerization and RNase H cleavage. All mutants demonstrated processivities of polymerization that were indistinguishable from wild-type enzyme under conditions in which deoxynucleoside triphosphates were not limiting. The V106A reverse transcriptase demonstrated a three-to fourfold slowing of both DNA 3-end-directed and RNA 5-end-directed RNase H cleavage relative to both wild-type and V179D enzymes, similar to what was observed for P236L in a previously published study (P. Gerondelis et al., J. Virol. 73:5803-5813, 1999). In contrast, the Y181C reverse transcriptase demonstrated a selective acceleration of the secondary RNase H cleavage step during both modes of RNase H cleavage. The relative replication fitness of these mutants in H9 cells was assessed in parallel infections as well as in growth competition experiments. Of the NNRTI-resistant mutants, V179D was more fit than Y181C, and both of these mutants were more fit than V106A, which demonstrated the greatest reduction in RNase H cleavage. These findings, in combination with results from previous work, suggest that abnormalities in RNase H cleavage are a common characteristic of HIV-1 mutants resistant to NNRTIs and that combined reductions in the rates of DNA 3-end-and RNA 5-end-directed cleavages are associated with significant reductions in the replication fitness of HIV-1.Infection with human immunodeficiency virus (HIV) is the cause of AIDS and affects over 30 million people worldwide (64). The primary targets of therapy for HIV infection include the viral protease and reverse transcriptase (RT). HIV type 1 (HIV-1) RT is a heterodimer consisting of 66-and 51-kDa subunits (p66 and p51, respectively) (3). p66 contains both the polymerase and the RNase H active sites of the enzyme (34, 37, 39). The RNase H domain is present in the carboxy-terminal third of p66. Although p51 is derived from p66 by proteolytic cleavage, it assumes a very different tertiary structure and does not contain a catalytic site (37, 39). The function of p51 is not known, but it may play a role in binding the tRNA 3Lys -template complex (3, 39) and in maintaining the structural integrity of the heterodimer (1).RNase H cleavage is essential for HIV-1 replication (61; for a review see reference 11). Two modes of RNase H cleavage have been described (Fig. 1). "Polymerase-dependent" cleavage is thought to occur in concert with DNA polymerization to degrade the genomic RNA during minus strand DNA synthesis (26,46). The position of the primary DNA 3Ј-end-directed cleavage occurs 15 to 18 nucleotides (nt) from the recessed 3Ј end of the DNA (26, 33); we have referred to this mode of cleavage as DNA 3Ј-end-directed RNase H cleavage. A second mode of RNase H cleavage occurs independently of DNA polymerization. The position ...

The reliable detection, sizing, and sorting of viruses and nanoparticles is important for biosensing, environmental monitoring, and quality control. Here we introduce an optical detection scheme for the real-time and label-free detection and recognition of single viruses and larger proteins. The method makes use of nanofluidic channels in combination with optical interferometry. Elastically scattered light from single viruses traversing a stationary laser focus is detected with a differential heterodyne interferometer and the resulting signal allows single viruses to be characterized individually. Heterodyne detection eliminates phase variations due to different particle trajectories, thus improving the recognition accuracy as compared to standard optical interferometry. We demonstrate the practicality of our approach by resolving nanoparticles of various sizes, and detecting and recognizing different species of human viruses from a mixture. The detection system can be readily integrated into larger nanofluidic architectures for practical applications.

SummaryMany biological processes and systems can be described by a set of differential equation (DE) models. However, literature in statistical inference for DE models is very sparse. We propose statistical estimation, model selection, and multimodel averaging methods for HIV viral fitness experiments in vitro that can be described by a set of nonlinear ordinary differential equations (ODE). The parameter identifiability of the ODE models is also addressed. We apply the proposed methods and techniques to experimental data of viral fitness for HIV-1 mutant 103N. We expect that the proposed modeling and inference approaches for the DE models can be widely used for a variety of biomedical studies.

We evaluated the replication efficiency of the HIV reverse transcriptase (RT) mutants K103N, G190A, and G190S, which confer resistance to the non-nucleoside RT inhibitor efavirenz, using growth competition assays in cell culture. In the absence of efavirenz, the fitness hierarchy was G190S < G190A < K103N < wild-type. The fitness reduction of G190S relative to K103N was less evident at high efavirenz concentrations, although K103N still replicated more efficiently. Efficiency of RNase H cleavage and RNA-dependent DNA synthesis from tRNA(Lys, 3) correlated with relative fitness, in biochemical studies of mutant RTs. Presteady state and steady state polymerization assays using DNA primers detected no abnormalities. This work is consistent with previous studies demonstrating that initiation of viral DNA synthesis is reduced in mutants with slowed RNase H cleavage, and suggests that both abnormalities contribute to the replication defect of these mutants. It also suggests that high concentrations of efavirenz are unlikely to favor the selection of G190S clinically.

Growth competition assays have been developed to quantify the relative fitnesses of human immunodeficiency virus (HIV-1) mutants. In this article we develop mathematical models to describe viral/cellular dynamic interactions in the assay experiment, from which new competitive fitness indices or parameters are defined. These indices include the log fitness ratio (LFR), the log relative fitness (LRF), and the production rate ratio (PRR). From the population genetics perspective, we clarify the confusion and correct the inconsistency in the definition of relative fitness in the literature of HIV-1 viral fitness. The LFR and LRF are easier to estimate from the experimental data than the PRR, which was misleadingly defined as the relative fitness in recent HIV-1 research literature. Calculation and estimation methods based on two data points and multiple data points were proposed and were carefully studied. In particular, we suggest using both standard linear regression (method of least squares) and a measurement error model approach for more-accurate estimates of competitive fitness parameters from multiple data points. The developed methodologies are generally applicable to any growth competition assays. A user-friendly computational tool also has been developed and is publicly available on the World Wide Web at http://www.urmc.rochester.edu/bstools/vfitness/virusfitness.htm.

Growth competition assays have been developed to quantify the relative fitness of HIV-1 mutants. In this article, we develop mathematical models to describe viral/cellular dynamic interactions in the assay system from which the competitive fitness indices or parameters are defined. In our previous HIV-viral fitness experiments, the concentration of uninfected target cells was assumed to be constant (Wu et al., 2006). But this may not be true in some experiments. In addition, dual infection may frequently occur in viral fitness experiments and may not be ignorable. Here, we relax these two assumptions and extend our earlier viral fitness model (Wu et al., 2006). The resulting models then become nonlinear ODE systems for which closed-form solutions are not achievable. In the new model, the viral relative fitness is a function of time since it depends on the target cell concentration. First, we studied the structure identifiability of the nonlinear ODE models. The identifiability analysis showed that all parameters in the proposed models are identifiable from the flow-cytometry-based experimental data that we collected. We then employed a global optimization approach (the differential evolution algorithm) to directly estimate the kinetic parameters as well as the relative fitness index in the nonlinear ODE models using nonlinear least square regression based on the experimental data. Practical identifiability was investigated via Monte Carlo simulations.

SUMMARY The relative fitness of a variant, according to population genetics theory, is that variant's relative contribution to successive generations. Most drug-resistant human immunodeficiency virus type 1 (HIV-1) variants have reduced replication fitness, but at least some of these deficits can be compensated for by the accumulation of second-site mutations. HIV-1 replication fitness also appears to influence the likelihood of a drug-resistant mutant emerging during treatment failure and is postulated to influence clinical outcomes. A variety of assays are available to measure HIV-1 replication fitness in cell culture; however, there is no agreement regarding which assays best correlate with clinical outcomes. A major limitation is that there is no high-throughput assay that incorporates an internal reference strain as a control and utilizes intact virus isolates. Some retrospective studies have demonstrated statistically significant correlations between HIV-1 replication fitness and clinical outcomes in some patient populations. However, different studies disagree as to which clinical outcomes are most closely associated with fitness. This may be in part due to assay design, sample size limitations, and differences in patient populations. In addition, the strength of the correlations between fitness and clinical outcomes is modest, suggesting that, at present, it would be difficult to utilize these assays for clinical management.

We designed a cell culture-based system to test the hypothesis that recombination events during HIV-1 replication would be more frequent near the dimerization initiation sequence (DIS). A 459-bp region spanning the DIS through the 5'-end of gag was sequenced and analyzed to determine the frequency and distribution of crossover sites. We found a strong preference for recombination events occurring within a 112-nt-long region encompassing the gag AUG (64% of crossovers occurred in this region, compared to 10-14% in surrounding regions with similar lengths). Surprisingly, the region immediately surrounding the DIS was not a preferred site of recombination. Analysis of recombination events using RNA templates transcribed in vitro revealed a preference for crossover sites at the start of the gag coding region, similar to that observed in cell culture. This recombinogenic region was unusually G-rich and promoted extensive pausing by RT in vitro. Template features that induce RT pausing very likely contribute to the observed template switching events in gag during minus-strand synthesis. The region in gag that was a preferred site for recombination also had an approximately 2-fold higher mutation frequency compared to the rest of the region sequenced, but mutations were no more common in recombinant compared to non-recombinant clones, suggesting that recombination events were not mutagenic.