Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Immunologic and virologic characterization of an ART-treated HIV-I patients cohort with long-term control of viremia

Background Long-term treatment of primary HIV-1 infection (PHI) may allow the immune reconstitution of responses lost during the acute viremic phase and decrease of peripheral reservoirs. This in turn may represent the best setting for the use of therapeutic vaccines in order to lower the viral set-point or control of viral rebound upon ART discontinuation. Methods We investigated a cohort of 16 patients who started ART at PHI, with treatment duration of ≥4 years and persistent aviremia (<50 HIV-1 copies/ml). The cohort was characterized in terms of viral subtype, cell-associated RNA, proviral DNA and HLA genotype. Secretion of IFN-γ, IL-2 and TNF-α by CD8 T-cells was analysed by polychromatic flowcytometry using a panel of 192 HIV-1-derived epitopes. Results This cohort is highly homogenous in terms of viral subtype: 81% clade B. We identified 44 epitope-specific responses: all patients had detectable responses to >1 epitope and the mean number of responding epitopes per patient was 3. The mean frequency of cytokines-secreting CD8 T-cells was 0.32%. CD8 T-cells secreting simultaneously IFN-γ, IL-2 and TNF-α made up for about 40% of the response and cells secreting at least 2 cytokines for about 80%, consistent with a highly polyfunctional CD8 T-cell profile. There was no difference in term of polyfunctionality when HLA restriction, or recognized viral regions and epitopes were considered. Proviral DNA was detectable in all patients but at low levels (mean = 108 copies/1 million PBMCs) while cell-associated mRNA was not detectable in 19% of patients (mean = 11 copies/1 million PBMCs when detectable). Conclusion Patients with sustained virological suppression after initiation of ART at PHI show polyfunctional CD8 T-cell and low levels of proviral DNA with an absence of residual replication in a substantial percentage of patients. The use of therapeutic vaccines in this population may promote low level of rebound viremia or control of viral replication upon ART cessation.

Host Genes Important to HIV Replication and Evolution.

Recent years have seen a significant increase in understanding of the host genetic and genomic determinants of susceptibility to HIV-1 infection and disease progression, driven in large part by candidate gene studies, genome-wide association studies, genome-wide transcriptome analyses, and large-scale in vitro genome screens. These studies have identified common variants in some host loci that clearly influence disease progression, characterized the scale and dynamics of gene and protein expression changes in response to infection, and provided the first comprehensive catalogs of genes and pathways involved in viral replication. Experimental models of AIDS and studies in natural hosts of primate lentiviruses have complemented and in some cases extended these findings. As the relevant technology continues to progress, the expectation is that such studies will increase in depth (e.g., to include host whole exome and whole genome sequencing) and in breadth (in particular, by integrating multiple data types).

Deletion of the viral anti-apoptotic gene F1L in the HIV/AIDS vaccine candidate MVA-C enhances immune responses against HIV-1 antigens.

Vaccinia virus (VACV) encodes an anti-apoptotic Bcl-2-like protein F1 that acts as an inhibitor of caspase-9 and of the Bak/Bax checkpoint but the role of this gene in immune responses is not known. Because dendritic cells that have phagocytosed apoptotic infected cells cross-present viral antigens to cytotoxic T cells inducing an antigen-specific immunity, we hypothesized that deletion of the viral anti-apoptotic F1L gene might have a profound effect on the capacity of poxvirus vectors to activate specific immune responses to virus-expressed recombinant antigens. This has been tested in a mouse model with an F1L deletion mutant of the HIV/AIDS vaccine candidate MVA-C that expresses Env and Gag-Pol-Nef antigens (MVA-C-ΔF1L). The viral gene F1L is not required for virus replication in cultured cells and its deletion in MVA-C induces extensive apoptosis and expression of immunomodulatory genes in infected cells. Analysis of the immune responses induced in BALB/c mice after DNA prime/MVA boost revealed that, in comparison with parental MVA-C, the mutant MVA-C-ΔF1L improves the magnitude of the HIV-1-specific CD8 T cell adaptive immune responses and impacts on the CD8 T cell memory phase by enhancing the magnitude of the response, reducing the contraction phase and changing the memory differentiation pattern. These findings reveal the immunomodulatory role of F1L and that the loss of this gene is a valid strategy for the optimization of MVA as vaccine vector.

The Individualized Genetic Barrier Predicts Treatment Response in a Large Cohort of HIV-1 Infected Patients.

