Polyelectrolyte Decomplexation via Addition of Salt: Charge Correlation Driven Zipper

Ion-Specific Effects on Ion and Polyelectrolyte Solvation

Ion-specific effects on aqueous solvation of monovalent counter ions, Na+ ${^+ }$ , K+ ${^+ }$ , Cl- ${^- }$ , and Br- ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl- ${^- }$ and Br- ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

Nonmonotonic electrophoretic mobility of rodlike polyelectrolytes by multivalent coions in added salt

It is well established that when multivalent counterions or salts are added to a solution of highly charged polyelectrolytes (PEs), correlation effects can cause charge inversion of the PE, leading to electrophoretic mobility (EM) reversal. In this work, we use coarse-grained molecular-dynamics simulations to unravel the less understood effect of coion valency on EM reversal for rigid DNA-like PEs. We find that EM reversal induced by multivalent counterions is suppressed with increasing coion valency in the salt added and eventually vanishes. Further, we find that EM is enhanced at fixed low salt concentrations for salts with monovalent counterions when multivalent coions with increasing valency are introduced. However, increasing the salt concentration causes a crossover that leads to EM reversal which is enhanced by increasing coion valency at high salt concentration. Remarkably, this multivalent coion-induced EM reversal persists even for low values of PE linear charge densities where multivalent counterions alone cannot induce EM reversal. These results facilitate tuning PE-PE interactions and self-assembly with both coion and counterion valencies.

A Perspective on the Glass Transition and the Dynamics of Polyelectrolyte Multilayers and Complexes

Polyelectrolyte multilayers (PEMs) or polyelectrolyte complexes (PECs), formed by layer-by-layer assembly or the mixing of oppositely charged polyelectrolytes (PEs) in aqueous solution, respectively, have potential applications in health, energy, and the environment. PEMs and PECs are very tunable because their structure and properties are influenced by factors such as pH, ionic strength, salt type, humidity, and temperature. Therefore, it is increasingly important to understand how these factors affect PECs and PEMs on a molecular level. In this Feature Article, we summarize our contributions to the field in the development of approaches to quantify the swelling, thermal properties, and dynamic mechanical properties of PEMs and PECs. First, the role of water as a plasticizer and in the glass-transition temperature (Tg) in both strong poly(diallyldimethylammonium)/poly(sodium 4-styrenesulfonate) (PDADMA/PSS) and weak poly(allylamine hydrochloride)/poly(acrylic acid) (PAH/PAA) systems is presented. Then, factors influencing the dynamics of PECs and PEMs are discussed. We also reflect on the swelling of PEMs in response to different salts and solvent additives. Last, the nature of water's microenvironment in PEMs/PECs is discussed. A special emphasis is placed on experimental techniques, along with molecular simulations. Taken together, this review presents an outlook and offers recommendations for future research directions, such as studying the additional effects of hydrogen-bonding hydrophobic interactions.

Theoretical and computational analysis of the electrophoretic polymer mobility inversion induced by charge correlations

Electrophoretic (EP) mobility reversal is commonly observed for strongly charged macromolecules in multivalent salt solutions. This curious effect takes place, e.g., when a charged polymer, such as DNA, adsorbs excess counterions so that the counterion-dressed surface charge reverses its sign, leading to the inversion of the polymer drift driven by an external electric field. In order to characterize this seemingly counterintuitive phenomenon that cannot be captured by electrostatic mean-field theories, we adapt here a previously developed strong-coupling-dressed Poisson-Boltzmann approach to the cylindrical geometry of the polyelectrolyte-salt system. Within the framework of this formalism, we derive an analytical polymer mobility formula dressed by charge correlations. In qualitative agreement with polymer transport experiments, this mobility formula predicts that the increment of the monovalent salt, the decrease of the multivalent counterion valency, and the increase of the dielectric permittivity of the background solvent suppress charge correlations and increase the multivalent bulk counterion concentration required for EP mobility reversal. These results are corroborated by coarse-grained molecular dynamics simulations showing how multivalent counterions induce mobility inversion at dilute concentrations and suppress the inversion effect at large concentrations. This re-entrant behavior, previously observed in the aggregation of like-charged polymer solutions, calls for verification by polymer transport experiments.

