Persistence Lengths and Structure Factors of Wormlike Polymers under Confinement

Unravelled proteins form blobs during translocation across nanopores

The electroosmotic-driven transport of unravelled proteins across nanopores is an important biological process that is now under investigation for the rapid analysis and sequencing of proteins. For this approach to work, however, it is crucial that the polymer is threaded in single file. Here we found that, contrary to the electrophoretic transport of charged polymers such as DNA, during polypeptide translocation blob-like structures typically form inside nanopores. Comparisons between different nanopore sizes, shapes and surface chemistries showed that under electroosmotic-dominated regimes single-file transport of polypeptides can be achieved using nanopores that simultaneously have an entry and an internal diameter that is smaller than the persistence length of the polymer, have a uniform non-sticky ( i . e . non-aromatic) nanopore inner surface, and using moderate translocation velocities.

Lysozyme stability in various deep eutectic solvents using molecular dynamics simulations

The ability of neat deep eutectic solvents (DESs) to influence protein structure and function has gained considerable interest due to the unstable nature of enzymes or therapeutic proteins, which are often exposed to thermal, chemical, or mechanical stresses when handled at an industrial scale. In this study, we simulated a model globular protein, lysozyme, in water and six choline chloride-based DES using molecular dynamics simulations, to investigate the structural changes in various solvent environments, giving insights into the overall stability of lysozyme. Root mean square deviation (RMSD) and root mean square fluctuations (RMSF) of the C-α backbone indicated that most DESs induced a less flexible and rigid lysozyme structure compared to water. The radius of gyration and end-to-end distance calculations pointed towards higher structural compactness in reline and levuline, while the structure of lysozyme considerably expanded in oxaline. Protein-solvent interactions were further analysed by hydrogen bonding interactions and radial distribution functions (RDF), which indicated a higher degree of lysozyme-hydrogen bond donor (HBD) interactions compared to lysozyme-choline hydrogen bonding. Surface area analysis revealed an overall % increase in total positive, negative, donor, and acceptor surface areas in malicine and oxaline compared to water and other DESs, indicating the exposure of a larger number of residues to interactions with the solvent. Reline, levuline, and polyol-based DESs comparatively stabilized lysozyme, even though changes in the secondary/tertiary structures were observed.Communicated by Ramaswamy H. Sarma.

The organization of double-stranded RNA in the chikungunya virus replication organelle

Alphaviruses are mosquito-borne, positive-sense single-stranded RNA viruses. Amongst the alphaviruses, chikungunya virus is notable as a large source of human illness, especially in tropical and subtropical regions. When they invade a cell, alphaviruses generate dedicated organelles for viral genome replication, so-called spherules. Spherules form as outward-facing buds at the plasma membrane, and it has recently been shown that the thin membrane neck that connects this membrane bud with the cytoplasm is guarded by a two-megadalton protein complex that contains all the enzymatic functions necessary for RNA replication. The lumen of the spherules contains a single copy of the negative-strand template RNA, present in a duplex with newly synthesized positive-sense RNA. Less is known about the organization of this double-stranded RNA as compared to the protein components of the spherule. Here, we analyzed cryo-electron tomograms of chikungunya virus spherules in terms of the organization of the double-stranded RNA replication intermediate. We find that the double-stranded RNA has a shortened apparent persistence length as compared to unconstrained double-stranded RNA. Around half of the genome is present in either of five conformations identified by subtomogram classification, each representing a relatively straight segment of ~25-32 nm. Finally, the RNA occupies the spherule lumen at a homogeneous density, but has a preferred orientation to be perpendicular to a vector pointing from the membrane neck towards the spherule center. Taken together, this analysis lays another piece of the puzzle of the highly coordinated alphavirus genome replication.

Constitutive Model of Erythrocyte Membranes with Distributions of Spectrin Orientations and Lengths

We present an analytical hyperelastic constitutive model of the red blood cell (erythrocyte) membrane based on recently improved characterizations of density and microscopic structure of its spectrin network from proteomics and cryo-electron tomography. The model includes distributions of both orientations and natural lengths of spectrin and updated copy numbers of proteins. By applying finite deformation to the spectrin network, we obtain the total free energy and stresses in terms of invariants of shear and area deformation. We generalize an expression of the initial shear modulus, which is independent of the number of molecular orientations within the network and also derive a simplified version of the model. We apply the model and its simplified version to analyze micropipette aspiration computationally and analytically and explore the effect of local cytoskeletal density change. We also explore the discrepancies among shear modulus values measured using different experimental techniques reported in the literature. We find that the model exhibits hardening behavior and can explain many of these discrepancies. Moreover, we find that the distribution of natural lengths plays a crucial role in the hardening behavior when the correct copy numbers of proteins are used. The initial shear modulus values we obtain using our current model (5.9-15.6 pN/μm) are close to the early estimates (6-9 pN/μm). This new, to our knowledge, constitutive model establishes a direct connection between the molecular structure of spectrin networks and constitutive laws and also defines a new picture of a much denser spectrin network than assumed in prior studies.