The success of combination antiretroviral therapy is limited by the evolutionary escape dynamics of HIV-1. We used Isotonic Conjunctive Bayesian Networks (I-CBNs), a class of probabilistic graphical models, to describe this process. We employed partial order constraints among viral resistance mutations, which give rise to a limited set of mutational pathways, and we modeled phenotypic drug resistance as monotonically increasing along any escape pathway. Using this model, the individualized genetic barrier (IGB) to each drug is derived as the probability of the virus not acquiring additional mutations that confer resistance. Drug-specific IGBs were combined to obtain the IGB to an entire regimen, which quantifies the virus' genetic potential for developing drug resistance under combination therapy. The IGB was tested as a predictor of therapeutic outcome using between 2,185 and 2,631 treatment change episodes of subtype B infected patients from the Swiss HIV Cohort Study Database, a large observational cohort. Using logistic regression, significant univariate predictors included most of the 18 drugs and single-drug IGBs, the IGB to the entire regimen, the expert rules-based genotypic susceptibility score (GSS), several individual mutations, and the peak viral load before treatment change. In the multivariate analysis, the only genotype-derived variables that remained significantly associated with virological success were GSS and, with 10-fold stronger association, IGB to regimen. When predicting suppression of viral load below 400 cps/ml, IGB outperformed GSS and also improved GSS-containing predictors significantly, but the difference was not significant for suppression below 50 cps/ml. Thus, the IGB to regimen is a novel data-derived predictor of treatment outcome that has potential to improve the interpretation of genotypic drug resistance tests.

Host genome influences on susceptibility to HIV-1

In vitro and in vivo analyses identified a significant component of heritability in cellular or host susceptibility to HIV-1. The bases for susceptibility can be traced to genetic differences (inter-species) resulting from evolutionary adaptation to exogenous (and endogenous) retroviral infections, and to intra-species and inter-individual (human) differences associated with genetic variation. We have completed large scale evolutionary analysis of genes involved in HIV life cycle and pathogenesis, as well as participating and conducting genome-wide association studies, linkage analysis, and transcriptome analysis. These studies allowed a better understanding of the influence of common human variants in HIV-1 susceptibility and define a number of experimental challenges in the filed: understanding of the role of rare and private mutations in susceptibility, and the development of better tools for the integration of data from large-scale studies.

Pharmacogenetics-based population pharmacokinetic analysis of etravirine in HIV-1 infected individuals.

OBJECTIVES: Etravirine (ETV) is metabolized by cytochrome P450 (CYP) 3A, 2C9, and 2C19. Metabolites are glucuronidated by uridine diphosphate glucuronosyltransferases (UGT). To identify the potential impact of genetic and non-genetic factors involved in ETV metabolism, we carried out a two-step pharmacogenetics-based population pharmacokinetic study in HIV-1 infected individuals. MATERIALS AND METHODS: The study population included 144 individuals contributing 289 ETV plasma concentrations and four individuals contributing 23 ETV plasma concentrations collected in a rich sampling design. Genetic variants [n=125 single-nucleotide polymorphisms (SNPs)] in 34 genes with a predicted role in ETV metabolism were selected. A first step population pharmacokinetic model included non-genetic and known genetic factors (seven SNPs in CYP2C, one SNP in CYP3A5) as covariates. Post-hoc individual ETV clearance (CL) was used in a second (discovery) step, in which the effect of the remaining 98 SNPs in CYP3A, P450 cytochrome oxidoreductase (POR), nuclear receptor genes, and UGTs was investigated. RESULTS: A one-compartment model with zero-order absorption best characterized ETV pharmacokinetics. The average ETV CL was 41 (l/h) (CV 51.1%), the volume of distribution was 1325 l, and the mean absorption time was 1.2 h. The administration of darunavir/ritonavir or tenofovir was the only non-genetic covariate influencing ETV CL significantly, resulting in a 40% [95% confidence interval (CI): 13-69%] and a 42% (95% CI: 17-68%) increase in ETV CL, respectively. Carriers of rs4244285 (CYP2C19*2) had 23% (8-38%) lower ETV CL. Co-administered antiretroviral agents and genetic factors explained 16% of the variance in ETV concentrations. None of the SNPs in the discovery step influenced ETV CL. CONCLUSION: ETV concentrations are highly variable, and co-administered antiretroviral agents and genetic factors explained only a modest part of the interindividual variability in ETV elimination. Opposing effects of interacting drugs effectively abrogate genetic influences on ETV CL, and vice-versa.

Molecular cloning of a novel human I-mfa domain-containing protein that differently regulates human T-cell leukemia virus type I and HIV-1 expression.

