Molecular recognition of CXCR4 by a dual tropic HIV-1 gp120 V3 loop

Highly Accurate and Efficient Deep Learning Paradigm for Full-Atom Protein Loop Modeling with KarmaLoop

Protein loop modeling is a challenging yet highly nontrivial task in protein structure prediction. Despite recent progress, existing methods including knowledge-based, ab initio, hybrid, and deep learning (DL) methods fall substantially short of either atomic accuracy or computational efficiency. To overcome these limitations, we present KarmaLoop, a novel paradigm that distinguishes itself as the first DL method centered on full-atom (encompassing both backbone and side-chain heavy atoms) protein loop modeling. Our results demonstrate that KarmaLoop considerably outperforms conventional and DL-based methods of loop modeling in terms of both accuracy and efficiency, with the average RMSDs of 1.77 and 1.95 Å for the CASP13+14 and CASP15 benchmark datasets, respectively, and manifests at least 2 orders of magnitude speedup in general compared with other methods. Consequently, our comprehensive evaluations indicate that KarmaLoop provides a state-of-the-art DL solution for protein loop modeling, with the potential to hasten the advancement of protein engineering, antibody-antigen recognition, and drug design.

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

Distinct genetic clusters in HIV-1 CRF01_AE-infected patients induced variable degrees of CD4+ T-cell loss

CRF01_AE strains have been shown to form multiple transmission clusters in China, and some clusters have disparate pathogenicity in Chinese men who have sex with men. This study focused on other CRF01_AE clusters prevalent in heterosexual populations. The CD4+ T-cell counts from both cross-section data in National HIV Molecular Epidemiology Survey and seropositive cohort data were used to evaluate the pathogenicity of the CRF01_AE clusters and other HIV-1 sub-types. Their mechanisms of pathogenicity were evaluated by co-receptor tropisms, predicted by genotyping and confirmed with virus isolate phenotyping, as well as inflammation parameters. Our research elucidated that individuals infected with CRF01_AE clusters 1 and 2 exhibited significantly lower baseline CD4+ T-cell counts and greater CD4+ T-cell loss in cohort follow-up, compared with other HIV-1 sub-types and CRF01_AE clusters. The increased pathogenesis of cluster 1 or 2 was associated with higher CXCR4 tropisms, higher inflammation/immune activation, and increased pyroptosis. The protein structure modeling analysis revealed that the envelope V3 loop of clusters 1 and 2 viruses is favorable for CXCR4 co-receptor usage. Imbedded with the most mutating reverse transcriptase, HIV-1 is one of the most variable viruses. CRF01_AE clusters 1 and 2 have been found to have evolved into more virulent strains in regions with predominant heterosexual infections. The virulent strains increased the pressure for early diagnosis and treatment in HIV patients. To save more lives, HIV-1 surveillance systems should be upgraded from serology and genotyping to phenotyping, which could support precision interventions for those infected by virulent viruses.

IMPORTANCE

Retroviruses swiftly adapt, employing error-prone enzymes for genetic and phenotypic evolution, optimizing survival strategies, and enhancing virulence levels. HIV-1 CRF01_AE has persistently undergone adaptive selection, and cluster 1 and 2 infections display lower counts and fast loss of CD4+ T cells than other HIV-1 sub-types and CRF01_AE clusters. Its mechanisms are associated with increased CXCR4 tropism due to an envelope structure change favoring a tropism shift from CCR5 to CXCR4, thereby shaping viral phenotype features and impacting pathogenicity. This underscores the significance of consistently monitoring HIV-1 genetic evolution and phenotypic transfer to see whether selection bias across risk groups alters the delicate balance of transmissible versus toxic trade-offs, since virulent strains such as CRF01_AE clusters 1 and 2 could seriously compromise the efficacy of antiviral treatment. Only through such early warning and diagnostic services can precise antiviral treatments be administered to those infected with more virulent HIV-1 strains.

CXCR4 Is a Potential Target for Anti-HIV Gene Therapy

The human immunodeficiency virus (HIV) epidemic is a global issue. The estimated number of people with HIV is 39,000,000 to date. Antiviral therapy is the primary approach to treat the infection. However, it does not allow for a complete elimination of the pathogen. The advances in modern gene therapy methods open up new possibilities of effective therapy. One of these areas of possibility is the development of technologies to prevent virus penetration into the cell. Currently, a number of technologies aimed at either the prevention of virus binding to the CCR5 coreceptor or its knockout are undergoing various stages of clinical trials. Since HIV can also utilize the CXCR4 coreceptor, technologies to modify this receptor are also required. Standard knockout of CXCR4 is impossible due to its physiological significance. This review presents an analysis of interactions between individual amino acids in CXCR4 and physiological ligands and HIV gp120. It also discusses potential targets for gene therapy approaches aimed at modifying the coreceptor.

Single-domain antibodies applied as antiviral immunotherapeutics

Viral infections have been the cause of high mortality rates throughout different periods in history. Over the last two decades, outbreaks caused by zoonotic diseases and transmitted by arboviruses have had a significant impact on human health. The emergence of viral infections in different parts of the world encourages the search for new inputs to fight pathologies of viral origin. Antibodies represent the predominant class of new drugs developed in recent years and approved for the treatment of various human diseases, including cancer, autoimmune and infectious diseases. A promising group of antibodies are single-domain antibodies derived from camelid heavy chain immunoglobulins, or VHHs, are biomolecules with nanometric dimensions and unique pharmaceutical and biophysical properties that can be used in the diagnosis and immunotherapy of viral infections. For viral neutralization to occur, VHHs can act in different stages of the viral cycle, including the actual inhibition of infection, to hindering viral replication or assembly. This review article addresses advances involving the use of VHHs in therapeutic propositions aimed to battle different viruses that affect human health.

Spermine and spermidine bind CXCR4 and inhibit CXCR4- but not CCR5-tropic HIV-1 infection

