Human immunodeficiency virus reverse transcriptase displays dramatically higher fidelity under physiological magnesium conditions in vitro

Mechanistic Interplay between HIV-1 Reverse Transcriptase Enzyme Kinetics and Host SAMHD1 Protein: Viral Myeloid-Cell Tropism and Genomic Mutagenesis

Human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) has been the primary interest among studies on antiviral discovery, viral replication kinetics, drug resistance, and viral evolution. Following infection and entry into target cells, the HIV-1 core disassembles, and the viral RT concomitantly converts the viral RNA into double-stranded proviral DNA, which is integrated into the host genome. The successful completion of the viral life cycle highly depends on the enzymatic DNA polymerase activity of RT. Furthermore, HIV-1 RT has long been known as an error-prone DNA polymerase due to its lack of proofreading exonuclease properties. Indeed, the low fidelity of HIV-1 RT has been considered as one of the key factors in the uniquely high rate of mutagenesis of HIV-1, which leads to efficient viral escape from immune and therapeutic antiviral selective pressures. Interestingly, a series of studies on the replication kinetics of HIV-1 in non-dividing myeloid cells and myeloid specific host restriction factor, SAM domain, and HD domain-containing protein, SAMHD1, suggest that the myeloid cell tropism and high rate of mutagenesis of HIV-1 are mechanistically connected. Here, we review not only HIV-1 RT as a key antiviral target, but also potential evolutionary and mechanistic crosstalk among the unique enzymatic features of HIV-1 RT, the replication kinetics of HIV-1, cell tropism, viral genetic mutation, and host SAMHD1 protein.

Deep Sequencing Analysis of Individual HIV-1 Proviruses Reveals Frequent Asymmetric Long Terminal Repeats

Effective strategies to eliminate human immunodeficiency virus type 1 (HIV-1) reservoirs are likely to require more thorough characterizations of proviruses that persist on antiretroviral therapy (ART). The rarity of infected CD4⁺ T-cells and related technical challenges have limited the characterization of integrated proviruses. Current approaches using next-generation sequencing can be inefficient and limited sequencing depth can make it difficult to link proviral sequences to their respective integration sites. Here, we report on an efficient method by which HIV-1 proviruses and their sites of integration are amplified and sequenced. Across five HIV-1-positive individuals on clinically effective ART, a median of 41.2% (n = 88 of 209) of amplifications yielded near-full-length proviruses and their 5'-host-virus junctions containing a median of 430 bp (range, 18 to 1,363 bp) of flanking host sequence. Unexpectedly, 29.5% (n = 26 of 88) of the sequenced proviruses had structural asymmetries between the 5' and 3' long terminal repeats (LTRs), commonly in the form of major 3' deletions. Sequence-intact proviruses were detected in 3 of 5 donors, and infected CD4⁺ T-cell clones were detected in 4 of 5 donors. The accuracy of the method was validated by amplifying and sequencing full-length proviruses and flanking host sequences directly from peripheral blood mononuclear cell DNA. The individual proviral sequencing assay (IPSA) described here can provide an accurate, in-depth, and longitudinal characterization of HIV-1 proviruses that persist on ART, which is important for targeting proviruses for elimination and assessing the impact of interventions designed to eradicate HIV-1. IMPORTANCE The integration of human immunodeficiency virus type 1 (HIV-1) into chromosomal DNA establishes the long-term persistence of HIV-1 as proviruses despite effective antiretroviral therapy (ART). Characterizing proviruses is difficult because of their rarity in individuals on long-term suppressive ART, their highly polymorphic sequences and genetic structures, and the need for efficient amplification and sequencing of the provirus and its integration site. Here, we describe a novel, integrated, two-step method (individual proviral sequencing assay [IPSA]) that amplifies the host-virus junction and the full-length provirus except for the last 69 bp of the 3' long terminal repeat (LTR). Using this method, we identified the integration sites of proviruses, including those that are sequence intact and replication competent or defective. Importantly, this new method identified previously unreported asymmetries between LTRs that have implications for how proviruses are detected and quantified. The IPSA method reported is unaffected by LTR asymmetries, permitting a more accurate and comprehensive characterization of the proviral landscape.

Genomic signatures of protease and reverse transcriptase genes from HIV-1 subtype C isolated from first-line ART patients in India

