Inma Rodríguez Cantalapiedra

Modeling of spatiotemporal patterns in bacterial colonies.

A diffusion-reaction model for the growth of bacterial colonies is presented. The often observed cooperative behavior developed by bacteria which increases their motility in adverse growth conditions is here introduced as a nonlinear diffusion term. The presence of this mechanism depends on a response which can present hysteresis. By changing only the concentrations of agar and initial nutrient, numerical integration of the proposed model reproduces the different patterns shown by Bacillus subtilis OG-01.

Dependency of calcium alternans on Ryanodine Receptor Refractoriness

Background Rapid pacing rates induce alternations in the cytosolic calcium concentration caused by fluctuations in calcium released from the sarcoplasmic reticulum (SR). However, the relationship between calcium alternans and refractoriness of the SR calcium release channel (RyR2) remains elusive. Methodology/Principal Findings To investigate how ryanodine receptor (RyR2) refractoriness modulates calcium handling on a beat-to-beat basis using a numerical rabbit cardiomyocyte model. We used a mathematical rabbit cardiomyocyte model to study the beat-to-beat calcium response as a function of RyR2 activation and inactivation. Bi-dimensional maps were constructed depicting the beat-to-beat response. When alternans was observed, a novel numerical clamping protocol was used to determine whether alternans was caused by oscillations in SR calcium loading or by RyR2 refractoriness. Using this protocol, we identified regions of RyR2 gating parameters where SR calcium loading or RyR2 refractoriness underlie the induction of calcium alternans, and we found that at the onset of alternans both mechanisms contribute. At low inactivation rates of the RyR2, calcium alternans was caused by alternation in SR calcium loading, while at low activation rates it was caused by alternation in the level of available RyR2s. Conclusions/Significance We have mapped cardiomyocyte beat-to-beat responses as a function of RyR2 activation and inactivation, identifying domains where SR calcium load or RyR2 refractoriness underlie the induction of calcium alternans. A corollary of this work is that RyR2 refractoriness due to slow recovery from inactivation can be the cause of calcium alternans even when alternation in SR calcium load is present.

Are SR Ca content fluctuations or SR refractoriness the key to atrial cardiac alternans?: insights from a human atrial model

Despite the important role of electromechanical alternans in cardiac arrhythmogenesis, its molecular origin is not well understood. The appearance of calcium alternans has often been associated to fluctuations in the sarcoplasmic reticulum (SR) Ca loading. However, cytosolic calcium alternans observed without concurrent oscillations in the SR Ca content suggests an alternative mechanism related to a dysfunction in the dynamics of the ryanodine receptor (RyR2). We have investigated the effect of SR release refractoriness in the appearance of alternans, using a mathematical model of a single human atrial cell, based on the model by Nygren et al. (30), where we modified the dynamics of the RyR2 and of SR Ca release. The genesis of calcium alternans was studied stimulating the cell for different periods and values of the RyR2 recovery time from inactivation. At fast rates cytosolic calcium alternans were obtained without concurrent SR Ca content fluctuations. A transition from regular response to alternans was also observed, changing the recovery time from inactivation of the RyR2. This transition was found to be hysteretic, so for a given set of parameters different responses were observed. We then studied the relevance of RyR2 refractoriness for the generation of alternans, reproducing the same protocols as in recent experiments. In particular, restitution of Ca release during alternans was studied with a S1S2 protocol, obtaining a different response if the S2 stimulation was given after a long or a short release. We show that the experimental results can be explained by RyR2 refractoriness, arising from a slow RyR2 recovery from inactivation, stressing the role of the RyR2 in the genesis of alternans.

Chaotic motion of space charge wave fronts in semiconductors under time-independent voltage bias

A standard drift-diffusion model of space charge wave propagation in semiconductors has been studied numerically and analytically under dc voltage bias. For sufficiently long samples, appropriate contact resistivity, and applied voltage-such that the sample is biased in a regime of negative differential resistance-we find chaos in the propagation of nonlinear fronts (charge monopoles of alternating sign) of electric field. The chaos is always low dimensional, but has a complex spatial structure; this behavior can be interpreted using a finite-dimensional asymptotic model in which the front (charge monopole) positions and the electrical current are the only dynamical variables.

