Single-exponential Fits Research Articles

Molecular Properties of Rhodopsin and Rod Function

Signal transduction in rod cells begins with photon absorption by rhodopsin and leads to the generation of an electrical response. The response profile is determined by the molecular properties of the phototransduction components. To examine how the molecular properties of rhodopsin correlate with the rod-response profile, we have generated a knock-in mouse with rhodopsin replaced by its E122Q mutant, which exhibits properties different from those of wild-type (WT) rhodopsin. Knock-in mouse rods with E122Q rhodopsin exhibited a photosensitivity about 70% of WT. Correspondingly, their single-photon response had an amplitude about 80% of WT, and a rate of decline from peak about 1.3 times of WT. The overall 30% lower photosensitivity of mutant rods can be explained by a lower pigment photosensitivity (0.9) and the smaller single-photon response (0.8). The slower decline of the response, however, did not correlate with the 10-fold shorter lifetime of the meta-II state of E122Q rhodopsin. This shorter lifetime became evident in the recovery phase of rod cells only when arrestin was absent. Simulation analysis of the photoresponse profile indicated that the slower decline and the smaller amplitude of the single-photon response can both be explained by the shift in the meta-I/meta-II equilibrium of E122Q rhodopsin toward meta-I. The difference in meta-III lifetime between WT and E122Q mutant became obvious in the recovery phase of the dark current after moderate photobleaching of rod cells. Thus, the present study clearly reveals how the molecular properties of rhodopsin affect the amplitude, shape, and kinetics of the rod response.

The Synthesis, Structural Characterization, and Receptor Specificity of the α-Conotoxin Vc1.1

The alpha-conotoxin Vc1.1 is a small disulfide-bonded peptide currently in development as a treatment for neuropathic pain. This study describes the synthesis, determination of the disulfide connectivity, and the determination of the three-dimensional structure of Vc1.1 using NMR spectroscopy. Vc1.1 was shown to inhibit nicotine-evoked membrane currents in isolated bovine chromaffin cells in a concentration-dependent manner and preferentially targets peripheral nicotinic acetylcholine receptor (nAChR) subtypes over central subtypes. Specifically, Vc1.1 is selective for alpha3-containing nAChR subtypes. The three-dimensional structure of Vc1.1 comprises a small alpha-helix spanning residues Pro6 to Asp11 and is braced by the I-III, II-IV disulfide connectivity seen in other alpha-conotoxins. A comparison of the structure of Vc1.1 with other alpha-conotoxins, taken together with nAChR selectivity data, suggests that the conserved proline at position 6 is important for binding, whereas a number of residues in the C-terminal portion of the peptide contribute toward the selectivity. The structure reported here should open new opportunities for further development of Vc1.1 or analogues as analgesic agents.

Kinetic Analysis of Steps in the Repair of Damaged DNA by Human O6-Alkylguanine-DNA Alkyltransferase

Rates of individual steps in the removal of alkyl groups from O6-methyl (Me) and -benzyl (Bz) guanine in oligonucleotides by human O6-alkylguanine DNA alkyltransferase (AGT) were estimated using rapid reaction kinetic methods. The overall reaction yields hyperbolic plots of rate versus AGT concentration for O6-MeG but linear plots for the O6-BzG reaction, which is approximately 100-fold faster. The binding of AGT and DNA (double-stranded 30-mer/36-mer complex) appears to be diffusion-limited. The rate of dissociation of the complex is approximately 25-fold slower (approximately 1 s(-1)) for DNA containing O6-MeG or O6-BzG than unmodified DNA. The fluorescent dC-analog 6-methylpyrrolo[2,3-d]pyrimidine-2(3H) one deoxyribonucleoside (pyrrolo dC), which pairs with G, was positioned opposite G, O6-MeG, or O6-BzG and used as a probe of the rate of base flipping. A rapid increase of fluorescence (k approximately 200 s(-1)) was observed with O6-MeG and O6-BzG and AGT but not with a Gly mutation at Arg128, which has been implicated in base flipping with crystal structures. Only weak and slower fluorescence changes were observed with G:pyrrolo dC or T:2-aminopurine pairs. These rate estimates were used in a kinetic model in which AGT binds and scans DNA rapidly, flips O6-alkylG residues, transfers the alkyl group in a chemical step that is rate-limiting in the case of O6-MeG but not O6-BzG, and releases the dealkylated DNA. The results explain the overall patterns of rates of alkyl group removal versus AGT concentration and the effects of the mutations, as well as the greater affinity of AGT for DNA with O6-alkylG lesions.

