Storage and retrieval of individual genomes

A repetitive sequence collection is one where portions of a base sequence of length n are repeated many times with small variations, forming a collection of total length N. Examples of such collections are version control data and genome sequences of individuals, where the differences can be expressed by lists of basic edit operations. Flexible and efficient data analysis on a such typically huge collection is plausible using suffix trees. However, suffix tree occupies O(N log N) bits, which very soon inhibits in-memory analyses. Recent advances in full-text self-indexing reduce the space of suffix tree to O(N log σ) bits, where σ is the alphabet size. In practice, the space reduction is more than 10-fold, for example on suffix tree of Human Genome. However, this reduction factor remains constant when more sequences are added to the collection. We develop a new family of self-indexes suited for the repetitive sequence collection setting. Their expected space requirement depends only on the length n of the base sequence and the number s of variations in its repeated copies. That is, the space reduction factor is no longer constant, but depends on N / n. We believe the structures developed in this work will provide a fundamental basis for storage and retrieval of individual genomes as they become available due to rapid progress in the sequencing technologies.

Boxicity and cubicity of asteroidal triple free graphs

The boxicity of a graph G, denoted box(G), is the least integer d such that G is the intersection graph of a family of d-dimensional (axis-parallel) boxes. The cubicity, denoted cub(G), is the least dsuch that G is the intersection graph of a family of d-dimensional unit cubes. An independent set of three vertices is an asteroidal triple if any two are joined by a path avoiding the neighbourhood of the third. A graph is asteroidal triple free (AT-free) if it has no asteroidal triple. The claw number psi(G) is the number of edges in the largest star that is an induced subgraph of G. For an AT-free graph G with chromatic number chi(G) and claw number psi(G), we show that box(G) <= chi(C) and that this bound is sharp. We also show that cub(G) <= box(G)([log(2) psi(G)] + 2) <= chi(G)([log(2) psi(G)] + 2). If G is an AT-free graph having girth at least 5, then box(G) <= 2, and therefore cub(G) <= 2 [log(2) psi(G)] + 4. (c) 2010 Elsevier B.V. All rights reserved.

Energy landscape view of phase transitions and slow dynamics in thermotropic liquid crystals

Thermotropic liquid crystals are known to display rich phase behavior on temperature variation. Although the nematic phase is orientationally ordered but translationally disordered, a smectic phase is characterized by the appearance of a partial translational order in addition to a further increase in orientational order. In an attempt to understand the interplay between orientational and translational order in the mesophases that thermotropic liquid crystals typically exhibit upon cooling from the high-temperature isotropic phase, we investigate the potential energy landscapes of a family of model liquid crystalline systems. The configurations of the system corresponding to the local potential energy minima, known as the inherent structures, are determined from computer simulations across the mesophases. We find that the depth of the potential energy minima explored by the system along an isochor grows through the nematic phase as temperature drops in contrast to its insensitivity to temperature in the isotropic and smectic phases. The onset of the growth of the orientational order in the parent phase is found to induce a translational order, resulting in a smectic-like layer in the underlying inherent structures; the inherent structures, surprisingly, never seem to sustain orientational order alone if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. The Arrhenius temperature dependence of the orientational relaxation time breaks down near the isotropic-nematic transition. We find that this breakdown occurs at a temperature below which the system explores increasingly deeper potential energy minima.

Biochemical characterization of the intracellular domain of the human guanylyl cyclase C receptor provides evidence for a catalytically active homotrimer

Guanylyl cyclase C (GCC) is the receptor for the family of guanylin peptides and bacterial heat-stable enterotoxins (ST). The receptor is composed of an extracellular, ligand-binding domain and an intracellular domain with a region of homology to protein kinases and a guanylyl cyclase catalytic domain. We have expressed the entire intracellular domain of GCC in insect cells and purified the recombinant protein, GCC-IDbac, to study its catalytic activity and regulation. Kinetic properties of the purified protein were similar to that of full-length GCC, and high activity was observed when MnGTP was used as the substrate. Nonionic detergents, which stimulate the guanylyl cyclase activity of membrane-associated GCC, did not appreciably increase the activity of GCC-IDbac, indicating that activation of the receptor by Lubrol involved conformational changes that required the transmembrane and/or the extracellular domain. The guanylyl cyclase activity of GCC-IDbac was inhibited by Zn2+, at concentrations shown to inhibit adenylyl cyclase, suggesting a structural homology between the two enzymes. Covalent crosslinking of GCC-IDbac indicated that the protein could associate as a dimer, but a large fraction was present as a trimer. Gel filtration analysis also showed that the major fraction of the protein eluted at a molecular size of a trimer, suggesting that the dimer detected by cross-linking represented subtle differences in the juxtaposition of the individual polypeptide chains. We therefore provide evidence that the trimeric state of GCC is catalytically active, and sequences required to generate the trimer are present in the intracellular domain of GCC.

