Orientational Correlation Function and Persistence Lengths of Flexible Polyelectrolytes

Exploring Mechanisms of Lipid Nanoparticle-Mucus Interactions in Healthy and Cystic Fibrosis Conditions

Mucus forms the first defense line of human lungs, and as such hampers the efficient delivery of therapeutics to the underlying epithelium. This holds particularly true for genetic cargo such as CRISPR-based gene editing tools which cannot readily surmount the mucosal barrier. While lipid nanoparticles (LNPs) emerge as versatile non-viral gene delivery systems that can help overcome the delivery challenge, many knowledge gaps remain, especially for diseased states such as cystic fibrosis (CF). This study provides fundamental insights into Cas9 mRNA or ribonucleoprotein-loaded LNP-mucus interactions in healthy and diseased states by assessing the impact of the genetic cargo, mucin sialylation, mucin concentration, ionic strength, pH, and polyethylene glycol (PEG) concentration and nature on LNP diffusivity leveraging experimental approaches and Brownian dynamics (BD) simulations. Taken together, this study identifies key mucus and LNP characteristics that are critical to enabling a rational LNP design for transmucosal delivery.

Dilute polyelectrolyte solutions: recent progress and open questions

Polyelectrolytes are a class of polymers possessing ionic groups on their repeating units. Since counterions can dissociate from the polymer backbone, polyelectrolyte chains are strongly influenced by electrostatic interactions. As a result, the physical properties of polyelectrolyte solutions are significantly different from those of electrically neutral polymers. The aim of this article is to highlight key results and some outstanding questions in the polyelectrolyte research from recent literature. We focus on the influence of electrostatics on conformational and hydrodynamic properties of polyelectrolyte chains. A compilation of experimental results from the literature reveals significant disparities with theoretical predictions. We also discuss a new class of polyelectrolytes called poly(ionic liquid)s that exhibit unique physical properties in comparison to ordinary polyelectrolytes. We conclude this review by listing some key research challenges in order to fully understand the conformation and dynamics of polyelectrolytes in solutions.

Length-dependent Intramolecular Coil-to-Globule Transition in Poly(ADP-ribose) Induced by Cations

Poly(ADP-ribose) (PAR), as part of a post-translational modification, serves as a flexible scaffold for noncovalent protein binding. Such binding is influenced by PAR chain length through a mechanism yet to be elucidated. Structural insights have been elusive, partly due to the difficulties associated with synthesizing PAR chains of defined lengths. Here, we employ an integrated approach combining molecular dynamics (MD) simulations with small-angle X-ray scattering (SAXS) experiments, enabling us to identify highly heterogeneous ensembles of PAR conformers at two different, physiologically relevant lengths: PAR 15 and PAR 22 . Our findings reveal that numerous factors including backbone conformation, base stacking, and chain length contribute to determining the structural ensembles. We also observe length-dependent compaction of PAR upon the addition of small amounts of Mg 2+ ions, with the 22-mer exhibiting ADP-ribose bundles formed through local intramolecular coil-to-globule transitions. This study illuminates how such bundling could be instrumental in deciphering the length-dependent action of PAR.

GRAPHICAL ABSTRACT

Interplay of electrolyte concentration and molecular weight of polyDADMAC on cellulose surface adsorption

Cationic polyelectrolytes (PEs) are commonly used additives in manufacturing of cellulose based products such as regenerated fibers and paper to tailor their product properties. Here we are studying the adsorption of poly(diallyldimethylammonium chloride), PD, on cellulose, using in situ surface plasmon resonance spectroscopy (SPR) measurements. We employ model surfaces from regenerated cellulose xanthate (CX) and trimethylsilyl cellulose (TMSC), mimicking industrially relevant regenerated cellulose substrates. The effects of the PDs molecular weight were strongly depending on the ionic strength and type of electrolyte (NaCl vs CaCl₂). Without electrolytes, the adsorption was monolayer-type, i.e. independent of molecular weight. At moderate ionic strength, adsorption increased due to more pronounced PE coiling, while at high ionic strength electrostatic shielding strongly reduced adsorption of PDs. Results exhibited pronounced differences for the chosen substrates (cellulose regenerated from xanthate (CX_reg) vs. regenerated from trimethylsilyl cellulose, TMSC_reg). Consistently higher adsorbed amounts of the PD were determined on CX_reg surfaces compared TMSC. This can be attributed to a more negative zeta potential, a higher AFM roughness and a higher degree of swelling (investigated by QCM-D) of the CX_reg substrates.

