Oscillatory genetic circuits: new designs and mathematical models

Introduction and motivation: A wide variety of organisms have developed in-ternal biomolecular clocks in order to adapt to cyclic changes of the environment. Clock operation involves genetic networks. These genetic networks have to be mod¬eled in order to understand the underlying mechanism of oscillations and to design new synthetic cellular clocks. This doctoral thesis has resulted in two contributions to the fields of genetic clocks and systems and synthetic biology, generally. The first contribution is a new genetic circuit model that exhibits an oscillatory behav¬ior through catalytic RNA molecules. The second and major contribution is a new genetic circuit model demonstrating that a repressor molecule acting on the positive feedback of a self-activating gene produces reliable oscillations. First contribution: A new model of a synthetic genetic oscillator based on a typical two-gene motif with one positive and one negative feedback loop is pre¬sented. The originality is that the repressor is a catalytic RNA molecule rather than a protein or a non-catalytic RNA molecule. This catalytic RNA is a ribozyme that acts post-transcriptionally by binding to and cleaving target mRNA molecules. This genetic clock involves just two genes, a mRNA and an activator protein, apart from the ribozyme. Parameter values that produce a circadian period in both determin¬istic and stochastic simulations have been chosen as an example of clock operation. The effects of the stochastic fluctuations are quantified by a period histogram and autocorrelation function. The conclusion is that catalytic RNA molecules can act as repressor proteins and simplify the design of genetic oscillators. Second and major contribution: It is demonstrated that a self-activating gene in conjunction with a simple negative interaction can easily produce robust matically validated. This model is comprised of two clearly distinct parts. The first is a positive feedback created by a protein that binds to the promoter of its own gene and activates the transcription. The second is a negative interaction in which a repressor molecule prevents this protein from binding to its promoter. A stochastic study shows that the system is robust to noise. A deterministic study identifies that the oscillator dynamics are mainly driven by two types of biomolecules: the protein, and the complex formed by the repressor and this protein. The main conclusion of this study is that a simple and usual negative interaction, such as degradation, se¬questration or inhibition, acting on the positive transcriptional feedback of a single gene is a sufficient condition to produce reliable oscillations. One gene is enough and the positive transcriptional feedback signal does not need to activate a second repressor gene. At the genetic level, this means that an explicit negative feedback loop is not necessary. Unlike many genetic oscillators, this model needs neither cooperative binding reactions nor the formation of protein multimers. Applications and future research directions: Recently, RNA molecules have been found to play many new catalytic roles. The first oscillatory genetic model proposed in this thesis uses ribozymes as repressor molecules. This could provide new synthetic biology design principles and a better understanding of cel¬lular clocks regulated by RNA molecules. The second genetic model proposed here involves only a repression acting on a self-activating gene and produces robust oscil¬lations. Unlike current two-gene oscillators, this model surprisingly does not require a second repressor gene. This result could help to clarify the design principles of cellular clocks and constitute a new efficient tool for engineering synthetic genetic oscillators. Possible follow-on research directions are: validate models in vivo and in vitro, research the potential of second model as a genetic memory, investigate new genetic oscillators regulated by non-coding RNAs and design a biosensor of positive feedbacks in genetic networks based on the operation of the second model Resumen Introduccion y motivacion: Una amplia variedad de organismos han desarro-llado relojes biomoleculares internos con el fin de adaptarse a los cambios ciclicos del entorno. El funcionamiento de estos relojes involucra redes geneticas. El mo delado de estas redes geneticas es esencial tanto para entender los mecanismos que producen las oscilaciones como para diseiiar nuevos circuitos sinteticos en celulas. Esta tesis doctoral ha dado lugar a dos contribuciones dentro de los campos de los circuitos geneticos en particular, y biologia de sistemas y sintetica en general. La primera contribucion es un nuevo modelo de circuito genetico que muestra un comportamiento oscilatorio usando moleculas de ARN cataliticas. La segunda y principal contribucion es un nuevo modelo de circuito genetico que demuestra que una molecula represora actuando sobre el lazo de un gen auto-activado produce oscilaciones robustas. Primera contribucion: Es un nuevo modelo de oscilador genetico sintetico basado en una tipica red genetica compuesta por dos genes con dos lazos de retroa-limentacion, uno positivo y otro negativo. La novedad de este modelo es que el represor es una molecula de ARN catalftica, en lugar de una protefna o una molecula de ARN no-catalitica. Este ARN catalitico es una ribozima que actua despues de la transcription genetica uniendose y cortando moleculas de ARN mensajero (ARNm). Este reloj genetico involucra solo dos genes, un ARNm y una proteina activadora, aparte de la ribozima. Como ejemplo de funcionamiento, se han escogido valores de los parametros que producen oscilaciones con periodo circadiano (24 horas) tanto en simulaciones deterministas como estocasticas. El efecto de las fluctuaciones es-tocasticas ha sido cuantificado mediante un histograma del periodo y la función de auto-correlacion. La conclusion es que las moleculas de ARN con propiedades cataliticas pueden jugar el misnio papel que las protemas represoras, y por lo tanto, simplificar el diseno de los osciladores geneticos. Segunda y principal contribucion: Es un nuevo modelo de oscilador genetico que demuestra que un gen auto-activado junto con una simple interaction negativa puede producir oscilaciones robustas. Este modelo ha sido estudiado y validado matematicamente. El modelo esta compuesto de dos partes bien diferenciadas. La primera parte es un lazo de retroalimentacion positiva creado por una proteina que se une al promotor de su propio gen activando la transcription. La segunda parte es una interaction negativa en la que una molecula represora evita la union de la proteina con el promotor. Un estudio estocastico muestra que el sistema es robusto al ruido. Un estudio determinista muestra que la dinamica del sistema es debida principalmente a dos tipos de biomoleculas: la proteina, y el complejo formado por el represor y esta proteina. La conclusion principal de este estudio es que una simple y usual interaction negativa, tal como una degradation, un secuestro o una inhibition, actuando sobre el lazo de retroalimentacion positiva de un solo gen es una condition suficiente para producir oscilaciones robustas. Un gen es suficiente y el lazo de retroalimentacion positiva no necesita activar a un segundo gen represor, tal y como ocurre en los relojes actuales con dos genes. Esto significa que a nivel genetico un lazo de retroalimentacion negativa no es necesario de forma explicita. Ademas, este modelo no necesita reacciones cooperativas ni la formation de multimeros proteicos, al contrario que en muchos osciladores geneticos. Aplicaciones y futuras lineas de investigacion: En los liltimos anos, se han descubierto muchas moleculas de ARN con capacidad catalitica. El primer modelo de oscilador genetico propuesto en esta tesis usa ribozimas como moleculas repre¬soras. Esto podria proporcionar nuevos principios de diseno en biologia sintetica y una mejor comprension de los relojes celulares regulados por moleculas de ARN. El segundo modelo de oscilador genetico propuesto aqui involucra solo una represion actuando sobre un gen auto-activado y produce oscilaciones robustas. Sorprendente-mente, un segundo gen represor no es necesario al contrario que en los bien conocidos osciladores con dos genes. Este resultado podria ayudar a clarificar los principios de diseno de los relojes celulares naturales y constituir una nueva y eficiente he-rramienta para crear osciladores geneticos sinteticos. Algunas de las futuras lineas de investigation abiertas tras esta tesis son: (1) la validation in vivo e in vitro de ambos modelos, (2) el estudio del potential del segundo modelo como circuito base para la construction de una memoria genetica, (3) el estudio de nuevos osciladores geneticos regulados por ARN no codificante y, por ultimo, (4) el rediseno del se¬gundo modelo de oscilador genetico para su uso como biosensor capaz de detectar genes auto-activados en redes geneticas.

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.

