Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genome-scale modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

54 records · Page 3

Genome-scale Design and Engineering of Non-model Yeast Organisms for Production of Biofuels and Bioproducts

The overall goal of this project was to develop genome-scale design and engineering tools for two non-model yeast organisms including Rhodotorula toruloides and Issatchenkia orientalis to produce high-levels of fatty acids-derived products and organic acids, respectively. The project was performed between 9/15/2017 and 9/14/2024 (the last two-years were no-cost extensions). The team consisted of Huimin Zhao (Lead PI) and Christopher Rao (Co-PI) from the University of Illinois at Urbana-Champaign (UIUC), Costas Maranas (Co-PI) from the Pennsylvania State University, Joshua Rabinowitz (Co-PI) and Martin Wuhr (Co-PI) from Princeton University, and Yasuo Yoshikuni (Co-PI) from the DOE Joint Genome Institute. The team has made great progress in both tool development and fundamental understanding of these two non-model yeasts. In total, there were 40 research publications (one of them is still under review) and one patent application as well as numerous oral presentations.

60 APPLIED LIFE SCIENCES↗

Systems-Level Modeling for CRISPR-Based Metabolic Engineering

The CRISPR-Cas system has enabled the development of sophisticated, multigene metabolic engineering programs through the use of guide RNA-directed activation or repression of target genes. To optimize biosynthetic pathways in microbial systems, we need improved models to inform design and implementation of transcriptional programs. Recent progress has resulted in new modeling approaches for identifying gene targets and predicting the efficacy of guide RNA targeting. Genome-scale and flux balance models have successfully been applied to identify targets for improving biosynthetic production yields using combinatorial CRISPR-interference (CRISPRi) programs. Here, the advent of new approaches for tunable and dynamic CRISPR activation (CRISPRa) promises to further advance these engineering capabilities. Once appropriate targets are identified, guide RNA prediction models can lead to increased efficacy in gene targeting. Developing improved models and incorporating approaches from machine learning may be able to overcome current limitations and greatly expand the capabilities of CRISPR-Cas9 tools for metabolic engineering.

59 BASIC BIOLOGICAL SCIENCES↗

Metabolic flux and resource balance in the oleaginous yeast Rhodotorula toruloides

The yeast Rhodotorula toruloides is a promising bioproduction organism due to its high lipid yields and ability to grow on cheap and abundant substrates. Quantitative, systems-level assessment of its metabolic activity is accordingly merited. Resource-balance analysis (RBA) models capture not only reaction stoichiometry but also enzyme requirements for catalysis, providing valuable tools for understanding metabolic trade-offs and optimizing metabolic engineering strategies. Here, in this work, we present systems-level measurements of R. toruloides metabolic flux based on isotope tracing and metabolic flux analysis. In combination with new proteomic measurements, these flux data are used to parameterize a genome-scale resource balance model rtRBA. We find that S. cerevisiae and R. toruloides grow at nearly indistinguishable rates using similar biosynthetic but dramatically different central metabolic programs. R. toruloides consumes one-fifth as much glucose, which it metabolizes primarily via the pentose phosphate pathway and TCA cycle unlike primarily glycolysis in S. cerevisiae . Overall, across these two divergent yeasts, protein abundances aligned more closely than metabolic flux. Resource balance modeling of these metabolic programs predicts superior theoretical yields but lower productivities in R. toruloides than S. cerevisiae for industrial chemicals, highlighting the value of rapid glucose uptake for productivity but respiratory metabolism for yields.

60 APPLIED LIFE SCIENCES↗

Coupling flux balance analysis with reactive transport modeling through machine learning for rapid and stable simulation of microbial metabolic switching

Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.

59 BASIC BIOLOGICAL SCIENCES↗

Deciphering metabolic differentiation during Bacillus subtilis sporulation

