Effect of disease-selective predation on prey infected by contact and external sources

We propose and analyze a simple mathematical model for susceptible prey (S)–infected prey (I)–predator (P) interaction, where the susceptible prey population (S) is infected directly from external sources as well as through contact with infected class (I) and the predator completely avoids consuming the infected prey. The model is analyzed to obtain different thresholds of the key parameters under which the system exhibits stability around the biologically feasible equilibria. Through numerical simulations we display the effects of external infection and the infection through contact on the system dynamics in the absence as well as in the presence of the predator. We compare the system dynamics when infection occurs only through contact, with that when it occurs through contact and external sources. Our analysis demonstrates that under a disease-selective predation, stability and oscillations of the system is determined by two key parameters: the external infection rate and the force of infection through contact. Due to the introduction of external infection, the predator and the prey population show limit-cycle oscillations over a range parametric values. We suggest that while predicting the dynamics of such an eco-epidemiological system, the modes of infection and the infection rates might be carefully investigated.

Arenavirus budding resulting from viral-protein-associated cell membrane curvature

Viral replication occurs within cells, with release (and onward infection) primarily achieved through two alternative mechanisms: lysis, in which virions emerge as the infected cell dies and bursts open; or budding, in which virions emerge gradually from a still living cell by appropriating a small part of the cell membrane. Virus budding is a poorly understood process that challenges current models of vesicle formation. Here, a plausible mechanism for arenavirus budding is presented, building on recent evidence that viral proteins embed in the inner lipid layer of the cell membrane. Experimental results confirm that viral protein is associated with increased membrane curvature, whereas a mathematical model is used to show that localized increases in curvature alone are sufficient to generate viral buds. The magnitude of the protein-induced curvature is calculated from the size of the amphipathic region hypothetically removed from the inner membrane as a result of translation, with a change in membrane stiffness estimated from observed differences in virion deformation as a result of protein depletion. Numerical results are based on experimental data and estimates for three arenaviruses, but the mechanisms described are more broadly applicable. The hypothesized mechanism is shown to be sufficient to generate spontaneous budding that matches well both qualitatively and quantitatively with experimental observations.

Modeling the summertime evolution of sea-ice melt ponds

1] We present a mathematical model describing the summer melting of sea ice. We simulate the evolution of melt ponds and determine area coverage and total surface ablation. The model predictions are tested for sensitivity to the melt rate of unponded ice, enhanced melt rate beneath the melt ponds, vertical seepage, and horizontal permeability. The model is initialized with surface topographies derived from laser altimetry corresponding to first-year sea ice and multiyear sea ice. We predict that there are large differences in the depth of melt ponds and the area of coverage between the two types of ice. We also find that the vertical seepage rate and the melt rate of unponded ice are important in determining the total surface ablation and area covered by melt ponds.

Multiple stationary solutions of an irradiated slab

A mathematical model describing the heat budget of an irradiated medium is introduced. The one-dimensional form of the equations and boundary conditions are presented and analysed. Heat transport at one face of the slab occurs by absorption (and reflection) of an incoming beam of short-wave radiation with a fraction of this radiation penetrating into the body of the slab, a diffusive heat flux in the slab and a prescribed incoming heat flux term. The other face of the slab is immersed in its own melt and is considered to be a free surface. Here, temperature continuity is prescribed and evolution of the surface is determined by a Stefan condition. These boundary conditions are flexible enough to describe a range of situations such as a laser shining on an opaque medium, or the natural environment of polar sea ice or lake ice. A two-stream radiation model is used which replaces the simple Beer’s law of radiation attenuation frequently used for semi-infinite domains. The stationary solutions of the governing equations are sought and it is found that there exists two possible stationary solutions for a given set of boundary conditions and a range of parameter choices. It is found that the existence of two stationary solutions is a direct result of the model of radiation absorption, due to its effect on the albedo of the medium. A linear stability analysis and numerical calculations indicate that where two stationary solutions exist, the solution corresponding to a larger thickness is always stable and the solution corresponding to a smaller thickness is unstable. Numerical simulations reveal that when there are two solutions, if the slab is thinner than the smaller stationary thickness it will melt completely, whereas if the slab is thicker than the smaller stationary thickness it will evolve toward the larger stationary thickness. These results indicate that other mechanisms (e.g. wave-induced agglomeration of crystals) are necessary to grow a slab from zero initial thickness in the parameter regime that yields two stationary solutions.

