Flux balance analysis of biological systems: applications and challenges

Systems level modelling and simulations of biological processes are proving to be invaluable in obtaining a quantitative and dynamic perspective of various aspects of cellular function. In particular, constraint-based analyses of metabolic networks have gained considerable popularity for simulating cellular metabolism, of which flux balance analysis (FBA), is most widely used. Unlike mechanistic simulations that depend on accurate kinetic data, which are scarcely available, FBA is based on the principle of conservation of mass in a network, which utilizes the stoichiometric matrix and a biologically relevant objective function to identify optimal reaction flux distributions. FBA has been used to analyse genome-scale reconstructions of several organisms; it has also been used to analyse the effect of perturbations, such as gene deletions or drug inhibitions in silico. This article reviews the usefulness of FBA as a tool for gaining biological insights, advances in methodology enabling integration of regulatory information and thermodynamic constraints, and finally addresses the challenges that lie ahead. Various use scenarios and biological insights obtained from FBA, and applications in fields such metabolic engineering and drug target identification, are also discussed. Genome-scale constraint-based models have an immense potential for building and testing hypotheses, as well as to guide experimentation.

Functional evaluation of genetic variants associated with endometriosis near GREB1

STUDY QUESTION: Do DNA variants in the growth regulation by estrogen in breast cancer 1 (GREB1) region regulate endometrial GREB1 expression and increase the risk of developing endometriosis in women? SUMMARY ANSWER: We identified new single nucleotide polymorphisms (SNPs) with strong association with endometriosis at the GREB1 locus although we did not detect altered GREB1 expression in endometriosis patients with defined genotypes. WHAT IS ALREADY KNOWN: Genome-wide association studies have identified the GREB1 region on chromosome 2p25.1 for increasing endometriosis risk. The differential expression of GREB1 has also been reported by others in association with endometriosis disease phenotype. STUDY DESIGN, SIZE, DURATION: Fine mapping studies comprehensively evaluated SNPs within the GREB1 region in a large-scale data set (>2500 cases and >4000 controls). Publicly available bioinformatics tools were employed to functionally annotate SNPs showing the strongest association signal with endometriosis risk. Endometrial GREB1 mRNA and protein expression was studied with respect to phases of the menstrual cycle (n = 2-45 per cycle stage) and expression quantitative trait loci (eQTL) analysis for significant SNPs were undertaken for GREB1 [mRNA (n = 94) and protein (n = 44) in endometrium]. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants in this study are females who provided blood and/or endometrial tissue samples in a hospital setting. The key SNPs were genotyped using Sequenom MassARRAY. The functional roles and regulatory annotations for identified SNPs are predicted by various publicly available bioinformatics tools. Endometrial GREB1 expression work employed qRT-PCR, western blotting and immunohistochemistry studies. MAIN RESULTS AND THE ROLE OF CHANCE: Fine mapping results identified a number of SNPs showing stronger association (0.004 < P < 0.032) with endometriosis risk than the original GWAS SNP (rs13394619) (P = 0.034). Some of these SNPs were predicted to have functional roles, for example, interaction with transcription factor motifs. The haplotype (a combination of alleles) formed by the risk alleles from two common SNPs showed significant association (P = 0.026) with endometriosis and epistasis analysis showed no evidence for interaction between the two SNPs, suggesting an additive effect of SNPs on endometriosis risk. In normal human endometrium, GREB1 protein expression was altered depending on the cycle stage (significantly different in late proliferative versus late secretory, P < 0.05) and cell type (glandular epithelium, not stromal cells). However, GREB1 expression in endometriosis cases versus controls and eQTL analyses did not reveal any significant changes. LIMITATIONS, REASONS FOR CAUTION: In silico prediction tools are generally based on cell lines different to our tissue and disease of interest. Functional annotations drawn from these analyses should be considered with this limitation in mind. We identified cell-specific and hormone-specific changes in GREB1 protein expression. The lack of a significant difference observed following our GREB1 expression studies may be the result of moderate power on mixed cell populations in the endometrial tissue samples. WIDER IMPLICATIONS OF THE FINDINGS: This study further implicates the GREB1 region on chromosome 2p25.1 and the GREB1 gene with involvement in endometriosis risk. More detailed functional studies are required to determine the role of the novel GREB1 transcripts in endometriosis pathophysiology. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants APP1012245, APP1026033, APP1049472 and APP1046880. There are no competing interests.

