Modeling the temporal dynamics of master regulators and CtrA proteolysis in Caulobacter crescentus cell cycle.

The cell cycle of Caulobacter crescentus involves the polar morphogenesis and an asymmetric cell division driven by precise interactions and regulations of proteins, which makes Caulobacter an ideal model organism for investigating bacterial cell development and differentiation. The abundance of molecular data accumulated on Caulobacter motivates system biologists to analyze the complex regulatory network of cell cycle via quantitative modeling. In this paper, We propose a comprehensive model to accurately characterize the underlying mechanisms of cell cycle regulation based on the study of: a) chromosome replication and methylation; b) interactive pathways of five master regulatory proteins including DnaA, GcrA, CcrM, CtrA, and SciP, as well as novel consideration of their corresponding mRNAs; c) cell cycle-dependent proteolysis of CtrA through hierarchical protease complexes. The temporal dynamics of our simulation results are able to closely replicate an extensive set of experimental observations and capture the main phenotype of seven mutant strains of Caulobacter crescentus. Collectively, the proposed model can be used to predict phenotypes of other mutant cases, especially for nonviable strains which are hard to cultivate and observe. Moreover, the module of cyclic proteolysis is an efficient tool to study the metabolism of proteins with similar mechanisms.

Cell cycle control and environmental response by second messengers in Caulobacter crescentus.

BACKGROUND: Second messengers, c-di-GMP and (p)ppGpp, are vital regulatory molecules in bacteria, influencing cellular processes such as biofilm formation, transcription, virulence, quorum sensing, and proliferation. While c-di-GMP and (p)ppGpp are both synthesized from GTP molecules, they play antagonistic roles in regulating the cell cycle. In C. crescentus, c-di-GMP works as a major regulator of pole morphogenesis and cell development. It inhibits cell motility and promotes S-phase entry by inhibiting the activity of the master regulator, CtrA. Intracellular (p)ppGpp accumulates under starvation, which helps bacteria to survive under stressful conditions through regulating nucleotide levels and halting proliferation. (p)ppGpp responds to nitrogen levels through RelA-SpoT homolog enzymes, detecting glutamine concentration using a nitrogen phosphotransferase system (PTS Ntr). This work relates the guanine nucleotide-based second messenger regulatory network with the bacterial PTS Ntr system and investigates how bacteria respond to nutrient availability. RESULTS: We propose a mathematical model for the dynamics of c-di-GMP and (p)ppGpp in C. crescentus and analyze how the guanine nucleotide-based second messenger system responds to certain environmental changes communicated through the PTS Ntr system. Our mathematical model consists of seven ODEs describing the dynamics of nucleotides and PTS Ntr enzymes. Our simulations are consistent with experimental observations and suggest, among other predictions, that SpoT can effectively decrease c-di-GMP levels in response to nitrogen starvation just as well as it increases (p)ppGpp levels. Thus, the activity of SpoT (or its homologues in other bacterial species) can likely influence the cell cycle by influencing both c-di-GMP and (p)ppGpp. CONCLUSIONS: In this work, we integrate current knowledge and experimental observations from the literature to formulate a novel mathematical model. We analyze the model and demonstrate how the PTS Ntr system influences (p)ppGpp, c-di-GMP, GMP and GTP concentrations. While this model does not consider all aspects of PTS Ntr signaling, such as cross-talk with the carbon PTS system, here we present our first effort to develop a model of nutrient signaling in C. crescentus.

Feature selection of gene expression data for Cancer classification using double RBF-kernels.

BACKGROUND: Using knowledge-based interpretation to analyze omics data can not only obtain essential information regarding various biological processes, but also reflect the current physiological status of cells and tissue. The major challenge to analyze gene expression data, with a large number of genes and small samples, is to extract disease-related information from a massive amount of redundant data and noise. Gene selection, eliminating redundant and irrelevant genes, has been a key step to address this problem. RESULTS: The modified method was tested on four benchmark datasets with either two-class phenotypes or multiclass phenotypes, outperforming previous methods, with relatively higher accuracy, true positive rate, false positive rate and reduced runtime. CONCLUSIONS: This paper proposes an effective feature selection method, combining double RBF-kernels with weighted analysis, to extract feature genes from gene expression data, by exploring its nonlinear mapping ability.

Computational prediction of therapeutic peptides based on graph index.

