Incorporation of enzyme concentrations into FBA and identification of optimal metabolic pathways

<p>Abstract</p> <p>Background</p> <p>In the present article, we propose a method for determining optimal metabolic pathways in terms of the level of concentration of the enzymes catalyzing various reactions in the entire metabolic network. The method, first of all, generates data on reaction fluxes in a pathway based on steady state condition. A set of constraints is formulated incorporating weighting coefficients corresponding to concentration of enzymes catalyzing reactions in the pathway. Finally, the rate of yield of the target metabolite, starting with a given substrate, is maximized in order to identify an optimal pathway through these weighting coefficients.</p> <p>Results</p> <p>The effectiveness of the present method is demonstrated on two synthetic systems existing in the literature, two pentose phosphate, two glycolytic pathways, core carbon metabolism and a large network of carotenoid biosynthesis pathway of various organisms belonging to different phylogeny. A comparative study with the existing extreme pathway analysis also forms a part of this investigation. Biological relevance and validation of the results are provided. Finally, the impact of the method on metabolic engineering is explained with a few examples.</p> <p>Conclusions</p> <p>The method may be viewed as determining an optimal set of enzymes that is required to get an optimal metabolic pathway. Although it is a simple one, it has been able to identify a carotenoid biosynthesis pathway and the optimal pathway of core carbon metabolic network that is closer to some earlier investigations than that obtained by the extreme pathway analysis. Moreover, the present method has identified correctly optimal pathways for pentose phosphate and glycolytic pathways. It has been mentioned using some examples how the method can suitably be used in the context of metabolic engineering.</p

Characterisation of feline renal cortical fibroblast cultures and their transcriptional response to transforming growth factor beta 1

Chronic kidney disease (CKD) is common in geriatric cats, and the most prevalent pathology is chronic tubulointerstitial inflammation and fibrosis. The cell type predominantly responsible for the production of extra-cellular matrix in renal fibrosis is the myofibroblast, and fibroblast to myofibroblast differentiation is probably a crucial event. The cytokine TGF-β1 is reportedly the most important regulator of myofibroblastic differentiation in other species. The aim of this study was to isolate and characterise renal fibroblasts from cadaverous kidney tissue of cats with and without CKD, and to investigate the transcriptional response to TGF-β1

Predominant Functional Expression of Kv1.3 by Activated Microglia of the Hippocampus after Status epilepticus

BACKGROUND:Growing evidence indicates that the functional state of microglial cells differs according to the pathological conditions that trigger their activation. In particular, activated microglial cells can express sets of Kv subunits which sustain delayed rectifying potassium currents (Kdr) and modulate differently microglia proliferation and ability to release mediators. We recently reported that hippocampal microglia is in a particular activation state after a status epilepticus (SE) and the present study aimed at identifying which of the Kv channels are functionally expressed by microglia in this model. METHODOLOGY/PRINCIPAL FINDINGS:SE was induced by systemic injection of kainate in CX3CR1(eGFP/+) mice and whole cell recordings of fluorescent microglia were performed in acute hippocampal slices prepared 48 h after SE. Microglia expressed Kdr currents which were characterized by a potential of half-maximal activation near -25 mV, prominent steady-state and cumulative inactivations. Kdr currents were almost abolished by the broad spectrum antagonist 4-Aminopyridine (1 mM). In contrast, tetraethylammonium (TEA) at a concentration of 1 mM, known to block Kv3.1, Kv1.1 and 1.2 subunits, only weakly reduced Kdr currents. However, at a concentration of 5 mM which should also affect Kv1.3 and 1.6, TEA inhibited about 30% of the Kdr conductance. Alpha-dendrotoxin, which selectively inhibits Kv1.1, 1.2 and 1.6, reduced only weakly Kdr currents, indicating that channels formed by homomeric assemblies of these subunits are not important contributors of Kdr currents. Finally, agitoxin-2 and margatoxin strongly inhibited the current. CONCLUSIONS/SIGNIFICANCE:These results indicate that Kv1.3 containing channels predominantly determined Kdr currents in activated microglia after SE

The Exopolysaccharide Matrix Modulates the Interaction between 3D Architecture and Virulence of a Mixed-Species Oral Biofilm

