Coevolution Analysis of HIV-1 Envelope Glycoprotein Complex

Selection among site-dependent structurally constrained substitution models of protein evolution by approximate Bayesian computation

MOTIVATION

The selection among substitution models of molecular evolution is fundamental for obtaining accurate phylogenetic inferences. At the protein level, evolutionary analyses are traditionally based on empirical substitution models but these models make unrealistic assumptions and are being surpassed by structurally constrained substitution (SCS) models. The SCS models often consider site-dependent evolution, a process that provides realism but complicates their implementation into likelihood functions that are commonly used for substitution model selection.

RESULTS

We present a method to perform selection among site-dependent SCS models, also among empirical and site-dependent SCS models, based on the approximate Bayesian computation (ABC) approach and its implementation into the computational framework ProteinModelerABC. The framework implements ABC with and without regression adjustments and includes diverse empirical and site-dependent SCS models of protein evolution. Using extensive simulated data, we found that it provides selection among SCS and empirical models with acceptable accuracy. As illustrative examples, we applied the framework to analyze a variety of protein families observing that SCS models fit them better than the corresponding best-fitting empirical substitution models.

AVAILABILITY AND IMPLEMENTATION

ProteinModelerABC is freely available from https://github.com/DavidFerreiro/ProteinModelerABC, can run in parallel and includes a graphical user interface. The framework is distributed with detailed documentation and ready-to-use examples.

Regulation of epitope exposure in the gp41 membrane-proximal external region through interactions at the apex of HIV-1 Env

Glycoprotein Env of human immunodeficiency virus type 1 (HIV-1) mediates viral entry through membrane fusion. Composed of gp120 and gp41 subunits arranged as a trimer-of-heterodimers, Env adopts a metastable, highly dynamic conformation on the virion surface. This structural plasticity limits the temporospatial exposure of many highly conserved, neutralizing epitopes, contributing to the difficulty in developing effective HIV-1 vaccines. Here, we employed antibody neutralization of HIV-1 infectivity to investigate how inter- and intra-gp120 interactions mediated by variable loops V1/V2 and V3 at the Env apex regulate accessibility of the gp41 membrane-proximal external region (MPER) at the Env base. Swapping the V3 loop from Env_SF162 into the Env_HXB2 background shifted MPER exposure from the prefusogenic state to a functional intermediate conformation that was distinct from the prehairpin-intermediate state sensitive to gp41-targeted fusion inhibitors. The V3-loop swap had a profound impact on global protein dynamics, biasing the equilibrium to a closed conformation resistant to most anti-gp120 antibodies, stabilizing the protein to both cold- and soluble CD4-induced Env inactivation, and increasing the CD4 requirements for viral entry. Further dissection of the Env_HXB2 V3 loop revealed that residue 306 uniquely modulated epitope exposure and trimer stability. The R306S substitution substantially decreased sensitivity to antibodies targeting the gp41 MPER and, surprisingly, the gp120 V3-loop crown (residues 312–315), but had only modest effects on exposure of intervening gp120 epitopes. Furthermore, the point mutation reduced soluble CD4-induced inactivation, but had no impact on cold inactivation. The residue appeared to exert its effects by electrostatically modifying the strength of intra-subunit interactions between the V1/V2 and V3 loops. The distinct patterns of neutralization and stability pointed to a novel prefusogenic Env conformation along the receptor activation pathway and suggested that apical Env-regulation of gp41 MPER exposure can be decoupled from much of the dynamics of gp120 subunits.

Surface glycoprotein Env is the main target for neutralizing antibodies elicited by HIV-1 vaccines. Env spontaneously fluctuates among different structures, limiting exposure of many attractive antibody-binding epitopes and, thereby, confounding vaccine development. To characterize these fluctuations, we examined how exposure of the MPER epitope found at the base of Env is regulated by interactions of the V3 loop located in the apex. Starting with an extremely flexible Env with a readily-exposed MPER, we identified two alterations that substantially restricted antibody access to the epitope. The first, a wholesale swap of V3 loops between HIV-1 strains, energetically stabilized Env in a closed structure that restricted access to antibodies throughout the protein. The second, a point mutation that altered V3-loop charge, specifically destabilized the MPER-exposed conformation but had minimal impact on antibody access to Env regions in between the apex and base. The results indicate that MPER exposure is not explicitly tied to the dynamics of Env regions between the apex and base and suggest a new structural fluctuation during Env activation.

Analysis of sequence diversity and selection pressure in HIV-1 clade C gp41 from India

Evaluation of viral diversity is critical for the rational design of treatment modalities against Human immunodeficiency virus (HIV). Predominated by HIV-1 clade C (HIV-1C), the epidemic in India represents the third largest population infected with HIV-1 globally. Glycoprotein 41 (gp41) is critical for viral replication and is a target for the design of therapeutic strategies. However, documentation of viral diversity of gp41 gene in infected individuals from India remains limited. Present study employed high throughput sequencing to examine variation in gp41 amplicons generated from blood derived viruses in 24 HIV-1C infected individuals from Mumbai, India. Sequence diversity profiles were documented in different functional domains of gp41. Furthermore, through a meta-analysis approach, all reported gp41 sequences from India (N = 70) were compared with those from South Africa (N = 126), country with the largest HIV epidemic globally, also predominated by HIV-1C. A total of 44 positions displayed statistically significant differential (p < 0.05) Shannon entropy in the two regions. This comparison also identified 11 codon sites undergoing distinct selection, 8 of which remained differentially selected in an extended comparison of data from Asia (N = 137) and Africa(N = 383). Assessment of correlated mutation networks associated with differentially selected residues revealed common as well as distinct interaction networks. Furthermore, codon usage analysis revealed 17 differentially selected codons (Mann-Whitney test, p < 0.001) in Asia and Africa. Dissimilar trends in GC content across codon positions were also observed. In depth understanding of these divergent evolutionary signatures through extended analysis with larger data-sets would assist development of effective interventions being considered for HIV-1C.

