MolArt: a molecular structure annotation and visualization tool

Visualization of automatically combined disease maps and pathway diagrams for rare diseases

Introduction: Investigation of molecular mechanisms of human disorders, especially rare diseases, require exploration of various knowledge repositories for building precise hypotheses and complex data interpretation. Recently, increasingly more resources offer diagrammatic representation of such mechanisms, including disease-dedicated schematics in pathway databases and disease maps. However, collection of knowledge across them is challenging, especially for research projects with limited manpower. Methods: In this article we present an automated workflow for construction of maps of molecular mechanisms for rare diseases. The workflow requires a standardized definition of a disease using Orphanet or HPO identifiers to collect relevant genes and variants, and to assemble a functional, visual repository of related mechanisms, including data overlays. The diagrams composing the final map are unified to a common systems biology format from CellDesigner SBML, GPML and SBML+layout+render. The constructed resource contains disease-relevant genes and variants as data overlays for immediate visual exploration, including embedded genetic variant browser and protein structure viewer. Results: We demonstrate the functionality of our workflow on two examples of rare diseases: Kawasaki disease and retinitis pigmentosa. Two maps are constructed based on their corresponding identifiers. Moreover, for the retinitis pigmentosa use-case, we include a list of differentially expressed genes to demonstrate how to tailor the workflow using omics datasets. Discussion: In summary, our work allows for an ad-hoc construction of molecular diagrams combined from different sources, preserving their layout and graphical style, but integrating them into a single resource. This allows to reduce time consuming tasks of prototyping of a molecular disease map, enabling visual exploration, hypothesis building, data visualization and further refinement. The code of the workflow is open and accessible at https://gitlab.lcsb.uni.lu/minerva/automap/.

The Metal-binding Protein Atlas (MbPA): an integrated database for curating metalloproteins in all aspects

Metal-binding proteins are essential for the vital activities and engage in their roles by acting in concert with metal cations. MbPA (The Metal-binding Protein Atlas) is the most comprehensive resource up to now dedicated to curating metal-binding proteins. Currently, it contains 106373 entries and 440187 sites related to 54 metals and 8169 species. Users can view all metal-binding proteins and species-specific proteins in MbPA. There are also metal-proteomics data that quantitatively describes protein expression in different tissues and organs. By analyzing the data of the amino acid residues at the metal-binding site, it is found that about 80% of the metal ions tend to bind to cysteine, aspartic acid, glutamic acid, and histidine. Moreover, we use Diversity Measure to confirm that the diversity of metal-binding is specific in different area of periodic table, and further elucidate the binding modes of 19 transition metals on 20 amino acids. In addition, MbPA also embraces 6855 potential pathogenic mutations related to metalloprotein. The resource is freely available at http://bioinfor.imu.edu.cn/mbpa.

Bioinformatics Tools and Knowledgebases to Assist Generating Targeted Assays for Plasma Proteomics

In targeted proteomics experiments, selecting the appropriate proteotypic peptides as surrogate for the target protein is a crucial pre-acquisition step. This step is largely a bioinformatics exercise that involves integrating information on the peptides and proteins and using various software tools and knowledgebases. We present here a few resources that automate and simplify the selection process to a great degree. These tools and knowledgebases were developed primarily to streamline targeted proteomics assay development and include PeptidePicker, PeptidePickerDB, MRMAssayDB, MouseQuaPro, and PeptideTracker. We have used these tools to develop and document thousands of targeted proteomics assays, many of them for plasma proteins with focus on human and mouse. An important aspect in all these resources is the integrative approach on which they are based. Using these tools in the first steps of designing a singleplexed or multiplexed targeted proteomic experiment can reduce the necessary experimental steps tremendously. All the tools and knowledgebases we describe here are Web-based and freely accessible so scientists can query the information conveniently from the browser. This chapter provides an overview of these software tools and knowledgebases, their content, and how to use them for targeted plasma proteomics. We further demonstrate how to use them with the results of the HUPO Human Plasma Proteome Project to produce a new database of 3.8 k targeted assays for known human plasma proteins. Upon experimental validation, these assays should help in the further quantitative characterizing of the plasma proteome.

Delineation of functionally essential protein regions for 242 neurodevelopmental genes