Semen is an important vector for sexual HIV-1 transmission. Although CXCR4-tropic (X4) HIV-1 may be present in semen, almost exclusively CCR5-tropic (R5) HIV-1 causes systemic infection after sexual intercourse. To identify factors that may limit sexual X4-HIV-1 transmission, we generated a seminal fluid-derived compound library and screened it for antiviral agents. We identified four adjacent fractions that blocked X4-HIV-1 but not R5-HIV-1 and found that they all contained spermine and spermidine, abundant polyamines in semen. We showed that spermine, which is present in semen at concentrations up to 14 mM, binds CXCR4 and selectively inhibits cell-free and cell-associated X4-HIV-1 infection of cell lines and primary target cells at micromolar concentrations. Our findings suggest that seminal spermine restricts sexual X4-HIV-1 transmission.

Profiling of proteome changes in plasma of HIV-infected patients receiving antiretroviral therapy

PURPOSE

Antiretroviral therapy (ART) prevents human immunodeficiency virus (HIV)-1 onward transmission and disease progression, leading to excellent prognosis in people living with HIV-1 (PWH). However, side effects, complications, and impaired immune reconstitution persist in some patients treated with ART. We aimed to profile proteome changes in plasma before and after ART to identify the molecular pathways altered by ART.

EXPERIMENTAL DESIGN

Quantitative proteomics analysis based on tandem mass tag (TMT) labeling was used to profile proteome changes of paired plasma samples from HIV-1 patients before receiving ART and after ART treatment.

RESULTS

A total of 1398 protein groups (PGs) were identified, in which 18 proteins were downregulated and 50 were upregulated in plasma from ART treated patients. Based on Ingenuity Pathway analysis (IPA), gap junction signaling and actin cytoskeleton signaling were enriched among upregulated proteins, while downregulated proteins were mainly participated in IL-15 signaling pathway. Patients with the low level of CSF1R and the high levels of MINPP1 and TGM3 showed better CD4⁺ T-cell recovery.

CONCLUSIONS

The present study provided plasma proteome changes after ART to elucidate the underlying mechanistic pathways in response to ART, and also identified potential targets to prompt immune reconstitution.

Unique profile of predominant CCR5-tropic in CRF07_BC HIV-1 infections and discovery of an unusual CXCR4-tropic strain

CRF07_BC is one of the most prevalent HIV-1 strains in China, which contributes over one-third of the virus transmissions in the country. In general, CRF07_BC is associated with slower disease progression, while the underlying mechanisms remain unclear. Our study focused on envelope proteins (Env) and its V3 loop which determine viral binding to co-receptors during infection of cells. We studied a large dataset of 3,937 env sequences in China and found that CRF07_BC had a unique profile of predominantly single CCR5 tropism compared with CCR5 and CXCR4 dual tropisms in other HIV-1 subtypes. The percentages of the CXCR4-tropic virus in B (3.7%) and CRF01_AE (10.4%) infection are much higher than that of CRF07_BC (0.1%), which is supported by median false-positive rates (FPRs) of 69.8%, 25.5%, and 13.4% for CRF07_BC, B, and CRF01_AE respectively, with a cutoff FPR for CXCR4-tropic at 2%. In this study, we identified the first pure CXCR4-tropic virus from one CRF07_BC-infected patient with an extremely low CD4⁺T cell count (7 cells/mm³). Structural analysis found that the V3 region of this virus has the characteristic 7T and 25R and a substitution of conserved “GPGQ” crown motif for “GPGH”. This study provided compelling evidence that CRF07_BC has the ability to evolve into CXCR4 strains. Our study also lay down the groundwork for studies on tropism switch, which were commonly done for other HIV-1 subtypes, for the long-delayed CRF07_BC.

Immune reconstruction effectiveness of combination antiretroviral therapy for HIV-1 CRF01_AE cluster 1 and 2 infected individuals

There are great disparities of the results in immune reconstruction (IR) of the HIV-1 infected patients during combined antiretroviral therapy (cART), due to both host polymorphisms and viral genetic subtypes. Identifying these factors and elucidating their impact on the IR could help to improve the efficacy. To study the factors influencing the IR, we conducted a 15-year retrospective cohort study of HIV-1 infected individuals under cART. The trend of CD4⁺ count changes was evaluated by the generalized estimating equations. Cox proportional model and propensity score matching were used to identify variables that affect the possibility of achieving IR. The tropism characteristics of virus were compared using the coreceptor binding model. In addition to baseline CD4⁺ counts and age implications, CRF01_AE cluster 1 was associated with a poorer probability of achieving IR than infection with cluster 2 (aHR, 1.39; 95%CI, 1.02-1.90) and other subtypes (aHR, 1.83; 95%CI, 1.31-2.56). The mean time from cART initiation to achieve IR was much longer in patients infected by CRF01_AE cluster 1 than other subtypes/sub-clusters (P < 0.001). In-depth analysis indicated that a higher proportion of CXCR4 viruses were found in CRF01_AE clusters 1 and 2 (P < 0.05), and showed tendency to favour CXCR4 binding to V3 signatures. This study indicated the immune restoration impairment found in patients were associated with HIV-1 CRF01_AE cluster 1, which was attributed to the high proportion of CXCR4-tropic viruses. To improve the effectiveness of cART, more efforts should be made in the early identification of HIV-1 subtype/sub-cluster and monitoring of virus phenotypes.

Computational Modulation of the V3 Region of Glycoprotein gp125 of HIV-2

HIV-2 infection is frequently neglected in HIV/AIDS campaigns. However, a special emphasis must be given to HIV-2 as an untreated infection that also leads to AIDS and death, and for which the efficacy of most available drugs is limited against HIV-2. HIV envelope glycoproteins mediate binding to the receptor CD4 and co-receptors at the surface of the target cell, enabling fusion with the cell membrane and viral entry. Here, we developed and optimized a computer-assisted drug design approach of an important HIV-2 glycoprotein that allows us to explore and gain further insights at the molecular level into protein structures and interactions crucial for the inhibition of HIV-2 cell entry. The 3D structure of a key HIV-2ROD gp125 region was generated by a homology modeling campaign. To disclose the importance of the main structural features and compare them with experimental results, 3D-models of six mutants were also generated. These mutations revealed the selective impact on the behavior of the protein. Furthermore, molecular dynamics simulations were performed to optimize the models, and the dynamic behavior was tackled to account for structure flexibility and interactions network formation. Structurally, the mutations studied lead to a loss of aromatic features, which is very important for the establishment of π-π interactions and could induce a structural preference by a specific coreceptor. These new insights into the structure-function relationship of HIV-2 gp125 V3 and surrounding regions will help in the design of better models and the design of new small molecules capable to inhibit the attachment and binding of HIV with host cells.

