Email updates

Keep up to date with the latest news and content from BMC Medical Genomics and BioMed Central.

Open Access Research article

Host sequence motifs shared by HIV predict response to antiretroviral therapy

William Dampier1, Perry Evans2, Lyle Ungar2 and Aydin Tozeren1*

Author affiliations

1 Center for Integrated Bioinformatics, Drexel University, Bossone Research Center 711, 3120 Market Street, Philadelphia, PA 19104, USA

2 Genomics and Computational Biology and Department of Computer and Information Science, University of Pennsylvania, Levine Hall, 3330 Walnut Street, Philadelphia, PA 19104, USA

For all author emails, please log on.

Citation and License

BMC Medical Genomics 2009, 2:47  doi:10.1186/1755-8794-2-47


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1755-8794/2/47


Received:4 December 2008
Accepted:23 July 2009
Published:23 July 2009

© 2009 Dampier et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The HIV viral genome mutates at a high rate and poses a significant long term health risk even in the presence of combination antiretroviral therapy. Current methods for predicting a patient's response to therapy rely on site-directed mutagenesis experiments and in vitro resistance assays. In this bioinformatics study we treat response to antiretroviral therapy as a two-body problem: response to therapy is considered to be a function of both the host and pathogen proteomes. We set out to identify potential responders based on the presence or absence of host protein and DNA motifs on the HIV proteome.

Results

An alignment of thousands of HIV-1 sequences attested to extensive variation in nucleotide sequence but also showed conservation of eukaryotic short linear motifs on the protein coding regions. The reduction in viral load of patients in the Stanford HIV Drug Resistance Database exhibited a bimodal distribution after 24 weeks of antiretroviral therapy, with 2,000 copies/ml cutoff. Similarly, patients allocated into responder/non-responder categories based on consistent viral load reduction during a 24 week period showed clear separation. In both cases of phenotype identification, a set of features composed of short linear motifs in the reverse transcriptase region of HIV sequence accurately predicted a patient's response to therapy. Motifs that overlap resistance sites were highly predictive of responder identification in single drug regimens but these features lost importance in defining responders in multi-drug therapies.

Conclusion

HIV sequence mutates in a way that preferentially preserves peptide sequence motifs that are also found in the human proteome. The presence and absence of such motifs at specific regions of the HIV sequence is highly predictive of response to therapy. Some of these predictive motifs overlap with known HIV-1 resistance sites. These motifs are well established in bioinformatics databases and hence do not require identification via in vitro mutation experiments.

Background

Human Immunodeficiency Virus (HIV) is a single stranded RNA virus that contains nine genes coding for fifteen proteins [1,2]. HIV has a powerful effect on the human immune system due to its ability to hijack hundreds of human proteins in continued infection [3]. HIV's POL gene codes for three important enzymes that are essential to the life cycle of the virus: the protein reverse transcriptase (RT) is common to all retroviruses and transcribes the viral RNA into double stranded DNA [1]. The RT enzyme has no proofreading ability [4] which explains the high mutation rate observed with in vitro experiments for the HIV virus [5]. POL also encodes the integrase protein which fuses the viral DNA produced by RT into the host genome [4]. The third enzyme coded by POL, protease (PR), is an enzyme that cleaves the multiple proteins coded by HIV's GAG and POL genes into separate functional units [1]. Mutations at the active sites of these three enzymes or inhibition of enzyme activity by drugs disrupt HIV's ability to replicate in host cells and thus block the infection cycle [6].

Most of the drugs that are currently used for controlling HIV infection target the three viral enzymes coded by the HIV POL gene. Antiretroviral drugs such as zidovudine (AZT), lamivudine (3TC), emtricitabine (FTC), zalcitabine (ddC), stavudine (D4T), didanosine (DDI) and nevirapine (NVP) target RT [7] whereas antiretroviral drugs such as indinavir (IDV), nelfinavir (NFV), and atazanavir (ATV) were designed as PR inhibitors [8]. Clinicians also use a set of entry and integrase inhibitors in HIV treatment [9-11]. When antiretroviral drug are used one at a time, eventually a drug resistant viral phenotype will emerge [12]. Viral loads (VL) from in vitro cultures of HIV infected immune cells have diminishing growth rates in the presence of antiretroviral therapy but eventually a resistant viral phenotype emerges [13]. The resistance conferring mutations in the viral genome have been extensively documented and these mutations have been correlated to response to therapy [13-16]. Combination of antiretroviral drugs has the advantage of targeting multiple stages of the viral life cycle. The multi-target Highly Active Antiretroviral therapies (HAART) exert a high level of evolutionary pressure on the virus by effectively requiring multiple simultaneous mutations to produce resistant strains [17-19]. As a result, the virus takes much longer time to develop resistance to several drugs at the same time [20].