Coassembly of Two GluR6 Kainate Receptor Splice Variants within a Functional Protein Complex

Kainate receptors (KAR) are composed of several distinct subunits and splice variants, but the functional relevance of this diversity remains largely unclear. Here we show that two splice variants of the GluR6 subunit, GluR6a and GluR6b, which differ in their C-terminal domains, do not show distinct functional properties, but coassemble as heteromers in vitro and in vivo. Using a proteomic approach combining affinity purification and MALDI-TOF mass spectrometry, we found that GluR6a and GluR6b interact with two distinct subsets of cytosolic proteins mainly involved in Ca(2+) regulation of channel function and intracellular trafficking. Guided by these results, we provide evidence that the regulation of native KAR function by NMDA receptors depends on the heteromerization of GluR6a and GluR6b and interaction of calcineurin with GluR6b. Thus, GluR6a and GluR6b bring in close proximity two separate subsets of interacting proteins that contribute to the fine regulation of KAR trafficking and function.

Chronic morphine treatment alters N-methyl-D-aspartate receptors in freshly isolated neurons from nucleus accumbens.

Although there is now evidence of a role for N-methyl-D-aspartate (NMDA) receptors in nucleus accumbens (NAcc) neurons in the effects of chronic opiate treatment, the cellular and molecular mechanisms underlying this phenomenon are still unclear. Therefore, we studied the effects of chronic morphine on the pharmacological and biophysical properties of NMDA receptors in freshly isolated medium spiny neurons from NAcc. We found that chronic morphine treatment did not alter the affinity for NMDA receptor agonists such as glutamate, homoquinolinic acid, and NMDA, but decreased the affinity of glycine, the allosteric NMDA receptor coagonist, from 2.24 +/- 0.15 microM to 5.1 +/- 1.45 microM. Chronic morphine treatment also altered the affinity of two noncompetitive NMDA receptor antagonists, 7-chloro-kynurenic acid and ifenprodil. However, morphine had no effect on a third antagonist, D-(-)-2-amino-5-phosphonopentanoic acid. Single-exponential fits of desensitized NMDA current tails gave tau values ranging from 0.5 to 4 s in neurons from both control and morphine-treated rats. However, a shift to the left of the distribution of tau values after morphine treatment revealed that NMDA current desensitization rate was accelerated in a majority of NAcc neurons. Taken together with our recent molecular studies, our data are consistent with a shift away from NMDA receptor subunit (NR) NR2B and 2C function toward increased NR2A subunit expression or function after chronic morphine, a process that could alter excitability and integrative properties and may represent a neuroadaptation of NAcc medium spiny neurons underlying morphine dependence.

Presynaptic Capacitance Measurements and Ca2+ Uncaging Reveal Submillisecond Exocytosis Kinetics and Characterize the Ca2+ Sensitivity of Vesicle Pool Depletion at a Fast CNS Synapse

The intracellular Ca2+ sensitivity of synaptic vesicle fusion is an important determinant of transmitter release probability, but it is unknown for most CNS synapses. We combined whole-cell membrane capacitance measurements and Ca2+ uncaging at the large calyx of Held nerve terminals to determine the Ca2+ sensitivity of synaptic vesicle fusion at a glutamatergic CNS synapse, independent of recording EPSCs. Capacitance increases measured 30-50 msec after elevating the intracellular Ca2+ concentration ([Ca2+]i) by Ca2+ uncaging were half-maximal at approximately 5 microm [Ca2+]i. At 10 microm [Ca2+]i, capacitance increases reached maximal values (256 +/- 125 fF; mean +/- SD), indicating the depletion of an average pool of approximately 4000 readily releasable vesicles. Vesicle pool depletion was confirmed in cross-depletion experiments, in which capacitance responses were measured after Ca2+ uncaging, or after combined stimuli of prolonged presynaptic depolarizations and Ca2+ uncaging. To analyze the Ca2+-dependent rates of vesicle pool depletion, the capacitance rise after Ca2+ uncaging was fitted with single- or double-exponential functions. The fast time constants of double-exponential fits, and the time constants of single-exponential fits were 2-3 msec at 10-15 microm [Ca2+]i and reached submillisecond values at 30 microm [Ca2+]i. These results suggest that three to five readily releasable vesicles can fuse within <1 msec at each active zone of a calyx of Held, given that [Ca2+]i rises sufficiently high. Submillisecond kinetics of exocytosis are reached at significantly lower [Ca2+]i than at ribbon-type sensory synapses previously investigated by capacitance measurements.

