Effect of beta-adrenergic blockade on dynamic electrical restitution in vivo

The slope of the action potential duration (APD) restitution curve may be a significant determinant of the propensity to develop ventricular fibrillation, with steeper slopes associated with a more arrhythmogenic substrate. We hypothesized that one mechanism by which beta-blockers reduce sudden cardiac death is by flattening the APD restitution curve. Therefore, we investigated whether infusion of esmolol modulates the APD restitution curve in vivo. In 10 Yorkshire pigs, dynamic APD restitution curves were determined from measurements of APD at 90% repolarization with a monophasic action potential catheter positioned against the right ventricular septum during right ventricular apical pacing in the basal state and during infusion of esmolol. APD restitution curves were fitted to the three-parameter (a, b, c) exponential equation, APD = a.[1 - e((-b.DI))] + c, where DI is the diastolic interval. Esmolol decreased the maximal APD slope, 0.68 +/- 0.14 vs. 0.94 +/- 0.24 (baseline), P = 0.002, and flattened the APD restitution curve at shorter DIs, 75 and 100 ms (P < 0.05). To compare the slopes of the APD restitution curves at similar steady states, slopes were also computed at points of intersection between the restitution curve and the lines representing pacing at a fixed cycle length (CL) of 200, 225, 250, 275, and 300 ms using the relationship CL = APD + DI. Esmolol decreased APD restitution slopes at CLs 200-275 ms (P < 0.05). Esmolol flattens the cardiac APD restitution curve in vivo, particularly at shorter CLs and DIs. This may represent a novel mechanism by which beta-blockers prevent sudden cardiac death.     

Alternans Resonance and Propagation Block during Supernormal Conduction in Cardiac Tissue with Decreased [K]o

Cardiac restitution is an important factor in arrhythmogenesis. Steep positive action potential duration and conduction velocity (CV) restitution slopes promote alternans and reentrant arrhythmias. We examined the consequences of supernormal conduction (characterized by a negative CV restitution slope) on patterns of conduction and alternans in strands of Luo-Rudy model cells and in cultured cardiac cell strands. Interbeat intervals (IBIs) were analyzed as a function of distance during S1S2 protocols and during pacing at alternating cycle lengths. Supernormal conduction was induced by decreasing [K⁺]_o. In control [K⁺]_o simulations, S1S2 intervals converged toward basic cycle length with a length constant determined by both CV and the CV restitution slope. During alternant pacing, the amplitude of IBI alternans converged with a shorter length constant, determined also by the action potential duration restitution slope. In contrast, during supernormal conduction, S1S2 intervals and the amplitude of alternans diverged. This amplification (resonance) led to phase-locked or more complex alternans patterns, and then to distal conduction block. The convergence/divergence of IBIs was verified in the cultured strands, in which naturally occurring tissue heterogeneities resulted in prominent discontinuities of the spatial IBI profiles. We conclude that supernormal conduction potentiates alternans and spatial analysis of IBIs represents a powerful method to locate tissue heterogeneities.     

Dynamics of action potential head-tail interaction during reentry in cardiac tissue: ionic mechanisms

In a sufficiently short reentry pathway, the excitation wave front (head) propagates into tissue that is partially refractory (tail) from the previous action potential (AP). We incorporate a detailed mathematical model of the ventricular myocyte into a one-dimensional closed pathway to investigate the effects of head-tail interaction and ion accumulation on the dynamics of reentry. The results were the following: 1) a high degree of head-tail interaction produces oscillations in several AP properties; 2) Ca(2+)-transient oscillations are in phase with AP duration oscillations and are often of greater magnitude; 3) as the wave front propagates around the pathway, AP properties undergo periodic spatial oscillations that produce complicated temporal oscillations at a single site; 4) depending on the degree of head-tail interaction, intracellular [Na(+)] accumulation during reentry either stabilizes or destabilizes reentry; and 5) elevated extracellular [K(+)] destabilizes reentry by prolonging the tail of postrepolarization refractoriness.     

The pinwheel experiment revisited: effects of cellular electrophysiological properties on vulnerability to cardiac reentry

