Functionalized Fullerene Targeting Human Voltage-Gated Sodium Channel, hNav1.7.

Mutations of hNav1.7 that cause its activities to be enhanced contribute to severe neuropathic pain. Only a small number of hNav1.7 specific inhibitors have been identified, most of which interact with the voltage-sensing domain of the voltage-activated sodium ion channel. In our previous computational study, we demonstrated that a [Lys6]-C84 fullerene binds tightly (affinity of 46 nM) to NavAb, the voltage-gated sodium channel from the bacterium Arcobacter butzleri. Here, we extend this work and, using molecular dynamics simulations, demonstrate that the same [Lys6]-C84 fullerene binds strongly (2.7 nM) to the pore of a modeled human sodium ion channel hNav1.7. In contrast, the fullerene binds only weakly to a mutated model of hNav1.7 (I1399D) (14.5 mM) and a model of the skeletal muscle hNav1.4 (3.7 mM). Comparison of one representative sequence from each of the nine human sodium channel isoforms shows that only hNav1.7 possesses residues that are critical for binding the fullerene derivative and blocking the channel pore.

The Bound Structures of 17ß-Estradiol-Binding Aptamers.

DNA aptamers can exhibit high affinity and selectivity towards their targets, but the aptamer-target complex structures are rarely available from crystallography and often difficult to elucidate. This is particularly true of small molecule targets, including 17ß-estradiol (E2), which is becoming one of the most widely encountered endocrine-disrupting chemicals in the environment. Using molecular dynamics simulations, we demonstrate that E2 binds to a thymine loop region common to all E2-specific aptamers in the literature. Analyzing these structures allows us to design new E2 binding sequences. As well as illuminating the essential sequence and structural factors for generating specificity for E2, we demonstrate the effectiveness of molecular dynamics simulations for aptamer science.

Molecular Mechanism of Binding between 17ß-Estradiol and DNA.

Although 17ß-estradiol (E2) is a natural molecule involved in the endocrine system, its widespread use in various applications has resulted in its accumulation in the environment and its classification as an endocrine-disrupting molecule. These molecules can interfere with the hormonal system, and have been linked to various adverse effects such as the proliferation of breast cancer. It has been proposed that E2 could contribute to breast cancer by the induction of DNA damage. Mass spectrometry has demonstrated that E2 can bind to DNA but the mechanism by which E2 interacts with DNA has yet to be elucidated. Using all-atom molecular dynamics simulations, we demonstrate that E2 intercalates (inserts between two successive DNA base pairs) in DNA at the location specific to estrogen receptor binding, known as the estrogen response element (ERE), and to other random sequences of DNA. Our results suggest that excess E2 has the potential to disrupt processes in the body which rely on binding to DNA, such as the binding of the estrogen receptor to the ERE and the activity of enzymes that bind DNA, and could lead to DNA damage.

Interaction of Boron Nitride Nanosheets with Model Cell Membranes.

Boron nitride nanomaterials have attracted attention for biomedical applications, due to their improved biocompatibility when compared with carbon nanomaterials. Recently, graphene and graphene oxide nanosheets have been shown, both experimentally and computationally, to destructively extract phospholipids from Escherichia coli. Boron nitride nanosheets (BNNSs) have exciting potential biological and environmental applications, for example the ability to remove oil from water. These applications are likely to increase the exposure of prokaryotes and eukaryotes to BNNSs. Yet, despite their promise, the interaction between BNNSs and cell membranes has not yet been investigated. Here, all-atom molecular dynamics simulations were used to demonstrate that BNNSs are spontaneously attracted to the polar headgroups of the lipid bilayer. The BNNSs do not passively cross the lipid bilayer, most likely due to the large forces experienced by the BNNSs. This study provides insight into the interaction of BNNSs with cell membranes and may aid our understanding of their improved biocompatibility.

Insertion mechanism and stability of boron nitride nanotubes in lipid bilayers.

We provide insight into the interaction of boron nitride nanotubes (BNNTs) with cell membranes to better understand their improved biocompatibility compared to carbon nanotubes (CNTs). Contrary to CNTs, no computational studies exist investigating the insertion mechanism and stability of BNNTs in membranes. Our molecular dynamics simulations demonstrate that BNNTs are spontaneously attracted to lipid bilayers and are stable once inserted. They insert via a lipid-mediated, passive insertion mechanism. BNNTs demonstrate similar characteristics to more biocompatible functionalized CNTs.

Binding of fullerenes and nanotubes to MscL.

