Estimation of the left ventricular relaxation time constant tau requires consideration of the pressure asymptote

The left ventricular isovolumic pressure decay, obtained by cardiac catheterization, is widely characterized by the time constant tau of the exponential regression p(t)=Pomega+(P0-Pomega)exp(-t/tau). However, several authors prefer to prefix Pomega=0 instead of coestimating the pressure asymptote empirically; others present tau values estimated by both methods that often lead to discordant results and interpretation of lusitropic changes. The present study aims to clarify the relations between the tau estimates from both methods and to decide for the more reliable estimate. The effect of presetting a zero asymptote on the tau estimate was investigated mathematically and empirically, based on left ventricular pressure decay data from isolated ejecting rat and guinea pig hearts at different preload and during spontaneous decrease of cardiac function. Estimating tau with preset Pomega=0 always yields smaller values than the regression with empirically estimated asymptote if the latter is negative and vice versa. The sequences of tau estimates from both methods can therefore proceed in reverse direction if tau and Pomega change in opposite directions between the measurements. This is exemplified by data obtained during an increasing preload in spontaneously depressed isolated hearts. The estimation of the time constant of isovolumic pressure fall with a preset zero asymptote is heavily biased and cannot be used for comparing the lusitropic state of the heart in hemodynamic conditions with considerably altered pressure asymptotes.     

Physical determinants of left ventricular isovolumic pressure decline: model prediction with in vivo validation

The rapid decline in pressure during isovolumic relaxation (IVR) is traditionally fit algebraically via two empiric indexes: tau, the time constant of IVR, or tau(L), a logistic time constant. Although these indexes are used for in vivo diastolic function characterization of the same physiological process, their characterization of IVR in the pressure phase plane is strikingly different, and no smooth and continuous transformation between them exists. To avoid the parametric discontinuity between tau and tau(L) and more fully characterize isovolumic relaxation in mechanistic terms, we modeled ventricular IVR kinematically, employing a traditional, lumped relaxation (resistive) and a novel elastic parameter. The model predicts IVR pressure as a function of time as the solution of d(2)P/dt(2) + (1/micro)dP/dt + E(k)P = 0, where micro (ms) is a relaxation rate (resistance) similar to tau or tau(L) and E(k) (1/s(2)) is an elastic (stiffness) parameter (per unit mass). Validation involved analysis of 310 beats (10 consecutive beats for 31 subjects). This model fit the IVR data as well as or better than tau or tau(L) in all cases (average root mean squared error for dP/dt vs. t: 29 mmHg/s for model and 35 and 65 mmHg/s for tau and tau(L), respectively). The solution naturally encompasses tau and tau(L) as parametric limits, and good correlation between tau and 1/microE(k) (tau = 1.15/microE(k) - 11.85; r(2) = 0.96) indicates that isovolumic pressure decline is determined jointly by elastic (E(k)) and resistive (1/mu) parameters. We conclude that pressure decline during IVR is incompletely characterized by resistance (i.e., tau and tau(L)) alone but is determined jointly by elastic (E(k)) and resistive (1/micro) mechanisms.     

Dissociation between myocardial relaxation and diastolic stiffness in the stunned heart: its prevention by ischemic preconditioning

The effects of myocardial stunning and ischemic preconditioning on left-ventricular developed pressure and end-diastolic pressure (diastolic stiffness) as well as on coronary-perfusion pressure were examined in isolated isovolumic rabbit hearts. The isovolumic relaxation was evaluated, and the time constant of pressure decay during the isovolumic period was calculated. Our experimental protocol comprised: 1) myocardial stunning-global ischemia (15 min) followed by reperfusion (30 min); 2) myocardial stunning-global ischemia (20 min) followed by reperfusion (30 min); and 3) ischemic preconditioning — a single cycle of brief global ischemia and reperfusion (5 min each), before a second ischemic period, of 20-min duration. There was no effect upon systolic and diastolic parameters when 15 and 20 minutes of ischemia were evaluated. In both stunned groups the left ventricular developed pressure first recovered to near control values, but then stabilized at only 60% of the control values. Whereas the isovolumic relaxation time constant was increased after 5 min of reperfusion, and return to control values at late reperfusion, the end diastolic pressure remained elevated during the entire period. Values of dP/dV calculated at common pressure levels, were used as a second index of diastolic stiffness. They were increased after stunning, as also was the coronary perfusion pressure. When the heart was preconditioned with a single episode of ischemia, the systolic and diastolic alterations were completely abolished. We thus concluded that diastolic abnormalities incurred by myocardial stunning consist in both an increase in diastolic stiffness and an early impairment of isovolumic relaxation. The increase in stiffness cannot result from incomplete relaxation since these two parameters become temporally dissociated during the reperfusion period.     

Diastolic ventricular-vascular stiffness and relaxation relation: elucidation of coupling via pressure phase plane-derived indexes

Because systole and diastole are coupled and systolic ventricular-vascular coupling has been characterized, we hypothesize that diastolic ventricular-vascular coupling (DVVC) exists and can be characterized in terms of relaxation and stiffness. To characterize and elucidate DVVC mechanisms, we introduce time derivative of pressure (dP/dt) vs. time-varying pressure [P(t)] (pressure phase plane, PPP)-derived analogs of ventricular and vascular "stiffness" and relaxation parameters. Although volume change (dV) = 0 during isovolumic periods, and time-varying left ventricular (LV) stiffness, typically expressed as change in pressure per unit change in volume (dP/dV), is undefined, our formulation allows determination of a PPP-derived stiffness analog during isovolumic contraction and relaxation. Similarly, an aortic stiffness analog is also derivable from the PPP. LV relaxation was characterized via tau, the time constant of isovolumic relaxation, and vascular (aortic pressure decay) relaxation was characterized in terms of its equivalent (windkessel) exponential decay time constant kappa. The results show that PPP-derived systolic and diastolic ventricular and vascular stiffness are strongly coupled [K(Ao)(+)=1.71(K(LV)(+)) +154, r=0.86; K(Ao)(-)=0.677(K(LV)(-))-5.53, r=0.86]. In support of the DVVC hypothesis, a strong linear correlation between relaxation (rate of pressure decay) indexes kappa and tau (kappa = 9.89tau - 90.3, r = 0.81) was also observed. The correlations observed underscore the role of long-term, steady-state DVVC as a diastolic function determinant. Awareness of the PPP-derived DVVC parameters provides insight into mechanisms and facilitates quantification of arterial stiffening and associated increase in diastolic chamber stiffness. The PPP method provides a tool for quantitative assessment and determination of the functional coupling of the vasculature to diastolic function.     

Spatio-temporal attributes of left ventricular pressure decay rate during isovolumic relaxation

Global left ventricular (LV) isovolumic relaxation rate has been characterized: 1) via the time constant of isovolumic relaxation τ or 2) via the logistic time constant τ(L). An alternate kinematic method, characterizes isovolumic relaxation (IVR) in accordance with Newton's Second Law. The model's parameters, stiffness E(k), and damping/relaxation μ result from best fit of model-predicted pressure to in vivo data. All three models (exponential, logistic, and kinematic) characterize global relaxation in terms of pressure decay rates. However, IVR is inhomogeneous and anisotropic. Apical and basal LV wall segments untwist at different times and rates, and transmural strain and strain rates differ due to the helically variable pitch of myocytes and sheets. Accordingly, we hypothesized that the exponential model (τ) or kinematic model (μ and E(k)) parameters will elucidate the spatiotemporal variation of IVR rate. Left ventricular pressures in 20 subjects were recorded using a high-fidelity, multipressure transducer (3 cm apart) catheter. Simultaneous, dual-channel pressure data was plotted in the pressure phase-plane (dP/dt vs. P) and τ, μ, and E(k) were computed in 1631 beats (average: 82 beats per subject). Tau differed significantly between the two channels (P < 0.05) in 16 of 20 subjects, whereas μ and E(k) differed significantly (P < 0.05) in all 20 subjects. These results show that quantifying the relaxation rate from data recorded at a single location has limitations. Moreover, kinematic model based analysis allows characterization of restoring (recoil) forces and resistive (crossbridge uncoupling) forces during IVR and their spatio-temporal dependence, thereby elucidating the relative roles of stiffness vs. relaxation as IVR rate determinants.     

Indexes of diastolic RV function: load dependence and changes after chronic RV pressure overload in lambs

