Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structured RNA”

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.

At least 487 records · Page 27

Improved Biofuel Production through Discovery and Engineering of Terpene Metabolism in Switchgrass

Project Objectives - Of the myriad specialized metabolites that plants deploy to adapt to environmental challenges, terpenes form the largest group. In many major crops, unique terpene blends serve as key stress defenses that directly impact plant fitness and yield. In addition, terpenes, such as bisabolene and pinene, are used for producing renewable biofuels. Essential to advancing a broader use of terpenes for biofuel feedstock engineering is a system-wide knowledge of the diverse biosynthetic machinery and defensive potential of often species-specific terpene blends. The proposed project would merge genome-wide enzyme discovery with comparative –omics, protein structural and plant microbiome studies to define the biosynthesis and stress-defensive functions of the switchgrass (Panicum virgatum) terpene network. These insights would be combined with developing and applying non-transgenic genome editing tools to design plants with desirable terpene blends for higher productivity and biofuel production on marginal lands. As a dedicated lignocellulosic feedstock for U.S. biofuel production with high net energy yield, stress tolerance, and available genome resources, switchgrass is well-suited for devising new avenues for biofuel production. Project Description – The diversity of plant terpene defenses is governed by species-specific families of terpene synthase (TPS) and cytochrome P450 monooxygenase (P450) enzymes. Mining of the switchgrass genome (genotype Alamo) identified ~100 TPS and P450 candidate genes, and combinatorial biochemical analysis of synthesized TPSs and P450s revealed more than a dozen enzymes with common and novel activities. In addition, several identified terpene metabolites and the corresponding transcripts were up-regulated in response to abiotic stressors. These findings demonstrate a unique switchgrass terpene network with probable importance to abiotic stress tolerance, thus providing a large chemical portfolio for optimizing crop resistance, yield, and biofuel composition. Leveraging these preliminary data, we propose to generate a genome-wide map of the switchgrass terpene metabolic network through multi-gene co-expression analyses that allow the efficient cross-validation of TPS and P450 functions. Key enzymes would further be applied to structure-function studies via X-ray protein crystallography, homology modeling and site-directed mutagenesis to gain mechanistic insight into the catalytic specificity of switchgrass terpene metabolism and provide gene and amino acid targets for genome editing. In tandem with terpene pathway discovery, system-wide metabolomics, transcriptomics and proteomics studies in switchgrass accessions of contrasting drought tolerance would define the role of switchgrass terpene metabolism in conferring abiotic stress resilience. Metabolic changes would be assessed in a combined approach of targeted (terpenes) and untargeted metabolite profiling using a high-resolution LC-MS/MS approach, differential gene expression analyses through multiplexed Illumina RNA sequencing, and quantitative analysis of high-priority pathway enzymes using multiple reaction monitoring (MRM). Drawing on these insights, knock-down/out mutants of stress-associated pathway nodes would be generated by optimizing transient virus-induced gene silencing (VIGS) and CRISPR/Cas9 systems under control of the Tobacco Rattle Virus (TRV). The resulting mutant lines would then be analyzed for stress susceptibility and the impact on the root microbiome to define gene functions in planta. Knowledge of terpene pathways, enzyme mechanisms and bioactivities would be applied to enhance switchgrass stress resilience and to tailor-make terpene blends for biofuel production. Here, TRV-enabled CRISPR/Cas9 genome editing, including allele-specific knock-out of redundant genes, engineering of enzyme specificity via structure-guided point mutations, and overexpression of terpene genes relevant to stress-protection or biofuel production, would be used to increase metabolic flux toward desired pathways. Broader Impacts - Integrating the system-wide discovery, mechanistic analysis and non-transgenic genome engineering of the switchgrass terpene network aligns the required steps to unlock the chemical potential of this important metabolite class to generate crops that are more resistant to stress and provide advanced biofuel production in light of rising climate pressures as foreseeable challenges for bioenergy crop cultivation. The proposed project would further offer interdisciplinary student training through active involvement in the project and integration of research concepts and outcomes into newly-developed graduate and undergraduate courses on Plant Biotechnology.

09 BIOMASS FUELS↗

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗

Nonhomologous end-joining uses distinct mechanisms to repair each strand of a double strand break

