Fact-checked by Grok 2 weeks ago

Phase response curve

A phase response curve (PRC) is a graphical tool that quantifies the transient shift in the cycle period of an oscillator—such as a biological rhythm—induced by a perturbation applied at varying phases of its oscillation, plotting the phase advance or delay against the timing of the stimulus.^[1] These curves are fundamental in analyzing how external or internal stimuli entrain or reset oscillatory dynamics, with the phase typically normalized from 0 to 1 over the unperturbed period T_0, and the shift calculated as \Delta \phi = (T_0 - T_1)/T_0, where T_1 is the perturbed cycle length.^[2] PRCs are classified into types based on their shape and the oscillator's bifurcation structure: Type 0 PRCs exhibit large shifts with a slope near -1, enabling strong resetting; Type 1 PRCs show smaller, primarily monotonic shifts with a slope near 0; and Type II PRCs are biphasic, allowing both advances and delays.^[1] In circadian biology, PRCs describe how zeitgebers like light or melatonin adjust the ~24-hour internal clock to environmental cycles, with phase delays typically occurring in the early subjective night and advances in the late subjective night, often featuring a "dead zone" of minimal response during the day.^[2] This framework, pioneered by Arthur Winfree in the 1970s, underpins models of entrainment, where repeated stimuli gradually align the oscillator to an external period, and is experimentally derived through protocols involving stimulus application followed by rhythm assays over multiple cycles.^[2] For instance, human PRCs to light exposure reveal peak sensitivity around the body's core temperature minimum, informing chronotherapy for sleep disorders.^[2] In neuroscience, PRCs characterize the response of neuronal oscillators to perturbations like synaptic inputs, measuring shifts in spike timing to predict synchronization in neural networks.^[3] Type I neuronal PRCs, common in integrators like Hodgkin-Huxley models near saddle-node bifurcations, produce only advances and hinder synchronization under excitatory coupling, while Type II PRCs, seen in resonator neurons near Hopf bifurcations, enable phase locking and oscillatory patterns in circuits.^[3] Methods to compute neuronal PRCs include direct perturbation trials, adjoint-based linearization for infinitesimal responses, and spike-triggered averages, requiring datasets of hundreds to thousands of spikes for accuracy.^[3] Beyond biology, PRCs apply to engineering and physics for controlling coupled oscillators, such as in laser arrays or power grids.^[1]

Fundamentals

Definition and Principles

A phase response curve (PRC) is a function that quantifies the transient phase shift—either an advance or delay—in the periodic cycle of a biological oscillator resulting from a perturbation of fixed duration and intensity delivered at different phases of the oscillation.^[1] This shift is measured as the difference between the old phase (before the perturbation) and the new phase (after the perturbation returns to the limit cycle), providing a map of how the timing of stimuli influences the oscillator's rhythm.^[2] The phase itself serves as a normalized time variable within the oscillator's cycle, commonly scaled from 0 to 1 (or equivalently 0 to $2\pi) to represent progress along the limit cycle, with phase 0 typically marking a reference point such as the onset or peak of the oscillation.^[1] Advances occur when the new phase exceeds the old phase (shifting events earlier relative to the unperturbed cycle), while delays happen when the new phase is less than the old phase (postponing events).^[2] PRCs are categorized into types based on their shape and the strength of the stimulus. Type 0 PRCs, elicited by strong perturbations, produce large, nonlinear phase shifts that can approach a full cycle in magnitude, often including dead zones of negligible response and a discontinuous profile near a slope of -1 in the phase transition curve.^[2] For weak stimuli, Type 1 PRCs yield small shifts with a continuous, nearly flat slope near 0, meaning the magnitude varies little with the phase of application. For weak perturbations, an additional distinction exists between Type 1 (monotonic) and Type II (biphasic) PRCs, reflecting different underlying bifurcation structures.^[2]^[1] This classification, originally developed in the context of limit-cycle oscillators, highlights how PRC shape reflects the underlying dynamics, such as the geometry of the oscillator's state space.^[1] Key principles of PRCs involve isochrons, which are surfaces of constant asymptotic phase in the oscillator's phase space; perturbations map the system to a new isochron, determining the resulting phase reset.^[1] For weak perturbations, the infinitesimal PRC approximates the first-order phase response as the integral of the perturbation weighted by the phase sensitivity, enabling analysis of subtle interactions like synchronization in ensembles of oscillators.^[4] To illustrate, consider a simple sinusoidal oscillator where the rhythm follows a smooth wave; a brief perturbation early in the upward swing might delay the peak by compressing the trajectory toward the limit cycle, while one near the downward swing could advance it by accelerating the return, demonstrating how PRCs reveal timing sensitivities without changing the cycle's intrinsic frequency.^[1] These concepts underpin applications in diverse biological systems, including circadian rhythms and neuronal firing patterns.^[2]

Mathematical Formulation

