CCTeta: A novel soluble guanylyl cyclase-interacting protein

Nitric oxide (NO) transduces most of its biological effects through activation of the heterodimeric enzyme, soluble guanylyl cyclase (sGC). Activation of sGC results in the production of 3′,5 ′-cyclic guanosine monophosphate (cGMP) from 5′ -guanosine triphosphate (GTP). In this thesis, we demonstrate a novel protein interaction between CCT (chaperonin containing t-complex polypeptide) subunit η and the α1β1 isoform of sGC. Using the yeast-two-hybrid system, CCTη was found to interact with the N-terminal portion of β1 subunit of sGC. This interaction was then confirmed in vitro with a co-immunoprecipitation from mouse brain. The interaction between these two proteins was further supported by a co-localization of the proteins within rat brain. Using the yeast-two-hybrid system, CCTη was found to bind to the N-terminal portion of sGC. In vitro assays with purified CCTη and Sf9 lysate expressing sGC resulted in a 33% inhibition of sodium nitroprusside (SNP)-stimulated sGC activity. The same assays were then performed using BAY41-2272, an NO-independent allosteric sGC activator, and CCTη had no effect on this activity. Furthermore, CCTη had no effect on the activity of αβCys105 sGC a constitutively active mutant that lacks a heme group. Of note is the fact that the full-length CCTη-expressing bacterial lysate inhibited the activity of sGC-expressing Sf9 lysate by 48% compared with GST alone. This indicates that the amino terminal 94 amino acids of CCTη are important to the inhibition of sGC activity. Lastly, a 45% inhibition of sGC activity by CCTη was seen in vivo in BE2 cells stably transfected with CCTη and treated with SNP. The fact that the inhibition of sGC was more pronounced with bacterial lysate expressing CCTη versus the purified CCTη implies that some factor in the bacterial lysate enhances the inhibitory effect of CCTη. Because the level of inhibition seen in bacterial lysate and in vivo experiments is similar, might imply that the factor that aids in CCTη effect on sGC is conserved. Together, these data suggest that CCTη is a novel type of sGC inhibitor that inhibits sGC by modifying the binding of NO to the heme group or the subsequent conformational changes induced by NO binding. ^

Elucidation of the role of 53BP1 in DNA double strand break (DSB) repair and tumor suppression

The protein p53 binding protein one (53BP1) was discovered in a yeast two-hybrid screen that used the DNA binding domain of p53 as bait. Cloning of full-length 53BP1 showed that this protein contains several protein domains which help make up the protein, which include two tandem BRCT domains and a amino-terminal serine/glutamine cluster domain (SCD). These are two protein domains are often seen in factors that are involved in the cellular response to DNA damage and control of cell cycle checkpoints and we hypothesize that 53BP1 is involved in the cellular response to DNA damage. In support of this hypothesis we observe that 53BP1 is phosphorylated and undergoes a dramatic nuclear re-localization in response to DNA damaging agents. 53BP1 also interacts with several factors that are important in the cellular response to DNA damage, such as the BRCA1 tumor suppressor, ATM and Rad3 related (ATR), and the phosphorylated version of the histone variant H2AX. Mice deficient in 53BP1 display increased sensitivity ionizing radiation (IR), a DNA damaging agent that introduces DNA double strand breaks (DSBs). In addition, 53BP1-deficient mice do not properly undergo the process of class switch recombination (CSR). We also observe that when a defect in 53BP1 is combined with a defect in p53; the resulting mice have an increased rate of formation of spontaneous tumors, notably the formation of B and T lineage lymphomas. The T lineage tumors arise by two distinct mechanisms: one driven by defects in cell cycle regulation and a second driven by defects in the ability to repair DNA DSBs. The B lineage tumors arise by the inability to repair DNA damage and over-expression of the oncogene c-myc. ^ With these observations, we conclude that not only does 53BP1 function in the cellular response to DNA damage, but it also works in concert with p53 to suppress tumor formation. ^