Identification of novel molecular determinants of co-receptor usage in HIV-1 subtype F V3 envelope sequences

HIV-1 determinants of coreceptor usage within the gp120 V3 loop have been broadly studied over the past years. This information has led to the development of state-of the-art bioinformatic tools that are useful to predict co-receptor usage based on the V3 loop sequence mainly of subtypes B, C and A. However, these methods show a poor performance for subtype F V3 loops, which are found in an increasing number of HIV-1 strains worldwide. In the present work we investigated determinants of viral tropisms in the understudied subtype F by looking at genotypic and structural information of coreceptor:V3 loop interactions in a novel group of 40 subtype F V3 loops obtained from HIV-1 strains phenotypically characterized either as syncytium inducing or non-syncytium inducing by the MT-2 assay. We provide novel information about estimated interactions energies between a set of V3 loops with known tropism in subtype F, that allowed us to improve predictions of the coreceptor usage for this subtype. Understanding genetic and structural features underlying HIV coreceptor usage across different subtypes is relevant for the rational design of preventive and therapeutic strategies aimed at limiting the HIV-1 epidemic worldwide.

Effects of sunlight exposure and vitamin D supplementation on HIV patients

Unlike many vitamins derived predominantly from food sources, vitamin D is produced endogenously in the skin upon exposure to sunlight. Ethnicity, skin pigmentation, socioeconomic status, geographic location, climate and sunscreen; all of these factors contribute to the amount of insolation for any given individual. Insufficient insolation creates the prerequisites for vitamin D deficiency. This is particularly true in HIV-infected individuals, who are highly vulnerable to vitamin D insufficiency/deficiency, as it plays a huge role in the musculoskeletal and cardiovascular systems. Antiretroviral therapy may also be a factor in vitamin D deficiency. Today, as the issues of preventing common skeletal and non-skeletal diseases with HIV-infected people are becoming highly relevant, the maintenance of vitamin D levels through exposure to sunlight or supplementation appears to be an effective and safe solution. This review focuses on studies concerning the potential role of vitamin D supplementation through adequate sunlight exposure or dietary intake in HIV-infected people. The biology and epidemiology of HIV infection, as well as the issues related to vitamin D deficiency, its status on immune function, the effect of vitamin D against HIV disease progression and other health aspects of this vitamin, are briefly explained.

Computational evolution of an RNA-binding protein towards enhanced oxidized-RNA binding

The oxidation of RNA has been implicated in the development of many diseases. Among the four ribonucleotides, guanosine is the most susceptible to oxidation, resulting in the formation of 8-oxo-7,8-dihydroguanosine (8-oxoG). Despite the limited knowledge about how cells regulate the detrimental effects of oxidized RNA, cellular factors involved in its regulation have begun to be identified. One of these factors is polynucleotide phosphorylase (PNPase), a multifunctional enzyme implicated in RNA turnover. In the present study, we have examined the interaction of PNPase with 8-oxoG in atomic detail to provide insights into the mechanism of 8-oxoG discrimination. We hypothesized that PNPase subunits cooperate to form a binding site using the dynamic SFF loop within the central channel of the PNPase homotrimer. We evolved this site using a novel approach that initially screened mutants from a library of beneficial mutations and assessed their interactions using multi-nanosecond Molecular Dynamics simulations. We found that evolving this single site resulted in a fold change increase in 8-oxoG affinity between 1.2 and 1.5 and/or selectivity between 1.5 and 1.9. In addition to the improvement in 8-oxoG binding, complementation of K12 Δpnp with plasmids expressing mutant PNPases caused increased cell tolerance to H₂O₂. This observation provides a clear link between molecular discrimination of RNA oxidation and cell survival. Moreover, this study provides a framework for the manipulation of modified-RNA protein readers, which has potential application in synthetic biology and epitranscriptomics.

Interactions between Curcumin Derivatives and Amyloid-β Fibrils: Insights from Molecular Dynamics Simulations

The aggregation of amyloid-β (Aβ) peptides into senile plaques is a hallmark of Alzheimer's disease (AD) and is hypothesized to be the primary cause of AD related neurodegeneration. Previous studies have shown the ability of curcumin to both inhibit the aggregation of Aβ peptides into oligomers or fibrils and reduce amyloids in vivo. Despite the promise of curcumin and its derivatives to serve as diagnostic, preventative, and potentially therapeutic AD molecules, the mechanism by which curcumin and its derivatives bind to and inhibit Aβ fibrils' formation remains elusive. Here, we investigated curcumin and a set of curcumin derivatives in complex with a hexamer peptide model of the Aβ1-42 fibril using nearly exhaustive docking, followed by multi-ns molecular dynamics simulations, to provide atomistic-detail insights into the molecules' binding and inhibitory properties. In the vast majority of the simulations, curcumin and its derivatives remain firmly bound in complex with the fibril through primarily three different principle binding modes, in which the molecules interact with residue domain 17LVFFA21, in line with previous experiments. In a small subset of these simulations, the molecules partly dissociate the outermost peptide of the Aβ1-42 fibril by disrupting β-sheets within the residue domain 12VHHQKLVFF20. A comparison between binding modes leading or not leading to partial dissociation of the outermost peptide suggests that the latter is attributed to a few subtle key structural and energetic interaction-based differences. Interestingly, partial dissociation appears to be either an outcome of high affinity interactions or a cause leading to high affinity interactions between the molecules and the fibril, which could partly serve as a compensation for the energy loss in the fibril due to partial dissociation. In conjunction with this, we suggest a potential inhibition mechanism of Αβ1-42 aggregation by the molecules, where the partially dissociated 16KLVFF20 domain of the outermost peptide could either remain unstructured or wrap around to form intramolecular interactions with the same peptide's 29GAIIG33 domain, while the molecules could additionally act as a patch against the external edge of the second outermost peptide's 16KLVFF20 domain. Thereby, individually or concurrently, these could prohibit fibril elongation.

Designer Amyloid Cell-Penetrating Peptides for Potential Use as Gene Transfer Vehicles