HAART therapies often reduce viral replication to undetectable levels. They decrease morbidity and mortality rates but nonetheless can be ineffective in some individuals [21,22]. Search for new antiretroviral drugs with different target sites along the HIV sequence is ongoing. Targeting the virus itself may not be enough, however, to block the progress of infection. One may also have to consider the set of host proteins playing crucial roles in viral replication as targets for therapy. Recently, researchers have identified sets of human proteins that interact with HIV proteins [23-25] and another set of host proteins required for HIV infection through a functional genomic screen (([26-28], but the modes of interaction of these host proteins with specific HIV proteins are yet to be fully explored. Nevertheless, the ability of HIV-1 viral proteins to bind within the host cell network is likely to play a critical role in disease progression [29]. It is possible that this new focus on host proteins interacting with HIV will lead to new therapies targeting host cells required for HIV infection [30].

In this study, we first cluster patients into responder and non-responder categories based on viral load response to antiretroviral therapy. We then used stepwise logistic regression to differentiate responders and non-responders using linear sequence motifs common to host and viral genomes as features. We focused on viral load in the responder/non-responder classification because recent studies indicate that CD4 cell count monitoring does not accurately identify individuals with virologic failure among patients taking antiviral therapy [31]. A novel aspect of our study is the recognition of bimodality [32] in the viral load reduction in antiretroviral therapy in patient data stored in the Stanford HIV Drug Resistance Database [33] both at eight weeks and twenty four weeks after the beginning of the therapy. In total, we used three different methods for assigning responder phenotype based on viral load. Multiple models of phenotype classification allowed us to identify the role of phenotype selection in determining significant features associated with drug response.

Another novel feature of our study is the treatment of drug response as a two body problem, namely that response to drugs is assumed to be affected by both the viral and host genotypes. We sought to identify linear motifs on the HIV sequence that are also found in the host and are functionally annotated: host transcription factor binding sequence motifs [34], miRNA binding sequence motifs on the nucleotide sequence [35] and eukaryotic linear motifs [36] on the protein amino acid sequences. The motivation to use such features in predicting responder/non-responder categories comes from the observed phenomena of the virus hijacking host cell apparatus for its self replication [37]. Another important motivation is to find a feature set based solely on viral sequence and not requiring a priori information obtained via virus-specific in vitro cell assays. This type of a feature set is attractive, as it can be used to explore the drug response of viruses to antiviral therapy in the absence of extensive data on resistance mutations. Previous research on quantitative prediction of patient response to antiretroviral drugs in HIV infection [38-43] has employed similar and even more advanced machine learning algorithms than used here, but has not made explicit use of biologically meaningful linear motifs.

Results

Responder/Non-responder classification

Clinical annotation of more than 2,000 RT sequence samples in the Stanford HIV Drug Resistance database contained measurements of VL at six time points during the course of twenty-four week therapy. The drugs used in various single and combination therapies as well as the numbers of HIV-1 individuals taking the therapy are shown in Table 1. As described in the methods section, the first classification method for responders and non-responders, SD or Standard Datenum [39], was based on the fold-change of the entire patient database between the 0 and 8 week time points. The SD method classifies patients as responders if their viral load decreases by 100-fold over this time period. All other patients are labelled as non-responders. As shown in Figure 1A, this led to binomial distribution with clear peaks identified for responders and non-responders. The second method for phenotype classification, Incremental Reduction (IR), is based on patients having a reduction of viral load in four out of six weeks. Figure 1B shows the sub populations of responders and non-responders for this classification as a function of VL at three different instances in the clinical trial. It is clear from the figure that responders move towards zero VL whereas non-responders are much less mobile in this setting. The third method for phenotype classification (BM) was based on the observation that viral load reduction after 24 weeks of therapy exhibited a bimodal distribution (Figure 1C). This method used a cutoff of 2,000 copies/mL to differentiate between responders and non-responders. Subpopulations corresponding to each drug regimen shown in Table 1 also exhibited similar bimodal distributions.

Table 1. Therapy Classification

thumbnailFigure 1. Responder Classifications. A graphical representation of the three phenotype classification methods: Standard Datenum (SD), Incremental Reduction (IR) and Bimodal classification (BM). Figure 1A: SD, A histogram showing the log10 change in viral load of all patients in the database. Patients labelled as "responders" are marked in pink and non-responders in "blue". Figure 1B: IR, Three scatter plots representing the viral load vs. CD4 counts for all patients in the database after 8, 12, and 24 weeks of therapy. Patients which decreased in viral load in 75% of their visits are labelled as "responders" and marked in pink; those that did not are labelled as "non-responders" and marked in blue. Figure 1C: BM, A histogram of the change in viral load after 24 weeks of therapy. Those patients that decreased by more than 2000 copies/ml were labelled as "responders" and are marked in pink; those that did not were labelled as "non-responders" and are marked in blue.

The overlap between these three methods is shown in the Venn diagram in Figure 2. More than half of the responders from each method are also declared responders by the other two methods. However, 244 of the 925 patients labelled as responders by the SD method at eight weeks are not considered responders after 24 weeks by the BM. This suggests that after a strong initial response to therapy, some patients regress between 8th and 24th week of intervention with antiretroviral drugs. We used these three clinically relevant phenotype classification methods to identify sequence motifs associated with the responder group in each classification.

thumbnailFigure 2. Venn Diagram. Venn diagram showing the intersection between responder sets corresponding to SD, IR, and BM classification.

Conserved linear motifs along HIV and their correlation with response to antiretroviral drugs

Our results show that the HIV sequence, although highly variant in nucleotide sequence, expresses eukaryotic linear motifs (ELMs) that are largely conserved over hundreds of subtype B and subtype C sequences, as shown in Figure 3. The motifs recognized in globular domain regions are not shown as they are less likely to be instrumental in the interactions of HIV-1 proteins with host targets. The figure illustrates the presence of ELMs at high density along the flexible, domain-free regions of the HIV proteins. ELMs found on HIV proteins are largely conserved in frequency of appearance in eukaryotic proteomes (unpublished observations) and as such these motifs are good candidates in feature selection for predicting response to antiviral drugs.

thumbnailFigure 3. Feature Annotation. Annotation of a short linear motifs (Eukareotytic Linear Motifs, miRNAs binding sites, human transcription factor binding sites) along the viral sequence for 100 subtype C and 500 subtype B sequences. The colour code is as follows: homology Islands (green), human miRNA binding-sites (blue), human TF sites (silver), cleavage ELMs (red), ligation ELMs (purple), modification ELMs (brown), and export ELMs (pink). The clinically annotated sequence region is shown in the black box.

We used Step-Wise Logistic Regression (SWLR) to classify patients into responder or non-responder categories based on the presence or absence of ELMs, miRNA binding sites, TF binding sites, and resistance sites, collectively referred to as features. SWLR employs an iterative algorithm to determine which features should be included in the final logistic regression model [44]. In brief, the algorithm starts with an initial group of features and fits a logistic regression model. It then discards any features with a near zero coefficient and determines which of the excluded features may have a non-zero coefficient if added to the model. This process repeats until it converges to a solution; In our experience this occurs within 100 iterations.

We used SWLR in 500 iterations of training and testing at equal proportions for all responder/non-responder samples shown in Table 1. The resulting Receiver Operator Characteristics (ROC) curves for IR classification for the therapy regimens presented in Table 1 are shown in Figure 4. These ROC curves show high prediction accuracy of responders with the features used in the model. The area under the ROC curve (AUC) is an indicator of the combined sensitivity (ability to detect true positives) and specificity (ability to detect true negatives) of the model. As shown in Figure 4, random mixing of the responder and non-responder populations by 20% drastically reduced AUC for all drug regimens. Random mixing by 50% resulted in AUC values nearly equal to 0.5 as would be expected for randomly selected populations. These results confirm the utility of the selected features for predicting responder/non-responder identity using logistic regression.

thumbnailFigure 4. ROC Curves. Receiver Operator Characteristic (ROC) curves determined by the stepwise-logistic regression (SWLR) for the therapy regimens presented in Table 1 using the IR classification. The BOLD blue shows the average ROC curve over 500 iterations. The solid black line indicates the prediction ability with 20% shuffling of the responder v non-responder categories. The dashed line indicates the corresponding averages of completely shuffled responder vs. non-responder categories.

The AUC values for all three phenotype classification methods are shown in Table 1. Note that AUC values for BM and IR phenotype classifications are similar and point to high accuracy of prediction of outcome with these classification methods. The SD method, on the other hand, gave AUC values that were somewhat smaller than the other two methods. It is possible that the feature set used in our SD analysis is not optimal for predicting responders after eight weeks of therapy.

Regression Coefficients

The average number of regression coefficients (features) found significant over 500 training/testing iterations ranged from five to ten, depending on the drug regimens presented in Table 1. These features corresponeded to two specific resistance sites (RS)s and ELMs. In a set of control SWLR computations, we used other motifs such as human transcription binding site motifs and miRNA binding motifs on the RT sequence, but none of them were found to be significant in regression. Shown in Figure 5 are regression coefficients with absolute values greater than 0.5 for the three phenotype classifications: SD (Figure 5A), IR (Figure 5B), and BM (Figure 5C). Note that the two resistance sites on the figure are highly predictive of outcome in single drug regimens such as AZT and DDI targeting RT along with the ELMs that overlap this part of the sequence. Mutation RS V108 is a strong indicator of poor response to AZT, DDI, 3TC, and AZT, 3TC combination at 8 weeks (SD classification) whereas RS M36 has a negative effect on a larger spectrum of drug combinations (Figure 5A). These two resistance sites are the only ones that emerged in the set of features that are highly correlated with response to antiretroviral drugs. However, the regression does not lose accuracy when resistance sites are excluded from the features used in the analysis (data not shown). In this restricted set the significance of ELMs overlapping the resistance sites increases to compensate for the deletion, confirming the important role this sequence region plays in signalling resistance to some of the antiretrovirals targeting RT. Our findings point to resistance sites (or overlapping ELMs) having strong correlation to response to single antiretroviral therapies, but response to HAART therapies are correlated strongly with functional host protein motifs that are also expressed by the RT.

thumbnailFigure 5. SWLR Feature Regression Coefficients. Heatmaps indicating the average of the SWLR regression coefficient for the motifs used in the classification. Blue colour in the ruler bar indicates that presence of an ELM motif creates greater likelihood of being in the responder category (R ELM) whereas red indicates greater likelihood of being in the non-responder category (NR ELM). Top Panel: SD; Middle Panel: IR, Bottom Panel: BM.

One of the most consistent predictors of positive outcome across therapy regimens is the presence of ELM-Lig-SH3-3 (Figure 5A). This is the motif recognized by the SH3 domains of host proteins with a non-canonical class II recognition capacity [45]. The SH3 domain is a protein-protein interaction module commonly found in intracellular signalling and adaptor proteins. The SH3 domains of multiple endocytic proteins have been recently implicated in binding ubiquitin, which serves as a signal for diverse cellular processes including protein destruction [45].

The two resistance sites and the ELMs that overlap them continue to be predictors of negative outcome in terms of response to subsets of antiretroviral therapies in phenotype classification based on incremental reduction of the VL (Figure 5B, IR Classification). In this case, the consistent positive predictor is the motif ELM-Lig-MAPK-1. MAPK interacting molecules that carry this docking motif help to regulate specific interaction in the MAPK cascade [46,47]. It is feasible that human MAPK is recognizing the ELM on these RT proteins, decreasing their efficacy through phosphorylation or other inhibition methods.

Figure 5C, showing the BM classification method, reveals the resistance site M36 as a consistent indicator for negative response and ELM-Lig-SH2-STAT 5 as a strong indicator for positive response to antiretroviral therapy. This ELM is a motif recognized by proteins that have a significant impact on innate immunity during sepsis [48]. The innate immune system provides immediate defence against infection and serves as the first line of host defence during infection [49]. Recent research point to the depletion of white blood cells associated with innate immunity and their recovery under HAART [50].

Among the host proteins that have been documented to interact with the HIV RT protein, those that have at least one of the ELMs shown in Figure 5 are presented in Table 2. The table contains 33 host proteins with varying functions closely related to the immune response and signalling. The most common gene ontology categories [51] and KEGG pathways [52] among these proteins include adenyl ribonucleotide binding, phosphorylation, cell death, and apoptosis and pathways such as natural killer cell mediated cytotoxicity and the MAPK signalling pathway (Table 3). Our present knowledge of the grammar of protein interactions between the host and the virus does not allow us to draw definitive models of the network of interactions that differentiates responders from non-responders in HAART therapies. Nonetheless, the results presented above provide a start towards constructing a plausible mechanism of how viral and host genotypes affect response to antiretroviral therapies

Table 2. Interacting Proteins

Table 3. Biological Context

Discussion

The deadly course of HIV infection eventually leading to AIDS and associated opportunistic infections has been altered for a majority of individuals under HAART therapies thanks to combination antiretroviral therapies. These therapies have also reduced viral load dramatically in most patients, rendering them much less effective in transmitting the virus to others [53]. Research has focused on discovering new drugs targeting HIV proteins as well as on identifying host proteins necessary for viral growth as further possible targets for drugs. However, the interaction between the viral and host genotypes jointly affecting an individual's response to antiretroviral drugs has not been fully explored.

In this study we hypothesized that those host sequence motifs that are involved in protein-protein and protein/DNA/RNA interactions and also found in viral genomes are features that could play important roles in determining HIV-1 disease progression. Our prediction technique determines whether a particular therapy regimen is complementary to the sequence profile of each patient. Our thinking is motivated by the accumulating experimental evidence that viruses utilize motifs found in the host genome and proteins for integrating into host cell molecular networks and hijacking their function for viral replication [54,55]. Using linear sequence motifs shared by both the host and the virus provides an approach for investigating the plausible mechanisms of host virus interactions and suggesting those that may be altered by antiretroviral drugs.

We have used known resistance sites and host motifs found on HIV reverse transcriptase as features for differentiating responders from non-responders (or weak responders) in stepwise logistic regression for 16 different combinations of antiretroviral drug regimens containing at least one drug against HIV reverse transcriptase. Responder phenotype was defined multiple ways to gain insights into drug response at 8 weeks (SD phenotype classification) and 24 weeks (BM phenotype classification) after the beginning of the therapy and somewhere in between (IR Phenotype classification). Host motifs that appear to be highly relevant to viral replication such as the transcription site binding motifs [56,57] and miRNA binding site sequence motifs [58,59] could not be included into the analysis because these motifs are not contained within the RT region. Two resistance sites on HIV RT were found to be indicators of negative outcome, especially for regimens consisting of antiretrovirals targeting RT, but their influence was lower in HAART therapies. For the HAART therapy cases, the ELMS that contained these two resistance sites could be deleted from the model without sacrificing prediction accuracy. On the other hand, a number of ELMs were strongly correlated with positive outcome at different stages of antiretroviral therapy. These ELMs were associated with binding events leading to phosphorylation, ubiquination and the innate immune response.

Our approach to relate HIV sequence motifs to the course of infection does not require a priori information about how the HIV sequence would mutate in the presence of antiretroviral drugs. We were able to make accurate predictions without the resistance site information available in the literature. The input to our machine learning algorithm is simply the HIV sequence. We use publicly available bioinformatics tools to annotate these sequences with host motifs relevant to outcome. We then identify the motifs on the sequence that differentiate between responders and non-responders. These motifs can then be linked to specific viral host protein interactions and the pathways of these interactions. The promise of our approach will be fully explored with the availability of clinically annotated HIV whole genome sequences obtained at different time points during HAART therapy.

Conclusion

Linear binding motifs found in both the host and viral proteomes constitute a set of features highly predictive of response to therapy involving different combinations of antiretroviral drugs. Stepwise logistic regression as used here utilizes only the HIV-1 sequence and does not require annotations of resistance sites specific to various antiretroviral drugs. This study emphasizes finding sequence motifs which facilitate binding between viral and host proteins. This binding may allow the hijacking of host protein binding sites from their usual binding partners and thus alter the signalling pathways of the host cell. Our study points to competitive binding of HIV proteins to host proteins using motifs found in the host as the mechanism of interplay between the host and pathogen genotypes in dictating response to therapy. Our method is applicable to other viral infections where the viral sequence is known but resistance sites to antiviral therapies have not yet been documented.

Methods

Data sources for HIV1 sequences and clinical phenotype assignment

This study utilizes sequence and clinical data from two distinct sources. All whole genome HIV-1 sequences were downloaded from the Los Alamos HIV Sequence Database http://www.hiv.lanl.gov/ webcite in order to get a motif expression map of the whole genome. As of 9/1/2006, this dataset consisted of 1,112 subtype B and 922 subtype C whole genome sequences, along with a smaller number of samples from other subtypes. This dataset also contained five reference sequences each for alignment of subtypes B and C.

We used data from the Stanford HIV Drug Resistance Database [33] in order to investigate the clinical relevance of host protein and DNA motifs on the RT region of the HIV-1 sequence. The Stanford database curates clinical information from drug trials on large HIV cohorts and associates them with the sequence coding the protein targeted by the drug. As of 11/15/2008, the database contained few PR region sequences. However, the dataset contained 2,019 RT sequences annotated with clinical parameters such as CD4 counts, VLs and the specific antiretroviral therapy as shown in Table 1. Each patient in this subset had at least 1 sequence fragment from RT, had 4 or more CD4 and VL measurements at 0, 2, 4, 8, 12, and 24 weeks during the course of a constant therapy regimen.

Phenotype Classification

We focused on VL in the responder/non-responder classification [31] and examined the patient population using three methods of responder/non-responder classification: Standard Datenum (SD), Incremental Reduction (IR) and BiModal Classification (BM). The Standard Datenum method labels patients as responders if their VL decreases by 100-fold over 8 weeks of therapy [39]. The reduction in VL over the 24 week period logged by the Stanford HIV Drug Resistance Database exhibited a bimodal distribution for the patient population. Parameters of this distribution were obtained using the expectance maximization method described in [32] and indicated that a reduction of 2000 copies/mL in viral load would accurately split the responder and non-responder distributions. We refer to this method as bimodal classification. The third method we used was designed to avoid potential noise issues that could arise from relying the VL measurement on a single clinical visit [60]. The phenotype classification according to incremental reduction of the VL is such that if a patient's VL decreases between at least four visit pairs, then those patients are labelled as responders.

Linear Motifs on HIV Genome and Proteome and Resistance Sites

Our classification method uses the presence and absence of short linear motifs on the HIV genome. These motifs can be grouped into three basic types: eukaryotic linear motifs (ELMs), nucleotide-based motifs and a priori-based resistance mutations. In order to evaluate the relative positions of nucleotide motifs and protein motifs on the same platform, we annotated the protein motifs back to their corresponding nucleotide positions. This could create some ambiguity since HIV has multiple overlapping reading frames. However, our clinical dataset only contained sequences from the RT region. We used a local BLASTx query [61] on a database of HIV-1 subtype B and C reference samples to translate the nucleotide fragments into their corresponding protein sequences (see Additional file 1). This ensured the proper translation even if the start and stop codons were missing from the sequence.

Additional file 1. Source Code. All python and MATLAB code required to produce the figures and tables shown in the manuscript. Documentation and unit-tests are provided to facilitate their usage.

Format: ZIP Size: 40KB Download fileOpen Data

The first feature group consisted of ELM ligation sites and subcellular targeting sequences. These were identified on HIV-1 protein amino acid sequences using the ELM webserver tool [36]. The webserver tool filters out ELMs that fall into the globular regions proteins due to their predicted location within the 3D structure of the protein [36]. The second feature group consisted of HIV-1 sequence motifs that corresponded to annotated human transcription factor (TF) binding site motifs and miRNA binding sites. We used the MATCH™ web server [62] to annotate the TF binding sites on HIV-1 sequences with the public version of the TRANSFAC ® database as of 11/14/08 [34]. We required a core similarity of 0.75 and a global similarity of 0.70 in parameter assignment and chose among alternatives the method that minimized false negatives [62]. For the annotation of miRNA binding sites, recognition sequences for human miRNA were obtained from a human miRNA database [35]. As of 11/14/08 this database contained 417 experimentally verified human miRNA binding recognition sequences. The HIV sequences were scanned using the RNAhybrid program [63] and the background parameters of the extreme value distribution were created from 1,000 random sequences with dinucleotide distributions identical to our compiled HIV-1 sequence database [63]. Any binding site which had a p < 0.01 was annotated as a potential miRNA binding site. The third group of features consisted of resistance mutation sites on HIV sequence [64]. In order to capture the known HIV-1 therapy resistance mutation sites on the amino acid sequence of RT, we created regular expressions similar to ELMs which identify the known resistance conferring mutations (RSs) from the Stanford HIV-1 Resistance Database [33].

Predicting Therapy Outcome

We used stepwise logistic regression (SWLR) to assess the potential of the extracted short linear sequence features along the RT sequence in differentiating between responders and non responders [44]. SWLR was implemented in the MATLAB™ 2007b Statistics toolbox [65] (see Additional file 1). This regression method employs an iterative algorithm to determining the features that should be included in a predictive model. We used p-value < 0.01 as an entrance cutoff and p-value > 0.1 as a removal cutoff. In our study the algorithm converged to a final solution within 50–200 iterations.

SWLR algorithm was applied to differentiate responders from non-responders in three different assignments of the phenotype for 500 iterations of 2-fold cross validation. Since the efficiency of the SWLR algorithm is sensitive to the class composition of the training data [44] we ensured that each training set consisted of roughly 50% responders and 50% non-responders. After each set of training we determined the specificity and sensitivity of our classifier on the independent testing data and plotted the receiver operator characteristics (ROC) curve for each iteration in our scheme. The area under the ROC curve (AUC) represents the likelihood that one can identify a responder accurately using the method. This procedure was performed independently for each therapy regime under consideration and for the whole population shown in Table 1.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

WD wrote the manuscript with AT and performed the analysis with PE. LU provided technical insight. All authors approved the final manuscript.

Acknowledgements

We would like to acknowledge Dr. Fatah Kashanchi and Dr. Louis Mansky for their critical review of our manuscript and input towards its improvement. The study was supported by the National Institute of Health (NIH) grant #232240 and by the National Science Foundation (NSF) grant # 235327. Additional support came from a Drexel University Calhoun Fellowship (WD) and from NIH training grant T32 HG000046 (PE).

References

  1. Frankel AD, Young JAT: HIV-1: Fifteen Proteins and an RNA.

    Annual Review of Biochemistry 1998, 67(1):1-25. PubMed Abstract | Publisher Full Text OpenURL

  2. Los Alamos HIV-1 Sequence Database [http://www.hiv.lanl.gov/] webcite

  3. Grant RM, Hecht FM, Warmerdam M, Liu L, Liegler T, Petropoulos CJ, Hellmann NS, Chesney M, Busch MP, Kahn JO: Time trends in primary HIV-1 drug resistance among recently infected persons.

    Jama 2002, 288(2):181-188. PubMed Abstract | Publisher Full Text OpenURL

  4. Kati WM, Johnson KA, Jerva LF, Anderson KS: Mechanism and fidelity of HIV reverse transcriptase.

    Journal of Biological Chemistry 1992, 267(36):25988-25997. PubMed Abstract | Publisher Full Text OpenURL

  5. Kuhner MK, Yamato J, Felsenstein J: Estimating effective population size and mutation rate from sequence data using Metropolis-Hastings sampling.

    Genetics 1995, 140(4):1421-1430. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Brass AL, Dykxhoorn DM, Benita Y, Yan N, Engelman A, Xavier RJ, Lieberman J, Elledge SJ: Identification of Host Proteins Required for HIV Infection Through a Functional Genomic Screen.

    AAAS 2008, 319:921. OpenURL

  7. De Clercq E: HIV inhibitors targeted at the reverse transcriptase.

    AIDS research and human retroviruses 1992, 8(2):119-134. Publisher Full Text OpenURL

  8. Deeks SG, Smith M, Holodniy M, Kahn JO: HIV-1 protease inhibitors. A review for clinicians.

    JAMA 1997, 277(2):145-153. PubMed Abstract | Publisher Full Text OpenURL

  9. Pommier Y, Pilon AA, Bajaj K, Mazumder A, Neamati N: HIV-1 integrase as a target for antiviral drugs.

    Antiviral chemistry & chemotherapy 1997, 8(6):463-483. OpenURL

  10. Nair V: REVIEW HIV integrase as a target for antiviral chemotherapy.

    Review of Medical Virology 2002, 12:179-193. Publisher Full Text OpenURL

  11. Pommier Y, Johnson AA, Marchand C: Integrase inhibitors to treat HIV/Aids.

    Nature 2005, 4(3):236-248. Publisher Full Text OpenURL

  12. Rambaut A, Posada D, Crandall KA, Holmes EC: The causes and consequences of HIV evolution.

    Nature Reviews Genetics 2004, 5(1):52-61. PubMed Abstract | Publisher Full Text OpenURL

  13. Johnson VA, Brun-Vezinet F, Clotet B, Gunthard HF, Kuritzkes DR, Pillay D, Schapiro JM, Richman DD: Update of the drug resistance mutations in HIV-1: 2007.

    Top HIV Medicine 2007, 15(4):119-125. OpenURL

  14. D'Aquila RT, Schapiro JM, Brun-Vézinet F, Clotet B, Md PD, Conway B, Demeter LM, Grant RM, Johnson VA, Kuritzkes DR: Drug Resistance Mutations in HIV-1.

    Top HIV Medicine 2002., 10(5) OpenURL

  15. Johnson VA, Brun-Vézinet F, Clotet B, Conway B, Md RTD, Demeter LM, Kuritzkes DR, Pillay D, Schapiro JM, Telenti A: Update of the Drug Resistance Mutations in HIV-1: 2004.

    Top HIV Medicine 2004, 292:119-24. OpenURL

  16. Johnson VA, Brun-Vezinet F, Clotet B, Kuritzkes DR, Pillay D, Schapiro JM, Richman DD: Update of the drug resistance mutations in HIV-1: Fall 2006.

    Top HIV Medicine 2006, 14(3):125-130. OpenURL

  17. Katz MH, Schwarcz SK, Kellogg TA, Klausner JD, Dilley JW, Gibson S, McFarland W: Impact of highly active antiretroviral treatment on HIV seroincidence among men who have sex with men: San Francisco.

    American journal of public health 2002, 92(3):388-394. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Mansky LM: The mutation rate of human immunodeficiency virus type 1 is influenced by the vpr gene.

    Virology 1996, 222(2):391-400. PubMed Abstract | Publisher Full Text OpenURL

  19. Mocroft A, Gill MJ, Davidson W, Phillips AN: Predictors of a viral response and subsequent virological treatment failure in patients with HIV starting a protease inhibitor.

    AIDS (London, England) 1998, 12(16):2161. PubMed Abstract | Publisher Full Text OpenURL

  20. Deeks SG: Treatment of antiretroviral-drug-resistant HIV-1 infection.

    The Lancet 2003, 362(9400):2002-2011. Publisher Full Text OpenURL

  21. Lucas GM, Chaisson RE, Moore RD: Highly active antiretroviral therapy in a large urban clinic: risk factors for virologic failure and adverse drug reactions.

    Annals of internal medicine 1999, 131(2):81-87. PubMed Abstract | Publisher Full Text OpenURL

  22. Scheer S, Chu PL, Klausner JD, Katz MH, Schwarcz SK: Effect of highly active antiretroviral therapy on diagnoses of sexually transmitted diseases in people with AIDS.

    Lancet 2001, 357(9254):432-435. PubMed Abstract | Publisher Full Text OpenURL

  23. Pinney JW, Dickerson JE, Fu W, Sanders-Beer BE, Ptak RG, Robertson DL: HIV-host interactions: a map of viral perturbation of the host system.

    AIDS (London, England) 2009, 23(5):549-554. PubMed Abstract | Publisher Full Text OpenURL

  24. Fu W, Sanders-Beer BE, Katz KS, Maglott DR, Pruitt KD, Ptak RG: Human immunodeficiency virus type 1, human protein interaction database at NCBI.

    Nucleic acids research 2009, (37 Database):D417. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Ptak RG, Fu W, Sanders-Beer BE, Dickerson JE, Pinney JW, Robertson DL, Rozanov MN, Katz KS, Maglott DR, Pruitt KD: Cataloguing the HIV-1 human protein interaction network.

    AIDS Research and Human Retroviruses 2008 2008, 24(12):1497-1502. Publisher Full Text OpenURL

  26. König R, Zhou Y, Elleder D, Diamond TL, Bonamy GMC, Irelan JT, Chiang C, Tu BP, De Jesus PD, Lilley CE: Global analysis of host-pathogen interactions that regulate early-stage HIV-1 replication.

    Cell 2008, 135(1):49-60. PubMed Abstract | Publisher Full Text OpenURL

  27. Brass AL, Dykxhoorn DM, Benita Y, Yan N, Engelman A, Xavier RJ, Lieberman J, Elledge SJ: Identification of host proteins required for HIV infection through a functional genomic screen.

    Science (New York, NY) 2008, 319(5865):921-926. OpenURL

  28. Fellay J, Shianna KV, Ge D, Colombo S, Ledergerber B, Weale M, Zhang K, Gumbs C, Castagna A, Cossarizza A: A whole-genome association study of major determinants for host control of HIV-1.

    Science (New York, NY) 2007, 317(5840):944. OpenURL

  29. Stauber RH, Pavlakis GN: Intracellular Trafficking and Interactions of the HIV-1 Tat Protein.

    Virology 1998, 252(1):126-136. PubMed Abstract | Publisher Full Text OpenURL

  30. Connor RI, Sheridan KE, Ceradini D, Choe S, Landau NR: Change in Coreceptor Use Correlates with Disease Progression in HIV-1-Infected Individuals.

    Journal of Experimental Medicine 1997, 185(4):621-628. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Moore DM, Awor A, Downing R, Kaplan J, Montaner JS, Hancock J, Were W, Mermin J: CD4+ T-Cell Count Monitoring Does Not Accurately Identify HIV-Infected Adults With Virologic Failure Receiving Antiretroviral Therapy.

    Journal of Acquired Immune Deficiency Syndromes 2008, 48(5):477-484. OpenURL

  32. Ertel A, Tozeren A: Switch-like genes populate cell communication pathways and are enriched for extracellular proteins.

    BMC Genomics 2008, 9:3. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  33. Rhee SY, Gonzales MJ, Kantor R, Betts BJ, Ravela J, Shafer RW: Human immunodeficiency virus reverse transcriptase and protease sequence database.

    Nucleic acids research 2003, 31(1):298-303. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, et al.: TRANSFAC: transcriptional regulation, from patterns to profiles.

    Nucleic acids research 2003, 31(1):374-378. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression.

    Nucleic acids research 2008, (36 Database):D149-153. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Puntervoll P, Linding R, Gemund C, Chabanis-Davidson S, Mattingsdal M, Cameron S, Martin DM, Ausiello G, Brannetti B, Costantini A, et al.: ELM server: A new resource for investigating short functional sites in modular eukaryotic proteins.

    Nucleic acids research 2003, 31(13):3625-3630. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Kadaveru K, Vyas J, Schiller MR: Viral infection and human disease – insights from minimotifs.

    Front Biosci 2008, 13:6455-6471. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Larder B, Wang D, Revell A, Montaner J, Harrigan R, De Wolf F, Lange J, Wegner S, Ruiz L, Pérez-Elías MJ: The development of artificial neural networks to predict virological response to combination HIV therapy.

    Antiviral therapy 2007, 12(1):15. PubMed Abstract OpenURL

  39. Rosen-Zvi M, Altmann A, Prosperi M, Aharoni E, Neuvirth H, Sonnerborg A, Schulter E, Struck D, Peres Y, Incardona F: Selecting anti-HIV therapies based on a variety of genomic and clinical factors.

    Bioinformatics (Oxford, England) 2008, 24(13):i399. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Nanni L, Lumini A: MppS: An ensemble of support vector machine based on multiple physicochemical properties of amino acids.

    Neurocomputing 2006, 69(13–15):1688-1690. Publisher Full Text OpenURL

  41. Beerenwinkel N, Daumer M, Oette M, Korn K, Hoffmann D, Kaiser R, Lengauer T, Selbig J, Walter H: Geno2pheno: Estimating phenotypic drug resistance from HIV-1 genotypes.

    Nucleic acids research 2003, 31(13):3850-3855. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Beerenwinkel N, Lengauer T, Daumer M, Kaiser R, Walter H, Korn K, Hoffmann D, Selbig J: Methods for optimizing antiviral combination therapies.

    Bioinformatics (Oxford, England) 2003, 19(Suppl 1):i16-25. PubMed Abstract | Publisher Full Text OpenURL

  43. Vermeiren H, Van Craenenbroeck E, Alen P, Bacheler L, Picchio G, Lecocq P: Prediction of HIV-1 drug susceptibility phenotype from the viral genotype using linear regression modeling.

    Journal of Virological Methods 2007, 145(1):47-55. PubMed Abstract | Publisher Full Text OpenURL

  44. Draper NR, Smith H: Applied Regression Analysis. New York; Wiley-Interscience; 1967.

  45. He Y, Hicke L, Radhakrishnan I: Structural basis for ubiquitin recognition by SH3 domains.

    Journal of molecular biology 2007, 373(1):190-196. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Biondi RM, Nebreda AR: Signalling specificity of Ser/Thr protein kinases through docking-site-mediated interactions.

    The Biochemical journal 2003, 372(Pt 1):1-13. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Jacobs D, Glossip D, Xing H, Muslin AJ, Kornfeld K: Multiple docking sites on substrate proteins form a modular system that mediates recognition by ERK MAP kinase.

    Genes & development 1999, 13(2):163-175. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Matsukawa A: STAT proteins in innate immunity during sepsis: lessons from gene knockout mice.

    Acta medica Okayama 2007, 61(5):239-245. PubMed Abstract | Publisher Full Text OpenURL

  49. Levy JA: The importance of the innate immune system in controlling HIV infection and disease.

    Trends in Immunology 2001, 22(6):312-316. PubMed Abstract | Publisher Full Text OpenURL

  50. Li D, Xu XN: NKT cells in HIV-1 infection.

    Cell research 2008, 18(8):817-822. PubMed Abstract | Publisher Full Text OpenURL

  51. Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, et al.: The Gene Ontology (GO) database and informatics resource.

    Nucleic acids research 2004, (32 Database):D258-261. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes.

    Nucleic acids research 2000, 28(1):27-30. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Castilla J, Jorge del Romero MD, Hernando V, Marincovich B, García S, Rodríguez C: Effectiveness of Highly Active Antiretroviral Therapy in Reducing Heterosexual Transmission of HIV.

    Journal of Acquired Immune Deficiency Syndromes 2005, 40(1):96. PubMed Abstract | Publisher Full Text OpenURL

  54. Garber ME, Wei P, KewalRamani VN, Mayall TP, Herrmann CH, Rice AP, Littman DR, Jones KA: The interaction between HIV-1 Tat and human cyclin T1 requires zinc and a critical cysteine residue that is not conserved in the murine CycT1 protein.

    1998, 12(22):3512-3527. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Longo F, Marchetti MA, Castagnoli L, Battaglia PA, Gigliani F: A Novel Approach to Protein-Protein Interaction: Complex Formation between the P53 Tumor Suppressor and the HIV Tat Proteins.

    Biochemical and Biophysical Research Communications 1995, 206(1):326-334. Publisher Full Text OpenURL

  56. Van Lint C, Amella CA, Emiliani S, John M, Jie T, Verdin E: Transcription factor binding sites downstream of the human immunodeficiency virus type 1 transcription start site are important for virus infectivity.

    The Journal of Virology 1997, 71(8):6113-6127. OpenURL

  57. Rockman MV, Hahn MW, Soranzo N, Goldstein DB, Wray GA: Positive Selection on a Human-Specific Transcription Factor Binding Site Regulating IL4 Expression.

    Current Biology 2003, 13(23):2118-2123. PubMed Abstract | Publisher Full Text OpenURL

  58. Hariharan M, Scaria V, Pillai B, Brahmachari SK: Targets for human encoded microRNAs in HIV genes.

    Biochemical and Biophysical Research Communications 2005, 337(4):1214-1218. Publisher Full Text OpenURL

  59. Huang J, Wang F, Argyris E, Chen K, Liang Z, Tian H, Huang W, Squires K, Verlinghieri G, Zhang H: Cellular microRNAs contribute to HIV-1 latency in resting primary CD4 T lymphocytes.

    Nature Medicine 2007, 13:1241-1247. PubMed Abstract | Publisher Full Text OpenURL

  60. Mulder J, McKinney N, Christopherson C, Sninsky J, Greenfield L, Kwok S: Rapid and simple PCR assay for quantitation of human immunodeficiency virus type 1 RNA in plasma: application to acute retroviral infection.

    Journal of clinical microbiology 1994, 32(2):292-300. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

    Nucleic acids research 1997, 25(17):3389-3402. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Kel AE, Gossling E, Reuter I, Cheremushkin E, Kel-Margoulis OV, Wingender E: MATCH: A tool for searching transcription factor binding sites in DNA sequences.

    Nucleic acids research 2003, 31(13):3576-3579. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes.

    RNA (New York, NY) 2004, 10(10):1507-1517. OpenURL

  64. Shafer RW, Schapiro JM: HIV-1 drug resistance mutations: an updated framework for the second decade of HAART.

    AIDS reviews 2008, 10(2):67-84. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  65. MATLAB 2007b [http://www.mathworks.com] webcite

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1755-8794/2/47/prepub