Disentangling steric and electrostatic factors in nanoscale transport through confined space

Controlling DNA Fragments Translocation across Nanopores with the Synergic Use of Site-Directed Mutagenesis, pH-Dependent Charge Tuning, and Electroosmotic Flow

Biological and solid-state nanopores are at the core of transformative techniques and nanodevices, democratizing the examination of matter and biochemical reactions at the single-molecule level, with low cost, portability, and simplicity in operation. One of the crucial hurdles in such endeavors is the fast analyte translocation, which limits characterization, and a rich number of strategies have been explored over the years to overcome this. Here, by site-directed mutagenesis on the α-hemolysin protein nanopore (α-HL), sought to replace selected amino acids with glycine, electrostatic binding sites were induced on the nanopore's vestibule and constriction region and achieved in the most favorable case a 20-fold increase in the translocation time of short single-stranded DNA (ssDNA) at neutral pH, with respect to the wild-type (WT) nanopore. We demonstrated an efficient tool of controlling the ssDNA translocation time, via the interplay between the nanopore-ssDNA surface electrostatic interactions and electroosmotic flow, all mediated by the pH-dependent ionization of amino acids lining the nanopore's translocation pathway. Our data also reveal the nonmonotonic, pH-induced alteration of ssDNA average translocation time. Unlike mildly acidic conditions (pH ∼ 4.7), at a pH ∼ 2.8 maintained symmetrically or asymmetrically across the WT α-HL, we evidenced the manifestation of a dominant electroosmotic flow, determining the speeding up of the ssDNA translocation across the nanopore by counteracting the ssDNA-nanopore attractive electrostatic interactions. We envision potential applications of the presented approach by enabling easy-to-use, real-time detection of short ssDNA sequences, without the need for complex biochemical modifications to the nanopore to mitigate the fast translocation of such sequences.

How molecular interactions tune the characteristic time of nanocomposite colloidal sensors

HYPOTHESIS

Mass transport critically controls the performance of colloidal metal-polymer sensors. We hypothesize that molecular-level pair interactions, such as electric, steric, and specific binding effects, govern the mass transport and, in return, the characteristic time of these sensors.

THEORY

Here we present a simple theory guided by experimental data to examine the sensing performance of two usually encountered archetypal metal-polymer sensors, namely (1) core-shell and (2) yolk-shell architectures. For this purpose, we use the static reactive density functional theory framework, determining how (i) charge, (ii) size, and (iii) non-covalent binding factors modulate the characteristic time.

FINDINGS

We show how an interplay between diffusivity and partitioning governs the sensing time of the sensors, where an anti-correlation cancellation between them renders the time non-trivial. Our study demonstrates that the convoluted substrate-hydrogel shell interaction controls the characteristic time of these colloidal sensors, especially when the sensors are in a collapsed state. Notably, the substrates with a high dipole moment tend to equilibrate greatly, but undesirably, at the shell-solution interface. With this, we encourage the formation of a metastable sorption state.

Laser-based temperature control to study the roles of entropy and enthalpy in polymer-nanopore interactions

Single-molecule approaches for probing the free energy of confinement for polymers in a nanopore environment are critical for the development of nanopore biosensors. We developed a laser-based nanopore heating approach to monitor the free energy profiles of such a single-molecule sensor. Using this approach, we measure the free energy profiles of two distinct polymers, polyethylene glycol and water-soluble peptides, as they interact with the nanopore sensor. Polyethylene glycol demonstrates a retention mechanism dominated by entropy with little sign of interaction with the pore, while peptides show an enthalpic mechanism, which can be attributed to physisorption to the nanopore (e.g., hydrogen bonding). To manipulate the energetics, we introduced thiolate-capped gold clusters [Au25(SG)18] into the pore, which increases the charge and leads to additional electrostatic interactions that help dissect the contribution that enthalpy and entropy make in this modified environment. These observations provide a benchmark for optimization of single-molecule nanopore sensors.

Accurate modeling of a biological nanopore with an extended continuum framework

