Next Article in Journal
Sexual Dimorphism in Kisspeptin Signaling
Next Article in Special Issue
CaMKII Inhibition Attenuates Distinct Gain-of-Function Effects Produced by Mutant Nav1.6 Channels and Reduces Neuronal Excitability
Previous Article in Journal
Investigation of Colonic Regeneration via Precise Damage Application Using Femtosecond Laser-Based Nanosurgery
Previous Article in Special Issue
KCNQ2 Selectivity Filter Mutations Cause Kv7.2 M-Current Dysfunction and Configuration Changes Manifesting as Epileptic Encephalopathies and Autistic Spectrum Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GABA Receptors Can Depolarize the Neuronal Membrane Potential via Quantum Tunneling of Chloride Ions: A Quantum Mathematical Study

by
Sager Nawafleh
1,*,
Abdallah Barjas Qaswal
2,*,
Aiman Suleiman
3,
Obada Alali
4,
Fuad Mohammed Zayed
2,
Mohammad Abu Orabi Al-Adwan
2 and
Mo’ath Bani Ali
2
1
Department of Anesthesia and Intensive Care Unit, The Hashemite University, Zarqa 13115, Jordan
2
School of Medicine, The University of Jordan, Amman 11942, Jordan
3
Department of Anesthesia, Intensive Care and Pain Management, Harvard Medical School, Beth Israel Deaconess Medical Center, Boston, MA 02215, USA
4
Department of Anesthesia and Intensive Care, Alabdali Clemenceau Hospital, Amman 11190, Jordan
*
Authors to whom correspondence should be addressed.
Cells 2022, 11(7), 1145; https://doi.org/10.3390/cells11071145
Submission received: 28 February 2022 / Revised: 22 March 2022 / Accepted: 24 March 2022 / Published: 28 March 2022

Abstract

:
GABA (gamma-aminobutyric acid) receptors represent the major inhibitory receptors in the nervous system and their inhibitory effects are mediated by the influx of chloride ions that tends to hyperpolarize the resting membrane potential. However, GABA receptors can depolarize the resting membrane potential and thus can also show excitatory effects in neurons. The major mechanism behind this depolarization is mainly attributed to the accumulation of chloride ions in the intracellular compartment. This accumulation leads to increase in the intracellular chloride concentration and depolarize the Nernst potential of chloride ions. When the membrane potential is relatively hyperpolarized, this will result in a chloride efflux instead of influx trying to reach their depolarized equilibrium potential. Here, we propose different mechanism based on a major consequence of quantum mechanics, which is quantum tunneling. The quantum tunneling model of ions is applied on GABA receptors and their corresponding chloride ions to show how chloride ions can depolarize the resting membrane potential. The quantum model states that intracellular chloride ions have higher quantum tunneling probability than extracellular chloride ions. This is attributed to the discrepancy in the kinetic energy between them. At physiological parameters, the quantum tunneling is negligible to the degree that chloride ions cannot depolarize the membrane potential. Under certain conditions such as early neuronal development, gain-of-function mutations, stroke and trauma that can lower the energy barrier of the closed gate of GABA receptors, the quantum tunneling is enhanced so that the chloride ions can depolarize the resting membrane potential. The major unique feature of the quantum tunneling mechanism is that the net efflux of chloride ions is attained without the need for intracellular accumulation of chloride ions as long as the energy barrier of the gate is reduced but still higher than the kinetic energy of the chloride ion as a condition for quantum tunneling to take place.

Graphical Abstract

1. Introduction

GABA (gamma-aminobutyric acid) receptors are pentameric ligand-gated receptors and are considered to be inhibitory receptors because they hyperpolarize the resting membrane potential by the influx of chloride ions. However, GABA receptors can depolarize the membrane potential if there is an intracellular chloride accumulation that increases the intracellular concentration and makes its Nernst potential less negative. Accordingly, when the chloride’s Nernst potential becomes less negative than the resting membrane potential, the opening of GABA receptors will result in chloride efflux instead of chloride influx. This chloride efflux will result in membrane depolarization instead of hyperpolarization [1,2].
The action of depolarization mediated by GABA receptors occurs physiologically during early development or pathologically in certain conditions such as trauma, stroke, and epilepsy [1,3,4,5,6]. The main mechanism behind this depolarization is the increase in the intracellular chloride ions, which makes the Nernst potential less negative [7,8]. The molecular mechanism behind this alteration in the homeostasis of chloride ions concentrations is the abnormal expression of NKCC1 transporter, which loads the chloride ions inside neurons, and KCC2 transporter, which extrudes chloride ions outside the neurons [4,5,6,8]. The higher expression of NKCC1 and the lower expression of KCC2 leads to increased intracellular chloride ions concentration, which results in membrane depolarization mediated by chloride ions efflux [4,5,6,8].
However, because the intracellular chloride concentrations cannot be reliably measured depending on the expression levels of chloride transporters and the estimation of the intracellular concentration is not accurate since the in-vitro studies do not take into consideration the neuronal activity, the actual intracellular chloride concentration may be not that high to induce membrane depolarization in the corresponding conditions. Hence, an alternative mechanism may be required to explain the depolarization action mediated by GABA receptors [8].
The model of quantum tunneling of ions has been studied previously in hypothetical biological frameworks and was found to generate solid mathematical solutions on atomic and molecular levels in nascent fields of neurophysiology and pharmacokinetics [9,10,11].
In our paper, we present a different explanation of the depolarization induced by GABA receptors based on the model of quantum tunneling of chloride ions through the closed gate. The quantum model will be able to explain the depolarization action of the GABA receptors without the requirement of intracellular chloride accumulation. This will provide an alternative mechanism of depolarization in case the increase in the intracellular chloride concentration is not sufficient to depolarize the membrane potential. Moreover, specific features of the quantum model will be presented to exhibit its distinction over the classical model. Moreover, as the action of depolarization in the context of GABA receptor’s function can result in either excitation or inhibition according to a recent study [8], we do not aim to explore this idea in the context of the quantum model, but we aim to focus on the action of depolarization itself. However, we will propose certain mechanisms by which the quantum model can be exploited in such a way that can be tested experimentally and be used to explain the physiological and pathophysiological effects of depolarization induced by GABA receptors.

2. Materials and Methods

The Mathematical Model of Chloride Ions Quantum Tunneling