Regulation of viral genome expression is the result of complex cooperation between viral proteins and host cell factors. We report here the characterization of a novel cellular factor sharing homology with the specific cysteine-rich C-terminal domain of the basic helix-loop-helix repressor protein I-mfa. The synthesis of this new factor, called HIC for Human I-mfa domain-Containing protein, is controlled at the translational level by two different codons, an ATG and an upstream non-ATG translational initiator, allowing the production of two protein isoforms, p32 and p40, respectively. We show that the HIC protein isoforms present different subcellular localizations, p32 being mainly distributed throughout the cytoplasm, whereas p40 is targeted to the nucleolus. Moreover, in trying to understand the function of HIC, we have found that both isoforms stimulate in T-cells the expression of a luciferase reporter gene driven by the human T-cell leukemia virus type I-long terminal repeat in the presence of the viral transactivator Tax. We demonstrate by mutagenesis that the I-mfa-like domain of HIC is involved in this regulation. Finally, we also show that HIC is able to down-regulate the luciferase expression from the human immunodeficiency virus type 1-long terminal repeat induced by the viral transactivator Tat. From these results, we propose that HIC and I-mfa represent two members of a new family of proteins regulating gene expression and characterized by a particular cysteine-rich C-terminal domain.

The EuroSIDA study: Regional differences in the HIV-1 epidemic and treatment response to antiretroviral therapy among HIV-infected patients across Europe--a review of published results.

EuroSIDA is a pan-European observational study that follows 14,265 HIV-infected patients from 31 European countries, Israel and Argentina, of which 2,560 are patients from eastern Europe (EE). The study group has performed several analyses addressing regional differences in the HIV-epidemic across Europe, where all countries were divided into five regions: south, west central, north, east central Europe and EE. Significant regional differences in patients' characteristics and pattern of AIDS diagnoses were documented. More patients from EE were diagnosed with tuberculosis compared to other regions. Significantly fewer HIV-infected patients in EE, who fulfilled the criteria for starting combination antiretroviral therapy (cART), actually received cART as compared with other regions of Europe. Those, receiving cART in EE had a lower initial virologic response rate irrespectively of the regimen used, although it has improved within years. Besides, treatment failure was more common in this region. Thus, improvements in the clinical management of HIV patients in EE are urgently needed. Strategies include creating scientific collaborations for HIV clinicians as well as teaching clinicians about the most advanced HIV management at clinically oriented courses held in eastern Europe.

Susceptibility and adaptation to human TRIM5α alleles at positive selected sites in HIV-1 capsid.

Numerous in vitro studies attribute to human TRIM5α some modest anti-HIV-1 activity and human population studies suggest some differential effect of TRIM5α polymorphisms on disease progression. If the activity of TRIM5α were relevant in vivo, it could result in positive selection on the viral capsid. To address this issue, we identified 10 positively selected sites in HIV-1 capsid from multiple viral strains and generated 17 clade B viruses carrying a minor (i.e. low frequency) residue or an alanine at those positions. All recombinant viruses were susceptible to the modest effect of common human TRIM5α and allelic variants R136Q, and H419Y; H43Y and G249D TRIM5α were generally inactive. Increased sensitivity to TRIM5α was observed for some capsid variants, suggesting that minor residues are selected against in human populations. On the other hand, the modest potency of human TRIM5α does not translate in escape mutations in the viral capsid.

Long sequence duplications, repeats, and palindromes in HIV-1 gp120: length variation in V4 as the product of misalignment mechanism.

We have shown that indels in gp120 V4 are associated to the presence of duplicated and palindromic sequences, suggesting that they may be produced by strand-slippage misalignment mechanism. Indels in V4 involved region-specific duplications 9 to 15 bp long, and repeats of various lengths, associated to trinucleotides AAT. No duplications were found in V3 and C3. The frequency of palindromic sequences in individual genes was found to be significantly higher in gp120 (p < or = 3.00E-7), and significantly lower in Tat (p < or = 9.00E-7) than the average frequency calculated over the full genome. The finding of elements of misalignment in association with indels in V4 suggests that these mutations may occur in proviral DNA after integration of HIV into the host genome. It also implies that occurrence of large indels in gp120 is not random but is directed by the presence and distribution of elements of misalignment in the HIV genome.

Impact of soluble CD26 on treatment outcome and hepatitis C virus-specific T cells in chronic hepatitis C virus genotype 1 infection.