Analytical and numerical solutions describing the inward solidification of a binary melt

We present a mathematical model describing the inward solidification of a slab, a circular cylinder and a sphere of binary melt kept below its equilibrium freezing temperature. The thermal and physical properties of the melt and solid are assumed to be identical. An asymptotic method, valid in the limit of large Stefan number is used to decompose the moving boundary problem for a pure substance into a hierarchy of fixed-domain diffusion problems. Approximate, analytical solutions are derived for the inward solidification of a slab and a sphere of a binary melt which are compared with numerical solutions of the unapproximated system. The solutions are found to agree within the appropriate asymptotic regime of large Stefan number and small time. Numerical solutions are used to demonstrate the dependence of the solidification process upon the level of impurity and other parameters. We conclude with a discussion of the solutions obtained, their stability and possible extensions and refinements of our study.

Modelling negative feedback networks for activating transcription factor 3 predicts a dominant role for miRNAs in immediate early gene regulation

Activating transcription factor 3 (Atf3) is rapidly and transiently upregulated in numerous systems, and is associated with various disease states. Atf3 is required for negative feedback regulation of other genes, but is itself subject to negative feedback regulation possibly by autorepression. In cardiomyocytes, Atf3 and Egr1 mRNAs are upregulated via ERK1/2 signalling and Atf3 suppresses Egr1 expression. We previously developed a mathematical model for the Atf3-Egr1 system. Here, we adjusted and extended the model to explore mechanisms of Atf3 feedback regulation. Introduction of an autorepressive loop for Atf3 tuned down its expression and inhibition of Egr1 was lost, demonstrating that negative feedback regulation of Atf3 by Atf3 itself is implausible in this context. Experimentally, signals downstream from ERK1/2 suppress Atf3 expression. Mathematical modelling indicated that this cannot occur by phosphorylation of pre-existing inhibitory transcriptional regulators because the time delay is too short. De novo synthesis of an inhibitory transcription factor (ITF) with a high affinity for the Atf3 promoter could suppress Atf3 expression, but (as with the Atf3 autorepression loop) inhibition of Egr1 was lost. Developing the model to include newly-synthesised miRNAs very efficiently terminated Atf3 protein expression and, with a 4-fold increase in the rate of degradation of mRNA from the mRNA/miRNA complex, profiles for Atf3 mRNA, Atf3 protein and Egr1 mRNA approximated to the experimental data. Combining the ITF model with that of the miRNA did not improve the profiles suggesting that miRNAs are likely to play a dominant role in switching off Atf3 expression post-induction.

Modeling migratory grazing of zooplankton on toxic and non-toxic phytoplankton

Migratory grazing of zooplankton between non-toxic phytoplankton (NTP) and toxic phytoplankton (TPP) is a realistic phenomena unexplored so far. The present article is a first step in this direction. A mathematical model of NTP–TPP-zooplankton with constant and variable zooplankton migration is proposed and analyzed. The asymptotic dynamics of the model system around the biologically feasible equilibria is explored through local stability analysis. The dynamics of the proposed system is explored and displayed for different combination of migratory parameters and toxin inhibition parameters. Our analysis suggests that the migratory grazing of zooplankton has a significant role in determining the dynamic stability and oscillation of phytoplankton zooplankton systems.

Toxin-allelopathy among phytoplankton species prevents competitive exclusion

Toxic or allelopathic compounds liberated by toxin-producing phytoplankton (TPP) acts as a strong mediator in plankton dynamics. On an analysis of a set of phytoplankton biomass data that have been collected by our group in the northwest part of the Bay of Bengal, and by analysis of a three-component mathematical model under a constant as well as a stochastic environment, we explore the role of toxin-allelopathy in determining the dynamic behavior of the competing phytoplankton species. The overall results, based on analytical and numerical wings, demonstrate that toxin-allelopathy due to the TPP promotes a stable co-existence of those competitive phytoplankton that would otherwise exhibit competitive exclusion of the weak species. Our study suggests that TPP might be a potential candidate for maintaining the co-existence and diversity of competing phytoplankton species.