Investigation of nematic to smectic phase transition and dynamical properties of strongly confined semiflexible polymers using Langevin dynamics

We investigated the nematic to smectic phase transition for strongly confined semiflexible polymer solutions in slit-like confinements using GPU-accelerated Langevin dynamics. We characterized the phase transitions from the nematic to smectic phases for semi-flexible polymer solutions as the polymer density increased. The dependence for the lyotropic nematic to smectic transition can be collapsed by scaling exponents between 0.2 and 0.3. The smectic C phase is found for all the cases with the polymer orientation director tilted with respect to smectic layer lateral alignment. As the chain rigidity increases, the transition density decreases for systems in which the polymer persistence length (P) to slit height (H) ratios are 1.25, 2.5, 3.75, 5 and 25. We also characterized the polymer dynamics for the isotropic-nematic-smectic transitions. The overall polymer diffusivity decreased steadily as the polymer density increased. We observed anomalous polymer diffusion along the nematic director near the isotropic-nematic transition, similar to previously reported behavior for nematic-forming ellipsoids. Polymer diffusivity decreased sharply by two orders of magnitude upon the nematic-smectic transition.

Charging and Release Mechanisms of Flexible Macromolecules in Droplets

We study systematically the charging and release mechanisms of a flexible macromolecule, modeled by poly(ethylene glycol) (PEG), in a droplet by using molecular dynamics simulations. We compare how PEG is solvated and charged by sodium Na⁺ ions in a droplet of water (H₂O), acetonitrile (MeCN), and their mixtures. Initially, we examine the location and the conformation of the macromolecule in a droplet bearing no net charge. It is revealed that the presence of charge carriers do not affect the location of PEG in aqueous and MeCN droplets compared with that in the neutral droplets, but the location of the macromolecule and the droplet size do affect the PEG conformation. PEG is charged on the surface of a sodiated aqueous droplet that is found close to the Rayleigh limit. Its charging is coupled to the extrusion mechanism, where PEG segments leave the droplet once they coordinate a Na⁺ ion or in a correlated motion with Na⁺ ions. In contrast, as PEG resides in the interior of a MeCN droplet, it is sodiated inside the droplet. The compact macro-ion transitions through partially unwound states to an extended conformation, a process occurring during the final stage of desolvation and in the presence of only a handful of MeCN molecules. For charged H₂O/MeCN droplets, the sodiation of PEG is determined by the H₂O component, reflecting its slower evaporation and preference over MeCN for solvating Na⁺ ions. We use the simulation data to construct an analytical model that suggests that the droplet surface electric field may play a role in the macro-ion-droplet interactions that lead to the extrusion of the macro-ion. This study provides the first evidence of the effect of the surface electric field by using atomistic simulations. Graphical Abstract ᅟ.

Knot Energy, Complexity, and Mobility of Knotted Polymers

The Coulomb energy E C is defined by the energy required to charge a conductive object and scales inversely to the self-capacity C, a basic measure of object size and shape. It is known that C is minimized for a sphere for all objects having the same volume, and that C increases as the symmetry of an object is reduced at fixed volume. Mathematically similar energy functionals have been related to the average knot crossing number 〈m〉, a natural measure of knot complexity and, correspondingly, we find E C to be directly related to 〈m〉 of knotted DNA. To establish this relation, we employ molecular dynamics simulations to generate knotted polymeric configurations having different length and stiffness, and minimum knot crossing number values m for a wide class of knot types relevant to the real DNA. We then compute E C for all these knotted polymers using the program ZENO and find that the average Coulomb energy 〈E C〉 is directly proportional to 〈m〉. Finally, we calculate estimates of the ratio of the hydrodynamic radius, radius of gyration, and the intrinsic viscosity of semi-flexible knotted polymers in comparison to the linear polymeric chains since these ratios should be useful in characterizing knotted polymers experimentally.

Stretching of Tethered DNA in Nanoslits