Cell-penetrating peptides are used extensively to deliver molecules into cells due to their unique characteristics such as rapid internalization, charge, and non-cytotoxicity. Amyloid fibril biomaterials were reported as gene transfer or retroviral infection enhancers; no cell internalization of the peptides themselves is reported so far. In this study, we focus on two rationally and computationally designed peptides comprised of β-sheet cores derived from naturally occurring protein sequences and designed positively charged and aromatic residues exposed at key residue positions. The β-sheet cores bestow the designed peptides with the ability to self-assemble into amyloid fibrils. The introduction of positively charged and aromatic residues additionally promotes DNA condensation and cell internalization by the self-assembled material formed by the designed peptides. Our results demonstrate that these designer peptide fibrils can efficiently enter mammalian cells while carrying packaged luciferase-encoding plasmid DNA, and they can act as a protein expression enhancer. Interestingly, the peptides additionally exhibited strong antimicrobial activity against the enterobacterium Escherichia coli.

Molecular Mechanism for Attractant Signaling to DHMA by E. coli Tsr

The attractant chemotaxis response of Escherichia coli to norepinephrine requires that it be converted to 3,4-dihydroxymandelic acid (DHMA) by the monoamine oxidase TynA and the aromatic aldehyde dehydrogenase FeaB. DHMA is sensed by the serine chemoreceptor Tsr, and the attractant response requires that at least one subunit of the periplasmic domain of the Tsr homodimer (pTsr) has an intact serine-binding site. DHMA that is generated in vivo by E. coli is expected to be a racemic mixture of the (R) and (S) enantiomers, so it has been unclear whether one or both chiral forms are active. Here, we used a combination of state-of-the-art tools in molecular docking and simulations, including an in-house simulation-based docking protocol, to investigate the binding properties of (R)-DHMA and (S)-DHMA to E. coli pTsr. Our studies computationally predicted that (R)-DHMA should promote a stronger attractant response than (S)-DHMA because of a consistently greater-magnitude piston-like pushdown of the pTsr α-helix 4 toward the membrane upon binding of (R)-DHMA than upon binding of (S)-DHMA. This displacement is caused primarily by interaction of DHMA with Tsr residue Thr156, which has been shown by genetic studies to be critical for the attractant response to L-serine and DHMA. These findings led us to separate the two chiral species and test their effectiveness as chemoattractants. Both the tethered cell and motility migration coefficient assays validated the prediction that (R)-DHMA is a stronger attractant than (S)-DHMA. Our study demonstrates that refined computational docking and simulation studies combined with experiments can be used to investigate situations in which subtle differences between ligands may lead to diverse chemotactic responses.

Isoflavones as Ah Receptor Agonists in Colon-Derived Cell Lines: Structure-Activity Relationships

Many of the protective responses observed for flavonoids in the gastrointestinal track resemble aryl hydrocarbon receptor (AhR)-mediated effects. Therefore, we examined the structure-activity relationships of isoflavones and isomeric flavone and flavanones as AhR ligands on the basis of their induction of CYP1A1, CYP1B1, and UGT1A1 gene expression in colon cancer Caco2 cells and young adult mouse colonocyte (YAMC) cells. Caco2 cells were significantly more Ah-responsive than YAMC cells, and this was due, in part, to flavonoid-induced cytotoxicity in the latter cell lines. The structure-activity relationships for the flavonoids were complex and both response and cell context specific; however, there was significant variability in the AhR activities of the isomeric substituted isoflavones and flavones. For example, 4',5,7-trihydroxyisoflavone (genistein) was AhR-inactive whereas 4',5,7-trihydroxyflavone (apigenin) induced CYP1A1, CYP1B1, and UGT1A1 in Caco2 cells. In contrast, both 5,7-dihydroxy-4-methoxy substituted isoflavone (biochanin A) and flavone (acacetin) induced all three AhR-responsive genes; 4',5,7-trimethoxyisoflavone was a potent AhR agonist, and the isomeric flavone was AhR-inactive. These results coupled with simulation studies modeling flavonoid interaction within the AhR binding pocket demonstrate that the orientation of the substituted phenyl ring at C-2 (flavones) or C-3 (isoflavones) on the common 4-H-chromen-4-one ring strongly influences the activities of isoflavones and flavones as AhR agonists.

Molecular Mechanism of HIV-1 Entry

HIV-1 envelope glycoprotein [Env; trimeric (gp160)₃ cleaved to (gp120/gp41)₃] attaches the virion to a susceptible cell and induces fusion of viral and cell membranes to initiate infection. It interacts with the primary receptor CD4 and coreceptor (e.g., chemokine receptor CCR5 or CXCR4) to allow viral entry by triggering large structural rearrangements and unleashing the fusogenic potential of gp41 to induce membrane fusion. Recent advances in structural biology of HIV-1 Env and its complexes with the cellular receptors have revealed molecular details of HIV-1 entry and yielded new mechanistic insights. In this review, I summarize our latest understanding of the HIV-1 membrane fusion process and discuss possible pathways for productive viral entry.

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.

Modeling interaction between gp120 HIV protein and CCR5 receptor

The study of the process of HIV entry into the host cell and the creation of biomimetic nanosystems that are able to selectively bind viral particles and proteins is a high priority research area for the development of novel diagnostic tools and treatment of HIV infection. Recently, we described multilayer nanoparticles (nanotraps) with heparin surface and cationic peptides comprising the N-terminal tail (Nt) and the second extracellular loop (ECL2) of CCR5 receptor, which could bind with high affinity some inflammatory chemokines, in particular, Rantes. Because of the similarity of the binding determinants in CCR5 structure, both for chemokines and gp120 HIV protein, here we expand this approach to the study of the interactions of these biomimetic nanosystems and their components with the peptide representing the V3 loop of the activated form of gp120. According to surface plasmon resonance results, a conformational rearrangement is involved in the process of V3 and CCR5 fragments binding. As in the case of Rantes, ECL2 peptide showed much higher affinity to V3 peptide than Nt (K_D  = 3.72 × 10^-8 and 1.10 × 10^-6  M, respectively). Heparin-covered nanoparticles bearing CCR5 peptides effectively bound V3 as well. The presence of both heparin and the peptides in the structure of the nanotraps was shown to be crucial for the interaction with the V3 loop. Thus, short cationic peptides ECL2 and Nt proved to be excellent candidates for the design of CCR5 receptor mimetics.

