|
|
||||||||
RENAL HEMODYNAMICS AND CARDIORENAL INTEGRATION
1Department of Molecular Pharmacology, Physiology, and Biotechnology, Brown University, Providence, Rhode Island; 2Department of Physics, Danish Technical University, Lyngby, Denmark; 3Department of Physics, Saratov State University, Saratov, Russia; 4Department of Physiology & Biophysics, University of South Florida, Tampa, Florida; and 5Department of Medical Physiology, Panum Institut, University of Copenhagen, Copenhagen N, Denmark
Submitted 9 August 2004 ; accepted in final form 19 January 2005
| ABSTRACT |
|---|
|
|
|---|
tubuloglomerular feedback; myogenic mechanism; nonlinear interactions; amplitude modulation
Simulation results with the model predicted that TGF modulates both the amplitude and the frequency of the myogenic oscillation, a result that has not been reported in the experimental literature. Frequency modulation is a variation in time of the frequency of an oscillation, and the two forms of modulation can arise because of nonlinear interactions with another oscillator. Standard spectral techniques assume that the signal being measured is stationary, a condition violated when there is frequency or amplitude modulation. We turned to wavelet transforms to measure the two forms of modulation (3). Wavelets are mathematical functions that represent a class of single waves that are soliton-like. The process of convolving them with time series, known as wavelet transforms, can provide separate measures of amplitude and frequency as continuous functions of time. We applied wavelet transforms to previously published results of blood flow in single renal efferent arterioles (25) and to the simulation results. The analysis confirmed the presence of amplitude and frequency modulation in the experimental results; the modulating frequency was the same as the frequency of the TGF oscillation. The results of these analyses are also presented in this paper. The comparison of simulation and experimental results is consistent with coupling of TGF to voltage-gated Ca channels.
The modulation measured from experimental results is more variable than the simulation predicts. We used arterial blood pressure records obtained by telemetry from normotensive rats as an input to the model. The result indicates that fluctuations of arterial pressure contribute to the variability of frequency and amplitude modulation.
| METHODS |
|---|
|
|
|---|
In the second study, Sprague-Dawley and spontaneously hypertensive rats rats were anesthetized with pentobarbital to permit the implantation of a blood pressure transducer, radio transmitter, and battery. A needle filled with a proprietary anticoagulant gel was inserted into the abdominal aorta, and the rest of the device was placed in the abdominal cavity. The incision was closed, and the rat was allowed to recover. After a suitable postoperative period, the rat was placed in a cage with a radio receiver nearby, and the blood pressure was recorded while the rat was free to move about the cage. Food and water were provided to the rat on the same schedule as to all other rats in the colony. The signal was passed through an analog low-pass filter with a cutoff frequency of 3 Hz, which was sufficient to remove the fluctuation due to the heart rate. Details are found in Ref. 10. Only results from normotensive rats are used in this paper.
Wavelet transforms. Preliminary model results with membrane-based coupling sites gave evidence of amplitude and frequency modulation of the faster oscillation by the slower one. Frequency modulation is an example of nonlinear interaction; amplitude modulation can arise from either a linear or nonlinear interaction. The equations we derived to model the interaction are nonlinear (18). Frequency and amplitude modulation can easily be seen and measured in the model results run at constant arterial pressure, but not in model results run with variable arterial pressure (18), or in results of experimental studies (see below). A uniform method was therefore needed for the analysis of all results. Standard spectral methods can give some information on the modulation process, but they are unable to separate the two forms if they coexist. We therefore made use of wavelet analysis, a class of analytic techniques that are well designed for such problems. The method provides an opportunity to study the temporal evolution of signals with different rhythmic components.
The wavelet transform of a signal x(t) is performed by convolving a wavelet function with the signal, as follows:
![]() | (1) |
where
![]() |
is a function constructed from a "mother" wavelet
(u) using scaling and translation. The latter wavelet is a function that, in general, can have an arbitrary shape provided it is soliton-like; that is, is localized over a time interval and has a zero average in time. The Tx (a,t) are the wavelet coefficients, a is a time scale parameter, t is time, and the asterisk refers to the complex conjugate. In the analysis of various rhythmic components the so-called Morlet function is often applied. Its simplified expression appropriate to the study of low-frequency dynamics has the following form:
![]() | (2) |
The wavelet represents a harmonic oscillation with frequency f = 1/a (angular frequency 2
f) and with amplitude modified in time by a Gaussian factor that describes how the wave arises and decays. Hereafter, we consider the timefrequency representation of the wavelet transform. With this, the coefficients Tx(f,t) measure the spectral contribution near the frequency f at time t of the observed signal. In addition to the wavelet transform coefficients Tx(f,t), we estimated the energy density, Ex(f,t), from
![]() | (3) |
where a = 1/f is a time scale and C is a constant chosen dependent on the form of
. The result is a surface in three-dimensional space Ex(f,t), with well-defined ridges corresponding to different oscillatory modes in the signal. Sections of this surface at fixed moments of t = t0 correspond to the local energy spectrum. At each sampled moment of time, the maxima of the local spectra provide frequencies and amplitudes of the rhythmic components. The instantaneous frequencies and amplitudes are then analyzed for entrainment and modulation of those various rhythmic components.
There is a complication with the application of wavelet methods that is specific to the renal blood flow data. All methods suppose that a spectrum is calculated within some data window. To obtain proper results, the spectral components, such as the frequencies being estimated, should be nearly constant within this window. In the case of wavelet analysis, this means that during 56 oscillations of the fast mode, the corresponding frequency should not change significantly. In our experience, this condition is usually satisfied if the modulated frequency is about 30 times the modulating frequency, but the ratio of the two oscillations for nephron data is 5:1. In the data we wish to analyze, there will be significant violation of any assumption of constant frequency within even the smallest useful sampling window. The limited resolution that would be imposed could produce spurious peaks in the modulation spectrum.
To solve the problem created by the relatively small ratio between the two frequencies, we estimated the modulation properties of the signal by introducing the following procedure:
1. From the original signal x(t), extract the time dependencies of the instantaneous frequencies ffast and fslow related to the fast and slow oscillatory modes, respectively, using wavelet analysis. These frequencies are identified using the energy spectrum, equation (3), by taking the energy maxima in the frequency bands 0.020.05 Hz for fslow(t), and 0.10.2 Hz for ffast(t). The two frequencies are identified for each time point, ti, and a new time series is constructed from the frequencies in the energy spectrum at each of the time points in the original time series.
2. ffast(t) is used now as the input signal for the next wavelet transform, that is, we take ffast(t) instead of x(t) in Eq. (1). Again, the wavelet coefficients and the energy density are estimated. The latter will contain information about all modes involved in the modulation process. Instead of the instantaneous frequency of the fast dynamics, we can take the instantaneous amplitude of this oscillatory mode and, therefore, it is possible to study the properties of amplitude modulation of the fast rhythm as well. The frequency resolution of the method is
0.01 Hz.
Model. The model and the numerical methods used to solve it are presented in detail in the accompanying paper (18). The model consists of three partial differential equations that describe pressure, flow, and NaCl concentration as functions of time and distance along the nephron, one initial valued ordinary differential equation that describes the rate of change of protein concentration in glomerular capillaries as a function of distance, and six initial valued ordinary differential equations that describe the rates of change of a number of intracellular variables in each of two afferent arteriolar segments. The input to the arteriolar models is the arterial pressure, and the model predicts the segmental resistance to flow, from which the GFR and renal plasma flow are calculated.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
Blood pressure fluctuations and frequency modulation. The simulation results of Fig. 1 show nearly perfect regularity in the amplitude and frequency modulation, whereas the experimental results in Fig. 3 show considerable variability in these measures. One possible difference in the conditions between the two is that the simulation was conducted with constant blood pressure, while the experiments were done with the fluctuations in arterial blood pressure that are normally present (10). Therefore, we used blood pressure records made by telemetry from conscious Sprague-Dawley rats to determine whether the blood pressure fluctuations might be responsible for the changing modulation. Fig. 6 contains the results. The arterial pressure record shown in Fig. 6A was used as an input to the same model as was used in the simulation of Fig. 1. The predicted nephron plasma flow rate shows prominent oscillations at the TGF frequency and higher-frequency smaller-amplitude oscillations as well. The amplitude of the TGF oscillations varies with the arterial pressure, higher pressures generating larger oscillations. Because TGF is coupled to voltage-gated Ca channels, the variation in the amplitude of TGF oscillations sets the stage for variation in the amplitude and frequency of the myogenic oscillation. Figure 6, C and D show the instantaneous amplitude and frequency calculated by wavelet transforms from the predicted nephron plasma flow rate. As can be seen, the amplitude and frequency vary in an irregular pattern, as is the case in the experiment shown in Fig. 3. The range of the amplitude modulation in Figs. 3 and 6 cannot be compared directly because, as explained in the legend to Fig. 3, the units are different. The frequency modulation in Fig. 6 varies over a greater range of frequencies than in Fig. 3. Note that the blood pressure in the experiment whose data are in Figs. 2 and 3 was lower than in the experiment of Fig. 6. In this model (18), as in experiments with rats (25), the amplitude of the TGF oscillation varies directly with arterial blood pressure. The greater amplitude of the frequency modulation in Fig. 6, compared with Fig. 3, could be attributed, at least in part, to the difference in blood pressures. System and measurement noise is another factor that could contribute to the difference in the modulation response. The data analyzed for Fig. 3 have all sources of noise for the measurement, whereas the data of Fig. 6 has only the variation in arterial blood pressure. Figure 7 compares the instantaneous frequency of the TGF oscillation with the frequency of the signal modulating myogenic amplitude and frequency. As can be seen in Fig. 7, both the TGF oscillation and the modulating signals occupy the same frequency band, a result consistent with the suggestion that the TGF signal is responsible for the modulation phenomena.
|
|
| DISCUSSION |
|---|
|
|
|---|
, the coupling coefficient, so that the result is not critically dependent on a particular value of this new coefficient. In addition, the model predicted that TGF modulates the frequency and the amplitude of the myogenic oscillation. This result has not been reported previously. To detect this phenomenon, it is necessary to measure the frequency and the amplitude of the oscillation as continuous functions of time. The methods that permit this kind of analysis include time frequency analysis (2) and wavelet transforms (3). We selected wavelet transforms because the method permits the complete separation of amplitude and frequency information and provides quantitative information. With this technique one convolves a wavelet function with the time series to be analyzed. Appropriate selection of coefficients permits the analysis of specific frequencies. This procedure was applied to the simulation results and to measured single-nephron blood flow rate. The analysis of the simulation results shows that the amplitude varies in synchrony with the TGF oscillation. This result means that the amplitude of the myogenic oscillation varies. Close examination of the simulation's time series, Fig. 1A, shows the variation. The magnitude of the myogenic oscillation will depend on Cai. In the model, TGF is coupled to voltage-gated Ca channels. If fewer Ca channels open because of the action of TGF, Ca will not rise as much, and the myogenic oscillation will be attenuated, and vice versa.
Figure 1 also shows the frequency of the myogenic oscillation in the simulation result. The oscillation switches back and forth between frequencies of 0.096 and 0.141 Hz. The clock that triggers the myogenic oscillation comprises a set of voltage and Ca sensitive K channels, voltage-gated Ca channels, and a leak pathway. The three pathways interact and form an autonomous limit cycle, the myogenic oscillator (5, 18). A change in any of the ionic currents will change the clock rate. The TGF signal does not have a switch-like shape, but the system that generates the clock signal is nonlinear, and it is therefore to be expected that the frequency will not trace the TGF signal exactly.
We have previously used wavelet transforms with tubular pressure measurements (24), which have less noise than the blood flow; the results of that study and of the analysis of single-nephron blood flow in this one are in accord. Both the amplitude and the frequency modulation result with measured single-nephron blood flow show considerable variation with time, in contrast to the simulation result. The question at issue is to what extent variation in TGF is responsible for this variability. We developed the double-wavelet method to answer this question. The two time series of Fig. 3 were themselves subjected to wavelet transforms to determine the frequency of the modulation effect. Fig. 4 shows the result. The TGF frequency is also shown and is not constant. It does, however, remain within a fairly narrow frequency band. The frequencies of the amplitude and the frequency of the myogenic oscillation also vary within this same band. The similarity of the frequency bands is consistent with the hypothesis that the TGF signal modulates both the amplitude and the frequency of the myogenic oscillation, at least in this one experiment. Results from 10 rats were obtained in (25). We analyzed all 10 records and calculated the instantaneous ratio of the modulation frequency to the TGF frequency. This procedure generated
50,000 data points. As can be seen in Fig. 5, the results cluster around a ratio of unity, which is also consistent with the hypothesis that TGF modulates the amplitude and the frequency of the myogenic oscillation.
We speculated that the naturally occurring fluctuations in mean arterial blood pressure could be responsible for at least part of the variability in TGF and therefore in the frequency and amplitude modulation. To test this suggestion, we used measured blood pressure records from conscious rats (10) as the input to the model. The arterial pressure in this record varied over nearly 60 mmHg. The model predicts that the amplitude of the TGF oscillation varies directly with the arterial pressure but that nephron plasma flow remains stable around a mean value of 140 nl/min. The amplitude of the myogenic oscillation varied at the frequency of TGF, but the largest changes occurred at times when the blood pressure changed values. Also, please note that the ordinate of the graph showing amplitude vs. time is larger than in Fig. 1. The results of Fig. 1 were obtained with a constant arterial pressure of 100 mmHg. The pressures in Fig. 6 were largely greater than 100 mmHg. The predicted amplitude of the myogenic oscillation varies with arterial pressure (see Fig. 10 in Ref. 18) and experimentally (25). The greater amplitude in Fig. 6 probably results from the higher pressures.
This numerical experiment also showed significant frequency modulation and considerable variability in that modulation. Recall that this is the same model that generated Fig. 1, with its regular switching of both amplitude and frequency. When a measured arterial pressure record is used as input to the model, both forms of modulation become more variable. We suggest therefore, that naturally occurring fluctuations in arterial pressure provide at least part of the cause of the variation in frequency and amplitude modulation.
We used voltage-gated Ca channels as the point in the cell that couples to TGF. We based this choice on the strict dependence of the myogenic mechanism on extracellular Ca, and the presence of voltage-gated Ca channels in the plasma membranes of the arteriolar cells. The behavior of the model matches experimental data on autoregulation and on system dynamics, as detailed above. There are other coupling points that have not been evaluated in this study but that may play an important role. There is experimental evidence to support the operation of these other mechanisms and no reason a priori to assume that the entire coupling burden is carried by a single mechanism. These include signals mediated through intracellular Ca stores; receptor-operated Ca channels, which also mediate actions on intracellular Ca stores (13, 19); and Ca-insensitive regulation, most likely an action on phosphorylation of myosin light chain (11). The phosphorylation mechanism is unlikely to produce significant frequency modulation because the action of this signaling pathway does not affect the constituents of the myogenic clock mechanism, but it could act in concert with other mechanisms. Store-operated Ca release is a very complex and important mechanism, but we are unable intuitively to predict its effect on system dynamics.
Nonlinear interactions between TGF and the myogenic mechanism are probably important signaling mechanisms in renal autoregulation. The two different forms of modulation we have presented, both of which arise from nonlinear interactions, solve two different problems. The first is time coordination and stabilization. The arterial pressure can change rapidly. Propagation of a pressure disturbance through the tubular system to the macula densa entails a delay of 1015 s (18). The myogenic mechanism can respond to the same signal in a fraction of this time. Unless there is some way to coordinate the two mechanisms, they could act in an uncoordinated manner and fail to provide adequate autoregulation. Amplitude modulation should permit the actions of the two mechanisms to act cooperatively. In human-made control systems, the introduction of a fast, low-ampltidude signal inside the feedback loop of a nonlinear controller has the effect of narrowing the nonlinear sector and stabilizing the system (26, 27). Because TGF is a nonlinear system, its interaction with the myogenic mechanism improves the stability of autoregulation. When the blood pressure increases, the kidney needs the afferent arteriolar resistance to increase, and vice versa. Without some form of coordination between the two mechanisms, it is to be expected that the system would fail to respond as the kidney needs it to and that the blood flow will fluctuate inappropriately. Amplitude modulation helps to prevent this.
Frequency modulation solves a second problem, which is the fidelity of signal transmission. Frequency modulation is well known in signal processing as a robust method for sending and receiving signals in a noisy environment. Arteriolar smooth muscle contraction is triggered by increasing Cai, but Ca has many more functions in these cells than controlling muscle contraction. It is no exaggeration to say that the Ca environment of these cells is noisy. Spatial inhomogeneity is one answer to this problem, and recent evidence points to localization of calmodulin with voltage-gated Ca channels, which are inhomogeneously distributed in the plasma membrane (21). Frequency modulation is another, complementary way to connect the signal and the target.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. Iliescu, R. Cazan, G. R. McLemore Jr., M. Venegas-Pont, and M. J. Ryan Renal blood flow and dynamic autoregulation in conscious mice Am J Physiol Renal Physiol, September 1, 2008; 295(3): F734 - F740. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Kleinstreuer, T. David, M. J. Plank, and Z. Endre Dynamic myogenic autoregulation in the rat kidney: a whole-organ model Am J Physiol Renal Physiol, June 1, 2008; 294(6): F1453 - F1464. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. V. Sosnovtseva, A. N. Pavlov, E. Mosekilde, K.-P. Yip, N.-H. Holstein-Rathlou, and D. J. Marsh Synchronization among mechanisms of renal autoregulation is reduced in hypertensive rats Am J Physiol Renal Physiol, November 1, 2007; 293(5): F1545 - F1555. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Balasubramanian, A. Ahmed, C.-M. Lo, J. S. K. Sham, and K.-P. Yip Integrin-mediated mechanotransduction in renal vascular smooth muscle cells: activation of calcium sparks Am J Physiol Regulatory Integrative Comp Physiol, October 1, 2007; 293(4): R1586 - R1594. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang, R. D. Loutzenhiser, and W. A. Cupples Frequency modulation of renal myogenic autoregulation by perfusion pressure Am J Physiol Regulatory Integrative Comp Physiol, September 1, 2007; 293(3): R1199 - R1204. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. N. Pavlov, V. A. Makarov, E. Mosekilde, and O. V. Sosnovtseva Application of wavelet-based tools to study the dynamics of biological processes Brief Bioinform, December 1, 2006; 7(4): 375 - 389. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. T. Layton, L. C. Moore, and H. E. Layton Multistability in tubuloglomerular feedback and spectral complexity in spontaneously hypertensive rats Am J Physiol Renal Physiol, July 1, 2006; 291(1): F79 - F97. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Raghavan, X. Chen, K.-P. Yip, D. J. Marsh, and K. H. Chon Interactions between TGF-dependent and myogenic oscillations in tubular pressure and whole kidney blood flow in both SDR and SHR Am J Physiol Renal Physiol, March 1, 2006; 290(3): F720 - F732. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |