Modeling Transport Through Synthetic Nanopores

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

DNA double helix, a tiny electromotor

Flowing fluid past chiral objects has been used for centuries to power rotary motion in man-made machines. By contrast, rotary motion in nanoscale biological or chemical systems is produced by biasing Brownian motion through cyclic chemical reactions. Here we show that a chiral biological molecule, a DNA or RNA duplex rotates unidirectionally at billions of revolutions per minute when an electric field is applied along the duplex, with the rotation direction being determined by the chirality of the duplex. The rotation is found to be powered by the drag force of the electro-osmotic flow, realizing the operating principle of a macroscopic turbine at the nanoscale. The resulting torques are sufficient to power rotation of nanoscale beads and rods, offering an engineering principle for constructing nanoscale systems powered by electric field.

Tunable Nanopore Arrays as the Basis for Ionic Circuits

There has been considerable interest in preparing ionic circuits capable of manipulating ionic and molecular transport in a solution. This direction of research is inspired by biological systems where multiple pores with different functionalities embedded in a cell membrane transmit external signals and underlie all physiological processes. In this manuscript, we describe the modeling of ion transport through small arrays of nanopores consisting of 3, 6, and 9 nanopores and an integrated gate electrode placed on the membrane surface next to one pore opening. We show that by tuning the gate voltage and strategically placing nanopores with nonlinear current-voltage characteristics, the local signal at the gate affects ionic transport through all nanopores in the array. Conditions were identified when the same gate voltage induced opposite rectification properties of neighboring nanopores. We also demonstrate that an ionic diode embedded in a nanopore array can modulate transport properties of neighboring pores even without a gate voltage. The results are explained by the role of concentration polarization and overlapping depletion zones on one side of the membrane. The modeling presented here is intended to become an inspiration to future experiments to create nanopore arrays that can transduce signals in space and time.

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

Resistive amplitude fingerprints during translocation of linear molecules through charged solid-state nanopores

We report the first analytical theory on the amplitude of resistive signals during molecular translocation through charged solid-state nanopores with variable cross-sectional area and piecewise-constant surface charge densities. By providing closed-form explicit algebraic expressions for the concentration profiles inside charged nanopores, this theory allows the prediction of baseline and translocation resistive signals without the need for numerical simulation of the electrokinetic phenomena. A transversely homogenized theory and an asymptotic expansion for weakly charged pores capture DC or quasi-static rectification due to field-induced intrapore concentration polarization (as a result of pore charge inhomogeneity or a translocating molecule). This theory, validated by simulations and experiments, is then used to explain why the amplitude of a single stranded DNA molecule can be twice as high as the amplitude of its double stranded counterpart. It also suggests designs for intrapore concentration polarization and volume exclusion effects that can produce biphasic and other amplitude fingerprints for high-throughput and yet discriminating molecular identification.

Boron nitride nanotubes as containers for targeted drug delivery of doxorubicin

Using molecular dynamics simulations, the adsorption and diffusion of doxorubicin drug molecules in boron nitride nanotubes are investigated. The interaction between doxorubicin and the nanotube is governed by van der Waals attraction. We find strong adsorption of doxorubicin to the wall for narrow nanotubes (radius of 9 Å). For larger radii (12 and 15 Å), the adsorption energy decreases, while the diffusion coefficient of doxorubicin increases. It does, however, not reach the values of pure water, as adsorption events still hinder the doxorubicin mobility. It is concluded that nanotubes wider than around 4 nm diameter can serve as efficient drug containers for targeted drug delivery of doxorubicin in cancer chemotherapy.

Effect of Octadecylamine Surfactant on DNA Interactions with Graphene Surfaces

Understanding of how to integrate DNA molecules with graphene materials is important for the development of biosensors and biomolecular logic circuits. For some of these applications, controlling DNA structural conformation on the graphene substrate is critically important and can be achieved through the use of self-assembled monolayers. Here, we performed all-atom molecular dynamics simulations to understand how various 1-octadecylamine (ODA) coatings of the graphene surface affect the conformation of double-stranded DNA (dsDNA) on the surface. The simulation results demonstrated that dsDNA structures become more stable as ODA concentration increases due to the formation of DNA-ODA hydrogen bonds and reduction of DNA-surface interactions, which aid in retaining internal DNA interactions. Specifically, the interaction of ODA molecules with DNA prevents nucleobases from forming π-π stacking interactions with the surface. Some dsDNA conformations, such as sharp kinks or unwinding, can occur more frequently in DNA with A-T sequences due to weaker pairing interactions than with G-C sequences. Furthermore, our results conclude that both DNA sequence and ODA concentration play an essential role in experimentally observed conformational changes of DNA on the graphene surface.