A mechanism of insulin-like growth factor binding protein 2-induced cell mobility in glioma

Insulin-like growth factor binding protein 2 (IGFBP2) is a protein known to be overexpressed in a majority of glioblastoma multiforme (GBM) tumors. While it is known the IGFBP2 is involved in promoting GBM tumor cell invasion, no mechanism exists for how the protein is involved in signal transduction pathways leading to enhanced cell invasion. ^ We follow up on preliminary microarray data on IGFBP2-overexpressing GBM cells and protein sequence analysis of IGFBP2 in generating the hypothesis that IGFBP2 interacts with integnn α5 in regulating cell mobility. Microarray data showing upregulation of integrin α5 by IGFBP2 is validated and evidence of protein-protein interaction between IGFBP2 and integrin α5 is found. The exact binding domain on IGFBP2 responsible for its interaction with integrin α5 is also determined, confirming our initial findings and reaffirming that the IGFBP2/integrin α5 interaction is specific. Disruption of this interaction resulted in attenuation of IGFBP2-enhanced cell mobility. Further, we found that cell mobility is only enhanced when IGFBP2 and integrin α5 are both overexpressed and able to interact with each other. ^ We also determined fibronectin to be a critical player in the activation of the IGFBP2/integrin α5 pathway. The activation of this pathway appears to be progressive and initiates once GBM cells have sufficiently established anchorage. ^

From protein microarray to biology: The discovery of epigenetic regulatory molecules

Chromatin, composed of repeating nucleosome units, is the genetic polymer of life. To aid in DNA compaction and organized storage, the double helix wraps around a core complex of histone proteins to form the nucleosome, and is therefore no longer freely accessible to cellular proteins for the processes of transcription, replication and DNA repair. Over the course of evolution, DNA-based applications have developed routes to access DNA bound up in chromatin, and further, have actually utilized the chromatin structure to create another level of complexity and information storage. The histone molecules that DNA surrounds have free-floating tails that extend out of the nucleosome. These tails are post-translationally modified to create docking sites for the proteins involved in transcription, replication and repair, thus providing one prominent way that specific genomic sequences are accessed and manipulated. Adding another degree of information storage, histone tail-modifications paint the genome in precise manners to influence a state of transcriptional activity or repression, to generate euchromatin, containing gene-dense regions, or heterochromatin, containing repeat sequences and low-density gene regions. The work presented here is the study of histone tail modifications, how they are written and how they are read, divided into two projects. Both begin with protein microarray experiments where we discover the protein domains that can bind modified histone tails, and how multiple tail modifications can influence this binding. Project one then looks deeper into the enzymes that lay down the tail modifications. Specifically, we studied histone-tail arginine methylation by PRMT6. We found that methylation of a specific histone residue by PRMT6, arginine 2 of H3, can antagonize the binding of protein domains to the H3 tail and therefore affect transcription of genes regulated by the H3-tail binding proteins. Project two focuses on a protein we identified to bind modified histone tails, PHF20, and was an endeavor to discover the biological role of this protein. Thus, in total, we are looking at a complete process: (1) histone tail modification by an enzyme (here, PRMT6), (2) how this and other modifications are bound by conserved protein domains, and (3) by using PHF20 as an example, the functional outcome of binding through investigating the biological role of a chromatin reader. ^

Identification of a conserved cluster in the RH domain of GRK critical for activation by GPCRs