Abstract The bacteriumBacillus subtilisundergoes asymmetric cell division during sporulation, producing a mother cell and a smaller forespore connected by the SpoIIQ-SpoIIIA (or Q-A) channel. The two cells differentiate metabolically, and the forespore becomes dependent on the mother cell for essential building blocks. Here, we investigate the metabolic interactions between mother cell and forespore using genome-scale metabolic and expression models as well as experiments. Our results indicate that nucleotides are synthesized in the mother cell and transported in the form of nucleoside di- or tri-phosphates to the forespore via the Q-A channel. However, if the Q-A channel is inactivated later in sporulation, then glycolytic enzymes can form an ATP and NADH shuttle, providing the forespore with energy and reducing power. Our integrated in silico and in vivo approach sheds light into the intricate metabolic interactions underlying cell differentiation inB. subtilis, and provides a foundation for future studies of metabolic differentiation.

Science & Technology - Other Topics↗

Comparative Performance Evaluation of Large Language Models for Extracting Molecular Interactions and Pathway Knowledge

Understanding the interactions and regulatory relationships among biomolecules is essential for deciphering complex biological systems and elucidating the mechanisms behind diverse biological functions. Traditionally, the collection of such molecular interaction data has relied on expert curation, a process that is both time-consuming and labor-intensive. To address these limitations, this study explores the use of large language models (LLMs) to automate the genome-scale extraction of molecular interaction knowledge. Here, we evaluate the performance of various LLMs on key biological tasks, including the identification of protein-protein interactions, detection of genes associated with pathways influenced by low-dose radiation, and inference of gene regulatory relationships. Our findings demonstrate that larger LLMs tend to perform better, particularly in extracting intricate gene and protein interactions. Despite their strengths, these models face challenges in recognizing functionally diverse gene groups and highly correlated regulatory relationships. Through a comprehensive analysis using established molecular interaction and pathway databases, we show that LLMs possess the potential to identify relevant biomolecules and predict their interactions, offering valuable insights and marking a significant step toward AI-driven biological knowledge discovery.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Application of functional genomics for domestication of novel non-model microbes

Abstract With the expansion of domesticated microbes producing biomaterials and chemicals to support a growing circular bioeconomy, the variety of waste and sustainable substrates that can support microbial growth and production will also continue to expand. The diversity of these microbes also requires a range of compatible genetic tools to engineer improved robustness and economic viability. As we still do not fully understand the function of many genes in even highly studied model microbes, engineering improved microbial performance requires introducing genome-scale genetic modifications followed by screening or selecting mutants that enhance growth under prohibitive conditions encountered during production. These approaches include adaptive laboratory evolution, random or directed mutagenesis, transposon-mediated gene disruption, or CRISPR interference (CRISPRi). Although any of these approaches may be applicable for identifying engineering targets, here we focus on using CRISPRi to reduce the time required to engineer more robust microbes for industrial applications. One-Sentence Summary The development of genome scale CRISPR-based libraries in new microbes enables discovery of genetic factors linked to desired traits for engineering more robust microbial systems.

59 BASIC BIOLOGICAL SCIENCES↗

Beyond Component Optimization: Systems Level Biodesign for Lanthanide Recovery

Global demand for lanthanides (Ln) is projected to rise sharply over the next decade, while geographically concentrated supply chains and the low concentrations and matrix complexity of secondary feedstocks limit the reach of conventional hydro- and pyrometallurgical separation. Engineered biological systems offer a selective, low-energy alternative, and component-level advances in Ln-binding proteins, AI-designed selective scaffolds, and cell-surface display platforms now rival synthetic chelators in affinity and selectivity. These components, however, remain functionally isolated. Currently, there are no engineered chassis coupling recognition, intracellular trafficking, accumulation, and controlled release into an end-to-end pipeline. Here, we outline how new biodesign strategies and chassis selection must move beyond bioleaching to encompass the full recovery pathway. Achieving this requires integrating AI/ML-guided design, genome-scale build tools, high-throughput phenotyping, and biophysical transport modeling within a Design–Build–Test–Learn cycle tuned to recognition, trafficking, accumulation, and release.

Biodesign↗

A constraint-based framework for exploring the impact of multireaction dependencies on metabolic functions