Adsorption and Diffusion of Cisplatin Molecules in Nanoporous Materials: A Molecular Dynamics Study

Using molecular dynamics simulations, the adsorption and diffusion of cisplatin drug molecules in nanopores is investigated for several inorganic materials. Three different materials are studied with widely-varying properties: metallic gold, covalent silicon, and silica. We found a strong influence of both the van der Waals and the electrostatic interaction on the adsorption behavior on the pore walls, which in turn influence the diffusion coefficients. While van der Waals forces generally lead to a reduction of the diffusion coefficient, the fluctuations in the electrostatic energy induced by orientation changes of the cisplatin molecule were found to help desorb the molecule from the wall.

Multiscale analysis of the effect of surface charge pattern on a nanopore's rectification and selectivity properties: From all-atom model to Poisson-Nernst-Planck

We report a multiscale modeling study for charged cylindrical nanopores using three modeling levels that include (1) an all-atom explicit-water model studied with molecular dynamics, and reduced models with implicit water containing (2) hard-sphere ions studied with the Local Equilibrium Monte Carlo simulation method (computing ionic correlations accurately), and (3) point ions studied with Poisson-Nernst-Planck theory (mean-field approximation). We show that reduced models are able to reproduce device functions (rectification and selectivity) for a wide variety of charge patterns, that is, reduced models are useful in understanding the mesoscale physics of the device (i.e., how the current is produced). We also analyze the relationship of the reduced implicit-water models with the explicit-water model and show that diffusion coefficients in the reduced models can be used as adjustable parameters with which the results of the explicit- and implicit-water models can be related. We find that the values of the diffusion coefficients are sensitive to the net charge of the pore but are relatively transferable to different voltages and charge patterns with the same total charge.

Rapid and Accurate Determination of Nanopore Ionic Current Using a Steric Exclusion Model

Nanopore sensing has emerged as a versatile approach to detection and identification of biomolecules. Presently, researchers rely on experience and intuition for choosing or modifying the nanopores to detect a target analyte. The field would greatly benefit from a computational method that could relate the atomic-scale geometry of the nanopores and analytes to the blockade nanopore currents they produce. Existing computational methods are either computationally too expensive to be used routinely in experimental laboratories or not sensitive enough to account for the atomic structure of the pore and the analytes. Here, we demonstrate a robust and inexpensive computational approach-the steric exclusion model (SEM) of nanopore conductance-that is orders of magnitude more efficient than all-atom MD and yet is sensitive enough to account for the atomic structure of the nanopore and the analyte. The method combines the computational efficiency of a finite element solver with the atomic precision of a nanopore conductance map to yield unprecedented speed and accuracy of ionic current prediction. We validate our SEM approach through comparison with the current blockades computed using the all-atom molecular dynamics method for a range of proteins confined to a solid-state nanopore, biological channels embedded in a lipid bilayer membranes, and blockade currents produced by DNA homopolymers in MspA. We illustrate potential applications of SEM by computing blockade currents produced by nucleosome proteins in a solid-state nanopore, individual amino acids in MspA, and by testing the effect of point mutations on amino acid distinguishability. We expect our SEM approach to become an integral part of future development of the nanopore sensing field.

Water-Compression Gating of Nanopore Transport

Electric field-driven motion of biomolecules is a process essential to many analytics methods, in particular, to nanopore sensing, where a transient reduction of nanopore ionic current indicates the passage of a biomolecule through the nanopore. However, before any molecule can be examined by a nanopore, the molecule must first enter the nanopore from the solution. Previously, the rate of capture by a nanopore was found to increase with the strength of the applied electric field. Here, we theoretically show that, in the case of narrow pores in graphene membranes, increasing the strength of the electric field can not only decrease the rate of capture, but also repel biomolecules from the nanopore. As the strong electric field polarizes water near and within the nanopore, the high gradient of the field also produces a strong dielectrophoretic force that compresses the water. The pressure difference caused by the sharp water density gradient produces a hydrostatic force that repels DNA or proteins from the nanopore, preventing, in certain conditions, their capture. We show that such local compression of fluid can regulate the transport of biomolecules through nanoscale passages in the absence of physical gates and sort proteins according to their phosphorylated states.

