Systems biology towards life in silico: mathematics of the control of living cells

Systems Biology is the science that aims to understand how biological function absent from macromolecules in isolation, arises when they are components of their system. Dedicated to the memory of Reinhart Heinrich, this paper discusses the origin and evolution of the new part of systems biology that relates to metabolic and signal-transduction pathways and extends mathematical biology so as to address postgenomic experimental reality. Various approaches to modeling the dynamics generated by metabolic and signal-transduction pathways are compared. The silicon cell approach aims to describe the intracellular network of interest precisely, by numerically integrating the precise rate equations that characterize the ways macromolecules’ interact with each other. The non-equilibrium thermodynamic or ‘lin–log’ approach approximates the enzyme rate equations in terms of linear functions of the logarithms of the concentrations. Biochemical Systems Analysis approximates in terms of power laws. Importantly all these approaches link system behavior to molecular interaction properties. The latter two do this less precisely but enable analytical solutions. By limiting the questions asked, to optimal flux patterns, or to control of fluxes and concentrations around the (patho)physiological state, Flux Balance Analysis and Metabolic/Hierarchical Control Analysis again enable analytical solutions. Both the silicon cell approach and Metabolic/Hierarchical Control Analysis are able to highlight where and how system function derives from molecular interactions. The latter approach has also discovered a set of fundamental principles underlying the control of biological systems. The new law that relates concentration control to control by time is illustrated for an important signal transduction pathway, i.e. nuclear hormone receptor signaling such as relevant to bone formation. It is envisaged that there is much more Mathematical Biology to be discovered in the area between molecules and Life.     

Analysis of a Critical Interaction within the Archaeal Box C/D Small Ribonucleoprotein Complex

In archaea and eukarya, box C/D ribonucleoprotein (RNP) complexes are responsible for 2′-O-methylation of tRNAs and rRNAs. The archaeal box C/D small RNP complex requires a small RNA component (sRNA) possessing Watson-Crick complementarity to the target RNA along with three proteins: L7Ae, Nop5p, and fibrillarin. Transfer of a methyl group from S-adenosylmethionine to the target RNA is performed by fibrillarin, which by itself has no affinity for the sRNA-target duplex. Instead, it is targeted to the site of methylation through association with Nop5p, which in turn binds to the L7Ae-sRNA complex. To understand how Nop5p serves as a bridge between the targeting and catalytic functions of the box C/D small RNP complex, we have employed alanine scanning to evaluate the interaction between the Pyrococcus horikoshii Nop5p domain and an L7Ae box C/D RNA complex. From these data, we were able to construct an isolated RNA-binding domain (Nop-RBD) that folds correctly as demonstrated by x-ray crystallography and binds to the L7Ae box C/D RNA complex with near wild type affinity. These data demonstrate that the Nop-RBD is an autonomously folding and functional module important for protein assembly in a number of complexes centered on the L7Ae-kinkturn RNP.Many biological RNAs require extensive modification to attain full functionality in the cell (1). Currently there are over 100 known RNA modification types ranging from small functional group substitutions to the addition of large multi-cyclic ring structures (2). Transfer RNA, one of many functional RNAs targeted for modification (3-6), possesses the greatest modification type diversity, many of which are important for proper biological function (7). Ribosomal RNA, on the other hand, contains predominantly two types of modified nucleotides: pseudouridine and 2′-O-methylribose (8). The crystal structures of the ribosome suggest that these modifications are important for proper folding (9, 10) and structural stabilization (11) in vivo as evidenced by their strong tendency to localize to regions associated with function (8, 12, 13). These roles have been verified biochemically in a number of cases (14), whereas newly emerging functional modifications are continually being investigated.Box C/D ribonucleoprotein (RNP)³ complexes serve as RNA-guided site-specific 2′-O-methyltransferases in both archaea and eukaryotes (15, 16) where they are referred to as small RNP complexes and small nucleolar RNPs, respectively. Target RNA pairs with the sRNA guide sequence and is methylated at the 2′-hydroxyl group of the nucleotide five bases upstream of either the D or D′ box motif of the sRNA (Fig. 1, star) (17, 18). In archaea, the internal C′ and D′ motifs generally conform to a box C/D consensus sequence (19), and each sRNA contains two guide regions ∼12 nucleotides in length (20). The bipartite architecture of the RNP potentially enables the complex to methylate two distinct RNA targets (21) and has been shown to be essential for site-specific methylation (22).Open in a separate windowFIGURE 1.Organization of the archaeal box C/D complex. The protein components of this RNP are L7Ae, Nop5p, and fibrillarin, which together bind a box C/D sRNA. The regions of the Box C/D sRNA corresponding to the conserved C, D, C′, and D′ boxes are labeled. The target RNA binds the sRNA through Watson-Crick pairing and is methylated by fibrillarin at the fifth nucleotide from the D/D′ boxes (star).In addition to the sRNA, the archaeal box C/D complex requires three proteins for activity (23): the ribosomal protein L7Ae (24, 25), fibrillarin, and the Nop56/Nop58 homolog Nop5p (Fig. 1). L7Ae binds to both box C/D and the C′/D′ motifs (26), which respectively comprise kink-turn (27) or k-loop structures (28), to initiate the assembly of the RNP (29, 30). Fibrillarin performs the methyl group transfer from the cofactor S-adenosylmethionine to the target RNA (31-33). For this to occur, the active site of fibrillarin must be positioned precisely over the specific 2′-hydroxyl group to be methylated. Although fibrillarin methylates this functional group in the context of a Watson-Crick base-paired helix (guide/target), it has little to no binding affinity for double-stranded RNA or for the L7Ae-sRNA complex (22, 26, 33, 34). Nop5p serves as an intermediary protein bringing fibrillarin to the complex through its association with both the L7Ae-sRNA complex and fibrillarin (22). Along with its role as an intermediary between fibrillarin and the L7Ae-sRNA complex, Nop5p possesses other functions not yet fully understood. For example, Nop5p self-dimerizes through a coiled-coil domain (35) that in most archaea and eukaryotic homologs includes a small insertion sequence of unknown function (36, 37). However, dimerization and fibrillarin binding have been shown to be mutually exclusive in Methanocaldococcus jannaschii Nop5p, potentially because of the presence of this insertion sequence (36). Thus, whether Nop5p is a monomer or a dimer in the active RNP is still under debate.In this study, we focus our attention on the Nop5p protein to investigate its interaction with a L7Ae box C/D RNA complex because both the fibrillarin-Nop5p and the L7Ae box C/D RNA interfaces are known from crystal structures (29, 35, 38). Individual residues on the surface of a monomeric form of Nop5p (referred to as mNop5p) (22) were mutated to alanine, and the effect on binding affinity for a L7Ae box C/D motif RNA complex was assessed through the use of electrophoretic mobility shift assays. These data reveal that residues important for binding cluster within the highly conserved NOP domain (39, 40). To demonstrate that this domain is solely responsible for the affinity of Nop5p for the preassembled L7Ae box C/D RNA complex, we expressed and purified it in isolation from the full Nop5p protein. The isolated Nop-RBD domain binds to the L7Ae box C/D RNA complex with nearly wild type affinity, demonstrating that the Nop-RBD is truly an autonomously folding and functional module. Comparison of our data with the crystal structure of the homologous spliceosomal hPrp31-15.5K protein-U4 snRNA complex (41) suggests the adoption of a similar mode of binding, further supporting a crucial role for the NOP domain in RNP complex assembly.     

Go contributes to olfactory reception in Drosophila melanogaster

Background  

Seven-transmembrane receptors typically mediate olfactory signal transduction by coupling to G-proteins. Although insect odorant receptors have seven transmembrane domains like G-protein coupled receptors, they have an inverted membrane topology and function as ligand-gated cation channels. Consequently, the involvement of cyclic nucleotides and G proteins in insect odor reception is controversial. Since the heterotrimeric G_oα subunit is expressed in Drosophila olfactory receptor neurons, we reasoned that G_o acts together with insect odorant receptor cation channels to mediate odor-induced physiological responses.     

Evidence for a Novel Marine Harmful Algal Bloom: Cyanotoxin (Microcystin) Transfer from Land to Sea Otters