External forces and confinement are two fundamental and complementary approaches for biopolymer stretching. By employing micro- and nanofluidics, we study the force-extension dynamics by simultaneously applying external forces and confinement to single-DNA molecules. In particular, we apply external electric fields to stretch single DNA molecules that are attached to microspheres anchored at a nanoslit entrance. Using this method, we measure the force-extension relation of tethered DNA and describe this relation with modified wormlike chain models. This allowed experimental validations of several theoretical predictions, including the increase in the global persistence length of confined DNA with increasing degree of confinement and the "confined Pincus" regime in slit confinement.

Comparison of a stripe and slab confinement for ring and linear macromolecules in nanochannel

The combined effects of the channel asymmetry and the closed chain topology on the chain extension, structure factor, and the orientation correlations were studied using coarse-grained molecular dynamics simulations for moderate chain lengths. These effects are related to applications in linearization experiments with a DNA molecule in nanofluidic devices. According to the aspect ratio, the channels are classified as a stripe or slabs. The chain segments do not have any freedom to move in the direction of the narrowest stripe size, being approximately the same size as the segment size. The chains of both ring and linear topologies are extended more in a stripe than in a slab; this effect is strengthened for a ring. For a ring in a stripe, the extension-confinement strength dependence leads to effective Flory exponents even larger than 3/4, which is characteristic for a self-avoiding two-dimensional chain. While the chain extension-confinement strength dependence for both topologies conforms to the de Gennes regime in a stripe, a linear chain undergoes gradual transition to the pseudoideal regime as the slab height increases in the slab-like confinement. For a confined circle, the onset of the pseudoideal regime is shifted to larger slab heights. The structure factor confirms the absence of the pseudoideal and extended de Gennes regime in a stripe and the transition from the extended to the pseudoideal regime of a circular and linear chain upon increasing the slab heights.

Density Functional Theory of Polymer Structure and Conformations

We present a density functional approach to quantitatively evaluate the microscopic conformations of polymer chains with consideration of the effects of chain stiffness, polymer concentration, and short chain molecules. For polystyrene (PS), poly(ethylene oxide) (PEO), and poly(methyl methacrylate) (PMMA) melts with low-polymerization degree, as chain length increases, they display different stretching ratios and show non-universal scaling exponents due to their different chain stiffnesses. In good solvent, increase of PS concentration induces the decline of gyration radius. For PS blends containing short (m1=1−100) and long (m=100) chains, the expansion of long chains becomes unobvious once m1 is larger than 40, which is also different to the scaling properties of ideal chain blends.

Confined semiflexible polymers suppress fluctuations of soft membrane tubes

We use Monte Carlo computer simulations to investigate tubular membrane structures with and without semiflexible polymers confined inside. At small values of membrane bending rigidity, empty fluid and non-fluid membrane tubes exhibit markedly different behavior, with fluid membranes adopting irregular, highly fluctuating shapes and non-fluid membranes maintaining extended tube-like structures. Fluid membranes, unlike non-fluid membranes, exhibit a local maximum in specific heat as their bending rigidity increases. The peak is coincident with a transition to extended tube-like structures. We further find that confining a semiflexible polymer within a fluid membrane tube reduces the specific heat of the membrane, which is a consequence of suppressed membrane shape fluctuations. Polymers with a sufficiently large persistence length can significantly deform the membrane tube, with long polymers leading to localized bulges in the membrane that accommodate regions in which the polymer forms loops. Analytical calculations of the energies of idealized polymer-membrane configurations provide additional insight into the formation of polymer-induced membrane deformations.

Ejecting phage DNA against cellular turgor pressure

We examine in vivo ejection of noncondensed DNA from tailed bacteriophages into bacteria. The ejection is dominantly governed by the physical conditions in the bacteria. The confinement of the DNA in the virus capsid only slightly helps the ejection, becoming completely irrelevant during its last stages. A simple calculation based on the premise of condensed DNA in the cell enables us to estimate the maximal bacterial turgor pressure against which the ejection can still be fully realized. The calculated pressure (~5 atm) shows that the ejection of DNA into Gram-negative bacteria could proceed spontaneously, i.e., without the need to invoke active mechanisms.

Force-Extension for DNA in a Nanoslit: Mapping between the 3D and 2D Limits