Abstract Metabolism operates under physico-chemical constraints that result in multireaction dependencies. Understanding how multireaction dependencies affect metabolic phenotypes remains challenging, hindering their biotechnological applications. Here, we propose the concept of a forcedly balanced complex that allows to efficiently determine the effects of specific multireaction dependencies on metabolic network functions in constrained-based models. Using this concept, we found that the fraction of multireaction dependencies induced by forcedly balanced complexes in genome-scale metabolic networks followed power law with exponential cut-off. We identified forcedly balanced complexes that are lethal in cancer but have little effect on growth in healthy tissue models. In addition, these forcedly balanced complexes are largely specific to models of particular cancer types. Therefore, multireaction dependencies resulting from forced balancing of complexes represent an innovative means to control cancers that, we argue, can be implemented via transporter engineering. The presented constraint-based approaches pave the way for using multireaction dependencies in metabolic engineering for diverse biotechnological applications.

Küken, Anika↗

Detectability of Varied Hybridization Scenarios Using Genome-Scale Hybrid Detection Methods

Hybridization events complicate the accurate reconstruction of phylogenies, as they lead to patterns of genetic heritability that are unexpected under traditional, bifurcating models of species trees. This phenomenon has led to the development of methods to infer these varied hybridization events, both methods that reconstruct networks directly, as well as summary methods that predict individual hybridization events from a subset of taxa. However, a lack of empirical comparisons between methods – especially those pertaining to large networks with varied hybridization scenarios – hinders their practical use. Here, we provide a comprehensive review of popular summary methods: TICR, MSCquartets, HyDe, Patterson’s D-Statistic (ABBA-BABA), D3, and Dp. TICR and MSCquartets are based on quartet concordance factors gathered from gene tree topologies and HyDe, Patterson’s D-Statistic, D3, and Dp use site pattern frequencies to identify hybridization events between sets of three taxa. We then use simulated data to address questions of method accuracy and ideal use scenarios by testing methods against complex networks which depict gene flow events that differ in depth (timing), quantity (single vs. multiple, overlapping hybridizations), and rate of gene flow (γ). We find that deeper or multiple hybridization events may introduce noise and weaken the signal of hybridization, leading to higher relative false negative rates across all methods. Despite some forms of hybridization eluding quartet-based detection methods, MSCquartets displays high precision in most scenarios. While HyDe results in high false negative rates when tested on hybridizations involving extinct or unsampled ghost lineages, HyDe is the only method able to identify the direction of hybridization, distinguishing the source parental lineages from recipient hybrid lineages. Lastly, we test the methods on a dataset of ultraconserved elements from the bee subfamily Nomiinae, finding possible hybridization events between clades which correspond to regions of poor support in the species tree estimated in a previous study.

Bjorner, Marianne B.↗

Competition and cooperation: The plasticity of bacterial interactions across environments

Bacteria live in diverse communities, forming complex networks of interacting species. A central question in bacterial ecology is whether species engage in cooperative or competitive interactions. But this question often neglects the role of the environment. Here, we use genome-scale metabolic networks from two different open-access collections (AGORA and CarveMe) to assess pairwise interactions of different microbes in varying environmental conditions (provision of different environmental compounds). By computationally simulating thousands of environments for 10,000 pairs of bacteria from each collection, we found that most pairs were able to both compete and cooperate depending on the availability of environmental resources. This modeling approach allowed us to determine commonalities between environments that could facilitate the potential for cooperation or competition between a pair of species. Namely, cooperative interactions, especially obligate, were most common in less diverse environments. Further, as compounds were removed from the environment, we found interactions tended to degrade towards obligacy. However, we also found that on average at least one compound could be removed from an environment to switch the interaction from competition to facultative cooperation or vice versa. Together our approach indicates a high degree of plasticity in microbial interactions in response to the availability of environmental resources.

Biochemistry & Molecular Biology↗

The Transcriptional Response of Diverse Saccharomyces Cerevisiae Strains to Simulated Microgravity