Abstract Nonhomologous end-joining repairs chromosomal double strand breaks, but it is unknown whether both strands are repaired by this pathway, and if one strand break’s repair path impacts the other. Here, we show that nonhomologous end-joining employs both of two a priori possible strategies. Strand breaks that can be directly ligated are joined near-simultaneously, with no effect of one strand break’s repair path on the other. More complex end structures require obligatorily ordered repair. The first strand to be repaired is used as template for repair of the opposite/second strand break, with the latter repair reaction occurring fastest when also coupled to nonhomologous end-joining. Enforced asymmetry in repair of each strand break can extend to the gap-filling polymerase employed, and whether the polymerases incorporate RNA or DNA. Our results resolve questions about pathway mechanism and identify a requirement for flexibility of the nonhomologous end-joining machinery for efficient repair of both strand breaks within diverse cellular double strand breaks.

Luthman, Adam J.↗

Self-assembly and condensation of intermolecular poly(UG) RNA quadruplexes

Abstract Poly(UG) or ‘pUG’ dinucleotide repeats are highly abundant sequences in eukaryotic RNAs. In Caenorhabditis elegans, pUGs are added to RNA 3′ ends to direct gene silencing within Mutator foci, a germ granule condensate. Here, we show that pUG RNAs efficiently self-assemble into gel condensates through quadruplex (G4) interactions. Short pUG sequences form right-handed intermolecular G4s (pUG G4s), while longer pUGs form left-handed intramolecular G4s (pUG folds). We determined a 1.05 Å crystal structure of an intermolecular pUG G4, which reveals an eight stranded G4 dimer involving 48 nucleotides, 7 different G and U quartet conformations, 7 coordinated potassium ions, 8 sodium ions and a buried water molecule. A comparison of the intermolecular pUG G4 and intramolecular pUG fold structures provides insights into the molecular basis for G4 handedness and illustrates how a simple dinucleotide repeat sequence can form complex structures with diverse topologies.

Biochemistry & Molecular Biology↗

A Faster Triphosphorylation Ribozyme

In support of the RNA world hypothesis, previous studies identified trimetaphosphate (Tmp) as a plausible energy source for RNA world organisms. In one of these studies, catalytic RNAs (ribozymes) that catalyze the triphosphorylation of RNA 5'-hydroxyl groups using Tmp were obtained by in vitro selection. One ribozyme (TPR1) was analyzed in more detail. TPR1 catalyzes the triphosphorylation reaction to a rate of 0.013 min -1 under selection conditions (50 mM Tmp, 100 mM MgCl 2 , 22°C). To identify a triphosphorylation ribozyme that catalyzes faster triphosphorylation, and possibly learn about its secondary structure TPR1 was subjected to a doped selection. The resulting ribozyme, TPR1e, contains seven mutations relative to TPR1, displays a previously unidentified duplex that constrains the ribozyme's structure, and reacts at a 24-fold faster rate than the parent ribozyme. Under optimal conditions (150 mM Tmp, 650 mM MgCl 2 , 40°C), the triphosphorylation rate of TRP1e reaches 6.8 min -1 .

Gregory F. Dolan↗

Origin and early evolution of photosynthesis

Photosynthesis was well-established on the earth at least 3.5 thousand million years ago, and it is widely believed that these ancient organisms had similar metabolic capabilities to modern cyanobacteria. This requires that development of two photosystems and the oxygen evolution capability occurred very early in the earth's history, and that a presumed phase of evolution involving non-oxygen evolving photosynthetic organisms took place even earlier. The evolutionary relationships of the reaction center complexes found in all the classes of currently existing organisms have been analyzed using sequence analysis and biophysical measurements. The results indicate that all reaction centers fall into two basic groups, those with pheophytin and a pair of quinones as early acceptors, and those with iron sulfur clusters as early acceptors. No simple linear branching evolutionary scheme can account for the distribution patterns of reaction centers in existing photosynthetic organisms, and lateral transfer of genetic information is considered as a likely possibility. Possible scenarios for the development of primitive reaction centers into the heterodimeric protein structures found in existing reaction centers and for the development of organisms with two linked photosystems are presented.

NASA Program Exobiology↗

Assembly of respiratory syncytial virus matrix protein lattice and its coordination with fusion glycoprotein trimers

Respiratory syncytial virus (RSV) is an enveloped, filamentous, negative-strand RNA virus that causes significant respiratory illness worldwide. RSV vaccines are available, however there is still significant need for research to support the development of vaccines and therapeutics against RSV and related Mononegavirales viruses. Individual virions vary in size, with an average diameter of ~130 nm and ranging from ~500 nm to over 10 µm in length. Though the general arrangement of structural proteins in virions is known, we use cryo-electron tomography and sub-tomogram averaging to determine the molecular organization of RSV structural proteins. We show that the peripheral membrane-associated RSV matrix (M) protein is arranged in a packed helical-like lattice of M-dimers. We report that RSV F glycoprotein is frequently observed as pairs of trimers oriented in an anti-parallel conformation to support potential interactions between trimers. Our sub-tomogram averages indicate the positioning of F-trimer pairs is correlated with the underlying M lattice. These results provide insight into RSV virion organization and may aid in the development of RSV vaccines and anti-viral targets.