In vivo analysis of end-plate noise of human extensor digitorum brevis muscle after intramuscularly injected botulinum toxin type A.

To evaluate changes occurring in the neuromuscular junction after injection with botulinum neurotoxin type A (BoTx), three healthy volunteers were injected with 10 U BoTx in the right extensor digitorum brevis muscle. In agreement with previous observations, amplitude of compound muscle action potential (CMAP) decreased to approximately 30% of the initial value at approximately day 8 and slowly returned to baseline values around day 250. Values of the acetylcholine receptor (AChR) open time were determined by spectral analysis of end-plate noise and from single exponential fits to the decay phase of individual miniature end-plate potentials (MEPPs). At baseline, the mean channel open times determined by end-plate noise analysis and the exponential fits were 1.1 +/- 0.2 ms and 1.20 +/- 0.04 ms, respectively. After BoTx injection, initially no end-plate noise could be recorded. From day 9 onwards, however, a gradual recurrence of end-plate noise was observed, with mean channel open times of approximately 2-5 ms, being maximal between days 20 to 140. In addition, the shape of many recorded MEPPs was different from the typical fast rising MEPPs observed at baseline. After day 80, end-plate noise gradually returned to normal and mean channel open times decreased slowly to baseline values. Our findings reflect the changed AChR characteristics of newly formed neuromuscular junctions, which are created after BoTx injection and gradually removed after restoration of the original neuromuscular junctions.

Nitrite binding to metmyoglobin and methemoglobin in comparison to nitric oxide binding

Nitrite binds reversibly to the ferriheme proteins metmyoglobin and methemoglobin in aqueous buffer solution at a physiological pH of 7.4. The spectral changes recorded for the formation of metMb(NO2-) differ significantly from those observed for the nitrosylation of metMb, which can be accounted for in terms of the different reaction products. Nitric oxide binding to metMb produces a nitrosyl product with Fe(II)-NO+ character, whereas the reaction with nitrite produces an Fe(III)-NO2- complex. The kinetics of the binding and release of nitrite by metMb and metHb were investigated by stopped-flow techniques at ambient and high pressure. The kinetic traces recorded for the reaction of nitrite with metMb exhibit excellent single-exponential fits, whereas nitrite binding to metHb is characterized by double-exponential kinetics which were assigned to the reactions of the alpha- and beta-chains of metHb with NO2-. The rate constants for the binding of nitrite to metMb and metHb were found to be much smaller than those reported for the binding of NO, such that nitrite impurities will not affect the latter reaction. The activation parameters (deltaH++,deltaS(ne),deltaV++) obtained from the temperature and pressure dependence of the reactions support the operation of a dissociative mechanism for the binding and release of nitrite, similar to that found for the binding and release of NO in metMb.

The Catalytic Tyr-9 of Glutathione S-Transferase A1-1 Controls the Dynamics of the C terminus

The glutathione S-transferase enzymes (GSTs) have a tyrosine or serine residue at their active site that hydrogen bonds to and stabilizes the thiolate anion of glutathione, GS(-). The importance of this hydrogen bond is obvious, in light of the enhanced nucleophilicity of GS(-) versus the protonated thiol. Several A-class GSTs contain a C-terminal segment that undergoes a ligand-dependent local folding reaction. Here, we demonstrate the effects of the Y9F substitution on binding affinity for glutathione conjugates and on rates of the order-disorder transition of the C terminus in rat GST A1-1. The equilibrium binding affinity of the glutathione conjugate, GS-NBD (NBD-Cl, 7-chloro-4-nitrobenzo-2-oxa-1, 3-diazole), was decreased from 4.09 microm to 0.641 microm upon substitution of Tyr-9 with Phe. This result was supported by isothermal titration calorimetry, with K(d) values of 1.51 microm and 0.391 microm for wild type and Y9F, respectively. The increase in binding affinity for the mutant is associated with dramatic decreases in rates for the C-terminal order-disorder transition, based on a stopped-flow kinetic analysis. The same effects were observed, qualitatively, for a second GSH conjugate, GS-ethacrynic acid. Apparently, the phenolic hydroxyl group of Tyr-9 is critical for orchestrating C-terminal dynamics and efficient product release, in addition to its role in lowering the pK(a) of GSH.

