A dynamic clamp protocol to artificially modify cell capacitance
Abstract
Dynamics of excitable cells and networks depend on the membrane time constant, set by membrane resistance and capacitance. Whereas pharmacological and genetic manipulations of ionic conductances of excitable membranes are routine in electrophysiology, experimental control over capacitance remains a challenge. Here, we present capacitance clamp, an approach that allows electrophysiologists to mimic a modified capacitance in biological neurons via an unconventional application of the dynamic clamp technique. We first demonstrate the feasibility to quantitatively modulate capacitance in a mathematical neuron model and then confirm the functionality of capacitance clamp in in vitro experiments in granule cells of rodent dentate gyrus with up to threefold virtual capacitance changes. Clamping of capacitance thus constitutes a novel technique to probe and decipher mechanisms of neuronal signaling in ways that were so far inaccessible to experimental electrophysiology.
Editor's evaluation
The manuscript introduces a new enhancement to the dynamic clamp technique, CapClamp that, analogous to the artificial conductances of standard Dynamic Clamp, allows the experimenter to adjust the somatic time constant by setting a new membrane artificial capacitance independent of any change in input resistance. The technique is shown to have application for studying temporal integration, energetic costs of spiking and bifurcations. The technique is rigorously tested in model and physiological application and is robust when sampling frequency of the feedback (clamp) loop is fast compared to the fastest electrical event in a neuron (usually action potentials), and for vertebrate neurons it should be 20KHz or faster and yet faster for fast spiking neurons.
https://doi.org/10.7554/eLife.75517.sa0Introduction
Membrane capacitance is a major biophysical parameter in neurons and other excitable cells, which determines how fast the membrane potential changes in response to a current (Golowasch et al., 2009; White and Hooper, 2013). How capacitance impacts electrical signaling and neuronal processing, however, can rarely be observed experimentally, because besides reduced values in myelinated axons (Hartline and Colman, 2007) most membranes appear to have a specific membrane capacitance in the range of 0.7–1.0 μF/cm2 (Gentet et al., 2000). The effects of capacitance changes can, therefore, so far only be compared via mathematical simulations, where capacitance is simple to control. Such modeling, for example, suggests that the reduced membrane capacitance observed in human pyramidal cells can serve to increase synaptic efficacy or propagation speed of action potentials (Eyal et al., 2016, but see BeaulieuLaroche et al., 2018). In contrast, experimental manipulation of capacitance remains challenging; in particular because changes in membrane area, thickness and lipid composition that affect capacitance might influence other membrane functions, such as the embedding of ion channels, with potentially unintended and uncontrolled consequences for electrical behavior. Here, we address this technical challenge by introducing capacitance clamp (CapClamp): an intracellular recording mode based on the dynamic clamp that emulates altered capacitance values in biological neurons (Robinson, 1994; Sharp et al., 1993). Via CapClamp, the voltage dynamics governed by the actual biophysics of a cell – active ion channels and synaptic inputs – can thus be flexibly probed under multiple ‘virtual’ capacitance conditions, which provides precise experimental control over this hitherto inaccessible parameter.
In addition to the analysis of biological capacitance adaptations, control over capacitance offers a distinct way to probe cellular electrical dynamics. Capacitance has a unique temporal role, because its direct effects are restricted to the membrane time constant whilst leaving the steady state IV function unaltered. In this way, capacitance differs from leak conductance, the other determinant of the time constant, as the latter also alters steadystate response amplitudes. For this reason, theoretical studies preferentially vary capacitance to investigate ion channel dynamics (Jaffe and Brenner, 2018; Franci et al., 2018) and qualitative switches (bifurcations) in neural excitability (Kirst et al., 2015; Hesse et al., 2017). Furthermore, effects of an altered capacitance can be informative about more complex, time scalerelated parameters like temperature or ion concentrations (Contreras et al., 2020). Such computational predictions, however, often rely on simplified neuron models, so a similar experimental control over capacitance would be desirable to test them in biological cells.
The proposed CapClamp alters capacitance in a virtual manner, combining the simplicity of computational control with the complex biophysics of a real neuron. It is inspired by the dynamic clamp technique, which has originally been developed to simulate the presence of additional conductances in a biological neuron relying on a fast feedback loop between intracellular recording and a computational model (Robinson, 1994; Sharp et al., 1993; Prinz et al., 2004; Economo et al., 2010). The precise control over these virtual conductances enables electrophysiological experiments that are more difficult or even impossible with traditional pharmacological or genetic means (Turrigiano et al., 1996; Svirskis et al., 2004; Prescott et al., 2008b; Hasenstaub et al., 2010; Szűcs et al., 2017; Pfeiffer et al., 2020). Here, we demonstrate how the dynamic clamp can be extended to enable manipulations of the apparent membrane capacitance by currents designed to speed up or slow down dynamics of the membrane potential. We derive a simple expression for these CapClamp currents, which can be applied in all excitable cells and only requires the experimenter to specify the original cell and the desired target capacitance. In an experiment based on a hardwareimplemented RC circuit, we verify that the CapClamp indeed correctly modifies the time constant. Via numerical simulations, we confirm that a clamped model neuron exhibits the same pronounced changes of firing and spike shape as a control cell with an altered capacitance. For an experimental demonstration, we clamp the nearsomatic capacitance of rat dentate gyrus granule cells and analyze how the induced local capacitance change affects their spiking behavior. Finally, we illustrate how the CapClamp can be used to probe signal integration and energy consumption of excitable cells in ways that so far were experimentally inaccessible.
Results
Capacitance clamp: A dynamic clamp protocol to mimic capacitance changes
Dynamic clamp relies on a fast feedback loop between an intracellular recording of a neuron and a computer that simulates virtual cellular or circuit components online. Originally, the dynamic clamp has been developed to study how a membrane conductance alters the neuron’s voltage dynamics (Sharp et al., 1993; Robinson, 1994). In each sampling interval (i.e. time interval between two voltage samplings), a digital model of the conductance receives the sampled membrane potential, updates the conductance state and sends the corresponding current value back to the amplifier. Given a sufficiently high update rate f_{dyn} (often ≥10 kHz), this current injected via the recording electrode makes the dynamics of the neuron appear as if the virtual channels represented by the conductance model were physically present in the membrane.
Whereas conductances gate ionic currents across the membrane, the capacitance determines how fast these currents can change the membrane potential. Every altered membrane property that results in a modified capacitance value, such as membrane area, thickness or lipid composition, affects this rate of change of the membrane potential (Figure 1A). To artificially mimic a modified capacitance, we therefore first asked whether a dynamic clamp protocol with its fast feedback loop between voltage sampling and current injection could adjust the ‘speed’ of a cell’s membrane potential (Figure 1B). Using the current balance equation, the basic mathematical description of membrane voltage dynamics, we derived a capacitance clamp (CapClamp) scheme with a simple expression for the clamping current ${I}_{\mathrm{dyn}}$ (see "Derivation of the CapClamp current" in Methods),
which only requires the experimenter to measure the cell capacitance ${C}_{c}$ in order to set a new target capacitance ${C}_{t}$. In every sampling interval $\mathrm{\Delta}t={{f}_{\mathrm{d}\mathrm{y}\mathrm{n}}}^{1}$, the CapClamp uses the measured cell capacitance value ${C}_{c}$ and the voltage derivative $\frac{{V}_{i}{V}_{i1}}{\mathrm{\Delta}t}$ to estimate the present membrane current and then increases ($C}_{t}<{C}_{c$) or decreases ($C}_{c}<{C}_{t$) the net current by insertion of a correction current in the next time bin. In this way, despite a physically unaltered capacitance, the membrane potential changes faster or, respectively, slower – as if the clamped cell actually had the different capacitance ${C}_{t}$ selected by the experimenter. In the following, we will demonstrate the CapClamp in simulated and experimental scenarios with increasing complexity ranging from a passive RC circuit up to biological neurons with a spatially extended morphology.
Clamping capacitance in a passive cell
The simplest scenario to apply the CapClamp is a single compartment passive cell, equivalent to an RC circuit. In the absence of active conductances, the effects of a capacitance change can be precisely formulated: the capacitance $C$ sets the membrane time constant $\tau =RC$, determining how fast the membrane potential changes in response to a current. Note that, in contrast to the resistance $R$, the change in capacitance leaves the voltage amplitude of the steadystate response unaltered. To quantitatively confirm the effects of clamping capacitance and the ability of an exclusively temporal control, we measured time constant and capacitance of a clamped RC circuit in experiment and analyzed the temporal filtering properties of a modeled clamped circuit using mathematical analysis.
To experimentally characterize a clamped passive cell, we implemented the CapClamp scheme in a dynamic clamp setup (see "Dynamic clamp setup" in Methods) and recorded voltage responses to current pulses from the simplest possible model cell, that is, a hardware implemented RC circuit, while clamping it at a range of target capacitances (Figure 1C). As expected for an RC circuit, the charging curve of the unclamped model cell was fit well by a single exponential, whose time constant ($\tau $ = 11.1ms) and voltage amplitude ($\mathrm{\Delta}V$ = −9.9 mV) allowed us to determine the circuit’s resistance $R$=99.4 MΩ and capacitance $C$=112.3 pF. This capacitance value was then used as the cell capacitance ${C}_{c}$ input for the CapClamp. Clamped at a decreased capacitance, the time constant shortened (${C}_{t}$ = 67.4 pF: $\tau $=6.6ms) and at an increased capacitance, it lengthened (${C}_{t}$ = 336.9 pF: $\tau $=33.0ms), but in both cases the steady state voltage amplitude remained the same. Accordingly, the measured capacitance of the clamped circuit confirmed the chosen target capacitance for the whole tested range from a 0.6 up to a 3fold change with respect to the original capacitance (e.g. ${C}_{t}$=67.4 pF: $C$=67.5 pF; ${C}_{t}$=336.9 pF: $C$=338.1 pF), whereas the measured resistance remained constant (Figure 1D).
As a consequence of the correctly transformed voltage response, the leak current in the clamped RC circuit also behaved as if the capacitance had changed. When the circuit was clamped, the leak current through the resistance, ${I}_{R}=\frac{V}{R}$, exhibited a shorter ($C}_{t}<{C}_{c$) or longer ($C}_{t}>{C}_{c$) transient until reaching steady state. Further, the charge $\mathrm{\Delta}Q({I}_{R},{I}_{\mathrm{ext}})$ deposited on the capacitance by the apparent ‘transmembrane’ current, the sum of leak and external stimulus current, reduced ($C}_{t}<{C}_{c$) or increased ($C}_{t}>{C}_{c$) to the extent expected for an altered capacitance (Figure 1C). In contrast, the overall deposited charge $\mathrm{\Delta}Q({I}_{R},{I}_{\mathrm{ext}},{I}_{\mathrm{dyn}})$, including the clamping current, attained the same steadystate amplitude in the clamped and the original circuit, reflecting that the physical capacitance did not change. For the simple RC circuit considered here, the distinction between the clamping current and the intrinsic ‘cellular’ currents might appear artificial, because all currents use the same charge carrier. In a biological neuron, however, this distinction becomes relevant, because the clamping currents through the recording electrode might rely on other charge carriers (depending on the used intracellular solution) than the cellular currents governed by multiple ion selective channel types.
For more complex stimuli than a simple current pulse, the temporal filtering properties of a clamped membrane determine how well the CapClamp mimics the chosen capacitance change. To generally assess these filtering properties, we analytically derived the frequencydependent impedance of a modeled clamped RC circuit using linear control theory (Figure 1—figure supplement 1b A, see "Impedance of a capacitanceclamped RC circuit" in Appendix 1). The derived impedance profiles confirmed the experimentally observed altered time constants. For example, an RC circuit clamped at an increased capacitance further attenuated nonzero frequencies reflecting its longer time constant. Overall, impedance amplitudes of a clamped RC and the corresponding target circuit fit well up to a tenth of the dynamic clamp frequency f_{dyn}, that is up to ≈2 kHz for a 20 kHz dynamic clamp system as used here (Figure 1—figure supplement 1B and C ). As high frequencies are heavily attenuated by the low pass filter of a cell’s membrane, these differences lead to relatively small deviations in the voltage responses. The mathematical analysis thus suggests that for a fast dynamic clamp system ( > 20 kHz), the CapClamp is expected to work well for most stimuli with time scales in the physiological range.
Simulation of the capclamp in a biophysical neuron model
In neurons with active spikegenerating conductances, capacitance changes impact neuronal firing via the interplay of the altered membrane time constant and the gating kinetics of the channels involved. As gating dynamics can be in the submillisecond range, for example for transient sodium channels, the CapClamp is expected to require a sufficiently high dynamic clamp frequency to accurately reproduce changes of spike shape or firing rate. To understand these requirements and lay the ground for investigations of capacitance changes in biological neurons, we simulated the CapClamp in a neuron model with biophysical channel dynamics and a singlecompartment morphology (see "Biophysical neuron model" in Methods). The simulation allowed us to compare the firing of the clamped neuron to the expected firing at this modified capacitance.
Specifically, we inspected the spiking responses to a depolarizing current for the original 150 pF, a decreased 90 pF and an increased 210 pF capacitance, for the latter two comparing clamped and expected dynamics (Figure 2A). Capacitance changes exerted a notable influence on both firing frequency and spike shape, which was mostly wellcaptured by the simulated CapClamp (Table 1). When the capacitance was decreased to 90 pF, spiking became faster and action potentials had a larger peak amplitude, a decreased duration and an increased afterhyperpolarization (AHP). When the capacitance was increased to 210 pF, the effects were opposite: spiking became slower and action potentials had a smaller peak amplitude, an increased duration and a reduced AHP. At decreased capacitances, spike amplitudes of the clamped neuron were larger than in the respective control simulation, a consequence of the limited tracking of the fast sodium current at the dynamic clamp frequency used (Figure 2B and C). Except for this brief overshoot, the CapClamp overall forced the membrane potential on the expected trajectory and correctly adjusted the resulting ionic currents and the gating variable dynamics of the active conductances. For example, at a reduced capacitance of 90 pF, sodium channels inactivated less during the fast rise of the AP and therefore the sodium current exhibited a second peak during AP repolarization (see the sodium inactivation variable $h$ at AP peak time in Figure 2B).
A subsequent comparison of simulated spiking for the whole range of tested target capacitances from 75 pF to 225 pF confirmed that the CapClamp reliably reproduced the main effects of a modified capacitance on spike shape (Figure 2D) and firing frequency (Figure 2E). Furthermore, the obtained frequencycurrent curves fit well with the theoretically expected reduction of excitability at higher capacitance: a decrease of gain proportional to $\frac{1}{C}$ and a constant rheobase current (see "Analytically expected effect of capacitance on the form of the fI curve" in Methods). A crucial factor for the CapClamp, especially for a good quantitative fit of the spike shape, is the dynamic clamp frequency – observable differences at a 20 kHz sampling frequency were strongly reduced for a sampling frequency of 100 kHz (Figure 2C and D). In this regard, the chosen neuron model is especially demanding because its rapid gating dynamics are fit to a fast spiking interneuron (Wang and Buzsáki, 1996). Taken together, our simulations show that capacitance impacts neuronal spiking from firing frequency to action potential shape and that the CapClamp is wellsuited to study these effects.
Experimental demonstration of the CapClamp in rat dentate gyrus granule cells
Biological neurons differ from the simple ‘cells’ considered so far, that is RC circuit and single compartment neuron model, in one major aspect: they can have complex morphologies, where the membrane potential varies between different compartments and membrane capacitance is distributed across the neuronal structure. As the CapClamp in contrast operates locally through the recording electrode, the emulated capacitance change is expected to be localized to the recorded compartment instead of affecting all compartments. To demonstrate such localized capacitance changes and study their effects on neuronal spiking, we applied the CapClamp in in vitro patchclamp recordings of rat dentate gyrus granule cells (DGGCs). Among morphologically complex cells, DGGCs appear wellsuited to test the CapClamp, because their morphological structure, consisting of a central soma and one to four primary apical dendrites as shown in Figure 3A (Rihn and Claiborne, 1990), translates to a relatively compact electrotonic structure (SchmidtHieber et al., 2007; Wybo et al., 2019).
Measurement of local nearsomatic capacitance
Most capacitance measurements aim to provide an accurate estimate of the global capacitance of a neuron (Golowasch et al., 2009; White and Hooper, 2013). To correctly infer the transmembrane and axial current, however, the CapClamp requires the local capacitance value of the compartment where the electrode is placed at. For the somatic DGGC recordings, we exploit that the current clamp step method – fitting charging curves via a sum of exponential terms – can also provide local capacitance information (Golowasch et al., 2009). DGGC charging curves consisted of a slow ($\tau}_{0$: 15.1 ± 4.8 ms, R_{0}: 127 ± 45 MΩ) and a fast ($\tau}_{1$: 0.77 ± 0.24 ms, R_{1}: 35 ± 15 MΩ) component. Such a response can be understood in terms of a two compartment circuit consisting of a near compartment, comprising the patched soma and its surrounding, coupled to a far, mostly dendritic, compartment as depicted in Figure 3A (for details on the mapping, see "Capacitance measurements" in Methods). Importantly, the slow and fast components can be mapped to the corresponding five circuit parameters: near capacitance ${C}_{n}$ (21.0 ± 9.4 pF), near resistance ${R}_{n}$ (854 ± 394 MΩ), coupling resistance ${R}_{a}$ (53 ± 20 MΩ), far capacitance ${C}_{f}$ (106 ± 33 pF), and far resistance ${R}_{f}$ (156 ± 60 MΩ) (Figure 3C). Accordingly, this nearsomatic capacitance ${C}_{n}$ represents the summed capacitance of the membrane area that is electrotonically close to the recording site and thus is the value that the CapClamp requires as input and should be able to modify.
Altered nearsomatic capacitance in DGGCs
To confirm the localized effect of the CapClamp, we repeated the above capacitance measurement while clamping DGGCs at values ranging from 0.6 to 3 times the original near capacitance. Figure 3B depicts how the charging of the membrane potential in an exemplary cell changed its shape in reaction to the clamp. Both slow and fast time constant lengthened with capacitances, whereas the associated resistances increased and decreased, respectively, such that their sum, the total input resistance (which is expected to be independent of capacitance), remained constant. These measured time constants and amplitudes matched the predicted ones for a two compartment circuit with a near capacitance at the chosen target values and all other circuit parameters at their original values. In a multicompartment simulation of a morphologically reconstructed DGGC, we could reproduce both the two compartment structure of DGGCs and the isolated modification of the near capacitance, further confirming the local control via the CapClamp.
Across 18 recorded cells, the CapClamp robustly altered DGGC charging curves and modified their charging time constants. Within the tested capacitance range, the slow time constant ${\tau}_{0}$ decreased by –0.8 (1.0 to 0.6) ms, median and interquartile range in parentheses, and increased up to 3.0 (2.4 to 3.9) ms, whereas the fast time constant ${\tau}_{1}$ changes ranged from –0.24 (0.29 to 0.20) ms up to 0.60 (0.36 to 0.86) ms (Figure 3D). To quantify how well these changes reflected an altered near capacitance, we evaluated the goodness of fit between the observed and expected time constants and resistances. In the majority of cells, Rsquared values were close to 1, indicating that the CapClamp induced the expected changes (${\tau}_{0}$: 0.87 (0.76 to 0.92), R_{0}: 0.77 (0.56 to 0.89), ${\tau}_{1}$: 0.76 (0.32 to 0.97), R_{1}: 0.85 (0.75 to 0.91)). The largest mismatches occurred for the fast time constant, especially at high capacitances, where the measured time constant was often shorter than predicted (Figure 3D). A small bias toward a shorter fast component is to be expected and also present in the multicompartment simulation, because this time constant was only about ten times longer than the sampling interval of 50 µs limiting its slowingdown by the CapClamp currents. Larger deviations of ${\tau}_{1}$ however could not be reproduced in numerical simulations and likely result from other error sources, such as the difficulty of fitting this small and short time constant in the presence of noise or imprecise estimates of the original near capacitance (see "Online measurement of capacitance" in Methods). Overall, in terms of circuit parameters, the capacitance measurements confirmed the targeted near capacitance change for 12 out of 18 cells within an average error of 10% (Figure 3E). In summary, the CapClamp achieved an isolated change of the nearsomatic capacitance in DGGCs and thereby allows to control the time constants of their passive voltage dynamics.
Nearsomatic capacitance governs action potential shape and firing frequency in dentate gyrus granule cells
In neurons such as the recorded DGGCs, where the axon directly emerges from the soma, the ability to clamp the nearsomatic capacitance provides control over the major capacitive load for the action potential generating site in the axon initial segment. Consequently, the CapClamp, although acting locally, is expected to impact action potential (AP) dynamics and excitability of a morphologically complex DGGC as demonstrated earlier for the simplified singlecompartment neuron model (Figure 2). To illustrate how the CapClamp can be applied to characterize neuronal firing, we compared spiking responses and fI curves across near capacitances ranging from 0.6 to 3 times the original value, corresponding to a range from 10 pF to 60 pF for the near and from 110 pF to 160 pF for the total (near and far) capacitance.
Clamping the nearsomatic capacitance in DGGCs, we observed pronounced changes in the spiking response to depolarizing current step, clearly visible in the raw voltage traces (Figure 4A). The most apparent change was an altered AP shape (Figure 4B) – a continuous reduction of AP peak amplitude (from 60 ± 10 mV at 0.6 ${\mathrm{C}}_{\mathrm{n}}$ to 22 ± 17 mV at 3 ${\mathrm{C}}_{\mathrm{n}}$ for 9 DGGCs) and a simultaneous broadening of AP width (from 0.78 ± 0.15 ms at 0.6 ${\mathrm{C}}_{\mathrm{n}}$ to 1.33 ± 0.48 ms at 3 ${\mathrm{C}}_{\mathrm{n}}$) with increasing capacitance (Figure 4C and D). In addition, fast afterhyperpolarization (fAHP) was diminished and disappeared in the majority of cells after increasing capacitance (fAHP in 8/9 cells at 0.6 ${\mathrm{C}}_{\mathrm{n}}$ and 2/9 at 3 ${\mathrm{C}}_{\mathrm{n}}$). Importantly, the observed disappearance of fAHP cannot be explained by increased capacitive filtering alone, as an increased capacitance would reduce the fAHP amplitude, but not abolish it. Thus, our data suggests that the somatic capacitive load in DGGCs is able to influence the AP generating currents.
To illustrate the interplay of capacitance and the AP generating currents, we compared the observed spikes with hypothetical ones obtained by assuming unaltered currents with respect to those at the original capacitance (see "Protocol 2: Analysis of fI curves and spike shapes" in Methods). Recorded and hypothetical spike shapes exhibited marked differences (Figure 4B). At 0.6fold decreased capacitances, for example, the recorded AP amplitude was significantly smaller than the hypothetical one (rec.: 60 ± 10 mV, hyp.: 94 ± 19 mV, onesided Wilcoxon signedrank $Z$=0, p<0.001), presumably reflecting a reduction of the driving force for the sodium current when the AP peak approaches the reversal potential of sodium. Furthermore, at threefold increased capacitance, as noted above the recorded spikes exhibited no fAHP in most cells while the hypothetical ones still did (fAHP rec: 2/9, hyp: 8/9) – potentially a result of a reduced activation of potassium channels due to lower AP amplitudes and/or earlier closing during the slowed AP repolarization. In contrast to driving force and gating dynamics, the channel kinetics, for example their activation curves, cannot be altered by capacitance. Correspondingly, the spike threshold, which reflects the voltage where sodium channels start to massively open, was not significantly correlated with near capacitance (Pearson correlation r = 0.10, p = 0.42). Taken together, our analysis indicates that an altered somatic capacitance affects both sodium and potassium currents underlying APs in DGGCs.
Nearsomatic capacitance also impacted DGGC excitability. With increasing capacitance, DGGCs became less excitable and firing frequencies significantly decreased (Figure 4D and F). From 0.6 to 3fold of the original near capacitance, the decrease was modest for low firing rates close to threshold (from 9.7 ±3.2 Hz to 7.8 ±3.9 Hz, Wilcoxon signedrank $Z$=45, p=0.002) and became more pronounced for high firing rates at the largest injected currents (from 23.3 ±6.4 Hz to 18,6 ±4.8 Hz, $Z$=45, p=0.002). In terms of the firing ratecurrent (fI) curves, the gain of the DGGCs significantly decreased with capacitance (from 1.82 ± 0.40 Hz/$\sqrt{\mathrm{pA}}$ at 0.6 ${C}_{n}$ to 1.48 ± 0.34 Hz/$\sqrt{\mathrm{pA}}$ at 3 ${C}_{n}$, $Z$=45, p=0.002), whereas the rheobase current remained relatively constant (from 185 ±82 pA at 0.6 ${C}_{n}$ to 184 ±77 pA at 3 ${C}_{n}$, twosided, $Z$=12, p=0.25). Across cells, the gain reduction obtained by linear regression was –0.10 (0.13 to 0.06) Hz/$\sqrt{\mathrm{pA}}$ per 10 pF near capacitance (median and interquartile range, significant slope in 8/9 cells, p < 0.1). Compared to the simulated neuron with a gain reduction of –0.22 Hz/$\sqrt{\mathrm{pA}}$ per 10 pF over the same capacitance range as in the DGGC experiments (see Figure 2D), DGGCs thus exhibit a weaker gain dependence on near capacitance. A biological factor for this reduced effect in the DGGCs is their overall smaller gain – reflecting the different set of ionic conductances compared to the WangBuzsáki model designed to mimic a fastspiking cortical interneuron. Assuming a scaling of gain with 1/C (as predicted theoretically for neurons with a continuous fI curve, see "Analytically expected effect of capacitance on the form of the fI curve" in Methods), the gain reduction is expected to be gain(${\mathrm{C}}_{\mathrm{c}}$)/${\mathrm{C}}_{\mathrm{c}}$, which is approximately twice as high for the model ($\approx \frac{3.8\frac{\mathrm{Hz}}{\sqrt{\mathrm{pA}}}}{150\mathrm{p}\mathrm{F}}=0.25$ Hz/$\sqrt{\mathrm{pA}}$ per 10 pF) compared to the average DGGC ($\approx \frac{1.8\frac{\mathrm{Hz}}{\sqrt{\mathrm{pA}}}}{127\mathrm{p}\mathrm{F}}=0.14$ Hz/$\sqrt{\mathrm{pA}}$ per 10 pF). A further technical factor for a weaker effect in the DGGCs is the local nature of the capacitance modification. Depending on the particular location and geometry of the axon initial segment, the influence of the clamped somatic compartment on AIS excitability can differ (Goethals and Brette, 2020). The altered excitability in the majority of DGGCs, however, demonstrates that clamping their near capacitance was sufficient to affect the capacitive load of their AIS. We conclude that a somatic capacitance clamp, altering perisomatic capacitance alone, is able to modify the inputoutput relationship of a real neuron.
Applications of the CapClamp
The CapClamp lends itself to either test hypotheses on the impact of capacitance or to exploit the control over the membrane time constant in order to to alter neuronal dynamics in informative ways. In the following, we briefly illustrate applications of the CapClamp from these two fields, applying the technique to experimentally explore effects of capacitance on temporal integration and energetic costs of spiking.
Temporal integration
A basic processing step in neuronal computation is temporal integration, the summation of timeseparated synaptic inputs (Krueppel et al., 2011; Athilingam et al., 2017). An upper limit for temporal integration, at least in the absence of dedicated active channels, is set by the membrane time constant $\tau =RC$, which is directly proportional to the cell’s capacitance. Hence, increasing the capacitance of a cell should make it a better integrator: if two brief inputs arrive separated by less than the membrane time constant, the cell will summate the responses so that the membrane potential after the second is higher than after the first one. Indeed, when we compared the responses of DGGCs clamped at decreased and increased nearsomatic capacitances to current pulse trains, increasing the capacitance allowed the cell to better “sum” 3ms pulses at an inter stimulus interval (ISI) of 5 ms as apparent by the stairlike voltage response with a higher ratio of last to first pulse response. (Figure 5A and B). At an ISI of 50ms, in contrast, neither capacitance was sufficient for temporal integration. The biological relevance of tailoring capacitance to temporal processing can, for example, be observed in auditory cells of the barnowl, which have no dendrites to reduce capacitive and resistive load and hence shorten their time constant such that they can perform submillisecond coincidence detection (Ashida et al., 2007).
Energy consumption during spiking
Action potentials are energetically expensive, because the contributing sodium and potassium ions need to be pumped back using ATP (Laughlin et al., 1998; Hasenstaub et al., 2010). The minimal amount of ionic charge required for an action potential is dictated by the capacitance as $Q=C\mathrm{\Delta}{V}_{AP}$, suggesting that a smaller capacitance is energetically favorable. In order to gauge how capacitance affects charge accumulation and energy consumption, we reexamined spike shapes for a fixed current input at different capacitances both in the simulated neuron and in the recorded DGGCs (Figure 5C and E). We found that despite a reduced amplitude at larger capacitances, these smaller spikes still required more depolarizing charge ${Q}_{+}=C\mathrm{\Delta}{V}_{AP}$ (Figure 5D and F). In the model, we tested whether this depolarizing charge provided a reliable indication of the sodium charge ${Q}_{{\mathrm{Na}}^{+}}$, which finally determines pump activity and energy consumption (Figure 5D). Due to the overlap of sodium and potassium currents, the sodium charge exceeded the net depolarizing charge, but as this overlap remained roughly constant, both charge measures increased by the same amount with capacitance. Taken together, in the tested model and the recorded DGGCs, energy consumption per action potential appears to be reduced at smaller capacitances. In line with this observation, it has been reported that perineuronal nets could decrease membrane capacitance of fastspiking interneurons, thereby facilitating highfrequency firing, while keeping energetic costs at bay (Tewari et al., 2018).
Discussion
The dynamic clamp is a valuable tool in intracellular recordings to examine the diverse roles of ionic conductances in excitable cells (Sharp et al., 1993; Prinz et al., 2004; Wilders, 2006; Economo et al., 2010). In this study, we introduced the capacitance clamp (CapClamp), an application of the dynamic clamp that allows electrophysiologists to mimic a modified membrane capacitance in a biological neuron. Via simulations of a biophysical neuron model, we confirmed that the CapClamp correctly captures how capacitance affects spike shapes and firing frequency. In recordings of rat dentate gyrus granule cells, we further verified that the CapClamp could accurately control the capacitance of the recorded somatic compartments. Moreover, we clamped this nearsomatic capacitance of DGGCs during spiking and found that, as predicted by our simulations, capacitance can modify the fI curve and alter the course of the spike generating currents. CapClamp can serve as a new probe to neuronal signaling and physiology. In the following, we highlight requirements for the CapClamp and discuss how this experimental control over capacitance can benefit the study of cellular electrical behavior.
Precise, flexible and local control over capacitance in all excitable cells
To our knowledge, the CapClamp is the first tool to experimentally study capacitance changes in a precise and flexible manner. The CapClamp owes its precision and flexibility to the virtual nature of the altered capacitance. In contrast, methods to physically modify the capacitance are affected by various undesired side effects. Dendritic pinching, decoupling dendrites from the soma, for instance greatly reduces membrane area and thereby capacitance, but also removes all dendritic conductances (Bekkers and Häusser, 2007). Capacitance alterations have also been reported after application of mefloquine, a drug binding to membrane phospholipids, but it also blocks gap junctions (Szoboszlay et al., 2016). A notable exception is the recent demonstration of engineered polymer synthesis in neuronal cell membranes, which alters their capacitance, but not their input resistance (Liu et al., 2020). In comparison, however, the CapClamp provides more accurate and dynamic control by allowing experimenters to test multiple selected capacitance values in a single cell (Figures 3 and 4), while being significantly simpler to implement.
The CapClamp can be applied in every excitable cell. Here, we focused on neurons, but the proposed clamping currents can also be used to study capacitance changes in other cells, including for example heart cells (Wilders, 2006; de Oliveira et al., 2015). In particular, no prior knowledge about the ionic or external currents in the clamped cell is required, so that the capacitance can be clamped during any experimental protocol (step current, ramp current, etc.) or during synaptic input. Furthermore, capacitance can be clamped in both electrotonically compact cells like oocytes (Ori et al., 2020) and noncompact cells like most neurons (Wybo et al., 2019), although in the latter case the CapClamp is limited locally to the capacitance of the recorded compartment (Figure 3). Consequently, the effects of clamping capacitance depend on the cell’s morphology and the recording site. The soma, for example, represents the major capacitive load for spike generation in vertebrate neurons, where the axon predominantly emerges close to the soma (Figure 4), but it is expected to exert less influence in neurons, where the axon comes out of the dendritic tree, a common feature of invertebrate neurons (Hesse and Schreiber, 2015), but also seen in mammalian neurons (Martina et al., 2000; Thome et al., 2014).
The major prerequisite to apply the CapClamp is a reliable capacitance measurement of the clamped compartment, which can be challenging, especially for electrotonically complex cells (Golowasch et al., 2009; White and Hooper, 2013). An imprecise capacitance estimate leads to erroneous clamping currents, which increase highfrequency noise for small errors and might even induce instabilities for larger errors. The measurement method presented for the recorded DGGCs, that is mapping the charging response to a two compartment circuit, could in principle be extended to cells with a larger number of compartments e.g. pyramidal cells (Edwards and Mulloney, 1984; Wybo et al., 2021). Yet, accurate multiexponential fitting is demanding and the assumption of uniform membrane properties underlying the mapping is a simplification, shown to be violated in some cells, such as GABAergic interneurons (Nörenberg et al., 2010). As an alternative, measurement protocols could be exploited that inherently yield local capacitance estimates, including fast voltage ramps (Golowasch et al., 2009) or sampling of voltage responses to fast fluctuating currents (Badel et al., 2008). Reliable capacitance measurements can further be used to compare measured and target capacitance of the clamped cell, which can serve as a first simple test to ensure the quality of the CapClamp.
A CapClamp on every rig
As a novel application of the established dynamic clamp technique, the CapClamp is an accessible and lowcost extension of a standard electrophysiology stack (Prinz et al., 2004; Economo et al., 2010). For an existing dynamic clamp setup, the sole requirement is to implement the calculation of the clamping currents (see Equation 5). Otherwise, multiple open source frameworks exist that only require a dedicated computer with a data acquisition card to enable the dynamic clamp in a conventional electrophysiology setup (Dorval et al., 2001; Benda et al., 2007; Kemenes et al., 2011; Linaro et al., 2015; Patel et al., 2017; Desai et al., 2017; Amaducci et al., 2019). To facilitate the usage of the technique, we provide code for the CapClamp scheme in the RELACS and RTXI frameworks (see "Data and software availability" in Appendix 1).
In CapClamp recordings, as in all dynamic clamp applications, a high sampling frequency and accurate voltage monitoring are key (Bettencourt et al., 2008). Whether a sampling frequency is sufficiently high can be tested by assuring that the observed voltage dynamics for example the spike amplitudes are invariant when the sampling frequency is decreased from the maximal possible value (Robinson, 1994). For the simulated fastspiking interneuron, we found a satisfactory clamp at a frequency of 20 kHz, which we expect to also be sufficient for most excitatory neurons, because they tend to have slower voltage dynamics (Hasenstaub et al., 2010). In our single electrode recordings, we focused on careful electrode compensation to avoid electrode artifacts in the recorded voltages which would lead to incorrectly estimated membrane currents and eventually instabilities. To improve voltage monitoring, future applications could either apply active electrode compensation (Brette et al., 2008; Bal and Destexhe, 2009) or resort to two electrode recordings, where current injection and voltage recordings are separated.
Modifying capacitance as a probe for cellular electrical dynamics
Via the CapClamp, experimenters can ask a question that was previously accessible only in theoretical work: What if capacitance was different? In contrast to the theoretical approach, the answers to this question do not have to rely on models of channel dynamics or other membrane properties, because the latter are provided by the biological cell itself (Sharp et al., 1993). Modifying capacitance with the CapClamp can serve either to investigate changes in this biophysical parameter or, more broadly, to alter the membrane time constant of a cell as a way to characterize its electrical dynamics.
Understanding the role of capacitance
The virtual capacitance changes induced by the CapClamp could serve to address two crucial questions about actual membrane biophysics: why capacitance appears to be biologically mostly constant (Gentet et al., 2000) and how exceptions to this rule can facilitate or deter neuronal function (Amzica and Neckelmann, 1999; Hartline and Colman, 2007; Eyal et al., 2016; Tewari et al., 2018). Capacitance is for example rarely tested for optimality  a common question in ion channel kinetics, which appear optimized for function and energy expenditure (Hasenstaub et al., 2010; Sengupta et al., 2010). Regarding energy consumption, our CapClamp experiments in DGGCs indicate that action potentials become energetically cheaper at lower capacitances (Figure 5E and F). Interestingly, reports of exceptional capacitance values mostly find reductions e.g. for myelinated axons (${C}_{m}\approx $ 0.05 uF/cm^{2} for a 10fold wrapped myelin sheath, see Castelfranco and Hartline, 2015) or human pyramidal cells (${C}_{m}\approx $ 0.5 uF/cm^{2}, see Eyal et al., 2016) suggesting that indeed the metabolic cost of AP generation could have been a contributing factor to capacitance adaptations. In addition, the recent hypothesis that perineuronal nets can reduce capacitance of interneurons in a similar way as myelination of axons suggests that capacitance adaptation could be more widespread in the brain than often assumed (Tewari et al., 2018). Moreover, understanding the role of capacitance can contribute to an improved understanding of infrared (Shapiro et al., 2017; CarvalhodeSouza et al., 2018) and ultrasonic (Krasovitski et al., 2011; Plaksin et al., 2014) stimulation of neural activity, whose effects are assumed to rely on rapid alteration of the capacitance.
Another application of the CapClamp might be to investigate changes of excitability associated with neuronal growth. During development, cell size can increase considerably, necessarily accompanied by a larger membrane capacitance (McComb et al., 2003). To maintain neural function, neurons need to compensate for this altered capacitance via a corresponding regulation of ionic conductances – a homeostatic process that is hypothesized to involve activitydependent channel expression (GorurShandilya et al., 2020). To disentangle the contributions of capacitance in this concurrent alteration with ionic conductances, the CapClamp could be combined with the ‘classic’ dynamic clamp, for example parallel changes of both capacitance and a leak conductance.
Altering the membrane time constant
A key contribution of the CapClamp is the isolated experimental control of the membrane time constant via changes in capacitance while leaving the ion channel conductances unaffected. In neuron models, monitoring response properties when changing the membrane time constant has been used to characterize a cell’s dynamical repertoire (Kirst et al., 2015; Hesse et al., 2017; Franci et al., 2018). As an experimental analogue, the CapClamp introduces this option for the characterization of biological neurons.
To optimally support neural processing, nerve cells exhibit qualitatively different response properties, which in some cases can be flexibly adapted to context. For example, neurons with class 2 excitability (marked by a jump of the fI curve to nonzero frequencies when exceeding threshold) can be switched to class 1 excitability (marked by a smooth transition with arbitrarily low frequencies) via neuromodulation (Stiefel et al., 2008; Stiefel et al., 2009), transforming them from resonators to integrators. These qualitative differences in response and processing properties can be characterized by bifurcation analysis (Izhikevich, 2006; Prescott et al., 2008a; Kirst et al., 2015; Hesse et al., 2017). Capacitance as a canonical parameter can induce transitions between excitability classes and the underlying bifurcation types, including the switch of neuronal dynamics from class 1 excitability (with regularly spiking neurons) to dynamics that include bistable firing with stochastic switches between spiking and rest (Hesse et al., 2017). Because computational properties can be expected to change with such qualitative switches in dynamics, it may be of interest to determine how close the dynamics of a given cell is to a transition. An estimate of this proximity to switches that can be obtained via the CapClamp by monitoring firing properties and qualitative changes thereof as a function of membrane capacitance. Dynamics in the vicinity of capacitanceinduced switches are likely to be also susceptible to switches induced by other parameters with similar temporal effects, such changes in temperature (Hesse et al., 2017) or ionic concentrations (Contreras et al., 2020). As such switches can involve regimes of exceptionally fast dynamics, for such measurements extra care should be given to ensure that the temporal resolution of the dynamic clamp is sufficiently high.
In addition to such qualitative changes of dynamics, the broad impact of the time constant (and therefore the capacitance) on firing frequency and spike shape could be applied for more quantitative studies of neuronal activity. On the one hand, observations of neural activity at different capacitances could for example be used to further constrain and improve fitting of conductancebased neuron models (Podlaski et al., 2017; Gouwens et al., 2018; Franci et al., 2018). On the other hand, it could serve to examine activitydependent physiological processes such as ion concentration dynamics (Contreras et al., 2020) or calcium controlled channel homeostasis (O’Leary et al., 2014; Temporal et al., 2014; Santin and Schulz, 2019).
Conclusion
Taken together, the presented CapClamp enables an accurate and flexible control over capacitance in biological neurons, a basic determinant of cellular excitability, that so far has been inaccessible in experiment. We expect that the CapClamp will, therefore, broaden and enrich the electrophysiological study of neurons and other excitable cells. With expanding techniques to sense and manipulate neural activity, the combination of modeling and targeted closedloop feedback that underlies the CapClamp (and more generally the dynamic clamp Chamorro et al., 2012) will further unlock experimental control over other previously inaccessible aspects of single neuron (Ullah and Schiff, 2009; Rivera et al., 2015; Harrigan et al., 2018) and network dynamics (Newman et al., 2015; Hocker and Park, 2019).
Materials and methods
Derivation of the capclamp current
Request a detailed protocolIn order to derive a dynamic clamp feedback scheme for the CapClamp, we compare the actual membrane potential dynamics at the original capacitance ${C}_{c}$ with the target dynamics at the chosen capacitance ${C}_{t}$. The actual dynamics of the cell, which for the moment is assumed to be isopotential, is given by the currentbalance equation of a single compartment
with capacitance ${C}_{c}$, membrane currents $I(V,t)$ (comprising all ionic and synaptic currents, as well as external stimuli) and the dynamic clamp current ${I}_{\mathrm{dyn}}(t)$. Note that ionic and synaptic contributions to the membrane currents $I(V,t)$ are voltagedependent, both with respect to driving force and gating dynamics, so that a voltage trajectory governed by a different capacitance also leads to a modified shape of the membrane currents. In the target dynamics, the dynamic clamp current is absent and the capacitance is modified to the desired value
Both membrane potential trajectories would coincide, if we chose a dynamic clamp current such that the righthand sides of actual (Equation 2) and target dynamics (Equation 3) become identical,
Generally, an exact model for the membrane currents $I(V,t)$ will not be available, as it would require knowledge about all active conductances and incoming synaptic inputs. Instead, the membrane current can be estimated from the stream of incoming voltage data using the discrete version of Equation 2
where $\mathrm{\Delta}t$ is the sampling interval. A prerequisite is the measurement of the cell capacitance ${C}_{c}$. Furthermore, for the estimation to be accurate, the samplin ginterval needs to be shorter than the fastest time scales of changes in the membrane currents for example sodium gating time constants. With this estimated membrane current, the complete expression for the CapClampcurrent reads
The above derivation assumes that the cell is isopotential. In the case of an electrotonically noncompact cell, the steps are identical, but the cell capacitance ${C}_{c}$ has to be replaced by the capacitance of the compartment where the recording electrode is located. Consequently, in a nonisopotential neuron, the mimicked capacitance modification is restricted to the compartment at the tip of the recording electrode  a constraint known as the space clamp that is shared by all clamping techniques (Prinz et al., 2004; BarYehuda and Korngreen, 2008).
The indexing above assumes a voltage sampling ${V}_{i}=V(i\mathrm{\Delta}t)$ and a quasiimmediate current injection ${I}_{\mathrm{d}\mathrm{y}\mathrm{n},i}={I}_{\mathrm{d}\mathrm{y}\mathrm{n}}(i\mathrm{\Delta}t)$ . However, sampling can take a nonnegligible amount of time, so that depending on the sampling system the currently available voltage actually represents the voltage from the previous cycle ${V}_{i}=V((i1)\mathrm{\Delta}t)$ . In this case, for a correct estimation of the membrane currents, the dynamic clamp current index has to be shifted correspondingly to ${I}_{\mathrm{d}\mathrm{y}\mathrm{n},i}=\frac{{C}_{c}{C}_{t}}{{C}_{t}}\left({C}_{c}\frac{{V}_{i}{V}_{i1}}{\mathrm{\Delta}t}{I}_{\mathrm{d}\mathrm{y}\mathrm{n},i2}\right)$.
Capacitance measurements
To apply the CapClamp, a prerequisite is to measure the capacitance of the recorded local compartment. Here, we use the current clamp protocol, which estimates the capacitance from the voltage response to a current step with amplitude ${I}_{ext}$,
where an ordering in terms of these time scales is assumed i.e. ${\tau}_{0}>\tau 1>\dots$. Depending on the morphology, this sum can have a large number of components (Major et al., 1993), but in practice often only two or three components can be reliably extracted. As described in Golowasch et al., 2009, the slowest component ${\tau}_{0}$ is the membrane time constant and allows to infer the total capacitance of a neuron by $C=\frac{{\tau}_{0}}{{R}_{0}}=\frac{{\tau}_{0}}{{v}_{0}}{I}_{ext}$. In the case of an isopotential cell, the membrane time constant is the only component in the charging curve and the total capacitance can be used for the CapClamp.
Measurement of near capacitance
Request a detailed protocolFor the case of two components ${\tau}_{0},{R}_{0}$ and ${\tau}_{1},{R}_{1}$ in the charging curve (Equation 6), an equivalent two compartment circuit can be identified comprising a near compartment with capacitance ${C}_{n}$ and resistance ${R}_{n}$ connected via a coupling resistance ${R}_{a}$ to a far compartment with capacitance ${C}_{f}$ and resistance ${R}_{f}$ (Golowasch et al., 2009). With the additional assumption of a uniform membrane time constant ${\tau}_{m}={R}_{n}{C}_{n}={R}_{f}{C}_{f}$, the fitted two components can be mapped to the values of these five circuit parameters, which in particular provides the near capacitance ${C}_{n}$ required for the CapClamp
When the capacitance is subsequently clamped to a kfold different value, ${C}_{t}=k{C}_{n}$, the uniformity assumption has to be correspondingly adjusted to ${R}_{n}{C}_{n}=k{R}_{f}{C}_{f}$ (see "Mapping between a charging curve with two components and a two compartment circuit" in Appendix 1).
CapClamp in dentate gyrus granule cells
Electrophysiology
Request a detailed protocolAcute brain slices were produced as described earlier (Booker et al., 2014). Briefly, rats were anesthetized (3% Isoflurane, Abbott, Wiesbaden, Germany) and then decapitated. Brains were removed quickly and transferred to carbogenated (95% O_{2} / 5% CO_{2}) icecold sucroseACSF containing (in $\mathrm{mM}$): 87 NaCl, 2.5 KCl, 25 NaHCO_{3}, 1.25 NaH_{2}PO_{4}, 25 glucose, 75 sucrose, 7 MgCl_{2}, 0.5 CaCl_{2}, 1 Napyruvate, 1 ascorbic acid. Horizontal brain slices of 300 µm thickness were cut using a Vibratome (VT1200 S, Leica, Wetzlar, Germany). Hippocampal tissue slices, were collected and placed in a submerged holding chamber filled with carbogenated sucrose ACSF at 3234 °C for 30 min and then at room temperature for 15 min before recording. Experiments were alternated between left and right hemisphere slices to prevent bias due to slice condition.
For recording, slices were transferred to a submerged chamber and superfused with prewarmed, carbogenated ACSF containing (in $\mathrm{mM}$): 125 NaCl, 2.5 KCl, 25 NaHCO_{3}, 1.25 NaH_{2}PO_{4}, 25 glucose, 1 MgCl_{2}, 2 CaCl_{2}, 1 Napyruvate, 1 ascorbic acid. The bath temperature was set to 3234 °C with a perfusion rate of 1213 ml/min. Slices were visualized using an upright microscope (AxioScope; Zeiss) equipped with infrared differential inference contrast optics and a digital camera (Retiga EX QImaging CCD, Teledyne Photometrics, AZ, USA). Granule cells from the DG were chosen based on their anatomical location within the cell body layer as well as their morphological appearance.
Wholecell patchclamp electrodes were produced from borosilicate glass capillaries (outer diameter , inner diameter 1 mm, Hilgenberg, Germany) using a horizontal puller (P97, Sutter Instruments, CA, USA) and filled with an intracellular solution consisting of (in $\mathrm{mM}$): Kgluconate 130, KCl 10, HEPES 10, EGTA 10, MgCl_{2} 2, Na_{2}ATP 2, Na_{2}GTP 0.3, Na_{2}Creatine 1 and 0.1% biocytin (adjusted to pH 7.3 and 315 mOsm), giving a series resistance of 2.54 MΩ. All recordings were performed with a SEC LX10 amplifier (npi electronic, Germany), filtered online at 20 kHz with the builtin Bessel filter, and digitized at 20 kHz (National Instruments, UK). Following breakthrough into wholecell configuration, we adjusted the bridge and capacitance compensation before switching to the dynamic clamp mode for recording. Cells were excluded if resting membrane potential was more depolarized than 45 mV. The liquid junction potential was not corrected.
Neuronal visualization and immunohistochemistry
Request a detailed protocolFollowing recording, selected cells were immersion fixed in 4% paraformaldehyde (PFA) in 0.1 M phosphate buffer (PB, pH 7.4) at 4 °C for 24–48 hr, slices were then transferred to fresh PB. Prior to immunohistochemical processing, slices were rinsed in PB, followed by PB buffered saline (PBS, 0.9% NaCl). Slices were then rinsed in PBS and incubated in a fluorescentconjugated streptavidin (Alexa Fluor647, 1:1000, Invitrogen, UK) in PBS solution containing 3% NGS, 0.1% TritonX100 and 0.05% NaN3 for 24 hr at 4 °C. Slices were rinsed in PBS and then desalted in PB before being mounted (FluoromountG, Southern Biotech) on 300µmthick metal spacers, coverslipped, sealed, and stored at 4 °C prior to imaging.
Confocal imaging and reconstruction
Request a detailed protocolDGGCs were imaged on a laser scanning confocal microscope (FV1000, Olympus, Japan). First, a low magnification (4 x, Olympus, Japan) overview image was taken to confirm the cellular type and localization to the DG, then high resolution zstacks were obtained with a 30x silicone oil immersion objective (N.A. 1.05, UPlanSApo, Olympus) over the whole extent of the cell (1 µm axial steps). Image stacks were stitched offline using the FIJI software package (https://imagej.net/software/fiji/imagej.net), then the cells were reconstructed and volume filled using Neutube (https://www.neutracing.com/neutracing.com) (Feng et al., 2015).
Dynamic clamp setup
Request a detailed protocolData acquisition and dynamic clamp loop were controlled by RELACS, V0.9.8, RRID:SCR_017280 using a dedicated computer with a Linuxbased real time operating system (https://www.rtai.org/rtai.org). The sampling frequency was set to 20 kHz and the recordings were performed in discontinuous current clamp with a duty cycle of 16.5 µs. We implemented a CapClamp procedure for RELACS that allows the user to online specify the measured capacitance ${C}_{c}$ and the desired target capacitance ${C}_{t}$ (for documentation and installation instruction, see "Data and software availability" in Appendix 1).
Online measurement of capacitance
Request a detailed protocolFor the online measurement of the local capacitance, DGGCs were subjected to twenty hyperpolarizing pulses of $200\text{}\mathrm{ms}$ length with $400\text{}\mathrm{ms}$ pauses and an amplitude chosen to produce a response of $5\text{}\mathrm{mV}$ in order to minimize interference from active ionic currents. Responses were averaged and the resulting mean trajectory was fit with a sum of exponentials using the LevenbergMarquardt method from the python library scipy (Virtanen et al., 2020). Fits were performed with one, two and three components and were compared via the Fstatistic (Bardsley et al., 1986). In all recorded DGGCs, the two component fit was significantly better than the one exponential fit ($p<0.05$, 18/18), whereas no cell exhibited a significant third component ($p<0.05$, 0/18). Finally, the extracted two components were mapped to a two compartment circuit as explained above and the near capacitance was then used in the subsequent CapClamp (Table 2).
An offline reexamination revealed that in several recorded cells the above fitting procedure yielded inaccurate estimates of the exponential components, e.g. very short fast components due to an artefactual voltage dip before pulse onset. To circumvent these problems, improved offline fits were performed for the artifactfree recharging at the pulse end ( see "Adapted fitting procedure of dentate gyrus charging curves" in Appendix 1). In 8/18 cells, the offline and the original online estimate of the near capacitance differed by less than 20%, but overall the offline measurement yielded higher capacitance values than originally used for the CapClamp (offline: $21.0\pm 9.4\text{}\mathrm{pF}\text{/}$, online: $14.9\pm 4.8\text{}\mathrm{pF}\text{/}$). In contrast to the online measurement, the offline procedure reported a better fit with three components for a subset of cells ($p<0.05$, 7/18), but for the analysis presented here the result of the two component fit is used in all cells.
Protocol 1: Verification of altered capacitance
Request a detailed protocolAfter online measurement of the capacitance, each DGGC was clamped at a range of capacitances from 60% to 300% of the original near capacitance. For each clamped capacitance, the above offline capacitance measurement protocol was repeated to see how the CapClamp altered the slow and fast components. These time scale and amplitude changes were then mapped to the corresponding two compartment circuit parameters to compare them to the target capacitance (see Measurement of near capacitance). Due to the difference between online and offline estimate of the original near capacitance, we corrected the original target capacitance to ${C}_{t}^{\mathrm{corr}}={C}_{c}^{\mathrm{off}}+\mathrm{\Delta}{C}_{t}$, which preserves the targeted capacitance change $\mathrm{\Delta}{C}_{t}={C}_{t}{C}_{c}^{\mathrm{on}}$. Equally, the clamping factors in the mapping were updated to $k=\frac{{C}_{t}^{\mathrm{corr}}}{{C}_{c}^{\mathrm{off}}}$.
Protocol 2: Analysis of fI curves and spike shapes
Request a detailed protocolIn a subset of cells, after measuring near capacitance, an fI curve was obtained for the original capacitance and for target capacitances in the above range. Current pulses were 1 s long and repeated three times, at amplitudes ranging from 90% to 200% of an estimated rheobase. This rheobase was estimated by the first occurrence of spiking in response to a ramp (length: 5 s, height: 250 pA). For a quantitative comparison, the resulting fI curves were fit by a squareroot function
which captured their type 1 firing with a continuous frequencycurrent relationship (Izhikevich, 2006, p. 168). Cells with more than 30% varying input resistance within the protocol and/or a nonmonotonically increasing fI curves were excluded from the analysis.
Spikes were detected as a minimum 10 mV elevation over the average depolarization during the pulse. For the mean action potential (AP) shape, varying spike forms from the initial (< 300 ms) part of the pulse were discarded. The extracted AP features were peak amplitude, threshold voltage and threshold time to peak (voltage derivative crossing 10 mV/ms), height (difference between peak and threshold), temporal width at half of the height and fast afterhyperpolarization (fAHP; a voltage dip of –0.5 mV or larger within 10 ms after the spike). For threshold and fAHP detection, the spike shape was filtered with a digital 4th order Butterworth filter with critical frequencies 3.3 kHz, respectively 1 kHz.
To detect, whether changes in capacitance affect the action potential generating currents, we compared the recorded spikes with hypothetical ones obtained by assuming unaltered currents with respect to the original near capacitance. For a target capacitance ${C}_{t}$, such a hypothetical spike would be a scaled version of the original spike,
where ${V}_{c}(t)$ is the spike form at the original cell capacitance ${C}_{c}$ and the initial time t_{0} was chosen to be ${t}_{spike}3\mathrm{m}\mathrm{s}$ short before onset of the spike generating currents. Changes in the measured spike shape compared to this hypothetical shape signal a change of the underlying currents.
Simulations of the CapClamp
Simulations of neuron models coupled to the CapClamp were implemented using the neuron simulator Brian2 (Stimberg et al., 2019) and the CapClamp was realized using the Brian2 provided NetworkOperation that updated the clamp current every sampling interval using Equation 5 with zero delay between voltage sampling and current injection (for links to the available code, see "Data and software availability" in Appendix 1).
Biophysical neuron model
Request a detailed protocolIn order to test the CapClamp in the presence of active ionic conductances, a WangBuzsáki (WB) neuron, a single compartment model of hippocampal interneurons, was used (Wang and Buzsáki, 1996). Gating dynamics and peak conductances of the transient sodium current and the delayed rectifier potassium current were modeled as described earlier (Hesse et al., 2017, Appendix A). The specific membrane capacitance was chosen as $C}_{m}=0.75\phantom{\rule{thinmathspace}{0ex}}\frac{\mu \mathrm{F}}{\mathrm{c}{\mathrm{m}}^{2}$ and the membrane area was set to $A=20000\phantom{\rule{thinmathspace}{0ex}}{\mu m}^{2}$, so that the original cell capacitance was 150 pF. When the capacitance is varied, the WB neuron undergoes a wellcharacterized series of bifurcations; in particular it exhibits a saddlenode loop (SNL) bifurcation at $C}_{m}=1.47\phantom{\rule{thinmathspace}{0ex}}\frac{\mu \mathrm{F}}{\mathrm{c}{\mathrm{m}}^{2}$ accompanied by an abrupt doubling of the firing rate (Hesse et al., 2017). For the demonstration of the CapClamp here, we decided to restrict the tested capacitances to the regime below this critical value, but we confirmed via additional simulations that the CapClamp continues to work beyond the bifurcation (data not shown).
Simulations were performed with the second order RungeKutta method, a time step of 1 $\mu \mathrm{s}$ and dynamic clamp loop frequencies up to 100 kHz. Analysis of spike shapes and fI curves was performed in the same way as for the recorded cells.
Analytically expected effect of capacitance on the form of the fI curve
Request a detailed protocolHow the form of the fI curve depends on capacitance can be analytically calculated for a singlecompartment conductancebased neuron model undergoing a saddlenode on a limit cycle bifurcation at spiking onset like the WB model considered here (Izhikevich, 2006), pp. 162–168; (Schleimer and Schreiber, 2018). In this case, the time between two spikes ${T}_{\mathrm{isi}}$ is dominated by the slow traversal T_{2} of the saddle node, which close to threshold is multiple times longer than the brief duration T_{1} of the spike and can be derived by considering solely local dynamics
where $a$ and $c$ parametrize the normal form of the dynamics around the saddle node and ${I}_{\mathrm{rheo}}={I}_{\mathrm{sn}}$ is the current value where the saddle node bifurcation occurs. Inverting the inter spike interval to get the frequency then gives the square root form of the fI curve (see Equation 8). Under the assumption of fast gating kinetics, the $\mathrm{gain}=\frac{\sqrt{ac}}{\pi}$ is expected to be proportional to the inverse of the capacitance $\frac{1}{C}$, because the relevant time scale for the local slow dynamics around the saddle node is the membrane time constant implying that the traversal duration scales as ${T}_{2}\propto \tau \propto C$. The rheobase current in contrast is expected to remain constant, because equilibrium points are independent of the time scales of the dynamics. Formally calculating the normal form parameters $a$ and $c$ confirms these expectations (see "Formal derivation of fI curve gain and rheobase dependence oncapacitance" in Appendix 1).
The $\frac{1}{C}$ dependence of the gain allows to estimate an expected gain reduction for small capacitance changes around the original capacitance ${C}_{c}$ by Taylor expansion
which we compare for both simulated neuron and DGGCs to the observed gain reduction.
Multicompartment model of a dentate gyrus granule cell
Request a detailed protocolFor a controlled test of the CapClamp in an electrotonically noncompact cell, a morphologically reconstruction of a recorded DGGC was used as the basis for a multicompartment simulation. Soma and the two dendritic trees had a total area of 14,126 $\mu {\mathrm{m}}^{2}$. The axon was removed for the simulation. Membrane properties were assumed to be uniform and chosen such that they reproduced the average values of the total capacitance and the membrane time constant observed in the experiments: $C}_{m}=\frac{\overline{{C}_{n}+{C}_{f}}}{A}\approx 0.9\phantom{\rule{thinmathspace}{0ex}}\frac{\mu \mathrm{F}}{\mathrm{c}{\mathrm{m}}^{2}$ and $R}_{m}=\frac{{\tau}_{0}}{{C}_{m}}\approx 16800\phantom{\rule{thinmathspace}{0ex}}\mathrm{\Omega}{\mathrm{c}\mathrm{m}}^{2$. The axial resistivity was chosen as $R}_{\mathrm{a}\mathrm{x}\mathrm{i}\mathrm{a}\mathrm{l}}=300\phantom{\rule{thinmathspace}{0ex}}\mathrm{\Omega}\mathrm{c}\mathrm{m$. Simulations were performed with exponential Euler integration, a time step of 10 $\mu \mathrm{s}$ and a dynamic clamp sampling frequency of 20 kHz. Capacitance measurement and clamp procedure were the same as in the recorded DGGCs (Table 2).
Appendix 1
Impedance of a capacitanceclamped RC circuit
The impedance of a cell captures its linear response to the whole range of input frequencies (see Figure 1—figure supplement 1b). In the following, we derive the impedance of a passive membrane, an RC circuit, with capacitance ${C}_{c}$ coupled to the CapClamp and compare it to the impedance of an RC circuit with the target capacitance ${C}_{t}$.
Analysis of the dynamic clamp via the Ztransform
In general, the dynamic clamp technique forms a digital filter, mapping the incoming sampled voltages to injected currents. For a sampling interval $\mathrm{\Delta}t$, a linear mapping such as the CapClamp has the form
where $N$ and $M$ determine history of voltage and current values, respectively, taken into account. For the CapClamp, the coefficients depend on cell capacitance ${C}_{c}$, target capacitance ${C}_{t}$ and the sampling interval (see Equation 5),
This linear mapping can be represented and analyzed using the Ztransform (Dorf and Bishop, 2010, Ch. 13),
where the transfer function follows from the properties of the Ztransform: linearity $\lambda {X}_{i}\stackrel{Z}{\to}\lambda \hat{X}(z)$ and delay transformation ${X}_{i1}\stackrel{Z}{\to}{z}^{1}\hat{X}(z)$ (Dorf and Bishop, 2010, Table 13.2),
If the cell also forms a linear system, like the RC circuit, the transfer function of the coupled system Appendix 1—figure 1 is given by Dorf and Bishop, 2010, Table 2.6.
where ${H}_{\mathrm{cell}}(z)$ is the Ztransform of the membrane filter, e.g. ${H}_{\mathrm{cell}}(z)={H}_{\mathrm{RC}}(z)$.
The transfer function of the coupled system ${H}_{\mathrm{cell}+\mathrm{dyn}}(z)$ can then be compared with the one of the target system ${H}_{\mathrm{target}}(z)$ (Appendix 1—figure 1). Additionally, the frequencydependent impedance can be retrieved from the transfer function by
Transfer function of the CapClamp
The Ztransform of the CapClamp filter can be read directly from the general form of the transfer function (Equation 12) and the CapClamp feedback coefficients (Equation 10),
Transfer function of the RC circuit
In an RC circuit, the dynamics of the voltage are
Thus, in a single time step $\mathrm{\Delta}t$, when the current is fixed, the voltage evolves as
where $\tau =RC$ is the time constant. Applying the Ztransform results in the transfer function
which is subsequently used as the cell’s transfer function ${H}_{\mathrm{cell}}(z)={H}_{\mathrm{RC}}(z)$.
Transfer function of the clamped RC circuit
Introducing $K=\frac{{C}_{c}{C}_{t}}{{C}_{t}}$ and ${h}_{c}=\frac{\mathrm{\Delta}t}{{\tau}_{c}}$, the RC circuit (Equation 16) and CapClamp (Equation 15) transfer functions can be combined using Equation 13 to get the transfer function of the combined system
In comparison, the transfer function of the target RC circuit reads
with ${h}_{t}=\frac{\mathrm{\Delta}t}{{\tau}_{t}}=\frac{\mathrm{\Delta}t}{R{C}_{t}}$ reflecting the different target capacitance.
Figure 1—figure supplement 1b compares the resulting impedances for decreased and increased capacitances. As discussed in the "Results" section , the impedance amplitudes fit well up to a tenth of the dynamic clamp frequency. A closer look at the transfer function explains the fit at low frequencies and the deviations at higher frequencies.
Input resistance is preserved
The input resistance is equal to the impedance at zero frequency, that is at $z={e}^{i2\pi 0}=1$, which for both coupled and target system is the original resistance,
Poles and zeros
For a further comparison, poles and zeros of the transfer functions are calculated. To simplify the expressions, it is assumed that the time constant of the original and target circuits are much larger than the sampling interval, that is ${h}_{c}\ll 1$ and ${h}_{t}\ll 1$.
Target circuit
The target circuit has no zero and a single pole located at
Capacitance clamped circuit
The coupled system has one zero at
The clamped circuit has two poles at
and
Comparison of poles
All poles and zeros for an RC circuit in its original state and clamped at decreased and increased capacitances are shown in Appendix 1—figure 1. The first pole of the clamped circuit coincides with the one of the target circuit: $p}_{c}^{(1)}=1\frac{{C}_{c}}{{C}_{c}}\frac{\mathrm{\Delta}t}{R{C}_{c}}=1{h}_{t}={p}_{t}^{(1)$. As these pole lies close to $z=1$, they determine the lower frequency response, which explains why the impedance amplitudes fit so well in this range.
In addition to moving the existent pole of the cell circuit to the one of the target circuit, the CapClamp creates an additional pole ${p}_{c}^{(2)}\approx \frac{{h}_{c}}{2}(\frac{{C}_{c}}{{C}_{t}}1)$ and a new zero ${r}_{c}^{(1)}=1\frac{{C}_{c}}{{C}_{t}}$. Thus, at an increased capacitance $C}_{t}>{C}_{c$, the new pole lies in the left half of the unit circle and thereby increases the impedance at higher frequencies. In contrast, at a decreased capacitance, the additional zero moves into the left half of the complex plane and thereby decreases the impedance at higher frequencies.
Stability
For the investigated RC circuit with $R$=100 MΩ and $C$ = 150 pF and a sampling interval of 50 us, both poles of the capacitance clamped system remain within the unit circle (Appendix 1—figure 1) for the tested range from 0.1 to 10 times the original capacitance. As the coupled system is naturally causal, this implies that the transfer function of the clamped circuit is stable for this range of target capacitances, i.e. there are no unstable oscillations.
Mapping between a charging curve with two components and a two compartment circuit
In the following, we explain how a charging curve of a cell with two components can be mapped to the parameters of a two compartment circuit, which we used to extract the local capacitance in the recorded dentate gyrus granule cells (see Figure 3). We first report the approach and results derived earlier (Golowasch et al., 2009) and then explain how to extend the mapping when the capacitance is clamped to a modified value.
Golowasch et al. derived expressions for the near capacitance and the other circuit parameters by comparing the impedance of a two compartment circuit in Figure 3A
with the impedance of a system whose response to a step currents is a sum of two exponentials
The comparison of these two impedances gives four equations linking the circuit parameters and the two components of the charging curve:
To solve this set of equations, they assume that the membrane time constant is the same in all compartments ${C}_{n}{R}_{n}={C}_{f}{R}_{f}={\tau}_{c}$. However in a clamped neuron, where the near capacitance is targeted to be modified to a kfold different value, this equation becomes
where $k=\frac{{C}_{n,clam.}}{{C}_{n,orig.}}$.
For the unclamped case, $k=1$, the mapping from the two components to the circuit parameters is
For the clamped case, $k\ne 1$, we used the python package sympy to solve the equations.
Adapted fitting procedure of dentate gyrus charging curves
The initial online capacitance measurement was based on fitting the charging curve at the beginning of the current pulse. Posterior analysis showed an artefactual voltage drop of –0.2 mV starting about 0.2ms before pulse onset (probably due to coupling of the DAQ measurement card and the motherboard of the dynamic clamp computer), which limited the reliability of the online fit for cells with a small fast component. As no such artifact was observed for the recharging at the end of the pulse, this part was used in an improved offline fit. Additional measures to improve the fit were: cut of the first 0.2ms after pulse end to minimize electrode artifacts, limiting the fit to the first 60 ms (3–4 times ${\tau}_{0}$) after the pulse to prioritize the early part of the charging curve and a switch to the python package lmfit for better evaluation of parameter confidence bounds (https://lmfit.github.io/lmfitpy/). Furthermore, the finite rise time of the current injection by the amplifier was taken into account by adapting the original form of the charging curve (Equation 6) to
where the current rise time of the amplifier ${\tau}_{a}$ ($87\pm 2\phantom{\rule{thinmathspace}{0ex}}\mu \mathrm{s}$) was obtained by fitting the recorded injected current for the current step command by a simple exponential. A comparison of the two exponential components and the resulting circuit parameters for the online and offline fitting procedures is show in Appendix 1—table 1.
For the charging curves under capacitance clamp, the fitting procedure for the charging curve with two exponentials was initialized with values as expected for the targeted capacitance change: mapping the fitting results of the unclamped response to a two compartment circuit, changing the near capacitance to the targeted value and finally mapping this altered circuit back to the expected time scale and amplitudes. This initialization improved the fits especially at increased near capacitances, where the amplitude of the fast component becomes smaller.
Formal derivation of fI curve gain and rheobase dependence on capacitance
To confirm the expectations of gain and rheobase dependence on capacitance in a single compartment neuron model (see "Analytically expected effect of capacitance on the form of the fI curve" in Methods), we here sketch the calculation of the normal form parameters $a$ and $c$ following Izhikevich, 2006, pp. 162–163. In principle, $a$ and $c$ can be calculated for arbitrary gating kinetics by projecting the dynamics on the center manifold (Schleimer and Schreiber, 2018). For the assumption of small gating time constants, however, they can be expressed in simpler terms using the steady state IV relation of the neuron divided by its membrane capacitance
taking the form $a=\frac{1}{2}\frac{{\mathrm{\partial}}^{2}\mathbf{I}(V,I)}{\mathrm{\partial}{V}^{2}}{}_{V={V}_{\mathrm{s}\mathrm{n}},I={I}_{\mathrm{s}\mathrm{n}}}$ and $c=\frac{\mathrm{\partial}\mathbf{I}(V,I)}{\mathrm{\partial}I}{}_{V={V}_{\mathrm{s}\mathrm{n}},I={I}_{\mathrm{s}\mathrm{n}}}$, where the saddle node voltage $V}_{\mathrm{s}\mathrm{n}$ and current $I}_{\mathrm{s}\mathrm{n}$ are given by the equations
In summary, both $a$ and $c$ are proportional to $\frac{1}{C}$ and the rheobase current ${I}_{\mathrm{rheo}}={I}_{\mathrm{sn}}={I}_{\mathrm{\infty}}({V}_{\mathrm{sn}})$ is independent of capacitance, thus confirming the expected scaling.
Data and software availability
Electrophysiological recordings of capacitance clamped dentate gyrus granule cells: Pfeiffer, P., & Tomás F. J. B. (2021). Capacitance clamp demonstration in rat dentate gyrus granule cells. https://doi.org/10.5281/zenodo.5552207
Project repository with capacitance clamp module for https://scicrunch.org/resolver/RRID:SCR_017280RELACS and custom analysis/simulation in python: Pfeiffer, P., Tomás, F. J. B., Wu, J., Schleimer, J.H., Vida, I., & Schreiber, S. (2021). Software for: A dynamic clamp protocol to artificially modify cell capacitance. https://doi.org/10.5281/zenodo.6322768
Capacitance clamp plugin for http://rtxi.org/RTXI, a realtime dataacquistion and control application for biological research that allows to extend a conventional electrophysiology setup for dynamic clamp experiments (Patel et al., 2017). Capacitance_clamp_rtxi_module: https://doi.org/10.5281/zenodo.5553946
Data availability
All data generated, analysis code as well as computational modelling code is uploaded on https://zenodo.org/, see article section Data and software availability.

ZenodoCapacitance Clamp Demonstration in Rat Dentate Gyrus Granule Cells.https://doi.org/10.5281/zenodo.5552207
References

Membrane capacitance of cortical neurons and glia during sleep oscillations and spikewave seizuresJournal of Neurophysiology 82:2731–2746.https://doi.org/10.1152/jn.1999.82.5.2731

Passive soma facilitates submillisecond coincidence detection in the owl’s auditory systemJournal of Neurophysiology 97:2267–2282.https://doi.org/10.1152/jn.00399.2006

Dynamic IV curves are reliable predictors of naturalistic pyramidalneuron voltage tracesJournal of Neurophysiology 99:656–666.https://doi.org/10.1152/jn.01107.2007

BookDynamicClamp: From Principles to ApplicationsNew York, NY: Springer US.https://doi.org/10.1007/9780387892795

Spaceclamp problems when voltage clamping neurons expressing voltagegated conductancesJournal of Neurophysiology 99:1127–1136.https://doi.org/10.1152/jn.01232.2007

From response to stimulus: adaptive sampling in sensory physiologyCurrent Opinion in Neurobiology 17:430–436.https://doi.org/10.1016/j.conb.2007.07.009

Effects of imperfect dynamic clamp: computational and experimental resultsJournal of Neuroscience Methods 169:282–289.https://doi.org/10.1016/j.jneumeth.2007.10.009

Wholecell patchclamp recordings from morphologically and neurochemicallyidentified hippocampal interneuronsJournal of Visualized Experiments e51706.https://doi.org/10.3791/51706

The evolution of vertebrate and invertebrate myelin: A theoretical computational studyJournal of Computational Neuroscience 38:521–538.https://doi.org/10.1007/s108270150552x

Increased cell membrane capacitance is the dominant mechanism of stretchdependent conduction slowing in the rabbit heart: A computational studyCellular and Molecular Bioengineering 8:237–246.https://doi.org/10.1007/s1219501503849

A Dynamic Clamp on Every RigENeuro 4:ENEURO.025017.2017.https://doi.org/10.1523/ENEURO.025017.2017

RealTime linux dynamic clamp: A fast and flexible way to construct virtual ion channels in living cellsAnnals of Biomedical Engineering 29:897–907.https://doi.org/10.1114/1.1408929

Dynamic clamp: alteration of response properties and creation of virtual realities in neurophysiologyThe Journal of Neuroscience 30:2407–2413.https://doi.org/10.1523/JNEUROSCI.595409.2010

Compartmental models of electrotonic structure and synaptic integration in an identified neuroneThe Journal of Physiology 348:89–113.https://doi.org/10.1113/jphysiol.1984.sp015101

Robust and tunable bursting requires slow positive feedbackJournal of Neurophysiology 119:1222–1234.https://doi.org/10.1152/jn.00804.2017

Direct measurement of specific membrane capacitance in neuronsBiophysical Journal 79:314–320.https://doi.org/10.1016/S00063495(00)76293X

Myopic control of neural dynamicsPLOS Computational Biology 15:e1006854.https://doi.org/10.1371/journal.pcbi.1006854

A computational model for how the fast afterhyperpolarization paradoxically increases gain in regularly firing neuronsJournal of Neurophysiology 119:1506–1520.https://doi.org/10.1152/jn.00385.2017

Dynamic clamp with StdpC softwareNature Protocols 6:405–417.https://doi.org/10.1038/nprot.2010.200

Realtime Electrophysiology: Using Closedloop Protocols to Probe Neuronal Dynamics and BeyondJournal of Visualized Experiments 1:e52320.https://doi.org/10.3791/52320

Genetically targeted chemical assembly of functional materials in living cells, tissues, and animalsScience (New York, N.Y.) 367:1372–1376.https://doi.org/10.1126/science.aay4866

Distal initiation and active propagation of action potentials in interneuron dendritesScience (New York, N.Y.) 287:295–300.https://doi.org/10.1126/science.287.5451.295

Electrophysiological differences in the CPG aerial respiratory behavior between juvenile and adult LymnaeaJournal of Neurophysiology 90:983–992.https://doi.org/10.1152/jn.00263.2003

Hard realtime closedloop electrophysiology with the RealTime eXperiment Interface (RTXIPLOS Computational Biology 13:e1005430.https://doi.org/10.1371/journal.pcbi.1005430

Biophysical basis for three distinct dynamical mechanisms of action potential initiationPLOS Computational Biology 4:e1000198.https://doi.org/10.1371/journal.pcbi.1000198

Pyramidal neurons switch from integrators in vitro to resonators under in vivolike conditionsJournal of Neurophysiology 100:3030–3042.https://doi.org/10.1152/jn.90634.2008

The dynamic clamp comes of ageTrends in Neurosciences 27:218–224.https://doi.org/10.1016/j.tins.2004.02.004

Dendritic growth and regression in rat dentate granule cells during late postnatal developmentBrain Research. Developmental Brain Research 54:115–124.https://doi.org/10.1016/01653806(90)900716

Towards a dynamic clamp for neurochemical modalitiesSensors (Basel, Switzerland) 15:10465–10480.https://doi.org/10.3390/s150510465

Conductance injectionTrends in Neurosciences 17:147–148.https://doi.org/10.1016/01662236(94)900884

Phaseresponse curves of ion channel gating kineticsMathematical Methods in the Applied Sciences 41:8844–8858.https://doi.org/10.1002/mma.5232

Subthreshold dendritic signal processing and coincidence detection in dentate gyrus granule cellsThe Journal of Neuroscience 27:8430–8441.https://doi.org/10.1523/JNEUROSCI.178707.2007

Action potential energy efficiency varies among neuron types in vertebrates and invertebratesPLOS Computational Biology 6:e1000840.https://doi.org/10.1371/journal.pcbi.1000840

The dynamic clamp: artificial conductances in biological neuronsTrends in Neurosciences 16:389–394.https://doi.org/10.1016/01662236(93)900046

The effects of cholinergic neuromodulation on neuronal phaseresponse curves of modeled cortical neuronsJournal of Computational Neuroscience 26:289–301.https://doi.org/10.1007/s1082700801119

Frequencydependent regulation of intrinsic excitability by voltageactivated membrane conductances, computational modeling and dynamic clampThe European Journal of Neuroscience 46:2429–2444.https://doi.org/10.1111/ejn.13708

Cellular shortterm memory from a slow potassium conductanceJournal of Neurophysiology 75:963–966.https://doi.org/10.1152/jn.1996.75.2.963

Tracking and control of neuronal HodgkinHuxley dynamicsPhysical Review. E, Statistical, Nonlinear, and Soft Matter Physics 79:1–4.https://doi.org/10.1103/PhysRevE.79.040901

Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network modelThe Journal of Neuroscience 16:6402–6413.

Contamination of currentclamp measurement of neuron capacitance by voltagedependent phenomenaJournal of Neurophysiology 110:257–268.https://doi.org/10.1152/jn.00993.2012

Dynamic clamp: A powerful tool in cardiac electrophysiologyThe Journal of Physiology 576:349–359.https://doi.org/10.1113/jphysiol.2006.115840

Electrical Compartmentalization in NeuronsCell Reports 26:1759–1773.https://doi.org/10.1016/j.celrep.2019.01.074
Decision letter

Ronald L CalabreseSenior and Reviewing Editor; Emory University, United States

Jorge GolowaschReviewer; New Jersey Institute of Technology, United States

Ronald L CalabreseReviewer; Emory University, United States
Our editorial process produces two outputs: i) public reviews designed to be posted alongside the preprint for the benefit of readers; ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.
Decision letter after peer review:
Thank you for submitting your article "A dynamic clamp protocol to artificially modify cell capacitance" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Ronald L Calabrese as the Senior and Reviewing Editor and Reviewer #3. The following individual involved in review of your submission has agreed to reveal their identity: Jorge Golowasch (Reviewer #1).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1) As discussed by Reviewer #2 in their public review, the changes in capacitance have a weak effect on excitability in real neurons. The authors should definitely apply the same comparison already performed for action potential, i.e., expected changes vs. real changes in FI curve following capacitance clamp. Maybe in DGGCs, not much effect is expected on gain, but this has to be clearly demonstrated. Otherwise, it raises concerns about the efficiency of capacitance clamp beyond the soma.
2) Please address all the concerns expressed in the Recommendations for the authors.
Reviewer #1 (Recommendations for the authors):
This is an excellent paper. The new dynamic clamp method described here to control membrane capacitance is based on sound theory, it is well described and tested. I have no major concerns. In fact, I appreciate the authors for developing this in such a timely manner for me, as I was planning to do something similar myself.
Reviewer #2 (Recommendations for the authors):
1) The authors demonstrate in the neuron model that manipulating capacitance not only affects action potential waveform but also significantly alters the excitability profile of the cell, modifying in particular the gain of the fI curve of the neuron. However, when tested on DGGCs, the effect of capacitance on action potential shape is very strong but the effect on excitability is very mild. The overall change in gain is close to 20% for a 5fold change in capacitance, while a similar change in capacitance induced a ~2fold in gain in the neuron model. While the authors say that the results in real neurons are similar to the ones obtained in the simulated neuron, the quantitative difference is large enough to contradict that statement. Moreover, this discrepancy questions the ability of the capacitance clamp to efficiently modify capacitance in real neurons. In fact, concerning action potential shape, the authors compare the effect of capacitance manipulation in real neurons with the expected effect (Figure 4B), but do not present this comparison for excitability measurements. It would have been very interesting to see whether the actual results significantly depart from the expected results.
2) Concerning the potential impact of capacitance clamp, and since most changes in capacitance in physiological contexts seem to be related to neuronal growth, it would have been really interesting to compare the impact of manipulating capacitance with the impact of manipulating concomitantly input resistance and capacitance, which is expected when neuronal size is changing, for instance during development. This type of comparison would also help underlining the significant contribution of changes in capacitance, which have been so far undermined. Testing these two manipulations in parallel would greatly help to disentangle the specific contributions of changes in membrane resistance and capacitance during neuronal growth, and would emphasize the value of the capacitance clamp tool.
Reviewer #3 (Recommendations for the authors):
The manuscript is very clearly written and wellfocused.
Lines 9394: "…change with respect to the original capacitance (e.g. Ct=67.4 pF: C=67.5 pF; Ct=336.9 pF: C=338.1 pF), whereas the… Why are these examples chosen and not a more uniform range?
Section 2.5.1: I found this very confusing. In an RC circuit the response sizes to a sequence of current pulses as measured by baseline to peak are identical. It is very confusing to say "…the cell's response to the second one should be higher than to the first one." Or to say "…as apparent by the larger step sizes in the stairlike voltage response and the finally higher ratio of last to first pulse response." Please rewrite.
Section 2.5.3: This section was the only part of the paper I found unconvincing. It is a foregone conclusion that the CapClamp can find the critical capacitance in the WangBuzsáki neuron, given Section 2.3. The failure of the technique to find the critical capacitance in a dentate gyrus granule cell is thus an ambiguous result. Is this a technical failure or a real result? I suggest deleting this section.
https://doi.org/10.7554/eLife.75517.sa1Author response
Essential revisions:
1) As discussed by Reviewer #2 in their public review, the changes in capacitance have a weak effect on excitability in real neurons. The authors should definitely apply the same comparison already performed for action potential, i.e., expected changes vs. real changes in FI curve following capacitance clamp. Maybe in DGGCs, not much effect is expected on gain, but this has to be clearly demonstrated. Otherwise, it raises concerns about the efficiency of capacitance clamp beyond the soma.
We thank the reviewers for pointing out this apparent discrepancy between model and experiment. The discrepancy, however, is not as strong as the text may have suggested. In the revised manuscript, we now provide a quantitative analysis of the effects on excitability in the simulated neuron models and the dentate granule cells (DGGCs) as well as a comparison to expected changes based on theoretical predictions for interspike interval duration (for details see the full answer to Recommendation 1 of reviewer 2 on page 6). In summary, this reanalysis shows that the discrepancy is actually smaller than suggested by the comparison in terms of relative capacitance changes. Capacitance effects on DGGC excitability in terms of gain reduction are reduced by a factor of 0.5 in comparison to those in the model. This remaining discrepancy is consistent with theoretical predictions – less excitable cells should also exhibit a smaller gain reduction, when capacitance is increased. The additional gain analysis is thus further evidence that in the recorded DGGCs the somatic capacitance clamp was sufficient to effectively change the capacitive load of the spike generating compartment. Thus, we are confident that the CapClamp also provides this control in other cell types with spike initiation at or close to the soma.
2) Please address all the concerns expressed in the Recommendations for the authors.
We thank all reviewers for the thorough and helpful comments and addressed all concerns expressed in the Recommendations.
Reviewer #2 (Recommendations for the authors):
1) The authors demonstrate in the neuron model that manipulating capacitance not only affects action potential waveform but also significantly alters the excitability profile of the cell, modifying in particular the gain of the fI curve of the neuron. However, when tested on DGGCs, the effect of capacitance on action potential shape is very strong but the effect on excitability is very mild. The overall change in gain is close to 20% for a 5fold change in capacitance, while a similar change in capacitance induced a ~2fold in gain in the neuron model. While the authors say that the results in real neurons are similar to the ones obtained in the simulated neuron, the quantitative difference is large enough to contradict that statement. Moreover, this discrepancy questions the ability of the capacitance clamp to efficiently modify capacitance in real neurons. In fact, concerning action potential shape, the authors compare the effect of capacitance manipulation in real neurons with the expected effect (Figure 4B), but do not present this comparison for excitability measurements. It would have been very interesting to see whether the actual results significantly depart from the expected results.
We thank the reviewer for pointing out the apparent discrepancy of gain dependence on capacitance between model and experiment. This discrepancy, however, is actually smaller than it appeared in the original manuscript. As Figure 4 on the DGGCs emphasized the near capacitance, which the CapClamp modifies, our presentation suggested comparing gains in terms of relative changes of near capacitance in the experiment and total capacitance in the model. These relative changes of the modified capacitance, however, are not the correct base for a comparison, because in the case of the DGGCs they do not take into account the large constant (far) dendritic capacitance.
In the revised manuscript, we therefore provide an explicit comparison between the model and experiment in terms of absolute capacitance changes, i.e., for the model C_{t}C_{cell} (Figure 2 E and page 7) and for the DGGCs C_{t}C_{near} (Figure 4 D and pages 1213). In other words, we report how much the fI curve gain changes per a 10 pF increase of capacitance. In addition, we use an analytical argument to discuss how much the gain is expected to change. Our analysis shows that the dependence of gain on capacitance in the DGGCs is weaker by a factor of ~0.5 in comparison to that of the model. Such a remaining difference is expected, because the DGGCs have an overall smaller gain (for the detailed argument see below). An additional factor might be the degree of electric coupling of somatic and axonal compartment, which can vary among cells. Yet, a significant gain reduction in the majority of DGGCs – within the expected range – shows that somatic capacitance control was sufficient to affect spike generation. Overall, we are confident that this reanalysis further confirms the ability of the capacitance clamp to effectively modify capacitance in real neurons.
Measured gain dependence on near capacitance in DGGCs: In the DGGCs, we explored a range of 60% to 300% of the original near capacitance. As the average near capacitance was ~20 pF, the total capacitance of ~120pF (near+far) thus changed in a range from ~110 pF to ~160 pF, i.e., roughly a 1.5 fold change and not a 5fold change as for the near capacitance. To clarify this smaller relative change with respect to the total capacitance, we mention this range explicitly in the revised manuscript (lines 206209) and show the total capacitance in the exemplary cell in Figure 4 A. As a measure of gain dependence, we performed a linear regression of gain versus absolute capacitance change for each recorded DGGC, which yielded a gain reduction of 0.10 (0.13 to 0.06) Hz/pA^{1/2} per 10 pF (median and interquartile range for 9 cells, see lines 239240). Gain reduction in DGGCs and the simulated neuron thus differ by a factor of ~0.5 (see next section) – a more similar effect size in experiment and model than suggested by the comparison in terms relative near capacitances.
Comparison of DGGCs gain dependence to the WangBuzsáki model: The revised manuscript explicitly shows how the gain of the WangBuzsaki (WB) neuron model depends on capacitance (see Figure 2 E and updated Table 1). At the original capacitance, the gain of the WB model is 3.8 Hz/pA^{1/2}, which is about twofold higher than the gain of the unclamped DGGCs 1.76 (1.72 to 1.86) Hz/pA^{1/2}. Across the large capacitance range tested in the model from 90 pF to 210 pF, the gain roughly follows a 1/C behavior with a steep slope for small capacitances and flatter one at large capacitances. The comparison to the experiment thus depends on the capacitance range over which this reduction is evaluated. One fair way of comparison is to pick a range around the original capacitance of 150 pF corresponding to the experimentally accessible manipulations (20 to +40 pF). In this range, the WB model gain reduction was 0.22 Hz/pA^{1/2} per 10 pF. The gain of the DGGCs thus is less sensitive (0.10 Hz/pA^{1/2} per 10 pF), but this reduction to 50% is a much milder difference and can be partially accounted for by the different excitability of the two neuron types (see lines 240250).
Expected dependence of gain on (near) capacitance in DGGCs: What is the origin of the smaller gain reduction in the DGGCs compared to the model neuron? One difference between WB model and DGGC lies in the composition of ionic currents that together with the baseline capacitance determine the original gain of a neuron. From this original gain, a rough expectation for the gain reduction can be formulated assuming a scaling of gain with 1/C – a good approximation for class 1 singlecompartment neurons under certain conditions (Izhikevich, 2006, see added Methods section “Analytically expected effect of capacitance on the form of the fI curve”, pages 2526). Then close to the original capacitance C_{c}, the gain reduction should be gain(C_{c})/C_{c}. As the model overall has a higher gain but approximately the same total capacitance as the DGGCs, it is therefore expected to also have a higher gain reduction, Concretely, according to this backontheenvelope calculation, the model gain reduction should be around 3.8 Hz/pA^{1/2}/ 150 pF = 0.25 Hz/pA^{1/2} per 10 pF, twice as high as the DGGC expectation 1.8 Hz/pA^{1/2} / 127 pF = 0.14 Hz/pA^{1/2} per 10 pF. Concluding, the observed discrepancy between model and experiment in gain reduction is consistent with their difference in overall excitability.
In addition to the biological differences, a further technical factor potentially decreasing the measured gain dependence in the DGGCs is the space clamp problem. Location and geometry of the axon initial segment (AIS), as well as the distribution of active channels, determine how much the capacitance of the somatic compartment impacts AIS excitability (Goethals, 2020, eLife). As these factors might vary among cells, they could contribute to the observed variability in the impact of somatic capacitance. The observed significant gain reductions in the majority of DGGCs, together with the effects on the spike shape, however, demonstrate that clamping capacitance of the soma can effectively modify the capacitive load of the AIS in a real neuron.
2) Concerning the potential impact of capacitance clamp, and since most changes in capacitance in physiological contexts seem to be related to neuronal growth, it would have been really interesting to compare the impact of manipulating capacitance with the impact of manipulating concomitantly input resistance and capacitance, which is expected when neuronal size is changing, for instance during development. This type of comparison would also help underlining the significant contribution of changes in capacitance, which have been so far undermined. Testing these two manipulations in parallel would greatly help to disentangle the specific contributions of changes in membrane resistance and capacitance during neuronal growth, and would emphasize the value of the capacitance clamp tool.
We thank the reviewer for this interesting idea. Neuronal growth provides a good example of a process that cannot be studied by classic dynamic clamp alone, because a modified cell size inevitably implies an altered capacitance as well as homeostatic regulation of channel numbers. As proposed by the reviewer, an experimenter can combine the capacitance clamp with virtual conductance injection via classic dynamic clamp to map neuronal firing in the space of capacitance and leak conductance, effectively simulating the passive electrical changes associated with different growth scenarios.
Such a combined experiment is unfortunately beyond the scope of this introductory study of the capacitance clamp. But we included a section in the Discussion to emphasize this possible combination of capacitance clamp and classic dynamic clamp in the context of neuronal growth (lines 366372).
Reviewer #3 (Recommendations for the authors):
The manuscript is very clearly written and wellfocused.
Lines 9394: "…change with respect to the original capacitance (e.g. Ct=67.4 pF: C=67.5 pF; Ct=336.9 pF: C=338.1 pF), whereas the… Why are these examples chosen and not a more uniform range?
Throughout the paper we used the range explored in the experiments from 60% to 300% of the capacitance to be modified and the indicated absolute capacitance values represent the corresponding limits of the tested range. As mentioned in the Discussion, this was the range where we had a reliable stable capacitance clamp in experiment.
Section 2.5.1: I found this very confusing. In an RC circuit the response sizes to a sequence of current pulses as measured by baseline to peak are identical. It is very confusing to say "…the cell's response to the second one should be higher than to the first one." Or to say "…as apparent by the larger step sizes in the stairlike voltage response and the finally higher ratio of last to first pulse response." Please rewrite.
It is correct that the response size is the same if the length of a single pulse is much longer than the membrane time constant, so that the circuit reaches its stationary state. For short pulses, however, as considered here in the context of temporal integration, the response size of the later pulses increases due to temporal summation (see Figure 5 A). We rewrote the section clarifying the context of brief input pulses accordingly (lines 259270).
Section 2.5.3: This section was the only part of the paper I found unconvincing. It is a foregone conclusion that the CapClamp can find the critical capacitance in the WangBuzsáki neuron, given Section 2.3. The failure of the technique to find the critical capacitance in a dentate gyrus granule cell is thus an ambiguous result. Is this a technical failure or a real result? I suggest deleting this section.
One of our original motivations to develop the capacitance clamp was to extend the experimental tools available to systematically scan real neurons for interesting bifurcations. As some bifurcations like the presented saddlenode loop bifurcation involve fast dynamics (switch to the strongly attracting manifold, see Hesse, 2017), we thus wanted to ensure that dynamics in this regime can still be mimicked with the CapClamp. We agree, however, with the comment of the Reviewer that the presented simulation results in Figure 5 G and H (original manuscript) are difficult to interpret biologically and the interpretation of the absence of the critical transition in experiments remains ambiguous. We therefore followed the suggestion to delete this section yet took the liberty to include a few sentences in the Discussion outlining the option to detect dynamical switches with CapClamp (see lines 378392).
https://doi.org/10.7554/eLife.75517.sa2Article and author information
Author details
Funding
Bundesministerium für Bildung und Forschung (01GQ1403)
 JanHendrik Schleimer
 Susanne Schreiber
Deutsche Forschungsgemeinschaft (GRK 1589/2)
 Paul Pfeiffer
 Federico José Barreda Tomás
Deutsche Forschungsgemeinschaft (EXC 257)
 Federico José Barreda Tomás
 Imre Vida
Deutsche Forschungsgemeinschaft (FOR 2134)
 Federico José Barreda Tomás
 Imre Vida
H2020 European Research Council (864243)
 Susanne Schreiber
Einstein Stiftung Berlin (EZ2014224)
 Jiameng Wu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Jan Benda and Lukas Sonnenberg for their dedicated support with dynamic clamp and fruitful discussions. We are grateful to Eve Marder and Ekaterina Morozova for being able to test the RTXI implementation of the CapClamp and apply the CapClamp in neurons of the crustacean stomatogastric ganglion. We thank Robert Gowers and Philipp Norton for valuable feedback on the manuscript.
This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No 864243). The article processing charge was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – 491192747 and the Open Access Publication Fund of HumboldtUniversität zu Berlin.
Ethics
All procedures and animal maintenance were performaed in accordance with institutional guidelines, the German Animal Welfare Act, the European Council Directive 86/609/EEC regarding the protection of animals, and guidelines from local authorities (Berlin, T0215/11).
Senior and Reviewing Editor
 Ronald L Calabrese, Emory University, United States