59 BASIC BIOLOGICAL SCIENCES↗

Elucidation of odd-chain dicarboxylate metabolism in Acinetobacter baylyi and application to polyethylene upcycling

Polyethylene (PE) is a versatile polymer, but its end-of-life management is challenging due to its recalcitrant structure. We present a promising approach combining chemical degradation and bio-upcycling to convert postconsumer PE waste into a value-added bioproduct. Specifically, PE was degraded into acetic acid and C 4 –C 7 dicarboxylic acids by nitric acid. We then elucidated the catabolic pathways for glutarate (C 5 ) and pimelate (C 7 ) in the nonmodel bacterium Acinetobacter baylyi ADP1 through RNA sequencing, phenotyping, and enzymatic assays. Whole-genome sequencing of evolved isolates also identified a crucial IclR family transcriptional regulator, DcaS, which acts as a repressor of dicarboxylate metabolism. The reverse-engineered strain exhibited enhanced substrate utilization compared to the wild-type strain. Using rational metabolic engineering, the PE deconstruction products were bioconverted into the valuable chemical lycopene, highlighting the potential of this microbial chassis to produce value-added bioproducts from postconsumer PE waste, thus promoting a circular economy for plastics.

metabolic engineering↗

Measurement of atom resolvability in cryo-EM maps with Q-scores

Cryogenic electron microscopy (cryo-EM) maps are now at the point where resolvability of individual atoms can be achieved. However, resolvability is not necessarily uniform throughout the map. We introduce a quantitative parameter to characterize the resolvability of individual atoms in cryo-EM maps, the map Q-score. Q-scores can be calculated for atoms in proteins, nucleic acids, water, ligands and other solvent atoms, using models fitted to or derived from cryo-EM maps. Q-scores can also be averaged to represent larger features such as entire residues and nucleotides. Averaged over entire models, Q-scores correlate very well with the estimated resolution of cryo-EM maps for both protein and RNA. Here, assuming the models they are calculated from are well fitted to the map, Q-scores can be used as a measure of resolvability in cryo-EM maps at various scales, from entire macromolecules down to individual atoms. Q-score analysis of multiple cryo-EM maps of the same proteins derived from different laboratories confirms the reproducibility of structural features from side chains down to water and ion atoms.

59 BASIC BIOLOGICAL SCIENCES↗

A distinct RNA recognition mechanism governs Np 4 decapping by RppH

Dinucleoside tetraphosphates, often described as alarmones because their cellular concentration increases in response to stress, have recently been shown to function in bacteria as precursors to nucleoside tetraphosphate (Np 4 ) RNA caps. Removal of this cap is critical for initiating 5' end-dependent degradation of those RNAs, potentially affecting bacterial adaptability to stress; however, the predominant Np 4 decapping enzyme in proteobacteria, ApaH, is inactivated by the very conditions of disulfide stress that enable Np 4 -capped RNAs to accumulate to high levels. Here, in this study, we show that, in Escherichia coli cells experiencing such stress, the RNA pyrophosphohydrolase RppH assumes a leading role in decapping those transcripts, preferring them as substrates over their triphosphorylated and diphosphorylated counterparts. Unexpectedly, this enzyme recognizes Np 4 -capped 5' ends by a mechanism distinct from the one it uses to recognize other 5' termini, resulting in a one-nucleotide shift in substrate specificity. The unique manner in which capped substrates of this kind bind to the active site of RppH positions the δ-phosphate, rather than the β-phosphate, for hydrolytic attack, generating triphosphorylated RNA as the primary product of decapping. Consequently, a second RppH-catalyzed deprotection step is required to produce the monophosphorylated 5' terminus needed to stimulate rapid RNA decay. The unconventional manner in which RppH recognizes Np 4 -capped 5' ends and its differential impact on the rates at which such termini are deprotected as a prelude to RNA degradation could have major consequences for reprogramming gene expression during disulfide stress.

59 BASIC BIOLOGICAL SCIENCES↗

Ecological physiology of Synechococcus sp. strain SH-94-5, a naturally occurring cyanobacterium deficient in nitrate assimilation

