UDock2: interactive real-time multi-body protein-protein docking software.

MOTIVATION: Protein-protein docking aims at predicting the geometry of protein interactions to gain insights into the mechanisms underlying these processes and develop new strategies for drug discovery. Interactive and user-oriented manipulation tools can support this task complementary to automated software. RESULTS: This article presents an interactive multi-body protein-protein docking software, UDock2, designed for research but also usable for teaching and popularization of science purposes due to its high usability. In UDock2, the users tackle the conformational space of protein interfaces using an intuitive real-time docking procedure with on-the-fly scoring. UDock2 integrates traditional computer graphics methods to facilitate the visualization and to provide better insight into protein surfaces, interfaces, and properties. AVAILABILITY AND IMPLEMENTATION: UDock2 is open-source, cross-platform (Windows and Linux), and available at http://udock.fr. The code can be accessed at https://gitlab.com/Udock/Udock2.

Comparative evaluation of shape retrieval methods on macromolecular surfaces: an application of computer vision methods in structural bioinformatics.

MOTIVATION: The investigation of the structure of biological systems at the molecular level gives insights about their functions and dynamics. Shape and surface of biomolecules are fundamental to molecular recognition events. Characterizing their geometry can lead to more adequate predictions of their interactions. In the present work, we assess the performance of reference shape retrieval methods from the computer vision community on protein shapes. RESULTS: Shape retrieval methods are efficient in identifying orthologous proteins and tracking large conformational changes. This work illustrates the interest for the protein surface shape as a higher-level representation of the protein structure that (i) abstracts the underlying protein sequence, structure or fold, (ii) allows the use of shape retrieval methods to screen large databases of protein structures to identify surficial homologs and possible interacting partners and (iii) opens an extension of the protein structure-function paradigm toward a protein structure-surface(s)-function paradigm. AVAILABILITYAND IMPLEMENTATION: All data are available online at http://datasetmachat.drugdesign.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Predicting Potential Endocrine Disrupting Chemicals Binding to Estrogen Receptor α (ERα) Using a Pipeline Combining Structure-Based and Ligand-Based in Silico Methods.

The estrogen receptors α (ERα) are transcription factors involved in several physiological processes belonging to the nuclear receptors (NRs) protein family. Besides the endogenous ligands, several other chemicals are able to bind to those receptors. Among them are endocrine disrupting chemicals (EDCs) that can trigger toxicological pathways. Many studies have focused on predicting EDCs based on their ability to bind NRs; mainly, estrogen receptors (ER), thyroid hormones receptors (TR), androgen receptors (AR), glucocorticoid receptors (GR), and peroxisome proliferator-activated receptors gamma (PPARγ). In this work, we suggest a pipeline designed for the prediction of ERα binding activity. The flagged compounds can be further explored using experimental techniques to assess their potential to be EDCs. The pipeline is a combination of structure based (docking and pharmacophore models) and ligand based (pharmacophore models) methods. The models have been constructed using the Environmental Protection Agency (EPA) data encompassing a large number of structurally diverse compounds. A validation step was then achieved using two external databases: the NR-DBIND (Nuclear Receptors DataBase Including Negative Data) and the EADB (Estrogenic Activity DataBase). Different combination protocols were explored. Results showed that the combination of models performed better than each model taken individually. The consensus protocol that reached values of 0.81 and 0.54 for sensitivity and specificity, respectively, was the best suited for our toxicological study. Insights and recommendations were drawn to alleviate the screening quality of other projects focusing on ERα binding predictions.

Predicting the affinity of Farnesoid X Receptor ligands through a hierarchical ranking protocol: a D3R Grand Challenge 2 case study.

The Drug Design Data Resource (D3R) Grand Challenges are blind contests organized to assess the state-of-the-art methods accuracy in predicting binding modes and relative binding free energies of experimentally validated ligands for a given target. The second stage of the D3R Grand Challenge 2 (GC2) was focused on ranking 102 compounds according to their predicted affinity for Farnesoid X Receptor. In this task, our workflow was ranked 5th out of the 77 submissions in the structure-based category. Our strategy consisted in (1) a combination of molecular docking using AutoDock 4.2 and manual edition of available structures for binding poses generation using SeeSAR, (2) the use of HYDE scoring for pose selection, and (3) a hierarchical ranking using HYDE and MM/GBSA. In this report, we detail our pose generation and ligands ranking protocols and provide guidelines to be used in a prospective computer aided drug design program.

Screening Explorer-An Interactive Tool for the Analysis of Screening Results.