Reviewers
 Jorge Golowasch, New Jersey Institute of Technology, United States
 Ronald L Calabrese, Emory University, United States
Publication history
 Received: November 12, 2021
 Preprint posted: November 13, 2021 (view preprint)
 Accepted: March 17, 2022
 Accepted Manuscript published: April 1, 2022 (version 1)
 Version of Record published: May 26, 2022 (version 2)
Copyright
© 2022, Pfeiffer et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 901
 Page views

 211
 Downloads

 0
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Neuroscience
One signature of the human brain is its ability to derive knowledge from language inputs, in addition to nonlinguistic sensory channels such as vision and touch. How does human language experience modulate the mechanism by which semantic knowledge is stored in the human brain? We investigated this question using a unique human model with varying amounts and qualities of early language exposure: early deaf adults who were born to hearing parents and had reduced early exposure and delayed acquisition of any natural human language (speech or sign), with early deaf adults who acquired sign language from birth as the control group that matches on nonlinguistic sensory experiences. Neural responses in a semantic judgment task with 90 written words that were familiar to both groups were measured using fMRI. The deaf group with reduced early language exposure, compared with the deaf control group, showed reduced semantic sensitivity, in both multivariate pattern (semantic structure encoding) and univariate (abstractness effect) analyses, in the left dorsal anterior temporal lobe (dATL). These results provide positive, causal evidence that language experience drives the neural semantic representation in the dATL, highlighting the roles of language in forming human neural semantic structures beyond nonverbal sensory experiences.

 Neuroscience
Across phyla, males often produce speciesspecific vocalizations to attract females. Although understanding the neural mechanisms underlying behavior has been challenging in vertebrates, we previously identified two anatomically distinct central pattern generators (CPGs) that drive the fast and slow clicks of male Xenopus laevis, using an ex vivo preparation that produces fictive vocalizations. Here, we extended this approach to four additional species, X. amieti, X. cliivi, X. petersii, and X. tropicalis, by developing ex vivo brain preparation from which fictive vocalizations are elicited in response to a chemical or electrical stimulus. We found that even though the courtship calls are speciesspecific, the CPGs used to generate clicks are conserved across species. The fast CPGs, which critically rely on reciprocal connections between the parabrachial nucleus and the nucleus ambiguus, are conserved among fastclick species, and slow CPGs are shared among slowclick species. In addition, our results suggest that testosterone plays a role in organizing fast CPGs in fastclick species, but not in slowclick species. Moreover, fast CPGs are not inherited by all species but monopolized by fastclick species. The results suggest that speciesspecific calls of the genus Xenopus have evolved by utilizing conserved slow and/or fast CPGs inherited by each species.