Interaction between two polyelectrolytes in monovalent aqueous salt solutions

We use the recently developed soft-potential-enhanced Poisson-Boltzmann (SPB) theory to study the interaction between two parallel polyelectrolytes (PEs) in monovalent ionic solutions in the weak-coupling regime. The SPB theory is fitted to ion distributions from coarse-grained molecular dynamics (MD) simulations and benchmarked against all-atom MD modelling for poly(diallyldimethylammonium) (PDADMA). We show that the SPB theory is able to accurately capture the interactions between two PEs at distances beyond the PE radius. For PDADMA positional correlations between the charged groups lead to locally asymmetric PE charge and ion distributions. This gives rise to small deviations from the SPB prediction that appear as short-range oscillations in the potential of mean force. Our results suggest that the SPB theory can be an efficient way to model interactions in chemically specific complex PE systems.

Modified Poisson–Boltzmann theory for polyelectrolytes in monovalent salt solutions with finite-size ions

We present a soft-potential-enhanced Poisson–Boltzmann (SPB) theory to efficiently capture ion distributions and electrostatic potential around rodlike charged macromolecules. The SPB model is calibrated with a coarse-grained particle-based model for polyelectrolytes (PEs) in monovalent salt solutions as well as compared to a full atomistic molecular dynamics simulation with the explicit solvent. We demonstrate that our modification enables the SPB theory to accurately predict monovalent ion distributions around a rodlike PE in a wide range of ion and charge distribution conditions in the weak-coupling regime. These include excess salt concentrations up to 1M and ion sizes ranging from small ions, such as Na+ or Cl−, to softer and larger ions with a size comparable to the PE diameter. The work provides a simple way to implement an enhancement that effectively captures the influence of ion size and species into the PB theory in the context of PEs in aqueous salt solutions.

Influence of temperature and salt on coacervation in an aqueous mixture of poly-L-lysine (PLL) and poly-(sodium styrene sulfonate) (PSSNa)

The mixture of poly-L-lysine (PLL) and long-chain PSSNa can lead to the formation of soluble complexes depending on pH, PLL concentration, ionic strength, and temperature. The influence of these stimuli was studied by zetametry, dynamic and ultra-small-angle light scattering, and turbidimetric measurements. First of all, we studied the stoichiometry of complexation, and then considered the influence of salt concentration and temperature on the behavior of the mixture at different pH values. These findings have allowed us to conclude that the polyelectrolyte-polypeptide stoichiometry is controlled by electrostatic interactions between opposite charges. At mass ratios between 1.8 and 2.3 and with net charges close to neutrality, unstable complexes were formed and flocculated due to the hydrophobic attraction leading to macroscopic phase separation. The linear charge density of the complex is also controlled by the ionic strength. Higher CaCl₂ concentrations reduce the complex stability and decrease the charge density, which leads to surface patch binding (SPB) at higher pH. Finally, the electrostatic interactions and strength of hydrogen bonds increased the stabilization of the complexes formed at temperatures lower than 45 °C. At temperatures higher than 45 °C, hydrophobic interactions became more dominant, causing a destabilization of the complexes.

pH-Induced Changes in Polypeptide Conformation: Force-Field Comparison with Experimental Validation

Microsecond-long all-atom molecular dynamics (MD) simulations, circular dichroism, laser Doppler velocimetry, and dynamic light-scattering techniques have been used to investigate pH-induced changes in the secondary structure, charge, and conformation of poly l-lysine (PLL) and poly l-glutamic acid (PGA). The employed combination of the experimental methods reveals for both PLL and PGA a narrow pH range at which they are charged enough to form stable colloidal suspensions, maintaining their α-helix content above 60%; an elevated charge state of the peptides required for colloidal stability promotes the peptide solvation as a random coil. To obtain a more microscopic view on the conformations and to verify the modeling performance, peptide secondary structure and conformations rising in MD simulations are also examined using three different force fields, i.e., OPLS-AA, CHARMM27, and AMBER99SB*-ILDNP. Ramachandran plots reveal that in the examined setup the α-helix content is systematically overestimated in CHARMM27, while OPLS-AA overestimates the β-sheet fraction at lower ionization degrees. At high ionization degrees, the OPLS-AA force-field-predicted secondary structure fractions match the experimentally measured distribution most closely. However, the pH-induced changes in PLL and PGA secondary structure are reasonably captured only by the AMBER99SB*-ILDNP force field, with the exception of the fully charged PGA in which the α-helix content is overestimated. The comparison to simulations results shows that the examined force fields involve significant deviations in their predictions for charged homopolypeptides. The detailed mapping of secondary structure dependency on pH for the polypeptides, especially finding the stable colloidal α-helical regime for both examined peptides, has significant potential for practical applications of the charged homopolypeptides. The findings raise attention especially to the pH fine tuning as an underappreciated control factor in surface modification and self-assembly.

