Search PubMedSearch

Biomedical subjects

W R Holmes

Publications and source records attributed to W R Holmes.

10 recordsLinked to original sources

Interpretation of time constant and electrotonic length estimates in multicylinder or branched neuronal structures.

1. We have investigated the theoretical and practical problems associated with the interpretation of time constants and the estimation of electrotonic length with equivalent cylinder formulas for neurons best represented as multiple cylinders or branched structures. Two analytic methods were used to compute the time constants and coefficients of passive voltage transients (and time constants of current transients under voltage clamp). One method, suitable for simple geometries, involves analytic solutions to boundary value problems. The other, suitable for neurons of any geometric complexity, is an algebraic approach based on compartmental models. Neither of these methods requires the simulation of transients. 2. We computed the time constants and coefficients of voltage transients for several hypothetical neurons and also for a spinal motoneuron whose morphology was characterized from serial reconstructions. These time constants and coefficients were used to generate voltage transients. Then exponential peeling, nonlinear regression, and transform methods were applied to these transients to test how well these procedures estimate the underlying time constants and coefficients. 3. For a serially reconstructed motoneuron with 732 compartments, we found that the theoretical and peeled tau 0 values were nearly equal, but the theoretical tau 1 was much larger than the peeled tau 1. The theoretical tau 1 could not be peeled because it was associated with a coefficient, C1, that had a very small value. In fact, there were 156 time constants between 1.0 and 6.0 ms, most of which had very small coefficients; none had a coefficient larger than 2% of the signal. The peeled value of tau 1 (called tau 1 peel) can be viewed as some sort of a weighted average of the time constants having the largest coefficients. 4. We studied simple hypothetical neurons to determine what interpretation could be applied to the multitude of theoretical time constants. We found that after tau 0, there was a group of time constants associated with eigenfunctions that were odd (or approximately odd) functions with respect to the soma. These time constants could be interpreted as "equalizing" time constants along particular paths between different pairs of dendritic terminals in the neuron. After this group of time constants, there was one that we call tau even because it was associated with an eigen-function that was approximately even with respect to the soma. This tau even could be interpreted as an equalizing time constant for charge equalization between proximal membrane (soma and proximal dendrites) and distal membrane (including all distal dendrites).4=

Animals

Electrotonic length estimates in neurons with dendritic tapering or somatic shunt.

1. Compartmental models were used to compute the time constants and coefficients of voltage and current transients in hypothetical neurons having tapering dendrites or soma shunt and in a serially reconstructed motoneuron with soma shunt. These time constants and coefficients were used in equivalent cylinder formulas for the electrotonic length, L, of a cell to assess the magnitude of the errors that result when the equivalent cylinder formulas are applied to neurons with dendritic tapering or soma shunt. 2. Of all the formulas for a cylinder (with sealed ends), the most commonly used formula, which we call L tau 0/tau 1 (the formula uses the current-clamp time constants tau 0 and tau 1), was the most robust estimator of L in structures that tapered linearly. When the diameter at the end of the cylinder was no less than 20% of the initial diameter, L tau 0/tau 1 underestimated the actual L by at most 10%. 3. The equivalent cylinder formulas for a cylinder were applied to neurons modeled as a cylinder with a shunted soma at one end. The formula for L based solely on voltage-clamp time constants gave an exact estimate of L. However, the second voltage-clamp time constant cannot be reliably obtained experimentally for neurons studied thus far. Of the remaining formulas, L tau 0/tau 1 was again the most robust estimator of L. This formula overestimated L with the size of the overestimates depending on beta, rho beta = 1, and the actual L of the cylinder, where beta is the soma shunt factor, and rho beta = 1 is the dendritic-to-somatic conductance ratio when beta = 1 (no shunt). When the actual L was 0.5 and the soma shunt was large, this formula overestimated L by two- to threefold, but when the actual L was 1.5, the overestimate was only 10-15% regardless of the size of the shunt. 4. In neurons modeled as two cylinders with soma shunt, the L tau 0/tau 1 value computed with the actual tau 0 and tau 1 values overestimated the average L by two to six times when soma shunt was large. However, the L tau 0/tau 1 estimates computed with tau 0 and tau 1 values estimated with the exponential fitting program DISCRETE from voltage transients computed for these neuron models were never this large because of inherent problems in estimating closely spaced time constants from data.(ABSTRACT TRUNCATED AT 400 WORDS)

Animals

Estimating the electrotonic structure of neurons with compartmental models.

1. A procedure based on compartmental modeling called the "constrained inverse computation" was developed for estimating the electrotonic structure of neurons. With the constrained inverse computation, a set of N electrotonic parameters are estimated iteratively with use of a Newton-Raphson algorithm given values of N parameters that can be measured or estimated from experimental data. 2. The constrained inverse computation is illustrated by several applications to the basic example of a neuron represented as one cylinder coupled to a soma. The number of unknown parameters estimated was different (ranging from 2 to 6) when different sets of constraints were chosen. The unknowns were chosen from the following: dendritic membrane resistivity Rmd, soma membrane resistivity Rms, intracellular resistivity Ri, membrane capacity Cm, dendritic membrane area AD, soma membrane area As, electrotonic length L, and resistivity-free length, rfl (rfl = 2l/d1/2 where l and d are length and diameter of the cylinder). The values of the unknown parameters were estimated from the values of an equal number of known parameters, which were chosen from the following: the time constants and coefficients of a voltage transient tau 0, tau 1, ..., C0, C1, ..., voltage-clamp time constants tau vc1, tau vc2, ..., and input resistance RN. Note that initially, morphological data were treated as unknown, rather than known. 3. When complete morphology was not known, parameters from voltage and current transients, combined with the input resistance were not sufficient to completely specify the electrotonic structure of the neuron. For a neuron represented as a cylinder coupled to a soma, there were an infinite number of combinations of Rmd, Rms, Ri, Cm, AS, AD, and L that could be fitted to the same voltage and current transients and input resistance. 4. One reason for the nonuniqueness when complete morphology was not specified is that the Ri estimate is intrinsically bound to the morphology. Ri enters the inverse computation only in the calculation of the electrotonic length of a compartment. The electrotonic length of a compartment is l[4 Ri/(dRmd)]1/2, where l and d are the length and diameter of the compartment. Without complete morphology, the inverse computation cannot distinguish between a change in d or l and a change in Ri. Even when morphology is known, the accuracy of the Ri estimate obtained by any fitting procedure is affected by systematic errors in length and diameter measurements (i.e., tissue shrinkage); the Ri estimate is inversely proportional to the length measurement and proportional to the square root of the diameter measurement.(ABSTRACT TRUNCATED AT 400 WORDS)

Animals