“Super-blooms” of cyanobacteria that produce potent and environmentally persistent biotoxins (microcystins) are an emerging global health issue in freshwater habitats. Monitoring of the marine environment for secondary impacts has been minimal, although microcystin-contaminated freshwater is known to be entering marine ecosystems. Here we confirm deaths of marine mammals from microcystin intoxication and provide evidence implicating land-sea flow with trophic transfer through marine invertebrates as the most likely route of exposure. This hypothesis was evaluated through environmental detection of potential freshwater and marine microcystin sources, sea otter necropsy with biochemical analysis of tissues and evaluation of bioaccumulation of freshwater microcystins by marine invertebrates. Ocean discharge of freshwater microcystins was confirmed for three nutrient-impaired rivers flowing into the Monterey Bay National Marine Sanctuary, and microcystin concentrations up to 2,900 ppm (2.9 million ppb) were detected in a freshwater lake and downstream tributaries to within 1 km of the ocean. Deaths of 21 southern sea otters, a federally listed threatened species, were linked to microcystin intoxication. Finally, farmed and free-living marine clams, mussels and oysters of species that are often consumed by sea otters and humans exhibited significant biomagnification (to 107 times ambient water levels) and slow depuration of freshwater cyanotoxins, suggesting a potentially serious environmental and public health threat that extends from the lowest trophic levels of nutrient-impaired freshwater habitat to apex marine predators. Microcystin-poisoned sea otters were commonly recovered near river mouths and harbors and contaminated marine bivalves were implicated as the most likely source of this potent hepatotoxin for wild otters. This is the first report of deaths of marine mammals due to cyanotoxins and confirms the existence of a novel class of marine “harmful algal bloom” in the Pacific coastal environment; that of hepatotoxic shellfish poisoning (HSP), suggesting that animals and humans are at risk from microcystin poisoning when consuming shellfish harvested at the land-sea interface.     

Catabolism of 4-hydroxy-2-trans-nonenal by THP1 monocytes/macrophages and inactivation of carboxylesterases by this lipid electrophile

Oxidative stress in cells and tissues leads to the formation of an assortment of lipid electrophiles, such as the quantitatively important 4-hydroxy-2-trans-nonenal (HNE). Although this cytotoxic aldehyde is atherogenic the mechanisms involved are unclear. We hypothesize that elevated HNE levels can directly inactivate esterase and lipase activities in macrophages via protein adduction, thus generating a biochemical lesion that accelerates foam cell formation and subsequent atherosclerosis. In the present study we examined the effects of HNE treatment on esterase and lipase activities in human THP1 monocytes/macrophages at various physiological scales (i.e., pure recombinant enzymes, cell lysate, and intact living cells). The hydrolytic activities of bacterial and human carboxylesterase enzymes (pnbCE and CES1, respectively) were inactivated by HNE in vitro in a time- and concentration-dependent manner. In addition, so were the hydrolytic activities of THP1 cell lysates and intact THP1 monocytes and macrophages. A single lysine residue (Lys105) in recombinant CES1 was modified by HNE via a Michael addition reaction, whereas the lone reduced cysteine residue (Cys389) was found unmodified. The lipolytic activity of cell lysates and intact cells was more sensitive to the inhibitory effects of HNE than the esterolytic activity. Moreover, immunoblotting analysis using HNE antibodies confirmed that several cellular proteins were adducted by HNE following treatment of intact THP1 monocytes, albeit at relatively high HNE concentrations (>50 μM). Unexpectedly, in contrast to CES1, the treatment of a recombinant human CES2 with HNE enhanced its enzymatic activity ∼3-fold compared to untreated enzyme. In addition, THP1 monocytes/macrophages can efficiently metabolize HNE, and glutathione conjugation of HNE is responsible for ∼43% of its catabolism. The functional importance of HNE-mediated inactivation of cellular hydrolytic enzymes with respect to atherogenesis remains obscure, although this study has taken a first step toward addressing this important issue by examining the potential of HNE to inhibit this biochemical activity in a human monocyte/macrophage cell line.     

The secret life of α-catenin: moonlighting in morphogenesis

Cadherin-based intercellular adhesions are important determinants of proper tissue architecture. These adhesions must be both stable and dynamic to maintain tissue integrity as cells undergo morphogenetic movements during development. The role of α-catenin in this process has been vigorously debated due to conflicting in vitro and in vivo evidence regarding its molecular mechanism of action. Recent data supports the classical view that α-catenin facilitates actin attachments at adherens junctions, but also suggests that α-catenin may act as a force transducer, and may have additional roles in the cytoplasm. These multiple functions for α-catenin converge on the regulation of adhesion and may help to explain its stable yet dynamic nature.     

Estimating occupancy dynamics in an anuran assemblage from Louisiana,USA

