Membrane current series monitoring: essential reduction of data points to finite number of stable parameters

In traditional studies of changes in cell membrane potential or trans-membrane currents a large part of the recorded data presents “a pure noise.” This noise results mainly from the random openings of membrane ionic channels. Different types of stationary or non-stationary noise analysis have been used in electrophysiological experiments for identification of channels kinetic states. But these methods have a limited power and often cannot answer to the main question of the experimental study: do external factors induce a significant change of channels kinetics? A new method suggested in the current study is based on the scaling properties of the beta-distribution function that allows reducing the series containing 200,000 and more data points to analysis of only 10–20 stable parameters. The following clusterization using the generalized Pearson correlation function allows taking into account the influence of an external factor and combine/separate different parameters of interest into a statistical cluster considering the influential parameter. This method which we call BRC (Beta distribution-Reduction-Clusterization) opens new possibilities in creation of a largely reduced database while extracting specific fingerprints of the long-term series. The BRC method was validated using patch clamp current recordings containing 250,000 data points obtained from the living cells and from open tip electrode. The numerical distinction between these two series in terms of the reduced parameters was obtained.

In traditional studies of changes in cell membrane potential or trans-membrane currents a large part of the recorded data presents "a pure noise." This noise results mainly from the random openings of membrane ionic channels. Different types of stationary or non-stationary noise analysis have been used in electrophysiological experiments for identification of channels kinetic states. But these methods have a limited power and often cannot answer to the main question of the experimental study: do external factors induce a significant change of channels kinetics? A new method suggested in the current study is based on the scaling properties of the beta-distribution function that allows reducing the series containing 200,000 and more data points to analysis of only 10-20 stable parameters. The following clusterization using the generalized Pearson correlation function allows taking into account the influence of an external factor and combine/separate different parameters of interest into a statistical cluster considering the influential parameter. This method which we call BRC (Beta distribution-Reduction-Clusterization) opens new possibilities in creation of a largely reduced database while extracting specific fingerprints of the long-term series. The BRC method was validated using patch clamp current recordings containing 250,000 data points obtained from the living cells and from open tip electrode. The numerical distinction between these two series in terms of the reduced parameters was obtained.

