Pred-O3, a web server to predict molecules, olfactory receptors and odor relationships.

The sense of smell is a biological process involving volatile molecules that interact with proteins called olfactory receptors to transmit a nervous message that allows the recognition of a perceived odor. However, the relationships between odorant molecules, olfactory receptors and odors (O3) are far from being well understood due to the combinatorial olfactory codes and large family of olfactory receptors. This is the reason why, based on 5802 odorant molecules and their annotations to 863 olfactory receptors (human) and 7029 odors and flavors annotations, a web server called Pred-O3 has been designed to provide insights into olfaction. Predictive models based on Artificial Intelligence have been developed allowing to suggest olfactory receptors and odors associated with a new molecule. In addition, based on the encoding of the odorant molecule's structure, physicochemical features related to odors and/or olfactory receptors are proposed. Finally, based on the structural models of the 98 olfactory receptors a systematic docking protocol can be applied and suggest if a molecule can bind or not to an olfactory receptor. Therefore, Pred-O3 is well suited to aid in the design of new odorant molecules and assist in fragrance research and sensory neuroscience. Pred-O3 is accessible at ' https://odor.rpbs.univ-paris-diderot.fr/'.

PEP-FOLD4: a pH-dependent force field for peptide structure prediction in aqueous solution.

Accurate and fast structure prediction of peptides of less 40 amino acids in aqueous solution has many biological applications, but their conformations are pH- and salt concentration-dependent. In this work, we present PEP-FOLD4 which goes one step beyond many machine-learning approaches, such as AlphaFold2, TrRosetta and RaptorX. Adding the Debye-Hueckel formalism for charged-charged side chain interactions to a Mie formalism for all intramolecular (backbone and side chain) interactions, PEP-FOLD4, based on a coarse-grained representation of the peptides, performs as well as machine-learning methods on well-structured peptides, but displays significant improvements for poly-charged peptides. PEP-FOLD4 is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD4. This server is free and there is no login requirement.

Exploring a Structural Data Mining Approach to Design Linkers for Head-to-Tail Peptide Cyclization.

Peptides have recently regained interest as therapeutic candidates, but their development remains confronted with several limitations including low bioavailability. Backbone head-to-tail cyclization, i.e., setting a covalent peptide bond linking the last amino acid with the first one, is one effective strategy of peptide-based drug design to stabilize the conformation of bioactive peptides while preserving peptide properties in terms of low toxicity, binding affinity, target selectivity, and preventing enzymatic degradation. Starting from an active peptide, it usually requires the design of a linker of a few amino acids to make it possible to cyclize the peptide, possibly preserving the conformation of the initial peptide and not affecting its activity. However, very little is known about the sequence-structure relationship requirements of designing linkers for peptide cyclization in a rational manner. Recently, we have shown that large-scale data-mining of available protein structures can lead to the precise identification of protein loop conformations, even from remote structural classes. Here, we transpose this approach to linkers, allowing head-to-tail peptide cyclization. First we show that given a linker sequence and the conformation of the linear peptide, it is possible to accurately predict the cyclized peptide conformation. Second, and more importantly, we show that it seems possible to elaborate on the information inferred from protein structures to propose effective candidate linker sequences constrained by length and amino acid composition, providing the first framework for the rational design of head-to-tail cyclization linkers. Finally, we illustrate this for two peptides using a limited set of amino-acids likely not to interfere with peptide function. For a linear peptide derived from Nrf2, the peptide cyclized starting from the experimental structure showed a 26-fold increase in the binding affinity. For urotensin II, a peptide already cyclized by a disulfide bond that exerts a broad array of biological activities, we were able, starting from models of the structure, to design a head-to-tail cyclized peptide, the first synthesized bicyclic 14-residue long urotensin II analogue, showing a retention of in vitro activity. Although preliminary, our results strongly suggest that such an approach has strong potential for cyclic peptide-based drug design.

