Equilibrium polyelectrolyte bundles with different multivalent counterion concentrations

DNA-protamine condensates under low salt conditions: molecular dynamics simulation with a simple coarse-grained model focusing on electrostatic interactions

Protamine, a small, strongly positively-charged protein, plays a key role in achieving chromatin condensation inside sperm cells and is also involved in the formulation of nanoparticles for gene therapy and packaging of mRNA-based vaccines against viral infection and cancer. The detailed mechanisms of such condensations are still poorly understood especially under low salt conditions where electrostatic interaction predominates. Our previous study, with a refined coarse-grained model in full consideration of the long-range electrostatic interactions, has demonstrated the crucial role of electrostatic interaction in protamine-controlled reversible DNA condensation. Therefore, we herein pay our attention only to the electrostatic interaction and devise a coarser-grained bead-spring model representing the right linear charge density on protamine and DNA chains but treating other short-range interactions as simply as possible, which would be suitable for real-scale simulations. Effective pair potential calculations and large-scale molecular dynamics simulations using this extremely simple model reproduce the phase behaviour of DNA in a wide range of protamine concentrations under low salt conditions, again revealing the importance of the electrostatic interaction in this process and providing a detailed nanoscale picture of bundle formation mediated by a charge disproportionation mechanism. Our simulations also show that protamine length alters DNA overcharging and in turn redissolution thresholds of DNA condensates, revealing the important role played by entropies and correlated fluctuations of condensing agents and thus offering an additional opportunity to design tailored nanoparticles for gene therapy. The control mechanism of DNA-protamine condensates will also provide a better microscopic picture of biomolecular condensates, i.e., membraneless organelles arising from liquid-liquid phase separation, that are emerging as key principles of intracellular organization. Such condensates controlled by post-translational modification of protamine, in particular phosphorylation, or by variations in protamine length from species to species may also be responsible for the chromatin-nucleoplasm patterning observed during spermatogenesis in several vertebrate and invertebrate species.

Manipulation of Confined Polyelectrolyte Conformations through Dielectric Mismatch

We demonstrate that a highly charged polyelectrolyte confined in a spherical cavity undergoes reversible transformations between amorphous conformations and a four-fold symmetry morphology as a function of dielectric mismatch between the media inside and outside the cavity. Surface polarization due to dielectric mismatch exhibits an extra "confinement" effect, which is most pronounced within a certain range of the cavity radius and the electrostatic strength between the monomers and counterions and multivalent counterions. For cavities with a charged surface, surface polarization leads to an increased amount of counterions adsorbed in the outer side, further compressing the confined polyelectrolyte into a four-fold symmetry morphology. The equilibrium conformation of the chain is dependent upon several key factors including the relative permittivities of the media inside and outside the cavity, multivalent counterion concentration, cavity radius relative to the chain length, and interface charge density. Our findings offer insights into the effects of dielectric mismatch in packaging and delivery of polyelectrolytes across media with different relative permittivities. Moreover, the reversible transformation of the polyelectrolyte conformations in response to environmental permittivity allows for potential applications in biosensing and medical monitoring.

Study of the conformation of polyelectrolyte aggregates using coarse-grained molecular dynamics simulations

The conformation of polyelectrolyte aggregates as a function of the backbone rigidity is investigated by coarse-grained molecular dynamics simulation. The polyelectrolyte is represented by a bead-spring chain with charged side chains. The simulations start from the uniform distributions of the polyelectrolytes, and the resultant polyelectrolyte conformation after a few microseconds exhibits spherical self-aggregates, clusters, or bending bundle-like aggregates, depending on the backbone rigidity. The interaggregate structures on a large scale are featured by the static structure factor (SSF). The simulated SSFs of the bending bundle-like aggregates are consistent with those of the small angle X-ray scattering (SAXS) measurement so we successfully assign the microscopic structures of polyelectrolytes to the SAXS measurement. The power-law of the SSFs for the bundle conditions is steeper than that of the conventional cylinder model. The present study finds that such discrepancy in the power-law results from the bending of the bundle-like aggregates. In addition, the relaxation behavior includes slow dynamics. The present study proposes that such slow dynamics results from diffusion-limited aggregation and from gliding processes to reduce local metastable folding within the aggregates.

Assembly of Triblock Amphiphilic Peptides into One-Dimensional Aggregates and Network Formation

