Diurnal rhythms (DR) in gene expression in human oral mucosa: Implications for gender differences in toxicity, response and survival and optimal timing of targeted therapy (Rx)

2507

Background: DR in processes relevant to oncology, including cell cycle, apoptosis, and angiogenesis, have been demonstrated (Nat Rev Cancer 3:350, Cancer Res 63:7277). In rodents, 10–20% of the genome has a 24-hour (h) rhythm in RNA expression (Cell 109:307). A molecular clock consisting of transcription/translation feedback loops of clock-genes controls this rhythmicity (Nat Rev Neurosci 4:649). We studied, for the first time, daily whole genome RNA expression patterning in healthy human volunteers. Methods: RNA samples extracted from oral mucosa biopsies (bx) obtained every 4h over 24h (6 bx) from 5 males (M) and 5 females (F) were subjected to microarray analysis (Affymetrix HG_U133_Plus2 chip, 54,679 transcripts). COSOPT, designed for circadian microarray time series analysis (Methods Enzymol 383: 149) was used to detect genes with significant rhythms (pMMC-Beta = 0.1). Expression patterns were visualized in GeneSpring GX7.3 (Agilent) and validated by real-time PCR and ANOVA. Results: There were 801 and 810 rhythmic genes in M and F respectively, with most genes peaking at 4AM or 4PM in M but at 6AM or 11AM in F. Only 90 rhythmic genes (including core clock-genes and clock controlled genes) were common to M and F, with over 700 genes only rhythmic in M and not in F and vice versa. The profiles of clock- genes and clock-controlled genes were inverted relative to the nocturnal rodent data. There were 75 and 67 rhythmic transcription factor genes in M and F respectively, with 28 common to both M and F. There were 71 rhythmic human cancer genes (Nat Rev Cancer 4:177) with significant gender differences. This group includes rhythmic gene products involved in signaling pathways currently targeted for cancer Rx. Conclusions: We show for the first time a significant gender specific DR in gene expression involving multiple genes of interest in oncology. This may contribute to the documented gender differences in toxicity, response and survival (J Clin Onc 24: 3562, NEJM 353:133), and can inform future trials of optimal timing of antisense Rx and other targeted Rx. The inverted DR in rodents vs. humans has implications for translating rodent data to human Rx trials.

No significant financial relationships to disclose.

Impact of pulsatility on the ensemble orderliness (approximate entropy) of neurohormone secretion

Regular patterns of neurohormone secretion are driven by underlying pulsatile and subordinate (feedback sensitive) dynamics. Measures of time-series orderliness, e.g., the approximate entropy (ApEn) statistic (Pincus SM. Proc Natl Acad Sci 88: 2297-2301, 1991), vividly discriminate pathological and physiological patterns of hormone release. To investigate how specific pulsatility features impact regularity estimates, we have examined the sensitivity of the ApEn metric to systematic variations in the frequency, amplitude, and half-life of simulated neurohormone pulse trains (Veldhuis JD, Carlson ML, and Johnson ML. Proc Natl Acad Sci 84: 7686-7690, 1987) and compared the impact of a high vs. low baseline luteinizing hormone (LH) pattern regularity state mimicking the normal female luteal phase and the young male, respectively. Shortening the interpulse interval length elevated ApEn in both pulsatility models, thereby signifying greater ensemble series irregularity. The frequency sensitivity of ApEn was robust to several complementary renditions of ApEn and to variations in experimental uncertainty, basal (nonpulsatile) LH secretion, and secretory burst amplitude. ApEn rose with increasing hormone half-life, especially in the face of low baseline variability emulated by midluteal LH secretion profiles. High variability of secretory burst amplitude, pulse duration, or interpeak intervals increased ApEn in the more orderly femalelike construct; in the highly irregular malelike LH pulse model, these variability changes had little effect on ApEn. In summary, the ensemble regularity statistic, ApEn, quantifies unequal pattern orderliness in neurohormone pulse trains with minimal dependence on mean pulse amplitude, interpulse baseline, or (subthreshold) sample uncertainty. Thus ApEn monitors changing secretory event frequency and interpulse variability with sensitivity to starting pattern regularity, providing a mechanistic linkage between model evolution and statistical change.

Modulation of dose intensity in aerodigestive tract cancers: strategies to reduce toxicity

Advances in diagnostic and therapeutic radiology and a better understanding of cell biology are being applied in practical ways to modulate treatment morbidity. Conformal radiotherapy targets the cancer precisely and can be combined with new systemically administered radiosensitizers. The successes of conventional chemoradiation programs support continued study of newer ways to deliver systemic radiosensitizing chemotherapy. However, chemoradiation creates a narrower therapeutic window compared to irradiation alone and increased treatment intensity, even with conformal chemoradiation techniques, can potentially result in frequent complications, detrimental treatment delays, and decreased quality of life. Treatment schedules employing a "best tolerated time" modelfor systemic administration of radiosensitizing chemotherapy, based on the concept of chronotolerance, offer attractive ways to address the challenging problem of normal tissue toxicity associated with conformal chemoradiation. This approach may be beneficial in the elderly and those medically unfit to tolerate traditional dose-intense combined-modality schedules. Further evaluation of this concept is warranted, based on existing data.

Seizures induce phase shifts of rat circadian rhythms

RATIONALE

Epileptic seizures may alter neuroendocrinological cycles. Light pulses induce phase shifts in circadian rhythms. Using hippocampal-kindled rats to ensure maximal clinical expression, we determined if seizures likewise induce phase shifts.

METHODS

We monitored the circadian rhythm of temperature (CRT) with intraperitoneal radiotelemetry in rats (n=21) isolated from time cues and light for 3-week trials. Seizures were triggered with hippocampal electrical stimulation at different circadian phases. Optimized, least-error phase shifts were calculated from preictal and postictal CRTs. Induced seizures were referenced to CRT (t(max)=00:00, 24-h circadian cycle).

RESULTS

Phase shifts (individual responses=57) differed across the circadian cycle. Rather than forming a clear phase-response curve, phase shifts were especially variable between 00:00 and 06:00 h.

CONCLUSIONS

This study demonstrates that electrically-induced seizures induce advances and delays in CRT in a phase-dependent fashion but in a pattern different from typical light-induced phase shifts. Disorders of circadian regulation may contribute to some of the altered endogenous cycles associated with epilepsy.

Decrease in beta-cell mass leads to impaired pulsatile insulin secretion, reduced postprandial hepatic insulin clearance, and relative hyperglucagonemia in the minipig

