SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Review

Microscopic modeling of charge transport in sensing proteins

Lino Reggiani12*, Jean-Francois Millithaler1 and Cecilia Pennetta2

Author Affiliations

1 Dipartimento di Ingegneria dell’Innovazione and CNISM, Università del Salento, Via Arnesano, Lecce, 73100, Italy

2 Dipartimento di Matematica e Fisica Ennio De Giorgi and CNISM, Università del Salento, Via Arnesano, Lecce, 73100, Italy

For all author emails, please log on.

Nanoscale Research Letters 2012, 7:340  doi:10.1186/1556-276X-7-340

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/340


Received:6 January 2012
Accepted:22 June 2012
Published:22 June 2012

© 2012 Reggiani et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Sensing proteins (receptors) are nanostructures that exhibit very complex behaviors (ions pumping, conformational change, reaction catalysis, etc). They are constituted by a specific sequence of amino acids within a codified spatial organization. The functioning of these macromolecules is intrinsically connected with their spatial structure, which modifications are normally associated with their biological function. With the advance of nanotechnology, the investigation of the electrical properties of receptors has emerged as a demanding issue. Beside the fundamental interest, the possibility to exploit the electrical properties for the development of bioelectronic devices of new generations has attracted major interest. From the experimental side, we investigate three complementary kinds of measurements: (1) current-voltage (I-V) measurements in nanometric layers sandwiched between macroscopic contacts, (2) I-V measurements within an AFM environment in nanometric monolayers deposited on a conducting substrate, and (3) electrochemical impedance spectroscopy measurements on appropriate monolayers of self-assembled samples. From the theoretical side, a microscopic interpretation of these experiments is still a challenging issue. This paper reviews recent theoretical results carried out within the European project, Bioelectronic Olfactory Neuron Device, which provides a first quantitative interpretation of charge transport experiments exploiting static and dynamic electrical properties of several receptors. To this purpose, we have developed an impedance network protein analogue (INPA) which considers the interaction between neighboring amino acids within a given radius as responsible of charge transfer throughout the protein. The conformational change, due to the sensing action produced by the capture of the ligand (photon, odour), induces a modification of the spatial structure and, thus, of the electrical properties of the receptor. By a scaling procedure, the electrical change of the receptor when passing from the native to the active state is used to interpret the macroscopic measurement obtained within different methods. The developed INPA model is found to be very promising for a better understanding of the role of receptor topology in the mechanism responsible of charge transfer. Present results point favorably to the development of a new generation of nano-biosensors within the lab-on-chip strategy.

Keywords:
Impedance network; Bacteriorhodopsin; Olfactory receptors; Conformational change; Nano-biosensors

Review

Receptors are proteins of relevant interest because of the fundamental role they play in living environments at a cellular level [1,2]. At present, the best known (and studied) are bacteriorhodopsin (bR), a light sensitive protein whose action is related to a proton pump, and some of the G protein-coupled receptors, a large class of proteins that are sensitive to the light or more generally to the capture of single or a few specific molecules (ligands). The activation of a receptor is driven by the capture of an external ligand [3] and starts with a conformational change followed by a biological chain of events finally detected by the brain. The intriguing question is whether this biological chain of detection could be replaced by an electrical chain monitored by the change of an electrical property of the single protein induced by the change of conformation. To answer this question, the investigation of the electrical properties of a protein is the essential step.