Genomic signatures of the protease and reverse transcriptase gene of HIV-1 from HIV infected North Indian patients who were under ART from 1 to ≤ 7 years were analyzed. The DNA from plasma samples of 9 patients and RNA from 57 patients were isolated and subjected to amplification for the protease and reverse transcriptase gene of HIV-1 subtype C. Then sequencing was carried out following the WHO dried blood spot protocol. The drug resistance mutation patterns were analyzed using the HIV Drug Resistance Database, Stanford University, USA. Lamivudine-associated drug-resistance mutations such as M184V/M184I, nevirapine-associated drug resistance mutations Y181C and H221Y, and efavirenz-associated drug resistance mutations M230I were observed in reverse transcriptase gene of archived DNA of two HIV-1 infected patients. No mutation was observed in the remaining 7 patients. Various computational tools and websites like viral epidemiological signature pattern analysis (VESPA), hyper mutation, SNAP version 2.1.1, and entropy were utilized for the analysis of the signature pattern of amino acids, hyper mutation, selection pressure, and Shannon entropy in the protease and reverse transcriptase gene sequences of the 9 archived DNA, 56 protease gene and 51 reverse transcriptase gene from the HIV-1 DNA amplified sequences of RNA. The HIV-1 Subtype-C (Gene bank accession number: AB023804) and first isolate HXB2 (Gene bank accession number: K03455.1) was taken as reference sequence. The signature amino acid sequences were identified in the protease and reverse transcriptase gene, no hyper mutation, highest entropy was marked in the amino acid positions and synonymous to non-synonymous nucleotide ratio was calculated in the protease and reverse transcriptase gene of 9 archived DNA sequences, 56 protease and 51 reverse transcriptase gene sequences of HIV-1 Subtype C isolates.

Processing of RNA Containing 8-Oxo-7,8-Dihydroguanosine (8-oxoG) by the Exoribonuclease Xrn-1

Understanding how oxidatively damaged RNA is handled intracellularly is of relevance due to the link between oxidized RNA and the progression/development of some diseases as well as aging. Among the ribonucleases responsible for the decay of modified (chemically or naturally) RNA is the exonuclease Xrn-1, a processive enzyme that catalyzes the hydrolysis of 5′-phosphorylated RNA in a 5′→3′ direction. We set out to explore the reactivity of this exonuclease towards oligonucleotides (ONs, 20-nt to 30-nt long) of RNA containing 8-oxo-7,8-dihydroguanosine (8-oxoG), obtained via solid-phase synthesis. The results show that Xrn-1 stalled at sites containing 8-oxoG, evidenced by the presence of a slower moving band (via electrophoretic analyses) than that observed for the canonical analogue. The observed fragment(s) were characterized via PAGE and MALDI-TOF to confirm that the oligonucleotide fragment(s) contained a 5′-phosphorylated 8-oxoG. Furthermore, the yields for this stalling varied from app. 5–30% with 8-oxoG located at different positions and in different sequences. To gain a better understanding of the decreased nuclease efficiency, we probed: 1) H-bonding and spatial constraints; 2) anti-syn conformational changes; 3) concentration of divalent cation; and 4) secondary structure. This was carried out by introducing methylated or brominated purines (m¹G, m^6,6A, or 8-BrG), probing varying [Mg²⁺], and using circular dichroism (CD) to explore the formation of structured RNA. It was determined that spatial constraints imposed by conformational changes around the glycosidic bond may be partially responsible for stalling, however, the results do not fully explain some of the observed higher stalling yields. We hypothesize that altered π-π stacking along with induced H-bonding interactions between 8-oxoG and residues within the binding site may also play a role in the decreased Xrn-1 efficiency. Overall, these observations suggest that other factors, yet to be discovered/established, are likely to contribute to the decay of oxidized RNA. In addition, Xrn-1 degraded RNA containing m¹G, and stalled mildly at sites where it encountered m^6,6A, or 8-BrG, which is of particular interest given that the former two are naturally occurring modifications.

Physiological magnesium concentrations increase fidelity of diverse reverse transcriptases from HIV-1, HIV-2, and foamy virus, but not MuLV or AMV