Most insulin is secreted in discrete pulses at an interval of approximately 6 min. Increased insulin secretion after meal ingestion is achieved through the mechanism of amplification of the burst mass. Conversely, in type 2 diabetes, insulin secretion is impaired as a consequence of decreased insulin pulse mass. beta-cell mass is reported to be deficient in type 2 diabetes. We tested the hypothesis that decreased beta-cell mass leads to decreased insulin pulse mass. Insulin secretion was examined before and after an approximately 60% decrease in beta-cell mass achieved by a single injection of alloxan in a porcine model. Alloxan injection resulted in stable diabetes (fasting plasma glucose 7.4 +/- 1.1 vs. 4.4 +/- 0.1 mmol/l; P < 0.01) with impaired insulin secretion in the fasting and fed states and during a hyperglycemic clamp (decreased by 54, 80, and 90%, respectively). Deconvolution analysis revealed a selective decrease in insulin pulse mass (by 54, 60, and 90%) with no change in pulse frequency. Rhythm analysis revealed no change in the periodicity of regular oscillations after alloxan administration in the fasting state but was unable to detect stable rhythms reliably after enteric or intravenous glucose stimulation. After alloxan administration, insulin secretion and insulin pulse mass (but not insulin pulse interval) decreased in relation to beta-cell mass. However, the decreased pulse mass (and pulse amplitude delivered to the liver) was associated with a decrease in hepatic insulin clearance, which partially offset the decreased insulin secretion. Despite hyperglycemia, postprandial glucagon concentrations were increased after alloxan administration (103.4 +/- 6.3 vs. 92.2 +/- 2.5 pg/ml; P < 0.01). We conclude that an alloxan-induced selective decrease in beta-cell mass leads to deficient insulin secretion by attenuating insulin pulse mass, and that the latter is associated with decreased hepatic insulin clearance and relative hyperglucagonemia, thereby emulating the pattern of islet dysfunction observed in type 2 diabetes.

The use of regularity as estimated by approximate entropy to distinguish saltatory growth

A nonlinear dynamics metric, approximate entropy (ApEn), is investigated as a diagnostic method for distinguishing between mathematical models, and the underlying mechanistic hypotheses that purport to describe the same time series experimental observations. ApEn measures the occurrence of pattern regularity within a time series, and is used here to investigate growth patterns in daily length growth. The notion investigated is that ApEn distributions for competing time series patterns expressed as mathematical formulations can be modelled by Monte Carlo and bootstrap methods and compared to the ApEn values for an original experimental data series. If the ApEn values for the different models do not overlap, then it is expected that ApEn can be utilized to distinguish these models and hypotheses, and to provide statistical assessment for the underlying biological patterns in experimental data. The conclusion is that the ApEn metric is successful as a time series diagnostic tool. It is a model-independent statistic that clearly differentiates saltatory growth from slowly varying continuous models of growth and serves to further document the saltatory nature of growth. This is a unique application of approximate entropy, illustrating the broad applicability of ApEn to biological time series, with the specific example of discriminating a saltatory growth process in longitudinal growth data. Future investigations of regularity in longitudinal time series in human biology with ApEn statistics are suggested.

A construct of interactive feedback control of the GH axis in the male

Growth hormone (GH) secretion is controlled by GH-releasing hormone (GHRH), the GH release-inhibiting hormone somatostatin (SRIF), and autofeedback connections. The ensemble network produces sexually dimorphic patterns of GH secretion. In an effort to formalize this system, we implemented a deterministically based autonomous feedback-driven construct of five principal dose-responsive regulatory interactions: GHRH drive of GH pituitary release, competitive inhibition of GH release by SRIF, GH autofeedback via SRIF with a time delay, delayed GH autonegative feedback on GHRH, and SRIF inhibition of GHRH secretion. This formulation engenders a malelike pattern of successive GH volleys due jointly to positive time-delayed feedback of GH on SRIF and negative feedback of SRIF on GH and GHRH. The multipeak volley is explicated as arising from a reciprocal interaction between GH and GHRH during periods of low SRIF secretion. The applicability of this formalism to neuroendocrine control is explored by initial parameter sensitivity analysis and is illustrated for selected feedback-dependent experimental paradigms. The present construct is not overparameterized and does not require an ad hoc pulse generator to achieve pulsatile GH output. Further evolution of interactive constructs could aid in exploring more complex feedback postulates that confer the vivid sexual dimorphism of female GH profiles.

Secretory process regularity monitors neuroendocrine feedback and feedforward signaling strength in humans

The present experiments examine the neuroregulatory hypothesis that the degree of sample-by-sample regularity of hormone output by an interlinked hypothalamopituitary target-organ system monitors the strength of feedback and/or feedforward signaling. To test this postulate and assess its generality, we implemented a total of nine thematically complementary perturbation experiments. In particular, we altered feedback or feedforward signaling selectively in two distinct neuroendocrine systems; namely, the growth hormone (GH) insulin-like growth factor type I (IGF-I) and the luteinizing hormone-testosterone axes. Four experimental paradigms comprised preferential reduction vs. enhancement of IGF-I or testosterone feedback signal strength; and, conversely, five others entailed selective attenuation vs. augmentation of GH-releasing hormone and gonadotropin-releasing hormone feedforward signal intensity. In these independent interventions, quantitation of subordinate (nonpulsatile) secretory pattern reproducibility via the approximate entropy statistic unmasked salient changes (P values typically <10(-3)) in the conditional regularity of serial hormone output with high consistency (96-100%). In particular, approximate entropy quantified degradation of secretory subpattern orderliness under either muted feedback restraint or heightened feedforward drive. Assuming valid interpretation of the biological constraints imposed, these experimental observations coincide with earlier reductionist mathematical predictions, wherein increased irregularity of coupled parameter output mirrors attenuated feedback and/or augmented feedforward coupling within an integrative system.

Orchestrated transcription of key pathways in Arabidopsis by the circadian clock

Like most organisms, plants have endogenous biological clocks that coordinate internal events with the external environment. We used high-density oligonucleotide microarrays to examine gene expression in Arabidopsis and found that 6% of the more than 8000 genes on the array exhibited circadian changes in steady-state messenger RNA levels. Clusters of circadian-regulated genes were found in pathways involved in plant responses to light and other key metabolic pathways. Computational analysis of cycling genes allowed the identification of a highly conserved promoter motif that we found to be required for circadian control of gene expression. Our study presents a comprehensive view of the temporal compartmentalization of physiological pathways by the circadian clock in a eukaryote.

Episodes of severe hypoglycemia in type 1 diabetes are preceded and followed within 48 hours by measurable disturbances in blood glucose