Latent TGF-β binding proteins -3 and -4 : transcriptional control and extracellular matrix targeting

Extracellular matrix (ECM) is a complex network of various proteins and proteoglycans which provides tissues with structural strength and resilience. By harvesting signaling molecules like growth factors ECM has the capacity to control cellular functions including proliferation, differentiation and cell survival. Latent transforming growth factor β (TGF-β) binding proteins (LTBPs) associate fibrillar structures of the ECM and mediate the efficient secretion and ECM deposition of latent TGF-β. The current work was conducted to determine the regulatory regions of LTBP-3 and -4 genes to gain insight into their tissue-specific expression which also has impact on TGF-β biology. Furthermore, the current research aimed at defining the ECM targeting of the N-terminal variants of LTBP-4 (LTBP-4S and -4L), which is required to understand their functions in tissues and to gain insight into conditions in which TGF-β is activated. To characterize the regulatory regions of LTBP-3 and -4 genes in silico and functional promoter analysis techniques were employed. It was found that the expression of LTBP-4S and -4L are under control of two independent promoters. This finding was in accordance with the observed expression patterns of LTBP-4S and -4L in human tissues. All promoter regions characterized in this study were TATAless, GC-rich and highly conserved between human and mouse species. Putative binding sites for Sp1 and GATA family of transcription factors were recognized in all of these regulatory regions. It is possible that these transcription factors control the basal expression of LTBP-3 and -4 genes. Smad binding element was found within the LTBP-3 and -4S promoter regions, but it was not present in LTBP-4L promoter. Although this element important for TGF-β signaling was present in LTBP-4S promoter, TGF-β did not induce its transcriptional activity. LTBP-3 promoter activity and mRNA expression instead were stimulated by TGF-β1 in osteosarcoma cells. It was found that the stimulatory effect of TGF-β was mediated by Smad and Erk MAPK signaling pathways. The current work explored the ECM targeting of LTBP-4S and identified binding partners of this protein. It was found that the N-terminal end of LTBP-4S possesses fibronectin (FN) binding sites which are critical for its ECM targeting. FN deficient fibroblasts incorporated LTBP-4S into their ECM only after addition of exogenous FN. Furthermore, LTBP-4S was found to have heparin binding regions, of which the C-terminal binding site mediated fibroblast adhesion. Soluble heparin prevented the ECM association of LTBP-4S in fibroblast cultures. In the current work it was observed that there are significant differences in the secretion, processing and ECM targeting of LTBP-4S and -4L. Interestingly, it was observed that most of the secreted LTBP-4L was associated with latent TGF-β1, whereas LTBP-4S was mainly secreted as a free form from CHO cells. This thesis provides information on transcriptional regulation of LTBP-3 and -4 genes, which is required for the deeper understanding of their tissue-specific functions. Further, the current work elucidates the structural variability of LTBPs, which appears to have impact on secretion and ECM targeting of TGF-β. These findings may advance understanding the abnormal activation of TGF-β which is associated with connective tissue disorders and cancer.

The intrinsic plasticity of farm businesses and their resilience to change. An Australian example.