Spaceflight imposes multiple stresses on biological systems resulting in genome-scale adaptations. Understanding these adaptations and their underlying molecular mechanisms is important to clarifying and reducing the risks associated with spaceflight. One such risk is infection by microbes present in spacecraft and their associated systems and inhabitants. This risk is compounded by results suggesting that some microbes may exhibit increased virulence after exposure to spaceflight conditions. The yeast, S. cerevisiae, is a powerful microbial model system, and its response to spaceflight has been studied for decades. However, to date, these studies have utilized common lab strains. Yet studies on trait variation in S. cerevisiae demonstrate that these lab strains are not representative of wild yeast and instead respond to environmental stimuli in an atypical manner. Thus, it is not clear how transferable these results are to the wild S. cerevisiae strains likely to be encountered during spaceflight. To determine if diverse S. cerevisiae strains exhibit a conserved response to simulated microgravity, we will utilize a collection of 100 S. cerevisiae strains isolated from clinical, environmental and industrial settings. We will place selected S. cerevisiae strains in simulated microgravity using a high-aspect rotating vessel (HARV) and document their transcriptional response by RNA-sequencing and quantify similarities and differences between strains. Our research will have a strong impact on the understanding of how genetic diversity of microorganisms effects their response to spaceflight, and will serve as a platform for further studies.

yeast↗

The Transcriptional Response of Diverse Saccharomyces Cerevisiae Strains to Simulated Microgravity

Spaceflight imposes multiple stresses on biological systems resulting in genome-scale adaptations. Understanding these adaptations and their underlying molecular mechanisms is important to clarifying and reducing the risks associated with spaceflight. One such risk is infection by microbes present in spacecraft and their associated systems and inhabitants. This risk is compounded by results suggesting that some microbes may exhibit increased virulence after exposure to spaceflight conditions. The yeast, S. cerevisiae, is a powerful microbial model system, and it's response to spaceflight has been studied for decades. However, to date, these studies have utilized common lab strains. Yet studies on trait variation in S. cerevisiae demonstrate that these lab strains are not representative of wild yeast and instead respond to environmental stimuli in an atypical manner. Thus, it is not clear how transferable these results are to the wild S. cerevisiae strains likely to be encountered during spaceflight. To determine if diverse S. cerevisiae strains exhibit a conserved response to simulated microgravity, we will utilize a collection of 100 S. cerevisiae strains isolated from clinical, environmental and industrial settings. We will place selected S. cerevisiae strains in simulated microgravity using a high-aspect rotating vessel (HARV) and document their transcriptional response by RNA-sequencing and quantify similarities and differences between strains. Our research will have a strong impact on the understanding of how genetic diversity of microorganisms effects their response to spaceflight, and will serve as a platform for further studies.

genetic diversity↗

The Transcriptional Response of Diverse Saccharomyces cerevisiae Strains to Simulated Microgravity

Spaceflight imposes multiple stresses on biological systems resulting in genome-scale adaptations. Understanding these adaptations and their underlying molecular mechanisms is important to clarifying and reducing the risks associated with spaceflight. One such risk is infection by microbes present in spacecraft and their associated systems and inhabitants. This risk is compounded by results suggesting that some microbes may exhibit increased virulence after exposure to spaceflight conditions. The yeast, S. cerevisiae, is a powerful microbial model system, and it's response to spaceflight has been studied for decades. However, to date, these studies have utilized common lab strains. Yet studies on trait variation in S. cerevisiae demonstrate that these lab strains are not representative of wild yeast and instead respond to environmental stimuli in an a typical manner. Thus, it is not clear how transferable these results are to the wild S. cerevisiae strains likely to be encountered during spaceflight. To determine if diverse S. cerevisiae strains exhibit a conserved response to simulated microgravity, we will utilize a collection of 100 S. cerevisiae strains isolated from clinical, environmental and industrial settings. We will place selected S. cerevisiae strains in simulated microgravity using a high-aspect rotating vessel (HARV) and document their transcriptional response by RNA-sequencing and quantify similarities and differences between strains. Our research will have a strong impact on the understanding of how genetic diversity of microorganisms effects their response to spaceflight, and will serve as a platform for further studies.

spaceflight↗

Revealing systematic changes in the transcriptome during the transition from exponential growth to stationary phase

