Coupled Systems of Differential-Algebraic and Kinetic Equations with Application to the Mathematical Modelling of Muscle Tissue

We consider a coupled system composed of a linear differential-algebraic equation (DAE) and a linear large-scale system of ordinary differential equations where the latter stands for the dynamics of numerous identical particles. Replacing the discrete particles by a kinetic equation for a particle density, we obtain in the mean-field limit the new class of partially kinetic systems. We investigate the influence of constraints on the kinetic theory of those systems and present necessary adjustments. We adapt the mean-field limit to the DAE model and show that index reduction and the mean-field limit commute. As a main result, we prove Dobrushin's stability estimate for linear systems. The estimate implies convergence of the mean-field limit and provides a rigorous link between the particle dynamics and their kinetic description. Our research is inspired by mathematical models for muscle tissue where the macroscopic behaviour is governed by the equations of continuum mechanics, often discretised by the finite element method, and the microscopic muscle contraction process is described by Huxley's sliding filament theory. The latter represents a kinetic equation that characterises the state of the actin-myosin bindings in the muscle filaments. Linear partially kinetic systems are a simplified version of such models, with focus on the constraints.Comment: 32 pages, 18 figure

Contrasting cellular uptake pathways for chlorido and iodido iminopyridine ruthenium arene anticancer complexes

The pathways involved in cellular uptake and accumulation of iminopyridine complexes of general formula [Ru(η6-p-cymene)(N,N-dimethyl-N′-[(E)-pyridine-2-ylmethylidene]benzene-1,4-diamine)X]PF6 bearing two different halido ligands X = Cl or I, have been explored. The ratio of passive/active cellular accumulation of Ru in A2780 human ovarian cancer cells is compared and contrasted with cisplatin. Also, saturation of cellular uptake, time-dependence of cellular influx/efflux equilibria, together with endocytotic pathways such as caveolae and facilitated diffusion are investigated and discussed. Temperature dependence studies of Ru accumulation in the A2780 cells show that in contrast to cisplatin (CDDP) and chlorido complex 1, which are taken up largely through active transport, the iodido complex 2 enters cells via passive transport. The cellular efflux of Ru is slow (ca. 25% retained after 72 h) and is partially inhibited by verapamil, implicating the P-gp protein in the efflux mechanism. Ouabain inhibition experiments suggest that the cellular uptake of these ruthenium complexes relies at least in part on facilitated diffusion, and in particular is dependent on the membrane potential. In addition the finding that depletion of cellular ATP with antimycin A had little effect on cellular Ru accumulation from iodido complex 2 is consistent with passive diffusion. In contrast, ATP depletion caused a major increase in cellular accumulation of ruthenium from chlorido complex 1

Spatio-Temporal Dependence of the Signaling Response in Immune-Receptor Trafficking Networks Regulated by Cell Density: A Theoretical Model

Cell signaling processes involve receptor trafficking through highly connected networks of interacting components. The binding of surface receptors to their specific ligands is a key factor for the control and triggering of signaling pathways. In most experimental systems, ligand concentration and cell density vary within a wide range of values. Dependence of the signal response on cell density is related with the extracellular volume available per cell. This dependence has previously been studied using non-spatial models which assume that signaling components are well mixed and uniformly distributed in a single compartment. In this paper, a mathematical model that shows the influence exerted by cell density on the spatio-temporal evolution of ligands, cell surface receptors, and intracellular signaling molecules is developed. To this end, partial differential equations were used to model ligand and receptor trafficking dynamics through the different domains of the whole system. This enabled us to analyze several interesting features involved with these systems, namely: a) how the perturbation caused by the signaling response propagates through the system; b) receptor internalization dynamics and how cell density affects the robustness of dose-response curves upon variation of the binding affinity; and c) that enhanced correlations between ligand input and system response are obtained under conditions that result in larger perturbations of the equilibrium . Finally, the results are compared with those obtained by considering that the above components are well mixed in a single compartment

Transcriptional Regulation: Effects of Promoter Proximal Pausing on Speed, Synchrony and Reliability

Recent whole genome polymerase binding assays in the Drosophila embryo have shown that a substantial proportion of uninduced genes have pre-assembled RNA polymerase-II transcription initiation complex (PIC) bound to their promoters. These constitute a subset of promoter proximally paused genes for which mRNA elongation instead of promoter access is regulated. This difference can be described as a rearrangement of the regulatory topology to control the downstream transcriptional process of elongation rather than the upstream transcriptional initiation event. It has been shown experimentally that genes with the former mode of regulation tend to induce faster and more synchronously, and that promoter-proximal pausing is observed mainly in metazoans, in accord with a posited impact on synchrony. However, it has not been shown whether or not it is the change in the regulated step per se that is causal. We investigate this question by proposing and analyzing a continuous-time Markov chain model of PIC assembly regulated at one of two steps: initial polymerase association with DNA, or release from a paused, transcribing state. Our analysis demonstrates that, over a wide range of physical parameters, increased speed and synchrony are functional consequences of elongation control. Further, we make new predictions about the effect of elongation regulation on the consistent control of total transcript number between cells. We also identify which elements in the transcription induction pathway are most sensitive to molecular noise and thus possibly the most evolutionarily constrained. Our methods produce symbolic expressions for quantities of interest with reasonable computational effort and they can be used to explore the interplay between interaction topology and molecular noise in a broader class of biochemical networks. We provide general-purpose code implementing these methods