Screening Explorer is a web-based application that allows for an intuitive evaluation of the results of screening experiments using complementary metrics in the field. The usual evaluation of screening results implies the separate generation and apprehension of the ROC, predictiveness, and enrichment curves and their global metrics. Similarly, partial metrics need to be calculated repeatedly for different fractions of a data set and there exists no handy tool that allows reading partial metrics simultaneously on different charts. For a deeper understanding of the results of screening experiments, we rendered their analysis straightforward by linking these metrics interactively in an interactive usable web-based application. We also implemented simple consensus scoring methods based on scores normalization, standardization (z-scores), and compounds ranking to evaluate the enrichments that can be expected through methods combination. Two demonstration data sets allow the users to easily apprehend the functions of this tool that can be applied to the analysis of virtual and experimental screening results. Screening Explorer is freely accessible at http://stats.drugdesign.fr .

Benchmarking Data Sets for the Evaluation of Virtual Ligand Screening Methods: Review and Perspectives.

Virtual screening methods are commonly used nowadays in drug discovery processes. However, to ensure their reliability, they have to be carefully evaluated. The evaluation of these methods is often realized in a retrospective way, notably by studying the enrichment of benchmarking data sets. To this purpose, numerous benchmarking data sets were developed over the years, and the resulting improvements led to the availability of high quality benchmarking data sets. However, some points still have to be considered in the selection of the active compounds, decoys, and protein structures to obtain optimal benchmarking data sets.

Evidence after imputation for a role of MICA variants in nonprogression and elite control of HIV type 1 infection.

Past genome-wide association studies (GWAS) involving individuals with AIDS have mainly identified associations in the HLA region. Using the latest software, we imputed 7 million single-nucleotide polymorphisms (SNPs)/indels of the 1000 Genomes Project from the GWAS-determined genotypes of individuals in the Genomics of Resistance to Immunodeficiency Virus AIDS nonprogression cohort and compared them with those of control cohorts. The strongest signals were in MICA, the gene encoding major histocompatibility class I polypeptide-related sequence A (P = 3.31 × 10(-12)), with a particular exonic deletion (P = 1.59 × 10(-8)) in full linkage disequilibrium with the reference HCP5 rs2395029 SNP. Haplotype analysis also revealed an additive effect between HLA-C, HLA-B, and MICA variants. These data suggest a role for MICA in progression and elite control of human immunodeficiency virus type 1 infection.

Synthesis and structure-activity relationship of non-peptidic antagonists of neuropilin-1 receptor.

Neuropilins (NRPs) are VEGF-A165 co-receptors over-expressed in tumor cells, and considered as targets in angiogenic-related pathologies. We previously identified compound 1, the first non-peptidic antagonist of the VEGF-A165/NRP binding, which exhibits in vivo anti-angiogenic and anti-tumor activities. We report here the synthesis and biological evaluations of new antagonists structurally-related to compound 1. Among these molecules, 4a, 4c and 4d show cytotoxic effects on HUVEC and MDA-MB-31 cells, and antagonize VEGF-A165/NRP-1 binding. This study confirmed our key structure-activity relationships hypothesis and paved the way to compound 1 'hit to lead' optimization.

Importance of the pharmacological profile of the bound ligand in enrichment on nuclear receptors: toward the use of experimentally validated decoy ligands.

The evaluation of virtual ligand screening methods is of major importance to ensure their reliability. Taking into account the agonist/antagonist pharmacological profile should improve the quality of the benchmarking data sets since ligand binding can induce conformational changes in the nuclear receptor structure and such changes may vary according to the agonist/antagonist ligand profile. We indeed found that splitting the agonist and antagonist ligands into two separate data sets for a given nuclear receptor target significantly enhances the quality of the evaluation. The pharmacological profile of the ligand bound in the binding site of the target structure was also found to be an additional critical parameter. We also illustrate that active compound data sets for a given pharmacological activity can be used as a set of experimentally validated decoy ligands for another pharmacological activity to ensure a reliable and challenging evaluation of virtual screening methods.

Efficient GPU computation of large protein Solvent-Excluded Surface.

The Solvent-Excluded Surface (SES) is an essential representation of molecules which is massively used in molecular modeling and drug discovery since it represents the interacting surface between molecules. Based on its properties, it supports the visualization of both large scale shapes and details of molecules. While several methods targeted its computation, the ability to process large molecular structures to address the introduction of big complex analysis while leveraging the massively parallel architecture of GPUs has remained a challenge. This is mostly caused by the need for consequent memory allocation or by the complexity of the parallelization of its processing. In this paper, we leverage the last theoretical advances made for the depiction of the SES to provide fast analytical computation with low impact on memory. We show that our method is able to compute the complete surface while handling large molecular complexes with competitive computation time costs compared to previous works.

