Mechanical activation opens a lipid-lined pore in OSCA ion channels.

OSCA/TMEM63 channels are the largest known family of mechanosensitive channels1-3, playing critical roles in plant4-7 and mammalian8,9 mechanotransduction. Here we determined 44 cryogenic electron microscopy structures of OSCA/TMEM63 channels in different environments to investigate the molecular basis of OSCA/TMEM63 channel mechanosensitivity. In nanodiscs, we mimicked increased membrane tension and observed a dilated pore with membrane access in one of the OSCA1.2 subunits. In liposomes, we captured the fully open structure of OSCA1.2 in the inside-in orientation, in which the pore shows a large lateral opening to the membrane. Unusually for ion channels, structural, functional and computational evidence supports the existence of a 'proteo-lipidic pore' in which lipids act as a wall of the ion permeation pathway. In the less tension-sensitive homologue OSCA3.1, we identified an 'interlocking' lipid tightly bound in the central cleft, keeping the channel closed. Mutation of the lipid-coordinating residues induced OSCA3.1 activation, revealing a conserved open conformation of OSCA channels. Our structures provide a global picture of the OSCA channel gating cycle, uncover the importance of bound lipids and show that each subunit can open independently. This expands both our understanding of channel-mediated mechanotransduction and channel pore formation, with important mechanistic implications for the TMEM16 and TMC protein families.

TLR7 gain-of-function genetic variation causes human lupus.

Although circumstantial evidence supports enhanced Toll-like receptor 7 (TLR7) signalling as a mechanism of human systemic autoimmune disease1-7, evidence of lupus-causing TLR7 gene variants is lacking. Here we describe human systemic lupus erythematosus caused by a TLR7 gain-of-function variant. TLR7 is a sensor of viral RNA8,9 and binds to guanosine10-12. We identified a de novo, previously undescribed missense TLR7Y264H variant in a child with severe lupus and additional variants in other patients with lupus. The TLR7Y264H variant selectively increased sensing of guanosine and 2',3'-cGMP10-12, and was sufficient to cause lupus when introduced into mice. We show that enhanced TLR7 signalling drives aberrant survival of B cell receptor (BCR)-activated B cells, and in a cell-intrinsic manner, accumulation of CD11c+ age-associated B cells and germinal centre B cells. Follicular and extrafollicular helper T cells were also increased but these phenotypes were cell-extrinsic. Deficiency of MyD88 (an adaptor protein downstream of TLR7) rescued autoimmunity, aberrant B cell survival, and all cellular and serological phenotypes. Despite prominent spontaneous germinal-centre formation in Tlr7Y264H mice, autoimmunity was not ameliorated by germinal-centre deficiency, suggesting an extrafollicular origin of pathogenic B cells. We establish the importance of TLR7 and guanosine-containing self-ligands for human lupus pathogenesis, which paves the way for therapeutic TLR7 or MyD88 inhibition.

Mechanistic basis for multidrug resistance and collateral drug sensitivity conferred to the malaria parasite by polymorphisms in PfMDR1 and PfCRT.

Polymorphisms in the Plasmodium falciparum multidrug resistance protein 1 (pfmdr1) gene and the Plasmodium falciparum chloroquine resistance transporter (pfcrt) gene alter the malaria parasite's susceptibility to most of the current antimalarial drugs. However, the precise mechanisms by which PfMDR1 contributes to multidrug resistance have not yet been fully elucidated, nor is it understood why polymorphisms in pfmdr1 and pfcrt that cause chloroquine resistance simultaneously increase the parasite's susceptibility to lumefantrine and mefloquine-a phenomenon known as collateral drug sensitivity. Here, we present a robust expression system for PfMDR1 in Xenopus oocytes that enables direct and high-resolution biochemical characterizations of the protein. We show that wild-type PfMDR1 transports diverse pharmacons, including lumefantrine, mefloquine, dihydroartemisinin, piperaquine, amodiaquine, methylene blue, and chloroquine (but not the antiviral drug amantadine). Field-derived mutant isoforms of PfMDR1 differ from the wild-type protein, and each other, in their capacities to transport these drugs, indicating that PfMDR1-induced changes in the distribution of drugs between the parasite's digestive vacuole (DV) and the cytosol are a key driver of both antimalarial resistance and the variability between multidrug resistance phenotypes. Of note, the PfMDR1 isoforms prevalent in chloroquine-resistant isolates exhibit reduced capacities for chloroquine, lumefantrine, and mefloquine transport. We observe the opposite relationship between chloroquine resistance-conferring mutations in PfCRT and drug transport activity. Using our established assays for characterizing PfCRT in the Xenopus oocyte system and in live parasite assays, we demonstrate that these PfCRT isoforms transport all 3 drugs, whereas wild-type PfCRT does not. We present a mechanistic model for collateral drug sensitivity in which mutant isoforms of PfMDR1 and PfCRT cause chloroquine, lumefantrine, and mefloquine to remain in the cytosol instead of sequestering within the DV. This change in drug distribution increases the access of lumefantrine and mefloquine to their primary targets (thought to be located outside of the DV), while simultaneously decreasing chloroquine's access to its target within the DV. The mechanistic insights presented here provide a basis for developing approaches that extend the useful life span of antimalarials by exploiting the opposing selection forces they exert upon PfCRT and PfMDR1.