Diastolic function is a major determinant of ventricular performance, especially when loading conditions are altered. We evaluated biventricular diastolic function in lambs and studied possible load dependence of diastolic parameters [minimum first derivative of pressure vs. time (dP/dt(min)) and time constant of isovolumic relaxation (tau)] in normal (n = 5) and chronic right ventricular (RV) pressure-overloaded (n = 5) hearts by using an adjustable band on the pulmonary artery (PAB). Pressure-volume relations were measured during preload reduction to obtain the end-diastolic pressure-volume relationship (EDPVR). In normal lambs, absolute dP/dt(min) and tau were lower in the RV than in the left ventricle whereas the chamber stiffness constant (b) was roughly the same. After PAB, RV tau and dP/dt(min) were significantly higher compared with control. The RV EDPVR indicated impaired diastolic function. During acute pressure reduction, both dP/dt(min) and tau showed a relationship with end-systolic pressure. These relationships could explain the increased dP/dt(min) but not the increased tau-value after banding. Therefore, the increased tau after banding reflects intrinsic myocardial changes. We conclude that after chronic RV pressure overload, RV early relaxation is prolonged and diastolic stiffness is increased, both indicative of impaired diastolic function.     

Stiffness and relaxation components of the exponential and logistic time constants may be used to derive a load-independent index of isovolumic pressure decay

In current practice, empirical parameters such as the monoexponential time constant tau or the logistic model time constant tauL are used to quantitate isovolumic relaxation. Previous work indicates that tau and tauL are load dependent. A load-independent index of isovolumic pressure decline (LIIIVPD) does not exist. In this study, we derive and validate a LIIIVPD. Recently, we have derived and validated a kinematic model of isovolumic pressure decay (IVPD), where IVPD is accurately predicted by the solution to an equation of motion parameterized by stiffness (Ek), relaxation (tauc), and pressure asymptote (Pinfinity) parameters. In this study, we use this kinematic model to predict, derive, and validate the load-independent index MLIIIVPD. We predict that the plot of lumped recoil effects [Ek.(P*max-Pinfinity)] versus resistance effects [tauc.(dP/dtmin)], defined by a set of load-varying IVPD contours, where P*max is maximum pressure and dP/dtmin is the minimum first derivative of pressure, yields a linear relation with a constant (i.e., load independent) slope MLIIIVPD. To validate the load independence, we analyzed an average of 107 IVPD contours in 25 subjects (2,669 beats total) undergoing diagnostic catheterization. For the group as a whole, we found the Ek.(P*max-Pinfinity) versus tauc.(dP/dtmin) relation to be highly linear, with the average slope MLIIIVPD=1.107+/-0.044 and the average r2=0.993+/-0.006. For all subjects, MLIIIVPD was found to be linearly correlated to the subject averaged tau (r2=0.65), tauL(r2=0.50), and dP/dtmin (r2=0.63), as well as to ejection fraction (r2=0.52). We conclude that MLIIIVPD is a LIIIVPD because it is load independent and correlates with conventional IVPD parameters. Further validation of MLIIIVPD in selected pathophysiological settings is warranted.     

Differences in the reaction of rat and guinea pig hearts to ischemia and reperfusion

The response of rat and guinea-pig hearts to ischemia and reperfusion has been studied in identical conditions. Total 15-min ischemia of isolated rat hearts at 36 degrees C induced an almost 3-fold rise in isovolumic left ventricular diastolic pressure as well as a fall in the developed pressure and heart rate. Guinea-pig hearts, in the same conditions, exhibited a more steep fall in heart rate, with no rise in diastolic pressure. With constant heart rate produced by electrical stimulation at 4 Hz, the difference between two groups remained unchanged, while a more rapid fall in developed pressure in guinea-pig hearts coincided with a more profound fall in extracellular pH and almost a 2-fold rise in extracellular K+ activity. Rapid elimination of K+ and H+ at the early stages of reperfusion was followed by fibrillation in the majority of guinea-pig hearts, while no fibrillation was observed in rat hearts.     

Differential right and left ventricular diastolic tolerance to acute afterload and NCX gene expression in Wistar rats

This study evaluated right ventricular (RV) and left ventricular (LV) diastolic tolerance to afterload and SERCA2a, phospholamban and sodium-calcium exchanger (NCX) gene expression in Wistar rats. Time constant tau and end diastolic pressure-dimension relation (EDPDR) were analyzed in response to progressive RV or LV afterload elevations, induced by beat-to-beat pulmonary trunk or aortic root constrictions, respectively. Afterload elevations decreased LV- tau, but increased RV-tau. Whereas LV- tau analyzed the major course of pressure fall, RV- tau only assessed the last fourth. Furthermore, RV afterload elevations progressively upward shifted RV EDPDR, whilst LV afterload elevations did not change LV-EDPDR. SERCA2a and phospholamban mRNA were similar in both ventricles. NCX-mRNA was almost 50 % lower in RV than in LV. Left ventricular afterload elevations, therefore, accelerated the pressure fall and did not induce diastolic dysfunction, indicating high LV diastolic tolerance to afterload. On the contrary, RV afterload elevations decelerated the late RV pressure fall and induced diastolic dysfunction, indicating small RV diastolic tolerance to afterload. These results support previous findings relating NCX with late Ca(2+) reuptake, late relaxation and diastolic dysfunction.     