Virulent biofilms are responsible for a range of infections, including oral diseases. All biofilms harbor a microbial-derived extracellular-matrix. The exopolysaccharides (EPS) formed on tooth-pellicle and bacterial surfaces provide binding sites for microorganisms; eventually the accumulated EPS enmeshes microbial cells. The metabolic activity of the bacteria within this matrix leads to acidification of the milieu. We explored the mechanisms through which the Streptococcus mutans-produced EPS-matrix modulates the three-dimensional (3D) architecture and the population shifts during morphogenesis of biofilms on a saliva-coated-apatitic surface using a mixed-bacterial species system. Concomitantly, we examined whether the matrix influences the development of pH-microenvironments within intact-biofilms using a novel 3D in situ pH-mapping technique. Data reveal that the production of the EPS-matrix helps to create spatial heterogeneities by forming an intricate network of exopolysaccharide-enmeshed bacterial-islets (microcolonies) through localized cell-to-matrix interactions. This complex 3D architecture creates compartmentalized acidic and EPS-rich microenvironments throughout the biofilm, which triggers the dominance of pathogenic S. mutans within a mixed-species system. The establishment of a 3D-matrix and EPS-enmeshed microcolonies were largely mediated by the S. mutans gtfB/gtfC genes, expression of which was enhanced in the presence of Actinomyces naeslundii and Streptococcus oralis. Acidic pockets were found only in the interiors of bacterial-islets that are protected by EPS, which impedes rapid neutralization by buffer (pH 7.0). As a result, regions of low pH (<5.5) were detected at specific locations along the surface of attachment. Resistance to chlorhexidine was enhanced in cells within EPS-microcolony complexes compared to those outside such structures within the biofilm. Our results illustrate the critical interaction between matrix architecture and pH heterogeneity in the 3D environment. The formation of structured acidic-microenvironments in close proximity to the apatite-surface is an essential factor associated with virulence in cariogenic-biofilms. These observations may have relevance beyond the mouth, as matrix is inherent to all biofilms

Parameter estimate of signal transduction pathways

BACKGROUND: The "inverse" problem is related to the determination of unknown causes on the bases of the observation of their effects. This is the opposite of the corresponding "direct" problem, which relates to the prediction of the effects generated by a complete description of some agencies. The solution of an inverse problem entails the construction of a mathematical model and takes the moves from a number of experimental data. In this respect, inverse problems are often ill-conditioned as the amount of experimental conditions available are often insufficient to unambiguously solve the mathematical model. Several approaches to solving inverse problems are possible, both computational and experimental, some of which are mentioned in this article. In this work, we will describe in details the attempt to solve an inverse problem which arose in the study of an intracellular signaling pathway. RESULTS: Using the Genetic Algorithm to find the sub-optimal solution to the optimization problem, we have estimated a set of unknown parameters describing a kinetic model of a signaling pathway in the neuronal cell. The model is composed of mass action ordinary differential equations, where the kinetic parameters describe protein-protein interactions, protein synthesis and degradation. The algorithm has been implemented on a parallel platform. Several potential solutions of the problem have been computed, each solution being a set of model parameters. A sub-set of parameters has been selected on the basis on their small coefficient of variation across the ensemble of solutions. CONCLUSION: Despite the lack of sufficiently reliable and homogeneous experimental data, the genetic algorithm approach has allowed to estimate the approximate value of a number of model parameters in a kinetic model of a signaling pathway: these parameters have been assessed to be relevant for the reproduction of the available experimental data

Expression Profiling of Major Histocompatibility and Natural Killer Complex Genes Reveals Candidates for Controlling Risk of Graft versus Host Disease

Background: The major histocompatibility complex (MHC) is the most important genomic region that contributes to the risk of graft versus host disease (GVHD) after haematopoietic stem cell transplantation. Matching of MHC class I and II genes is essential for the success of transplantation. However, the MHC contains additional genes that also contribute to the risk of developing acute GVHD. It is difficult to identify these genes by genetic association studies alone due to linkage disequilibrium in this region. Therefore, we aimed to identify MHC genes and other genes involved in the pathophysiology of GVHD by mRNA expression profiling. Methodology/Principal Findings: To reduce the complexity of the task, we used genetically well-defined rat inbred strains and a rat skin explant assay, an in-vitro-model of the graft versus host reaction (GVHR), to analyze the expression of MHC, natural killer complex (NKC), and other genes in cutaneous GVHR. We observed a statistically significant and strong up or down regulation of 11 MHC, 6 NKC, and 168 genes encoded in other genomic regions, i.e. 4.9%, 14.0%, and 2.6% of the tested genes respectively. The regulation of 7 selected MHC and 3 NKC genes was confirmed by quantitative real-time PCR and in independent skin explant assays. In addition, similar regulations of most of the selected genes were observed in GVHD-affected skin lesions of transplanted rats and in human skin explant assays. Conclusions/Significance: We identified rat and human MHC and NKC genes that are regulated during GVHR in skin explant assays and could therefore serve as biomarkers for GVHD. Several of the respective human genes, including HLA-DMB, C2, AIF1, SPR1, UBD, and OLR1, are polymorphic. These candidates may therefore contribute to the genetic risk of GVHD in patients