Neurodevelopmental disorders (NDDs), including severe paediatric epilepsy, autism and intellectual disabilities are heterogeneous conditions in which clinical genetic testing can often identify a pathogenic variant. For many of them, genetic therapies will be tested in this or the coming years in clinical trials. In contrast to first-generation symptomatic treatments, the new disease-modifying precision medicines require a genetic test-informed diagnosis before a patient can be enrolled in a clinical trial. However, even in 2022, most identified genetic variants in NDD genes are 'variants of uncertain significance'. To safely enrol patients in precision medicine clinical trials, it is important to increase our knowledge about which regions in NDD-associated proteins can 'tolerate' missense variants and which ones are 'essential' and will cause a NDD when mutated. In addition, knowledge about functionally indispensable regions in the 3D structure context of proteins can also provide insights into the molecular mechanisms of disease variants. We developed a novel consensus approach that overlays evolutionary, and population based genomic scores to identify 3D essential sites (Essential3D) on protein structures. After extensive benchmarking of AlphaFold predicted and experimentally solved protein structures, we generated the currently largest expert curated protein structure set for 242 NDDs and identified 14 377 Essential3D sites across 189 gene disorders associated proteins. We demonstrate that the consensus annotation of Essential3D sites improves prioritization of disease mutations over single annotations. The identified Essential3D sites were enriched for functional features such as intermembrane regions or active sites and discovered key inter-molecule interactions in protein complexes that were otherwise not annotated. Using the currently largest autism, developmental disorders, and epilepsies exome sequencing studies including >360 000 NDD patients and population controls, we found that missense variants at Essential3D sites are 8-fold enriched in patients. In summary, we developed a comprehensive protein structure set for 242 NDDs and identified 14 377 Essential3D sites in these. All data are available at https://es-ndd.broadinstitute.org for interactive visual inspection to enhance variant interpretation and development of mechanistic hypotheses for 242 NDDs genes. The provided resources will enhance clinical variant interpretation and in silico drug target development for NDD-associated genes and encoded proteins.

Exploration and comparison of molecular mechanisms across diseases using MINERVA Net

Protein function is often interpreted using molecular interaction diagrams, encoding roles a given protein plays in various molecular mechanisms. Information about disease-related mechanisms can be inferred from disease maps, knowledge repositories containing manually constructed systems biology diagrams. Disease maps hosted on the Molecular Interaction Network VisuAlization (MINERVA) Platform are individually accessible through a REST API interface of each instance, making it challenging to systematically explore their contents. To address this challenge, we introduce the MINERVA Net web service, a repository of open-access disease maps allowing users to publicly share minimal information about their maps. The MINERVA Net repository provides REST API endpoints of particular disease maps, which then can be individually queried for content. In this article, we describe the concept of MINERVA Net and illustrate its use by comparing proteins and their interactions in three different disease maps.

GWYRE: A resource for mapping variants onto experimental and modeled structures of human protein complexes

•

Structure of protein complexes is important for interpreting genetic variation.

•

Data on single amino acid variants is available from high-throughput sequencing.

•

Integrated modeling approach was applied to proteins and their complexes.

•

GWYRE resource incorporates predicted protein complexes with mapped mutations.

Rapid progress in structural modeling of proteins and their interactions is powered by advances in knowledge-based methodologies along with better understanding of physical principles of protein structure and function. The pool of structural data for modeling of proteins and protein–protein complexes is constantly increasing due to the rapid growth of protein interaction databases and Protein Data Bank. The GWYRE (Genome Wide PhYRE) project capitalizes on these developments by advancing and applying new powerful modeling methodologies to structural modeling of protein–protein interactions and genetic variation. The methods integrate knowledge-based tertiary structure prediction using Phyre2 and quaternary structure prediction using template-based docking by a full-structure alignment protocol to generate models for binary complexes. The predictions are incorporated in a comprehensive public resource for structural characterization of the human interactome and the location of human genetic variants. The GWYRE resource facilitates better understanding of principles of protein interaction and structure/function relationships. The resource is available at http://www.gwyre.org.

OUP accepted manuscript

Motivation

Mapping positional features from one-dimensional (1D) sequences onto three-dimensional (3D) structures of biological macromolecules is a powerful tool to show geometric patterns of biochemical annotations and provide a better understanding of the mechanisms underpinning protein and nucleic acid function at the atomic level.

Results