The model of quantum tunneling of ions has been applied on the voltage-gated channels [12,13]. Herein, we apply the quantum model on the closed GABA receptors and chloride ions. This is valid because the voltage-gated channel and the GABA receptor possess a closed hydrophobic gate that blocks the permeation of ions [14,15].
The quantum tunneling probability of chloride ions through the closed gate of GABA receptors can be calculated by the following equation [12,13]:
T Q = e 8 m ħ × 2 w 3 g ( g K E ) 3
where T Q is the probability of tunneling, m is the mass of chloride ion (5.9 × 10−26 kg), g is the height of the energy barrier of the closed gate in the GABA receptor or simply ‘the energy of the gate’, w is the length of the gate, and KE is the kinetic energy of chloride ion.
Equation (1) has been applied before on the closed gate of the voltage-gated channels, and in the present study, we are going to apply the model of quantum tunneling of chloride ions on the closed gate of the GABA receptor. The hydrophobic gating in voltage-gated channels and pentameric ligand-gated channels, including GABA receptors has been shown to be a reasonable mechanism that controls the flow of ions [14,15,16,17,18]. The hydrophobic gate forms an energy barrier higher than the kinetic energy of the passing ions [14,15,16,17,18]. This gate is formed by hydrophobic residues that come close to each other to form a narrow pore. In the previous studies, the closed gate was viewed as barrier similar to barrier shape of a regular electric field [12,13]. This means that the energy barrier takes the shape of a triangle, see Figure 1. The experimental evidence showed that the energy barrier for the ion while passing through a hydrophobic material, which includes a hydrophobic membrane or a hydrophobic gate in ion channels, changes in a shape similar to the shape of the barrier of a regular electric field or the shape of a triangle [17,18,19,20,21,22,23], which is represented by two possible triangle shapes in Figure 1 that are indicated by the green and red curves. We don’t say that the hydrophobic gate is a regular electric field, rather we say that the barrier shape of the gate can be viewed as the barrier shape of a regular electric field based on the experimental observations [17,18,19,20,21,22,23]. The curves in Figure 1 may not reflect the exact relationship between the energy barrier and the ion’s position through the gate, but they provide a reasonable approximation of the barrier shape based on the experimental results [17,18,19,20,21,22,23]. This approximation is rational to be made because there is no clear mathematical equation that can describe how the energy barrier changes with respect to the ion’s position; hence approximating the energy profile as a linear change will make it easier to obtain the mathematical relationship, as it will be explained in the following paragraphs.
Accordingly, the barrier energy of the closed hydrophobic gate can be described by the following equation according to the green curve in Figure 1:
U ( x ) = g w x
where g is the height of the triangle or the height of the energy barrier or simply ‘the energy of the gate’, w is the width of the barrier or the length of the gate, and x is the position of ion while passing through the closed gate.
The original equation that calculates the tunneling probability of a particle through a barrier of arbitrary shape [24]:
T Q = e 8 m ħ X 1 X 2 U ( x ) K E d x
The integral in Equation (3) can be solved as the following:
X 1 X 2 g w x K E d x = 2 w 3 g ( g w x 2 K E ) 3 2 w 3 g ( g w x 1 K E ) 3
x2 is at the end of the gate ( x 2 = w ) and x1 is where U(x1) = KE, hence Equation (4) becomes:
X 1 X 2 g w x K E d x = 2 w 3 g ( g K E ) 3
By substituting Equation (5) in Equation (3), we get Equation (1).
Our solved integral is based on the green curve in Figure 1, but if the integral is solved according to the red curve, we will get the same result as the following:
The left red line can be represented by the following equation U ( x ) = 2 g w x , while the right red line can be represented by the following equation U ( x ) = 2 g w ( x w ) and if the integral is taken through the forbidden region from x1′ to x2′:
x 1 w 2 2 g w x K E d x + w 2 x 2 2 g w ( x w ) K E d x = ( w 3 g ( 2 g w ( w 2 ) K E ) 3 w 3 g ( 2 g w ( x 1 ) K E ) 3 ) + ( w 3 g ( 2 g w ( x 2 w ) K E ) 3 w 3 g ( 2 g w ( w 2 w ) K E ) 3 )
x1′ and x2′ are where U(x) = KE, hence:
x 1 w 2 2 g w x K E d x + w 2 x 2 2 g w ( x w ) K E d x = w 3 g ( g K E ) 3 + w 3 g ( g K E ) 3 = 2 w 3 g ( g K E ) 3
Accordingly, both barrier shapes are equivalent to estimate the tunneling probability because the solved integrals in both cases are the same. For this reason, we can say that the energy barrier of the hydrophobic gate can be viewed as the barrier shape of a regular electric field based on the approximations made on the actual curves of the energy profiles of ions while passing through the closed hydrophobic gate and hydrophobic membranes [17,18,19,20,21,22,23]. So, in a conclusion, the barrier shapes in Figure 1 are consistent and reasonable for the following reasons: (1) the curves in Figure 1 are triangle barriers similar to the barrier shapes obtained based on the potential mean force calculations in these references [17,18,19,20,21,22,23], (2) up to authors’ knowledge, there is no clear mathematical equation that predicts how the barrier energy changes with respect to the ion’s position in a hydrophobic material.
As chloride ions are present in two compartments, which are intracellular and extracellular compartments, we expect different kinetic energy for the intracellular and extracellular chloride ions. As chloride ions are anions and the membrane potential of cells are negative inside with respect to outside, we expect that intracellular chloride ions will get a kinetic energy from the membrane potential as they go from the intracellular compartment to the extracellular compartment through the GABA receptor. However, the kinetic energy of the extracellular chloride ions will decrease as they go through the membrane potential. See Figure 2.
As this article [15] indicates that the closed gate is in the middle of the GABA receptor, other locations cannot be excluded. Therefore, we will make our investigations based on different locations. We will refer to the location of the gate by n and we will choose three values for n, which are n = 1, n = 2, and n = 3, see Figure 3. Accordingly, this means that the kinetic energy of the extracellular chloride ion will decrease by q V m n , while the kinetic energy of the intracellular chloride ion will increase by ( 1 1 n ) q V m .
Therefore, the kinetic energy of the intracellular chloride ions can be calculated by the following equation:
K E ( C l I ) = K E T h e r m a l + q V m ( 1 1 n )
In addition, the kinetic energy of the extracellular chloride ions can be calculated by the following equation:
K E ( C l E ) = K E T h e r m a l q V m n
where K E T h e r m a l is the average thermal kinetic energy, q is the charge of chloride ion (1.6 × 10−19 C), Vm is the cellular membrane potential, and n is the location of the gate.
If we take the average thermal energy K E T h e r m a l = 3 2 K B T = 0.64 × 10 20 J for extracellular ions, this will result in a negative kinetic energy. For example, if n = 1 and Vm = 0.07 V, then K E ( C l E ) = 0.48 × 10 20 J. Such a negative value can be used and substituted in Equation (1) to calculate the tunneling probability. However, to get a reasonable positive value of the kinetic energy of the extracellular ions, the following equation can be used:
K E ( C l E ) = 0 + ( 3 K B T q V m n ) 2 = 3 2 K B T q V m 2 n
Equation (10) gives the average kinetic energy of the extracellular chloride ions of the two values, which are 0 J and 2 K E T h e r m a l q V m n = 2 ( 3 2 K B T ) q V m n = 3 K B T q V m n . The first value 0 J means that the membrane potential decreased the thermal energy of the chloride ion to zero once it reached the closed gate, while we choose 2 K E T h e r m a l as an estimation of the highest kinetic energy possible with a reasonable fraction of chloride ions since increasing the thermal kinetic energy further will make the fraction of chloride ions much less to be reasonable to calculate the average kinetic energy. As we said, this approach is just to provide a more meaningful physical values of the kinetic energy and not to provide the exact values.
Accordingly, the kinetic energy of the intracellular and the extracellular chloride ions can be calculated by the following equations, respectively:
K E ( C l I ) = 3 2 K B T + q V m ( 1 1 n )
K E ( C l E ) = 3 2 K B T q V m 2 n
where KB = 1.38 × 10−23 J/K and T is the body temperature (310 K).
Equation (12) calculates the average kinetic energy of the extracellular chloride ions that can reach the closed gate. In such way, if n = 1 and Vm = 0.07 V, K E ( C l E ) = 3 2 K B T q V m 2 n = 0.08 × 10 20 J. This value of kinetic energy is reasonable because it is a positive value and it will be less than the energy of the intracellular chloride ions if the same variables are substituted in Equation (11), which is K E ( C l I ) = 0.64 × 10 20 J. Moreover, if Vm = 0 V, both intracellular and extracellular chloride ions have the same average kinetic energy 3 2 K B T . In conclusion, Equations (11) and (12) are simple, but they are consistent and reasonable because they describe how the membrane potential changes the kinetic energy of the intracellular and extracellular chloride ions with meaningful physical values, also in the absence of membrane potential, both of them have the same average kinetic energy. This is said because choosing values higher or less than 2KEThermal will result in a discrepancy between the intracellular and extracellular chloride ions when the case is Vm = 0 V, which is not consistent because both of them should have the same average thermal kinetic energy in such a case.
Accordingly, the quantum tunneling probability of the intracellular chloride ions can be calculated by the following equation:
T Q ( C l I ) = e 8 m ħ × 2 w 3 g ( g ( q V m ( 1 1 n ) + 3 2 K B T ) ) 3
while the quantum tunneling probability of the extracellular chloride ions can be calculated by the following equation:
T Q ( C l E ) = e 8 m ħ × 2 w 3 g ( g ( 3 2 K B T q V m 2 n ) ) 3
To simplify Equations (13) and (14), the following mathematical simplifications will be considered: we take 10−20 as a common factor from the square root and multiply the exponent by ( 10 20 ) 3 = 10 30 ( 10 20 ) 3 = 10 30 , we take L = w 10 10 and multiply the exponent by 10−10 and we take G = g 10 20 and divide the exponent by 10−20. Therefore, when we say for example that the energy of the gate G = 2 J, then the actual energy of the gate g = 2 × 10−20 J. Also, when we say that L = 1 m, then the actual gate length w = 1 × 10−10 m. Eventually, Equations (13) and (14) become:
T Q ( C l I ) = e M ( C l ) × L G ( G 16 V m ( 1 1 n ) 0.64 ) 3
T Q ( C l E ) = e M ( C l ) × L G ( G 0.64 + 8 V m n ) 3
where M(Cl) is just a factor that contains the mass of the ion, the reduced Planck constant and the multiplying and dividing factors that are the result of the simplifications made above.
Hence, M ( C l ) = 8 m ħ × 2 × 10 10 × 10 30 3 × 10 20 = 8 × 5.9 × 10 26 1.05 × 10 34 × 2 × 10 20 3 = 43.62 .
The quantum tunneling of ions can result in a quantum unitary conductance of the closed GABA receptors. To obtain the quantum unitary conductance of the closed GABA receptors, the following equation will be applied [24,25]:
C Q G A B A = q 2 h T Q ( C l )
where q is the charge of chloride ion (1.6 × 10−19 C), h is the Planck constant (6.6 × 10−34 Js), and TQ(Cl) is the quantum tunneling probability of chloride ions.
Eventually, the quantum membrane conductance at certain density D of GABA receptors can be calculated by the following equation [26]:
M C Q G A B A = D × C Q G A B A
To assess the influence of quantum membrane conductance induced by the quantum tunneling of chloride ions on the resting membrane conductance, the Goldman–Hodgkin–Katz (GHK) equation should be used. Before going through the effect of the quantum tunneling of chloride ions, we should first determine the membrane potential by sodium and potassium ions at the physiological parameters. To achieve this, the classical version of the GHK should be used [26]:
[ K ] E M C K + [ N a ] E M C N a = e F V m R T ( [ K ] I M C K + [ N a ] I M C N a )
where [K]E is the extracellular potassium concentration, [K]I is the intracellular potassium concentration, [Na]E is the extracellular sodium concentration, [Na]I is the intracellular sodium concentration, MCK is the leaky membrane conductance of potassium ions at the resting, MCNa is the leaky membrane conductance of sodium ions at the resting state, F is Faraday’s constant (96,485.33 C/mol), R is the gas constant (8.31 J/Kmol), T is the body temperature (310 K), and Vm is the membrane potential.
We added the minus sign in this term e F V m R T to get an absolute value of the membrane potential, which is negative inside with regard to outside. Hence, we will deal with positive values of membrane potential in the following sections but it is considered to be negative inside with regard to outside.
In our paper, we set the following values: [K]E = 4 mmol/L [26], [Na]E = 142 mmol/L [26], [K]I = 140 mmol/L [26], [Na]I = 14 mmol/L [26], MCK = 0.5 mS/cm2 [26], and MCNa = 0.022 mS/cm2 [26] in Equation (19) to get a membrane potential of 0.07 V. These values will be the reference values of investigation in the present study to show how the quantum tunneling of chloride ions can change the membrane potential value.
To assess the effect of the classical transport of chloride ions through the open GABA receptors, the following equation can be used [26]:
[ K ] E M C K + [ N a ] E M C N a + [ C l ] I M C C G A B A = e F V m R T ( [ K ] I M C K + [ N a ] I M C N a + [ C l ] E M C C G A B A )
where [Cl]I is the intracellular chloride concentration, [Cl]E is the extracellular chloride concentration, and MCC−GABA is the membrane conductance of chloride ions due to the classical opening of the GABA receptor.
On the other hand, to show how the quantum tunneling of chloride ions through the closed GABA receptors affect the resting membrane potential, the quantum version of the GHK equation is used [27]:
[ K ] E M C K + [ N a ] E M C N a + [ C l ] I M C Q G A B A ( I ) = e F V m R T ( [ K ] I M C K + [ N a ] I M C N a + [ C l ] E M C Q G A B A ( E ) )
where M C Q G A B A ( I ) is the quantum membrane conductance of the intracellular chloride ions and M C Q G A B A ( E ) is the quantum membrane conductance of the extracellular chloride ions.
The values of chloride concentration are [Cl]E = 120 mmol/L [26] and [Cl]I = 5 mmol/L [26].