In normal heart, ventricular fibrillation can be induced by a single properly timed strong electrical or mechanical stimulus. A mechanism first proposed by Winfree and coined the "pinwheel experiment" emphasizes the timing and strength of the stimulus in inducing figure-of-eight reentry. However, the effects of cellular electrophysiological properties on vulnerability to reentry in the pinwheel scenario have not been investigated. In this study, we extend Winfree's pinwheel experiment to show how the vulnerability to reentry is affected by the graded action potential responses induced by a strong premature stimulus, action potential duration (APD), and APD restitution in simulated monodomain homogeneous two-dimensional tissue. We find that a larger graded response, longer APD, or steeper APD restitution slope reduces the vulnerable window of reentry. Strong graded responses and long APD promote tip-tip interactions at long coupling intervals, causing the two initiated spiral wave tips to annihilate. Steep APD restitution promotes wave front-wave back interaction, causing conduction block in the central common pathway of figure-of-eight reentry. We derive an analytical treatment that shows good agreement with numerical simulation results.     

Relationship of length of uterus in prepubertal pigs and number of corpora lutea and fetuses at 30 days of gestation

Relationships of the length of the uterus at one reproductive stage to the length at other stages and the effect on potential litter size were determined. In Experiment 1, the length of the uterus was measured in situ at 20, 60, or 100 days of age at laparotomy with 20 gilts in each of the three age groups. Forty days after the initial measurement, the uterus was again measured in situ, gilts were ovariectomized and hysterectomized, and associations among uterine measurements at the two different stages were determined. Correlations between uterine length in situ and between uterine length and weight 40 days later were all greater than 0.75 (P < 0.001). The length of one uterine horn increased from 13.9 cm at 20 days of age to 36.7 cm at 140 days of age (P < 0.01). In Experiment 2, 66 gilts were unilaterally hysterectomized and ovariectomized (UHOX) at 150 days of age and the ovary was weighed. Length of the one horn was measured in 36 gilts. At 10 days after first estrus, the length of the remaining uterine horn was measured at laparotomy and corpora lutea were counted in 53 gilts. In 15 of the 53 gilts the remaining uterine horn was removed to obtain uterine weights. At the second or third estrus, 38 gilts were mated and at Day 30 of gestation, 31 gilts were pregnant. The gilts were killed, and the length of the uterus measured and corpora lutea (CL) and fetuses were counted. The length of one uterine horn at 150 days of age was 70 cm with a range of 47–110 cm. At 10 days after first estrus, length had increased to a mean of 141 cm with a range of 86–194 cm and at 30 days of gestation the mean was 244 cm with a range of 186–311 cm. There were 12.2 CL at first estrus, which was not different from 12.4 CL at the second or third estrus. The mean number of fetuses in one horn at 30 days of gestation was 9.5 with 77% prenatal survival. Length of uterine horn at 150 days of age was correlated with uterine horn length (r = 0.56, P < 0.001) at 10 days after first estrus and number of live fetuses (r=0.39, P<0.05) at 30 days of gestation. At Day 10 after first estrus, uterine length was not correlated with the number of CL, whereas at Day 30 of gestation, the number of CL and uterine length were correlated with the number of live fetuses in those gilts with below the mean number of live fetuses, but not in those gilts with above the mean of live fetuses. The number of live fetuses (r=0.66, P < 0.001) and fetal survival (r=0.63; P < 0.001) were correlated with uterine horn length in pregnant UHOX gilts. Length of the prepubertal uterus gives an indication of postpubertal length and the potential litter size in pigs.     

The Transfer Functions of Cardiac Tissue during Stochastic Pacing

