Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Personalized drug administrations using Support Vector Machine

The decision-making process regarding drug dose, regularly used in everyday medical practice, is critical to patients' health and recovery. It is a challenging process, especially for a drug with narrow therapeutic ranges, in which a medical doctor decides the quantity (dose amount) and frequency (dose interval) on the basis of a set of available patient features and doctor's clinical experience (a priori adaptation). Computer support in drug dose administration makes the prescription procedure faster, more accurate, objective, and less expensive, with a tendency to reduce the number of invasive procedures. This paper presents an advanced integrated Drug Administration Decision Support System (DADSS) to help clinicians/patients with the dose computing. Based on a support vector machine (SVM) algorithm, enhanced with the random sample consensus technique, this system is able to predict the drug concentration values and computes the ideal dose amount and dose interval for a new patient. With an extension to combine the SVM method and the explicit analytical model, the advanced integrated DADSS system is able to compute drug concentration-to-time curves for a patient under different conditions. A feedback loop is enabled to update the curve with a new measured concentration value to make it more personalized (a posteriori adaptation).

Risk factors of postictal generalized EEG suppression in generalized convulsive seizures.

OBJECTIVE: To identify the clinical determinants of occurrence of postictal generalized EEG suppression (PGES) after generalized convulsive seizures (GCS). METHODS: We reviewed the video-EEG recordings of 417 patients included in the REPO2MSE study, a multicenter prospective cohort study of patients with drug-resistant focal epilepsy. According to ictal semiology, we classified GCS into 3 types: tonic-clonic GCS with bilateral and symmetric tonic arm extension (type 1), clonic GCS without tonic arm extension or flexion (type 2), and GCS with unilateral or asymmetric tonic arm extension or flexion (type 3). Association between PGES and person-specific or seizure-specific variables was analyzed after correction for individual effects and the varying number of seizures. RESULTS: A total of 99 GCS in 69 patients were included. Occurrence of PGES was independently associated with GCS type (p < 0.001) and lack of early administration of oxygen (p < 0.001). Odds ratio (OR) for GCS type 1 in comparison with GCS type 2 was 66.0 (95% confidence interval [CI 5.4-801.6]). In GCS type 1, risk of PGES was significantly increased when the seizure occurred during sleep (OR 5.0, 95% CI 1.2-20.9) and when oxygen was not administered early (OR 13.4, 95% CI 3.2-55.9). CONCLUSION: The risk of PGES dramatically varied as a function of GCS semiologic characteristics. Whatever the type of GCS, occurrence of PGES was prevented by early administration of oxygen.

Expression of Staphylococcus aureus clumping factor A in Lactococcus lactis subsp. cremoris using a new shuttle vector.

Staphylococcus aureus harbors redundant adhesins mediating tissue colonization and infection. To evaluate their intrinsic role outside of the staphylococcal background, a system was designed to express them in Lactococcus lactis subsp. cremoris 1363. This bacterium is devoid of virulence factors and has a known genetic background. A new Escherichia coli-L. lactis shuttle and expression vector was constructed for this purpose. First, the high-copy-number lactococcal plasmid pIL253 was equipped with the oriColE1 origin, generating pOri253 that could replicate in E. coli. Second, the lactococcal promoters P23 or P59 were inserted at one end of the pOri253 multicloning site. Gene expression was assessed by a luciferase reporter system. The plasmid carrying P23 (named pOri23) expressed luciferase constitutively at a level 10,000 times greater than did the P59-containing plasmid. Transcription was absent in E. coli. The staphylococcal clumping factor A (clfA) gene was cloned into pOri23 and used as a model system. Lactococci carrying pOri23-clfA produced an unaltered and functional 130-kDa ClfA protein attached to their cell walls. This was indicated both by the presence of the protein in Western blots of solubilized cell walls and by the ability of ClfA-positive lactococci to clump in the presence of plasma. ClfA-positive lactococci had clumping titers (titer of 4,112) similar to those of S. aureus Newman in soluble fibrinogen and bound equally well to solid-phase fibrinogen. These experiments provide a new way to study individual staphylococcal pathogenic factors and might complement both classical knockout mutagenesis and modern in vivo expression technology and signature tag mutagenesis.

