Dynamic simulation of regulatory networks using SQUAD.

BACKGROUND: The ambition of most molecular biologists is the understanding of the intricate network of molecular interactions that control biological systems. As scientists uncover the components and the connectivity of these networks, it becomes possible to study their dynamical behavior as a whole and discover what is the specific role of each of their components. Since the behavior of a network is by no means intuitive, it becomes necessary to use computational models to understand its behavior and to be able to make predictions about it. Unfortunately, most current computational models describe small networks due to the scarcity of kinetic data available. To overcome this problem, we previously published a methodology to convert a signaling network into a dynamical system, even in the total absence of kinetic information. In this paper we present a software implementation of such methodology. RESULTS: We developed SQUAD, a software for the dynamic simulation of signaling networks using the standardized qualitative dynamical systems approach. SQUAD converts the network into a discrete dynamical system, and it uses a binary decision diagram algorithm to identify all the steady states of the system. Then, the software creates a continuous dynamical system and localizes its steady states which are located near the steady states of the discrete system. The software permits to make simulations on the continuous system, allowing for the modification of several parameters. Importantly, SQUAD includes a framework for perturbing networks in a manner similar to what is performed in experimental laboratory protocols, for example by activating receptors or knocking out molecular components. Using this software we have been able to successfully reproduce the behavior of the regulatory network implicated in T-helper cell differentiation. CONCLUSION: The simulation of regulatory networks aims at predicting the behavior of a whole system when subject to stimuli, such as drugs, or determine the role of specific components within the network. The predictions can then be used to interpret and/or drive laboratory experiments. SQUAD provides a user-friendly graphical interface, accessible to both computational and experimental biologists for the fast qualitative simulation of large regulatory networks for which kinetic data is not necessarily available.

Characterization of the Aspergillus nidulans biotin biosynthetic gene cluster and use of the bioDA gene as a new transformation marker.

The genes involved in the biosynthesis of biotin were identified in the hyphal fungus Aspergillus nidulans through homology searches and complementation of Escherichia coli biotin-auxotrophic mutants. Whereas the 7,8-diaminopelargonic acid synthase and dethiobiotin synthetase are encoded by distinct genes in bacteria and the yeast Saccharomyces cerevisiae, both activities are performed in A. nidulans by a single enzyme, encoded by the bifunctional gene bioDA. Such a bifunctional bioDA gene is a genetic feature common to numerous members of the ascomycete filamentous fungi and basidiomycetes, as well as in plants and oömycota. However, unlike in other eukaryota, the three bio genes contributing to the four enzymatic steps from pimeloyl-CoA to biotin are organized in a gene cluster in pezizomycotina. The A. nidulans auxotrophic mutants biA1, biA2 and biA3 were all found to have mutations in the 7,8-diaminopelargonic acid synthase domain of the bioDA gene. Although biotin auxotrophy is an inconvenient marker in classical genetic manipulations due to cross-feeding of biotin, transformation of the biA1 mutant with the bioDA gene from either A. nidulans or Aspergillus fumigatus led to the recovery of well-defined biotin-prototrophic colonies. The usefulness of bioDA gene as a novel and robust transformation marker was demonstrated in co-transformation experiments with a green fluorescent protein reporter, and in the efficient deletion of the laccase (yA) gene via homologous recombination in a mutant lacking non-homologous end-joining activity.

Leukemic cluster growth in culture is an independent risk factor for acute myeloid leukemia and short survival in patients with myelodysplastic syndrome