The force-extension relation for a semiflexible polymer confined in a nanoslit is investigated. Both the effective correlation length and force-extension relation change as the chain goes from 3D (large slit heights) to 2D (tight confinement). At low forces, correlations along the polymer give an effective dimensionality. The strong force limit can be interpolated with the weak force limit for two regimes: when confinement dominates over extensile force and vice versa. These interpolations give good agreement with simulations for all slit heights and forces. We thus generalize the Marko-Siggia force-extension relation for DNA and other semiflexible biopolymers in nanoconfinement.

Confined polyelectrolytes: The complexity of a simple system

The interaction between polyelectrolytes and counterions in confined situations and the mutual relationship between chain conformation and ion condensation is an important issue in several areas. In the biological field, it assumes particular relevance in the understanding of the packaging of nucleic acids, which is crucial in the design of gene delivery systems. In this work, a simple coarse-grained model is used to assess the cooperativity between conformational change and ion condensation in spherically confined backbones, with capsides permeable to the counterions. It is seen that the variation on the degree of condensation depends on counterion valence. For monovalent counterions, the degree of condensation passes through a minimum before increasing as the confining space diminishes. In contrast, for trivalent ions, the overall tendency is to decrease the degree of condensation as the confinement space also decreases. Most of the particles reside close to the spherical wall, even for systems in which the density is higher closer to the cavity center. This effect is more pronounced, when monovalent counterions are present. Additionally, there are clear variations in the charge along the concentric layers that cannot be totally ascribed to polyelectrolyte behavior, as shown by decoupling the chain into monomers. If both chain and counterions are confined, the formation of a counterion rich region immediately before the wall is observed. Spool and doughnut-like structures are formed for stiff chains, within a nontrivial evolution with increasing confinement.

DNA conformation in nanochannels: Monte Carlo simulation studies using a primitive DNA model

We have performed canonical ensemble Monte Carlo simulations of a primitive DNA model to study the conformation of 2.56 ~ 21.8 μm long DNA molecules confined in nanochannels at various ionic concentrations with the comparison of our previous experimental findings. In the model, the DNA molecule is represented as a chain of charged hard spheres connected by fixed bond length and the nanochannels as planar hard walls. System potentials consist of explicit electrostatic potential along with short-ranged hard-sphere and angle potentials. Our primitive model system provides valuable insight into the DNA conformation, which cannot be easily obtained from experiments or theories. First, the visualization and statistical analysis of DNA molecules in various channel dimensions and ionic strengths verified the formation of locally coiled structures such as backfolding or hairpin and their significance even in highly stretched states. Although the folding events mostly occur within the region of ~0.5 μm from both chain ends, significant portion of the events still take place in the middle region. Second, our study also showed that two controlling factors such as channel dimension and ionic strength widely used in stretching DNA molecules have different influence on the local DNA structure. Ionic strength changes local correlation between neighboring monomers by controlling the strength of electrostatic interaction (and thus the persistence length of DNA), which leads to more coiled local conformation. On the other hand, channel dimension controls the overall stretch by applying the geometric constraint to the non-local DNA conformation instead of directly affecting local correlation. Third, the molecular weight dependence of DNA stretch was observed especially in low stretch regime, which is mainly due to the fact that low stretch modes observed in short DNA molecules are not readily accessible to much longer DNA molecules, resulting in the increase in the stretch of longer DNA molecules.

Conformation Model of Back-Folding and Looping of a Single DNA Molecule Confined Inside a Nanochannel

Currently, no theory is available to describe the conformation of DNA confined in a channel when the nanochannel diameter is around the persistence length. Back-folded hairpins in the undulating wormlike chain conformation result in the formation of loops, which reduces the stretch of the molecule in the longitudinal direction of the channel. A cooperativity model is used to quantify the frequency and size of the loop domains. The predictions agree with results from the Monte Carlo simulation.

Intermediate filaments in small configuration spaces

Intermediate filaments play a key role in cell mechanics. Apart from their great importance from a biomedical point of view, they also act as a very suitable micrometer-sized model system for semiflexible polymers. We perform a statistical analysis of the thermal fluctuations of individual filaments confined in microchannels. The small channel width and the resulting deflections at the walls give rise to a reduction of the configuration space by about 2 orders of magnitude. This circumstance enables us to precisely measure the intrinsic persistence length of vimentin intermediate filaments and to show that they behave as ideal wormlike chains; we observe that small fluctuations in perpendicular planes decouple. Furthermore, the inclusion of results for confined actin filaments demonstrates that the Odijk confinement regime is valid over at least 1 order of magnitude in persistence length.

Effective stiffening of DNA due to nematic ordering causes DNA molecules packed in phage capsids to preferentially form torus knots

