Dynamic regulation of the tryptophan operon: a modeling study and comparison with experimental data.
Ontology highlight
ABSTRACT:
SUBMITTER: Santillan M
PROVIDER: S-EPMC29262 | biostudies-literature | 2001 Feb
REPOSITORIES: biostudies-literature
Publications
Dynamic regulation of the tryptophan operon: a modeling study and comparison with experimental data.
Proceedings of the National Academy of Sciences of the United States of America 20010201 4
A mathematical model for regulation of the tryptophan operon is presented. This model takes into account repression, feedback enzyme inhibition, and transcriptional attenuation. Special attention is given to model parameter estimation based on experimental data. The model's system of delay differential equations is numerically solved, and the results are compared with experimental data on the temporal evolution of enzyme activity in cultures of Escherichia coli after a nutritional shift (minimal ...[more]
Similar Datasets
Project description:BACKGROUND:Complete genomic sequences of closely related organisms, such as the chlamydiae, afford the opportunity to assess significant strain differences against a background of many shared characteristics. The chlamydiae are ubiquitous intracellular parasites that are important pathogens of humans and other organisms. Tryptophan limitation caused by production of interferon-gamma by the host and subsequent induction of indoleamine dioxygenase is a key aspect of the host-parasite interaction. It appears that the chlamydiae have learned to recognize tryptophan depletion as a signal for developmental remodeling. The consequent non-cultivable state of persistence can be increasingly equated to chronic disease conditions. RESULTS:The genes encoding enzymes of tryptophan biosynthesis were the focal point of this study. Chlamydophila psittaci was found to possess a compact operon containing PRPP synthase, kynureninase, and genes encoding all but the first step of tryptophan biosynthesis. All but one of the genes exhibited translational coupling. Other chlamydiae (Chlamydia trachomatis, C. muridarum and Chlamydophila pneumoniae) lack genes encoding PRPP synthase, kynureninase, and either lack tryptophan-pathway genes altogether or exhibit various stages of reductive loss. The origin of the genes comprising the trp operon does not seem to have been from lateral gene transfer. CONCLUSIONS:The factors that accommodate the transition of different chlamydial species to the persistent (chronic) state of pathogenesis include marked differences in strategies deployed to obtain tryptophan from host resources. C. psittaci appears to have a novel mechanism for intercepting an early intermediate of tryptophan catabolism and recycling it back to tryptophan. In effect, a host-parasite metabolic mosaic has evolved for tryptophan recycling.
| S-EPMC126876 | biostudies-literature
Project description:Secondary structure prediction is an important problem in RNA bioinformatics because knowledge of structure is critical to understanding the functions of RNA sequences. Significant improvements in prediction accuracy have recently been demonstrated though the incorporation of experimentally obtained structural information, for instance using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) mapping. However, such mapping data is currently available only for a limited number of RNA sequences. In this article, we present a method for extending the benefit of experimental mapping data in secondary structure prediction to homologous sequences. Specifically, we propose a method for integrating experimental mapping data into a comparative sequence analysis algorithm for secondary structure prediction of multiple homologs, whereby the mapping data benefits not only the prediction for the specific sequence that was mapped but also other homologs. The proposed method is realized by modifying the TurboFold II algorithm for prediction of RNA secondary structures to utilize basepairing probabilities guided by SHAPE experimental data when such data are available. The SHAPE-mapping-guided basepairing probabilities are obtained using the RSample method. Results demonstrate that the SHAPE mapping data for a sequence improves structure prediction accuracy of other homologous sequences beyond the accuracy obtained by sequence comparison alone (TurboFold II). The updated version of TurboFold II is freely available as part of the RNAstructure software package.
| S-EPMC5529333 | biostudies-literature
Project description:Dynamic causal modeling (DCM) is a methodological approach to study effective connectivity among brain regions. Based on a set of observations and a biophysical model of brain interactions, DCM uses a Bayesian framework to estimate the posterior distribution of the free parameters of the model (e.g. modulation of connectivity) and infer architectural properties of the most plausible model (i.e. model selection). When modeling electrophysiological event-related responses, the estimation of the model relies on the integration of the system of delay differential equations (DDEs) that describe the dynamics of the system. In this technical note, we compared two numerical schemes for the integration of DDEs. The first, and standard, scheme approximates the DDEs (more precisely, the state of the system, with respect to conduction delays among brain regions) using ordinary differential equations (ODEs) and solves it with a fixed step size. The second scheme uses a dedicated DDEs solver with adaptive step sizes to control error, making it theoretically more accurate. To highlight the effects of the approximation used by the first integration scheme in regard to parameter estimation and Bayesian model selection, we performed simulations of local field potentials using first, a simple model comprising 2 regions and second, a more complex model comprising 6 regions. In these simulations, the second integration scheme served as the standard to which the first one was compared. Then, the performances of the two integration schemes were directly compared by fitting a public mismatch negativity EEG dataset with different models. The simulations revealed that the use of the standard DCM integration scheme was acceptable for Bayesian model selection but underestimated the connectivity parameters and did not allow an accurate estimation of conduction delays. Fitting to empirical data showed that the models systematically obtained an increased accuracy when using the second integration scheme. We conclude that inference on connectivity strength and delay based on DCM for EEG/MEG requires an accurate integration scheme.
| S-EPMC5929904 | biostudies-other
Project description:Decreasing trends in measles mortality have been reported in recent years. However, such estimates of measles mortality have depended heavily on assumed regional measles case fatality risks (CFRs) and made little use of mortality data from low- and middle-income countries in general and India, the country with the highest measles burden globally, in particular.We constructed a dynamic model of measles transmission in India with parameters that were empirically inferred using spectral analysis from a time series of measles mortality extracted from the Million Death Study, an ongoing longitudinal study recording deaths across 2.4 million Indian households and attributing causes of death using verbal autopsy. The model was then used to estimate the measles CFR, the number of measles deaths, and the impact of vaccination in 2000-2015 among under-five children in India and in the states of Bihar and Uttar Pradesh (UP), two states with large populations and the highest numbers of measles deaths in India.We obtained the following estimated CFRs among under-five children for the year 2005: 0.63% (95% confidence interval (CI): 0.40-1.00%) for India as a whole, 0.62% (0.38-1.00%) for Bihar, and 1.19% (0.80-1.75%) for UP. During 2000-2015, we estimated that 607,000 (95% CI: 383,000-958,000) under-five deaths attributed to measles occurred in India as a whole. If no routine vaccination or supplemental immunization activities had occurred from 2000 to 2015, an additional 1.6 (1.0-2.6) million deaths for under-five children would have occurred across India.We developed a data- and model-driven estimation of the historical measles dynamics, CFR, and vaccination impact in India, extracting the periodicity of epidemics using spectral and coherence analysis, which allowed us to infer key parameters driving measles transmission dynamics and mortality.
| S-EPMC5550950 | biostudies-other
Project description:Large-scale observational data from citizen science efforts are becoming increasingly common in ecology, and researchers often choose between these and data from intensive local-scale studies for their analyses. This choice has potential trade-offs related to spatial scale, observer variance, and interannual variability. Here we explored this issue with phenology by comparing models built using data from the large-scale, citizen science USA National Phenology Network (USA-NPN) effort with models built using data from more intensive studies at Long Term Ecological Research (LTER) sites. We built statistical and process based phenology models for species common to each data set. From these models, we compared parameter estimates, estimates of phenological events, and out-of-sample errors between models derived from both USA-NPN and LTER data. We found that model parameter estimates for the same species were most similar between the two data sets when using simple models, but parameter estimates varied widely as model complexity increased. Despite this, estimates for the date of phenological events and out-of-sample errors were similar, regardless of the model chosen. Predictions for USA-NPN data had the lowest error when using models built from the USA-NPN data, while LTER predictions were best made using LTER-derived models, confirming that models perform best when applied at the same scale they were built. This difference in the cross-scale model comparison is likely due to variation in phenological requirements within species. Models using the USA-NPN data set can integrate parameters over a large spatial scale while those using an LTER data set can only estimate parameters for a single location. Accordingly, the choice of data set depends on the research question. Inferences about species-specific phenological requirements are best made with LTER data, and if USA-NPN or similar data are all that is available, then analyses should be limited to simple models. Large-scale predictive modeling is best done with the larger-scale USA-NPN data, which has high spatial representation and a large regional species pool. LTER data sets, on the other hand, have high site fidelity and thus characterize inter-annual variability extremely well. Future research aimed at forecasting phenology events for particular species over larger scales should develop models that integrate the strengths of both data sets.
| S-EPMC7378950 | biostudies-literature
Project description:Dynamic Causal Modeling (DCM) can be used to quantify cognitive function in individuals as effective connectivity. However, ambiguity among subjects in the number and location of discernible active regions prevents all candidate models from being compared in all subjects, precluding the use of DCM as an individual cognitive phenotyping tool. This paper proposes a solution to this problem by treating missing regions in the first-level analysis as missing data, and performing estimation of the time course associated with any missing region using one of four candidate methods: zero-filling, average-filling, noise-filling using a fixed stochastic process, or one estimated using expectation-maximization. The effect of this estimation scheme was analyzed by treating it as a preprocessing step to DCM and observing the resulting effects on model evidence. Simulation studies show that estimation using expectation-maximization yields the highest classification accuracy using a simple loss function and highest model evidence, relative to other methods. This result held for various dataset sizes and varying numbers of model choice. In real data, application to Go/No-Go and Simon tasks allowed computation of signals from the missing nodes and the consequent computation of model evidence in all subjects compared to 62 and 48 percent respectively if no preprocessing was performed. These results demonstrate the face validity of the preprocessing scheme and open the possibility of using single-subject DCM as an individual cognitive phenotyping tool.
| S-EPMC4082189 | biostudies-literature
Project description:RNA is a unique bio-macromolecule that can both record genetic information and perform biological functions in a variety of molecular processes, including transcription, splicing, translation, and even regulating protein function. RNAs adopt specific three-dimensional conformations to enable their functions. Experimental determination of high-resolution RNA structures using x-ray crystallography is both laborious and demands expertise, thus, hindering our comprehension of RNA structural biology. The computational modeling of RNA structure was a milestone in the birth of bioinformatics. Although computational modeling has been greatly improved over the last decade showing many successful cases, the accuracy of such computational modeling is not only length-dependent but also varies according to the complexity of the structure. To increase credibility, various experimental data were integrated into computational modeling. In this review, we summarize the experiments that can be integrated into RNA structure modeling as well as the computational methods based on these experimental data. We also demonstrate how computational modeling can help the experimental determination of RNA structure. We highlight the recent advances in computational modeling which can offer reliable structure models using high-throughput experimental data.
| S-EPMC7649352 | biostudies-literature
Project description:Deterministic dynamic causal modeling (DCM) for fMRI data is a sophisticated approach to analyse effective connectivity in terms of directed interactions between brain regions of interest. To date it is difficult to know if acquired fMRI data will yield precise estimation of DCM parameters. Focusing on parameter identifiability, an important prerequisite for research questions on directed connectivity, we present an approach inferring if parameters of an envisaged DCM are identifiable based on information from fMRI data. With the freely available "attention to motion" dataset, we investigate identifiability of two DCMs and show how different imaging specifications impact on identifiability. We used the profile likelihood, which has successfully been applied in systems biology, to assess the identifiability of parameters in a DCM with specified scanning parameters. Parameters are identifiable when minima of the profile likelihood as well as finite confidence intervals for the parameters exist. Intermediate epoch duration, shorter TR and longer session duration generally increased the information content in the data and thus improved identifiability. Irrespective of biological factors such as size and location of a region, attention should be paid to densely interconnected regions in a DCM, as those seem to be prone to non-identifiability. Our approach, available in the DCMident toolbox, enables to judge if the parameters of an envisaged DCM are sufficiently determined by underlying data without priors as opposed to primarily reflecting the Bayesian priors in a SPM-DCM. Assessments with the DCMident toolbox prior to a study will lead to improved identifiability of the parameters and thus might prevent suboptimal data acquisition. Thus, the toolbox can be used as a preprocessing step to provide immediate statements on parameter identifiability.
| S-EPMC4335185 | biostudies-other
Project description:BackgroundThe enormous database of microbial DNA generated from the Sargasso Sea metagenome provides a unique opportunity to locate genes participating in different biosynthetic pathways and to attempt to understand the relationship and evolution of those genes. In this article, an analysis of the Sargasso Sea metagenome is made with respect to the seven genes of the tryptophan pathway.ResultsAt least 5% of all the genes that are related to amino acid biosynthesis are tryptophan (trp) genes. Many contigs and scaffolds contain whole or split operons that are similar to previously analyzed trp gene organizations. Only two scaffolds discovered in this analysis possess a different operon organization of tryptophan pathway genes than those previously known. Many marine organisms lack an operon-type organization of these genes or have mini-operons containing only two trp genes. In addition, the trpB genes from this search reveal that the dichotomous division between trpB_1 and trpB_2 also occurs in organisms from the Sargasso Sea. One cluster was found to contain trpB sequences that were closely related to each other but distinct from most known trpB sequences.ConclusionThe data show that trp genes are widely dispersed within this metagenome. The novel organization of these genes and an unusual group of trpB_1 sequences that were found among some of these Sargasso Sea bacteria indicate that there is much to be discovered about both the reason for certain gene orders and the regulation of tryptophan biosynthesis in marine bacteria.
| S-EPMC2395257 | biostudies-literature
Project description:The malleability of intrinsically disordered proteins (IDPs) has generated great interest in understanding how their conformations respond to crowded cellular environments. Experiments can report gross properties such as fluorescence resonance energy transfer (FRET) efficiency but cannot resolve the conformational ensembles of IDPs and their interactions with macromolecular crowders. Computation can in principle provide the latter information but in practice has been hampered by the enormous expense for realistic modeling of IDPs and crowders and for sufficient conformational sampling. Here, taking advantage of a powerful method called FMAP (fast Fourier transform-based modeling of atomistic protein-crowder interactions), we computed how the conformational ensembles of three IDPs are modified in concentrated polyethylene glycol (PEG) 6000 solutions. We represented the IDPs at the all-atom level and the PEG molecules at a coarse-grained level and calculated the experimental observable, i.e., FRET efficiency. Whereas accounting for only steric repulsion of PEG led to overestimation of crowding effects, quantitative agreement with experimental data was obtained upon including mild IDP-PEG attraction. The present work demonstrates that realistic modeling of IDPs under crowded conditions for direct comparison with experiments is now achievable.
| S-EPMC6447477 | biostudies-literature