SPLINDID: a semi-parametric, model-based method for obtaining transcription rates and gene regulation parameters from genomic and proteomic expression profiles

Gene-Centric Model Approaches for Accurate Prediction of Pesticide Biodegradation in Soils

Pesticides are widely used in agriculture despite their negative impact on ecosystems and human health. Biogeochemical modeling facilitates the mechanistic understanding of microbial controls on pesticide turnover in soils. We propose to inform models of coupled microbial dynamics and pesticide turnover with measurements of the abundance and expression of functional genes. To assess the advantages of informing models with genetic data, we developed a novel "gene-centric" model and compared model variants of differing structural complexity against a standard biomass-based model. The models were calibrated and validated using data from two batch experiments in which the degradation of the pesticides dichlorophenoxyacetic acid (2,4-D) and 2-methyl-4-chlorophenoxyacetic acid (MCPA) were observed in soil. When calibrating against data on pesticide mineralization, the gene-centric and biomass-based models performed equally well. However, accounting for pesticide-triggered gene regulation allows improved performance in capturing microbial dynamics and in predicting pesticide mineralization. This novel modeling approach also reveals a hysteretic relationship between pesticide degradation rates and gene expression, implying that the biodegradation performance in soils cannot be directly assessed by measuring the expression of functional genes. Our gene-centric model provides an effective approach for exploiting molecular biology data to simulate pesticide degradation in soils.

Integrated Dissection of Cysteine Oxidative Post-translational Modification Proteome During Cardiac Hypertrophy

Cysteine oxidative modification of cellular proteins is crucial for many aspects of cardiac hypertrophy development. However, integrated dissection of multiple types of cysteine oxidative post-translational modifications (O-PTM) of proteomes in cardiac hypertrophy is currently missing. Here we developed a novel discovery platform that encompasses a customized biotin switch-based quantitative proteomics pipeline and an advanced analytic workflow to comprehensively profile the landscape of cysteine O-PTM in an ISO-induced cardiac hypertrophy mouse model. Specifically, we identified a total of 1655 proteins containing 3324 oxidized cysteine sites by at least one of the following three modifications: reversible cysteine O-PTM, cysteine sulfinylation (CysSO₂H), and cysteine sulfonylation (CysSO₃H). Analyzing the hypertrophy signatures that are reproducibly discovered from this computational workflow unveiled four biological processes with increased cysteine O-PTM. Among them, protein phosphorylation, creatine metabolism, and response to elevated Ca²⁺ pathways exhibited an elevation of cysteine O-PTM in early stages, whereas glucose metabolism enzymes were increasingly modified in later stages, illustrating a temporal regulatory map in cardiac hypertrophy. Our cysteine O-PTM platform depicts a dynamic and integrated landscape of the cysteine oxidative proteome, through the extracted molecular signatures, and provides critical mechanistic insights in cardiac hypertrophy. Data are available via ProteomeXchange with identifier PXD010336.

Reconstructing Dynamic Promoter Activity Profiles from Reporter Gene Data

Accurate characterization of promoter activity is important when designing expression systems for systems biology and metabolic engineering applications. Promoters that respond to changes in the environment enable the dynamic control of gene expression without the necessity of inducer compounds, for example. However, the dynamic nature of these processes poses challenges for estimating promoter activity. Most experimental approaches utilize reporter gene expression to estimate promoter activity. Typically the reporter gene encodes a fluorescent protein that is used to infer a constant promoter activity despite the fact that the observed output may be dynamic and is a number of steps away from the transcription process. In fact, some promoters that are often thought of as constitutive can show changes in activity when growth conditions change. For these reasons, we have developed a system of ordinary differential equations for estimating dynamic promoter activity for promoters that change their activity in response to the environment that is robust to noise and changes in growth rate. Our approach, inference of dynamic promoter activity (PromAct), improves on existing methods by more accurately inferring known promoter activity profiles. This method is also capable of estimating the correct scale of promoter activity and can be applied to quantitative data sets to estimate quantitative rates.

Analysis of gene coexpression by B-spline based CoD estimation

The gene coexpression study has emerged as a novel holistic approach for microarray data analysis. Different indices have been used in exploring coexpression relationship, but each is associated with certain pitfalls. The Pearson's correlation coefficient, for example, is not capable of uncovering nonlinear pattern and directionality of coexpression. Mutual information can detect nonlinearity but fails to show directionality. The coefficient of determination (CoD) is unique in exploring different patterns of gene coexpression, but so far only applied to discrete data and the conversion of continuous microarray data to the discrete format could lead to information loss. Here, we proposed an effective algorithm, CoexPro, for gene coexpression analysis. The new algorithm is based on B-spline approximation of coexpression between a pair of genes, followed by CoD estimation. The algorithm was justified by simulation studies and by functional semantic similarity analysis. The proposed algorithm is capable of uncovering both linear and a specific class of nonlinear relationships from continuous microarray data. It can also provide suggestions for possible directionality of coexpression to the researchers. The new algorithm presents a novel model for gene coexpression and will be a valuable tool for a variety of gene expression and network studies. The application of the algorithm was demonstrated by an analysis on ligand-receptor coexpression in cancerous and noncancerous cells. The software implementing the algorithm is available upon request to the authors.