Multi-scale support vector algorithms for hot spot detection and modelling

The algorithmic approach to data modelling has developed rapidly these last years, in particular methods based on data mining and machine learning have been used in a growing number of applications. These methods follow a data-driven methodology, aiming at providing the best possible generalization and predictive abilities instead of concentrating on the properties of the data model. One of the most successful groups of such methods is known as Support Vector algorithms. Following the fruitful developments in applying Support Vector algorithms to spatial data, this paper introduces a new extension of the traditional support vector regression (SVR) algorithm. This extension allows for the simultaneous modelling of environmental data at several spatial scales. The joint influence of environmental processes presenting different patterns at different scales is here learned automatically from data, providing the optimum mixture of short and large-scale models. The method is adaptive to the spatial scale of the data. With this advantage, it can provide efficient means to model local anomalies that may typically arise in situations at an early phase of an environmental emergency. However, the proposed approach still requires some prior knowledge on the possible existence of such short-scale patterns. This is a possible limitation of the method for its implementation in early warning systems. The purpose of this paper is to present the multi-scale SVR model and to illustrate its use with an application to the mapping of Cs137 activity given the measurements taken in the region of Briansk following the Chernobyl accident.

A drug administration decision support system

Drug delivery is one of the most common clinical routines in hospitals, and is critical to patients' health and recovery. It includes a decision making process in which a medical doctor decides the amount (dose) and frequency (dose interval) on the basis of a set of available patients' feature data and the doctor's clinical experience (a priori adaptation). This process can be computerized in order to make the prescription procedure in a fast, objective, inexpensive, non-invasive and accurate way. This paper proposes a Drug Administration Decision Support System (DADSS) to help clinicians/patients with the initial dose computing. The system is based on a Support Vector Machine (SVM) algorithm for estimation of the potential drug concentration in the blood of a patient, from which a best combination of dose and dose interval is selected at the level of a DSS. The addition of the RANdom SAmple Consensus (RANSAC) technique enhances the prediction accuracy by selecting inliers for SVM modeling. Experiments are performed for the drug imatinib case study which shows more than 40% improvement in the prediction accuracy compared with previous works. An important extension to the patient features' data is also proposed in this paper.

Prevalence and determinants of QT interval prolongation in medical inpatients.

BACKGROUND: QT interval prolongation carries an increased risk of torsade de pointes and death. AIM: We sought to determine the prevalence of QT prolongation in medical inpatients and to identify determinants of this condition. METHODS: We enrolled consecutive patients who were admitted to the internal medicine ward and who had an electrocardiogram performed within 24 h of admission. We collected information on baseline patient characteristics and the use of QT-prolonging drugs. Two blinded readers manually measured the QT intervals. QT intervals were corrected for heart rate using the traditional Bazett formula and the linear regression-based Framingham formula. We used logistic regression to identify patient characteristics and drugs that were independently associated with QTc prolongation. RESULTS: Of 537 inpatients, 22.3% had a prolonged QTc based on the Bazett formula. The adjusted odds for QTc prolongation based on the Bazett correction were significantly higher in patients who had liver disease (OR 2.9, 95% CI: 1.5-5.6), hypokalaemia (OR 3.3, 95% CI: 1.9-5.6) and who were taking ≥1 QT-prolonging drug at admission (OR 1.7, 95% CI: 1.1-2.6). Overall, 50.8% of patients with QTc prolongation received additional QT-prolonging drugs during hospitalisation. CONCLUSIONS: The prevalence of QTc prolongation was high among medical inpatients but depended on the method used to correct for heart rate. The use of QT-prolonging drugs, hypokalaemia and liver disease increased the risk of QTc prolongation. Many patients with QTc prolongation received additional QT-prolonging drugs during hospitalisation, further increasing the risk of torsade de pointes and death.

Highly inducible synthesis of heterologous proteins in epithelial cells carrying a glucocorticoid-responsive vector.