The restitution properties of cardiac action potential duration (APD) and conduction velocity (CV) are important factors in arrhythmogenesis. They determine alternans, wavebreak, and the patterns of reentrant arrhythmias. We developed a novel approach to characterize restitution using transfer functions. Transfer functions relate an input and an output quantity in terms of gain and phase shift in the complex frequency domain. We derived an analytical expression for the transfer function of interbeat intervals (IBIs) during conduction from one site (input) to another site downstream (output). Transfer functions can be efficiently obtained using a stochastic pacing protocol. Using simulations of conduction and extracellular mapping of strands of neonatal rat ventricular myocytes, we show that transfer functions permit the quantification of APD and CV restitution slopes when it is difficult to measure APD directly. We find that the normally positive CV restitution slope attenuates IBI variations. In contrast, a negative CV restitution slope (induced by decreasing extracellular [K⁺]) amplifies IBI variations with a maximum at the frequency of alternans. Hence, it potentiates alternans and renders conduction unstable, even in the absence of APD restitution. Thus, stochastic pacing and transfer function analysis represent a powerful strategy to evaluate restitution and the stability of conduction.     

Cardiac dynamics: a simplified model for action potential propagation

This paper analyzes a new semiphysiological ionic model, used recently to study reexitations and reentry in cardiac tissue [I.R. Cantalapiedra et al, PRE 82 011907 (2010)]. The aim of the model is to reproduce action potencial morphologies and restitution curves obtained, either from experimental data, or from more complex electrophysiological models. The model divides all ion currents into four groups according to their function, thus resulting into fast-slow and inward-outward currents. We show that this simplified model is flexible enough as to accurately capture the electrical properties of cardiac myocytes, having the advantage of being less computational demanding than detailed electrophysiological models. Under some conditions, it has been shown to be amenable to mathematical analysis. The model reproduces the action potential (AP) change with stimulation rate observed both experimentally and in realistic models of healthy human and guinea pig myocytes (TNNP and LRd models, respectively). When simulated in a cable it also gives the right dependence of the conduction velocity (CV) with stimulation rate. Besides reproducing correctly these restitution properties, it also gives a good fit for the morphology of the AP, including the notch typical of phase 1. Finally, we perform simulations in a realistic geometric model of the rabbit’s ventricles, finding a good qualitative agreement in AP propagation and the ECG. Thus, this simplified model represents an alternative to more complex models when studying instabilities in wave propagation.     

Dynamical effects of diffusive cell coupling on cardiac excitation and propagation: a simulation study

Cell coupling is considered to be important for cardiac action potential propagation and arrhythmogenesis. We carried out computer simulations to investigate the effects of stimulation strength and cell-to-cell coupling on action potential duration (APD) restitution, APD alternans, and stability of reentry in models of isolated cell, one-dimensional cable, and two-dimensional tissue. Phase I formulation of the Luo and Rudy action potential model was used. We found that stronger stimulation resulted in a shallower APD restitution curve and onset of APD alternans at a faster pacing rate. Reducing diffusive coupling between cells prolonged APD. Weaker diffusive currents along the direction of propagation steepened APD restitution and caused APD alternans to occur at a slower pacing rate in tissue. Diffusive current due to curvature changed APD but had little effect on APD restitution slope and onset of instability. Heterogeneous cell coupling caused APD inhomogeneities in space. Reduction in coupling strength either uniformly or randomly had little effect on the rotation period and stability of a reentry, but random cell decoupling slowed the rotation period and, thus, stabilized the reentry, preventing it from breaking up into multiple waves. Therefore, in addition to its effects on action potential conduction velocity, diffusive cell coupling also affects APD in a rate-dependent manner, causes electrophysiological heterogeneities, and thus modulates the dynamics of cardiac excitation. These effects are brought about by the modulation of ionic current activation and inactivation.     

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 700 ms. 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.     

Effects of Na(+) channel and cell coupling abnormalities on vulnerability to reentry: a simulation study