InterEvDock3: a combined template-based and free docking server with increased performance through explicit modeling of complex homologs and integration of covariation-based contact maps.

The InterEvDock3 protein docking server exploits the constraints of evolution by multiple means to generate structural models of protein assemblies. The server takes as input either several sequences or 3D structures of proteins known to interact. It returns a set of 10 consensus candidate complexes, together with interface predictions to guide further experimental validation interactively. Three key novelties were implemented in InterEvDock3 to help obtain more reliable models: users can (i) generate template-based structural models of assemblies using close and remote homologs of known 3D structure, detected through an automated search protocol, (ii) select the assembly models most consistent with contact maps from external methods that implement covariation-based contact prediction with or without deep learning and (iii) exploit a novel coevolution-based scoring scheme at atomic level, which leads to significantly higher free docking success rates. The performance of the server was validated on two large free docking benchmark databases, containing respectively 230 unbound targets (Weng dataset) and 812 models of unbound targets (PPI4DOCK dataset). Its effectiveness has also been proven on a number of challenging examples. The InterEvDock3 web interface is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/InterEvDock3/.

Glutathionylation primes soluble glyceraldehyde-3-phosphate dehydrogenase for late collapse into insoluble aggregates.

Protein aggregation is a complex physiological process, primarily determined by stress-related factors revealing the hidden aggregation propensity of proteins that otherwise are fully soluble. Here we report a mechanism by which glycolytic glyceraldehyde-3-phosphate dehydrogenase of Arabidopsis thaliana (AtGAPC1) is primed to form insoluble aggregates by the glutathionylation of its catalytic cysteine (Cys149). Following a lag phase, glutathionylated AtGAPC1 initiates a self-aggregation process resulting in the formation of branched chains of globular particles made of partially misfolded and totally inactive proteins. GSH molecules within AtGAPC1 active sites are suggested to provide the initial destabilizing signal. The following removal of glutathione by the formation of an intramolecular disulfide bond between Cys149 and Cys153 reinforces the aggregation process. Physiological reductases, thioredoxins and glutaredoxins, could not dissolve AtGAPC1 aggregates but could efficiently contrast their growth. Besides acting as a protective mechanism against overoxidation, S-glutathionylation of AtGAPC1 triggers an unexpected aggregation pathway with completely different and still unexplored physiological implications.

Molecular dynamics simulations reveal statistics and microscopic mechanisms of water permeation in membrane-embedded artificial water channel nanoconstructs.

Understanding water transport mechanisms at the nanoscale level remains a challenge for theoretical chemical physics. Major advances in chemical synthesis have allowed us to discover new artificial water channels, rivaling with or even surpassing water conductance and selectivity of natural protein channels. In order to interpret experimental features and understand microscopic determinants for performance improvements, numerical approaches based on all-atom molecular dynamics simulations and enhanced sampling methods have been proposed. In this study, we quantify the influence of microscopic observables, such as channel radius and hydrogen bond connectivity, and of meso-scale features, such as the size of self-assembly blocks, on the permeation rate of a self-assembled nanocrystal-like artificial water channel. Although the absolute permeation rate extrapolated from these simulations is overestimated by one order of magnitude compared to the experimental measurement, the detailed analysis of several observed conductive patterns in large assemblies opens new pathways to scalable membranes with enhanced water conductance for the future design.

DaReUS-Loop: a web server to model multiple loops in homology models.

Loop regions in protein structures often have crucial roles, and they are much more variable in sequence and structure than other regions. In homology modeling, this leads to larger deviations from the homologous templates, and loop modeling of homology models remains an open problem. To address this issue, we have previously developed the DaReUS-Loop protocol, leading to significant improvement over existing methods. Here, a DaReUS-Loop web server is presented, providing an automated platform for modeling or remodeling loops in the context of homology models. This is the first web server accepting a protein with up to 20 loop regions, and modeling them all in parallel. It also provides a prediction confidence level that corresponds to the expected accuracy of the loops. DaReUS-Loop facilitates the analysis of the results through its interactive graphical interface and is freely available at http://bioserv.rpbs.univ-paris-diderot.fr/services/DaReUS-Loop/.

