Universal principles of lineage architecture and stem cell identity in renewing tissues.

Adult tissues in multicellular organisms typically contain a variety of stem, progenitor and differentiated cell types arranged in a lineage hierarchy that regulates healthy tissue turnover. Lineage hierarchies in disparate tissues often exhibit common features, yet the general principles regulating their architecture are not known. Here, we provide a formal framework for understanding the relationship between cell molecular 'states' and cell 'types', based on the topology of admissible cell state trajectories. We show that a self-renewing cell type - if defined as suggested by this framework - must reside at the top of any homeostatic renewing lineage hierarchy, and only there. This architecture arises as a natural consequence of homeostasis, and indeed is the only possible way that lineage architectures can be constructed to support homeostasis in renewing tissues. Furthermore, under suitable feedback regulation, for example from the stem cell niche, we show that the property of 'stemness' is entirely determined by the cell environment, in accordance with the notion that stem cell identities are contextual and not determined by hard-wired, cell-intrinsic characteristics. This article has an associated 'The people behind the papers' interview.

Stratification of asthma by lipidomic profiling of induced sputum supernatant.

BACKGROUND: Asthma is a chronic respiratory disease with significant heterogeneity in its clinical presentation and pathobiology. There is need for improved understanding of respiratory lipid metabolism in asthma patients and its relation to observable clinical features. OBJECTIVE: We performed a comprehensive, prospective, cross-sectional analysis of the lipid composition of induced sputum supernatant obtained from asthma patients with a range of disease severities, as well as from healthy controls. METHODS: Induced sputum supernatant was collected from 211 adults with asthma and 41 healthy individuals enrolled onto the U-BIOPRED (Unbiased Biomarkers for the Prediction of Respiratory Disease Outcomes) study. Sputum lipidomes were characterized by semiquantitative shotgun mass spectrometry and clustered using topologic data analysis to identify lipid phenotypes. RESULTS: Shotgun lipidomics of induced sputum supernatant revealed a spectrum of 9 molecular phenotypes, highlighting not just significant differences between the sputum lipidomes of asthma patients and healthy controls, but also within the asthma patient population. Matching clinical, pathobiologic, proteomic, and transcriptomic data helped inform the underlying disease processes. Sputum lipid phenotypes with higher levels of nonendogenous, cell-derived lipids were associated with significantly worse asthma severity, worse lung function, and elevated granulocyte counts. CONCLUSION: We propose a novel mechanism of increased lipid loading in the epithelial lining fluid of asthma patients resulting from the secretion of extracellular vesicles by granulocytic inflammatory cells, which could reduce the ability of pulmonary surfactant to lower surface tension in asthmatic small airways, as well as compromise its role as an immune regulator.

Cryo-EM Map Anisotropy Can Be Attenuated by Map Post-Processing and a New Method for Its Estimation.

One of the most important challenges in cryogenic electron microscopy (cryo-EM) is the substantial number of samples that exhibit preferred orientations, which leads to an uneven coverage of the projection sphere. As a result, the overall quality of the reconstructed maps can be severely affected, as manifested by the presence of anisotropy in the map resolution. Several methods have been proposed to measure the directional resolution of maps in tandem with experimental protocols to address the problem of preferential orientations in cryo-EM. Following these works, in this manuscript we identified one potential limitation that may affect most of the existing methods and we proposed an alternative approach to evaluate the presence of preferential orientations in cryo-EM reconstructions. In addition, we also showed that some of the most recently proposed cryo-EM map post-processing algorithms can attenuate map anisotropy, thus offering alternative visualization opportunities for cases affected by moderate levels of preferential orientations.

An information-theoretic approach to single cell sequencing analysis.

BACKGROUND: Single-cell sequencing (sc-Seq) experiments are producing increasingly large data sets. However, large data sets do not necessarily contain large amounts of information. RESULTS: Here, we formally quantify the information obtained from a sc-Seq experiment and show that it corresponds to an intuitive notion of gene expression heterogeneity. We demonstrate a natural relation between our notion of heterogeneity and that of cell type, decomposing heterogeneity into that component attributable to differential expression between cell types (inter-cluster heterogeneity) and that remaining (intra-cluster heterogeneity). We test our definition of heterogeneity as the objective function of a clustering algorithm, and show that it is a useful descriptor for gene expression patterns associated with different cell types. CONCLUSIONS: Thus, our definition of gene heterogeneity leads to a biologically meaningful notion of cell type, as groups of cells that are statistically equivalent with respect to their patterns of gene expression. Our measure of heterogeneity, and its decomposition into inter- and intra-cluster, is non-parametric, intrinsic, unbiased, and requires no additional assumptions about expression patterns. Based on this theory, we develop an efficient method for the automatic unsupervised clustering of cells from sc-Seq data, and provide an R package implementation.