Reverse transcriptases (RTs) are typically assayed using optimized Mg²⁺ concentrations (~5-10 mM) several-fold higher than physiological cellular free Mg²⁺ (~0.5 mM). Recent analyses demonstrated that HIV-1, but not Moloney murine leukaemia (MuLV) or avain myeloblastosis (AMV) virus RTs has higher fidelity in low Mg²⁺. In the current report, lacZα-based α-complementation assays were used to measure the fidelity of several RTs including HIV-1 (subtype B and A/E), several drug-resistant HIV-1 derivatives, HIV-2, and prototype foamy virus (PFV), all which showed higher fidelity using physiological Mg²⁺, while MuLV and AMV RTs demonstrated equivalent fidelity in low and high Mg²⁺. In 0.5 mM Mg²⁺, all RTs demonstrated approximately equal fidelity, except for PFV which showed higher fidelity. A Next Generation Sequencing (NGS) approach that used barcoding to determine mutation profiles was used to examine the types of mutations made by HIV-1 RT (type B) in low (0.5 mM) and high (6 mM) Mg²⁺ on a lacZα template. Unlike α-complementation assays which are dependent on LacZα activity, the NGS assay scores mutations at all positions and of every type. Consistent with α-complementation assays, a ~four-fold increase in mutations was observed in high Mg²⁺. These findings help explain why HIV-1 RT displays lower fidelity in vitro (with high Mg²⁺ concentrations) than other RTs (e.g. MuLV and AMV), yet cellular fidelity for these viruses is comparable. Establishing in vitro conditions that accurately represent RT's activity in cells is pivotal to determining the contribution of RT and other factors to the mutation profile observed with HIV-1.

Witnessing the structural evolution of an RNA enzyme

An RNA polymerase ribozyme that has been the subject of extensive directed evolution efforts has attained the ability to synthesize complex functional RNAs, including a full-length copy of its own evolutionary ancestor. During the course of evolution, the catalytic core of the ribozyme has undergone a major structural rearrangement, resulting in a novel tertiary structural element that lies in close proximity to the active site. Through a combination of site-directed mutagenesis, structural probing, and deep sequencing analysis, the trajectory of evolution was seen to involve the progressive stabilization of the new structure, which provides the basis for improved catalytic activity of the ribozyme. Multiple paths to the new structure were explored by the evolving population, converging upon a common solution. Tertiary structural remodeling of RNA is known to occur in nature, as evidenced by the phylogenetic analysis of extant organisms, but this type of structural innovation had not previously been observed in an experimental setting. Despite prior speculation that the catalytic core of the ribozyme had become trapped in a narrow local fitness optimum, the evolving population has broken through to a new fitness locale, raising the possibility that further improvement of polymerase activity may be achievable.

HIV Reverse Transcriptase Pre-Steady-State Kinetic Analysis of Chain Terminators and Translocation Inhibitors Reveals Interactions between Magnesium and Nucleotide 3'-OH

Deoxythymidine triphosphate analogues with various 3' substituents in the sugar ring (-OH (dTTP)), -H, -N₃, -NH₂, -F, -O-CH₃, no group (2',3'-didehydro-2',3'-dideoxythymidine triphosphate (d4TTP)), and those retaining the 3'-OH but with 4' additions (4'-C-methyl, 4'-C-ethyl) or sugar ring modifications (d-carba dTTP) were evaluated using pre-steady-state kinetics in low (0.5 mM) and high (6 mM) Mg²⁺ with HIV reverse transcriptase (RT). Analogues showed diminished observed incorporation rate constants (k _obs) compared to dTTP ranging from about 2-fold (3'-H, -N₃, and d4TTP with high Mg²⁺) to >10-fold (3'-NH₂ and 3'-F with low Mg²⁺), while 3'-O-CH₃ dTTP incorporated much slower than other analogues. Illustrating the importance of interactions between Mg²⁺ and the 3'-OH, k _obs using 5 μM dTTP and 0.5 mM Mg²⁺ was only modestly slower (1.6-fold) than with 6 mM Mg²⁺, while analogues with 3' alterations incorporated 2.8-5.1-fold slower in 0.5 mM Mg²⁺. In contrast, 4'-C-methyl and d-carba dTTP, which retain the 3'-OH, were not significantly affected by Mg²⁺. Consistent with these results, analogues with 3' modifications were better inhibitors in 6 versus 0.5 mM Mg²⁺. Equilibrium dissociation constant (K _D) and maximum incorporation rate (k _pol) determinations for dTTP and analogues lacking a 3'-OH indicated that low Mg²⁺ caused a several-fold greater reduction in k _pol with the analogues but did not significantly affect K _D, results consistent with a role for 3'-OH/Mg²⁺ interactions in catalysis rather than nucleotide binding. Overall, results emphasize the importance of previously unreported interactions between Mg²⁺ and the 3'-OH of the incoming nucleotide and suggest that inhibitors with 3'-OH groups may have advantages in low free Mg²⁺ in physiological settings.