In patients with myelodysplastic syndrome (MDS) precursor cell cultures (colony-forming unit cells, CFU-C) can provide an insight into the growth potential of malignant myeloid cells. In a retrospective single-center study of 73 untreated MDS patients we assessed whether CFU-C growth patterns were of prognostic value in addition to established criteria. Abnormalities were classified as qualitative (i.e. leukemic cluster growth) or quantitative (i.e. strongly reduced/absent growth). Thirty-nine patients (53%) showed leukemic growth, 26 patients (36%) had strongly reduced/absent colony growth, and 12 patients showed both. In a univariate analysis the presence of leukemic growth was associated with strongly reduced survival (at 10 years 4 vs. 34%, p = 0.004), and a high incidence of transformation to AML (76 vs. 32%, p = 0.01). Multivariate analysis identified leukemic growth as a strong and independent predictor of early death (relative risk 2.12, p = 0.03) and transformation to AML (relative risk 2.63, p = 0.04). Quantitative abnormalities had no significant impact on the disease course. CFU- C assays have significant predictive value in addition to established prognostic factors in MDS. Leukemic growth identifies a subpopulation of MDS patients with poor prognosis.

Universal screening and decolonization for control of MRSA in nursing homes : a cluster randomized controlled study

OBJECTIVE The risk of carrying methicillin-resistant Staphylococcus aureus (MRSA) is higher among nursing home (NH) residents than in the general population. However, control strategies are not clearly defined in this setting. In this study, we compared the impact of standard precautions either alone (control) or combined with screening of residents and decolonization of carriers (intervention) to control MRSA in NHs. DESIGN Cluster randomized controlled trial SETTING NHs of the state of Vaud, Switzerland PARTICIPANTS Of 157 total NHs in Vaud, 104 (67%) participated in the study. INTERVENTION Standard precautions were enforced in all participating NHs, and residents underwent MRSA screening at baseline and 12 months thereafter. All carriers identified in intervention NHs, either at study entry or among newly admitted residents, underwent topical decolonization combined with environmental disinfection, except in cases of MRSA infection, MRSA bacteriuria, or deep skin ulcers. RESULTS NHs were randomly allocated to a control group (51 NHs, 2,412 residents) or an intervention group (53 NHs, 2,338 residents). Characteristics of NHs and residents were similar in both groups. The mean screening rates were 86% (range, 27%-100%) in control NHs and 87% (20%-100%) in intervention NHs. Prevalence of MRSA carriage averaged 8.9% in both control NHs (range, 0%-43%) and intervention NHs (range, 0%-38%) at baseline, and this rate significantly declined to 6.6% in control NHs and to 5.8% in intervention NHs after 12 months. However, the decline did not differ between groups (P=.66). CONCLUSION Universal screening followed by decolonization of carriers did not significantly reduce the prevalence of the MRSA carriage rate at 1 year compared with standard precautions

A Cluster Randomised Controlled Trial to Evaluate a Policy of Making Hip Protectors Available to Residents of Nursing Homes

Objectives: to evaluate the effectiveness of a policy of making hip protectors available to residents of nursing homes. Design: a cluster randomised controlled trial of the policy in nursing and residential homes, with the home as the unit of randomisation. Setting: 127 nursing and residential homes in the greater Belfast area of Northern Ireland. Participants: 40 homes in the intervention group (representing 1,366 occupied beds) and 87 homes in the control group (representing 2,751 occupied beds). Interventions: a policy of making hip protectors available free of charge to residents of nursing homes and supporting the implementation process by employing a nurse facilitator to encourage staff in the homes to promote their use, over a 72-week period. Main outcome measures: the rate of hip fractures in intervention and control homes, and the level of adherence to use of hip protectors. Results: there were 85 hip fractures in the intervention homes and 163 in the control homes. The mean fracture rate per 100 residents was 6.22 in the intervention homes and 5.92 in the control homes, giving an adjusted rate ratio for the intervention group compared to the control group of 1.05 (95% CI 0.77, 1.43, P = 0.76). Initial acceptance of the hip protectors was 37.2% (508/1,366) with adherence falling to 19.9% (272/1,366) at 72 weeks.Conclusions: making hip protectors available to residents of nursing and residential homes did not reduce the rate of hip fracture. This research does not support the introduction of a policy of providing hip protectors to residents of nursing homes.

Characterization of the hcnABC gene cluster encoding hydrogen cyanide synthase and anaerobic regulation by ANR in the strictly aerobic biocontrol agent Pseudomonas fluorescens CHA0.

