Kinetics of loop formation in polymer chains

We investigate the kinetics of loop formation in ideal flexible polymer chains (the Rouse model), and polymers in good and poor solvents. We show for the Rouse model, using a modification of the theory of Szabo, Schulten, and Schulten, that the time scale for cyclization is tau(c) approximately tau(0)N(2) (where tau(0) is a microscopic time scale and N is the number of monomers), provided the coupling between the relaxation dynamics of the end-to-end vector and the looping dynamics is taken into account. The resulting analytic expression fits the simulation results accurately when a, the capture radius for contact formation, exceeds b, the average distance between two connected beads. Simulations also show that when a < b, tau(c) approximately N(alpha)(tau), where 1.5 < alpha(tau) < or = 2 in the range 7 < N < 200 used in the simulations. By using a diffusion coefficient that is dependent on the length scales a and b (with a < b), which captures the two-stage mechanism by which looping occurs when a < b, we obtain an analytic expression for tauc that fits the simulation results well. The kinetics of contact formation between the ends of the chain are profoundly effected when interactions between monomers are taken into account. Remarkably, for N < 100, the values of tau(c) decrease by more than 2 orders of magnitude when the solvent quality changes from good to poor. Fits of the simulation data for tau(c) to a power law in N (tau(c) approximately N(alpha)(tau)) show that alpha(tau) varies from about 2.4 in a good solvent to about 1.0 in poor solvents. The effective exponent alpha(tau) decreases as the strength of the attractive monomer-monomer interactions increases. Loop formation in poor solvents, in which the polymer adopts dense, compact globular conformations, occurs by a reptation-like mechanism of the ends of the chain. The time for contact formation between beads that are interior to the chain in good solvents changes nonmonotonically as the loop length varies. In contrast, the variation in interior loop closure time is monotonic in poor solvents. The implications of our results for contact formation in polypeptide chains, RNA, and single-stranded DNA are briefly outlined.     

Persistence lengths and structure factors of wormlike polymers under confinement

The behavior of semiflexible chains modeling wormlike polymers such as DNA and actin in confined spaces was explored by coarse-grained Monte Carlo simulations. The persistence length P, mean end-to-end distance R2, mean radius of gyration Rg2, and the size ratio R2/Rg2 were computed for chains in slits, cylinders, and spheres. It was found that the intrinsic persistence length of a free chain undergoes on confinement substantial alteration into the apparent persistence length. The qualitative differences were found in trends of the apparent persistence lengths between slits and cylinders on one side and spheres on the other side. The quantities P, R2, Rg2, and R2/Rg2 display similar dependences upon squeezing the chains in nanopores. The above quantities change nonmonotonically with confinement in slits and cylinders, whereas they drop smoothly with decreasing radius of a sphere. For elongation of a chain in a cylinder, two regimes corresponding to strong and moderate confinements were found and compared to experiments and predictions of the blob and Odijk theories. In a spherical cavity, the toroidal chain structure with a hole in the center was detected under strong confinements. The scattering form factor S(q) computed for semiflexible confined chains revealed three regimes of behavior in a slit and a cylinder that matched up well with the scaling theory. The complex form of the function S(q) computed for a sphere was interpreted as a sign of the toroidal structure. A reasonable agreement was found between the simulations and measurements of DNA and actin filaments, confined in nano- and microfluidic channels and spherical droplets, pertaining to the changes of the persistence lengths, chain elongation, and toroidal structure formation.     

Influence of surface interactions on folding and forced unbinding of semiflexible chains

We have investigated the folding and forced unbinding transitions of adsorbed semiflexible polymer chains using theory and simulations. These processes describe, at an elementary level, a number of biologically relevant phenomena that include adhesive interactions between proteins and tethering of receptors to cell walls. The binding interface is modeled as a solid surface, and the wormlike chain (WLC) is used for the semiflexible chain (SC). Using Langevin simulations, in the overdamped limit we examine the ordering kinetics of racquet-like and toroidal structures in the presence of an attractive interaction between the surface and the polymer chain. For a range of interactions, temperature, and the persistence length, l(p), we obtained the monomer density distribution, n(x), (x is the perpendicular distance of a tagged chain end from the surface) for all of the relevant morphologies. There is a single peak in n(x) inside the range of attractive forces, b, for chains in the extended conformations, whereas in racquet and toroidal structures there is an additional peak at x approximately b. The simulated results for n(x) are in good agreement with theory. The formation of toroids on the surface appears to be a first-order transition as evidenced by the bimodal distribution in n(x). The theoretical result underestimates the simulated n(x) for x < b and follows n(x) closely for x >/= b; the calculated density agrees exactly with n(x) in the range x < b. The chain-surface interaction is probed by subjecting the surface structures to a pulling force, f. The average extension, x( f), as a function of f exhibits a sigmoidal profile with sharp all-or-none transition at the unfolding force threshold f = f(c) which increases for more structured states. Simulated x(f) compare well with the theoretical predictions. The critical force, f(c), is a function of l(s)/l(c) for a fixed temperature, where l(c) and l(s) are the length scales that express the strength of the intramolecular and SC-surface attraction, respectively. For a fixed l(s), f(c) increases as l(p) decreases.     

