Ion channel analysis device

Abstract

To provide an ion channel analysis device capable of obtaining a mathematical parameter for analyzing a state of an ion channel.SOLUTION: An ion channel analysis device 1 includes: a detection unit 20 that detects an ion current I flowing through an ion channel 10; and an analysis unit 30 that acquires first data 40 related to the ion current I detected by the detection unit 20. There is a regression model M that uses, as teacher data, the known first data 40 and a known mathematical parameter 70 corresponding to the known first data 40. The analysis unit 30 inputs the first data 40 related to the ion current I detected by the detection unit 20 to the regression model M, thereby outputting the mathematical parameter 70 corresponding to the input first data 40.SELECTED DRAWING: Figure 1

Description

Technical Field

[0001] The present disclosure relates to an ion channel analysis device.

Background Art

[0002] There are transmembrane proteins called ion channels in the cell membrane of a living body (for example, Patent Document 1). Ion channels serve as passageways for ions in the cell membrane. Ions carry electric charges. In the cell membrane, ions are passively permeated through ion channels. Ion channels open and close.

[0003] When an ion channel opens, the ion current flowing through the ion channel increases. When the ion channel closes, the ion current flowing through the ion channel decreases.

Prior Art Documents

Patent Documents

[0004]

Patent Document 1

Summary of the Invention

Problems to be Solved by the Invention

[0005] By the way, it is conceivable to analyze the state of an ion channel based on data obtained from the ion current flowing through the ion channel (for example, current waveforms, histograms, etc.).

[0006] In order to analyze the state of an ion channel based on data obtained from the ion current flowing through the ion channel, it is necessary to obtain mathematical parameters corresponding to the data, but it has been difficult to obtain such mathematical parameters.

[0007] The present disclosure has been made in view of such points, and an object thereof is to provide an ion channel analysis device capable of obtaining mathematical parameters for analyzing the state of an ion channel.

Means for Solving the Problems

[0008] The ion channel analysis device according to the present disclosure includes a detection unit that detects an ion current flowing through an ion channel, and an analysis unit that obtains first data related to the ion current detected by the detection unit. There is a regression model using the known first data and the known mathematical parameters corresponding to the known first data as teacher data. The analysis unit inputs the first data related to the ion current detected by the detection unit into the regression model, and outputs the mathematical parameters corresponding to the input first data.

[0009] According to such a configuration, the regression model uses the known first data and the known mathematical parameters corresponding to the known first data as teacher data. That is, the regression model learns the correlation between the known first data and the known mathematical parameters.

[0010] The analysis unit inputs the first data related to the ion current detected by the detection unit into the regression model that reflects the correlation between the known first data and the known mathematical parameters. Thereby, the analysis unit can output the mathematical parameters corresponding to the input first data.

[0011] By using the output mathematical parameters, the state of the ion channel can be analyzed.

[0012] As described above, it is possible to provide an ion channel analysis device capable of obtaining mathematical parameters for analyzing the state of an ion channel.

[0013] In one embodiment, the first data includes a current waveform indicating a temporal change in the ion current, and the regression model uses the known current waveform and the known mathematical parameters corresponding to the known current waveform as training data. The analysis unit inputs the current waveform related to the ion current detected by the detection unit into the regression model, and outputs the mathematical parameters corresponding to the input current waveform.

[0014] According to such a configuration, by applying the current waveform as the first data, the analysis unit outputs the mathematical parameters corresponding to the input current waveform. Thereby, mathematical parameters effective for analyzing the state of the ion channel can be obtained.

[0015] In one embodiment, the first data includes a histogram corresponding to a current waveform indicating a temporal change in the ion current and showing the frequency for each value of the ion current, and the regression model uses the known histogram and the known mathematical parameters corresponding to the known histogram as training data. The analysis unit inputs the histogram related to the ion current detected by the detection unit into the regression model, and outputs the mathematical parameters corresponding to the input histogram.

[0016] According to such a configuration, by applying the histogram as the first data, the analysis unit outputs the mathematical parameters corresponding to the input histogram. Thereby, similar to the case of the current waveform, mathematical parameters effective for analyzing the state of the ion channel can be obtained. In particular, when the histogram is applied as the first data, the number of items of the mathematical parameters required for analyzing the state of the ion channel can be reduced compared to the case where the current waveform is applied as the first data.

[0017] In one embodiment, the analysis unit analyzes the state of the ion channel based on the output mathematical parameters.

[0018] According to such a configuration, the operation of analyzing the state of the ion channel can be automated.

[0019] In one embodiment, the mathematical parameter includes information regarding the open / closed state of the ion channel, and the analysis unit analyzes the open / closed state of the ion channel based on the output mathematical parameter.

[0020] According to such a configuration, the open / closed state of the ion channel can be grasped.

[0021] In one embodiment, the mathematical parameter includes information regarding the opening probability as the open / closed state of the ion channel, and the analysis unit obtains the opening probability of the ion channel based on the output mathematical parameter.

[0022] According to such a configuration, the opening probability of the ion channel can be obtained.

[0023] In one embodiment, the ion current I includes a base current value and a peak current value larger than the base current value, and the mathematical parameter includes, as information regarding the opening probability, the base current value, the peak current value, a time constant from the base current value to the peak current value, a time constant from the peak current value to the base current value, a duration at the base current value, and a duration at the peak current value.

[0024] According to such a configuration, the opening probability of the ion channel can be accurately obtained.

[0025] In one embodiment, the ion channel analysis device is applied to the ion channel formed in the artificial cell membrane.

[0026] According to such a configuration, it is effective in conducting experiments using an artificial cell membrane or the like.

Advantages of the Invention

[0027] According to the present disclosure, it is possible to provide an ion channel analysis device capable of obtaining mathematical parameters for analyzing the state of an ion channel.

Brief Description of the Drawings

[0028]

Figure 1

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

Modes for Carrying Out the Invention

[0029] Hereinafter, embodiments of the present disclosure will be described in detail with reference to the drawings. The following description of the preferred embodiments is merely illustrative in nature and is in no way intended to limit the present disclosure, its applications, or its uses.

