Deciphering interaction fingerprints from protein molecular surfaces using geometric deep learning

Predicting interactions between proteins and other biomolecules solely based on structure remains a challenge in biology. A high-level representation of protein structure, the molecular surface, displays patterns of chemical and geometric features that fingerprint a protein’s modes of interactions with other biomolecules. We hypothesize that proteins participating in similar interactions may share common fingerprints, independent of their evolutionary history. Fingerprints may be difficult to grasp by visual analysis but could be learned from large-scale datasets. We present MaSIF (molecular surface interaction fingerprinting), a conceptual framework based on a geometric deep learning method to capture fingerprints that are important for specific biomolecular interactions. We showcase MaSIF with three prediction challenges: protein pocket-ligand prediction, protein–protein interaction site prediction and ultrafast scanning of protein surfaces for prediction of protein–protein complexes. We anticipate that our conceptual framework will lead to improvements in our understanding of protein function and design.

SPServer: split-statistical potentials for the analysis of protein structures and protein-protein interactions

Jaume Bonet Martinez

BackgroundStatistical potentials, also named knowledge-based potentials, are scoring functions derived from empirical data that can be used to evaluate the quality of protein folds and protein-protein interaction (PPI) structures. In previous works we decomposed the statistical potentials in different terms, named Split-Statistical Potentials, accounting for the type of amino acid pairs, their hydrophobicity, solvent accessibility and type of secondary structure. These potentials have been successfully used to identify near-native structures in protein structure prediction, rank protein docking poses, and predict PPI binding affinities.ResultsHere, we present the SPServer, a web server that applies the Split-Statistical Potentials to analyze protein folds and protein interfaces. SPServer provides global scores as well as residue/residue-pair profiles presented as score plots and maps. This level of detail allows users to: (1) identify potentially problematic regions on protein structures; (2) identify disrupting amino acid pairs in protein interfaces; and (3) compare and analyze the quality of tertiary and quaternary structural models. ConclusionsWhile there are many web servers that provide scoring functions to assess the quality of either protein folds or PPI structures, SPServer integrates both aspects in a unique easy-to-use web server. Moreover, the server permits to locally assess the quality of the structures and interfaces at a residue level and provides tools to compare the local assessment between structures.Server addresshttps://sbi.upf.edu/spserver/.

BMC2021

Peptide Structure and Radical Chemistry in the Gas Phase Revealed by Product Ion Abundance Analysis in Electron Capture and Transfer Dissociation Mass Spectrometry

Hisham Ben Hamidane