Fourier transform infrared spectroscopy investigation of water microenvironments in polyelectrolyte multilayers at varying temperatures

Polyelectrolyte multilayers (PEMs) are thin films formed by the alternating deposition of oppositely charged polyelectrolytes. Water plays an important role in influencing the physical properties of PEMs, as it can act both as a plasticizer and swelling agent. However, the way in which water molecules distribute around and hydrate ion pairs has not been fully quantified with respect to both temperature and ionic strength. Here, we examine the effects of temperature and ionic strength on the hydration microenvironments of fully immersed poly(diallyldimethylammonium)/polystyrene sulfonate (PDADMA/PSS) PEMs. This is accomplished by tracking the OD stretch peak using attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy at 0.25-1.5 M NaCl and 35-70 °C. The OD stretch peak is deconvoluted into three peaks: (1) high frequency water, which represents a tightly bound microenvironment, (2) low frequency water, which represents a loosely bound microenvironment, and (3) bulk water. In general, the majority of water absorbed into the PEM exists in a bound state, with little-to-no bulk water observed. Increasing temperature slightly reduces the amount of absorbed water, while addition of salt increases the amount of absorbed water. Finally, a van't Hoff analysis is applied to estimate the enthalpy (11-22 kJ mol-1) and entropy (48-79 kJ mol-1 K-1) of water exchanging from low to high frequency states.

Comparing water-mediated hydrogen-bonding in different polyelectrolyte complexes

All-atom molecular dynamics simulations are used to investigate the polyelectrolyte-specific influence of hydration and temperature on water diffusion in hydrated polyelectrolyte complexes (PECs). Two model PECs were compared: poly(allylamine hydrochloride) (PAH)-poly(sodium 4-styrenesulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA)-poly(acrylic acid) (PAA). The findings show that the strength of the hydrogen bonding i.e. polyelectrolyte water interaction has enormous influence on the water mobility, which has implications for PEC structure and properties. A 10-fold difference in the average water diffusion coefficient between PAH-PSS and PDADMA-PAA PECs at the same hydration level is observed. The vast majority of the water molecules hydrating the PDADMA-PAA PECs, for hydrations in the range of 26-38 wt%, are effectively immobilized, whereas for PAH-PSS PECs the amount of immobilized water decreases with hydration. This points to the polyelectrolyte-specific character of the PE-water hydrogen bonding relationship with temperature. PAA-water hydrogen bonds are found to be significantly less sensitive to temperature than for PSS-water. The polyelectrolyte-water interactions, investigated via radial distribution function, hydrogen bond distance and angle distributions, are connected with resulting structure of the PECs. The PDADMA-PAA and PAH-PSS PECs are prepared experimentally and the states of water at different hydration levels is determined using differential scanning calorimetry (DSC). Experiments confirm the differences between PDADMA-PAA and PAH-PSS PECs observed in the theoretical modelling. The results suggest that the initial predictions of the PEC's bonding with water can be based on simple molecular-level considerations.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.

Solution-Controlled Conformational Switching of an Anchored Wireframe DNA Nanostructure