The phase response curve (PRC) quantifies the phase shift \Delta \phi induced by a perturbation of magnitude \epsilon applied at phase \phi in a limit cycle oscillator, approximated for small \epsilon by the linear relation \Delta \phi = Z(\phi) \epsilon \mod 1, where Z(\phi) is the infinitesimal PRC function representing the sensitivity to infinitesimal perturbations.^[1] This formulation assumes weak perturbations where higher-order effects are negligible, and the phase \phi is normalized such that the unperturbed period corresponds to \phi \in [0, 1).^[5] For finite perturbations, the PRC deviates from linearity, requiring direct computation of \Delta \phi as the difference between perturbed and unperturbed return times to a reference phase, but the infinitesimal case serves as a foundational approximation for analytical tractability.^[1] In limit cycle oscillators governed by \dot{x} = f(x), where x is the state vector and f defines the dynamics on the stable limit cycle x_0(t), the infinitesimal PRC Z(\phi) is derived via the adjoint method. The adjoint equation in time coordinates is \frac{dZ}{dt} = - \left( \frac{\partial f}{\partial x}(x_0(t)) \right)^T Z(t), with periodic boundary conditions Z(0) = Z(T) (where T is the period) and normalization Z(t)^T \dot{x}_0(t) = 1 to ensure the phase advance is \Delta \phi = \int Z(t)^T g(t) \, dt for a perturbation g(t).^[5] Transforming to phase coordinates \phi = \theta(t), the equation becomes \frac{dZ}{d\phi} = - \left( \frac{\partial f}{\partial x}(x_0(\phi)) \right)^T Z(\phi) + \lambda, where \lambda is a Lagrange multiplier enforcing the normalization \int_0^1 Z(\phi)^T \frac{\partial x_0}{\partial \phi}(\phi) \, d\phi = 1.^[6] The phase resetting map iterates the dynamics under periodic forcing, given by \phi_{n+1} = \phi_n + \Delta \phi(\phi_n) \mod 1, which determines stable phase-locked states when the map has fixed points.^[1] For synchronization analysis, this map reveals Arnold tongues in the parameter space of forcing frequency and amplitude, regions where the oscillator locks to rational frequency ratios (e.g., 1:1 entrainment), with tongue width scaling as the coupling strength for weak interactions.^[1] Numerical computation of PRCs often employs direct simulation: solve the perturbed ODE \dot{x} = f(x) + \epsilon g(\phi) from initial phase \phi, measure the asymptotic phase shift by interpolation to the isochron or Poincaré section, and repeat across \phi \in [0,1) to construct Z(\phi).^[1] This method is versatile for arbitrary models but computationally intensive compared to the adjoint approach, which integrates the linear adjoint equation once along the limit cycle.^[5] In weakly coupled regimes, the infinitesimal PRC suffices for phase reduction to \dot{\phi} = \omega + \epsilon Z(\phi) \cdot h(\phi), but strong coupling necessitates higher-order PRCs incorporating quadratic and beyond terms for accurate phase shifts.^[1] PRC types distinguish behaviors: Type 1 PRCs produce only advances or delays (associated with saddle-node bifurcations), while Type II PRCs enable both, leading to complex dynamics like relative coordination in networks.^[1]

Circadian Rhythms

Light-Induced PRC

In mammalian circadian systems, light serves as the principal zeitgeber, inducing phase shifts in the rhythm of the suprachiasmatic nucleus (SCN), the central circadian pacemaker. These shifts follow a characteristic pattern: exposure during the early subjective night typically produces phase delays, while light in the late subjective night elicits phase advances, with minimal effects observed during the subjective day. This type 0 phase response curve (PRC) to bright light pulses aligns the endogenous clock with environmental day-night cycles, preventing drift in free-running conditions.^[7] Pioneering empirical studies in the 1980s established the human light PRC using controlled bright light exposures. In one seminal investigation, 3-hour pulses of approximately 10,000 lux administered at various circadian phases resulted in phase shifts of 1 to 2 hours, with delays up to 1.5 hours in the early evening and advances of similar magnitude in the early morning, relative to the core body temperature minimum. Similar patterns emerge in rodents, such as Syrian hamsters, where light pulses during early subjective night delay the rhythm by 0.5 to 1.5 hours, and late-night exposures advance it comparably. Across species, the PRC exhibits a biphasic structure with distinct advance and delay zones flanking a dead zone during the subjective day, though the absolute timing differs: nocturnal animals like rodents show delays shortly after activity onset, while diurnal species, such as the rodent Arvicanthis niloticus, display an inverted profile relative to their daytime activity, yet retain the core advance-delay morphology.^[7] The photic signal originates from intrinsically photosensitive retinal ganglion cells (ipRGCs), a subset of retinal ganglion cells that express the photopigment melanopsin, conferring sensitivity to short-wavelength blue light around 480 nm. These ipRGCs project directly to the SCN via the retinohypothalamic tract, conveying light information that modulates clock gene expression, such as Per1 and Per2, to drive entrainment. Unlike rods and cones, which contribute indirectly, melanopsin-mediated signaling is essential for sustained phase shifts, as demonstrated in melanopsin-knockout mice that exhibit diminished light-induced resetting. The magnitude of light-induced phase shifts depends on stimulus parameters, including intensity and duration. Phase responses scale logarithmically with light intensity, such that a tenfold increase from 100 to 1,000 lux can double the shift size, though responses saturate above ~10,000 lux; even low intensities (50-100 lux) elicit measurable effects if timed appropriately. Duration also influences efficacy: brief pulses (e.g., 30 minutes) produce smaller shifts than extended exposures (3-6 hours), but the relationship is nonlinear, with diminishing returns beyond 2-3 hours. Timing is referenced to internal circadian markers, such as dim light melatonin onset (DLMO) in humans, where light shortly after DLMO maximizes delays.^[9]^[9] In clinical applications, light-induced PRCs inform chronotherapy protocols to mitigate circadian misalignment from jet lag or shift work. For eastward jet lag (phase advance required), morning bright light exposure (2,000-10,000 lux for 1-2 hours) accelerates adaptation by 1-2 days per treatment, while westward travel benefits from evening light to induce delays. Shift workers employ scheduled bright light during night shifts combined with dark therapy (e.g., blue-blocking glasses) during the day to realign the SCN, reducing sleep inertia and improving alertness; such interventions can shift rhythms by up to 2 hours per week with consistent application.^[10]^[10]

Melatonin and Pharmacological PRCs

Melatonin administration induces phase shifts in human circadian rhythms according to a phase response curve (PRC) characterized by phase advances during the early subjective night and phase delays in the morning or late subjective day.^[11] This PRC is approximately 12 hours out of phase with the light-induced PRC, such that melatonin advances occur when light would delay and vice versa.^[12] Seminal studies in the 1990s, building on earlier 1980s work identifying melatonin's role in circadian regulation, demonstrated these effects using physiological doses of 0.5 to 5 mg administered orally under dim light conditions to avoid suppression of endogenous melatonin.^[11]^[13] The phase-shifting mechanism of melatonin involves binding to MT1 and MT2 receptors in the suprachiasmatic nucleus (SCN), the master circadian pacemaker.^[14] Activation of MT1 receptors acutely suppresses SCN neuronal firing, while MT2 receptors mediate the phase shifts by altering clock gene expression and neuronal synchronization.^[15] This contrasts with other pharmacological agents; for instance, benzodiazepines like triazolam produce phase shifts resembling non-photic stimuli, with advances during the subjective day and delays at night, likely via GABA_A receptor modulation in the SCN.^[16] Vasopressin antagonists, targeting V1a receptors in SCN efferent pathways, can attenuate or reverse certain phase shifts, highlighting vasopressin's role in consolidating circadian outputs but with less direct phase-shifting potency than melatonin.^[17] Optimal dosing for phase advances is 3-5 mg of fast-release melatonin taken at dusk, approximately 2-4 hours before the dim light melatonin onset (DLMO), to maximize advances of up to 1-2 hours while minimizing soporific effects during the day.^[18]^[19] This timing and dose are applied in treating circadian rhythm sleep-wake disorders (CRSWDs) like delayed sleep-wake phase disorder, with FDA-approved extended-release formulations (e.g., 2 mg) available for specific indications such as insomnia in older adults.^[20] Other pharmacological agents mimicking melatonin's effects include dual MT1/MT2 agonists like tasimelteon, which exhibit similar PRCs with phase advances in the biological evening, as shown in clinical trials from the 2000s onward.^[21] Historical trials of melatonin agonists in the 1990s laid the groundwork, leading to tasimelteon's FDA approval in 2014 for non-24-hour sleep-wake disorder in blind individuals, where it entrains rhythms with daily 20 mg doses taken 1 hour before bedtime.^[22] Serotonin-modulating agents with melatonin-like activity, such as agomelatine (a MT1/MT2 agonist and 5-HT2C antagonist), also produce comparable phase advances but are primarily used for depression with circadian components.^[23] Melatonin and its agonists generally have minimal side effects at therapeutic doses, including transient drowsiness or headache, and are considered safe for short-term use (1-2 months) in CRSWD therapy.^[20] However, interactions with light exposure can reduce efficacy by suppressing endogenous melatonin, and long-term use requires monitoring due to potential disruptions in natural circadian regulation or emerging concerns about cardiovascular risks.^[24] Additive phase-shifting effects with timed light exposure have been observed but require coordinated protocols.^[25]

Non-Photic Stimuli

Non-photic stimuli, such as physical activity, feeding schedules, social interactions, and temperature changes, serve as secondary zeitgebers that can induce phase shifts in circadian rhythms, though typically weaker than those elicited by light. These cues contribute to entrainment by influencing the suprachiasmatic nucleus (SCN) indirectly or by resynchronizing peripheral clocks, complementing photic signals in maintaining circadian alignment, particularly in scenarios where light exposure is limited.^[26] Exercise represents a prominent non-photic stimulus capable of generating phase response curves (PRCs) in humans, with phase advances observed from afternoon and early evening sessions, such as 1-hour treadmill workouts at moderate intensity. A comprehensive human exercise PRC demonstrated maximum phase advances of approximately 90 minutes when exercise occurred between 1:00 pm and 4:00 pm relative to habitual wake time, while evening exercise from 7:00 pm to 10:00 pm produced delays up to 70 minutes. Studies from the 1990s, including those involving nocturnal exercise, reported phase delays of about 1 hour in both young and older adults, with effects persisting across age groups.^[27]^[28]^[29] The mechanisms underlying exercise-induced phase shifts involve arousal-mediated sympathetic nervous system activation, which indirectly modulates SCN activity through geniculo-hypothalamic tract signaling or changes in core body temperature, rather than direct retinal input. In animal models, such as Syrian hamsters, voluntary wheel-running during the subjective day elicits robust phase advances of 1-2 hours, mimicking human responses and highlighting conserved pathways across species.^[30] Beyond exercise, scheduled feeding acts as a potent non-photic cue, particularly for peripheral oscillators like the liver clock, where time-restricted feeding can phase shift hepatic gene expression independently of the SCN-driven central clock. For instance, restricting food intake to a 4-6 hour window during the rest phase advances liver clock phases by up to 4-6 hours in rodents, decoupling peripheral rhythms from light cues. Social cues, such as interactions with conspecifics exhibiting strong circadian patterns, can enhance rhythmicity or induce small phase shifts (20-40 minutes) in rodents like deer mice, while temperature pulses (e.g., warming by 3-4°C during the subjective day) produce phase delays in Syrian hamsters via thermosensitive neurons in the SCN.^[31] Despite their efficacy, non-photic stimuli generally produce smaller phase shifts (30-60 minutes) compared to light, with high variability influenced by exercise intensity, duration, and individual chronotype—evening types exhibit larger advances from morning exercise than morning types. These limitations underscore their role as adjuncts rather than primary entrainers. In practical applications, combining non-photic stimuli like exercise with light exposure enhances overall entrainment, aiding adaptation for shift workers—where timed evening exercise promotes phase delays aligning rhythms to night schedules—and athletes, who benefit from optimized timing to mitigate jet lag or irregular training. Such integrated approaches have shown additive effects in laboratory simulations of shift work, improving sleep quality and performance.^[29]^[32]

Interactions and Modeling