3. Results

In this section, we will evaluate the values of quantum tunneling probability of chloride ions through the closed GABA receptors, the values of quantum conductance of GABA receptors, and the values of quantum membrane conductance of GABA receptors. Eventually, we will be able to investigate the effect of quantum tunneling of chloride ions on the resting membrane potential.
The work that GABA agonist performs to activate the GABA receptors is called the energy of activation. The energy of activation is a reasonable estimation for the energy barrier height of the closed gate of GABA receptors because agonists such as GABA induce conformational changes to open the closed gate, hence what prevents chloride ion from passing this closed gate is the energy required to open it. The activation energy required to open the closed GABA receptors by GABA ligand is 4.3 kcal/mol = 3 × 10−20 J [28]. However, a more accurate estimation of the barrier height should be made based on the calculations of potential mean force (PMF) while ion passes through the closed gate. Up to authors’ knowledge, there is no evidence on the exact value of the barrier height for the closed gate of GABA receptor, but as GABA receptors are similar to the other pentameric ligand-gated ion channels such as nACh receptors and 5-HT3 receptors, then the obtained value of 10KBT = 4.3 × 10−20 J from the closed nACh receptors can be used to estimate the barrier height of the closed gate for GABA receptor [17,18]. In both ways of estimation, the energy barrier is within the order of magnitude 10−20. Accordingly, we are going to choose the range (2 − 6) × 10−20 J for the barrier height g with an average value 4 × 10−20 J to evaluate its influence on the quantum tunneling of chloride ions. The most important aspect of this range is that its values are within the order of magnitude of 10−20, which is consistent with the experimental observations. Moreover, we will show later that the influential quantum tunneling, which affects the membrane potential, occurs at lower values and around 2 × 10−20 J because the quantum model explains the depolarization effect of GABA receptors based on the drop in the height of the energy barrier possibly due to either physiological or pathological conditions. Hence, any value higher than 2 × 10−20 J will not be significant to enhance the quantum tunneling of chloride ions. This indicates that the energy barrier made by the hydrophobic gate is high enough to block the classical permeation and the quantum tunneling at normal conditions in which the barrier is higher than 2 × 10−20 J.
Additionally, the closed gate is formed by the aggregation of amino-acid residues. Hence, it is expected that the length of the gate equals the length of a residue within the alpha-helix 1.5 × 10−10 m [29,30]. In our study, the length range (1 − 2) × 10−10 m is considered for the purposes of evaluation with an average value of 1.5 × 10−10 m. Moreover, we will show how the length of the gate can influence the quantum tunneling-induced membrane depolarization by setting values higher than 1.5 × 10−10 m in case each subunit of GABA receptor contributes by more than one residue to the hydrophobic gate.
The range of gate location n (1–3) is considered for the purposes of evaluation with an average value of 2 (at the middle of the GABA receptor) to account for the different possible locations of the gate in the GABA receptor. Moreover, the membrane potential range 0–0.07 V is considered for the purposes of investigation.

3.1. The Quantum Tunneling Probability of Chloride Ions through the Closed GABA Receptors

The quantum tunneling probability of intracellular chloride ions and extracellular chloride ions can be calculated by the following equations, respectively:
T Q ( C l I ) = e 43.62 × L G ( G 16 V m ( 1 1 n ) 0.64 ) 3
T Q ( C l E ) = e 43.62 × L G ( G 0.64 + 8 V m n ) 3  
Based on Equations (22) and (23), the relationship between the energy of the gate and the tunneling probability of chloride ions is represented at different values of gate length to show its influence on this relationship. See Figure 4.
Based on Equations (22) and (23), the relationship between the energy of the gate and the tunneling probability of chloride ions is represented at different values of membrane potential to show its influence on this relationship. See Figure 5.
Based on Equations (22) and (23), the relationship between the energy of the gate and the tunneling probability of chloride ions is represented at different values of gate location to show its influence on this relationship. See Figure 6.

3.2. The Quantum Unitary Conductance of Closed GABA Receptors Mediated by the Quantum Tunneling of Chloride Ions

The quantum unitary conductance of closed GABA receptors for the intracellular chloride ions and the extracellular chloride ions can be calculated by the following equations, respectively:
C Q G A B A ( I ) = 3.88 × 10 5 × e 43.62 × L G ( G 16 V m ( 1 1 n ) 0.64 ) 3
C Q G A B A ( E ) = 3.88 × 10 5 × 43.62 × L G ( G 0.64 + 8 V m n ) 3
where q 2 h = ( 1.6 × 10 19 ) 2 6.6 × 10 34 = 3.88 × 10 5 S. The unit of quantum unitary conductance used in the present study is Siemens (S).
Based on Equations (24) and (25), the relationship between the energy of the gate and the quantum unitary conductance is represented at different values of gate length. See Figure 7.
Based on Equations (24) and (25), the relationship between the energy of the gate and the quantum unitary conductance is represented at different values of membrane potential. See Figure 8.
Based on Equations (24) and (25), the relationship between the energy of the gate and the quantum unitary conductance is represented at different values of gate location. See Figure 9.

3.3. The Quantum Membrane Conductance from the Closed GABA Receptors

The quantum membrane conductance from the closed GABA receptors for the intracellular chloride ions and the extracellular chloride ions can be calculated by the following equations, respectively:
M C Q G A B A ( I ) = 3.88 × 10 2 × e 43.62 × L G ( G 16 V m ( 1 1 n ) 0.64 ) 3 × D
M C Q G A B A ( E ) = 3.88 × 10 2 × 43.62 × L G ( G 0.64 + 8 V m n ) 3 × D
where 3.88 × 10−5 S is converted to 3.88 × 10−2 mS by multiplying by 103 and D is the GABA receptors density (channels/cm2). Hence, the unit of the quantum membrane conductance used in the present study is mS/cm2. We will use three different values of GABA receptors density [31] to observe its influence on quantum conductance. These three values are: 10 receptors/µm2 (109 channels/cm2), 100 receptors/µm2 (1010 channels/cm2), and 1000 receptors/ µm2 (1011 channels/cm2).
Based on Equations (26) and (27), the relationship between the energy of the gate and the quantum membrane conductance is represented at different values of gate length. See Figure 10.
Based on Equations (26) and (27), the relationship between the energy of the gate and the quantum membrane conductance is represented at different values of membrane potential. See Figure 11.
Based on Equations (26) and (27), the relationship between the energy of the gate and the quantum membrane conductance is represented at different values of gate location. See Figure 12.
Based on Equations (26) and (27), the relationship between the energy of the gate and the quantum membrane conductance is represented at different values of receptors density. See Figure 13.