To this purpose, we mention recent experiments on the current-voltage (I-V) characterization of bR monolayers sandwiched between metallic contacts [4,5]. Measurements were carried out in dark and when the sample was irradiated with a green light to which bR is sensitive. At increasing voltages, the current was found to exhibit a superlinear behavior which is common to both the illumination conditions. Furthermore, for each voltage value, the current is found to be significantly enhanced up to about a factor of three by the presence of light. These results, besides confirming the very low electrical conductivity of bR, suggest an underlying tunneling mechanism of charge transport [6,7]. On the other hand (most interesting for technological applications), following the I-V response, it is possible to monitor the protein activation, i.e., the conformational change. By implementing an atomic force microscope (AFM) technique, I-V measurements on bR were extended to higher applied voltages (near to about 10 V) where a cross over from direct to injection (or Fowler Nordheim [6]) tunneling regime was evident [8,9]. On the same subject, recent electrical experiments have been performed on the rat olfactory receptor (OR) OR I7 and on the human OR 17-40 properly deposited on a functionalized gold substrate with the technique of molecular self-assembly. These measurements [10,11] showed the possibility to monitor the sensing action of the protein and, in turn, their conformational change by means of the modification of the electrochemical impedance spectrum in the presence of a controlled flux of specific odorants (heptanal, octanal). From a theoretical side, a microscopic interpretation of all these experiments is still a challenging issue.

The aim of the present work is to report on a microscopic model, recently developed by the authors [3,12-16], which is able to interpret the electrical properties of a single protein and their modifications due to the sensing action in terms of the change of the protein tertiary structure. By implementing an impedance network protein analogue (INPA), the theoretical model is henceforth called INPA model. To this scope, we further included in the model a sequential tunneling mechanism [14-16] able to account for the superlinear I-V characteristics experimentally evidenced in bR and similar proteins at increasing applied voltages. Then, the model is validated on available experiments and its predictability discussed in the perspective of further experimental confirmations.

Model and materials

The INPA model is briefly summarized below. Each protein is associated with a topological network (graph) consisting of nodes and links. Once the tertiary structure of the protein is given by the protein database (PDB) [17] or similarly, the graph reproduces its main features. To this purpose, the nodes are in correspondence with the amino acids, taken as single interacting centers, and their position coincides with the related Cα atom. Each couple of nodes is connected with a link when their distance is less than an assigned interacting radius, RC. In this way, the protein network analogue becomes a set of intercrossing spheres of radius RC. In principle, all the values of RC from zero to tens of nanometers (the protein size) are possible. Actually, this value depends on the kind of interaction under consideration [12,18-20]. In the present case, as interaction, we consider the transfer of charge. Accordingly, RC is taken in the range 5 to 15 Å, as will be detailed in the following. Finally, the network gives a skeleton of the protein topology in a fixed conformation.

Figure 1 shows how the full atomic representation of a protein transforms into the sketch of a possible graph. The graph turns into an impedance network when an elemental impedance is associated with each link. In the present case, following the experimental outcomes, elementary impedance consisting of a resistance R parallelly connected with a capacitance C as shown in Figure 2. The elemental impedance between the i-th and j-th nodes is taken as [21]:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/340/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/340/mathml/M1">View MathML</a>

(1)

thumbnailFigure 1. Schematic of the INPA model. The left-hand side shows the all-atom picture of a protein. The right-hand side is the partial sketch of the equivalent network obtained for a given value of the interaction radius, RC (see text).

thumbnailFigure 2. Schematic of the elemental impedance. The left-hand side is the selection of the interacting amino acids (nodes) within a distance smaller than the interaction radius, RC. The right-hand side is the schematic view of the elementary RC impedance replacing one link between the α-carbon atoms of the i-th and j-th amino acids when their distance li,j is smaller than RC.

where <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/340/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/340/mathml/M2">View MathML</a>, is the cross-sectional area between the spheres centered on the ij nodes, respectively; li,j is the distance between these centers; ρis the resistivity taken to be the same for every amino acid with the indicative value, ρ = 1010Ωm (here, the microscopic mechanism responsible of charge transfer is not specified); <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/340/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/340/mathml/M3">View MathML</a> is the imaginary unit; ε0 is the vacuum permittivity; and ωis the circular frequency of the applied harmonic voltage. The relative dielectric constant pertaining to the couple of ij amino acids, εi,j, is expressed in terms of the intrinsic polarizability of the ij amino acids [21].