MRI assessment of LV relaxation by untwisting rate: a new isovolumic phase measure of tau

Most noninvasive measures of diastolic function are made during left ventricular (LV) filling and are therefore subject to "pseudonormalization," because variation in left atrial (LA) pressure may confound the estimation of relaxation rate. Counterclockwise twist of the LV develops during ejection, but untwisting occurs rapidly during isovolumic relaxation, before mitral opening. We hypothesized that the rate of untwisting might reflect the process of relaxation independent of LA pressure. Recoil rate (RR), the velocity of LV untwisting, was measured by tagged magnetic resonance imaging and regressed against the relaxation time constant (tau), recorded by catheterization, in 10 dogs at baseline and after dobutamine, saline, esmolol, and methoxamine treatment. RR correlated closely (average r = -0.86) with tau and was unaffected by elevated LA pressure. Multiple regression showed that tau, but not LA or aortic pressure, was an independent predictor of RR (P < 0.0001, P = 0.99, and P = 0.18, respectively). The rate of recoil of torsion, determined wholly noninvasively, provides an isovolumic phase, preload-independent assessment of LV relaxation. Use of this novel parameter should allow the detailed study of diastolic function in states known to affect filling rates, such as aging, hypertension, and congestive heart failure.     

Half-logistic time constant: a more reliable lusitropic index than monoexponential time constant regardless of temperature in canine left ventricle

Temperature changes influence cardiac diastolic function. The monoexponential time constant (tauE), which is a conventional lusitropic index of the rate of left ventricular (LV) pressure fall, increases with cooling and decreases with warming. We have proposed that a half-logistic time constant (tauL) is a better lusitropic index than tauE at normothermia. In the present study, we investigated whether tauL can remain a superior measure as temperature varies. The isovolumic relaxation LV pressure curves from the minimum of the first time derivative of LV pressure (dP/dtmin) to the LV end-diastolic pressure were analyzed at 30, 33, 36, 38, and 40 degrees C in excised, cross-circulated canine hearts. tauL and tauE were evaluated by curve-fitting using the least squares method and applying the half-logistic equation, P(t) = PA/[1 + exp(t/tauL)] + PB, and the monoexponential equation, P(t) = P0exp(-t/tauE) + Pinfinity. Both tauL and tauE increased significantly with decreasing temperature and decreased with increasing temperature. The half-logistic correlation coefficient (r) values were significantly higher than the monoexponential r values at the 5 above-mentioned temperatures. This implies that the superiority of the goodness of the half-logistic fit is not temperature dependent. The half-logistic model characterizes the amplitude and time course of LV pressure fall more reliably than the monoexponential model. Hence, we concluded that tauL is a more useful lusitropic index regardless of temperature.     

Isovolumic pressure-to-early rapid filling decay rate relation: model-based derivation and validation via simultaneous catheterization echocardiography.

Transmitral Doppler echocardiography is the preferred method of noninvasive diastolic function assessment. Correlations between catheterization-based measures of isovolumic relaxation (IVR) and transmitral, early rapid filling (Doppler E-wave)-derived parameters have been observed, but no model-based, causal explanation has been offered. IVR has also been characterized in terms of its duration as IVR time (IVRT) and by tau, the time-constant of IVR, by approximating the terminal left ventricular IVR pressure contour as Pt= Pinfinity + P(o)e(-t/tau), where Pt is the continuity of pressure, Pinfinity and Po are constants, t is time, and tau is the time constant of IVR. To characterize the relation between IVR and early rapid filling more fully, simultaneous (micromanometric) left ventricular pressure and transmitral Doppler E-wave data from 25 subjects undergoing elective cardiac catheterization and having normal physiology were analyzed. The time constant tau was determined from the dP/dt vs. P (phase) plane and, simultaneous Doppler E-waves provided global indexes of chamber viscosity/relaxation (c), chamber stiffness (k), and load (xo). We hypothesize that temporal continuity of pressure decay at mitral valve opening and physiological constraints permit the algebraic derivation of linear relations relating 1/tau to both peak atrioventricular pressure gradient (kxo) and E-wave-derived viscosity/relaxation (c) but does not support a similar, causal (linear) relation between deceleration time and tau or IVRT. Both predicted linear relations were observed: kxo to 1/tau (r = 0.71) and viscosity/relaxation to 1/tau (r = 0.71). Similarly, as anticipated, only a weak linear correlation between deceleration time and IVRT or tau was observed (r = 0.41). The observed in vivo relationship provides insight into the isovolumic mechanism of relaxation and the changing-volume mechanism of early rapid filling via a link of the respective relaxation properties.     