The role of dynamic instabilities in the initiation of reentry in diseased (remodeled) hearts remains poorly explored. Using computer simulations, we studied the effects of altered Na(+) channel and cell coupling properties on the vulnerable window (VW) for reentry in simulated two-dimensional cardiac tissue with and without dynamic instabilities. We related the VW for reentry to effects on conduction velocity, action potential duration (APD), effective refractory period dispersion and restitution, and concordant and discordant APD alternans. We found the following: 1). reduced Na(+) current density and slowed recovery promoted postrepolarization refractoriness and enhanced concordant and discordant APD alternans, which increased the VW for reentry; 2). uniformly reduced cell coupling had little effect on cellular electrophysiological properties and the VW for reentry. However, randomly reduced cell coupling combined with decoupling promoted APD dispersion and alternans, which also increased the VW for reentry; 3). the combination of decreased Na(+) channel conductance, slowed Na(+) channel recovery, and cellular uncoupling synergistically increased the VW for reentry; and 4) the VW for reentry was greater when APD restitution slope was steep than when it was flat. In summary, altered Na(+) channel and cellular coupling properties increase vulnerability to reentrant arrhythmias. In remodeled hearts with altered Na(+) channel properties and cellular uncoupling, dynamic instabilities arising from electrical restitution exert important influences on the VW for reentry.     

Structural Heterogeneity Modulates Effective Refractory Period: A Mechanism of Focal Arrhythmia Initiation

Reductions in electrotonic loading around regions of structural and electrophysiological heterogeneity may facilitate capture of focal triggered activity, initiating reentrant arrhythmias. How electrotonic loading, refractoriness and capture of focal ectopics depend upon the intricate nature of physiological structural anatomy, as well as pathological tissue remodelling, however, is not well understood. In this study, we performed computational bidomain simulations with anatomically-detailed models representing the rabbit left ventricle. We used these models to quantify the relationship between local structural anatomy and spatial heterogeneity in action potential (AP) characteristics, electrotonic currents and effective refractory periods (ERPs) under pacing and restitution protocols. Regions surrounding vessel cavities, in addition to tissue surfaces, had significantly lower peak downstream electrotonic currents than well coupled myocardium ( vs A/cm²), with faster maximum AP upstroke velocities ( vs mV/ms), although noticeably very similar APDs ( vs ms) and AP restitution properties. Despite similarities in APDs, ERPs in regions of low electrotonic load in the vicinity of surfaces, intramural vessel cavities and endocardial structures were up to ms shorter compared to neighbouring well-coupled tissue, leading to regions of sharp ERP gradients. Consequently, focal extra-stimuli timed within this window of ERP heterogeneity between neighbouring regions readily induced uni-directional block, inducing reentry. Most effective induction sites were within channels of low ERPs between large vessels and epicardium. Significant differences in ERP driven by reductions in electrotonic loading due to fine-scale physiological structural heterogeneity provides an important mechanism of capture of focal activity and reentry induction. Application to pathological ventricles, particularly myocardial infarction, will have important implications in anti-arrhythmia therapy.     

Steeper restitution slopes across right ventricular endocardium in patients with cardiomyopathy at high risk of ventricular arrhythmias

Steep action potential duration (APD) restitution slopes (>1) and spatial APD restitution heterogeneity provide the substrate for ventricular fibrillation in computational models and experimental studies. Their relationship to ventricular arrhythmia vulnerability in human cardiomyopathy has not been defined. Patients with cardiomyopathy [left ventricular (LV) ejection fraction <40%] and no history of ventricular arrhythmias underwent risk stratification with programmed electrical stimulation or T wave alternans (TWA). Low-risk patients (n = 10) had no inducible ventricular tachycardia (VT) or negative TWA, while high-risk patients (n = 8) had inducible VT or positive TWA. Activation recovery interval (ARI) restitution slopes were measured simultaneously from 10 right ventricular (RV) endocardial sites during an S1-S2 pacing protocol. ARI restitution slope heterogeneity was defined as the coefficient of variation of slopes. Mean ARI restitution slope was significantly steeper in the high-risk group compared with the low-risk group [1.16 (SD 0.31) vs. 0.59 (SD 0.19), P = 0.0002]. The proportion of endocardial recording sites with a slope >1 was significantly larger in the high-risk patients [47% (SD 35) vs. 13% (SD 21), P = 0.022]. Spatial heterogeneity of ARI restitution slopes was similar between the two groups [29% (SD 16) vs. 39% (SD 34), P = 0.48]. There was an inverse linear relationship between the ARI restitution slope and the minimum diastolic interval (P < 0.001). In cardiomyopathic patients at high risk of ventricular arrhythmias, ARI restitution slopes along the RV endocardium are steeper, but restitution slope heterogeneity is similar compared with those at low risk. Steeper ARI restitution slopes may increase the propensity for ventricular arrhythmias in patients with impaired left ventricular function.     