ABSTRACT The composition of bacterial transcriptomes is determined by the transcriptional regulatory network (TRN). The TRN regulates the transition from one physiological state to another. Here, we use independent component analysis to monitor the composition of the transcriptome during the transition from the exponential growth phase to the stationary phase. With Escherichia coli K-12 MG1655 as a model strain, we trigger the transition using carbon, nitrogen, and sulfur starvation. We find that (i) the transition to the stationary phase accompanies common transcriptome changes, including increased stringent responses and reduced production of cellular building blocks and energy regardless of the limiting element; (ii) condition-specific changes are strongly associated with transcriptional regulators ( e.g. , Crp, NtrC, CysB, Cbl) responsible for metabolizing the limiting element; and (iii) the shortage of each limiting element differentially affects the production of amino acids and extracellular polymers. This study demonstrates how the combination of genome-scale datasets and new data analytics reveals the fundamental characteristics of a key transition in the life cycle of bacteria. IMPORTANCE Nutrient limitations are critical environmental perturbations in bacterial physiology. Despite its importance, a detailed understanding of how bacterial transcriptomes are adjusted has been limited. By utilizing independent component analysis (ICA) to decompose transcriptome data, this study reveals key regulatory events that enable bacteria to adapt to nutrient limitations. The findings not only highlight common responses, such as the stringent response, but also condition-specific regulatory shifts associated with carbon, nitrogen, and sulfur starvation. The insights gained from this work advance our knowledge of bacterial physiology, gene regulation, and metabolic adaptation.

Lim, Hyun Gyu (ORCID:0000000204692388)↗

Transporter annotations are holding up progress in metabolic modeling

Mechanistic, constraint-based models of microbial isolates or communities are a staple in the metabolic analysis toolbox, but predictions about microbe-microbe and microbe-environment interactions are only as good as the accuracy of transporter annotations. A number of hurdles stand in the way of comprehensive functional assignments for membrane transporters. These include general or non-specific substrate assignments, ambiguity in the localization, directionality and reversibility of a transporter, and the many-to-many mapping of substrates, transporters and genes. In this perspective, we summarize progress in both experimental and computational approaches used to determine the function of transporters and consider paths forward that integrate both. Investment in accurate, high-throughput functional characterization is needed to train the next-generation of predictive tools toward genome-scale metabolic network reconstructions that better predict phenotypes and interactions. More reliable predictions in this domain will benefit fields ranging from personalized medicine to metabolic engineering to microbial ecology.

Casey, John↗

Bacterial population-level trade-offs between drought tolerance and resource acquisition traits impact decomposition

Microbes drive fundamental ecosystem processes, such as decomposition. Environmental stressors are known to affect microbes, their fitness, and the ecosystem functions that they perform; yet, understanding the causal mechanisms behind this influence has been difficult. We used leaf litter on soil surface as a model in situ system to assess changes in bacterial genomic traits and decomposition rates for 18 months with drought as a stressor. We hypothesized that genome-scale trade-offs due to investment in stress tolerance traits under drought reduce the capacity for bacterial populations to carry out decomposition, and that these population-level trade-offs scale up to impact emergent community traits, thereby reducing decomposition rates. We observed drought tolerance mechanisms that were heightened in bacterial populations under drought, identified as higher gene copy numbers in metagenome-assembled genomes. A subset of populations under drought had reduced carbohydrate-active enzyme genes that suggested—as a trade-off—a decline in decomposition capabilities. These trade-offs were driven by community succession and taxonomic shifts as distinct patterns appeared in populations. We show that trait–trade-offs in bacterial populations under drought could scale up to reduce overall decomposition capabilities and litter decay rates. Using a trait-based approach to assess the population ecology of soil bacteria, we demonstrate genome-level trade-offs in response to drought with consequences for decomposition rates.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling the Activity of Single Genes