3.4. The Quantum Tunneling-Induced Membrane Depolarization via Quantum Tunneling of Chloride Ions through the Closed GABA Receptors

We substitute the values of ions concentration and their conductance values in Equation (21) to get the following equation:
5.12 + 0.19 × D × e 43.62 × L G ( G 16 V m ( 1 1 n ) 0.64 ) 3 = e 37.45 V m ( 70.31 + 4.66 × D × e 43.62 × L G ( G 0.64 + 8 V m n ) 3 )
Equation (28) enables us to assess the influence of the quantum tunneling of chloride ions according to several factors including: the energy of the gate, the length of the gate, and the location of the gate.
Based on Equation (28), we can investigate the influence of gate energy on the resting membrane potential at different values of gate length. See Figure 14.
Based on Equation (28), we can investigate the influence of the gate energy on the resting membrane potential at different values of gate location. See Figure 15.
Based on Equation (28), we can investigate the influence of the gate energy on the resting membrane potential at different values of receptors density. See Figure 16.

3.5. The Effect of the Intracellular Chloride Concentration on the Resting Membrane Potential Based on the Quantum and Classical Models

To investigate the influence of the intracellular chloride concentration in the context of the quantum tunneling of chloride ions, see Figure 17.
The graphs in Figure 17 are plotted according to Equation (28) by setting different values of intracellular chloride concentration.
To investigate the influence of the intracellular chloride concentration on the resting membrane potential via the classical opening of certain numbers of GABA receptors, the concentration of ions and their conductance values are substituted in Equation (20) to get the following equation:
5.12 + C C G A B A × D × [ C l ] I = e 37.45 V m ( 70.31 + C C G A B A × D × [ C l ] E )
where CC−GABA is the classical unitary conductance of open GABA receptor (30 × 10−12 S [32]) and D is the density of GABA receptors. Moreover, C C G A B A × D = M C C G A B A is the equation of the membrane conductance due to the opening of GABA receptors.
According to Equation (29), the relationship between the intracellular chloride concentration and the resting membrane potential can be investigated at different values of membrane conductance from the opening of GABA receptors MCC−GABA. See Figure 18.

4. Discussion

The present work introduces the mechanism of membrane depolarization by GABA receptors based on the quantum tunneling model. The basic idea behind this mechanism is that chloride ions can pass through the closed gate of GABA receptors via quantum tunneling. Quantum tunneling is a probabilistic quantum event; hence, the quantum tunneling probability should be calculated. To assess the influence of chloride ions quantum tunneling on the resting membrane conductance, the unitary quantum conductance of the GABA receptor and the quantum membrane conductance should be calculated. The following sections will focus on the numerical descriptions of these quantum electrical properties of GABA receptors and chloride ions to provide a better understanding of the depolarization action. The quantum tunneling model is proposed because it can explain the depolarization action of GABA receptors without the requirement of the intracellular chloride ions accumulation as in the classical interpretation, see Figure 19.
The increase in the intracellular chloride concentration may be overestimated because the measurement is based on the expression levels of chloride transporters and the experimental studies did not take into consideration the factor of neuronal activity, which affects the intracellular chloride concentration [8]. Therefore, the quantum model seems helpful to explain the depolarization effect in cases where the intracellular chloride concentration is normal.
Based on Figure 4, Figure 5 and Figure 6:
  • It is clear that the tunneling probability of intracellular chloride ions is higher than that of extracellular chloride ions. This is attributed to the difference in the kinetic energy when they reach the closed gate. The higher the kinetic energy, the higher the tunneling probability. As the membrane potential is negative inside with regard to outside, it is expected that the intracellular chloride ion acquires kinetic energy while passing in GABA receptor until hitting the closed gate, while the extracellular chloride ions will lose kinetic energy as they pass in GABA receptors until reaching the closed gate. See Figure 20.
The case here is the opposite if it is compared with the case of cations such as lithium, sodium, and potassium ions [9,13] because the extracellular cations will get kinetic energy higher than the intracellular cations. Interestingly, such discrepancy in terms of tunneling probability generates what is called ‘quantum gradient’ that favors the efflux of chloride ions, which results in membrane depolarization, as will be discussed later.
2.
The tunneling probability of both intracellular and extracellular chloride ions decreases as the energy of the gate and the length of the gate increase. On the other hand, the membrane potential has a differential effect on chloride ions. As the membrane potential increases, the tunneling probability of intracellular chloride ions increases, whereas the tunneling probability of extracellular chloride ions decreases as the membrane potential increases. Additionally, as the value of n increases, the membrane potential available for intracellular chloride ions increases and thus their tunneling probability increases. However, as the value of n increases, the membrane potential opposing the movement of extracellular chloride ions and the loss of their kinetic energy decrease and thus their tunneling probability increases.
The binding of GABA ligand to its receptor is required to open the closed gate by decreasing its energy barrier [28]. When the gate is open, the chloride ion can pass classically. The quantum tunneling model implies that chloride ions can pass through the closed gate at different levels of gate energy and without binding of GABA to its receptor. The values of tunneling probability at different levels of gate energy were clearly presented in the previous figures according to different factors. To determine the significance of this quantum tunneling, it is vital to consider the quantum unitary conductance of GABA receptor and to be compared with the classical conductance when the receptor is open.
Based on Figure 7, Figure 8 and Figure 9:
  • According to the quantum tunneling model, the quantum unitary conductance of the intracellular chloride ions is higher than the quantum conductance of the extracellular chloride ions. On the other hand, the classical model assigns the same unitary conductance for both extracellular and intracellular chloride ions.
  • If we take the classical unitary conductance of open GABA receptor to be 30 × 10−12 S [32], it seems obvious that the quantum unitary conductance is extremely small until reaching near G = 2 J, at which the quantum unitary conductance becomes comparable to the classical value. Therefore, it is expected that the quantum tunneling of chloride ions becomes influential when the energy of the gate of GABA receptors decreases to 2 J or less, as will be explained later in the paper. Based on the quantum tunneling model, it is predicted that GABA receptors have persistent basal activity in the absence of GABA ligands under certain conditions that can be physiological or pathological if they cause a decrease in the energy of the closed gate. The quantum tunneling effect will be extremely low at gate energy values higher than 2 J.
  • The quantum unitary conductance changes according to the energy of the gate, the length of the gate, the membrane potential, and the location of the gate as it was discussed regarding the tunneling probability.
Accordingly, the quantum membrane conductance due to quantum tunneling of chloride ions through the closed GABA receptors can be investigated. Thus, based on Figure 10, Figure 11, Figure 12 and Figure 13:
  • The quantum membrane conductance of intracellular chloride ions is higher than the quantum membrane conductance of the extracellular chloride ions. This discrepancy in the quantum conductance will result in a net efflux of chloride ions.
  • If the quantum conductance values are compared with the leaky classical conductance values of potassium and sodium ions at the resting state, which are 0.5 mS/cm2 and 0.022 mS/cm2, then it is clear that the values of quantum conductance for the intracellular ions become comparable to the classical ones once the energy of the closed gate drops to 2 J. This is said to provide an empirical prediction regarding the influence of quantum tunneling on the resting membrane potential. As the quantum conductance of the extracellular ions remain extremely low across all different values of gate energy and the quantum conductance of the intracellular ions becomes comparable to the leaky conductance of sodium and potassium ions around G = 2 J, then it is expected that chloride ions can depolarize the resting membrane potential via quantum tunneling through the closed GABA receptors.

4.1. The Influence of Chloride Ions Quantum Tunneling on the Resting Membrane Potential

4.1.1. The Influence of Quantum Tunneling of Chloride Ions at Different Values of Gate Length L

Equation (28) is used to investigate the influence of chloride ions quantum tunneling on the resting membrane potential. According to the quantum tunneling model, it is expected that chloride ions will depolarize the resting membrane potential because intracellular chloride ions have higher quantum membrane conductance if they are compared with the extracellular ions. This will result in a net chloride efflux that depolarizes the membrane potential. To make the depolarization effect of GABA receptors more understandable, a numerical description will be useful.
Based on Figure 14, Figure 15 and Figure 16, it is clear that the membrane potential remains unaffected across a certain range of gate energy until reaching a gate energy value called the energy curving point Gcur at which the action of depolarization begins. To calculate this value, we choose Vm = 0.086 V to be the beginning value of depolarization then we substitute it in Equation (28) and solve the equation for G to get Gcur based on the setting values in Figure 14, Figure 15 and Figure 16. By looking at the mathematical graphs, it is clear that the depolarization starts near the value of G = 2 J. The end of the graphs is at G = 1 J. Hence, the rate of depolarization can be calculated by the following equation:
R = 0.086 V m ( G = 1 ) G c u r 1  
We used the value G = 1 J and we did not use it before in the discussion of tunneling probability and quantum conductance to ensure that the energy barrier of the closed gate G is higher than the kinetic energy of the ion so that the principle of quantum tunneling can be applied; otherwise, the ion can pass classically through the gate. Moreover, Equation (28) is more dynamical equation if it is compared with the equations of tunneling probability and equations of quantum conductance because it take into consideration the interactions of two factors, which are the membrane potential and the energy of the gate, hence longer range of gate energy can be used.
To provide a numerical description of the depolarization action of GABA receptors at different values of gate length, see Table 1.
As the value of G c u r increases, this means that the GABA receptor is more sensitive to the drop in the energy of the gate to depolarize the resting membrane potential because it will be more able to depolarize the membrane potential at higher values of G.
Based on Table 1, we can conclude that as the gate length increases, the closed GABA receptor becomes less sensitive to depolarize the membrane potential. This is indicated by the decreasing values of Gcur. On the other hand, the rate of depolarization increases as the gate length increases to compensate for the decreasing sensitivity.
Our mathematical investigation in the present paper focuses on the range (1 − 2) × 10−10 m for the gate length. However, it will be useful to challenge the quantum model at higher values of the gate length. See Figure 21.
Figure 21 shows that the quantum tunneling-induced membrane depolarization is still evident and valid even at higher values of gate length in case there is more than one hydrophobic residue from each subunit of the GABA receptor that contributes to the hydrophobic gating. As the gate length L increases, the energy curving point decreases and the absolute value of membrane potential at G = 1 J increases (less depolarization). This means that as the gate length increases, the range of gate energy, in which the quantum tunneling can affect the membrane potential gets narrower. Herein, we stress on the idea that the kinetic energy of the ion is less than the energy of the gate across this range of influential quantum tunneling that can depolarize the membrane potential.