This study quantifies blood glucose (BG) disturbances occurring before and after episodes of severe hypoglycemia (SH). For 6-8 months, 85 individuals with type 1 diabetes and a history of SH (age, 44+/-10 yr; 41 women and 44 men; duration of diabetes, 26+/-11 yr; hemoglobin A1c, 7.7+/-1.1%) used Lifescan One Touch BG meters for self-monitoring three to five times daily and recorded the date and time of SH episodes in diaries. For each subject, the timing of SH episodes was located in the temporal stream of SMBG readings recorded by the meter, and characteristics, including the Low BG index (LBGI), were computed in 24-h increments. In the 24-h period before the SH episode LBGI rose (P < 0.001), average BG was lower (P = 0.001), and BG variance increased (P = 0.001). In the 24 h after SH, LBGI and BGvariance remained elevated (P < 0.001), but average BG returned to baseline. These disturbances disappeared in 48 h. On the basis of LBGI we identified subjects at low, moderate, and high risk of SH, who reported, on the average, 1.7, 3.4, and 7.4 SH episodes (P < 0.005) during the study. In addition, we designed an algorithm that predicted 50% of all SH episodes that occurred in this subject group. We conclude that episodes of SH are preceded and followed by quantifiable BG disturbances, which could be used to devise warnings of imminent SH.

Nuclear scintigraphic assessment of intestinal dysfunction after combined treatment with 9-amino-20(S)-camptothecin (9-AC) and irradiation

PURPOSE

The camptothecins (CPTs) are potent radiosensitizers of malignant tumors in vivo. The extent of normal tissue damage after combined CPT and radiation treatment is unknown. In this article, a jejunal absorption assay with (99m)Tc- pertechnetate (Na[(99m)TcO(4)]) was used to assess C3H/Kam mice given total body irradiation (TBI) of 4 Gy, 6 Gy, and 8 Gy, 2 mg/kg single intramuscular injection of 9-AC or a combination of 2 mg/kg 9-AC + 4 Gy TBI. We also correlated the absorption data with morphologic changes in the jejunal mucosa.

MATERIALS AND METHODS

((99m)TcO(4))(-) absorption from the intestinal lumen into the circulation was studied with dynamic gamma-scintigraphy combined with a multichannel analyzer to record the radiometry data in a time-dependent fashion. Jejunal cross sections were scored for the number of cells per villus and the percentage of apoptotic and mitotic cells in the crypt compartment. The jejunal microcolony assay was used to quantify jejunal crypt survival.

RESULTS

A dose-dependent decrease in the absorption function was observed 3.5 days following TBI. The mean absorption rate was reduced to 89 +/- 16% of control in response to a sublethal 4 Gy TBI and dropped to 47. 5 (9.8% in response to 8 Gy TBI. The mean rate of intestinal absorption was delayed by single sublethal 2 mg/kg 9-AC injection to 62 (11% in comparison with control values. The combination of a single 4 Gy TBI with a 9-AC treatment decreased the ((99m)TcO(4))(-) jejunal absorption in an additive fashion producing absorption lifetime values more than twofold longer than controls. The decrease in ((99m)TcO(4))(-) absorption at 3.5 days after irradiation, 9-AC treatment or the combination of the two agents correlates with the number of cells per villus and the percentage of apoptotic cells in the crypt compartment.

CONCLUSION

Dynamic enteroscintigraphy with (99m)Tc-pertechnetate is a sensitive functional assay for rapid evaluation of radiation and chemotherapy induced intestinal damage. Reduced intestinal absorptive function has a cellular basis and correlates directly with the numbers of cells lost per villus in a treatment-dependent manner.

Dual epileptic foci in a single patient express distinct temporal patterns dependent on limbic versus nonlimbic brain location

How timing information is transferred from the suprachiasmatic nucleus to other regions of the brain to mediate activity, either physiological or pathological, is largely unclear. A patient with medically refractory epilepsy and a well-documented, long-term seizure diary provided a unique means to demonstrate how susceptibility to chronobiological modulation varies with brain region. Evaluation for epilepsy surgery disclosed two independent epileptic foci, one limbic and the other nonlimbic. Seizures from both foci occurred periodically with a dominant period of 24 hours but were out of phase with each other. Temporal lobe seizures occurred maximally in the light portion of the daily light-dark cycle, and parietal lobe seizures occurred nocturnally and out of phase with limbic seizures. These data suggest that neuronal excitation and inhibition, depending on the anatomical system involved in epilepsy, may be differently affected by circadian modulation.

Effects of circadian regulation and rest-activity state on spontaneous seizures in a rat model of limbic epilepsy

PURPOSE

Circadian regulation via the suprachiasmatic nuclei and rest-activity state may influence expression of limbic seizures.

METHODS

Male rats (n = 14) were made epileptic by electrical stimulation of the hippocampus, causing limbic status epilepticus and subsequent seizures. We monitored seizures with intrahippocampal electrodes in 12-12-h light/dark (LD) cycles and in continuous dark (DD). We used radiotelemetry monitoring of activity to measure state and body temperature to determine circadian phase. Cosinor analysis and chi2 tests determined whether seizures occurred rhythmically when plotted by phase. State was defined as inactive or active in 10-min epochs based on whether activity count was below or above a cut-off value validated from video observation.

RESULTS

In LD, the peak seizure occurrence was 14:59 h after circadian temperature peak (95% confidence limit, 13:37-16:19). Phasic seizure occurrence persisted in DD for 14:05 (12:31-15:38), p < 0.0001, against uniform mean distribution. In LD, 14,787 epochs contained 1, 268 seizures; seizures preferentially occurred during inactive epochs (965 observed, 878 expected in proportion to the overall distribution of inactive versus active epochs; p < 0.001). In DD, 20, 664 epochs contained 1,609 seizures; seizures had no preferential occurrence by state (999 observed, 1,025 expected; p = 0.16).

CONCLUSIONS

Limbic seizures occurred with an endogenous circadian rhythm. Seizures preferentially struck during inactivity during entrainment to the light-dark cycle.

Nuclear scintigraphic assessment of radiation-induced intestinal dysfunction

Radiation-induced damage to the intestine can be measured by abnormalities in the absorption of various nutrients. Changes in intestinal absorption occur after irradiation because of loss of the intestinal absorptive surface and a consequent decrease in active transport. In our study, the jejunal absorption of (99m)Tc-pertechnetate, an actively transported gamma-ray emitter, was assessed in C3H/Kam mice given total-body irradiation with doses of 4, 6, 8 and 12.5 Gy and correlated with morphological changes in the intestinal epithelium. The absorption of (99m)Tc-pertechnetate from the intestinal lumen into the circulation was studied with a dynamic gamma-ray-scintigraphy assay combined with a multichannel analyzer to record the radiometry data automatically in a time-dependent regimen. The resulting radioactivity-time curves obtained for irradiated animals were compared to those for control animals. A dose-dependent decrease in absorptive function was observed 3.5 days after irradiation. The mean absorption rate was reduced to 78.8 +/- 9.3% of control levels in response to 4 Gy total-body irradiation (mean +/- SEM tracer absorption lifetime was 237 +/- 23 s compared to 187 +/- 12 s in nonirradiated controls) and to 28.3 +/- 3.7% in response to 12.5 Gy (660 +/- 76 s). The decrease in absorption of (99m)Tc-pertechnetate at 3.5 days after irradiation correlated strongly (P < 0.001) with TBI dose, with the number of cells per villus, and with the percentage of cells in the crypt compartment that were apoptotic or mitotic. A jejunal microcolony assay showed no loss of crypts and hence no measured dose-response effects after 4, 6 or 8 Gy TBI. These results show that dynamic enteroscintigraphy with sodium (99m)Tc-pertechnetate is a sensitive functional assay for rapid evaluation of radiation-induced intestinal damage in the clinically relevant dose range and has a cellular basis.

Hypothalamic neuronal loss and altered circadian rhythm of temperature in a rat model of mesial temporal lobe epilepsy

PURPOSE

Numerous dysfunctions in endogenous hypothalamic function have been associated with mesial temporal lobe epilepsy (MTLE). One endogenous activity is the circadian rhythm of temperature (CRT). In this study we examined whether hypothalamically mediated function is altered in the electrically induced, self-sustained, limbic status epilepticus model of MTLE. We then wished to determine whether there was a structural basis for regulatory alterations.

METHODS

We measured CRT with peritoneal temperature telemetry obtained in light-entrained (LD) and in free-running, constant-dark (DD) conditions. CRT from epileptic and controls of normal animals and kindled animals were quantized by fast Fourier transform-nonlinear least squares analysis to determine rhythmic complexity.

RESULTS

The circadian component of CRT was preserved in all animals. In DD, CRTs of epileptic animals were more complex than those of normal animals. CRT of kindled animals showed no increased complexity after electrically induced seizures. Neuronal density was decreased in regions of the anterior and posterior hypothalamus but not in the suprachiasmatic nuclei from the epileptic rats.

CONCLUSIONS

Alterations in CRT due to the epileptic state were independent of isolated seizures. Altered circadian thermoregulation in epileptic rats corresponded to regional hypothalamic neuronal loss. Structural changes of the hypothalamus may explain alterations in endogenous rhythms in MTLE.

Growth hormone (GH) receptor blockade with a PEG-modified GH (B2036-PEG) lowers serum insulin-like growth factor-I but does not acutely stimulate serum GH

B2036-PEG, a GH receptor (GH-R) antagonist, is an analog of GH that is PEG-modified to prolong its action. Nine mutations alter the binding properties of this molecule, preventing GH-R dimerization and GH action. A potential therapeutic role of B2036-PEG is to block GH action, e.g. in refractory acromegaly. A phase I, placebo-controlled, single rising-dose study was performed in 36 normal young men (ages, 18-37 yr; within 15% ideal body weight). Four groups received a single s.c. injection of either placebo (n = 3 in each group, total n = 12) or B2036-PEG (0.03, 0.1, 0.3, or 1.0 mg/kg; n = 6 each dose). B2036-PEG and GH concentrations were measured 0, 0.25, 0.5, 1, 3, 6, 9, 12, 24, 36, 48, 72, 96, 120, and 144 h after dosing. Serum insulin-like growth factor-I was measured before and 1-7 days after dosing. All doses were well tolerated, with no serious or severe adverse reactions. B2036-PEG, at 1.0 mg/kg, reduced insulin-like growth factor-I by 49 +/- 6% on day 5 (P < 0.001 vs. placebo). GH was measured by two independent methods: 1) modified Nichols chemiluminescence assay (empirically corrected for B2036-PEG cross-reactivity); and 2) direct GH two-site immunoassay, using monoclonal antibodies that did not react with B2036-PEG. There was good agreement between the two methods. GH did not change substantially at any B2036-PEG dose, suggesting that B2036-PEG does not interact with hypothalamic GH-Rs to block short-loop feedback. B2036-PEG may thus block peripheral GH action without enhancing its secretion.

