Discrete modeling for integration and analysis of large-scale signaling networks

Most biological processes are orchestrated by large-scale molecular networks which are described in large-scale model repositories and whose dynamics are extremely complex. An observed phenotype is a state of this system that results from control mechanisms whose identification is key to its understanding. The Biological Pathway Exchange (BioPAX) format is widely used to standardize the biological information relative to regulatory processes. However, few modeling approaches developed so far enable for computing the events that control a phenotype in large-scale networks.Here we developed an integrated approach to build large-scale dynamic networks from BioPAX knowledge databases in order to analyse trajectories and to identify sets of biological entities that control a phenotype. The Cadbiom approach relies on the guarded transitions formalism, a discrete modeling approach which models a system dynamics by taking into account competition and cooperation events in chains of reactions. The method can be applied to every BioPAX (large-scale) model thanks to a specific package which automatically generates Cadbiom models from BioPAX files.The Cadbiom framework was applied to the BioPAX version of two resources (PID, KEGG) of the Pathway Commons database and to the Atlas of Cancer Signalling Network (ACSN). As a case-study, it was used to characterize sets of biological entities implicated in the epithelial-mesenchymal transition. Our results highlight the similarities between the PID and ACSN resources in terms of biological content, and underline the heterogeneity of usage of the BioPAX semantics limiting the fusion of models that require curation. Causality analyses demonstrate the smart complementarity of the databases in terms of combinatorics of controllers that explain a phenotype. From a biological perspective, our results show the specificity of controllers for epithelial and mesenchymal phenotypes that are consistent with the literature and identify a novel signature for intermediate states.     

Lipid content and apoptosis of in vitro-produced bovine embryos as determinants of susceptibility to vitrification

The objective was to evaluate supplementation of fetal calf serum (FCS) and phenazine ethosulfate (PES), a metabolic regulator that inhibits fatty acid synthesis, in culture media during in vitro production (IVP) of bovine embryos. Taking oocyte fertilization (n = 4,320) as Day 0, four concentrations of FCS (0, 2.5, 5, and 10%) and three periods of exposure to PES (without addition—Control; after 60 h—PES Day 2.5 of embryo culture; and after 96 h—PES Day 4) were evaluated. Increasing FCS concentration in the culture media enhanced lipid accumulation (P < 0.05), increased apoptosis in fresh (2.5%: 19.1 ± 1.8 vs 10%: 28.4 ± 2.3, P < 0.05; mean ± SEM) and vitrified (2.5%: 42.8 ± 2.7 vs 10%: 69.2 ± 3.4, P < 0.05) blastocysts, and reduced blastocoele re-expansion after vitrification (2.5%: 81.6 ± 2.5 vs 10%: 67.3 ± 3.5, P < 0.05). The addition of PES in culture media, either from Days 2.5 or 4, reduced lipid accumulation (P < 0.05) and increased blastocoele re-expansion after vitrification (Control: 72.0 ± 3.0 vs PES Day 2.5: 79.9 ± 2.8 or PES Day 4: 86.2 ± 2.4, P < 0.05). However, just the use of PES from D4 reduced apoptosis in vitrified blastocysts (Control: 52.0 ± 3.0 vs PES Day 4: 39.2 ± 2.4, P < 0.05). Independent of FCS withdrawal or PES addition to culture media, the in vivo control group had lesser lipid accumulation, a lower apoptosis rate, and greater cryotolerance (P < 0.05). The increased lipid content was moderately correlated with apoptosis in vitrified blastocysts (r = 0.64, P = 0.01). In contrast, the increased apoptosis in fresh blastocysts was strongly correlated with apoptosis in vitrified blastocysts (r = 0.94, P < 0.0001). Therefore, using only 2.5% FCS and the addition of PES from Day 4, increased the survival of IVP embryos after vitrification. Moreover, embryo quality, represented by the fresh apoptosis rate, was better than lipid content for predicting embryo survival after vitrification.     