Peptide assembly plays a key role in both neurological diseases and development of novel biomaterials with well-defined nanostructures. Synthetic model peptides provide a unique platform to explore the role of intermolecular interactions in the assembly process. A triblock peptide architecture designed by the Hartgerink group is a versatile system which relies on Coulomb interactions, hydrogen bonding, and hydrophobicity to guide these peptides' assembly at three different length scales: β-sheets, double-wall ribbon-like aggregates, and finally a highly porous network structure which can support gels with ≤1% by weight peptide concentration. In this study, by using molecular dynamics simulations of a structure based implicit solvent coarse grained model, we analyzed this hierarchical assembly process. Parametrization of our CG model is based on multiple-state points from atomistic simulations, which enables this model to represent the conformational adaptability of the triblock peptide molecule based on the surrounding medium. Our results indicate that emergence of the double-wall β-sheet packing mechanism, proposed in light of the experimental evidence, strongly depends on the subtle balance of the intermolecular forces. We demonstrate that, even though backbone hydrogen bonding dominates the early nucleation stages, depending on the strength of the hydrophobic and Coulomb forces, alternative structures such as zero-dimensional aggregates with two β-sheets oriented orthogonally (which we refer to as a cross-packed structure) and β-sheets with misoriented hydrophobic side chains are also feasible. We discuss the implications of these competing structures for the three different length scales of assembly by systematically investigating the influence of density, counterion valency, and hydrophobicity.

A route to self-assemble suspended DNA nano-complexes

Highly charged polyelectrolytes can self-assemble in presence of condensing agents such as multivalent cations, amphiphilic molecules or proteins of opposite charge. Aside precipitation, the formation of soluble micro- and nano-particles has been reported in multiple systems. However a precise control of experimental conditions needed to achieve the desired structures has been so far hampered by the extreme sensitivity of the samples to formulation pathways. Herein we combine experiments and molecular modelling to investigate the detailed microscopic dynamics and the structure of self-assembled hexagonal bundles made of short dsDNA fragments complexed with small basic proteins. We suggest that inhomogeneous mixing conditions are required to form and stabilize charged self-assembled nano-aggregates in large excess of DNA. Our results should help re-interpreting puzzling behaviors reported for a large class of strongly charged polyelectrolyte systems.

Generic, phenomenological, on-the-fly renormalized repulsion model for self-limited organization of terminal supraparticle assemblies

Self-limited, or terminal, supraparticles have long received great interest because of their abundance in biological systems (DNA bundles and virus capsids) and their potential use in a host of applications ranging from photonics and catalysis to encapsulation for drug delivery. Moreover, soft, uniform colloidal aggregates are a promising candidate for quasicrystal and other hierarchical assemblies. In this work, we present a generic coarse-grained model that captures the formation of self-limited assemblies observed in various soft-matter systems including nanoparticles, colloids, and polyelectrolytes. Using molecular dynamics simulations, we demonstrate that the assembly process is self-limited when the repulsion between the particles is renormalized to balance their attraction during aggregation. The uniform finite-sized aggregates are further shown to be thermodynamically stable and tunable with a single dimensionless parameter. We find large aggregates self-organize internally into a core-shell morphology and exhibit anomalous uniformity when the constituent nanoparticles have a polydisperse size distribution.

Polyelectrolyte properties of filamentous biopolymers and their consequences in biological fluids

Anionic polyelectrolyte filaments are common in biological cells. DNA, RNA, the cytoskeletal filaments F-actin, microtubules, and intermediate filaments, and polysaccharides such as hyaluronan that form the pericellular matrix all have large net negative charge densities distributed over their surfaces. Several filamentous viruses with diameters and stiffnesses similar to those of cytoskeletal polymers also have similar negative charge densities. Extracellular protein filaments such collagen, fibrin and elastin, in contrast, have notably smaller charge densities and do not behave as highly charged polyelectrolytes in solution. This review summarizes data that demonstrate generic counterion-mediated effects on four structurally unrelated biopolymers of similar charge density: F-actin, vimentin, Pf1 virus, and DNA, and explores the possible biological and pathophysiological consequences of the polyelectrolyte properties of biological filaments.

Electrostatic repulsion-driven crystallization model arising from filament networks