By positioning the input and output electrical contacts on the first and last node, respectively for a given applied bias, the network is solved within a linear Kirchhoff scheme, and its global impedance spectrum, Z(ω), is calculated in the standard frequency range 0.1 to 105 Hz. By construction, besides the small-signal dynamic response, this network produces a parameter dependent static I-V characteristic determined as

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/340/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/340/mathml/M4">View MathML</a>

(2)

However, since the purpose is to monitor the impedance variation of the receptor due to the conformational change, the relative change of impedance, rather than its absolute value, is taken as the relevant output of the modeling. The spatial distribution of the obtained network is found to be, in general, highly irregular. Accordingly, it describes the protein topology only when the tree structure is dominant, i.e., when the value of RC is sufficiently small so that the links between nearest neighbors predominate. We remark that the tree structure might not suffice to reproduce the global functioning of the protein. Thus, it is necessary to find an optimal value of RC which can reproduce both the structure and the function of the protein [12,13,20]. From one hand, the RC value should be greater than the mean distance between contiguous Cα atoms, say 3.8 Å; otherwise, the network appears to be not connected, and the value of the single protein impedance tends to be infinite. From another hand, for very large values of RC, say 80 Å, each node will link all the others and the network is no longer descriptive of the protein topology. In other words, for values of RC above about 80 Å the value of the impedance of the single protein no longer depends on RC. Accordingly, the best agreement between structure and function was obtained with values of RC in the range 5 to 70 Å [12,13]. To determine nonlinear I-V characteristics, the model was implemented to include a sequential tunneling mechanism of charge transfer by means of a stochastic selection of low resistive links ruled by the voltage-dependent probability [14,15]:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/340/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/340/mathml/M5">View MathML</a>

(3)

where Vi,j is the potential drop between the i-th and j-th node; m is an electron effective mass (taken as the bare one); e is the electron charge; is the reduced Planck constant; and Φ is the barrier height. The model has been further implemented to include also injection tunneling [16].

Comparison between theory and experiments

As anticipated above, the INPA is able to interpret the behavior of the I-V characteristics carried out in bR with an electrode-bilayer-electrode structure and when the protein is illuminated or less by green light [4]. To this purpose, the model is applied to the PDB entries 2NTU (for the native state) and 2NTW (for the activated state of bR) [17]. The comparison between theory and experiments proceeds as follows. The theoretical model describes the I-V modifications of a single protein when it undergoes a conformational change. On the other hand, the experimental results are carried out on a macroscopic sample that contains a macroscopic number of proteins. To compare the results, we normalize the current value obtained for the single protein in its native state at 1 V to that of the experiment in dark. For the case of bR, the normalization corresponds to multiply the current of the single protein by a factor 108 to 109. The same normalization factor is used for the single protein in its activated state, which is then compared with the experimental value under green light.

Figure 3 reports the simulated results calculated for both the native and activated states in comparison with experiments. The trend evidenced by experiments is reproduced when considering an interacting radius of 6 Å and a single energy barrier between neighboring amino acids of 59 meV [14,15]. The agreement with experiments is improved when using a Gaussian distribution of barrier heights with average value of 69 meV and a standard deviation of 44 meV [14,15].

thumbnailFigure 3. Current-voltage characteristics on a large-area monolayer of bacteriorhodopsin. Measured I-V characteristics of the native and activated state of bR are compared with theoretical calculations carried out using the 2NTU and 2NTW structures for the native and activated states, respectively. (a) - Full squares (circles) refer to calculations of the native (activated) state) carried out with Φ = 53 meV the same for each barrier; continuous curves are guides to the eyes, dashed curves refer to experiments carried out on a monolayer of bR as schematically shown in the inset of the figure [4]. (b) - The same as Figure 3 with calculations carried out with an average Φ = 69 meV Gaussian distributed with a variance σ = 44 meV. All the calculated currents are normalized to the experimental value of current in dark at 1 V [14].