Glucose respiration and fermentation in Zygosaccharomyces bailii and Saccharomyces cerevisiae express different sensitivity patterns to ethanol and acetic acid

In the yeast Zygosaccharomyces bailii ISA 1307, respiration and fermentation ofglucose were exponentially inhibited by ethanol, both processes displaying similar sensitivity tothe alcohol. Moreover, the degree of inhibition on fermentation was of the same magnitude as thatreported for Saccharomyces cerevisiae. Acetic acid also inhibited these two metabolicprocesses in Z. bailii , with the kinetics of inhibition again being exponential. However,inhibition of fermentation was much less pronounced than in S. cerevisiae. The valuesestimated with Z. bailii for the minimum inhibitory concentration of acetic acid rangedfrom 100 to 240 mmol l⁻¹ total acetic acid compared with values of near zeroreported for S. cerevisiae. The inhibitory effects of acetic acid on Z. bailii were notsignificantly potentiated by ethanol.     

Role of histidine-50, glutamic acid-96, and histidine-137 in the ribonucleolytic mechanism of the ribotoxin alpha-sarcin

alpha-Sarcin is a ribotoxin secreted by the mold Aspergillus giganteus that degrades the ribosomal RNA by acting as a cyclizing ribonuclease. Three residues potentially involved in the mechanism of catalysis--histidine-50, glutamic acid-96, and histidine-137--were changed to glutamine. Three different single mutation variants (H50Q, E96Q, H137Q) as well as a double variant (H50/137Q) and a triple variant (H50/137Q/E96Q) were prepared and isolated to homogeneity. These variants were spectroscopically (circular dichroism, fluorescence emission, and proton nuclear magnetic resonance) characterized. According to these results, the three-dimensional structure of these variants of alpha-sarcin was preserved; only very minor local changes were detected. All the variants were inactive when assayed against either intact ribosomes or poly(A). The effect of pH on the ribonucleolytic activity of alpha-sarcin was evaluated against the ApA dinucleotide. This assay revealed that only the H50Q variant still retained its ability to cleave a phosphodiester bond, but it did so to a lesser extent than did wild-type alpha-sarcin. The results obtained are interpreted in terms of His137 and Glu96 as essential residues for the catalytic activity of alpha-sarcin (His137 as the general acid and Glu96 as the general base) and His50 stabilizing the transition state of the reaction catalyzed by alpha-sarcin.     

Characterization of the molecular properties of KtrC,a second RCK domain that regulates a Ktr channel in Bacillus subtilis

RCK (regulating conductance of K⁺) domains are common regulatory domains that control the activity of eukaryotic and prokaryotic K⁺ channels and transporters. In bacteria these domains play roles in osmoregulation, regulation of turgor and membrane potential and in pH homeostasis. Whole-genome sequencing unveiled RCK gene redundancy, however the biological role of this redundancy is not well understood. In Bacillus subtilis, there are two closely related RCK domain proteins (KtrA and KtrC) that regulate the activity of the Ktr cation channels. KtrA has been well characterized but little is known about KtrC. We have characterized the structural and biochemical proprieties of KtrC and conclude that KtrC binds ATP and ADP, just like KtrA. However, in solution KtrC exist in a dynamic equilibrium between octamers and non-octameric species that is dependent on the bound ligand, with ATP destabilizing the octameric ring relative to ADP. Accordingly, KtrC-ADP crystal structures reveal closed octameric rings similar to those in KtrA, while KtrC-ATP adopts an open assembly with RCK domains forming a super-helix. In addition, both KtrC-ATP and -ADP octamers are stabilized by the signaling molecule cyclic-di-AMP, which binds to KtrC with high affinity. In contrast, c-di-AMP binds with 100-fold lower affinity to KtrA. Despite these differences we show with an E. coli complementation assay that KtrC and KtrA are interchangeable and able to form functional transporters with both KtrB and KtrD. The distinctive properties of KtrC, in particular ligand-dependent assembly/disassembly, suggest that this protein has a specific physiological role that is distinct from KtrA.     

