Computational Design of a Carbon Nanotube Fluorofullerene Biosensor
Abstract
: Carbon nanotubes offer exciting opportunities for devising highly-sensitive detectors of specific molecules in biology and the environment. Detection limits as low as 10−11 M have already been achieved using nanotube-based sensors. We propose the design of a biosensor comprised of functionalized carbon nanotube pores embedded in a silicon-nitride or other membrane, fluorofullerene-Fragment antigen-binding (Fab fragment) conjugates, and polymer beads with complementary Fab fragments. We show by using molecular and stochastic dynamics that conduction through the (9, 9) exohydrogenated carbon nanotubes is 20 times larger than through the Ion Channel Switch ICS™ biosensor, and fluorofullerenes block the nanotube entrance with a dissociation constant as low as 37 pM. Under normal operating conditions and in the absence of analyte, fluorofullerenes block the nanotube pores and the polymer beads float around in the reservoir. When analyte is injected into the reservoir the Fab fragments attached to the fluorofullerene and polymer bead crosslink to the analyte. The drag of the much larger polymer bead then acts to pull the fluorofullerene from the nanotube entrance, thereby allowing the flow of monovalent cations across the membrane. Assuming a tight seal is formed between the two reservoirs, such a biosensor would be able to detect one channel opening and thus one molecule of analyte making it a highly sensitive detection design.1. Introduction
Nanotubes have already attracted attention for use as highly sensitive biosensors due to their unique electrical and mechanical properties. Functionalized carbon nanotubes that display distinct characteristics triggered by changes in their surrounding environment offer exciting opportunities for the detection of specific molecules, particularly in biological and environmental samples [1]. Already, nanotube-based biosensors have shown detection limits as low as 10−11 M (for reviews see references [2] and [3]). The majority of nanotube-based sensors exploit their electrical properties, but attention has been drawn to sensors which make use of the nanotube pore for ionic conduction [4–9]. Most reported carbon nanotube-based biosensors are enzyme-modified electrodes for detecting different biomolecules (such as glucose) and environmental pollutants (such as organo-phosphate pesticides) [3]. Owing to their sensitivity, these biosensors may lead to numerous technological advances in the medical industry.
The excellent electrical property of carbon nanotubes has led to their incorporation into numerous biosensor designs. Moreover, since the current flows along the nanotube surface, the electrical conductance of the nanotube is highly sensitive to chemical changes on its sidewalls [10]. In these devices the binding of a target molecule to the carbon nanotube surface is detected by a change in electrical signal, such as through a change in conductance [11]. These types of biosensors have been used to detect E. coli [11], or glucose molecules [12,13]. Moreover, they can also be used to study dynamic processes, such as the release of ATP from living cells [10]. Sotiropoulou and Chaniotakis [12] used an array of multi-walled carbon nanotubes grown on a platinum substrate. When glucose oxidase enzymes are immobilized to the nanotube ends a change in electrical current occurs. By measuring currents that flow across the nanowire, they demonstrate the detection of glucose to 0.19 mM, which is among the best reported values for glucose biosensors. Carbon nanotube forests have also been developed to detect various proteins, such as human serum albumin [14], prostate specific antigen [15] and horseradish peroxidase [16]. Antibodies are covalently attached to the carboxylated ends of the carbon nanotubes and the specific protein is detected by a change in electrical current. For example, these forests could detect 1 nM of human serum albumin [14]. Liu et al. [17] demonstrated the potential of single-wall carbon nanotubes to detect volatile organic compounds in human breath to diagnose lung cancer. Their results show that the resistance of the nanotubes increases when exposed to both polar and non-polar molecules. Similarly, resistance changes of carbon nanotubes have also been used to detect viruses [18] and the toxic gases NO2 and NH3 [19]. Recently, Park et al. [20] developed a biosensor which mimics the canine nose comprised of olfactory nanovesicle fused carbon nanotube transistors which can detect hexanal, an indicator of oxidation of food, at 1 fM concentration by monitoring conductance change of the nanotube. As well as their unique electrical properties, nanotubes are also widely investigated due to their remarkable mechanical properties. In such sensors, the nanotube acts as a cantilever which is sensitive to the molecule of interest [21]. Similar to the change in electrical signal generated by the immobilization of specific proteins onto the nanotube, mechanical deformations of the nanotube cantilever can also induce changes to electrical conductance [21].
One of the simplest uses of nanotubes for sensing applications is the detection of analytes due to a change in transmembrane current, such as proposed in this paper. Carbon nanotube membranes with nanotube diameters of approximately 7.5 nm in diameter have been used to detect streptavidin [4] and desthiobiotin [5] when the ionic flux is reduced by a factor of 15 and 24, respectively. Nednoor et al. [6] also modulate ionic flux through their array of aligned carbon nanotubes 7 nm in diameter embedded in a polystyrene matrix by phosphorylation. Phosphorylation and dephosphorylation of synthetic peptides at the nanotube entrance modulate ionic flux as a result of antibody attachment by reversing the charges at the nanotube entrance. Gold nanotubes have also been used to detect specific proteins [7–9]. Siwy et al. [9] have embedded conically shaped gold nanotubes with diameters 5–9 nm in a polymeric membrane. Proteins are detected when they occlude the nanotube entrance, thus preventing the flow of ionic currents. Similarly, Kobayashi and Martin [7] fabricated membranes constructed from gold nanotubes with inside diameters ranging from 1 to 4 nm and achieved detection limits of 10−11 M by measuring the change in transmembrane current when an analyte effectively plugs the nanotube entrance. A disadvantage with these size-based selectivity designs is that the analyte species must have molecular dimensions comparable to the inside diameter of the nanotubes; thus, in practice, it would be difficult to distinguish between molecules of a similar size. Moreover, only proteins with comparable diameters to the nanotube mouth can be detected. In contrast, Steinle et al. [8] mimicked ligand-gated ion channels where their gold nanotube embedded membrane is switched from the “off” state to the “on” state in response to a chemical stimulus. The “off” state is obtained by making the membranes hydrophobic so that water and ions cannot enter the pore. Then, by adding sufficiently high concentrations of hydrophobic ionic species, such as a drug, the membrane switches to the “on” state since the pores fill with water and electrolyte.
The Ion Channel Switch ICS™ biosensor [22] makes use of ionic flow through the bacterial channel gramicidin-A. Gramicidin-A is a polypeptide consisting of 15 amino acid residues which forms a narrow pore of approximately 2 Å in radius. It was one of the first antibiotics used clinically as it forms a channel across bacteria cell membranes and selectively conducts monovalent cations [23,24]. This biosensor has been used for the detection of low molecular weight drugs, proteins and microorganisms [25]. It is constructed from a gold electrode to which is tethered a lipid membrane containing gramicidin ion channels linked to antibodies. The switch has a high gain, where a single channel conducts up to a million ions per second, and can detect analytes from sub-picomolar to micromolar concentrations in less than 10 minutes [22]. There are 108 gramicidin channels per cm2, so that each electrode area of 0.03 cm2 contains 3 × 106 channels with approximately half of these being dimerized [26] or open so that the measured current is of the order of microamperes. These currently used biosensors respond to concentrations as low as 10 fM [26]. Work is continuing to miniaturize these ICS™ biosensors, but electrode sizes of 1 μm radius only contain a few ion channels and therefore measured current is of the order of picoamperes and measurement noise becomes a significant issue [26].
Here we illustrate a design, based on a similar principle to the ICS™ biosensor [22], of an ultra-sensitive detector, comprised of an array of functionalized carbon nanotubes and fluorinated fullerenes (fluorofullerenes). Under normal conditions all carbon nanotubes are blocked by the fluorofullerenes. The fluorofullerenes are attached to Fragment antigen-binding (Fab fragments) and therefore not only block the nanotube pore but also act as recognition elements. Upon detection of some analyte the fluorofullerenes are stripped from the nanotube entrance, allowing monovalent cations to conduct at a rate of 82 pA per nanotube. We demonstrate, using molecular and stochastic dynamics, the permeation characteristics of monovalent cations through the unblocked carbon nanotubes, construct the free energy of binding of two fluorofullerenes to the nanotube pore and calculate their respective binding affinity. Finally, we outline the construction and function of our proof-of-concept biosensor design.
2. Theoretical Section
Using molecular dynamics (MD) and distributional molecular dynamics (DMD) simulations we examine a (9, 9) single-walled carbon nanotube (CNT) with an effective radius of 4.53 Å and length of 36 Å embedded in a lipid bilayer separating two reservoirs. The nanotube is terminated with hydrogen atoms with a partial charge of 0.115e, and −0.115e on carbon atoms directly bonded to these hydrogen atoms [27]. In addition, two regions on the nanotube outside surface are hydrated (exohydrogenated). These two regions are located from z = −11.8 to −6.9 Å and from z = 6.9 to 11.8 Å. Further details of the methodology can be found in our earlier papers [28–30].
2.1. Molecular Dynamics
MD simulations are used to calculate the profile of potential of mean force (PMF) of (i) a potassium or chloride ion traversing through the CNT and (ii) two fluorofullerenes (C60F36 and C60F60) binding to the (9, 9) CNT. All MD simulations are performed using NAMD 2.8 and visualized using VMD 1.9 [31,32]. All simulations use the CHARMM27 force field [33,34] and TIP3P water model, a constant pressure (1 atm) and temperature 310 (K). Atom coordinates and topology file for the C60 fullerene were taken from the NAMD mailing list [35]. Parameters for the fluorine atoms were obtained from CHARMM36 force field [33,34,36] and Dunlap et al. [37]. The carbon and fluorine atoms of the fluorofullerenes are given a partial charge of 0.071e and −0.071e, respectively [38]. The nanotube is embedded in a 3-palmitoyl-2-oleoyl-D-glycero-1-phosphatidylethanolamine (POPE) lipid bilayer separating two reservoirs. The PMF for the potassium and chloride ions is conducted using umbrella sampling at an ionic concentration of 0 mM to avoid double-counting of the ion-ion interactions in subsequent DMD simulations in which the ion-ion interactions are explicitly dealt with using macroscopic electrostatics [28]. The effect of this is discussed in our earlier paper [30]. Therefore, the simulation contained either one potassium or one chloride ion.
The PMF for the binding of the fluorofullerene to the carbon nanotube is obtained using umbrella sampling. The fluorofullerenes were initially placed at the nanotube entrance and equilibrated for 1 ns. Using the equilibrated nanotube-fullerene structure we generate sampling windows by performing steered molecular dynamics. A force constant of 10 kcal/mol/Å2 is applied to pull the fluorofullerene away from the nanotube entrance. The channel central axis is used as the reaction coordinate. The pulling generates a continuous number of configurations along the permeation pathway so that we can construct an umbrella sampling window every 0.5 Å. During umbrella sampling the centre of mass of the fluorofullerene is confined to be within a cylinder of 8 Å centred on the channel axis and beyond this a harmonic constraint of 20 kcal/mol/Å2 is applied. In addition, a force constant of 10 kcal/mol/Å2 is applied in the z direction to constrain the centre of mass to the sampling window. The centre of mass coordinates of the fluorofullerene are saved every 0.5 ps. The weighted histogram analysis method [39,40] and centre of mass coordinates are used to construct the PMF along the z direction. Each sampling window is run for 5 ns. The PMF is shown to converge as the depth changes by less than 0.5 kT when simulations are run for a further 1 ns. The dissociation constant (Kd) in the unit of molar is estimated to be [41,42]:
2.2. Distributional Molecular Dynamics
DMD simulations [28–30] are used to determine the current-voltage profile of potassium ions traversing through the carbon nanotube. DMD allows us to reproduce the distribution of ion trajectories from MD as closely as possible but enables simulations that are two to three orders of magnitude faster than classical MD simulations.
We first determine the free energy profile and frictional and random force (friction kernel) at discrete segments within the channel using MD. Then, we incorporate these results into stochastic dynamics simulations based on the non-linear generalized Langevin equation to determine the current-voltage profile, and ion binding sites within the channel. The simulation space is separated into two regions. In the first region, the pore (nanotube) is divided into thin slices and the forces acting on each ion and diffusion coefficient are obtained from the free energy and friction kernel generated using MD so that certain properties of the ion trajectory are reproduced in DMD. Macroscopic electrostatics is used for ion-ion interactions within the pore. We use a dielectric constant of 2 for the membrane region and 60 for water within the nanotube. These values have been shown to provide excellent agreement with experimental data for biological ion channels which have pores of a similar size [43]. In the second region, the reservoir, normal Brownian dynamics and macroscopic electrostatics is performed [44].
3. Results and Discussion
3.1. Ion Conduction
The (9, 9) exohydrogenated CNT 36 Å in length is shown in Figure 1(A). Figure 1(B) shows the axial free energy encountered by potassium (K+) and chloride (Cl−) ions as they traverse the CNT. The free energy profile is derived from the three-dimensional PMF and is therefore the energy each ion encounters as it moves along the axis. Potassium ions encounter two energy wells of approximately 3 kT at the each end of the nanotube, located at ±9 Å. These wells are separated by a shallow energy barrier of 2.2 kT. In contrast, chloride ions encounter an insurmountable energy barrier of approximately 20 kT. The location of the energy wells for potassium ions align with the location of the exohydrogenated regions of the nanotube shown in Figure 1(A). In the absence of an applied potential, potassium ions tend to dwell in the centre of the two exohydrogenated regions of the nanotube (at ±9 Å), as shown in Figure 1(C). On average, 1.2 potassium ions are present in the 36 Å length nanotube. The nanotube length examined in this paper is much shorter than the smallest nanotubes produced experimentally (200 Å) [45]. To maintain low barriers for potassium ions it will be necessary to exohydrogenate repeating sections of the entire length of the long nanotube to generate an undulating profile. In the exohydrogenated sections of the nanotube an ion experiences a negatively lined pore as a result of the negative charges on the carbon atoms (−0.115e).
The current-voltage profile for potassium ions for the (9, 9) exohydrogenated CNT, obtained with symmetrical solutions of 500 mM in both reservoirs, is shown in Figure 2(A). The current-voltage relationship is linear, with a current of 82 pA at 200 mV, equivalent to 5.1 × 108 ions per second. No chloride conductance was observed for all investigated voltages and concentrations. The magnitude of current across the nanotube plotted against the ionic concentration of potassium ions in the reservoirs follows the Michaelis-Menten form Figure 2(B). For example, at an applied potential of 200 mV, an ionic concentration of 200 mM reduces the current to 52 pA. The current increases rapidly with an increasing ionic concentration initially and then saturates with a further increase in concentration. The current saturates at approximately 127 pA, and the concentration for the half-maximum current occurs at 277 mM.
A single channel in the ICS™ biosensor [22] has an experimentally measured current of 4 pA at 200 mV [23,24]. The carbon nanotube presented here demonstrates a 20-fold increase in current compared to the ICS™ biosensor. For the nanotube to function as a biosensor it is necessary to block the nanotube pore, preferably when no analyte is present so that we observe conduction when analyte is present. In the next section we investigate two blocker molecules constructed from fluorinated C60 fullerenes.
3.2. Blockage by Fluorofullerenes
The carbon nanotubes examined in the previous section are permanently open in that monovalent cations move through the pore under the influence of an applied potential. For these nanotubes to form part of a biosensor we want the pore to be closed normally and to open when a specific molecule we wish to detect is present, so that the biosensor is “switched on” upon detection. It is therefore necessary to determine a suitable blocking molecule which can non-covalently bind to the nanotube entrance with sufficient energy that its probability of reopening is low. Moreover, the binding strength must not be too strong such that the blocking molecule cannot be stripped away when the molecule of interest is detected.
We examine two fluorofullerenes, namely C60F36 and C60F60, as potential blocking molecules. We assume T symmetry for the C60F36 fluorofullerene taken from Taylor [46] which is the most stable of the three isomers [47,48]. Fluorinated fullerenes were first synthesized in 1991 by the reaction of molecular fluorine with C60 [49]. Fluorofullerene molecules from C60F2 to C60F102 have been detected in the gas phase [49]. C60F18, C60F36 and C60F48 are readily produced [38,46] and experimentalists have shown that only species with an even number of fluorine atoms are observed [38]. It is difficult to add more than 48 fluorine atoms to the C60 cage, but fullerenes with higher fluorination have been observed in small concentrations such as the C60F60 fullerene [46,49]. We chose the C60F36 and C60F60 fluorofullerenes so as to obtain a range as to the fluorofullerene's potential binding strength to the carbon nanotube. Ab initio calculations have shown that tube-like isomers of C60F60 are more stable than the cage isomer [50]. In this paper we only study the cage isomer.
The centre of mass of the fluorofullerenes are placed 3 Å from the nanotube entrance (22 Å) and after 1 ns equilibrium simulations the fluorofullerenes remain at the nanotube entrance. To verify the predicted binding of these fluorofullerenes to the nanotube we determine the PMF for the unbinding from the nanotube. As shown in Figure 3, the PMF reaches a minimum at 22 Å and 23 Å for the C60F36 and C60F60 fullerenes, respectively, where the nanotube centre of mass is at 0 Å. For the C60F60 fullerene the PMF increases rapidly from −25.6 kT at 23 Å to −7.3 kT at 24.5 Å. The C60F36 fullerene exhibits a shallower well of −11.9 kT. It is assumed that the properties for the window at 35 Å are similar to bulk, and therefore the PMF is set to zero at this point. Unlike the C60F36 fullerene, the C60F60 fullerene also exhibits another shallower minimum energy position at approximately 25.5 Å. At this position the C60F60 fullerene binds to the edge of the nanotube partially occluding the nanotube entrance. The fluorofullerenes non-covalently bind to the nanotube entrance through favourable electrostatic and van der Waals interactions.
The bound configuration of C60F36 and C60F60 demonstrate complete blockage of the nanotube pore. This is illustrated in Figure 4 for the C60F60 fullerene, a similar blockage is observed for the C60F36 fullerene. Using Equation (1) we obtain a dissociation constant, Kd of 87 μM and 37 pM for the C60F36 and C60F60 fullerenes, respectively.
3.3. Biosensor Construction and Function
There are three main components of our biosensor design: the carbon nanotube pores embedded in a membrane, a fluorofullerene with a Fab fragment attached to its surface, and a soluble polymer bead with complementary Fab fragments. The function of our biosensor relies on the removal of the fluorofullerene from the nanotube entrance in response to analyte. However, to reduce noise of the system it is important that we choose a fluorofullerene which exhibits a large binding energy such that the probability of channels spontaneously opening and closing is as small as possible. In this section we will describe the construction of our three components and how they fit together to detect analyte.
It is possible to embed the carbon nanotubes into a highly dense array enabling the biosensor the greatest chance of detecting the analyte. For example, our exohydrogenated (9, 9) carbon nanotubes could be embedded in a silicon-nitride membrane similar to the work of Holt et al. [51]. They obtained a pore density of 2.5 × 1011 cm−2 in their carbon nanotube-silicon nitride membrane. However, their membrane thickness is 2–3 μm. To make our theoretical calculations tractable computationally we use a short tube. For the nanotubes we examine in this paper to span this thickness we must increase their length by approximately 1,000 times. But a longer nanotube will essentially produce the same results once we exohydrogenate sections along its entire length. For instance, the PMF for potassium ions would then be an undulating profile along the entire length with multiple energy wells and barriers and a K+ would only ever face a 2.2 kT barrier as in Figure 1(A).
From our MD simulations the C60F60 fullerene exhibits the largest binding energy, with a depth of 25.6 kT (Figure 3). The probability of the nanotube pore opening and closing using this fluorofullerene will be very small, approximately 10−8, so that noise in the system will be minimal and we can simply calculate our baseline current and the mean duration a pore stays open once it opens spontaneously. Next, we must attach Fab fragment to the fluorofullerene. Taylor et al. [52] demonstrated the possibility of including oxygen atoms on prefabricated fluorofullerenes, where epoxides on the nanotube surface are produced by nucleophilic substitution of F by OH. Moreover, fluorofullerenols, fluorofullerenes with hydroxides attached, are readily formed from C60F36 by reaction with moist solvent and can be obtained in significant yield [47]. Using this methodology we propose to attach 2–3 hydroxyl groups on the C60F60 fullerene, such that it will become C60F58(OH)2. These hydroxyl groups can then be used to functionally attach Fab fragments. First, C60-PEG conjugation will be performed on the fluorofullerenols as previously reported [53,54] with biotin attached to the PEG chain. This biotin-PEG chain would then be coupled to a biotinylated antibody fragment (Fab) using streptavidin intermediates, similar to that described by Cornell et al. [22]. These fluorofullerene conjugates would be placed on one side of the carbon nanotube membrane.
We propose to utilize polymer beads, with radii of approximately 15 nm [55], which have been widely investigated for use in medicine and biology through to opto-electronics [55]. The surface of the soluble polymer beads will be functionalized with antibodies or complementary Fab fragments to that on the fluorofullerene with synthesis according to the procedures described by Uchida et al. [56], van Erp et al. [57], Sarobe et al. [58], and Kawaguchi [55]. For example, Sarobe et al. [58] covalently coupled immunoglobulin G (IgG) to chloromethylstyrene beads to create stable sensitized latex beads. Uchida et al. [56] immobilized both intact IgG or Fab fragments of antibovine serum albumin antibody to polymer microspheres using an oxidized dextran as a binder. Moreover, DNA has been immobilized on polymer beads using biotin-avidin binding [55]. The beads will be placed on the same side of the membrane as the fluorofullerene conjugates.
When no analyte is present in the top reservoir, the fluorofullerene conjugates will block the entrances to the nanotube pores and the polymer beads will float around Figure 5(A). When analyte is injected into the top reservoir the Fab fragments present on the polymer bead and fullerene conjugate will crosslink to the analyte. The flow caused by the injection of analyte will act to pull the fluorofullerene from the nanotube entrance due to the drag forces on the polymer bead. This will in turn unblock the pore, as illustrated in Figure 5(B), and a potassium current will be detected and thus the analyte. Using the present design it would not be possible to detect multiple analytes with distinct binding affinities, such as a low or high affinity analyte. The current flowing across the open nanotubes may be monitored with a conventional current-to-voltage amplifier, such as the one used routinely for measuring ionic channel currents. If finer resolutions are required, we can make use of the signal processing techniques based on the hidden Markov model [59] or fluctuation analyses [60,61]. To open a pore, the drag force FD on these polymer beads must be sufficient to pull the fluorofullerenes from the nanotube entrance. In other words, the drag force must be greater than the fluorofullerene-nanotube binding energy. If we assume laminar flow past a sphere, Stokes Law states that FD = 6πμaU, where a is the radius of the sphere (15 nm), U is the flow speed and μ is the viscosity. If we assume an injection speed [62] of 1.2 mL/min and viscosity of water we obtain a drag force of approximately 1.2 × 10−13 N for a polymer bead 15 nm in radius. This is many orders of magnitude larger than the binding force of the fluorofullerene which for the C60F60 fullerene is approximately 3.5 × 10−20 N. The binding affinity between the Fab fragments and the analyte must also be equal to or greater than the fluorofullerene-nanotube binding affinity to ensure that the fluorofullerenes can be pulled from the nanotube mouth. We determine dissociation constants of 87 μM and 37 pM for the C60F36 and C60F60 fullerenes, respectively. Friguet et al. [63] determined the dissociation constant of antigen-antibody equilibria in solution in the nM range. Therefore, to ensure the fluorofullerene can be pulled from the nanotube entrance the fullerene must be functionalized with somewhere between 36 and 60 fluorine atoms. Fullerenes can be fabricated with between 2 and 60 fluorine atoms [47] but the C60F48 fullerene, which is formed readily [46], may provide the best blocking characteristics.
As described in Section 3.1, a single nanotube channel facilitates the flow of 5.1 × 108 ions per second. This magnitude is significantly larger than through the Ion Channel Switch ICS™ biosensor where a single channel facilitates 2.5 × 107 ions per second [22]. Moreover, assuming a similar pore density to that observed by Holt et al. [51], our membrane would exhibit a high channel density. For example, assuming all channels are initially blocked by fluorofullerene conjugates, there would be 2.5 × 1011 available recognition sites in a device 1 cm2. In work to miniaturize the ICS™ biosensor to 1 μm radius, noise becomes a significant issue since at this size there are only a few ion channels and measured current is of the order of picoamperes [26]. Using the biosensor design in this paper a 1 μm2 membrane would contain 2,500 nanotube channels and the measured current per nanotube is an order of magnitude larger. Therefore, although the background noise of our proposed biosensor is not known, it is likely that the sensitivity of our biosensor will be an order of magnitude smaller than the ICS™ biosensor. Furthermore, at voltages above 150–200 mV the ICS™ biosensor will not function since the membrane will break down. In contrast, voltages as high as 1 V could be applied to the fluorofullerene-nanotube biosensor further increasing the current flow through the tube. Some of the channels in the membrane of the ICS™ biosensor can diffuse laterally within the plane of the membrane [22]. The antibodies on the mobile channels scan an area of 1 μm2 in less than 5 min, and the speed and sensitivity of the biosensor response is proportional to the number of binding sites accessible to each mobile channel [22]. In other words, the sensitivity to analyte is enhanced since each gramicidin channel is able to cross link with any of many sites at the sensor surface to which the molecule to be detected may bind. In contrast, the fluorofullerene-carbon nanotube biosensor proposed here would be constructed of nanotubes fixed into their membrane, thereby removing the sensitivity gained from in-plane diffusion. However, as mentioned our membrane would exhibit a high density of channels which would improve its sensitivity due to the large number of recognition sites.
4. Conclusions
Using both molecular and stochastic dynamics, we show that an exohydrogenated (9, 9) carbon nanotube with hydrogen terminated ends exclusively conducts monovalent cations with a current of 82 pA at 200 mV. In comparison gramicidin-A conducts potassium ions with an experimentally measured current of 4 pA at 200mV [23,24]. Furthermore, we demonstrate that it is possible to completely block the nanotube pore using fluorofullerenes, with dissociation constants of 87 μM and 37 pM for the C60F36 and C60F60 fullerenes, respectively.
We propose at the proof-of-concept level the design of a biosensor comprised of the functionalized carbon nanotube pores embedded in a membrane, a fluorofullerene with a Fab fragment attached to its surface, and a soluble polymer bead with complementary Fab fragments. Under normal operating conditions, when no analyte is present, the nanotubes are blocked by a C60F60 fluorofullerene-Fab conjugate and the polymer beads float around in the reservoir. When analyte is injected into the top reservoir the Fab fragments present on the polymer bead and fullerene conjugate crosslink to the analyte. Then, the drag force on the much larger polymer bead will act to pull the fluorofullerene away from the nanotube and unblock the pore and thus allow the flow of monovalent cations across the membrane.
The advantages of this carbon nanotube fluorofullerene biosensor compared to the ICS™ biosensor are: (i) a 20-fold increase in potassium ion conductance, (ii) the potential to activate a large number of channels since the nanotubes can be densely packed into membrane, and (iii) the potential to develop a more sensitive current amplifier since current can be measured directly. Furthermore, assuming a tight seal is formed between the two reservoirs such a biosensor would be able to detect the flow of currents caused by the opening of only a few nanotubes and thus a few molecules of analyte making it a highly sensitive detection design.
Much work needs to be done before such a device is realized. For example, as a proof-of-concept we have examined only one nanotube diameter and terminal functionalization. Future work will be necessary to study the effect of nanotube diameter, and carbon nanotube tip functionalization on the fullerene binding energy. It is now possible to fabricate hydrogenated carbon nanotubes using an atomic hydrogen beam, and these functionalized nanotubes have been shown to be stable at room temperature [64]. Unfortunately, lack of control of functionalization location remains a challenge. However, Raghuveer et al. [65] demonstrated the ability to site-selectively functionalize the outside surface of multi-walled carbon nanotubes using focused-ion-beam irradiation and subsequent mild chemical treatments. Further work is needed in this area to facilitate the creation of devices such as our proposed biosensor.
Acknowledgments
We thank Silvie Ngo and Rhys Hawkins for expert graphics preparation. This work was supported by the Australian Research Council through a Discovery Early Career Researcher Award. The work was supported by the NCI National Facility at the Australian National University.
References
- Menard-Moyon, C.; Kostarelos, K.; Prato, M.; Blanco, A. Functionalized carbon nanotubes for probing and modulating molecular functions. Chem. Biol. 2010, 17, 107–115. [Google Scholar]
- Kohli, P.; Wirtz, M.; Martin, C.R. Nanotube membrane based biosensors. Electroanalysis 2004, 16, 9–18. [Google Scholar]
- Wang, J.; Lin, Y. Functionalized carbon nanotubes and nanofibers for biosensing applications. Trends Anal. Chem. 2008, 27, 619–626. [Google Scholar]
- Hinds, B.J.; Chopra, N.; Rantell, T.; Andrews, R.; Gavalas, V.; Bachas, L.G. Aligned multiwalled carbon nanotube membranes. Science 2004, 303, 62–65. [Google Scholar]
- Nednoor, P.; Chopra, N.; Gavalas, V.; Bachas, L.G.; Hinds, B.J. Reversible biochemical switching of ionic transport through aligned carbon nanotube membranes. Chem. Mater. 2005, 17, 3595–3599. [Google Scholar]
- Nednoor, P.; Gavalas, V.G.; Chopra, N.; Hinds, B.J.; Bachas, L.G. Carbon nanotube based biomimetic membranes: Mimicking protein channels regulated by phosphorylation. J. Mater. Chem. 2007, 17, 1755–1757. [Google Scholar]
- Kobayashi, Y.; Martin, C.R. Highly sensitive methods for electroanalytical chemistry based on nanotubule membranes. Anal. Chem. 1999, 71, 3665–3672. [Google Scholar]
- Steinle, E.D.; Mitchell, D.T.; Wirtz, M.; Lee, S.B.; Young, V.Y.; Martin, C.R. Ion channel mimetic micropore and nanotube membrane sensors. Anal. Chem. 2002, 74, 2416–2422. [Google Scholar]
- Siwy, Z.; Trofin, L.; Kohli, P.; Baker, L.A.; Trautmann, C.; Martin, C.R. Protein biosensors based on biofunctionalized conical gold nanotubes. J. Am. Chem. Soc. 2005, 127, 5000–5001. [Google Scholar]
- Huang, Y.; Sudibya, H.G.; Fu, D.; Xue, R.; Dong, X.; Li, L.J.; Chen, P. Label-free detection of ATP release from living astrocytes with high temporal resolution using carbon nanotube network. Biosens. Bioelectron. 2009, 24, 2716–2720. [Google Scholar]
- Lee, J.O.; So, H.M. Single-walled carbon nanotube biosensors. In Nanomaterials for the Life Sciences; Kumar, C.S.S.R., Ed.; Wiley-VCH Verlag: Weinheim, Germany, 2011; Volume 9, pp. 217–237. [Google Scholar]
- Sotiropoulou, S.; Chaniotakis, N.A. Carbon nanotube array-based biosensor. Anal. Bioanal. Chem. 2003, 375, 103–105. [Google Scholar]
- Lin, Y.; Lu, F.; Tu, Y.; Ren, Z. Glucose biosensors based on carbon nanotube nanoelectrode ensembles. Nano Lett. 2004, 4, 191–195. [Google Scholar]
- Yu, X.; Kim, S.N.; Papadimitrakopoulos, F.; Rusling, J.F. Protein immunosensor using single-wall carbon nanotube forests with electrochemical detection of enzyme labels. Mol. BioSyst. 2005, 1, 70–78. [Google Scholar]
- Yu, X.; Munge, B.; Patel, V.; Jensen, G.; Bhirde, A.; Gong, J.D.; Kim, S.N.; Gillespie, J.; Gutkind, J.S.; Papadimitrakopoulos, F.; Rusling, J.F. Carbon nanotube amplification strategies for highly sensitive immunodetection of cancer biomarkers. J. Am. Chem. Soc. 2006, 128, 11199–11205. [Google Scholar]
- O'Connor, M.; Kim, S.N.; Killard, A.J.; Forster, R.J.; Smyth, M.R.; Papadimitrakopoulos, F.; Rusling, J.F. Mediated amperometric immunosensing using single walled carbon nanotube forests. Analyst 2004, 129, 1176–1180. [Google Scholar]
- Liu, F.L.; Xiao, P.; Fang, H.L.; Dai, H.F.; Qiao, L.; Zhang, Y.H. Single-walled carbon nanotube-based biosensors for the detection of volatile organic compounds of lung cancer. Physica. E 2011, 44, 367–372. [Google Scholar]
- Bhattacharya, M.; Hong, S.; Lee, D.; Cui, T.; Goyal, S.M. Carbon nanotube based sensors for the detection of viruses. Sens. Actuators B Chem. 2011, 155, 67–74. [Google Scholar]
- Kong, J.; Franklin, N.R.; Zhou, C.; Chapline, M.G.; Peng, S.; Cho, K.; Dai, H. Nanotube molecular wires as chemical sensors. Science 2000, 287, 622–625. [Google Scholar]
- Park, J.; Lim, J.H.; Jin, H.J.; Namgung, S.; Lee, S.H.; Park, T.H.; Hong, S. A bioelectronic sensor based on canine olfactory nanovesicle-carbon nanotube hybrid structures for the fast assessment of food quality. Analyst 2012, 137, 3249–3254. [Google Scholar]
- Arlett, J.L.; Myers, E.B.; Roukes, M.L. Comparative advantages of mechanical biosensors. Nat. Nanotechnol. 2011, 6, 203–215. [Google Scholar]
- Cornell, B.A.; Braach-Maksvytis, V.L.B.; King, L.G.; Osman, P.D.J.; Raguse, B.; Wieczorek, L.; Pace, R.J. A biosensor that uses ion-channel switches. Nature 1997, 387, 580–583. [Google Scholar]
- Edwards, S.; Corry, B.; Kuyucak, S.; Chung, S.H. Continuum electrostatics fails to describe ion permeation in the gramicidin channel. Biophys. J. 2002, 83, 1348–1360. [Google Scholar]
- Andersen, O.S.; Koeppe, R.E.; Roux, B. Gramicidin channels. IEEE Trans. Nanobiosci. 2005, 4, 10–20. [Google Scholar]
- Separovic, F.; Cornell, B.A. Gated ion channel-based biosensor device. In Biological Membrane Ion Channels; Chung, S.H., Anderson, O.S., Krishnamurthy, V.V., Eds.; Springer: New York, NY, USA, 2007; pp. 595–621. [Google Scholar]
- Krishnamurthy, V.; Monfared, S.M.; Cornell, B. Ion channel biosensors—Part II: Dynamic modeling, analysis, and statistical signal processing. IEEE Trans. Nanotechnol. 2010, 9, 313–321. [Google Scholar]
- Lu, D.; Li, Y.; Ravaiolo, U.; Schulten, K. Empirical nanotube model for biological applications. J. Phys. Chem. B 2005, 109, 11461–11467. [Google Scholar]
- Gordon, D.; Krishnamurthy, V.; Chung, S.H. Generalized Langevin models of molecular dynamics simulations with applications to ion channels. J. Chem. Phys. 2009, 131, 134102. [Google Scholar]
- Hilder, T.A.; Gordon, D.; Chung, S.H. Synthetic chloride-selective carbon nanotubes examined by using molecular and stochastic dynamics. Biophys. J. 2010, 99, 1734–1742. [Google Scholar]
- Hilder, T.A.; Chung, S.H. Carbon nanotube as a gramicidin analogue. Chem. Phys. Lett. 2011, 501, 423–426. [Google Scholar]
- Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kale, L.; Schulten, K.J. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar]
- Humphrey, W.; Dalke, A.; Schulten, K.J. VMD: Visual molecular dynamics. J. Mol. Graphics 1996, 14, 33–38. [Google Scholar]
- Brooks, R.R.; Bruccoleri, R.E.; Olafson, B.D.; States, D.J.; Swaminathan, S.; Karplus, M. CHARMM: A program for macromolecular energy, minimization, and dynamics calculations. J. Comput. Chem. 1983, 4, 187–217. [Google Scholar]
- MacKerell, A.D., Jr.; Brooks, B.; Brooks, C.L., III; Nilsson, L.; Roux, B.; Won, Y.; Karplus, M. CHARMM: The energy function and its parameterization. In Encyclopedia of Computational Chemistry; Schleyer, P.V.R., Ed.; John Wiley & Sons: Chichester, UK, 1998; Volume 1, pp. 271–277. [Google Scholar]
- Fullerene Topology. Available online: http://www.ks.uiuc.edu/Research/namd/mailing_list/namd-l/11066.html (accessed on 26 March 2012).
- Brooks, B.R.; Brooks, C.L., III; Mackerell, A.D.; Nilsson, L.; Petrella, R.J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The biomolecular simulation program. J. Comput. Chem. 2009, 30, 1545–1614. [Google Scholar]
- Dunlap, B.I.; Brenner, D.W.; Mintmire, J.W.; Mowrey, R.C.; White, C.T. Geometric and electronic structures of C60H60, C60F60, and C60H36. J. Phys. Chem. 1991, 95, 5763–5768. [Google Scholar]
- Sheka, E.F. Stepwise computational synthesis of fullerene C60 derivatives. Fluorinated fullerenes C60F2k. J. Exp. Theor. Phys. 2010, 111, 397–414. [Google Scholar]
- Kumar, S.; Rosenberg, J.M.; Bouzida, D.; Swendsen, R.H.; Kollman, P.A. Multidimensional free-energy calculations using the weighted histogram analysis method. J. Comput. Chem. 1995, 16, 1339–1350. [Google Scholar]
- WHAM: The Weighted Histogram Analysis Method. Available online: http://membrane.urmc.rochester.edu/content/wham (accessed on 1 November 2011).
- Allen, T.W.; Andersen, O.S.; Roux, B. Energetics of ion conduction through the gramicidin channel. Proc. Natl. Acad. Sci. USA 2004, 101, 117–122. [Google Scholar]
- Chen, R.; Robinson, A.; Gordon, D.; Chung, S.H. Modeling the binding of three toxins to the voltage gated potassium channel (Kv1.3). Biophys. J. 2011, 101, 2652–2660. [Google Scholar]
- Ng, J.A.; Vora, T.; Krishnamurthy, V.; Chung, S.H. Estimating the dielectric constant of the channel protein and pore. Euro. Biophys. J. 2008, 37, 213–222. [Google Scholar]
- Hoyles, M.; Kuyucak, S.; Chung, S.H. Computer simulation of ion conductance in membrane channels. Phys. Rev. E 1998, 58, 3654–3661. [Google Scholar]
- Gu, Z.; Peng, H.; Hauge, R.H.; Smalley, R.E.; Margrave, J.L. Cutting single-wall carbon nanotubes through fluorination. Nano Lett. 2002, 2, 1009–1013. [Google Scholar]
- Taylor, R. Lecture Notes on Fullerene Chemistry; Imperial College Press: London, UK, 1999. [Google Scholar]
- Taylor, R. Why fluorinate fullerenes? J. Fluorine Chem. 2004, 125, 359–368. [Google Scholar]
- Bulusheva, L.G.; Okotrub, A.V.; Shnitov, V.V.; Bryzgalov, V.V.; Boltalina, O.V.; Goldt, I.V.; Vyalikh, D.V. Electronic structure of C60F36 studied by quantum-chemical modeling of experimental photoemission and X-ray absorption spectra. J. Chem. Phys. 2009, 130, 014704. [Google Scholar]
- Boltalina, O.V.; Galeva, N.A. Direct fluorination of fullerenes. Russ. Chem. Rev. 2000, 69, 609–621. [Google Scholar]
- Jia, J.; Wu, H.S.; Xu, X.H.; Zhang, X.M.; Jiao, H. Fused five-membered rings determine the stability of C60F60. J. Am. Chem. Soc. 2008, 130, 3985–3988. [Google Scholar]
- Holt, J.K.; Park, H.G.; Wang, Y.; Stadermann, M.; Artyukhin, A.B.; Grigoropoulos, C.P.; Noy, A.; Bakajin, O. Fast mass transport through sub-2-nanometer carbon nanotubes. Science 2006, 312, 1034–1037. [Google Scholar]
- Taylor, R.; Langley, G.J.; Brisdon, A.K.; Holloway, J.H.; Hope, E.G.; Kroto, H.W.; Walton, D.R.M. Highly oxygenated derivatives of fluorinated C60, and the mode of fragmentation of the fluorinated cage under electron impact ionization conditions. J. Chem. Soc. Chem. Commun. 1993, 10, 875–878. [Google Scholar]
- Tabata, Y.; Murakami, Y.; Ikada, Y. Photodynamic effect of polyethylene glycol-modified fullerene on tumor. Jpn. J. Cancer Res. 1997, 88, 1108–1116. [Google Scholar]
- Liu, J.; Ohta, S.; Sonoda, A.; Yamada, M.; Yamamoto, M.; Nitta, N.; Murata, K.; Tabata, Y. Preparation of PEG-conjugated fullerene containing Gd3+ ions for photodynamic therapy. J. Contr. Release 2007, 117, 104–110. [Google Scholar]
- Kawaguchi, H. Functional polymer microspheres. Prog. Polym. Sci. 2000, 25, 1171–1210. [Google Scholar]
- Uchida, T.; Hosaka, S.; Murao, Y. Complement activation by polymer binding IgG. Biomaterials 1984, 5, 281–283. [Google Scholar]
- Van Erp, R.; Linders, Y.E.M.; van Sommeren, A.P.G.; Gribnau, T.C.J. Characterization of monoclonal antibodies physically absorbed onto polystyrene latex particles. J. Immunol. Meth. 1992, 152, 191–199. [Google Scholar]
- Sarobe, J.; Miraballes, I.; Molina, J.A.; Forcada, J.; Hidalgo-Alvarez, R. Nephelometric assay of immunoglobulin G chemically bound to chloromethyl polystyrene beads. Polym. Adv. Technol. 1996, 7, 749–753. [Google Scholar]
- Chung, S.H.; Moore, J.B.; Xia, L.; Premkumar, L.S.; Gage, P.W. Characterization of single channel currents using digital signal processing techniques based on hidden Markov models. Phil. Trans. Roy. Soc. London B 1990, 329, 256–285. [Google Scholar]
- Chung, S.H.; Pulford, G.W. Fluctuation analysis of patch-clamp or whole-cell recordings containing many single channels. J. Neurosci. Meth. 1993, 50, 369–384. [Google Scholar]
- Pulford, G.W.; Kennedy, R.A.; Chung, S.H. Identification of individual channel kinetics from recordings containing many identical channels. Signal Proc. 1995, 43, 207–221. [Google Scholar]
- Liu, Y.; Roberts, D.C.S.; Morgan, D. Sensitization of the reinforcing effects of self-administered cocaine in rats: Effects of dose and intravenous injection speed. Euro. J. Neurosci. 2005, 22, 195–200. [Google Scholar]
- Friguet, B.; Chaffotte, A.F.; Djavadi-Ohaniance, L.; Goldberg, M.E. Measurements of the true affinity constant in solution of antigen-antibody complexes by enzyme-linked immunosorbent assay. J. Immunol. Meth. 1985, 77, 305–319. [Google Scholar]
- Nikitin, A.; Ogasawara, H.; Mann, D.; Denecke, R.; Zhang, Z.; Dai, H.; Cho, K.J.; Nilson, A. Hydrogenation of single-walled carbon nanotubes. Phys. Rev. Lett. 2005, 95, 225507. [Google Scholar]
- Raghuveer, M.S.; Kumar, A.; Frederick, M.J.; Louie, G.P.; Ganesan, P.G.; Ramanath, G. Site-selective functionalization of nano- or macroscopic segments of carbon nanotubes. Adv. Mater. 2006, 18, 547–552. [Google Scholar]
© 2012 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).
Share and Cite
Hilder, T.A.; Pace, R.J.; Chung, S.-H. Computational Design of a Carbon Nanotube Fluorofullerene Biosensor. Sensors 2012, 12, 13720-13735. https://doi.org/10.3390/s121013720
Hilder TA, Pace RJ, Chung S-H. Computational Design of a Carbon Nanotube Fluorofullerene Biosensor. Sensors. 2012; 12(10):13720-13735. https://doi.org/10.3390/s121013720
Chicago/Turabian StyleHilder, Tamsyn A., Ron J. Pace, and Shin-Ho Chung. 2012. "Computational Design of a Carbon Nanotube Fluorofullerene Biosensor" Sensors 12, no. 10: 13720-13735. https://doi.org/10.3390/s121013720
APA StyleHilder, T. A., Pace, R. J., & Chung, S. -H. (2012). Computational Design of a Carbon Nanotube Fluorofullerene Biosensor. Sensors, 12(10), 13720-13735. https://doi.org/10.3390/s121013720