Texture perception determined by doy protein isolate and inulin addition in potato puree: links with mechanical and microstructural features

This study evaluated the effect of adding soy protein isolate (SPI) and long-chain perception, trained and untrained panel inulin (INL) blends with 10 different SPI : INL ratios on the textural, rheological and 17 microstructural properties of freshly made and frozen/thawed potato puree. All the potato puree samples were subjected to a sensory texture pro?le analysis and a 21 trained panel rated the intensity of six descriptors, while an untrained panel did the same on six selected frozen/thawed products. The main SPI : INL ratio effect remained signi?cant for all the descriptors evaluated, when the analysis of variance was applied considering the untrained assessors as random effects. However, only trained panel scores for creaminess corresponded well with untrained assessor. Rheological ?ow index values were linked with variations in perceived consistency, and geometric and surface textural attributes were explained by structural features such as the presence of INL crystallites and SPI coarse strands.

Identification of thaumatin-like protein and aspartyl protease as new major allergens in lettuce (Lactuca sativa)

Scope: Today, about 2–8% of the population of Western countries exhibits some type of food allergy whose impact ranges from localized symptoms confined to the oral mucosa to severe anaphylactic reactions. Consumed worldwide, lettuce is a Compositae family vegetable that can elicit allergic reactions. To date, however, only one lipid transfer protein has been described in allergic reaction to lettuce. The aim of this study was to identify potential new allergens involved in lettuce allergy. Methods and results: Sera from 42 Spanish lettuce-allergic patients were obtained from pa-tients recruited at the outpatient clinic. IgE-binding proteins were detected by SDS-PAGE and immunoblotting. Molecular characterization of IgE-binding bands was performed by MS. Thaumatin was purified using the Agilent 3100 OFFGEL system. The IgE-binding bands recognized in the sera of more than 50% of patients were identified as lipid transfer protein (9 kDa), a thaumatin-like protein (26 kDa), and an aspartyl protease (35 and 45 kDa). ELISA inhibition studies were performed to confirm the IgE reactivity of the purified allergen. Conclusion: Two new major lettuce allergens—a thaumatin-like protein and an aspartyl protease—have been identified and characterized. These allergens may be used to improve both diagnosis and treatment of lettuce-allergic patients.

