Realistic Single Cell Modeling – from Experiment to Simulation

This tutorial gives a step by step explanation of how one uses experimental data to construct a biologically realistic multicompartmental model. Special emphasis is given on the many ways that this process can be imprecise. The tutorial is intended for both experimentalists who want to get into computer modeling and for computer scientists who use abstract neural network models but are curious about biological realistic modeling. The tutorial is not dependent on the use of a specific simulation engine, but rather covers the kind of data needed for constructing a model, how they are used, and potential pitfalls in the process.

Color and Scent: How Single Genes Influence Pollinator Attraction

A major function of angiosperm flowers is the recruitment of animal pollinators that serve to transfer pollen among conspecific plants. Distinct sets of floral characteristics, called pollination syndromes, are correlated with visitation by specific groups of pollinators. Switches among pollination syndromes have occurred in many plant families. Such switches must have involved coordinated changes in multiple traits and multiple genes. Two well-studied floral traits affecting pollinator attraction are petal color and scent production. We review current knowledge about the biosynthetic pathways for floral color and scent production and their interaction at the genetic and biochemical levels. A key question in the field concerns the genes that underlie natural variation in color and scent and how such genes affect pollinator preference, reproductive isolation, and ultimately speciation.

Two-dimensional assembly and local redox-activity of molecular hybrid structures in an electrochemical environment

The self-assembly and redox-properties of two viologen derivatives, N-hexyl-N-(6-thiohexyl)-4,4-bipyridinium bromide (HS-6V6-H) and N,N-bis(6-thiohexyl)-4,4-bipyridinium bromide (HS-6V6-SH), immobilized on Au(111)-(1x1) macro-electrodes were investigated by cyclic voltammetry, surface enhanced infrared spectroscopy (SEIRAS) and in situ scanning tunneling microscopy (STM). Depending on the assembly conditions one could distinguish three different types of adlayers for both viologens: a low coverage disordered and an ordered striped phase of flat oriented molecules as well as a high coverage monolayer composed of tilted viologen moieties. Both molecules, HS-6V6-H and HS-6V6-SH, were successfully immobilized on Au(poly) nano-electrodes, which gave a well-defined redox-response in the lower pA–current range. An in situ STM configuration was employed to explore electron transport properties of single molecule junctions Au(T)|HS-6V6-SH(HS-6V6-H)|Au(S). The observed sigmoidal potential dependence, measured at variable substrate potential ES and at constant bias voltage (ET–ES), was attributed to electronic structure changes of the viologen moiety during the one-electron reduction/re-oxidation process V2+ V+. Tunneling experiments in asymmetric, STM-based junctions Au(T)-S-6V6-H|Au(S) revealed current (iT)–voltage (ET) curves with a maximum located at the equilibrium potential of the redox-process V2+ V+. The experimental iT–ET characteristics of the HS-6V6-H–modified tunneling junction were tentatively attributed to a sequential two-step electron transfer mechanism.

Electrochemical characterization of self-assembled Ferrocene-terminated alkanethiol monolayers on low-index gold single crystal electrodes

We present a voltammetric and in situ STM study of 11-ferrocenyl-1-undecanethiol (FcC11) assembled on low-index single crystal and polycrystalline gold electrodes. The influence of electrode surface structure as well as of structure defects in the self-assembled FcC11 monolayers on the electrochemical response during the oxidation and reduction of the terminal ferrocene group is explored. The nature of the redox peaks is discussed in detail. We identified the coexistence of disordered FcC11 regions with 2D patches of “locally ordered” FcC11 species. We demonstrate that close-packed domains are preferentially formed at atomically flat terraces. Increasing the defect density of the substrate surface leads to a decreasing amount of locally ordered FcC11 molecules.

Reduction in total plasma ghrelin levels following catecholamine depletion: relation to bulimic and depressive symptoms

There is increasing preclinical and clinical evidence of the important role played by the gastric peptide hormone ghrelin in the pathogenesis of symptoms of depression and eating disorders. To investigate the role of ghrelin and its considered counterpart, peptide tyrosine tyrosine (PYY), in the development of bulimic and depressive symptoms induced by catecholamine depletion, we administered the tyrosine hydroxylase inhibitor alpha-methyl-paratyrosine (AMPT) in a randomized, double-blind, placebo-controlled crossover, single-site experimental trial to 29 healthy controls and 20 subjects with fully recovered bulimia nervosa (rBN). We found a decrease between peprandial and postprandial plasma ghrelin levels (p < 0.0001) and a postprandial rise in plasma PYY levels (p < 0.0001) in both conditions in the entire study population. Plasma ghrelin levels decreased in the entire study population after treatment with AMPT compared to placebo (p < 0.006). AMPT-induced changes in plasma ghrelin levels were negatively correlated with AMPT-induced depressive symptoms (p < 0.004). Plasma ghrelin and plasma PYY levels were also negatively correlated (p < 0.05). We did not observe a difference in ghrelin or PYY response to catecholamine depletion between rBN subjects and healthy controls, and there was no correlation between plasma ghrelin and PYY levels and bulimic symptoms induced by catecholamine depletion. These findings suggest a relationship between catecholamines and ghrelin with depressive symptoms.