4.1.2. The Influence of Quantum Tunneling of Chloride Ions at Different Values of Gate Location n

To provide a numerical description of the depolarization of GABA receptors at different values of gate location n, see Table 2.
Based on Table 2, we can conclude that as the value of gate location n increases, the closed GABA receptor becomes more sensitive to depolarize the resting membrane potential indicated by the increasing values of Gcur. On the other hand, the rate of depolarization decreases as the n value increases.

4.1.3. The Influence of Quantum Tunneling of Chloride Ions on the Resting Membrane Potential at Different Values of Channels Density D

To provide a numerical description of the depolarization by GABA receptors at different values of receptors density D, see Table 3.
Based on Table 3, we can conclude that as the receptors density increases, the closed GABA receptor becomes more sensitive to depolarize the resting membrane potential, which is indicated by the increasing values of Gcur. On the other hand, the rate of depolarization decreases as the density of receptors increases.

4.2. The Comparison between the Classical and Quantum Models of GABA Receptors in Terms of the Intracellular Chloride Concentration

According to Figure 17, it is obvious that the depolarization action mediated by the quantum tunneling of chloride ions is affected minimally by increasing the intracellular chloride concentration. Interestingly, even at intracellular chloride concentration of 5 mmol/L, closed GABA receptors can depolarize the resting membrane potential starting from G near 2 J until G = 1 J with a specific rate of depolarization. At the gate energy values higher than 2 J, the depolarization action via quantum tunneling is lost even at high intracellular chloride concentration because the tunneling probability will be extremely low to get significant quantum conductance to change the membrane potential, as we discussed before.
On the other hand, according to Figure 18, the classical model predicts that depolarization action starts around the intracellular chloride concentration of 9 mmol/L. At intracellular concentrations less than 9 mmol/L, hyperpolarization is induced by GABA receptors opening because the absolute value of Nernst potential for chloride ions is higher than the absolute value of the initial membrane potential, which is 70 mV. Therefore, a net chloride influx occurs and manifests as hyperpolarization.
Moreover, the quantum model is mainly concerned with the energy of the closed gate more than the chloride concentrations, while the classical model mainly depends on the intracellular concentration to explain the depolarization action. Therefore, the same intracellular chloride concentration can have different values of depolarized membrane potential according to the value of gate energy between 2 J and 1 J based on the quantum tunneling model.

4.3. Possible Proposed Mechanisms behind the Drop in the Energy of the Closed Gate of GABA Receptors and Potential Implications and Perspectives

The major aim of the present article is to explain the depolarization action of GABA receptors according to the drop in the energy of the gate. This, as presented before, enables us to comprehend this effect under normal intracellular chloride concentration. In this section, we will propose certain mechanisms that can be responsible for the drop in the energy of the gate. These mechanisms may be used to provide more concrete approaches to test the quantum tunneling model of chloride ions.
The classical view of the GABA receptor function relies on the ability of the agonist, such as the GABA ligand, to open the closed gate by performing a work equal to the activation energy [28]. This gate is said to be open because its energy barrier is less than the kinetic energy of chloride ions, hence it allows the passage of both intracellular and extracellular chloride ions by the same permeability. The agonist decreases the energy barrier of the closed gate by dilating the pore and increasing the hydration probability [17,18]. Furthermore, the net flow is determined by the resting membrane potential and the Nernst potential of chloride ions, which is modulated by the intracellular chloride ions. On the other hand, the quantum tunneling model determines the function of GABA receptor according to the energy of the closed gate.
Accordingly, if certain physiological and pathological mechanisms decrease the energy of the closed gate, the quantum tunneling of chloride ions can be observed. However, the decrease in the energy of the gate must be made so that the gate is still closed and the energy barrier of the gate is higher than the kinetic energy of the chloride ion. According to the presented results, the critical energy value at which the depolarization effect starts is around 2 J. The major advantage of the quantum model over the classical model is that the explanation of the depolarization action of GABA receptors does not require the action of the intracellular accumulation of chloride ions. Therefore, the quantum model will be useful to explain the events of depolarization in which there is no chloride ions accumulation especially that the previously conducted studies might overestimate the values of the intracellular chloride concentration and did not take into consideration the neuronal activity that can influence the chloride concentrations. Thus, an alternative mechanism is required to explain such conditions of depolarization.
Our aim by listing the following mechanisms is to show how the quantum tunneling model can be useful and how it can be exploited to be put under experimental investigation. Here we list a few possible mechanisms with potential implications:
  • We propose that during early development, the binding of GABA ligand to its receptor does not yield the full activation energy needed to open the closed gate; hence, the gate remains closed but with less energy barrier than before. Accordingly, the quantum tunneling of the intracellular chloride ions will be augmented more significantly than the tunneling of the extracellular chloride ions. This yields a quantum gradient that produces a net chloride efflux that leads to membrane depolarization, which has been observed during early neuronal development [2,3,4].
  • The gain-of-function mutations are well-established genetic conditions that affect ion channels such as ligand-gated channels and voltage-gated channels [33,34]. If these gain-of-function mutations decrease the gate energy of the closed GABA receptor, this will enhance the quantum tunneling of chloride ions to depolarize the membrane potential. Interestingly, this quantum tunneling-induced membrane depolarization occurs even without GABA binding to the receptor because the drop of the gate energy had occurred. Therefore, these mutations may be responsible for generating outward anion current that depolarizes membrane potential, which explains the increased tendency of seizures in patients with such mutations [33,34].
  • Acquired conditions such as stroke and trauma are implicated as pathological events that affect the function of GABA receptors [5,6]. As such conditions can decrease the energy required to open the voltage-gated channels [35,36], this can be extrapolated and applied on GABA receptors because they share the same hydrophobic gating mechanism [14,17]. Accordingly, these conditions can decrease the energy of the closed gate in GABA receptors, which leads to enhancing the quantum tunneling of chloride ions.
  • Since the depolarization mediated by GABA receptors can result in excitatory effects or inhibitory effects [8], this idea can be used to explain the mechanism of action of certain classes of agents and medications. These agents include the anesthetic agents that are used to shut down the consciousness so that certain surgery can be done. The classical view is that these agents activate the GABA receptors to induce hyperpolarizing currents and thus decrease neuronal activity [37,38]. In the perspective of the quantum tunneling model, these agents can hypothetically be viewed as quantum tunneling enhancers by lowering the values of gate energy facilitating the chloride efflux, which leads to quantum tunneling-induced membrane depolarization. This depolarization may inhibit neuronal activity, which explains the anesthetic effects of these agents. Moreover, the quantum tunneling model predicts the ability of chloride ions to depolarize the membrane potential by a large degree, reaching less than 10 mV (negative inside with regard to outside), as is clearly evident in Figure 14, Figure 15, Figure 16 and Figure 17. This large depolarization significantly increases the fraction of the inactivated sodium channels blocking the neuronal activity even though that membrane potential is depolarized, and this can be served to explain the inhibitory effects of anesthetic agents on consciousness.
  • Interestingly, as the quantum tunneling model predicts the ability of GABA receptor to permeate chloride ions without the binding of GABA ligand at certain range of gate energy, this comes consistent with the spontaneous activity of GABA receptors without the requirement of GABA binding, which has been observed and mentioned in the literature [39].
These proposed mechanisms require experimental verification to prove the scientific validity of the quantum tunneling model of chloride ions. Hence, the present study represents the first stage of the investigation process towards better understanding of the biological function of GABA receptors from the quantum mechanical perspective and towards more established quantum biology and quantum medicine.

5. Conclusions

