1. Introduction
The electrolyte-oxide-semiconductor field-effect transistor (EOSFET) has been introduced as an ion sensitive sensor (ISFET), where the interaction of ions with the gate oxide induces a changed source-drain current at a constant gate voltage applied to a reference electrode in the bath electrolyte [
1]. A more specific interface-induced response can be achieved for biological analytes when the gate oxide is modified by suitable coatings [
2]. However, response of an EOSFET may also be induced by ionic currents in the electrolyte, which give rise to a changed electrical potential beyond the electrical double layer of the oxide. The present paper describes the concepts of current-induced transistor sensorics for electroactive cells with a focus on planar cell-substrate junctions. The situation is illustrated in
Figure 1 [
3]: the adherent cell membrane and the semiconductor chip are separated by a thin cleft of extracellular electrolyte. Ionic and capacitive currents through the membrane give rise to a current flow along the cleft. That ionic current induces a voltage
with respect to the bath and modulates the electronic source-drain current of a transistor at constant gate voltage.
Three issues are addressed in particular: (i) the cell-substrate junction is described as a planar capacitive core-coat conductor, which allows relating the cellular currents with changes of the electrical potential and changes of ion concentrations in the junction; (ii) in certain situations, the current-induced response is superposed by the interface-induced response of an ion-sensitive gate due to changed ion concentrations; and (iii) the inhomogeneous distribution of ion conductances in the cell membrane is crucial for a current-induced transistor response with intact electrogenic cells.
This paper is a scholarly overview of concepts and does not intend to provide a comprehensive review of the field, though some inconsistencies in the literature are addressed. The examples for various kinds of transistor recordings with electroactive cells are chosen from studies in my own lab.
The structure of the paper is as follows: First, the physics of the EOSFET is considered in order to identify the current-induced and the interface-induced transistor response with electrogenic cells. Then, the theory of a planar capacitive core-coat conductor is presented, which relates cellular currents with changes of the electrical potential and of ion concentrations in a cell-substrate junction. The dynamics of a core-coat conductor is described, and the reduction to a one-compartment model is derived for constant and for changing ion concentrations. On that basis, the transistor recording of voltage-activated ion channels is considered for an intracellular voltage-clamp. In this context, the role of an ion-sensitive gate oxide and the issue of inhomogeneous membrane conductances are pointed out. Current-induced transistor sensorics with intact cells is described for neuronal action potentials recorded at soma and axon, and for a receptor-cell-transistor biosensor.
Several technical issues have to be envisaged for cell-based transistor sensorics: (i) The gate must be smaller than a cell, and far smaller than a traditional transistor sensorics gate. As a consequence, the noise power, which is reciprocal to the gate area, becomes a problem; (ii) Current-induced extracellular voltages are in the range of a few microvolts to a few millivolts, lower than usual interface-induced signals. The signal-to-noise ratio is a critical issue; (iii) The transistor signals vary due to a variability of the cells, to the variable interactions of cells and substrate, and to a variable position of the transistor in the cell-substrate junction. These problems are not addressed here.
2. Transistor Sensorics with Cellular Currents
First some general aspects of the EOSFET are considered to avoid misunderstandings. On that basis, the interface-induced and current-induced responses of transistor recording are identified for electrogenic cells.
Gate voltage and electrical potential. The source–drain current of an EOSFET is determined by band bending, the drop of the electrical potential in the semiconductor
from its surface at the gate oxide to the bulk [
4]. Band bending is controlled by a gate voltage
that is applied between the wires joined to a reference electrode (usually Ag/AgCl) in the bulk electrolyte and to the backside metal contact of the bulk semiconductor. (Practically, a voltage
is applied to the semiconductor with the electrode on ground potential.) We consider the relation of
with the difference
of the electrical potentials in the bulk electrolyte and in the bulk semiconductor, and then the relation of that difference with the band bending
.
The metal wires (identical materials) are in electronical equilibrium with the metal contacts (different materials) of the electrolyte and of the semiconductor, respectively, and these contacts are in electronical equilibrium with the redox solution of the reference electrode and with the bulk semiconductor, respectively. These equilibrations are described as alignments of the electrochemical potential of electrons (Fermi energies)
(
elementary charge) in the different phases with components of an electronic chemical potential and the electrical potential. (For the quasi-Fermi energy of a redox solution see [
5].) As a consequence, the gate voltage controls the difference
of the electronic electrochemical potentials in the redox solution and in the semiconductor [
6]. When we take into account the electrical equilibration of the redox solution and the bulk electrolyte through a salt-bridge with
, the relation between the gate voltage and the difference of the electrical potentials in the bulk electrolyte and the bulk semiconductor is expressed by Equation (1).
The electronic chemical potentials in the redox solution and in the semiconductor are important for the value of the flatband and the threshold voltage of an EOSFET (see
Appendix A), but irrelevant for changes of the gate voltage around a working voltage of the field-effect transistor above the threshold of strong inversion with
.
The overall difference
of the electrical potentials can be partitioned into several components according to Equation (2a) with the potential drops
and
in the electrolyte and in the semiconductor (both from the surface at the oxide to the bulk), with the potential jumps
and
across the interfaces oxide-to-electrolyte and oxide-to-semiconductor, and with the potential drop
across the oxide (electrolyte-to-semiconductor). The latter is split into contributions of a fixed charge
at the oxide/semiconductor interface and of the space charge in the semiconductor according to Equation (2b) for a p-channel transistor in depletion approximation (gate capacitance
, permittivity
, donor concentration
) [
4].
The fabrication of a chip determines the potential drop
as well as the potential jump
. The conditions of the electrolyte determine the potential jump
(Si=O and Si-OH groups, oriented H
2O, double layer of Si-O
− groups and adsorbed ions) [
7]. The term quasi-Helmholtz layer is used in analogy to a metal/electrolyte interface with
. The conditions of the electrolyte also determine the potential drop
with a component
in a diffuse electrical double layer (Gouy-Chapman) and a component
in the free electrolyte beyond the double layer, which is induced by a current flow in the electrolyte.
vanishes in equilibrium for the electroneutral electrolyte. Considering Equations (1) and (2a), a changed gate voltage
is partitioned into a changed current-induced drop of the electrical potential in the free electrolyte, a changed potential drop across the oxide/electrolyte interface (with Gouy-Chapman and quasi-Helmholtz layer), and a changed band bending (with the changed potential drop
due the space charge) according to Equation (3).
Transistor and sensor. Equation (3) can be read in two ways: (i) In a transistor mode, the conditions in the electrolyte are constant. A changed gate voltage
induces a changed band bending with a changed source-drain current above the threshold of strong inversion; (ii) In a sensor mode, the gate voltage is held constant. A changed electrical potential in the free electrolyte or across the Gouy-Chapman and quasi-Helmholtz layers gives rise to a changed band bending. The resulting change of the source-drain current can be expressed as an apparent change of the gate voltage with Equation (4).
The interface-induced response accounts for the classical ISFET sensorics, whereas the current-induced response is a typical issue of the transistor sensorics with electrogenic cells.
Interpretation. The goal of transistor sensorics with electrogenic cells is an interpretation of a signal in terms of the generating cellular currents. According to the core-coat conductor theory of a cell-substrate junction, as considered in the next section, cellular currents induce changes of the electrical potential in the junction and changes of the ion concentrations, which give rise to an interface-induced response .
The interpretation of a signal
is involved if the interface-induced response plays a role. In a first step, the interface-induced response must be identified in the total transistor signal. This task is difficult because the peculiar (slow) dynamics of ion concentrations affects the interface-induced response as well as the current-induced response. In a second step, a calibration in terms of concentration changes
is required. In principle, such a relation is obtained by titrating a cell-free transistor with the ion concentration
in the bath electrolyte yielding an ion sensitivity of the gate per decade
according to Equation (5) for a limited concentration range.
However, the chemical nature of the gate surface is affected by cell adhesion. With an unknown ion sensitivity in the junction, an assignment of a changing concentration is not possible.
An unequivocal interpretation of transistor signals with electroactive cells is achieved if the interface–induced response plays no role with . Two situations can be envisaged: (i) the activation of cellular ion currents is so short that changes of the ion concentrations in the junction can be neglected; and (ii) the gate is chemically passivated such that changing ion concentrations are not able to induce an interface-induced response.
Finally, it may be noted that the sensor concept as sketched here for the EOSFET is also applicable for an electrolyte-oxide-metal-oxide-semiconductor (EOMOS) transistor and –
mutatis mutandis – for other potentiometric electrodes (see
Appendix B).
3. Planar Capacitive Core-Coat Conductor
Ionic and capacitive currents through an adherent cell membrane induce changes of the electrical potential and of the ion concentrations in a cell-substrate junction. These relations are described in terms of a planar capacitive core-coat conductor.
The thin film of extracellular electrolyte in a cell-substrate junction is insulated from the intracellular electrolyte by the plasma membrane and from the semiconductor by an oxide layer (
Figure 1b). For cultured mammalian cells, the distance of the two insulating coats is
as measured by fluorescence interference contrast (FLIC) microscopy [
8]. It is far larger than the Debye length of the electrolyte with
and far smaller than the width of the adhesion area with a radius in the range of
. The condition
is similar to a relation for neuronal axons, where a thin thread of intracellular electrolyte is insulated from the bath by the cylindrical plasma membrane with
for Debye length, radius, and length.
The electrical features of an axon are characterized by three properties: (i) the cytoplasmatic core as well as the bath are electroneutral beyond thin space charge layers at the membrane; (ii) the electrical charge in these space charge layers determines the voltage across the membrane by a capacitor relation; and (iii) the profiles of the ion concentrations and of the electrical potential across the narrow core are neglected as compared with the profiles along the core. As a consequence, axons are described as 1D linear capacitive core-coat conductors with ionic currents through the coat and along the core, usually for constant intracellular ion concentration [
9] but occasionally also for changing ion concentrations [
10]. With the same rationale, cell-substrate junctions have been described as 2D planar capacitive core-coat conductors for constant [
11] and for changing ion concentrations [
12,
13].
Balance of ion flux. The ion flux through the adherent membrane in a cell-substrate junction gives rise to a changing (in time) ion concentration in the core and to a changing (in space) ion flux along the core. The flux along the core is driven by gradients (in two dimensions) of the electrical potential
(drift mobility
) and of the concentration
(diffusion coefficient
, thermal energy
). The flux balance is described by Equation (6) with the ionic charge number
and the electrical current
per unit area of the adherent membrane [
14].
Balance of charge. Changing ion concentrations give rise to a changing electrical charge in the junction. A net charge
in the core electrolyte would violate the electroneutrality condition. That charge is located in thin space charge layers along membrane and semiconductor oxide, respectively, with counter-charges induced across the two coats [
15]. A changing charge is related with changing differences
and
of the electrical potential in the junction with respect to cell and semiconductor by a differential capacitor relation. With a capacitance
of the membrane and an average capacitance
of the substrate (differing from the local capacitance
of the gate oxide), we obtain Equation (7). (A different approach has been chosen in [
16], where the net charge is assigned to a space charge in the electrolyte film. See
Appendix C).
Current balance and limit of constant ion concentrations. The balance of electrical charge (Equation (7)) can be combined with the dynamics of the ion concentrations (Equation (6)). We obtain Equation (8) for the balance of the total current, when the electrical potentials in junction, cell, and semiconductor are expressed as voltages with respect to the bulk electrolyte with
,
, and
defining a junction capacitance
(Distinguish capacitance
and ion concentrations
!).
For small changes of the ion concentrations, the diffusion term in Equation (8) can be cancelled and the concentrations can be excluded from the spatial derivative of the drift term. With a conductivity
of the electrolyte in the junction, we obtain Equation (9).
This “2D cable equation” for a planar capacitive core-coat conductor with constant ion concentrations is homologous to the common 1D cable equation for neuronal axons. Two peculiarities, however, must be noted: (i) the voltage in the junction depends on the voltage changes in the two different adjacent media of cell and substrate; and (ii) at its periphery, the core conductor is open to the bath electrolyte.
Transistor response. The capacitive and ionic currents through the adherent membrane of a cell-substrate junction give rise to a coupled dynamics of the extracellular ion concentrations
and of the extracellular electrical potential
, as described by Equations (6) and (7). The apparent change of the gate voltage
for an EOSFET in the junction (Equation (4)) is caused by a current-induced response due to a changing voltage in the junction and by an interface-induced response due to changing ion concentrations according to Equation (10).
4. Dynamics of Planar Capacitive Core-Coat Conductor
The time constants for changes of the voltage and of the ion concentrations in a cell-substrate junction are crucial for the response of a transistor to cellular currents. We consider a circular junction (radius , area ) with a homogeneous sheet resistance .
Voltage dynamics. The dynamics of the junction voltage
for constant ion concentrations is described by Equation (11a) with a radial coordinate
(at a constant bias voltage of the semiconductor). The ionic membrane currents depend on the conductances per unit area
according to Equation (11b) , which are driven by the transmembrane voltage and the reversal voltage
(determined by the difference of the logarithmic concentrations across the membrane).
First, we apply a voltage step
to a cell without membrane conductances for a boundary condition
and an initial condition
. The sharp pulse of capacitive current induces a junction voltage
, which decays according to Equation (12) with a relaxation of modes (Bessel functions
,
, zeros
of
) [
17] (p. 199).
To reveal the response to an ion current, we apply a voltage step to a cell with a single ion conductance
starting from
. Besides the capacitive response, an increasing junction voltage is induced according to Equation (13) for
with a stationary parabolic profile [
17] (p. 204).
The complete response is described by a superposition of the Equations (12) and (13). The dynamics is characterized by the time constant of the lowest mode.
Typical parameters for cultured mammalian cells (HEK293, hippocampal neuron) are
,
,
,
, and
. The time-dependent profiles
upon a voltage step
are plotted in
Figure 2 for the capacitive current and for the ionic current [
18]. The capacitive current induces an initial voltage of 77
mV. The ionic current gives rise to a stationary parabolic voltage profile with an amplitude of 2.5
mV. The dynamics is dominated by the time constant
.
Quasi-
stationary response. The electrical time constant of a cell-substrate junction is in the microsecond range and far shorter than the time scale of the natural dynamics of capacitive and ionic currents of cells in the range of 0.1–1
ms. As a consequence, the junction voltage is almost relaxed upon stimulation. The 2D cable equation (Equation (9)) can be reduced to Equation (14) for a homogeneous junction and small voltages
.
In this quasi-stationary situation, the planar capacitive core-coat conductor appears as a sheet conductor that is driven from the cell by a 2D (capacitive and ionic) current-source density. There is a certain analogy to a volume conductor with a 3D current-source density [
19,
20]. Considering Equation (13), the quasi-stationary voltage profile is parabolic for a circular junction according to Equation (15).
Concentration dynamics. The concentration dynamics of ions in a planar core-coat conductor is governed by drift and diffusion according to Equation (6). The drift term is determined by the electrical potential, which depends on the concentration dynamics of all ions (Equation (7)). When we disregard the drift, we obtain Equation (16) for the pure, independent ion diffusion in a homogeneous circular junction.
A step of ion current
induces a concentration change according to Equation (17) for a boundary condition
and an initial condition
with the concentration
in the bulk electrolyte. The dynamics is characterized by the time constant
of the lowest mode.
For
(with
and
),
,
, and
, a stationary concentration change up to 18.6
mM is attained with a time constant
. The dynamics of ion concentrations is far slower than the dynamics of the voltage. It is illustrated by
Figure 2b, if the time axis is rescaled by a factor
. A significant concentration change appears only if the duration of the ion current is comparable with the time constant
.
Electrical characterization of cell-substrate junction. The capacitive transient upon stimulation by a voltage-step (Equation (12),
Figure 2a) can be used to determine the sheet resistance of a cell-substrate junction. More reliable, however, is the application of an ac voltage with an amplitude
(angular frequency
). The frequency-dependent voltage profile
in a homogeneous circular junction without membrane conductance can be derived from Equation (11) according to Equation (18) with the modified Bessel function
. This approach has been used for leech neurons with an array of transistors in the cell–substrate junction [
21].
A related approach is the application of an ac voltage
to the capacitor structure of a chip. The response
across the adherent membrane has been probed with a voltage-sensitive fluorescent dye [
22,
23]. The relation of stimulus and response is similar to Equation (18) but more involved due to the coupled response of junction and cell. An average sheet resistance of
has been obtained for HEK293 cells on a silica surface coated with fibronectin and of
for hippocampal rat neurons on poly–lysine [
24].
A particularly elegant method is the recording of the thermal voltage fluctuations in a cell-substrate junction [
25]. The autocorrelation function of
, averaged over a circular sensor with a radius
in the center of a circular junction, is expressed by Equation (19) [
26].
Practically, the power spectral density (PSD) of the voltage fluctuations is probed, and the data are fitted with a theoretical PSD, which is obtained by Fourier transformation of the autocorrelation function. For low frequencies, it is given by Equation (20) [
26].
The PSD for rat neurons on poly–lysine has been probed with an EOSFET at low frequency with
[
25]. According to Equation (20), with
, a sheet resistance of
is obtained, which is similar to the resistance evaluated from the ac stimulation.
The PSD for a wide range of frequencies has been probed for neurons from
Lymnaea stagnalis using an array of EOMOS transistors, as shown in
Figure 3 [
26]. The data are fitted with a sheet resistance
and a junction capacitance (whole cell) of
. The unusually large sheet resistance, which indicates a particularly close adhesion, was due to a special preparation of the dissociated snail neurons. The perfect fit of the data over five orders of magnitude of frequencies confirms that the cell-substrate junction is indeed a capacitive core-coat conductor.
5. One-Compartment Model
The typical issue of a cell-semiconductor assembly is the electrical coupling of the ionic and the electronic system across the planar junction. The profiles of the electrical potential and of the ion concentrations are not so important. For that reason, the planar core-coat conductor can be drastically reduced: drift and diffusion of ions along the differential area elements of the junction are replaced by global drift and diffusion from a single compartment to the bath electrolyte. On the one hand, this simple approach provides a surprisingly good approximation for the dynamics of the planar core-coat conductor, and, on the other hand, it is adequate if a cell-substrate junction is not sufficiently well defined with respect to the profiles of resistance, capacitance, and membrane conductance and to the position of a transistor. The “one-compartment model” is derived from the planar core-coat conductor for constant ion concentrations as well as for changing concentrations.
One-
compartment model for constant ion concentrations. We start with the 2D cable equation (Equation (9)) for a homogeneous circular junction (Equation (11)). We replace the divergence of the electrical current by a global current from a junction compartment to the bath according to Ohm’s law (formally by a substitution
) and obtain Equation (21a).
The Ohmic conductance per unit area as defined by Equation (21b) combines three crucial parameters of a cell-substrate junction, the distance , the area , and the resistivity . The scaling factor is obtained by comparing the one-compartment model with a circular core-coat conductor as follows.
When a voltage step
is applied to a cell without membrane conductance, the capacitive current induces a junction voltage, which decays according to Equation (22a). Starting from
for a single ion conductance
, the voltage step gives rise to an increasing junction voltage according to Equation (22b) for
with a stationary voltage
.
When the time constant is compared with the time constant of the circular core-coat conductor (Equation (12)), we obtain . When the stationary voltage is compared with the stationary parabolic voltage profile (Equation (13)), the scaling factor is or , respectively, if the maximum voltage or the average voltage is considered.
The one-compartment model for constant ion concentrations corresponds to a simple equivalent electrical circuit with a single node assigned to the cell-substrate junction and an Ohmic seal conductance [
27]. The formal equivalent circuit, however, does not provide a relation of the seal conductance with the structural parameters of a cell-substrate junction as it is obtained when the model is derived from the planar core-coat conductor [
28].
One-
compartment model with changing ion concentrations. A one-compartment model can be also defined for changing ion concentrations [
13]. We start with the drift-diffusion equation (Equation (6)) for a homogeneous junction. The divergence of the ion flux is replaced by the global flux from a junction compartment to the bath with Ohm’s law for the drift and with Fick’s first law for the diffusion. We obtain Equation (23a), where the ion conductivity
and reversal voltage
for the adherent membrane depend on the ion concentration according to Equations (23b) [
29].
The voltage
is related with the total electrical charge according to the differential capacitor relation of Equation (24a) (at a constant bias voltage of the semiconductor). For a constant membrane voltage, the integral capacitor relation of Equation (24b) can be used considering an initial condition
and the electroneutrality condition
for the bath electrolyte.
6. Transistor Recording of Ion Conductances
Voltage-activated membrane conductances are crucial for the transistor sensorics of electrogenic cells. As examples, we consider the inactivating sodium conductance and the non-inactivating potassium conductance of the classical Hodgkin-Huxley (HH) model [
30]. They are described by maximum conductances
and
and by the open probabilities of independent gates
for activation and
for inactivation according to Equation (25a). The open probabilities obey rate equations (Equation (25b)) with the rate constants
and
, which depend on the transmembrane voltage.
Quasi-stationary response. The time constant for the dynamics of the gates at common transmembrane voltages is in a range of 0.1–1 ms for the activation of the Na conductance and in a range of 1–10 ms for the inactivation of the Na conductance and the activation of the K conductance. Thus, the dynamics is far slower than the electrical relaxation of the cell-substrate junction but faster than the ion diffusion in the junction. The relation has two consequences: (i) the junction voltage follows the membrane current in a quasi–stationary state; and (ii) the change of the extracellular ion concentrations due to the inactivating Na conductance is small, whereas it may be large for the non-inactivating K conductance if it is held open by a constant intracellular voltage.
To elucidate the transistor response of ion channels, the voltage in an adherent cell is controlled by a whole-cell patch-pipette. With this approach a transistor signal
can be compared with the total current through the cell membrane
such that a current-induced and an interface-induced responses can be distinguished. As a guide, various numerical simulations of the capacitive core-coat conductor are performed with the HH conductances (rate constants at 18.5 °C) using the one-compartment model with changing ion concentrations according to Equations (23) and (24a) with an area
, a distance
, and tabulated diffusion coefficients. A scaling factor
ηJ = 5.78π takes into account the partial spatial averaging by a transistor [
31].
Voltage-clamp with constant ion concentrations. Almost constant ion concentrations in the junction are expected for the Na conductance due to the fast inactivation and the high extracellular Na concentration. This assumption is confirmed by a simulation with (extracellular electrolyte 50 mM·NaCl, 2 mM·KCl, intracellular electrolyte 12 mM Na ions).
Considering Equation (21), the quasi-stationary junction voltage is described by Equation (26a) for
at a constant intracellular voltage. The current through the whole membrane with an area
and an average conductance
is expressed by Equation (26b).
If the dynamics of the conductance is the same in the whole membrane, the voltage
and the current
exhibit the same waveform and are proportional according to Equation (27).
Nav1.4 channel. Transistor recordings under voltage–clamp have been reported for the sodium channel Na
v1.4 in HEK293 cells at an extracellular Na concentration of 50
mM [
32]. The Na conductance was activated by depolarizing steps of the intracellular voltage. The transistor response
and the membrane current
reflected both the typical waveform of a transient sodium current with an activation phase and an inactivation phase. The apparent change of the gate voltage is plotted in
Figure 4 versus the total membrane current. There is a perfect proportionality, which is fitted with a relation
. The result indicates that the transistor signal is indeed due to a pure current-induced response with
. By ac stimulation and recording the response of transistor and current, a parameter
has been measured. Considering Equation (27), a depleted conductance in the adherent membrane with
was estimated.
Voltage-
clamp with changing ion concentration. Changing ion concentrations are expected under voltage-clamp for the non-inactivating K conductance at the low physiological K concentration. The situation is illustrated by a simulation with
(extracellular electrolyte 5
mM·KCl, 135
mM·NaCl, intracellular electrolyte 40
mM K ions). The cell is depolarized by 140
mV for 15
ms starting from the reversal voltage of potassium. The results are plotted in
Figure 5. We have to distinguish three levels: (i) the current-induced response for constant ion concentrations; (ii) the effect of the ion concentrations on the current-induced response; and (iii) the additional effect of the ion concentrations on the interface-induced response.
The junction voltage for constant ion concentrations is plotted in
Figure 5d as trace 0. Upon depolarization, a stationary level of
is reached within a millisecond. The dynamics is determined by the activation of the potassium conductance (
Figure 5b). The stationary voltage disappears suddenly upon repolarization to the reversal voltage.
The junction voltage with changing ion concentrations is shown in
Figure 5d as trace 1. Within a millisecond, it reaches almost the level of
. It relaxes within several milliseconds to a stationary voltage of
. This decay corresponds to the dynamics of the ion concentrations (
Figure 5c). Upon repolarization, the voltage drops suddenly to a negative level and relaxes in a fast phase due to the dynamics of the ion conductance and in a slow phase due to the dynamics of the concentrations.
The K concentration in the junction increases from 5
mM to 15
mM (
Figure 5c). The logarithmic change
may give rise to an interface-induced response (Equation (10)). The total transistor signal
is plotted in
Figure 5e for an insensitive gate (trace 1) with
and for three ion sensitivities
. For a low sensitivity of 1
mV (trace 2), the interface-induced response compensates partially the ion-dependent dynamics of the current-induced response upon depolarization and repolarization. For a sensitivity of 3
mV (trace 3), the transistor signal exhibits two positive components upon depolarization, a fast response, which is determined by the dynamics of the conductance, and a slow response, which is determined by the dynamics of the K concentration. Upon repolarization, the transistor signal drops suddenly due to the disappearance of the membrane current. Then, the interface-induced response relaxes with the ion concentrations. For an ion sensitivity of 12
mV (trace 4), the current-induced and the interface-induced response are fused to a smooth, large signal up to
. The interface-induced component is six times larger than the current–induced response. Upon repolarization, the signal drops suddenly due to the disappearance of the membrane current. Then, the interface-induced response relaxes with
.
Transistor recordings of non–inactivating K conductances under voltage clamp have been reported for hippocampal neurons [
33] and for HEK293 cells with the recombinant channels Maxi-K [
12,
34], Kv1.3 [
13,
35], and EAG [
16].
Kv1.3 channel. The Kv1.3 channel has been studied on chips that were cleaned with an alkaline detergent [
13] and with an acidic detergent [
35]. Recordings of the “acidic experiment” are plotted in
Figure 6a. Upon depolarization, an effective change of the gate voltage up to about 2.5
mV appears within 2
ms. It decays partially during depolarization. Upon repolarization, the signal drops suddenly and a transient negative signal appears. Recordings of the “alkaline experiment” are displayed in
Figure 6b, wherein a smooth increase of the transistor signal up to 7.5
mV is observed. Upon repolarization, the signal drops suddenly by 2
mV, and the remaining amplitude relaxes slowly. (The conditions of the two experiments differ with respect to channel expression, cell size, and transistor quality; the results can be compared only with respect to the dynamics and the relative amplitudes of the fast and slow components).
For the same ion channel, the waveforms of the transistor signals are rather different. When we compare the highest traces of
Figure 6a,b with the simulation of
Figure 5e, it is apparent that the waveform of the “acidic experiment” resembles trace 1 for a pure current
–induced response with
, whereas the waveform of the “alkaline experiment” resembles trace 4 with a superposed interface-induced response for an ion sensitivity 12
mV of the gate. The different ion sensitivities must be related with the cleaning procedures. In fact, after cleaning with the alkaline detergent, the chip surface was highly hydrophilic, whereas it was extremely hydrophobic after cleaning with the acidic detergent due to the deposition (selfassembly) of a cationic amphiphile, which blocks the negative charges of silica.
It appears that rather complex transistor signals can be described in terms of a capacitive core-coat conductor with changing ion concentrations. However, two problems must be considered.
(i) The core-coat conductor theory is not able explain the transistor recordings in terms of the cellular currents with a changing extracellular voltage and changing ion concentrations as long as the ion sensitivity of the gate in the cell-substrate junction is unknown. A calibration of a cell-free transistor is not adequate. For example, the sensitivity
that matches the “alkaline experiment” is distinctly lower than the sensitivity
that was obtained by titration [
13]. The discrepancy must be assigned to a modification of the chip surface by cell adhesion such that negative oxide charges are partially compensated by adsorbed poly-lysine and extracellular matrix protein. (Such a compensation of negative surface charges by a cationic polymer had been observed for a monomolecular lipid film using a fluorescent probe [
36]).
(ii) The pattern of the transistor signals (
Figure 6) is matched quite nicely by the computations (
Figure 5e). However, the interface-induced response is slower than the computed dynamics of the ion concentrations. A low ion mobility in the electrolyte film of the cell-substrate junction can be excluded because the resistivity in the junction is the same as in the bulk electrolyte [
24]. The effect may be due to a slow kinetics of ion binding [
13,
37].
EAG channel. Current-induced and interface-induced signals have been also distinguished in transistor recordings of the EAG potassium channel under voltage-clamp [
16]. The current-induced response was low around 0.2
mV due to a low membrane current. Also the interface-induced response was low due to a low ion sensitivity of 0.7
mV of the gate (maybe caused by a cleaning procedure with an incubation in 20% H
2SO
4). The data were explained in terms of the one-compartment model with constant ion concentrations for the current-induced response and by a Poisson-Nernst-Planck approach to compute the changed ion concentrations, which gave rise to the interface-induced response (see
Appendix C).
Maxi-K channel with high extracellular K concentration. In a first report on transistor recording of a recombinant ion channel (Maxi-K) under voltage clamp [
34], a high extracellular potassium concentration of 100
mM has been applied. The intention was to suppress the interface-induced response, which had been observed when the transistor was cleaned with the alkaline detergent [
12]. The results are displayed in
Figure 7. Upon depolarization, the noisy signals with low time resolution indicate stationary transistor signaly up to 1
mV, which are proportional to the total membrane current with
.
A simulation (extracellular electrolyte 100
mM·KCl, 35
mM·NaCl, intracellular electrolyte 140
mM K ions) with an ion sensitivity of 12
mV for the gate oxide yields a biphasic waveform of the transistor signal that resembles trace 3 in
Figure 5e with a fast current-induced response
and a stationary signal
that is reached within several milliseconds. Thus the high extracellular K concentration lowers drastically the interface-induced component (
cf. Figure 5e), but does not abolish it completely. On the basis of the simulation, a relation
is estimated for the current-induced component of the Maxi-K recording. With a weighting factor
(measured by ac stimulation), an enhancement of the Maxi-K conductance in the adherent membrane of
is obtained using Equation (27) (in comparison to a value of 3.4 estimated in [
34]).
K conductance in hippocampal neurons. A first transistor recording under voltage-clamp with voltage-activated potassium conductances has been reported for cultured hippocampal neurons [
33]. The membrane current and the transistor signal are plotted in
Figure 8. (The chips were cleaned with a mild dish detergent). Upon depolarization the membrane current exhibits an inactivating A-type potassium conductance and a non-inactivating K-type potassium conductance. The transistor response is quite different. The apparent change of the gate voltage increases rapidly by about 80
µV and further more slowly to a stationary level of 120
µV. Upon deactivation, the voltage drops suddenly by about 80
µV with a subsequent slow relaxation. The fast and slow signal were not distinguished in [
33].
It appears that the transistor does not record an A
-type conductance in the junction at all, whereas it responds to the non
–inactivating K
-type conductance with a biphasic waveform that resembles again trace 3 in
Figure 5e. Attributing the fast signal to a current-induced response
and considering the stationary membrane current
, we obtain a relation
. With a weighting factor
(obtained from a capacitive transient), an enhancement of the K-type conductance in the adherent membrane
is obtained using Equation (27) (in comparison to a value of 5.7 estimated in [
33]), whereas the A-type conductance in the adherent membrane is low with
.
7. Transistor Recording of Neuronal Action Potentials
An action potential (AP) in neurons arises from the dynamics of the membrane voltage and of voltage-activated conductances. In the Hodgkin-Huxley model, the AP is an almost monophasic positive transient
of the membrane voltage with a fast rise due to a Na inward current, a slower decay due to a K outward current with an inactivating Na current, and a minor negative tail [
30]. Two issues must be considered when transistors are applied for extracellular recording of APs in an intact cell: (i) the dynamics of the voltage in the closed cell must be taken into account; and (ii) a local compensation of the ionic and capacitive currents in the adherent membrane must be avoided.
Current compensation. Let us consider a free neuron with a homogeneous membrane in a homogeneous electrolyte. The global current balance for the whole cell (membrane area
) with homogeneous ion conductances
is described by Equation (28).
The global current balance implies a local compensation of the ionic and capacitive currents in every area element of the membrane. There is no current flow to the surrounding electrolyte and no extracellular voltage that could be recorded. To obtain an extracellular voltage, the membrane must be inhomogeneous with an enhanced current in one region and a depleted current in another region. For a homogeneous capacitance, that inhomogeneity is created by an inhomogeneous distribution of ion conductances, which is given by the intrinsic structure of a cell or induced by the interaction with the substrate.
Current balance for an adherent cell. We consider a two
–domain model of the cell membrane with ion conductances
in the free membrane (area
) and ion conductances
in the adherent membrane (area
). Considering that the conductances are activated for a short time during a neuronal AP, the effect of changing extracellular ion concentrations is disregarded. We use the one-compartment model for constant ion concentrations according to Equation (21) with a quasi-stationary dynamics of the junction voltage. For small voltages
, we obtain Equation (29a). The current balance of the closed cell is expressed by Equation (29b).
Two simple models are considered for inhomogeneous ion conductances with in the adherent and free membrane.
Leak conductance in soma-
substrate junction. In the first model the cell-substrate junction is formed by part of a neuronal soma with low voltage-activated conductances. The interaction with the substrate is assumed to induce a leak conductance
with
. The junction voltage is determined by the leak current and the capacitive current through the adherent membrane according to Equation (28a). We obtain Equation (30) when we give up the condition
.
For a low leak conductance, the junction voltage is dominated by the capacitive current and proportional to the first derivative of the intracellular voltage. If an almost monophasic AP is created by voltage-activated conductances in the free membrane (e.g., in the axon hillock), a biphasic positive-negative transient is induced (A-type response). If the leak conductance in the junction is large, the voltage transient is proportional to the intracellular voltage, and the response to an AP is monophasic positive (B-type response).
In fact, two kinds of transistor signals had been observed in the original transistor recordings of action potentials for adhering leech neurons [
38] with a small response
resembling
and a large response
resembling
. The interpretation of these recordings in terms of an A-type and B-type response of the leak model was given in a subsequent paper [
39]. It was illustrated by a test experiment.
Defined Gaussian transients
with
were applied to leech neurons using a whole-cell patch-pipette. A-type and B-type responses of the transistor signal
were observed for different cells, as shown in
Figure 9 [
40]. The two waveforms can be parametrized according to Equation (30) for the A-type response with
,
and for the B-type response with an enhanced leak conductance
(responsible for a changed waveform) and a reduced junction conductance
(responsible for a larger amplitude).
Reversible transitions between A-type and B-type recordings of action potentials have been observed when neurons were pressed toward the substrate and when the force was released again [
41,
42]. These transitions were connected with enhanced/reduced leak conductance and a reduced/enhanced junction conductance, respectively. The mechanism is not clear.
Depleted/enhanced voltage-activated conductances. In the second model all voltage-activated conductances in a neuronal membrane are depleted or enhanced in the cell-substrate junction by the same factor
with
. When all ion currents are eliminated from Equations (28a) and (28b), considering Equation (21b), the striking Equation (31) is obtained as a relation between the extracellular and the intracellular voltage [
43].
For complete depletion (), the junction voltage is determined by the capacitive current through the adhesion area . For a homogeneous membrane (), there is no response at all due to local current compensation. For strong enhancement (), the junction voltage is determined by the ionic current through , which appears in Equation (31) as a negative capacitive current through . The first limit corresponds to a soma–substrate junction without leak conductance as considered above (biphasic positive-negative response to an AP). The second limit corresponds to a junction where the local ion currents dominate with a biphasic negative-positive response to an AP.
Dissociated leech neurons have been attached with their axon stump to an array of EOSFETs [
44], and an AP was elicited with an impaled micropipette. A biphasic negative-positive transistor signal was observed, as shown in
Figure 10. The current-induced response
is induced by the large local ion currents in the axon stump. Due to the global current balance in the cell, the signal reflects the inverted capacitive current through the soma membrane. The asymmetry of the biphasic
with a deeper trough corresponds to the asymmetry of an almost monophasic
.
Various waveforms of transistor recordings of APs were reported for cultured hippocampal rat neurons [
3,
45] as well as for dissociated invertebrate neurons from
Lymnaea [
46] and
Aplysia [
47]. They were assigned to a selective enhancement and depletion of different voltage-activated conductances in the adherent membrane.
8. Transistor Recording of Propagating Action Potential
When an action potential propagates along an axon, the different components of the membrane current (capacitive, sodium, potassium) are separated in space. As a consequence, local membrane currents exist, even for a structurally homogeneous axon. An extracellular voltage is induced, which can be probed by transistor recording.
Propagating AP. The membrane current per unit area of an axon is balanced by the changing current along the core according to Equation (32) (cylinder radius
, core resistivity
).
The coupling of the membrane voltage with the voltage
–activated ion conductances gives rise to a wave
that propagates at a velocity
with an invariant waveform. The voltage profile at an arbitrary time is plotted in
Figure 11 (top) as computed with the original Hodgkin-Huxley parameters
,
,
, and
(rate constants at 18.5 °C) for a radius
, typical for an unmyelinated mammalian axon. The velocity is
. The wave is almost monophasic with a steep rise and a smooth drop: The amplitude is 90
mV with a full width at half maximum of 260
µm in space and 300
µs in time. The voltage wave is connected with a current wave
, which is proportional to the curvature of the intracellular voltage in space as well as in time (because the voltage obeys a wave equation) according to Equation (33).
This current wave induces a drop of the electrical potential in the surrounding. For a cable-substrate junction, the voltage with respect to the bulk electrolyte is partitioned into a voltage drop in the junction and a voltage in the surrounding electrolyte with . We consider separately a cable in a homogeneous electrolyte and a cable-substrate junction.
Homogeneous electrolyte. The quasi–stationary response of the electrical potential in a volume conductor (resistivity
) to a current source density
is described by a 3D Poisson equation
[
19]. For a propagating AP, the contributions of all currents along the cylinder are superposed. The extracellular voltage
at a radial distance
is given by Equation (34) when the membrane current per unit length
is placed along the cylinder axis [
48]. (A logarithmic divergence of the integral is avoided due to the finite width of the current wave). A reaction of the extracellular voltage is neglected with
.
The current
is obtained from the Hodgkin-Huxley model for
. Equation (34) is evaluated with a resistivity
typical for a culture medium. The profile of the extracellular voltage at the surface of the axon
is plotted in
Figure 11 (center). The wave is almost biphasic positive-negative and distinctly narrower than the monophasic wave of the intracellular voltage. Peak and trough are separated by 120
µm in space and 140
µs in time. The amplitudes are 6
µV and −9
µV [
49].
Cable-substrate junction. When a cylinder is attached to a planar substrate, a lane-shaped junction is formed with a distance
that increases from a minimum
at the midline towards the cylinder radius
. We replace that profile by a planar junction with a distance
and match its resistance by choosing a half-width
, where the actual distance would be
. This planar junction is driven by the current
. The quasi-stationary response of the electrical potential is determined by a 2D Poisson equation
(Equation (14)). The curvature along the lane is neglected because the membrane current changes smoothly along the axon. The voltage drop
with a boundary condition
is expressed by the stationary parabolic response across the lane-junction according to Equation (35).
For an adhering axon, we choose
as the smallest distance reported for a cell membrane on a silicon chip [
50]. The half-width of the junction is
for
. With
, the sheet resistance is
. The profile
at an arbitrary time is plotted in
Figure 11 (bottom). The response along the lane is almost biphasic positive-negative with peak and trough separated by 120
µm in space and 140
µs in time. The amplitudes in the center of the junction are 1.3
µV and −1.9
µV.
The total voltage in the cable-substrate junction is obtained by superposing the voltage drop in the junction and the voltage in the surrounding electrolyte. Considering that the current of an adhering cable flows into one half-space, the amplitudes of the voltage in the electrolyte are 12 µV and −18 µV. The contribution of the electrolyte is by an order of magnitude larger than the contribution of the cable-substrate junction with 1.3 µV and −1.9 µV.
Transistor recording of axon in retina. The extracellular voltage of a propagating AP has been studied by transistor recording of axons in the rabbit retina [
51]. The retina was mounted with ganglion cells down on an array of EOMOS transistors. Maps of the transistor response
are displayed in
Figure 12. The velocity of the AP is 1.5
m/s. The transistor response is biphasic positive-negative with peak and trough separated by 120
µm in space and 80
µs in time. The amplitudes are about 500
µV and −675
µV.
The asymmetric biphasic waveform and the separation of peak and trough correspond to the computation shown in
Figure 11. The agreement indicates that the transistor response is due to a current-induced response with
. The velocity may be matched by choosing a larger radius or lower intracellular resistivity according to the relation
with
[
30]. Striking is the large amplitude of the recordings. Two issues must be considered: (i) for an axon near the substrate, even without close adhesion, the whole membrane current flows into a half-space with the mounted retina; and (ii) the retina has a resistivity of about 1250 Ω
cm [
52]. Thus, the result of the computation shown in
Figure 11 must be multiplied by a factor 25, yielding amplitudes of 150
µV and −240
µV. The remaining difference with respect to the experiment may due to a larger axon radius or higher membrane conductances.
Transistor recording with adhering axon in culture. The extracellular voltage of a propagating AP has been studied for adhering axons of cultured mammalian neurons using an array of metal electrodes on a CMOS chip (see
Appendix B) [
53]. The amplitudes of the recordings were in the range of 10
µV. Peak and trough of biphasic signals were separated by about 150
µs. These values are in agreement with the computation shown in
Figure 11.
The propagation of an AP along an adhering mammalian axon has also been probed with an array of nanowire transistors on a planar substrate [
54]. Almost monophasic positive transients of the nanowire current were reported corresponding to an extracellular voltage in the range of −85
mV [
55] or −6.6
mV, respectively [
56]. Even the lower estimate is four orders of magnitude larger than computed for the voltage drop in a cable-substrate junction and three orders of magnitude larger than computed for the voltage in the surrounding electrolyte. It is difficult to see how the parameters could be changed to avoid that discrepancy. Moreover, the waveform of the reported transients, which resemble the inverted intracellular voltage of a local AP, are difficult to reconcile with the waveform of the extracellular voltage of a propagating AP [
57].
9. Receptor-Cell-Transistor Sensor
Transistor recording with intact cells is an approach to develop biosensors with cellular receptors that bind agonists or pharmaceuticals. The receptor may activate an ion conductance in the plasma membrane, which gives rise to a current-induced response of the transistor. Considering that ligand-activated conductances remain open for a long time due to slow inactivation, a change of the extracellular ion concentrations is expected. An evaluation in terms of a current-induced response requires a gate oxide with a weak ion sensitivity. First, the concept is considered with the one-compartment model for constant ion concentrations (quasi-stationary response, small voltages), and then the prototype with a serotonin receptor is described.
Receptor-activated conductance. We assume that the receptor-activated conductance with a high reversal voltage
dominates the membrane. According to the two-domain model of Equation (29), the current balances in the junction and in the cell are expressed by Equation (36), with a receptor-activated conductance
in the adherent membrane and an average conductance
.
When the conductance is activated, the intracellular voltage relaxes rapidly towards the reversal voltage according to Equation (36b). According to Equation (36a), there is no stationary junction voltage, only a short transient, which is difficult to detect.
Supplementary K conductance. A stationary driving voltage
is established if a voltage-activated potassium conductance
is added with
. The junction voltage results from a superposition of the ion currents of both conductances according to Equation (37a). According to Equation (37b). The intracellular voltage approaches a resting voltage
.
If the potassium current is eliminated, we obtain Equation (38), where the junction voltage is determined by the driving voltage
and the inhomogeneity of the receptor-activated conductance.
On this basis, the transistor sensorics is described as follows: a cell stays at its low natural resting voltage with a deactivated receptor conductance and a deactivated K conductance. The agonist is added, and current flows into the cell through the receptor-activated conductance. The intracellular voltage rises such that the K conductance is activated. K current flows out of the cell, and a resting voltage is established. The superposition of the two ion currents through the adherent membrane yields a junction voltage for a depletion and a junction voltage for an enhancement . Finally, the receptor-activated conductance is inactivated, the intracellular voltage drops, the K conductance is deactivated, and the natural resting voltage is reestablished.
Inhomogeneity induced by diffusion. It is not necessary that the inhomogeneity of a receptor-activated conductance is structurally implemented. When an agonist is added to the bath, it first reaches the free membrane of an adherent cell and, with a delay, the adherent membrane due to diffusion of the agonist along the cell-substrate junction. The first phase with an activated free membrane is equivalent to . The second phase with an activated adherent membrane and an inactivated conductance in the free membrane is equivalent to . Prerequisite of such an effect is a matched dynamics of diffusion and inactivation.
The delayed activation of the adherent membrane by diffusion has been demonstrated for the serotonin receptor 5-HT3A in a HEK293 cell under voltage-clamp [
58]. A transistor with a weak ion sensitivity was used as described above. Serotonin (100
µM) was added with a double-barreled pipette. The apparent change of the gate voltage
is plotted in
Figure 13 versus the total membrane current
. Three phases can be distinguished: (i) the membrane inward current increases for 20
ms without response of the transistor due to an activation of the free membrane; (ii) the inward current decreases and the apparent change of the gate voltage increases for 70
ms due to inactivation of the free membrane and activation of the adherent membrane; and (iii) both the membrane current and the transistor signal disappear due to complete inactivation.
Transistor sensorics with serotonin receptor. A complete biosensor with the 5-HT3A receptor has been implemented with a supplementary Kv1.3 channel in HEK293 cells on a transistor with a low ion sensitivity of the gate [
35]. A result is shown in
Figure 14. When serotonin is applied, a positive transistor signal of about 600
µV appears. Within about 100
ms, it changes to a negative signal of about −450
µV, which relaxes slowly. The whole signal is abolished if the specific antagonist tropisetron is added.
The biphasic waveform of the transistor record is expected for a current-induced response with . The two phases correspond to a switched inhomogeneity of the membrane from depletion to enhancement of the receptor-activated conductance in the junction, as caused by diffusion and inactivation. The asymmetry of the positive and negative wave and the small amplitudes are due to the strong overlap of the two phases as the adherent membrane becomes activated at a time when the free membrane is only partially inactivated.