[0030] <First Embodiment> (Ion Channel Analysis Device) The ion channel analysis device 1 according to the first embodiment will be described. The ion channel analysis device 1 is applied to the ion channel 10. The ion channel analysis device 1 analyzes the state of the ion channel 10. In this example, the ion channel 10 is formed in the artificial cell membrane 5.

[0031] FIG. 1 shows an ion channel 10 formed in an artificial cell membrane 5. The artificial cell membrane 5 mimics the cell membrane of a living body and is artificially formed. When a water droplet 4 is dropped into an organic solvent (oil) in which lipid molecules 2 are dispersed, a monolayer of lipid molecules 2 is spontaneously formed around the water droplet 4. When two water droplets 4 are brought into contact, the monolayers of lipid molecules 2 overlap at the contact portion between the two water droplets 4. An artificial cell membrane 5 composed of a bilayer of lipid molecules 2 is formed between the two water droplets 4. One of the two water droplets 4 mimics the inside of a living cell, and the other of the two water droplets 4 mimics the outside of a living cell.

[0032] The protein 3 that serves as the origin of the ion channel 10 is dispersed in the water droplet 4. When the protein 3 adheres to the artificial cell membrane 5, the ion channel 10 is formed.

[0033] The ion channel 10 is a passage that penetrates the artificial cell membrane 5. The ion channel 10 communicates the two water droplets 4 with each other. Ions passively pass through the ion channel 10 between the two water droplets 4. The ion channel 10 is a passage for ions in the artificial cell membrane 5. Ions carry a charge. Examples of ions include potassium ions and sodium ions.

[0034] The ion channel 10 opens and closes. When the ion channel 10 opens, the ion current I flowing through the ion channel 10 between the two water droplets 4 increases. When the ion channel 10 closes, the ion current I flowing through the ion channel 10 between the two water droplets 4 decreases.

[0035] The ion channel analyzer 1 includes a detection unit 20 and an analysis unit 30. The detection unit 20 is a known ammeter. The detection unit 20 detects the ion current I flowing through the ion channel 10. The detection unit 20 is connected in series to both ends (one water droplet 4 side and the other water droplet 4 side) of the ion channel 10. Further, the detection unit 20 has a function of applying a command voltage to the ion channel 10.

[0036] The analysis unit 30 is connected to the detection unit 20. The analysis unit 30 is built into the ion channel analysis device 1 main body. The analysis unit 30 is, for example, a computer. The analysis unit 30 includes, for example, a processor mounted on a substrate and a memory device that stores software for operating the processor. The analysis unit 30 performs the arithmetic processing described later.

[0037] The analysis unit 30 obtains the first data 40 related to the ion current I detected by the detection unit 20. The analysis unit 30 analyzes the state of the ion channel 10 based on the first data 40 related to the ion current I. In particular, in this example, the analysis unit 30 analyzes the open / closed state of the ion channel 10.

[0038] (Current waveform) The first data 40 includes a current waveform 50. The current waveform 50 shows the time change of the ion current I. FIG. 2 shows the current waveform 50 of the ion current I. In FIG. 2, the horizontal axis represents time t, and the vertical axis represents the value of the ion current I. The unit of time is, for example, [s]. The unit of the ion current I is, for example, [A]. The current waveform 50 is a collection of sets of data of time t and data of the ion current I.

[0039] In the current waveform 50, the ion current I includes a leakage current value IL as a base current value and a peak current value IP. The leakage current value IL corresponds to the value of the ion current I when the ion channel 10 is closed. Even when the ion channel 10 is closed, the ion current I does not completely become zero, but leaks through the artificial cell membrane 5 and flows a little through the ion channel 10. Therefore, the leakage current value IL is slightly larger than zero.

[0040] The magnitude of the leakage current value IL in the ion current I is determined in proportion to the magnitude of the command voltage applied to the ion channel 10. If the applied command voltage is a negative potential, the leakage current value IL in the ion current I also becomes a negative value.

[0041] The peak current value IP is greater than the leakage current value IL in absolute value. The peak current value IP corresponds to the value of the ion current I when the ion channel 10 is open. When the ion channel 10 is open, more ion current I flows through the ion channel 10. Therefore, the peak current value IP is greater than the leakage current value IL in absolute value. According to FIG. 2, the ion current I repeats the leakage current value IL and the peak current value IP. The leakage current value IL is not completely constant and fluctuates randomly within the range of the noise σ. The peak current value IP is not completely constant and fluctuates randomly within the range of the noise σ.

[0042] When the ion channel 10 changes from closed to open, that is, when the ion current I reaches the peak current value IP from the leakage current value IL, it takes a time of the first time constant τo. When the ion channel 10 changes from open to closed, that is, when the ion current I reaches the leakage current value IL from the peak current value IP, it takes a time of the second time constant τc.

[0043] The ion current I persists at the leakage current value IL for the closed duration Tc. The closed durations Tc at the leakage current values IL of the 1st to 4th valleys are denoted as Tc1, Tc2, Tc3, and Tc4. The total closed duration Tc in the target period is the sum of Tc1 to Tc4. That is, Tc = Tc1 + Tc2 + Tc3 + Tc4.

[0044] The ion current I persists at the peak current value IP for the open duration To. The open durations To at the peak current values IP of the 1st to 4th peaks are denoted as To1, To2, To3, and To4. The total open duration To in the target period is the sum of To1 to To4. That is, To = To1 + To2 + To3 + To4.

[0045] As for the open / closed state of the ion channel 10, there is an open probability Po. The open probability Po is the ratio of the open duration To during which the ion channel 10 is open to the target period. The open probability Po is approximately obtained by To / (Tc + To) (Po = To / (Tc + To)).

[0046] As the open / closed state of ion channel 10, there is a closed probability Pc. The closed probability Pc is the ratio of the closed duration Tc during which ion channel 10 is closed to the target period. The closed probability Pc is roughly obtained by Tc / (Tc + To) (Pc = Tc / (Tc + To)).

