Quantitative Biology

Computational prediction and functional analysis of arsenic-binding proteins in human cells Abstract Background Arsenic has a broad anti-cancer ability against hematologic malignancies and solid tumors. To systematically understand the biological functions of arsenic, we need to identify arsenic-binding proteins in human cells. However, due to lack of effective theoretical tools and experimental methods, only a few arsenic-binding proteins have been identified. Methods Based on the crystal structure of ArsM, we generated a single mutation free energy profile for arsenic binding using free energy perturbation methods. Multiple validations provide an indication that our computational model has the ability to predict arsenic-binding proteins with desirable accuracy. We subsequently apply this computational model to scan the entire human genome to identify all the potential arsenic-binding proteins. Results: The computationally predicted arsenic-binding proteins show a wide range of biological functions, especially in the signaling transduction pathways. In the signaling transduction pathways, arsenic directly binds to the key factors (e.g., Notch receptors, Notch ligands, Wnt family proteins, TGF-beta, and their interacting proteins) and results in significant inhibitions on their enzymatic activities, further having a crucial impact on the related signaling pathways. Conclusions Arsenic has a significant impact on signaling transduction in cells. Arsenic binding to proteins can lead to dysfunctions of the target proteins, having crucial impacts on both signaling pathway and gene transcription. We hope that the computationally predicted arsenic-binding proteins and the functional analysis can provide a novel insight into the biological functions of arsenic, revealing a mechanism for the broad anti-cancer of arsenic.

Current challenges and solutions of de novo assembly Abstract Background Next-generation sequencing (NGS) technologies have fostered an unprecedented proliferation of high-throughput sequencing projects and a concomitant development of novel algorithms for the assembly of short reads. However, numerous technical or computational challenges in de novo assembly still remain, although many new ideas and solutions have been suggested to tackle the challenges in both experimental and computational settings. Results In this review, we first briefly introduce some of the major challenges faced by NGS sequence assembly. Then, we analyze the characteristics of various sequencing platforms and their impact on assembly results. After that, we classify de novo assemblers according to their frameworks (overlap graph-based, de Bruijn graph-based and string graph-based), and introduce the characteristics of each assembly tool and their adaptation scene. Next, we introduce in detail the solutions to the main challenges of de novo assembly of next generation sequencing data, single-cell sequencing data and single molecule sequencing data. At last, we discuss the application of SMS long reads in solving problems encountered in NGS assembly. Conclusions This review not only gives an overview of the latest methods and developments in assembly algorithms, but also provides guidelines to determine the optimal assembly algorithm for a given input sequencing data type.

Predicting enhancer-promoter interaction from genomic sequence with deep neural networks Abstract Background In the human genome, distal enhancers are involved in regulating target genes through proximal promoters by forming enhancer-promoter interactions. Although recently developed high-throughput experimental approaches have allowed us to recognize potential enhancer-promoter interactions genome-wide, it is still largely unclear to what extent the sequence-level information encoded in our genome help guide such interactions. Methods Here we report a new computational method (named "SPEID") using deep learning models to predict enhancer-promoter interactions based on sequence-based features only, when the locations of putative enhancers and promoters in a particular cell type are given. Results Our results across six different cell types demonstrate that SPEID is effective in predicting enhancer-promoter interactions as compared to state-of-the-art methods that only use information from a single cell type. As a proof-of-principle, we also applied SPEID to identify somatic non-coding mutations in melanoma samples that may have reduced enhancer-promoter interactions in tumor genomes. Conclusions This work demonstrates that deep learning models can help reveal that sequence-based features alone are sufficient to reliably predict enhancer-promoter interactions genome-wide.

QB : Embracing the future of quantitative understanding and engineering of life

High-throughput quantification of microbial birth and death dynamics using fluorescence microscopy Abstract Background Microbes live in dynamic environments where nutrient concentrations fluctuate. Quantifying fitness in terms of birth rate and death rate in a wide range of environments is critical for understanding microbial evolution and ecology. Methods Here, using high-throughput time-lapse microscopy, we have quantified how Saccharomyces cerevisiae mutants incapable of synthesizing an essential metabolite (auxotrophs) grow or die in various concentrations of the required metabolite.We establish that cells normally expressing fluorescent proteins lose fluorescence upon death and that the total fluorescence in an imaging frame is proportional to the number of live cells even when cells form multiple layers. We validate our microscopy approach of measuring birth and death rates using flow cytometry, cell counting, and chemostat culturing. Results For lysine-requiring cells, very low concentrations of lysine are not detectably consumed and do not support cell birth, but delay the onset of death phase and reduce the death rate compared to no lysine. In contrast, in low hypoxanthine, hypoxanthine-requiring cells can produce new cells, yet also die faster than in the absence of hypoxanthine. For both strains, birth rates under various metabolite concentrations are better described by the sigmoidal-shaped Moser model than the well-known Monod model, while death rates can vary with metabolite concentration and time. Conclusions Our work reveals how time-lapse microscopy can be used to discover non-intuitive microbial birth and death dynamics and to quantify growth rates in many environments.

A survey of web resources and tools for the study of TCM network pharmacology Abstract Background Traditional Chinese medicine (TCM) treats diseases in a holistic manner, while TCM formulae are multi-component, multi-target agents at the molecular level. Thus there are many parallels between the key ideas of TCM pharmacology and network pharmacology. These years, TCM network pharmacology has developed as an interdisciplinary of TCM science and network pharmacology, which studies the mechanism of TCM at the molecular level and in the context of biological networks. It provides a new research paradigm that can use modern biomedical science to interpret the mechanism of TCM, which is promising to accelerate the modernization and internationalization of TCM. Results In this paper we introduce state-of-the-art free data sources, web servers and softwares that can be used in the TCM network pharmacology, including databases of TCM, drug targets and diseases, web servers for the prediction of drug targets, and tools for network and functional analysis. Conclusions This review could help experimental pharmacologists make better use of the existing data and methods in their study of TCM.

