A Pycellerator Tutorial.

We present a tutorial on using Pycellerator for biomolecular simulations. Models are described in human readable (and editable) text files (UTF8 or ASCII) containing collections of reactions, assignments, initial conditions, function definitions, and rate constants. These models are then converted into a Python program that can optionally solve the system, e.g., as a system of differential equations using ODEINT, or be run by another program. The input language implements an extended version of the Cellerator arrow notation, including mass action, Hill functions, S-Systems, MWC, and reactions with user-defined kinetic laws. Simple flux balance analysis is also implemented. We will demonstrate the implementation and analysis of progressively more complex models, starting from simple mass action through indexed cascades. Pycellerator can be used as a library that is integrated into other programs, run as a command line program, or in iPython notebooks. It is implemented in Python 2.7 and available under an open source GPL license

Modeling the organization of the WUSCHEL expression domain in the shoot apical meristem

Motivation: The above-ground tissues of higher plants are generated from a small region of cells situated at the plant apex called the shoot apical meristem. An important genetic control circuit modulating the size of the Arabidopsis thaliana meristem is a feed-back network between the CLAVATA3 and WUSCHEL genes. Although the expression patterns for these genes do not overlap, WUSCHEL activity is both necessary and sufficient (when expressed ectopically) for the induction of CLAVATA3 expression. However, upregulation of CLAVATA3 in conjunction with the receptor kinase CLAVATA1 results in the downregulation of WUSCHEL. Despite much work, experimental data for this network are incomplete and additional hypotheses are needed to explain the spatial locations and dynamics of these expression domains. Predictive mathematical models describing the system should provide a useful tool for investigating and discriminating among possible hypotheses, by determining which hypotheses best explain observed gene expression dynamics. Results: We are developing a method using in vivo live confocal microscopy to capture quantitative gene expression data and create templates for computational models. We present two models accounting for the organization of the WUSCHEL expression domain. Our preferred model uses a reaction-diffusion mechanism in which an activator induces WUSCHEL expression. This model is able to organize the WUSCHEL expression domain. In addition, the model predicts the dynamical reorganization seen in experiments where cells, including the WUSCHEL domain, are ablated, and it also predicts the spatial expansion of the WUSCHEL domain resulting from removal of the CLAVATA3 signal

A Dynamic Model of Interactions of Ca^(2+), Calmodulin, and Catalytic Subunits of Ca^(2+)/Calmodulin-Dependent Protein Kinase II

During the acquisition of memories, influx of Ca^(2+) into the postsynaptic spine through the pores of activated N-methyl-D-aspartate-type glutamate receptors triggers processes that change the strength of excitatory synapses. The pattern of Ca^(2+) influx during the first few seconds of activity is interpreted within the Ca^(2+)-dependent signaling network such that synaptic strength is eventually either potentiated or depressed. Many of the critical signaling enzymes that control synaptic plasticity, including Ca^(2+)/calmodulin-dependent protein kinase II (CaMKII), are regulated by calmodulin, a small protein that can bind up to 4 Ca^(2+) ions. As a first step toward clarifying how the Ca^(2+)-signaling network decides between potentiation or depression, we have created a kinetic model of the interactions of Ca^(2+), calmodulin, and CaMKII that represents our best understanding of the dynamics of these interactions under conditions that resemble those in a postsynaptic spine. We constrained parameters of the model from data in the literature, or from our own measurements, and then predicted time courses of activation and autophosphorylation of CaMKII under a variety of conditions. Simulations showed that species of calmodulin with fewer than four bound Ca^(2+) play a significant role in activation of CaMKII in the physiological regime, supporting the notion that processing ofCa^(2+) signals in a spine involves competition among target enzymes for binding to unsaturated species of CaM in an environment in which the concentration of Ca^(2+) is fluctuating rapidly. Indeed, we showed that dependence of activation on the frequency of Ca^(2+) transients arises from the kinetics of interaction of fluctuating Ca^(2+) with calmodulin/CaMKII complexes. We used parameter sensitivity analysis to identify which parameters will be most beneficial to measure more carefully to improve the accuracy of predictions. This model provides a quantitative base from which to build more complex dynamic models of postsynaptic signal transduction during learning