Observation that DNA molecules in bacteriophage capsids preferentially form torus type of knots provided a sensitive gauge to evaluate various models of DNA arrangement in phage heads. Only models resulting in a preponderance of torus knots could be considered as close to reality. Recent studies revealed that experimentally observed enrichment of torus knots can be qualitatively reproduced in numerical simulations that include a potential inducing nematic arrangement of tightly packed DNA molecules within phage capsids. Here, we investigate what aspects of the nematic arrangement are crucial for inducing formation of torus knots. Our results indicate that the effective stiffening of DNA by the nematic arrangement not only promotes knotting in general but is also the decisive factor in promoting formation of DNA torus knots in phage capsids.

Critical polyelectrolyte adsorption under confinement: planar slit, cylindrical pore, and spherical cavity

We explore the properties of adsorption of flexible polyelectrolyte chains in confined spaces between the oppositely charged surfaces in three basic geometries. A method of approximate uniformly valid solutions for the Green function equation for the eigenfunctions of polymer density distributions is developed to rationalize the critical adsorption conditions. The same approach was implemented in our recent study for the "inverse" problem of polyelectrolyte adsorption onto a planar surface, and on the outer surface of rod-like and spherical obstacles. For the three adsorption geometries investigated, the theory yields simple scaling relations for the minimal surface charge density that triggers the chain adsorption, as a function of the Debye screening length and surface curvature. The encapsulation of polyelectrolytes is governed by interplay of the electrostatic attraction energy toward the adsorbing surface and entropic repulsion of the chain squeezed into a thin slit or small cavities. Under the conditions of surface-mediated confinement, substantially larger polymer linear charge densities are required to adsorb a polyelectrolyte inside a charged spherical cavity, relative to a cylindrical pore and to a planar slit (at the same interfacial surface charge density). Possible biological implications are discussed briefly in the end.

Brownian dynamics simulation of nucleocytoplasmic transport: a coarse-grained model for the functional state of the nuclear pore complex

The nuclear pore complex (NPC) regulates molecular traffic across the nuclear envelope (NE). Selective transport happens on the order of milliseconds and the length scale of tens of nanometers; however, the transport mechanism remains elusive. Central to the transport process is the hydrophobic interactions between karyopherins (kaps) and Phe-Gly (FG) repeat domains. Taking into account the polymeric nature of FG-repeats grafted on the elastic structure of the NPC, and the kap-FG hydrophobic affinity, we have established a coarse-grained model of the NPC structure that mimics nucleocytoplasmic transport. To establish a foundation for future works, the methodology and biophysical rationale behind the model is explained in details. The model predicts that the first-passage time of a 15 nm cargo-complex is about 2.6±0.13 ms with an inverse Gaussian distribution for statistically adequate number of independent Brownian dynamics simulations. Moreover, the cargo-complex is primarily attached to the channel wall where it interacts with the FG-layer as it passes through the central channel. The kap-FG hydrophobic interaction is highly dynamic and fast, which ensures an efficient translocation through the NPC. Further, almost all eight hydrophobic binding spots on kap-β are occupied simultaneously during transport. Finally, as opposed to intact NPCs, cytoplasmic filaments-deficient NPCs show a high degree of permeability to inert cargos, implying the defining role of cytoplasmic filaments in the selectivity barrier.

Perforating and spanning the nuclear envelope (NE), the nuclear pore complex (NPC) is a supramolecular assembly that regulates all traffic between the nucleus and cytoplasm. As the unique gateway to the nucleus, NPC selectively facilitates the transport of large cargo while offering a relatively unobstructed pathway for small molecules and ions. Despite the high throughput of about 1000 translocations per NPC per second, the NPC strictly controls the passage of individual cargos. However, the dynamic mechanism of nucleocytoplasmic transport is poorly understood. It is too difficult to experiment on the transport mechanism within the confined geometry of this tiny pore in vivo. Currently, only computational techniques can elucidate the detailed events happening at this tiny pore with a refined spatiotemporal resolution to account for transient bonds. Based on experimental data regarding the NPC structure and nucleocytoplasmic transport, we have established a coarse-grained model of the functional state of the NPC. The model mimics nucleocytoplasmic transport and allows us to directly observe the processes happening within the pore from a biophysical perspective. The first-passage time of a single cargo-complex is found to be about 2.6 ms. Furthermore, kap-FG hydrophobic bonds are highly dynamic and short-lived, ensuring efficient transport.