High affinity CXCR4 inhibitors generated by linking low affinity peptides

G-protein coupled receptors (GPCRs) are implicated in many diseases and attractive targets for drug discovery. Peptide fragments derived from protein ligands of GPCRs are commonly used as probes of GPCR function and as leads for drug development. However, these peptide fragments lack the structural integrity of their parent full-length protein ligands and often show low receptor affinity, which limits their research and therapeutic values. It remains a challenge to efficiently generate high affinity peptide inhibitors of GPCRs. We have investigated a combinational approach involving the synthetic covalent linkage of two low affinity peptide fragments to determine if the strategy can yield high affinity GPCR inhibitors. We examined this design approach using the chemokine receptor CXCR4 as a model of GPCR system. Here, we provide a proof of concept demonstration by designing and synthesizing two peptides, AR5 and AR6, that combine a peptide fragment derived from two viral ligands of CXCR4, vMIP-II and HIV-1 envelope glycoprotein gp120. AR5 and AR6 display nanomolar binding affinity, in contrast to the weak micromolar CXCR4 binding of each peptide fragment alone, and inhibit HIV-1 entry via CXCR4. Further studies were carried out for the representative peptide AR6 using western blotting and site-directed mutagenesis in conjunction with molecular dynamic simulation and binding free energy calculation to determine how the peptide interacts with CXCR4 and inhibits its downstream signaling. These results demonstrate that this combinational approach is effective for generating nanomolar active inhibitors of CXCR4 and may be applicable to other GPCRs.

Structure-based lead optimization to improve antiviral potency and ADMET properties of phenyl-1H-pyrrole-carboxamide entry inhibitors targeted to HIV-1 gp120

We are continuing our concerted effort to optimize our first lead entry antagonist, NBD-11021, which targets the Phe43 cavity of the HIV-1 envelope glycoprotein gp120, to improve antiviral potency and ADMET properties. In this report, we present a structure-based approach that helped us to generate working hypotheses to modify further a recently reported advanced lead entry antagonist, NBD-14107, which showed significant improvement in antiviral potency when tested in a single-cycle assay against a large panel of Env-pseudotyped viruses. We report here the synthesis of twenty-nine new compounds and evaluation of their antiviral activity in a single-cycle and multi-cycle assay to derive a comprehensive structure-activity relationship (SAR). We have selected three inhibitors with the high selectivity index for testing against a large panel of 55 Env-pseudotyped viruses representing a diverse set of clinical isolates of different subtypes. The antiviral activity of one of these potent inhibitors, 55 (NBD-14189), against some clinical isolates was as low as 63 nM. We determined the sensitivity of CD4-binding site mutated-pseudoviruses to these inhibitors to confirm that they target HIV-1 gp120. Furthermore, we assessed their ADMET properties and compared them to the clinical candidate attachment inhibitor, BMS-626529. The ADMET data indicate that some of these new inhibitors have comparable ADMET properties to BMS-626529 and can be optimized further to potential clinical candidates.

Co-evolution techniques are reshaping the way we do structural bioinformatics

Co-evolution techniques were originally conceived to assist in protein structure prediction by inferring pairs of residues that share spatial proximity. However, the functional relationships that can be extrapolated from co-evolution have also proven to be useful in a wide array of structural bioinformatics applications. These techniques are a powerful way to extract structural and functional information in a sequence-rich world.

A comprehensive identification-evidence based alternative for HIV/AIDS treatment with HAART in the healthcare industries

BACKGROUND AND OBJECTIVE

The HIV/AIDS-related issue has given rise to a priority concern in which potential new therapies are increasingly highlighted to lessen the negative impact of highly active anti-retroviral therapy (HAART) in the healthcare industry. With the motivation of "medical applications," this study focuses on the main advanced feature selection techniques and classification approaches that reflect a new architecture, and a trial to build a hybrid model for interested parties.

METHODS

This study first uses an integrated linear-nonlinear feature selection technique to identify the determinants influencing HAART medication and utilizes organizations of different condition-attributes to generate a hybrid model based on a rough set classifier to study evolving HIV/AIDS research in order to improve classification performance.

RESULTS

The proposed model makes use of a real data set from Taiwan's specialist medical center. The experimental results show that the proposed model yields a satisfactory result that is superior to the listed methods, and the core condition-attributes PVL, CD4, Code, Age, Year, PLT, and Sex were identified in the HIV/AIDS data set. In addition, the decision rule set created can be referenced as a knowledge-based healthcare service system as the best of evidence-based practices in the workflow of current clinical diagnosis.

CONCLUSIONS

This study highlights the importance of these key factors and provides the rationale that the proposed model is an effective alternative to analyzing sustained HAART medication in follow-up studies of HIV/AIDS treatment in practice.