Despite the broad success of biological nanopores as powerful instruments for the analysis of proteins and nucleic acids at the single-molecule level, a fast simulation methodology to accurately model their nanofluidic properties is currently unavailable. This limits the rational engineering of nanopore traits and makes the unambiguous interpretation of experimental results challenging. Here, we present a continuum approach that can faithfully reproduce the experimentally measured ionic conductance of the biological nanopore Cytolysin A (ClyA) over a wide range of ionic strengths and bias potentials. Our model consists of the extended Poisson-Nernst-Planck and Navier-Stokes (ePNP-NS) equations and a computationally efficient 2D-axisymmetric representation for the geometry and charge distribution of the nanopore. Importantly, the ePNP-NS equations achieve this accuracy by self-consistently considering the finite size of the ions and the influence of both the ionic strength and the nanoscopic scale of the pore on the local properties of the electrolyte. These comprise the mobility and diffusivity of the ions, and the density, viscosity and relative permittivity of the solvent. Crucially, by applying our methodology to ClyA, a biological nanopore used for single-molecule enzymology studies, we could directly quantify several nanofluidic characteristics difficult to determine experimentally. These include the ion selectivity, the ion concentration distributions, the electrostatic potential landscape, the magnitude of the electro-osmotic flow field, and the internal pressure distribution. Hence, this work provides a means to obtain fundamental new insights into the nanofluidic properties of biological nanopores and paves the way towards their rational engineering.

Ensemble-Based Steered Molecular Dynamics Predicts Relative Residence Time of A2A Receptor Binders

Drug-target residence time, the length of time for which a small molecule stays bound to its receptor target, has increasingly become a key property for optimization in drug discovery programs. However, its in silico prediction has proven difficult. Here we describe a method, using atomistic ensemble-based steered molecular dynamics (SMD), to observe the dissociation of ligands from their target G protein-coupled receptor in a time scale suitable for drug discovery. These dissociation simulations accurately, precisely, and reproducibly identify ligand-residue interactions and quantify the change in ligand energy values for both protein and water. The method has been applied to 17 ligands of the A_2A adenosine receptor, all with published experimental kinetic binding data. The residues that interact with the ligand as it dissociates are known experimentally to have an effect on binding affinities and residence times. There is a good correlation ( R² = 0.79) between the computationally calculated change in water-ligand interaction energy and experimentally determined residence time. Our results indicate that ensemble-based SMD is a rapid, novel, and accurate semi-empirical method for the determination of drug-target relative residence time.

Synthetically Diversified Protein Nanopores: Resolving Click Reaction Mechanisms

Nanopores are emerging as a powerful tool for the investigation of nanoscale processes at the single-molecule level. Here, we demonstrate the methionine-selective synthetic diversification of α-hemolysin (α-HL) protein nanopores and their exploitation as a platform for investigating reaction mechanisms. A wide range of functionalities, including azides, alkynes, nucleotides, and single-stranded DNA, were incorporated into individual pores in a divergent fashion. The ion currents flowing through the modified pores were used to observe the trajectory of a range of azide-alkyne click reactions and revealed several short-lived intermediates in Cu(I)-catalyzed azide-alkyne [3 + 2] cycloadditions (CuAAC) at the single-molecule level. Analysis of ion-current fluctuations enabled the populations of species involved in rapidly exchanging equilibria to be determined, facilitating the resolution of several transient intermediates in the CuAAC reaction mechanism. The versatile pore-modification chemistry offers a useful approach for enabling future physical organic investigations of reaction mechanisms at the single-molecule level.

Information Dynamics of a Nonlinear Stochastic Nanopore System

Nanopores have become a subject of interest in the scientific community due to their potential uses in nanometer-scale laboratory and research applications, including infectious disease diagnostics and DNA sequencing. Additionally, they display behavioral similarity to molecular and cellular scale physiological processes. Recent advances in information theory have made it possible to probe the information dynamics of nonlinear stochastic dynamical systems, such as autonomously fluctuating nanopore systems, which has enhanced our understanding of the physical systems they model. We present the results of local (LER) and specific entropy rate (SER) computations from a simulation study of an autonomously fluctuating nanopore system. We learn that both metrics show increases that correspond to fluctuations in the nanopore current, indicating fundamental changes in information generation surrounding these fluctuations.

Discrimination of supramolecular chirality using a protein nanopore

Supramolecular chirality may emerge from self-assembly processes to yield architectures that differ only in the topological arrangement of their constituent parts. Since the properties of the resulting enantiomeric assemblies are identical, purification and characterisation can be challenging. Here, we have examined the hypothesis that the intrinsic chirality of a protein nanopore can be exploited to detect supramolecular chirality. Transient blockages in the ion current flowing through a single membrane-spanning α-haemolysin nanopore were shown to discriminate between M4L6 tetrahedral coordination cages of opposing chiralities. The single-molecule nature of the approach facilitated direct access to the rates of association and dissociation with the nanopore, which allowed the concentrations of the enantiomeric supramolecular assemblies to be determined in situ. Thus, we have established that a protein nanopore can be used to discriminate the chiral topologies of supramolecular assemblies, even when they are too large to fully enter the nanopore.