Avoiding Drug Resistance in HIV Reverse Transcriptase

HIV reverse transcriptase (RT) is an enzyme that plays a major role in the replication cycle of HIV and has been a key target of anti-HIV drug development efforts. Because of the high genetic diversity of the virus, mutations in RT can impart resistance to various RT inhibitors. As the prevalence of drug resistance mutations is on the rise, it is necessary to design strategies that will lead to drugs less susceptible to resistance. Here we provide an in-depth review of HIV reverse transcriptase, current RT inhibitors, novel RT inhibitors, and mechanisms of drug resistance. We also present novel strategies that can be useful to overcome RT's ability to escape therapies through drug resistance. While resistance may not be completely avoidable, designing drugs based on the strategies and principles discussed in this review could decrease the prevalence of drug resistance.

Translesion synthesis by AMV, HIV, and MMLVreverse transcriptases using RNA templates containing inosine, guanosine, and their 8-oxo-7,8-dihydropurine derivatives

Inosine is ubiquitous and essential in many biological processes, including RNA-editing. In addition, oxidative stress on RNA has been a topic of increasing interest due, in part, to its potential role in the development/progression of disease. In this work we probed the ability of three reverse transcriptases (RTs) to catalyze the synthesis of cDNA in the presence of RNA templates containing inosine (I), 8-oxo-7,8-dihydroinosine (8oxo-I), guanosine (G), or 8-oxo-7,8-dihydroguanosine (8-oxoG), and explored the impact that these purine derivatives have as a function of position. To this end, we used 29-mers of RNA (as template) containing the modifications at position-18 and reverse transcribed DNA using 17-mers, 18-mers, or 19-mers (as primers). Generally reactivity of the viral RTs, AMV / HIV / MMLV, towards cDNA synthesis was similar for templates containing G or I as well as for those with 8-oxoG or 8-oxoI. Notable differences are: 1) the use of 18-mers of DNA (to explore cDNA synthesis past the lesion/modification) led to inhibition of DNA elongation in cases where a G:dA wobble pair was present, while the presence of I, 8-oxoI, or 8-oxoG led to full synthesis of the corresponding cDNA, with the latter two displaying a more efficient process; 2) HIV RT is more sensitive to modified base pairs in the vicinity of cDNA synthesis; and 3) the presence of a modification two positions away from transcription initiation has an adverse impact on the overall process. Steady-state kinetics were established using AMV RT to determine substrate specificities towards canonical dNTPs (N = G, C, T, A). Overall we found evidence that RNA templates containing inosine are likely to incorporate dC > dT > > dA, where reactivity in the presence of dA was found to be pH dependent (process abolished at pH 7.3); and that the absence of the C2-exocyclic amine, as displayed with templates containing 8-oxoI, leads to increased selectivity towards incorporation of dA over dC. The data will be useful in assessing the impact that the presence of inosine and/or oxidatively generated lesions have on viral processes and adds to previous reports where I codes exclusively like G. Similar results were obtained upon comparison of AMV and MMLV RTs.

Reverse Transcriptase: From Transcriptomics to Genome Editing

Reverse transcriptases (RTs) are enzymes that can generate a complementary strand of DNA (cDNA) from RNA. Coupled with PCR, RTs have been widely used to detect RNAs and to clone expressed genes. Classical retroviral RTs have been improved by protein engineering. These enzymes and newly characterized RTs are key elements in the development of next-generation sequencing techniques that are now being applied to the study of transcriptomics. In addition, engineered RTs fused to a CRISPR/Cas9 nickase have recently shown great potential as tools to manipulate eukaryotic genomes. In this review, we discuss the properties and uses of wild type and engineered RTs in biotechnological applications, from conventional RT-PCR to recently introduced prime editing.

Decreased Km to dNTPs is an essential M-MuLV reverse transcriptase adoption required to perform efficient cDNA synthesis in One-Step RT-PCR assay

Personalized medicine and advanced diagnostic tools based on RNA analysis are focusing on fast and direct One-Step RT-PCR assays. First strand complementary DNA (cDNA) synthesized by the reverse transcriptase (RT) is exponentially amplified in the end-point or real-time PCR. Even a minor discrepancy in PCR conditions would result in big deviations during the data analysis. Thus, One-Step RT-PCR composition is typically based on the PCR buffer. In this study, we have used compartmentalized ribosome display technique for in vitro evolution of the Moloney Murine Leukemia Virus reverse transcriptase (M-MuLV RT) that would be able to perform efficient full-length cDNA synthesis in PCR buffer optimized for Thermus aquaticus DNA polymerase. The most frequent mutations found in a selected library were analyzed. Aside from the mutations, which switch off RNase H activity of RT and are beneficial for the full-length cDNA synthesis, we have identified several mutations in the active center of the enzyme (Q221R and V223A/M), which result in 4-5-fold decrease of Km for dNTPs (<0.2 mM). The selected mutations are in surprising agreement with the natural evolution process because they transformed the active center from the oncoretroviral M-MuLV RT-type to the lenitiviral enzyme-type. We believe that this was the major and essential phenotypic adjustment required to perform fast and efficient cDNA synthesis in PCR buffer at 0.2-mM concentration of each dNTP.

Capsid-CPSF6 Interaction Licenses Nuclear HIV-1 Trafficking to Sites of Viral DNA Integration

HIV-1 integration into the host genome favors actively transcribed genes. Prior work indicated that the nuclear periphery provides the architectural basis for integration site selection, with viral capsid-binding host cofactor CPSF6 and viral integrase-binding cofactor LEDGF/p75 contributing to selection of individual sites. Here, by investigating the early phase of infection, we determine that HIV-1 traffics throughout the nucleus for integration. CPSF6-capsid interactions allow the virus to bypass peripheral heterochromatin and penetrate the nuclear structure for integration. Loss of interaction with CPSF6 dramatically alters virus localization toward the nuclear periphery and integration into transcriptionally repressed lamina-associated heterochromatin, while loss of LEDGF/p75 does not significantly affect intranuclear HIV-1 localization. Thus, CPSF6 serves as a master regulator of HIV-1 intranuclear localization by trafficking viral preintegration complexes away from heterochromatin at the periphery toward gene-dense chromosomal regions within the nuclear interior.

Transcriptional inaccuracy threshold attenuates differences in RNA-dependent DNA synthesis fidelity between retroviral reverse transcriptases

In M13mp2 lacZα forward mutation assays measuring intrinsic fidelity of DNA-dependent DNA synthesis, wild-type human immunodeficiency virus type 1 (HIV-1) RTs of group M/subtype B previously showed >10-fold higher error rates than murine leukaemia virus (MLV) and avian myeloblastosis virus (AMV) RTs. An adapted version of the assay was used to obtain error rates of RNA-dependent DNA synthesis for several RTs, including wild-type HIV-1_BH10, HIV-1_ESP49, AMV and MLV RTs, and the high-fidelity mutants of HIV-1_ESP49 RT K65R and K65R/V75I. Our results showed that there were less than two-fold differences in fidelity between the studied RTs with error rates ranging within 2.5 × 10^-5 and 3.5 × 10^-5. These results were consistent with the existence of a transcriptional inaccuracy threshold, generated by the RNA polymerase while synthesizing the RNA template used in the assay. A modest but consistent reduction of the inaccuracy threshold was achieved by lowering the pH and Mg²⁺ concentration of the transcription reaction. Despite assay limitations, we conclude that HIV-1_BH10 and HIV-1_ESP49 RTs are less accurate when copying DNA templates than RNA templates. Analysis of the RNA-dependent mutational spectra revealed a higher tendency to introduce large deletions at the initiation of reverse transcription by all HIV-1 RTs except the double-mutant K65R/V75I.

Diaryl ethers with carboxymethoxyphenacyl motif as potent HIV-1 reverse transcriptase inhibitors with improved solubility

