Quantification of Deaminase Activity-Dependent and -Independent Restriction of HIV-1 Replication Mediated by APOBEC3F and APOBEC3G through Experimental-Mathematical Investigation
- Tomoko Kobayashia,
- Yoshiki Koizumib,
- Junko S. Takeuchia,
- Naoko Misawaa,
- Yuichi Kimuraa,
- Satoru Moritac,
- Kazuyuki Aiharad,e,
- Yoshio Koyanagia,
- Shingo Iwamif and
- Kei Satoa
- aLaboratory of Viral Pathogenesis, Institute for Virus Research, Kyoto University, Kyoto, Japan
- bSchool of Medicine, College of Medical, Pharmaceutical and Health Sciences, Kanazawa University, Ishikawa, Japan
- cDepartment of Mathematical and Systems Engineering, Shizuoka University, Shizuoka, Japan
- dInstitute of Industrial Science, The University of Tokyo, Tokyo, Japan
- eGraduate School of Information Science and Technology, The University of Tokyo, Tokyo, Japan
- fDepartment of Biology, Faculty of Sciences, Kyushu University, Fukuoka, Japan
- S. R. Ross, Editor
ABSTRACT
APOBEC3F and APOBEC3G cytidine deaminases potently inhibit human immunodeficiency virus type 1 (HIV-1) replication by enzymatically inserting G-to-A mutations in viral DNA and/or impairing viral reverse transcription independently of their deaminase activity. Through experimental and mathematical investigation, here we quantitatively demonstrate that 99.3% of the antiviral effect of APOBEC3G is dependent on its deaminase activity, whereas 30.2% of the antiviral effect of APOBEC3F is attributed to deaminase-independent ability. This is the first report quantitatively elucidating how APOBEC3F and APOBEC3G differ in their anti-HIV-1 modes.
TEXT
Human immunodeficiency virus type 1 (HIV-1), the causative agent of AIDS, hijacks lines of cellular proteins for its replication (reviewed in reference 1). On the other hand, accumulating evidence has revealed that human cells intrinsically possess anti-HIV-1 proteins (reviewed in reference 2). Human apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) proteins are cellular cytidine deaminases, and certain APOBEC3 proteins, particularly APOBEC3F (3, 4) and APOBEC3G (5), are the intrinsic restriction factors against HIV-1. APOBEC3F and APOBEC3G are incorporated into assembling HIV-1 particles and brought into the newly infected cells. During reverse transcription (RT) in the infected cells, these APOBEC3 proteins enzymatically convert C in the viral minus-strand DNA to U, resulting in G-to-A mutation in the nascent plus-strand DNA. These mutations can become nonsynonymous or even lethal, which severely debilitates subsequent viral replication.
To antagonize these APOBEC3F- and APOBEC3G-mediated antiviral effect, viral infectivity factor (Vif), an accessory protein of HIV-1, degrades APOBEC3F and APOBEC3G through a ubiquitin/proteasome-dependent pathway (reviewed in references 6 and 7). Previous studies have directly demonstrated that Vif is a prerequisite for HIV-1 replication in vitro (e.g., human primary CD4+ T lymphocytes and macrophages) (8, 9) and in vivo (e.g., humanized mouse models) (10, 11). Moreover, G-to-A mutations have been clearly observed in the proviral DNA of HIV-1-infected patients, although the frequencies of G-to-A mutations seem to vary among individuals (12–23).
It has been reported that the incorporated APOBEC3F and APOBEC3G potently impair the HIV-1 RT process even prior to the insertion of G-to-A mutations (24–26). Moreover, this inhibition is independent of these APOBEC3s' deaminase activity (26). These findings indicate that APOBEC3F and APOBEC3G potently inhibit HIV-1 replication through at least two distinct modes: (i) deaminase activity-dependent G-to-A mutations of viral DNA and (ii) impairment of viral RT independent of deaminase activity (Fig. 1). However, how much of the inhibition of HIV-1 replication by APOBEC3 proteins is attributed to each of these two modes has not been quantitatively revealed.
Schematic diagram of the HIV-1 life cycle and the antiviral effect of APOBEC3 proteins. (Left) Incoming WT APOBEC3 proteins (purple) cause (i) the inhibition of the viral RT process independent of their deaminase activity (blue arrow) and (ii) deaminase activity-dependent C-to-U mutations in minus-strand viral DNA (vDNA) resulting in G-to-A mutations in plus-strand vDNA (red arrow), both of which lead to the abrogation of viral replication; (Right) on the other hand, incoming CI APOBEC3 proteins (blue) cause only the inhibition of the viral RT process independent of their deaminase activity (blue arrow).
In order to quantitatively determine the relative extent of deaminase activity-dependent and -independent antiviral activity of APOBEC3 proteins, we performed a single-round infection assay using a cell culture system. 293T cells and TZM-bl cells (obtained through the NIH AIDS Research and Reference Reagent Program) (27) were maintained in Dulbecco's modified Eagle medium (Sigma) containing 10% fetal calf serum (FCS) and antibiotics. Flag-tagged wild-type (WT) APOBEC3F (GenBank accession number NM_145298.5) and WT APOBEC3G (GenBank accession number NM_021822.3) expression plasmids are based on pcDNA3.1 (Life Technologies).
Regarding the deaminase activity of APOBEC3 proteins, it has been reported that the glutamate (E) in the catalytic site of APOBEC3 protein is involved in the proton shuttling during deamination reaction (28). Also, it has been already revealed that APOBEC3F E251 and APOBEC3G E259 are the catalytic glutamates based on mutagenesis experiments (APOBEC3F E251Q and APOBEC3G E259Q, respectively) (26, 29) and their crystal structures (30, 31). Upon these findings, we constructed catalytically inactive (CI) mutants, Flag-tagged APOBEC3F E251Q and APOBEC3G E259Q expression plasmids. These plasmids were constructed by using the GENEART site-directed mutagenesis system (Life Technologies), with WT plasmids as the template and the following primers: APOBEC3F E251Q forward, 5′-ACC CAT TGT CAT GCA CAG AGG TGC TTC CTC TCT-3′; APOBEC3F E251Q reverse, 5′-AGA GAG GAA GCA CCT CTG TGC ATG ACA ATG GGT-3′; APOBEC3G E259Q forward, 5′-GAA GGC CGC CAT GCA CAG CTG TGC TTC CTG GAC-3′; APOBEC3G E259Q reverse, 5′-GTC CAG GAA GCA CAG CTG TGC ATG GCG GCC TTC-3′. The sequences of the constructed plasmids were confirmed by using an ABI Prism 3130xl genetic analyzer (Applied Biosystems), and the data were analyzed by Sequencher v5.1 software (Hitachi).
To investigate the antiviral effect of WT and CI APOBEC3 proteins, 1 μg of a vif-deficient HIV-1-producing plasmid (pNL4-3Δvif; based on HIV-1 strain NL4-3 [GenBank accession number M19921.2]) (32) was cotransfected with the Flag-tagged plasmid expressing either WT APOBEC3F, CI APOBEC3F, WT APOBEC3G, or CI APOBEC3G at 5 different doses (0, 11, 33, 100, or 300 ng) into 293T cells by using Lipofectamine 2000 (Life Technologies) according to the manufacturer's protocol. At 48 h posttransfection, the culture supernatant was harvested, centrifuged, and then filtered through a 0.45-μm-pore-size filter (Millipore) to produce virus solution. To detect a viral protein (p24CA) and virion-incorporated Flag-tagged APOBEC3 proteins, SDS-PAGE/Western blotting was performed as previously described (11, 33–36). Briefly, 1 ml of the virus solution was ultracentrifuged at 100,000 × g for 1 h at 4°C using a TL-100 instrument (Beckman). The pellet was lysed with 1× SDS buffer, and the levels of viral protein (p24CA) and virion-incorporated Flag-tagged APOBEC3 proteins were detected by SDS-PAGE/Western blotting. For the Western blotting, anti-p24CA (ViroStat) and anti-Flag (clone M2; Sigma) antibodies were used. To measure the infectivity of virus solution, the TZM-bl assay was performed as previously described (11, 33–36). Briefly, 100 μl of the virus solution was inoculated into TZM-bl cells, and the β-galactosidase activity was measured by using the Galacto-Star mammalian reporter gene assay system (Roche) and a 1420 ALBOSX multilabel counter instrument (PerkinElmer) according to the manufacturer's protocol. To validate that the CI mutants used in this study (APOBEC3F E251Q and APOBEC3G E259Q) were defective in their deaminase activities, 100 μl of the virus solutions was treated with DNase I (TaKaRa) at 37°C for 2 h for the removal of the transfected plasmids and then was inoculated into TZM-bl cells. DNA was extracted by using the DNeasy blood and tissue kit (Qiagen), and the viral DNA was amplified by PCR using the following primers: forward, 5′-GTT TGG AAA GGA CCA GCA AA-3′; reverse, 5′-GCC CAA GTA TCC CCG TAA GT-3′. PCR was performed by using PrimeSTAR GXL DNA polymerase (TaKaRa) according to the manufacturer's protocol in the following sequence of conditions: (i) 98°C for 2 min; (ii) 98°C for 10 s; (iii) 55°C for 15 s; (iv) 68°C for 45 s, repeat steps ii to iv for 30 cycles; (v) 68°C for 5 min. The resultant 774 bp of the viral DNA fragment was cloned into pCR-Blunt II-TOPO vector (Life Technologies) by using the Zero Blunt TOPO PCR cloning kit (Life Technologies) according to the manufacturer's protocol. The sequencing PCR was performed by using M13 sequencing primers and an ABI Prism 3130xl genetic analyzer (Applied Biosystems), and the data were analyzed by Sequencher v5.1 software (Hitachi) and Hypermut v2.0 (http://www.hiv.lanl.gov/content/sequence/HYPERMUT/hypermut.html) (37).
As shown in Fig. 2A, APOBEC3 expression did not affect the level of released HIV-1 particles. In addition, APOBEC3 proteins were incorporated into the released virions in dose-dependent manners, and the levels of WT APOBEC3 proteins in the released virions were comparable to those of CI mutants (Fig. 2A). Next, we assessed the infectivity of the harvested viruses by TZM-bl assay (11, 33–36) and found that both WT and CI APOBEC3 proteins attenuated HIV-1 infectivity in dose-dependent manners (Fig. 2B). Consistent with previous reports (38–40), the anti-HIV-1 effect of WT APOBEC3G was greater than that of WT APOBEC3F (Fig. 2B). On the other hand, it was of interest that the antiviral effect of CI APOBEC3G (APOBEC3G E259Q) was smaller than that of CI APOBEC3F (APOBEC3F E251Q) (Fig. 2B). We then confirmed that the CI mutants did not insert G-to-A mutations at all, while WT APOBEC3 proteins elicited G-to-A mutations in dose-dependent manners (see Fig. 4A and B). These results strongly suggest that WT APOBEC3 proteins exert both G-to-A mutation-mediated (deaminase activity-dependent) and RT inhibition-mediated (i.e., deaminase activity-independent) anti-HIV-1 effects, whereas CI mutants exert only the latter effect (Fig. 1).
Antiviral ability of APOBEC3F and APOBEC3G proteins. (A) Western blotting of the released HIV-1 particles. The level of HIV-1 particles released into the culture supernatant was analyzed by Western blotting, and representative results are shown. (B) Quantification of antiviral activity of APOBEC3F and APOBEC3G by experimental and mathematical analyses. The infectivity of released HIV-1 particles was quantified by TZM-bl assay. The assay was performed in triplicate. The averages of experimental data are shown as dots with standard deviations represented by bars, and a graph of the best-fit model is superimposed.
In order to quantitatively estimate the contribution of these two mechanisms of action (i.e., RT inhibition and G-to-A mutation)
in restricting HIV-1 replication, we built a mathematical model from the experimental data. Because RT inhibition and G-to-A
mutation are independently caused by WT APOBEC3 proteins, their effect as a whole is represented as the union of two independent
events. Here, we estimated the combination of each antiviral effect by Bliss independence, which is a major index for the
evaluation of drug combination effects (41, 42). We define x as the amount of transfected APOBEC3 expression plasmid (ng), and fRT(x) and fMu(x) are the unaffected fractions of viral infectivity by RT inhibition and G-to-A mutation, respectively. Then, the total effect
of the two antiviral activities, fWT(x), is assumed to be the product of the individual effect of the two antiviral activities as follows:
(1)
We analyzed the experimental data of viral infectivity at the 4 different doses (Fig. 2A), using the phenomenological model, equation 1. Furthermore, we assume that RT inhibition and G-to-A mutation are dose-dependently described by simple Hill functions of
x, and then we have
(2) and
(3) respectively. Here, hRT and hMu are the Hill coefficients, and KRT and KMu are the amounts of transfected APOBEC3 expression plasmid (ng) required to decrease the viral infectivity by 50%. Using equations 2 and 3, the viral infectivities affected by WT and CI APOBEC3 are calculated by
(4) and
(5) respectively. To estimate the parameters hRT, hMu, KRT, and KMu, we fitted the above models fWT(x) and fCI(x) (i.e., equations 4 and 5) to 9 sets of the corresponding experimental measurement of HIV-1 infectivity by using a nonlinear least-squares method,
which minimizes the sum of squared residuals (SSR) as follows:
(6) where EWT(xi) and ECI(xi) are the experimental data of the infectivity of viruses released from the cells cotransfected with different doses of WT
or CI APOBEC3 expression plasmid xi (i = 1, 2, 3, 4, and 5 correspond to xi = 0, 11, 33, 100, and 300 ng). The principle axis method in Mathematica 8.0 (Wolfram) was used for searching the optimal parameters.
The experiments shown in Fig. 2B were performed in triplicate for each WT and CI dose. By pairing these two sets of three values, we obtained 3 × 3 = 9 combinations
between WT and CI values. For each combination, we estimated h and K values by fitting (i.e., we obtained 9 h values and 9 K values for APOBEC3F and APOBEC3G, respectively). To calculate the confidence intervals for the estimated parameters, we used
a nonparametric bootstrap method with 20,000 resamples of 9 KMu values of APOBEC3F and APOBEC3G, respectively. Using the mean values of the estimated parameters, fMu(x) predicts the anti-HIV-1 effect of APOBEC3F and APOBEC3G solely mediated by G-to-A mutations (Fig. 2B, represented in red lines). Intriguingly, these analyses revealed that the G-to-A mutation-mediated anti-HIV-1 effect of
APOBEC3G (hMu of APOBEC3G, 1.93 ± 0.18; KMu of APOBEC3G, 5.38 ± 0.49; 95% confidence interval, 5.09 to 5.69) was significantly greater than that of APOBEC3F (hMu of APOBEC3F, 1.41 ± 0.39; KMu of APOBEC3F, 40.8 ± 21.3; 95% confidence interval, 28.96 to 54.95) (Fig. 2B; KMus between APOBEC3F and APOBEC3G, P = 0.038 by the nonparametric bootstrap method). In addition, sequencing analysis revealed that the hypermutation index (38, 43), which represents the level of G-to-A mutation in viral DNA, of APOBEC3G was higher than that of APOBEC3F (see Fig. 4B). Because the levels of APOBEC3F and APOBEC3G in the released virions were comparable (Fig. 2A), these results suggest that APOBEC3G more frequently inserts G-to-A mutations and, therefore, more potently attenuates viral
infection than APOBEC3F does.
Furthermore, to quantitatively estimate the relative contribution of G-to-A mutation and RT inhibition to the anti-HIV-1 effect,
we defined the relative contribution rates of G-to-A mutation and RT inhibition to the total antiviral effects, RCMu(x) and RCRT(x), as follows:
(7)
(8) where fRT(x) − fWT(x) means the amplification of antiviral effects by adding the G-to-A mutation to the condition of preexisting RT inhibition,
and fRT(x) − fWT(x) + fMu(x) − fWT(x) is the sum of each amplification of the antiviral effects by G-to-A mutation and RT inhibition (Fig. 3A). In other words, RCMu(x) represents the contribution of the G-to-A mutation to the total antiviral effects relative to that of the RT inhibition.
Calculating RCMu(x) and RCRT(x) in WT APOBEC3F and WT APOBEC3G with equations 7 and 8 shows that RCRT(x) in WT APOBEC3F (30.2% ± 4.18%; 95% confidence interval, 17.07% to 45.38%) was greater than that in WT APOBEC3G (0.66% ±
0.69%; 95% confidence interval, 0.42% to 0.93%) (Fig. 3B; P = 0.043 by the nonparametric bootstrap method).
Quantification of mutation-dependent and -independent antiviral effects by APOBEC3F and APOBEC3G. (A) A Venn diagram of the relative contribution rates of deaminase activity-dependent G-to-A mutations [Mutation, RCMu(x)] and deaminase activity-independent RT inhibition [RT inhibition, RCRT(x)] to the total antiviral effect of APOBEC3 proteins. (B) Relative contribution rate of the different modes of antiviral activities with APOBEC3F and APOBEC3G. RCMu(x) (red area) and RCRT(x) (blue area) of APOBEC3F and APOBEC3G are estimated as described in the text. RCRT(x) in APOBEC3F, 11 ng, 26.5% ± 36.7%; 33 ng, 27.9% ± 30.6%; 100 ng, 30.6% ± 20.2%; 300 ng, 36.0% ± 12.9%. RCRT(x) in APOBEC3G, 11 ng, 1.64% ± 1.02%; 33 ng, 0.64% ± 0.41%; 100 ng, 0.26% ± 0.17%; 300 ng, 0.11% ± 0.08%.
Summarizing these mathematical analyses (Fig. 3B), the relative contribution of the deaminase activity-dependent antiviral effect of APOBEC3G [RCMu(x) = 99.3% ± 0.69%] was 150-fold greater than its deaminase-independent (i.e., RT inhibition-dependent) antiviral effect [RCRT(x) = 0.66% ± 0.69%]. This indicates that APOBEC3G strongly depends on its catalytic activity in inhibiting HIV-1 replication.
On the other hand, the relative contribution of deaminase activity to the antiviral effect of APOBEC3F [RCMu(x) = 69.8% ± 4.18%] was only 2.3-fold greater than that of its RT inhibition [RCRT(x) = 30.2% ± 4.18%]. These findings suggest that APOBEC3G depends more on its deaminase-dependent antiviral effect than APOBEC3F
does. In this regard, it is known that APOBEC3F induces mainly GA-to-AA mutations, while APOBEC3G induces more of GG-to-AG mutations (italics and underlining indicate the mutation) (7, 44). Consistent with previous observations (7, 44), we observed that WT APOBEC3F and WT APOBEC3G inserted GA-to-AA and GG-to-AG mutations, respectively, in dose-dependent manners (Fig. 4B). For a detailed assessment of the impact of G-to-A mutations caused by APOBEC3F and APOBEC3G, we performed in silico analyses. If there are L lethal mutation sites (i.e., TGG [tryptophan] to TAG [stop codon]) in a full-genome proviral DNA sequence, which is constructed by F sites in total, we can calculate the proportion S of proviral DNA with at least one stop codon after the proviral DNA suffers from N times of APOBEC3-mediated mutations per full-length proviral DNA as follows:
(9)
Mutation signatures and in silico analysis. (A and B) Summary of viral DNA sequence analyzed in this study. (A) Raw sequence data (the sequences are available upon request). (B) Summarized table. a, The number in parenthesis represents the amount of transfected APOBEC3-expressing plasmid; b, the numbers of analyzed nucleotides and clones (in parenthesis) are shown; c, hypermutation index (38, 43) is calculated as follows: (no. of G-to-A mutations − no. of A-to-G mutations)/sequence length (bp) × 1,000; d, the number and the percentage (in parenthesis) of dinucleotide contexts of G-to-A mutation sites in total G-to-A mutations are shown. ND, not determined, probably due to a high level of G-to-A hypermutations mediated by APOBEC3G in the primer-binding region of template DNA that resulted in the failure of the PCR process under this condition (53, 54). (C) In silico analyses. The effects of G-to-A mutations mediated by APOBEC3F (GA to AA, red line) or APOBEC3G (GG to AG, blue line) on the proportion of proviral DNA with stop codons (S, y axis) were estimated with equation 9 as described in the text.
To perform this in silico simulation, we obtained 1,596 full-genome HIV-1 sequences registered in the Los Alamos HIV sequence database (http://www.hiv.lanl.gov). In the 1,596 data sets (containing 13,799,872 bases), the following motifs were picked up by using the seqinR package for the R statistical computing environment v2.15.2 (http://www.r-project.org): the number of APOBEC3F-targeted motifs (FAPOBEC3F, the number of GA sites), 1,134,377; the number of APOBEC3F-targeted motifs directly leading to lethal mutation (LAPOBEC3F, the number of TGGA with TGG as a codon, leading to TGAA), 59,788; the number of APOBEC3G-targeted motifs (FAPOBEC3G, the number of GG sites), 929,607; the number of APOBEC3G-targeted motifs directly leading to lethal mutation (LAPOBEC3G, the number of TGG with TGG as a codon, leading to TAG), 151,234. Calculating SAPOBEC3F and SAPOBEC3G with equation 9 shows that APOBEC3G-mediated mutations (GG to AG) quickly become lethal (Fig. 4C; 3.90 mutations with 50% being lethal, 13.0 mutations with 90% being lethal). On the other hand, the APOBEC3F-mediated mutation (GA to AA) less frequently became lethal (Fig. 4C, 12.80 mutations with 50% being lethal, 42.53 mutations with 90% being lethal). Taken together, these findings well explain the differences in the profiles of antiviral activities between these APOBEC3 proteins, which were observed in our experimental and mathematical investigations (Fig. 2B and 3). They strongly suggest that APOBEC3G impairs HIV-1 replication mostly dependent on its deaminase activity. In addition, APOBEC3G is more effective than APOBEC3F because of its (i) higher G-to-A mutation frequency than APOBEC3F (indicated by hypermutation index in Fig. 4B) and (ii) mutation signature readily resulting in lethal mutation (GG-to-AG mutations; Fig. 4C). On the other hand, APOBEC3F suppresses HIV-1 replication with the combination of G-to-A mutations and RT inhibition. In our in silico simulation (Fig. 4C), we assumed that only termination mutations (L) result in the abrogation of viral replication. However, there are other mutations which can result in the failure of viral replication, such as (i) long-terminal repeat, (ii) polypurine tracts, (iii) tRNA-binding regions, (iv) initiation codons (i.e., ATG), and (v) nonsynonymous mutations in viral genes. Therefore, as previously described (44), it should be noted that our simulation (Fig. 4C) may underestimate the antiviral effect of G-to-A mutations caused by APOBEC3.
Although we observed the dose-dependent increase of G-to-A mutations by WT APOBEC3F and APOBEC3G (Fig. 4B), we detected some amplicons without G-to-A mutations (Fig. 4A). This raises a possibility that a portion of viruses released from the cells transfected with APOBEC3-expressing plasmid may not package APOBEC3 proteins. On the other hand, it is known that the viral clones suffered from APOBEC3-mediated G-to-A mutations are less efficiently amplified by PCR than the clones without any mutations, because the primer-binding sites can be also mutated (44). Therefore, the nonmutated amplicons detected (Fig. 4A) may be due to a technical bias on PCR rather than the presence of virions without packaged APOBEC3 proteins. Moreover, a biochemical analysis has previously estimated that approximately 4 to 9 APOBEC3 molecules are packaged per vif-deficient HIV-1 particle (45). However, it seems technically difficult to quantify the absolute number of APOBEC3 molecules per virion. To fully reveal this issue, further investigations will be needed.
APOBEC3-associated mutation patterns (i.e., GG-to-AG and/or GA-to-AA mutations) have been observed in HIV-1-infected CD4+ T lymphocytes (38, 39, 46) and macrophages (38) as well as in infected humanized mouse models (10, 11) and individuals (14, 22), strongly suggesting that some APOBEC3 proteins, at least APOBEC3F and APOBEC3G, are closely associated with HIV-1 infection, replication, and pathogenesis. In addition, APOBEC3F and APOBEC3G are highly associated with single-strand nucleic acids, including viral RNA and DNA, which is a prerequisite for their packaging into the virions and exerting their antiviral activity in the newly infected cells (47–49). APOBEC3F and APOBEC3G catalyze C-to-U mutations by directly binding to minus-strand viral DNA, and the binding of APOBEC3 to minus-strand viral DNA results in the inhibition of the viral RT process, which is independent of APOBEC3's catalytic activity (47, 48). Therefore, it is technically unable to prepare the APOBEC3 mutant causing G-to-A mutations without RT inhibition and to experimentally elucidate the sole effect of deaminase activity of APOBEC3 proteins on HIV-1 replication. By applying mathematical modeling to the experimental data, here we quantitatively estimated the sole effect of G-to-A mutations by APOBEC3 proteins, which cannot be determined through conventional experimental methods (Fig. 2B; the sole effects of G-to-A mutations by APOBEC3F and APOBEC3G are represented in red lines). Moreover, since the combination of the two different antiviral mechanisms has some synergistic or additive effects to inhibit viral replication, the combination of deaminase activity-dependent and -independent effects cannot be estimated by linearly summing the experimental measurement of their individual antiviral effect. Similar problems occur in evaluating the nonlinear effects of two-drug combinations, and, therefore, mathematical models such as Bliss independence are useful and necessary for estimating the combination effects by using single-drug effects. In this study, we applied Bliss independence to evaluate the antiviral effects of G-to-A mutation and RT inhibition in consideration of the nonlinear effects of two different antiviral mechanisms of action.
So far, our experimental and mathematical approach has quantitatively revealed the replication dynamics of retroviruses (50, 51) and enteroviruses (52). To the best of our knowledge, here we quantitatively determined the antiviral effect of cellular proteins, human APOBEC3, using the experimental data for the first time. The synergistic experimental and mathematical strategy is a powerful approach to quantitatively investigate the dynamics of virus infections in a way that is impossible by conventional experimental strategies.
ACKNOWLEDGMENTS
We thank Keiko Okano (Research Administration Office, Kyoto University) for proofreading the manuscript.
This study was supported in part by grants from the following: the Aihara Innovative Mathematical Modeling Project, the Japan Society for the Promotion of Science (JSPS), through the Funding Program for World-Leading Innovative R&D on Science and Technology (FIRST Program), initiated by the Council for Science and Technology Policy of Japan (to Kazuyuki Aihara, Shingo Iwami, and Kei Sato); Takeda Science Foundation (to Kei Sato); Sumitomo Foundation Research Grant (to Kei Sato); Senshin Medical Research Foundation (to Kei Sato); Kyushu University Interdisciplinary Programs in Education and Projects in Research Development (to Shingo Iwami); Grant-in-Aid for Young Scientists B25800092 (to Shingo Iwami) from JSPS; Grants-in-Aid for Scientific Research B24390112 and S22220007 (to Yoshio Koyanagi) from JSPS; Grant-in-Aid for Scientific Research on Innovative Areas 24115008 (to Yoshio Koyanagi) from the Ministry of Education, Culture, Sports, Science and Technology of Japan; and Research on HIV/AIDS (to Yoshio Koyanagi) from the Ministry of Health, Labor and Welfare of Japan.
FOOTNOTES
- Received 8 January 2014.
- Accepted 3 March 2014.
- Address correspondence to Kei Sato, ksato{at}virus.kyoto-u.ac.jp, or Shingo Iwami, siwami{at}kyushu-u.org.
-
Tomoko Kobayashi and Yoshiki Koizumi contributed equally to this study.
-
Published ahead of print 12 March 2014
REFERENCES
- Copyright © 2014, American Society for Microbiology. All Rights Reserved.