Nav1.2 channel mutations preventing fast inactivation lead to SCN2A encephalopathy.

SCN2A gene-related early-infantile developmental and epileptic encephalopathy (EI-DEE) is a rare and severe disorder that manifests in early infancy. SCN2A mutations affecting the fast inactivation gating mechanism can result in altered voltage dependence and incomplete inactivation of the encoded neuronal Nav1.2 channel and lead to abnormal neuronal excitability. In this study, we evaluated clinical data of seven missense Nav1.2 variants associated with DEE and performed molecular dynamics simulations, patch-clamp electrophysiology, and dynamic clamp real-time neuronal modelling to elucidate the molecular and neuron-scale phenotypic consequences of the mutations. The N1662D mutation almost completely prevented fast inactivation without affecting activation. The comparison of wild-type and N1662D channel structures suggested that the ambifunctional hydrogen bond formation between residues N1662 and Q1494 is essential for fast inactivation. Fast inactivation could also be prevented with engineered Q1494A or Q1494L Nav1.2 channel variants, whereas Q1494E or Q1494 K variants resulted in incomplete inactivation and persistent current. Molecular dynamics simulations revealed a reduced affinity of the hydrophobic IFM-motif to its receptor site with N1662D and Q1494L variants relative to wild-type. These results demonstrate that the interactions between N1662 and Q1494 underpin the stability and the orientation of the inactivation gate and are essential for the development of fast inactivation. Six DEE-associated Nav1.2 variants, with mutations mapped to channel segments known to be implicated in fast inactivation were also evaluated. Remarkably, the L1657P variant also prevented fast inactivation and produced biophysical characteristics that were similar to those of N1662D, whereas the M1501 V, M1501T, F1651C, P1658S, and A1659 V variants resulted in biophysical properties that were consistent with gain-of-function and enhanced action potential firing of hybrid neurons in dynamic action potential clamp experiments. Paradoxically, low densities of N1662D or L1657P currents potentiated action potential firing, whereas increased densities resulted in sustained depolarization. Our results provide novel structural insights into the molecular mechanism of Nav1.2 channel fast inactivation and inform treatment strategies for SCN2A-related EI-DEE. The contribution of non-inactivating Nav1.2 channels to neuronal excitability may constitute a distinct cellular mechanism in the pathogenesis of SCN2A-related DEE.

Lipid redistribution in the highly curved footprint of Piezo1.

We investigate the effects on the distribution of lipids in the plasma membrane that are caused by the insertion of a protein, Piezo1, that significantly distorts the membrane toward the cytosol. From coarse-grained molecular dynamics simulations, we find that the major effects occur in the outer, extracellular, leaflet. The mol fraction of cholesterol increases significantly in the curved region of the membrane close to Piezo1, while those of phosphatidylcholine and of sphingomyelin decrease. In the inner leaflet, mol fractions of cholesterol and of phosphatidylethanolamine decrease slightly as the protein is approached, while that of phosphatidylserine increases slightly. The mol fraction of phosphatidylcholine decreases markedly as the protein is approached. Most of these results are understood in the context of a theoretical model that utilizes two elements: 1) a coupling between the leaflets' actual curvatures and their compositionally dependent spontaneous curvatures and 2) the dependence of the spontaneous curvatures not only on the mol fractions of the phospholipids, but also on the effect that cholesterol has on the spontaneous curvatures of the phospholipids.

New insights from modelling studies and molecular dynamics simulations of the DIS5-S6 extracellular linker of the skeletal muscle sodium channel NaV1.4.