Competing effects of toxin-producing phytoplankton on overall plankton populations in the Bay of Bengal

The coexistence of a large number of phytoplankton species on a seemingly limited variety of resources is a classical problem in ecology, known as ‘the paradox of the plankton’. Strong fluctuations in species abundance due to the external factors or competitive interactions leading to oscillations, chaos and short-term equilibria have been cited so far to explain multi-species coexistence and biodiversity of phytoplankton. However, none of the explanations has been universally accepted. The qualitative view and statistical analysis of our field data establish two distinct roles of toxin-producing phytoplankton (TPP): toxin allelopathy weakens the interspecific competition among phytoplankton groups and the inhibition due to ingestion of toxic substances reduces the abundance of the grazer zooplankton. Structuring the overall plankton population as a combination of nontoxic phytoplankton (NTP), toxic phytoplankton, and zooplankton, here we offer a novel solution to the plankton paradox governed by the activity of TPP. We demonstrate our findings through qualitative analysis of our sample data followed by analysis of a mathematical model.

Bond graph models of DC-DC converters operating in both CCM and DCM

In this paper, Bond Graphs are employed to develop a novel mathematical model of conventional switched-mode DC-DC converters valid for both continuous and discontinuous conduction modes. A unique causality bond graph model of hybrid models is suggested with the operation of the switch and the diode to be represented by a Modulated Transformer with a binary input and a resistor with fixed conductance causality. The operation of the diode is controlled using an if-then function within the model. The extracted hybrid model is implemented on a Boost and Buck converter with their operations to change from CCM to DCM and to return to CCM. The vector fields of the models show validity in a wide operation area and comparison with the simulation of the converters using PSPICE reveals high accuracy of the proposed model, with the Normalised Root Means Square Error and the Maximum Absolute Error remaining adequately low. The model is also experimentally tested on a Buck topology.

Regulation of early steps of GPVI signal transduction by phosphatases: a systems biology approach

We present a data-driven mathematical model of a key initiating step in platelet activation, a central process in the prevention of bleeding following Injury. In vascular disease, this process is activated inappropriately and causes thrombosis, heart attacks and stroke. The collagen receptor GPVI is the primary trigger for platelet activation at sites of injury. Understanding the complex molecular mechanisms initiated by this receptor is important for development of more effective antithrombotic medicines. In this work we developed a series of nonlinear ordinary differential equation models that are direct representations of biological hypotheses surrounding the initial steps in GPVI-stimulated signal transduction. At each stage model simulations were compared to our own quantitative, high-temporal experimental data that guides further experimental design, data collection and model refinement. Much is known about the linear forward reactions within platelet signalling pathways but knowledge of the roles of putative reverse reactions are poorly understood. An initial model, that includes a simple constitutively active phosphatase, was unable to explain experimental data. Model revisions, incorporating a complex pathway of interactions (and specifically the phosphatase TULA-2), provided a good description of the experimental data both based on observations of phosphorylation in samples from one donor and in those of a wider population. Our model was used to investigate the levels of proteins involved in regulating the pathway and the effect of low GPVI levels that have been associated with disease. Results indicate a clear separation in healthy and GPVI deficient states in respect of the signalling cascade dynamics associated with Syk tyrosine phosphorylation and activation. Our approach reveals the central importance of this negative feedback pathway that results in the temporal regulation of a specific class of protein tyrosine phosphatases in controlling the rate, and therefore extent, of GPVI-stimulated platelet activation.

Body mass index values for newborns according to gestational age