The secondary metabolite hydrogen cyanide (HCN) is produced by Pseudomonas fluorescens from glycine, essentially under microaerophilic conditions. The genetic basis of HCN synthesis in P. fluorescens CHA0 was investigated. The contiguous structural genes hcnABC encoding HCN synthase were expressed from the T7 promoter in Escherichia coli, resulting in HCN production in this bacterium. Analysis of the nucleotide sequence of the hcnABC genes showed that each HCN synthase subunit was similar to known enzymes involved in hydrogen transfer, i.e., to formate dehydrogenase (for HcnA) or amino acid oxidases (for HcnB and HcnC). These similarities and the presence of flavin adenine dinucleotide- or NAD(P)-binding motifs in HcnB and HcnC suggest that HCN synthase may act as a dehydrogenase in the reaction leading from glycine to HCN and CO2. The hcnA promoter was mapped by primer extension; the -40 sequence (TTGGC ... ATCAA) resembled the consensus FNR (fumarate and nitrate reductase regulator) binding sequence (TTGAT ... ATCAA). The gene encoding the FNR-like protein ANR (anaerobic regulator) was cloned from P. fluorescens CHA0 and sequenced. ANR of strain CHA0 was most similar to ANR of P. aeruginosa and CydR of Azotobacter vinelandii. An anr mutant of P. fluorescens (CHA21) produced little HCN and was unable to express an hcnA-lacZ translational fusion, whereas in wild-type strain CHA0, microaerophilic conditions strongly favored the expression of the hcnA-lacZ fusion. Mutant CHA21 as well as an hcn deletion mutant were impaired in their capacity to suppress black root rot of tobacco, a disease caused by Thielaviopsis basicola, under gnotobiotic conditions. This effect was most pronounced in water-saturated artificial soil, where the anr mutant had lost about 30% of disease suppression ability, compared with wild-type strain CHA0. These results show that the anaerobic regulator ANR is required for cyanide synthesis in the strictly aerobic strain CHA0 and suggest that ANR-mediated cyanogenesis contributes to the suppression of black root rot.

Production of full-length cDNA sequences by sequencing and analysis of expressed sequence tags from Schistosoma mansoni

The number of sequences generated by genome projects has increased exponentially, but gene characterization has not followed at the same rate. Sequencing and analysis of full-length cDNAs is an important step in gene characterization that has been used nowadays by several research groups. In this work, we have selected Schistosoma mansoni clones for full-length sequencing, using an algorithm that investigates the presence of the initial methionine in the parasite sequence based on the positions of alignment start between two sequences. BLAST searches to produce such alignments have been performed using parasite expressed sequence tags produced by Minas Gerais Genome Network against sequences from the database Eukaryotic Cluster of Orthologous Groups (KOG). This procedure has allowed the selection of clones representing 398 proteins which have not been deposited as S. mansoni complete CDS in any public database. Dedicated sequencing of 96 of such clones with reads from both 5' and 3' ends has been performed. These reads have been assembled using PHRAP, resulting in the production of 33 full-length sequences that represent novel S. mansoni proteins. These results shall contribute to construct a more complete view of the biology of this important parasite.

Leukemic cluster growth in culture is an independent risk factor for acute myeloid leukemia and short survival in patients with myelodysplastic syndrome.

In patients with myelodysplastic syndrome (MDS) precursor cell cultures (colony-forming unit cells, CFU-C) can provide an insight into the growth potential of malignant myeloid cells. In a retrospective single-center study of 73 untreated MDS patients we assessed whether CFU-C growth patterns were of prognostic value in addition to established criteria. Abnormalities were classified as qualitative (i.e. leukemic cluster growth) or quantitative (i.e. strongly reduced/absent growth). Thirty-nine patients (53%) showed leukemic growth, 26 patients (36%) had strongly reduced/absent colony growth, and 12 patients showed both. In a univariate analysis the presence of leukemic growth was associated with strongly reduced survival (at 10 years 4 vs. 34%, p = 0.004), and a high incidence of transformation to AML (76 vs. 32%, p = 0.01). Multivariate analysis identified leukemic growth as a strong and independent predictor of early death (relative risk 2.12, p = 0.03) and transformation to AML (relative risk 2.63, p = 0.04). Quantitative abnormalities had no significant impact on the disease course. CFU-C assays have a significant predictive value in addition to established prognostic factors in MDS. Leukemic growth identifies a subpopulation of MDS patients with poor prognosis.

