SIMULIT : un modèle de simulation pour l'analyse et la planification de l'activité hospitalière

Le modèle développé à l'Institut universitaire de médecine sociale et préventive de Lausanne utilise un programme informatique pour simuler les mouvements d'entrées et de sorties des hôpitaux de soins généraux. Cette simulation se fonde sur les données récoltées de routine dans les hôpitaux; elle tient notamment compte de certaines variations journalières et saisonnières, du nombre d'entrées, ainsi que du "Case-Mix" de l'hôpital, c'est-à-dire de la répartition des cas selon les groupes cliniques et l'âge des patients.

Intravenous streptomycin dosing regimen in a patient undergoing hemodialysis: Plasma level monitoring and pharmacokinetic simulation

Introduction: Streptomycin, as other aminoglycosides, exhibits concentration-dependent bacterial killing but has a narrow therapeutic window. It is primarily eliminated unchanged by the kidneys. Data and dosing information to achieve a safe regimen in patients with chronic renal failure undergoing hemodialysis (HD) are scarce. Although main adverse reactions are related to prolonged, elevated serum concentrations, literature recommendation is to administer streptomycin after each HD. Patients (or Materials) and Methods: We report the case of a patient with end-stage renal failure, undergoing HD, who was successfully treated with streptomycin for gentamicin-resistant Enterococcus faecalis bacteremia with prosthetic arteriovenous fistula infection. Streptomycin was administered intravenously 7.5 mg/kg, 3 hours before each dialysis (3 times a week) during 6 weeks in combination with amoxicillin. Streptomycin plasma levels were monitored with repeated blood sampling before, after, and between HD sessions. A 2-compartment model was used to reconstruct the concentration time profile over days on and off HD. Results: Streptomycin trough plasma-concentration was 2.8 mg/L. It peaked to 21.4 mg/L 30 minutes after intravenous administration, decreased to 18.2 mg/L immediately before HD, and dropped to 4.5 mg/L at the end of a 4-hour HD session. Plasma level increased again to 5.7 mg/L 2 hours after the end of HD and was 2.8 mg/L 48 hours later, before the next administration and HD. The pharmacokinetics of streptomycin was best described with a 2-compartment model. The computer simulation fitted fairly well to the observed concentrations during or between HD sessions. Redistribution between the 2 compartments after the end of HD reproduced the rebound of plasma concentrations after HD. No significant toxicity was observed during treatment. The outcome of the infection was favorable, and no sign of relapse was observed after a follow-up of 3 months. Conclusion: Streptomycin administration of 7.5 mg/kg 3 hours before HD sessions in a patient with end-stage renal failure resulted in an effective and safe dosing regimen. Monitoring plasma levels along with pharmacokinetic simulation document the suitability of this dosing scheme, which should replace current dosage recommendations for streptomycin in HD.

Total shoulder arthroplasty: downward inclination of the glenoid component to balance supraspinatus deficiency.

HYPOTHESIS: Supraspinatus deficiency associated with total shoulder arthroplasty (TSA) provokes eccentric loading and may induce loosening of the glenoid component. A downward inclination of the glenoid component has been proposed to balance supraspinatus deficiency. METHODS: This hypothesis was assessed by a numeric musculoskeletal model of the glenohumeral joint during active abduction. Three cases were compared: TSA with normal muscular function, TSA with supraspinatus deficiency, and TSA with supraspinatus deficiency and downward inclination of the glenoid. RESULTS: Supraspinatus deficiency increased humeral migration and eccentric loading. A downward inclination of the glenoid partly balanced the loss of stability, but this potential advantage was counterbalanced by an important stress increase within the glenoid cement. The additional subchondral bone reaming required to incline the glenoid component indeed reduced the bone support, increasing cement deformation and stress. CONCLUSION: Glenoid inclination should not be obtained at the expense of subchondral bone support.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

The effects of physiologically plausible connectivity structure on local and global dynamics in large scale brain models.