In the CryoEM-structure of the hSkMNaV1.4 ion channel (PDB:6AGF), the 59-residue DIS5-S6 linker peptide was omitted due to absence of electron density. This peptide is intriguing - comprised of unique sequence and found only in mammalian skeletal muscle sodium ion channels. To probe potential physiological and evolutionary significance, we constructed an homology model of the complete hSkMNaV1.4 channel. Rather than a flexible random coil potentiating drift across the channel, the linker folds into a compact configuration through self-assembling secondary structural elements. Analogous sequences from 48 mammalian organisms show hypervariability with between 40% and 100% sequence similarity. To investigate structural implications, sequences from 14 representative organisms were additionally modelled. All showed highly conserved N-and C-terminal residues closely superimposed, suggesting a critical functional role. An optimally located asparagine residue within the conserved region was investigated for N-linked glycosylation and MD simulations carried out. Results suggest a complex glycan added at this site in the linker may form electrostatic interactions with the DIV voltage sensing domain and be mechanistically involved in channel gating. The relationship of unique sequence, compact configuration, potential glycosylation and MD simulations are discussed relative to SkMNaV1.4 structure and function.

Characterizing fenestration size in sodium channel subtypes and their accessibility to inhibitors.

Voltage-gated sodium channels (Nav) underlie the electrical activity of nerve and muscle cells. Humans have nine different subtypes of these channels, which are the target of small-molecule inhibitors commonly used to treat a range of conditions. Structural studies have identified four lateral fenestrations within the Nav pore module that have been shown to influence Nav pore blocker access during resting-state inhibition. However, the structural differences among the nine subtypes are still unclear. In particular, the dimensions of the four individual fenestrations across the Nav subtypes and their differential accessibility to pore blockers is yet to be characterized. To address this, we applied classical molecular dynamics simulations to study the recently published structures of Nav1.1, Nav1.2, Nav1.4, Nav1.5, and Nav1.7. Although there is significant variability in the bottleneck sizes of the Nav fenestrations, the subtypes follow a common pattern, with wider DI-II and DIII-IV fenestrations, a more restricted DII-III fenestration, and the most restricted DI-IV fenestration. We further identify the key bottleneck residues in each fenestration and show that the motions of aromatic residue sidechains govern the bottleneck radii. Well-tempered metadynamics simulations of Nav1.4 and Nav1.5 in the presence of the pore blocker lidocaine also support the DI-II fenestration being the most likely access route for drugs. Our computational results provide a foundation for future in vitro experiments examining the route of drug access to sodium channels. Understanding the fenestrations and their accessibility to drugs is critical for future analyses of diseases mutations across different sodium channel subtypes, with the potential to inform pharmacological development of resting-state inhibitors and subtype-selective drug design.

Modeling thermodiffusion in aqueous sodium chloride solutions-Which water model is best?

Thermodiffusion is the migration of a species due to a temperature gradient and is the driving phenomenon in many applications ranging from early cancer detection to uranium enrichment. Molecular dynamics (MD) simulations can be a useful tool for exploring the rather complex thermodiffusive behavior of species, such as proteins and ions. However, current MD models of thermodiffusion in aqueous ionic solutions struggle to quantitatively predict the Soret coefficient, which indicates the magnitude and direction of species migration under a temperature gradient. In this work, we aim to improve the accuracy of MD thermodiffusion models by assessing how well different water models can recreate thermodiffusion in a benchmark aqueous NaCl solution. We tested four of the best available rigid non-polarizable water models (TIP3P-FB, TIP4P-FB, OPC3, and OPC) and the commonly used TIP3P and SPC/E water models for their ability to predict the inversion temperature and Soret coefficient in 0.5, 2, and 4M aqueous NaCl solutions. Each water model predicted a noticeably different ion distribution yielding different inversion temperatures and magnitudes of the Soret coefficient. By comparing the modeled Soret coefficients to published experimental values, we determine TIP3P-FB to be the water model that best recreates thermodiffusion in aqueous NaCl solutions. Our findings can aid future works in selecting the most accurate rigid non-polarizable water model, including water and ion parameters for investigating thermodiffusion through MD simulations.

Differences in local anaesthetic and antiepileptic binding in the inactivated state of human sodium channel Nav1.4.