A Monte Carlo simulation study of branched polymers

Monte Carlo simulations are presented for the static properties of highly branched polymer molecules. The molecules consist of a semiflexible backbone of hard-sphere monomers with semiflexible side chains, also composed of hard-sphere monomers, attached to either every backbone bead or every other backbone bead. The conformational properties and structure factor of this model are investigated as a function of the stiffness of the backbone and side chains. The average conformations of the side chains are similar to self-avoiding random walks. The simulations show that there is a stiffening of the backbone as degree of crowding is increased, for example, if the branch spacing is decreased or side chain length is increased. The persistence length of the backbone is relatively insensitive to the stiffness of the side chains over the range investigated. The simulations reproduce most of the qualitative features of the structure factor observed in experiment, although the magnitude of the stiffening of the backbone is smaller than in experiment.     

Mesoscopic properties of semiflexible amyloid fibrils

We have determined the contour length, persistence length, bending rigidity, and critical percolation concentration for semiflexible amyloid fibrils formed from the globular proteins beta-lactoglobulin, bovine serum albumin, and ovalbumin. The persistence length was estimated using an adjusted random contact model for highly charged semiflexible chains. We have found contour lengths in the range of 50 nm to 10 microm and persistence lengths in the range of 16 nm to 1.6 microm. This wide range of contour and persistence lengths and the ease of preparation of these amyloid fibrils make them ideal model systems for the study of semiflexible polymers.     

Monte Carlo simulation of the self-assembly and phase behavior of semiflexible equilibrium polymers

Grand canonical Monte Carlo simulations of a simple model semiflexible equilibrium polymer system, consisting of hard sphere monomers reversibly self-assembling into chains of arbitrary length, have been performed using a novel sampling method to add or remove multiple monomers during a single MC move. Systems with two different persistence lengths and a range of bond association constants have been studied. We find first-order lyotropic phase transitions between isotropic and nematic phases near the concentrations predicted by a statistical thermodynamic theory, but with significantly narrower coexistence regions. A possible contribution to the discrepancy between theory and simulation is that the length distribution of chains in the nematic phase is bi-exponential, differing from the simple exponential distribution found in the isotropic phase and predicted from a mean-field treatment of the nematic. The additional short length-scale characterizing the distribution appears to arise from the lower orientational order of short chains. The dependence of this length-scale on chemical potential, bond association constant, and total monomer concentration has been examined.     

Viscoelastic properties of an isolated polymer chain with arbitrary flexibility,chain length and hydrodynamic interactions from 4‐dimensional Dirac propagator

Kholodenko's theory of semiflexible polymer chains, the conformation and properties of which are obtained from the Dirac propagator, shows applicability to dilute solutions of semiflexible polymers of arbitrary persistence and contour lengths by calculating the static scattering function and the squared end‐to‐end distance of the polymer chain. In the present work, the theory is extended and applied to obtain the intrinsic viscosity with consideration of hydrodynamic interactions. The intrinsic viscosity formula is derived as function of chain length and persistence length. The hydrodynamic interactions are also taken into account following the Kirkwood and Riseman scheme. From this calculation, we obtain the general expression for the intrinsic viscosity and diffusion coefficients covering the whole range of chain flexibilities without confusion with the excluded volume effects. Calculated limiting values of hydrodynamical observables are in complete agreement with those known for random coils and rigid rods.     

Structure of flexible and semiflexible polyelectrolyte chains in confined spaces of slit micro/nanochannels