Cluster analysis of behavioural and event-related potentials during a contingent negative variation paradigm in remitting-relapsing and benign forms of multiple sclerosis

Background: Event-related potentials (ERPs) may be used as a highly sensitive way of detecting subtle degrees of cognitive dysfunction. On the other hand, impairment of cognitive skills is increasingly recognised as a hallmark of patients suffering from multiple sclerosis (MS). We sought to determine the psychophysiological pattern of information processing among MS patients with the relapsing-remitting form of the disease and low physical disability considered as two subtypes: 'typical relapsing-remitting' (RRMS) and 'benign MS' (BMS). Furthermore, we subjected our data to a cluster analysis to determine whether MS patients and healthy controls could be differentiated in terms of their psychophysiological profile.Methods: We investigated MS patients with RRMS and BMS subtypes using event-related potentials (ERPs) acquired in the context of a Posner visual-spatial cueing paradigm. Specifically, our study aimed to assess ERP brain activity in response preparation (contingent negative variation -CNV) and stimuli processing in MS patients. Latency and amplitude of different ERP components (P1, eN1, N1, P2, N2, P3 and late negativity -LN) as well as behavioural responses (reaction time -RT; correct responses -CRs; and number of errors) were analyzed and then subjected to cluster analysis. Results: Both MS groups showed delayed behavioural responses and enhanced latency for long-latency ERP components (P2, N2, P3) as well as relatively preserved ERP amplitude, but BMS patients obtained more important performance deficits (lower CRs and higher RTs) and abnormalities related to the latency (N1, P3) and amplitude of ERPs (eCNV, eN1, LN). However, RRMS patients also demonstrated abnormally high amplitudes related to the preparation performance period of CNV (cCNV) and post-processing phase (LN). Cluster analyses revealed that RRMS patients appear to make up a relatively homogeneous group with moderate deficits mainly related to ERP latencies, whereas BMS patients appear to make up a rather more heterogeneous group with more severe information processing and attentional deficits. Conclusions: Our findings are suggestive of a slowing of information processing for MS patients that may be a consequence of demyelination and axonal degeneration, which also seems to occur in MS patients that show little or no progression in the physical severity of the disease over time.

Cl gene cluster encoding several small nucleolar RNAs: a comparison amongst trypanosomatids

Small nucleolar RNAs (snoRNAs) are small non-coding RNAs that modify RNA molecules such as rRNA and snRNA by guiding 2'-O-ribose methylation (C/D box snoRNA family) and pseudouridylation reactions (H/ACA snoRNA family). H/ACA snoRNAs are also involved in trans-splicing in trypanosomatids. The aims of this work were to characterise the Cl gene cluster that encodes several snoRNAs in Trypanosoma rangeli and compare it with clusters from Trypanosoma cruzi, Trypanosoma brucei, Leishmania major, Leishmania infantum, Leishmania braziliensis and Leptomonas collosoma. The T. rangeli Cl gene cluster is an 801 base pair (bp) repeat sequence that encodes three C/D (Cl1, Cl2 and Cl4) and three H/ACA (Cl3, Cl5 and Cl6) snoRNAs. In contrast to T. brucei, the Cl3 and Cl5 homologues have not been annotated in the Leishmania or T. cruzi genome projects (http//:www.genedb.org). Of note, snoRNA transcribed regions have a high degree of sequence identity among all species and share gene synteny. Collectively, these findings suggest that the Cl cluster could constitute an interesting target for therapeutic (gene silencing) or diagnostic intervention strategies (PCR-derived tools).

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).