Imbibition, but not release from stratification, sets the circadian clock in Arabidopsis seedlings

Circadian rhythms in the abundance of the CAT2 catalase mRNA were not seen in etiolated seedlings but developed upon illumination. These circadian oscillations were preceded by a rapid and transient induction of CAT2 mRNA abundance that varied strikingly according to the timing (circadian phase) of the onset of illumination. This variation oscillated with a circadian periodicity of approximately 28 hr, indicating that the circadian oscillator is running in etiolated seedlings and regulates (gates) the induction of CAT2 by light. Moreover, because we assayed populations of seedlings, we infer that the individual clocks among populations of etiolated seedlings were synchronized before the onset of illumination. What developmental or environmental signals synchronized the clocks among seedlings? Varying the phase of the onset of illumination relative to release from stratification failed to affect the acute induction of CAT2, indicating that the temperature step from 4 to 22 degrees C associated with release from stratification did not reset the circadian clock. However, the acute induction of CAT2 mRNA varied with time after imbibition, demonstrating that imbibition provides a signal capable of resetting the circadian clock and of synchronizing the clocks among populations of seedlings.

Circadian rhythms of locomotor activity in zebrafish

As part of an effort to characterize the circadian system of the zebrafish, we examined the circadian regulation of locomotor activity in adult males and females. Gross locomotor activity was measured using infrared movement detectors. The effects of light, dark, and temperature on the amplitude, phase, and free-running periods of locomotor rhythms were determined. When zebrafish were maintained in a 12-h light:12 h dark cycle at 25 degrees C, 86% of the fish were most active during the light phase of the cycle. The phases of free-running rhythms measured after transfer of fish from light cycles to constant conditions indicate that this diurnal activity profile reflects entrained circadian rhythmicity. When animals were maintained in constant conditions, the proportion that expressed significant circadian rhythmicity depended on ambient temperature. At 21 degrees C, 73% of the animals were rhythmic in constant darkness, and 65% were rhythmic in constant light. Fewer (28-59%) were rhythmic at 18 degrees, 25 degrees, and 28.5 degrees C. The free-running period of rhythmic animals was not affected by temperature within this range. The average period was shorter in constant light (LL; 12 lx) than in constant darkness (DD) in all but one experiment, and the difference was statistically significant for animals held at 21 degrees C. These data indicate that zebrafish locomotor activity is regulated by a circadian clock that is temperature compensated. Because rhythmicity is most robust at 21 degrees C, this would be the optimal temperature for future studies of the physiological basis of zebrafish behavioral rhythms.