Modulation of DNA conformation in electrolytic nanodroplets

The behavior of deoxyribonucleic acid (DNA) molecules in confinement is of profound importance in various bioengineering and medical applications. In the present study, all-atom molecular dynamics simulation is utilized to investigate the transition of the double-strand DNA (dsDNA) conformation in the electrolytic nanodroplet. Three typical conformations, i.e., C-shaped, folded S-shaped, and double C-shaped, are observed for different droplet sizes and ionic concentrations. To reveal the physics underlying this phenomenon, the characteristics of the dsDNA molecules, such as the overcharging intensity, the end-to-end distance, the radius of gyration, etc. are analyzed in detail based on the numerical results. It is found that the transition can be ascribed to the buckling of the polymer molecules under the compression due to the confinement of the nanodroplet, and it can be modulated by the ionic concentration in the electrolyte. Generally, nanoscale confinement dominates dsDNA behavior over the electrostatic effects in smaller nanodroplets, while the latter becomes more important for larger nanodroplets. This competition results in the persistence length increasing with the nanodroplet radii. Based on these discussions, a non-dimensional elasto-capillary number μ is proposed to classify the dsDNA conformations into three regions.

Actin crosslinker competition and sorting drive emergent GUV size-dependent actin network architecture

The proteins that make up the actin cytoskeleton can self-assemble into a variety of structures. In vitro experiments and coarse-grained simulations have shown that the actin crosslinking proteins α-actinin and fascin segregate into distinct domains in single actin bundles with a molecular size-dependent competition-based mechanism. Here, by encapsulating actin, α-actinin, and fascin in giant unilamellar vesicles (GUVs), we show that physical confinement can cause these proteins to form much more complex structures, including rings and asters at GUV peripheries and centers; the prevalence of different structures depends on GUV size. Strikingly, we found that α-actinin and fascin self-sort into separate domains in the aster structures with actin bundles whose apparent stiffness depends on the ratio of the relative concentrations of α-actinin and fascin. The observed boundary-imposed effect on protein sorting may be a general mechanism for creating emergent structures in biopolymer networks with multiple crosslinkers.

Actin crosslinker competition and sorting drive emergent GUV size-dependent actin network architecture

The proteins that make up the actin cytoskeleton can self-assemble into a variety of structures. In vitro experiments and coarse-grained simulations have shown that the actin crosslinking proteins α-actinin and fascin segregate into distinct domains in single actin bundles with a molecular size-dependent competition-based mechanism. Here, by encapsulating actin, α-actinin, and fascin in giant unilamellar vesicles (GUVs), we show that physical confinement can cause these proteins to form much more complex structures, including rings and asters at GUV peripheries and centers; the prevalence of different structures depends on GUV size. Strikingly, we found that α-actinin and fascin self-sort into separate domains in the aster structures with actin bundles whose apparent stiffness depends on the ratio of the relative concentrations of α-actinin and fascin. The observed boundary-imposed effect on protein sorting may be a general mechanism for creating emergent structures in biopolymer networks with multiple crosslinkers.

By encapsulating proteins in giant unilamellar vesicles, Bashirzadeh et al find that actin crosslinkers, α-actinin and fascin, can self-assemble with actin into complex structures that depend on the degree of confinement. Further analysis and modeling show that α-actinin and fascin sort to separate domains of these structures. These insights may be generalizable to other biopolymer networks containing crosslinkers.

Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics

Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.

Ion Transport and the True Transference Number in Nonaqueous Polyelectrolyte Solutions for Lithium Ion Batteries

Nonaqueous polyelectrolyte solutions have been recently proposed as high Li+ transference number electrolytes for lithium ion batteries. However, the atomistic phenomena governing ion diffusion and migration in polyelectrolytes are poorly understood, particularly in nonaqueous solvents. Here, the structural and transport properties of a model polyelectrolyte solution, poly(allyl glycidyl ether-lithium sulfonate) in dimethyl sulfoxide, are studied using all-atom molecular dynamics simulations. We find that the static structural analysis of Li+ ion pairing is insufficient to fully explain the overall conductivity trend, necessitating a dynamic analysis of the diffusion mechanism, in which we observe a shift from largely vehicular transport to more structural diffusion as the Li+ concentration increases. Furthermore, we demonstrate that despite the significantly higher diffusion coefficient of the lithium ion, the negatively charged polyion is responsible for the majority of the solution conductivity at all concentrations, corresponding to Li+ transference numbers much lower than previously estimated experimentally. We quantify the ion-ion correlations unique to polyelectrolyte systems that are responsible for this surprising behavior. These results highlight the need to reconsider the approximations typically made for transport in polyelectrolyte solutions.