Fragment Merging Using a Graph Database Samples Different Catalogue Space than Similarity Search.

Fragment merging is a promising approach to progressing fragments directly to on-scale potency: each designed compound incorporates the structural motifs of overlapping fragments in a way that ensures compounds recapitulate multiple high-quality interactions. Searching commercial catalogues provides one useful way to quickly and cheaply identify such merges and circumvents the challenge of synthetic accessibility, provided they can be readily identified. Here, we demonstrate that the Fragment Network, a graph database that provides a novel way to explore the chemical space surrounding fragment hits, is well-suited to this challenge. We use an iteration of the database containing >120 million catalogue compounds to find fragment merges for four crystallographic screening campaigns and contrast the results with a traditional fingerprint-based similarity search. The two approaches identify complementary sets of merges that recapitulate the observed fragment-protein interactions but lie in different regions of chemical space. We further show our methodology is an effective route to achieving on-scale potency by retrospective analyses for two different targets; in analyses of public COVID Moonshot and Mycobacterium tuberculosis EthR inhibitors, potential inhibitors with micromolar IC50 values were identified. This work demonstrates the use of the Fragment Network to increase the yield of fragment merges beyond that of a classical catalogue search.

3DBionotes COVID-19 edition.

SUMMARY: The web platform 3DBionotes-WS integrates multiple web services and an interactive web viewer to provide a unified environment in which biological annotations can be analyzed in their structural context. Since the COVID-19 outbreak, new structural data from many viral proteins have been provided at a very fast pace. This effort includes many cryogenic electron microscopy (cryo-EM) studies, together with more traditional ones (X-rays, NMR), using several modeling approaches and complemented with structural predictions. At the same time, a plethora of new genomics and interactomics information (including fragment screening and structure-based virtual screening efforts) have been made available from different servers. In this context, we have developed 3DBionotes-COVID-19 as an answer to: (i) the need to explore multiomics data in a unified context with a special focus on structural information and (ii) the drive to incorporate quality measurements, especially in the form of advanced validation metrics for cryo-EM. AVAILABILITY AND IMPLEMENTATION: https://3dbionotes.cnb.csic.es/ws/covid19. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Advances in Xmipp for Cryo-Electron Microscopy: From Xmipp to Scipion.

Xmipp is an open-source software package consisting of multiple programs for processing data originating from electron microscopy and electron tomography, designed and managed by the Biocomputing Unit of the Spanish National Center for Biotechnology, although with contributions from many other developers over the world. During its 25 years of existence, Xmipp underwent multiple changes and updates. While there were many publications related to new programs and functionality added to Xmipp, there is no single publication on the Xmipp as a package since 2013. In this article, we give an overview of the changes and new work since 2013, describe technologies and techniques used during the development, and take a peek at the future of the package.

MicrographCleaner: A python package for cryo-EM micrograph cleaning using deep learning.

Cryo-EM Single Particle Analysis workflows require tens of thousands of high-quality particle projections to unveil the three-dimensional structure of macromolecules. Conventional methods for automatic particle picking tend to suffer from high false-positive rates, hampering the reconstruction process. One common cause of this problem is the presence of carbon and different types of high-contrast contaminations. In order to overcome this limitation, we have developed MicrographCleaner, a deep learning package designed to discriminate, in an automated fashion, between regions of micrographs which are suitable for particle picking, and those which are not. MicrographCleaner implements a U-net-like deep learning model trained on a manually curated dataset compiled from over five hundred micrographs. The benchmarking, carried out on approximately one hundred independent micrographs, shows that MicrographCleaner is a very efficient approach for micrograph preprocessing. MicrographCleaner (micrograph_cleaner_em) package is available at PyPI and Anaconda Cloud and also as a Scipion/Xmipp protocol. Source code is available at https://github.com/rsanchezgarc/micrograph_cleaner_em.