Temporal distribution of partial seizures: comparison of an animal model with human partial epilepsy

Seizures do not often strike randomly but may occur in circadian patterns. We compared daily times of partial seizures determined by continuous electroencephalography among patients with mesial temporal lobe epilepsy (MTLE; n = 64), those with extratemporal lobe (XTLE; n = 26) or lesional temporal lobe epilepsy (LTLE; n = 8), and a rat model similar to MTLE in which rats become epileptic after electrically induced limbic status epilepticus (postlimbic status [PLS]; n = 20). Rats were maintained on a 12-hour light/dark cycle with lights on at 0700 hours. The distributions of seizures were fitted by cosinor analysis to determine time of peak seizure incidence +/- 95% confidence interval (95% CI). The mean fraction +/- SD of seizures recorded during light was 63 +/- 17% in PLS animals and 60 +/- 21% in humans. Peak incidence of seizures for PLS rats (547 seizures) was 1645 (95% CI = 1448,1830) and for MTLE subjects (774 seizures) was 1500 (95% CI = 1324,1636). Seizures from XTLE (465 seizures) and LTLE (48 seizures) did not fit a cosinor model and occurred no more frequently during light than dark. In conclusion, limbic seizures in humans and PLS rats occur more often during light than dark and have similar cosinor daily distributions. The chronological similarity between human MTLE and PLS rat epilepsy suggests that limbic seizure occurrence has a relation to the circadian regulatory system.

Statistically accurate estimation of hormone concentrations and associated uncertainties: methodology, validation, and applications

We describe a data reduction procedure to assign statistically accurate estimates of unknown hormone concentrations, with associated uncertainties, based on experimental uncertainties in sample replicates and the fitted calibration curve. Three mathematical calibration curve functions are considered. The one providing optimal statistical characterization of reference calibrators is chosen for unknown evaluation. Experimental error is addressed by assigning and propagating uncertainty estimates for each measured response (including zero-dose responses) by an empirically determined discrete uncertainty profile and by propagating calibration curve uncertainty. Discrete uncertainty profiles account for both response precision (replicability) and accuracy (deviation from predicted calibration curves) without relying on assumed theoretical response variance-assay response relations. The validity of assigning variable response weighting by this procedure was assessed by Monte Carlo simulations based on chemiluminescence growth hormone calibration curves. Much-improved accuracy and estimated precision are achieved for unknown hormone concentrations, particularly extremely low concentrations, by using this variable response weighting procedure.

Long-term monitoring of circadian rhythms in c-fos gene expression from suprachiasmatic nucleus cultures

BACKGROUND

The AP-1 family of transcription factors has been implicated in the control of the expression of many genes in response to environmental signals. Previous studies have provided temporal profiles for c-fos expression by taking measurements from many animals at several points in time, but these studies provide limited information about dynamic changes in expression. Here, we have devised a method of continuously measuring c-fos expression.

RESULTS

A transgenic mouse line expressing the human c-fos promoter linked to the firefly luciferase reporter gene (fos/luc) was generated to continuously monitor c-fos gene expression. A second transgenic mouse line expressing luciferase under the control of the cytomegalovirus promoter (CMV/luc) served as a control. Luminescence originating from identifiable brain regions was imaged from fos/luc brain slice cultures. Expression of the fos/luc transgene accurately reflected transcriptional responses of the endogenous c-fos gene. Dynamic changes in fos/luc expression in suprachiasmatic nuclei (SCN) explant cultures were monitored continuously, and luminescence showed almost 24 hour rhythms lasting up to five circadian cycles. In contrast, bioluminescence monitored from CMV/luc SCN explant cultures was not rhythmic.

CONCLUSION

The fos/luc transgenic mouse will be useful for long-term, non-invasive monitoring of c-fos transcriptional responses to the changing cellular environment. Circadian rhythms in c-fos expression can be monitored non-invasively in real time from the SCN, clearly demonstrating that c-fos transcription is regulated by the circadian clock.

Aging alters feed-forward and feedback linkages between LH and testosterone in healthy men

To discern the effect of aging on coordinate luteinizing hormone (LH) and testosterone secretion, we sampled healthy older men (age 62-74 yr, n = 11) and young controls (age 21-34 yr, n = 13) every 2.5 min overnight. Deconvolution analysis and cross-correlation were used to relate serum LH concentrations to calculated testosterone secretion rates (feed-forward stimulation), as well as serum testosterone concentrations to computed LH secretion rates (feedback inhibition). Despite statistically similar mean serum LH and testosterone concentrations in the young and older men, older individuals had diminished feed-forward stimulation of LH concentrations on calculated testosterone secretion rates, as well as delayed feedback inhibition of testosterone concentrations on computed LH secretion rates.

Effects of synergistic signaling by phytochrome A and cryptochrome1 on circadian clock-regulated catalase expression

Persistent oscillation in constant conditions is a defining characteristic of circadian rhythms. However, in plants transferred into extended dark conditions, circadian rhythms in mRNA abundance commonly damp in amplitude over two or three cycles to a steady state level of relatively constant, low mRNA abundance. In Arabidopsis, catalase CAT3 mRNA oscillations damp rapidly in extended dark conditions, but unlike catalase CAT2 and the chlorophyll a/b binding protein gene CAB, in which the circadian oscillations damp to low steady state mRNA abundance, CAT3 mRNA oscillations damp to high steady state levels of mRNA abundance. Mutational disruption of either phytochrome- or cryptochrome-mediated light perception prevents damping of the oscillations in CAT3 mRNA abundance and reveals strong circadian oscillations that persist for multiple cycles in extended dark conditions. Damping of CAT3 mRNA oscillations specifically requires phytochrome A but not phytochrome B and also requires the cryptochrome1 blue light receptor. Therefore, we conclude that synergistic signaling mediated through both phytochrome A and cryptochrome1 is required for damping of circadian CAT3 mRNA oscillations in extended dark conditions.

Quantitative analysis of Drosophila period gene transcription in living animals

To determine the in vivo regulatory pattern of the clock gene period (per), the authors recently developed transgenic Drosophila carrying a luciferase cDNA fused to the promoter region of per. They have now carried out noninvasive, high time-resolution experiments allowing high-throughput monitoring of circadian bioluminescence rhythms in individual living adults for several days. This immediately solved several problems (resulting directly from individual asynchrony within a population) that have accompanied previous biochemical experiments in which groups of animals were sacrificed at each time point. Furthermore, the authors have developed numerical analysis methods for automatically determining rhythmicity associated with bioluminescence records from single flies. This has revealed some features of per gene transcription that were previously unappreciated and provides a general strategy for the analysis of rhythmic time series in the study of molecular rhythms.