This paper examines the idea that plasticity in farm management introduces resilience to change and allows farm businesses to perform when operating in highly variable environments. We also argue for the need to develop and apply more integrative assessments of farm performance that combine the use of modelling tools with deliberative processes involving farmers and researchers in a co-learning process, to more effectively identify and implement more productive and resilient farm businesses. In a plastic farming system, farm management is highly contingent on environmental conditions. In plastic farming systems farm managers constantly vary crops and inputs based on the availability of limited and variable resources (e.g. land, water, finances, labour, machinery, etc.), and signals from its operating environment (e.g. climate, markets), with the objective of maximising a number of, often competing, objectives (e.g. maximise profits, minimise risks, etc.). In contrast in more rigid farming systems farm management is more calendar driven and relatively fixed sequences of crops are regularly followed over time and across the farm. Here we describe the application of a whole farm simulation model to (i) compare, in silico, the sensitivity of two farming systems designs of contrasting levels of plasticity, operating in two contrasting environments, when exposed to a stressor in the form of climate change scenarios;(ii) investigate the presence of interactions and feedbacks at the field and farm levels capable of modifying the intensity and direction of the responses to climate signals; and (iii) discuss the need for the development and application of more integrative assessments in the analysis of impacts and adaptation options to climate change. In both environments, the more plastic farm management strategy had higher median profits and was less risky for the baseline and less intensive climate change scenarios (2030). However, for the more severe climate change scenarios (2070), the benefit of plastic strategies tended to disappear. These results suggest that, to a point, farming systems having higher levels of plasticity would enable farmers to more effectively respond to climate shifts, thus ensuring the economic viability of the farm business. Though, as the intensity of the stress increases (e.g. 2070 climate change scenario) more significant changes in the farming system might be required to adapt. We also found that in the case studies analysed here, most of the impacts from the climate change scenarios on farm profit and economic risk originated from important reductions in cropping intensity and changes in crop mix rather than from changes in the yields of individual crops. Changes in cropping intensity and crop mix were explained by the combination of reductions in the number of sowing opportunities around critical times in the cropping calendar, and to operational constraints at the whole farm level i.e. limited work capacity in an environment having fewer and more concentrated sowing opportunities. This indicates that indirect impacts from shifts in climate on farm operations can be more important than direct impacts from climate on the yield of individual crops. The results suggest that due to the complexity of farm businesses, impact assessments and opportunities for adaptation to climate change might also need to be pursued at higher integration levels than the crop or the field. We conclude that plasticity can be a desirable characteristic in farming systems operating in highly variable environments, and that integrated whole farm systems analyses of impacts and adaptation to climate change are required to identify important interactions between farm management decision rules, availability of resources, and farmer's preference.

The Pseudomonas syringae-derived HrpA pilins - molecular characterization and biotechnological application of the transcripts

Many Gram-negative bacteria pathogenic to plants and animals possess type III secretion systems that are used to cause disease. Effector proteins are injected into host cells using the type III secretion machineries. Despite vigorous studies, the nature of the secretion signal for type III secreted proteins still remains elusive. Both mRNA and proteinaceous signals have been proposed. Findings on coupling of translation to secretion by the type III secretion systems are also still contradictory. This study dealt with the secretion signal of HrpA from Pseudomonas syringae pathovar tomato. HrpA is the major component of the type III secretion system-associated Hrp pilus and a substrate for the type III secretion systems. The secretion signal was shown to reside in the first 15 codons or amino acids, a location typical for type III secretion signals. Translation of HrpA in the absence of a functional type III secretion system was established, but it does not exclude the possibility of coupling of translation to secretion when the secretion apparatus is present. The hrpA transcripts from various unrelated plant pathogenic bacteria were shown to be extremely stable. The biological relevance of this observation is unknown, but possible explanations include the high prevalence of HrpA protein, an mRNA secretion signal or timing of secretion. The hrpA mRNAs are stable over a wide range of temperatures, in the absence of translating ribosomes and even in the heterologous host Escherichia coli. The untranslated regions (UTRs) of hrpA transcripts from at least 20 pathovars of Pseudomonas syringae are highly homologous, whilst their coding regions exhibit low similarity. The stable nature of hrpA messenger RNAs is likely to be due to the folding of their 5 and 3 UTRs. In silico the UTRs seem to form stem-loop structures, the hairpin structures in the 3 UTRs being rich in guanidine and cytosine residues. The stable nature of the hrpA transcript redirected the studies to the stabilization of heterologous transcripts and to the use of stable messenger RNAs in recombinant protein production. Fragments of the hrpA transcript can be used to confer stability on heterologous transcripts from several sources of bacterial and eukaryotic origin, and to elevate the levels of production of the corresponding recombinant proteins several folds. hrpA transcript stabilizing elements can be used for improving the yields of recombinant proteins even in Escherichia coli, one of the most commonly used industrial protein production hosts.

Mechanisms of KCC2 upregulation during development