Disparate impact on CD4 T cell count by two distinct HIV-1 phylogenetic clusters from the same clade

HIV-1 evolved into various genetic subtypes and circulating recombinant forms (CRFs) in the global epidemic. The same subtype or CRF is usually considered to have similar phenotype. Being one of the world's major CRFs, CRF01_AE infection was reported to associate with higher prevalence of CXCR4 (X4) viruses and faster CD4 decline. However, the underlying mechanisms remain unclear. We identified eight phylogenetic clusters of CRF01_AE in China and hypothesized that they may have different phenotypes. In the National HIV Molecular Epidemiology Survey, we discovered that people infected by CRF01_AE cluster 4 had significantly lower CD4 counts (391 vs. 470, P < 0.0001) and higher prevalence of X4-using viruses (17.1% vs. 4.4%, P < 0.0001) compared with those infected by cluster 5. In an MSM cohort, X4-using viruses were only isolated from seroconvertors in cluster 4, which was associated with low a CD4 count within the first year of infection (141 vs. 440, P = 0.003). Using a coreceptor binding model, we identified unique V3 signatures in cluster 4 that favor CXCR4 use. We demonstrate that the HIV-1 phenotype and pathogenicity can be determined at the phylogenetic cluster level in the same subtype. Since its initial spread to humans from chimpanzees, estimated to be the first half of the 20th century, HIV-1 continues to undergo rapid evolution in larger and more diverse populations. The divergent phenotype evolution of two major CRF01_AE clusters highlights the importance of monitoring the genetic evolution and phenotypic shift of HIV-1 to provide early warning of the appearance of more pathogenic strains.

Structural basis of coreceptor recognition by HIV-1 envelope spike

HIV-1 envelope glycoprotein (Env), which consists of trimeric (gp160)₃ cleaved to (gp120 and gp41)₃, interacts with the primary receptor CD4 and a coreceptor (such as chemokine receptor CCR5) to fuse viral and target-cell membranes. The gp120-coreceptor interaction has previously been proposed as the most crucial trigger for unleashing the fusogenic potential of gp41. Here we report a cryo-electron microscopy structure of a full-length gp120 in complex with soluble CD4 and unmodified human CCR5, at 3.9 Å resolution. The V3 loop of gp120 inserts into the chemokine-binding pocket formed by seven transmembrane helices of CCR5, and the N terminus of CCR5 contacts the CD4-induced bridging sheet of gp120. CCR5 induces no obvious allosteric changes in gp120 that can propagate to gp41; it does bring the Env trimer close to the target membrane. The N terminus of gp120, which is gripped by gp41 in the pre-fusion or CD4-bound Env, flips back in the CCR5-bound conformation and may irreversibly destabilize gp41 to initiate fusion. The coreceptor probably functions by stabilizing and anchoring the CD4-induced conformation of Env near the cell membrane. These results advance our understanding of HIV-1 entry into host cells and may guide the development of vaccines and therapeutic agents.

Elucidating the multi-targeted anti-amyloid activity and enhanced islet amyloid polypeptide binding of β-wrapins

β-wrapins are engineered binding proteins stabilizing the β-hairpin conformations of amyloidogenic proteins islet amyloid polypeptide (IAPP), amyloid-β, and α-synuclein, thus inhibiting their amyloid propensity. Here, we use computational and experimental methods to investigate the molecular recognition of IAPP by β-wrapins. We show that the multi-targeted, IAPP, amyloid-β, and α-synuclein, binding properties of β-wrapins originate mainly from optimized interactions between β-wrapin residues and sets of residues in the three amyloidogenic proteins with similar physicochemical properties. Our results suggest that IAPP is a comparatively promiscuous β-wrapin target, probably due to the low number of charged residues in the IAPP β-hairpin motif. The sub-micromolar affinity of β-wrapin HI18, specifically selected against IAPP, is achieved in part by salt-bridge formation between HI18 residue Glu10 and the IAPP N-terminal residue Lys1, both located in the flexible N-termini of the interacting proteins. Our findings provide insights towards developing novel protein-based single- or multi-targeted therapeutics.

Virtual Screening of Chemical Compounds for Discovery of Complement C3 Ligands

The complement system is our first line of defense against foreign pathogens, but when it is not properly regulated, complement is implicated in the pathology of several autoimmune and inflammatory disorders. Compstatin is a peptidic complement inhibitor that acts by blocking the cleavage of complement protein C3 to the proinflammatory fragment C3a and opsonin fragment C3b. In this study, we aim to identify druglike small-molecule complement inhibitors with physicochemical, geometric, and binding properties similar to those of compstatin. We employed two approaches using various high-throughput virtual screening methods, which incorporate molecular dynamics (MD) simulations, pharmacophore model design, energy calculations, and molecular docking and scoring. We have generated a library of 274 chemical compounds with computationally predicted binding affinities for the compstatin binding site of C3. We have tested subsets of these chemical compounds experimentally for complement inhibitory activity, using hemolytic assays, and for binding affinity, using microscale thermophoresis. As a result, although none of the compounds showed inhibitory activity, compound 29 was identified to exhibit weak competitive binding against a potent compstatin analogue, therefore validating our computational approaches. Additional docking and MD simulation studies suggest that compound 29 interacts with C3 residues, which have been shown to be important in binding of compstatin to the C3c fragment of C3. Compound 29 is amenable to physicochemical optimization to acquire inhibitory properties. Additionally, it is possible that some of the untested compounds will demonstrate binding and inhibition in future experimental studies.

A novel amyloid designable scaffold and potential inhibitor inspired by GAIIG of amyloid beta and the HIV-1 V3 loop

The GAIIG sequence, common to the amyloid beta peptide (residues 29-33) and to the HIV-1 gp120 (residues 24-28 in a typical V3 loop), self-assembles into amyloid fibrils, as suggested by theory and the experiments presented here. The longer YATGAIIGNII sequence from the V3 loop also self-assembles into amyloid fibrils, of which the first three and the last two residues are outside the amyloid GAIIG core. We postulate that this sequence, with suitably selected modifications at the flexible positions, can serve as a designable scaffold for novel amyloid-based materials. Moreover, we report the single crystal X-ray structure of the beta-breaker peptide GAIPIG at 1.05 Å resolution. The structural information provided in this study could serve as the basis for structure-based design of potential inhibitors of amyloid formation.