The crystallization of bundles in filament networks interacting via long-range repulsions in confinement is described by a phenomenological model. The model demonstrates the formation of the hexagonal crystalline order via the interplay of the confinement potential and the filament-filament repulsion. Two distinct crystallization mechanisms in the short- and large- screening length regimes are discussed, and the phase diagram is obtained. Simulation of large bundles predicts the existence of topological defects within the bundled filaments. This electrostatic repulsion-driven crystallization model arising from studying filament networks can even find a more general context extending to charged colloidal systems.

Solvation effects on like-charge attraction

We present results of molecular dynamics simulations of the electrostatic interaction between two parallel charged rods in the presence of divalent counterions. Such polyelectrolytes have been considered as a simple model for understanding electrostatic interactions in highly charged biomolecules such as DNA. Since there are correlations between the free charge carriers, the phenomenon of like charge attraction appears for specific parameters. We explore the role of solvation effects and the resulting deviations from Coulomb's law on the nanoscale on this peculiar phenomenon. The behavior of the force between the charged rods in a simulation with atomistic representation of water molecules is completely different from a model in which water is modeled as a continuum dielectric. By calculating counterion-rodion pair correlation functions, we find that the presence of water molecules changes the structure of the counterion cloud and results in both qualitative and quantitative changes of the force between highly charged polyelectrolytes.

Ensemble equivalence for counterion condensation on a two-dimensional charged disk

We study the counterion condensation on a two-dimensional charged disk in the limit of infinite dilution, and compare the energy-temperature relation obtained from the canonical free energy and microcanonical entropy. The microcanonical entropy is piecewise linear in energy, and is shown to be concave for all energies. As a result, even though the interactions are long-ranged, the energy-temperature relation and hence the counterion condensation transition points are identical in both the ensembles.

Counterion condensation theory of attraction between like charges in the absence of multivalent counterions

There is abundant experimental evidence suggesting the existence of attractive interactions among identically charged polyelectrolytes in ordinary salt solutions. The presence of multivalent counterions is not required. We review the relevant literature in detail and conclude that it merits more attention than it has received. We discuss also some recent observations of a low ionic strength attraction of negatively charged DNA to the region of a negatively charged glass nanoslit where the floor of the nanoslit meets the walls, again in the absence of multivalent ions. On the theoretical side, it has become clear that purely electrostatic interactions require the presence of multivalent counterions if they are to generate like-charge attraction. Any theory of like-charge attraction in the absence of multivalent counterions must therefore contain a non-electrostatic component. We point out that counterion condensation theory, which has predicted like-charge polyelectrolyte attraction in an intermediate range of distances in ordinary 1:1 salt conditions, contains both electrostatic and non-electrostatic elements. The non-electrostatic component of the theory is the modeling constraint that the counterions fall into two explicit populations, condensed and uncondensed. As reviewed in the paper, this physically motivated constraint is supported by strong experimental evidence. We proceed to offer an explanation of the nanoslit observations by showing in an idealized model that the line of intersection of two intersecting planes is a virtual polyelectrolyte. Since we have previously developed a counterion condensation theory of attraction of two like-charged polyelectrolytes, our suggestion is that the DNA is attracted to the virtual polyelectrolytes that may be located in the nanoslit where floor meets walls. We present the detailed calculations needed to document this suggestion: an extension of previous theory to the case of polyelectrolytes with like but not identical charges; the demonstration of counterion condensation on a plane with bare charge density greater than an explicitly exhibited critical value; a calculation of the free energy of the plane; a calculation of the interaction of a line charge polyelectrolyte with a like-charged plane; and the detailed demonstration that the line of intersection of two planes is a virtual polyelectrolyte.

Computer simulations of the growth of synthetic peptide fibres

We present a coarse-grained computer model designed to study the growth of fibres in a synthetic self-assembling peptide system. The system consists of two 28 residue α-helical sequences, denoted AB and CD, in which the interactions between the half peptides, A, B, C and D, may be tuned individually to promote different types of growth behaviour. In the model, AB and CD are represented by double ended rods, with interaction sites distributed along their lengths. Monte Carlo simulations are performed to follow fibre growth. It is found that lateral and longitudinal growth of the fibre are governed by different mechanisms--the former is diffusion limited with a very small activation energy for the addition of units, whereas the latter occurs via a process of secondary nucleation at the fibre ends. As a result, longitudinal growth generally proceeds more slowly than lateral growth. Furthermore, it is shown that the aspect ratio of the growing fibre may be controlled by adjusting the temperature and the relative strengths of the interactions. The predictions of the model are discussed in the context of published data from real peptide systems.