Rotational relaxation of perylene in n-alcohols and n-alkanes studied by two-photon-induced anisotropy decay

Two-photon-induced anisotropy decays have been measured for perylene in a series of n-alcohols and n-alkanes and in 2-propanol and cyclohexane. Perylene fluorescence was excited at 570 nm and detected by time-correlated single-photon counting. The two-photon-induced anisotropy decays r 1( t) and r 2( t) were measured for two-photon excitation with linearly and circularly polarized excitation pulses. Rotational correlation times determined from single-exponential fits to anisotropy decays were shorter for linearly polarized excitation than for circularly polarized excitation in all cases with the possible exception of 2-propanol at 20°C. This result demonstrates the existence of at least two underlying rotational diffusion times, which are weighted differently in the decays of the anisotropic distributions generated with linear and circular polarization as predicted by the theory of two-photon-induced anisotropy decay [J. Chem. Phys. 101 (1994) 10283]. The viscosity dependencies of the rotational correlation time follow slip (for alkanes) or sub-slip (for alcohols) boundary conditions. A solvent dependence in the ratio of rotational correlation times determined with linearly and circularly polarized excitation is shown to indicate non-hydrodynamic behavior and can be described as a solvent dependence in the relative values of the principal diffusion coefficients, i.e. the effective rotor shape. The results are consistent with more anisotropic diffusion in alcohols than in alkanes, with faster spinning about the axis normal to the plane of the molecule relative to tumbling about in-plane axes. The initial ( t=0) anisotropy values r 1(0) and r 2(0) together with the two-photon polarization parameter Ω were used to analyze the two-photon tensor governing the transition.

Measurement of mass diffusion coefficients using nonexponential forced Rayleigh scattering signals

Recent reports on mass-diffusion forced Rayleigh scattering (FRS) experiments have emphasized that the detected signal arises as a difference between two exponentially decaying fields diffracted from “complementary” ground-state and photoproduct population gratings. A mass-diffusion coefficient has nevertheless often been extracted by forcing a single-exponential fit to the data, especially in cases where the decay appears to be monotonic. In this paper, we use simulations and experiments to evaluate the accuracy of single-exponential fits for FRS profiles, and we propose a useful alternative method for obtaining a meaningful rate constant in cases in which the error in the single-exponential analysis is large. We begin by noting from the complementary grating model that (1) severe deviations from single-exponential decay can occur for an arbitrarily small (but nonzero) difference in the ground-state/photoproduct rate constants, and (2) the first cumulant of a FRS decay—in contrast to that of a dynamic light scattering profile—does not (in general) represent a physically useful decay rate. These statements apply to both monotonic and nonmonotonic decays. We then show that a combination of the first two FRS cumulants provides a physically useful mean rate constant. Finally, to address these issues experimentally, we have reexamined the diffusion of methyl red (MR) through 2-propanol at room temperature, a system previously analyzed using single-exponential fits. The new experiments, carried out at higher sensitivity than the previous studies, show that the MR/2-propanol signal is nonmonotonic. The geometric-mean diffusion coefficient obtained from the curvature of the local maximum is compared to the diffusion coefficient inferred from single-exponential fits, and it is found that the latter is larger by nearly a factor of 2. The results reported here should prove important in improving the accuracy of the FRS technique.

Time-Resolved Fluorescence Study of Conformational Dynamics in Opioid Peptides

Rotational correlation times have been determined from fluorescence anisotropy decays of the tyrosyl residue in the opioid pentapeptides DPDPE (Tyr-d-Pen-Gly-Phe-d-Pen), DPDPE(SH)2, and [Leu5]-enkephalin, revealing internal peptide motions. Fluorescence decays were measured by time-correlated single-photon counting. For all three peptides, the fluorescence emission is characterized by three-exponential intensity decays with amplitudes that are consistent with ground-state populations of rotamers of the tyrosyl side chain. Rotational correlation times in water determined from single-exponential fits are 80−130 ps, in good agreement with molecular dynamics simulations [Wang, Y.; Kuczera, K. J. Phys. Chem. 1996, 100, 2555−2563]. Internal peptide motions were studied by measurement of the rotational correlation times in solutions of 50% propylene glycol in water over the temperature range from 258 to 313 K. Two distinct temperature regions were observed. In the low-temperature regime the thermal viscosity coe...