Coarse-grained and atomic resolution biomolecular docking with the ATTRACT approach.

The ATTRACT protein-protein docking program has been employed to predict protein-protein complex structures in CAPRI rounds 38-45. For 11 out of 16 targets acceptable or better quality solutions have been submitted (~70%). It includes also several cases of peptide-protein docking and the successful prediction of the geometry of carbohydrate-protein interactions. The option of combining rigid body minimization and simultaneous optimization in collective degrees of freedom based on elastic network modes was employed and systematically evaluated. Application to a large benchmark set indicates a modest improvement in docking performance compared to rigid docking. Possible further improvements of the docking approach in particular at the scoring and the flexible refinement steps are discussed.

String method solution of the gating pathways for a pentameric ligand-gated ion channel.

Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of ß-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.

Water permeation across artificial I-quartet membrane channels: from structure to disorder.

Artificial water channels (AWCs) have been designed for water transport across membranes with the aim to mimic the high water permeability observed for biological systems such as aquaporins (∼108-109 water molecules per s per channel), as well as their selectivity to reject ion permeation at the same time. Recent works on designed self-assembling alkylureido-ethylimidazole compounds forming imidazole-quartet channels (I-quartets), have shown both high water permeability and total ionic-rejection. I-quartets are thus promising candidates for further development of AWCs. However, the molecular mechanism of water permeation as well as I-quartet organization and stability in a membrane environment need to be fully understood to guide their optimal design. Here, we use a wide range of all-atom molecular dynamics (MD) simulations and their analysis to understand the structure/activity relationships of the I-quartet channels. Four different types with varying alkyl chain length or chirality have been studied in a complex fully hydrated lipid bilayer environment at both microsecond and nanosecond scale. Microsecond simulations show two distinct behaviors; (i) two out of four systems maintain chiral dipolar oriented water wires, but also undergo a strong reorganization of the crystal shape, (ii) the two other I-quartet channels completely lose the initial organization, nonetheless keeping a water transport activity. Short MD simulations with higher time resolution were conducted to characterize the dynamic properties of water molecules in these model channels and provided a detailed hypothesis on the molecular mechanism of water permeation. The ordered confined water was characterized with quantitative measures of hydrogen-bond life-time and single particle dynamics, showing variability among I-quartet channels. We will further discuss the underlying assumptions, currently based on self-aggregation simulations and crystal patches embedded in lipid bilayer simulations and attempt to describe possible alternative approaches to computationally capture the water permeation mechanism and the self-assembly process of these AWCs.

Allosteric and hyperekplexic mutant phenotypes investigated on an α1 glycine receptor transmembrane structure.

The glycine receptor (GlyR) is a pentameric ligand-gated ion channel (pLGIC) mediating inhibitory transmission in the nervous system. Its transmembrane domain (TMD) is the target of allosteric modulators such as general anesthetics and ethanol and is a major locus for hyperekplexic congenital mutations altering the allosteric transitions of activation or desensitization. We previously showed that the TMD of the human α1GlyR could be fused to the extracellular domain of GLIC, a bacterial pLGIC, to form a functional chimera called Lily. Here, we overexpress Lily in Schneider 2 insect cells and solve its structure by X-ray crystallography at 3.5 Å resolution. The TMD of the α1GlyR adopts a closed-channel conformation involving a single ring of hydrophobic residues at the center of the pore. Electrophysiological recordings show that the phenotypes of key allosteric mutations of the α1GlyR, scattered all along the pore, are qualitatively preserved in this chimera, including those that confer decreased sensitivity to agonists, constitutive activity, decreased activation kinetics, or increased desensitization kinetics. Combined structural and functional data indicate a pore-opening mechanism for the α1GlyR, suggesting a structural explanation for the effect of some key hyperekplexic allosteric mutations. The first X-ray structure of the TMD of the α1GlyR solved here using GLIC as a scaffold paves the way for mechanistic investigation and design of allosteric modulators of a human receptor.

