Towards a computational reconstruction of the electrodynamics of premature and full term human labour

We apply virtual tissue engineering to the full term human uterus with a view to reconstruction of the spatiotemporal patterns of electrical activity of the myometrium that control mechanical activity via intracellular calcium. The three-dimensional geometry of the gravid uterus has been reconstructed from segmented in vivo magnetic resonance imaging as well as ex vivo diffusion tensor magnetic resonance imaging to resolve fine scale tissue architecture. A late-pregnancy uterine smooth muscle cell model is constructed and bursting analysed using continuation algorithms. These cell models are incorporated into partial differential equation models for tissue synchronisation and propagation. The ultimate objective is to develop a quantitative and predictive understanding of the mechanisms that initiate and regulate labour.     

Simultaneous modelling of metabolic, genetic and product-interaction networks

The creation of cell models from annotated genome information, as well as additional data from other databases, requires both a format and medium for its distribution. Standards are described for the representation of the data in the form of Document Type Definitions (DTDs) for XML files. Separate DTDs are detailed for genetic, metabolic and gene product-interaction networks, which can be used to hold information on individual subsystems, or which may be combined to create a whole cell DTD. In the execution of this work, a fifth DTD was also created for a metabolite thesaurus, which allows incorporation of metabolite synonyms and generic nomenclature data into the models. A gene-regulation classification scheme was also created, to facilitate incorporation of gene regulatory information in an efficient manner. The work is described with particular reference to the metabolic network of Escherichia coli, which contains 808 individual enzymes. The assignment of confidence levels to these data, through the use of Gene Ontology evidence codes, is highlighted. In silico investigations may now be performed using the mathematical simulation workbench, DBsolve, which incorporates the facility to introduce data directly from XML.     

Temperature dependence of the epidermal growth factor receptor signaling network can be accounted for by a kinetic model.

Stimulation of isolated hepatocytes with epidermal growth factor (EGF) causes rapid tyrosine phosphorylation of the EGF receptor (EGFR) and adapter/target proteins, which was monitored with 1 and 2 s resolution at 37, 20, and 4 degrees C. The temporal responses detected for multiple signaling proteins involve both transient and sustained phosphorylation patterns, which change dramatically at low temperatures. To account quantitatively for complex responses, we employed a mechanistic kinetic model of the EGFR pathway, formulated in molecular terms as cascades of protein interactions and phosphorylation and dephosphorylation reactions. Assuming differential temperature dependencies for different reaction groups, such as SH2 and PTB domain-mediated interactions, the EGFR kinase, and the phosphatases, good quantitative agreement was obtained between computer-simulated and measured responses. The kinetic model demonstrates that, for each protein-protein interaction, the dissociation rate constant, k(off), strongly decreases at low temperatures, whereas this decline may or may not be accompanied by a large decrease in the k(on) value. Temperature-induced changes in the maximal activities of the reactions catalyzed by the EGFR kinase were moderate, compared to such changes in the V(max) of the phosphatases. However, strong changes in both the V(max) and K(m) for phosphatases resulted in moderate changes in the V(max)/K(m) ratio, comparable to the corresponding changes in EGFR kinase activity, with a single exception for the receptor phosphatase at 4 degrees C. The model suggests a significant decrease in the rates of the EGF receptor dimerization and its dephosphorylation at 4 degrees C, which can be related to the phase transition in the membrane lipids. A combination of high-resolution experimental monitoring and molecular level kinetic modeling made it possible to quantitatively account for the temperature dependence of the integrative signaling responses.     

Trigger factor assists the refolding of heterodimeric but not monomeric luciferases

The refolding of thermally inactivated protein by ATP-independent trigger factor (TF) and ATP-dependent DnaKJE chaperones was comparatively analyzed. Heterodimeric (αβ) bacterial luciferases of Aliivibrio fischeri, Photobacterium leiognathi, and Vibrio harveyi as well as monomeric luciferases of Vibrio harveyi and Luciola mingrelica (firefly) were used as substrates. In the presence of TF, thermally inactivated heterodimeric bacterial luciferases refold, while monomeric luciferases do not refold. These observations were made both in vivo (Escherichia coli ΔdnaKJ containing plasmids with tig gene) and in vitro (purified TF). Unlike TF, the DnaKJE chaperone system refolds both monomeric and heterodimeric luciferases with equal efficiency.     