As therapeutic peptides have been taken into consideration in disease therapy in recent years, many biologists spent time and labor to verify various functional peptides from a large number of peptide sequences. In order to reduce the workload and increase the efficiency of identification of functional proteins, we propose a sequence-based model, q-FP (functional peptide prediction based on the q-Wiener Index), capable of recognizing potentially functional proteins. We extract three types of features by mixing graphic representation and statistical indices based on the q-Wiener index and physicochemical properties of amino acids. Our support-vector-machine-based model achieves an accuracy of 96.71%, 93.34%, 98.40%, and 91.40% for anticancer, virulent, and allergenic proteins datasets, respectively, by using 5-fold cross validation.

Protein sequence analysis by incorporating modified chaos game and physicochemical properties into Chou's general pseudo amino acid composition.

In this contribution we introduced a novel graphical method to compare protein sequences. By mapping a protein sequence into 3D space based on codons and physicochemical properties of 20 amino acids, we are able to get a unique P-vector from the 3D curve. This approach is consistent with wobble theory of amino acids. We compute the distance between sequences by their P-vectors to measure similarities/dissimilarities among protein sequences. Finally, we use our method to analyze four datasets and get better results compared with previous approaches.

The complete chloroplast genome of Arachis lutescens Krapov. & Rigoni (Fabaceae).

Arachis lutescens Krapov. & Rigoni 1958 is an important species due to their potentially extensive applications for cultivated peanut breeding. The whole chloroplast genome of A. lutescens was successfully assembled and annotated for the first time. The complete chloroplast genome of A. lutescens is a typically circular structure of 156,398 bp with a GC content of 36.3%. It comprises a large single-copy (LSC) region of 85,950 bp, a small single-copy (SSC) region of 18,800 bp, and two inverted repeat regions (IRs) of 25,824 bp, each. The plastome of A. lutescens contains a total of 125 genes, including 81 protein-coding genes, 36 tRNAs, and eight rRNAs. The phylogenetic analysis strongly supports the close relationship between A. lutescens and cultivated peanut clades. This study contributes to our understanding of the molecular characteristics and evolutionary relationships of this plant species.

Turing-pattern model of scaffolding proteins that establish spatial asymmetry during the cell cycle of Caulobacter crescentus.

The crescent-shaped bacterium Caulobacter crescentus divides asymmetrically into a sessile (stalked) cell and a motile (flagellated) cell. This dimorphic cell division cycle is driven by the asymmetric appearance of scaffolding proteins at the cell's stalk and flagellum poles. The scaffolding proteins recruit enzyme complexes that phosphorylate and degrade a master transcription factor, CtrA, and the abundance and phosphorylation state of CtrA control the onset of DNA synthesis and the differentiation of stalked and flagellated cell types. In this study, we use a Turing-pattern mechanism to simulate the spatiotemporal dynamics of scaffolding proteins in Caulobacter and how they influence the abundance and intracellular distribution of CtrA∼P. Our mathematical model captures crucial features of wild-type and mutant strains and predicts the distributions of CtrA∼P and signaling proteins in mutant strains. Our model accounts for Caulobacter polar morphogenesis and shows how spatial localization and phosphosignaling cooperate to establish asymmetry during the cell cycle.

Uncovering the dynamic effects of DEX treatment on lung cancer by integrating bioinformatic inference and multiscale modeling of scRNA-seq and proteomics data.

Lung cancer is one of the leading causes of cancer-related death, with a five-year survival rate of 18%. It is a priority for us to understand the underlying mechanisms affecting lung cancer therapeutics' implementation and effectiveness. In this study, we combine the power of Bioinformatics and Systems Biology to comprehensively uncover functional and signaling pathways of drug treatment using bioinformatics inference and multiscale modeling of both scRNA-seq data and proteomics data. Based on a time series of lung adenocarcinoma derived A549 cells after DEX treatment, we first identified the differentially expressed genes (DEGs) in those lung cancer cells. Through the interrogation of regulatory network of those DEGs, we identified key hub genes including TGFß, MYC, and SMAD3 varied underlie DEX treatment. Further gene set enrichment analysis revealed the TGFß signaling pathway as the top enriched term. Those genes involved in the TGFß pathway and their crosstalk with the ERBB pathway presented a strong survival prognosis in clinical lung cancer samples. With the basis of biological validation and literature-based curation, a multiscale model of tumor regulation centered on both TGFß-induced and ERBB-amplified signaling pathways was developed to characterize the dynamic effects of DEX therapy on lung cancer cells. Our simulation results were well matched to available data of SMAD2, FOXO3, TGFß1, and TGFßR1 over the time course. Moreover, we provided predictions of different doses to illustrate the trend and therapeutic potential of DEX treatment. The innovative and cross-disciplinary approach can be further applied to other computational studies in tumorigenesis and oncotherapy. We released the approach as a user-friendly tool named BIMM (Bioinformatic Inference and Multiscale Modeling), with all the key features available at https://github.com/chenm19/BIMM.