Kinetics of lysozyme refolding: structural characterization of a non-specifically collapsed state using time-resolved X-ray scattering

We report time-resolved small angle X-ray scattering (SAXS) studies of the structural characteristics of the collapsed state of lysozyme from hen-egg white (HEL) obtained on initiating refolding by rapidly changing solvent conditions from 8 M to 1.1 M urea at pH 2.9. At this reduced pH the lifetime, of about one second, of the non-specifically collapsed ensemble is considerably prolonged relative to its value at pH 5.2. The SAXS studies are combined with time resolved measurements of tryptophan fluorescence and of the rate of formation of native molecules using interrupted refolding experiments. We observe large burst phase changes in intrinsic tryptophan fluorescence and in the radius of gyration ( R g) which is reduced from 22 Å in the fully unfolded state to ∼19 to 20 Å. Subsequent decrease of the R g to the value for native lysozyme (15 Å) follows the time course of formation of native molecules. Single exponential fits to the singular value decomposition (SVD) components of the SAXS data allow reconstruction of the SAXS profile at early time points of refolding. The results of this analysis suggest a globular shape of the collapsed state. A similar fit to the forward scattering amplitude, I(0), suggests that the collapsed state has a solvent accessible surface area which is considerably increased relative to that of the native protein. These results show directly that the non-specifically collapsed state formed during the burst phase in lysozyme refolding indeed represents a molecular compaction and a change in shape from a fully denatured random coil state (albeit restricted by disulfide bonds) to an ensemble of globular conformations which, however, have not yet formed a solvent-protected hydrophobic core.

Decay of calcium transients after electrical stimulation in rat fast- and slow-twitch skeletal muscle fibres.

1. Calcium transients were calculated from fura-2 fluorescence signals (corrected for kinetic delays in the Ca(2+)-fura-2 reaction) from single rat skeletal muscle fibres, either fully dissociated from the fast-twitch flexor digitorum brevis (FDB) muscle or in small bundles from the slow-twitch soleus muscle. Fibres or bundles were embedded in agarose gel to inhibit movement and stimulated by single or trains of 1-2 ms electrical pulses (100 Hz, 2-400 ms train duration). 2. The rate constant of decay of [Ca2+] determined from single-exponential fits to the final decay phase of [Ca2+] after a single action potential was considerably faster in FDB fibres than in soleus fibres. As the stimulation duration increased, the rate constant of [Ca2+] decay decreased for both the FDB and soleus fibres, but the effect was greater in FDB than in soleus fibres. 3. Using the magnitude of the decline in the rate constant of [Ca2+] decay with increasing stimulation duration as an index of relative contribution of the saturable Ca2+ binding sites on parvalbumin, subpopulations termed 'high', 'medium' and 'low', referring to estimated parvalbumin content, were determined within each group of FDB and soleus fibres. In fibres assigned to the 'high' and 'medium' groups, parvalbumin was the major contributor (50-73%) to the [Ca2+] decay rate constant after a single action potential. In fibres in the 'low' group, parvalbumin contributed only 0-28% to the rate constant of [Ca2+] decay. 4. Fluorescence recordings using mag-fura-2, a lower-affinity Ca2+ indicator expected to be in equilibrium with myoplasmic Ca2+, gave similar values for both the [Ca2+] decay rate constant after a single action potential and the decrease in this rate constant with increased stimulation duration, as found for the fura-2 [Ca2+] transients from FDB and soleus fibres. Thus, the observed differences in decay rate of Ca2+ were not introduced by kinetic correction of the fura-2 recordings, but are attributed to differences in the Ca2+ binding and transport properties of fast- and slow-twitch mammalian fibres.

Accommodation of mouse DRG growth cones to electrically induced collapse: Kinetic analysis of calcium transients and set‐point theory