Gamma-aminobutyric acid (GABA) acting through ionotropic GABAA receptors plays a crucial role in the activity of the central nervous system (CNS). It triggers Ca2+ rise providing trophic support in developing neurons and conducts fast inhibitory function in mature neuronal networks. There is a developmental change in the GABAA reversal potential towards more negative levels during the first two postnatal weeks in rodent hippocampus. This change provides the basis for mature GABAergic activity and is attributable to the developmental expression of the neuron-specific potassium chloride cotransporter 2 (KCC2). In this work we have studied the mechanisms responsible for the control of KCC2 developmental expression. As a model system we used hippocampal dissociated cultures plated from embryonic day (E) 17 mice embryos before the onset of KCC2 expression. We showed that KCC2 was significantly up-regulated during the first two weeks of culture development. Interestingly, the level of KCC2 upregulation was not altered by chronic pharmacological blockage of action potentials as well as GABAergic and glutamatergic synaptic transmission. By in silico analysis of the proximal KCC2 promoter region we identified 10 candidate transcription factor binding sites that are highly conserved in mammalian KCC2 genes. One of these transcription factors, namely early growth response factor 4 (Egr4), had similar developmental profile as KCC2 and considerably increased the activity of mouse KCC2 gene in neuronal cells. Next we investigated the involvement of neurotrophic factors in regulation of Egr4 and KCC2 expression. We found that in immature hippocampal cultures Egr4 and KCC2 levels were strongly up-regulated by brain derived neurotrophic factor (BDNF)and neurturin. The effect of neurotrophic factors was dependent on the activation of a mitogen activated protein kinase (MAPK) signal transduction pathway. Intact Egr4-binding site in proximal KCC2 promoter was required for BDNF-induced KCC2 transcription. In vitro data were confirmed by several in vivo experiments where we detected an upregulation of KCC2 protein levels after intrahippocampal administration of BDNF or neurturin. Importantly, a MAPK-dependent rise in Egr4 and KCC2 expression levels was also observed after a period of kainic acid-induced seizure activity in neonatal rats suggesting that neuronal activity might be involved in Egr4-mediated regulation of KCC2 expression. Finally we demonstrated that the mammalian KCC2 gene (alias Slc12a5) generated two neuron-specific isoforms by using alternative promoters and first exons. A novel isoform of KCC2, termed KCC2a, differed from the previously known KCC2b isoform by 40 unique N-terminal amino acid residues. KCC2a expression was restricted to CNS,remained relatively constant during postnatal development, and contributed 20 50% of total KCC2 mRNA expression in the neonatal mouse brainstem and spinal cord. In summary, our data provide insight into the complex regulation of KCC2 expression during early postnatal development. Although basal KCC2 expression seems to be intrinsically regulated, it can be further augmented by neurotrophic factors or by enhanced activity triggering MAPK phosphorylation and Egr4 induction. Additional KCC2a isoform, regulated by another promoter, provides basal KCC2 level in neonatal brainstem and spinal cord required for survival of KCC2b knockout mice.

Characterisation, cloning and production of industrially interesting enzymes: gluconolactone oxidase of Penicillium cyaneo-fulvum and gluconate 5-dehydrogenase of Gluconobacter suboxydans