Pharmacodynamics simulation of HOEC by a computational model of arachidonic acid metabolic network Abstract Background Arachidonic acid (AA) metabolic network is activated in the most inflammatory related diseases, and small-molecular drugs targeting AA network are increasingly available. However, side effects of above mentioned drugs have always been the biggest obstacle. (+)-2-(1-hydroxyl-4-oxocyclohexyl) ethyl caffeate (HOEC), a natural product acted as an inhibitor of 5-lipoxygenase (5-LOX) and 15-LOX in vitro, exhibited weaker therapeutic effect in high dose than that in low dose to collagen induced arthritis (CIA) rats. In this study, we tried to elucidate the potential regulatory mechanism by using quantitative pharmacology. Methods First, we generated an experimental data set by monitoring the dynamics of AA metabolites' concentration in A23187 stimulated and different doses of HOEC co-incubated RAW264.7. Then we constructed a dynamic model of A23187-stimulated AA metabolic model to evaluate how a model-based simulation of AA metabolic data assists to find the most suitable treatment dose by predicting the pharmacodynamics of HOEC. Results Compared to the experimental data, the model could simulate the inhibitory effect of HOEC on 5-LOX and 15-LOX, and reproduced the increase of the metabolic flux in the cyclooxygenase (COX) pathway. However, a concomitant, early-stage of stimulation-related decrease of prostaglandins (PGs) production in HOEC incubated RAW264.7 cells was not simulated in the model. Conclusion Using the model, we predict that higher dose of HOEC disrupts the flux balance in COX and LOX of the AA network, and increased COX flux can interfere the curative effects of LOX inhibitor on resolution of inflammation which is crucial for the efficient and safe drug design.

Time-scale separation and stochasticity conspire to impact phenotypic dynamics in the canonical and inverted Bacillus subtilis core genetic regulation circuits Abstract Background In this work, we study two seemingly unrelated aspects of core genetic nonlinear dynamical control of the competence phenotype in Bacillus subtilis, a common Gram-positive bacterium living in the soil. Methods We focus on hitherto unchartered aspects of the dynamics by exploring the effect of time-scale separation between transcription and translation and, as well, the effect of intrinsic molecular stochasticity. We consider these aspects of regulatory control as two possible evolutionary handles. Results Hence, using theory and computations, we study how the onset of oscillations breaks the excitability-based competence phenotype in two topologically close evolutionary-competing circuits: the canonical "wild-type" regulation circuit selected by Evolution and the corresponding indirect-feedback inverted circuit that failed to be selected by Evolution, as was shown elsewhere, due to dynamical reasons. Conclusions Relying on in-silico perturbation of the living state, we show that the canonical core genetic regulation of excitability-based competence is more robust against switching to phenotype-breaking oscillations than the inverted feedback organism. We show how this is due to time-scale separation and stochasticity.

Computational methods and applications for quantitative systems pharmacology Abstract Background Quantitative systems pharmacology (QSP) is an emerging discipline that integrates diverse data to quantitatively explore the interactions between drugs and multi-scale systems including small compounds, nucleic acids, proteins, pathways, cells, organs and disease processes. Results Various computational methods such as ADME/Tevaluation, molecular modeling, logical modeling, network modeling, pathway analysis, multi-scale systems pharmacology platforms and virtual patient for QSP have been developed. We reviewed the major progresses and broad applications in medical guidance, drug discovery and exploration of pharmacodynamic material basis and mechanism of traditional Chinese medicine. Conclusion QSP has significant achievements in recent years and is a promising approach for quantitative evaluation of drug efficacy and systematic exploration of mechanisms of action of drugs.

Insights into the antineoplastic mechanism of Chelidonium majus via systems pharmacology approach Abstract Background The antineoplastic activity of Chelidonium majus has been reported, but its mechanism of action (MoA) is unsuspected. The emerging theory of systems pharmacology may be a useful approach to analyze the complicated MoA of this multi-ingredient traditional Chinese medicine (TCM). Methods We collected the ingredients and related compound-target interactions of C. majus from several databases. The bSDTNBI (balanced substructure-drug-target network-based inference) method was applied to predict each ingredient's targets. Pathway enrichment analysis was subsequently conducted to illustrate the potential MoA, and prognostic genes were identified to predict the certain types of cancers that C. majus might be beneficial in treatment. Bioassays and literature survey were used to validate the in silico results. Results Systems pharmacology analysis demonstrated that C. majus exerted experimental or putative interactions with 18 cancer-associated pathways, and might specifically act on 13 types of cancers. Chelidonine, sanguinarine, chelerythrine, berberine, and coptisine, which are the predominant components of C. majus, may suppress the cancer genes by regulating cell cycle, inducing cell apoptosis and inhibiting proliferation. Conclusions The antineoplastic MoA of C. majus was investigated by systems pharmacology approach. C. majus exhibited promising pharmacological effect against cancer, and may consequently be useful material in further drug development. The alkaloids are the key components in C. majus that exhibit anticancer activity.

Alexandros Sfakianakis
Anapafseos 5 . Agios Nikolaos
Crete.Greece.72100
2841026182

6948891480

alsfakia@gmail.com