Induction of ventricular arrhythmias following mechanical impact: A simulation study in 3D

Commotio cordis, mechanical induction of heart rhythm disturbances, including sudden cardiac death, in the absence of corresponding structural damage, has been reported with increasing frequency in young individuals participating in sporting activities. Recently, the electrophysiological changes during c. cordis have been attributed to mechano-electric feedback, and particularly, to the recruitment of stretch-activated ion channels. The underlying mechanisms, however, by which a mechanical impact results in ventricular fibrillation, remain unknown. This study employs a 3D realistic model of rabbit ventricular geometry and fiber orientation to elucidate the electrophysiological mechanisms involved in arrhythmia induction following acute mechanical stimulation of the heart. Impact effects are modeled through stretch-activated ion channel activation in a 3D region of the ventricles representing the impact profile. Both cation-nonselective and potassium-selective stretch-activated ion channels are recruited upon mechanical impact. The impact is administered at various coupling intervals following pacing at the apex. To aid in the interpretation of results, the effect of mechanical stimulation on single cell action potentials is also examined. The results demonstrate that the region of impact is characterized by different types of cellular responses, including generation of a new action potential, shortening, or lengthening of action potential duration. The impact induces sustained reentry only when (1) a new activation is elicited by mechanical stimulation (caused by activation of cation-nonselective stretch-activated ion channels), and (2) upon return to the original region of impact, this activation does not encounter an extension of action potential duration (prevented by activation of potassium-selective stretch-activated ion channels).     

Toward prediction of the local onset of alternans in the heart

A beat-to-beat variation in the cardiac action potential duration is a phenomenon known as alternans. Alternans has been linked to ventricular fibrillation, and thus the ability to predict the onset of alternans could be clinically beneficial. Theoretically, it has been proposed that the slope of a restitution curve, which relates the duration of the action potential to the preceding diastolic interval, can predict the onset of alternans. Experimentally, however, this hypothesis has not been consistently proven, mainly because of the intrinsic complexity of the dynamics of cardiac tissue. It was recently shown that the restitution portrait, which combines several restitution curves simultaneously, is associated with the onset of alternans in isolated myocytes. Our main purpose in this study was to determine whether the restitution portrait is correlated with the onset of alternans in the heart, where the dynamics include a spatial complexity. We performed optical mapping experiments in isolated Langendorff-perfused rabbit hearts in which alternans was induced by periodic pacing at different frequencies, and identified the local onset of alternans, B^onset. We identified two regions of the heart: the area that exhibited alternans at B^onset (1:1_alt) and the area that did not (1:1). We constructed two-dimensional restitution portraits for the epicardial surface of the heart and measured the spatial distribution of three different slopes (the dynamic restitution slope, , and two local S1-S2 slopes, S₁₂ and ) separately for these two regions. We found that the S₁₂ and slopes differed significantly between the 1:1_alt and 1:1 regions just before the onset of alternans, and slopes were statistically similar. In addition, we found that the slopes of the dynamic restitution curve S_dyn were also statistically similar between these two regions. On the other hand, the quantitative values of all slopes were significantly different from the theoretically predicted value of one. These results demonstrate that the slopes measured in the restitution portrait correlate with the onset of alternans in the heart.     

Effects of Na(+) and K(+) channel blockade on vulnerability to and termination of fibrillation in simulated normal cardiac tissue