3DBIONOTES v3.0: crossing molecular and structural biology data with genomic variations.

MOTIVATION: Many diseases are associated to single nucleotide polymorphisms that affect critical regions of proteins as binding sites or post translational modifications. Therefore, analysing genomic variants with structural and molecular biology data is a powerful framework in order to elucidate the potential causes of such diseases. RESULTS: A new version of our web framework 3DBIONOTES is presented. This version offers new tools to analyse and visualize protein annotations and genomic variants, including a contingency analysis of variants and amino acid features by means of a Fisher exact test, the integration of a gene annotation viewer to highlight protein features on gene sequences and a protein-protein interaction viewer to display protein annotations at network level. AVAILABILITY AND IMPLEMENTATION: The web server is available at https://3dbionotes.cnb.csic.es. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. CONTACT: Spanish National Institute for Bioinformatics (INB ELIXIR-ES) and Biocomputing Unit, National Centre of Biotechnology (CSIC)/Instruct Image Processing Centre, C/ Darwin nº 3, Campus of Cantoblanco, 28049 Madrid, Spain.

BIPSPI: a method for the prediction of partner-specific protein-protein interfaces.

Motivation: Protein-Protein Interactions (PPI) are essentials for most cellular processes and thus, unveiling how proteins interact is a crucial question that can be better understood by identifying which residues are responsible for the interaction. Computational approaches are orders of magnitude cheaper and faster than experimental ones, leading to proliferation of multiple methods aimed to predict which residues belong to the interface of an interaction. Results: We present BIPSPI, a new machine learning-based method for the prediction of partner-specific PPI sites. Contrary to most binding site prediction methods, the proposed approach takes into account a pair of interacting proteins rather than a single one in order to predict partner-specific binding sites. BIPSPI has been trained employing sequence-based and structural features from both protein partners of each complex compiled in the Protein-Protein Docking Benchmark version 5.0 and in an additional set independently compiled. Also, a version trained only on sequences has been developed. The performance of our approach has been assessed by a leave-one-out cross-validation over different benchmarks, outperforming state-of-the-art methods. Availability and implementation: BIPSPI web server is freely available at http://bipspi.cnb.csic.es. BIPSPI code is available at https://github.com/bioinsilico/BIPSPI. Docker image is available at https://hub.docker.com/r/bioinsilico/bipspi/. Supplementary information: Supplementary data are available at Bioinformatics online.

Validation of electron microscopy initial models via small angle X-ray scattering curves.

MOTIVATION: Cryo electron microscopy (EM) is currently one of the main tools to reveal the structural information of biological macromolecules. The re-construction of three-dimensional (3D) maps is typically carried out following an iterative process that requires an initial estimation of the 3D map to be refined in subsequent steps. Therefore, its determination is key in the quality of the final results, and there are cases in which it is still an open issue in single particle analysis (SPA). Small angle X-ray scattering (SAXS) is a well-known technique applied to structural biology. It is useful from small nanostructures up to macromolecular ensembles for its ability to obtain low resolution information of the biological sample measuring its X-ray scattering curve. These curves, together with further analysis, are able to yield information on the sizes, shapes and structures of the analyzed particles. RESULTS: In this paper, we show how the low resolution structural information revealed by SAXS is very useful for the validation of EM initial 3D models in SPA, helping the following refinement process to obtain more accurate 3D structures. For this purpose, we approximate the initial map by pseudo-atoms and predict the SAXS curve expected for this pseudo-atomic structure. The match between the predicted and experimental SAXS curves is considered as a good sign of the correctness of the EM initial map. AVAILABILITY AND IMPLEMENTATION: The algorithm is freely available as part of the Scipion 1.2 software at http://scipion.i2pc.es/.

3DBIONOTES v2.0: a web server for the automatic annotation of macromolecular structures.

MOTIVATION: Complementing structural information with biochemical and biomedical annotations is a powerful approach to explore the biological function of macromolecular complexes. However, currently the compilation of annotations and structural data is a feature only available for those structures that have been released as entries to the Protein Data Bank. RESULTS: To help researchers in assessing the consistency between structures and biological annotations for structural models not deposited in databases, we present 3DBIONOTES v2.0, a web application designed for the automatic annotation of biochemical and biomedical information onto macromolecular structural models determined by any experimental or computational technique. AVAILABILITY AND IMPLEMENTATION: The web server is available at http://3dbionotes-ws.cnb.csic.es. CONTACT: jsegura@cnb.csic.es. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Enantioselective Synthesis of 2-Amino-1,1-diarylalkanes Bearing a Carbocyclic Ring Substituted Indole through Asymmetric Catalytic Reaction of Hydroxyindoles with Nitroalkenes.