Self-assembled DNA origami nanostructures have a high degree of programmable spatial control that enables nanoscale molecular manipulations. A surface-tethered, flexible DNA nanomesh is reported herein which spontaneously undergoes sharp, dynamic conformational transitions under physiological conditions. The transitions occur between two major macrostates: a spread state dominated by the interaction between the DNA nanomesh and the BSA/streptavidin surface and a surface-avoiding contracted state. Due to a slow rate of stochastic transition events on the order of tens of minutes, the dynamic conformations of individual structures can be detected in situ with DNA PAINT microscopy. Time series localization data with automated imaging processing to track the dynamically changing radial distribution of structural markers are combined. Conformational distributions of tethered structures in buffers with elevated pH exhibit a calcium-dependent domination of the spread state. This is likely due to electrostatic interactions between the structures and immobilized surface proteins (BSA and streptavidin). An interaction is observed in solution under similar buffer conditions with dynamic light scattering. Exchanging between solutions that promote one or the other state leads to in situ sample-wide transitions between the states. The technique herein can be a useful tool for dynamic control and observation of nanoscale interactions and spatial relationships.

Hydration and Temperature Response of Water Mobility in Poly(diallyldimethylammonium)-Poly(sodium 4-styrenesulfonate) Complexes

The combination of all-atom molecular dynamics simulations with differential scanning calorimetry (DSC) has been exploited to investigate the influence of temperature and hydration on the water distribution and mobility in poly(diallyldimethylammonium) (PDADMA) and poly(sodium 4-styrenesulfonate) (PSS) complexes. The findings show that the vast majority of the water molecules hydrating the polyelectrolyte complexes (PECs) with 18-30 wt % hydration are effectively immobilized due to the strong interactions between the PE charge groups and water. Temperature and hydration were found to decrease similarly the fraction of strongly bound water. Additionally, at low hydration or at low temperatures, water motions become dominantly local vibrations and rotations instead of translational motion; translation dominance is recovered in a similar fashion by increase of both temperature and hydration. DSC experiments corroborate the simulation findings by showing that nonfreezing, bound water dominates in hydrated PECs at comparable hydrations. Our results raise attention to water as an equal variable to temperature in the design and engineering of stimuli-responsive polyelectrolyte materials and provide mechanistic explanation for the similarity.

Ability of the Poisson-Boltzmann equation to capture molecular dynamics predicted ion distribution around polyelectrolytes

Here, we examine polyelectrolyte (PE) and ion chemistry specificity in ion condensation via all-atom molecular dynamics (MD) simulations and assess the ability of the Poisson-Boltzmann (PB) equation to describe the ion distribution predicted by the MD simulations. The PB model enables the extraction of parameters characterizing ion condensation. We find that the modified PB equation which contains the effective PE radius and the energy of the ion-specific interaction as empirical fitting parameters describes ion distribution accurately at large distances but close to the PE, especially when strongly localized charge or specific ion binding sites are present, the mean field description of PB fails. However, the PB model captures the MD predicted ion condensation in terms of the Manning radius and fraction of condensed counterions for all the examined PEs and ion species. We show that the condensed ion layer thickness in our MD simulations collapses on a single master curve for all the examined simple, monovalent ions (Na⁺, Br⁺, Cs⁺, Cl^-, and Br^-) and PEs when plotted against the Manning parameter (and consequently the PE line charge density). The significance of this finding is that, contrary to the Manning radius extracted from the mean field PB model, the condensed layer thickness in the all atom detail MD modelling does not depend on the PE chemistry or counterion type. Furthermore, the fraction of condensed counterions in the MD simulations exceeds the PB theory prediction. The findings contribute toward understanding and modelling ion distribution around PEs and other charged macromolecules in aqueous solutions, such as DNA, functionalized nanotubes, and viruses.

Repulsion between oppositely charged rod-shaped macromolecules: Role of overcharging and ionic confinement

The interaction between two oppositely charged rod-shaped macro-ions in a micro-ion solution is investigated via Monte Carlo simulations of the primitive model. The focus is on the asymmetry in rod and/or ion charge, i.e., conditions where oppositely charged objects can repel one another. For equally and oppositely charged rods with asymmetric z:1 micro-ions, repulsion may be induced by overcharging one of the rods with the z valent ions. For asymmetrically charged rods in a symmetric z:z micro-ion solution, a repulsive interaction-at separation of the order of one ion diameter-can arise via an unbalanced osmotic pressure contribution from the ionic atmosphere in the inter-rod space, and an attractive interaction-at a smaller separation-may occur due to a "squeezing out" of the micro-ions from the space between the rods (with a consequent gain in entropy). The thermodynamics of each mechanism is investigated in terms of rod charge and size and micro-ion valence, size, and concentration. Our findings contribute to the understanding of the complex role of charge asymmetry on the interaction of, for example, oppositely charged polyelectrolytes, functionalized nanotubes, and rod-like biomolecules, e.g., viruses.