Structural and magnetic properties of Sr2Fe1+xMo1-xO6 (-1 <= x <= 0.25)

We have synthesized the solid solution Sr2Fe1+xMo1-xO6 with -1 <= x <= 0.25, the composition x=0 corresponding to the well-known double-perovskite system Sr2FeMoO6. We report structural and magnetic properties of the above system, exhibiting systematic variations across the series. These results restrict the range of models that can explain magnetism in this family of compounds, providing an understanding of the magnetic structure.

A note on the equivalence of shock manifold equations

The shock manifold equation is a first order nonlinear partial differential equation, which describes the kinematics of a shockfront in an ideal gas with constant specific heats. However, it was found that there was more than one of these shock manifold equations, and the shock surface could be embedded in a one parameter family of surfaces, obtained as a solution of any of these shock manifold equations. Associated with each shock manifold equation is a set of characteristic curves called lsquoshock raysrsquo. This paper investigates the nature of various associated shock ray equations.

High-temperature superconductivity in bismuth-copper oxides of the type BimMnCupOx (M=Mg, Ca, Sr, Ba, Bi)

Several oxides of the Bi m M n Cu p O x family (m=2, 3;n=2, 3, 4;p=1, 2, 3 and M=alkaline earth or Bi), possessing structures similar to the Aurivillius family of oxides, show highT c superconductivity.

High-yield bacterial expression and structural characterization of recombinant human insulin-like growth factor binding protein-2

The diverse biological activities of the insulin-like growth factors (IGF-1 and IGF-2) are mediated by the IGF-1 receptor (IGF-1R). These actions are modulated by a family of six IGF-binding proteins (ICFBP-1-6; 22-31 kDa) that via high affinity binding to the IGFs (K-D similar to 300-700 pM) both protect the IGFs in the circulation and attenuate IGF action by blocking their receptor access In recent years, IGFBPs have been implicated in a variety of cancers However, the structural basis of their interaction with IGFs and/or other proteins is not completely understood A critical challenge in the structural characterization of full-length IGFBPs has been the difficulty in expressing these proteins at levels suitable for NMR/X-ray crystallography analysis Here we describe the high-yield expression of full-length recombinant human IGFBP-2 (rhIGFBP-2) in Eschericha coli Using a single step purification protocol, rhIGFBP-2 was obtained with >95% purity and structurally characterized using NMR spectroscopy. The protein was found to exist as a monomer at the high concentrations required for structural studies and to exist in a single conformation exhibiting a unique intra-molecular disulfide-bonding pattern The protein retained full biologic activity. This study represents the first high-yield expression of wild-type recombinant human IGFBP-2 in E coli and first structural characterization of a full-length IGFBP (C) 2010 Elsevier Inc. All rights reserved

Discovery of oxidative enzymes for food engineering : Tyrosinase and sulfhydryl oxidase