PLO3S: Protein LOcal Surficial Similarity Screening.

The study of protein molecular surfaces enables to better understand and predict protein interactions. Different methods have been developed in computer vision to compare surfaces that can be applied to protein molecular surfaces. The present work proposes a method using the Wave Kernel Signature: Protein LOcal Surficial Similarity Screening (PLO3S). The descriptor of the PLO3S method is a local surface shape descriptor projected on a unit sphere mapped onto a 2D plane and called Surface Wave Interpolated Maps (SWIM). PLO3S allows to rapidly compare protein surface shapes through local comparisons to filter large protein surfaces datasets in protein structures virtual screening protocols.

Multiple structures for virtual ligand screening: defining binding site properties-based criteria to optimize the selection of the query.

Structure based virtual ligand screening (SBVLS) methods are widely used in drug discovery programs. When several structures of the target are available, protocols based either on single structure docking or on ensemble docking can be used. The performance of the methods depends on the structure(s) used as a reference, whose choice requires retrospective enrichment studies on benchmarking databases which consume additional resources. In the present study, we have identified several trends in the properties of the binding sites of the structures that led to the optimal performance in retrospective SBVLS tests whatever the docking program used (Surflex-dock or ICM). By assessing their hydrophobicity and comparing their volume and opening, we show that the selection of optimal structures should be possible with no requirement of prior retrospective enrichment studies. If the mean binding site volume is lower than 350 A(3), the structure with the smaller volume should be preferred. In the other cases, the structure with the largest binding site should be preferred. These optimal structures may be either selected for a single structure docking strategy or an ensemble docking strategy. When constructing an ensemble, the opening of the site might be an interesting criterion additionaly to its volume as the most closed structures should not be preferred in the large systems. These "binding site properties-based" guidelines could be helpful to optimize future prospective drug discovery protocols when several structures of the target are available.

Multicohort genomewide association study reveals a new signal of protection against HIV-1 acquisition.

BACKGROUND: To date, only mutations in CCR5 have been shown to confer resistance to human immunodeficiency virus type 1 (HIV-1) infection, and these explain only a small fraction of the observed variability in HIV susceptibility. METHODS: We performed a meta-analysis between 2 independent European genomewide association studies, each comparing HIV-1 seropositive cases with normal population controls known to be HIV uninfected, to identify single-nucleotide polymorphisms (SNPs) associated with the HIV-1 acquisition phenotype. SNPs exhibiting P < 10(-5) in this first stage underwent second-stage analysis in 2 independent US cohorts of European descent. RESULTS: After the first stage, a single highly significant association was revealed for the chromosome 8 rs6996198 with HIV-1 acquisition and was replicated in both second-stage cohorts. Across the 4 groups, the rs6996198-T allele was consistently associated with a significant reduced risk of HIV-1 infection, and the global meta-analysis reached genomewide significance: P(combined) = 7.76 × 10(-8). CONCLUSIONS: We provide strong evidence of association for a common variant with HIV-1 acquisition in populations of European ancestry. This protective signal against HIV-1 infection is the first identified outside the CCR5 nexus. First clues point to a potential functional role for a nearby candidate gene, CYP7B1, but this locus warrants further investigation.

Review of in silico studies dedicated to the nuclear receptor family: Therapeutic prospects and toxicological concerns.

Being in the center of both therapeutic and toxicological concerns, NRs are widely studied for drug discovery application but also to unravel the potential toxicity of environmental compounds such as pesticides, cosmetics or additives. High throughput screening campaigns (HTS) are largely used to detect compounds able to interact with this protein family for both therapeutic and toxicological purposes. These methods lead to a large amount of data requiring the use of computational approaches for a robust and correct analysis and interpretation. The output data can be used to build predictive models to forecast the behavior of new chemicals based on their in vitro activities. This atrticle is a review of the studies published in the last decade and dedicated to NR ligands in silico prediction for both therapeutic and toxicological purposes. Over 100 articles concerning 14 NR subfamilies were carefully read and analyzed in order to retrieve the most commonly used computational methods to develop predictive models, to retrieve the databases deployed in the model building process and to pinpoint some of the limitations they faced.

Surface-based protein domains retrieval methods from a SHREC2021 challenge.