The quantum tunneling model of chloride ions can explain the depolarization effect of GABA receptors at normal intracellular chloride ions without the requirement of the intracellular chloride ions accumulations. The quantum tunneling-induced membrane depolarization has been found to start at around G = 2 J.

Author Contributions

Conceptualization, A.B.Q.; methodology, A.B.Q. and S.N.; software, A.B.Q.; validation, A.B.Q., A.S. and S.N.; formal analysis, A.B.Q.; investigation, S.N., A.B.Q. and A.S.; resources, S.N., A.B.Q. and A.S.; data curation, A.B.Q., writing—original draft preparation, A.B.Q.; writing—review and editing, A.B.Q., S.N., A.S., O.A., F.M.Z., M.A.O.A.-A. and M.B.A.; visualization, A.S. and A.B.Q.; supervision, S.N.; project administration, S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Spitzer, N.C. How GABA generates depolarization. J. Physiol. 2010, 588, 757. [Google Scholar] [CrossRef] [PubMed]
  2. Ben-Ari, Y.E.; Cherubini, E.N.; Corradetti, R.E.; Gaiarsa, J. Giant synaptic potentials in immature rat CA3 hippocampal neurons. J. Physiol. 1989, 416, 303–325. [Google Scholar] [CrossRef] [PubMed]
  3. Kasyanov, A.M.; Safiulina, V.F.; Voronin, L.L.; Cherubini, E. GABA-mediated giant depolarizing potentials as coincidence detectors for enhancing synaptic efficacy in the developing hippocampus. Proc. Natl. Acad. Sci. USA 2004, 101, 3967–3972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Fritschy, J.M. Epilepsy, E/I balance and GABAA receptor plasticity. Front. Mol. Neurosci. 2008, 1, 5. [Google Scholar] [CrossRef] [Green Version]
  5. Jaenisch, N.; Witte, O.W.; Frahm, C. Downregulation of potassium chloride cotransporter KCC2 after transient focal cerebral ischemia. Stroke 2010, 41, e151–e159. [Google Scholar] [CrossRef]
  6. Dzhala, V.; Valeeva, G.; Glykys, J.; Khazipov, R.; Staley, K. Traumatic alterations in GABA signaling disrupt hippocampal network activity in the developing brain. J. Neurosci. 2012, 32, 4017–4031. [Google Scholar] [CrossRef]
  7. Herbison, A.E.; Moenter, S.M. Depolarising and hyperpolarising actions of GABAA receptor activation on GnRH neurons: Towards an emerging consensus. J. Neuroendocrinol. 2011, 23, 557. [Google Scholar] [CrossRef] [Green Version]
  8. Kilb, W. When Are Depolarizing GABAergic Responses Excitatory? Front. Mol. Neurosci. 2021, 14, 747835. [Google Scholar] [CrossRef]
  9. Khreesha, L.; Qaswal, A.B.; Al Omari, B.; Albliwi, M.A.; Ababneh, O.; Albanna, A.; Abunab’ah, A.; Iswaid, M.; Alarood, S.; Guzu, H.; et al. Quantum Tunneling-Induced Membrane Depolarization Can Explain the Cellular Effects Mediated by Lithium: Mathematical Modeling and Hypothesis. Membranes 2021, 11, 851. [Google Scholar] [CrossRef]
  10. Alrabayah, M.; Qaswal, A.B.; Suleiman, A.; Khreesha, L. Role of Potassium Ions Quantum Tunneling in the Pathophysiology of Phantom Limb Pain. Brain Sci. 2020, 10, 241. [Google Scholar] [CrossRef] [Green Version]
  11. Qaswal, A.B.; Suleiman, A.; Guzu, H.; Harb, T.; Atiyat, B. The potential role of lithium as an antiviral agent against SARS-CoV-2 via membrane depolarization: Review and hypothesis. Sci. Pharm. 2021, 89, 11. [Google Scholar] [CrossRef]
  12. Qaswal, A.B. Quantum tunneling of ions through the closed voltage-gated channels of the biological membrane: A mathematical model and implications. Quantum Rep. 2019, 1, 19. [Google Scholar] [CrossRef] [Green Version]
  13. Qaswal, A.B.; Ababneh, O.; Khreesha, L.; Al-Ani, A.; Suleihat, A.; Abbad, M. Mathematical Modeling of Ion Quantum Tunneling Reveals Novel Properties of Voltage-Gated Channels and Quantum Aspects of Their Pathophysiology in Excitability-Related Disorders. Pathophysiology 2021, 28, 10. [Google Scholar] [CrossRef]
  14. Aryal, P.; Sansom, M.S.; Tucker, S.J. Hydrophobic gating in ion channels. J. Mol. Biol. 2015, 427, 121–130. [Google Scholar] [CrossRef] [Green Version]
  15. Bali, M.; Akabas, M.H. The location of a closed channel gate in the GABAA receptor channel. J. Gen. Physiol. 2007, 129, 145–159. [Google Scholar] [CrossRef] [Green Version]
  16. Naffaa, M.M.; Hung, S.; Chebib, M.; Johnston, G.A.; Hanrahan, J.R. GABA-ρ receptors: Distinctive functions and molecular pharmacology. Br. J. Pharmacol. 2017, 174, 1881–1894. [Google Scholar] [CrossRef] [Green Version]
  17. Rao, S.; Klesse, G.; Lynch, C.I.; Tucker, S.J.; Sansom, M.S. Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions. J. Phys. Chem. B 2021, 125, 981–994. [Google Scholar] [CrossRef]
  18. Beckstein, O.; Sansom, M.S. A hydrophobic gate in an ion channel: The closed state of the nicotinic acetylcholine receptor. Phys. Biol. 2006, 3, 147. [Google Scholar] [CrossRef]
  19. Tepper, H.L.; Voth, G.A. Mechanisms of passive ion permeation through lipid bilayers: Insights from simulations. J. Phys. Chem. B 2006, 110, 21327–21337. [Google Scholar] [CrossRef] [Green Version]
  20. Khavrutskii, I.V.; Gorfe, A.A.; Lu, B.; McCammon, J.A. Free energy for the permeation of Na+ and Cl ions and their ion-pair through a zwitterionic dimyristoyl phosphatidylcholine lipid bilayer by umbrella integration with harmonic fourier beads. J. Am. Chem. Soc. 2009, 131, 1706–1716. [Google Scholar] [CrossRef] [Green Version]
  21. Jia, Z.; Yazdani, M.; Zhang, G.; Cui, J.; Chen, J. Hydrophobic gating in BK channels. Nat. Commun. 2018, 9, 3408. [Google Scholar] [CrossRef] [PubMed]
  22. Xin, Y.; Huang, A.; Hu, Q.; Shi, H.; Wang, M.; Xiao, Z.; Zheng, X.; Di, Z.; Chu, P.K. Barrier Reduction of Lithium Ion Tunneling through Graphene with Hybrid Defects: First-Principles Calculations. Adv. Theory Simul. 2018, 1, 1700009. [Google Scholar] [CrossRef]
  23. Poltavsky, I.; Zheng, L.; Mortazavi, M.; Tkatchenko, A. Quantum tunneling of thermal protons through pristine graphene. J. Chem. Phys. 2018, 148, 204707. [Google Scholar] [CrossRef] [PubMed]
  24. Serway, R.A.; Moses, C.J.; Moyer, C.A. Modern Physics; Cengage Learning: Boston, MA, USA, 2004. [Google Scholar]
  25. Chen, F.; Hihath, J.; Huang, Z.; Li, X.; Tao, N.J. Measurement of single-molecule conductance. Annu. Rev. Phys. Chem. 2007, 58, 535–564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Bertil, H.; Bertil, H. Ion Channels of Excitable Membranes, 3rd ed.; Sinauer Associates: Sunderland, MA, USA, 2007; Volume 2001. [Google Scholar]
  27. Qaswal, A.B. Quantum Electrochemical Equilibrium: Quantum Version of the Goldman–Hodgkin–Katz Equation. Quantum Rep. 2020, 2, 17. [Google Scholar] [CrossRef]
  28. Jones, M.V.; Jonas, P.; Sahara, Y.; Westbrook, G.L. Microscopic kinetics and energetics distinguish GABAA receptor agonists from antagonists. Biophys. J. 2001, 81, 2660–2670. [Google Scholar] [CrossRef] [Green Version]
  29. Moran, L.A.; Horton, R.A.; Scrimgeour, K.G.; Perry, M.D. Principles of Biochemistry; Pearson: London, UK, 2014. [Google Scholar]
  30. Lodish, H.; Berk, A.; Kaiser, C.A.; Krieger, M.; Scott, M.P.; Bretscher, A.; Ploegh, H.; Matsudaira, P. Molecular Cell Biology; Macmillan: London, UK, 2008. [Google Scholar]
  31. Martín-Belmonte, A.; Aguado, C.; Alfaro-Ruíz, R.; Moreno-Martínez, A.E.; De La Ossa, L.; Martínez-Hernández, J.; Buisson, A.; Shigemoto, R.; Fukazawa, Y.; Luján, R. Density of GABAB receptors is reduced in granule cells of the hippocampus in a mouse model of Alzheimer’s disease. Int. J. Mol. Sci. 2020, 21, 2459. [Google Scholar] [CrossRef] [Green Version]
  32. Eghbali, M.; Birnir, B.; Gage, P.W. Conductance of GABAA channels activated by pentobarbitone in hippocampal neurons from newborn rats. J. Physiol. 2003, 552, 13–22. [Google Scholar] [CrossRef]
  33. Ahring, P.K.; Liao, V.W.; Gardella, E.; Johannesen, K.M.; Krey, I.; Selmer, K.K.; Stadheim, B.F.; Davis, H.; Peinhardt, C.; Koko, M.; et al. Gain-of-function variants in GABRD reveal a novel pathway for neurodevelopmental disorders and epilepsy. Brain J. Neurol. 2021, 144, awab391. [Google Scholar] [CrossRef]
  34. Crunelli, V.; Di Giovanni, G.; Parri, H.R.; Errington, A.C. Gain-of-function of thalamic extrasynaptic GABA-A receptors in typical absence seizures. In Extrasynaptic GABAA Receptors; Springer: New York, NY, USA, 2014; pp. 223–237. [Google Scholar]
  35. Morris, C.E. Voltage-gated channel mechanosensitivity: Fact or friction? Front. Physiol. 2011, 2, 25. [Google Scholar] [CrossRef] [Green Version]
  36. Reeves, D.; Ursell, T.; Sens, P.; Kondev, J.; Phillips, R. Membrane mechanics as a probe of ion-channel gating mechanisms. Phys. Rev. E 2008, 78, 041901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wolfe, R.C. Inhaled anesthetic agents. J. PeriAnesthesia Nurs. 2020, 35, 441–444. [Google Scholar] [CrossRef] [PubMed]
  38. Son, Y. Molecular mechanisms of general anesthesia. Korean J. Anesthesiol. 2010, 59, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. O’Neill, N.; Sylantyev, S. The functional role of spontaneously opening GABAA receptors in neural transmission. Front. Mol. Neurosci. 2019, 12, 72. [Google Scholar] [CrossRef]