Enzymes offer many advantages in industrial processes, such as high specificity, mild treatment conditions and low energy requirements. Therefore, the industry has exploited them in many sectors including food processing. Enzymes can modify food properties by acting on small molecules or on polymers such as carbohydrates or proteins. Crosslinking enzymes such as tyrosinases and sulfhydryl oxidases catalyse the formation of novel covalent bonds between specific residues in proteins and/or peptides, thus forming or modifying the protein network of food. In this study, novel secreted fungal proteins with sequence features typical of tyrosinases and sulfhydryl oxidases were iden-tified through a genome mining study. Representatives of both of these enzyme families were selected for heterologous produc-tion in the filamentous fungus Trichoderma reesei and biochemical characterisation. Firstly, a novel family of putative tyrosinases carrying a shorter sequence than the previously characterised tyrosinases was discovered. These proteins lacked the whole linker and C-terminal domain that possibly play a role in cofactor incorporation, folding or protein activity. One of these proteins, AoCO4 from Aspergillus oryzae, was produced in T. reesei with a production level of about 1.5 g/l. The enzyme AoCO4 was correctly folded and bound the copper cofactors with a type-3 copper centre. However, the enzyme had only a low level of activity with the phenolic substrates tested. Highest activity was obtained with 4-tert-butylcatechol. Since tyrosine was not a substrate for AoCO4, the enzyme was classified as catechol oxidase. Secondly, the genome analysis for secreted proteins with sequence features typical of flavin-dependent sulfhydryl oxidases pinpointed two previously uncharacterised proteins AoSOX1 and AoSOX2 from A. oryzae. These two novel sulfhydryl oxidases were produced in T. reesei with production levels of 70 and 180 mg/l, respectively, in shake flask cultivations. AoSOX1 and AoSOX2 were FAD-dependent enzymes with a dimeric tertiary structure and they both showed activity on small sulfhydryl compounds such as glutathione and dithiothreitol, and were drastically inhibited by zinc sulphate. AoSOX2 showed good stabil-ity to thermal and chemical denaturation, being superior to AoSOX1 in this respect. Thirdly, the suitability of AoSOX1 as a possible baking improver was elucidated. The effect of AoSOX1, alone and in combi-nation with the widely used improver ascorbic acid was tested on yeasted wheat dough, both fresh and frozen, and on fresh water-flour dough. In all cases, AoSOX1 had no effect on the fermentation properties of fresh yeasted dough. AoSOX1 nega-tively affected the fermentation properties of frozen doughs and accelerated the damaging effects of the frozen storage, i.e. giving a softer dough with poorer gas retention abilities than the control. In combination with ascorbic acid, AoSOX1 gave harder doughs. In accordance, rheological studies in yeast-free dough showed that the presence of only AoSOX1 resulted in weaker and more extensible dough whereas a dough with opposite properties was obtained if ascorbic acid was also used. Doughs containing ascorbic acid and increasing amounts of AoSOX1 were harder in a dose-dependent manner. Sulfhydryl oxidase AoSOX1 had an enhancing effect on the dough hardening mechanism of ascorbic acid. This was ascribed mainly to the produc-tion of hydrogen peroxide in the SOX reaction which is able to convert the ascorbic acid to the actual improver dehydroascorbic acid. In addition, AoSOX1 could possibly oxidise the free glutathione in the dough and thus prevent the loss of dough strength caused by the spontaneous reduction of the disulfide bonds constituting the dough protein network. Sulfhydryl oxidase AoSOX1 is therefore able to enhance the action of ascorbic acid in wheat dough and could potentially be applied in wheat dough baking.

The function of NR3B and NR4A orphan nuclear receptors in osteoblasts

Nuclear receptors (NRs) comprise a large family of proteins that mediate the effects of small lipophilic molecules such as steroid hormones. In addition, there are a group of NRs which lack identified natural ligands and are referred as orphan NRs. In this thesis, the function of two such orphan NR families, the NR3B (ERRα, ERRβ and ERRγ) and the NR4A family (NGFI-B, Nurr1 and Nor1), was studied. NR3B and NR4A receptors regulate many biological processes such as energy metabolism and carcinogenesis. In addition, NR3B and NR4A receptors are expressed in bone. Therefore, the signaling and function of NR3B and NR4A orphan nuclear receptors was studied specifically in osteoblasts. NR4A receptors were found to be regulated by NR3B receptors and the Wnt/β-catenin signaling pathway as ERRα, ERRγ and β-catenin repressed the transcriptional activity of NR4A receptors in U2-OS cells. NGFI-B was found to repress the transcriptional activity of ERRγ in HeLa cells. The phytoestrogen equol was identified as a new agonist for ERRγ and ERRβ in PC-3, U2-OS, and SaOS-2 cells. Equol increased the transcriptional activity of ERRγ by increasing ERRγ co-activator binding and by inducing a conformational change in the ligand binding pocket of ERRγ. The growth inhibitory effect of equol on PC-3 prostate cancer cells was decreased by blocking ERRγ expression by siRNA. Therefore, ERRγ could mediate some of the beneficial health effects of equol. The Wnt/β-catenin signaling pathway is important for the differentiation and function of osteoblasts. NR3B and NR4A receptors were found to repress the transcriptional activity mediated by β-catenin in U2-OS cells. The mesenchymal stem cells (MSCs) isolated from ERRα knockout (KO) mice showed diminished proliferation and osteoblastic differentiation compared to the wild-type cells. The overexpression of ERRα in osteoblastic MC3T3-E1 cell line increased their mineralization. Bone sialoprotein (BSP) was shown to be a direct target gene for ERRα and ERRγ as the BSP promoter was activated by ERRα or ERRγ and PGC-1α in HeLa cells. The adipogenic differentiation of ERRα KO MSCs was also decreased and they expressed less adipogenic marker genes. In conclusion, the studies described in this thesis demonstrated that the transcriptional activity of NR3B and NR4A receptors can be regulated by other orphan NRs and signaling pathways in osteoblasts. NR3B receptors can also be regulated by ligands and a new agonist, equol, was identified for ERRβ and ERRγ. New roles for NR3B and NR4A were also identified as they were shown to converge with the Wnt signaling pathway in osteoblasts, ERRγ was shown to mediate the growth inhibitory effect of equol in prostate cancer cells, and ERRα was shown to regulate positively MSC proliferation, osteoblastic differentiation and adipogenesis.

