Opportunities and Challenges in Building a Spatiotemporal Multi-scale Model of the Human Pancreatic ß Cell.

The construction of a predictive model of an entire eukaryotic cell that describes its dynamic structure from atomic to cellular scales is a grand challenge at the intersection of biology, chemistry, physics, and computer science. Having such a model will open new dimensions in biological research and accelerate healthcare advancements. Developing the necessary experimental and modeling methods presents abundant opportunities for a community effort to realize this goal. Here, we present a vision for creation of a spatiotemporal multi-scale model of the pancreatic ß-cell, a relevant target for understanding and modulating the pathogenesis of diabetes.

A Regulatory Circuit Controlling the Dynamics of NFκB cRel Transitions B Cells from Proliferation to Plasma Cell Differentiation.

Humoral immunity depends on efficient activation of B cells and their subsequent differentiation into antibody-secreting cells (ASCs). The transcription factor NFκB cRel is critical for B cell proliferation, but incorporating its known regulatory interactions into a mathematical model of the ASC differentiation circuit prevented ASC generation in simulations. Indeed, experimental ectopic cRel expression blocked ASC differentiation by inhibiting the transcription factor Blimp1, and in wild-type (WT) cells cRel was dynamically repressed during ASC differentiation by Blimp1 binding the Rel locus. Including this bi-stable circuit of mutual cRel-Blimp1 antagonism into a multi-scale model revealed that dynamic repression of cRel controls the switch from B cell proliferation to ASC generation phases and hence the respective cell population dynamics. Our studies provide a mechanistic explanation of how dysregulation of this bi-stable circuit might result in pathologic B cell population phenotypes and thus offer new avenues for diagnostic stratification and treatment.

scAMAC: self-supervised clustering of scRNA-seq data based on adaptive multi-scale autoencoder.

Cluster assignment is vital to analyzing single-cell RNA sequencing (scRNA-seq) data to understand high-level biological processes. Deep learning-based clustering methods have recently been widely used in scRNA-seq data analysis. However, existing deep models often overlook the interconnections and interactions among network layers, leading to the loss of structural information within the network layers. Herein, we develop a new self-supervised clustering method based on an adaptive multi-scale autoencoder, called scAMAC. The self-supervised clustering network utilizes the Multi-Scale Attention mechanism to fuse the feature information from the encoder, hidden and decoder layers of the multi-scale autoencoder, which enables the exploration of cellular correlations within the same scale and captures deep features across different scales. The self-supervised clustering network calculates the membership matrix using the fused latent features and optimizes the clustering network based on the membership matrix. scAMAC employs an adaptive feedback mechanism to supervise the parameter updates of the multi-scale autoencoder, obtaining a more effective representation of cell features. scAMAC not only enables cell clustering but also performs data reconstruction through the decoding layer. Through extensive experiments, we demonstrate that scAMAC is superior to several advanced clustering and imputation methods in both data clustering and reconstruction. In addition, scAMAC is beneficial for downstream analysis, such as cell trajectory inference. Our scAMAC model codes are freely available at https://github.com/yancy2024/scAMAC.

MUSCLE: multi-view and multi-scale attentional feature fusion for microRNA-disease associations prediction.