Circadian rhythms in mouse suprachiasmatic nucleus explants on multimicroelectrode plates

The suprachiasmatic nucleus (SCN) of the mammalian hypothalamus functions as a circadian pacemaker. This study used multimicroelectrode plates to measure extracellular action potential activity simultaneously from multiple sites within the cultured mouse SCN. Neurons within the isolated mouse SCN expressed a circadian rhythm in spontaneous firing rate for weeks in culture.

Conditional circadian dysfunction of the Arabidopsis early-flowering 3 mutant

Photoperiodic responses, such as the daylength-dependent control of reproductive development, are associated with a circadian biological clock. The photoperiod-insensitive early-flowering 3 (elf3) mutant of Arabidopsis thaliana lacks rhythmicity in two distinct circadian-regulated processes. This defect was apparent only when plants were assayed under constant light conditions. elf3 mutants retain rhythmicity in constant dark and anticipate light/dark transitions under most light/dark regimes. The conditional arrhythmic phenotype suggests that the circadian pacemaker is intact in darkness in elf3 mutant plants, but the transduction of light signals to the circadian clock is impaired.

Novel features of drosophila period Transcription revealed by real-time luciferase reporting

The rapid turnover of luciferase and the sensitive, non-invasive nature of its assay make this reporter gene uniquely situated for temporal gene expression studies. To determine the in vivo regulatory pattern of the Drosophila clock gene period (per), we generated transgenic strains carrying a luciferase cDNA fused to the promoter region of the per gene. This has allowed us to monitor circadian rhythms of bioluminescence from pacemaker cells within the head for several days in individual living adults. These high time-resolution experiments permitted neuronal per transcription and opens the door to vastly simplified experiments in general chronobiology and studies of temporally regulated transcription in a wide range of experimental systems.

Thermodynamics of interaction of the fusion-inhibiting peptide Z-D-Phe-L-Phe-Gly with dioleoylphosphatidylcholine vesicles: direct calorimetric determination

The binding of the fusion-inhibiting peptide Z-D-Phe-L-Phe-Gly to unilamellar lipid vesicles of dioleoylphosphatidylcholine (DOPC) was investigated by isothermal titration calorimetry (ITC). The peptide Z-D-Phe-L-Phe-Gly is known to inhibit fusion of myxo- and paramyxoviruses with cells as well as cell-cell and vesicle-vesicle fusion in model systems. Calorimetric titrations conducted over a range of temperatures permitted characterization of the thermodynamics of the interaction of Z-D-Phe-L-Phe-Gly with model DOPC lipid membranes. Simultaneous global analysis of 15 ITC binding curves acquired at four different temperatures allowed determination of the equilibrium site association constant (K), stoichiometry of binding (n), binding enthalpy change (delta H), and heat capacity change of binding (delta Cp) in a single set of experiments. The binding affinity and enthalpy change per mole of DOPC bound at 25 degrees C was log K = 2.463 +/- 0.075 and delta H = -1.07 +/- 0.12 kcal/mol DOPC while the binding heat capacity change per mole of DOPC bound was delta Cp = -20.3 +/- 2.8 cal/(K.mol DOPC) with a temperature dependence (from 10-45 degrees C) of d(delta Cp)/dT = 0.37 +/- 0.18 cal/(K2.mol DOPC). A temperature-independent binding stoichiometry was determined to be n = 5.56 +/- 0.33 DOPC molecules per Z-D-Phe-L-Phe-Gly. A comparison of these results with previous peptide-lipid binding studies is discussed as is their relevance to a current model of the interaction of fusion-inhibiting peptides with phospholipid membranes.

The regulation of circadian period by phototransduction pathways in Arabidopsis

Transgenic Arabidopsis plants expressing a luciferase gene fused to a circadian-regulated promoter exhibited robust rhythms in bioluminescence. The cyclic luminescence has a 24.7-hour period in white light but 30- to 36-hour periods under constant darkness. Either red or blue light shortened the period of the wild type to 25 hours. A phytochrome-deficient mutation lengthened the period in continuous red light but had little effect in continuous blue light, whereas seedlings carrying mutations that activate light-dependent pathways in darkness maintained shorter periods in constant darkness. These results suggest that both phytochrome- and blue light-responsive photoreceptor pathways control the period of the circadian clock.

Enhanced sensitivity growth hormone (GH) chemiluminescence assay reveals lower postglucose nadir GH concentrations in men than women

Modifications were made to a commercially available human (h) GH chemiluminescence assay (Nichols Luma Tag hGH assay), which improved its sensitivity to 0.002 micrograms/L. The results of this assay had a high correlation with those of the Nichols hGH immunoradiometric assay (IRMA; r = 0.91; P < 0.001). The addition of recombinant hGH-binding protein (0.1-10 nmol/L) to standards and serum samples caused a dose-responsive reduction in measured GH in both the chemiluminescence assay and the IRMA; at physiological concentrations of hGH-binding protein, a 10-20% reduction was observed. Fifteen normal young adults (nine men and six women) underwent a standard 100-g oral glucose tolerance test, and plasma GH was measured from 30 min before until 5 h after glucose ingestion. GH was measurable in all samples with the chemiluminescence assay, but fell below the sensitivity of the IRMA in 59% of the samples. There was no difference between baseline or peak glucose levels in male and female subjects, but serum GH concentrations (mean +/- SD) measured by the enhanced sensitivity chemiluminescence assay were lower in male than female subjects at both baseline (0.12 +/- 0.08 vs. 2.3 +/- 2.3 micrograms/L; P < 0.01) and the postglucose GH nadir (0.029 +/- 0.014 vs. 0.25 +/- 0.23 micrograms/L; P < 0.01). The high correlation between baseline and nadir GH (r = 0.82; P < 0.001) and the equivalent fractional decline in mean GH levels in men and women after glucose administration (67 +/- 17% vs. 84 +/- 8%; P = 0.06) suggest that the lower GH levels in men after glucose treatment are due to lower baseline values and not to a greater suppressive effect of glucose.

Two-dimensional differential scanning calorimetry: simultaneous resolution of intrinsic protein structural energetics and ligand binding interactions by global linkage analysis

A general theoretical development for the design and analysis of two-dimensional thermal stability surfaces of proteins is presented. The surfaces are generated from multiple excess heat capacity profiles (<delta Cp> vs T) obtained at varying concentrations of an interacting ligand. The energetics of both the intrinsic protein stability and the protein-ligand interaction are simultaneously resolved by employing statistical thermodynamic models in global linkage analysis. This formalism allows resolution of the intrinsic protein folding-unfolding parameters (enthalpy, entropy, and heat capacity changes) as well as the ligand interaction parameters (binding stoichiometry, enthalpy, entropy, and heat capacity changes). The theory has been applied to the case of ribonuclease A and its interaction with cytidine-2'-monophosphate. The accuracy of the thermodynamic parameters obtained by this approach compares within error with those parameters that can be obtained by direct measurements.