Steady-State Stability of Current Mode Active-Clamp ZVS DC-DC Converters

Active-clamp dc-dc converters are pulsewidth-modulated converters having two switches featuring zero-voltage switching at frequencies beyond 100 kHz. Generalized equivalent circuits valid for steady-state and dynamic performance have been proposed for the family of active-clamp converters. The active-clamp converter is analyzed for its dynamic behavior under current control in this paper. The steady-state stability analysis is presented. On account of the lossless damping inherent in the active-clamp converters, it appears that the stability region in the current-controlled active-clamp converters get extended for duty ratios, a little greater than 0.5, unlike in conventional hard-switched converters. The conventional graphical approach fails to assess the stability of current-controlled active-clamp converters due to the coupling between the filter inductor current and resonant inductor current. An analysis that takes into account the presence of the resonant elements is presented to establish the condition for stability. This method correctly predicts the stability of the current-controlled active-clamp converters. A simple expression for the maximum duty cycle for subharmonic free operation is obtained. The results are verified experimentally.

Classification of Protein Kinases on the Basis of Both Kinase and Non-Kinase Regions

Background: Protein phosphorylation is a generic way to regulate signal transduction pathways in all kingdoms of life. In many organisms, it is achieved by the large family of Ser/Thr/Tyr protein kinases which are traditionally classified into groups and subfamilies on the basis of the amino acid sequence of their catalytic domains. Many protein kinases are multidomain in nature but the diversity of the accessory domains and their organization are usually not taken into account while classifying kinases into groups or subfamilies. Methodology: Here, we present an approach which considers amino acid sequences of complete gene products, in order to suggest refinements in sets of pre-classified sequences. The strategy is based on alignment-free similarity scores and iterative Area Under the Curve (AUC) computation. Similarity scores are computed by detecting common patterns between two sequences and scoring them using a substitution matrix, with a consistent normalization scheme. This allows us to handle full-length sequences, and implicitly takes into account domain diversity and domain shuffling. We quantitatively validate our approach on a subset of 212 human protein kinases. We then employ it on the complete repertoire of human protein kinases and suggest few qualitative refinements in the subfamily assignment stored in the KinG database, which is based on catalytic domains only. Based on our new measure, we delineate 37 cases of potential hybrid kinases: sequences for which classical classification based entirely on catalytic domains is inconsistent with the full-length similarity scores computed here, which implicitly consider multi-domain nature and regions outside the catalytic kinase domain. We also provide some examples of hybrid kinases of the protozoan parasite Entamoeba histolytica. Conclusions: The implicit consideration of multi-domain architectures is a valuable inclusion to complement other classification schemes. The proposed algorithm may also be employed to classify other families of enzymes with multidomain architecture.

Modelling initial survival and growth of radiata pine in New Zealand.

A sensitive framework has been developed for modelling young radiata pine survival, its growth and its size class distribution, from time of planting to age 5 or 6 years. The data and analysis refer to the Central North Island region of New Zealand. The survival function is derived from a Weibull probability density function, to reflect diminishing mortality with the passage of time in young stands. An anamorphic family of trends was used, as very little between-tree competition can be expected in young stands. An exponential height function was found to fit best the lower portion of its sigmoid form. The most appropriate basal area/ha exponential function included an allometric adjustment which resulted in compatible mean height and basal area/ha models. Each of these equations successfully represented the effects of several establishment practices by making coefficients linear functions of site factors, management activities and their interactions. Height and diameter distribution modelling techniques that ensured compatibility with stand values were employed to represent the effects of management practices on crop variation. Model parameters for this research were estimated using data from site preparation experiments in the region and were tested with some independent data sets.