Modeling transitions in body composition: the approach to steady state for anthropometric measures and physiological functions in the Minnesota human starvation study

Background

This study evaluated whether the changes in several anthropometric and functional measures during caloric restriction combined with walking and treadmill exercise would fit a simple model of approach to steady state (a plateau) that can be solved using spreadsheet software (Microsoft Excel^®). We hypothesized that transitions in waist girth and several body compartments would fit a simple exponential model that approaches a stable steady-state.

Methods

The model (an equation) was applied to outcomes reported in the Minnesota starvation experiment using Microsoft Excel's Solver^® function to derive rate parameters (k) and projected steady state values. However, data for most end-points were available only at t = 0, 12 and 24 weeks of caloric restriction. Therefore, we derived 2 new equations that enable model solutions to be calculated from 3 equally spaced data points.

Results

For the group of male subjects in the Minnesota study, body mass declined with a first order rate constant of about 0.079 wk^-1. The fractional rate of loss of fat free mass, which includes components that remained almost constant during starvation, was 0.064 wk^-1, compared to a rate of loss of fat mass of 0.103 wk^-1. The rate of loss of abdominal fat, as exemplified by the change in the waist girth, was 0.213 wk^-1.

On average, 0.77 kg was lost per cm of waist girth. Other girths showed rates of loss between 0.085 and 0.131 wk^-1. Resting energy expenditure (REE) declined at 0.131 wk^-1. Changes in heart volume, hand strength, work capacity and N excretion showed rates of loss in the same range. The group of 32 subjects was close to steady state or had already reached steady state for the variables under consideration at the end of semi-starvation.

Conclusion

When energy intake is changed to new, relatively constant levels, while physical activity is maintained, changes in several anthropometric and physiological measures can be modeled as an exponential approach to steady state using software that is widely available. The 3 point method for parameter estimation provides a criterion for testing whether change in a variable can be usefully modelled with exponential kinetics within the time range for which data are available.

Simulated maximum likelihood method for estimating kinetic rates in gene expression

MOTIVATION

Kinetic rate in gene expression is a key measurement of the stability of gene products and gives important information for the reconstruction of genetic regulatory networks. Recent developments in experimental technologies have made it possible to measure the numbers of transcripts and protein molecules in single cells. Although estimation methods based on deterministic models have been proposed aimed at evaluating kinetic rates from experimental observations, these methods cannot tackle noise in gene expression that may arise from discrete processes of gene expression, small numbers of mRNA transcript, fluctuations in the activity of transcriptional factors and variability in the experimental environment.

RESULTS

In this paper, we develop effective methods for estimating kinetic rates in genetic regulatory networks. The simulated maximum likelihood method is used to evaluate parameters in stochastic models described by either stochastic differential equations or discrete biochemical reactions. Different types of non-parametric density functions are used to measure the transitional probability of experimental observations. For stochastic models described by biochemical reactions, we propose to use the simulated frequency distribution to evaluate the transitional density based on the discrete nature of stochastic simulations. The genetic optimization algorithm is used as an efficient tool to search for optimal reaction rates. Numerical results indicate that the proposed methods can give robust estimations of kinetic rates with good accuracy.

Application of SPLINDID, a semiparametric, model-based method for pharmacogenomic modeling of mRNA dynamics

PURPOSE

This study was conducted to evaluate the applicability of SPLINDID, a semiparametric, model-based approach for obtaining transcription rates from the pharmacodynamics of mRNA expression.

METHODS

A nonparametric exponential cubic spline function was used to obtain the transcription rate profile and the dynamics of mRNA expression was fitted using compartmental approaches. The transcription rate profile and mRNA degradation parameter was estimated using maximum likelihood method of ADAPT II software.

RESULTS

Data sets containing noise for mRNA levels were simulated for four diverse pharmaceutically relevant conditions: receptor nonlinearity, a model in which the variant mRNAs differing in mRNA degradation constants were transcribed and for a minimal model of the cell cycle. SPLINDID was able to fit the data sets and accurately recapitulate the transcription rate profiles normalized to the mRNA degradation rate constants. The model was also challenged using experimental data containing time profiles of cell-cycle-regulated genes.

CONCLUSIONS

The SPLINDID approach is flexible in capturing complicated/complex mRNA profiles that are encountered in many experimental data sets.