An asymmetric catalytic reaction of hydroxyindoles with nitroalkenes leading to the Friedel-Crafts alkylation in the carbocyclic ring of indole is presented. The method is based on the activating/directing effects of the hydroxy group situated in the carbocyclic ring of the indole providing nitroalkylated indoles functionalizated at the C-4, C-5, and C-7 positions with high yield, regio-, and enantioselectivity. The optically enriched nitroalkanes were transformed efficiently in optically enriched 2-amino-1,1-diarylalkanes bearing a carbocyclic ring substituted indole.

3DCONS-DB: A Database of Position-Specific Scoring Matrices in Protein Structures.

Many studies have used position-specific scoring matrices (PSSM) profiles to characterize residues in protein structures and to predict a broad range of protein features. Moreover, PSSM profiles of Protein Data Bank (PDB) entries have been recalculated in many works for different purposes. Although the computational cost of calculating a single PSSM profile is affordable, many statistical studies or machine learning-based methods used thousands of profiles to achieve their goals, thereby leading to a substantial increase of the computational cost. In this work we present a new database compiling PSSM profiles for the proteins of the PDB. Currently, the database contains 333,532 protein chain profiles involving 123,135 different PDB entries.

3DIANA: 3D Domain Interaction Analysis: A Toolbox for Quaternary Structure Modeling.

Electron microscopy (EM) is experiencing a revolution with the advent of a new generation of Direct Electron Detectors, enabling a broad range of large and flexible structures to be resolved well below 1 nm resolution. Although EM techniques are evolving to the point of directly obtaining structural data at near-atomic resolution, for many molecules the attainable resolution might not be enough to propose high-resolution structural models. However, accessing information on atomic coordinates is a necessary step toward a deeper understanding of the molecular mechanisms that allow proteins to perform specific tasks. For that reason, methods for the integration of EM three-dimensional maps with x-ray and NMR structural data are being developed, a modeling task that is normally referred to as fitting, resulting in the so called hybrid models. In this work, we present a novel application-3DIANA-specially targeted to those cases in which the EM map resolution is medium or low and additional experimental structural information is scarce or even lacking. In this way, 3DIANA statistically evaluates proposed/potential contacts between protein domains, presents a complete catalog of both structurally resolved and predicted interacting regions involving these domains and, finally, suggests structural templates to model the interaction between them. The evaluation of the proposed interactions is computed with DIMERO, a new method that scores physical binding sites based on the topology of protein interaction networks, which has recently shown the capability to increase by 200% the number of domain-domain interactions predicted in interactomes as compared to previous approaches. The new application displays the information at a sequence and structural level and is accessible through a web browser or as a Chimera plugin at http://3diana.cnb.csic.es.

The transition to synchronization of networked systems.

We study the synchronization properties of a generic networked dynamical system, and show that, under a suitable approximation, the transition to synchronization can be predicted with the only help of eigenvalues and eigenvectors of the graph Laplacian matrix. The transition comes out to be made of a well defined sequence of events, each of which corresponds to a specific clustered state. The network's nodes involved in each of the clusters can be identified, and the value of the coupling strength at which the events are taking place can be approximately ascertained. Finally, we present large-scale simulations which show the accuracy of the approximation made, and of our predictions in describing the synchronization transition of both synthetic and real-world large size networks, and we even report that the observed sequence of clusters is preserved in heterogeneous networks made of slightly non-identical systems.

MetaFun: unveiling sex-based differences in multiple transcriptomic studies through comprehensive functional meta-analysis.