Electrical stimulation causes growth cones of mouse dorsal root ganglion neurons to collapse. During chronic stimulation, however, growth cones resume motility. In addition, these growth cones are now resistant to the collapsing effects of subsequent stimulation, a process we term accommodation. We compared the kinetics of electrically induced Ca2+ transients in naive and accommodated growth cones in order to determine whether the accommodation process results from a change in the Ca2+ transient, or a change in the Ca2+ sensitivity of the growth cones. Three kinetics were determined: (1) the initial increase to peak Ca2+ levels produced by 10 Hz stimulation; (2) recovery from peak Ca2+ levels during stimulus trains lasting 15 min; and (3) clearing of Ca2+ from growth cones after terminating the stimulus. These kinetics were analyzed using single exponential fits to changes in fura-2 fluorescence ratios. The electrically evoked increase in Ca2+ was significantly slower in accommodated growth cones (tau = 6.0 s) compared to naive growth cones (tau = 1.4 s). Despite the slower increase of [Ca2+]i in accommodated growth cones, peak [Ca2+]i was similar to that reached in naive growth cones, and the steady-state Ca2+ level was significantly elevated after chronic stimulation. Thus, accommodated growth cones maintained outgrowth at [Ca2+]i that caused collapse initially. Time course experiments show that accommodation is a slow process (t 1/2 = about 3 h). Accommodation did not induce measurable changes in the rates of Ca2+ homeostasis during or after stimulus trains. The kinetics of Ca2+ recovery during (tau = 90 s) and after 15 min of stimulation (tau = 8.5 s) was not significantly different in accommodated versus naive growth cones. Rates of 45Ca2+ efflux were also similar in both types of growth cones. These results suggest two regulatory processes contributing to growth cone motility during chronic stimulation: (1) recovery of [Ca2+]i to levels permissive to neurite outgrowth, and (2) an increase in the range of optimal [Ca2+]i for growth cone motility. These adaptive responses of mammalian growth cones to chronic stimulation could be involved in the modulation of CNS development by electrical activity of neurons.

Multiple kinetic components of sodium channel inactivation in rabbit Schwann cells.

1. We have studied the kinetics of inactivation of sodium currents evoked in Schwann cells from neonatal and adult rabbits with patch-clamp recording techniques. The decay both of whole-cell currents and of ensemble currents obtained from outside-out patches was reasonably well-described by single exponential fits which gave values for the time constant, tau h, similar to those found for such currents in nerve (about 0.5 ms at 0 mV). Although inclusion of an additional exponential component usually improved the fits to the decay of these currents, the relative amplitude of the slower component (time constant 3-6 ms) was always small. 2. The time course of recovery from steady-state inactivation clearly consisted of two exponential components. At -120 mV, recovery from steady-state inactivation at -50 mV consisted of a fast component with a time constant of 2.2 +/- 0.2 ms and a much slower component with a time constant of 1.2 +/- 0.2 s (n = 9). The relative amplitude of the slow component (expressed as a fraction of the sodium current when inactivation was removed completely) was 0.56 +/- 0.03. The corresponding amplitudes of the slow component when similar experiments were done from holding potentials of 0 and -70 mV were 0.80 +/- 0.04 and 0.36 +/- 0.03, respectively (n = 5 and 8). 3. The onset of steady-state inactivation also followed a bi-exponential time course. The time constant of the slower component was similar at each potential examined (0, -50 and -70 mV), being 6-7 s. The relative amplitude of the slow component of onset depended on membrane potential, and it was similar (at each potential examined) to the corresponding amplitude of the slow component of recovery. 4. The inclusion of 40 mM-iodate ions in the pipette solution slowed the decay of whole-cell sodium currents, as did the extracellular application of venom (10 micrograms ml-1) from the scorpion Leiurus quinquestriatus. Prolonged exposure of the cells to Leiurus venom appeared to increase the steady-state amount of slow inactivation. 5. Records of single-channel sodium currents tended to cluster into records which either did, or did not, contain openings. This apparently non-random behaviour depended on membrane potential, and on the frequency at which the test steps were repeated, in the way expected if it resulted from slow inactivation.(ABSTRACT TRUNCATED AT 400 WORDS)

Delayed end ejection increases isovolumic ventricular relaxation rate in isolated perfused canine hearts.