Synechococcus sp. strain SH-94-5 is a nitrate assimilation-deficient cyanobacterium which was isolated from an ammonium-replete hot spring in central Oregon. While this clone could grow on ammonium and some forms of organic nitrogen as sole nitrogen sources, it could not grow on either nitrate or nitrite, even under conditions favoring passive diffusion. It was determined that this clone does not express functional nitrate reductase or nitrite reductase and that the lack of activity of either enzyme is not due to inactivation of the cyanobacterial nitrogen control protein NtcA. A few other naturally occurring cyanobacterial strains are also nitrate assimilation deficient, and phylogenetic analyses indicated that the ability to utilize nitrate has been independently lost at least four times during the evolutionary history of the cyanobacteria. This phenotype is associated with the presence of environmental ammonium, a negative regulator of nitrate assimilation gene expression, which may indicate that natural selection to maintain functional copies of nitrate assimilation genes has been relaxed in these habitats. These results suggest how the evolutionary fates of conditionally expressed genes might differ between environments and thereby effect ecological divergence and biogeographical structure in the microbial world.

Nitrates/metabolism↗

Phylogenetic mapping of bacterial morphology

The availability of a meaningful molecular phylogeny for bacteria provides a context for examining the historical significance of various developments in bacterial evolution. Herein, the classical morphological descriptions of selected members of the domain Bacteria are mapped upon the genealogical ancestry deduced from comparison of small-subunit rRNA sequences. For the species examined in this study, a distinct pattern emerges which indicates that the coccus shape has arisen and accumulated independently multiple times in separate lineages and typically survived as a persistent end-state morphology. At least two other morphologies persist but have evolved only once. This study demonstrates that although bacterial morphology is not useful in defining bacterial phylogeny, it is remarkably consistent with that phylogeny once it is known. An examination of the experimental evidence available for morphogenesis as well as microbial fossil evidence corroborates these findings. It is proposed that the accumulation of persistent morphologies is a result of the biophysical properties of peptidoglycan and their genetic control, and that an evolved body-plan strategy based on peptidoglycan may have been a fate-sealing step in the evolution of Bacteria. More generally, this study illustrates that significant evolutionary insights can be obtained by examining biological and biochemical data in the context of a reliable phylogenetic structure.

NASA Discipline Exobiology↗

Engineered mRNA–ribosome fusions for facile biosynthesis of selenoproteins

Ribosomes are often used in synthetic biology as a tool to produce desired proteins with enhanced properties or entirely new functions. However, repurposing ribosomes for producing designer proteins is challenging due to the limited number of engineering solutions available to alter the natural activity of these enzymes. In this study, we advance ribosome engineering by describing a novel strategy based on functional fusions of ribosomal RNA (rRNA) with messenger RNA (mRNA). Specifically, we create an mRNA–ribosome fusion called RiboU, where the 16S rRNA is covalently attached to selenocysteine insertion sequence (SECIS), a regulatory RNA element found in mRNAs encoding selenoproteins. When SECIS sequences are present in natural mRNAs, they instruct ribosomes to decode UGA codons as selenocysteine (Sec, U) codons instead of interpreting them as stop codons. This enables ribosomes to insert Sec into the growing polypeptide chain at the appropriate site. Our work demonstrates that the SECIS sequence maintains its functionality even when inserted into the ribosome structure. As a result, the engineered ribosomes RiboU interpret UAG codons as Sec codons, allowing easy and site-specific insertion of Sec in a protein of interest with no further modification to the natural machinery of protein synthesis. To validate this approach, we use RiboU ribosomes to produce three functional target selenoproteins in Escherichia coli by site-specifically inserting Sec into the proteins’ active sites. Overall, our work demonstrates the feasibility of creating functional mRNA–rRNA fusions as a strategy for ribosome engineering, providing a novel tool for producing Sec-containing proteins in live bacterial cells.

59 BASIC BIOLOGICAL SCIENCES↗

Evolution of major metabolic innovations in the Precambrian

A combination of information on the metabolic capabilities of prokaryotes with a composite phylogenetic tree depicting an overview of prokaryote evolution based on the sequences of bacterial ferredoxin, 2Fe-2S ferredoxin, 5S ribosomal RNA, and c-type cytochromes shows three zones of major metabolic innovation in the Precambrian. The middle of these, which reflects the genesis of oxygen-releasing photosynthesis and aerobic respiration, links metabolic innovations of the anaerobic stem on the one hand and, on the other, proliferation of aerobic bacteria and the symbiotic associations leading to the eukaryotes. Those pathways where information on the structure of the enzymes is known are especially considered. Halobacterium and Thermoplasma (archaebacteria) do not belong to a totally independent line on the basis of the composite tree but branch from the eukaryote cytoplasmic line.