Shear-solvo defect annihilation of diblock copolymer thin films over a large area

Achieving defect-free block copolymer (BCP) nanopatterns with a long-ranged orientation over a large area remains a persistent challenge, impeding the successful and widespread application of BCP self-assembly. Here, we demonstrate a new experimental strategy for defect annihilation while conserving structural order and enhancing uniformity of nanopatterns. Sequential shear alignment and solvent vapor annealing generate perfectly aligned nanopatterns with a low defect density over centimeter-scale areas, outperforming previous single or sequential combinations of annealing. The enhanced order quality and pattern uniformity were characterized in unprecedented detail via scattering analysis and incorporating new mathematical indices using elaborate image processing algorithms. In addition, using an advanced sampling method combined with a coarse-grained molecular simulation, we found that domain swelling is the driving force for enhanced defect annihilation. The superior quality of large-scale nanopatterns was further confirmed with diffraction and optical properties after metallized patterns, suggesting strong potential for application in optoelectrical devices.

Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Double-stranded (ds) RNAs play essential roles in many processes of cell metabolism. The knowledge of three-dimensional (3D) structure, stability, and flexibility of dsRNAs in salt solutions is important for understanding their biological functions. In this work, we further developed our previously proposed coarse-grained model to predict 3D structure, stability, and flexibility for dsRNAs in monovalent and divalent ion solutions through involving an implicit structure-based electrostatic potential. The model can make reliable predictions for 3D structures of extensive dsRNAs with/without bulge/internal loops from their sequences, and the involvement of the structure-based electrostatic potential and corresponding ion condition can improve the predictions for 3D structures of dsRNAs in ion solutions. Furthermore, the model can make good predictions for thermal stability for extensive dsRNAs over the wide range of monovalent/divalent ion concentrations, and our analyses show that the thermally unfolding pathway of dsRNA is generally dependent on its length as well as its sequence. In addition, the model was employed to examine the salt-dependent flexibility of a dsRNA helix, and the calculated salt-dependent persistence lengths are in good accordance with experiments.

Correlation anisotropy and stiffness of DNA molecules confined in nanochannels

The anisotropy of orientational correlations in DNA molecules confined in cylindrical channels is explored by Monte Carlo simulations using a coarse-grained model of double-stranded (ds) DNA. We find that the correlation function ⟨C(s)⟩_⊥ in the transverse (confined) dimension exhibits a region of negative values in the whole range of channel sizes. Such a clear-cut sign of the opposite orientation of chain segments represents a microscopic validation of the Odijk deflection mechanism in narrow channels. At moderate-to-weak confinement, the negative ⟨C(s)⟩_⊥ correlations imply a preference of DNA segments for transverse looping. The inclination for looping can explain a reduction of stiffness as well as the enhanced knotting of confined DNA relative to that detected earlier in bulk at some channel sizes. Furthermore, it is shown that the orientational persistence length P_or fails to convey the apparent stiffness of DNA molecules in channels. Instead, correlation lengths P_∥ and P_⊥ in the axial and transverse directions, respectively, encompass the channel-induced modifications of DNA stiffness.

Surface Forces between Highly Charged Cationic Polyelectrolytes Adsorbed to Silica: How Control of pH and the Adsorbed Amount Determines the Net Surface Charge

Atomic force microscopy (AFM) and quartz crystal microbalance with dissipation (QCM-D) were employed to investigate the pH-dependent adsorption of poly(diallyldimethylammonium chloride) (polyDADMAC) to silica surfaces as well as the surface forces between these layers. It was found that polyDADMAC adopted a relatively flat conformation when adsorbed to a silica surface and that the adsorbed amount increased with increasing pH. From the surface force measurements it is evident that the surface undergoes a charge reversal upon saturation with polyDADMAC, at the three different investigated pH values, and that some degree of charge regulation of the silica surface takes place during the adsorption process. Finally, the overcharging phenomenon is discussed in terms of a geometrical mismatch due to the different average spacing between the surface charges on the silica surface and the size of the DADMAC monomer.