In search of new non-nucleoside reverse transcriptase inhibitors (NNRTIs) with improved solubility, two series of novel diaryl ethers with phenacyl moiety were designed and evaluated for their HIV-1 reverse transcriptase inhibition potentials. All compounds exhibited good to excellent results with IC₅₀ at low micromolar to submicromolar concentrations. Two most active compounds (7e and 7 g) exhibit inhibitory potency comparable or even better than that of nevirapine and rilpivirine. Furthermore, SupT1 and CD4⁺ cell infectivity assays for the most promising (7e) have confirmed its strong antiviral potential while docking studies indicate a novel binding interactions responsible for high activity.

Mechanisms and Factors That Drive Extensive Human Immunodeficiency Virus Type-1 Hypervariability: An Overview

The extensive hypervariability of human immunodeficiency virus type-1 (HIV-1) populations represents a major barrier against the success of currently available antiretroviral therapy. Moreover, it is still the most important obstacle that faces the development of an effective preventive vaccine against this infectious virus. Indeed, several factors can drive such hypervariability within and between HIV-1 patients. These factors include: first, the very low fidelity nature of HIV-1 reverse transcriptase; second, the extremely high HIV-1 replication rate; and third, the high genomic recombination rate that the virus has. All these factors together with the APOBEC3 proteins family and the immune and antiviral drugs pressures drive the extensive hypervariability of HIV-1 populations. Studying these factors and the mechanisms that drive such hypervariability will provide valuable insights that may guide the development of effective therapeutic and preventive strategies against HIV-1 infection in the near future. To this end, in this review, we summarized recent advances in this area of HIV-1 research.

What do docking and QSAR tell us about the design of HIV-1 reverse transcriptase nonnucleoside inhibitors?

Despite vigorous studies, effective nonnucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) are still in demand, not only due to toxicity and detrimental side effects of currently used drugs but also because of the emergence of multidrug-resistant viral strains. In this contribution, we present results of docking of 47 inhibitors to 107 allosteric centers of HIV-1 reverse transcriptase. Based on the average binding scores, we have constructed QSAR equations to elucidate directions of further developments in the inhibitor design that come from this structural data.

Next-generation sequencing-based analysis of reverse transcriptase fidelity

In this study, we devised a simple and rapid method to analyze fidelity of reverse transcriptase (RT) using next-generation sequencing (NGS). The method comprises a cDNA synthesis reaction from standard RNA with a primer containing a tag of 14 randomized bases and the RT to be tested, PCR using high-fidelity DNA polymerase, and NGS. By comparing the sequence of each read with the reference sequence, mutations were identified. The mutation can be identified to be due to an error introduced by either cDNA synthesis, PCR, or NGS based on whether the sequence reads with the same tag contain the same mutation or not. The error rates in cDNA synthesis with Moloney murine leukemia virus (MMLV) RT thermostable variant MM4 or the recently developed 16-tuple variant of family B DNA polymerase with RT activity, RTX, from Thermococcus kodakarensis, were 0.75-1.0 × 10^-4 errors/base, while that in the reaction with the wild-type human immunodeficiency virus type 1 (HIV-1) RT was 2.6 × 10^-4 errors/base. Overall, our method could precisely evaluate the fidelity of various RTs with different reaction conditions in a high-throughput manner without the use of expensive optics and troublesome adaptor ligation.

Physiological Mg2+ Conditions Significantly Alter the Inhibition of HIV-1 and HIV-2 Reverse Transcriptases by Nucleoside and Non-Nucleoside Inhibitors in Vitro