In situ dehydration behavior of zeolite-like pentagonite: A single-crystal X-ray study

The structural modifications upon heating of pentagonite, Ca(VO)(Si4O10)·4H2O (space group Ccm21, a=10.3708(2), b=14.0643(2), c=8.97810(10) Å, V=1309.53(3) Å3) were investigated by in situ temperature dependent single-crystal X-ray structure refinements. Diffraction data of a sample from Poona district (India) have been measured in steps of 25 up to 250 °C and in steps of 50 °C between 250 and 400 °C. Pentagonite has a porous framework structure made up by layers of silicate tetrahedra connected by V4+O5 square pyramids. Ca and H2O molecules are extraframework occupants. Room temperature diffraction data allowed refinement of H positions. The hydrogen-bond system links the extraframework occupants to the silicate layers and also interconnects the H2O molecules located inside the channels. Ca is seven-fold coordinated forming four bonds to O of the tetrahedral framework and three bonds to extraframework H2O. The H2O molecule at O9 showing a high displacement parameter is not bonded to Ca. The dehydration in pentagonite proceeds in three steps. At 100 °C the H2O molecule at O8 was released while O9 moved towards Ca. As a consequence the displacement parameter of H2O at O9 halved compared to that at room temperature. The unit-cell volume decreased to 1287.33(3) Å3 leading to a formula with 3H2O per formula unit (pfu). Ca remained seven-fold coordinated. At 175 °C Ca(VO)(Si4O10)·3H2O transformed into a new phase with 1H2O molecule pfu characterized by doubling of the c axis and the monoclinic space group Pn. Severe bending of specific TOT angles led to contraction of the porous three-dimensional framework. In addition, H2O at O9 was expelled while H2O at O7 approached a position in the center of the channel. The normalized volume decreased to 1069.44(9) Å3. The Ca coordination reduced from seven- to six-fold. At 225 °C a new anhydrous phase with space group Pna21 but without doubling of c had formed. Release of H2O at O7 caused additional contraction of TOT angles and volume reduction (V=1036.31(9) Å3). Ca adopted five-fold coordination. During heating excursion up to 400 °C this anhydrous phase remained preserved. Between room temperature and 225 °C the unit-cell volume decreased by 21% due to dehydration. The dehydration steps compare well with the thermo-gravimetric data reported in the literature.

Animal Guts as Ideal Chemical Reactors: Maximizing Absorption Rates

I solved equations that describe coupled hydrolysis in and absorption from a continuously stirred tank reactor (CSTR), a plug flow reactor (PFR), and a batch reactor (BR) for the rate of ingestion and/or the throughput time that maximizes the rate of absorption (=gross rate of gain from digestion). Predictions are that foods requiring a single hydrolytic step (e.g., disaccharides) yield ingestion rates that vary inversely with the concentration of food substrate ingested, whereas foods that require multiple hydrolytic and absorptive reactions proceeding in parallel (e.g., proteins) yield maximal ingestion rates at intermediate substrate concentrations. Counterintuitively, then, animals acting to maximize their absorption rates should show compensatory ingestion (more rapid feeding on food of lower concentration), except for the lower range of diet quality fur complex diets and except for animals that show purely linear (passive) uptake. At their respective maxima in absorption rates, the PFR and BR yield only modestly higher rates of gain than the CSTR but do so at substantially lower rates of ingestion. All three ideal reactors show milder than linear reduction in rate of absorption when throughput or holding time in the gut is increased (e.g., by scarcity or predation hazard); higher efficiency of hydrolysis and extraction offset lower intake. Hence adding feeding costs and hazards of predation is likely to slow ingestion rates and raise absorption efficiencies substantially over the cost-free optima found here.

Ultrathin strut biodegradable polymer sirolimus-eluting stent versus durable polymer everolimus-eluting stent for percutaneous coronary revascularisation (BIOSCIENCE): a randomised, single-blind, non-inferiority trial