INTRODUCTION
During electrophysiological studies it is common to record rather long tracks of signals. These signals are registered as temporal variations of cell membrane potential or trans-membrane currents induced by the opening of some ligand-or voltage-gated or even chaotic ionic channels. Usually the principal aim of such a study is the registration of some macroscopic signals-evoked or spontaneous-and the change of parameters of these signals characterizes the total effect of some actions that are located in the experimental object. But a large part of the record forms a so-called "empty track" containing a "pure noise" only. It is well known that this noise reflects mainly the result of random openings of transmembrane ionic channels. Different types of stationary or non-stationary noise analysis have been used for identification of these channels' states (Neher and Sakmann, 1976;Sigworth, 1980Sigworth, , 1985Sigworth, , 1986Läuger, 1985;Traynelisa and Jaramilloa, 1998;Alvarez et al., 2002;Venkataramanan and Sigworth, 2002).
Unfortunately, these methods have not come into widespread use among physiologists since they often cannot answer the main question of the study: If this drug or this change of environment state induces the reliable change of channels condition or not?
Thus, there is an urgent task to develop a special language that can be compact and reliable in order to describe accurately very long current streams (long-time series) with hidden signals and noise in terms of a finite and statistically understandable set of reduced parameters. In this paper we want to show how to develop this special language based on an example of the analysis of signals recorded in rat's spinal cord slices. Besides this problem we want to show how to detect the presence of the biological object inside the experimental set. For this purpose we also recorded data representing the dependence of the current vs. time when the biological object is absent. Examples of currents recorded in a living cell and with empty electrodes are shown in Figure 1. It is well noticeable that these two signals are apparently very similar. Even though generally distinguishable by an experienced observer the reliability of these differences cannot be numerically evaluated without some special analytic methods.
To the authors' best knowledge one method is basically suitable for quantitative analysis of the different long-time series. This method was introduced by Peng et al. (1994) and nowadays it is known as detrended fluctuation analysis (DFA). It was well described in literature by their creators (Ossadnik et al., 1994;Peng et al., 1995) and found its application in analysis of biomedical (Penzel et al., 2003;Jospin et al., 2007;Burr et al., 2008) and other (Hausdorff et al., 1995(Hausdorff et al., , 1996 data. But it is necessary to note that the DFA algorithm works well only for certain types of non-stationary time series (especially having slowly varying trends), it is not designed to handle all possible non-stationarities in real-world data. This algorithm was not free also from uncontrollable errors that are associated with approximate fitting of detrended fluctuations by the segments of straight lines or by the parabolic or high order polynomials (Kantelhardt et al., 2001). The final straight line with power-law exponent α DFA is obtained as a slope in a double-log scale as a result of the fitting procedure and contains the fitting error that depends also on the type of segmentation of the initial series considered. These uncontrollable errors (usually they are not properly analyzed in the literature) can lead to different results in calculation of the desired value of the α DFA and other associated fitting parameters in analysis of the same long-time series.
A technique, called scale-dependent Lyapunov exponent (SDLE, see Gao et al., 2006Gao et al., , 2012bGao et al., , 2013Hu et al., 2010), provides a more comprehensive characterization of complex time series. Some of DFA's limitations have been overcome recently as well by using a new method called adaptive fractal analysis (AFA, see Gao et al., 2010Gao et al., , 2011Gao et al., , 2012aRiley et al., 2012;Kuznetsov et al., 2013). AFA has been shown to be able to determine global trends, remove noise, perform fractal analysis and multiscale decomposition and present data as a curve. However, new tools could be developed specifically designed to show and estimate even mild differences between two long time series.
Thus, it would be desirable to have a new method with "high resolution" (10-20 significant parameters) to distinguish more accurately the experimental data and effect of treatments. In this paper we demonstrate such method based on some invariant properties of the beta-distribution function; furthermore this method admits a procedure that controls the error in each stage of its application. From our point of view the effectiveness of new approach is based on the monotone behavior of the primary fitting parameters that admit the secondary fit. This peculiarity allows compressing initial fitting parameters with the help of the secondary fit and present initial data set in more compact form.
The four fitting parameters (A, B, α, β) of beta-distribution can be interpreted and used for quantitative reading of fluctuations arising on different scales of the long-time series considered. In previous papers (Nigmatullin, 2010;Nigmatullin et al., 2012) based on the principle of the strong correlation of random sequences it was shown that the cumulative (integral) curve obtained from the sequence of the ranged amplitudes (SRA) can be described with high accuracy by means of the beta-distribution function. In other words, any detrended random sequence being transformed to the SRA (when all amplitudes of the initial sequence are sorted out and located in the descending order y 1 > y 2 > . . . > y N ) after elimination of its mean value and subsequent integration, forms a bell-like curve J(x) that can be fit (with controllable relative error) by the function: (1) Here the limiting values x 0 < x N define the ends of the location interval of the random sequence considered. In many cases the parameters x 0 , x N are known. Other quantitative parameters (A, B, α, β) should be found from the fitting procedure of the function J(x) to the curve Jb(x). The power-law exponents (α, β) reflect the fractal properties of the random sequence considered and the presence of the memory that is expressed in the behavior of the corresponding SRAs. The criterion for the verification of the presence of memory in two random sequences which are compared is as follows. If one SRA being plotted with respect to another one forms a curve close to a straight line then these two random curves are defined as a having a relative memory and can be considered as being strongly correlated. This important property allows transforming any segment of a random sequence to a beta-distribution function and "read" this segment in terms of four unknown fitting parameters (A, B, α, β). Such transformation from 30 to 50 or even more initial points belonging to a random sequence can be read in terms of these four parameters only. This allows us to suggest a new type of spectroscopy based on some scaling properties of the beta-distribution. This transformation is called Fluctuation Spectroscopy based on Beta-Distribution (FSBD). In general we suggest a method which we call BRC (Beta distribution-Reduction-Clusterization). The basic problem that is solved in this paper by using the BRC method can be formulated as follows: Is it possible to suggest a reliable method with controllable error that has a wide range of applicability and which has a flexible small set (10-20) of statistically understandable parameters for quantitative characterization of the differences between long-time series?
Interneurons were voltage-clamped at −65 mV in the wholecell configuration after obtaining GV seals (usually not less than 2 GV) by means of a patch-clamp amplifier (Axopatch 200B; Molecular Devices, Sunnyvale, CA, USA). Compensation of capacitance (Cm) and series resistance (Rs) was achieved with the inbuilt circuitry of the amplifier. Series resistance was compensated by 40 ÷ 70% and did not change appreciably from the beginning to the end of the experiments, indicating stable recording conditions. The tracks used for comparison were recorded by the immersion of filled patch-pipettes in artificial cerebrospinal fluid; the patch-pipettes were voltage-clamped at −65 mV too.
Then all data were sampled at 10 kHz and stored on-line with a PC using the pClamp 10.0/Clampex 10.0 software package (Molecular Devices).