Multi-drug resistance is becoming an increasing problem in the treatment of bacterial infections and diseases. The mechanosensitive channel of large conductance (MscL) is highly conserved among prokaryotes. Evidence suggests that a pharmacological agent that can affect the gating of, or block the current through, MscL has significant potential as a new class of antimicrobial compound capable of targeting a range of pathogenic bacteria with minimal side-effects to infected patients. Using molecular dynamics we examine the binding of fullerenes and nanotubes to MscL and demonstrate that both are stable within the MscL pore. We predict that fullerenes will attenuate the flow of ions through MscL by reducing the pore volume available to water and ions, but nanotubes will prevent pore closure resulting in a permanently open pore. Moreover, we confirm experimentally that it is possible to attenuate the flow of ions through MscL using a C60-γ cyclodextrin complex.

What have we learnt about the mechanisms of rapid water transport, ion rejection and selectivity in nanopores from molecular simulation?

Nanopores have demonstrated an extraordinary ability to allow water molecules to pass through their interiors at rates far exceeding expectations based on continuum theory. Moreover, simulation studies suggest that particular nanoscale pores have the potential to discriminate between water and salts as well as to distinguish between a range of different ion types. Some of the unusual features of transport in these nanopores have been elucidated with molecular dynamics simulation, specifically the spontaneous filling and rapid transport of water, the rejection of ions and the selection between ions. The main focus of this review, however, is the physical mechanisms which act to produce such remarkable behaviour at this scale, drawing on the many studies that have been conducted in the last decade. Since molecular dynamics simulations allow the motion of individual atoms to be followed over time, they have the potential to provide fundamental insight into the reasons why transport in nanoscale pores differs from expectations based on macroscopic theory. Gaining an understanding of the mechanisms of transport in these tiny pores should guide future experiments in this area aimed at developing novel technologies and improving existing membrane separation techniques.

Multi-ion versus single-ion conduction mechanisms can yield current rectification in biological ion channels.

There is clear evidence that the net magnitude of negative charge at the intracellular end of inwardly rectifying potassium channels helps to generate an asymmetry in the magnitude of the current that will pass in each direction. However, a complete understanding of the physical mechanism that links these charges to current rectification has yet to be obtained. Using Brownian dynamics, we compare the conduction mechanism and binding sites in rectifying and non-rectifying channel models. We find that in our models, rectification is a consequence of asymmetry in the hydrophobicity and charge of the pore lining. As a consequence, inward conduction can occur by a multi-ion conduction mechanism. However, outward conduction is restricted, since there are fewer ions at the intracellular entrance and outwardly moving ions must cross the pore on their own. We pose the question as to whether the same mechanism could be at play in inwardly rectifying potassium channels.

Designing a C84 fullerene as a specific voltage-gated sodium channel blocker.

Fullerene derivatives demonstrate considerable potential for numerous biological applications, such as the effective inhibition of HIV protease. Recently, they were identified for their ability to indiscriminately block biological ion channels. A fullerene derivative which specifically blocks a particular ion channel could lead to a new set of drug leads for the treatment of various ion channel-related diseases. Here, we demonstrate their extraordinary potential by designing a fullerene which mimics some of the functions of µ-conotoxin, a peptide derived from cone snail venom which potently binds to the bacterial voltage-gated sodium channel (NavAb). We show, using molecular dynamics simulations, that the C84 fullerene with six lysine derivatives uniformly attached to its surface is selective to NavAb over a voltage-gated potassium channel (Kv1.3). The side chain of one of the lysine residues protrudes into the selectivity filter of the channel, while the methionine residues located just outside of the channel form hydrophobic contacts with the carbon atoms of the fullerene. The modified C84 fullerene strongly binds to the NavAb channel with an affinity of 46 nM but binds weakly to Kv1.3 with an affinity of 3 mM. This potent blocker of NavAb may serve as a structural template from which potent compounds can be designed for the targeting of mammalian Nav channels. There is a genuine need to target mammalian Nav channels as a form of treatment of various diseases which have been linked to their malfunction, such as epilepsy and chronic pain.

Conduction and block of inward rectifier K+ channels: predicted structure of a potent blocker of Kir2.1.