Molecular Modeling of Chemoreceptor:Ligand Interactions

Docking algorithms have been widely used to elucidate ligand:receptor interactions that are important in biological function. Here, we introduce an in-house developed docking-refinement protocol that combines the following innovative features. (1) The use of multiple short molecular dynamics (MD) docking simulations, with residues within the binding pocket of the receptor unconstrained, so that the binding modes of the ligand in the binding pocket may be exhaustively examined. (2) The initial positioning of the ligand within the binding pocket based on complementary shape, and the use of both harmonic and quartic spherical potentials to constrain the ligand in the binding pocket during multiple short docking simulations. (3) The selection of the most probable binding modes generated by the short docking simulations using interaction energy calculations, as well as the subsequent application of all-atom MD simulations and physical-chemistry based free energy calculations to elucidate the most favorable binding mode of the ligand in complex with the receptor. In this chapter, we provide step-by-step instructions on how to computationally investigate the binding of small-molecule ligands to protein receptors by examining as control and test cases, respectively, the binding of L-serine and R-3,4-dihydroxymandelic acid (R-DHMA) to the Escherichia coli chemoreceptor Tsr. Similar computational strategies can be used for the molecular modeling of a series of ligand:protein receptor interactions.

Protein structural disorder of the envelope V3 loop contributes to the switch in human immunodeficiency virus type 1 cell tropism

Human immunodeficiency virus type 1 (HIV-1) envelope gp120 is partly an intrinsically disordered (unstructured/disordered) protein as it contains regions that do not fold into well-defined protein structures. These disordered regions play important roles in HIV’s life cycle, particularly, V3 loop-dependent cell entry, which determines how the virus uses two coreceptors on immune cells, the chemokine receptors CCR5 (R5), CXCR4 (X4) or both (R5X4 virus). Most infecting HIV-1 variants utilise CCR5, while a switch to CXCR4-use occurs in the majority of infections. Why does this ‘rewiring’ event occur in HIV-1 infected patients? As changes in the charge of the V3 loop are associated with this receptor switch and it has been suggested that charged residues promote structure disorder, we hypothesise that the intrinsic disorder of the V3 loop is permissive to sequence variation thus contributing to the switch in cell tropism. To test this we use three independent data sets of gp120 to analyse V3 loop disorder. We find that the V3 loop of X4 virus has significantly higher intrinsic disorder tendency than R5 and R5X4 virus, while R5X4 virus has the lowest. These results indicate that structural disorder plays an important role in HIV-1 cell tropism and CXCR4 binding. We discuss the potential evolutionary mechanisms leading to the fixation of disorder promoting mutations and the adaptive potential of protein structural disorder in viral host adaptation.

HIV-1 subtype CRF01_AE and B differ in utilization of low levels of CCR5, Maraviroc susceptibility and potential N-glycosylation sites

HIV subtypes not only predominate in different geographical regions but also differ in key phenotypic characteristics. To determine if genotypic and/or phenotypic differences in the Envelope (Env) glycoprotein can explain subtype related differences, we cloned 37 full length Envs from Subtype B and AE HIV infected individuals from Singapore. Our data demonstrates that CRF01_AE Envs have lower Potential N Glycosylation Sites and higher risk of ×4 development. Phenotypically, CRF01_AE were less infectious than subtype B Envs in cells expressing low levels of CCR5. Moreover, the Maraviroc IC₅₀ was higher for subtype B Envs and correlated with infectivity in low CCR5 expressing cells as well as PNGS. Specifically, the glycosylation site N301 in the V3 loop was seen less frequently in AE subtype and CXCR4 topic viruses. CRF01_AE differs from B subtype in terms of CCR5 usage and Maraviroc susceptibility which may have implications for HIV pathogenesis and virus evolution.

Structural Basis for Apelin Control of the Human Apelin Receptor

Apelin receptor (APJR) is a key regulator of human cardiovascular function and is activated by two different endogenous peptide ligands, apelin and Elabela, each with different isoforms diversified by length and amino acid sequence. Here we report the 2.6-Å resolution crystal structure of human APJR in complex with a designed 17-amino-acid apelin mimetic peptide agonist. The structure reveals that the peptide agonist adopts a lactam constrained curved two-site ligand binding mode. Combined with mutation analysis and molecular dynamics simulations with apelin-13 binding to the wild-type APJR, this structure provides a mechanistic understanding of apelin recognition and binding specificity. Comparison of this structure with that of other peptide receptors suggests that endogenous peptide ligands with a high degree of conformational flexibility may bind and modulate the receptors via a similar two-site binding mechanism.

Position 22 of the V3 loop is associated with HIV infectivity

Human immunodeficiency virus subtype 1B (HIV-1B) binds to the CD4 receptor and co-receptor CCR5 or CXCR4 to enter T lymphocytes. The amino acid sequence of the HIV envelope glycoprotein V3 region determines the co-receptor tropism, thereby influencing the infectivity of the virus. Our research group previously found that the amino acid at position 22 of the V3 region may affect the infectivity of the virus, and in this study, we tested this hypothesis. We constructed pseudoviruses by changing the amino acids at position 22 of the V3 region in CCR5-tropic and CXCR4-tropic viruses and tested their infectivity. When the amino acid at V3 position 22 was altered in the CCR5- and CXCR4-tropic viruses, their ability to infect cells decreased to 20.6% and 17.14%, respectively. Therefore, we propose that residue 22 in the V3 region of subtype HIV-1B significantly influences the infectivity of the virus.

Conformational properties of the third variable loop of HIV-1AD8 envelope glycoprotein in the liganded conditions