Peptide and protein ion abundances in mass spectrometry (MS) and tandem mass spectrometry (MS/MS) may provide orthogonal, or complementary, information to their mass measurements. Collision induced dissociation (CID), the most commonly used MS/MS technique, demonstrates high variability of fragmentation pathways and their relative product abundances as a function of experimental parameters. Moreover, a reliable prediction of relative abundances of CID products is possible only for relatively small systems, but not for larger biomolecules. Recently developed radical-induced fragmentation techniques, e.g., electron capture and transfer dissociation (ECD/ETD), appear to demonstrate more repeatable and reproducible fragmentation patterns. However, it has been reported that ECD/ETD product ion abundance distribution for peptides is rather flat and extracting additional peptide sequence information is yet to be done. In the present work, an attempt is made to understand the peptide's ECD/ETD product ion abundance formation fundamentals and to reveal the possible implications of this source of information for improving peptide and protein structural characterization. Specifically, we have shown that ECD/ETD of peptides is sensitive to the chemical environment in terms of peptide primary structure (amino acid nature) and possibly higher order structure (hydrogen bond connectivity). After demonstrating the analytical validity of the ECD/ETD product ion abundance methodology, its parameter dependence and its applicability to different experimental setups, results will be presented and rationalized along two main axes of comprehension: gas phase peptide radical chemistry and structure. Whereas fundamental aspects of ECD/ETD (mechanism, thermodynamics) can be addressed with the former, additional ECD/ETD yield modulations were observed and presumably relate to structural properties of ions in the gas phase. ECD specificity will be demonstrated on the class of amphipathic peptides, where periodicity is revealed in product ion abundance and correlated to peptide primary and secondary structure. To decipher between both contributions, sequence variations of amphipathic peptides were probed by ECD and correlated to gas phase structures suggested by ion mobility mass spectrometry and supported by molecular dynamics computations. After probing amino acid side chain influence, unraveling the chemical environment specificity of ECD/ETD led us to investigate peptide backbone modifications. Peptides containing beta amino acids revealed two main effects in ECD and ETD. We demonstrate that individual substitutions yield c and z type fragments only for amino acids that provide sufficient z• radical stabilization, whereas multiple substitutions shift the preferential fragmentation pathway to the otherwise minor a-type ion formation pathway indicating possible change in peptide protonation. Finally, understanding the partitioning of radical ions between N- and C-terminal fragments led us to investigate ion internal energy dependence. The extent of hydrogen atom transfer, the [c+z] complex lifetime and secondary fragmentation pathways (w-type ion formation) are directly modulated by ion internal energy as demonstrated by decoupling ETD from charge reduced CID. To summarize, ECD/ETD product ion abundance methodology developed in the thesis advances biomolecule radical chemistry understanding and provides a potentially useful methodology for peptide and protein structural characterization.

EPFL2011

Molecular Modeling of Bacterial Nanomachineries

Matteo Thomas Degiacomi

Proteins have the ability to assemble in multimeric states to perform their specific biological function. Unfortunately, characterizing experimentally these structures at atomistic resolution is usually difficult. For this reason, in silico methodologies aiming at predicting how multiple protein copies arrange to forma multimeric complex would be desirable. We present Parallel OptimizationWorkbench (POW), a swarm intelligence based optimization framework able to deal, in principle, with any optimization problem. We show that POW can be applied to biologically relevant problems such as prediction of protein assemblies and the parameterization of a Coarse-Grained force field for proteins. By combining POW optimizations, Molecular Dynamics simulations, Poisson-Boltzmann calculations and a variety of experiments, we subsequently study two bacterial nanomachieries: Aeromonas hydrophila's pore-forming toxin aerolysin, and Yersinia enterocolitica injectisome. These structures are challenging both for their size, and for the timescales involved in their functioning. Aerolysin is a pore-forming toxin secreted as an hydrophilic monomer. By means of large conformational changes, the protein heptamerizes on the target cell's surface, and finally inserts β-barrel into its lipid bilayer, causing cell death. The main hurdle in the study of this structure is the complexity of the mode of action, which spans timescales currently unreachable by classical molecular dynamics. We show that aerolysin C-terminal region has the dual role of preventing premature oligomerization and helping the folding of tertiary structure, qualifying therefore as an intramolecular chaperone. We study the transmembrane β-barrel properties and compare them with those of the homologous protein α-hemolysin. We show that aerolysin's barrel is more rigid than α-hemolysin's, and should be anion selective. We present models for aerolysin heptamer both in prepore and, for the first time, in membrane-inserted conformation. Our results are validated experimentally, and are consistent with known biochemical and structural data. The injectisome is an example of a type III secretion system. Its most striking feature is probably its size: hundreds of proteins assemble in a unique structure spanning the Gram-negative bacterial double membrane, and protruding outside the cell as a needle for tenth of nanometers. Obtaining an atomistic representation of this massive structure, and therefore some insights about its mode of action, is one of the greatest challenges. We show that the final length of injectisome's needle is determined by the secondary structure content of a ruler protein located inside its cavity during assembly. Using POW, we also produce the first model for Yersinia injectisome's basal body, highlighting the flexibility of this region in adapting between the inner and outer membranes. As a whole, this work demonstrates that a synergy of dry and wet experiments can provide precious insights into macromolecular structure and function.