MicroRNAs (miRNAs) synergize with various biomolecules in human cells resulting in diverse functions in regulating a wide range of biological processes. Predicting potential disease-associated miRNAs as valuable biomarkers contributes to the treatment of human diseases. However, few previous methods take a holistic perspective and only concentrate on isolated miRNA and disease objects, thereby ignoring that human cells are responsible for multiple relationships. In this work, we first constructed a multi-view graph based on the relationships between miRNAs and various biomolecules, and then utilized graph attention neural network to learn the graph topology features of miRNAs and diseases for each view. Next, we added an attention mechanism again, and developed a multi-scale feature fusion module, aiming to determine the optimal fusion results for the multi-view topology features of miRNAs and diseases. In addition, the prior attribute knowledge of miRNAs and diseases was simultaneously added to achieve better prediction results and solve the cold start problem. Finally, the learned miRNA and disease representations were then concatenated and fed into a multi-layer perceptron for end-to-end training and predicting potential miRNA-disease associations. To assess the efficacy of our model (called MUSCLE), we performed 5- and 10-fold cross-validation (CV), which got average the Area under ROC curves of 0.966${\pm }$0.0102 and 0.973${\pm }$0.0135, respectively, outperforming most current state-of-the-art models. We then examined the impact of crucial parameters on prediction performance and performed ablation experiments on the feature combination and model architecture. Furthermore, the case studies about colon cancer, lung cancer and breast cancer also fully demonstrate the good inductive capability of MUSCLE. Our data and code are free available at a public GitHub repository: https://github.com/zht-code/MUSCLE.git.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane.

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. Therefore, development of new anti-TB drugs remains one of the highest public health priorities. Mtb has evolved a complex cell envelope that represents a formidable barrier to antibiotics. The Mtb cell envelop consists of four distinct layers enriched for Mtb specific lipids and glycans. Although the outer membrane, comprised of mycolic acid esters, has been extensively studied, less is known about the plasma membrane, which also plays a critical role in impacting antibiotic efficacy. The Mtb plasma membrane has a unique lipid composition, with mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs) comprising more than 50% of the lipids. However, the role of PIMs in the structure and function of the membrane remains elusive. Here, we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to shape the biophysical properties of mycobacterial plasma membranes. We assess both symmetric and asymmetric assemblies of the Mtb plasma membrane and compare this with residue distributions of Mtb integral membrane protein structures. To further validate the model, we tested known anti-TB drugs and demonstrated that our models agree with experimental results. Thus, our work sheds new light on the organization of the mycobacterial plasma membrane. This paves the way for future studies on antibiotic development and understanding Mtb membrane protein function.

Phenotypic plasticity: A missing element in the theory of vegetation pattern formation.

Regular spatial patterns of vegetation are a common sight in drylands. Their formation is a population-level response to water stress that increases water availability for the few via partial plant mortality. At the individual level, plants can also adapt to water stress by changing their phenotype. Phenotypic plasticity of individual plants and spatial patterning of plant populations have extensively been studied independently, but the likely interplay between the two robust mechanisms has remained unexplored. In this paper, we incorporate phenotypic plasticity into a multi-level theory of vegetation pattern formation and use a fascinating ecological phenomenon, the Namibian "fairy circles," to demonstrate the need for such a theory. We show that phenotypic changes in the root structure of plants, coupled with pattern-forming feedback within soil layers, can resolve two puzzles that the current theory fails to explain: observations of multi-scale patterns and the absence of theoretically predicted large-scale stripe and spot patterns along the rainfall gradient. Importantly, we find that multi-level responses to stress unveil a wide variety of more effective stress-relaxation pathways, compared to single-level responses, implying a previously underestimated resilience of dryland ecosystems.

Airborne disease transmission during indoor gatherings over multiple time scales: Modeling framework and policy implications.

Indoor superspreading events are significant drivers of transmission of respiratory diseases. In this work, we study the dynamics of airborne transmission in consecutive meetings of individuals in enclosed spaces. In contrast to the usual pairwise-interaction models of infection where effective contacts transmit the disease, we focus on group interactions where individuals with distinct health states meet simultaneously. Specifically, the disease is transmitted by infected individuals exhaling droplets (contributing to the viral load in the closed space) and susceptible ones inhaling the contaminated air. We propose a modeling framework that couples the fast dynamics of the viral load attained over meetings in enclosed spaces and the slow dynamics of disease progression at the population level. Our modeling framework incorporates the multiple time scales involved in different setups in which indoor events may happen, from single-time events to events hosting multiple meetings per day, over many days. We present theoretical and numerical results of trade-offs between the room characteristics (ventilation system efficiency and air mass) and the group's behavioral and composition characteristics (group size, mask compliance, testing, meeting time, and break times), that inform indoor policies to achieve disease control in closed environments through different pathways. Our results emphasize the impact of break times, mask-wearing, and testing on facilitating the conditions to achieve disease control. We study scenarios of different break times, mask compliance, and testing. We also derive policy guidelines to contain the infection rate under a certain threshold.