Figure 4 compares the numerical and experimental data for the measurements carried out on bR within an AFM environment. The position of the contacts is then changed in order to reproduce the experiments obtained at different indentation of the AFM tip. At increasing depths of the extended contact, the net effect leads to a reduction in the number of amino acids involved in the electrical transport. As a consequence, we expect higher currents and a shift to lower potential values for the crossover between direct and injection tunneling regimes. Numerical calculations reported in Figure 4 confirm these expectations [16]. In particular, the flat contact that simulates the tip indentation was placed at the depths of 0.50 nm and 1.0 nm, respectively, from the top of the protein. For these depths, at increasing tip indentations, the experimental data exhibit a quantitative agreement with the modeling in the region of high voltages where the injection tunneling prevails. Even if the experimental indentation of the tip is greater for about a factor of three with respect to that of the simulations, we consider the agreement between theory and experiments to be satisfactory in view of the simplifications needed to convert the single protein calculation into the macroscopic value measured [16]. However, at low voltages, the results of the simulations underestimate for up to an order of magnitude the values of experiments. The disagreement at low voltages is overcome by assuming the existence of a leakage contribution, probably associated with the complexity of the contact regions constituted by trimers and lipids, with respect to the used model [22]. Accordingly, the single resistance associated with each link is replaced by the parallel of two resistances, one pertaining to the protein and the other to a more realistic modeling of the contact regions. Interestingly enough, the microscopic interpretation is obtained with a length-independent electron effective mass, contrary to the case of [8,9] where, to fit experiments, the carrier effective mass should increase over one order of magnitude at increasing the tip indentation. Both these features are consequences of the sequential tunneling mechanism [14-16] that replaces the single tunneling mechanisms of the metal-insulator-metal model previously used [8,9].

thumbnailFigure 4. Current voltage characteristics within an AFM technique on a monolayer of bacteriorhodopsin. Measured I-V characteristics of the native state of bR are compared with theoretical calculations carried out using the 2NTU structure at the reported different indentations of the AFM tip. The AFM technique is sketched in the inset. The thick continuous lines show experimental data [8] and the symbols numerical simulations obtained with extended contacts including a leakage current [16].

The impedance spectrum of a single olfactory receptor is investigated over a wide range of frequencies, and the results are reported by the associated Nyquist plot. This plot draws the negative imaginary part versus the real part of the global impedance within a given frequency range (typically from 1 mHz to 100 kHz as in experiments [23]). Figure 5 reports the Nyquist plots of human OR 17-40 with the impedance normalized to the static value of the native state. Symbols pertain to experiments, with crosses referring to the absence of a specific odorant, full (empty) squares pertain to the presence of the specific odorant helional (heptanal) at the concentrations reported in the figure. Curves pertain to theoretical results where the single protein is taken to be representative of the entire sample, and with continuous (dashed) lines referring to the native (active) state [3]. The agreement between theory and experiments is found to be satisfactory from a qualitative point of view and, apart for a underestimation of the maximum experimental value, acceptable from a quantitative point of view. We remark that the near-ideal semicircle shape of the experimental Nyquist plot is well reproduced, thus, confirming that the network impedance model behaves closely to a single RC circuit as expected by the presence of a rather uniform distribution of time constants associated with the different values of the resistance and capacitance of the links [13]. The scarcity of experimental data and the lack of a well-certified knowledge for the tertiary structures of the proteins under investigation lead us to consider these results as a first but significant step towards a microscopic modeling of the electrical properties of this olfactory receptor.

thumbnailFigure 5. EIS investigation on human OR 17-40. Nyquist plot for the human OR 17-40 in the absence and in the presence of specific ligands (heptanal and helional molecules). The impedances are normalized to the static value of the native state Zreal(0). In the upper figure, symbols report experiments carried out with a standard three electrode electrolytic cell sketched on the right-hand side part of the figure (with WE, working electrode; CE, counter electrode; RE, reference electrode). Here, crosses refer to the absence of an odorant and empty (full) squares to the presence of an heptanal (helional) odorants at concentration of 10−10M. In the lower part of the figure, lines report theoretical results with the continuous curve referring to the native state configuration with RC = 70 Å and the dashed line referring to the active state configuration with RC = 46 Å [3].