Objective: The combination of twho anthropometric parameters has been more appropriate to assess body composition and proportions in children, with special attention to the Body Mass Index (BMI), as it relates weight and length. However the BMI values for the neonatal period have not been determined yet. This study shows the BMI for newborns at different gestational ages represented in a normal smoothed percentile curve. Methods: Retrospective study including 2,406 appropriate for gestational age newborns following the Alexander et al curve (1996) from 29 to 42 weeks of gestational age. Weight and lenght were measured following standard procedures. For the construction aof a normal smoothed percentile curve, the 3(rd) 5(th), 10(th), 25(th), 5(th), 75(th), 90(th) and 95(th) percentiles were determined and a statistical procedure based on the mathematical model ""sinosuoidal fit"" was applied to establish a curve that estimates biological growth parameters. Results: The Body Mass Index values for gestational age in all percentiles shows a steady increase up to 38 weeks, levels off up to the 40(th) week, followed by a slight decrease to the 42(nd) week in both genders. Conclusion: The results show a direct correlation between gestational age and Body Mass Index for both genders in the nine percentiles, and can provide a useful reference to assess intra-uterine proportional growth.

DC traction load flow including AC distribution network

In the first part some information and characterisation about an AC distribution network that feeds traction substations and their possible influences on the DC traction load flow are presented. Those influences are investigated and mathematically modelled. To corroborate the mathematical model, an example is presented and their results are confronted with real measurements.

Dissociation of circadian and light inhibition of melatonin release through forced desynchronization in the rat

Pineal melatonin release exhibits a circadian rhythm with a tight nocturnal pattern. Melatonin synthesis is regulated by the master circadian clock within the hypothalamic suprachiasmatic nucleus (SCN) and is also directly inhibited by light. The SCN is necessary for both circadian regulation and light inhibition of melatonin synthesis and thus it has been difficult to isolate these two regulatory limbs to define the output pathways by which the SCN conveys circadian and light phase information to the pineal. A 22-h light-dark (LD) cycle forced desynchrony protocol leads to the stable dissociation of rhythmic clock gene expression within the ventrolateral SCN (vlSCN) and the dorsomedial SCN (dmSCN). In the present study, we have used this protocol to assess the pattern of melatonin release under forced desynchronization of these SCN subregions. In light of our reported patterns of clock gene expression in the forced desynchronized rat, we propose that the vlSCN oscillator entrains to the 22-h LD cycle whereas the dmSCN shows relative coordination to the light-entrained vlSCN, and that this dual-oscillator configuration accounts for the pattern of melatonin release. We present a simple mathematical model in which the relative coordination of a single oscillator within the dmSCN to a single light-entrained oscillator within the vlSCN faithfully portrays the circadian phase, duration and amplitude of melatonin release under forced desynchronization. Our results underscore the importance of the SCN`s subregional organization to both photic input processing and rhythmic output control.

Molecular characterization and fumonisin production by Fusarium verticillioides isolated from corn grains of different geographic origins in Brazil

Gibberella moniliformis is most commonly associated with maize worldwide and produces high levels of fumonisins, some of the most agriculturally important mycotoxins. Studies demonstrate that molecular methods can be helpful for a rapid identification of Fusarium species and their levels of toxin production. The purpose of this research was to apply molecular methods (AFLP, TEF-1 alpha partial gene sequencing and PCR based on MAT alleles) for the identification of Fusarium species isolated from Brazilian corn and to verify if real time RT-PCR technique based on FUM1 and FUM19 genes is appropriated to estimate fumonisins B(1) and B(2) production levels. Among the isolated strains, 96 were identified as Fusarium verricillioides, and four as other Fusarium species. Concordant phylogenies were obtained by AFLP and TEF-1 alpha sequencing, permitting the classification of the different species into distinct clades. Concerning MAT alleles, 70% of the F. verricillioides isolates carried the MAT-1 and 30% MAT-2. A significant correlation was observed between the expression of the genes and toxin production r=0.95 and r=0.79 (correlation of FUM1 with FB(1) and FB(2), respectively, P < 0.0001): r=0.93 and r =0.78 (correlation of FUM19 with FB(1) and FB(2). respectively, P < 0.0001). Molecular methods used in this study were found to be useful for the rapid identification of Fusarium species. The high and significant correlation between FUM1 and FUM19 expression and fumonisins production suggests that real time RT-PCR is suitable for studies considering the influence of abiotic and biotic factors on expression of these genes. This is the first report concerning the expression of fumonisin biosynthetic genes in Fusarium strains isolated from Brazilian agricultural commodity. (c) 2010 Elsevier B.V. All rights reserved.