Scale-Dependent Stiffness and Internal Tension of a Model Brush Polymer

Bottle-brush polymers exhibit closely grafted side chains that interact by steric repulsion, thereby causing stiffening of the main polymer chain. We use single-molecule elasticity measurements of model brush polymers to quantify this effect. We find that stiffening is only significant on long length scales, with the main chain retaining flexibility on short scales. From the elasticity data, we extract an estimate of the internal tension generated by side-chain repulsion; this estimate is consistent with the predictions of blob-based scaling theories.

Visualizing single-stranded nucleic acids in solution

Single-stranded nucleic acids (ssNAs) are ubiquitous in many key cellular functions. Their flexibility limits both the number of high-resolution structures available, leaving only a small number of protein-ssNA crystal structures, while forcing solution investigations to report ensemble averages. A description of the conformational distributions of ssNAs is essential to more fully characterize biologically relevant interactions. We combine small angle X-ray scattering (SAXS) with ensemble-optimization methods (EOM) to dynamically build and refine sets of ssNA structures. By constructing candidate chains in representative dinucleotide steps and refining the models against SAXS data, a broad array of structures can be obtained to match varying solution conditions and strand sequences. In addition to the distribution of large scale structural parameters, this approach reveals, for the first time, intricate details of the phosphate backbone and underlying strand conformations. Such information on unperturbed strands will critically inform a detailed understanding of an array of problems including protein-ssNA binding, RNA folding and the polymer nature of NAs. In addition, this scheme, which couples EOM selection with an iteratively refining pool to give confidence in the underlying structures, is likely extendable to the study of other flexible systems.

Scaling analysis of random walks with persistence lengths: Application to self-avoiding walks

We develop an approach for performing scaling analysis of N-step random walks (RWs). The mean square end-to-end distance, 〈R[over ⃗]_{N}^{2}〉, is written in terms of inner persistence lengths (IPLs), which we define by the ensemble averages of dot products between the walker's position and displacement vectors, at the jth step. For RW models statistically invariant under orthogonal transformations, we analytically introduce a relation between 〈R[over ⃗]_{N}^{2}〉 and the persistence length, λ_{N}, which is defined as the mean end-to-end vector projection in the first step direction. For self-avoiding walks (SAWs) on 2D and 3D lattices we introduce a series expansion for λ_{N}, and by Monte Carlo simulations we find that λ_{∞} is equal to a constant; the scaling corrections for λ_{N} can be second- and higher-order corrections to scaling for 〈R[over ⃗]_{N}^{2}〉. Building SAWs with typically 100 steps, we estimate the exponents ν_{0} and Δ_{1} from the IPL behavior as function of j. The obtained results are in excellent agreement with those in the literature. This shows that only an ensemble of paths with the same length is sufficient for determining the scaling behavior of 〈R[over ⃗]_{N}^{2}〉, being that the whole information needed is contained in the inner part of the paths.

Salt contribution to the flexibility of single-stranded nucleic acid offinite length

Nucleic acids are negatively charged macromolecules and their structure properties are strongly coupled to metal ions in solutions. In this article, the salt effects on the flexibility of single-stranded (ss) nucleic acid chain ranging from 12 to 120 nucleotides are investigated systematically by the coarse-grained Monte Carlo simulations where the salt ions are considered explicitly and the ss chain is modeled with the virtual-bond structural model. Our calculations show that, the increase of ion concentration causes the structural collapse of ss chain and multivalent ions are much more efficient in causing such collapse, and both trivalent/small divalent ions can induce more compact state than a random relaxation state. We found that monovalent, divalent, and trivalent ions can all overcharge ss chain, and the dominating source for such overcharging changes from ion-exclusion-volume effect to ion Coulomb correlations. In addition, the predicted Na(+) and Mg(2+)-dependent persistence length l(p)'s of ss nucleic acid are in accordance with the available experimental data, and through systematic calculations, we obtained the empirical formulas for l(p) as a function of [Na(+)], [Mg(2+)] and chain length.

Counter-ion distribution around flexible polyelectrolytes having different molecular architecture

