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.

Identification of early metabolic markers of various degrees of ischemia in the mouse brain by localized H-1 magnetic resonance spectroscopy

Objectives: Magnetic resonance (MR) imaging and spectroscopy (MRS) allow the establishment of the anatomical evolution and neurochemical profiles of ischemic lesions. The aim of the present study was to identify markers of reversible and irreversible damage by comparing the effects of 10-mins middle cerebral artery occlusion (MCAO), mimicking a transient ischemic attack, with the effects of 30-mins MCAO, inducing a striatal lesion. Methods: ICR-CD1 mice were subjected to 10-mins (n = 11) or 30-mins (n = 9) endoluminal MCAO by filament technique at 0 h. The regional cerebral blood flow (CBF) was monitored in all animals by laser- Doppler flowmetry with a flexible probe fixed on the skull with < 20% of baseline CBF during ischemia and > 70% during reperfusion. All MR studies were carried out in a horizontal 14.1T magnet. Fast spin echo images with T2-weighted parameters were acquired to localize the volume of interest and evaluate the lesion size. Immediately after adjustment of field inhomogeneities, localized 1H MRS was applied to obtain the neurochemical profile from the striatum (6 to 8 microliters). Six animals (sham group) underwent nearly identical procedures without MCAO. Results: The 10-mins MCAO induced no MR- or histologically detectable lesion in most of the mice and a small lesion in some of them. We thus had two groups with the same duration of ischemia but a different outcome, which could be compared to sham-operated mice and more severe ischemic mice (30-mins MCAO). Lactate increase, a hallmark of ischemic insult, was only detected significantly after 30-mins MCAO, whereas at 3 h post ischemia, glutamine was increased in all ischemic mice independently of duration and outcome. In contrast, glutamate, and even more so, N-acetyl-aspartate, decreased only in those mice exhibiting visible lesions on T2-weighted images at 24 h. Conclusions: These results suggest that an increased glutamine/glutamate ratio is a sensitive marker indicating the presence of an excitotoxic insult. Glutamate and NAA, on the other hand, appear to predict permanent neuronal damage. In conclusion, as early as 3 h post ischemia, it is possible to identify early metabolic markers manifesting the presence of a mild ischemic insult as well as the lesion outcome at 24 h.

Dynamic graphics of parametrically linked multivariate methods used in compositional data analysis

Many multivariate methods that are apparently distinct can be linked by introducing oneor more parameters in their definition. Methods that can be linked in this way arecorrespondence analysis, unweighted or weighted logratio analysis (the latter alsoknown as "spectral mapping"), nonsymmetric correspondence analysis, principalcomponent analysis (with and without logarithmic transformation of the data) andmultidimensional scaling. In this presentation I will show how several of thesemethods, which are frequently used in compositional data analysis, may be linkedthrough parametrizations such as power transformations, linear transformations andconvex linear combinations. Since the methods of interest here all lead to visual mapsof data, a "movie" can be made where where the linking parameter is allowed to vary insmall steps: the results are recalculated "frame by frame" and one can see the smoothchange from one method to another. Several of these "movies" will be shown, giving adeeper insight into the similarities and differences between these methods

Identification of reference genes for quantitative RT-PCR in ascending aortic aneurysms.

Hypertension and congenital aortic valve malformations are frequent causes of ascending aortic aneurysms. The molecular mechanisms of aneurysm formation under these circumstances are not well understood. Reference genes for gene activity studies in aortic tissue that are not influenced by aortic valve morphology and its hemodynamic consequences, aortic dilatation, hypertension, or antihypertensive medication are not available so far. This study determines genes in ascending aortic tissue that are independent of these parameters. Tissue specimens from dilated and undilated ascending aortas were obtained from 60 patients (age ≤70 years) with different morphologies of the aortic valve (tricuspid undilated n = 24, dilated n = 11; bicuspid undilated n = 6, dilated n = 15; unicuspid dilated n = 4). Of the studied individuals, 36 had hypertension, and 31 received ACE inhibitors or AT1 receptor antagonists. The specimens were obtained intraoperatively from the wall of the ascending aorta. We analyzed the expression levels of 32 candidate reference genes by quantitative RT-PCR (RT-qPCR). Differential expression levels were assessed by parametric statistics. The expression analysis of these 32 genes by RT-qPCR showed that EIF2B1, ELF1, and PPIA remained constant in their expression levels in the different specimen groups, thus being insensitive to aortic valve morphology, aortic dilatation, hypertension, and medication with ACE inhibitors or AT1 receptor antagonists. Unlike many other commonly used reference genes, the genes EIF2B1, ELF1, and PPIA are neither confounded by aortic comorbidities nor by antihypertensive medication and therefore are most suitable for gene expression analysis of ascending aortic tissue.

Kinematics and dynamic complexity of postural transitions in frail elderly subjects.

The aim of this study was to propose a methodology allowing a detailed characterization of body sit-to-stand/stand-to-sit postural transition. Parameters characterizing the kinematics of the trunk movement during sit-to-stand (Si-St) postural transition were calculated using one initial sensor system fixed on the trunk and a data logger. Dynamic complexity of these postural transitions was estimated by fractal dimension of acceleration-angular velocity plot. We concluded that this method provides a simple and accurate tool for monitoring frail elderly and to objectively evaluate the efficacy of a rehabilitation program.