Figure 1. The figure represents a schematic diagram of two possible shapes of the energy barrier of the gate with respect to the ion’s position x through the gate. The green curve represents a triangle energy barrier in which the energy increases linearly till the end of the gate reaching the maximum energy value g similar to the energy barrier of the regular electric field, while the red curve represents a triangle barrier that shows a symmetry at the center similar to what observed experimentally in the hydrophobic materials. The red curve describes a linear increase in the energy of the gate till reaching the center where it reaches the maximum value g and then it starts to decrease by the same rate till the end of the gate. x1 is where the kinetic energy of the ion equals the barrier energy at the green curve, while x1′ and x2′ are the ion’s positions where the kinetic energy of the ion equals the barrier energy on the red curve.
Figure 1. The figure represents a schematic diagram of two possible shapes of the energy barrier of the gate with respect to the ion’s position x through the gate. The green curve represents a triangle energy barrier in which the energy increases linearly till the end of the gate reaching the maximum energy value g similar to the energy barrier of the regular electric field, while the red curve represents a triangle barrier that shows a symmetry at the center similar to what observed experimentally in the hydrophobic materials. The red curve describes a linear increase in the energy of the gate till reaching the center where it reaches the maximum value g and then it starts to decrease by the same rate till the end of the gate. x1 is where the kinetic energy of the ion equals the barrier energy at the green curve, while x1′ and x2′ are the ion’s positions where the kinetic energy of the ion equals the barrier energy on the red curve.
Cells 11 01145 g001
Figure 2. The figure represents a schematic diagram of chloride ion movement in the GABA receptor. The intracellular chloride ion moves from the negative side of the membrane towards the positive side, and this increases its kinetic energy until hitting the closed gate in the center of the GABA receptor, while the extracellular chloride ion moves from the positive side towards the negative side and this decreases its kinetic energy until hitting the closed gate.
Figure 2. The figure represents a schematic diagram of chloride ion movement in the GABA receptor. The intracellular chloride ion moves from the negative side of the membrane towards the positive side, and this increases its kinetic energy until hitting the closed gate in the center of the GABA receptor, while the extracellular chloride ion moves from the positive side towards the negative side and this decreases its kinetic energy until hitting the closed gate.
Cells 11 01145 g002
Figure 3. The figure represents different locations of the closed gate in the GABA receptor. The gate location affects the kinetic energy of the chloride ion once it reaches the closed gate. As the n value increases, the membrane potential available to enhance the kinetic energy of the intracellular anion increases, hence its kinetic energy increases. On the other hand, as the n value increases, the membrane potential available to dampen the kinetic energy of the extracellular anion decreases, hence the decrease in its kinetic energy is reduced.
Figure 3. The figure represents different locations of the closed gate in the GABA receptor. The gate location affects the kinetic energy of the chloride ion once it reaches the closed gate. As the n value increases, the membrane potential available to enhance the kinetic energy of the intracellular anion increases, hence its kinetic energy increases. On the other hand, as the n value increases, the membrane potential available to dampen the kinetic energy of the extracellular anion decreases, hence the decrease in its kinetic energy is reduced.
Cells 11 01145 g003
Figure 4. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the tunneling probability of the intracellular chloride ions log 10 ( T Q ) I and the tunneling probability of the extracellular chloride ions log 10 ( T Q ) E , respectively. The relationships are investigated at different values of gate length and according to the setting values above each figure.
Figure 4. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the tunneling probability of the intracellular chloride ions log 10 ( T Q ) I and the tunneling probability of the extracellular chloride ions log 10 ( T Q ) E , respectively. The relationships are investigated at different values of gate length and according to the setting values above each figure.
Cells 11 01145 g004
Figure 5. (a,b): The figure represents the relationship between the energy of gate G and the common logarithm of the tunneling probability of the intracellular chloride ions log 10 ( T Q ) I and the tunneling probability of the extracellular chloride ions log 10 ( T Q ) E , respectively. The relationships are investigated at different values of membrane potential and according to the setting values indicated above each graph.
Figure 5. (a,b): The figure represents the relationship between the energy of gate G and the common logarithm of the tunneling probability of the intracellular chloride ions log 10 ( T Q ) I and the tunneling probability of the extracellular chloride ions log 10 ( T Q ) E , respectively. The relationships are investigated at different values of membrane potential and according to the setting values indicated above each graph.
Cells 11 01145 g005
Figure 6. (a,b): The figure represents the relationship between the energy of gate G and the common logarithm of the tunneling probability of the intracellular chloride ions log 10 ( T Q ) I and the tunneling probability of the extracellular chloride ions log 10 ( T Q ) E , respectively. The relationships are investigated at different values of gate location and according to the setting values indicated above each graph.
Figure 6. (a,b): The figure represents the relationship between the energy of gate G and the common logarithm of the tunneling probability of the intracellular chloride ions log 10 ( T Q ) I and the tunneling probability of the extracellular chloride ions log 10 ( T Q ) E , respectively. The relationships are investigated at different values of gate location and according to the setting values indicated above each graph.
Cells 11 01145 g006
Figure 7. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum unitary conductance of the intracellular chloride ions log 10 ( C Q G A B A ) I and the quantum unitary conductance of the extracellular chloride ions log 10 ( C Q G A B A ) E , respectively. The relationships are investigated at different values of gate length and according to the setting values above each graph.
Figure 7. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum unitary conductance of the intracellular chloride ions log 10 ( C Q G A B A ) I and the quantum unitary conductance of the extracellular chloride ions log 10 ( C Q G A B A ) E , respectively. The relationships are investigated at different values of gate length and according to the setting values above each graph.
Cells 11 01145 g007
Figure 8. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum unitary conductance of the intracellular chloride ions log 10 ( C Q G A B A ) I and the quantum unitary conductance of the extracellular chloride ions log 10 ( C Q G A B A ) E , respectively. The relationships are investigated at different values of membrane potential and according to the setting values above each graph.
Figure 8. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum unitary conductance of the intracellular chloride ions log 10 ( C Q G A B A ) I and the quantum unitary conductance of the extracellular chloride ions log 10 ( C Q G A B A ) E , respectively. The relationships are investigated at different values of membrane potential and according to the setting values above each graph.
Cells 11 01145 g008
Figure 9. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum unitary conductance of the intracellular chloride ions log 10 ( C Q G A B A ) I and the quantum unitary conductance of the extracellular chloride ions log 10 ( C Q G A B A ) E , respectively. The relationships are investigated at different values of gate location and according to the setting values above each graph.
Figure 9. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum unitary conductance of the intracellular chloride ions log 10 ( C Q G A B A ) I and the quantum unitary conductance of the extracellular chloride ions log 10 ( C Q G A B A ) E , respectively. The relationships are investigated at different values of gate location and according to the setting values above each graph.
Cells 11 01145 g009
Figure 10. (a,b): The figure represents the relationship between the energy of the gate and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of gate length and according to the setting values above each graph.
Figure 10. (a,b): The figure represents the relationship between the energy of the gate and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of gate length and according to the setting values above each graph.
Cells 11 01145 g010
Figure 11. (a,b): The figure represents the relationship between the energy of the gate and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of membrane potential and according to the setting values above each graph.
Figure 11. (a,b): The figure represents the relationship between the energy of the gate and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of membrane potential and according to the setting values above each graph.
Cells 11 01145 g011
Figure 12. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of gate location and according to the setting values above each graph.
Figure 12. (a,b): The figure represents the relationship between the energy of the gate G and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of gate location and according to the setting values above each graph.
Cells 11 01145 g012
Figure 13. (a,b): The figure represents the relationship between the energy of gate G and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of receptors density and according to the setting values above each graph.
Figure 13. (a,b): The figure represents the relationship between the energy of gate G and the common logarithm of the quantum membrane conductance of the intracellular chloride ions log 10 ( M C Q G A B A ) I and the quantum membrane conductance of the extracellular chloride ions log 10 ( M C Q G A B A ) E , respectively. The relationships are investigated at different values of receptors density and according to the setting values above each graph.
Cells 11 01145 g013
Figure 14. The figure represents the influence of the energy of gate G on the resting membrane potential according to different values of gate length and according to the setting values above the figure. The observed depolarization effect is due to quantum tunneling of chloride ions through GABA receptors.
Figure 14. The figure represents the influence of the energy of gate G on the resting membrane potential according to different values of gate length and according to the setting values above the figure. The observed depolarization effect is due to quantum tunneling of chloride ions through GABA receptors.
Cells 11 01145 g014
Figure 15. The figure represents the influence of the energy of gate G on the resting membrane potential according to different values of gate location and according to the setting values above the figure.
Figure 15. The figure represents the influence of the energy of gate G on the resting membrane potential according to different values of gate location and according to the setting values above the figure.
Cells 11 01145 g015
Figure 16. The figure represents the influence of the energy of gate G on the resting membrane potential according to different values of receptors density and according to the setting values above the figure.
Figure 16. The figure represents the influence of the energy of gate G on the resting membrane potential according to different values of receptors density and according to the setting values above the figure.
Cells 11 01145 g016
Figure 17. The figure represents the relationship between the energy of gate G and the resting membrane potential under the influence of the quantum tunneling of chloride ions at different intracellular chloride concentrations.
Figure 17. The figure represents the relationship between the energy of gate G and the resting membrane potential under the influence of the quantum tunneling of chloride ions at different intracellular chloride concentrations.
Cells 11 01145 g017
Figure 18. The figure represents the relationship between the intracellular chloride ions and the resting membrane potential under the influence of the classical transport of chloride ions through open GABA receptors at different values of membrane conductance of chloride ions.
Figure 18. The figure represents the relationship between the intracellular chloride ions and the resting membrane potential under the influence of the classical transport of chloride ions through open GABA receptors at different values of membrane conductance of chloride ions.
Cells 11 01145 g018
Figure 19. The figure represents the classical view that explains the depolarization action of GABA receptor. (a): The decrease in the intracellular chloride concentration will favor the influx of chloride ions that tends to hyperpolarize the membrane potential when the GABA receptor opens. (b): The increase in the concentration of the intracellular chloride ions will favor the efflux of chloride ions that tends to depolarize the resting membrane potential.
Figure 19. The figure represents the classical view that explains the depolarization action of GABA receptor. (a): The decrease in the intracellular chloride concentration will favor the influx of chloride ions that tends to hyperpolarize the membrane potential when the GABA receptor opens. (b): The increase in the concentration of the intracellular chloride ions will favor the efflux of chloride ions that tends to depolarize the resting membrane potential.
Cells 11 01145 g019
Figure 20. The figure represents the quantum mechanical view of the depolarization action of GABA receptor. (a): The intracellular chloride ion with higher kinetic energy, which is represented by a higher quantum wave frequency, has a higher quantum tunneling probability through the closed gate, which is represented by a higher wave amplitude after passing the gate. (b): The extracellular chloride ion with a lower kinetic energy, which is represented by a lower quantum wave frequency, has a lower quantum tunneling probability through the closed gate, which is represented by a lower wave amplitude after passing the gate. (c): Due to the discrepancy between the intracellular and extracellular chloride ions in terms of the tunneling probability, a quantum gradient favors the efflux of chloride ions, which will depolarize the membrane potential. The quantum efflux gradient is represented by a dotted line.
Figure 20. The figure represents the quantum mechanical view of the depolarization action of GABA receptor. (a): The intracellular chloride ion with higher kinetic energy, which is represented by a higher quantum wave frequency, has a higher quantum tunneling probability through the closed gate, which is represented by a higher wave amplitude after passing the gate. (b): The extracellular chloride ion with a lower kinetic energy, which is represented by a lower quantum wave frequency, has a lower quantum tunneling probability through the closed gate, which is represented by a lower wave amplitude after passing the gate. (c): Due to the discrepancy between the intracellular and extracellular chloride ions in terms of the tunneling probability, a quantum gradient favors the efflux of chloride ions, which will depolarize the membrane potential. The quantum efflux gradient is represented by a dotted line.
Cells 11 01145 g020
Figure 21. The figure represents the relationship between the energy of gate G and the membrane potential at gate length L values higher than 2 m to show its influence on the depolarization effect.
Figure 21. The figure represents the relationship between the energy of gate G and the membrane potential at gate length L values higher than 2 m to show its influence on the depolarization effect.
Cells 11 01145 g021
Table 1. The table represents the curving energy point, the membrane potential at G = 1 J, and the average rate of depolarization at different values of gate length L and according to the other setting values in Figure 14.
Table 1. The table represents the curving energy point, the membrane potential at G = 1 J, and the average rate of depolarization at different values of gate length L and according to the other setting values in Figure 14.
L = 1 m
IonGcurVm at G = 1 JR
Cl2.35 J0.0064 V0.046 V/J
L = 1.5 m
IonGcurVm at G = 1 JR
Cl1.98 J0.0044 V0.066 V/J
L = 2 m
IonGcurVm at G = 1 JR
Cl1.8 J0.0039 V0.081 V/J
Table 2. The table represents the curving energy point, the membrane potential at G = 1 J, and the average rate of depolarization at different values of gate location n and according to the other setting values in Figure 15.
Table 2. The table represents the curving energy point, the membrane potential at G = 1 J, and the average rate of depolarization at different values of gate location n and according to the other setting values in Figure 15.
n = 1
IonGcurVm at G = 1 JR
Cl1.2 J 0.0062 V0.31 V/J
n = 2
IonGcurVm at G = 1 JR
Cl1.98 J 0.0044 V0.066 V/J
n = 3
IonGcurVm at G = 1 JR
Cl2.23 J 0.004 V0.053 V/J
Table 3. The table represents the curving energy point, the membrane potential at G = 1 J, and the average rate of depolarization at different values of channels density D and according to the other setting values in Figure 16.
Table 3. The table represents the curving energy point, the membrane potential at G = 1 J, and the average rate of depolarization at different values of channels density D and according to the other setting values in Figure 16.
D = 109 channels/cm2
IonGcurVm at G = 1 JR
Cl1.91 J 0.0044 V0.071 V/J
D = 1010 channels/cm2
IonGcurVm at G = 1 JR
Cl1.98 J 0.0044 V0.066 V/J
D = 1011 channels/cm2
IonGcurVm at G = 1 JR
Cl2.05 J 0.0044 V0.062 V/J
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nawafleh, S.; Qaswal, A.B.; Suleiman, A.; Alali, O.; Zayed, F.M.; Al-Adwan, M.A.O.; Bani Ali, M. GABA Receptors Can Depolarize the Neuronal Membrane Potential via Quantum Tunneling of Chloride Ions: A Quantum Mathematical Study. Cells 2022, 11, 1145. https://doi.org/10.3390/cells11071145