[0047] The sum of the open probability Po and the closed probability Pc is 1 (Po + Pc = 1). If one of the open probability Po and the closed probability Pc is determined, the other is also determined. Obtaining the open probability Po and obtaining the closed probability Pc are essentially the same. Hereinafter, the open probability Po will be described.

[0048] In the target period, in addition to the closed duration Tc and the open duration To, there are a first time constant τo and a second time constant τc. When accurately obtaining the open probability Po, the first time constant τo and the second time constant τc should not be assigned to either the closed duration Tc or the open duration To.

[0049] For example, an appropriate threshold value IA can be set for the ion current I. It can be considered to assign the time when the ion current I is smaller than the threshold value IA to the closed duration Tc, and the time when the ion current I is larger than the threshold value IA to the open duration To.

[0050] In the case of a human, for the given current waveform 50, the threshold value IA can be set visually each time between the leakage current value IL and the peak current value IP. Thereby, in the case of a human, an accurate open probability Po considering the first time constant τo and the second time constant τc can be obtained.

[0051] However, it is difficult for a computer to automatically and accurately obtain the opening probability Po based on the current waveform 50. In the first place, even if the current waveform 50 is given, the computer cannot recognize the leakage current value IL and the peak current value IP. For this reason, for the given current waveform 50, the computer cannot always appropriately set the threshold value IA in the middle between the leakage current value IL and the peak current value IP. Also, the influence of the noise σ makes it difficult for the computer to recognize the leakage current value IL and the peak current value IP.

[0052] It may be considered to fix the threshold value IA in advance at an appropriate predetermined value, but the form of the current waveform 50 varies each time depending on the type of the ion channel 10 and the magnitude of the command voltage applied to the ion channel 10. For example, both the leakage current value IL and the peak current value IP increase or decrease (the entire current waveform 50 moves upward or downward). For this reason, if the threshold value IA is fixed in advance at an appropriate predetermined value, both the leakage current value IL and the peak current value IP in the current waveform 50 may become larger or smaller than the threshold value IA, and the accurate opening probability Po cannot be obtained.

[0053] Thus, it is difficult for a computer to automatically and accurately obtain the opening probability Po based on the current waveform 50.

[0054] (Histogram) The first data 40 includes a histogram 60. The histogram 60 corresponds to the current waveform 50. Specifically, the histogram 60 is converted from the current waveform 50. The histogram 60 shows the frequency F for each value of the ion current I. FIG. 3 shows the histogram 60 of the ion current I. In FIG. 2, the horizontal axis represents the ion current I, and the vertical axis represents the frequency F. The unit of the frequency F is dimensionless [-]. The frequency F is also a probability.

[0055] The histogram 60 is originally a bar graph, but by increasing the number of data points, it becomes a curve as shown in FIG. 3. In this example, the histogram 60 is composed of a Gaussian distribution. The Gaussian distribution is also called a normal distribution.

[0056] In the histogram 60, a closed-side peak portion 61 corresponding to the closed duration Tc (leak current IL) is formed on the side with a small ion current I (left side), and an open-side peak portion 62 corresponding to the open duration To (peak current IP) is formed on the side with a large ion current I (right side). In the histogram 60, a skirt portion 63 corresponding to the first time constant τo, the second time constant τc, and the noise σ is formed between the closed-side peak portion 61 and the open-side peak portion 62.

[0057] The closed-side peak portion 61 has a closed area A1. The open-side peak portion 62 has an open area A2. The skirt portion 63 has a slight skirt area A3.

[0058] The opening probability Po of the ion channel 10 is approximately obtained by A2 / (A1 + A2) (Po = A2 / (A1 + A2)). The closing probability Pc of the ion channel 10 is approximately obtained by A1 / (A1 + A2) (Pc = A1 / (A1 + A2)).

[0059] As described above, obtaining the opening probability Po and obtaining the closing probability Pc are essentially the same. Therefore, the opening probability Po will be described below.

[0060] In accurately obtaining the opening probability Po, the skirt area A3 must be distributed to either the closed area A1 or the open area A2. For example, an appropriate threshold value IA is set for the ion current I. It is conceivable to distribute the area when the ion current I is smaller than the threshold value IA to the closed area A1, and distribute the area when the ion current I is larger than the threshold value IA to the open area A2.

[0061] For a human, for the given histogram 60, the closed-side peak portion 61 and the open-side peak portion 62 can be visually identified, and the threshold value IA can be set each time in the middle (skirt portion 63) of the two. Thereby, for a human, an accurate opening probability Po considering the skirt area A3 of the skirt portion 63 can be obtained.

[0062] However, it is difficult for a computer to automatically and accurately obtain the opening probability Po based on the histogram 60. In the first place, even if the histogram 60 is given, the computer cannot distinguish between the closed-side peak portion 61 and the open-side peak portion 62. For this reason, the computer cannot appropriately set the threshold value IA at the middle (the skirt portion 63) between the closed-side peak portion 61 and the open-side peak portion 62 each time.

[0063] It may be considered to fix the threshold value IA in advance at an appropriate predetermined value. However, the mode of the histogram 60 varies each time depending on the type of the ion channel 10 and the magnitude of the command voltage applied to the ion channel 10. For example, when both the leak current value IL and the peak current value IP increase or decrease, the entire histogram 60 moves to the right or left. For this reason, when the threshold value IA is fixed in advance at an appropriate predetermined value, in the histogram 60, both the closed-side peak portion 61 and the open-side peak portion 62 may become larger or smaller than the threshold value IA (may be located on the right or left), and the accurate opening probability Po cannot be obtained.

[0064] Thus, it is difficult for a computer to automatically and accurately obtain the opening probability Po based on the histogram 60.

[0065] (Regression Model) FIG. 4 shows the regression model M according to the first embodiment. The analysis unit 30 of the ion channel analyzer 1 uses the regression model M. The regression model M is in the external server 6. Specifically, the regression model M is stored in the external server 6. The analysis unit 30 accesses the regression model M in the external server 6, inputs data to the regression model M, and receives the data output by the regression model M.