A glucocorticoid-responsive vector is described which allows for the highly inducible expression of complementary DNAs (cDNAs) in stably transfected mammalian cell lines. This vector, pLK-neo, composed of a variant mouse mammary tumor virus long terminal repeat promoter, containing a hormone regulatory element, a Geneticin resistance-encoding gene in a simian virus 40 transcription unit, and a polylinker insertion site for heterologous cDNAs, was used to express the polymeric immunoglobulin (poly-Ig) receptor and the thymocyte marker, Thy-1, in Madin-Darby canine kidney (MDCK) cells and in murine fibroblast L cells. A high level of poly-Ig receptor or Thy-1 mRNA accumulation was observed in MDCK cells in response to dexamethasone with a parallel ten- to 200-fold increase in protein synthesis depending on the recombinant protein and the transfected cell clone.

Birth records from Swiss married couples analyzed over of the past 35 years reveal an aging of first-time mothers by 5.1 years while the interpregnancy interval has shortened

RésuméIntroduction : Le travail de thèse est un article publié dans le journal « Fertility & Sterility » s'intitulant « Birth records from Swiss married couples analyzed over the past 35 years reveal an aging of first-time mothers by 5.1 years while the interpregnancy interval has shortened ».Méthodes : Les données concernant l'âge auquel les femmes mariées donnaient naissance ont été obtenues de l'Office fédéral de la statistique. Nous avons examiné la période allant de 1969 à 2006. Cet intervalle de temps choisi nous a permis de prendre en considération un total de 2'716'370 naissances. L'âge moyen des parturientes à la naissance de leur 1er, 2ème et 3ème enfant a été calculé pour chaque année à l'aide du logiciel Excel. Grâce à ces données on a pu obtenir les intervalles inter-gestationnels théoriques entre le 1er et 2ème enfant, ainsi qu'entre le 2ème et 3ème enfant.Résultats : Nous pouvons constater que l'intervalle inter-gestationnel théorique entre la première et la deuxième grossesse était de 23.2 mois en 1969 et passait à 13.0 mois en 2006. L'intervalle compris entre la deuxième et la troisième grossesse passait de 22.4 mois en 1969 à 7.9 mois en 2006. Notre analyse suggère donc que c'est le facteur social qui exerce un effet plus important que le facteur biologique sur les intervalles inter-gestationnels, car ces intervalles diminuaient malgré l'augmentation de l'âge des mères à la naissance.

Substitution of (R,S)-methadone by (R)-methadone: Impact on QTc interval.

Methadone is administered as a chiral mixture of (R,S)-methadone. The opioid effect is mainly mediated by (R)-methadone, whereas (S)-methadone blocks the human ether-à-go-go-related gene (hERG) voltage-gated potassium channel more potently, which can cause drug-induced long QT syndrome, leading to potentially lethal ventricular tachyarrhythmias. To investigate whether substitution of (R,S)-methadone by (R)-methadone could reduce the corrected QT (QTc) interval, (R,S)-methadone was replaced by (R)-methadone (half-dose) in 39 opioid-dependent patients receiving maintenance treatment for 14 days. (R)-methadone was then replaced by the initial dose of (R,S)-methadone for 14 days (n = 29). Trough (R)-methadone and (S)-methadone plasma levels and electrocardiogram measurements were taken. The Fridericia-corrected QT (QTcF) interval decreased when (R,S)-methadone was replaced by a half-dose of (R)-methadone; the median (interquartile range [IQR]) values were 423 (398-440) milliseconds (ms) and 412 (395-431) ms (P = .06) at days 0 and 14, respectively. Using a univariate mixed-effect linear model, the QTcF value decreased by a mean of -3.9 ms (95% confidence interval [CI], -7.7 to -0.2) per week (P = .04). The QTcF value increased when (R)-methadone was replaced by the initial dose of (R,S)-methadone for 14 days; median (IQR) values were 424 (398-436) ms and 424 (412-443) ms (P = .01) at days 14 and 28, respectively. The univariate model showed that the QTcF value increased by a mean of 4.7 ms (95% CI, 1.3-8.1) per week (P = .006). Substitution of (R,S)-methadone by (R)-methadone reduces the QTc interval value. A safer cardiac profile of (R)-methadone is in agreement with previous in vitro and pharmacogenetic studies. If the present results are confirmed by larger studies, (R)-methadone should be prescribed instead of (R,S)-methadone to reduce the risk of cardiac toxic effects and sudden death.