The work covered in this thesis is focused on the development of technology for bioconversion of glucose into D-erythorbic acid (D-EA) and 5-ketogluconic acid (5-KGA). The task was to show on proof-of-concept level the functionality of the enzymatic conversion or one-step bioconversion of glucose to these acids. The feasibility of both studies to be further developed for production processes was also evaluated. The glucose - D-EA bioconversion study was based on the use of a cloned gene encoding a D-EA forming soluble flavoprotein, D-gluconolactone oxidase (GLO). GLO was purified from Penicillium cyaneo-fulvum and partially sequenced. The peptide sequences obtained were used to isolate a cDNA clone encoding the enzyme. The cloned gene (GenBank accession no. AY576053) is homologous to the other known eukaryotic lactone oxidases and also to some putative prokaryotic lactone oxidases. Analysis of the deduced protein sequence of GLO indicated the presence of a typical secretion signal sequence at the N-terminus of the enzyme. No other targeting/anchoring signals were found, suggesting that GLO is the first known lactone oxidase that is secreted rather than targeted to the membranes of the endoplasmic reticulum or mitochondria. Experimental evidence supports this analysis, as near complete secretion of GLO was observed in two different yeast expression systems. Highest expression levels of GLO were obtained using Pichia pastoris as an expression host. Recombinant GLO was characterised and the suitability of purified GLO for the production of D-EA was studied. Immobilised GLO was found to be rapidly inactivated during D-EA production. The feasibility of in vivo glucose - D-EA conversion using a P. pastoris strain co-expressing the genes of GLO and glucose oxidase (GOD, E.C. 1.1.3.4) of A. niger was demonstrated. The glucose - 5-KGA bioconversion study followed a similar strategy to that used in the D-EA production research. The rationale was based on the use of a cloned gene encoding a membrane-bound pyrroloquinoline quinone (PQQ)-dependent gluconate 5-dehydrogenase (GA 5-DH). GA 5-DH was purified to homogeneity from the only source of this enzyme known in literature, Gluconobacter suboxydans, and partially sequenced. Using the amino acid sequence information, the GA 5-DH gene was cloned from a genomic library of G. suboxydans. The cloned gene was sequenced (GenBank accession no. AJ577472) and found to be an operon of two adjacent genes encoding two subunits of GA 5-DH. It turned out that GA 5-DH is a rather close homologue of a sorbitol dehydrogenase from another G. suboxydans strain. It was also found that GA 5-DH has significant polyol dehydrogenase activity. The G. suboxydans GA 5-DH gene was poorly expressed in E. coli. Under optimised conditions maximum expression levels of GA 5-DH did not exceed the levels found in wild-type G. suboxydans. Attempts to increase expression levels resulted in repression of growth and extensive cell lysis. However, the expression levels were sufficient to demonstrate the possibility of bioconversion of glucose and gluconate into 5-KGA using recombinant strains of E. coli. An uncharacterised homologue of GA 5-DH was identified in Xanthomonas campestris using in silico screening. This enzyme encoded by chromosomal locus NP_636946 was found by a sequencing project of X. campestris and named as a hypothetical glucose dehydrogenase. The gene encoding this uncharacterised enzyme was cloned, expressed in E. coli and found to encode a gluconate/polyol dehydrogenase without glucose dehydrogenase activity. Moreover, the X. campestris GA 5-DH gene was expressed in E. coli at nearly 30 times higher levels than the G. suboxydans GA 5-DH gene. Good expressability of the X. campestris GA-5DH gene makes it a valuable tool not only for 5-KGA production in the tartaric acid (TA) bioprocess, but possibly also for other bioprocesses (e.g. oxidation of sorbitol into L-sorbose). In addition to glucose - 5-KGA bioconversion, a preliminary study of the feasibility of enzymatic conversion of 5-KGA into TA was carried out. Here, the efficacy of the first step of a prospective two-step conversion route including a transketolase and a dehydrogenase was confirmed. It was found that transketolase convert 5-KGA into TA semialdehyde. A candidate for the second step was suggested to be succinic dehydrogenase, but this was not tested. The analysis of the two subprojects indicated that bioconversion of glucose to TA using X. campestris GA 5-DH should be prioritised first and the process development efforts in future should be focused on development of more efficient GA 5-DH production strains by screening a more suitable production host and by protein engineering.

Computational methods for the prediction of the structure and interactions of coiled-coil peptides

The past several years have seen significant advances in the development of computational methods for the prediction of the structure and interactions of coiled-coil peptides. These methods are generally based on pairwise correlations of amino acids, helical propensity, thermal melts and the energetics of sidechain interactions, as well as statistical patterns based on Hidden Markov Model (HMM) and Support Vector Machine (SVM) techniques. These methods are complemented by a number of public databases that contain sequences, motifs, domains and other details of coiled-coil structures identified by various algorithms. Some of these computational methods have been developed to make predictions of coiled-coil structure on the basis of sequence information; however, structural predictions of the oligomerisation state of these peptides still remains largely an open question due to the dynamic behaviour of these molecules. This review focuses on existing in silico methods for the prediction of coiled-coil peptides of functional importance using sequence and/or three-dimensional structural data.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Characterization of sequence and structural features of the Candida krusei Enolase

The incidence of human infections by the fungal pathogen Candida species has been increasing in recent years. Enolase is an essential protein in fungal metabolism. Sequence data is available for human and a number of medically important fungal species. An understanding of the structural and functional features of fungal enolases may provide the structural basis for their use as a target for the development of new anti-fungal drugs. We have obtained the sequence of the enolase of Candida krusei (C. krusei), as it is a significant medically important fungal pathogen. We have then used multiple sequence alignments with various enolase isoforms in order to identify C. krusei specific amino acid residues. The phylogenetic tree of enolases shows that the C. krusei enolase assembles on the tree with the fungal genes. Importantly, C. krusei lacks four amino acids in the active site compared to human enolase, as revealed by multiple sequence alignments. These differences in the substrate binding site may be exploited for the design of new anti-fungal drugs to selectively block this enzyme. The lack of the important amino acids in the active site also indicates that C. krusei enolase might have evolved as a member of a mechanistically diverse enolase superfamily catalying somewhat different reactions.