Barnabas, J.↗

Carbonaceous Meteorites Contain a Wide Range of Extraterrestrial Nucleobases

All terrestrial organisms depend on nucleic acids (RNA and DNA), which use pyrimidine and purine nucleobases to encode genetic information. Carbon-rich meteorites may have been important sources of organic compounds required for the emergence of life on the early Earth; however, the origin and formation of nuc1eobases in meteorites has been debated for over 50 y. So far, the few nuc1eobases reported in meteorites are biologically common and lacked the structural diversity typical of other indigenous meteoritic organics. Here, we investigated the abundance and distribution of nucleobases and nucleobase analogs in formic acid extracts of 12 different meteorites by liquid chromatography-mass spectrometry. The Murchison and Lonewolf Nunataks 94102 meteorites contained a diverse suite of nucleobases, which included three unusual and terrestrially rare nucleobase analogs; purine, 2,6-diminopurine, and 6,8-diaminopurine. In a parallel experiment, we found an identical suite of nucleobases and nucleobase analogs generated in reactions of ammonium cyanide. Additionally, these nucleobase analoge were not detected above our parts-per-billion detection limits in any of the procedural blanks, control samples, a terrestrial soil sample, and an Antarctic ice sample. Our results demonstrate that the purines detected in meteorites are consistent with products of ammonium cyanide chemistry, which provides a plausible mechanism for their synthesis in the asteroid parent bodies, and strongly supports an extraterrestrial origin. The discovery of new nucleobase analogs in meteorites also expands the prebiotic molecular inventory available for constructing the first genetic molecules.

Callahan, Michael P.↗

Antisense targeting of decoy exons can reduce intron retention and increase protein expression in human erythroblasts

The decoy exon model has been proposed to regulate a subset of intron retention (IR) events involving predominantly larger introns (>1 kb). Splicing reporter studies have shown that decoy splice sites are essential for activity, suggesting that decoys act by engaging intron-terminal splice sites and competing with cross-intron interactions required for intron excision. The decoy model predicts that antisense oligonucleotides may be able to block decoy splice sites in endogenous pre-mRNA, thereby reducing IR and increasing productive gene expression. Indeed, we now demonstrate that targeting a decoy 5' splice site in the O-GlcNAc transferase (OGT) gene reduced IR from ~80% to ~20% in primary human erythroblasts, accompanied by increases in spliced OGT RNA and OGT protein expression. The remaining OGT IR was refractory to antisense treatment and might be mediated by independent mechanism(s). In contrast, other retained introns were strongly dependent on decoy function, since antisense targeting of decoy 5' splice sites greatly reduced (SNRNP70) or nearly eliminated (SF3B1) IR in two widely expressed splicing factors, and also greatly reduced IR in transcripts encoding the erythroid-specific structural protein, α-spectrin (SPTA1). These results show that modulating decoy exon function can dramatically alter IR and suggest that dynamic regulation of decoy exons could be a mechanism to fine-tune gene expression post-transcriptionally in many cell types.

59 BASIC BIOLOGICAL SCIENCES↗

Improved genome editing by an engineered CRISPR-Cas12a

CRISPR-Cas12a is an RNA-guided, programmable genome editing enzyme found within bacterial adaptive immune pathways. Unlike CRISPR-Cas9, Cas12a uses only a single catalytic site to both cleave target double-stranded DNA (dsDNA) (cis-activity) and indiscriminately degrade single-stranded DNA (ssDNA) (trans-activity). To investigate how the relative potency of cis- versus trans-DNase activity affects Cas12a-mediated genome editing, we first used structure-guided engineering to generate variants of Lachnospiraceae bacterium Cas12a that selectively disrupt trans-activity. The resulting engineered mutant with the biggest differential between cis - and trans -DNase activity in vitro showed minimal genome editing activity in human cells, motivating a second set of experiments using directed evolution to generate additional mutants with robust genome editing activity. Notably, these engineered and evolved mutants had enhanced ability to induce homology-directed repair (HDR) editing by 2–18-fold compared to wild-type Cas12a when using HDR donors containing mismatches with crRNA at the PAM-distal region. Finally, a site-specific reversion mutation produced improved Cas12a (iCas12a) variants with superior genome editing efficiency at genomic sites that are difficult to edit using wild-type Cas12a. This strategy establishes a pipeline for creating improved genome editing tools by combining structural insights with randomization and selection. The available structures of other CRISPR-Cas enzymes will enable this strategy to be applied to improve the efficacy of other genome-editing proteins.

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↗