[0066] The regression model M is a supervised machine learning model. Specific regression models M include, for example, Ridge regression, GDBT (Gradient Boosting Decision Tree), Multi-Layer Perceptron (MLP), Convolutional Neural Network (CNN) model, Transformer model, and the like.

[0067] The regression model M learns as follows. The regression model M uses the known first data 40 and the known mathematical parameters 70 as the training data. The known mathematical parameters 70 correspond to the known first data 40.

[0068] During the learning of the regression model M, the known first data 40 is given as the training data to the input layer of the regression model M, and the known mathematical parameters 70 are given as the training data to the output layer of the regression model M.

[0069] For example, the regression model M uses the known current waveform 50 and the known mathematical parameters 70 as the training data. At this time, the known mathematical parameters 70 correspond to the known current waveform 50. Also, the known current waveform 50 is given as the training data to the input layer of the regression model M, and the known mathematical parameters 70 are given as the training data to the output layer of the regression model M.

[0070] Alternatively, the regression model M uses the known histogram 60 and the known mathematical parameters 70 as the training data. At this time, the known mathematical parameters 70 correspond to the known histogram 60. Also, the known histogram 60 is given as the training data to the input layer of the regression model M, and the known mathematical parameters 70 are given as the training data to the output layer of the regression model M.

[0071] The mathematical parameter 70 is an input variable corresponding to the first data 40 (current waveform 50 and histogram 60). The known first data 40 (current waveform 50 and histogram 60) is generated by numerical simulation based on the known mathematical parameter 70. While changing the conditions of the mathematical parameter 70, a number of mathematical parameters 70 are set. For each mathematical parameter 70 with different conditions, the corresponding first data 40 (current waveform 50 and histogram 60) is generated. As a result, a number of sets of the first data 40 (current waveform 50 and histogram 60) and the mathematical parameter 70 are prepared according to the conditions.

[0072] The mathematical parameter 70 includes information regarding the open / closed state of the ion channel 10. The open / closed state of the ion channel 10 refers to the state regarding the presence or absence of opening / closing of the ion channel 10, the degree of opening, and the degree of closing. Specifically, the mathematical parameter 70 includes information regarding the open probability Po as the open / closed state of the ion channel 10.

[0073] More specifically, as information regarding the open probability Po of the ion channel 10, the mathematical parameter 70 includes the leak current value IL, the peak current value IP, the first time constant τo from the leak current value IL to reaching the peak current value IP, the second time constant τc from the peak current value IP to reaching the leak current value IL, the closed duration Tc at the leak current value IL, and the open duration To at the peak current value IP.

[0074] Also, the mathematical parameter 70 may include the open probability Po itself as information regarding the open probability Po of the ion channel 10.

[0075] Note that, regardless of whether the known current waveform 50 or the known histogram 60 is used as the teacher data given to the input layer, it is preferable to give these known mathematical parameters 70 (items not enclosed in parentheses in FIG. 4) to the output layer.

[0076] The mathematical parameter 70 may include the number of opening and closing events D (the number of repetitions of the opening and closing of the ion channel 10 during the target period), the opening and closing timing E (the time when the ion channel 10 opens and closes), the closed duration Tc for each opening and closing event (such as Tc1 to Tc4), and the open duration To for each opening and closing event (such as To1 to To4) as information regarding the opening probability Po of the ion channel 10.

[0077] When using the known current waveform 50 as the teacher data given to the input layer, it is preferable to give these known mathematical parameters 70 (the items enclosed in parentheses in FIG. 4) to the output layer. However, when using the known histogram 60 as the teacher data given to the input layer, it is not essential to give these known mathematical parameters 70 (the items enclosed in parentheses in FIG. 4) to the output layer.

[0078] The analysis unit 30 executes the regression model M as follows. The analysis unit 30 inputs the first data 40 related to the ion current I detected by the detection unit 20 into the regression model M, and outputs the mathematical parameter 70 corresponding to the input first data 40. The analysis unit 30 accesses the regression model M in the external server 6 and inputs the first data 40 related to the detected ion current I into the regression model M. The analysis unit 30 acquires the mathematical parameter 70 output by the regression model M in the external server 6 (corresponding to the first data 40). The analysis unit 30 outputs the acquired mathematical parameter 70.

[0079] As the mathematical parameter 70 corresponding to the first data 40, information regarding the opening probability Po (open / closed state) of the ion channel 10 is output.

[0080] For example, the analysis unit 30 inputs the current waveform 50 related to the ion current I detected by the detection unit 20 into the regression model M, and outputs the mathematical parameter 70 corresponding to the input current waveform 50.

[0081] As mathematical parameters 70 corresponding to the current waveform 50, a leakage current value IL, a peak current value IP, a first time constant τo, a second time constant τc, a closed duration Tc, an open duration To, and an open probability Po itself are output. Further, as mathematical parameters 70 corresponding to the current waveform 50, the number of opening and closing events D, the opening and closing timing E, the closed duration Tc for each opening and closing event (for example, Tc1 to Tc4, etc.), and the open duration To for each opening and closing event (for example, To1 to To4, etc.) are output.

[0082] Alternatively, the analysis unit 30 inputs the histogram 60 related to the ion current I detected by the detection unit 20 into the regression model M, and outputs mathematical parameters 70 corresponding to the input histogram 60.

[0083] As mathematical parameters 70 corresponding to the histogram 60, a leakage current value IL, a peak current value IP, a first time constant τo, a second time constant τc, a closed duration Tc, an open duration To, and an open probability Po itself are output.

[0084] The analysis unit 30 analyzes the state of the ion channel 10 based on the output mathematical parameters 70. Specifically, the analysis unit 30 analyzes the open / closed state of the ion channel 10 based on the output mathematical parameters 70.

[0085] More specifically, the analysis unit 30 obtains the open probability Po of the ion channel 10 based on the output mathematical parameters 70. Since specific numerical values are given as the mathematical parameters 70, the analysis unit 30 can easily obtain the open probability Po.