Potential of mean force and transient states in polyelectrolyte pair complexation

The pair association between two polyelectrolytes (PEs) of the same size but opposite charge is systematically studied in terms of the potential of mean force (PMF) along their center-of-mass reaction coordinate via coarse-grained, implicit-solvent, explicit-salt computer simulations. The focus is set on the onset and the intermediate transient stages of complexation. At conditions above the counterion-condensation threshold, the PE association process exhibits a distinct sliding-rod-like behavior where the polymer chains approach each other by first stretching out at a critical distance close to their contour length, then "shaking hand" and sliding along each other in a parallel fashion, before eventually folding into a neutral complex. The essential part of the PMF for highly charged PEs can be very well described by a simple theory based on sliding charged "Debye-Hückel" rods with renormalized charges in addition to an explicit entropy contribution owing to the release of condensed counterions. Interestingly, at the onset of complex formation, the mean force between the PE chains is found to be discontinuous, reflecting a bimodal structural behavior that arises from the coexistence of interconnected-rod and isolated-coil states. These two microstates of the PE complex are balanced by subtle counterion release effects and separated by a free-energy barrier due to unfavorable stretching entropy.

Development of the modern theory of polymeric complex coacervation

Oppositely charged polymers can undergo the process of complex coacervation, which refers to a liquid-liquid phase separation driven by electrostatic attraction. These materials have demonstrated considerable promise as the basis for complex, self-assembled materials. In this review, we provide a broad overview of the theoretical tools used to understand the physical properties of polymeric coacervates. In particular, we discuss historic theories (Voorn-Overbeek, Random Phase Approximation), and then describe recent developments in the field (Field Theoretic, Counterion Release, Molecular Simulation, and Polymer Reference Interaction Site Model methods). We provide context for these methods, and map out the patchwork of theoretical models that are used to describe a diverse array of coacervate systems. We use this review of the literature to clarify a number of important theoretical challenges remaining in our physical understanding of complex coacervation.

Atomic-Scale Molecular Dynamics Simulations of DNA-Polycation Complexes: Two Distinct Binding Patterns

Synthetic cationic polymers represent a promising class of delivery vectors for gene therapy. Here, we employ atomistic molecular dynamics simulations to gain insight into the structure and properties of complexes of DNA with four linear polycations: polyethylenimine (PEI), poly-l-lysine (PLL), polyvinylamine (PVA), and polyallylamine (PAA). These polycations differ in their polymer geometries, protonation states, and hydrophobicities of their backbone chains. Overall, our results demonstrate for the first time the existence of two distinct patterns of binding of DNA with polycations. For PEI, PLL, and PAA, the complex is stabilized by the electrostatic attraction between protonated amine groups of the polycation and phosphate groups of DNA. In contrast, PVA demonstrates an alternative binding pattern as it gets embedded into the DNA major groove. It is likely that both the polymer topology and affinity of the backbone chain of PVA to the DNA groove are responsible for such behavior. The differences in binding patterns can have important biomedical implications: embedding PVA into a DNA groove makes it less sensitive to changes in the aqueous environment (pH level, ionic strength, etc.) and could therefore hinder the intracellular release of genetic material from a delivery vector, leading to lower transfection activity.

First Steps in the Aggregation Process of Copolymers Based on Thymine Monomers: Characterization by Molecular Dynamics Simulations and Atomic Force Microscopy