We present a new library designed to display fully customizable interactive views between 1D positional features of protein and/or nucleic acid sequences and their 3D structures as isolated chains or components of macromolecular assemblies.

Availability and implementation

https://github.com/rcsb/rcsb-saguaro-3d.

Supplementary information

Supplementary data are available at Bioinformatics online.

Humans and machines in biomedical knowledge curation: hypertrophic cardiomyopathy molecular mechanisms' representation

Background

Biomedical knowledge is dispersed in scientific literature and is growing constantly. Curation is the extraction of knowledge from unstructured data into a computable form and could be done manually or automatically. Hypertrophic cardiomyopathy (HCM) is the most common inherited cardiac disease, with genotype–phenotype associations still incompletely understood. We compared human- and machine-curated HCM molecular mechanisms’ models and examined the performance of different machine approaches for that task.

Results

We created six models representing HCM molecular mechanisms using different approaches and made them publicly available, analyzed them as networks, and tried to explain the models’ differences by the analysis of factors that affect the quality of machine-curated models (query constraints and reading systems’ performance). A result of this work is also the Interactive HCM map, the only publicly available knowledge resource dedicated to HCM. Sizes and topological parameters of the networks differed notably, and a low consensus was found in terms of centrality measures between networks. Consensus about the most important nodes was achieved only with respect to one element (calcium). Models with a reduced level of noise were generated and cooperatively working elements were detected. REACH and TRIPS reading systems showed much higher accuracy than Sparser, but at the cost of extraction performance. TRIPS proved to be the best single reading system for text segments about HCM, in terms of the compromise between accuracy and extraction performance.

Conclusions

Different approaches in curation can produce models of the same disease with diverse characteristics, and they give rise to utterly different conclusions in subsequent analysis. The final purpose of the model should direct the choice of curation techniques. Manual curation represents the gold standard for information extraction in biomedical research and is most suitable when only high-quality elements for models are required. Automated curation provides more substance, but high level of noise is expected. Different curation strategies can reduce the level of human input needed. Biomedical knowledge would benefit overwhelmingly, especially as to its rapid growth, if computers were to be able to assist in analysis on a larger scale.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13040-021-00279-2.

AVIA 3.0: interactive portal for genomic variant and sample level analysis

SUMMARY

The Annotation, Visualization and Impact Analysis (AVIA) is a web application combining multiple features to annotate and visualize genomic variant data. Users can investigate functional significance of their genetic alterations across samples, genes and pathways. Version 3.0 of AVIA offers filtering options through interactive charts and by linking disease relevant data sources. Newly incorporated services include gene, variant and sample level reporting, literature and functional correlations among impacted genes, comparative analysis across samples and against data sources such as TCGA and ClinVar, and cohort building. Sample and data management is now feasible through the application, which allows greater flexibility with sharing, reannotating and organizing data. Most importantly, AVIA's utility stems from its convenience for allowing users to upload and explore results without any a priori knowledge or the need to install, update and maintain software or databases. Together, these enhancements strengthen AVIA as a comprehensive, user-driven variant analysis portal.

AVAILABILITYAND IMPLEMENTATION

AVIA is accessible online at https://avia-abcc.ncifcrf.gov.

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

GPCRsignal: webserver for analysis of the interface between G-protein-coupled receptors and their effector proteins by dynamics and mutations