MARPPI: boosting prediction of protein-protein interactions with multi-scale architecture residual network.

Protein-protein interactions (PPIs) are a major component of the cellular biochemical reaction network. Rich sequence information and machine learning techniques reduce the dependence of exploring PPIs on wet experiments, which are costly and time-consuming. This paper proposes a PPI prediction model, multi-scale architecture residual network for PPIs (MARPPI), based on dual-channel and multi-feature. Multi-feature leverages Res2vec to obtain the association information between residues, and utilizes pseudo amino acid composition, autocorrelation descriptors and multivariate mutual information to achieve the amino acid composition and order information, physicochemical properties and information entropy, respectively. Dual channel utilizes multi-scale architecture improved ResNet network which extracts protein sequence features to reduce protein feature loss. Compared with other advanced methods, MARPPI achieves 96.03%, 99.01% and 91.80% accuracy in the intraspecific datasets of Saccharomyces cerevisiae, Human and Helicobacter pylori, respectively. The accuracy on the two interspecific datasets of Human-Bacillus anthracis and Human-Yersinia pestis is 97.29%, and 95.30%, respectively. In addition, results on specific datasets of disease (neurodegenerative and metabolic disorders) demonstrate the ability to detect hidden interactions. To better illustrate the performance of MARPPI, evaluations on independent datasets and PPIs network suggest that MARPPI can be used to predict cross-species interactions. The above shows that MARPPI can be regarded as a concise, efficient and accurate tool for PPI datasets.

MMSMAPlus: a multi-view multi-scale multi-attention embedding model for protein function prediction.

Protein is the most important component in organisms and plays an indispensable role in life activities. In recent years, a large number of intelligent methods have been proposed to predict protein function. These methods obtain different types of protein information, including sequence, structure and interaction network. Among them, protein sequences have gained significant attention where methods are investigated to extract the information from different views of features. However, how to fully exploit the views for effective protein sequence analysis remains a challenge. In this regard, we propose a multi-view, multi-scale and multi-attention deep neural model (MMSMA) for protein function prediction. First, MMSMA extracts multi-view features from protein sequences, including one-hot encoding features, evolutionary information features, deep semantic features and overlapping property features based on physiochemistry. Second, a specific multi-scale multi-attention deep network model (MSMA) is built for each view to realize the deep feature learning and preliminary classification. In MSMA, both multi-scale local patterns and long-range dependence from protein sequences can be captured. Third, a multi-view adaptive decision mechanism is developed to make a comprehensive decision based on the classification results of all the views. To further improve the prediction performance, an extended version of MMSMA, MMSMAPlus, is proposed to integrate homology-based protein prediction under the framework of multi-view deep neural model. Experimental results show that the MMSMAPlus has promising performance and is significantly superior to the state-of-the-art methods. The source code can be found at https://github.com/wzy-2020/MMSMAPlus.

Exploring pharmacological active ingredients of traditional Chinese medicine by pharmacotranscriptomic map in ITCM.