One of the most critical aspects of G Protein Coupled Receptors (GPCRs) regulation is their rapid and acute desensitization following agonist stimulation. Phosphorylation of these receptors by GPCR kinases (GRK) is a major mechanism of desensitization. Considerable evidence from studies of rhodopsin kinase and GRK2 suggests there is an allosteric docking site for the receptor distinct from the GRK catalytic site. While the agonist-activated GPCR appears crucial for GRK activation, the molecular details of this interaction remain unclear. Recent studies suggested an important role for the N- and C-termini and domains in the small lobe of the kinase domain in allosteric activation; however, neither the mechanism of action of that site nor the RH domain contributions have been elucidated. To search for the allosteric site, we first indentified evolutionarily conserved sites within the RH and kinase domains presumably deterministic of protein function employing evolutionary trace (ET) methodology and crystal structures of GRK6. Focusing on a conserved cluster centered on helices 3, 9, and 10 in the RH domain, key residues of GRK5 and 6 were targeted for mutagenesis and functional assays. We found that a number of double mutations within helices 3, 9, and 10 and the N-terminus markedly reduced (50–90%) the constitutive phosphorylation of the β-2 Adrenergic Receptor (β2AR) in intact cells and phosphorylation of light-activated rhodopsin (Rho*) in vitro as compared to wild type (WT) GRK5 or 6. Based on these results, we designed peptide mimetics of GRK5 helix 9 both computationally and through chemical modifications with the goal of both confirming the importance of helix 9 and developing a useful inhibitor to disrupt the GPCR-GRK interaction. Several peptides were found to block Rho* phosphorylation by GRK5 including the native helix 9 sequence, Peptide Builder designed-peptide preserving only the key ET residues, and chemically locked helices. Most peptidomimetics showed inhibition of GRK5 activity greater than 80 % with an IC50 of ∼ 30 µM. Alanine scanning of helix 9 has further revealed both essential and non-essential residues for inhibition. Importantly, substitution of Arg 169 by an alanine in the native helix 9-based peptide gave an almost complete inhibition at 30 µM with an IC50 of ∼ 10 µM. In summary we report a previously unrecognized crucial role for the RH domain of GRK5 and 6, and the subsequent identification of a lead peptide inhibitor of protein-protein interaction with potential for specific blockade of GPCR desensitization. ^

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Analysis of VirB4, a membrane-associated ATPase subunit essential for T -complex transport in Agrobacterium tumefaciens

Agrobacterium tumefaciens is a plant pathogen with the unique ability to export oncogenic DNA-protein complexes (T-complexes) to susceptible plant cells and cause crown gall tumors. Delivery of the T-complexes across the bacterial membranes requires eleven VirB proteins and VirD4, which are postulated to form a transmembrane transporter. This thesis examines the subcellular localization and oligomeric structure of the 87-kDa VirB4 protein, which is one of three essential ATPases proposed to energize T-complex transport and/or assembly. Results of subcellular localization studies showed that VirB4 is tightly associated with the cytoplasmic membrane, suggesting that it is a membrane-spanning protein. The membrane topology of VirB4 was determined by using a nested deletion strategy to generate random fusions between virB4 and the periplasmically-active alkaline phosphatase, $\sp\prime phoA$. Analysis of PhoA and complementary $\beta$-galactosidase reporter fusions identified two putative periplasmically-exposed regions in VirB4. A periplasmic exposure of one of these regions was further confirmed by protease susceptibility assays using A. tumefaciens spheroplasts. To gain insight into the structure of the transporter, the topological configurations of other VirB proteins were also examined. Results from hydropathy analyses, subcellular localization, protease susceptibility, and PhoA reporter fusion studies support a model that all of the VirB proteins localize at one or both of the bacterial membranes. Immunoprecipitation and Co$\sp{2+}$ affinity chromatography studies demonstrated that native VirB4 (87-kDa) and a functional N-terminally tagged HIS-VirB4 derivative (89-kDa) interact and that the interaction is independent of other VirB proteins. A $\lambda$ cI repressor fusion assay supplied further evidence for VirB4 dimer formation. A VirB4 dimerization domain was localized to the N-terminal third of the protein, as judged by: (i) transdominance of an allele that codes for this region of VirB4; (ii) co-retention of a His-tagged N-terminal truncation derivative and native VirB4 on Co$\sp{2+}$ affinity columns; and (iii) dimer formation of the N-terminal third of VirB4 fused to the cI repressor protein. Taken together, these findings are consistent with a model that VirB4 is topologically configured as an integral cytoplasmic membrane protein with two periplasmic domains and that VirB4 assembles as homodimers via an N-terminal dimerization domain. Dimer formation is postulated to be essential for stabilization of VirB4 monomers during T-complex transporter assembly. ^