Voltage-gated sodium channels play a vital role in nerve and muscle cells, enabling them to encode and transmit electrical signals. Currently, there exist several classes of drugs that aim to inhibit these channels for therapeutic purposes, including local anesthetics, antiepileptics and antiarrhythmics. However, sodium-channel-inhibiting drugs lack subtype specificity; instead, they inhibit all sodium channels in the human body. Improving understanding of the mechanisms of binding of existing nonselective drugs is important in providing insight into how subtype-selective drugs could be developed. This study used molecular dynamics simulations to investigate the binding of the antiepileptics carbamazepine and lamotrigine and the local anesthetic lidocaine in neutral and charged states to the recently resolved human Nav1.4 channel. Replica exchange solute tempering was used to enable greater sampling of each compound within the pore. It was found that all four compounds show similarities in their binding sites within the pore. However, the positions of the carbamazepine and lamotrigine did not occlude the center of the pore but preferentially bound to homologous domain DII and DIII. The charged and neutral forms of lidocaine positioned themselves more centrally in the pore, with more common interactions with DIV. The best localized binding site was for charged lidocaine, whose aromatic moiety interacted with Y1593, whereas the amine projected toward the selectivity filter. Comparisons with our previous simulations and published structures highlight potential differences between tonic and use-dependent block related to conformational changes occurring in the pore.

Binding of azobenzene and p-diaminoazobenzene to the human voltage-gated sodium channel Nav1.4.

The activity of voltage-gated ion channels can be controlled by the binding of photoswitches inside their internal cavity and subsequent light irradiation. We investigated the binding of azobenzene and p-diaminoazobenzene to the human Nav1.4 channel in the inactivated state by means of Gaussian accelerated molecular dynamics simulations and free-energy computations. Three stable binding pockets were identified for each of the two photoswitches. In all the cases, the binding is controlled by the balance between the favorable hydrophobic interactions of the ligands with the nonpolar residues of the protein and the unfavorable polar solvation energy. In addition, electrostatic interactions between the ligand and the polar aminoacids are also relevant for p-diaminoazobenzene due to the presence of the amino groups on the benzene moieties. These groups participate in hydrogen bonding in the most favorable binding pocket and in long-range electrostatic interactions in the other pockets. The thermodinamically preferred binding sites found for both photoswitches are close to the selectivity filter of the channel. Therefore, it is very likely that the binding of these ligands will induce alterations in the ion conduction through the channel.

Protonation state of inhibitors determines interaction sites within voltage-gated sodium channels.

Voltage-gated sodium channels are essential for carrying electrical signals throughout the body, and mutations in these proteins are responsible for a variety of disorders, including epilepsy and pain syndromes. As such, they are the target of a number of drugs used for reducing pain or combatting arrhythmias and seizures. However, these drugs affect all sodium channel subtypes found in the body. Designing compounds to target select sodium channel subtypes will provide a new therapeutic pathway and would maximize treatment efficacy while minimizing side effects. Here, we examine the binding preferences of nine compounds known to be sodium channel pore blockers in molecular dynamics simulations. We use the approach of replica exchange solute tempering (REST) to gain a more complete understanding of the inhibitors' behavior inside the pore of NavMs, a bacterial sodium channel, and NavPas, a eukaryotic sodium channel. Using these simulations, we are able to show that both charged and neutral compounds partition into the bilayer, but neutral forms more readily cross it. We show that there are two possible binding sites for the compounds: (i) a site on helix 6, which has been previously determined by many experimental and computational studies, and (ii) an additional site, occupied by protonated compounds in which the positively charged part of the drug is attracted into the selectivity filter. Distinguishing distinct binding poses for neutral and charged compounds is essential for understanding the nature of pore block and will aid the design of subtype-selective sodium channel inhibitors.

Piezo1 Forms Specific, Functionally Important Interactions with Phosphoinositides and Cholesterol.

Touch, hearing, and blood pressure regulation require mechanically gated ion channels that convert mechanical stimuli into electrical currents. One such channel is Piezo1, which plays a key role in the transduction of mechanical stimuli in humans and is implicated in diseases, such as xerocytosis and lymphatic dysplasia. There is building evidence that suggests Piezo1 can be regulated by the membrane environment, with the activity of the channel determined by the local concentration of lipids, such as cholesterol and phosphoinositides. To better understand the interaction of Piezo1 with its environment, we conduct simulations of the protein in a complex mammalian bilayer containing more than 60 different lipid types together with electrophysiology and mutagenesis experiments. We find that the protein alters its local membrane composition, enriching specific lipids and forming essential binding sites for phosphoinositides and cholesterol that are functionally relevant and often related to Piezo1-mediated pathologies. We also identify a number of key structural connections between the propeller and pore domains located close to lipid-binding sites.