Dysfunction of Kir2.1, thought to be the major component of inward currents, I(K1), in the heart, has been linked to various channelopathies, such as short Q-T syndrome. Unfortunately, currently no known blockers of Kir2.x channels exist. In contrast, Kir1.1b, predominantly expressed in the kidney, is potently blocked by an oxidation-resistant mutant of the honey bee toxin tertiapin (tertiapin-Q). Using various computational tools, we show that both channels are closed by a hydrophobic gating mechanism and inward rectification occurs in the absence of divalent cations and polyamines. We then demonstrate that tertiapin-Q binds to the external vestibule of Kir1.1b and Kir2.1 with K(d) values of 11.6 nM and 131 µM, respectively. We find that a single mutation of tertiapin-Q increases the binding affinity for Kir2.1 by 5 orders of magnitude (K(d) = 0.7 nM). This potent blocker of Kir2.1 may serve as a structural template from which potent compounds for the treatment of various diseases mediated by this channel subfamily, such as cardiac arrhythmia, can be developed.

Conductance properties of the inwardly rectifying channel, Kir3.2: molecular and Brownian dynamics study.

Using the recently unveiled crystal structure, and molecular and Brownian dynamics simulations, we elucidate several conductance properties of the inwardly rectifying potassium channel, Kir3.2, which is implicated in cardiac and neurological disorders. We show that the pore is closed by a hydrophobic gating mechanism similar to that observed in Kv1.2. Once open, potassium ions move into, but not out of, the cell. The asymmetrical current-voltage relationship arises from the lack of negatively charged residues at the narrow intracellular mouth of the channel. When four phenylalanine residues guarding the intracellular gate are mutated to glutamate residues, the channel no longer shows inward rectification. Inward rectification is restored in the mutant Kir3.2 when it becomes blocked by intracellular Mg(2+). Tertiapin, a polypeptide toxin isolated from the honey bee, is known to block several subtypes of the inwardly rectifying channels with differing affinities. We identify critical residues in the toxin and Kir3.2 for the formation of the stable complex. A lysine residue of tertiapin protrudes into the selectivity filter of Kir3.2, while two other basic residues of the toxin form hydrogen bonds with acidic residues located just outside the channel entrance. The depth of the potential of mean force encountered by tertiapin is -16.1kT, thus indicating that the channel will be half-blocked by 0.4µM of the toxin.

Computational design of a carbon nanotube fluorofullerene biosensor.

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(TM) 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.

Computational modeling of transport in synthetic nanotubes.

Synthetic nanotubes that have the ability to broadly mimic the function of biological ion channels have extraordinary potential for various applications, from ultrasensitive biosensors to efficient water purification devices. As a result of their immense potential, the design and fabrication of such synthetic nanotubes is rapidly gaining momentum. We briefly review recent theoretical and experimental studies on nanoscale cylindrical hollow tubes constructed from carbon, boron, and nitrogen atoms that are able to selectively transport water molecules, cations (positively charged ions), or anions (negatively charged ions) similar to various biological ion channels. FROM THE CLINICAL EDITOR: This review discusses the current status of synthetic nanotube research, including recent theoretical and experimental studies on nanoscale cylindrical hollow tubes constructed from carbon, boron, and nitrogen atoms that are able to selectively transport water molecules, cations or anions similar to biological ion channels.

Synthetic cation-selective nanotube: permeant cations chaperoned by anions.

The ability to design ion-selective, synthetic nanotubes which mimic biological ion channels may have significant implications for the future treatment of bacteria, diseases, and as ultrasensitive biosensors. We present the design of a synthetic nanotube made from carbon atoms that selectively allows monovalent cations to move across and rejects all anions. The cation-selective nanotube mimics some of the salient properties of biological ion channels. Before practical nanodevices are successfully fabricated it is vital that proof-of-concept computational studies are performed. With this in mind we use molecular and stochastic dynamics simulations to characterize the dynamics of ion permeation across a single-walled (10, 10), 36 Å long, carbon nanotube terminated with carboxylic acid with an effective radius of 5.08 Å. Although cations encounter a high energy barrier of 7 kT, its height is drastically reduced by a chloride ion in the nanotube. The presence of a chloride ion near the pore entrance thus enables a cation to enter the pore and, once in the pore, it is chaperoned by the resident counterion across the narrow pore. The moment the chaperoned cation transits the pore, the counterion moves back to the entrance to ferry another ion. The synthetic nanotube has a high sodium conductance of 124 pS and shows linear current-voltage and current-concentration profiles. The cation-anion selectivity ratio ranges from 8 to 25, depending on the ionic concentrations in the reservoirs.

Synthetic chloride-selective carbon nanotubes examined by using molecular and stochastic dynamics.