Proteins are essential to nearly all cellular mechanism and the effectors of the cells activities. As such, they often interact through their surface with other proteins or other cellular ligands such as ions or organic molecules. The evolution generates plenty of different proteins, with unique abilities, but also proteins with related functions hence similar 3D surface properties (shape, physico-chemical properties, …). The protein surfaces are therefore of primary importance for their activity. In the present work, we assess the ability of different methods to detect such similarities based on the geometry of the protein surfaces (described as 3D meshes), using either their shape only, or their shape and the electrostatic potential (a biologically relevant property of proteins surface). Five different groups participated in this contest using the shape-only dataset, and one group extended its pre-existing method to handle the electrostatic potential. Our comparative study reveals both the ability of the methods to detect related proteins and their difficulties to distinguish between highly related proteins. Our study allows also to analyze the putative influence of electrostatic information in addition to the one of protein shapes alone. Finally, the discussion permits to expose the results with respect to ones obtained in the previous contests for the extended method. The source codes of each presented method have been made available online.

Computationally driven discovery of SARS-CoV-2 Mpro inhibitors: from design to experimental validation.

We report a fast-track computationally driven discovery of new SARS-CoV-2 main protease (Mpro) inhibitors whose potency ranges from mM for the initial non-covalent ligands to sub-µM for the final covalent compound (IC50 = 830 ± 50 nM). The project extensively relied on high-resolution all-atom molecular dynamics simulations and absolute binding free energy calculations performed using the polarizable AMOEBA force field. The study is complemented by extensive adaptive sampling simulations that are used to rationalize the different ligand binding poses through the explicit reconstruction of the ligand-protein conformation space. Machine learning predictions are also performed to predict selected compound properties. While simulations extensively use high performance computing to strongly reduce the time-to-solution, they were systematically coupled to nuclear magnetic resonance experiments to drive synthesis and for in vitro characterization of compounds. Such a study highlights the power of in silico strategies that rely on structure-based approaches for drug design and allows the protein conformational multiplicity problem to be addressed. The proposed fluorinated tetrahydroquinolines open routes for further optimization of Mpro inhibitors towards low nM affinities.

Multiple-cohort genetic association study reveals CXCR6 as a new chemokine receptor involved in long-term nonprogression to AIDS.

BACKGROUND: The compilation of previous genomewide association studies of AIDS shows a major polymorphism in the HCP5 gene associated with both control of the viral load and long-term nonprogression (LTNP) to AIDS. METHODS: To look for genetic variants that affect LTNP without necessary control of the viral load, we reanalyzed the genomewide data of the unique LTNP Genomics of Resistance to Immunodeficiency Virus (GRIV) cohort by excluding "elite controller" patients, who were controlling the viral load at very low levels (<100 copies/mL). RESULTS: The rs2234358 polymorphism in the CXCR6 gene was the strongest signal (P=2.5 x 10(-7); odds ratio, 1.85) obtained for the genomewide association study comparing the 186 GRIV LTNPs who were not elite controllers with 697 uninfected control subjects. This association was replicated in 3 additional independent European studies, reaching genomewide significance of P(combined)=9.7 x 10(-10). This association with LTNP is independent of the CCR2-CCR5 locus and the HCP5 polymorphisms. CONCLUSIONS: The statistical significance, the replication, and the magnitude of the association demonstrate that CXCR6 is likely involved in the molecular etiology of AIDS and, in particular, in LTNP, emphasizing the power of extreme-phenotype cohorts. CXCR6 is a chemokine receptor that is known as a minor coreceptor in human immunodeficiency virus type 1 infection but could participate in disease progression through its role as a mediator of inflammation.

Neuropilin 1 and Neuropilin 2 gene invalidation or pharmacological inhibition reveals their relevance for the treatment of metastatic renal cell carcinoma.

