Evolutionarily conserved protein motifs drive interactions between the plant nucleoskeleton and nuclear pores.

The nucleoskeleton forms a filamentous meshwork under the nuclear envelope and contributes to the regulation of nuclear shape and gene expression. To understand how the Arabidopsis (Arabidopsis thaliana) nucleoskeleton physically connects to the nuclear periphery in plants, we investigated the Arabidopsis nucleoskeleton protein KAKU4 and sought for functional regions responsible for its localization at the nuclear periphery. We identified 3 conserved peptide motifs within the N-terminal region of KAKU4, which are required for intermolecular interactions of KAKU4 with itself, interaction with the nucleoskeleton protein CROWDED NUCLEI (CRWN), localization at the nuclear periphery, and nuclear elongation in differentiated tissues. Unexpectedly, we find these motifs to be present also in NUP82 and NUP136, 2 plant-specific nucleoporins from the nuclear pore basket. We further show that NUP82, NUP136, and KAKU4 have a common evolutionary history predating nonvascular land plants with KAKU4 mainly localizing outside the nuclear pore suggesting its divergence from an ancient nucleoporin into a new nucleoskeleton component. Finally, we demonstrate that both NUP82 and NUP136, through their shared N-terminal motifs, interact with CRWN and KAKU4 proteins revealing the existence of a physical continuum between the nuclear pore and the nucleoskeleton in plants.

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.

DockSurf: A Molecular Modeling Software for the Prediction of Protein/Surface Adhesion.

The elucidation of structural interfaces between proteins and inorganic surfaces is a crucial aspect of bionanotechnology development. Despite its significance, the interfacial structures between proteins and metallic surfaces are yet to be fully understood, and the lack of experimental investigation has impeded the development of many devices. To overcome this limitation, we suggest considering the generation of protein/surface structures as a molecular docking problem with a homogenous plan as the target. To this extent, we propose a new software, DockSurf, which aims to quickly propose reliable protein/surface structures. Our approach considers the conformational exploration with Euler's angles, which provide a cartography instead of a unique structure. Interaction energies were derived from quantum mechanics computations for a set of small molecules that describe protein atom types and implemented in a Derjaguin, Landau, Verwey, and Overbeek potential for the consideration of large systems such as proteins. The validation of DockSurf software was conducted with molecular dynamics for corona proteins with gold surfaces and provided enthusiastic results. This software is implemented in the RPBS platform to facilitate widespread access to the scientific community.

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/.

Proteo3Dnet: a web server for the integration of structural information with interactomics data.

Proteo3Dnet is a web server dedicated to the analysis of mass spectrometry interactomics experiments. Given a flat list of proteins, its aim is to organize it in terms of structural interactions to provide a clearer overview of the data. This is achieved using three means: (i) the search for interologs with resolved structure available in the protein data bank, including cross-species remote homology search, (ii) the search for possibly weaker interactions mediated through Short Linear Motifs as predicted by ELM-a unique feature of Proteo3Dnet, (iii) the search for protein-protein interactions physically validated in the BioGRID database. The server then compiles this information and returns a graph of the identified interactions and details about the different searches. The graph can be interactively explored to understand the way the core complexes identified could interact. It can also suggest undetected partners to the experimentalists, or specific cases of conditionally exclusive binding. The interest of Proteo3Dnet, previously demonstrated for the difficult cases of the proteasome and pragmin complexes data is, here, illustrated in the context of yeast precursors to the small ribosomal subunits and the smaller interactome of 14-3-3zeta frequent interactors. The Proteo3Dnet web server is accessible at http://bioserv.rpbs.univ-paris-diderot.fr/services/Proteo3Dnet/.

The CH1α domain of mucosal gp41 IgA contributes to antibody specificity and antiviral functions in HIV-1 highly exposed Sero-Negative individuals.

The antibody molecule comprises a variable domain conferring antigen specificity and affinity distinct from the heavy chain constant (CH) domains dictating effector functions. We here interrogate this paradigm by evaluating the unique influence of the CH1α domain on epitope specificity and functions using two mucosal gp41-specific Fab-IgAs (FabA) derived from HIV-1 highly-exposed but persistently seronegative individuals (HESN). These HESN develop selectively affinity-matured HIV-1-specific mucosal IgA that target the gp41 viral envelope and might provide protection although by unclear mechanisms. Isotype-switching FabAs into Fab-IgGs (FabGs) results in a >10-fold loss in affinity for HIV-1 clade A, B, and C gp41, together with reduced neutralization of HIV-1 cross-clade. The FabA conformational epitopes map selectively on gp41 in 6-Helix bundle and pre-fusion conformations cross-clade, unlike FabGs. Finally, we designed in silico, a 12 amino-acid peptide recapitulating one FabA conformational epitope that inhibits the FabA binding to gp41 cross-clade and its neutralizing activity. Altogether, our results reveal that the CH1α domain shapes the antibody paratope through an allosteric effect, thereby strengthening the antibody specificity and functional activities. Further, they clarify the mechanisms by which these HESN IgAs might confer protection against HIV-1-sexual acquisition. The IgA-specific epitope we characterized by reverse vaccinology could help designing a mucosal HIV-1 vaccine.