SCALING PROPERTIES OF THE BETA-DISTRIBUTION AND DESCRIPTION OF THE TREATMENT PROCEDURE
In this section we want to demonstrate the scaling properties of Expression (1). We subject x, x 0 and x N in Expression (1) to the following scaling transformations, keeping the power-law exponents α and β invariable: which gives the following beta transformation: where A = A · ξ (α + β) . This is the accurate mathematical result that follows from the scaling transformation of the initial coordinates.
In order to have a simple criterion for comparison of the two beta-distributions let us calculate the values of two extreme points x,x belonging to the functions Jb(x) and Jb(x ) respectively.
From Expressions (3) it follows that for the scaling transformation (2) the heightsH,H of the extreme points of the two belllike distributions at the fixed values of the power-law exponents α and β and parameter B should coincide with each other. Besides this criterion it is necessary to take into account the scaling relationship between the heightsH,H . If two power-law exponents α and β are subjected to the scaling transformation at the fixed value of the length = x N − x 0 : then simple manipulations lead to the second scaling relationship: Here the amplitudes A and A are defined by relationships (1) and (2), respectively. The consideration of the scaling properties of the beta-distribution allows one to suggest the following two steps.
Step 1. This step includes the formation of the sequence of the range amplitudes (SRA) when all amplitudes located on the fixed Step 2. Numerical integration of the SRA with respect to its mean value and subsequent fit to the function (1). Figure 2 illustrates this transformation which is realized after application of these two steps.
Each sub-segment having equal length is transformed to its SRA (Figure 2A) in Step 1, and the integration of the SRAs with respect to its subtracted mean value gives finally the desired bell-like curve that can be fit to Expression (1) in Step 2. Mathematically these two steps correspondingly are expressed as: Here the integer index j (j = 1, 2, . . . , N) numerates the number of data points in the fixed segment = x N -x 0 containing initially 30-50 data points. trapezoid method. As a result of calculation of Expression (6) we obtain the desired bell-like curve J(x j ). Figure 2B shows the quality of the fitting of the bell-like curve obtained to the beta-distribution. In order to have the value of the relative error: to be limited to a few percentages (2-5)% we should choose the length of the minimal segment min of the initial series containing initially 30-50 data points. In Expression (7) the value N defines the number of data points that enters in the segment of the length . Thus, the first reduction criterion should be written as: Here the scaling parameter ξ has the same meaning as in Expression (2). This requirement allows one to consider the long-time series containing the total number of data points (j = 1, 2, . . . , N total ) in terms of the reduced parameters of the beta-distribution (A, B, α, β) depending on parameter k. Further it is convenient to rewrite condition (8) in the following form changing the numeration of the current parameter k: where [in comparison with (8)] the value 1 should coincide with the minimal value 30 < min < 50 giving the condition for finding the limiting value of K (the total number of segments is equaled to K + 1). In the opposite case, the value K + 1 should give the maximal length coinciding with the value N total . As a result of this reduction procedure one can transform N total data points to 4.(K + 1) parameters. But this step is not sufficient. If the functions A k , B k , (α + β) k have monotonic behavior one can realize further reduction to the primary set of the fitting parameters describing these functions. Now it is necessary to explain why the sum of the parameters (α + β) is selected instead of considering each-power law exponent separately. This selection is based on the comparison of these exponents with the single power-law exponent α DFA figuring as the basic parameter in the DFA. It is easy to see that relationship α + β = 1 with α ≈ β ≈ 0.5 (for this case betadistribution looks like a semicircle) corresponds to a distribution with the absence of power-law correlations in the time series. From another side it gives for α DFA = 0.5. Comparison with these two power-law exponents leads us to the following approximate expression: One can notice also that Expression (10) does not contradict other well-known power-law exponents (Hausdorff et al., 1995;Burr et al., 2008) We want to stress here that this requirement is approximate and can serve as an indication for division of long-time series with fractal structure (because it does not contradict with wellknown inequalities) known before from series with self-similar structure. The left-hand inequality follows from the requirement β f > 0 and does not contradict with numerical results obtained in other papers (Penzel et al., 2003;Jospin et al., 2007;Burr et al., 2008). We should also note that the equality (α + β) = 2 corresponds to a uniform amplitude distribution. The uniform distribution leads to the degeneration of the corresponding SRA to a straight line (Nigmatullin, 2010). The beta-distribution in this case is described by a parabolic curve. If one of the power-law exponent (say α → 0) then the position of extreme pointx → x 0 . Because of normalization w 1 + w 2 = 1 β → 1. This statement is valid also in the opposite case when α → 1, β → 0. So, the last relationship (11) can be considered as a specific fractal test in our further calculations. Here we should also note that in practical applications the existence of the interval 0 < α + β < 1 and inequality α + β > 2 also are possible. For the first case, for small values of α and β the beta-distribution degenerates to a rectanglelike curve. In the second case the values of the derivatives on the ends (x 0 , x N ) of the beta-distribution have zero values. These two cases correspond to degeneration of the fractal properties of the time-series analyzed. The verification of relationship (11) on the Weierstrass-Mandelbrot function that represents itself the selfaffine function (see its definition in Feder, 1988) confirms the relationship (11). So, for practical purposes it is useful to work with the combination of (α + β).
The statistical and geometrical meaning of other parameters entering to (1) can be explained as follows. The value of the amplitude A together with the height H of the beta-distribution is associated with intensity of the fluctuations analyzed. As one can see from Figure 3A the angle of the SRA slope counted off from zero point (after elimination of its mean value) is proportional to the height of the corresponding fluctuation that is expressed in the form of a beta-distribution in Figure 3B. If this angle approaches the vertical axis, the height of the distribution becomes large. In the opposite case when this angle tends to zero the height of the distribution is small. See Figure 3B where the first 14 beta-distributions are shown. The measure of asymmetry can be connected with parameters B and the values of weight factors w 1,2 that are defined by Expression (3). The value w 1 = 0.5 corresponds to the complete symmetry of the distribution in the horizontal direction. Any shift of this parameter to the left-(w 1 < 0.5) or to the right-hand side (w 1 > 0.5) reflects the horizontal asymmetry of the distribution. A small asymmetry of this distribution in vertical direction is controlled by the parameter B.
Step 3. After selection of the scaling parameter ξ and the limiting value K from Expression (9) one can obtain a family of bell-like curves that can be fitted to Expression (1). The calculated fitting parameters A k , α k , β k , B k , k = 1, 2, . . . , K + 1 from Expression (1) are obtained. The set of these bell-like curves and the corresponding fitting parameters forms the total fluctuation spectrum based on the beta-distribution (FSBD). Each part of this FSBD contains the corresponding beta-distribution: Step 4. In order to subject them to the scale-invariant properties described above it is necessary to average this family of distributions and consider only one weighted distribution: located in the given interval k . Here the parameter NBd k coincides with number of beta-distributions calculated for the given k. Figure 4 shows the averaged beta-distribution obtained for the cell number 3. If N total = 250,000 then from condition (9) at the given 1 = 32 and ξ = 2 we obtain that K = 13. So, the total number of beta-distributions NBd 1 = N total / 1 = 8333. The first 14 distributions belonging to this family is shown in Figure 3B.
Step 5. Further calculations are reduced to the analysis of the functional dependencies A k , α k , β k , B k , k = 1, 2, . . . , K + 1 with respect to the variable k. We define them as the primary fitting parameters characterizing the averaged distribution (13). Further analysis shows that the amplitude A k has monotonic behavior and can be described by a simple exponential behavior: Preliminary calculations show that this monotonic behavior is conserved for the long-time series without any trend. The presence of trend distorts this behavior. This dependence follows after substitution of Expression (9) in relationship (2) for the amplitudes. The perfect fit of this monotone curve is shown in Figure 5A. Other dependencies are not so simple but nevertheless they can be identified from simple power-law and exponential hypothesis with the help of the eigencoordinates (ECs) method (Baleanu et al., 2011;Ciurea et al., 2011). The dependences <(α + β) k >≡ S k (αβ) and <B k > have also monotonic character and can be fitted by means of two simple functions: These functions are shown, respectively, in Figures 5B,C. So, finally we obtain 10 fitting parameters that can be combined with 9 parameters figuring in Expressions (14) and (15) K+1 . The behavior of this weight factor is shown in Figure 5D. These ten parameters can be used as the primary set of the fitting parameters for creation of a specific "fingerprint" of the long-time series considered. The idea of clusterization of these parameters is discussed in Results Section. Further analysis shows that the distribution of the heights and mean values of the SRAs obtained for the family of distributions at 1 also forms two other different beta-distributions. These distributions are important also for clusterization purposes because initially the information about the secondary distribution of the heights of the initially formed beta-distributions family and mean values of the corresponding SRA were not taken into account. The distributions of the heights and mean values together with their beta-distributions are shown in Figures 6, 7, correspondingly. After fitting of these two distributions one can obtain in addition 5 significant parameters characterizing each beta-distribution separately.
A mn , (α + β) mn , w 1,mn , max (Bd mn ), mean(SRA mn ) . (16) These ten additional parameters we define as the secondary fitting parameters. The statistical meaning of these parameters are the following. The parameters A H, mn characterize the amplitudes of beta-distributions referring, correspondingly, to the heights (H) and mean values (mn). The sum (α + β) H, mn contains the information about their power-law exponents, w 1, H, mn gives the information about their asymmetry, max(Bd H , Bd mn ) signifies their heights, and the fifth parameter SRA H, mn contains information about the mean values of these two additional distributions.
From our point of view, these 20 (10 primary and 10 secondary) significant parameters [figuring in Expressions (14)-(16)] combined together can completely characterize the behavior of fluctuations associated with the long-time series analyzed and containing N total = 2.5.10 5 ÷ 10 6 and even more data points.