Multiscale modeling of a rectifying bipolar nanopore: explicit-water versus implicit-water simulations

In a multiscale modeling approach, we present computer simulation results for a rectifying bipolar nanopore at two modeling levels. In an all-atom model, we use explicit water to simulate ion transport directly with the molecular dynamics technique. In a reduced model, we use implicit water and apply the Local Equilibrium Monte Carlo method together with the Nernst-Planck transport equation. This hybrid method makes the fast calculation of ion transport possible at the price of lost details. We show that the implicit-water model is an appropriate representation of the explicit-water model when we look at the system at the device (i.e., input vs. output) level. The two models produce qualitatively similar behavior of the electrical current for different voltages and model parameters. Looking at the details of concentration and potential profiles, we find profound differences between the two models. These differences, however, do not influence the basic behavior of the model as a device because they do not influence the z-dependence of the concentration profiles which are the main determinants of current. These results then address an old paradox: how do reduced models, whose assumptions should break down in a nanoscale device, predict experimental data? Our simulations show that reduced models can still capture the overall device physics correctly, even though they get some important aspects of the molecular-scale physics quite wrong; reduced models work because they include the physics that is necessary from the point of view of device function. Therefore, reduced models can suffice for general device understanding and device design, but more detailed models might be needed for molecular level understanding.

Modulation of Molecular Flux Using a Graphene Nanopore Capacitor

Modulation of ionic current flowing through nanoscale pores is one of the fundamental biological processes. Inspired by nature, nanopores in synthetic solid-state membranes are being developed to enable rapid analysis of biological macromolecules and to serve as elements of nanofludic circuits. Here, we theoretically investigate ion and water transport through a graphene-insulator-graphene membrane containing a single, electrolyte-filled nanopore. By means of all-atom molecular dynamics simulations, we show that the charge state of such a graphene nanopore capacitor can regulate both the selectivity and the magnitude of the nanopore ionic current. At a fixed transmembrane bias, the ionic current can be switched from being carried by an equal mixture of cations and anions to being carried almost exclusively by either cationic or anionic species, depending on the sign of the charge assigned to both plates of the capacitor. Assigning the plates of the capacitor opposite sign charges can either increase the nanopore current or reduce it substantially, depending on the polarity of the bias driving the transmembrane current. Facilitated by the changes of the nanopore surface charge, such ionic current modulations are found to occur despite the physical dimensions of the nanopore being an order of magnitude larger than the screening length of the electrolyte. The ionic current rectification is accompanied by a pronounced electro-osmotic effect that can transport neutral molecules such as proteins and drugs across the solid-state membrane and thereby serve as an interface between electronic and chemical signals.

Graphene Nanopores for Protein Sequencing

An inexpensive, reliable method for protein sequencing is essential to unraveling the biological mechanisms governing cellular behavior and disease. Current protein sequencing methods suffer from limitations associated with the size of proteins that can be sequenced, the time, and the cost of the sequencing procedures. Here, we report the results of all-atom molecular dynamics simulations that investigated the feasibility of using graphene nanopores for protein sequencing. We focus our study on the biologically significant phenylalanine-glycine repeat peptides (FG-nups)-parts of the nuclear pore transport machinery. Surprisingly, we found FG-nups to behave similarly to single stranded DNA: the peptides adhere to graphene and exhibit step-wise translocation when subject to a transmembrane bias or a hydrostatic pressure gradient. Reducing the peptide's charge density or increasing the peptide's hydrophobicity was found to decrease the translocation speed. Yet, unidirectional and stepwise translocation driven by a transmembrane bias was observed even when the ratio of charged to hydrophobic amino acids was as low as 1:8. The nanopore transport of the peptides was found to produce stepwise modulations of the nanopore ionic current correlated with the type of amino acids present in the nanopore, suggesting that protein sequencing by measuring ionic current blockades may be possible.

What controls open-pore and residual currents in the first sensing zone of alpha-hemolysin nanopore? Combined experimental and theoretical study