Impact of climatic variation on populations of pine processionary moth Thaumetopoea pityocampa in a core area of its distribution

• 1 There is growing appreciation of climatic effects on insect population dynamics at the margins of distribution limits. Climatic effects might be less important and/or involve different drivers and processes near the centre of distributions.
• 2 We evaluated the effects of interannual variation in temperatures, radiation and precipitation on populations of pine processionary moth Thaumetopoea pityocampa in Central–South Portugal, a low altitude area with a relatively mild Mediterranean climate near the centre of the north–south range of the species.
• 3 We tested for effects of climate on mortality of young larvae, growth rates, final larval mass and fecundity.
• 4 Results indicated high mortality of early instars associated with low minimum and maximum daily temperatures and low precipitation. Low minimum temperatures were further associated with high parasitism by the larval parasitoid Phryxe caudata (Rondani) (Diptera, Tachinidae). Furthermore, larval growth rates were higher with high solar radiation during December and January, which was itself negatively related to precipitation and air temperature. Slow larval growth rates led to lower final mass at pupation in the spring, and smaller egg masses and smaller initial colony sizes during the next autumn.
• 5 Thus climatic factors, and temperature in particular, apparently contributed to population dynamics of T. pityocampa in the core of its distribution, as well as at its northern limits. The most specific climatic parameters of importance, however, and the connections between climate, physiology and insect demographics in the core area were clearly different from northern areas.

     

Carbon-13 nuclear magnetic resonance analysis of [1-13C]glucose metabolism in Crithidia fasciculata. Evidence of CO2 fixation by phosphoenolpyruvate carboxykinase

The non-invasive technique of 13C nuclear magnetic resonance was applied to study glucose metabolism in vivo in the insect parasite Crithidia fasciculata. It was found that under anaerobic conditions [1-13C]glucose underwent a glycolytic pathway whose main metabolic products were identified as [2-13C]ethanol, [2-13C]succinate and [1,3-13C2]glycerol. These metabolites were excreted by C. fasciculata into the incubation medium, while in the cells [3-13C]phosphoenolpyruvate was also detected in addition to the aforementioned compounds. The C3 acid is apparently the acceptor of the primary CO2 fixation reaction, which leads in Trypanosomatids to the synthesis of succinate. By addition of sodium bicarbonate to the incubation mixture L-[3-13C]malate was detected among the excretion products, while the ethanol:succinate ratio of 2.0 in the absence of bicarbonate changed to a ratio of 0.6 in the presence of the latter. This was due to a shift of the balance between carboxylation of phosphoenolpyruvate, leading to succinate, and pyruvate decarboxylation leading to ethanol. The addition of 25% 2H2O to the incubation mixture led to the formation of [2-13C, 2-2H]ethanol derived from the prior incorporation of 2H+ into pyruvate in the reactions mediated by either pyruvate kinase or malic enzyme. However, no 2H+ incorporation into L-malate was detected, excluding the possibility that the latter was formed by carboxylation of pyruvate, and lending support to the idea that L-malate results from the carboxylation of phosphoenolpyruvate to oxaloacetate by phosphoenolpyruvate carboxykinase. The formation of [2-13C, 2-2H]-succinate under the same conditions reflected the uptake of 2H+ during the reduction of fumarate. When the incubations were carried out in the presence of 100% 2H2O, several [1-13C, 1-2H]ethanol species were detected, as well as [2-13C, 2-2H]malate and [1,3-13C2, 1,3-2H2]glycerol. The former deuterated compounds reflect the existence of NAD2H species when the incubations were carried out in 100% 2H2O, while the incorporation of 2H+ into [1,3-13C2]glycerol must be attributed to the phosphoglucose-isomerase-mediated reaction during glycolysis.