BACKGROUND: While sex-based differences in various health scenarios have been thoroughly acknowledged in the literature, we lack sufficient tools and methods that allow for an in-depth analysis of sex as a variable in biomedical research. To fill this knowledge gap, we created MetaFun as an easy-to-use web-based tool to meta-analyze multiple transcriptomic datasets with a sex-based perspective to gain major statistical power and biological soundness. DESCRIPTION: MetaFun is a complete suite that allows the analysis of transcriptomics data and the exploration of the results at all levels, performing single-dataset exploratory analysis, differential gene expression, gene set functional enrichment, and finally, combining results in a functional meta-analysis. Which biological processes, molecular functions or cellular components are altered in a common pattern in different transcriptomic studies when comparing male and female patients? This and other biological questions of interest can be answered with the use of MetaFun. This tool is available at https://bioinfo.cipf.es/metafun while additional help can be found at https://gitlab.com/ubb-cipf/metafunweb/-/wikis/Summary . CONCLUSIONS: Overall, Metafun is the first open-access web-based tool to identify consensus biological functions across multiple transcriptomic datasets, helping to elucidate sex differences in numerous diseases. Its use will facilitate the generation of novel biological knowledge that can be used in the research and application of Personalized Medicine considering the sex of patients.

Persistent homology analysis distinguishes pathological bone microstructure in non-linear microscopy images.

We present a topological method for the detection and quantification of bone microstructure from non-linear microscopy images. Specifically, we analyse second harmonic generation (SHG) and two photon excited autofluorescence (TPaF) images of bone tissue which capture the distribution of matrix (fibrillar collagen) structure and autofluorescent molecules, respectively. Using persistent homology statistics with a signed Euclidean distance transform filtration on binary patches of images, we are able to quantify the number, size, distribution, and crowding of holes within and across samples imaged at the microscale. We apply our methodology to a previously characterized murine model of skeletal pathology whereby vascular endothelial growth factor expression was deleted in osteocalcin-expressing cells (OcnVEGFKO) presenting increased cortical porosity, compared to wild type (WT) littermate controls. We show significant differences in topological statistics between the OcnVEGFKO and WT groups and, when classifying the males, or females respectively, into OcnVEGFKO or WT groups, we obtain high prediction accuracies of 98.7% (74.2%) and 77.8% (65.8%) respectively for SHG (TPaF) images. The persistence statistics that we use are fully interpretable, can highlight regions of abnormality within an image and identify features at different spatial scales.

Multidisciplinary ecosystem to study lifecourse determinants and prevention of early-onset burdensome multimorbidity (MELD-B) - protocol for a research collaboration.

Background: Most people living with multiple long-term condition multimorbidity (MLTC-M) are under 65 (defined as 'early onset'). Earlier and greater accrual of long-term conditions (LTCs) may be influenced by the timing and nature of exposure to key risk factors, wider determinants or other LTCs at different life stages. We have established a research collaboration titled 'MELD-B' to understand how wider determinants, sentinel conditions (the first LTC in the lifecourse) and LTC accrual sequence affect risk of early-onset, burdensome MLTC-M, and to inform prevention interventions. Aim: Our aim is to identify critical periods in the lifecourse for prevention of early-onset, burdensome MLTC-M, identified through the analysis of birth cohorts and electronic health records, including artificial intelligence (AI)-enhanced analyses. Design: We will develop deeper understanding of 'burdensomeness' and 'complexity' through a qualitative evidence synthesis and a consensus study. Using safe data environments for analyses across large, representative routine healthcare datasets and birth cohorts, we will apply AI methods to identify early-onset, burdensome MLTC-M clusters and sentinel conditions, develop semi-supervised learning to match individuals across datasets, identify determinants of burdensome clusters, and model trajectories of LTC and burden accrual. We will characterise early-life (under 18 years) risk factors for early-onset, burdensome MLTC-M and sentinel conditions. Finally, using AI and causal inference modelling, we will model potential 'preventable moments', defined as time periods in the life course where there is an opportunity for intervention on risk factors and early determinants to prevent the development of MLTC-M. Patient and public involvement is integrated throughout.

Geometry and symmetry in biochemical reaction systems.

Complex systems of intracellular biochemical reactions have a central role in regulating cell identities and functions. Biochemical reaction systems are typically studied using the language and tools of graph theory. However, graph representations only describe pairwise interactions between molecular species and so are not well suited to modelling complex sets of reactions that may involve numerous reactants and/or products. Here, we make use of a recently developed hypergraph theory of chemical reactions that naturally allows for higher-order interactions to explore the geometry and quantify functional redundancy in biochemical reactions systems. Our results constitute a general theory of automorphisms for oriented hypergraphs and describe the effect of automorphism group structure on hypergraph Laplacian spectra.