Atomistic molecular dynamic simulations were performed to study the structure of isolated VBT-VBA (vinylbenzylthymine-vinylbenzyltriethylammonium chloride) copolymer chains in water at different monomeric species ratios (1:1 and 1:4). The geometric parameters of the structure that the copolymers form in equilibrium together with the basic interactions that stabilize them were determined. Atomic force microscopy (AFM) measurements of dried diluted concentrations of the two copolymers onto highly oriented pyrolytic graphite (HOPG) substrates were carried out to study their aggregation arrangement. The experiments show that both copolymers arrange in fiber-like structures. Comparing the diameters predicted by the simulation results and those obtained by AFM, it can be concluded that individual copolymers arrange in bunches of two chains, stabilized by contra-ions-copolymer interactions for the 1:1 copolymerization ratio at the ionic strength of our samples. In contrast, for the 1:4 system the individual copolymer chains do not aggregate in bunches. These results remark the relevance of the copolymerization ratio and ionic strength of the solvent in the mesoscopic structure of these materials.

The influence of ionic strength and mixing ratio on the colloidal stability of PDAC/PSS polyelectrolyte complexes

Polyelectrolyte complexes (PECs) form by mixing polycation and polyanion solutions together, and have been explored for a variety of applications. One challenge for PEC processing and application is that under certain conditions the as-formed PECs aggregate and precipitate out of suspension over the course of minutes to days. This aggregation is governed by several factors such as electrostatic repulsion, van der Waals attractions, and hydrophobic interactions. In this work, we explore the boundary between colloidally stable and unstable complexes as it is influenced by polycation/polyanion mixing ratio and ionic strength. The polymers examined are poly(diallyldimethylammonium chloride) (PDAC) and poly(sodium 4-styrenesulfonate) (PSS). Physical properties such as turbidity, hydrodynamic size, and zeta potential are investigated upon complex formation. We also perform detailed molecular dynamics simulations to examine the structure and effective charge distribution of the PECs at varying mixing ratios and salt concentrations to support the experimental findings. The results suggest that the colloidally stable/unstable boundary possibly marks the screening effects from added salt, resulting in weakly charged complexes that aggregate. At higher salt concentrations, the complexes initially form and then gradually dissolve into solution.

Chemistry specificity of DNA-polycation complex salt response: a simulation study of DNA, polylysine and polyethyleneimine

In this work, the chemistry specific stability determining factors of DNA-polycation complexes are examined by performing all-atom molecular dynamics simulations. To this end, we conduct a systematic variation of polycation line charge through polyethyleneimine (PEI) protonation and polycation chemistry via comparison with poly-l-lysine (PLL). Our simulations show that increasing line charge of the polycation alone does not lead to more salt tolerant complexes. Instead, the effective charge compensation by the polycation correlates with the increased stability of the complex against additional salt. The salt stability of PEI-DNA complexes also links to the proton sponge property of weak polycations, commonly assumed to be behind the effectivity of PEI as a gene delivery vector. Examination of the complexes reveals the mechanism behind this behaviour; more Cl(-) ions are attracted by the protonated complexes but, in contrast to the common depiction of the proton sponge behaviour, the ion influx does not cause swelling of the complex structure itself. However, PEI protonation leads to release of PEI while DNA remains tightly bound to the complex. Jointly, these findings shed light on the stability determining factors of DNA-polycation complexes, raise charge distribution as an important stability determining contributor, and indicate that the effectivity of PEI in gene delivery is likely to result from the freed PEI facilitating gene transfection.

Controlling polymer translocation and ion transport via charge correlations

We develop a correlation-corrected transport theory in order to predict ionic and polymer transport properties of membrane nanopores under physical conditions where mean-field electrostatics breaks down. The experimentally observed low KCl conductivity of open α-hemolysin pores is quantitatively explained by the presence of surface polarization effects. Upon the penetration of a DNA molecule into the pore, these polarization forces combined with the electroneutrality of DNA sets a lower boundary for the ionic current, explaining the weak salt dependence of blocked pore conductivities at dilute ion concentrations. The addition of multivalent counterions to the solution results in the reversal of the polymer charge and the direction of the electroosmotic flow. With trivalent spermidine or quadrivalent spermine molecules, the charge inversion is strong enough to stop the translocation of the polymer and to reverse its motion. This mechanism can be used efficiently in translocation experiments in order to improve the accuracy of DNA sequencing by minimizing the translocation velocity of the polymer.