The conformational dynamics of the HIV-1 envelope glycoprotein gp120 and gp41 (Env) remains poorly understood. Here we examined how the V3 loop conformation is regulated in the liganded state using a panel of recombinant HIV-1NL4-3 clones bearing HIV-1AD8 Env by two experimental approaches, one adopting a monoclonal neutralizing antibody KD-247 (suvizumab) that recognizes the tip of the V3 loop, and the other assessing the function of the V3 loop. A significant positive correlation of the Env-KD-247 binding was detected between the liganded and unliganded conditions. Namely, the mutation D163G located in the V2 loop, which enhances viral susceptibility to KD-247 by 59.4-fold, had little effect on the sCD4-induced increment of the virus-KD-247 binding. By contrast, a virus with the S370N mutation in the C3 region increased the virus-KD-247 binding by 91.4-fold, although it did not influence the KD-247-mediated neutralization. Co-receptor usage and the susceptibility to CCR5 inhibitor Maraviroc were unaffected by D163G and S370N mutations. Collectively, these data suggest that the conformation of the liganded V3-loop of HIV-1AD8 Env is still under regulation of other Env domains aside from the V3 loop, including V2 and C3. Our results give an insight into the structural properties of HIV-1 Env and viral resistance to entry inhibitors by non-V3 loop mutations.

Highly Accurate Structure-Based Prediction of HIV-1 Coreceptor Usage Suggests Intermolecular Interactions Driving Tropism

HIV-1 entry into host cells is mediated by interactions between the V3-loop of viral glycoprotein gp120 and chemokine receptor CCR5 or CXCR4, collectively known as HIV-1 coreceptors. Accurate genotypic prediction of coreceptor usage is of significant clinical interest and determination of the factors driving tropism has been the focus of extensive study. We have developed a method based on nonlinear support vector machines to elucidate the interacting residue pairs driving coreceptor usage and provide highly accurate coreceptor usage predictions. Our models utilize centroid-centroid interaction energies from computationally derived structures of the V3-loop:coreceptor complexes as primary features, while additional features based on established rules regarding V3-loop sequences are also investigated. We tested our method on 2455 V3-loop sequences of various lengths and subtypes, and produce a median area under the receiver operator curve of 0.977 based on 500 runs of 10-fold cross validation. Our study is the first to elucidate a small set of specific interacting residue pairs between the V3-loop and coreceptors capable of predicting coreceptor usage with high accuracy across major HIV-1 subtypes. The developed method has been implemented as a web tool named CRUSH, CoReceptor USage prediction for HIV-1, which is available at http://ares.tamu.edu/CRUSH/.

Use of network model to explore dynamic and allosteric properties of three GPCR homodimers

We used an elastic network model and protein structure network to study three class A GPCR homodimers.

Progress toward rationally designed small-molecule peptide and peptidomimetic CXCR4 antagonists

Over the last 5 years, X-ray structures of CXCR4 in complex with three different ligands (the small-molecule antagonist IT1t, the polypeptide antagonist CVX15 and the viral chemokine antagonist vMIP-II) have been released. In addition to the inherent scientific value of these specific X-ray structures, they provide a reliable structural foundation for studies of the molecular interactions between CXCR4 and its key peptide ligands (CXCL12 and HIV-1 gp120), and serve as valuable templates for further development of small-molecule CXCR4 antagonists with therapeutic potential. We here review recent computational studies of the molecular interactions between CXCR4 and its peptide ligands – based on the X-ray structures of CXCR4 – and the current status of small-molecule peptide and peptidomimetic CXCR4 antagonists.

Chemokine and chemokine receptor structure and interactions: implications for therapeutic strategies

The control of cell migration by chemokines involves interactions with two types of receptors: seven transmembrane chemokine-type G protein-coupled receptors and cell surface or extracellular matrix-associated glycosaminoglycans. Coordinated interaction of chemokines with both types of receptors is required for directional migration of cells in numerous physiological and pathological processes. Accumulated structural information, culminating most recently in the structure of a chemokine receptor in complex with a chemokine, has led to a view where chemokine oligomers bind to glycosaminoglycans through epitopes formed when chemokine subunits come together, while chemokine monomers bind to receptors in a pseudo two-step mechanism of receptor activation. Exploitation of this structural knowledge has and will continue to provide important information for therapeutic strategies, as described in this review.

Function-oriented development of CXCR4 antagonists as selective human immunodeficiency virus (HIV)-1 entry inhibitors

Motivated by the pivotal role of CXCR4 as an HIV entry co-receptor, we herein report a de novo hit-to-lead effort on the identification of subnanomolar purine-based CXCR4 antagonists against HIV-1 infection. Compound 24, with an EC50 of 0.5 nM against HIV-1 entry into host cells and an IC50 of 16.4 nM for inhibition of radioligand stromal-derived factor-1α (SDF-1α) binding to CXCR4, was also found to be highly selective against closely related chemokine receptors. We rationalized that compound 24 complementarily interacted with the critical CXCR4 residues that are essential for binding to HIV-1 gp120 V3 loop and subsequent viral entry. Compound 24 showed a 130-fold increase in anti-HIV activity compared to that of the marketed CXCR4 antagonist, AMD3100 (Plerixafor), whereas both compounds exhibited similar potency in mobilization of CXCR4(+)/CD34(+) stem cells at a high dose. Our study offers insight into the design of anti-HIV therapeutics devoid of major interference with SDF-1α function.

The impact of HIV-1 on neurogenesis: implications for HAND

HIV-1 infection, in addition to its destructive effects on the immune system, plays a role in the development of neurocognitive deficits. Indeed up to 50% of long-term HIV infected patients suffer from HIV-associated neurocognitive disorders (HAND). These deficits have been well characterized and defined clinically according to a number of cognitive parameters. HAND is often accompanied by atrophy of the brain including inhibition of neurogenesis, especially in the hippocampus. Many mechanisms have been proposed as contributing factors to HAND including induction of oxidative stress in the central nervous system (CNS), chronic microglial-mediated neuroinflammation, amyloid-beta (Aβ) deposition, hyperphosphorylated tau protein, and toxic effects of combination antiretroviral therapy (cART). In these review we focus solely on recent experimental evidence suggesting that disturbance by HIV-1 results in impairment of neurogenesis as one contributing factor to HAND. Impaired neurogenesis has been linked to cognitive deficits and other neurodegenerative disorders. This article will highlight recently identified pathological mechanisms which potentially contribute to the development of impaired neurogenesis by HIV-1 or HIV-1-associated proteins from both animal and human studies.

Insights into the mechanism of C5aR inhibition by PMX53 via implicit solvent molecular dynamics simulations and docking