Na(+) and K(+) channel-blocking drugs have anti- and proarrhythmic effects. Their effects during fibrillation, however, remain poorly understood. We used computer simulation of a two-dimensional (2-D) structurally normal tissue model with phase I of the Luo-Rudy action potential model to study the effects of Na(+) and K(+) channel blockade on vulnerability to and termination of reentry in simulated multiple-wavelet and mother rotor fibrillation. Our main findings are as follows: 1) Na(+) channel blockade decreased, whereas K(+) channel blockade increased, the vulnerable window of reentry in heterogeneous 2-D tissue because of opposing effects on dynamical wave instability. 2) Na(+) channel blockade increased the cycle length of reentry more than it increased refractoriness. In multiple-wavelet fibrillation, Na(+) channel blockade first increased and then decreased the average duration or transient time () of fibrillation. In mother rotor fibrillation, Na(+) channel blockade caused peripheral fibrillatory conduction block to resolve and the mother rotor to drift, leading to self-termination or sustained tachycardia. 3) K(+) channel blockade increased dynamical instability by steepening action potential duration restitution. In multiple-wavelet fibrillation, this effect shortened because of enhanced wave instability. In mother rotor fibrillation, this effect converted mother rotor fibrillation to multiple-wavelet fibrillation, which then could self-terminate. Our findings help illuminate, from a theoretical perspective, the possible underlying mechanisms of termination of different types of fibrillation by antiarrhythmic drugs.     

Spatial heterogeneity of the restitution portrait in rabbit epicardium

Spatial heterogeneity of repolarization can provide a substrate for reentry to occur in myocardium. This heterogeneity may result from spatial differences in action potential duration (APD) restitution. The restitution portrait (RP) measures many aspects of rate-dependent restitution: the dynamic restitution curve (RC), S1-S2 RC, and short-term memory response. We used the RP to characterize epicardial patterns of spatial heterogeneity of restitution that were repeatable across animals. New Zealand White rabbit ventricles were paced from the epicardial apex, midventricle, or base, and optical action potentials were recorded from the same three regions. A perturbed downsweep pacing protocol was applied that measured the RP over a range of cycle lengths from 1,000 to 140 ms. The time constant of short-term memory measured close to the stimulus was dependent on location. In the midventricle the mean time constant was 19.1 +/- 1.1 s, but it was 39% longer at the apex (P < 0.01) and 23% longer at the base (P = 0.03). The S1-S2 RC slope was dependent on pacing site (P = 0.015), with steeper slope when pacing from the apex than from the base. There were no significant repeatable spatial patterns in steady-state APD at all cycle lengths or in dynamic RC slope. These results indicate that transient patterns of epicardial heterogeneity of APD may occur after a change in pacing rate. Thus it may affect cardiac electrical stability at the onset of a tachycardia or during a series of ectopic beats. Differences in restitution with respect to pacing site suggest that vulnerability may be affected by the location of reentry or ectopic foci.     

Effects of simulated ischemia on spiral wave stability

Regional hyperkalemia during acute myocardial ischemia is a major factor promoting electrophysiological abnormalities leading to ventricular fibrillation (VF). However, steep action potential duration restitution, recently proposed to be a major determinant of VF, is typically decreased rather than increased by hyperkalemia and acute ischemia. To investigate this apparent contradiction, we simulated the effects of regional hyperkalemia and other ischemic components (anoxia and acidosis) on the stability of spiral wave reentry in simulated two-dimensional cardiac tissue by use of the Luo-Rudy ventricular action potential model. We found that the hyperkalemic "ischemic" area promotes wavebreak in the surrounding normal tissue by accelerating the rate of spiral wave reentry, even after the depolarized ischemic area itself has become unexcitable. Furthermore, wavebreak and fibrillation can be prevented if the dynamical instability of the normal tissue is reduced significantly by targeting electrical restitution properties, suggesting a novel therapeutic approach.     

High-frequency oscillations and circadian rhythm of the membrane potential in Spinach leaves

The microelectrode technique was used to follow oscillations in membrane potential in mesophyll cells of spinach (Spinacia oleracea L.) during exposure do different photoperiodic conditions. Both high-frequency oscillations and circadian variations were observed. The circadian rhythm was imposed on the period of high-frequency oscillation during short days as well as in continuous light: The free-running period was 25.2 h. The average period of high-frequency oscillation increased from 7.64 min in the dark to 19.95 min in the light within several minutes after dark to light transition. This period length coincides with the established period length for oscillations in the redox potential in the chloroplast suspensions of spinach.Abbreviations CL continuous light - SD short day - MP membrane potential