Thermomechanical relaxation and different water states in cottonseed protein derived bioplastics

Thermomechanical relaxation events and different water states in cottonseed protein bioplastics are presented whilst investigating the effects of aldehyde cross-linking agents. Thermomechanical relaxation of cottonseed protein bioplastics associated with protein denaturation, moisture absorption and broad glass transitions (Tg) were observed from DSC and DMA measurements. It was shown that variation of the aldehyde influences the storage modulus at very low temperature (below Tg). From measurements of the water fusion point, enthalpy, vaporisation, and weight loss, three water states in the water-absorbed bioplastics are suggested; namely strongly-bound-to-polymer, weakly-bound-to-polymer and bulk-like water. The water content and unreacted cross-linking agents are influential factors in controlling formation of the different water states, whilst the selection of different aldehydes was found to be negligible. These results could be valuable for adjusting the thermomechanical relaxations of protein based bioplastics, and tailoring their properties in wet environments.

Specific association of the gene product of PKD2 with the TRPC1 channel

The function(s) of the genes (PKD1 and PKD2) responsible for the majority of cases of autosomal dominant polycystic kidney disease is unknown. While PKD1 encodes a large integral membrane protein containing several structural motifs found in known proteins involved in cell–cell or cell–matrix interactions, PKD2 has homology to PKD1 and the major subunit of the voltage-activated Ca2+ channels. We now describe sequence homology between PKD2 and various members of the mammalian transient receptor potential channel (TRPC) proteins, thought to be activated by G protein-coupled receptor activation and/or depletion of internal Ca2+ stores. We show that PKD2 can directly associate with TRPC1 but not TRPC3 in transfected cells and in vitro. This association is mediated by two distinct domains in PKD2. One domain involves a minimal region of 73 amino acids in the C-terminal cytoplasmic tail of PKD2 shown previously to constitute an interacting domain with PKD1. However, distinct residues within this region mediate specific interactions with TRPC1 or PKD1. The C-terminal domain is sufficient but not necessary for the PKD2–TRPC1 association. A more N-terminal domain located within transmembrane segments S2 and S5, including a putative pore helical region between S5 and S6, is also responsible for the association. Given the ability of the TRPC to form functional homo- and heteromultimeric complexes, these data provide evidence that PKD2 may be functionally related to TRPC proteins and suggest a possible role of PKD2 in modulating Ca2+ entry in response to G protein-coupled receptor activation and/or store depletion.