Effective monitoring programs are designed to track changes in the distribution, occurrence, and abundance of species. We developed an extension of Royle and Kéry's (2007) single species model to estimate simultaneously temporal changes in probabilities of detection, occupancy, colonization, extinction, and species turnover using data on calling anuran amphibians, collected from 2002 to 2006 in the Lower Mississippi Alluvial Valley of Louisiana, USA. During our 5-year study, estimates of occurrence probabilities declined for all 12 species detected. These declines occurred primarily in conjunction with variation in estimates of local extinction probabilities (cajun chorus frog [Pseudacris fouquettei], spring peeper [P. crucifer], northern cricket frog [Acris crepitans], Cope's gray treefrog [Hyla chrysoscelis], green treefrog [H. cinerea], squirrel treefrog [H. squirella], southern leopard frog [Lithobates sphenocephalus], bronze frog [L. clamitans], American bullfrog [L. catesbeianus], and Fowler's toad [Anaxyrus fowleri]). For 2 species (eastern narrow-mouthed toad [Gastrophryne carolinensis] and Gulf Coast toad [Incilius nebulifer]), declines in occupancy appeared to be a consequence of both increased local extinction and decreased colonization events. The eastern narrow-mouthed toad experienced a 2.5-fold increase in estimates of occupancy in 2004, possibly because of the high amount of rainfall received during that year, along with a decrease in extinction and increase in colonization of new sites between 2003 and 2004. Our model can be incorporated into monitoring programs to estimate simultaneously the occupancy dynamics for multiple species that show similar responses to ecological conditions. It will likely be an important asset for those monitoring programs that employ the same methods to sample assemblages of ecologically similar species, including those that are rare. By combining information from multiple species to decrease the variance on estimates of individual species, our results are advantageous compared to single-species models. This feature enables managers and researchers to use an entire community, rather than just one species, as an ecological indicator in monitoring programs. © 2011 The Wildlife Society.     

Ab initio protein structure prediction

Steady progress has been made in the field of ab initio protein folding. A variety of methods now allow the prediction of low-resolution structures of small proteins or protein fragments up to approximately 100 amino acid residues in length. Such low-resolution structures may be sufficient for the functional annotation of protein sequences on a genome-wide scale. Although no consistently reliable algorithm is currently available, the essential challenges to developing a general theory or approach to protein structure prediction are better understood. The energy landscapes resulting from the structure prediction algorithms are only partially funneled to the native state of the protein. This review focuses on two areas of recent advances in ab initio structure prediction-improvements in the energy functions and strategies to search the caldera region of the energy landscapes.     

Brief debrisoquin administration to assess central dopaminergic function in children

Central dopaminergic (DA) function in children was assessed by monitoring plasma-free homovanillic acid (pHVA) levels after brief (18 hour) administration with debrisoquin sulfate, a peripherally active antihypertensive agent that blocks peripheral, but not central, HVA production. Brief debrisoquin administration resulted in marked reductions in pHVA in each of six patients studied. In five of the six patients, post-debrisoquin pHVA levels remained relatively stable over the six-hour period of observation. No significant cardiovascular or behavioral side effects of debrisoquin were observed. The brief debrisoquin administration method appears to be a safe, simple, and potentially valid peripheral technique for evaluating aspects of central dopaminergic function in children with neuropsychiatric disorders. Additional work is needed to further establish this method's validity and reliability.     

The role of lysyl oxidase and collagen crosslinking during sea urchin development

Lysyl oxidase, the only enzyme involved in collagen crosslinking, is shown to be present in embryos of the sea urchin Strongylocentrotus purpuratus. The enzyme specific activity increases over six-fold during development, showing the greatest rise during gastrulation and prism larva formation. The enzyme is inhibited by the specific inhibitor, beta-aminoproprionitrile (BAPN). Continuous BAPN treatment of S. purpuratus and Lytechinus pictus embryos from late cleavage stages onward increases the amount of noncrosslinked collagen present in prism larvae. When BAPN is added at the 128- or 256-cell stage it causes developmental arrest at the mesenchyme blastula stage. Embryos can be maintained in the arrested state for at least 96 h and will resume normal development and morphogenesis following BAPN removal. If BAPN is added after the mesenchyme blastula stage, it has little adverse effect on development; consequently nonspecific toxic effects of the drug are unlikely. The results suggest that lysyl oxidase and collagen crosslinking play a vital role in primary mesenchyme migration, gastrulation, and morphogenesis during sea urchin development and indicate that BAPN may be very useful in studying the extracellular matrix-cell interactions at the cellular and molecular level.