The electrophoretic transport of single-stranded DNA through biological nanopores such as alpha-hemolysin (αHL) is a promising and cost-effective technology with the potential to revolutionize genomics. The rational design of pores with the controlled polymer translocation rates and high contrast between different nucleotides could improve significantly nanopore sequencing applications. Here, we apply a combination of theoretical and experimental methods in an attempt to elucidate several selective modifications in the pore which were proposed to be central for the effective discrimination between purines and pyrimidines. Our nanopore test set includes the wild type αHL and six mutants (E111N/M113X/K147N) in which the cross-section and chemical functionality of the first constriction zone of the pore are modified. Electrophysiological recordings were combined with all-atom Molecular Dynamics simulations (MD) and a recently developed Brownian Dynamics (BROMOC) protocol to investigate residual ion currents and pore-DNA interactions for two homo-polymers e.g. poly(dA)40 or poly(dC)40 blocking the pore. The calculated residual currents and contrast in the poly(dA)40/poly(dC)40 blocked pore are in qualitative agreement with the experimental recordings. We showed that a simple structural metric allows rationalization of key elements in the emergent contrast between purines and pyrimidines in the modified αHL mutants. The shape of the pore and its capacity for hydrogen bonding to a translocated polynucleotide are two essential parameters for contrast optimization. To further probe the impact of these two factors in the ssDNA sensing, we eliminated the effect of the primary constriction using serine substitutions (i.e. E111S/M113S/T145S/K147S) and increased the hydrophobic volume of the central residue in the secondary constriction (L135I). This pore modification sharply increased the contrast between Adenine (A) and Cytosine (C).

Numerical investigation of molecular nano-array in potential-energy profile for a single dsDNA

A Rigorous numerical investigation on dsDNA translocation in quasi-2-dimensional nano-array filter is performed using Molecular Dynamics (MD) method. Various dsDNA molecules with different sizes are chosen in order to model Ogston sieving in a nano-array filter. The radius of gyration of dsDNA molecule is less than the characteristic length of the shallow region in nano-array. The dsDNA molecule is assumed to be in the 0.05M NaCl electrolyte. MD shows acceptable results for potential-energy profile for nano-array filter. According to the MD outcomes, the dsDNA electrophoretic mobility decreases almost linearly with dsDNA size and show the same trend as Ogston sieving for gel electrophoresis. In addition, different shapes for nano-array filter are studied for a unique dsDNA molecule. It is concluded that steeping the nano-array wall can cause the retardation of dsDNA translocation and decreases dsDNA electrophoretic mobility.

Slowing DNA Transport Using Graphene-DNA Interactions

Slowing down DNA translocation speed in a nanopore is essential to ensuring reliable resolution of individual bases. Thin membrane materials enhance spatial resolution but simultaneously reduce the temporal resolution as the molecules translocate far too quickly. In this study, the effect of exposed graphene layers on the transport dynamics of both single (ssDNA) and double-stranded DNA (dsDNA) through nanopores is examined. Nanopore devices with various combinations of graphene and Al2O3 dielectric layers in stacked membrane structures are fabricated. Slow translocations of ssDNA in nanopores drilled in membranes with layers of graphene are reported. The increased hydrophobic interactions between the ssDNA and the graphene layers could explain this phenomenon. Further confirmation of the hydrophobic origins of these interactions is obtained through reporting significantly faster translocations of dsDNA through these graphene layered membranes. Molecular dynamics simulations confirm the preferential interactions of DNA with the graphene layers as compared to the dielectric layer verifying the experimental findings. Based on our findings, we propose that the integration of multiple stacked graphene layers could slow down DNA enough to enable the identification of nucleobases.

Microsecond simulations of DNA and ion transport in nanopores with novel ion-ion and ion-nucleotides effective potentials

We developed a novel scheme based on the grand-canonical Monte Carlo/Brownian dynamics simulations and have extended it to studies of ion currents across three nanopores with the potential for single-stranded DNA (ssDNA) sequencing: solid-state nanopore Si₃N₄, α-hemolysin, and E111N/M113Y/K147N mutant. To describe nucleotide-specific ion dynamics compatible with ssDNA coarse-grained model, we used the inverse Monte Carlo protocol, which maps the relevant ion-nucleotide distribution functions from all-atom molecular dynamics (MD) simulations. Combined with the previously developed simulation platform for Brownian dynamics simulations of ion transport, it allows for microsecond- and millisecond-long simulations of ssDNA dynamics in the nanopore with a conductance computation accuracy that equals or exceeds that of all-atom MD simulations. In spite of the simplifications, the protocol produces results that agree with the results of previous studies on ion conductance across open channels and provide direct correlations with experimentally measured blockade currents and ion conductances that have been estimated from all-atom MD simulations.

Conformational transitions and stop-and-go nanopore transport of single-stranded DNA on charged graphene

