Benchmarking AlphaFold-Generated Structures of Chemokine-Chemokine Receptor Complexes.

AlphaFold and AlphaFold-Multimer have become two essential tools for the modeling of unknown structures of proteins and protein complexes. In this work, we extensively benchmarked the quality of chemokine-chemokine receptor structures generated by AlphaFold-Multimer against experimentally determined structures. Our analysis considered both the global quality of the model, as well as key structural features for chemokine recognition. To study the effects of template and multiple sequence alignment parameters on the results, a new prediction pipeline called LIT-AlphaFold (https://github.com/LIT-CCM-lab/LIT-AlphaFold) was developed, allowing extensive input customization. AlphaFold-Multimer correctly predicted differences in chemokine binding orientation and accurately reproduced the unique binding orientation of the CXCL12-ACKR3 complex. Further, the predictions of the full receptor N-terminus provided insights into a putative chemokine recognition site 0.5. The accuracy of chemokine N-terminus binding mode prediction varied between complexes, but the confidence score permitted the distinguishing of residues that were very likely well positioned. Finally, we generated a high-confidence model of the unsolved CXCL12-CXCR4 complex, which agreed with experimental mutagenesis and cross-linking data.

Comparing transmembrane protein structures with ATOLL.

SUMMARY: The 3D structure of transmembrane helices plays a key role in the function of membrane proteins. While visual inspection can usually discern the distinctive features of a helix bundle, simply translating them into a 2D diagram can be difficult. ATOLL (Aligned Transmembrane dOmains Layout fLattening) projects the helix bundle onto the lipid bilayer plane, thereby facilitating the comparison of different structures of the same membrane protein or structures of different membrane proteins. AVAILABILITY AND IMPLEMENTATION: ATOLL is a program written in Python3. The source code is freely available on the web at https://github.com/LIT-CCM-lab/ATOLL. ATOLL is implemented into a web server (https://atoll.drugdesign.unistra.fr/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Inhibition of the NAD salvage pathway in schistosomes impairs metabolism, reproduction, and parasite survival.

NAD, a key co-enzyme required for cell metabolism, is synthesized via two pathways in most organisms. Since schistosomes apparently lack enzymes required for de novo NAD biosynthesis, we evaluated whether these parasites, which infect >200 million people worldwide, maintain NAD homeostasis via the NAD salvage biosynthetic pathway. We found that intracellular NAD levels decline in schistosomes treated with drugs that block production of nicotinamide or nicotinamide mononucleotide-known NAD precursors in the non-deamidating salvage pathway. Moreover, in vitro inhibition of the NAD salvage pathway in schistosomes impaired egg production, disrupted the outer membranes of both immature and mature parasites and caused loss of mobility and death. Inhibiting the NAD salvage pathway in schistosome-infected mice significantly decreased NAD levels in adult parasites, which correlated with reduced egg production, fewer liver granulomas and parasite death. Thus, schistosomes, unlike their mammalian hosts, appear limited to one metabolic pathway to maintain NAD-dependent metabolic processes.

Unsupervised Classification of G-Protein Coupled Receptors and Their Conformational States Using IChem Intramolecular Interaction Patterns.

Over the past decade, the ever-growing structural information on G-protein coupled receptors (GPCRs) has revealed the three-dimensional (3D) characteristics of a receptor structure that is competent for G-protein binding. Structural markers are now commonly used to distinguish GPCR functional states, especially when analyzing molecular dynamics simulations. In particular, the position of the sixth helix within the seven transmembrane domains (TMs) is directly related to the coupling of the G-protein. Here, we show that the structural pattern defined by transmembrane intramolecular interactions (hydrogen bonds excluding backbone/backbone interactions, ionic bonds and aromatic interactions) is suitable for comparison of GPCR 3D structures and unsupervised distinction of the receptor states. First, we analyze a microsecond long molecular dynamic simulation of the human ß2-adrenergic receptor (ADRB2). Clustering of the 3D structures by pattern similarity identifies stable states which match the conformational classes defined by structural markers. Furthermore, the method directly spots the few state-specific interactions. Transforming pattern into graph, we extend the method to the comparison of different GPCRs. Clustering all GPCR experimentally determined structures by clique relative size first separates receptors, then their conformational states, thereby suggesting that the interaction patterns are specific of the receptor sequence and that the interaction signatures of conformational states are not shared across distant homologues.

Local Interaction Density (LID), a Fast and Efficient Tool to Prioritize Docking Poses.

Ligand docking at a protein site can be improved by prioritizing poses by similarity to validated binding modes found in the crystal structures of ligand/protein complexes. The interactions formed in the predicted model are searched in each of the reference 3D structures, taken individually. We propose to merge the information provided by all references, creating a single representation of all known binding modes. The method is called LID, an acronym for Local Interaction Density. LID was benchmarked in a pose prediction exercise on 19 proteins and 1382 ligands using PLANTS as docking software. It was also tested in a virtual screening challenge on eight proteins, with a dataset of 140,000 compounds from DUD-E and PubChem. LID significantly improved the performance of the docking program in both pose prediction and virtual screening. The gain is comparable to that obtained with a rescoring approach based on the individual comparison of reference binding modes (the GRIM method). Importantly, LID is effective with a small number of references. LID calculation time is negligible compared to the docking time.

Structural Searching of Biosynthetic Enzymes to Predict Protein Targets of Natural Products.

Recently, we have demonstrated that site comparison methodology using flavonoid biosynthetic enzymes as the query could automatically identify structural features common to different flavonoid-binding proteins, allowing for the identification of flavonoid targets such as protein kinases. With the aim of further validating the hypothesis that biosynthetic enzymes and therapeutic targets can contain a similar natural product imprint, we collected a set of 159 crystallographic structures representing 38 natural product biosynthetic enzymes by searching the Protein Databank. Each enzyme structure was used as a query to screen a repository of approximately 10 000 ligandable sites by active site similarity. We report a full analysis of the screening results and highlight three retrospective examples where the natural product validates the method, thereby revealing novel structural relationships between natural product biosynthetic enzymes and putative protein targets of the natural product. From a prospective perspective, our work provides a list of up to 64 potential novel targets for 25 well-characterized natural products.

Do Fragments and Crystallization Additives Bind Similarly to Drug-like Ligands?

The success of fragment-based drug design (FBDD) hinges upon the optimization of low-molecular-weight compounds (MW < 300 Da) with weak binding affinities to lead compounds with high affinity and selectivity. Usually, structural information from fragment-protein complexes is used to develop ideas about the binding mode of similar but drug-like molecules. In this regard, crystallization additives such as cryoprotectants or buffer components, which are highly abundant in crystal structures, are frequently ignored. Thus, the aim of this study was to investigate the information present in protein complexes with fragments as well as those with additives and how they relate to the binding modes of their drug-like counterparts. We present a thorough analysis of the binding modes of crystallographic additives, fragments, and drug-like ligands bound to four diverse targets of wide interest in drug discovery and highly represented in the Protein Data Bank: cyclin-dependent kinase 2, ß-secretase 1, carbonic anhydrase 2, and trypsin. We identified a total of 630 unique molecules bound to the catalytic binding sites, among them 31 additives, 222 fragments, and 377 drug-like ligands. In general, we observed that, independent of the target, protein-fragment interaction patterns are highly similar to those of drug-like ligands and mostly cover the residues crucial for binding. Crystallographic additives are also able to show conserved binding modes and recover the residues important for binding in some of the cases. Moreover, we show evidence that the information from fragments and drug-like ligands can be applied to rescore docking poses in order to improve the prediction of binding modes.

sc-PDB: a 3D-database of ligandable binding sites--10 years on.

The sc-PDB database (available at http://bioinfo-pharma.u-strasbg.fr/scPDB/) is a comprehensive and up-to-date selection of ligandable binding sites of the Protein Data Bank. Sites are defined from complexes between a protein and a pharmacological ligand. The database provides the all-atom description of the protein, its ligand, their binding site and their binding mode. Currently, the sc-PDB archive registers 9283 binding sites from 3678 unique proteins and 5608 unique ligands. The sc-PDB database was publicly launched in 2004 with the aim of providing structure files suitable for computational approaches to drug design, such as docking. During the last 10 years we have improved and standardized the processes for (i) identifying binding sites, (ii) correcting structures, (iii) annotating protein function and ligand properties and (iv) characterizing their binding mode. This paper presents the latest enhancements in the database, specifically pertaining to the representation of molecular interaction and to the similarity between ligand/protein binding patterns. The new website puts emphasis in pictorial analysis of data.

Probing the catalytic mechanism of bovine CD38/NAD+ glycohydrolase by site directed mutagenesis of key active site residues.

Bovine CD38/NAD(+) glycohydrolase catalyzes the hydrolysis of NAD(+) to nicotinamide and ADP-ribose and the formation of cyclic ADP-ribose via a stepwise reaction mechanism. Our recent crystallographic study of its Michaelis complex and covalently-trapped intermediates provided insights into the modalities of substrate binding and the molecular mechanism of bCD38. The aim of the present work was to determine the precise role of key conserved active site residues (Trp118, Glu138, Asp147, Trp181 and Glu218) by focusing mainly on the cleavage of the nicotinamide-ribosyl bond. We analyzed the kinetic parameters of mutants of these residues which reside within the bCD38 subdomain in the vicinity of the scissile bond of bound NAD(+). To address the reaction mechanism we also performed chemical rescue experiments with neutral (methanol) and ionic (azide, formate) nucleophiles. The crucial role of Glu218, which orients the substrate for cleavage by interacting with the N-ribosyl 2'-OH group of NAD(+), was highlighted. This contribution to catalysis accounts for almost half of the reaction energy barrier. Other contributions can be ascribed notably to Glu138 and Asp147 via ground-state destabilization and desolvation in the vicinity of the scissile bond. Key interactions with Trp118 and Trp181 were also proven to stabilize the ribooxocarbenium ion-like transition state. Altogether we propose that, as an alternative to a covalent acylal reaction intermediate with Glu218, catalysis by bCD38 proceeds through the formation of a discrete and transient ribooxocarbenium intermediate which is stabilized within the active site mostly by electrostatic interactions.

A single-residue change in the HIV-1 V3 loop associated with maraviroc resistance impairs CCR5 binding affinity while increasing replicative capacity.

BACKGROUND: Maraviroc (MVC) is an allosteric CCR5 inhibitor used against HIV-1 infection. While MVC-resistant viruses have been identified in patients, it still remains incompletely known how they adjust their CD4 and CCR5 binding properties to resist MVC inhibition while preserving their replicative capacity. It is thought that they maintain high efficiency of receptor binding. To date however, information about the binding affinities to receptors for inhibitor-resistant HIV-1 remains limited. RESULTS: Here, we show by means of viral envelope (gp120) binding experiments and virus-cell fusion kinetics that a MVC-resistant virus (MVC-Res) that had emerged as a dominant viral quasispecies in a patient displays reduced affinities for CD4 and CCR5 either free or bound to MVC, as compared to its MVC-sensitive counterpart isolated before MVC therapy. An alanine insertion within the GPG motif (G310_P311insA) of the MVC-resistant gp120 V3 loop is responsible for the decreased CCR5 binding affinity, while impaired binding to CD4 is due to sequence changes outside V3. Molecular dynamics simulations of gp120 binding to CCR5 further emphasize that the Ala insertion alters the structure of the V3 tip and weakens interaction with CCR5 ECL2. Paradoxically, infection experiments on cells expressing high levels of CCR5 also showed that Ala allows MVC-Res to use CCR5 efficiently, thereby improving viral fusion and replication efficiencies. Actually, although we found that the V3 loop of MVC-Res is required for high levels of MVC resistance, other regions outside V3 are sufficient to confer a moderate level of resistance. These sequence changes outside V3, however, come with a replication cost, which is compensated for by the Ala insertion in V3. CONCLUSION: These results indicate that changes in the V3 loop of MVC-resistant viruses can augment the efficiency of CCR5-dependent steps of viral entry other than gp120 binding, thereby compensating for their decreased affinity for entry receptors and improving their fusion and replication efficiencies. This study thus sheds light on unsuspected mechanisms whereby MVC-resistant HIV-1 could emerge and grow in treated patients.

Similarity between Flavonoid Biosynthetic Enzymes and Flavonoid Protein Targets Captured by Three-Dimensional Computing Approach.

Natural products are made by nature through interaction with biosynthetic enzymes. They also exert their effect as drugs by interaction with proteins. To address the question "Do biosynthetic enzymes and therapeutic targets share common mechanisms for the molecular recognition of natural products?", we compared the active site of five flavonoid biosynthetic enzymes to 8077 ligandable binding sites in the Protein Data Bank using two three-dimensional-based methods (SiteAlign and Shaper). Virtual screenings efficiently retrieved known flavonoid targets, in particular protein kinases. A consistent performance obtained for variable site descriptions (presence/absence of water, variable boundaries, or small structural changes) indicated that the methods are robust and thus well suited for the identification of potential target proteins of natural products. Finally, our results suggested that flavonoid binding is not primarily driven by shape, but rather by the recognition of common anchoring points.

Schistosoma mansoni NAD(+) catabolizing enzyme: identification of key residues in catalysis.

Schistosoma mansoni NAD(+) catabolizing enzyme (SmNACE), a distant homolog of mammalian CD38, shows significant structural and functional analogy to the members of the CD38/ADP-ribosyl cyclase family. The hallmark of SmNACE is the lack of ADP-ribosyl cyclase activity that might be ascribed to subtle changes in its active site. To better characterize the residues of the active site we determined the kinetic parameters of nine mutants encompassing three acidic residues: (i) the putative catalytic residue Glu202 and (ii) two acidic residues within the 'signature' region (the conserved Glu124 and the downstream Asp133), (iii) Ser169, a strictly conserved polar residue and (iv) two aromatic residues (His103 and Trp165). We established the very important role of Glu202 and of the hydrophobic domains overwhelmingly in the efficiency of the nicotinamide-ribosyl bond cleavage step. We also demonstrated that in sharp contrast with mammalian CD38, the 'signature' Glu124 is as critical as Glu202 for catalysis by the parasite enzyme. The different environments of the two Glu residues in the crystal structure of CD38 and in the homology model of SmNACE could explain such functional discrepancies. Mutagenesis data and 3D structures also indicated the importance of aromatic residues, especially His103, in the stabilization of the reaction intermediate as well as in the selection of its conformation suitable for cyclization to cyclic ADP-ribose. Finally, we showed that inhibition of SmNACE by the natural product cyanidin requires the integrity of Glu202 and Glu124, but not of His103 and Trp165, hence suggesting different recognition modes for substrate and inhibitor.

Structural analysis of neomycin B and kanamycin A binding Aminoglycosides Modifying Enzymes (AME) and bacterial ribosomal RNA.

Aminoglycosides are crucial antibiotics facing challenges from bacterial resistance. This study addresses the importance of aminoglycoside modifying enzymes in the context of escalating resistance. Drawing upon over two decades of structural data in the Protein Data Bank, we focused on two key antibiotics, neomycin B and kanamycin A, to explore how the aminoglycoside structure is exploited by this family of enzymes. A systematic comparison across diverse enzymes and the RNA A-site target identified common characteristics in the recognition mode, while assessing the adaptability of neomycin B and kanamycin A in various environments.

Exploration of the orthosteric/allosteric interface in human M1 muscarinic receptors by bitopic fluorescent ligands.

Bitopic binding properties apply to a variety of muscarinic compounds that span and simultaneously bind to both the orthosteric and allosteric receptor sites. We provide evidence that fluorescent pirenzepine derivatives, with the M1 antagonist fused to the boron-dipyrromethene [Bodipy (558/568)] fluorophore via spacers of varying lengths, exhibit orthosteric/allosteric binding properties at muscarinic M1 receptors. This behavior was inferred from a combination of functional, radioligand, and fluorescence resonance energy transfer binding experiments performed under equilibrium and kinetic conditions on enhanced green fluorescent protein-fused M1 receptors. Although displaying a common orthosteric component, the fluorescent compounds inherit bitopic properties from a linker-guided positioning of their Bodipy moiety within the M1 allosteric vestibule. Depending on linker length, the fluorophore is allowed to reach neighboring allosteric domains, overlapping or not with the classic gallamine site, but distinct from the allosteric indolocarbazole "WIN" site. Site-directed mutagenesis, as well as molecular modeling and ligand docking studies based on recently solved muscarinic receptor structures, further support the definition of two groups of Bodipy-pirenzepine derivatives exhibiting distinct allosteric binding poses. Thus, the linker may dictate pharmacological outcomes for bitopic molecules that are hardly predictable from the properties of individual orthosteric and allosteric building blocks. Our findings also demonstrate that the fusion of a fluorophore to an orthosteric ligand is not neutral, as it may confer, unless carefully controlled, unexpected properties to the resultant fluorescent tracer. Altogether, this study illustrates the importance of a "multifacet" experimental approach to unravel and validate bitopic ligand binding mechanisms.

Modeling the allosteric modulation of CCR5 function by Maraviroc.

Maraviroc is a non-peptidic, low molecular weight CC chemokine receptor 5 (CCR5) ligand that has recently been marketed for the treatment of HIV infected individuals. This review discusses recent molecular modeling studies of CCR5 by homology to CXC chemokine receptor 4, their contribution to the understanding of the allosteric mode of action of the inhibitor and their potential for the development of future drugs with improved efficiency and preservation of CCR5 biological functions.

Structural Insights into Molecular Recognition and Receptor Activation in Chemokine-Chemokine Receptor Complexes.

The chemokine system is a key player in the functioning of the immune system and a sought-after target for drug candidates. The number of experimental structures of chemokines in complex with chemokine receptors has increased rapidly over the past few years, providing essential information for rational development of chemokine receptor ligands. Here, we perform a comparative analysis of all chemokine-chemokine receptor structures, with the aim of characterizing the molecular recognition processes and highlighting the relationships between chemokine structures and functional processes. The structures show conserved interaction patterns between the chemokine core and the receptor N-terminus, while interactions near ECL2 display subfamily-specific features. Detailed analyses of the interactions of the chemokine N-terminal domain in the 7TM cavities reveal activation mechanisms for CCR5, CCR2, and CXCR2 and a mechanism for biased agonism in CCR1.

Predicting the duration of action of ß2-adrenergic receptor agonists: Ligand and structure-based approaches.

Agonists of the ß2 adrenergic receptor (ADRB2) are an important class of medications used for the treatment of respiratory diseases. They can be classified as short acting (SABA) or long acting (LABA), with each class playing a different role in patient management. In this work we explored both ligand-based and structure-based high-throughput approaches to classify ß2-agonists based on their duration of action. A completely in-silico prediction pipeline using an AlphaFold generated structure was used for structure-based modelling. Our analysis identified the ligands' 3D structure and lipophilicity as the most relevant features for the prediction of the duration of action. Interaction-based methods were also able to select ligands with the desired duration of action, incorporating the bias directly in the structure-based drug discovery pipeline without the need for further processing.

Mining the Protein Data Bank to inspire fragment library design.

The fragment approach has emerged as a method of choice for drug design, as it allows difficult therapeutic targets to be addressed. Success lies in the choice of the screened chemical library and the biophysical screening method, and also in the quality of the selected fragment and structural information used to develop a drug-like ligand. It has recently been proposed that promiscuous compounds, i.e., those that bind to several proteins, present an advantage for the fragment approach because they are likely to give frequent hits in screening. In this study, we searched the Protein Data Bank for fragments with multiple binding modes and targeting different sites. We identified 203 fragments represented by 90 scaffolds, some of which are not or hardly present in commercial fragment libraries. By contrast to other available fragment libraries, the studied set is enriched in fragments with a marked three-dimensional character (download at 10.5281/zenodo.7554649).

Allosteric model of maraviroc binding to CC chemokine receptor 5 (CCR5).

Maraviroc is a nonpeptidic small molecule human immunodeficiency virus type 1 (HIV-1) entry inhibitor that has just entered the therapeutic arsenal for the treatment of patients. We recently demonstrated that maraviroc binding to the HIV-1 coreceptor, CC chemokine receptor 5 (CCR5), prevents it from binding the chemokine CCL3 and the viral envelope glycoprotein gp120 by an allosteric mechanism. However, incomplete knowledge of ligand-binding sites and the lack of CCR5 crystal structures have hampered an in-depth molecular understanding of how the inhibitor works. Here, we addressed these issues by combining site-directed mutagenesis (SDM) with homology modeling and docking. Six crystal structures of G-protein-coupled receptors were compared for their suitability for CCR5 modeling. All CCR5 models had equally good geometry, but that built from the recently reported dimeric structure of the other HIV-1 coreceptor CXCR4 bound to the peptide CVX15 (Protein Data Bank code 3OE0) best agreed with the SDM data and discriminated CCR5 from non-CCR5 binders in a virtual screening approach. SDM and automated docking predicted that maraviroc inserts deeply in CCR5 transmembrane cavity where it can occupy three different binding sites, whereas CCL3 and gp120 lie on distinct yet overlapped regions of the CCR5 extracellular loop 2. Data suggesting that the transmembrane cavity remains accessible for maraviroc in CCL3-bound and gp120-bound CCR5 help explain our previous observation that the inhibitor enhances dissociation of preformed ligand-CCR5 complexes. Finally, we identified residues in the predicted CCR5 dimer interface that are mandatory for gp120 binding, suggesting that receptor dimerization might represent a target for new CCR5 entry inhibitors.

New insights into the mechanisms whereby low molecular weight CCR5 ligands inhibit HIV-1 infection.

CC chemokine receptor 5 (CCR5) is a G-protein-coupled receptor for the chemokines CCL3, -4, and -5 and a coreceptor for entry of R5-tropic strains of human immunodeficiency virus type 1 (HIV-1) into CD4(+) T-cells. We investigated the mechanisms whereby nonpeptidic, low molecular weight CCR5 ligands block HIV-1 entry and infection. Displacement binding assays and dissociation kinetics demonstrated that two of these molecules, i.e. TAK779 and maraviroc (MVC), inhibit CCL3 and the HIV-1 envelope glycoprotein gp120 binding to CCR5 by a noncompetitive and allosteric mechanism, supporting the view that they bind to regions of CCR5 distinct from the gp120- and CCL3-binding sites. We observed that TAK779 and MVC are full and weak inverse agonists for CCR5, respectively, indicating that they stabilize distinct CCR5 conformations with impaired abilities to activate G-proteins. Dissociation of [(125)I]CCL3 from CCR5 was accelerated by TAK779, to a lesser extent by MVC, and by GTP analogs, suggesting that inverse agonism contributes to allosteric inhibition of the chemokine binding to CCR5. TAK779 and MVC also promote dissociation of [(35)S]gp120 from CCR5 with an efficiency that correlates with their ability to act as inverse agonists. Displacement experiments revealed that affinities of MVC and TAK779 for the [(35)S]gp120-binding receptors are in the same range (IC(50) ∼6.4 versus 22 nm), although we found that MVC is 100-fold more potent than TAK779 for inhibiting HIV infection. This suggests that allosteric CCR5 inhibitors not only act by blocking gp120 binding but also alter distinct steps of CCR5 usage in the course of HIV infection.