Understanding the behavior of a polyelectrolyte in confined spaces has direct relevance in design and manipulation of microfluidic devices, as well as transport in living organisms. In this paper, a coarse-grained model of anionic semiflexible polyelectrolyte is applied, and its structure and dynamics are fully examined with Brownian dynamics (BD) simulations both in bulk solution and under confinement between two negatively charged parallel plates. The modeling is based on the nonlinear bead-spring discretization of a continuous chain with additional long-range electrostatic, Lennard-Jones, and hydrodynamic interactions between pairs of beads. The authors also consider the steric and electrostatic interactions between the bead and the confining wall. Relevant model parameters are determined from experimental rheology data on the anionic polysaccharide xanthan reported previously. For comparison, both flexible and semiflexible models are developed accompanying zero and finite intrinsic persistence lengths, respectively. The conformational changes of the polyelectrolyte chain induced by confinements and their dependence on the screening effect of the electrolyte solution are faithfully characterized with BD simulations. Depending on the intrinsic rigidity and the medium ionic strength, the polyelectrolyte can be classified as flexible, semiflexible, or rigid. Confined flexible and semiflexible chains exhibit a nonmonotonic variation in size, as measured by the radius of gyration and end-to-end distance, with changing slit width. For the semiflexible chain, this is coupled to the variations in long-range bond vector correlation. The rigid chain, realized at low ionic strength, does not have minima in size but exhibits a sigmoidal transition. The size of confined semiflexible and rigid polyelectrolytes can be well described by the wormlike chain model once the electrostatic effects are taken into account by the persistence length measured at long length scale.     

Light scattering measurements of anisotropic viscoelastic coefficients of a main-chain polymer nematic liquid crystal

We report the first systematic measurements of elastic coefficients and viscosities that clearly demonstrate the distinction between rigid and semiflexible behaviour of a main-chain polymer nematic liquid crystal, a solution of poly-γ-benzyl glutamate (PBG) in mixed organic solvents. Quasi-elastic Rayleigh scattering studies show a crossover which occurs at a molecular chain length near the persistence length of PBG as the chain length increases. The results are in qualitatively good agreement with recent theoretical predictions for semiflexible chains. It is seen that bending distortions of an individual polymer play an important role in the fundamental nature of nematic elasticity and viscosity.     

刚柔相嵌液晶高分子的赝二级相变温度

讨论了刚柔相嵌液晶高分子的向列相－各向同性相转变与其分子结构的关系，给出了该一级相变的赝二级相变温度Ｔ与这类液晶高分子的液晶基元和间隔基的长度，柔顺性（相关长度）以及它们之间的相互作用的关系的解析表示式，分析了液晶基元与间隔基连接处的表观弯曲（接口效应）对Ｔ的影响，文中的结论与实验相符。     

刚柔相嵌液晶高分子的赝二级相变温度

讨论了刚柔相嵌液晶高分子的向列相一各向同性相转变与其分子结构的关系．给出了该一级相变的赝二级相变温度Ｔ＊与这类液晶高分子的液晶基元和间隔基的长度、柔顺性（相关长度）以及它们之间的相互作用的关系的解析表示式．分析了液晶基元与间隔基连接处的表观弯曲（接口效应）对Ｔ＊的影响．文中的结论与实验相符．     

Phase behavior of semiflexible polymer chains

Monte Carlo simulations were performed on semiflexible polymer chains with the goal of delineating their isotropic-nematic (IN) and gas-liquid coexistence envelopes. The chain monomers are spherical beads that interact via a square-well potential with all other beads. Bonded beads are connected by strings chosen so that bond length varies between 1.01sigma and 1.05sigma (where sigma is the hard sphere diameter). The stiffness of the molecules is controlled via a potential between beads separated by two bonds; this potential restricts the distance between these beads to be between 2.02sigma and 2.1sigma. The vapor-liquid coexistence and IN coexistence curves are obtained using computer simulations. An IN transition is found for 10相似文献   

Ordered Toroid Structures of Nanoparticles in Self-attractive Semi exible Polymer/Nanoparticle Composites

By employing dynamic Monte Carlo simulations, we investigate a coil-to-toroid transition of self-attractive semiflexible polymers and the spatial distributions of nanoparticles in selfattractive semiflexible polymer/nanoparticle composites. The conformation of self-attractive semiflexible polymers depends on bending energy and self-attractive interactions between monomers in polymer chains. A three-stage process of toroid formation for self-attractive semiflexible chains is shown: several isolated toroids, a loose toroid structure, and a compact toroid structure. Utilizing the compact toroid conformations of self-attractive semiflexible chains, we can control effectively the spatial distributions of nanoparticles in self-attractive semiflexible polymer nanocomposites, and an unconventional toroid structure of nanoparticles is observed.     