Is there an efficient trap or collection method for sampling Anopheles darlingi and other malaria vectors that can describe the essential parameters affecting transmission dynamics as effectively as human landing catches? - A Review

Distribution, abundance, feeding behaviour, host preference, parity status and human-biting and infection rates are among the medical entomological parameters evaluated when determining the vector capacity of mosquito species. To evaluate these parameters, mosquitoes must be collected using an appropriate method. Malaria is primarily transmitted by anthropophilic and synanthropic anophelines. Thus, collection methods must result in the identification of the anthropophilic species and efficiently evaluate the parameters involved in malaria transmission dynamics. Consequently, human landing catches would be the most appropriate method if not for their inherent risk. The choice of alternative anopheline collection methods, such as traps, must consider their effectiveness in reproducing the efficiency of human attraction. Collection methods lure mosquitoes by using a mixture of olfactory, visual and thermal cues. Here, we reviewed, classified and compared the efficiency of anopheline collection methods, with an emphasis on Neotropical anthropophilic species, especially Anopheles darlingi, in distinct malaria epidemiological conditions in Brazil.

Fault diagnosis by refining the parameter uncertainty space of nonlinear dynamic systems

This paper deals with fault detection and isolation problems for nonlinear dynamic systems. Both problems are stated as constraint satisfaction problems (CSP) and solved using consistency techniques. The main contribution is the isolation method based on consistency techniques and uncertainty space refining of interval parameters. The major advantage of this method is that the isolation speed is fast even taking into account uncertainty in parameters, measurements, and model errors. Interval calculations bring independence from the assumption of monotony considered by several approaches for fault isolation which are based on observers. An application to a well known alcoholic fermentation process model is presented

Constraint satisfaction techniques under uncertain conditions for fault diagnosis in nonlinear dynamic systems

The speed of fault isolation is crucial for the design and reconfiguration of fault tolerant control (FTC). In this paper the fault isolation problem is stated as a constraint satisfaction problem (CSP) and solved using constraint propagation techniques. The proposed method is based on constraint satisfaction techniques and uncertainty space refining of interval parameters. In comparison with other approaches based on adaptive observers, the major advantage of the presented method is that the isolation speed is fast even taking into account uncertainty in parameters, measurements and model errors and without the monotonicity assumption. In order to illustrate the proposed approach, a case study of a nonlinear dynamic system is presented

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Developing methodologies to allow the robust prediction dynamic context dominated by perturbations

El principal objectiu del projecte era desenvolupar millores conceptuals i metodològiques que permetessin una millor predicció dels canvis en la distribució de les espècies (a una escala de paisatge) derivats de canvis ambientals en un context dominat per pertorbacions. En un primer estudi, vàrem comparar l'eficàcia de diferents models dinàmics per a predir la distribució de l'hortolà (Emberiza hortulana). Els nostres resultats indiquen que un model híbrid que combini canvis en la qualitat de l'hàbitat, derivats de canvis en el paisatge, amb un model poblacional espacialment explícit és una aproximació adequada per abordar canvis en la distribució d'espècies en contextos de dinàmica ambiental elevada i una capacitat de dispersió limitada de l'espècie objectiu. En un segon estudi abordarem la calibració mitjançant dades de seguiment de models de distribució dinàmics per a 12 espècies amb preferència per hàbitats oberts. Entre les conclusions extretes destaquem: (1) la necessitat de que les dades de seguiment abarquin aquelles àrees on es produeixen els canvis de qualitat; (2) el biaix que es produeix en la estimació dels paràmetres del model d'ocupació quan la hipòtesi de canvi de paisatge o el model de qualitat d'hàbitat són incorrectes. En el darrer treball estudiarem el possible impacte en 67 espècies d’ocells de diferents règims d’incendis, definits a partir de combinacions de nivells de canvi climàtic (portant a un augment esperat de la mida i freqüència d’incendis forestals), i eficiència d’extinció per part dels bombers. Segons els resultats dels nostres models, la combinació de factors antropogènics del regim d’incendis, tals com l’abandonament rural i l’extinció, poden ser més determinants per als canvis de distribució que els efectes derivats del canvi climàtic. Els productes generats inclouen tres publicacions científiques, una pàgina web amb resultats del projecte i una llibreria per a l'entorn estadístic R.

Assembling genes from predicted exons in linear time with dynamic programming

