Analysis of Physicochemical and Structural Properties Determining HIV-1 Coreceptor Usage

The relationship of HIV tropism with disease progression and the recent development of CCR5-blocking drugs underscore the importance of monitoring virus coreceptor usage. As an alternative to costly phenotypic assays, computational methods aim at predicting virus tropism based on the sequence and structure of the V3 loop of the virus gp120 protein. Here we present a numerical descriptor of the V3 loop encoding its physicochemical and structural properties. The descriptor allows for structure-based prediction of HIV tropism and identification of properties of the V3 loop that are crucial for coreceptor usage. Use of the proposed descriptor for prediction results in a statistically significant improvement over the prediction based solely on V3 sequence with 3 percentage points improvement in AUC and 7 percentage points in sensitivity at the specificity of the 11/25 rule (95%). We additionally assessed the predictive power of the new method on clinically derived ‘bulk’ sequence data and obtained a statistically significant improvement in AUC of 3 percentage points over sequence-based prediction. Furthermore, we demonstrated the capacity of our method to predict therapy outcome by applying it to 53 samples from patients undergoing Maraviroc therapy. The analysis of structural features of the loop informative of tropism indicates the importance of two loop regions and their physicochemical properties. The regions are located on opposite strands of the loop stem and the respective features are predominantly charge-, hydrophobicity- and structure-related. These regions are in close proximity in the bound conformation of the loop potentially forming a site determinant for the coreceptor binding. The method is available via server under http://structure.bioinf.mpi-inf.mpg.de/. Author Summary: Human Immunodeficiency Virus (HIV) requires one of the chemokine coreceptors CCR5 or CXCR4 for entry into the host cell. The capacity of the virus to use one or both of these coreceptors is termed tropism. Monitoring HIV tropism is of high importance due to the relationship of the emergence of CXCR4-tropic virus with the progression of immunodeficiency and for patient treatment with the recently developed CCR5 antagonists. Computational methods for predicting HIV tropism are based on sequence and on structure of the third variable region (V3 loop) of the viral gp120 protein — the major determinant of the HIV tropism. Limitations of the existing methods include the limited insights they provide into the biochemical determinants of coreceptor usage, high computational load of the structure-based methods and low prediction accuracy on clinically derived patient samples. Here we propose a numerical descriptor of the V3 loop encoding the physicochemical and structural properties of the loop. The new descriptor allows for server-based prediction of viral tropism with accuracy comparable to that of established sequence-based methods both on clonal and clinically derived patient data as well as for the interpretation of the properties of the loop relevant for tropism. The server is available under http://structure.bioinf.mpi-inf.mpg.de/.