[0086] For example, the analysis unit 30 obtains the open probability Po based on the mathematical parameters 70 corresponding to the histogram 60 (see FIG. 3). In the mathematical parameters 70 corresponding to the histogram 60, the leakage current value IL and the peak current value IP are given as specific numerical values. The analysis unit 30 calculates the intermediate value between the leakage current value IL and the peak current value IP, and sets this intermediate value as the threshold value IA.

[0087] In the histogram 60, the analysis unit 30 distributes the area of the portion on the side (left side) of the ion current I smaller than the threshold value IA to the closed area A1 of the closed-side peak portion 61, and distributes the area of the portion on the side (right side) of the ion current I larger than the threshold value IA to the open area A2 of the open-side peak portion 62. At this time, the skirt area A3 of the skirt portion 63 is distributed to the closed area A1 and the open area A2. Thereby, the analysis unit 30 can easily obtain the opening probability Po by Po = A2 / (A1 + A2).

[0088] Also, in the mathematical parameter 70 corresponding to the histogram 60, the opening probability Po itself is given by the regression model M. The analysis unit 30 may directly apply the opening probability Po as the mathematical parameter 70 given by the regression model M.

[0089] For example, by comparing the opening probability Po obtained using the above-described threshold value IA with the opening probability Po as the mathematical parameter 70 given by the regression model M, the validity (reliability) of the mathematical parameter 70 obtained from the regression model M may be confirmed.

[0090] Alternatively, the analysis unit 30 may obtain the opening probability Po based on the mathematical parameter 70 corresponding to the current waveform 50 (see FIG. 2). Note that when based on the mathematical parameter 70 corresponding to the current waveform 50, the number of items of the mathematical parameter 70 increases (see FIG. 4) compared to the case of being based on the mathematical parameter 70 corresponding to the histogram 60, and thus the calculation becomes complicated.

[0091] Also when based on the mathematical parameter 70 corresponding to the current waveform 50, similar to the case of being based on the mathematical parameter 70 corresponding to the histogram 60, the opening probability Po as the mathematical parameter 70 given by the regression model M may be directly applied.

[0092] (Operation and Effect of the First Embodiment) In the ion channel analysis device 1 according to the present embodiment, the detection unit 20 detects the ion current I flowing through the ion channel 10, and the analysis unit 30 obtains the first data 40 detected by the detection unit 20.

[0093] The regression model M uses the known first data 40 and the known mathematical parameters 70 corresponding to the known first data 40 as training data. That is, the regression model M is made to learn the correlation between the known first data 40 and the known mathematical parameters 70.

[0094] The analysis unit 30 inputs the first data 40 related to the ion current I detected by the detection unit 20 into the regression model M that reflects the correlation between the known first data 40 and the known mathematical parameters 70. As a result, the analysis unit 30 can output the mathematical parameters 70 corresponding to (correlated with) the input first data 40.

[0095] Based on the output mathematical parameters 70, the state of the ion channel 10 can be analyzed.

[0096] As described above, it is possible to provide the ion channel analysis device 1 capable of obtaining the mathematical parameters 70 for analyzing the state of the ion channel 10.

[0097] The regression model M is not just a data set that simply collects the known first data 40 and the known mathematical parameters 70, but a learning model. Since the regression model M is a learning model, even if the first data 40 input into the regression model M at the time of execution is in a complex form (for example, the noise σ of the current waveform 50 is large as shown in FIG. 2, or the tail part 63 of the histogram 60 is larger than zero as shown in FIG. 3), the regression model M with improved accuracy through learning outputs reliable and appropriate mathematical parameters 70.

[0098] By applying the current waveform 50 as the first data 40, the analysis unit 30 outputs mathematical parameters 70 corresponding to the input current waveform 50. Thereby, mathematical parameters 70 effective for analyzing the state of the ion channel 10 can be obtained.

[0099] By applying the histogram 60 as the first data 40, the analysis unit 30 outputs mathematical parameters 70 corresponding to the input histogram 60. Thereby, similar to the case of the current waveform 50, mathematical parameters 70 effective for analyzing the state of the ion channel 10 can be obtained.

[0100] In particular, when the histogram 60 is applied as the first data 40, the number of items of the mathematical parameters 70 required for analyzing the state of the ion channel 10 can be reduced compared to the case where the current waveform 50 is applied as the first data 40 (see FIG. 4). Thereby, it becomes easier to analyze the state of the ion channel 10 based on the mathematical parameters 70.

[0101] By the analysis unit 30 itself analyzing the state of the ion channel 10 based on the mathematical parameters 70, the operation of analyzing the state of the ion channel 10 can be automated.

[0102] Since the mathematical parameters 70 include information regarding the open / closed state of the ion channel 10, the open / closed state of the ion channel 10 can be grasped.

[0103] In particular, since the mathematical parameters 70 include information regarding the open probability Po as the open / closed state of the ion channel 10, the open probability Po of the ion channel 10 can be obtained.

[0104] The mathematical parameters 70 include, as information regarding the open probability Po, the leakage current value IL, the peak current value IP, the first time constant τo, the second time constant τc, the closed duration Tc, and the open duration To. Thereby, based on the mathematical parameters 70, the open probability Po of the ion channel 10 can be accurately obtained.

[0105] Since the ion channel analyzer 1 is applied to the ion channel 10 formed in the artificial cell membrane 5, it is effective for conducting experiments (such as those related to pharmaceuticals) using the artificial cell membrane 5.

[0106] (Modification of the First Embodiment) In the above embodiment, the threshold value IA is used to obtain the opening probability Po of the ion channel 10. However, the present invention is not limited to this, and the opening probability Po may be obtained by other methods.

[0107] In the above embodiment, the opening probability Po of the ion channel 10 is obtained based on the mathematical parameter 70. However, the present invention is not limited to this. The closing probability Pc may be obtained based on the mathematical parameter 70 (obtaining the opening probability Po and the closing probability Pc are essentially the same), or other opening and closing states of the ion channel 10 (such as the degree of opening or closing of the ion channel 10) may be analyzed, or other states of the ion channel 10 (such as the conductance that can be calculated from the applied voltage and the observed ion current in the ion channel 10) may be analyzed.