In a number of programs for gene structure prediction in higher eukaryotic genomic sequences, exon prediction is decoupled from gene assembly: a large pool of candidate exons is predicted and scored from features located in the query DNA sequence, and candidate genes are assembled from such a pool as sequences of nonoverlapping frame-compatible exons. Genes are scored as a function of the scores of the assembled exons, and the highest scoring candidate gene is assumed to be the most likely gene encoded by the query DNA sequence. Considering additive gene scoring functions, currently available algorithms to determine such a highest scoring candidate gene run in time proportional to the square of the number of predicted exons. Here, we present an algorithm whose running time grows only linearly with the size of the set of predicted exons. Polynomial algorithms rely on the fact that, while scanning the set of predicted exons, the highest scoring gene ending in a given exon can be obtained by appending the exon to the highest scoring among the highest scoring genes ending at each compatible preceding exon. The algorithm here relies on the simple fact that such highest scoring gene can be stored and updated. This requires scanning the set of predicted exons simultaneously by increasing acceptor and donor position. On the other hand, the algorithm described here does not assume an underlying gene structure model. Indeed, the definition of valid gene structures is externally defined in the so-called Gene Model. The Gene Model specifies simply which gene features are allowed immediately upstream which other gene features in valid gene structures. This allows for great flexibility in formulating the gene identification problem. In particular it allows for multiple-gene two-strand predictions and for considering gene features other than coding exons (such as promoter elements) in valid gene structures.

Chroma binary similarity and local alignment applied to cover song identification

We present a new technique for audio signal comparison based on tonal subsequence alignment and its application to detect cover versions (i.e., different performances of the same underlying musical piece). Cover song identification is a task whose popularity has increased in the Music Information Retrieval (MIR) community along in the past, as it provides a direct and objective way to evaluate music similarity algorithms.This article first presents a series of experiments carried outwith two state-of-the-art methods for cover song identification.We have studied several components of these (such as chroma resolution and similarity, transposition, beat tracking or Dynamic Time Warping constraints), in order to discover which characteristics would be desirable for a competitive cover song identifier. After analyzing many cross-validated results, the importance of these characteristics is discussed, and the best-performing ones are finally applied to the newly proposed method. Multipleevaluations of this one confirm a large increase in identificationaccuracy when comparing it with alternative state-of-the-artapproaches.

Influence of a lifestyle intervention in preschool children on physiological and psychological parameters (Ballabeina): study design of a cluster randomized controlled trial.

Childhood obesity and physical inactivity are increasing dramatically worldwide. Children of low socioeconomic status and/or children of migrant background are especially at risk. In general, the overall effectiveness of school-based programs on health-related outcomes has been disappointing. A special gap exists for younger children and in high risk groups. This paper describes the rationale, design, curriculum, and evaluation of a multicenter preschool randomized intervention study conducted in areas with a high migrant population in two out of 26 Swiss cantons. Twenty preschool classes in the German (canton St. Gallen) and another 20 in the French (canton Vaud) part of Switzerland were separately selected and randomized to an intervention and a control arm by the use of opaque envelopes. The multidisciplinary lifestyle intervention aimed to increase physical activity and sleep duration, to reinforce healthy nutrition and eating behaviour, and to reduce media use. According to the ecological model, it included children, their parents and the teachers. The regular teachers performed the majority of the intervention and were supported by a local health promoter. The intervention included physical activity lessons, adaptation of the built infrastructure; promotion of regional extracurricular physical activity; playful lessons about nutrition, media use and sleep, funny homework cards and information materials for teachers and parents. It lasted one school year. Baseline and post-intervention evaluations were performed in both arms. Primary outcome measures included BMI and aerobic fitness (20 m shuttle run test). Secondary outcomes included total (skinfolds, bioelectrical impedance) and central (waist circumference) body fat, motor abilities (obstacle course, static and dynamic balance), physical activity and sleep duration (accelerometry and questionnaires), nutritional behaviour and food intake, media use, quality of life and signs of hyperactivity (questionnaires), attention and spatial working memory ability (two validated tests). Researchers were blinded to group allocation. The purpose of this paper is to outline the design of a school-based multicenter cluster randomized, controlled trial aiming to reduce body mass index and to increase aerobic fitness in preschool children in culturally different parts of Switzerland with a high migrant population. Trial Registration: (clinicaltrials.gov) NCT00674544.

Back to square one: Identification issues in DSGE models

We investigate identifiability issues in DSGE models and their consequences for parameter estimation and model evaluation when the objective function measures the distance between estimated and model impulse responses. We show that observational equivalence, partial and weak identification problems are widespread, that they lead to biased estimates, unreliable t-statistics and may induce investigators to select false models. We examine whether different objective functions affect identification and study how small samples interact with parameters and shock identification. We provide diagnostics and tests to detect identification failures and apply them to a state-of-the-art model.

Dynamic graphics of parametrically linked multivariate methods used in compositional data analysis

Many multivariate methods that are apparently distinct can be linked by introducing oneor more parameters in their definition. Methods that can be linked in this way arecorrespondence analysis, unweighted or weighted logratio analysis (the latter alsoknown as "spectral mapping"), nonsymmetric correspondence analysis, principalcomponent analysis (with and without logarithmic transformation of the data) andmultidimensional scaling. In this presentation I will show how several of thesemethods, which are frequently used in compositional data analysis, may be linkedthrough parametrizations such as power transformations, linear transformations andconvex linear combinations. Since the methods of interest here all lead to visual mapsof data, a "movie" can be made where where the linking parameter is allowed to vary insmall steps: the results are recalculated "frame by frame" and one can see the smoothchange from one method to another. Several of these "movies" will be shown, giving adeeper insight into the similarities and differences between these methods.