We explore the monovalent counter-ion distribution around flexible highly-charged polyelectrolytes with different molecular architectures (linear chains, stars, and unknotted and trefoil rings) using molecular dynamics simulations that include an explicit solvent that interacts with the polyelectrolyte. In particular, we find that the molecular topology influences the fraction of counter-ions transiently associating with the polyelectrolyte on a scale of the order of the chain segments, forming a "condensed" counter-ion interfacial layer. As with the hydrogen bonding of water to proteins and other polymers, the persistence time of these interfacial "bound" counter-ions is relatively short, O(1 ps), and we characterize the fluctuations in the number of the counter-ions populating the interfacial layer. We also find that the counter-ions are distributed in a non-uniform fashion on the polyelectrolyte backbone, forming dynamical clusters whose form and average size is sensitive to molecular architecture. In addition, we find that the residual bound counter-ions, not located in either the interfacial layer or the bulk solution, form a diffuse ionic cloud around the polyelectrolyte due to the uncompensated polyelectrolyte charge along the backbone. Generally charge valence strongly influences the extent of the diffuse counter-ion cloud, but in the case of monovalent counter-ions, we find that the size of the diffuse counter-ion cloud nearly coincides with the polyelectrolyte radius of gyration, independent of molecular topology.

Coarse-grained simulation study of sequence effects on DNA hybridization in a concentrated environment

A novel coarse-grained model is developed to elucidate thermodynamics and kinetic mechanisms of DNA self-assembly. It accounts for sequence and solvent conditions to capture key experimental results such as sequence-dependent thermal property and salt-dependent persistence length of ssDNA and dsDNA. Moreover, constant-temperature simulations on two single strands of a homogeneous sequence show two main mechanisms of hybridization: a slow slithering mechanism and a one-order faster zippering mechanism. Furthermore, large-scale simulations at a high DNA strand concentration demonstrate that DNA self-assembly is a robust and enthalpically driven process in which the formation of double helices is deciphered to occur via multiple self-assembly pathways including the strand displacement mechanism. However, sequence plays an important role in shifting the majority of one pathway over the others and controlling size distribution of self-assembled aggregates. This study yields a complex picture on the role of sequence on programmable self-assembly and demonstrates a promising simulation tool that is suitable for studies in DNA nanotechnology.

Measuring the differential stoichiometry and energetics of ligand binding to macromolecules by single-molecule force spectroscopy: an extended theory

Many chemical techniques exist for measuring the stoichiometry of ligand binding to a macromolecule; however, these techniques are often specific to certain ligands or require the presumption of specific binding models. Here, we further develop a previously reported, general, thermodynamic method for extracting the change in number of ligands bound to a macromolecule as that macromolecule undergoes a conformational transition driven by mechanical stretching, for example, by magnetic tweezers or optical trapping. We extend the theory of this method to consider systems with many ligands, experiments conducted in different thermodynamic ensembles (e.g., constant-force, constant-extension), and experiments in which the system is not at equilibrium. Further, we show that analysis of the same single-molecule mechanical manipulation data yields a measure of the differential free energy of stabilization due to ligand binding, that is, the free energy contribution by which ligand binding favors one conformation of the macromolecule over another. We interpret an existing data set measuring ion binding to RNA and DNA in terms of this free energy.

The snakelike chain character of unstructured RNA

In the absence of base-pairing and tertiary structure, ribonucleic acid (RNA) assumes a random-walk conformation, modulated by the electrostatic self-repulsion of the charged, flexible backbone. This behavior is often modeled as a Kratky-Porod "wormlike chain" (WLC) with a Barrat-Joanny scale-dependent persistence length. In this study we report measurements of the end-to-end extension of poly(U) RNA under 0.1 to 10 pN applied force and observe two distinct elastic-response regimes: a low-force, power-law regime characteristic of a chain of swollen blobs on long length scales and a high-force, salt-valence-dependent regime consistent with ion-stabilized crumpling on short length scales. This short-scale structure is additionally supported by force- and salt-dependent quantification of the RNA ion atmosphere composition, which shows that ions are liberated under stretching; the number of ions liberated increases with increasing bulk salt concentration. Both this result and the observation of two elastic-response regimes directly contradict the WLC model, which predicts a single elastic regime across all forces and, when accounting for scale-dependent persistence length, the opposite trend in ion release with salt concentration. We conclude that RNA is better described as a "snakelike chain," characterized by smooth bending on long length scales and ion-stabilized crumpling on short length scales. In monovalent salt, these two regimes are separated by a characteristic length that scales with the Debye screening length, highlighting the determining importance of electrostatics in RNA conformation.

Branched-linear polyion complexes investigated by Monte Carlo simulations