In circadian systems, multiple stimuli can interact to produce net phase shifts that exceed those of individual cues, often through additive effects where phase response curves (PRCs) superimpose. For instance, combining morning bright light exposure with early evening exogenous melatonin administration in humans results in larger phase advances of the dim light melatonin onset (DLMO) compared to either stimulus alone, with empirical evidence from controlled trials showing advances of approximately 68 minutes for the combination versus 37-42 minutes for single treatments. This superposition has been demonstrated in 2000s human studies, such as those using 0.5 mg evening melatonin paired with bright light to yield additive shifts of up to 1.5 hours in circadian phase.^[33]^[34]^[35] Modeling these interactions typically involves representing phase shifts as vectors in polar coordinates, where the circadian cycle is treated as a 24-hour circle, allowing the net effect of multiple zeitgebers to be computed via vector addition to predict overall entrainment. This approach assumes linearity, enabling the summation of individual PRC contributions to forecast combined outcomes like enhanced advances from light and melatonin. However, limitations arise from the nonlinear dynamics inherent in biological oscillators, where strong stimuli or interactions may saturate responses or produce non-additive effects, deviating from simple superposition and requiring more complex models to capture saturation or interactions.^[2]^[36] The circadian system exhibits hierarchical entrainment, with the suprachiasmatic nucleus (SCN) serving as the core clock that synchronizes peripheral oscillators in tissues like the liver through humoral and neural signals. While the SCN rapidly adjusts to primary cues such as light-dark cycles, peripheral clocks respond more slowly to secondary zeitgebers like feeding, leading to transient desynchrony during shifts; conflicting cues, such as nighttime meals opposing SCN-driven rest signals, exacerbate this misalignment and are linked to metabolic disorders like obesity.^[37] Computational tools facilitate modeling of these multi-stimulus interactions by incorporating multiple PRCs into simulations for predicting entrainment and guiding interventions. For example, the Circadian Performance Simulation Software (CPSS), developed for aerospace applications, integrates PRCs for light, melatonin, and sleep to simulate personalized chronotherapy outcomes, optimizing timing to minimize desynchrony. Similarly, open-source packages like the circadian library enable analysis of combined PRC effects in rhythmic data, supporting tailored strategies for rhythm adjustment.^[38]^[39] Case studies in jet lag protocols highlight practical applications, integrating light, melatonin, and sleep timing to accelerate adaptation across time zones. Pre-flight strategies, such as gradually advancing sleep schedules with morning bright light (e.g., 3000 lux for 1-2 hours) and 3-5 mg evening melatonin, have reduced eastward jet lag symptoms by 1-2 days in field trials. NASA astronaut research during space shuttle missions demonstrates similar outcomes, where combined light exposure and melatonin dosing shifted circadian rhythms by up to 9 hours to match orbital schedules, though persistent sleep fragmentation underscored the need for integrated sleep hygiene to fully mitigate desynchrony.^[40]^[41]

Neuronal Systems

PRC in Single Neurons

In single neurons, phase response curves (PRCs) characterize how brief perturbations, such as synaptic inputs or current injections, alter the timing of action potentials or bursts in oscillating excitable cells. Depolarizing inputs typically cause phase advances by accelerating the neuron toward the next firing event, while hyperpolarizing inputs induce phase delays by slowing the approach to threshold.^[42] These responses are often described by type 1 PRCs, which exhibit only advances for depolarizing stimuli and only delays for hyperpolarizing ones, a pattern prevalent in neurons involved in central pattern generators (CPGs).^[42] Experimental determination of neuronal PRCs involves isolating the oscillator and applying controlled perturbations at various phases of the cycle. Early measurements drew from general oscillator theory but adapted for neural systems using intracellular recordings in reduced preparations, such as brain slices or isolated ganglia.^[43] Common techniques include the direct pulse method, where brief current pulses (e.g., 3 ms duration, 40–110 pA) are injected at intervals corresponding to different phases during periodic firing (7–12 Hz), and the phase shift is quantified as the difference in interspike intervals.^[42] For more naturalistic inputs, dynamic clamp simulates synaptic conductances (e.g., 10–500 nS) in real-time on voltage-clamped neurons, as demonstrated in slice preparations from rodent hippocampus or isolated crustacean ganglia.^[44] Continuous low-amplitude noise injections, analyzed via spike-triggered averaging, provide estimates of infinitesimal PRCs, minimizing bias from finite perturbation strength.^[43] The shape of a neuron's PRC arises from underlying ionic conductances that govern excitability across the oscillation cycle. In advance zones, where perturbations near the firing phase elicit strong shifts, activation of voltage-gated Na⁺ channels facilitates rapid depolarization, amplifying the effect of excitatory inputs.^[45] Conversely, delay zones, often during the interspike interval, involve K⁺ conductances (e.g., delayed rectifier or Ca²⁺-activated types) that promote repolarization and hyperpolarization, enhancing inhibitory perturbations.^[45] Adaptations of the Hodgkin-Huxley model incorporate these dynamics to simulate PRCs, revealing how conductance ratios (e.g., Na⁺ vs. K⁺) determine the transition from type 1 to type 2 profiles under varying current injections.^[44] A representative example is found in pacemaker neurons of the crustacean stomatogastric ganglion, such as the pyloric dilator (PD) cells, which generate rhythmic bursts. Brief inhibitory pulses advance the burst phase when delivered early in the cycle but delay it when applied late, while excitatory pulses show the opposite pattern, demonstrating phase-dependent resets critical for rhythm stability.^[44] These PRCs, measured via dynamic clamp in isolated preparations, highlight how synaptic inputs of varying strength and duration (e.g., 200–1000 ms) can fully reset the oscillator, with shifts scaling nonlinearly for stronger stimuli.^[44] Compared to circadian PRCs, neuronal versions operate on millisecond timescales driven by ion channel kinetics, rather than hours-long cycles mediated by transcriptional feedback. Neuronal oscillators also exhibit heightened sensitivity to weak perturbations (e.g., subthreshold currents <10 pA), enabling precise phase control in real-time neural computation, unlike the robustness of circadian systems to minor noise.^[42]

PRC in Neural Networks