Euclideanization of Majorana and Weyl Fermions

A continuous procedure is presented for euclideanization of Majorana and Weyl fermions without doubling their degrees of freedom. The Euclidean theory so obtained is SO(4) invariant and Osterwalder-Schrader (OS) positive. This enables us to define a one-complex parameter family of the N=1 supersymmetric Yang-Mills (SSYM) theories which interpolate between the Minkowski and a Euclidean SSYM theory. The interpolating action, and hence the Euclidean action, manifests all the continous symmetries of the original Minkowski space theory.

Regulation and Function of GATA Transcription Factors in Adrenocortical Tumors and Granulosa Cells

Transcription factors play a key role in tumor development, in which dysfunction of genes regulating tissue growth and differentiation is a central phenomenon. The GATA family of transcription factors consists of six members that bind to a consensus DNA sequence (A/T)GATA(A/G) in gene promoters and enhancers. The two GATA factors expressed in the adrenal cortex are GATA-4 and GATA-6. In both mice and humans, GATA-4 can be detected only during the fetal period, whereas GATA-6 expression is abundant both throughout development and in the adult. It is already established that GATA factors are important in both normal development and tumorigenesis of several endocrine organs, and expression of GATA-4 and GATA-6 is detected in adrenocortical tumors. The aim of this study was to elucidate the function of these factors in adrenocortical tumor growth. In embryonal development, the adrenocortical cells arise and differentiate from a common pool with gonadal steroidogenic cells, the urogenital ridge. As the adult adrenal cortex undergoes constant renewal, it is hypothesized that undifferentiated adrenocortical progenitor cells reside adjacent to the adrenal capsule and give rise to daughter cells that differentiate and migrate centripetally. A diverse array of hormones controls the differentiation, growth and survival of steroidogenic cells in the adrenal gland and the gonads. Factors such as luteinizing hormone and inhibins, traditionally associated with gonadal steroidogenic cells, can also influence the function of adrenocortical cells in physiological and pathophysiological states. Certain inbred strains of mice develop subcapsular adrenocortical tumors in response to gonadectomy. In this study, we found that these tumors express GATA-4, normally absent from the adult adrenal cortex, while GATA-6 expression is downregulated. Gonadal markers such as luteinizing hormone receptor, anti-Müllerian hormone and P450c17 are also expressed in the neoplastic cells, and the tumors produce gonadal hormones. The tumor cells have lost the expression of melanocortin-2 receptor and the CYP enzymes necessary for the synthesis of corticosterone and aldosterone. By way of xenograft studies utilizing NU/J nude mice, we confirmed that chronic gonadotropin elevation is sufficient to induce adrenocortical tumorigenesis in susceptible inbred strains. Collectively, these studies suggest that subcapsular adrenocortical progenitor cells can, under certain conditions, adopt a gonadal fate. We studied the molecular mechanisms involved in gene regulation in endocrine cells in order to elucidate the role of GATA factors in endocrine tissues. Ovarian granulosa cells express both GATA-4 and GATA-6, and the TGF-β signaling pathway is active in these cells. Inhibin-α is both a target gene for, and an atypical or antagonistic member of the TGF-β growth factor superfamily. In this study, we show that GATA-4 is required for TGF-β-mediated inhibin-α promoter activation in granulosa cells, and that GATA-4 physically interacts with Smad3, a TGF-β downstream protein. Apart from the regulation of steroidogenesis and other events in normal tissues, TGF-β signaling is implicated in tumors of multiple organs, including the adrenal cortex. Another signaling pathway found often to be aberrantly active in adrenocortical tumors is the Wnt pathway. As both of these pathways regulate the expression of inhibin-α, a transcriptional target for GATA-4 and GATA-6, we wanted to investigate whether GATA factors are associated with the components of these signaling cascades in human adrenocortical tumors. We found that the expression of Wnt co-receptors LRP5 and LRP6, Smad3, GATA-6 and SF-1 was diminished in adrenocortical carcinomas with poor outcome. All of these factors drive inhibin-α expression, and their expression in adrenocortical tumors correlated with that of inhibin-α. The results support a tumor suppressor role previously suggested for inhibin-α in the mouse adrenal cortex, and offer putative pathways associated with adrenocortical tumor aggressiveness. Unraveling the role of GATA factors and associated molecules in human and mouse adrenocortical tumors could ultimately contribute to the development of diagnostic tools and future therapies for these diseases.