Regulations of SOX9 activity during chondrogenesis

Sox9 is a transcription factor required for chondrocyte differentiation and cartilage formation. In an effort to identify SOX9 interacting protein(s), we screened a chondrocyte cDNA library with a modified yeast two-hybrid method, Son of Sevenless (SOS) recruitment system (SRS). The catalytic subunit of cyclic AMP-dependent protein kinase A (PKA-Cα) and a new long form of c-Maf transcription factor (Lc-Maf) were found to interact specifically with SOX9. We showed here that two PKA phosphorylation consensus sites of SOX9 could be phosphorylated by PKA in vitro as well as in vivo. PKA phosphorylation of SOX9 increases its DNA binding and transcriptional activities on a Col2a1 chondrocyte-specific enhancer. Mutations of these two PKA phosphorylation sites markedly decreased the activation of SOX9 by PKA. ^ To test whether parathyroid hormone-related peptide (PTHrP) signaling results in SOX9 phosphorylation, we generated a phosphospecific antibody that specifically recognizes SOX9 that is phosphorylated at serine 181 (S 181) one of the two consensus PKA phosphorylation sites. Addition of PTHrP to COS7 cells cotransfected with SOX9 and PTH/PTHrP receptor strongly increased phosphorylation of SOX9 at S181; this phosphorylation was blocked by a PKA-specific inhibitor. In similar experiments we showed that PTHrP increased the activity of a SOX9-dependent Col2a1 enhancer. This increase in activity was abolished when a SOX9 mutant was used containing serine-to-alanine substitution in the two consensus PKA phosphorylation sites of SOX9. Using our phosphospecific SOX9 antibody we showed by immunohistochemistry of mouse embryos that Sox9 phosphorylated at S181 was localized almost exclusively in the pre-hypertrophic zone of the growth plate, an area corresponding to the major site of expression of PTH/PTHrP receptor. In contrast, no phosphorylation of Sox9 at S181 was detected in growth plates of PTH/PTHrP receptor null mutant mice. Sox9, regardless of phosphorylation state, was present in all chondrocytes of both genotypes except in hypertrophic chondrocytes. Thus, Sox9 is a target of PTHrP signaling and the PTHrP-dependent phosphorylation of SOX9 enhances its transcriptional activity. ^ In order to investigate the in vivo function of Sox9 phosphorylation by PKA, we are generating a mouse model of mutant Sox9 harboring point mutations in two PKA phosphorylation sites. Preliminary results indicated that heterozygous mice containing half amount of mutant Sox9 that can not be phosphorylated by PKA have normal skeletal phenotype and homozygous mice are being generated. ^ Lc-Maf encodes an extra ten amino acids at the carboxyl terminus of c-Maf and contains a completely different 3′ untranslated region. The interaction between SOX9 and Lc-Maf was further confirmed by co-immunoprecipitation and GST-pull down assays, which mapped the interacting domains of SOX9 to HMG DNA binding domain and that of Lc-Maf to basic leusine zipper motif. In situ hybridizations showed that RNA of Lc-Maf coexpressed with those of Sox9 and Col2a1 in areas of mesenchymal condensation during the early stages of mouse embryo development. A DNA binding site of Lc-Maf was identified at the 5′ part of a 48-bp Col2a1 enhancer element near the HMG binding site of SOX9. Lc-Maf and SOX9 synergistically activated a luciferase reporter plasmid containing a Col2al enhancer and increased the transcription of endogenous Col2a1 gene. In summary, Lc-Maf is the first identified SOX9-interating protein during chondrogenesis and may be an important activator of Col2a1 gene. ^