In neural networks, phase response curves (PRCs) from individual neurons serve as building blocks to predict collective dynamics, where synaptic interactions propagate phase shifts across the population to generate synchronized rhythms or desynchronized states. Synaptic coupling modulates these shifts depending on the timing and type of input: excitatory synapses typically induce phase advances by depolarizing the postsynaptic neuron during its rising phase, accelerating the approach to the spike threshold, while inhibitory synapses cause phase delays by hyperpolarizing the neuron, postponing the spike.^[46]^[47] These effects are particularly pronounced in type I PRCs common to many cortical neurons, where the PRC is positive for excitatory perturbations early in the cycle and negative for inhibitory ones later.^[48] To model these interactions in large populations, phase reduction techniques derive effective coupling functions from individual PRCs convolved with synaptic kinetics, extending the classic Kuramoto model to account for pulse-coupled neural oscillators and heterogeneous PRC shapes. In such extensions, the interaction term incorporates the PRC to capture how spike timing influences phase evolution, enabling analysis of synchronization transitions in excitatory-inhibitory (E-I) balanced networks.^[49]^[50] For instance, weak coupling approximations yield a mean-field description where stable synchrony emerges when the first Fourier mode of the PRC aligns with the coupling delay.^[51] Synchronization patterns in these networks often manifest as in-phase locking, as seen in theta (4-8 Hz) rhythms of the hippocampus, where mutual excitatory coupling via AMPA synapses drives coherent phase advances to sustain population bursts. Similarly, gamma oscillations (30-100 Hz) in cortical circuits arise from E-I loops, with inhibitory PRCs enforcing precise timing for clustered firing.^[52]^[53] Conversely, desynchronization can occur with mismatched PRCs across neurons, such as in heterogeneous populations where variability in intrinsic frequencies or PRC shapes leads to phase diffusion and anti-phase clustering, disrupting coherent activity.^[54]^[52] Pathological implications of PRC-mediated synchronization include epilepsy, where abnormal resets from hypersynchronous excitatory inputs amplify seizure-like activity through positive feedback in PRC shapes, as demonstrated in models of cortical networks.^[55] In Parkinson's disease, excessive synchronization in basal ganglia circuits contributes to tremor, with properties of inhibitory PRCs in the globus pallidus influencing clustered oscillations and synchronization stability.^[55]^[56] Computational modeling leverages PRC-based phase reduction for stability analysis in large networks, reducing multidimensional dynamics to a low-dimensional phase equation whose Jacobian eigenvalues determine the robustness of synchronous states against noise or heterogeneity. This approach reveals that splay states or cluster synchrony are stable when the PRC's odd Fourier components dominate, providing insights into rhythm robustness in brain-scale simulations.^[51]^[57]^[58] Experimental evidence from cortical slices confirms PRC-mediated entrainment, where periodic optogenetic or electrical stimuli applied to ferret visual cortex or rat entorhinal slices induce phase locking to external rhythms, with measured PRCs predicting the strength of synchronization based on input timing relative to endogenous oscillations.^[3]^[59]

Historical Context

Early Conceptualization

The concept of synchronization in oscillating systems predates biological applications, tracing back to Christiaan Huygens' 17th-century observation that two pendulum clocks suspended from the same beam would gradually synchronize their swings, either in phase or antiphase, due to mechanical coupling through the supporting structure.^[60] This non-biological precursor laid a foundational idea for understanding how perturbations could adjust the timing of rhythms, though it was not initially linked to living systems. The first explicit application to biological rhythms emerged in the 1940s with studies on tidal influences, where Frank A. Brown Jr. examined how environmental cycles like tides entrained activity patterns in intertidal organisms such as fiddler crabs, suggesting phase adjustments to periodic stimuli without formalizing the response as a curve.^[61] In the mid-1950s, Colin S. Pittendrigh advanced these ideas through experimental work on the fruit fly Drosophila pseudoobscura, focusing on the eclosion rhythm—the timed emergence of adults from pupae. Pittendrigh and his collaborator Virginia G. Bruce demonstrated that brief light pulses applied at different phases of the circadian cycle induced phase-dependent resets, with delays occurring early in the subjective night and advances later, revealing the rhythm's sensitivity to zeitgebers like light.^[62] This work established the phase-dependent nature of entrainment in biological clocks, using population-level eclosion data to infer an underlying oscillator model where stimuli shifted the rhythm's phase rather than merely suppressing or advancing it uniformly.^[63] Building on Pittendrigh's empirical foundation, Arthur T. Winfree provided a mathematical formalization of phase resetting in the 1960s, introducing the "phase reset" concept to describe how perturbations alter an oscillator's trajectory on a limit cycle. In his seminal 1967 paper, Winfree modeled populations of coupled oscillators, showing that phase shifts depend on the timing and strength of the stimulus, and summarized these ideas in his 1980 book The Geometry of Biological Time, which synthesized decades of rhythm research into geometric interpretations of phase space. Winfree's framework emphasized the isochronal structure of oscillators, where points near a singularity could lead to complete resets, providing a theoretical basis for predicting entrainment in circadian systems.^[64] Early applications extended to neuronal systems in the 1960s, with David H. Perkel and colleagues studying phase shifts in identified neurons of the sea slug Aplysia californica. By injecting current pulses into bursting neurons like R15, they observed that perturbations at specific phases induced measurable advances or delays in the ongoing oscillation, demonstrating phase-dependent responses analogous to those in circadian rhythms.^[65] These intracellular experiments highlighted how synaptic or electrical inputs could reset neuronal firing patterns, bridging single-cell electrophysiology with broader oscillator theory. Early conceptualizations of phase responses were limited by their focus on binary outcomes—simple advances or delays—rather than the full curvature of a response function, often overlooking the continuous, stimulus-strength-dependent nature of shifts and the potential for type 0 resetting near critical phases.^[66] This simplified view sufficed for initial entrainment models but constrained deeper analyses of complex interactions until later refinements.

Key Developments and Applications