Synthetic channels, such as nanotubes, offer the possibility of ion-selective nanoscale pores which can broadly mimic the functions of various biological ion channels, and may one day be used as antimicrobial agents, or for treatment of cystic fibrosis. We have designed a carbon nanotube that is selectively permeable to anions. The virtual nanotubes are constructed from a hexagonal array of carbon atoms (graphene) rolled up to form a tubular structure, with an effective radius of 4.53 Å and length of 34 Å. The pore ends are terminated with polar carbonyl groups. The nanotube thus formed is embedded in a lipid bilayer and a reservoir containing ionic solutions is added at each end of the pore. The conductance properties of these synthetic channels are then examined with molecular and stochastic dynamics simulations. Profiles of the potential of mean force at 0 mM reveal that a cation moving across the pore encounters an insurmountable free energy barrier of ∼25 kT in height. In contrast, for anions, there are two energy wells of ∼12 kT near each end of the tube, separated by a central free energy barrier of 4 kT. The conductance of the pore, with symmetrical 500 mM solutions in the reservoirs, is 72 pS at 100 mV. The current saturates with an increasing ionic concentration, obeying a Michaelis-Menten relationship. The pore is normally occupied by two ions, and the rate-limiting step in conduction is the time taken for the resident ion near the exit gate to move out of the energy well.

Boron nitride nanotubes selectively permeable to cations or anions.

Biological ion channels in membranes are selectively permeable to specific ionic species. They maintain the resting membrane potential, generate propagated action potentials, and control a wide variety of cell functions. Here it is demonstrated theoretically that boron nitride nanotubes have the ability to carry out some of the important functions of biological ion channels. Boron nitride nanotubes with radii of 4.83 and 5.52 A embedded in a silicon nitride membrane are selectively permeable to cations and anions, respectively. They broadly mimic some of the permeation characteristics of gramicidin and chloride channels. Using distributional molecular dynamics, which is a combination of molecular and stochastic dynamics simulations, the properties of these engineered nanotubes are characterized, such as the free energy encountered by charged particles, the water-ion structure within the pore, and the current-voltage and current-concentration profiles. These engineered nanotubes have potential applications as sensitive biosensors, antibiotics, or filtration devices.

Salt rejection and water transport through boron nitride nanotubes.

Nanotube-based water-purification devices have the potential to transform the field of desalination and demineralization through their ability to remove salts and heavy metals without significantly affecting the fast flow of water molecules. Boron nitride nanotubes have shown superior water flow properties compared to carbon nanotubes, and are thus expected to provide a more efficient water purification device. Using molecular dynamics simulations it is shown that a (5, 5) boron nitride nanotube embedded in a silicon nitride membrane can, in principle, obtain 100% salt rejection at concentrations as high as 1 M owing to a high energy barrier while still allowing water molecules to flow at a rate as high as 10.7 water molecules per nanosecond (or 0.9268 L m(-2) h(-1)). Furthermore, ions continue to be rejected under the influence of high hydrostatic pressures up to 612 MPa. When the nanotube radius is increased to 4.14 A the tube becomes cation-selective, and at 5.52 A the tube becomes anion-selective.

Maximum velocity for a single water molecule entering a carbon nanotube.

Carbon nanotubes, despite their hydrophobic nature, rapidly fill with water and allow super fast fluid flow through their interior due to the almost frictionless nanotube surface. The question arises as to whether it is possible to maximize the uptake (suction energy) of water and thus generate the highest possible fluid flow. In this paper, we outline the concepts of an acceptance condition and the suction energy and subsequently examine the suction characteristics of a single water molecule entering a carbon nanotube. In particular, we find that for the hydrogen atoms oriented towards the tube end, the radius of the carbon nanotube must be at least 3.464 A (or 0.3464 nm) for acceptance of a water molecule, and that a radius of 3.95 A provides the maximum uptake or suction energy.

Modeling the loading and unloading of drugs into nanotubes.

One of the most promising applications of nanotechnology is that of drug delivery, and in particular the targeted delivery of drugs using nanotubes. Functionalized nanotubes might be able to target specific cells, become ingested, and then release their contents in response to a chemical trigger. This will have significant implications for the future treatment of patients, particularly those suffering from cancer, for whom presently the nonspecific nature of chemotherapy often kills healthy normal cells. Research to date has largely been through experiments investigating toxicity, biocompatibility, solubility, functionalization, and cellular uptake. More recently, the loading and unloading of molecular cargo has gained momentum from both experimental and theoretical investigations. This Review focuses on the loading and unloading of molecular cargo and highlights recent theoretical investigations, which to date have received very little attention in the review literature. The development of nanotube drug-delivery capsules is of vital concern for the improvement of medical treatment, and mathematical modeling tends to facilitate such development and provides a quicker route to applications of the technology. This Review highlights the latest progress in terms of theoretical investigations and provides a focus for the development of the next generation of medical therapeutics.