Electrolyte Gradient-Based Modulation of Molecular Transport through Nanoporous Gold Membranes

Nanopores, and nanoporous materials in general, are interesting for applications in chemical and biomolecular transport as pore sizes are on the same scale as the dimension of many (bio)chemical species. Many studies have focused on either single pores or small arrays of cylindrical pores, which are convenient in terms of their amenability toward computational modeling of transport phenomenon. However, the limited overall porosity may inhibit transport flux as well as the eventual implementation of these materials as active separation elements. Inspired by its relatively high porosity, we have explored nanoporous gold (NPG) as a membrane across which small molecular species can be transported. NPG offers a random, bicontinuous pore geometry, while also being inherently conductive and readily amenable to surface modification-attributes that may be enabling in the pursuit of size- and charge-based approaches to molecular separations. NPG was fabricated via a free-corrosion process whereby immersion of Au-containing alloys in concentrated nitric acid preferentially dissolves the less noble metals (e.g., Ni, Cu). Average pore diameters of 50 ± 20 nm were obtained as verified under scanning electron microscopy. NPG membranes were sandwiched between two reservoirs, and the selective transport of chemical species across the membrane in the presence of an ionic strength gradient was investigated. The flux of small molecules were monitored by UV-vis absorption spectrometry and found to be dependent upon the direction and magnitude of the ionic strength gradient. Moreover, transport trends underscored the effects of surface charge in a confined environment, considering that the pore diameters were on the same scale as the electrical double layer experienced by molecules transiting the membrane. Under such conditions, the transport of anions and cations through NPG was found to depend on an induced electric field as well as ion advection. Further electrical and surface chemical modulations of transport are expected to engender increased membrane functionality.

Temperature Effect on Ionic Current and ssDNA Transport through Nanopores

We have investigated the role of electrostatic interactions in the transport of nucleic acids and ions through nanopores. The passage of DNA through nanopores has so far been conjectured to involve a free-energy barrier for entry, followed by a downhill translocation where the driving voltage accelerates the polymer. We have tested the validity of this conjecture by using two toxins, α-hemolysin and aerolysin, which differ in their shape, size, and charge. The characteristic timescales in each toxin as a function of temperature show that the entry barrier is ∼15 kBT and the translocation barrier is ∼35 kBT, although the electrical force in the latter step is much stronger. Resolution of this fact, using a theoretical model, reveals that the attraction between DNA and the charges inside the barrel of the pore is the most dominant factor in determining the translocation speed and not merely the driving electrochemical potential gradient.

DNA modulates solvent isotope effects in a nanopore

Here we investigate the modulation of solvent isotope effects by the entry of DNA molecules into individual α-haemolysin nanopores. Solvent isotope effects in D2O versus H2O were enhanced (kH/kD ≈ 1.6) compared to the bulk (kH/kD ≈ 1.2), except when the pore was most blocked (kH/kD ≤ 1.1).

Comparative analysis of nucleotide translocation through protein nanopores using steered molecular dynamics and an adaptive biasing force

The translocation of nucleotide molecules across biological and synthetic nanopores has attracted attention as a next generation technique for sequencing DNA. Computer simulations have the ability to provide atomistic-level insight into important states and processes, delivering a means to develop a fundamental understanding of the translocation event, for example, by extracting the free energy of the process. Even with current supercomputing facilities, the simulation of many-atom systems in fine detail is limited to shorter timescales than the real events they attempt to recreate. This imposes the need for enhanced simulation techniques that expand the scope of investigation in a given timeframe. There are numerous free energy calculation and translocation methodologies available, and it is by no means clear which method is best applied to a particular problem. This article explores the use of two popular free energy calculation methodologies in a nucleotide-nanopore translocation system, using the α-hemolysin nanopore. The first uses constant velocity-steered molecular dynamics (cv-SMD) in conjunction with Jarzynski's equality. The second applies an adaptive biasing force (ABF), which has not previously been applied to the nucleotide-nanpore system. The purpose of this study is to provide a comprehensive comparison of these methodologies, allowing for a detailed comparative assessment of the scientific merits, the computational cost, and the statistical quality of the data obtained from each technique. We find that the ABF method produces results that are closer to experimental measurements than those from cv-SMD, whereas the net errors are smaller for the same computational cost.