A Replica Exchange Molecular Dynamics Simulation of a Single Polyethylene Chain: Temperature Dependence of Structural Properties and Chain Conformational Study at the Equilibrium Melting Temperature

The conformational properties of a finite length polyethylene chain were explored over a wide range of temperatures using a replica exchange molecular dynamics simulation providing high quality simulation data representative for the equilibrium behavior of the chain molecule. The radial distribution function (RDF) and the structure factor S(q) of the chain as a function of temperature are analyzed in detail. The different characteristic peaks in the RDF and S(q) were assigned to specific distances in the chain and structural changes occurring with the temperature. In S(q), a peak characteristic for the order in the solid state was found and used to determine the equilibrium melting temperature. A detailed scaling analysis of the structure factor covering the full q range was performed according to the work of Hammouda. In the Θ region, a quantitative analysis of the full structure factor was done using the equivalent Kuhn chain, which enabled us to assign the Θ region of our chain and to demonstrate, in our particular case, the failure of the Gaussian chain approach. The chain conformational properties at the equilibrium melting temperature are discussed using conformational distribution functions, using the largest principal component of the radius of gyration and shape parameters as order parameters. We demonstrate that for the system studied here, the Landau free energy expression based on this conformational distribution information leads to erroneous conclusions concerning the thermodynamic transition behavior. Finally, we focus on the instantaneous conformational properties at the equilibrium melting temperature and give a detailed analysis of the conformational shapes using different shape parameters and a simulation snapshot. We show that the chain does not only take the lamellar rod-like and globular conformational shapes, typical of the solid and liquid states, but can also explore many other conformational states, including the toroidal conformational state. It is the first demonstration that a flexible molecule like PE can also take a toroidal conformational state, which is normally linked to stiffer chains.

Micelle-like architecture of the monomer ensemble of Alzheimer's amyloid-β peptide in aqueous solution and its implications for Aβ aggregation

Aggregation of amyloid-β (Aβ) peptide, a 39- to 43-residue fragment of the amyloid precursor protein, is associated with Alzheimer's disease, the most common form of dementia in the elderly population. Several experimental studies have tried to characterize the atomic details of amyloid fibrils, which are the final product of Aβ aggregation. Much less is known about species forming during the early stages of aggregation, in particular about the monomeric state of the Aβ peptide that may be viewed as the product of the very first step in the hypothesized amyloid cascade. Here, the equilibrium ensembles of monomeric Aβ alloforms Aβ(1-40) and Aβ(1-42) are investigated by Monte Carlo simulations using an atomistic force field and implicit solvent model that have been shown previously to correctly reproduce the ensemble properties of other intrinsically disordered polypeptides. Our simulation results indicate that at physiological temperatures, both alloforms of Aβ assume a largely collapsed globular structure. Conformations feature a fluid hydrophobic core formed, on average, by contacts both within and between the two segments comprising residues 12-21 and 24-40/42, respectively. Furthermore, the 11 N-terminal residues are completely unstructured, and all charged side chains, in particular those of Glu22 and Asp23, remain exposed to solvent. Taken together, these observations indicate a micelle-like† architecture at the monomer level whose implications for oligomerization, as well as fibril formation and elongation, are discussed. We establish quantitatively the intrinsic disorder of Aβ and find the propensity to form regular secondary structure to be low but sequence specific. In the presence of a global and unspecific bias for backbone conformations to populate the β-basin, the β-sheet propensity along the sequence is consistent with the arrangement of the monomer within the fibril, as derived from solid-state NMR data. These observations indicate that the primary sequence partially encodes fibril structure, but that fibril elongation must be thought of as a templated assembly step.

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.

Energy-driven asymmetric partitioning of a semiflexible polymer between interconnected cavities

The distribution of a semiflexible chain in the volume of two interconnected spherical cavities of equal size has been investigated by using Monte Carlo simulations. The chain possessed an extension exceeding that of the cavity, leading to large probabilities of translocated states despite the entropic penalty of passing the narrow passage. Furthermore, an asymmetric state with unequal subchain lengths in the two cavities was more favorable than the symmetric state. The preference for the asymmetric state is driven by the bending energy. Basically, in the symmetric state both subchains are forced to be bent, whereas in the asymmetric case only one of the subchains must bend, leading to an overall smaller bending penalty and overall smaller free energy of the asymmetric state. These results are in contrast to the entropy-controlled partitioning of polymers into confinement and the symmetric translocation state appearing for flexible polymers.