The concept of phase response curves (PRCs) emerged in the mid-20th century as a tool to quantify how external stimuli perturb endogenous oscillators, with Jürgen Aschoff's 1965 review introducing response curves to describe phase-dependent shifts in circadian rhythms of birds and mammals exposed to light pulses.^[67] Colin Pittendrigh advanced this framework in the 1960s and 1970s by demonstrating PRCs in Drosophila and rodents, showing how light elicits delays early in the subjective night and advances late, which informed early models of entrainment to environmental cycles.^[68] A pivotal development occurred in 1976 when Pittendrigh and Serge Daan analyzed PRC variability across nocturnal rodents, revealing that PRC shape correlates with ecological adaptations, such as period lability, and proposed dual pacemaker models (morning and evening) to explain entrainment stability.^[69] Arthur Winfree's 1980 monograph, The Geometry of Biological Time, formalized PRC theory mathematically, distinguishing Type 1 PRCs (weak stimuli causing proportional small shifts) from Type 0 PRCs (strong stimuli enabling full-cycle resets), providing a geometric basis for understanding singularity points where oscillators lose amplitude. This work influenced subsequent experimental designs, as seen in compilations like Carl H. Johnson's 1990 PRC Atlas, which cataloged over 200 PRCs from diverse organisms and stimuli, facilitating comparative studies on oscillator robustness.^[70] In neuronal systems, PRCs gained traction in the 1990s for modeling synchronization, with G. Bard Ermentrout and Nancy Kopell's analyses showing how PRC shapes (e.g., Type II with advances and delays) predict network synchrony in weakly coupled oscillators, applied to central pattern generators in locomotion.^[51] By the 2000s, PRCs extended to human applications, such as Shawn D. Youngstedt's 2019 study establishing exercise-induced PRCs with delays up to 1 hour when timed in the evening, informing non-pharmacological interventions.^[26] Key applications of PRCs include chronotherapy, where timing light or melatonin pulses based on individual PRCs treats circadian sleep disorders; for instance, morning light advances phases by up to 2 hours in delayed sleep phase syndrome, improving alignment with social schedules.^[71] In oncology, PRC-guided drug scheduling enhances efficacy, as circadian variations in cell proliferation amplify chemotherapy impacts when aligned with tumor PRCs, reducing toxicity in clinical trials.^[72] For jet lag and shift work, PRC models predict optimal exposure times, recommending light strategies to shift phases by about 1 hour per day.^[73] In neuroscience, PRCs underpin deep brain stimulation for Parkinson's disease, where phase-specific pulses desynchronize aberrant beta oscillations, improving motor symptoms; simulations using PRCs aid in predicting synchronization in heterogeneous networks.^[51] By the 2020s, optogenetics and computational models have refined PRC applications in neuroscience and medicine.^[74] These developments highlight PRCs' role in bridging single-oscillator dynamics to systemic therapies, with ongoing refinements via optogenetics revealing stimulus-intensity effects on PRC steepness.^[3]

References