Computer simulation of bottle-brush polymers with flexible backbone: good solvent versus theta solvent conditions

By molecular dynamics simulation of a coarse-grained bead-spring-type model for a cylindrical molecular brush with a backbone chain of N(b) effective monomers to which with grafting density σ side chains with N effective monomers are tethered, several characteristic length scales are studied for variable solvent quality. Side chain lengths are in the range 5 ≤ N ≤ 40, backbone chain lengths are in the range 50 ≤ N(b) ≤ 200, and we perform a comparison to results for the bond fluctuation model on the simple cubic lattice (for which much longer chains are accessible, N(b) ≤ 1027, and which corresponds to an athermal, very good, solvent). We obtain linear dimensions of the side chains and the backbone chain and discuss their N-dependence in terms of power laws and the associated effective exponents. We show that even at the theta point the side chains are considerably stretched, their linear dimension depending on the solvent quality only weakly. Effective persistence lengths are extracted both from the orientational correlations and from the backbone end-to-end distance; it is shown that different measures of the persistence length (which would all agree for Gaussian chains) are not mutually consistent with each other and depend distinctly both on N(b) and the solvent quality. A brief discussion of pertinent experiments is given.     

Distribution functions for filaments under tension

We develop a biased Monte Carlo simulation technique to measure the distribution functions of the extension and the end-to-end distance of fluctuating filaments stretched by external force. The method is applicable for arbitrary ratio of the persistence length to the contour length and for arbitrary forces, and also for the case of steric constraints, such as an external wall. The fundamental idea underlying the algorithm is to account explicitly for the length-scale dependence of the effective elastic moduli. We find that orientational fluctuations and wall effects produce non-Gaussian distributions for nearly rigid filaments in the small to intermediate force regime. The simulation results are tested against analytic expressions for the force-extension curves, both in the semiflexible and nearly stiff limits.     

Influence of intrachain hydrogen bonding on the cyclization of a polystyrene chain

Cyclization of a polystyrene chain (M_n = 10,600; M_w/M_n = 1.09) both ends labeled with 4-(1-pyrenyl)butanoamide groups was studied in cyclohexane between 25 and 95°C. The amide groups (peptide bonds) at both ends can form an intrachain hydrogen bond between the amide hydrogen at one chain end and the carbonyl oxygen at the other. The presence of two sets of conformers, random coils, and chains cyclized through hydrogen bonding, complicates the data analysis. The pyrene excimer kinetics of this polymer is well described by a model composed of two monomers (hydrogen bonded and nonbonded chains) and one excimer, in equilibrium. The cyclization rate constant for hydrogen-bonded chains is larger than the one for nonhydrogen-bonded chains. The pyrene excimer binding energy (ca. 1.6 kcal/mol) is lower than the published value for nonhydrogen-bonded chains (～ 9 kcal/mol), suggesting that intrachain hydrogen bonding hinders the stabilization of the excimer. © 1994 John Wiley & Sons, Inc.     

Role of chain stiffness on the conformation of single polyelectrolytes in salt solutions

Conformation of single polyelectrolytes in tetravalent salt solutions is investigated under the framework of a coarse-grained model, using Langevin dynamics simulations. The chain size, studied by the radius of gyration, shows three different variational behaviors with salt concentration, depending on the chain stiffness. According to the size variations, polyelectrolytes of fixed chain length are classified into three categories: (1) flexible chain, for which the variation shows a curve similar to a tilted L, (2) semiflexible chain, whose curve resembles U, and (3) rigid chain, for which the curve is a straight line. The wormlike chain model with persistence length predicted by the Odijk-Skolnick-Fixman theory is found to be able to qualitatively describe the end-to-end distance at low salt concentration not only for semiflexible and rigid chains but also for flexible chain. In a low salt region, a flexible polyelectrolyte extends more significantly than a semiflexible chain, in reference of the size of their uncharged counterparts, and in a high salt region, regardless of chain stiffness, a chain attains a dimension comparable to that of its neutral polymer. The chain stiffness influences both the local and the global chain structures. A flexible chain exhibits a zigzagged local structure in the presence of salt ions, and the condensed structure is a disordered, random globule. A semiflexible chain is locally smooth, and the condensed structure is orderly packed, taking a form such as hairpin or toroid. Moreover, the chain stiffness can also affect the nature of the coil-globule transition. The transition occurred in a discrete manner for semiflexible chain, whereas it occurred in a continuous way for flexible chain. This discrete feature happened not only at low salt concentration when a semiflexible chain collapsed but also at high salt concentration when the collapsed chain is reexpanded. At the end, the effects of chain stiffness and salt concentration on the conformation of single polyelectrolytes are summarized in a schematic state diagram.     