Conclusions

The microscopic modeling of the electrical properties of a sensing protein is briefly reviewed with the objective of providing a unifying theoretical interpretation of available experiments. In this framework, protein tertiary structure plays a primary role. Indeed, by determining the network structure, its conformational change induces the network change and, consequently, the change of the electrical properties of the protein. The possibility to monitor the conformational change by means of electrical measurements is identified as the key point for a better knowledge of charge transport in biological materials and for exploiting new frontiers of application of sensing proteins. The qualitative and quantitative agreements between the numerical results and experiments pose the INPA, implemented to account for a sequential tunneling mechanism of charge transfer, as a physical plausible model to investigate the electrical properties in other proteins pertaining to or less to the receptor family.

Competing interests

The authors declare that they have no competing interests.

Authors’contributions

LR participated in the design of the study and drafted the manuscript. JFM participated in the design of the study and performed the statistical analysis. CP conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Authors’ information

LR is full professor in Physics of Matter at the Engineering Faculty of Salento University. JM received a PhD in micro-electronics from Montpellier University and has a post doc position at the Dipartimento di Ingegneria dell’ Innovazione of the Salento University. CP is an associate professor in Physics of Matter at the Science Faculty of Salento University.

Acknowledgements

The collaboration of Dr. Eleonora Alfinito is deeply acknowledged. This research is carried out within the Bioelectronic Olfactory Neuron Device (BOND) Project sponsored by the Commission of the European Community (CEC) within the 7th Program (grant agreement: 228685-2).