Reverse transcriptases (RTs) are typically assayed in vitro with 5-10 mM Mg²⁺, whereas the free Mg²⁺ concentration in cells is much lower. Artificially high Mg²⁺ concentrations used in vitro can misrepresent different properties of human immunodeficiency virus (HIV) RT, including fidelity, catalysis, pausing, and RNase H activity. Here, we analyzed nucleoside (NRTIs) and non-nucleoside RT inhibitors (NNRTIs) in primer extension assays at different concentrations of free Mg²⁺. At low concentrations of Mg²⁺, NRTIs and dideoxynucleotides (AZTTP, ddCTP, ddGTP, and 3TCTP) inhibited HIV-1 and HIV-2 RT synthesis less efficiently than they did with large amounts of Mg²⁺, whereas inhibition by the "translocation-defective RT inhibitor" EFdA (4'-ethynyl-2-fluoro-2'-deoxyadenosine) was unaffected by Mg²⁺ concentrations. Steady-state kinetic analyses revealed that the reduced level of inhibition at low Mg²⁺ concentrations resulted from a 3-9-fold (depending on the particular nucleotide and inhibitor) less efficient incorporation (based on k_cat/K_m) of these NRTIs under this condition compared to incorporation of natural dNTPs. In contrast, EFdATP was incorporated with an efficiency similar to that of its analogue dATP at low Mg²⁺ concentrations. Unlike NRTIs, NNRTIs (nevirapine, efavirenz, and rilviripine), were approximately 4-fold (based on IC₅₀ values) more effective at low than at high Mg²⁺ concentrations. Drug-resistant HIV-1 RT mutants also displayed the Mg²⁺-dependent difference in susceptibility to NRTIs and NNRTIs. In summary, analyzing the efficiency of inhibitors under more physiologically relevant low-Mg²⁺ conditions yielded results dramatically different from those from measurements using commonly employed high-Mg²⁺ in vitro conditions. These results also emphasize differences in Mg²⁺ sensitivity between the translocation inhibitor EFdATP and other NRTIs.

Ultrasensitive single-genome sequencing: accurate, targeted, next generation sequencing of HIV-1 RNA

BACKGROUND

Although next generation sequencing (NGS) offers the potential for studying virus populations in unprecedented depth, PCR error, amplification bias and recombination during library construction have limited its use to population sequencing and measurements of unlinked allele frequencies. Here we report a method, termed ultrasensitive Single-Genome Sequencing (uSGS), for NGS library construction and analysis that eliminates PCR errors and recombinants, and generates single-genome sequences of the same quality as the "gold-standard" of HIV-1 single-genome sequencing assay but with more than 100-fold greater depth.

RESULTS

Primer ID tagged cDNA was synthesized from mixtures of cloned BH10 wild-type and mutant HIV-1 transcripts containing ten drug resistance mutations. First, the resultant cDNA was divided and NGS libraries were generated in parallel using two methods: uSGS and a method applying long PCR primers to attach the NGS adaptors (LP-PCR-1). Second, cDNA was divided and NGS libraries were generated in parallel comparing 3 methods: uSGS and 2 methods adapted from more recent reports using variations of the long PCR primers to attach the adaptors (LP-PCR-2 and LP-PCR-3). Consistently, the uSGS method amplified a greater proportion of cDNAs, averaging 30% compared to 13% for LP-PCR-1, 21% for LP-PCR-2 and 14% for LP-PCR-3. Most importantly, when the uSGS sequences were binned according to their primer IDs, 94% of the bins did not contain PCR recombinant sequences versus only 55, 75 and 65% for LP-PCR-1, 2 and 3, respectively. Finally, when uSGS was applied to plasma samples from HIV-1 infected donors, both frequent and rare variants were detected in each sample and neighbor-joining trees revealed clusters of genomes driven by the linkage of these mutations, showing the lack of PCR recombinants in the datasets.

CONCLUSIONS

The uSGS assay can be used for accurate detection of rare variants and for identifying linkage of rare alleles associated with HIV-1 drug resistance. In addition, the method allows accurate in-depth analyses of the complex genetic relationships of viral populations in vivo.

Resolution of Specific Nucleotide Mismatches by Wild-Type and AZT-Resistant Reverse Transcriptases during HIV-1 Replication

A key contributor to HIV-1 genetic variation is reverse transcriptase errors. Some mutations result because reverse transcriptase (RT) lacks 3' to 5' proofreading exonuclease and can extend mismatches. However, RT also excises terminal nucleotides to a limited extent, and this activity contributes to AZT resistance. Because HIV-1 mismatch resolution has been studied in vitro but only indirectly during replication, we developed a novel system to study mismatched base pair resolution during HIV-1 replication in cultured cells using vectors that force template switching at defined locations. These vectors generated mismatched reverse transcription intermediates, with proviral products diagnostic of mismatch resolution mechanisms. Outcomes for wild-type (WT) RT and an AZT-resistant (AZT(R)) RT containing a thymidine analog mutation set-D67N, K70R, D215F, and K219Q-were compared. AZT(R) RT did not excise terminal nucleotides more frequently than WT, and for the majority of tested mismatches, both WT and AZT(R) RTs extended mismatches in more than 90% of proviruses. However, striking enzyme-specific differences were observed for one mispair, with WT RT preferentially resolving dC-rC pairs either by excising the mismatched base or switching templates prematurely, while AZT(R) RT primarily misaligned the primer strand, causing deletions via dislocation mutagenesis. Overall, the results confirmed HIV-1 RT's high capacity for mismatch extension during virus replication and revealed dramatic differences in aberrant intermediate resolution repertoires between WT and AZT(R) RTs on one mismatched replication intermediate. Correlating mismatch extension frequencies observed here with reported viral mutation rates suggests a complex interplay of nucleotide discrimination and mismatch extension drives HIV-1 mutagenesis.