BACKGROUND: Interferon and ribavirin therapy for chronic hepatitis C virus (HCV) infection yields sustained virological response (SVR) rates of 50-80%. Several factors such as non-1 genotype, beneficial IL28B genetic variants, low baseline IP-10, and the functionality of HCV-specific T cells predict SVR. With the pending introduction of new therapies for HCV entailing very rapid clearance of plasma HCV RNA, the importance of baseline biomarkers likely will increase in order to tailor therapy. CD26 (DPPIV) truncates the chemokine IP-10 into a shorter antagonistic form, and this truncation of IP-10 has been suggested to influence treatment outcome in patients with chronic HCV infection patients. In addition, previous reports have shown CD26 to be a co-stimulator for T cells. The aim of the present study was to assess the utility of CD26 as a biomarker for treatment outcome in chronic hepatitis C and to define its association with HCV-specific T cells. METHODS: Baseline plasma from 153 genotype 1 and 58 genotype 2/3 infected patients enrolled in an international multicenter phase III trial (DITTO-HCV) and 36 genotype 1 infected patients participating in a Swedish trial (TTG1) were evaluated regarding baseline soluble CD26 (sCD26) and the functionality of HCV-specific CD8(+) T cells. RESULTS: Genotype 1 infected patients achieving SVR in the DITTO (P = 0.002) and the TTG1 (P = 0.02) studies had lower pretreatment sCD26 concentrations compared with non-SVR patients. Sixty-five percent of patients with sCD26 concentrations below 600 ng/mL achieved SVR compared with 39% of the patients with sCD26 exceeding 600 ng/mL (P = 0.01). Patients with sCD26 concentrations below 600 ng/mL had significantly higher frequencies of HCV-specific CD8(+) T cells (P = 0.02). CONCLUSIONS: Low baseline systemic concentrations of sCD26 predict favorable treatment outcome in chronic HCV infection and may be associated with higher blood counts of HCV-specific CD8(+) T cells.

Minor protease inhibitor mutations at baseline do not increase the risk for a virological failure in HIV-1 subtype B infected patients.

BACKGROUND: Minor protease inhibitor (PI) mutations often exist as polymorphisms in HIV-1 sequences from treatment-naïve patients. Previous studies showed that their presence impairs the antiretroviral treatment (ART) response. Evaluating these findings in a larger cohort is essential. METHODS: To study the impact of minor PI mutations on time to viral suppression and time to virological failure, we included patients from the Swiss HIV Cohort Study infected with HIV-1 subtype B who started first-line ART with a PI and two nucleoside reverse transcriptase inhibitors. Cox regression models were performed to compare the outcomes among patients with 0 and ≥ 1 minor PI mutation. Models were adjusted for baseline HIV-1 RNA, CD4 cell count, sex, transmission category, age, ethnicity, year of ART start, the presence of nucleoside reverse transcriptase inhibitor mutations, and stratified for the administered PIs. RESULTS: We included 1199 patients of whom 944 (78.7%) received a boosted PI. Minor PI mutations associated with the administered PI were common: 41.7%, 16.1%, 4.7% and 1.9% had 1, 2, 3 or ≥ 4 mutations, respectively. The time to viral suppression was similar between patients with 0 (reference) and ≥ 1 minor PI mutation (multivariable hazard ratio (HR): 1.1 [95% confidence interval (CI): 1.0-1.3], P = .196). The time to virological failure was also similar (multivariable HR:.9 [95% CI:.5-1.6], P = .765). In addition, the impact of each single minor PI mutation was analyzed separately: none was significantly associated with the treatment outcome. CONCLUSIONS: The presence of minor PI mutations at baseline has no effect on the therapy outcome in HIV infected individuals.

Preclinical evaluation of the immunogenicity of C-type HIV-1-based DNA and NYVAC vaccines in the Balb/C mouse model.

As part of a European initiative (EuroVacc), we report the design, construction, and immunogenicity of two HIV-1 vaccine candidates based on a clade C virus strain (CN54) representing the current major epidemic in Asia and parts of Africa. Open reading frames encoding an artificial 160-kDa GagPolNef (GPN) polyprotein and the external glycoprotein gp120 were fully RNA and codon optimized. A DNA vaccine (DNA-GPN and DNA-gp120, referred to as DNA-C), and a replication-deficient vaccinia virus encoding both reading frames (NYVAC-C), were assessed regarding immunogenicity in Balb/C mice. The intramuscular administration of both plasmid DNA constructs, followed by two booster DNA immunizations, induced substantial T-cell responses against both antigens as well as Env-specific antibodies. Whereas low doses of NYVAC-C failed to induce specific CTL or antibodies, high doses generated cellular as well as humoral immune responses, but these did not reach the levels seen following DNA vaccination. The most potent immune responses were detectable using prime:boost protocols, regardless of whether DNA-C or NYVAC-C was used as the priming or boosting agent. These preclinical findings revealed the immunogenic response triggered by DNA-C and its enhancement by combining it with NYVAC-C, thus complementing the macaque preclinical and human phase I clinical studies of EuroVacc.