Control over interactions with biomolecules holds the key to applications of graphene in biotechnology. One such application is nanopore sequencing, where a DNA molecule is electrophoretically driven through a graphene nanopore. Here we investigate how interactions of single-stranded DNA and a graphene membrane can be controlled by electrically biasing the membrane. The results of our molecular dynamics simulations suggest that electric charge on graphene can force a DNA homopolymer to adopt a range of strikingly different conformations. The conformational response is sensitive to even very subtle nucleotide modifications, such as DNA methylation. The speed of DNA motion through a graphene nanopore is strongly affected by the graphene charge: a positive charge accelerates the motion, whereas a negative charge arrests it. As a possible application of the effect, we demonstrate stop-and-go transport of DNA controlled by the charge of graphene. Such on-demand transport of DNA is essential for realizing nanopore sequencing.

Structure and response to flow of the glycocalyx layer

The glycocalyx is a sugar-rich layer located at the luminal part of the endothelial cells. It is involved in key metabolic processes and its malfunction is related to several diseases. To understand the function of the glycocalyx, a molecular level characterization is necessary. In this article, we present large-scale molecular-dynamics simulations that provide a comprehensive description of the structure and dynamics of the glycocalyx. We introduce the most detailed, to-date, all-atom glycocalyx model, composed of lipid bilayer, proteoglycan dimers, and heparan sulfate chains with realistic sequences. Our results reveal the folding of proteoglycan ectodomain and the extended conformation of heparan sulfate chains. Furthermore, we study the glycocalyx response under shear flow and its role as a flypaper for binding fibroblast growth factors (FGFs), which are involved in diverse functions related to cellular differentiation, including angiogenesis, morphogenesis, and wound healing. The simulations show that the glycocalyx increases the effective concentration of FGFs, leading to FGF oligomerization, and acts as a lever to transfer mechanical stimulus into the cytoplasmic side of endothelial cells.

Slowing down single-molecule trafficking through a protein nanopore reveals intermediates for peptide translocation

The microscopic details of how peptides translocate one at a time through nanopores are crucial determinants for transport through membrane pores and important in developing nano-technologies. To date, the translocation process has been too fast relative to the resolution of the single molecule techniques that sought to detect its milestones. Using pH-tuned single-molecule electrophysiology and molecular dynamics simulations, we demonstrate how peptide passage through the α-hemolysin protein can be sufficiently slowed down to observe intermediate single-peptide sub-states associated to distinct structural milestones along the pore, and how to control residence time, direction and the sequence of spatio-temporal state-to-state dynamics of a single peptide. Molecular dynamics simulations of peptide translocation reveal the time- dependent ordering of intermediate structures of the translocating peptide inside the pore at atomic resolution. Calculations of the expected current ratios of the different pore-blocking microstates and their time sequencing are in accord with the recorded current traces.

Ion transport through a graphene nanopore

Molecular dynamics simulation is utilized to investigate the ionic transport of NaCl in solution through a graphene nanopore under an applied electric field. Results show the formation of concentration polarization layers in the vicinity of the graphene sheet. The nonuniformity of the ion distribution gives rise to an electric pressure which drives vortical motions in the fluid if the electric field is sufficiently strong to overcome the influence of viscosity and thermal fluctuations. The relative importance of hydrodynamic transport and thermal fluctuations in determining the pore conductivity is investigated. A second important effect that is observed is the mass transport of water through the nanopore, with an average velocity proportional to the applied voltage and independent of the pore diameter. The flux arises as a consequence of the asymmetry in the ion distribution which can be attributed to differing mobilities of the sodium and chlorine ions and to the polarity of water molecules. The accumulation of liquid molecules in the vicinity of the nanopore due to re-orientation of the water dipoles by the local electric field is seen to result in a local increase in the liquid density. Results confirm that the electric conductance is proportional to the nanopore diameter for the parameter regimes that we simulated. The occurrence of fluid vortices is found to result in an increase in the effective electrical conductance.

Optimization of the molecular dynamics method for simulations of DNA and ion transport through biological nanopores

Molecular dynamics (MD) simulations have become a standard method for the rational design and interpretation of experimental studies of DNA translocation through nanopores. The MD method, however, offers a multitude of algorithms, parameters, and other protocol choices that can affect the accuracy of the resulting data as well as computational efficiency. In this chapter, we examine the most popular choices offered by the MD method, seeking an optimal set of parameters that enable the most computationally efficient and accurate simulations of DNA and ion transport through biological nanopores. In particular, we examine the influence of short-range cutoff, integration timestep and force field parameters on the temperature and concentration dependence of bulk ion conductivity, ion pairing, ion solvation energy, DNA structure, DNA-ion interactions, and the ionic current through a nanopore.