BACKGROUND Refinements in stent design affecting strut thickness, surface polymer, and drug release have improved clinical outcomes of drug-eluting stents. We aimed to compare the safety and efficacy of a novel, ultrathin strut cobalt-chromium stent releasing sirolimus from a biodegradable polymer with a thin strut durable polymer everolimus-eluting stent. METHODS We did a randomised, single-blind, non-inferiority trial with minimum exclusion criteria at nine hospitals in Switzerland. We randomly assigned (1:1) patients aged 18 years or older with chronic stable coronary artery disease or acute coronary syndromes undergoing percutaneous coronary intervention to treatment with biodegradable polymer sirolimus-eluting stents or durable polymer everolimus-eluting stents. Randomisation was via a central web-based system and stratified by centre and presence of ST segment elevation myocardial infarction. Patients and outcome assessors were masked to treatment allocation, but treating physicians were not. The primary endpoint, target lesion failure, was a composite of cardiac death, target vessel myocardial infarction, and clinically-indicated target lesion revascularisation at 12 months. A margin of 3·5% was defined for non-inferiority of the biodegradable polymer sirolimus-eluting stent compared with the durable polymer everolimus-eluting stent. Analysis was by intention to treat. The trial is registered with ClinicalTrials.gov, number NCT01443104. FINDINGS Between Feb 24, 2012, and May 22, 2013, we randomly assigned 2119 patients with 3139 lesions to treatment with sirolimus-eluting stents (1063 patients, 1594 lesions) or everolimus-eluting stents (1056 patients, 1545 lesions). 407 (19%) patients presented with ST-segment elevation myocardial infarction. Target lesion failure with biodegradable polymer sirolimus-eluting stents (69 cases; 6·5%) was non-inferior to durable polymer everolimus-eluting stents (70 cases; 6·6%) at 12 months (absolute risk difference -0·14%, upper limit of one-sided 95% CI 1·97%, p for non-inferiority <0·0004). No significant differences were noted in rates of definite stent thrombosis (9 [0·9%] vs 4 [0·4%], rate ratio [RR] 2·26, 95% CI 0·70-7·33, p=0·16). In pre-specified stratified analyses of the primary endpoint, biodegradable polymer sirolimus-eluting stents were associated with improved outcome compared with durable polymer everolimus-eluting stents in the subgroup of patients with ST-segment elevation myocardial infarction (7 [3·3%] vs 17 [8·7%], RR 0·38, 95% CI 0·16-0·91, p=0·024, p for interaction=0·014). INTERPRETATION In a patient population with minimum exclusion criteria and high adherence to dual antiplatelet therapy, biodegradable polymer sirolimus-eluting stents were non-inferior to durable polymer everolimus-eluting stents for the combined safety and efficacy outcome target lesion failure at 12 months. The noted benefit in the subgroup of patients with ST-segment elevation myocardial infarction needs further study. FUNDING Clinical Trials Unit, University of Bern, and Biotronik, Bülach, Switzerland.

Fast parallel kriging-based stepwise uncertainty reduction with application to the identification of an excursion set

Stepwise uncertainty reduction (SUR) strategies aim at constructing a sequence of points for evaluating a function  f in such a way that the residual uncertainty about a quantity of interest progressively decreases to zero. Using such strategies in the framework of Gaussian process modeling has been shown to be efficient for estimating the volume of excursion of f above a fixed threshold. However, SUR strategies remain cumbersome to use in practice because of their high computational complexity, and the fact that they deliver a single point at each iteration. In this article we introduce several multipoint sampling criteria, allowing the selection of batches of points at which f can be evaluated in parallel. Such criteria are of particular interest when f is costly to evaluate and several CPUs are simultaneously available. We also manage to drastically reduce the computational cost of these strategies through the use of closed form formulas. We illustrate their performances in various numerical experiments, including a nuclear safety test case. Basic notions about kriging, auxiliary problems, complexity calculations, R code, and data are available online as supplementary materials.

Forward treatment planning for modulated electron radiotherapy (MERT) employing Monte Carlo methods