Functionally relevant large scale brain dynamics operates within the framework imposed by anatomical connectivity and time delays due to finite transmission speeds. To gain insight on the reliability and comparability of large scale brain network simulations, we investigate the effects of variations in the anatomical connectivity. Two different sets of detailed global connectivity structures are explored, the first extracted from the CoCoMac database and rescaled to the spatial extent of the human brain, the second derived from white-matter tractography applied to diffusion spectrum imaging (DSI) for a human subject. We use the combination of graph theoretical measures of the connection matrices and numerical simulations to explicate the importance of both connectivity strength and delays in shaping dynamic behaviour. Our results demonstrate that the brain dynamics derived from the CoCoMac database are more complex and biologically more realistic than the one based on the DSI database. We propose that the reason for this difference is the absence of directed weights in the DSI connectivity matrix.

The activation process of the alpha1B-adrenergic receptor: potential role of protonation and hydrophobicity of a highly conserved aspartate.

In this study, a quantitative approach was used to investigate the role of D142, which belongs to the highly conserved E/DRY sequence, in the activation process of the alpha1B-adrenergic receptor (alpha1B-AR). Experimental and computer-simulated mutagenesis were performed by substituting all possible natural amino acids at the D142 site. The resulting congeneric set of proteins together with the finding that all the receptor mutants show various levels of constitutive (agonist-independent) activity enabled us to quantitatively analyze the relationships between structural/dynamic features and the extent of constitutive activity. Our results suggest that the hydrophobic/hydrophilic character of D142, which could be regulated by protonation/deprotonation of this residue, is an important modulator of the transition between the inactive (R) and active (R*) state of the alpha1B-AR. Our study represents an example of quantitative structure-activity relationship analysis of the activation process of a G protein-coupled receptor.

Simulations in evolution. II. Relative fitness and the propagation of mutants.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous article presented a histogram model [1] consisting in populations of individuals whose number changed under the influence of variation and/or fitness, the total population remaining constant. Individuals are classified into bins, and the content of each bin is calculated generation after generation by an Excel spreadsheet. Here, we apply the histogram model to a stable population with fitness F(1)=1.00 in which one or two fitter mutants emerge. In a first scenario, a single mutant emerged in the population whose fitness was greater than 1.00. The simulations ended when the original population was reduced to a single individual. The histogram model was validated by excellent agreement between its predictions and those of a classical continuous function (Eqn. 1) which predicts the number of generations needed for a favorable mutation to spread throughout a population. But in contrast to Eqn. 1, our histogram model is adaptable to more complex scenarios, as demonstrated here. In the second and third scenarios, the original population was present at time zero together with two mutants which differed from the original population by two higher and distinct fitness values. In the fourth scenario, the large original population was present at time zero together with one fitter mutant. After a number of generations, when the mutant offspring had multiplied, a second mutant was introduced whose fitness was even greater. The histogram model also allows Shannon entropy (SE) to be monitored continuously as the information content of the total population decreases or increases. The results of these simulations illustrate, in a graphically didactic manner, the influence of natural selection, operating through relative fitness, in the emergence and dominance of a fitter mutant.

Postoperative Cognitive Dysfunction (POCD) and Inflammatory Markers in Elderly Patients