References

  1. American Institute of Physics: Olfaction and electronic nose. In Proceeding 13 International Symposium: April 15-17 2009; Brescia, Italy. Edited by Pardo M, Sberveglieri G. AIP Press; 2009:115-118. OpenURL

  2. Bioelectronic Olfactory Neuron Device (BOND): Collaborative project FP7-NMP-2008-SMALL-2, GA number 228685. [ http://bondproject.org/ webcite]

  3. Alfinito E, Millithaler JF, Reggiani L, Zine N, Jaffrezic-Renault N: Topological and electrical properties of 17-40 human olfactory receptor for the realization of a nanobiosensor.

    RCS Adv 2011, 1:123-127. OpenURL

  4. Jin YD, Friedman N, Sheves M, He T, Cahen D: Bacteriorhodopsin (bR) as an electronic conduction medium: current transport through bR-containing monolayers.

    PNAS 2006, 103:8601-8606. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Jin YD, Friedman N, Sheves M, Cahen D: Bacteriorhodopsin-mnolayer-based planar metal-insulator-metal junctions via mimetic vesicle fusion: preparation, characterization and bio-optoelectronic characteristics.

    Adv Funct Mater 2007, 17:1417-1428. Publisher Full Text OpenURL

  6. Simmons JG: Generalized formula for the elctron tunnel effect between similar electrodes separated by a thin insulating film.

    J Appl Phys 1963, 34:1793-1803. Publisher Full Text OpenURL

  7. Wang W, Lee T, Reed WA: Electron tunnelling in self-assembled monolayers.

    Rep Prog Phys 2005, 68:523-544. Publisher Full Text OpenURL

  8. Casuso I, Fumagalli L, Samitier J, Padrós E, Reggiani L, Akimov V, Gomila G: Electron transport through supported biomembranes at the nanoscale by conductive atomic force microscopy.

    Nanotechnology 2007, 18:465503-465511. PubMed Abstract | Publisher Full Text OpenURL

  9. Casuso I, Fumagalli L, Samitier J, Padrós E, Reggiani L, Akimov V, Gomila G: Nanoscale electrical conductivity of the purple membrane monolayer.

    Phys Rev E 2007, 76:041919-041924. OpenURL

  10. Hou Y, Helali S, Zhang A, Jaffrezic-Renault N, Martelet C, Minic J, Gorojankina T, Persuy M-A, Pajot-Augy E, Salesse R, Bessueille F, Samitier J, Errachid A, Akimov V, Reggiani L, Pennetta C, Alfinito E: Immobilization of rhodopsin on a self-assembled multilayer and its specific detection by electrochemical impedance spectroscopy.

    Biosens Bioelectron 2006, 21:1393-1402. PubMed Abstract | Publisher Full Text OpenURL

  11. Hou Y, Jaffrezic-Renault N, Martelet C, Zhang A, Minic J, Gorojankina T, Persuy M-A, Pajot-Augy E, Salesse R, Akimov V, Reggiani L, Pennetta C, Alfinito E, Ruiz O, Gomila G, Samitier J, Errachid A: A novel detection strategy for odorant molecules based on controlled bioengineering of rat olfactory receptor I7.

    Biosens Bioelectron 2007, 22:1550-1555. PubMed Abstract | Publisher Full Text OpenURL

  12. Alfinito E, Pennetta C, Reggiani L: A network model to correlate conformational change and the impedance spectrum of single proteins.

    Nanotechnology 2008, 19:065202-065214. PubMed Abstract | Publisher Full Text OpenURL

  13. Alfinito E, Pennetta C, Reggiani L: Topological change and impedance spectrum of rat olfactory I7: a comparative analysis with bovine rhodopsin and bacteriorhodopsin.

    J Appl Phys 2009, 105:084703-084708. Publisher Full Text OpenURL

  14. Alfinito E, Reggiani L: Detecting conformational change by current transport in bacteriorhodopsin.

    Europhys Lett 2009, 85:68002-68008. Publisher Full Text OpenURL

  15. Alfinito E, Reggiani L: The role of topology in electrical properties of bR and rat-olfactory receptor.

    Phys Rev E 2010, 81:032902-032905. OpenURL

  16. Alfinito E, Millithaler J-F, Reggiani L: Charge transport in purple membrane monolayers: a sequential tunneling approach.

    Phys Rev E 2011, 83:042902-042905. OpenURL

  17. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne1 PE: The protein deata bank.

    Nucleic Acids Res 2000, 28:235-242. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Tirion MM: Large amplitude elastic motions in proteins from a single-parameter, atomic analysis.

    Phys Rev Lett 1996, 77:1905-1908. PubMed Abstract | Publisher Full Text OpenURL

  19. Micheletti C, Carloni P, Maritan A: Accurate and efficient description of protein vibrational dynamics: comparing molecular dynamics and Gaussian models.

    Proteins 2004, 55:635-647. PubMed Abstract | Publisher Full Text OpenURL

  20. Juanico B, Sanejouand YH, Piazza F, De Los Rios P: Discrete breathers in nonlinear network models of proteins.

    Phys Rev Lett 2007, 99:238104-238107. PubMed Abstract | Publisher Full Text OpenURL

  21. Pennetta C, Akimov V, Alfinito E, Reggiani L, Gorojankina T, Minic J, Pajot E, Persuy MA, Salesse R, Casuso I, Errachid A, Gomila G, Ruiz O, Samitier J, Hou Y, Jaffrezic N, Ferrari G, Fumagalli L, Sampietro M: Towards the realization of nanobiosensors based on G-protein-coupled receptors. In Nanodevices for the Life Sciences (Nanotechnology of the Life Science). Volume 4. Edited by Kumar CSSR. Weinheim: Wiley-VCH; 2006:217-240. OpenURL

  22. Hianik T: Biological membranes and membrane mimics. In Bioelectrochemistry: Fundamentals, Experimental Techniques and Applications. Edited by Bartlett PN. England: Wiley J & Sons; 2008:87-148. OpenURL

  23. Benilova IV, Vidic M, Pajot-Augy E, Soldatkin AP, Martelet C, Jaffrezic-Renault N: Electrochemical study of human olfactory receptor OR 17-40 stimulation by odorants in solution.

    Mater Sci Eng C 2008, 28:633-639. Publisher Full Text OpenURL