We sought to determine the ejection variables that are principally responsible for increases in isovolumic ventricular relaxation rate observed with increases in stroke volume. In nine isolated canine hearts, left ventricular ejection was controlled by patterns specially designed to isolate the ejection parameters most critical to isovolumic relaxation rate. When stroke volume was augmented by increases in end-diastolic volume (EDV) with end-systolic volume (ESV) held constant, isovolumic ventricular relaxation rate was unchanged, as gauged by the time constants of single-exponential fits to decaying pressure. In contrast, when ESV was decreased with EDV held constant, isovolumic relaxation time constants decreased significantly, from approximately 100 to 70 msec (protocol I). The important difference in these two situations might have been that the time of end ejection was delayed in the case with faster isovolumic relaxation. To rule out other parameters that may have influenced isovolumic relaxation, ejection velocity was varied in another protocol (protocol II) by either delays in time of the onset of ejection or advances in end-ejection time, always with constant ESV and EDV. Here isovolumic relaxation was progressively slowed as end ejection occurred earlier, whereas isovolumic relaxation rate was insensitive to changes in the onset of ejection, consistent with the unique importance of end ejection to isovolumic relaxation. In fact, our analysis reveals the remarkable finding that changes in isovolumic relaxation time constant produced by either protocol I or protocol II could be related quantitatively to end ejection by a single curve.(ABSTRACT TRUNCATED AT 250 WORDS)

Resolvability of fluorescence lifetime distributions using phase fluorometry

The analysis of the fluorescence decay using discrete exponential components assumes that a small number of species is present. In the absence of a definite kinetic model or when a large number of species is present, the exponential analysis underestimates the uncertainty of the recovered lifetime values. A different approach to determine the lifetime of a population of molecules is the use of probability density functions and lifetime distributions. Fluorescence decay data from continuous distributions of exponentially decaying components were generated. Different magnitudes of error were added to the data to simulate experimental conditions. The resolvability of the distributional model was studied by fitting the simulated data to one and two exponentials. The maximum width of symmetric distributions (uniform, gaussian, and lorentzian), which cannot be distinguished from single and double exponential fits for statistical errors of 1 and 0.1%, were determined. The width limits are determined by the statistical error of the data. It is also shown that, in the frequency domain, the discrete exponential analysis does not uniformly weights all the components of a distribution. This systematic error is less important when probability and distribution functions are used to recover the decay. Finally, it is shown that real lifetime distributions can be proved using multimodal probability density functions. In the companion paper that follows we propose a physical approach, which provides lifetime distribution functions for the tryptophan decay in proteins. In the third companion paper (Alcala, J.R., E. Gratton, and F.J. Prendergast, 1987, Biophys. J., in press) we use the distribution functions obtained to fit data from the fluorescence decay of single tryptophan proteins.

Atom–atom potentials via electron gas theory

The Gordon–Kim electron gas model is used to derive a set of short range repulsive potentials for homoatomic pairs containing atoms up to Kr. These potentials correlate extremely well with experimental van der Waals and nonbonded radii. Single exponential fits to the repulsive potentials are given, and are combined with an approximate dispersion energy term derived from experimental and theoretical atomic dipole polarizabilities and C6 constants to form a set of internally consistent atom–atom potentials of the exp-6 form.

Possibilities and limitations of the time-correlated single photon counting technique: a comparative study of correction methods for the wavelength dependence of the instrument response function

One of the principal reasons for poor fits in time-correlated single photon counting fluorescence decay measurements is the wavelength dependence of the instrument response function. Over the years several correction methods have been developed to account for or correct for this effect. The most recent and widely used procedures are critically compared in this paper. A comparison of the channel shift technique, the pseudoscatterer technique, the excitation pulse-shape mimic technique and the delta function convolution method demonstrates that the latter method is the superior one. This method requires the measurement of the fluorescence decay of a reference compound with single exponential decay kinetics under identical conditions as used for the sample. A modified functional form is used to describe the sample decay law. Since no experimental data are altered in any way, a correct residual analysis is possible. Simulations show that fluorescence lifetimes can be recovered accurately, as long as the reference decay time is sufficiently different from that of the sample. Biexponential decays can be resolved successfully when the two decay times are well separated and different from that of the reference. The difference between the standard normal variates of chi-square for single and double exponential fits can be used to distinguish between mono- and bi-exponential decays. The usefulness of the delta function convolution method is demonstrated by fluorescence decay measurements of a series of different samples with very short decay times. Fluorescence lifetimes as short as 10 ps could be resolved accurately and reliably.