Mismatched Primer Extension Assays

Steady state kinetic assays have been a reliable way to estimate fidelity of several polymerases (Menendez-Arias, 2009; Rezende and Prasad, 2004; Svarovskaia et al., 2003). The ability to analyze the extension of primers with specific mismatches at the 3' end is a major strength of the mismatched primer extension assays. Recently, we used the mismatched primer extension assays to show that the fidelity of HIV RT increases dramatically when concentration of Mg²⁺ is reduced to a physiologically relevant concentration (~0.25 mM) (Achuthan et al., 2014). Here, we describe in detail how to perform the mismatched primer extension assay to measure the standard extension efficiency using human immunodeficiency virus reverse transcriptase (HIV RT) at 2 mM Mg²⁺. The relative fidelity of the polymerase can then be estimated using the standard extension efficiency. The assay described here is based on the method published in Mendelman et al. (1990).

Primer Extension Reactions for the PCR- based α- complementation Assay

The PCR- based- α- complementation assay is an effective technique to measure the fidelity of polymerases, especially RNA-dependent RNA polymerases (RDRP) and Reverse Transcriptases (RT). It has been successfully employed to determine the fidelity of the poliovirus polymerase 3D-pol (DeStefano, 2010) as well as the human immunodeficiency virus Reverse Transcriptase (HIV RT) (Achuthan et al., 2014). A major advantage of the assay is that since the PCR step is involved, even the low yield of products obtained after two rounds of low yield of RNA synthesis (for RDRP) or reverse transcription (for RT) can be measured using the assay. The assay also mimics the reverse transcription process, since both RNA- and DNA- directed RT synthesis steps are performed. We recently used this assay to show that the HIV RT, at physiologically relevant magnesium concentration, has accuracy in the same range as other reverse transcriptases (Achuthan et al., 2014). Here, we describe in detail how to prepare the inserts using the primer extension reactions. The prepared inserts are then processed further in the PCR- based- α- complementation assay.

Alternative divalent cations (Zn²⁺, Co²⁺, and Mn²⁺) are not mutagenic at conditions optimal for HIV-1 reverse transcriptase activity

BACKGROUND

Fidelity of DNA polymerases can be influenced by cation co-factors. Physiologically, Mg(2+) is used as a co-factor by HIV reverse transcriptase (RT) to perform catalysis; however, alternative cations including Mn(2+), Co(2+), and Zn(2+) can also support catalysis. Although Zn(2+) supports DNA synthesis, it inhibits HIV RT by significantly modifying RT catalysis. Zn(2+) is currently being investigated as a component of novel treatment options against HIV and we wanted to investigate the fidelity of RT with Zn(2+).

METHODS

We used PCR-based and plasmid-based alpha complementation assays as well as steady-state misinsertion and misincorporation assays to examine the fidelity of RT with Mn(2+), Co(2+), and Zn(2+).

RESULTS

The fidelity of DNA synthesis by HIV-1 RT was approximately 2.5 fold greater in Zn(2+) when compared to Mg(2+) at cation conditions optimized for nucleotide catalysis. Consistent with this, RT extended primers with mismatched 3' nucleotides poorly and inserted incorrect nucleotides less efficiently using Zn(2+) than Mg(2+). In agreement with previous literature, we observed that Mn(2+) and Co(2+) dramatically decreased the fidelity of RT at highly elevated concentrations (6 mM). However, surprisingly, the fidelity of HIV RT with Mn(2+) and Co(2+) remained similar to Mg(2+) at lower concentrations that are optimal for catalysis.

CONCLUSION

This study shows that Zn(2+), at optimal extension conditions, increases the fidelity of HIV-1 RT and challenges the notion that alternative cations capable of supporting polymerase catalysis are inherently mutagenic.