Complexes formed by one charged and branched copolymer with an oppositely charged and linear polyion have been investigated by Monte Carlo simulations. A coarse-grained description has been used, in which the main chain of the branched polyion and the linear polyion possess the same absolute charge and charge density. The spatial extension and other structural properties, such as bond-angle orientational correlation function, asphericity, and scaling analysis of formed complexes, at varying branching density and side-chain length of the branched polyion, have been explored. In particular, the balance between cohesive Coulomb attraction and side-chain repulsions resulted in two main structures of a polyion complex. These structures are (i) a globular polyion core surrounded by side chains appearing at low branching density and (ii) an extended polyion core with side chains still being expelled at high branching density. The globule-to-extended transition occurred at a crossover branching density being practically independent of the side chain length.

Double-exponential decay of orientational correlations in semiflexible polyelectrolytes

In this paper we revisited the problem of persistence length of polyelectrolytes. We performed a series of Molecular Dynamics simulations using the Debye-Hückel approximation for electrostatics to test several equations which go beyond the classical description of Odijk, Skolnick and Fixman (OSF). The data confirm earlier observations that in the limit of large contour separations the decay of orientational correlations can be described by a single-exponential function and the decay length can be described by the OSF relation. However, at short countour separations the behaviour is more complex. Recent equations which introduce more complicated expressions and an additional length scale could describe the results very well on both the short and the long length scale. The equation of Manghi and Netz when used without adjustable parameters could capture the qualitative trend but deviated in a quantitative comparison. Better quantitative agreement within the estimated error could be obtained using three equations with one adjustable parameter: 1) the equation of Manghi and Netz; 2) the equation proposed by us in this paper; 3) the equation proposed by Cannavacciuolo and Pedersen. Two characteristic length scales can be identified in the data: the intrinsic or bare persistence length and the electrostatic persistence length. All three equations use a single parameter to describe a smooth crossover from the short-range behaviour dominated by the intrinsic stiffness of the chain to the long-range OSF-like behaviour.

Stretching semiflexible polymer chains: evidence for the importance of excluded volume effects from Monte Carlo simulation

Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled by self-avoiding walks on the simple cubic lattice (d = 3 dimensions) and square lattice (d = 2 dimensions), varying chain stiffness by an energy penalty ε(b) for chain bending. In the absence of excluded volume interactions, the persistence length l(p) of the polymers would then simply be l(p) = l(b)(2d - 2)(-1)q(b) (-1) with q(b) = exp(-ε(b)/k(B)T), the bond length l(b) being the lattice spacing, and k(B)T is the thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method (PERM), both q(b) and the chain length N are varied over a wide range (0.005 ≤ q(b) ≤ 1, N ≤ 50,000), and also a stretching force f is applied to one chain end (fixing the other end at the origin). In the absence of this force, in d = 2 a single crossover from rod-like behavior (for contour lengths less than l(p)) to swollen coils occurs, invalidating the Kratky-Porod model, while in d = 3 a double crossover occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that are swollen due to the excluded volume interaction. If the stretching force is applied, excluded volume interactions matter for the force versus extension relation irrespective of chain stiffness in d = 2, while theories based on the Kratky-Porod model are found to work in d = 3 for stiff chains in an intermediate regime of chain extensions. While for q(b) ≪ 1 in this model a persistence length can be estimated from the initial decay of bond-orientational correlations, it is argued that this is not possible for more complex wormlike chains (e.g., bottle-brush polymers). Consequences for the proper interpretation of experiments are briefly discussed.

Adsorption of PEO-PPO-PEO triblock copolymers with end-capped cationic chains of poly(2-dimethylaminoethyl methacrylate)