An ultrasoft X-ray multi-microbeam irradiation system for studies of DNA damage responses by fixed- and live-cell fluorescence microscopy

Localized induction of DNA damage is a valuable tool for studying cellular DNA damage responses. In recent decades, methods have been developed to generate DNA damage using radiation of various types, including photons and charged particles. Here we describe a simple ultrasoft X-ray multi-microbeam system for high dose-rate, localized induction of DNA strand breaks in cells at spatially and geometrically adjustable sites. Our system can be combined with fixed- and live-cell microscopy to study responses of cells to DNA damage

Phosphoproteomics-Based Modeling Defines the Regulatory Mechanism Underlying Aberrant EGFR Signaling

BACKGROUND: Mutation of the epidermal growth factor receptor (EGFR) results in a discordant cell signaling, leading to the development of various diseases. However, the mechanism underlying the alteration of downstream signaling due to such mutation has not yet been completely understood at the system level. Here, we report a phosphoproteomics-based methodology for characterizing the regulatory mechanism underlying aberrant EGFR signaling using computational network modeling. METHODOLOGY/PRINCIPAL FINDINGS: Our phosphoproteomic analysis of the mutation at tyrosine 992 (Y992), one of the multifunctional docking sites of EGFR, revealed network-wide effects of the mutation on EGF signaling in a time-resolved manner. Computational modeling based on the temporal activation profiles enabled us to not only rediscover already-known protein interactions with Y992 and internalization property of mutated EGFR but also further gain model-driven insights into the effect of cellular content and the regulation of EGFR degradation. Our kinetic model also suggested critical reactions facilitating the reconstruction of the diverse effects of the mutation on phosphoproteome dynamics. CONCLUSIONS/SIGNIFICANCE: Our integrative approach provided a mechanistic description of the disorders of mutated EGFR signaling networks, which could facilitate the development of a systematic strategy toward controlling disease-related cell signaling

A domain-based approach to predict protein-protein interactions

<p>Abstract</p> <p>Background</p> <p>Knowing which proteins exist in a certain organism or cell type and how these proteins interact with each other are necessary for the understanding of biological processes at the whole cell level. The determination of the protein-protein interaction (PPI) networks has been the subject of extensive research. Despite the development of reasonably successful methods, serious technical difficulties still exist. In this paper we present DomainGA, a quantitative computational approach that uses the information about the domain-domain interactions to predict the interactions between proteins.</p> <p>Results</p> <p>DomainGA is a multi-parameter optimization method in which the available PPI information is used to derive a quantitative scoring scheme for the domain-domain pairs. Obtained domain interaction scores are then used to predict whether a pair of proteins interacts. Using the yeast PPI data and a series of tests, we show the robustness and insensitivity of the DomainGA method to the selection of the parameter sets, score ranges, and detection rules. Our DomainGA method achieves very high explanation ratios for the positive and negative PPIs in yeast. Based on our cross-verification tests on human PPIs, comparison of the optimized scores with the structurally observed domain interactions obtained from the iPFAM database, and sensitivity and specificity analysis; we conclude that our DomainGA method shows great promise to be applicable across multiple organisms.</p> <p>Conclusion</p> <p>We envision the DomainGA as a first step of a multiple tier approach to constructing organism specific PPIs. As it is based on fundamental structural information, the DomainGA approach can be used to create potential PPIs and the accuracy of the constructed interaction template can be further improved using complementary methods. Explanation ratios obtained in the reported test case studies clearly show that the false prediction rates of the template networks constructed using the DomainGA scores are reasonably low, and the erroneous predictions can be filtered further using supplementary approaches such as those based on literature search or other prediction methods.</p

Rapid and Accurate Prediction and Scoring of Water Molecules in Protein Binding Sites

Water plays a critical role in ligand-protein interactions. However, it is still challenging to predict accurately not only where water molecules prefer to bind, but also which of those water molecules might be displaceable. The latter is often seen as a route to optimizing affinity of potential drug candidates. Using a protocol we call WaterDock, we show that the freely available AutoDock Vina tool can be used to predict accurately the binding sites of water molecules. WaterDock was validated using data from X-ray crystallography, neutron diffraction and molecular dynamics simulations and correctly predicted 97% of the water molecules in the test set. In addition, we combined data-mining, heuristic and machine learning techniques to develop probabilistic water molecule classifiers. When applied to WaterDock predictions in the Astex Diverse Set of protein ligand complexes, we could identify whether a water molecule was conserved or displaced to an accuracy of 75%. A second model predicted whether water molecules were displaced by polar groups or by non-polar groups to an accuracy of 80%. These results should prove useful for anyone wishing to undertake rational design of new compounds where the displacement of water molecules is being considered as a route to improved affinity

EMSL Pore Scale Modeling Challenge/Workshop

Report covers the background for the workshop, objectives, important research directions, necessary capabilities and overall recommendations