PatchSearch: a web server for off-target protein identification.

The large number of proteins found in the human body implies that a drug may interact with many proteins, called off-target proteins, besides its intended target. The PatchSearch web server provides an automated workflow that allows users to identify structurally conserved binding sites at the protein surfaces in a set of user-supplied protein structures. Thus, this web server may help to detect potential off-target protein. It takes as input a protein complexed with a ligand and identifies within user-defined or predefined collections of protein structures, those having a binding site compatible with this ligand in terms of geometry and physicochemical properties. It is based on a non-sequential local alignment of the patch over the entire protein surface. Then the PatchSearch web server proposes a ligand binding mode for the potential off-target, as well as an estimated affinity calculated by the Vinardo scoring function. This novel tool is able to efficiently detects potential interactions of ligands with distant off-target proteins. Furthermore, by facilitating the discovery of unexpected off-targets, PatchSearch could contribute to the repurposing of existing drugs. The server is freely available at http://bioserv.rpbs.univ-paris-diderot.fr/services/PatchSearch.

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/.

Probing Protein Interaction Networks by Combining MS-Based Proteomics and Structural Data Integration.

Protein-protein interactions play a major role in the molecular machinery of life, and various techniques such as AP-MS are dedicated to their identification. However, those techniques return lists of proteins devoid of organizational structure, not detailing which proteins interact with which others. Proposing a hierarchical view of the interactions between the members of the flat list becomes highly tedious for large data sets when done by hand. To help hierarchize this data, we introduce a new bioinformatics protocol that integrates information of the multimeric protein 3D structures available in the Protein Data Bank using remote homology detection, as well as information related to Short Linear Motifs and interaction data from the BioGRID. We illustrate on two unrelated use-cases of different complexity how our approach can be useful to decipher the network of interactions hidden in the list of input proteins, and how it provides added value compared to state-of-the-art resources such as Interactome3D or STRING. Particularly, we show the added value of using homology detection to distinguish between orthologs and paralogs, and to distinguish between core obligate and more facultative interactions. We also demonstrate the potential of considering interactions occurring through Short Linear Motifs.

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.

InterEvDock2: an expanded server for protein docking using evolutionary and biological information from homology models and multimeric inputs.

Computational protein docking is a powerful strategy to predict structures of protein-protein interactions and provides crucial insights for the functional characterization of macromolecular cross-talks. We previously developed InterEvDock, a server for ab initio protein docking based on rigid-body sampling followed by consensus scoring using physics-based and statistical potentials, including the InterEvScore function specifically developed to incorporate co-evolutionary information in docking. InterEvDock2 is a major evolution of InterEvDock which allows users to submit input sequences - not only structures - and multimeric inputs and to specify constraints for the pairwise docking process based on previous knowledge about the interaction. For this purpose, we added modules in InterEvDock2 for automatic template search and comparative modeling of the input proteins. The InterEvDock2 pipeline was benchmarked on 812 complexes for which unbound homology models of the two partners and co-evolutionary information are available in the PPI4DOCK database. InterEvDock2 identified a correct model among the top 10 consensus in 29% of these cases (compared to 15-24% for individual scoring functions) and at least one correct interface residue among 10 predicted in 91% of these cases. InterEvDock2 is thus a unique protein docking server, designed to be useful for the experimental biology community. The InterEvDock2 web interface is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/InterEvDock2/.

Probing the quality control mechanism of the Escherichia coli twin-arginine translocase with folding variants of a de novo-designed heme protein.