[1]
Phase response curve - Scholarpedia
Dec 30, 2006 · A phase response curve (PRC) tabulates the transient change in the cycle period of an oscillator induced by a perturbation as a function of the ...Definition and Methodology · Advanced PRC Methodology
[2]
Phase Response Curves - Biological Clocks, Rhythms, and ... - NCBI
Phase response curves (PRCs) measure how a stimulus affects a biological clock's phase, plotting the change in phase (ϕnew - ϕold) against the initial phase (ϕ ...
[3]
A comparison of methods to determine neuronal phase-response ...
The phase-response curve (PRC) of a regularly firing neuron quantifies the shift in the next spike time as a function of the timing of a small perturbation ...
[4]
[PDF] Type I Membranes, Phase Resetting Curves, and Synchrony
In Rinzel and Ermentrout (1989) we review the classification of ex- citable membranes by Hodgkin (1948) in terms of their dynamics as a current is injected.
[5]
[PDF] Stochastic Neuroscience
which is called the adjoint for the limit cycle; it is also called the infinitesimal. PRC. The adjoint is the unique periodic solution to: dZ dt= −DXF(X0(t)) ...
[6]
[PDF] Infinitesimal Phase Response Curves and Adjoint Variational Equation
In this document, I am trying to convince the reader about a systematic approach to introduce Infinitesimal. Phase Response Curves (iPRCs) for small ...
[7]
Bright Light Induction of Strong (Type 0) Resetting of the Human ...
The response of the human circadian pacemaker to light was measured in 45 resetting trials. Each trial consisted of an initial endogenous circadian phase ...
[8]
and dark-pulse phase response curves for diurnal Octodon degus ...
Two different PRCs were generated in response to 20-min light pulses. The majority of animals produced significant phase delays between circadian time (CT) 0 ...Missing: species differences biphasic
[9]
Effect of Light on Human Circadian Physiology - PMC
These studies have revealed, how the timing, intensity, duration, and wavelength of light affect the human biological clock.Figure 1 · Figure 2 · Figure 3
[10]
Circadian Rhythm Sleep Disorders: Part I, Basic Principles, Shift ...
Although bright light exposure (3000-10,000 lux) has been shown to produce robust phase shifts, even modest intensities (50-600 lux) can produce substantial ...
[11]
Melatonin shifts human circadian rhythms according to a phase ...
A physiological dose of orally administered melatonin shifts circadian rhythms in humans according to a phase-response curve (PRC) that is nearly opposite in ...Missing: 1980s | Show results with:1980s
[12]
Phase shifting the human circadian clock using melatonin
Recently, we have described a phase response curve (PRC) to melatonin in humans that is about 12 h out of phase with the PRC to light. Exogenous melatonin ...
[13]
[PDF] Melatonin Shifts Human Circadian Rhythms According to a Phase ...
Summary: A physiological dose of orally administered melatonin shifts circadian rhythms in humans according to a phase-response curve (PRC) that is nearly.
[14]
MT1 and MT2 Melatonin Receptors: A Therapeutic Perspective - PMC
Endogenous melatonin released from the pineal gland at night may feedback onto the SCN and activate MT1 and MT2 receptors to phase shift local and overt ...
[15]
Melatonin receptors: Role on sleep and circadian rhythm regulation
Two G-protein coupled melatonin receptors, the MT1 and MT2, inhibit neuronal activity and phase shift circadian firing rhythms in the SCN, respectively. Recent ...
[16]
Use of benzodiazepines to manipulate the circadian clock regulating ...
Triazolam has also been found to facilitate the rate of reentrainment of the activity rhythm following an 8-hour advance or delay in the light-dark cycle.
[17]
[PDF] Vasopressin Resets the Central Circadian Clock in a Manner ...
Conclusions: Collectively, this work indicates that AVP signaling resets SCN molecular rhythms in conjunction with VIP signaling and in a manner influenced ...
[18]
[PDF] Clinical Practice Guideline for the Treatment of Intrinsic Circadian ...
A fast-release formulation of melatonin was. 251 utilized, with dosages ranging from 3-5 mg, taken between 18:00-19:00 (no circadian-. 252 based timing), for 4 ...
[19]
Human Phase Response Curves to Three Days of Daily Melatonin
Maximum advances occurred when 0.5 mg melatonin was taken in the afternoon, 2–4 h before the DLMO, or 9–11 h before sleep midpoint.
[20]
Melatonin - Mayo Clinic
However, evidence suggests that melatonin supplements promote sleep and are safe for short-term use. Melatonin can be used to treat delayed sleep phase and ...
[21]
Comparative Review of Approved Melatonin Agonists for the ...
Currently, only one melatonin receptor agonist, tasimelteon, is approved for the treatment of a CRSWD: non–24‐hour sleep‐wake disorder (or non‐24).
[22]
Tasimelteon for non-24-hour sleep–wake disorder in totally blind ...
Aug 4, 2015 · We assessed safety and efficacy (in terms of circadian entrainment and maintenance) of once-daily tasimelteon, a novel dual-melatonin receptor agonist.
[23]
Comparative Review of Approved Melatonin Agonists for the ...
Aug 8, 2016 · The human phase response curve (PRC) to melatonin is about 12 hours out of phase with the PRC to light. Chronobiol Int 1998; 15: 71–83 ...
[24]
Current Insights into the Risks of Using Melatonin as a Treatment for ...
Jan 12, 2023 · Melatonin appears to have a favorable safety profile in this population, however there is a dearth of evidence regarding the safety of prolonged use.
[25]
Clinical Implications of the Melatonin Phase Response Curve - PMC
In this issue of JCEM, the Burgess et al. report (1) on the melatonin phase response curve (PRC) has many clinical implications worthy of further discussion ...
[26]
Human circadian phase–response curves for exercise - Youngstedt
Feb 19, 2019 · The phase–response curve describing the magnitude and direction of circadian rhythm phase shifts, depending on the time of the zeigeber (time ...Missing: biological | Show results with:biological
[27]
Human circadian phase–response curves for exercise - PMC
Exercise causes circadian phase shifts. Phase advances occur at 7:00 am and 1:00 pm to 4:00 pm, and delays from 7:00 pm to 10:00 pm.
[28]
Circadian phase-shifting effects of nocturnal exercise in older ...
Feb 6, 2003 · Circadian phase-shifting effects of nocturnal exercise in older compared with young adults. Erin K. Baehr,1,2 Charmane I. Eastman,3 William ...
[29]
Intermittent bright light and exercise to entrain human circadian ...
Bright light can phase shift human circadian rhythms, and recent studies have suggested that exercise can also produce phase shifts in humans.
[30]
Age-related circadian disorganization caused by sympathetic ...
Jan 5, 2017 · Activation of the sympathetic nervous system by stress/exercise is an important signaling pathway for peripheral clock phase adaptation. To ...
[31]
Social influences on mammalian circadian rhythms: animal and ...
Mar 15, 2007 · ... circadian responsiveness to social cues in the diurnal rodent ... social stimuli do not phase shift the circadian temperature rhythm in rats.
[32]
Move the night way: how can physical activity facilitate adaptation to ...
Mar 2, 2024 · Strategically timed physical activity could help shift workers adapt their circadian rhythms to night shift work by eliciting a phase delay.
[33]
Combination of Light and Melatonin Time Cues for Phase ... - NIH
Morning bright light combined with early evening exogenous melatonin induced a greater phase advance of the DLMO than either treatment alone.
[34]
https://doi.org/10.1111/j.1600-079X.2004.00117.x
[35]
https://doi.org/10.1210/jc.2005-1009
[36]
(PDF) Chapter 1 Phase Response Curves - ResearchGate
Oct 21, 2014 · Phase response curves (PRCs) are widely used in circadian clocks, neuroscience, and heart physiology. They quantify the response of an oscillator to pulse-like ...
[37]
A minimal model of peripheral clocks reveals differential circadian re ...
The mammalian circadian system consists of a hierarchical network of oscillators, where the central clock in the suprachiasmatic nucleus (SCN) is thought to ...
[38]
Mathematical Modeling of Circadian/Performance Countermeasures
The mathematical modeling and the available Circadian Performance Simulation Software (CPSS) can be used to simulate and quantitatively evaluate different ...Missing: multiple | Show results with:multiple
[39]
circadian by Arcascope - GitHub Pages
A computational package for the simulation and analysis of circadian rhythms. Install circadian can be installed via pip.
[40]
Advancing Circadian Rhythms Before Eastward Flight: A Strategy to ...
A gradually advancing sleep schedule with intermittent morning bright light can be used to advance circadian rhythms before eastward flight.
[41]
[PDF] Sleep, Circadian Rhythms, and Performance During Space Shuttle ...
Ground-based research has shown that administration of melatonin when sleep is attempted at cir- cadian phases at which the body's own melatonin is absent.
[42]
Hippocampal CA1 pyramidal neurons exhibit type 1 phase ... - NIH
We introduce a novel, bias-correction method to estimate the infinitesimal phase-resetting curve (iPRC) and confirm type 1 excitability in hippocampal pyramidal ...Missing: seminal | Show results with:seminal
[43]
A Comparison of Methods to Determine Neuronal Phase-Response ...
The phase-response curve (PRC) is an important tool to determine the excitability type of single neurons which reveals consequences for their synchronizing ...Missing: seminal | Show results with:seminal
[44]
The Functional Consequences of Changes in the Strength and ...
The phase responses of the biological neuron show the same behavior as the simulated PRCs of the model bursting pacemaker: for progressively longer stimulus ...
[45]
https://doi.org/10.1007/978-1-4614-0965-9_10
[46]
Effect of Phase Response Curve Shape and Synaptic Driving Force ...
Jul 1, 2017 · Abstract. The role of the phase response curve (PRC) shape on the synchrony of synaptically coupled oscillating neurons is examined.<|control11|><|separator|>
[47]
Effect of phase response curve shape and synaptic driving force on ...
The role of the phase response curve (PRC) shape on the synchrony of synaptically coupled oscillating neurons is examined.
[48]
Phase response curve for a type I excitatory and inhibitory neuron...
Phase response curve for a type I excitatory and inhibitory neuron representing the phase change in the spiking response of the neuron to a depolarizing pulse, ...<|separator|>
[49]
Kuramoto Model for Excitation-Inhibition-Based Oscillations
Jun 13, 2018 · The Kuramoto model (KM) is a theoretical paradigm for investigating the emergence of rhythmic activity in large populations of oscillators.Missing: extensions PRCs
[50]
Kuramoto model for populations of quadratic integrate-and-fire ...
Jan 3, 2022 · In this paper, we aim to theoretically substantiate the use of the KM for neuronal modeling by providing a mathematical link between the ...Missing: PRCs | Show results with:PRCs
[51]
Phase-response curves and synchronized neural networks - PMC
We review the principal assumptions underlying the application of phase-response curves (PRCs) to synchronization in neuronal networks.
[52]
Phase-Resetting Curves Determine Synchronization, Phase Locking ...
Apr 22, 2009 · The analytical method directly calculates the phasic relationships and stability of oscillatory modes that are exhibited by the network from the ...
[53]
Emergence of global synchronization in directed excitatory networks ...
Feb 24, 2020 · The phase-response curve (PRC) is an experimentally obtainable measure of cellular properties that quantifies the shift in the next spike time ...
[54]
Synchronization of Excitatory Neurons with Strongly Heterogeneous ...
Nov 26, 2007 · However, the PRCs recorded from various brain areas, which include the hippocampus [12] , the entorhinal cortex [13] , the somatosensory cortexMissing: mismatched | Show results with:mismatched
[55]
Desynchronization boost by non-uniform coordinated reset ... - NIH
May 17, 2013 · Desynchronization of abnormal neural synchrony is theoretically compelling because of the complex dynamical mechanisms involved. We here present ...
[56]
Response of Neuronal Populations to Phase-Locked Stimulation
Apr 9, 2025 · This study validates a mathematical model of coupled oscillators in predicting the response of neural activity to stimulation for the first time.
[57]
Deep brain stimulation-guided optogenetic rescue of parkinsonian ...
May 13, 2020 · Both DBS and SST opto-activation increased the correlation coefficients when compared to parkinsonian condition, with smoother response profiles ...
[58]
Predicting the Existence and Stability of Phase-Locked Mode in ...
Aug 1, 2017 · Predicting phase-locked modes in large neural networks usually requires as a first step a complexity reduction to manageable subnetworks of ...
[59]
A framework for macroscopic phase-resetting curves for generalised ...
Aug 1, 2022 · Inferring the PRC and developing a systematic phase reduction theory for large-scale brain rhythms remains an outstanding challenge. Here we ...
[60]
Optogenetic stimulation effectively enhances intrinsically generated ...
Oct 22, 2013 · Our study shows that mild stimulation protocols that do not enforce particular activity patterns onto the network can be highly effective ...
[61]
Researchers prove Huygens was right about pendulum ... - Phys.org
Mar 29, 2016 · In 1665 Christiaan Huygens discovered that two pendulum clocks, hung from the same wooden structure, will always oscillate in synchronicity.
[62]
COMPENSATION FOR TEMPERATURE IN THE METABOLISM AND ...
Dehnel and Earl Segal, who themselves have contributed the original work from this laboratory. REFERENCES.<|separator|>
[63]
https://www.pnas.org/doi/pdf/10.1073/pnas.58.5.1862
[64]
Circadian systems, II. The oscillation in the individual Drosophila pupa
Solid circles indicate observed eclosion peaks (medians); populations illuminated at ct 20 exploit mainly the advanced eighth gate, but a few flies (small ...
[65]
Jet Wimp*
(Adapted from The Geometry of Biological Time (GBT).) see Figure 1. More generally, the term phase singularity re- fers to a locus of points where the phase is ...
[66]
Biological oscillators can be stopped—Topological study of a phase ...
Phase response curves, which are measured in these experiments, can be classified into two types. The one is the curve with one mapping degree (Type 1), and the ...
[67]
Phase Response Curves
A discontinuous fluence response curve means that some process "downstream" from the photopigment's absorption of light is converting the initially continuous ...
[68]
Can phase response curves of various treatments of circadian ...
Aschoff, 1965. J. Aschoff. Response curves in circadian periodicity. J. Aschoff (Ed.), Circadian Clocks, North Holland Publ. Co, Amsterdam (1965), pp. 95-111.
[69]
A Functional analysis of circadian pacemakers in nocturnal rodents
Phase response curves for 15′ bright light pulses of four species of nocturnal rodents are described. All show delay phase shifts early in the subjective night.
[70]
A functional analysis of circadian pacemakers in nocturnal rodents
Daan, S., Pittendrigh, C. S.: A functional analysis of circadian pacemakers in nocturnal rodents. II. The variability of phase response curves. J. comp ...
[71]
PRC Atlas - Carl Johnson Laboratory | Vanderbilt University
A Phase Response Curve (PRC) is a plot of phase-shifts as a function of circadian phase of a stimulus such as light pulses, temperature pulses, or pulses of ...
[72]
Phase Response Curve (PRC) - Circadian Sleep Disorders Network
Oct 14, 2021 · The Phase Response Curve (PRC) is a curve describing the relationship between a stimulus, like light, and a shift in circadian rhythm.
[73]
An Overview of the Circadian Clock in the Frame of Chronotherapy
This review discusses the links between an altered circadian clock and the rise of pathologies. We then sum up the proofs of concept advocating for the ...
[74]
Circadian Rhythms, Disease and Chronotherapy - PMC
Exploiting such rhythm synchronizers can enhance body clock functions, thereby improving circadian rhythms and perhaps reducing susceptibility to diseases that ...