GPCRsignal (https://gpcrsignal.biomodellab.eu/) is a webserver devoted to signaling complexes of G-protein–coupled receptors (GPCRs). The recent improvement in cryo-electron microscopy resulted in the determination of a large number of high-resolution structures of GPCRs bound to their effector proteins: G proteins or arrestins. Analyzing the interfaces between receptor and an effector protein is of high importance since a selection of proper G protein or specific conformation of arrestin leads to changes of signaling that can significantly affect action of drugs. GPCRsignal provides a possibility of running molecular dynamics simulations of all currently available GPCR-effector protein complexes for curated structures: wild-type, with crystal/cryo-EM mutations, or with mutations introduced by the user. The simulations are performed in an implicit water-membrane environment, so they are rather fast. User can run several simulations to obtain statistically valid results. The simulations can be analyzed separately using dynamic FlarePlots for particular types of interactions. One can also compare groups of simulations in Interaction frequency analysis as HeatMaps and also in interaction frequency difference analysis as sticks, linking the interacting residues, of different color and size proportional to differences in contact frequencies.

PredictProtein - Predicting Protein Structure and Function for 29 Years

Since 1992 PredictProtein (https://predictprotein.org) is a one-stop online resource for protein sequence analysis with its main site hosted at the Luxembourg Centre for Systems Biomedicine (LCSB) and queried monthly by over 3,000 users in 2020. PredictProtein was the first Internet server for protein predictions. It pioneered combining evolutionary information and machine learning. Given a protein sequence as input, the server outputs multiple sequence alignments, predictions of protein structure in 1D and 2D (secondary structure, solvent accessibility, transmembrane segments, disordered regions, protein flexibility, and disulfide bridges) and predictions of protein function (functional effects of sequence variation or point mutations, Gene Ontology (GO) terms, subcellular localization, and protein-, RNA-, and DNA binding). PredictProtein's infrastructure has moved to the LCSB increasing throughput; the use of MMseqs2 sequence search reduced runtime five-fold (apparently without lowering performance of prediction methods); user interface elements improved usability, and new prediction methods were added. PredictProtein recently included predictions from deep learning embeddings (GO and secondary structure) and a method for the prediction of proteins and residues binding DNA, RNA, or other proteins. PredictProtein.org aspires to provide reliable predictions to computational and experimental biologists alike. All scripts and methods are freely available for offline execution in high-throughput settings.

Significant non-existence of sequences in genomes and proteomes

Minimal absent words (MAWs) are minimal-length oligomers absent from a genome or proteome. Although some artificially synthesized MAWs have deleterious effects, there is still a lack of a strategy for the classification of non-occurring sequences as potentially malicious or benign. In this work, by using Markovian models with multiple-testing correction, we reveal significant absent oligomers, which are statistically expected to exist. This suggests that their absence is due to negative selection. We survey genomes and proteomes covering the diversity of life and find thousands of significant absent sequences. Common significant MAWs are often mono- or dinucleotide tracts, or palindromic. Significant viral MAWs are often restriction sites and may indicate unknown restriction motifs. Surprisingly, significant mammal genome MAWs are often present, but rare, in other mammals, suggesting that they are suppressed but not completely forbidden. Significant human MAWs are frequently present in prokaryotes, suggesting immune function, but rarely present in human viruses, indicating viral mimicry of the host. More than one-fourth of human proteins are one substitution away from containing a significant MAW, with the majority of replacements being predicted harmful. We provide a web-based, interactive database of significant MAWs across genomes and proteomes.

An Update on MRMAssayDB: A Comprehensive Resource for Targeted Proteomics Assays in the Community

Precise multiplexed quantification of proteins in biological samples can be achieved by targeted proteomics using multiple or parallel reaction monitoring (MRM/PRM). Combined with internal standards, the method achieves very good repeatability and reproducibility enabling excellent protein quantification and allowing longitudinal and cohort studies. A laborious part of performing such experiments lies in the preparation steps dedicated to the development and validation of individual protein assays. Several public repositories host information on targeted proteomics assays, including NCI’s Clinical Proteomic Tumor Analysis Consortium assay portals, PeptideAtlas SRM Experiment Library, SRMAtlas, PanoramaWeb, and PeptideTracker, with all offering varying levels of details. We introduced MRMAssayDB in 2018 as an integrated resource for targeted proteomics assays. The Web-based application maps and links the assays from the repositories, includes comprehensive up-to-date protein and sequence annotations, and provides multiple visualization options on the peptide and protein level. We have extended MRMAssayDB with more assays and extensive annotations. Currently it contains >828 000 assays covering >51 000 proteins from 94 organisms, of which >17 000 proteins are present in >2400 biological pathways, and >48 000 mapping to >21 000 Gene Ontology terms. This is an increase of about four times the number of assays since introduction. We have expanded annotations of interaction, biological pathways, and disease associations. A newly added visualization module for coupled molecular structural annotation browsing allows the user to interactively examine peptide sequence and any known PTMs and disease mutations, and map all to available protein 3D structures. Because of its integrative approach, MRMAssayDB enables a holistic view of suitable proteotypic peptides and commonly used transitions in empirical data. Availability: http://mrmassaydb.proteincentre.com.

Visualizing protein structures - tools and trends

Molecular visualization is fundamental in the current scientific literature, textbooks and dissemination materials. It provides an essential support for presenting results, reasoning on and formulating hypotheses related to molecular structure. Tools for visual exploration of structural data have become easily accessible on a broad variety of platforms thanks to advanced software tools that render a great service to the scientific community. These tools are often developed across disciplines bridging computer science, biology and chemistry. This mini-review was written as a short and compact overview for scientists who need to visualize protein structures and want to make an informed decision which tool they should use. Here, we first describe a few 'Swiss Army knives' geared towards protein visualization for everyday use with an existing large user base, then focus on more specialized tools for peculiar needs that are not yet as broadly known. Our selection is by no means exhaustive, but reflects a diverse snapshot of scenarios that we consider informative for the reader. We end with an account of future trends and perspectives.

MISCAST: MIssense variant to protein StruCture Analysis web SuiTe

Human genome sequencing efforts have greatly expanded, and a plethora of missense variants identified both in patients and in the general population is now publicly accessible. Interpretation of the molecular-level effect of missense variants, however, remains challenging and requires a particular investigation of amino acid substitutions in the context of protein structure and function. Answers to questions like 'Is a variant perturbing a site involved in key macromolecular interactions and/or cellular signaling?', or 'Is a variant changing an amino acid located at the protein core or part of a cluster of known pathogenic mutations in 3D?' are crucial. Motivated by these needs, we developed MISCAST (missense variant to protein structure analysis web suite; http://miscast.broadinstitute.org/). MISCAST is an interactive and user-friendly web server to visualize and analyze missense variants in protein sequence and structure space. Additionally, a comprehensive set of protein structural and functional features have been aggregated in MISCAST from multiple databases, and displayed on structures alongside the variants to provide users with the biological context of the variant location in an integrated platform. We further made the annotated data and protein structures readily downloadable from MISCAST to foster advanced offline analysis of missense variants by a wide biological community.

The Atlas of Inflammation Resolution (AIR)

Acute inflammation is a protective reaction by the immune system in response to invading pathogens or tissue damage. Ideally, the response should be localized, self-limited, and returning to homeostasis. If not resolved, acute inflammation can result in organ pathologies leading to chronic inflammatory phenotypes. Acute inflammation and inflammation resolution are complex coordinated processes, involving a number of cell types, interacting in space and time. The biomolecular complexity and the fact that several biomedical fields are involved, make a multi- and interdisciplinary approach necessary. The Atlas of Inflammation Resolution (AIR) is a web-based resource capturing an essential part of the state-of-the-art in acute inflammation and inflammation resolution research. The AIR provides an interface for users to search thousands of interactions, arranged in inter-connected multi-layers of process diagrams, covering a wide range of clinically relevant phenotypes. By mapping experimental data onto the AIR, it can be used to elucidate drug action as well as molecular mechanisms underlying different disease phenotypes. For the visualization and exploration of information, the AIR uses the Minerva platform, which is a well-established tool for the presentation of disease maps. The molecular details of the AIR are encoded using international standards. The AIR was created as a freely accessible resource, supporting research and education in the fields of acute inflammation and inflammation resolution. The AIR connects research communities, facilitates clinical decision making, and supports research scientists in the formulation and validation of hypotheses. The AIR is accessible through https://air.bio.informatik.uni-rostock.de.

WITHDRAWN: The Atlas of Inflammation Resolution (AIR)

The Publisher regrets that this article is an accidental duplication of an article that has already been published, https://doi.org/10.1016/j.mam.2020.100894. The duplicate article has therefore been withdrawn. The full Elsevier Policy on Article Withdrawal can be found at https://www.elsevier.com/about/our-business/policies/article-withdrawal.

Closing the gap between formats for storing layout information in systems biology

The understanding of complex biological networks often relies on both a dedicated layout and a topology. Currently, there are three major competing layout-aware systems biology formats, but there are no software tools or software libraries supporting all of them. This complicates the management of molecular network layouts and hinders their reuse and extension. In this paper, we present a high-level overview of the layout formats in systems biology, focusing on their commonalities and differences, review their support in existing software tools, libraries and repositories and finally introduce a new conversion module within the MINERVA platform. The module is available via a REST API and offers, besides the ability to convert between layout-aware systems biology formats, the possibility to export layouts into several graphical formats. The module enables conversion of very large networks with thousands of elements, such as disease maps or metabolic reconstructions, rendering it widely applicable in systems biology.