The central dogma of molecular biology states that information is stored in DNA, transcribed to messenger RNA (mRNA) and then translated into proteins. This picture is significantly augmentated when we consider the action of certain proteins in regulating transcription. These transcription factors provide a feedback pathway by which genes can regulate one another's expression as mRNA and then as protein. To review: DNA, RNA and proteins have different functions. DNA is the molecular storehouse of genetic information. When cells divide, the DNA is replicated, so that each daughter cell maintains the same genetic information as the mother cell. RNA acts as a go-between from DNA to proteins. Only a single copy of DNA is present, but multiple copies of the same piece of RNA may be present, allowing cells to make huge amounts of protein. In eukaryotes (organisms with a nucleus), DNA is found in the nucleus only. RNA is copied in the nucleus then translocates(moves) outside the nucleus, where it is transcribed into proteins. Along the way, the RNA may be spliced, i.e., may have pieces cut out. RNA then attaches to ribosomes and is translated to proteins. Proteins are the machinery of the cell other than DNA and RNA, all the complex molecules of the cell are proteins. Proteins are specialized machines, each of which fulfills its own task, which may be transporting oxygen, catalyzing reactions, or responding to extracellular signals, just to name a few. One of the more interesting functions a protein may have is binding directly or indirectly to DNA to perform transcriptional regulation, thus forming a closed feedback loop of gene regulation. The structure of DNA and the central dogma were understood in the 50s; in the early 80s it became possible to make arbitrary modifications to DNA and use cellular machinery to transcribe and translate the resulting genes; more recently, genomes (i.e., the complete DNA sequence) of many organisms have been sequenced. This large-scale sequencing began with simple organisms, viruses and bacteria, progressed to eukaryotes such as yeast, and more recently (1998) progressed to a multi-cellular animal, the nematode Caenorhabditis elegans. Sequencers have now moved on to the fruit fly Drosophila melanogaster, whose sequence is slated for completion by the end of 1999. The human genome project is expected to determine the complete sequence of all 3 billion bases of human DNA within the next five years. In the wake of genome-scale sequencing, further instrumentation is being developed to assay gene expression and function on a comparably large scale. Much of the work in computational biology focuses on computational tools used in sequencing, finding genes that are related to a particular gene, finding which parts of the DNA code for proteins and which do not, understanding what proteins will be formed from a given length of DNA, predicting how the proteins will fold from a one-dimensional structure into a three dimensional structure, and so on. Much less computational work has been done regarding the function of proteins. One reason for this is that different proteins function very differently, and so work on protein function is very specific to certain classes of proteins. There are, for example, proteins such enzymes that catalyze various intracellular reactions, receptors that respond to extracellular signals and ion channels that regulate the flow of charged particles into and out of the cell. In this chapter, we will consider a particular class of proteins called transcription factors(TFs), which are responsible for regulating when a certain gene is expressed in a certain cell, which cells it is express in, and how much is expressed. Understanding these processes will involve developing a deeper understanding of transcription, translation, and the cellular processes that control those processes. All of these elements fall under the aegis of gene regulation or more narrowly transcriptional regulation. Some of the key questions in gene regulation are: What genes are expressed in a certain cell at a certain time? How does gene expression differ from cell to cell in a multicellular organism? Which proteins act as transcription factors, i.e., are important in regulating gene expression? From questions like these, we hope to understand which genes are important for various macroscopic processes. Nearly all of the cells of a multicellular organism contain the same DNA. Yet this same genetic information yields a large number of different cell types. The fundamental difference between a neuron and a liver cell, for example, is which genes are expressed. Thus understanding gene regulation is an important step in understanding development. Furthermore, understanding the usual genes that are expressed in cells may give important clues about various diseases. Some diseases, such as sickle cell anemia and cystic fibrosis, are caused by defects in single, non-regulatory genes; others, such as certain cancers, are caused when the cellular control circuitry malfunctions - an understanding of these diseases will involve pathways of multiple interacting gene products. There are numerous challenges in the area of understanding and modeling gene regulation. First and foremost, biologists would like to develop a deeper understanding of the processes involved, including which genes and families of genes are important, how they interact, etc. From a computation point of view, there has been embarrassingly little work done. In this chapter there are many areas in which we can phrase meaningful, non-trivial computational questions, but questions that have not been addressed. Some of these are purely computational (what is a good algorithm for dealing with a model of type X) and others are more mathematical (given a system with certain characteristics, what sort of model can one use? How does one find biochemical parameters from system-level behavior using as few experiments as possible?). In addition to biological and algorithmic problems, there is also the ever-present issue of theoretical biology - what general principles can be derived from these systems, what can one do with models other than just simulate time-courses, what can be deduced about a class of systems without knowing all the details? The fundamental challenge to computationalists and theorists is to add value to the biology - to use models, modeling techniques and algorithms to understand the biology in new ways.

Mjolsness, Eric↗