[0108] In the above embodiment, the analysis unit 30 analyzes the state of the ion channel 10 based on the output mathematical parameter 70. However, the present invention is not limited to this. The analysis unit 30 may perform operations until the mathematical parameter 70 is output, and then a person may analyze the state of the ion channel 10 based on the mathematical parameter 70.

[0109] In the above embodiment, the current waveform 50 and the histogram 60 are exemplified as the first data 40. However, the present invention is not limited to this. For example, as the first data 40, conductance (unit: Siemens) obtained by dividing the ion current I by the command voltage may be used.

[0110] The regression model M may be stored in the memory of the analysis unit 30 instead of being stored in the external server 6.

[0111] <Second Embodiment> (Mixture Gaussian Model) The ion channel analysis device 1 according to the second embodiment will be described. In the following description, components having the same configuration as those in the above embodiment may be denoted by the same reference numerals, and detailed descriptions thereof may be omitted. FIG. 5 shows a first analysis target histogram 80 and a reference histogram 90 according to the second embodiment.

[0112] The detection unit 20 detects the ion current I flowing through the ion channel 10. The analysis unit 30 obtains a first analysis target histogram 80 based on the current waveform 50 related to the ion current I detected by the detection unit 20. The current waveform 50 shows the temporal change of the ion current I (see FIG. 2).

[0113] The first analysis target histogram 80 is the same as the histogram 60 according to the first embodiment (see FIG. 3). However, although details will be described later, the first analysis target histogram 80 has a smaller number of data compared to the second analysis target histogram 85 described later.

[0114] The first analysis target histogram 80 corresponds to the current waveform 50. Specifically, the first analysis target histogram 80 is converted from the current waveform 50. The first analysis target histogram 80 shows the frequency F for each value of the ion current I.

[0115] The first analysis target histogram 80 is composed of a Gaussian distribution (normal distribution). In the first analysis target histogram 80, a closed-side peak 81 corresponding to the closed duration Tc (leakage current IL) is formed on the smaller side (left side) of the ion current I, and an open-side peak 82 corresponding to the open duration To (peak current IP) is formed on the larger side (right side) of the ion current I. In the first analysis target histogram 80, a skirt portion 83 corresponding to the first time constant τo, the second time constant τc, and the noise σ is formed between the closed-side peak 81 and the open-side peak 82.

[0116] When the histogram is composed of a Gaussian distribution as in the first analysis target histogram 80, the number of peaks such as the closed-side peak 81 and the open-side peak 82 in the first analysis target histogram 80 is referred to as the number of clusters N. The number of clusters N of the first analysis target histogram 80 is 2 (N = 2).

[0117] In order to analyze the state (especially the open / closed state) of the ion channel 10, it is necessary to determine the number of clusters N of the histogram. For example, when the number of clusters N is 2, it can be seen that there are a cluster corresponding to when the ion channel 10 is closed (for example, the closed-side peak 81) and a cluster corresponding to when the ion channel 10 is open (for example, the open-side peak 82). If the number of clusters N is known, the state (especially the open / closed state) of the ion channel 10 can be known.

[0118] As can be seen from FIG. 5, in the case of a human, by visually observing the first analysis target histogram 80, the closed-side peak 81 and the open-side peak 82 can be distinguished, so it can be easily determined that the number of clusters N of the first analysis target histogram 80 is 2.

[0119] However, even if the first analysis target histogram 80 is given to the computer, the computer cannot distinguish between the closed-side peak 81 and the open-side peak 82. For this reason, it is difficult for the computer to determine the number of clusters N of the first analysis target histogram 80, and some kind of ingenuity is required.

[0120] The analysis unit 30 of the ion channel analysis device 1 uses the mixture Gaussian model G. The mixture Gaussian model G is in the external server 6. Specifically, the mixture Gaussian model G is stored in the external server 6. The analysis unit 30 accesses the mixture Gaussian model G in the external server 6, inputs data to the mixture Gaussian model G, and receives data output by the mixture Gaussian model G.

[0121] The mixture Gaussian model G is an unsupervised machine learning model. The mixture Gaussian model G is one of the clustering methods, and it is a method of analyzing by considering that the data is generated from a plurality of Gaussian distributions (normal distributions).

[0122] The mixture Gaussian model G in the external server 6 has a plurality of reference histograms 90. The plurality of reference histograms 90 are generated by the mixture Gaussian model G. The reference histogram 90 is composed of Gaussian distributions for each number N of clusters by the mixture Gaussian model G. In other words, the reference histogram 90 is generated by the mixture Gaussian model G for each number N of clusters.

[0123] In this example, the mixture Gaussian model G generates five reference histograms 90 composed of Gaussian distributions with the number N of clusters ranging from 1 to 5. When there is only one ion channel 10, two numbers N of clusters are sufficient, but when two or more ion channels 10 may exist, three or more numbers N of clusters are required.

[0124] The analysis unit 30 inputs the first analysis target histogram 80 into the mixture Gaussian model G in the external server 6. The mixture Gaussian model G generates and outputs a plurality of reference histograms 90 corresponding to the input first analysis target histogram 80. The analysis unit 30 acquires the plurality of reference histograms 90 output by the mixture Gaussian model G.

[0125] The analysis unit 30 compares the first analysis target histogram 80 with the reference histograms 90 for each number N of clusters. Specifically, the analysis unit 30 compares the first analysis target histogram 80 with the reference histogram 90 with the number N of clusters being 1 (N = 1), the reference histogram 90 with the number N of clusters being 2 (N = 2), the reference histogram 90 with the number N of clusters being 3 (N = 3), the reference histogram 90 with the number N of clusters being 4 (N = 4), and the reference histogram 90 with the number N of clusters being 5 (N = 5) in sequence.

[0126] The analysis unit 30 determines the number N of clusters when the first analysis target histogram 80 is closest to the reference histogram 90 as the number N of clusters of the first analysis target histogram 80. For example, the analysis unit 30 obtains the degree of coincidence between the first analysis target histogram 80 and the reference histogram 90 for each number N of clusters, and determines the number N of clusters when the degree of coincidence is the highest (the degree of matching is the highest, the difference is the smallest, etc.).