We study the adsorption of a symmetric triblock copolymer of ethylene oxide, EO, and propylene oxide, PO, end-capped with quarternized poly(2-dimethylaminoethyl methacrylate), DMAEMA (DMAEMA(24)-EO(132)PO(50)EO(132)-DMAEMA(24)). Light scattering and tensiometry are used to measure the relative size of the associated structures and surface excess at the air-liquid interface. The adsorbed amount, the amount of coupled water, and the viscoelasticity of the adsorbed polymer layer are measured on hydrophobic and hydrophilic surfaces (polypropylene, cellulose, and silica) by using quartz crystal microgravimetry (QCM) and surface plasmon resonance (SPR) at different ionic strengths and temperatures. The results of the experiments are compared with those obtained after adsorption of the uncharged precursor copolymer, without the cationic end-caps (EO(132)PO(50)EO(132)). DMAEMA(24)-EO(132)PO(50)EO(132)-DMAEMA(24) possesses higher affinity with the negatively charged silica and cellulose surfaces while the uncharged copolymer adsorbs to a larger extent on polypropylene surfaces. In this latter case, adsorption increases with increasing solution ionic strength and temperature. Adsorption of EO(132)PO(50)EO(132) on silica surfaces has little effect on the water contact angle (WCA), while adsorption of DMAEMA(24)-EO(132)PO(50)EO(132)-DMAEMA(24) increases the WCA of silica to 32°, indicating a large density of exposed PPO blocks upon adsorption. After adsorption of EO(132)PO(50)EO(132) and DMAEMA(24)-EO(132)PO(50)EO(132)-DMAEMA(24) on PP, the WCA is reduced by ≈14° and ≈28°, respectively, due to the exposed hydrophilic EO and highly water-soluble DMAEMA segments on the surfaces. The extent of surface coverage at saturation at the polypropylene/liquid interfaces (≈31 and 40 nm(2)/molecule obtained by QCM and SPR, respectively) is much lower, as expected, when compared with results obtained at the air/liquid interface, where a tighter packing is observed. The percentage of water coupled to the adsorbed cationic polymer decreases with solution ionic strength. Overall, these observations are ascribed to the effects of electrostatic screening, polymer hydrodynamic size, and solvency.

Persistence length of DNA molecules confined in nanochannels

The influence of confinement on the persistence length of dsDNA molecules under a high ionic strength environment was explored by coarse-grained Monte Carlo simulations in channels of different profiles. It was found that under confinement three definitions of the persistence length of DNA molecules were not equivalent and represented different properties. In case of the global quantities, the projection and the WLC persistence lengths, the apparent values up to several hundred nanometres are observed for DNA confined in narrow channels. The orientational correlation function cos theta(s) of confined DNA shows a complex pattern, distinctive for semiflexible polymers. At weak and moderate confinements the function cos theta(s) suggests an unexpected increase in the apparent DNA flexibility. The orientational persistence length computed from the initial slope of the function cos theta(s) mirrors only short-scale correlations and gives the value close to the intrinsic persistence length of DNA. The simulation data of direct relevance to experimental studies of DNA in microfluidic devices are compared with analytical theories for stiff chains.

Conformational characteristics of single flexible polyelectrolyte chain

The behaviour of a flexible anionic chain of 150 univalent and negatively charged beads connected by a harmonic-like potential with each other in the presence of an equal number of positive and free counterions, is studied in molecular dynamics simulations with Langevin thermostat in a wide range of temperatures. Simulations were carried out for several values of the bending parameter, corresponding to fully flexible polyion, moderately and strongly stiff polyion as well as for the case when bend conformation is preferable to the straight one. We have found that in all cases three regimes can be distinguished, which can be characterized as "random coil", observed at high temperatures; "extended conformation" observed at moderate temperatures (of the order of 1 in reduced units), and compact "globular conformation" attained at low temperatures. While the transition between high-temperature random and extended conformations is gradual, the transition from the extended coil to the globular state, taking place at a temperature of about 0.2 in reduced units, is of abrupt character resembling a phase transition.

Swelling of biological and semiflexible polyelectrolytes

We have developed a theoretical model of swelling of semiflexible (biological) polyelectrolytes in salt solutions. Our approach is based on separation of length scales which allowed us to split a chain's electrostatic energy into two parts that describe local and remote electrostatic interactions along the polymer backbone. The local part takes into account interactions between charged monomers that are separated by distances along the polymer backbone shorter than the chain's persistence length. These electrostatic interactions renormalize chain persistence length. The second part includes electrostatic interactions between remote charged pairs along the polymer backbone located at distances larger than the chain persistence length. These interactions are responsible for chain swelling. In the framework of this approach we calculated effective chain persistence length and chain size as a function of the Debye screening length, chain degree of ionization, bare persistence length and chain degree of polymerization. Our crossover expression for the effective chain's persistence length is in good quantitative agreement with the experimental data on DNA. We have been able to fit experimental datasets by using two adjustable parameters: DNA ionization degree (α = 0.15-0.17) and a bare persistence length (l(p) = 40-44 nm).

A mesoscale model of DNA and its renaturation