Asymptotic analysis of the Gunn effect with realistic boundary conditions

A general asymptotic analysis of the Gunn effect in n-type GaAs under general boundary conditions for metal-semiconductor contacts is presented. Depending on the parameter values in the boundary condition of the injecting contact, different types of waves mediate the Gunn effect. The periodic current oscillation typical of the Gunn effect may be caused by moving charge-monopole accumulation or depletion layers, or by low- or high-field charge-dipole solitary waves. A new instability caused by multiple shedding of ~low-field! dipole waves is found. In all cases the shape of the current oscillation is described in detail: we show the direct relationship between its major features ~maxima, minima, plateaus, etc.! and several critical currents ~which depend on the values of the contact parameters!. Our results open the possibility of measuring contact parameters from the analysis of the shape of the current oscillation. @S1063-651X~97!04708-9#

Reexcitation mechanisms in epicardial tissue: Role of Ito density heterogeneities and INa inactivation kinetics

Dispersion of action potential repolarization is known to be an important arrhythmogenic factor in cardiopathies such as Brugada syndrome. In this work, we analyze the effect of a variation in sodium current (I(Na)) inactivation and a heterogeneous rise of transient outward current (I(to)) in the probability of reentry in epicardial tissue. We use the Luo-Rudy model of epicardial ventricular action potential to study wave propagation in a one-dimensional fiber. Spatial dispersion in repolarization is introduced by splitting the fiber into zones with different strength of I(to). We then analyze the pro-arrhythmic effect of a variation in the relaxation time and steady-state of the sodium channel fast inactivating gate h. We quantify the probability of reentry measuring the percentage of reexcitations that occurs in 200 beats. We find that, for high stimulation rates, this percentage is negligible, but increases notably for pacing periods above 700ms. Surprisingly, with decreasing I(Na) inactivation time, the percentage of reexcitations does not grow monotonically, but presents vulnerable windows, separated by values of the I(Na) inactivation speed-up where reexcitation does not occur. By increasing the strength of L-type calcium current I(CaL) above a certain threshold, reexcitation disappears. Finally, we show the formation of reentry in stimulated two-dimensional epicardial tissue with modified I(Na) kinetics and I(to) heterogeneity. Thus, we confirm that while I(to) dispersion is necessary for phase-2 reentry, altered sodium inactivation kinetics influences the probability of reexcitation in a highly nonlinear fashion.

Semiclassical theory of shot noise in ballistic n ¿ -i-n ¿ semiconductor structures: Relevance of Pauli and long-range Coulomb correlations

We work out a semiclassical theory of shot noise in ballistic

Universality of the Gunn effect: Self-sustained oscillations mediated by solitary waves

{n}^{+}\ensuremath{-}i\ensuremath{-}{n}^{+}

Cardiac dynamics: a simplified model for action potential propagation

semiconductor structures aiming at studying two fundamental physical correlations coming from Pauli exclusion principle and long-range Coulomb interaction. The theory provides a unifying scheme which, in addition to the current-voltage characteristics, describes the suppression of shot noise due to Pauli and Coulomb correlations in the whole range of system parameters and applied bias. The whole scenario is summarized by a phase diagram in the plane of two dimensionless variables related to the sample length and contact chemical potential. Here different regions of physical interest can be identified where only Coulomb or only Pauli correlations are active, or where both are present with different relevance. The predictions of the theory are proven to be fully corroborated by Monte Carlo simulations.

Stationary states and phase diagram for a model of the Gunn effect under realistic boundary conditions

The Gunn effect consists of time-periodic oscillations of the current flowing through an external purely resistive circuit mediated by solitary wave dynamics of the electric field on an attached appropriate semiconductor. By means of an asymptotic analysis, it is argued that Gunn-like behavior occurs in specific classes of model equations. As an illustration, an example related to the constrained Cahn-Allen equation is analyzed.