Mechanisms underlying ischemic diastolic dysfunction: relation between rigor,calcium homeostasis,and relaxation rate

Increased diastolic chamber stiffness (upward arrow DCS) during ischemia may result from increased diastolic calcium, rigor, or reduced velocity of relaxation. We tested these potential mechanisms during severe ischemia in isolated red blood cell-perfused isovolumic rabbit hearts. Ischemia (coronary flow reduced 83%) reduced left ventricular (LV) contractility by 70%, which then remained stable. DCS progressively increased. When LV end-diastolic pressure had increased 5 mmHg, myofilament calcium responsiveness was altered with 50 mmol/l NH(4)Cl or 10 mmol/l butanedione monoxime. These affected contractility (i.e., a calcium-mediated force) but not upward arrow DCS. Second, quick length changes reversed upward arrow DCS, supporting a rigor mechanism. Third, ischemia increased the time constant of isovolumic pressure decline from 47 +/- 3 to 58 +/- 3 ms (P < 0.02) but concomitantly abbreviated the contraction-relaxation cycle, i.e., pressure dissipation occurred earlier without diastolic tetanization. Finally, to assess any link between rate of relaxation and upward arrow DCS, hearts were exposed to 10 mmol/l calcium. Calcium doubled contractility and accelerated relaxation velocity, but without affecting upward arrow DCS. Thus upward arrow DCS developed during ischemia despite severely reduced contractility via a rigor (and not calcium mediated) mechanism. Calcium resequestration capacity was preserved, and reduced relaxation velocity was not linked to upward arrow DCS.     

Load dependence of ventricular performance explained by model of calcium-myofilament interactions

Although a simple concept of load-independent behavior of the intact heart evolved from early studies of isolated, intact blood-perfused hearts, more recent studies showed that, as in isolated muscle, the mode of contraction (isovolumic vs. ejection) impacts on end-systolic elastance. The purpose of the present study was to test whether a four-state model of myofilament interactions with length-dependent rate constants could explain the complex contractile behavior of the intact, ejecting heart. Studies were performed in isolated, blood-perfused canine hearts with intracellular calcium transients measured by macroinjected aequorin. Measured calcium transients were used as the driving function for the model, and length-dependent rate constants yielding the highest concordance between measured and model-predicted midwall stress at different isovolumic volumes were determined. These length-dependent rate constants successfully predicted contractile behavior on ejecting contractions. This, along with additional model analysis, suggests that length-dependent changes in calcium binding affinity may not be an important factor contributing to load-dependent contractile performance in the intact heart under physiological conditions.     

Phase plane analysis of left ventricular hemodynamics.

We sought to extract additional physiological information from the time-dependent left ventricular (LV) pressure contour and thereby gain new insights into ventricular function. We used phase plane analysis to characterize high-fidelity pressure data in selected subjects undergoing elective cardiac catheterization. The standard hemodynamic indexes of LV systolic and diastolic function derived from the time-dependent LV pressure contour could be easily obtained using the phase plane method. Additional novel attributes of the phase plane pressure loop, such as phase plane pressure loop area, graphical representation of the isovolumic relaxation time constant, and quantitative measures of beat-to-beat systolic-diastolic coupling were characterized. The asymmetry between the pressures at which maximum isovolumic pressure rise and pressure fall occur, as well as their load dependence, were also easily quantitated. These results indicate that the phase plane method provides a novel window for physiological discovery and has theoretical and applied advantages in quantitative ventricular function characterization.     

Positive inotropic and negative lusitropic effects of endothelin receptor agonism in vivo