PURPOSE This paper describes the development of a forward planning process for modulated electron radiotherapy (MERT). The approach is based on a previously developed electron beam model used to calculate dose distributions of electron beams shaped by a photon multi leaf collimator (pMLC). METHODS As the electron beam model has already been implemented into the Swiss Monte Carlo Plan environment, the Eclipse treatment planning system (Varian Medical Systems, Palo Alto, CA) can be included in the planning process for MERT. In a first step, CT data are imported into Eclipse and a pMLC shaped electron beam is set up. This initial electron beam is then divided into segments, with the electron energy in each segment chosen according to the distal depth of the planning target volume (PTV) in beam direction. In order to improve the homogeneity of the dose distribution in the PTV, a feathering process (Gaussian edge feathering) is launched, which results in a number of feathered segments. For each of these segments a dose calculation is performed employing the in-house developed electron beam model along with the macro Monte Carlo dose calculation algorithm. Finally, an automated weight optimization of all segments is carried out and the total dose distribution is read back into Eclipse for display and evaluation. One academic and two clinical situations are investigated for possible benefits of MERT treatment compared to standard treatments performed in our clinics and treatment with a bolus electron conformal (BolusECT) method. RESULTS The MERT treatment plan of the academic case was superior to the standard single segment electron treatment plan in terms of organs at risk (OAR) sparing. Further, a comparison between an unfeathered and a feathered MERT plan showed better PTV coverage and homogeneity for the feathered plan, with V95% increased from 90% to 96% and V107% decreased from 8% to nearly 0%. For a clinical breast boost irradiation, the MERT plan led to a similar homogeneity in the PTV compared to the standard treatment plan while the mean body dose was lower for the MERT plan. Regarding the second clinical case, a whole breast treatment, MERT resulted in a reduction of the lung volume receiving more than 45% of the prescribed dose when compared to the standard plan. On the other hand, the MERT plan leads to a larger low-dose lung volume and a degraded dose homogeneity in the PTV. For the clinical cases evaluated in this work, treatment plans using the BolusECT technique resulted in a more homogenous PTV and CTV coverage but higher doses to the OARs than the MERT plans. CONCLUSIONS MERT treatments were successfully planned for phantom and clinical cases, applying a newly developed intuitive and efficient forward planning strategy that employs a MC based electron beam model for pMLC shaped electron beams. It is shown that MERT can lead to a dose reduction in OARs compared to other methods. The process of feathering MERT segments results in an improvement of the dose homogeneity in the PTV.

Effect of Salivary Contamination and Decontamination on Bond Strength of Two One-Step Self-Etching Adhesives to Dentin of Primary and Permanent Teeth

Purpose: To evaluate the effects of human saliva contamination and two decontamination procedures at different stages of the bonding procedure on the bond strength of two one-step self-etching adhesives to primary and permanent dentin. Materials and Methods: Extracted human primary and permanent molars (210 of each) were ground to mid-coronal dentin. The dentin specimens were randomly divided into 7 groups (n = 15/group/molar type) for each adhesive (Xeno V+ and Scotchbond Universal): no saliva contamination (control); saliva contamination before or after light curing of the adhesives followed by air drying, rinsing with water spray/air drying, or by rinsing with water spray/air drying/reapplication of the adhesives. Resin composite (Filtek Z250) was applied on the treated dentin surfaces. The specimens were stored at 37°C and 100% humidity for 24 h. After storage, shear bond strength (SBS) was measured and data analyzed with nonparametric ANOVA followed by exact Wilcoxon rank sum tests. Results: Xeno V+ generated significantly higher SBS than Scotchbond Universal when no saliva contamination occurred. Saliva contamination reduced SBS of Xeno V+, with the reduction being more pronounced when contamination occurred before light curing than after. In both situations, decontamination involving reapplication of the adhesive restored SBS. Saliva contamination had no significant effect on Scotchbond Universal. There were no differences in SBS between primary and permanent teeth. Conclusion: Rinsing with water and air drying followed by reapplication of the adhesive restored bond strength to saliva-contaminated dentin.

Probing the Electrocatalytic Oxygen Reduction Reaction Reactivity of Immobilized Multicopper Oxidase CueO

The bioelectrocatalytic (oxygen reduction reaction, ORR) properties of the multicopper oxidase CueO immobilized on gold electrodes were investigated. Macroscopic electrochemical techniques were combined with in situ scanning tunneling microscopy (STM) and surface-enhanced Raman spectroscopy at the ensemble and at the single-molecule level. Self-assembled monolayer of mercaptopropionic acid, cysteamine, and p-aminothiophenol were chosen as redox mediators. The highest ORR activity was observed for the protein attached to amino-terminated adlayers. In situ STM experiments revealed that the presence of oxygen causes distinct structure and electronic changes in the metallic centers of the enzyme, which determine the rate of intramolecular electron transfer and, consequently, affect the rate of electron tunneling through the protein. Complementary Raman spectroscopy experiments provided access for monitoring structural changes in the redox state of the type 1 copper center of the immobilized enzyme during the CueO-catalyzed oxygen reduction cycle. These results unequivocally demonstrate the existence of a direct electronic communication between the electrode substrate and the type 1 copper center.