Signaling pathways in the transcriptional and post-translational regulation of the human glutathione S -transferase P1 gene

The human glutathione S-transferase P1 (GSTP1) protein is an endogenous inhibitor of c-jun N-terminal kinases (JNKs) and an important phase II detoxification enzyme. ^ Recent identification of a cAMP response element (CRE) in the 5 ′-region of the human GSTP1 gene and several putative phosphorylation sites for the Ser/Thr protein kinases, including, cAMP-dependent protein kinases (PKAs), protein kinases C (PKCs), and JNKs in the GSTP1 protein raised the possibility that signaling pathways may play an important role in the transcriptional and post-translational regulation of GSTP1 gene. This study examined (a) whether the signaling pathway mediated by CAMP, via the GSTP1 CRE, is involved in the transcriptional regulation of the GSTP1 gene, (b) whether signaling pathways mediated by the Ser/Thr protein kinases (PKAs, PKCs, and JNKs) induce post-translational modification, viz. phosphorylation of the GSTP1 protein, and (c) whether such phosphorylation of the GSTP1 protein alters its functions in metabolism and in JNK signaling. ^ The first major finding in this study is the establishment of the human GSTP1 gene as a novel CAMP responsive gene in which transcription is activated via an interaction between PKA activated CRE binding protein-1 (CREB-1) and the CRE in the 5′-regulatory region. ^ The second major finding in this study is the observation that the GSTP1 protein undergoes phosphorylation and functionally activated by second messenger-activated protein kinases, PKA and PKC, in tumor cells with activated signaling pathways. Following phosphorylation by PKA or PKC, the catalytic activity of the GSTP1 protein was significantly enhanced, as indicated by a decrease in its Km (2- to 3.6-fold) and an increase in Kcat/ Km (1.6- to 2.5-fold) for glutathione. Given the frequent over-expression of GSTP1 and the aberrant PKA/PKC signaling cascade observed in tumors, these findings suggest that phosphorylation of GSTP1 may contribute to the malignant progression and drug-resistant phenotype of these tumors. ^ The third major finding in this study is that the GSTP1 protein, an inhibitor of JNKs, undergoes significant phosphorylation in tumor cells with activated JNK signaling pathway and in those under oxidative stress. Following phosphorylation by JNK, the ability of GSTP1 to inhibit JNK downstream function, i.e. c-jun phosphorylation, was significantly enhanced, suggesting a feedback mechanism of regulation of JNK-mediated cellular signaling. (Abstract shortened by UMI.) ^

Increased glutamine in leaves of poplar transgenic with pine GS1a caused greater anthranilate synthetase a-subunit (ASA1) transcript and protein abundances: an auxin-related mechanism for enhanced growth in GS transgenics?

The initial reaction in the pathway leading to the production of indole-3-acetic acid (IAA) in plants is the reaction between chorismate and glutamine to produce anthranilate, catalysed by the enzyme anthranilate synthase (ASA; EC 4.1.3.27). Compared with non-transgenic controls, leaves of transgenic poplar with ectopic expression of the pine cytosolic glutamine synthetase (GS1a; EC 6.3.1.2) produced significantly greater glutamine and significantly enhanced ASA a-subunit (ASA1) transcript and protein (approximately 130% and 120% higher than in the untransformed controls, respectively). Similarly, tobacco leaves fed with 30 mM glutamine and 2 mM chorismate showed enhanced ASA1 transcript and protein (175% and 90% higher than controls, respectively). Furthermore, free IAA was significantly elevated both in leaves of GS1a transgenic poplar and in tobacco leaves fed with 30 mM glutamine and 2 mM chorismate. These results indicated that enhanced cellular glutamine may account for the enhanced growth in GS transgenic poplars through the regulation of auxin biosynthesis

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.