AMA Style

Nawafleh S, Qaswal AB, Suleiman A, Alali O, Zayed FM, Al-Adwan MAO, Bani Ali M. GABA Receptors Can Depolarize the Neuronal Membrane Potential via Quantum Tunneling of Chloride Ions: A Quantum Mathematical Study. Cells. 2022; 11(7):1145. https://doi.org/10.3390/cells11071145

Chicago/Turabian Style

Nawafleh, Sager, Abdallah Barjas Qaswal, Aiman Suleiman, Obada Alali, Fuad Mohammed Zayed, Mohammad Abu Orabi Al-Adwan, and Mo’ath Bani Ali. 2022. "GABA Receptors Can Depolarize the Neuronal Membrane Potential via Quantum Tunneling of Chloride Ions: A Quantum Mathematical Study" Cells 11, no. 7: 1145. https://doi.org/10.3390/cells11071145

APA Style

Nawafleh, S., Qaswal, A. B., Suleiman, A., Alali, O., Zayed, F. M., Al-Adwan, M. A. O., & Bani Ali, M. (2022). GABA Receptors Can Depolarize the Neuronal Membrane Potential via Quantum Tunneling of Chloride Ions: A Quantum Mathematical Study. Cells, 11(7), 1145. https://doi.org/10.3390/cells11071145

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop