A stochastic theory of phase transitions in human hand movement

The order parameter equation for the relative phase of correlated hand movements, derived in a previous paper by Haken et al. (1985), is extended to a time-dependent stochastic differential equation. Its solutions are determined close to stationary points and for the transition region. Remarkably good agreement between this theory and recent experiments done by Kelso and Scholz (1985) is found, and new predictions are offered.     

Exponential equilibration of genetic circuits using entropy methods

Journal of Mathematical Biology - We analyse a continuum model for genetic circuits based on a partial integro-differential equation initially proposed in Friedman et al. (Phys Rev Lett...     

Microstructural analysis of deformation-induced hypoxic damage in skeletal muscle

Deep pressure ulcers are caused by sustained mechanical loading and involve skeletal muscle tissue injury. The exact underlying mechanisms are unclear, and the prevalence is high. Our hypothesis is that the aetiology is dominated by cellular deformation (Bouten et al. in Ann Biomed Eng 29:153-163, 2001; Breuls et al. in Ann Biomed Eng 31:1357-1364, 2003; Stekelenburg et al. in J App Physiol 100(6):1946-1954, 2006) and deformation-induced ischaemia. The experimental observation that mechanical compression induced a pattern of interspersed healthy and dead cells in skeletal muscle (Stekelenburg et al. in J App Physiol 100(6):1946-1954, 2006) strongly suggests to take into account the muscle microstructure in studying damage development. The present paper describes a computational model for deformation-induced hypoxic damage in skeletal muscle tissue. Dead cells stop consuming oxygen and are assumed to decrease in stiffness due to loss of structure. The questions addressed are if these two consequences of cell death influence the development of cell injury in the remaining cells. The results show that weakening of dead cells indeed affects the damage accumulation in other cells. Further, the fact that cells stop consuming oxygen after they have died, delays cell death of other cells.     

Genome-Wide Analysis of Gene Expression in Stationary Phase and Genetic Characterization of Stationary-Phase-Dependent Halocin Gene Expression in the Haloarchaeon Haloferax mediterranei

The stationary phase of microbial growth is a very complex state regulated by various environmental and physiological factors.An intensive study of stationary phase could promote a comprehensive understanding of the complete life cycle of microorganisms,and may provide important insights into their adaptation to harsh and nutrient-depleted conditions.Although the underlying mechanisms have been well-studied in bacteria and yeasts (Herman,2002;Navarro Llorens et al.,2010),less is known about this growth phase in archaea yet.The haloarchaeon Haloferax mediterranei has served as a good model for studying haloarchaeal physiology and metabolism for several decades because of its accelerated growth,remarkable metabolic ability and genomic stability (Han et al.,2012).During stationary phase,H.mediterranei can produce halocin H4 (Cheung et al.,1997),synthesize gas vesicles (J(a)ger et al.,2002),secrete extracellular polysaccharide (Antón et al.,1988) and accumulate poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV)(Cai et al.,2012).Due to these specific features,we selected H.mediterranei as a model system to investigate the archaeal gene expression and regulation during the stationary phase.     

A unified framework for unraveling the functional interaction structure of a biomolecular network based on stimulus-response experimental data

We propose a unified framework for the identification of functional interaction structures of biomolecular networks in a way that leads to a new experimental design procedure. In developing our approach, we have built upon previous work. Thus we begin by pointing out some of the restrictions associated with existing structure identification methods and point out how these restrictions may be eased. In particular, existing methods use specific forms of experimental algebraic equations with which to identify the functional interaction structure of a biomolecular network. In our work, we employ an extended form of these experimental algebraic equations which, while retaining their merits, also overcome some of their disadvantages. Experimental data are required in order to estimate the coefficients of the experimental algebraic equation set associated with the structure identification task. However, experimentalists are rarely provided with guidance on which parameters to perturb, and to what extent, to perturb them. When a model of network dynamics is required then there is also the vexed question of sample rate and sample time selection to be resolved. Supplying some answers to these questions is the main motivation of this paper. The approach is based on stationary and/or temporal data obtained from parameter perturbations, and unifies the previous approaches of Kholodenko et al. (PNAS 99 (2002) 12841-12846) and Sontag et al. (Bioinformatics 20 (2004) 1877-1886). By way of demonstration, we apply our unified approach to a network model which cannot be properly identified by existing methods. Finally, we propose an experiment design methodology, which is not limited by the amount of parameter perturbations, and illustrate its use with an in numero example.     

Modelling Circulating Tumour Cells for Personalised Survival Prediction in Metastatic Breast Cancer