Inverse bifurcation analysis: application to simple gene systems

BACKGROUND: Bifurcation analysis has proven to be a powerful method for understanding the qualitative behavior of gene regulatory networks. In addition to the more traditional forward problem of determining the mapping from parameter space to the space of model behavior, the inverse problem of determining model parameters to result in certain desired properties of the bifurcation diagram provides an attractive methodology for addressing important biological problems. These include understanding how the robustness of qualitative behavior arises from system design as well as providing a way to engineer biological networks with qualitative properties. RESULTS: We demonstrate that certain inverse bifurcation problems of biological interest may be cast as optimization problems involving minimal distances of reference parameter sets to bifurcation manifolds. This formulation allows for an iterative solution procedure based on performing a sequence of eigen-system computations and one-parameter continuations of solutions, the latter being a standard capability in existing numerical bifurcation software. As applications of the proposed method, we show that the problem of maximizing regions of a given qualitative behavior as well as the reverse engineering of bistable gene switches can be modelled and efficiently solved

Hadronic Mass Moments in Inclusive Semileptonic B Meson Decays

We have measured the first and second moments of the hadronic mass-squared distribution in B -> X_c l nu, for P(lepton) > 1.5 GeV/c. We find <M_X^2 - M_D[Bar]^2> = 0.251 +- 0.066 GeV^2, )^2 > = 0.576 +- 0.170 GeV^4, where M_D[Bar] is the spin-averaged D meson mass. From that first moment and the first moment of the photon energy spectrum in b -> s gamma, we find the HQET parameter lambda_1 (MS[Bar], to order 1/M^3 and beta_0 alpha_s^2) to be -0.24 +- 0.11 GeV^2. Using these first moments and the B semileptonic width, and assuming parton-hadron duality, we obtain |V_cb| = 0.0404 +- 0.0013.Comment: 11 pages postscript, also available through http://w4.lns.cornell.edu/public/CLNS, submitted to PR

Observation of the Ωc0\Omega_{c}^{0} Charmed Baryon at CLEO

The CLEO experiment at the CESR collider has used 13.7 fb$^{-1}$ of data to search for the production of the $\Omega_c^0$ (css-ground state) in $e^{+}e^{-}$ collisions at $\sqrt{s} \simeq 10.6$ {\rm GeV}. The modes used to study the $\Omega_c^0$ are $\Omega^- \pi^+$, $\Omega^- \pi^+ \pi^0$, $\Xi^- K^- pi^+ \pi^+$, $\Xi^0 K^- pi^+$, and $\Omega^- \pi^+ \pi^- \pi^+$. We observe a signal of 40.4$\pm$9.0(stat) events at a mass of 2694.6$\pm$2.6(stat)$\pm$1.9(syst) {\rm MeV/$c^2$}, for all modes combined.Comment: 10 pages postscript, also available through http://w4.lns.cornell.edu/public/CLN

SBMLsqueezer: A CellDesigner plug-in to generate kinetic rate equations for biochemical networks

<p>Abstract</p> <p>Background</p> <p>The development of complex biochemical models has been facilitated through the standardization of machine-readable representations like SBML (Systems Biology Markup Language). This effort is accompanied by the ongoing development of the human-readable diagrammatic representation SBGN (Systems Biology Graphical Notation). The graphical SBML editor CellDesigner allows direct translation of SBGN into SBML, and vice versa. For the assignment of kinetic rate laws, however, this process is not straightforward, as it often requires manual assembly and specific knowledge of kinetic equations.</p> <p>Results</p> <p>SBMLsqueezer facilitates exactly this modeling step via automated equation generation, overcoming the highly error-prone and cumbersome process of manually assigning kinetic equations. For each reaction the kinetic equation is derived from the stoichiometry, the participating species (e.g., proteins, mRNA or simple molecules) as well as the regulatory relations (activation, inhibition or other modulations) of the SBGN diagram. Such information allows distinctions between, for example, translation, phosphorylation or state transitions. The types of kinetics considered are numerous, for instance generalized mass-action, Hill, convenience and several Michaelis-Menten-based kinetics, each including activation and inhibition. These kinetics allow SBMLsqueezer to cover metabolic, gene regulatory, signal transduction and mixed networks. Whenever multiple kinetics are applicable to one reaction, parameter settings allow for user-defined specifications. After invoking SBMLsqueezer, the kinetic formulas are generated and assigned to the model, which can then be simulated in CellDesigner or with external ODE solvers. Furthermore, the equations can be exported to SBML, LaTeX or plain text format.</p> <p>Conclusion</p> <p>SBMLsqueezer considers the annotation of all participating reactants, products and regulators when generating rate laws for reactions. Thus, for each reaction, only applicable kinetic formulas are considered. This modeling scheme creates kinetics in accordance with the diagrammatic representation. In contrast most previously published tools have relied on the stoichiometry and generic modulators of a reaction, thus ignoring and potentially conflicting with the information expressed through the process diagram. Additional material and the source code can be found at the project homepage (URL found in the Availability and requirements section).</p