Optical imaging of functional domains in the cortex of the awake and behaving monkey

As demonstrated by anatomical and physiological studies, the cerebral cortex consists of groups of cortical modules, each comprising populations of neurons with similar functional properties. This functional modularity exists in both sensory and association neocortices. However, the role of such cortical modules in perceptual and cognitive behavior is unknown. To aid in the examination of this issue we have applied the high spatial resolution optical imaging methodology to the study of awake, behaving animals. In this paper, we report the optical imaging of orientation domains and blob structures, approximately 100–200 μm in size, in visual cortex of the awake and behaving monkey. By overcoming the spatial limitations of other existing imaging methods, optical imaging will permit the study of a wide variety of cortical functions at the columnar level, including motor and cognitive functions traditionally studied with positron-emission tomography or functional MRI techniques.

Structure and function of the Bacillus hybrid enzyme GluXyn-1: Native-like jellyroll fold preserved after insertion of autonomous globular domain

The 1,3–1,4-β-glucanase from Bacillus macerans (wtGLU) and the 1,4-β-xylanase from Bacillus subtilis (wtXYN) are both single-domain jellyroll proteins catalyzing similar enzymatic reactions. In the fusion protein GluXyn-1, the two proteins are joined by insertion of the entire XYN domain into a surface loop of cpMAC-57, a circularly permuted variant of wtGLU. GluXyn-1 was generated by protein engineering methods, produced in Escherichia coli and shown to fold spontaneously and have both enzymatic activities at wild-type level. The crystal structure of GluXyn-1 was determined at 2.1 Å resolution and refined to R = 17.7% and R(free) = 22.4%. It shows nearly ideal, native-like folding of both protein domains and a small, but significant hinge bending between the domains. The active sites are independent and accessible explaining the observed enzymatic activity. Because in GluXyn-1 the complete XYN domain is inserted into the compact folding unit of GLU, the wild-type-like activity and tertiary structure of the latter proves that the folding process of GLU does not depend on intramolecular interactions that are short-ranged in the sequence. Insertion fusions of the GluXyn-1 type may prove to be an easy route toward more stable bifunctional proteins in which the two parts are more closely associated than in linear end-to-end protein fusions.

Resistance gene N-mediated de novo synthesis and activation of a tobacco mitogen-activated protein kinase by tobacco mosaic virus infection

Salicylic acid-induced protein kinase (SIPK) and wounding-induced protein kinase (WIPK), two distinct members of the mitogen-activated protein (MAP) kinase family, are activated in tobacco resisting infection by tobacco mosaic virus (TMV). WIPK activation by TMV depends on the disease-resistance gene N because infection of susceptible tobacco not carrying the N gene failed to activate WIPK. Activation of WIPK required not only posttranslational phosphorylation but also a preceding rise in its mRNA and de novo synthesis of WIPK protein. The induction by TMV of WIPK mRNA and protein also occurred systemically. Its activation at the mRNA, protein, and enzyme levels was independent of salicylic acid. The regulation of WIPK at multiple levels by an N gene-mediated signal(s) suggests that this MAP kinase may be an important component upstream of salicylic acid in the signal-transduction pathway(s) leading to local and systemic resistance to TMV.