End-monomer Dynamics in Semiflexible Polymers

Spurred by an experimental controversy in the literature, we investigate the end-monomer dynamics of semiflexible polymers through Brownian hydrodynamic simulations and dynamic mean-field theory. Precise experimental observations over the last few years of end-monomer dynamics in the diffusion of double-stranded DNA have given conflicting results: one study indicated an unexpected Rouse-like scaling of the mean squared displacement (MSD) ?r(2)(t)? ~ t(1/2) at intermediate times, corresponding to fluctuations at length scales larger than the persistence length but smaller than the coil size; another study claimed the more conventional Zimm scaling ?r(2)(t)? ~ t(2/3) in the same time range. Using hydrodynamic simulations, analytical and scaling theories, we find a novel intermediate dynamical regime where the effective local exponent of the end-monomer MSD, α(t) = d log?r(2)(t)?/d log t, drops below the Zimm value of 2/3 for sufficiently long chains. The deviation from the Zimm prediction increases with chain length, though it does not reach the Rouse limit of 1/2. The qualitative features of this intermediate regime, found in simulations and in an improved mean-field theory for semiflexible polymers, in particular the variation of α(t) with chain and persistence lengths, can be reproduced through a heuristic scaling argument. Anomalously low values of the effective exponent α are explained by hydrodynamic effects related to the slow crossover from dynamics on length scales smaller than the persistence length to dynamics on larger length scales.     

End-to-end vs interior loop formation kinetics in unfolded polypeptide chains

The conformational search for favorable intramolecular interactions during protein folding is limited by intrachain diffusion processes. Recent studies on the dynamics of loop formation in unfolded polypeptide chains have focused on loops involving residues near the chain ends. During protein folding, however, most contacts are formed between residues in the interior of the chain. We compared the kinetics of end-to-end loop formation (type I loops) to the formation of end-to-interior (type II loops) and interior-to-interior loops (type III loops) using triplet-triplet energy transfer from xanthone to naphthylalanine. The results show that formation of type II and type III loops is slower compared to type I loops of the same size and amino acid sequence. The rate constant for type II loop formation decreases with increasing overall chain dimensions up to a limiting value, at which loop formation is about 2.5-fold slower for type II loops compared to type I loops. Comparing type II loops of different loop size and amino acid sequence shows that the ratio of loop dimension over total chain dimension determines the rate constant for loop formation. Formation of type III loops is 1.7-fold slower than formation of type II loops, indicating that local chain motions are strongly coupled to motions of other chain segments which leads to faster dynamics toward the chain ends. Our results show that differences in the kinetics of formation of type I, type II, and type III loops are mainly caused by differences in internal flexibility at the different positions in the polypeptide chain. Interactions of the polypeptide chain with the solvent contribute to the kinetics of loop formation, which are strongly viscosity-dependent. However, the observed differences in the kinetics of formation of type I, type II, and type III loops are not due to the increased number of peptide-solvent interactions in type II and type III loops compared to type I loops as indicated by identical viscosity dependencies for the kinetics of formation of the different types of loops.     

Theoretical analysis of polydispersity in the nematic phase of self-assembled semiflexible chains

In recent simulations of semiflexible equilibrium polymers [X. Lu and J. T. Kindt, J. Chem. Phys. 120, 10328 (2004)] a roughly biexponential distribution of chain lengths was observed in the nematic phase. In this study, we show that a theory representing the nematic phase as an equilibrium mixture of randomly oriented chains below a critical length Lc and chains above Lc obeying a Gaussian orientational distribution reproduces this distribution qualitatively. The agreement between predicted and simulated phase boundaries is improved, with a narrower isotropic-nematic coexistence region, compared with a monoexponential chain distribution. We find, furthermore, that the critical length Lc scales as the inverse of monomer concentration, irrespective of the bond strength of the chains, but that at the phase boundary, the fraction of monomers found in disordered chains peaks at a certain bond strength.