Structural basis for ion permeation mechanism in pentameric ligand-gated ion channels.

To understand the molecular mechanism of ion permeation in pentameric ligand-gated ion channels (pLGIC), we solved the structure of an open form of GLIC, a prokaryotic pLGIC, at 2.4 Å. Anomalous diffraction data were used to place bound anions and cations. This reveals ordered water molecules at the level of two rings of hydroxylated residues (named Ser6' and Thr2') that contribute to the ion selectivity filter. Two water pentagons are observed, a self-stabilized ice-like water pentagon and a second wider water pentagon, with one sodium ion between them. Single-channel electrophysiology shows that the side-chain hydroxyl of Ser6' is crucial for ion translocation. Simulations and electrostatics calculations complemented the description of hydration in the pore and suggest that the water pentagons observed in the crystal are important for the ion to cross hydrophobic constriction barriers. Simulations that pull a cation through the pore reveal that residue Ser6' actively contributes to ion translocation by reorienting its side chain when the ion is going through the pore. Generalization of these findings to the pLGIC family is proposed.

Salt-Excluding Artificial Water Channels Exhibiting Enhanced Dipolar Water and Proton Translocation.

Aquaporins (AQPs) are biological water channels known for fast water transport (∼10(8)-10(9) molecules/s/channel) with ion exclusion. Few synthetic channels have been designed to mimic this high water permeability, and none reject ions at a significant level. Selective water translocation has previously been shown to depend on water-wires spanning the AQP pore that reverse their orientation, combined with correlated channel motions. No quantitative correlation between the dipolar orientation of the water-wires and their effects on water and proton translocation has been reported. Here, we use complementary X-ray structural data, bilayer transport experiments, and molecular dynamics (MD) simulations to gain key insights and quantify transport. We report artificial imidazole-quartet water channels with 2.6 Špores, similar to AQP channels, that encapsulate oriented dipolar water-wires in a confined chiral conduit. These channels are able to transport ∼10(6) water molecules/s, which is within 2 orders of magnitude of AQPs' rates, and reject all ions except protons. The proton conductance is high (∼5 H(+)/s/channel) and approximately half that of the M2 proton channel at neutral pH. Chirality is a key feature influencing channel efficiency.

Structural basis for alcohol modulation of a pentameric ligand-gated ion channel.

Despite its long history of use and abuse in human culture, the molecular basis for alcohol action in the brain is poorly understood. The recent determination of the atomic-scale structure of GLIC, a prokaryotic member of the pentameric ligand-gated ion channel (pLGIC) family, provides a unique opportunity to characterize the structural basis for modulation of these channels, many of which are alcohol targets in brain. We observed that GLIC recapitulates bimodal modulation by n-alcohols, similar to some eukaryotic pLGICs: methanol and ethanol weakly potentiated proton-activated currents in GLIC, whereas n-alcohols larger than ethanol inhibited them. Mapping of residues important to alcohol modulation of ionotropic receptors for glycine, γ-aminobutyric acid, and acetylcholine onto GLIC revealed their proximity to transmembrane cavities that may accommodate one or more alcohol molecules. Site-directed mutations in the pore-lining M2 helix allowed the identification of four residues that influence alcohol potentiation, with the direction of their effects reflecting α-helical structure. At one of the potentiation-enhancing residues, decreased side chain volume converted GLIC into a highly ethanol-sensitive channel, comparable to its eukaryotic relatives. Covalent labeling of M2 positions with an alcohol analog, a methanethiosulfonate reagent, further implicated residues at the extracellular end of the helix in alcohol binding. Molecular dynamics simulations elucidated the structural consequences of a potentiation-enhancing mutation and suggested a structural mechanism for alcohol potentiation via interaction with a transmembrane cavity previously termed the "linking tunnel." These results provide a unique structural model for independent potentiating and inhibitory interactions of n-alcohols with a pLGIC family member.