Protein kinase mutants of human ATR increase sensitivity to UV and ionizing radiation and abrogate cell cycle checkpoint control

In fission yeast, the rad3 gene product plays a critical role in sensing DNA structure defects and activating damage response pathways. A structural homologue of rad3 in humans (ATR) has been identified based on sequence similarity in the protein kinase domain. General information regarding ATR expression, protein kinase activity, and cellular localization is known, but its function in human cells remains undetermined. In the current study, the ATR protein was examined by gel filtration of protein extracts and was found to exist predominantly as part of a large protein complex. A kinase-inactivated form of the ATR gene was prepared by site-directed mutagenesis and was used in transfection experiments to probe the function of this complex. Introduction of this kinase-dead ATR into a normal fibroblast cell line, an ATM-deficient fibroblast line derived from a patient with ataxia–telangiectasia, or a p53 mutant cell line all resulted in significant losses in cell viability. Clones expressing the kinase-dead ATR displayed increased sensitivity to x-rays and UV and a loss of checkpoint control. We conclude that ATR functions as a critical part of a protein complex that mediates responses to ionizing and UV radiation in human cells. These responses include effects on cell viability and cell cycle checkpoint control.

Evidence of insulin-stimulated phosphorylation and activation of the mammalian target of rapamycin mediated by a protein kinase B signaling pathway

The effects of insulin on the mammalian target of rapamycin, mTOR, were investigated in 3T3-L1 adipocytes. mTOR protein kinase activity was measured in immune complex assays with recombinant PHAS-I as substrate. Insulin-stimulated kinase activity was clearly observed when immunoprecipitations were conducted with the mTOR antibody, mTAb2. Insulin also increased by severalfold the 32P content of mTOR that was determined after purifying the protein from 32P-labeled adipocytes with rapamycin⋅FKBP12 agarose beads. Insulin affected neither the amount of mTOR immunoprecipitated nor the amount of mTOR detected by immunoblotting with mTAb2. However, the hormone markedly decreased the reactivity of mTOR with mTAb1, an antibody that activates the mTOR protein kinase. The effects of insulin on increasing mTOR protein kinase activity and on decreasing mTAb1 reactivity were abolished by incubating mTOR with protein phosphatase 1. Interestingly, the epitope for mTAb1 is located near the COOH terminus of mTOR in a 20-amino acid region that includes consensus sites for phosphorylation by protein kinase B (PKB). Experiments were performed in MER-Akt cells to investigate the role of PKB in controlling mTOR. These cells express a PKB-mutant estrogen receptor fusion protein that is activated when the cells are exposed to 4-hydroxytamoxifen. Activating PKB with 4-hydroxytamoxifen mimicked insulin by decreasing mTOR reactivity with mTAb1 and by increasing the PHAS-I kinase activity of mTOR. Our findings support the conclusion that insulin activates mTOR by promoting phosphorylation of the protein via a signaling pathway that contains PKB.

Crystal structure of a multifunctional 2-Cys peroxiredoxin heme-binding protein 23 kDa/proliferation-associated gene product

Heme-binding protein 23 kDa (HBP23), a rat isoform of human proliferation-associated gene product (PAG), is a member of the peroxiredoxin family of peroxidases, having two conserved cysteine residues. Recent biochemical studies have shown that HBP23/PAG is an oxidative stress-induced and proliferation-coupled multifunctional protein that exhibits specific bindings to c-Abl protein tyrosine kinase and heme, as well as a peroxidase activity. A 2.6-Å resolution crystal structure of rat HBP23 in oxidized form revealed an unusual dimer structure in which the active residue Cys-52 forms a disulfide bond with conserved Cys-173 from another subunit by C-terminal tail swapping. The active site is largely hydrophobic with partially exposed Cys-173, suggesting a reduction mechanism of oxidized HBP23 by thioredoxin. Thus, the unusual cysteine disulfide bond is involved in peroxidation catalysis by using thioredoxin as the source of reducing equivalents. The structure also provides a clue to possible interaction surfaces for c-Abl and heme. Several significant structural differences have been found from a 1-Cys peroxiredoxin, ORF6, which lacks the C-terminal conserved cysteine corresponding to Cys-173 of HBP23.