Protein transport across the cytoplasmic membrane of bacterial cells is mediated by either the general secretion (Sec) system or the twin-arginine translocase (Tat). The Tat machinery exports folded and cofactor-containing proteins from the cytoplasm to the periplasm by using the transmembrane proton motive force as a source of energy. The Tat apparatus apparently senses the folded state of its protein substrates, a quality-control mechanism that prevents premature export of nascent unfolded or misfolded polypeptides, but its mechanistic basis has not yet been determined. Here, we investigated the innate ability of the model Escherichia coli Tat system to recognize and translocate de novo-designed protein substrates with experimentally determined differences in the extent of folding. Water-soluble, four-helix bundle maquette proteins were engineered to bind two, one, or no heme b cofactors, resulting in a concomitant reduction in the extent of their folding, assessed with temperature-dependent CD spectroscopy and one-dimensional 1H NMR spectroscopy. Fusion of the archetypal N-terminal Tat signal peptide of the E. coli trimethylamine-N-oxide (TMAO) reductase (TorA) to the N terminus of the protein maquettes was sufficient for the Tat system to recognize them as substrates. The clear correlation between the level of Tat-dependent export and the degree of heme b-induced folding of the maquette protein suggested that the membrane-bound Tat machinery can sense the extent of folding and conformational flexibility of its substrates. We propose that these artificial proteins are ideal substrates for future investigations of the Tat system's quality-control mechanism.

The pepATTRACT web server for blind, large-scale peptide-protein docking.

Peptide-protein interactions are ubiquitous in the cell and form an important part of the interactome. Computational docking methods can complement experimental characterization of these complexes, but current protocols are not applicable on the proteome scale. pepATTRACT is a novel docking protocol that is fully blind, i.e. it does not require any information about the binding site. In various stages of its development, pepATTRACT has participated in CAPRI, making successful predictions for five out of seven protein-peptide targets. Its performance is similar or better than state-of-the-art local docking protocols that do require binding site information. Here we present a novel web server that carries out the rigid-body stage of pepATTRACT. On the peptiDB benchmark, the web server generates a correct model in the top 50 in 34% of the cases. Compared to the full pepATTRACT protocol, this leads to some loss of performance, but the computation time is reduced from ∼18 h to ∼10 min. Combined with the fact that it is fully blind, this makes the web server well-suited for large-scale in silico protein-peptide docking experiments. The rigid-body pepATTRACT server is freely available at http://bioserv.rpbs.univ-paris-diderot.fr/services/pepATTRACT.

A Free Web-Based Protocol to Assist Structure-Based Virtual Screening Experiments.

Chemical biology and drug discovery are complex and costly processes. In silico screening approaches play a key role in the identification and optimization of original bioactive molecules and increase the performance of modern chemical biology and drug discovery endeavors. Here, we describe a free web-based protocol dedicated to small-molecule virtual screening that includes three major steps: ADME-Tox filtering (via the web service FAF-Drugs4), docking-based virtual screening (via the web service MTiOpenScreen), and molecular mechanics optimization (via the web service AMMOS2 [Automatic Molecular Mechanics Optimization for in silico Screening]). The online tools FAF-Drugs4, MTiOpenScreen, and AMMOS2 are implemented in the freely accessible RPBS (Ressource Parisienne en Bioinformatique Structurale) platform. The proposed protocol allows users to screen thousands of small molecules and to download the top 1500 docked molecules that can be further processed online. Users can then decide to purchase a small list of compounds for in vitro validation. To demonstrate the potential of this online-based protocol, we performed virtual screening experiments of 4574 approved drugs against three cancer targets. The results were analyzed in the light of published drugs that have already been repositioned on these targets. We show that our protocol is able to identify active drugs within the top-ranked compounds. The web-based protocol is user-friendly and can successfully guide the identification of new promising molecules for chemical biology and drug discovery purposes.

InterEvDock: a docking server to predict the structure of protein-protein interactions using evolutionary information.

The structural modeling of protein-protein interactions is key in understanding how cell machineries cross-talk with each other. Molecular docking simulations provide efficient means to explore how two unbound protein structures interact. InterEvDock is a server for protein docking based on a free rigid-body docking strategy. A systematic rigid-body docking search is performed using the FRODOCK program and the resulting models are re-scored with InterEvScore and SOAP-PP statistical potentials. The InterEvScore potential was specifically designed to integrate co-evolutionary information in the docking process. InterEvDock server is thus particularly well suited in case homologous sequences are available for both binding partners. The server returns 10 structures of the most likely consensus models together with 10 predicted residues most likely involved in the interface. In 91% of all complexes tested in the benchmark, at least one residue out of the 10 predicted is involved in the interface, providing useful guidelines for mutagenesis. InterEvDock is able to identify a correct model among the top10 models for 49% of the rigid-body cases with evolutionary information, making it a unique and efficient tool to explore structural interactomes under an evolutionary perspective. The InterEvDock web interface is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/InterEvDock/.

PEP-FOLD3: faster de novo structure prediction for linear peptides in solution and in complex.

Structure determination of linear peptides of 5-50 amino acids in aqueous solution and interacting with proteins is a key aspect in structural biology. PEP-FOLD3 is a novel computational framework, that allows both (i) de novo free or biased prediction for linear peptides between 5 and 50 amino acids, and (ii) the generation of native-like conformations of peptides interacting with a protein when the interaction site is known in advance. PEP-FOLD3 is fast, and usually returns solutions in a few minutes. Testing PEP-FOLD3 on 56 peptides in aqueous solution led to experimental-like conformations for 80% of the targets. Using a benchmark of 61 peptide-protein targets starting from the unbound form of the protein receptor, PEP-FOLD3 was able to generate peptide poses deviating on average by 3.3Å from the experimental conformation and return a native-like pose in the first 10 clusters for 52% of the targets. PEP-FOLD3 is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD3.

BCSearch: fast structural fragment mining over large collections of protein structures.

Resources to mine the large amount of protein structures available today are necessary to better understand how amino acid variations are compatible with conformation preservation, to assist protein design, engineering and, further, the development of biologic therapeutic compounds. BCSearch is a versatile service to efficiently mine large collections of protein structures. It relies on a new approach based on a Binet-Cauchy kernel that is more discriminative than the widely used root mean square deviation criterion. It has statistics independent of size even for short fragments, and is fast. The systematic mining of large collections of structures such as the complete SCOPe protein structural classification or comprehensive subsets of the Protein Data Bank can be performed in few minutes. Based on this new score, we propose four innovative applications: BCFragSearch and BCMirrorSearch, respectively, search for fragments similar and anti-similar to a query and return information on the diversity of the sequences of the hits. BCLoopSearch identifies candidate fragments of fixed size matching the flanks of a gaped structure. BCSpecificitySearch analyzes a complete protein structure and returns information about sites having few similar fragments. BCSearch is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/BCSearch.

MTiOpenScreen: a web server for structure-based virtual screening.

Open screening endeavors play and will play a key role to facilitate the identification of new bioactive compounds in order to foster innovation and to improve the effectiveness of chemical biology and drug discovery processes. In this line, we developed the new web server MTiOpenScreen dedicated to small molecule docking and virtual screening. It includes two services, MTiAutoDock and MTiOpenScreen, allowing performing docking into a user-defined binding site or blind docking using AutoDock 4.2 and automated virtual screening with AutoDock Vina. MTiOpenScreen provides valuable starting collections for screening, two in-house prepared drug-like chemical libraries containing 150 000 PubChem compounds: the Diverse-lib containing diverse molecules and the iPPI-lib enriched in molecules likely to inhibit protein-protein interactions. In addition, MTiOpenScreen offers users the possibility to screen up to 5000 small molecules selected outside our two libraries. The predicted binding poses and energies of up to 1000 top ranked ligands can be downloaded. In this way, MTiOpenScreen enables researchers to apply virtual screening using different chemical libraries on traditional or more challenging protein targets such as protein-protein interactions. The MTiOpenScreen web server is free and open to all users at http://bioserv.rpbs.univ-paris-diderot.fr/services/MTiOpenScreen/.

Missense mutations in SLC26A8, encoding a sperm-specific activator of CFTR, are associated with human asthenozoospermia.

The cystic fibrosis transmembrane conductance regulator (CFTR) is present in mature sperm and is required for sperm motility and capacitation. Both these processes are controlled by ions fluxes and are essential for fertilization. We have shown that SLC26A8, a sperm-specific member of the SLC26 family of anion exchangers, associates with the CFTR channel and strongly stimulates its activity. This suggests that the two proteins cooperate to regulate the anion fluxes required for correct sperm motility and capacitation. Here, we report on three heterozygous SLC26A8 missense mutations identified in a cohort of 146 men presenting with asthenozoospermia: c.260G>A (p.Arg87Gln), c.2434G>A (p.Glu812Lys), and c.2860C>T (p.Arg954Cys). These mutations were not present in 121 controls matched for ethnicity, and statistical analysis on a control population of 8,600 individuals (from dbSNP and 1000 Genomes) showed them to be associated with asthenozoospermia with a power > 95%. By cotransfecting Chinese hamster ovary (CHO)-K1 cells with SLC26A8 variants and CFTR, we showed that the physical interaction between the two proteins was partly conserved but that the capacity to activate CFTR-dependent anion transport was completely abolished for all mutants. Biochemical studies revealed the presence of much smaller amounts of protein for all variants, but these amounts were restored to wild-type levels upon treatment with the proteasome inhibitor MG132. Immunocytochemistry also showed the amounts of SLC26A8 in sperm to be abnormally small in individuals carrying the mutations. These mutations might therefore impair formation of the SLC26A8-CFTR complex, principally by affecting SLC26A8 stability, consistent with an impairment of CFTR-dependent sperm-activation events in affected individuals.