Role of sn-1-saturated,sn-2-polyunsaturated phospholipids in control of membrane receptor conformational equilibrium: effects of cholesterol and acyl chain unsaturation on the metarhodopsin I in equilibrium with metarhodopsin II equilibrium

The effect of phospholipid bilayer acyl chain packing free volume on the equilibrium concentration of the form of photolyzed rhodopsin which initiates visual signal transduction, metarhodopsin II (meta II), is examined in reconstituted systems formed from the saturated phospholipid dimyristoylphosphatidylcholine (DMPC) and in the polyunsaturated phospholipid sn-1-palmitoyl-sn-2-arachidonoylphosphatidylcholine (PAPC) with and without 30 mol% cholesterol. The extent of meta II formation is determined from both flash photolysis measurements and rapidly acquired absorbance spectra. Equilibrium and dynamic properties of the lipid bilayer are characterized by the dynamic fluorescence properties of 1,6-diphenyl-1,3,5-hexatriene (DPH). DPH orientational properties are characterized by fv, a parameter which reflects the volume available for probe reorientation in the bilayer, relative to that available in an unhindered, isotropic environment [Straume, M., & Litman, B. J. (1987) Biochemistry 26, 5121-5126]. The metarhodopsin I in equilibrium with meta II equilibrium constant, Keq has a linear relationship with fv for rhodopsin in PAPC vesicles with and without cholesterol as well as for rhodopsin in DMPC vesicles, and these two correlation lines have different slopes. The correlations between Keq and fv in PAPC and DMPC systems are compared with a similar correlation in the native rod outer segment disk membrane and one reported previously in an egg phosphatidylcholine (egg PC) system [Mitchell, D. C., Straume, M., Miller, J. L., & Litman, B. J. (1990) Biochemistry 29, 9143-9149].(ABSTRACT TRUNCATED AT 250 WORDS)

Modulation of metarhodopsin formation by cholesterol-induced ordering of bilayer lipids

The effect of lipid ordering on the kinetics and extent of metarhodopsin II (meta II) formation was evaluated in bovine rhodopsin which had been reconstituted into phosphatidylcholine vesicles containing 0, 15, and 30 mol% cholesterol. The rate of establishment of the dynamic equilibrium between metarhodopsin I (meta I) and the two kinetically distinguished forms of meta II in the branched meta II model [meta IIfast and meta IIslow; Straume, M., Mitchell, D. C., Miller, J. L., & Litman, B. J. (1990) Biochemistry (preceding paper in this issue)] is derived from kinetic measurements of rhodopsin photolysis in these vesicle systems at several temperatures. Values of the meta I in equilibrium with meta IItotal equilibrium constant, Keq, are calculated from the derived model-dependent rate constants, and are shown to be equivalent to those derived from rapidly acquired absorbance spectra. The presence of 30 mol% cholesterol reduces Keq by approximately 50% between 10 and 37 degrees C. Analysis of the model-dependent parameters in terms of delta H and delta S reveals that cholesterol raises the free energy of meta IIslow, relative to meta I, by increasing delta H whereas it raises the relative free energy of meta IIfast by making delta S meta IIfast relative to meta I less positive. The reduction in Keq by both temperature and cholesterol is found to be directly correlated with a parameter that reflects the free volume available for molecular motion in the hydrophobic core of the bilayer [Straume, M., & Litman, B. J. (1988) Biochemistry 27, 7723-7733].(ABSTRACT TRUNCATED AT 250 WORDS)

Interconversion of metarhodopsins I and II: a branched photointermediate decay model

Flash photolysis experiments designed to monitor the establishment of the metarhodopsin I to metarhodopsin II equilibrium are interpreted according to a branched model in which two spectrally indistinguishable but kinetically distinguishable forms of metarhodopsin II are postulated to exist in equilibrium with a common pool of metarhodopsin I. This interpretation arises from the consistent requirement for at least three exponentials for a valid description of the observed growth of absorbance at 380 nm following bleaching of bovine rhodopsin in rod outer segment disk membranes. Analysis of the 380-nm transient absorbance data permitted direct determination of the five physically interpretable individual rate constants of the model. This analysis represents a more explicit interpretation of kinetic data than that employed in earlier experiments of this kind, which involved estimating only apparent rates and apparent amplitudes of discrete multiexponential functions. The 380-nm absorbance contributions of all relevant species contributing to the observed dynamic absorbance change were accounted for simultaneously during nonlinear least-squares estimation of the model rate parameters. Analysis of deconvoluted equilibrium spectra acquired from samples identical with those used in the kinetics experiments confirmed the metarhodopsin I-metarhodopsin II equilibrium constants, Keq, derived from the dynamic analyses. It is shown that Keq varies from 1.28 at 10 degrees C to 7.3 at 37 degrees C and that approximately 90% of the metarhodopsin II present is in the form of metarhodopsin IIslow over the temperature range 10-37 degrees C. A physical interpretation of this decay model is discussed in the context of a distribution of metarhodopsin II structural and energetic states.

Resolvability of free energy changes for oxygen binding and subunit association by human hemoglobin

Probability distributions of the free energy changes for oxygen binding, subunit association, and quaternary enhancement by human hemoglobin were obtained from Monte Carlo simulations performed on two independent sets of variable protein concentration equilibrium oxygen-binding data. Uncertainties in unliganded and fully liganded dimer to tetramer association free energy changes (0 delta G'2 and 4 delta G'2) were accounted for in the simulations. Distributions of the dimer to tetramer association free energy changes for forming singly and triply liganded tetramers (1 delta G'2 and 3 delta G'2) are well defined and quite symmetric, whereas that for forming doubly liganded tetramers (2 delta G'2) is poorly defined and highly asymmetric. The distribution of the dimer stepwise oxygen-binding free-energy change (delta g'2i) is well defined and quite symmetric as are those of the tetramer stepwise oxygen-binding free-energy changes for binding the first and last oxygens to tetramers (delta g'41 and delta g'44). Distributions of the intermediate tetramer stepwise oxygen-binding free-energy changes (delta g'42 and delta g'43) are poorly defined and highly asymmetric, but are compensatory in that their sum (delta g'4[2 + 3]) is again well defined and nearly symmetric. Distributions of the free energy changes corresponding to the tetramer product Adair oxygen binding constants (delta G'4i) are well defined and quite symmetric for i = 1, 3, 4 but not for i = 2. The distribution of delta g'44 - delta g'2i (the quaternary enhancement free energy change) is relatively narrow, nearly symmetric, and confined to the negative free-energy domain. This suggests that the quaternary enhancement free energy change (a) may be resolved with good confidence from this data and (b) is finite and negative under the conditions of these experiments. Our results also suggest two different four-state combinatorial switch models that provide accurate characterization of hemoglobin's functional behavior.