A mesoscale model of DNA is presented (3SPN.1), extending the scheme previously developed by our group. Each nucleotide is mapped onto three interaction sites. Solvent is accounted for implicitly through a medium-effective dielectric constant and electrostatic interactions are treated at the level of Debye-Hückel theory. The force field includes a weak, solvent-induced attraction, which helps mediate the renaturation of DNA. Model parameterization is accomplished through replica exchange molecular dynamics simulations of short oligonucleotide sequences over a range of composition and chain length. The model describes the melting temperature of DNA as a function of composition as well as ionic strength, and is consistent with heat capacity profiles from experiments. The dependence of persistence length on ionic strength is also captured by the force field. The proposed model is used to examine the renaturation of DNA. It is found that a typical renaturation event occurs through a nucleation step, whereby an interplay between repulsive electrostatic interactions and colloidal-like attractions allows the system to undergo a series of rearrangements before complete molecular reassociation occurs.

Structural transitions of encapsidated polyelectrolytes

Conformations and structural transitions of polyelectrolytes strictly confined onto a spherical 2D surface have been investigated by scaling descriptions based on physical arguments concerning polyelectrolyte adsorption onto planar surface and liquid crystals as well as by Monte Carlo simulations using a bead-spring model with short-range and electrostatic repulsions. In case of the electrostatic screened regime, a disordered-ordered (spiral) transition at increasing persistence length of the chain was found. It was predicted that the transition occurred when the persistence length is comparable with the mean spacing between adjacent strands of the ordered chain. The presence of a non-screened electrostatic repulsion led to a more complex behavior with i) a re-entrant order-disorder transition and ii) a tennis ball texture as an additional smectic/nematic structure. The various competing structures given by the theory were recovered by the Monte Carlo simulations, which also indicated that the tennis ball texture was favored over the spiral structure by the long-range interactions for semi-flexible chains.

Polymer properties of polythymine as revealed by translational diffusion

Biopolymers, such as single-stranded DNA (ssDNA), are often described as semiflexible polymers or wormlike chains. We investigated the length dependence of diffusional properties of homogeneous ssDNA (polythymine) with up to 100 nucleotides using fluorescence correlation spectroscopy. We found that the hydrodynamic radius Rh scales according to a power law, with an exponent between 0.5 and 0.7 depending on ionic strength I. With Rh being proportional to the square root of the persistence length Lp, we found that Lp approximately Im, with m=-0.22+/-0.01 for polythymine with 100 residues. For comparison, we performed molecular dynamics (MD) simulations with a force field that accounts for short-range interactions in vacuum, and determined the characteristic polymer properties end-to-end distance R, radius of gyration S, and persistence length Lp of various labeled and nonlabeled polythymine derivatives. We found excellent agreement for the length dependence of simulated S and experimental Rh measured at 100 mM NaCl, revealing that electrostatic interactions are completely shielded in aqueous solution at such ionic strength. MD simulations further showed that polythymine with >approximately 30 residues can be described as a semiflexible polymer with negligible influence of the fluorescent label; and that static flexibility is limited by geometrical and steric constraints as expressed by an intrinsic persistence length of approximately 1.7 nm. These results provide a benchmark for theories and MD simulations describing the influence of electrostatic interactions on polyelectrolyte properties, and thus help to develop a complete and accurate description of ssDNA.

Diblock polyampholytes grafted onto spherical particles: effect of stiffness, charge density, and grafting density

The structure of spherical brushes formed by symmetric diblock polyampholytes end-grafted onto small spherical particles in aqueous solution is examined within the framework of the so-called primitive model using Monte Carlo simulations. The properties of the two blocks are identical except for the sign of their charges. Three different chain flexibilities corresponding to flexible, semiflexible, and stiff blocks are considered at various polyampholyte linear charge densities and grafting densities. The link between the two blocks is flexible at all conditions, and the grafted segments are laterally mobile. Radial and lateral spatial distribution functions of different types and single-chain properties are analyzed. The brush structure strongly depends on the chain flexibility. With flexible chains, a disordered polyelectrolyte complex is formed at the surface of the particle, the complex becoming more compact at increasing linear charge density. With stiff blocks, the inner blocks are radially oriented. At low linear charged density, the outer blocks are orientationally disordered, whereas at increasing electrostatic interaction the two blocks of a polyampholyte are parallel and close to each other, leading to an ordered structure referred to as a polyampholyte star. As the grafting density is increased, the brush thickness responds differently for flexible and nonflexible chains, depending on a different balance between electrostatic interactions and excluded volume effects.