Ductal carcinoma is one of the most common cancers among women, and the main cause of death is the formation of metastases. The development of metastases is caused by cancer cells that migrate from the primary tumour site (the mammary duct) through the blood vessels and extravasating they initiate metastasis. Here, we propose a multi-compartment model which mimics the dynamics of tumoural cells in the mammary duct, in the circulatory system and in the bone. Through a branching process model, we describe the relation between the survival times and the four markers mainly involved in metastatic breast cancer (EPCAM, CD47, CD44 and MET). In particular, the model takes into account the gene expression profile of circulating tumour cells to predict personalised survival probability. We also include the administration of drugs as bisphosphonates, which reduce the formation of circulating tumour cells and their survival in the blood vessels, in order to analyse the dynamic changes induced by the therapy.We analyse the effects of circulating tumour cells on the progression of the disease providing a quantitative measure of the cell driver mutations needed for invading the bone tissue. Our model allows to design intervention scenarios that alter the patient-specific survival probability by modifying the populations of circulating tumour cells and it could be extended to other cancer metastasis dynamics.     

Charge transport in DNA model with vibrational and rotational coupling motions

The dynamics of the Peyrard-Bishop model for vibrational motion of DNA dynamics, which has been extended by taking into account the rotational motion for the nucleotides (Silva et al., J. Biol. Phys. 34, 511–519, 2018) is studied. We report on the presence of the modulational instability (MI) of a plane wave for charge migration in DNA and the generation of soliton-like excitations in DNA nucleotides. We show that the original differential-difference equation for the DNA dynamics can be reduced in the continuum approximation to a set of three coupled nonlinear equations. The linear stability analysis of continuous wave solutions of the coupled systems is performed and the growth rate of instability is found numerically. Numerical simulations show the validity of the analytical approach with the generation of wave packets provided that the wave numbers fall in the instability domain.     

Effects of rapid buffers on Ca2+ diffusion and Ca2+ oscillations.

Based on realistic mechanisms of Ca2+ buffering that include both stationary and mobile buffers, we derive and investigate models of Ca2+ diffusion in the presence of rapid buffers. We obtain a single transport equation for Ca2+ that contains the effects caused by both stationary and mobile buffers. For stationary buffers alone, we obtain an expression for the effective diffusion constant of Ca2+ that depends on local Ca2+ concentrations. Mobile buffers, such as fura-2, BAPTA, or small endogenous proteins, give rise to a transport equation that is no longer strictly diffusive. Calculations are presented to show that these effects can modify greatly the manner and rate at which Ca2+ diffuses in cells, and we compare these results with recent measurements by Allbritton et al. (1992). As a prelude to work on Ca2+ waves, we use a simplified version of our model of the activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane to illustrate the way in which Ca2+ buffering can affect both the amplitude and existence of Ca2+ oscillations.     

The kinetics of monospermic and polyspermic fertilization in free-spawning marine invertebrates

The equation of Vogel et al. (1982) is widely used in fertilization studies of free-spawning marine invertebrates to predict the percentage of viable eggs that will be fertilized at any specified levels of gamete concentration and contact time. Here, the random collision model that underlies the Vogel et al. equation is extended to distinguish between monospermic and polyspermic fertilization, and separate equations for the percentages of monospermic and polyspermic fertilization are obtained. These equations provide an explanation for empirical observations which have shown a decreased percentage of successful egg development at high sperm concentrations. Comparison is made with an earlier heuristic attempt (Styan, 1998) to predict the extent of polyspermic fertilization, and it is found that this earlier method can underestimate the percentage of polyspermic fertilization by up to 10 percent. Moreover, the approach used here retains the flexibility to model changes in sperm concentration due to dispersal mechanisms, and is able to model different mechanisms for the block to polyspermy.     

Effect of replacement of fish meal by meat and bone meal and poultry by-product meal in diets on the growth and immune response of Macrobrachium nipponense

The potential use of poultry by-product meal (PBM) and meat and bone meal (MBM) as alternative dietary protein sources for juvenile Macrobrachium nipponense was studied by a 70-day growth trial. Triplicate groups of M. nipponense (initial body weight: 0.37 g) were fed at 20.7-22.4 degrees C on each of the five isoenergetic and isonitrogenous diets (protein content about 38%) with different replacement of fish meal by MBM or PBM. The control diet used white fish meal as the sole protein source, the other four diets were prepared with 15% or 50% fish meal protein substituted by either MBM (MBM(15), MBM(50)) or PBM (PBM(15), PBM(50)). The results showed that replacement of fish meal by MBM in diets did not affect growth performance of M. nipponense (P > 0.05), while specific growth rate in PBM(15) was significantly higher than that in other groups (P < 0.05). Survival rates of shrimp fed with MBM(15) diet were significantly higher than that in other groups (P < 0.05). No significant differences in immunological parameters, including total haemocyte count (THC), phenoloxidase activity (PO) and respiratory burst (O(2)(-)), were observed between the shrimps that were fed five experimental diets, and all determined immunological parameters in control groups were slightly higher than those in replacement groups. In conclusion, either MBM or PBM investigated could replace up to 50% fish meal protein in diets for M. nipponense.     