Equilibrium and dynamic bilayer structural properties of unsaturated acyl chain phosphatidylcholine-cholesterol-rhodopsin recombinant vesicles and rod outer segment disk membranes as determined from higher order analysis of fluorescence anisotropy decay

Limited-frequency phase-modulation fluorometry of diphenylhexatriene (DPH) and trimethylammonium-diphenylhexatriene (TMA-DPH) was used to characterize the equilibrium and dynamic lipid structural properties of (1) reconstituted egg phosphatidylcholine (egg PC)-rhodopsin vesicles varying in rhodopsin content from 0 to approximately 1 mol %, (2) reconstituted PC-cholesterol-rhodopsin vesicles containing approximately 1 mol % rhodopsin and 0, approximately 15, or approximately 30 mol % cholesterol with egg PC, DOPC (di-18:1-PC), or PAPC (16:0,20:4-PC) as the phospholipid constituent, and (3) native bovine rod outer segment disk membranes. Experiments were conducted at 37, 25, 15, and 5 degrees C. Fluorescence lifetime analysis was performed by fitting the data to a constrained, discrete, biexponential model. Rotational depolarization properties were considered by a model requiring a single rotational diffusion coefficient and capable of producing orthogonal, bimodal orientational distributions for DPH and unimodal distributions for TMA-DPH [Straume, M., & Litman, B. J. (1987) Biochemistry 26, 5113-5120]. Unbleached rhodopsin reduced mean fluorophore lifetimes in proportion to the amount of protein present in PC vesicles as a result of probe-to-retinal energy transfer by (1) redistributing the relative lifetime contributions in favor of the short lifetime population and (2) reducing the lifetimes of each derived population. Lifetimes were increased by cholesterol and by reduction of the temperature, but the relative proportions of derived short- and long-lifetime populations were not affected. TMA-DPH lifetimes were more sensitive (in a relative manner) than were those of DPH. These observations are interpreted in terms of cholesterol and reduced temperature each inhibiting water penetrability into these bilayers, with a greater effect occurring in the headgroup and interfacial regions (probed by TMA-DPH) than in the hydrophobic bilayer interior (probed by DPH). Diunsaturated DOPC-rhodopsin recombinants were more resistant to temperature-dependent lifetime changes than were mixed-chain egg PC or PAPC vesicles. This suggests less favorable interaction of rhodopsin with diunsaturated PCs than with mixed-chain PCs. Lifetimes in disk membranes exhibited this same temperature dependence although DPH in disks had lifetimes longer than those seen in recombinant vesicles. TMA-DPH lifetimes in disks were more similar to those observed in cholesterol-containing recombinants. It would therefore appear that the large proportion of small, charged (at pH 7) phosphatidylethanolamine and phosphatidylserine headgroups present in disks reduces water penetrability into the d

Three-state combinatorial switch models as applied to the binding of oxygen by human hemoglobin

We have generated a series of all 6561 unique, discrete three-state combinatorial switch models to describe the partitioning of the cooperative oxygen-binding free change among the 10 variously ligated forms of human hemoglobin tetramers. These models were inspired by the experimental observation of Smith and Ackers that the cooperative free energy of the intersubunit contact regions of the 10 possible ligated forms of human hemoglobin tetramers can be represented by a particular distribution of three distinct energy levels [Smith, F. R., & Ackers, G. K. (1985) Proc. Natl. Acad. Sci. U.S.A. 82, 5347-5351]. A statistical thermodynamic formulation accounting for both dimer-tetramer equilibria and ligand binding properties of hemoglobin solutions as a function of oxygen and protein concentrations was utilized to exhaustively test these thermodynamic models. In this series of models each of the 10 ligated forms of the hemoglobin tetramer can exist in one, and only one, of three possible energy levels; i.e., each ligated form was assumed to be associated with a discrete energy state. This series of models includes all possible ways that the 10 ligation states of hemoglobin can be distributed into three distinct cooperative energy levels. The mathematical models, as presented here, do not permit equilibria between energy states to exist for any of the 10 unique ligated forms of hemoglobin tetramers. These models were analyzed by nonlinear least-squares estimation of the free energy parameters characteristic of this statistical thermodynamic development.(ABSTRACT TRUNCATED AT 250 WORDS)

Influence of cholesterol on equilibrium and dynamic bilayer structure of unsaturated acyl chain phosphatidylcholine vesicles as determined from higher order analysis of fluorescence anisotropy decay

The influence of cholesterol on equilibrium and dynamic bilayer structure in minimally to highly unsaturated phosphatidylcholine (PC) vesicles has been examined by characterization of the dynamic fluorescence properties of 1,6-diphenyl-1,3,5-hexatriene (DPH) and 1-[4-(trimethylammonio)phenyl]-6-phenyl-1,3,5-hexatriene (TMA-DPH). Large, unilamellar egg PC, palmitoyloleoyl-PC (POPC), dioleoyl-PC (DOPC), palmitoylarachidonoyl-PC (PAPC), and palmitoyldocosahexaenoyl-PC (P-22:6-PC) vesicles containing no cholesterol or approximately 15 or 30 mol % cholesterol have been examined. Equilibrium and dynamic DPH orientational properties were analyzed according to an orthogonal, bimodal orientational distribution function [Straume, M., & Litman, B.J. (1987) Biochemistry (preceding paper in this issue)]. The same mathematical formalism was applied to TMA-DPH except that probe orientational probability was permitted only in the distribution peak aligned parallel to the bilayer normal. TMA-DPH fluorescence lifetimes were consistently increased by incorporation of cholesterol into these vesicles. Greater acyl chain unsaturation and increasing temperature each promoted reduction of lifetimes in the presence or absence of cholesterol. DPH lifetimes were much less sensitive than those of TMA-DPH to changes in composition or temperature. This behavior is consistent with reduced water penetrability into liquid-crystalline bilayers as cholesterol content is increased and as acyl chain unsaturation and temperature are reduced. Cholesterol also induces substantial equilibrium ordering of the bilayer both at the hydrophobic core and at the bilayer-water interface. DPH orientational distributions were shifted in favor of alignment parallel to the acyl side chains. The distributions of both probes were narrowed in response to incorporation of cholesterol.(ABSTRACT TRUNCATED AT 250 WORDS)