With the emergence of high-throughput technologies, computational screening based on gene expression profiles has become one of the most effective methods for drug discovery. More importantly, profile-based approaches remarkably enhance novel drug-disease pair discovery without relying on drug- or disease-specific prior knowledge, which has been widely used in modern medicine. However, profile-based systematic screening of active ingredients of traditional Chinese medicine (TCM) has been scarcely performed due to inadequate pharmacotranscriptomic data. Here, we develop the largest-to-date online TCM active ingredients-based pharmacotranscriptomic platform integrated traditional Chinese medicine (ITCM) for the effective screening of active ingredients. First, we performed unified high-throughput experiments and constructed the largest data repository of 496 representative active ingredients, which was five times larger than the previous one built by our team. The transcriptome-based multi-scale analysis was also performed to elucidate their mechanism. Then, we developed six state-of-art signature search methods to screen active ingredients and determine the optimal signature size for all methods. Moreover, we integrated them into a screening strategy, TCM-Query, to identify the potential active ingredients for the special disease. In addition, we also comprehensively collected the TCM-related resource by literature mining. Finally, we applied ITCM to an active ingredient bavachinin, and two diseases, including prostate cancer and COVID-19, to demonstrate the power of drug discovery. ITCM was aimed to comprehensively explore the active ingredients of TCM and boost studies of pharmacological action and drug discovery. ITCM is available at http://itcm.biotcm.net.

Detection of PatIent-Level distances from single cell genomics and pathomics data with Optimal Transport (PILOT).

Although clinical applications represent the next challenge in single-cell genomics and digital pathology, we still lack computational methods to analyze single-cell or pathomics data to find sample-level trajectories or clusters associated with diseases. This remains challenging as single-cell/pathomics data are multi-scale, i.e., a sample is represented by clusters of cells/structures, and samples cannot be easily compared with each other. Here we propose PatIent Level analysis with Optimal Transport (PILOT). PILOT uses optimal transport to compute the Wasserstein distance between two individual single-cell samples. This allows us to perform unsupervised analysis at the sample level and uncover trajectories or cellular clusters associated with disease progression. We evaluate PILOT and competing approaches in single-cell genomics or pathomics studies involving various human diseases with up to 600 samples/patients and millions of cells or tissue structures. Our results demonstrate that PILOT detects disease-associated samples from large and complex single-cell or pathomics data. Moreover, PILOT provides a statistical approach to find changes in cell populations, gene expression, and tissue structures related to the trajectories or clusters supporting interpretation of predictions.

MRUNet-3D: A multi-stride residual 3D UNet for lung nodule segmentation.

Obtaining an accurate segmentation of the pulmonary nodules in computed tomography (CT) images is challenging. This is due to: (1) the heterogeneous nature of the lung nodules; (2) comparable visual characteristics between the nodules and their surroundings. A robust multi-scale feature extraction mechanism that can effectively obtain multi-scale representations at a granular level can improve segmentation accuracy. As the most commonly used network in lung nodule segmentation, UNet, its variants, and other image segmentation methods lack this robust feature extraction mechanism. In this study, we propose a multi-stride residual 3D UNet (MRUNet-3D) to improve the segmentation accuracy of lung nodules in CT images. It incorporates a multi-slide Res2Net block (MSR), which replaces the simple sequence of convolution layers in each encoder stage to effectively extract multi-scale features at a granular level from different receptive fields and resolutions while conserving the strengths of 3D UNet. The proposed method has been extensively evaluated on the publicly available LUNA16 dataset. Experimental results show that it achieves competitive segmentation performance with an average dice similarity coefficient of 83.47 % and an average surface distance of 0.35 mm on the dataset. More notably, our method has proven to be robust to the heterogeneity of lung nodules. It has also proven to perform better at segmenting small lung nodules. Ablation studies have shown that the proposed MSR and RFIA modules are fundamental to improving the performance of the proposed model.

MSCAN: multi-scale self- and cross-attention network for RNA methylation site prediction.