CLUSTERIZATION OF FINAL PARAMETERS BASED ON THE GENERALIZED PEARSON CORRELATION FUNCTION
For clusterization purposes one can suggest more accurate selection of similar sequences based on internal correlations. For this aim we introduce the generalized Pearson correlation function (GPCF) (Nigmatullin, 2010;Nigmatullin et al., 2012).
where expression: determines the generalized mean value (GMV)-function of the K-th order. Here the generalized mean value (GMV) function determines the mean value for all range of the moments (see Expression (19) below). The set of parameters (s 1 ,s 2 ,. . .,s K ) determines the type of the random sequence compared. The GPCF p determined by Expression (17)  , (B) nrm j (y) = y j max ( y j ) , y j = y j − min (y j ). (20b) j = 1, 2, . . . , N, 0 < nrm(y) < 1.
Here, as it was done above, the set y j defines an initial random sequence that can contain a trend or can be compared with another trendless sequence. The symbol | . . . | and index j (j = 1.2,. . .,N) determine the absolute value and number of the measured points, correspondingly. The second case (B) in [20(b)] corresponds to the case when the initial sequence is positive. If the limits mn and mx in (20) have opposite signs and accept sufficiently large values, then the GPCF function has two plateaus (equaled unit at small numbers of mn (i.e., GPCF mn = 1) and another limiting value GPCF mx depends on the degree of internal correlation between two random sequences compared. This right-hand limit (defined as Lm) is located between two values: The appearance of two plateaus implies that all information about possible correlations is complete and further increasing of the limiting numbers (mx, mn) figuring in (19) is useless. Numerous   figure) can be used as the statistically significant parameters for characterizing of the long-time series considered. So, in the results of this complete analysis one can obtain 20 statistically significant parameters that can be used for the detailed classification of the long-time series containing 2.5.10 5 ÷ 10 6 data points.
tests showed that the high degree of correlations between two random sequences is achieved when Lm = 1, while the lowest correlations are observed when Lm = M. This empirical observation, having a general character for all random sequences, allows us to introduce new correlation parameter CC (complete correlation)-factor, which is determined as: We would like to stress here that this factor is determined on the total set of the fractional moments located between exp(mn) and exp(mx). As it was mentioned above, in practical calculations for many cases it is sufficient to put mn = −15 and mx = +15. The CC-factor accepts the unit values when the degree of correlation is high while the case Lm = M corresponds to the lowest (remnant) degree of correlations that can be observed between the compared random sequences. In addition, we want to stress also the following fact. This CC-factor does not depend on the amplitudes of the random sequences. The pair random sequences compared should be normalized to the interval: 0 ≤ y j ≤ 1. It reflects the internal structure of correlations of the compared random sequences based presumably on the similarity of their probability distribution functions that are not known in many cases. Recent example related to application of the statistics of the fractional moments was considered in paper (Nigmatullin et al., 2012). So, the CC-factor (22) can be used for clusterization of the significant parameters based on the following idea. For a set of significant parameters referring to one qualitative factor one can calculate the limits of CC-factor: Here the low correlation limit cf min is determined by the sampling volume and conditions of experiment that should be almost the same for two qualitative factors compared (control/influence of another qualitative factor).

