Structure-based identification of novel scaffolds as potential HIV-1 entry inhibitors involving CCR5

C-C chemokine receptor 5 (CCR5), which is part of the chemokine receptor family, is a member of the G protein-coupled receptor superfamily. The interactions of CCR5 with HIV-1 during viral entry position it as an effective therapeutic target for designing potent antiviral therapies. The small-molecule Maraviroc was approved by the FDA as a CCR5 drug in 2007, while clinical trials failure has characterised many of the other CCR5 inhibitors. Thus, the continual identification of potential CCR5 inhibitors is, therefore, warranted. In this study, a structure-based discovery approach has been utilised to screen and retrieved novel potential CCR5 inhibitors from the Asinex antiviral compound (∼ 8,722) database. Explicit lipid-bilayer molecular dynamics simulation, in silico physicochemical and pharmacokinetic analyses, were further performed for the top compounds. A total of 23 structurally diverse compounds with binding scores higher than Maraviroc were selected. Subsequent molecular dynamics (MD) simulations analysis of the top four compounds LAS 51495192, BDB 26405401, BDB 26419079, and LAS 34154543, maintained stability at the CCR5 binding site. Furthermore, these compounds made pertinent interactions with CCR5 residues critical for the HIV-1 gp120-V3 loop binding such as Trp86, Tyr89, Phe109, Tyr108, Glu283 and Tyr251. Additionally, the predicted in silico physicochemical and pharmacokinetic descriptors of the selected compounds were within the acceptable range for drug-likeness. The results suggest positive indications that the identified molecules may represent promising CCR5 entry inhibitors. Further structural optimisations and biochemical testing of the proposed compounds may assist in the discovery of effective HIV-1 therapy.Communicated by Ramaswamy H. Sarma.

Mind the Gap—Deciphering GPCR Pharmacology Using 3D Pharmacophores and Artificial Intelligence

G protein-coupled receptors (GPCRs) are amongst the most pharmaceutically relevant and well-studied protein targets, yet unanswered questions in the field leave significant gaps in our understanding of their nuanced structure and function. Three-dimensional pharmacophore models are powerful computational tools in in silico drug discovery, presenting myriad opportunities for the integration of GPCR structural biology and cheminformatics. This review highlights success stories in the application of 3D pharmacophore modeling to de novo drug design, the discovery of biased and allosteric ligands, scaffold hopping, QSAR analysis, hit-to-lead optimization, GPCR de-orphanization, mechanistic understanding of GPCR pharmacology and the elucidation of ligand–receptor interactions. Furthermore, advances in the incorporation of dynamics and machine learning are highlighted. The review will analyze challenges in the field of GPCR drug discovery, detailing how 3D pharmacophore modeling can be used to address them. Finally, we will present opportunities afforded by 3D pharmacophore modeling in the advancement of our understanding and targeting of GPCRs.

Screening and bioavailability evaluation of anti-oxidative selenium-containing peptides from soybeans based on specific structures

Our previous study has evaluated the antioxidant capacity and identified the sequences of soybean selenium-containing peptides. Herein, pharmacophore screening, gastrointestinal simulation and in vivo pharmacokinetics were performed to predict the potentials of selenium-containing peptides in terms of antioxidant activity, safety and bioavailability. A pharmacophore model with 6 structure features was constructed for virtual screening to determine the potential activities of 85 selenium sequences from soybean peptides. Strong reversing effects (p < 0.05) of the targeted sequences were observed in tumor necrosis factor-α (TNF-α)-induced inflammatory cytokines and adhesion factors burst in EA·hy926/Caco-2 co-culture cell models. Ser-Phe-Gln-SeMet (SFQSeM), a promising peptide selected from both virtual screening and cell models, was proved to be stable in the gastrointestinal tract and could be transported across the Caco-2 monolayer via the paracellular pathway. Additionally, SFQSeM showed a long residence time (89.42 ± 1.34 min) and half-life (81.60 ± 11.88 min) after consumption, and it induced lower liver alanine/aspartate transaminase (ALT/AST) and serum nitric oxide (NO) levels compared to Na2SeO3 and SeMet (p < 0.05). The potency of SFQSeM against oxidative stress as well as its oral bioavailability and low risk highlight its potential utility as an effective Se nutritional supplement.

Allosteric mechanism of an oximino-piperidino-piperidine antagonist for the CCR5 chemokine receptor

The first step for the HIV-1 virus infecting host cell is bound with the CCR5 chemokine receptor. A set of allosteric inhibitors of oximino-piperidino-piperidine antagonists for CCR5 chemokine receptor was discovered. However, the allosteric mechanism of these inhibitors is still unsolved. Therefore, residue-level dynamics correlation network combining with on molecular dynamics simulation was used to investigate the allosteric mechanism. The dynamics correlation network of bound CCR5 is significantly different from that of free CCR5. The community of the most active complex suggests that the allosteric information can freely transfer from the allosteric site to the effector site of the second extracellular loop, while the information transfers bottleneck for the less active one. Here, a hypothesis was proposed that "binding-induced allosteric mechanism" was used to reveal the allosteric regulation of antagonists and the network perturbation confirmed it. Finally, the shortest path algorithm was used to identify the possible allosteric pathway with Gly173-Lys171-Thr177-Tyr89-LIG which was evaluated by the network perturbation of key residue. Furthermore, the efficiency of allostery for the most active system is the highest among these antagonist complexes. The strategy targeting the allosteric pathway can be used to design novel inhibitors of HIV-1 virus.

Using accelerated molecular dynamics simulation to shed light on the mechanism of activation/deactivation upon mutations for CCR5

In this work, accelerated molecular dynamics (aMD) simulations were used to study different effects of G286F and R126 mutations on the activity of CCR5. Potential of Mean Force (PMF) results indicate that there are stable inactive-like states and active-like ones existing in the conformation space of the wild type (WT), confirming that CCR5 could possess to some extent constitutive activity. But the R126N mutation could constrain CCR5 in the inactive state through influencing the TXP motif and limiting the movements of TM5 and TM6. In contrast, the G286F mutation promotes the activity of the receptor by increasing the distance of TM2–TM6 and the flexibility of the intracellular part of TM5 and changing the H-bonding in the TXP motif. The observations from the cross correlation analysis further show that the R126N mutation dramatically reduces the motion correlations between TMs, which should partly contribute to the deactivation of CCR5. Compared with the WT system, TM6 and TM7 in the G286F mutant are loosely correlated with other regions, which should be conducive to drive the movement of TM6 and TM7 toward the active conformation. In addition, the result from the protein structure network (PSN) analysis reveals that the shortest pathways connecting the extracellular and the intracellular domains are highly conserved in the three systems despite the different mutations, in which the hydrogen bond plays a pivotal role. However, the G286F mutation shortens the lifetime of the pathway with respect to the R126N mutation, which may be associated with the different activities of the two mutants. The pathway connecting the ligand-binding site and the G-protein region reveals that the allosteric communication between TM6 and TM7 is enhanced by the R126N mutation while the G286F mutation induces the activation of the G-protein pocket by arousing more residues in the NPxxY region to participate in the pathway.

In this work, accelerated molecular dynamics (aMD) simulations were used to study different effects of G286F and R126 mutations on the activity of CCR5.