BACKGROUND: Epi-transcriptome regulation through post-transcriptional RNA modifications is essential for all RNA types. Precise recognition of RNA modifications is critical for understanding their functions and regulatory mechanisms. However, wet experimental methods are often costly and time-consuming, limiting their wide range of applications. Therefore, recent research has focused on developing computational methods, particularly deep learning (DL). Bidirectional long short-term memory (BiLSTM), convolutional neural network (CNN), and the transformer have demonstrated achievements in modification site prediction. However, BiLSTM cannot achieve parallel computation, leading to a long training time, CNN cannot learn the dependencies of the long distance of the sequence, and the Transformer lacks information interaction with sequences at different scales. This insight underscores the necessity for continued research and development in natural language processing (NLP) and DL to devise an enhanced prediction framework that can effectively address the challenges presented. RESULTS: This study presents a multi-scale self- and cross-attention network (MSCAN) to identify the RNA methylation site using an NLP and DL way. Experiment results on twelve RNA modification sites (m6A, m1A, m5C, m5U, m6Am, m7G, Ψ, I, Am, Cm, Gm, and Um) reveal that the area under the receiver operating characteristic of MSCAN obtains respectively 98.34%, 85.41%, 97.29%, 96.74%, 99.04%, 79.94%, 76.22%, 65.69%, 92.92%, 92.03%, 95.77%, 89.66%, which is better than the state-of-the-art prediction model. This indicates that the model has strong generalization capabilities. Furthermore, MSCAN reveals a strong association among different types of RNA modifications from an experimental perspective. A user-friendly web server for predicting twelve widely occurring human RNA modification sites (m6A, m1A, m5C, m5U, m6Am, m7G, Ψ, I, Am, Cm, Gm, and Um) is available at http://47.242.23.141/MSCAN/index.php . CONCLUSIONS: A predictor framework has been developed through binary classification to predict RNA methylation sites.

Light3DHS: A lightweight 3D hippocampus segmentation method using multiscale convolution attention and vision transformer.

The morphological analysis and volume measurement of the hippocampus are crucial to the study of many brain diseases. Therefore, an accurate hippocampal segmentation method is beneficial for the development of clinical research in brain diseases. U-Net and its variants have become prevalent in hippocampus segmentation of Magnetic Resonance Imaging (MRI) due to their effectiveness, and the architecture based on Transformer has also received some attention. However, some existing methods focus too much on the shape and volume of the hippocampus rather than its spatial information, and the extracted information is independent of each other, ignoring the correlation between local and global features. In addition, many methods cannot be effectively applied to practical medical image segmentation due to many parameters and high computational complexity. To this end, we combined the advantages of CNNs and ViTs (Vision Transformer) and proposed a simple and lightweight model: Light3DHS for the segmentation of the 3D hippocampus. In order to obtain richer local contextual features, the encoder first utilizes a multi-scale convolutional attention module (MCA) to learn the spatial information of the hippocampus. Considering the importance of local features and global semantics for 3D segmentation, we used a lightweight ViT to learn high-level features of scale invariance and further fuse local-to-global representation. To evaluate the effectiveness of encoder feature representation, we designed three decoders of different complexity to generate segmentation maps. Experiments on three common hippocampal datasets demonstrate that the network achieves more accurate hippocampus segmentation with fewer parameters. Light3DHS performs better than other state-of-the-art algorithms.

The heritability and structural correlates of resting-state fMRI complexity.

The complexity of fMRI signals quantifies temporal dynamics of spontaneous neural activity, which has been increasingly recognized as providing important insights into cognitive functions and psychiatric disorders. However, its heritability and structural underpinnings are not well understood. Here, we utilize multi-scale sample entropy to extract resting-state fMRI complexity in a large healthy adult sample from the Human Connectome Project. We show that fMRI complexity at multiple time scales is heritable in broad brain regions. Heritability estimates are modest and regionally variable. We relate fMRI complexity to brain structure including surface area, cortical myelination, cortical thickness, subcortical volumes, and total brain volume. We find that surface area is negatively correlated with fine-scale complexity and positively correlated with coarse-scale complexity in most cortical regions, especially the association cortex. Most of these correlations are related to common genetic and environmental effects. We also find positive correlations between cortical myelination and fMRI complexity at fine scales and negative correlations at coarse scales in the prefrontal cortex, lateral temporal lobe, precuneus, lateral parietal cortex, and cingulate cortex, with these correlations mainly attributed to common environmental effects. We detect few significant associations between fMRI complexity and cortical thickness. Despite the non-significant association with total brain volume, fMRI complexity exhibits significant correlations with subcortical volumes in the hippocampus, cerebellum, putamen, and pallidum at certain scales. Collectively, our work establishes the genetic basis and structural correlates of resting-state fMRI complexity across multiple scales, supporting its potential application as an endophenotype for psychiatric disorders.