Purification and characterization of a CENP-B homologue protein that binds to the centromeric K-type repeat DNA of Schizosaccharomyces pombe

We have purified and characterized a novel 60-kDa protein that binds to centromeric K-type repeat DNA from Schizosaccharomyces pombe. This protein was initially purified by its ability to bind to the autonomously replicating sequence 3002 DNA. Cloning of the gene encoding this protein revealed that it possesses significant homology to the mammalian centromere DNA-binding protein CENP-B and S. pombe Abp1, and this gene was designated as cbh+ (CENP-B homologue). Cbh protein specifically interacts in vitro with the K-type repeat DNA, which is essential for centromere function. The Cbh-binding consensus sequence was determined by DNase I footprinting assays as PyPuATATPyPuTA, featuring an inverted repeat of the first four nucleotides. Based on its binding activity to centromeric DNA and homology to centromere proteins, we suggest that this protein may be a functional homologue of the mammalian CENP-B in S. pombe.

Primary structure and tissue distribution of two novel proline-rich γ-carboxyglutamic acid proteins

Two human cDNAs that encode novel vitamin K-dependent proteins have been cloned and sequenced. The predicted amino acid sequences suggest that both are single-pass transmembrane proteins with amino-terminal γ-carboxyglutamic acid-containing domains preceded by the typical propeptide sequences required for posttranslational γ-carboxylation of glutamic acid residues. The polypeptides, with deduced molecular masses of 23 and 17 kDa, are proline-rich within their putative cytoplasmic domains and contain several copies of the sequences PPXY and PXXP, motifs found in a variety of signaling and cytoskeletal proteins. Accordingly, these two proteins have been called proline-rich Gla proteins (PRGP1 and PRGP2). Unlike the γ-carboxyglutamic acid domain-containing proteins of the blood coagulation cascade, the two PRGPs are expressed in a variety of extrahepatic tissues, with PRGP1 and PRGP2 most abundantly expressed in the spinal cord and thyroid, respectively, among those tissues tested. Thus, these observations suggest a novel physiological role for these two new members of the vitamin K-dependent family of proteins.

Structural mimicry of a native protein by a minimized binding domain

The affinity between molecules depends both on the nature and presentation of the contacts. Here, we observe coupling of functional and structural elements when a protein binding domain is evolved to a smaller functional mimic. Previously, a 38-residue form of the 59-residue B-domain of protein A, termed Z38, was selected by phage display. Z38 contains 13 mutations and binds IgG only 10-fold weaker than the native B-domain. We present the solution structure of Z38 and show that it adopts a tertiary structure remarkably similar to that observed for the first two helices of B-domain in the B-domain/Fc complex [Deisenhofer, J. (1981) Biochemistry 20, 2361–2370], although it is significantly less stable. Based on this structure, we have improved on Z38 by designing a 34-residue disulfide-bonded variant (Z34C) that has dramatically enhanced stability and binds IgG with 9-fold higher affinity. The improved stability of Z34C led to NMR spectra with much greater chemical shift dispersion, resulting in a more precisely determined structure. Z34C, like Z38, has a structure virtually identical to the equivalent region from native protein A domains. The well-defined hydrophobic core of Z34C reveals key structural features that have evolved in this small, functional domain. Thus, the stabilized two-helix peptide, about half the size and having one-third of the remaining residues altered, accurately mimics both the structure and function of the native domain.