Molecular modeling of Bt Cry1Ac (DI–DII)–ASAL (Allium sativum lectin)–fusion protein and its interaction with aminopeptidase N (APN) receptor of Manduca sexta

Genetic engineering of Bacillus thuringiensis (Bt) Cry proteins has resulted in the synthesis of various novel toxin proteins with enhanced insecticidal activity and specificity towards different insect pests. In this study, a fusion protein consisting of the DI–DII domains of Cry1Ac and garlic lectin (ASAL) has been designed in silico by replacing the DIII domain of Cry1Ac with ASAL. The binding interface between the DI–DII domains of Cry1Ac and lectin has been identified using protein–protein docking studies. Free energy of binding calculations and interaction profiles between the Cry1Ac and lectin domains confirmed the stability of fusion protein. A total of 18 hydrogen bonds was observed in the DI–DII–lectin fusion protein compared to 11 hydrogen bonds in the Cry1Ac (DI–DII–DIII) protein. Molecular mechanics/Poisson–Boltzmann (generalized-Born) surface area [MM/PB (GB) SA] methods were used for predicting free energy of interactions of the fusion proteins. Protein–protein docking studies based on the number of hydrogen bonds, hydrophobic interactions, aromatic–aromatic, aromatic–sulphur, cation–pi interactions and binding energy of Cry1Ac/fusion proteins with the aminopeptidase N (APN) of Manduca sexta rationalised the higher binding affinity of the fusion protein with the APN receptor compared to that of the Cry1Ac–APN complex, as predicted by ZDOCK, Rosetta and ClusPro analysis. The molecular binding interface between the fusion protein and the APN receptor is well packed, analogously to that of the Cry1Ac–APN complex. These findings offer scope for the design and development of customized fusion molecules for improved pest management in crop plants.

Structure-Based Design of DevR Inhibitor Active against Nonreplicating Mycobacterium tuberculosis

Antitubercular treatment is directed against actively replicating organisms. There is an urgent need to develop drugs targeting persistent subpopulations of Mycobacterium tuberculosis. The DevR response regulator is believed to play a key role in bacterial dormancy adaptation during hypoxia. We developed a homology-based model of DevR and used it for the rational design of inhibitors. A phenylcoumarin derivative (compound 10) identified by in silico pharmacophore-based screening of 2.5 million compounds employing protocols with some novel features including a water-based pharmacophore query, was characterized further. Compound 10 inhibited DevR binding to target DNA, down-regulated dormancy genes transcription, and drastically reduced survival of hypoxic but not nutrient-starved dormant bacteria or actively growing organ ` isms. Our findings suggest that compound 10 ``locks'' DevR in an inactive conformation that is unable to bind cognate DNA and induce the dormancy regulon. These results provide proof-of-concept for DevR as a novel target to develop molecules with sterilizing activity against tubercle bacilli.

Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition

The rapid increase in genome sequence information has necessitated the annotation of their functional elements, particularly those occurring in the non-coding regions, in the genomic context. Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. In this analysis, we demonstrate that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. We have therefore obtained a set of free energy threshold values, for genomic DNA with varying GC content and used them as generic criteria for predicting promoter regions in several microbial genomes, using an in-house developed tool `PromPredict'. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (50.8% GC) and B. subtilis (43.5% GC) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10-8) level: 14q24.2 (rs227425, P-value 3.98 × 10-13, SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10-9, CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n 5 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis. © The Author 2013. Published by Oxford University Press.

Functional characterization of lysophosphatidic acid phosphatase from Arabidopsis thaliana

Lysophosphatidic acid (LPA) acts as a signaling molecule that regulates diverse cellular processes and it can rapidly be metabolized by phosphatase and acyltransferase LPA phosphatase gene has not been identified and characterized in plants so far The BLAST search revealed that the At3g03520 is similar to phospholipase family. and distantly related to bacterial phosphatases The conserved motif. (J)4XXXNXSFD, was identified in both At3g03520 like phospholipases and acid phosphatases In silico expression analysis of At3g03520 revealed a high expression during phosphate starvation and abiotic stresses. This gene was overexpressed in Escherichia coli and shown to posses LPA specific phosphatase activity These results Suggest that this gene possibly plays a role in signal transduction and storage lipid synthesis.