BACKGROUND: Despite the improvement of relapse-free survival mediated by anti-angiogenic drugs like sunitinib (Sutent®), or by combinations of anti-angiogenic drugs with immunotherapy, metastatic clear cell Renal Cell Carcinoma (mccRCC) remain incurable. Hence, new relevant treatments are urgently needed. The VEGFs coreceptors, Neuropilins 1, 2 (NRP1, 2) are expressed on several tumor cells including ccRCC. We analyzed the role of the VEGFs/NRPs signaling in ccRCC aggressiveness and evaluated the relevance to target this pathway. METHODS: We correlated the NRP1, 2 levels to patients' survival using online available data base. Human and mouse ccRCC cells were knocked-out for the NRP1 and NRP2 genes by a CRISPR/Cas9 method. The number of metabolically active cells was evaluated by XTT assays. Migration ability was determined by wound closure experiments and invasion ability by using Boyden chamber coated with collagen. Production of VEGFA and VEGFC was evaluated by ELISA. Experimental ccRCC were generated in immuno-competent/deficient mice. The effects of a competitive inhibitor of NRP1, 2, NRPa-308, was tested in vitro and in vivo with the above-mentioned tests and on experimental ccRCC. NRPa-308 docking was performed on both NRPs. RESULTS: Knock-out of the NRP1 and NRP2 genes inhibited cell metabolism and migration and stimulated the expression of VEGFA or VEGFC, respectively. NRPa-308 presented a higher affinity for NRP2 than for NRP1. It decreased cell metabolism and migration/invasion more efficiently than sunitinib and the commercially available NRP inhibitor EG00229. NRPa-308 presented a robust inhibition of experimental ccRCC growth in immunocompetent and immunodeficient mice. Such inhibition was associated with decreased expression of several pro-tumoral factors. Analysis of the TCGA database showed that the NRP2 pathway, more than the NRP1 pathway correlates with tumor aggressiveness only in metastatic patients. CONCLUSIONS: Our study strongly suggests that inhibiting NRPs is a relevant treatment for mccRCC patients in therapeutic impasses and NRPa-308 represents a relevant hit.

High-resolution mining of the SARS-CoV-2 main protease conformational space: supercomputer-driven unsupervised adaptive sampling.

We provide an unsupervised adaptive sampling strategy capable of producing µs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs). The global exploration problem is decomposed into a set of separate MD trajectories that can be restarted within a selective process to achieve sufficient phase-space sampling. Accurate statistical properties can be obtained through reweighting. Within this highly parallel setup, the Tinker-HP package can be powered by an arbitrary large number of GPUs on supercomputers, reducing exploration time from years to days. This approach is used to tackle the urgent modeling problem of the SARS-CoV-2 Main Protease (Mpro) producing more than 38 µs of all-atom simulations of its apo (ligand-free) dimer using the high-resolution AMOEBA PFF. The first 15.14 µs simulation (physiological pH) is compared to available non-PFF long-timescale simulation data. A detailed clustering analysis exhibits striking differences between FFs, with AMOEBA showing a richer conformational space. Focusing on key structural markers related to the oxyanion hole stability, we observe an asymmetry between protomers. One of them appears less structured resembling the experimentally inactive monomer for which a 6 µs simulation was performed as a basis for comparison. Results highlight the plasticity of the Mpro active site. The C-terminal end of its less structured protomer is shown to oscillate between several states, being able to interact with the other protomer, potentially modulating its activity. Active and distal site volumes are found to be larger in the most active protomer within our AMOEBA simulations compared to non-PFFs as additional cryptic pockets are uncovered. A second 17 µs AMOEBA simulation is performed with protonated His172 residues mimicking lower pH. Data show the protonation impact on the destructuring of the oxyanion loop. We finally analyze the solvation patterns around key histidine residues. The confined AMOEBA polarizable water molecules are able to explore a wide range of dipole moments, going beyond bulk values, leading to a water molecule count consistent with experimental data. Results suggest that the use of PFFs could be critical in drug discovery to accurately model the complexity of the molecular interactions structuring Mpro.

Comparative evaluation of 3D virtual ligand screening methods: impact of the molecular alignment on enrichment.

In the early stage of drug discovery programs, when the structure of a complex involving a target and a small molecule is available, structure-based virtual ligand screening methods are generally preferred. However, ligand-based strategies like shape-similarity search methods can also be applied. Shape-similarity search methods consist in exploring a pseudo-binding-site derived from the known small molecule used as a reference. Several of these methods use conformational sampling algorithms which are also shared by corresponding docking methods: for example Surflex-dock/Surflex-sim, FlexX/FlexS, ICM, and OMEGA-FRED/OMEGA-ROCS. Using 11 systems issued from the challenging "own" subsets of the Directory of Useful Decoys (DUD-own), we evaluated and compared the performance of the above-cited programs in terms of molecular alignment accuracy, enrichment in active compounds, and enrichment in different chemotypes (scaffold-hopping). Since molecular alignment is a crucial aspect of performance for the different methods, we have assessed its impact on enrichment. We have also illustrated the paradox of retrieving active compounds with good scores even if they are inaccurately positioned. Finally, we have highlighted possible positive aspects of using shape-based approaches in drug-discovery protocols when the structure of the target in complex with a small molecule is known.