Linkage and association studies of QTL for nuclear families by mixed models

The transmission disequilibrium test (TDT) has been utilized to test the linkage and association between a genetic trait locus and a marker. Spielman et al. (1993) introduced TDT to test linkage between a qualitative trait and a marker in the presence of association. In the presence of linkage, TDT can be applied to test for association for fine mapping (Martin et al., 1997; Spielman and Ewens, 1996). In recent years, extensive research has been carried out on the TDT between a quantitative trait and a marker locus (Allison, 1997; Fan et al., 2002; George et al., 1999; Rabinowitz, 1997; Xiong et al., 1998; Zhu and Elston, 2000, 2001). The original TDT for both qualitative and quantitative traits requires unrelated offspring of heterozygous parents for analysis, and much research has been carried out to extend it to fit for different settings. For nuclear families with multiple offspring, one approach is to treat each child independently for analysis. Obviously, this may not be a valid method since offspring of one family are related to each other. Another approach is to select one offspring randomly from each family for analysis. However, with this method much information may be lost. Martin et al. (1997, 2000) constructed useful statistical tests to analyse the data for qualitative traits. In this paper, we propose to use mixed models to analyse sample data of nuclear families with multiple offspring for quantitative traits according to the models in Amos (1994). The method uses data of all offspring by taking into account their trait mean and variance-covariance structures, which contain all the effects of major gene locus, polygenic loci and environment. A test statistic based on mixed models is shown to be more powerful than the test statistic proposed by George et al. (1999) under moderate disequilibrium for nuclear families. Moreover, it has higher power than the TDT statistic which is constructed by randomly choosing a single offspring from each nuclear family.     

Mathematical Model of Bone Remodeling Captures the Antiresorptive and Anabolic Actions of Various Therapies

A better understanding of the molecular pathways regulating the bone remodeling process should help in the development of new antiresorptive regulators and anabolic regulators, that is, regulators of bone resorption and of bone formation. Understanding the mechanisms by which parathyroid hormone (PTH) influences bone formation and how it switches from anabolic to catabolic action is important for treating osteoporosis (Poole and Reeve in Curr Opin Pharmacol 5:612–617, 2005). In this paper we describe a mathematical model of bone remodeling that incorporates, extends, and integrates several models of particular aspects of this biochemical system (Cabal et al. in J Bone Miner Res 28(8):1830–1836, 2013; Lemaire et al. in J Theor Biol 229:293–309, 2004; Peterson and Riggs in Bone 46:49–63, 2010; Raposo et al. in J Clin Endocrinol Metab 87(9):4330–4340, 2002; Ross et al. in J Disc Cont Dyn Sys Series B 17(6):2185–2200, 2012). We plan to use this model as a bone homeostasis platform to develop anabolic and antiresorptive compounds. The model will allow us to test hypotheses about the dynamics of compounds and to test the potential benefits of combination therapies. At the core of the model is the idealized account of osteoclast and osteoblast signaling given by Lemaire et al. (J Theor Biol 229:293–309, 2004). We have relaxed some of their assumptions about the roles of osteoprotegerin, transforming growth factor \(\upbeta \), and receptor activator of nuclear factor \(\upkappa \)B ligand; we have devised more detailed models of the interactions of these species. We have incorporated a model of the effect of calcium sensing receptor antagonists on remodeling (Cabal et al. in J Bone Miner Res 28(8):1830–1836, 2013). We have also incorporated a basic model of the effects of vitamin D on calcium homeostasis. We have included a simple model of the mechanism proposed by Bellido et al. (2003), Ross et al. (J Disc Cont Dyn Sys Series B 17(6):2185–2200, 2012), of the influence of PTH on osteoblast apoptosis, a mechanism that accounts for the anabolic response to pulsatile PTH administration. Finally, we have devised a simple model of the administration and effects of bisphosphonates. The biomarkers in the model are procollagen type 1 amino-terminal propeptide and C-terminal telopeptide. Bone mineral density is the model’s principal endpoint.     

Multilocus linkage tests based on affected relative pairs

For complex diseases, recent interest has focused on methods that take into account joint effects at interacting loci. Conditioning on effects of disease loci at known locations can lead to increased power to detect effects at other loci. Moreover, use of joint models allows investigation of the etiologic mechanisms that may be involved in the disease. Here we present a method for simultaneous analysis of the joint genetic effects at several loci that uses affected relative pairs. The method is a generalization of the two-locus LOD-score analysis for affected sib pairs proposed by Cordell et al. We derive expressions for the relative risk, lambdaR, to a relative of an affected individual, in terms of the additive and epistatic components of variance at an arbitrary number of disease loci, and we show how these can be used to fit a likelihood model to the identity-by-descent sharing among pairs of affected relatives in extended pedigrees. We implement the method by use of a stepwise strategy in which, given evidence of linkage to disease at m-1 locations on the genome, we calculate the conditional likelihood curve across the genome for an mth disease locus, using multipoint methods similar to those proposed by Kruglyak et al. We evaluate the properties of our method by use of simulated data and present an application to real data from families with insulin-dependent diabetes mellitus.     