Electrophysiological brain imaging based on simulation-driven deep learning in the context of epilepsy.

Identifying the location, the spatial extent and the electrical activity of distributed brain sources in the context of epilepsy through ElectroEncephaloGraphy (EEG) recordings is a challenging task because of the highly ill-posed nature of the underlying Electrophysiological Source Imaging (ESI) problem. To guarantee a unique solution, most existing ESI methods pay more attention to solve this inverse problem by imposing physiological constraints. This paper proposes an efficient ESI approach based on simulation-driven deep learning. Epileptic High-resolution 256-channels scalp EEG (Hr-EEG) signals are simulated in a realistic manner to train the proposed patient-specific model. More particularly, a computational neural mass model developed in our team is used to generate the temporal dynamics of the activity of each dipole while the forward problem is solved using a patient-specific three-shell realistic head model and the boundary element method. A Temporal Convolutional Network (TCN) is considered in the proposed model to capture local spatial patterns. To enable the model to observe the EEG signals from different scale levels, the multi-scale strategy is leveraged to capture the overall features and fine-grain features by adjusting the convolutional kernel size. Then, the Long Short-Term Memory (LSTM) is used to extract temporal dependencies among the computed spatial features. The performance of the proposed method is evaluated through three different scenarios of realistic synthetic interictal Hr-EEG data as well as on real interictal Hr-EEG data acquired in three patients with drug-resistant partial epilepsy, during their presurgical evaluation. A performance comparison study is also conducted with two other deep learning-based methods and four classical ESI techniques. The proposed model achieved a Dipole Localization Error (DLE) of 1.39 and Normalized Hamming Distance (NHD) of 0.28 in the case of one patch with SNR of 10 dB. In the case of two uncorrelated patches with an SNR of 10 dB, obtained DLE and NHD were respectively 1.50 and 0.28. Even in the more challenging scenario of two correlated patches with an SNR of 10 dB, the proposed approach still achieved a DLE of 3.74 and an NHD of 0.43. The results obtained on simulated data demonstrate that the proposed method outperforms the existing methods for different signal-to-noise and source configurations. The good behavior of the proposed method is also confirmed on real interictal EEG data. The robustness with respect to noise makes it a promising and alternative tool to localize epileptic brain areas and to reconstruct their electrical activities from EEG signals.

Multi-scale brain attributes contribute to the distribution of diffuse glioma subtypes.

Gliomas are primary brain tumors and are among the most malignant types. Adult-type diffuse gliomas can be classified based on their histological and molecular signatures as IDH-wildtype glioblastoma, IDH-mutant astrocytoma, and IDH-mutant and 1p/19q-codeleted oligodendroglioma. Recent studies have shown that each subtype of glioma has its own specific distribution pattern. However, the mechanisms underlying the specific distributions of glioma subtypes are not entirely clear despite partial explanations such as cell origin. To investigate the impact of multi-scale brain attributes on glioma distribution, we constructed cumulative frequency maps for diffuse glioma subtypes based on T1w structural images and evaluated the spatial correlation between tumor frequency and diverse brain attributes, including postmortem gene expression, functional connectivity metrics, cerebral perfusion, glucose metabolism, and neurotransmitter signaling. Regression models were constructed to evaluate the contribution of these factors to the anatomic distribution of different glioma subtypes. Our findings revealed that the three different subtypes of gliomas had distinct distribution patterns, showing spatial preferences toward different brain environmental attributes. Glioblastomas were especially likely to occur in regions enriched with synapse-related pathways and diverse neurotransmitter receptors. Astrocytomas and oligodendrogliomas preferentially occurred in areas enriched with genes associated with neutrophil-mediated immune responses. The functional network characteristics and neurotransmitter distribution also contributed to oligodendroglioma distribution. Our results suggest that different brain transcriptomic, neurotransmitter, and connectomic attributes are the factors that determine the specific distributions of glioma subtypes. These findings highlight the importance of bridging diverse scales of biological organization when studying neurological dysfunction.