PROCESSING OF THE LONG-TIME MEMBRANE CURRENT SERIES
In previous Section we described in details (S1-S5) basic steps of treatment of an arbitrary long-time series. Here we want to make some general remarks related to this procedure. If the long-time series considered contains the clearly expressed but random trend then its random behavior can disturb the monotonic behavior of the primary 9 parameters figuring in the fitting functions (15) and (16). In this cases one can recommend to apply the POLS (procedure of the optimal linear smoothing) described in papers (Baleanu et al., 2011;Ciurea et al., 2011;Nigmatullin et al., 2012) or simple numeric differentiation. These two procedures help to suppress the hidden random trend and obtain the monotonic behavior for the 9 parameters figuring in (15) and (16). In the shown figures we used the scaling factor ξ = 2. For the rational values of ξ from the interval (1, 2) Expression (9) can be modified as: So, numerical calculations realized at ξ = 1.5 show that results are not changed essentially, only the integer variable k in Expressions (15) and (16) is replaced as k → μk. We think that this method has a wide range of its applicability and these two modifications can be taken into account in order to express the long-range time series in terms of 20 significant parameters. In similar manner as it was treated the membrane currents for the randomly taken interneuron-3 one can treat other long-time series related to other (1, 2, 4, 5, 6, 7) interneurons. Besides, in order to differentiate these random sequences recorded without presence of a biological object we treated in the same manner 6 random sequences corresponding to empty electrode. The next problem is associated with the finding of criterion of clusterization that helps to combine these "control" membrane currents to one strongly-correlated cluster based on the values of the significant parameters. For each cell these parameters are collected in Table 1. For 6 files corresponding to pure solute (without presence of the cell) the results are collected in Table 2. How to differentiate these 20 quantitative parameters (in this case a qualitative factor is associated with the presence/absence of a biological cell) from each other? The simplest classification can be related to calculation of the mean value and standard deviation of the calculated significant parameter in each row. But more effective scheme of clusterization based on the statistics of the fractional moments and the usage of the complete correlation factor is considered in the next section.
For the clusterization of final parameters we have used new correlation parameter CC described in "Materials and Methods" section Expression (22). The calculation of the CC-factor (in our case it is based on a set of membrane currents associated with 3 "control" measurements for each chosen cell from the total set of currents representing other 7 biologic cells) which is considered as the complex correlation matrix (see Table 3) having minimal dimension (7 × 7) leads to the minimal value cf min = 0.9238. The result is not changed essentially if one calculates numerically the corresponding integrals with respect to their normalized significant parameters and then considers their CCfactors. The tendency of the strong correlations between columns of Table 1 is conserved, only the boundary of the correlation interval is slightly increased achieving the value Jcf min = 0.9736. So, using the method of clusterization based on the statistics of the fractional moments and Expression (22) one can say that all "control" currents measured for the sampling 7 × 7 = 49 form the strongly-correlated cluster with limits [0.9238, 1] for the initial set of significant parameters (20 parameters for each sampling) and [0.9736, 1] (for the corresponding integrals that are obtained by direct trapezoid method from the normalized significant parameters). In accordance with this method of clusterization one can make the following conclusion: if any another series having 20 significant parameters will give the CC-factor located in the interval [0.9238, 1] then it can be considered as the "friend" file belonging to this cluster, in the opposite case it can be considered as a "strange" file. For more reliable identification the saying above can be referred to the integrated columns formed from 20 normalized significant parameters. In the same manner we treated the files corresponding to the electrode currents recorded in normal saline solution without presence of biological object. The 20 desired parameters for 6 files are collected in Table 2. Their correlation matrix presented by Table 4 form another cluster. But attempt to combine the currents corresponding to the living o cells with currents corresponding to empty electrodes located in saline solution is unsuccessful. If we compare the correlation matrix of Table 5 with the previous ones (Tables 3, 4) then one can notice that the last matrix is uncorrelated (all elements are close to zero). It means that the presence of the biologic cell completely changes the statistical structure of the current and from qualitative point of view the long-time random sequences of currents recorded for both cases (presence/absence of biological cell) are different.
So, new clusterization method helps to express quantitatively the internal factor as the presence/absence of the living cell (compare this statement with series shown on Figure 1 where the corresponding currents look similar to each other). Definitely, more accurate measurements are needed in order to differentiate from many mixed factors that form a time-series for biological and non-biological objects a specific predominant factor that plays an essential role in this differentiation. But this problem merits a separate research.