Dehydration of the zeolite merlinoite from the Khibiny massif, Russia: an in situ temperature-dependent single-crystal X-ray study

Dehydration behaviour of the zeolite merlinoite, NaK11[Al12Si20O64]·15H2O, from the Khibiny massif (Russia) was studied by means of single-crystal X-ray diffraction conjoined with step-wise heating to 225 C. At room temperature merlinoite has the space group Immm with a = 14.0312(5), b = 14.2675(6), c = 10.0874(4) Å, and V = 2019.40(14) Å3. At 75 °C the merlinoite structure undergoes pronounced dehydration accompanied by a phase transition to a structure that has the space group P42/nmc and remains consistent at elevated temperature. A fully dehydrated phase occurs at 200 °C (at 225 °C: a = 13.341(4), b = 13.341(4), c = 9.707(4) Å, V = 1727.7(12) Å3). Dehydration-induced framework distortion and symmetry were found to be different from those observed for synthetic potassium merlinoite with the K11.5[Al11.5Si20.5O64]·15H2O composition.

In-vitro activity of taurolidine on single species and a multispecies population associated with periodontitis.

The antimicrobial activity of taurolidine was compared with minocycline against microbial species associated with periodontitis (four single strains and a 12-species mixture). Minimal inhibitory concentrations (MICs) and minimal bactericidal concentrations (MBCs), killing as well as activities on established and forming single-species biofilms and a 12-species biofilm were determined. The MICs of taurolidine against single species were always 0.31 mg/ml, the MBCs were 0.64 mg/ml. The used mixed microbiota was less sensitive to taurolidine, MIC and the MBC was 2.5 mg/ml. The strains and the mixture were completely killed by 2.5 mg/ml taurolidine, whereas 256 μg/ml minocycline reduced the bacterial counts of the mixture by 5 log10 colony forming units (cfu). Coating the surface with 10 mg/ml taurolidine or 256 μg/ml minocycline prevented completely biofilm formation of Porphyromonas gingivalis ATCC 33277 but not of Aggregatibacter actinomycetemcomitans Y4 and the mixture. On 4.5 d old biofilms, taurolidine acted concentration dependent with a reduction by 5 log10 cfu (P. gingivalis ATCC 33277) and 7 log10 cfu (A. actinomycetemcomitans Y4) when applying 10 mg/ml. Minocycline decreased the cfu counts by 1-2 log10 cfu independent of the used concentration. The reduction of the cfu counts in the 4.5 d old multi-species biofilms was about 3 log10 cfu after application of any minocycline concentration and after using 10 mg/ml taurolidine. Taurolidine is active against species associated with periodontitis, even within biofilms. Nevertheless a complete elimination of complex biofilms by taurolidine seems to be impossible and underlines the importance of a mechanical removal of biofilms prior to application of taurolidine.

Probing the effect of minor groove interactions on the catalytic efficiency of DNAzymes 8–17 and 10–23

DNAzymes (Dz) 8–17 and 10–23 are two widely studied and well-characterized RNA-cleaving DNA catalysts. In an effort to further improve the understanding of the fragile interactions and dynamics of the enzymatic mechanism, this study examines the catalytic efficiency of minimally modified DNAzymes. Five single mutants of Dz8–17 and Dz10–23 were prepared by replacing the adenine residues in the corresponding catalytic cores with 3-deazaadenine units. Kinetic assays were used to assess the effect on the catalytic activity and thereby identify the importance of hydrogen bonding that arises from the N3 atoms. The results suggest that modifications at A15 and A15.0 of Dz8–17 have a significant influence and show a reduction in catalytic activity. Modification at each location in Dz10–23 results in a decrease of the observed rate constants, with A12 appearing to be the most affected with a reduction of ∼80% of kobs and ∼25% of the maximal cleavage rate compared to the wild-type DNAzyme. On the other hand, modification of A12 in Dz8–17 showed an ∼130% increase in kobs, thus unraveling a new potential site for the introduction of chemical modifications. A pH-profile analysis showed that the chemical cleavage step is rate-determining, regardless of the presence and/or location of the mutation. These findings point towards the importance of the N3-nitrogens of certain adenine nucleotides located within the catalytic cores of the DNAzymes for efficient catalytic activity and further suggest that they might directly partake in maintaining the appropriate tertiary structure. Therefore, it appears that minor groove interactions constitute an important feature of DNAzymes as well as ribozymes.