FRET-Modulated Multihybrid Nanoparticles for Brightness-Equalized Single-Wavelength Barcoding.

Semiconductor quantum dots (QDs) are the most versatile fluorophores for Förster resonance energy transfer (FRET) because they can function as both donors and acceptors for a multitude of fluorophores. However, a complete understanding of multidonor-multiacceptor FRET networks on QDs and their full employment into advanced fluorescence sensing and imaging have not been accomplished. Here, we provide a holistic photophysical analysis of such multidonor-QD-multiacceptor FRET systems using time-resolved and steady-state photoluminescence (PL) spectroscopy and Monte Carlo simulations. Multiple terbium complex (Tb) donors (1-191 units) and Cy5.5 dye acceptors (1-60 units) were attached to a central QD, and the entire range of combinations of FRET pathways was investigated by Tb, QD, and Cy5.5 PL. Experimental and simulation results were in excellent agreement and could disentangle the distinct contributions of hetero-FRET, homo-FRET, and dye dimerization. The FRET efficiency was independent of the number of Tb donors and dependent on the number of Cy5.5 acceptors, which could be used to independently adapt the PL intensity by the number of Tb donors and the PL lifetime by the number of Cy5.5 acceptors. We used this unique tuning capability to prepare Tb-QD-Cy5.5 conjugates with distinct QD PL lifetimes but similar QD PL intensities. These brightness-equalized multihybrid FRET nanoparticles were applied to optical barcoding via three time-gated PL intensity detection windows, which resulted in simple RGB ratios. Direct applicability was demonstrated by an efficient RGB distinction of different nanoparticle-encoded microbeads within the same field of view with both single-wavelength excitation and detection on a standard fluorescence microscope.

A computational assessment of the permeability and salt rejection of carbon nanotube membranes and their application to water desalination.

Membranes made from nanomaterials such as nanotubes and graphene have been suggested to have a range of applications in water filtration and desalination, but determining their suitability for these purposes requires an accurate assessment of the properties of these novel materials. In this study, we use molecular dynamics simulations to determine the permeability and salt rejection capabilities for membranes incorporating carbon nanotubes (CNTs) at a range of pore sizes, pressures and concentrations. We include the influence of osmotic gradients and concentration build up and simulate at realistic pressures to improve the reliability of estimated membrane transport properties. We find that salt rejection is highly dependent on the applied hydrostatic pressure, meaning high rejection can be achieved with wider tubes than previously thought; while membrane permeability depends on salt concentration. The ideal size of the CNTs for desalination applications yielding high permeability and high salt rejection is found to be around 1.1 nm diameter. While there are limited energy gains to be achieved in using ultra-permeable CNT membranes in desalination by reverse osmosis, such membranes may allow for smaller plants to be built as is required when size or weight must be minimized. There are diminishing returns in further increasing membrane permeability, so efforts should focus on the fabrication of membranes containing narrow or functionalized CNTs that yield the desired rejection or selection properties rather than trying to optimize pore densities.

A functional correlate of severity in alternating hemiplegia of childhood.

OBJECTIVE: Mutations in ATP1A3, the gene that encodes the α3 subunit of the Na(+)/K(+) ATPase, are the primary cause of alternating hemiplegia of childhood (AHC). Correlations between different mutations and AHC severity were recently reported, with E815K identified in severe and D801N and G947R in milder cases. This study aims to explore the molecular pathological mechanisms in AHC and to identify functional correlates for mutations associated with different levels of disease severity. METHODS: Human wild type ATP1A3, and E815K, D801N and G947R mutants were expressed in Xenopus laevis oocytes and Na(+)/K(+) ATPase function measured. Structural homology models of the human α3 subunit containing AHC mutations were created. RESULTS: The AHC mutations examined all showed similar levels of reduction in forward cycling. Wild type forward cycling was reduced by coexpression with any mutant, indicating dominant negative interactions. Proton transport was measured and found to be selectively impaired only in E815K. Homology modeling showed that D801 and G947 lie within or near known cation binding sites while E815 is more distal. Despite its effect on proton transport, E815K was also distant from the proposed proton transport route. INTERPRETATION: Loss of forward cycling and dominant negativity are common and likely necessary pathomechanisms for AHC. In addition, loss of proton transport correlated with severity of AHC. D801N and G947R are likely to directly disrupt normal Na(+)/K(+) binding while E815K may disrupt forward cycling and proton transport via allosteric mechanisms yet to be elucidated.