[0127] In this example, the first analysis target histogram 80 is closest to the reference histogram 90 with the number of clusters N being 2 (N = 2). The analysis unit 30 determines that the number of clusters N of the first analysis target histogram 80 is 2 (N = 2).

[0128] Furthermore, the analysis unit 30 determines whether the number of clusters N of the first analysis target histogram 80 is 2 or more (N ≥ 2?).

[0129] When the number of clusters N of the first analysis target histogram 80 is 2 or more (N ≥ 2), the analysis unit 30 determines that the ion channel 10 is open and closed. This is because when the number of clusters is 2 or more, it can be seen that there are at least a cluster corresponding to when the ion channel 10 is closed (for example, the closed-side peak 81) and a cluster corresponding to when the ion channel 10 is open (for example, the open-side peak 82). Note that "the ion channel 10 is open and closed" means that the ion channel 10 has two states, an open state and a closed state.

[0130] When the number of clusters N of the first analysis target histogram 80 is less than 2 (N < 2), that is, when the number of clusters N of the first analysis target histogram 80 is 1 (N = 1), the analysis unit 30 determines that the ion channel 10 is not open and closed (the ion channel 10 is always open or always closed).

[0131] After determining the number of clusters N (= 2) of the first analysis target histogram 80, the analysis unit 30 obtains a second analysis target histogram 85 based on the current waveform 50. Specifically, the analysis unit 30 converts from the current waveform 50 to the second analysis target histogram 85. The second analysis target histogram 85 is the same as the histogram 60 according to the first embodiment (see FIG. 3).

[0132] The second analysis target histogram 85 corresponds to the current waveform 50. The second analysis target histogram 85 corresponds to the first analysis target histogram 80. The second analysis target histogram 85 has a larger number of data than the first analysis target histogram 80. In other words, the first analysis target histogram 80 has a smaller number of data than the second analysis target histogram 85.

[0133] Based on the second analysis target histogram 85 and the determined number of clusters N (= 2), the analysis unit 30 outputs the Gaussian distribution mathematical parameters 100 corresponding to the second analysis target histogram 85. The Gaussian distribution mathematical parameters 100 are information regarding the Gaussian distribution that constitutes the second analysis target histogram 85, and are, for example, the mean, variance, standard deviation, specific gravity, etc.

[0134] Based on the output Gaussian distribution mathematical parameters 100, the analysis unit 30 may analyze the state of the ion channel 10 (for example, the opening probability Po, etc.).

[0135] FIG. 6 is a flowchart of the determination of the number of clusters N according to the second embodiment. Starting from the start, in the first step S1, the detection unit 20 detects the ion current I flowing through the ion channel 10.

[0136] In the second step S2, based on the current waveform 50 related to the ion current I detected by the detection unit 20, the analysis unit 30 obtains the first analysis target histogram 80. The first analysis target histogram 80 has a smaller number of data than the second analysis target histogram 85.

[0137] In the third step S3, the analysis unit 30 inputs the first analysis target histogram 80 into the mixture Gaussian model G in the external server 6.

[0138] In the fourth step S4, the mixture Gaussian model G outputs a plurality of reference histograms 90 corresponding to the first analysis target histogram 80 for each number of clusters N.

[0139] In the fifth step S5, the analysis unit 30 acquires a plurality of reference histograms 90 for each number of clusters N output by the mixture Gaussian model G. Then, the analysis unit 30 compares the first analysis target histogram 80 with the reference histograms 90 for each number of clusters N.

[0140] In the sixth step S6, the analysis unit 30 determines the number of clusters N (= 1 to 5) when the first analysis target histogram 80 is closest to the reference histogram 90 as the number of clusters N of the first analysis target histogram 80. In this example, the number of clusters N of the first analysis target histogram 80 is determined to be 2 (N = 2).

[0141] In the seventh step S7, the analysis unit 30 determines whether the number of clusters N of the first analysis target histogram 80 is 2 or more (N ≧ 2?). If it is determined that the number of clusters N of the first analysis target histogram 80 is 2 or more (N ≧ 2), the process proceeds to the eighth step S8. If it is determined that the number of clusters N of the first analysis target histogram 80 is less than 2 (N < 2), that is, if it is determined that the number of clusters N of the first analysis target histogram 80 is 1 (N = 1), the process proceeds to the ninth step S9.

[0142] In the eighth step S8 (the number of clusters N is 2 or more (N ≧ 2)), the analysis unit 30 determines that the ion channel 10 is open and closed. When the eighth step S8 ends, the process proceeds to the tenth step S10.

[0143] In the ninth step S9 (the number of clusters N is less than 2 (N < 2)), the analysis unit 30 determines that the ion channel 10 is not open and closed (the ion channel 10 is always open or always closed). When the ninth step S9 ends, the process proceeds to the tenth step S10.

[0144] In the tenth step S10, the analysis unit 30 obtains a second analysis target histogram 85 with a larger number of data than the first analysis target histogram 80 based on the current waveform 50.

[0145] In the 11th step S11, the analysis unit 30 outputs a Gaussian distribution mathematical parameter 100 corresponding to the second analysis target histogram 85 based on the second analysis target histogram 85 and the determined number of clusters N (= 2). And then it reaches the end.

[0146] (Operation and Effect of the Second Embodiment) According to the ion channel analysis device 1 according to the present embodiment, by comparing the first analysis target histogram 80 with the reference histogram 90 for each number of clusters N and obtaining the number of clusters N when the first analysis target histogram 80 is closest to the reference histogram 90, the number of clusters N of the first analysis target histogram 80 can be determined.

[0147] Here, as described above, the number of clusters N is a factor effective for analyzing the state (especially the open / closed state) of the ion channel 10.

[0148] As described above, it is possible to provide the ion channel analysis device 1 capable of obtaining the number of clusters N for analyzing the state of the ion channel 10.