The endothelin (ET) system is involved in the regulation of myocardial function in health as well as in several diseases, such as congestive heart failure, myocardial infarction, and septic myocardial depression. Conflicting results have been reported regarding the acute contractile properties of ET-1. We therefore investigated the effects of intracoronary infusions of ET-1 and of the selective ET(B) receptor-selective agonist sarafotoxin 6c with increasing doses in anesthetized pigs. Myocardial effects were measured through analysis of the left ventricular pressure-volume relationship. ET-1 elicited increases in the myocardial contractile status (end-systolic elastance value of 0.94 +/- 0.11 to 1.48 +/- 0.23 and preload recruitable stroke work value of 68.7 +/- 4.7 to 83.4 +/- 7.2) that appear to be mediated through ET(A) receptors, whereas impairment in left ventricular isovolumic relaxation (tau = 41.5 +/- 1.4 to 58.1 +/- 5.0 and t(1/2) = 23.0 +/- 0.7 to 30.9 +/- 2.6, where tau is the time constant for pressure decay and t(1/2) is the half-time for pressure decay) was ET(B) receptor dependent. In addition, intravenous administration of ET-1 impaired ventricular relaxation but had no effect on contractility. Intracoronary sarafotoxin 6c administration caused impairments in left ventricular relaxation (tau from 43.3 +/- 1.8 to 54.4 +/- 3.4) as well as coronary vasoconstriction. In conclusion, ET-1 elicits positive inotropic and negative lusitropic myocardial effects in a pig model, possibly resulting from ET(A) and ET(B) receptor activation, respectively.     

Changes in rate of relaxation of sniffs with diaphragmatic fatigue in humans

The rate of relaxation of the diaphragm after stimulated (4 subjects) and voluntary (8 subjects) contractions was compared in normal young men. Stimulated contractions were induced by supramaximal unilateral phrenic nerve stimulation and voluntary contractions by short, sharp sniffs of varying tensions against an occluded airway. The rate of relaxation of the diaphragm was calculated from the rate of decline of transdiaphragmatic pressure (Pdi). In both conditions the maximum relaxation rate (MRR) was proportional to the peak transdiaphragmatic pressure (Pdi), whereas the time constant (tau) of the later exponential decline in Pdi was independent of Pdi. The mean +/- SE rate constant of relaxation (MRR/Pdi) was 0.0078 +/- 0.0002 ms-1 and the mean tau was 57 +/- 3.8 ms for stimulated contractions. The rate of relaxation after sniffs was not different, and it was not affected by either the lung volume at which occluded sniffs were performed (in the range of residual volume to functional residual capacity + 1 liter) or by the relative contribution gastric pressure made to Pdi. After diaphragmatic fatigue was induced by inspiring against a high alinear resistance there was a decrease in relaxation rate. In the 1st min postfatigue MRR/Pdi decreased (0.0063 +/- 0.0003 ms-1; P less than 0.005) and tau increased (83 +/- 5 ms; P less than 0.005). Both values returned to prefatigue levels within 5 min of the end of the studies. We conclude that the sniff may prove to be clinically useful in the detection of diaphragmatic fatigue.     

A paradigm for quantifying ventricular contraction

The left ventricle may be described as a time, volume and flow dependent pressure generator. First, isovolumic pressure is measured at various end-diastolic volumes. Next, pressure is adjusted to account for small changes accompanying ejection, denoted the ejection effect. The resulting analytical function can describe pressure generation and ventricular outflow of the ventricle under a wide range of contractile and vascular conditions. This paradigm is unique in separating isovolumic from ejecting ventricular properties, as well as ventricular from vascular conditions.     

Analysis of fluorescence anisotropy decays by a least square method

A method of fluorescence anisotropy decay analysis is described in this work. The transient anisotropy r(ex)(t) measured in a photocounting pulsefluorimeter is fitted by a non linear least square procedure to the ratio of convolutions of the apparatus response function g(t) by sums of appropriate exponential functions. This method takes rigorously into account the apparatus response function and is applicable to any shape of the later as well as to any values of fluorescence decay times and correlation times. The performances of the method have been tested with data simulated from measured response functions corresponding to an air lamp and a high pressure nitrogen lamp. The statistical standard errors of the anisotropy deca parameters have been found to be smaller than the standard errors previously calculated for the moment method. A systematic error delta in the fluorescence decay time entailed an error deltatheta in the correlation time such as Deltatheta/theta < deltatau/tau. By this method, good fitting of experimental data have been achieved very conveniently and accurately.