DISCUSSION
It is well known that cellular membrane is the element which largely provides cell functioning. Cell membrane has so many functions that it is difficult even to list-anyone can find them all in each textbook on cell biology. In general the membrane provides all interaction of the cell with the external environment including the perception of the effect of active substances. Withal Frontiers in Computational Neuroscience www.frontiersin.org September 2014 | Volume 8 | Article 120 | 9     the membrane comprises a lot of elements which produce socalled "membrane noise"-rather small variations of membrane potential or trans-membrane current; mainly they are different types of ion channels, transporters and pumps. There are many active substances affecting the operation of these elements so the action of these substances actually can be detected by analyzing the membrane noise. But even if some substance does not affect channels, transporters or pumps directly its action often can be detected by noise analysis too. For example if the substance affects G protein-coupled receptors or state of membrane lipids-in many cases it leads to the changes in the functioning of ion channels (Tillman and Cascio, 2003;Inanobe and Kurachi, 2014) and, accordingly, to the noise changes. So the analysis of the long-time series of noise can help to detect the action of many substances when we cannot detect this action differently. For analysis of the long-time series we applied new BRC method based on the beta-distribution function. Four parameters of the beta-distribution function can be used for description of the local fluctuations and the averaged beta-distributions can be applied for quantitative reading of series containing large number of data points. The fluctuation spectroscopy based on beta distribution allows realizing the essential reduction (2.5-10).10 5 data points to 20 quantitative parameters only [see Expressions (14)-(16)] that contain the basic information calculated from three basic beta-distributions: (1) distribution over different segments (scales), (2) the secondary beta-distributions over their heights and (3) distributions over mean values. This reduction becomes possible thanks to the invariant properties that are expressed by formulae (3) and (5). We suppose that this approach can be applied successfully for the unified additional analysis of fluctuations of different long-time series that present the results of monitoring of biological, medical and other data reflecting the results of response of the complex system considered with respect to some external factor. In particular, this BRC method is applicable to testing the action of antagonist of receptor and ion channels when the modification based on different type of interaction (with binding site or with the open channel with different kinetics). In such experiments in order to understand the Frontiers in Computational Neuroscience www.frontiersin.org September 2014 | Volume 8 | Article 120 | 11 mechanism of action of some new substances we only need to compare the FSBD parameter changes caused by this substance with typical changes stored in the database.

FUNDING
This work was partially supported (Andrei I. Skorinkin) by RF grant "Leading Scientific School" and RFBR grant.