Introduction: Particularly in elderly patients, the brain responds to a systemic inflammatory response with an increased production of inflammatory mediators. This has hypothetically been linked to the development of postoperative cognitive dysfunction (POCD). Methods: We investigated 31 patients aged >65 yrs undergoing elective major surgery under standardized general anaesthesia (thiopental, sevoflurane, fentanyl, atracurium). Cognitive function was measured preoperatively and 7 days postoperatively using the extended version of the Consortium to Establish a Registry for Alzheimer's Disease - Neuropsychological Assessment Battery (CERAD-NAB, validated German version) for which we developed a diagnostic cut-off in healthy elderly volunteers. Systemic C-reactive protein (CRP) and interleukin 6 (IL-6) were measured preoperatively, 2 days postoperatively, and 7 days postoperatively. Values for CRP, IL-6, operative characteristics and hospital length of stay in patients with POCD and without POCD were compared using the Mann- Whitney U test and are shown as median [range]. Results: Fourteen patients (45%) developed POCD. Values for CRP were not statistically different in patients with POCD and without POCD but tended to be higher in patients with POCD 2 days postoperatively. Patients with POCD had significantly higher IL-6 values on postoperative days 2 and 7 (table 1). These patients also had a significantly longer duration of anaesthesia (305 [195-620] vs.190 [150-560] min, p = 0.034), larger intraoperative blood loss (425 [0-1600] vs. 100 [0-1500] ml, p = 0.018) and longer hospital stays (15 [8-45] vs. 8 [4-40] days, p = 0.008). Table 1 POCD (n = 14) No POCD (n = 17) p value CRP (mg/dl) preop. 4.0 [1.0-245] 4.2 [0.3-36.2] 0.6 2 days postop. 223 [20-318] 98 [4.5-384] 0.07 7 days postop. 58 [15-147] 44 [11-148] 0.2 IL-6 (U/ml) preop. 2[2-28.1] 2 [2-7.3] 0.8 2 days postop. 56 [17-315] 20 [2-123] 0.009 7 days postop. 9[2-77] 4 [2-16] 0.03 Interpretation: In this small group of patients, high IL-6 values postoperatively were associated with POCD supporting a role for systemic inflammation in the development of POCD. In patients with POCD, duration of anaesthesia was significantly longer, and intraoperative blood losses were larger. These risk factors will need to be confirmed in a larger group of patients. The difference in length of stay may be indicative of postoperative complications, which have been linked to POCD earlier.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

Analysis of the role of mobility-lifetime products in the performance of amorphous silicon p-i-n solar cells

An analytical model of an amorphous silicon p-i-n solar cell is presented to describe its photovoltaic behavior under short-circuit conditions. It has been developed from the analysis of numerical simulation results. These results reproduce the experimental illumination dependence of short-circuit resistance, which is the reciprocal slope of the I(V) curve at the short-circuit point. The recombination rate profiles show that recombination in the regions of charged defects near the p-i and i-n interfaces should not be overlooked. Based on the interpretation of the numerical solutions, we deduce analytical expressions for the recombination current and short-circuit resistance. These expressions are given as a function of an effective ¿¿ product, which depends on the intensity of illumination. We also study the effect of surface recombination with simple expressions that describe its influence on current loss and short-circuit resistance.

Algebraic decay of velocity fluctuations in a confined fluid

Computer simulations of a colloidal particle suspended in a fluid confined by rigid walls show that, at long times, the velocity correlation function decays with a negative algebraic tail. The exponent depends on the confining geometry, rather than the spatial dimensionality. We can account for the tail by using a simple mode-coupling theory which exploits the fact that the sound wave generated by a moving particle becomes diffusive.

The switch from conventional to atypical antipsychotic treatment should not be based exclusively on the presence of cognitive deficits. A pilot study in individuals with schizophrenia.

Background: Atypical antipsychotics provide better control of the negative and affective symptoms of schizophrenia when compared with conventional neuroleptics; nevertheless, their heightened ability to improve cognitive dysfunction remains a matter of debate. This study aimed to examine the changes in cognition associated with long-term antipsychotic treatment and to evaluate the effect of the type of antipsychotic (conventional versus novel antipsychotic drugs) on cognitive performance over time. Methods: In this naturalistic study, we used a comprehensive neuropsychological battery of tests to assess a sample of schizophrenia patients taking either conventional (n = 13) or novel antipsychotics (n = 26) at baseline and at two years after. Results: Continuous antipsychotic treatment regardless of class was associated with improvement on verbal fluency, executive functions, and visual and verbal memory. Patients taking atypical antipsychotics did not show greater cognitive enhancement over two years than patients taking conventional antipsychotics. Conclusions Although long-term antipsychotic treatment slightly improved cognitive function, the switch from conventional to atypical antipsychotic treatment should not be based exclusively on the presence of these cognitive deficits.

Dosage optimization of treatments using population pharmacokinetic modeling and simulation.

Pharmacokinetic variability in drug levels represent for some drugs a major determinant of treatment success, since sub-therapeutic concentrations might lead to toxic reactions, treatment discontinuation or inefficacy. This is true for most antiretroviral drugs, which exhibit high inter-patient variability in their pharmacokinetics that has been partially explained by some genetic and non-genetic factors. The population pharmacokinetic approach represents a very useful tool for the description of the dose-concentration relationship, the quantification of variability in the target population of patients and the identification of influencing factors. It can thus be used to make predictions and dosage adjustment optimization based on Bayesian therapeutic drug monitoring (TDM). This approach has been used to characterize the pharmacokinetics of nevirapine (NVP) in 137 HIV-positive patients followed within the frame of a TDM program. Among tested covariates, body weight, co-administration of a cytochrome (CYP) 3A4 inducer or boosted atazanavir as well as elevated aspartate transaminases showed an effect on NVP elimination. In addition, genetic polymorphism in the CYP2B6 was associated with reduced NVP clearance. Altogether, these factors could explain 26% in NVP variability. Model-based simulations were used to compare the adequacy of different dosage regimens in relation to the therapeutic target associated with treatment efficacy. In conclusion, the population approach is very useful to characterize the pharmacokinetic profile of drugs in a population of interest. The quantification and the identification of the sources of variability is a rational approach to making optimal dosage decision for certain drugs administered chronically.

Spatial behavior in groups: an agent-based approach

We present an agent-based model with the aim of studying how macro-level dynamics of spatial distances among interacting individuals in a closed space emerge from micro-level dyadic and local interactions. Our agents moved on a lattice (referred to as a room) using a model implemented in a computer program called P-Space in order to minimize their dissatisfaction, defined as a function of the discrepancy between the real distance and the ideal, or desired, distance between agents. Ideal distances evolved in accordance with the agent's personal and social space, which changed throughout the dynamics of the interactions among the agents. In the first set of simulations we studied the effects of the parameters of the function that generated ideal distances, and in a second set we explored how group macrolevel behavior depended on model parameters and other variables. We learned that certain parameter values yielded consistent patterns in the agents' personal and social spaces, which in turn led to avoidance and approaching behaviors in the agents. We also found that the spatial behavior of the group of agents as a whole was influenced by the values of the model parameters, as well as by other variables such as the number of agents. Our work demonstrates that the bottom-up approach is a useful way of explaining macro-level spatial behavior. The proposed model is also shown to be a powerful tool for simulating the spatial behavior of groups of interacting individuals.

Comparative functional analysis of two fibroblast growth factor receptor 1 (FGFR1) mutations affecting the same residue (R254W and R254Q) in isolated hypogonadotropic hypogonadism (IHH).

FGFR1 mutations have been identified in both Kallmann syndrome and normosmic HH (nIHH). To date, few mutations in the FGFR1 gene have been structurally or functionally characterized in vitro to identify molecular mechanisms that contribute to the disease pathogenesis. We attempted to define the in vitro functionality of two FGFR1 mutants (R254W and R254Q), resulting from two different amino acid substitutions of the same residue, and to correlate the in vitro findings to the patient phenotypes. Two unrelated GnRH deficient probands were found to harbor mutations in FGFR1 (R254W and R254Q). Mutant signaling activity and expression levels were evaluated in vitro and compared to a wild type (WT) receptor. Signaling activity was determined by a FGF2/FGFR1 dependent transcription reporter assay. Receptor total expression levels were assessed by Western blot and cell surface expression was measured by a radiolabeled antibody binding assay. The R254W maximal receptor signaling capacity was reduced by 45% (p<0.01) while R254Q activity was not different from WT. However, both mutants displayed diminished total protein expression levels (40 and 30% reduction relative to WT, respectively), while protein maturation was unaffected. Accordingly, cell surface expression levels of the mutant receptors were also significantly reduced (35% p<0.01 and 15% p<0.05, respectively). The p.R254W and p.R254Q are both loss-of-function mutations as demonstrated by their reduced overall and cell surface expression levels suggesting a deleterious effect on receptor folding and stability. It appears that a tryptophan substitution at R254 is more disruptive to receptor structure than the more conserved glutamine substitution. No clear correlation between the severity of in vitro loss-of-function and phenotypic presentation could be assigned.