Dynamics of fluorescent dyes attached to G-quadruplex DNA and their effect on FRET experiments.

FRET spectroscopy is a promising approach for investigating the dynamics of G-quadruplex DNA folds and improving the targeting of G-quadruplexes by potential anticancer compounds. To better interpret such experiments, classical and replica-exchange molecular dynamics simulations and fluorescence-lifetime measurements are used to understand the behavior of a range of Cy3-based dyes attached to the 3' end of G-quadruplex DNA. The simulations revealed that the dyes interact extensively with the G-quadruplex. Identification of preferred dye positions relative to the G-quadruplex in the simulations allows the impact of dye-DNA interactions on FRET results to be determined. All the dyes show significant deviations from the common approximation of being freely rotating and not interacting with the host, but one of the Cy3 dye analogues is slightly closer to this case.

Locating the route of entry and binding sites of benzocaine and phenytoin in a bacterial voltage gated sodium channel.

Sodium channel blockers are used to control electrical excitability in cells as a treatment for epileptic seizures and cardiac arrhythmia, and to provide short term control of pain. Development of the next generation of drugs that can selectively target one of the nine types of voltage-gated sodium channel expressed in the body requires a much better understanding of how current channel blockers work. Here we make use of the recently determined crystal structure of the bacterial voltage gated sodium channel NavAb in molecular dynamics simulations to elucidate the position at which the sodium channel blocking drugs benzocaine and phenytoin bind to the protein as well as to understand how these drugs find their way into resting channels. We show that both drugs have two likely binding sites in the pore characterised by nonspecific, hydrophobic interactions: one just above the activation gate, and one at the entrance to the the lateral lipid filled fenestrations. Three independent methods find the same sites and all suggest that binding to the activation gate is slightly more favourable than at the fenestration. Both drugs are found to be able to pass through the fenestrations into the lipid with only small energy barriers, suggesting that this can represent the long posited hydrophobic entrance route for neutral drugs. Our simulations highlight the importance of a number of residues in directing drugs into and through the fenestration, and in forming the drug binding sites.

Binding and channeling of alternative substrates in the enzyme DmpFG: a molecular dynamics study.

DmpFG is a bifunctional enzyme comprised of an aldolase subunit, DmpG, and a dehydrogenase subunit, DmpF. The aldehyde intermediate produced by the aldolase is channeled directly through a buried molecular channel in the protein structure from the aldolase to the dehydrogenase active site. In this study, we have investigated the binding of a series of progressively larger substrates to the aldolase, DmpG, using molecular dynamics. All substrates investigated are easily accommodated within the active site, binding with free energy values comparable to the physiological substrate 4-hydroxy-2-ketovalerate. Subsequently, umbrella sampling was utilized to obtain free energy surfaces for the aldehyde intermediates (which would be generated from the aldolase reaction on each of these substrates) to move through the channel to the dehydrogenase DmpF. Small substrates were channeled with limited barriers in an energetically feasible process. We show that the barriers preventing bulky intermediates such as benzaldehyde from moving through the wild-type protein can be removed by selective mutation of channel-lining residues, demonstrating the potential for tailoring this enzyme to allow its use for the synthesis of specific chemical products. Furthermore, positions of transient escape routes in this flexible channel were determined.

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.

An entropic mechanism of generating selective ion binding in macromolecules.

Several mechanisms have been proposed to explain how ion channels and transporters distinguish between similar ions, a process crucial for maintaining proper cell function. Of these, three can be broadly classed as mechanisms involving specific positional constraints on the ion coordinating ligands which arise through: a "rigid cavity", a 'strained cavity' and 'reduced ligand fluctuations'. Each operates in subtly different ways yet can produce markedly different influences on ion selectivity. Here we expand upon preliminary investigations into the reduced ligand fluctuation mechanism of ion selectivity by simulating how a series of model systems respond to a decrease in ligand thermal fluctuations while simultaneously maintaining optimal ion-ligand binding distances. Simple abstract-ligand models, as well as simple models based upon the ion binding sites in two amino acid transporters, show that limiting ligand fluctuations can create ion selectivity between Li(+), Na(+) and K(+) even when there is no strain associated with the molecular framework accommodating the different ions. Reducing the fluctuations in the position of the coordinating ligands contributes to selectivity toward the smaller of two ions as a consequence of entropic differences.