[0149] Without relying on human visual inspection, the number of clusters N of the first analysis target histogram 80 can be grasped.

[0150] Based on whether the number of clusters N of the first analysis target histogram 80 is 2 or more, it is possible to easily determine whether the ion channel 10 is open or closed.

[0151] A reference histogram 90 with high accuracy can be generated by the mixture Gaussian model G.

[0152] By using the first analysis target histogram 80 with a small number of data, the calculation for determining the number of clusters N can be speeded up. On the other hand, by using the second analysis target histogram 85 with a large number of data, a Gaussian distribution mathematical parameter 100 with high accuracy can be obtained.

[0153] Since the ion channel analyzer 1 is applied to the ion channel 10 formed in the artificial cell membrane 5, it is effective in conducting experiments using the artificial cell membrane 5 and the like.

[0154] (Modification of the Second Embodiment) As described above, when the number of clusters N of the first analysis target histogram 80 is less than 2 (N < 2), that is, when the number of clusters N of the first analysis target histogram 80 is 1 (N = 1), it is determined that the ion channel 10 is not opening and closing, that is, the ion channel 10 is either always open or always closed.

[0155] In order to determine whether the ion channel 10 is always open or always closed, a threshold value IA related to the ion current I is set for the first analysis target histogram 80, and it may be determined whether the ion channel 10 is always open or always closed depending on whether the ion current I is greater than or less than the threshold value IA. In this case, if the ion current I is greater than the threshold value IA, it is determined that the ion channel 10 is always open, and if the ion current I is less than the threshold value IA, it is determined that the ion channel 10 is always closed.

[0156] In the above embodiment, the analysis unit 30 determines whether the ion channel 10 is opening and closing based on whether the number of clusters N of the first analysis target histogram 80 is 2 or more, but it is not limited thereto. The analysis unit 30 may stop at determining the number of clusters N. Based on the number of clusters N determined by the analysis unit 30, a person may analyze the state of the ion channel 10.

[0157] The mixture Gaussian model G may be stored in the memory of the analysis unit 30 instead of being stored in the external server 6.

[0158] <Other Embodiments> As described above, the present disclosure has been described with reference to preferred embodiments, but such descriptions are not limiting matters, and of course, various modifications, substitutions, or combinations are possible.

[0159] The analysis unit 30 may be provided outside the main body of the ion channel analysis device 1 instead of being built into the main body of the ion channel analysis device 1.

[0160] In the above embodiment, the ion channel analysis device 1 is applied to the artificial cell membrane 5, but it is not limited thereto, and for example, it may be applied to a biological cell membrane or the like.

[0161] In the above embodiment, the case where there is one ion channel 10 is illustrated, but it is not limited thereto. There may be a plurality of ion channels 10. In this case, the current waveform 50 may be a superposition of the currents I flowing through the plurality of ion channels 10.

Industrial Applicability

[0162] Since the present disclosure can be applied to the ion channel analysis device 1, it is extremely useful and has high industrial applicability.

Explanation of Signs

[0163] I Ion current IL Leak current value (base current value) IP Peak current value IA Threshold value t Time τo First time constant (time constant) τc Second time constant (time constant) To Open duration (duration) Tc Closed duration (duration) σ Noise Po Open probability Pc Closed probability F Frequency A1 Closed area A2 Open area A3 Skirt area D Number of opening and closing events E Opening and closing timing M Regression model G Mixture Gaussian model N Number of clusters 1 Ion channel analysis device 2 lipid molecules 3 organic solvents 4 water droplets 5 artificial cell membranes 6 external servers 10 ion channels 20 detection unit 30 analysis unit 40 First data 50 current waveform 60 histogram 61 closed-side peak 62 open-side peak 63 base 70 mathematical parameters 80 First histogram to be analyzed 81 closed-side peak 82 open-side peak 83 base 85 Second histogram to be analyzed 90 reference histogram 100 Gaussian distribution mathematical parameters (mathematical parameters)

Claims

1. A detection unit that detects an ion current flowing through an ion channel, and an analysis unit that obtains first data related to the ion current detected by the detection unit. The ion channel analysis device is provided with: There is a regression model that uses the known first data and known mathematical parameters corresponding to the known first data as training data, and the analysis unit outputs the mathematical parameters corresponding to the input first data by inputting the first data related to the ion current detected by the detection unit into the regression model.

2. The first data includes a current waveform indicating a temporal change in the ion current, the regression model uses the known current waveform and known mathematical parameters corresponding to the known current waveform as training data, and the analysis unit outputs the mathematical parameters corresponding to the input current waveform by inputting the current waveform related to the ion current detected by the detection unit into the regression model. The ion channel analysis device according to claim 1.

3. The first data includes a histogram corresponding to a current waveform indicating a temporal change in the ion current and showing the frequency for each value of the ion current, the regression model uses the known histogram and known mathematical parameters corresponding to the known histogram as training data, and the analysis unit outputs the mathematical parameters corresponding to the input histogram by inputting the histogram related to the ion current detected by the detection unit into the regression model. The ion channel analysis device according to claim 1 or 2.

4. The analysis unit analyzes the state of the ion channel based on the output mathematical parameters. The ion channel analysis device according to claim 1 or 2.

5. The mathematical parameter includes information regarding the open / closed state of the ion channel. The ion channel analysis device according to claim 1 or 2, wherein the analysis unit analyzes the open / closed state of the ion channel based on the output mathematical parameter.

6. The mathematical parameter includes information regarding the opening probability as the open / closed state of the ion channel. The ion channel analysis device according to claim 5, wherein the analysis unit obtains the opening probability of the ion channel based on the output mathematical parameter.

7. The ion current I includes a base current value and a peak current value greater than the base current value. The mathematical parameter includes, as information regarding the opening probability, the base current value, the peak current value, a time constant from the base current value to the peak current value, a time constant from the peak current value to the base current value, a duration at the base current value, and a duration at the peak current value. The ion channel analysis device according to claim 6.

8. The ion channel analysis device according to claim 1 or 2, applied to the ion channel formed in the artificial cell membrane.

Images