Relations Between Bacterial Biomass and Carbon Cycle in Marine Sediments: An Early Diagenetic Model

A new model for early diagenetic processes has been developed through a new formula explicitly accounting for microbial population dynamics. Following a mechanistic approach based on enzymatic reactions, a new model has been proposed for oxic mineralisation and denitrification. It incorporates the dynamics of bacterial metabolism. We find a general formula for inhibition processes of which some other mathematical expressions are particular cases. Moreover a fast numerical algorithm has been developed. It allows us to perform simulations of different diagenetic models in non-steady states. We use this algorithm to compare our model to a classical one (Soetaert et al., 1996). Dynamical evolutions of a perturbation of particulate organic carbon (POC) input are studied for both models. The results are very similar for stationary cases. But with variable inputs, the bacterial biomass dynamics brings about noticeable differences, and these are discussed.     

Open loop optogenetic control of simulated cortical epileptiform activity

We present a model for the use of open loop optogenetic control to inhibit epileptiform activity in a meso scale model of the human cortex. The meso scale cortical model first developed by Liley et al. (2001) is extended to two dimensions and the nature of the seizure waves is studied. We adapt to the meso scale a 4 state functional model of Channelrhodopsin-2 (ChR2) ion channels. The effects of pulsed and constant illumination on the conductance of these ion channels is presented. The inhibitory cell population is targeted for the application of open loop control. Seizure waves are successfully suppressed and the inherent properties of the optogenetic channels ensures charge balance in the cortex, protecting it from damage.     

Intrinsic binding effects in mixed counterionic polyelectrolyte systems: extension of condensation theory and comparison with voltammetric data

Voltammetric speciation data for the potassium/zinc/polymethacrylate system, recently obtained for various charge densities of the polyelectrolyte (Díaz-Cruz et al., Anal. Chim. Acta, 264 (1992) 163) and for different concentrations of monovalent counterions (van den Hoop and van Leeuwen, Anal. Chim. Acta, 273 (1993) 275), are compared with theoretical predictions computed according to a new thermodynamic model developed by Paoletti et al. (Biophys. Chem., 41 (1991) 73) and recently extended by Benegas and Paoletti (in preparation). The model allows: (i) the simultaneous condensation of both monovalent and divalent counterions and (ii) can account for a certain specific affinity of the polyelectrolyte for one type of the counterion over the other. For various charge densities of the polyelectrolyte, experimentally obtained speciation data for the K/Zn/PMA system agree well with theoretical predictions by considering an extra reduced molar affinity energy of -4RT for the Zn(2+) polyelectrolyte binding. The agreement between experimental and theoretical values for the distribution of Zn(2+) ions over the free and bound state becomes less perfect for relatively high concentrations of monovalent counterions.     

Density functional theory for the selective adsorption of small molecules on a surface modified with polymer brushes

A density functional method based on weighted density approximation is extended to study the selective adsorption of small molecules on a surface modified with end-grafted square-well chains. The excess part of the Helmholtz free energy functional is divided into two components: the hard sphere repulsion and the square-well attraction. The equation of state for hard sphere chain fluids developed by Liu et al. is used to calculate the repulsive part of the excess Helmholtz free energy functional, and the equation of state for square-well chain fluid with variable range developed by Li et al. is employed to calculate the attractive part. With this theoretical model, we examine the physical properties of the grafted polymer and the selective adsorption of small molecules on the modified surface.     

The time for oxygen release in photosynthesis: reconciliation of flash polarography with other measurement techniques

The time for oxygen release in photosynthesis has been reported to be 30–130 ms when measured by flash polarography under low polarization voltages (Plijter et al. 1988), in opposition to 1–3 ms with light modulated oxygen polarography (Jolio et al. 1966), with the detection of produced oxygen in a flowing sample (Etienne 1968) or with photoacoustic detection of oxygen evolution (Canaani et al. 1988). However, we show here that flash polarographic measurements require properly cleaned electrodes, a precise polarization voltage, as well as a short polarization time of the electrodes. When these criteria were met, an oxygen release in less than 2 ms could be measured by flash polarography under low polarization voltages, in accordance with the other techniques. But under high polarization voltages, the interpretation of the polarographic response to oxygen production must take into account the diffusion of oxygen, the capacitance of the platinum electrode and the oxygen release time. We present a model of the electrode response taking into account these factors; by interpreting the response of the electrodes with this model, we found an oxygen release time of 1.7 ms. These evidences support strongly a short oxygen release time of 1–3 ms.