Spatial proximity between the VPAC1 receptor and the amino terminus of agonist and antagonist peptides reveals distinct sites of interaction.

Vasoactive intestinal peptide (VIP) plays a major role in pathophysiology. Our previous studies demonstrated that the VIP sequence 6-28 interacts with the N-terminal ectodomain (N-ted) of its receptor, VPAC1. Probes for VIP and receptor antagonist PG97-269 were synthesized with a photolabile residue/Bpa at various positions and used to explore spatial proximity with VPAC1. PG97-269 probes with Bpa at position 0, 6, and 24 behaved as high-affinity receptor antagonists (K(i)=12, 9, and 7 nM, respectively). Photolabeling experiments revealed that the [Bpa(0)]-VIP probe was in physical contact with VPAC1 Q(135), while [Bpa(0)]-PG97-269 was covalently bound to G(62) residue of N-ted, indicating different binding sites. In contrast, photolabeling with [Bpa(6)]- and [Bpa(24)]-PG97-269 showed that the distal domains of PG97-269 interacted with N-ted, as we previously showed for VIP. Substitution with alanine of the K(143), T(144), and T(147) residues located in the first transmembrane domain of VPAC1 induced a loss of receptor affinity (IC(50)=1035, 874, and 2070 nM, respectively), and pharmacological studies using VIP2-28 indicated that these three residues play an important role in VPAC1 interaction with the first histidine residue of VIP. These data demonstrate that VIP and PG97-269 bind to distinct domains of VPAC1.

Molecular mechanism for the dual alcohol modulation of Cys-loop receptors.

Cys-loop receptors constitute a superfamily of pentameric ligand-gated ion channels (pLGICs), including receptors for acetylcholine, serotonin, glycine and γ-aminobutyric acid. Several bacterial homologues have been identified that are excellent models for understanding allosteric binding of alcohols and anesthetics in human Cys-loop receptors. Recently, we showed that a single point mutation on a prokaryotic homologue (GLIC) could transform it from a channel weakly potentiated by ethanol into a highly ethanol-sensitive channel. Here, we have employed molecular simulations to study ethanol binding to GLIC, and to elucidate the role of the ethanol-enhancing mutation in GLIC modulation. By performing 1-µs simulations with and without ethanol on wild-type and mutated GLIC, we observed spontaneous binding in both intra-subunit and inter-subunit transmembrane cavities. In contrast to the glycine receptor GlyR, in which we previously observed ethanol binding primarily in an inter-subunit cavity, ethanol primarily occupied an intra-subunit cavity in wild-type GLIC. However, the highly ethanol-sensitive GLIC mutation significantly enhanced ethanol binding in the inter-subunit cavity. These results demonstrate dramatic effects of the F(14')A mutation on the distribution of ligands, and are consistent with a two-site model of pLGIC inhibition and potentiation.

Class-B GPCR activation: is ligand helix-capping the key?

The class B family of G-protein-coupled receptors (GPCRs) regulates essential physiological functions such as exocrine and endocrine secretions, feeding behaviour, metabolism, growth, and neuro- and immuno-modulations. These receptors are activated by endogenous peptide hormones including secretin, glucagon, vasoactive intestinal peptide, corticotropin-releasing factor and parathyroid hormone. We have identified a common structural motif that is encoded in all class B GPCR-ligand N-terminal sequences. We propose that this local structure, a helix N-capping motif, is formed upon receptor binding and constitutes a key element underlying class B GPCR activation. The folded backbone conformation imposed by the capping structure could serve as a template for a rational design of drugs targeting class B GPCRs in several diseases.