Modeling Pressure-Driven Transport of Proteins through a Nanochannel

Reducing the size of a nanofluidic channel not only creates new opportunities for high-precision manipulation of biological macromolecules, but also makes the performance of the entire nanofluidic system more susceptible to undesirable interactions between the transported biomolecules and the walls of the channel. In this manuscript, we report molecular dynamics simulations of a pressure-driven flow through a silica nanochannel that characterized, with atomic resolution, adsorption of a model protein to its surface. Although the simulated adsorption of the proteins was found to be nonspecific, it had a dramatic effect on the rate of the protein transport. To determine the relative strength of the protein-silica interactions in different adsorbed states, we simulated flow-induced desorption of the proteins from the silica surface. Our analysis of the protein conformations in the adsorbed states did not reveal any simple dependence of the adsorption strength on the size and composition of the protein-silica contact, suggesting that the heterogeneity of the silica surface may be a important factor.

Modeling nanopores for sequencing DNA

Using nanopores to sequence DNA rapidly and at a low cost has the potential to radically transform the field of genomic research. However, despite all the exciting developments in the field, sequencing DNA using a nanopore has yet to be demonstrated. Among the many problems that hinder development of the nanopore sequencing methods is the inability of current experimental techniques to visualize DNA conformations in a nanopore and directly relate the microscopic state of the system to the measured signal. We have recently shown that such tasks could be accomplished through computation. This chapter provides step-by-step instructions of how to build atomic scale models of biological and solid-state nanopore systems, use the molecular dynamics method to simulate the electric field-driven transport of ions and DNA through the nanopores, and analyze the results of such computational experiments.

Computational microscopy of the role of protonable surface residues in nanoprecipitation oscillations

A novel phenomenon has recently been reported in polymeric nanopores. This phenomenon, so-called nanoprecipitation, is characterized by the transient formation of precipitates in the nanopore lumen, producing a sequence of low and high conductance states in the ionic current through the pore. By means of all-atom molecular dynamics simulations, we studied nanoprecipitation for polyethylene terephthalate nanopore immersed in electrolytic solution containing calcium phosphate, covering a total simulation time of 1.24 micros. Our results suggest that protonable surface residues at the nanopore surface, namely carboxyl groups, trigger the formation of precipitates that strongly adhere to the surface, blocking the pore and producing the low conductance state. On the basis of the simulations, we propose a mechanism for the formation of the high conductance state; the mechanism involves detachment of the precipitate from the surface due to reprotonation of carboxyl groups and subsequent translocation of the precipitate out of the pore.

Deciphering ionic current signatures of DNA transport through a nanopore

Within just a decade from the pioneering work demonstrating the utility of nanopores for molecular sensing, nanopores have emerged as versatile systems for single-molecule manipulation and analysis. In a typical setup, a gradient of the electrostatic potential captures charged solutes from the solution and forces them to move through a single nanopore, across an otherwise impermeable membrane. The ionic current blockades resulting from the presence of a solute in a nanopore can reveal the type of the solute, for example, the nucleotide makeup of a DNA strand. Despite great success, the microscopic mechanisms underlying the functionality of such stochastic sensors remain largely unknown, as it is not currently possible to characterize the microscopic conformations of single biomolecules directly in a nanopore and thereby unequivocally establish the causal relationship between the observables and the microscopic events. Such a relationship can be determined using molecular dynamics-a computational method that can accurately predict the time evolution of a molecular system starting from a given microscopic state. This article describes recent applications of this method to the process of DNA transport through biological and synthetic nanopores.

Functionalization of Artificial Freestanding Composite Nanomembranes

Artificial nanomembranes may be defined as synthetic freestanding structures with a thickness below 100 nm and a very large aspect ratio, of at least a few orders of magnitude. Being quasi-2D, they exhibit a host of unusual properties useful for various applications in energy harvesting, sensing, optics, plasmonics, biomedicine, etc. We review the main approaches to nanomembrane functionalization through nanocompositing, which ensures performance far superior to that of simple nanomembranes. These approaches include lamination (stacking of nanometer-thin strata of different materials), introduction of nanoparticle fillers into the nanomembrane scaffold, nanomembrane surface sculpting and modification through patterning (including formation of nanohole arrays and introduction of ion channels similar in function to those in biological nanomembranes). We also present some of our original results related to functionalization of metal matrix composite nanomembranes.