Scale-Dependent Rheology of Synovial Fluid Lubricating Macromolecules.

Synovial fluid (SF) is the complex biofluid that facilitates the exceptional lubrication of articular cartilage in joints. Its primary lubricating macromolecules, the linear polysaccharide hyaluronic acid (HA) and the mucin-like glycoprotein proteoglycan 4 (PRG4 or lubricin), interact synergistically to reduce boundary friction. However, the precise manner in which these molecules influence the rheological properties of SF remains unclear. This study aimed to elucidate this by employing confocal microscopy and multiscale rheometry to examine the microstructure and rheology of solutions containing recombinant human PRG4 (rhPRG4) and HA. Contrary to previous assumptions of an extensive HA-rhPRG4 network, it is discovered that rhPRG4 primarily forms stiff, gel-like aggregates. The properties of these aggregates, including their size and stiffness, are found to be influenced by the viscoelastic characteristics of the surrounding HA matrix. Consequently, the rheology of this system is not governed by a single length scale, but instead responds as a disordered, hierarchical network with solid-like rhPRG4 aggregates distributed throughout the continuous HA phase. These findings provide new insights into the biomechanical function of PRG4 in cartilage lubrication and may have implications in the development of HA-based therapies for joint diseases like osteoarthritis.

TMNP: a transcriptome-based multi-scale network pharmacology platform for herbal medicine.

One of the most difficult problems that hinder the development and application of herbal medicine is how to illuminate the global effects of herbs on the human body. Currently, the chemo-centric network pharmacology methodology regards herbs as a mixture of chemical ingredients and constructs the 'herb-compound-target-disease' connections based on bioinformatics methods, to explore the pharmacological effects of herbal medicine. However, this approach is severely affected by the complexity of the herbal composition. Alternatively, gene-expression profiles induced by herbal treatment reflect the overall biological effects of herbs and are suitable for studying the global effects of herbal medicine. Here, we develop an online transcriptome-based multi-scale network pharmacology platform (TMNP) for exploring the global effects of herbal medicine. Firstly, we build specific functional gene signatures for different biological scales from molecular to higher tissue levels. Then, specific algorithms are designed to measure the correlations of transcriptional profiles and types of gene signatures. Finally, TMNP uses pharmacotranscriptomics of herbal medicine as input and builds associations between herbs and different biological scales to explore the multi-scale effects of herb medicine. We applied TMNP to a single herb Astragalus membranaceus and Xuesaitong injection to demonstrate the power to reveal the multi-scale effects of herbal medicine. TMNP integrating herbal medicine and multiple biological scales into the same framework, will greatly extend the conventional network pharmacology model centering on the chemical components, and provide a window for systematically observing the complex interactions between herbal medicine and the human body. TMNP is available at http://www.bcxnfz.top/TMNP.

A multi-step and multi-scale bioinformatic protocol to investigate potential SARS-CoV-2 vaccine targets.

The COVID-19 pandemic has highlighted the need to come out with quick interventional solutions that can now be obtained through the application of different bioinformatics software to actively improve the success rate. Technological advances in fields such as computer modeling and simulation are enriching the discovery, development, assessment and monitoring for better prevention, diagnosis, treatment and scientific evidence generation of specific therapeutic strategies. The combined use of both molecular prediction tools and computer simulation in the development or regulatory evaluation of a medical intervention, are making the difference to better predict the efficacy and safety of new vaccines. An integrated bioinformatics pipeline that merges the prediction power of different software that act at different scales for evaluating the elicited response of human immune system against every pathogen is proposed. As a working example, we applied this problem solving protocol to predict the cross-reactivity of pre-existing vaccination interventions against SARS-CoV-2.