Background

The complement protein C5a acts by primarily binding and activating the G-protein coupled C5a receptor C5aR (CD88), and is implicated in many inflammatory diseases. The cyclic hexapeptide PMX53 (sequence Ace-Phe-[Orn-Pro-dCha-Trp-Arg]) is a full C5aR antagonist of nanomolar potency, and is widely used to study C5aR function in disease.

Results

We construct for the first time molecular models for the C5aR:PMX53 complex without the a priori use of experimental constraints, via a computational framework of molecular dynamics (MD) simulations, docking, conformational clustering and free energy filtering. The models agree with experimental data, and are used to propose important intermolecular interactions contributing to binding, and to develop a hypothesis for the mechanism of PMX53 antagonism.

Conclusion

This work forms the basis for the design of improved C5aR antagonists, as well as for atomic-detail mechanistic studies of complement activation and function. Our computational framework can be widely used to develop GPCR-ligand structural models in membrane environments, peptidomimetics and other chemical compounds with potential clinical use.

Elucidating a key anti-HIV-1 and cancer-associated axis: the structure of CCL5 (Rantes) in complex with CCR5

CCL5 (RANTES) is an inflammatory chemokine which binds to chemokine receptor CCR5 and induces signaling. The CCL5:CCR5 associated chemotactic signaling is of critical biological importance and is a potential HIV-1 therapeutic axis. Several studies provided growing evidence for the expression of CCL5 and CCR5 in non-hematological malignancies. Therefore, the delineation of the CCL5:CCR5 complex structure can pave the way for novel CCR5-targeted drugs. We employed a computational protocol which is primarily based on free energy calculations and molecular dynamics simulations, and report, what is to our knowledge, the first computationally derived CCL5:CCR5 complex structure which is in excellent agreement with experimental findings and clarifies the functional role of CCL5 and CCR5 residues which are associated with binding and signaling. A wealth of polar and non-polar interactions contributes to the tight CCL5:CCR5 binding. The structure of an HIV-1 gp120 V3 loop in complex with CCR5 has recently been derived through a similar computational protocol. A comparison between the CCL5 : CCR5 and the HIV-1 gp120 V3 loop : CCR5 complex structures depicts that both the chemokine and the virus primarily interact with the same CCR5 residues. The present work provides insights into the blocking mechanism of HIV-1 by CCL5.

Practical Aspects of Free-Energy Calculations: A Review

Free-energy calculations in the framework of classical molecular dynamics simulations are nowadays used in a wide range of research areas including solvation thermodynamics, molecular recognition, and protein folding. The basic components of a free-energy calculation, that is, a suitable model Hamiltonian, a sampling protocol, and an estimator for the free energy, are independent of the specific application. However, the attention that one has to pay to these components depends considerably on the specific application. Here, we review six different areas of application and discuss the relative importance of the three main components to provide the reader with an organigram and to make nonexperts aware of the many pitfalls present in free energy calculations.

Molecular recognition of CCR5 by an HIV-1 gp120 V3 loop

The binding of protein HIV-1 gp120 to coreceptors CCR5 or CXCR4 is a key step of the HIV-1 entry to the host cell, and is predominantly mediated through the V3 loop fragment of HIV-1 gp120. In the present work, we delineate the molecular recognition of chemokine receptor CCR5 by a dual tropic HIV-1 gp120 V3 loop, using a comprehensive set of computational tools predominantly based on molecular dynamics simulations and free energy calculations. We report, what is to our knowledge, the first complete HIV-1 gp120 V3 loop : CCR5 complex structure, which includes the whole V3 loop and the N-terminus of CCR5, and exhibits exceptional agreement with previous experimental findings. The computationally derived structure sheds light into the functional role of HIV-1 gp120 V3 loop and CCR5 residues associated with the HIV-1 coreceptor activity, and provides insights into the HIV-1 coreceptor selectivity and the blocking mechanism of HIV-1 gp120 by maraviroc. By comparing the binding of the specific dual tropic HIV-1 gp120 V3 loop with CCR5 and CXCR4, we observe that the HIV-1 gp120 V3 loop residues 13-21, which include the tip, share nearly identical structural and energetic properties in complex with both coreceptors. This result paves the way for the design of dual CCR5/CXCR4 targeted peptides as novel potential anti-AIDS therapeutics.

Elucidating a key component of cancer metastasis: CXCL12 (SDF-1α) binding to CXCR4

The chemotactic signaling induced by the binding of chemokine CXCL12 (SDF-1α) to chemokine receptor CXCR4 is of significant biological importance and is a potential therapeutic axis against HIV-1. However, as CXCR4 is overexpressed in certain cancer cells, the CXCL12:CXCR4 signaling is involved in tumor metastasis, progression, angiogenesis, and survival. Motivated by the pivotal role of the CXCL12:CXCR4 axis in cancer, we employed a comprehensive set of computational tools, predominantly based on free energy calculations and molecular dynamics simulations, to obtain insights into the molecular recognition of CXCR4 by CXCL12. We report, what is to our knowledge, the first computationally derived CXCL12:CXCR4 complex structure which is in remarkable agreement with experimental findings and sheds light into the functional role of CXCL12 and CXCR4 residues which are associated with binding and signaling. Our results reveal that the CXCL12 N-terminal domain is firmly bound within the CXCR4 transmembrane domain, and the central 24–50 residue domain of CXCL12 interacts with the upper N-terminal domain of CXCR4. The stability of the CXCL12:CXCR4 complex structure is attributed to an abundance of nonpolar and polar intermolecular interactions, including salt bridges formed between positively charged CXCL12 residues and negatively charged CXCR4 residues. The success of the computational protocol can mainly be attributed to the nearly exhaustive docking conformational search, as well as the heterogeneous dielectric implicit water-membrane-water model used to simulate and select the optimum conformations. We also recently utilized this protocol to elucidate the binding of an HIV-1 gp120 V3 loop in complex with CXCR4, and a comparison between the molecular recognition of CXCR4 by CXCL12 and the HIV-1 gp120 V3 loop shows that both CXCL12 and the HIV-1 gp120 V3 loop share the same CXCR4 binding pocket, as they mostly interact with the same CXCR4 residues.