The use of pure and impure placebo interventions in primary care - a qualitative approach

Background: Placebos play an important role in clinical trials and several surveys have shown that they are also common in daily practice. Previous research focused primarily on the frequency of placebo use in outpatient care. Our aim was to explore physicians' views on the use of placebos in daily practice, whereby distinction was made between pure placebos (substances with no pharmacological effect, e.g. sugar pills) and impure placebos (substances with pharmacological effect but not on the condition being treated, e.g. antibiotics in viral infections or vitamins). Methods: We performed semi-structured interviews with a sample of twelve primary care physicians (PCPs). The interview addressed individual definitions of a placebo, attitudes towards placebos and the participants' reasons for prescribing them. The interviews were transcribed and analysed using qualitative content analysis. Results: The definition of a placebo given by the majority of the PCPs in our study was one which actually only describes pure placebos. This definition, combined with the fact that most impure placebos were not regarded as placebos at all, means that most of the participating PCPs were not aware of the extent to which placebos are used in daily practice. The PCPs stated that they use placebos (both pure and impure) mainly in the case of non-severe diseases for which there was often no satisfactory somatic explanation. According to the PCPs, cases like this are often treated by complementary and alternative therapies and these, too, are associated with placebo effects. However, all PCPs felt that the ethical aspects of such treatment were unclear and they were unsure as to how to communicate the use of placebos to their patients. Most of them would appreciate ethical guidelines on how to deal with this issue. Conclusions: Many PCPs seem to be unaware that some of the drugs they prescribe are classified as impure placebos. Perceptions of effectiveness and doubts about the legal and ethical aspects of the use of placebos by PCPs may discourage their application. Dissemination of guidelines and consensus papers may be one approach, but it has to be acknowledged that the topic itself is in conflict with the PCPs' perception of themselves as professional and reliable physicians

On the effects of mechanical stress of biological membranes in modeling of swelling dynamics of biological systems

We highlight mechanical stretching and bending of membranes and the importance of membrane deformations in the analysis of swelling dynamics of biological systems, including cells and subcellular organelles. Membrane deformation upon swelling generates tensile stress and internal pressure, contributing to volume changes in biological systems. Therefore, in addition to physical (internal/external) and chemical factors, mechanical properties of the membranes should be considered in modeling analysis of cellular swelling. Here we describe an approach that considers mechanical properties of the membranes in the analysis of swelling dynamics of biological systems. This approach includes membrane bending and stretching deformations into the model, producing a more realistic description of swelling. We also discuss the effects of membrane stretching on swelling dynamics. We report that additional pressure generated by membrane bending is negligible, compared to pressures generated by membrane stretching, when both membrane surface area and volume are variable parameters. Note that bending deformations are reversible, while stretching deformation may be irreversible, leading to membrane disruption when they exceed a certain threshold level. Therefore, bending deformations need only be considered in reversible physiological swelling, whereas stretching deformations should also be considered in pathological irreversible swelling. Thus, the currently proposed approach may be used to develop a detailed biophysical model describing the transition from physiological to pathological swelling mode.National Aeronautics & Space Administration (NASA):80NSSC19M0049; PR Space Grant (NASA):NNX15AI11Hinfo:eu-repo/semantics/publishedVersio