The metabolic pathway collection: an update.

The Metabolic Pathway Collection from EMP is an extraction of data from the larger Enzymes and Metabolic Pathways database (EMP). This extraction has been made publicly available in the hope that others will find it useful for a variety of purposes. The original release in October 1995 contained 1814 distinct pathways. The current collection contains 2180. Metabolic reconstructions for the first completely sequenced organisms-Haemophilus influenzae,Mycoplasma genitalium,Saccharomyces cerevisiaeandMethanococcus janaschii-are all included in the current release. All of the pathways in the collections are available as ASCII files in the form generated by the main curator, Evgeni Selkov. In addition, we are offering a more structured encoding of a subset of the collection; our initial release of this subcollection includes all of the pathways inMycoplasma genitalium, and we ultimately intend to offer the entire collection in this form as well.     

PpCas9 from Pasteurella pneumotropica — a compact Type II-C Cas9 ortholog active in human cells

CRISPR-Cas defense systems opened up the field of genome editing due to the ease with which effector Cas nucleases can be programmed with guide RNAs to access desirable genomic sites. Type II-A SpCas9 from Streptococcus pyogenes was the first Cas9 nuclease used for genome editing and it remains the most popular enzyme of its class. Nevertheless, SpCas9 has some drawbacks including a relatively large size and restriction to targets flanked by an ‘NGG’ PAM sequence. The more compact Type II-C Cas9 orthologs can help to overcome the size limitation of SpCas9. Yet, only a few Type II-C nucleases were fully characterized to date. Here, we characterized two Cas9 II-C orthologs, DfCas9 from Defluviimonas sp.20V17 and PpCas9 from Pasteurella pneumotropica. Both DfCas9 and PpCas9 cleave DNA in vitro and have novel PAM requirements. Unlike DfCas9, the PpCas9 nuclease is active in human cells. This small nuclease requires an ‘NNNNRTT’ PAM orthogonal to that of SpCas9 and thus potentially can broaden the range of Cas9 applications in biomedicine and biotechnology.     

DBSolve Optimum: a software package for kinetic modeling which allows dynamic visualization of simulation results

Background  

Systems biology research and applications require creation, validation, extensive usage of mathematical models and visualization of simulation results by end-users. Our goal is to develop novel method for visualization of simulation results and implement it in simulation software package equipped with the sophisticated mathematical and computational techniques for model development, verification and parameter fitting.     

Kinetic Modeling of Energy Metabolism and Superoxide Generation in Hepatocyte Mitochondria

Direct nonenzymatic oxidation of semiquinone by oxygen is one of the main sources of superoxide radicals in mitochondria. Using all the known data on hepatocyte mitochondria, we have revealed the correlation between the rate of superoxide generation by the bc ₁complex and the transmembrane potential (). Assuming that the main electrogenic stage of the Qcycle is the electron transfer between the cytochrome bhemes, then the rate of superoxide generation sharply increases when grows from 150 to 180 mV. However, this interrelation is ambiguous. Indeed, the increase of the generation rate with the growth of the potential can occur faster when succinate dehydrogenase is inhibited by malonate than when external ADP is exhausted. When the potential is changed by adding phosphate or potassium (K⁺), the rate of production remains constant, although the comparison of the rates at the same reveals the effect of phosphate or potassium. It turned out that the rate of generation is a function of rather than any of its components. Phosphate and K⁺have practically no influence on , since the change in is compensated by pH. The rate of superoxide generation by the bc ₁complex is a multiple function of the electron-transfer activity of enzymes, the processes determining the membrane potential (e.g., loading), and the oxygen concentration. The kinetic model proposed in this work may serve to understand how the superoxide production is regulated.