Tractography dissection variability: What happens when 42 groups dissect 14 white matter bundles on the same dataset?

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process

Future perspectives in melanoma research: meeting report from the "Melanoma Bridge";: Napoli, December 3rd-6th 2014.

The fourth "Melanoma Bridge Meeting" took place in Naples, December 3-6th, 2014. The four topics discussed at this meeting were: Molecular and Immunological Advances, Combination Therapies, News in Immunotherapy, and Tumor Microenvironment and Biomarkers. Until recently systemic therapy for metastatic melanoma patients was ineffective, but recent advances in tumor biology and immunology have led to the development of new targeted and immunotherapeutic agents that prolong progression-free survival (PFS) and overall survival (OS). New therapies, such as mitogen-activated protein kinase (MAPK) pathway inhibitors as well as other signaling pathway inhibitors, are being tested in patients with metastatic melanoma either as monotherapy or in combination, and all have yielded promising results. These include inhibitors of receptor tyrosine kinases (BRAF, MEK, and VEGFR), the phosphatidylinositol 3 kinase (PI3K) pathway [PI3K, AKT, mammalian target of rapamycin (mTOR)], activators of apoptotic pathway, and the cell cycle inhibitors (CDK4/6). Various locoregional interventions including radiotherapy and surgery are still valid approaches in treatment of advanced melanoma that can be integrated with novel therapies. Intrinsic, adaptive and acquired resistance occur with targeted therapy such as BRAF inhibitors, where most responses are short-lived. Given that the reactivation of the MAPK pathway through several distinct mechanisms is responsible for the majority of acquired resistance, it is logical to combine BRAF inhibitors with inhibitors of targets downstream in the MAPK pathway. For example, combination of BRAF/MEK inhibitors (e.g., dabrafenib/trametinib) have been demonstrated to improve survival compared to monotherapy. Application of novel technologies such sequencing have proven useful as a tool for identification of MAPK pathway-alternative resistance mechanism and designing other combinatorial therapies such as those between BRAF and AKT inhibitors. Improved survival rates have also been observed with immune-targeted therapy for patients with metastatic melanoma. Immune-modulating antibodies came to the forefront with anti-CTLA-4, programmed cell death-1 (PD-1) and PD-1 ligand 1 (PD-L1) pathway blocking antibodies that result in durable responses in a subset of melanoma patients. Agents targeting other immune inhibitory (e.g., Tim-3) or immune stimulating (e.g., CD137) receptors and other approaches such as adoptive cell transfer demonstrate clinical benefit in patients with melanoma as well. These agents are being studied in combination with targeted therapies in attempt to produce longer-term responses than those more typically seen with targeted therapy. Other combinations with cytotoxic chemotherapy and inhibitors of angiogenesis are changing the evolving landscape of therapeutic options and are being evaluated to prevent or delay resistance and to further improve survival rates for this patient population. This meeting's specific focus was on advances in combination of targeted therapy and immunotherapy. Both combination targeted therapy approaches and different immunotherapies were discussed. Similarly to the previous meetings, the importance of biomarkers for clinical application as markers for diagnosis, prognosis and prediction of treatment response was an integral part of the meeting. The overall emphasis on biomarkers supports novel concepts toward integrating biomarkers into contemporary clinical management of patients with melanoma across the entire spectrum of disease stage. Translation of the knowledge gained from the biology of tumor microenvironment across different tumors represents a bridge to impact on prognosis and response to therapy in melanoma

Seismic evidence for a lower-mantle origin of the Iceland plume

Iceland, one of the most thoroughly investigated hotspots1,2,3, is generally accepted to be the manifestation of an upwelling mantle plume4. Yet whether the plume originates from the lower mantle or from a convective instability at a thermal boundary layer between the upper and lower mantle near 660 km depth5,6 remains unconstrained. Tomographic inversions of body-wave delay times show that low seismic velocities extend to at least 400 km depth beneath central Iceland7,8, but cannot resolve structure at greater depth. Here we report lateral variations in the depths of compressional-to-shear wave conversions at the two seismic discontinuities marking the top and bottom of the mantle transition zone beneath Iceland. We find that the transition zone is 20 km thinner than in the average Earth9 beneath central and southern Iceland, but is of normal thickness beneath surrounding areas, a result indicative of a hot and narrow plume originating from the lower mantle.This work was supported by the US NSFPeer Reviewe

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline, and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams