Mesoscopic description of solvent effects on polymer dynamics

Solvent effects on polymer dynamics and structure are investigated using a mesoscopic solvent model that accounts for hydrodynamic interactions among the polymer beads. The simulation method combines molecular dynamics of the polymer chain, interacting with the solvent molecules through intermolecular forces, with mesoscopic multiparticle collision dynamics for the solvent molecules. Changes in the intermolecular forces between the polymer beads and mesoscopic solvent molecules are used to vary the solvent conditions from those for good to poor solvents. Polymer collapse and expansion dynamics following changes in solvent conditions are studied for homopolymer and block copolymer solutions. The frictional properties of polymers are also investigated.     

Foldamer dynamics expressed via Markov state models. II. State space decomposition

The structural landscape of poly-phenylacetylene (pPA), otherwise known as m-phenylene ethynylene oligomers, has been shown to consist of a very diverse set of conformations, including helices, turns, and knots. Defining a state space decomposition to classify these conformations into easily identifiable states is an important step in understanding the dynamics in relation to Markov state models. We define the state decomposition of pPA oligomers in terms of the sequence of discretized dihedral angles between adjacent phenyl rings along the oligomer backbone. Furthermore, we derive in mathematical detail an approach to further reduce the number of states by grouping symmetrically equivalent states into a single parent state. A more challenging problem requires a formal definition for knotted states in the structural landscape. Assuming that the oligomer chain can only cross the ideal helix path once, we propose a technique to define a knotted state derived from a helical state determined by the position along the helical nucleus where the chain crosses the ideal helix path. Several examples of helical states and knotted states from the pPA 12-mer illustrate the principles outlined in this article.     

Molecular dynamics simulations of helix-forming, amine-functionalized m-poly(phenyleneethynylene)s

We present here the results of all-atom and united-atom molecular dynamics (MD) simulations that were used to examine the folding behavior of an amine-functionalized m-poly(phenyleneethynylene) (m-PPE) oligomer in aqueous environment. The parallelized GROMACS MD simulation code and OPLS force field were used for multiple MD simulations of m-PPE oligomers containing 24 phenyl rings in extended, coiled and helix conformations separately in water to determine the minimum energy conformation of the oligomer in aqueous solvent and what interactions are most important in determining this structure. Simulation results showed that the helix is the preferred minimum energy conformation of a single oligomer in water and that Lennard-Jones interactions are the dominant forces for the stabilization of the helix. In addition, these solvophobic interactions are strong enough to maintain the helix conformation at temperatures up to 523 K.     

Compactization of rigid-chain amphiphilic macromolecules with local helical structure

Conformational characteristics of amphiphilic macromolecules with secondary local helical structuring are studied by the method of molecular dynamics for different properties of a helix (bending angles between neighboring vectors of the bond and internal rotation angle) and different rigidities of its fixation. Extended helices with high distances between helical turns and dense helices in which neighboring turns directly adjoin each other are studied. As the quality of a solvent deteriorates, extended helices experience a well-pronounced coil-globule transition, whose amplitude increases with an increase in chain rigidity, while the dimensions of dense helices gradually change. In a poor solvent, extended helices formed “collagen-like” structures, flexible chains of dense helices produce hairpin structures, and rigid macromolecules of dense helices form rodlike globules with an almost ideal local helical order. Independently of helix parameters, a deterioration in solvent quality leads to stabilization of the local secondary structure.     

AFM characterization of nanoscale ribbons grown from a series of oligo(p-phenyleneethynylene) derivatives

The key to optimizing the properties of molecular scale wires lies in understanding and controlling the solid-state morphologies. This paper examines the influence of oligomer chain length, solvent, and concentration on the formation of nanoscale ribbons on mica substrates from solutions of oligo(p-phenyleneethynylene)s (OPEs) with hexyloxy side chains and thioacetyl end groups. The OPEs are of different molecular chain lengths, in which the numbers ofp-dihexyloxyphenyleneethynylene repeat units, n, are 1, 3, 5, and 7, respectively, with their two ends capped with 4-thioacetylphenyl alligator groups. The atomic force microscope (AFM) is employed to investigate the thin film morphology and study the self-assembled organizations. Solvent and concentration are found to exert a strong influence on thin film morphology. Under suitable conditions, OPEs with 7 p-dihexyloxyphenyleneethynylene repeat units are driven to form micrometer-long nanoribbons, oriented preferably along the 3-fold symmetry axes of the mica substrate. The cross section of the nanoribbons is composed of 7 molecules as evaluated by AFM characterization. On the other hand, oligomers with shorter chain lengths (n = 1, 3, and 5) produce thin films featuring globular nanoaggregates, chains consisting of elongated grains, and rods, respectively. Plausible reasons for the variation in thin film morphology are discussed, based on the results obtained from investigation of oligomer chain length, solvent, and concentration effects. A subtle balance among molecular size and physicochemical properties of solute molecules, solvent molecules, and substrate is crucial for the formation of desired structures. Among them, oligomer chain length plays a key role in thin film morphology, and the critical number of repeat units in OPE/poly(p-phenyleneethynylene) molecules for the formation of nanoribbon structures with a molecular cross section is supposed to be 8 or 9.     

On the modeling of aggregates of an optically active regioregular polythiophene

The conformational properties of the optically active regioregular poly[(R)-3-(4-(4-ethyl-2-oxazolin-2-yl) phenyl) thiophene] (PEOPT) were explored by molecular dynamics on a single chain using several solvents of increasing polarity. Furthermore, their aggregate formation was studied over a wide range of temperatures using a replica exchange molecular dynamics simulation providing simulation data representative of the equilibrium behaviour of their aggregates. Results show a clear tendency of PEOPT to keep a syn-gauche conformation between continuous backbone thiophene rings favouring a bent chain structure in solvent. After studying their aggregation behaviour in acetonitrile, a strong tendency to pack stabilizing structures that reinforce the chirality of the polymer, in concordance with experimental data, was found. Two different aggregated structures were observed depending on oligomer length, a self-assembled helical aggregate based on stacked octamers and a bent double helix aggregate in large oligomers.     

Hydrogen bond-stabilized helix formation of a m-phenylene ethynylene oligomer

[reaction: see text] Incorporation of a single hydrogen bonded beta-turn mimic in the backbone of a m-phenylene ethynylene oligomer is shown to affect the thermodynamic properties of the folding reaction. Oligomers 1 and 2 both undergo solvophobic helix formation, but hydrogen bonded oligomer 1 was found to form a more stable helix with a higher tolerance to solvent denaturation than isomeric, non-hydrogen bonded oligomer 2.     

Coarse grained model of entangled polymer melts

A coarse graining procedure aimed at reproducing both the chain structure and dynamics in melts of linear monodisperse polymers is presented. The reference system is a bead-spring-type representation of the melt. The level of coarse graining is selected equal to the number of beads in the entanglement segment, Ne. The coarse model is still discrete and contains blobs each representing Ne consecutive beads in the fine scale model. The mapping is defined by the following conditions: the probability of given state of the coarse system is equal to that of all fine system states compatible with the respective coarse state, the dissipation per coarse grained object is similar in the two systems, constraints to the motion of a representative chain exist in the fine phase space, and the coarse phase space is adjusted such to represent them. Specifically, the chain inner blobs are constrained to move along the backbone of the coarse grained chain, while the end blobs move in the three-dimensional embedding space. The end blobs continuously redefine the diffusion path for the inner blobs. The input parameters governing the dynamics of the coarse grained system are calibrated based on the fine scale model behavior. Although the coarse model cannot reproduce the whole thermodynamics of the fine system, it ensures that the pair and end-to-end distribution functions, the rate of relaxation of segmental and end-to-end vectors, the Rouse modes, and the diffusion dynamics are properly represented.     

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.     

Universal consequences of the presence of excluded volume interactions in dilute polymer solutions undergoing shear flow

The role of solvent quality in determining the universal material properties of dilute polymer solutions undergoing steady simple shear flow is examined. A bead-spring chain representation of the polymer molecule is used, and the influence of solvent molecules on polymer conformations is modelled by a narrow Gaussian excluded volume potential that acts pairwise between the beads of the chain. Brownian dynamics simulations data, acquired for chains of finite length, and extrapolated to the limit of infinite chain length, are shown to be model independent. This feature of the narrow Gaussian potential, which leads to results identical to a delta-function repulsive potential, enables the prediction of both universal crossover scaling functions and asymptotic behavior in the excluded volume limit. Universal viscometric functions, obtained by this procedure, are found to exhibit increased shear thinning with increasing solvent quality. In the excluded volume limit, they are found to obey power law scaling with the characteristic shear rate beta, in close agreement with previously obtained renormalization group results. The presence of excluded volume interactions is also shown to lead to a weakening of the alignment of the polymer chain with the flow direction.     

Chain conformations and phase separation in polymer solutions with varying solvent quality

Molecular dynamics simulations are used to investigate the conformations of a single polymer chain, represented by the Kremer-Grest bead-spring model, in a solution with a Lennard-Jones liquid as the solvent when the interaction strength between the polymer and solvent is varied. Results show that when the polymer-solvent interaction is unfavorable, the chain collapses as one would expect in a poor solvent. For more attractive polymer-solvent interactions, the solvent quality improves and the chain is increasingly solvated and exhibits ideal and then swollen conformations. However, as the polymer-solvent interaction strength is increased further to be more than about twice the strength of the polymer-polymer and solvent-solvent interactions, the chain exhibits an unexpected collapsing behavior. Correspondingly, for strong polymer-solvent attractions, phase separation is observed in the solutions of multiple chains. These results indicate that the solvent becomes effectively poor again at very attractive polymer-solvent interactions. Nonetheless, the mechanism of chain collapsing and phase separation in this limit differs from the case with a poor solvent rendered by unfavorable polymer-solvent interactions. In the latter, the solvent is excluded from the domain of the collapsed chains while in the former, the solvent is still present in the pervaded volume of a collapsed chain or in the polymer-rich domain that phase separates from the pure solvent. In the limit of strong polymer-solvent attractions, the solvent behaves as a glue to stick monomers together, causing a single chain to collapse and multiple chains to aggregate and phase separate.     

Conformational Aspects and Intramolecular Phase Separation of Alternating Copolymacromonomers: A Computer Simulation Study

Summary: Using the bond fluctuation model (BFM), simulations have been performed on molecular bottle brushes with two chemically different types of side chains. In the first part of this work, rigidity of the backbone is imposed. The influence of solvent quality and side chain length on the intramolecular phase separation of side chain material, leading to Janus cylinders, has been investigated and compared to theoretical predictions. In the second part of this work, the restriction of rigidity for the backbone is relaxed. In a poor solvent, the side chain material collapses into a globular state. Several globules containing each one type of side chain material are connected together by backbone material. When imposing different solvent conditions for both types of side chains, a bending of the backbone is found as predicted by theory and observed in very recent experiments.

Comparison between typical snapshots in θ solvent. The side chain length is 25 beads. Dark and light coloured side chains repel each other.     

The mechanism of the self-assembly of associating DNA molecules under shear flow: Brownian dynamics simulation

A shear flow induces the assembly of DNAs with the sticky spots. In order to strictly interpret the mechanism of shear-induced DNA assembly, Brownian dynamics simulations with the bead-spring model were carried out for these molecules at various ranges of the Weissenberg numbers (We). We calculate a formation time and analyze the radial distribution function of end beads and the probability distribution of fractional extension at the formation time to understand the mechanism of shear-induced assembly. At low Weissenberg number the formation time, which is defined as an elapsed time until a multimer forms for the first time, decreases rapidly, reaching a plateau at We = 1000. A shear flow changes the radial distribution of end beads, which is almost the same regardless of the Weissenberg number. A shear flow deforms and stretches the molecules and generates different distributions between end beads with a stickly spot. The fractional extension progresses rapidly in shear flow from a Gaussian-like distribution to a uniform distribution. The progress of the distribution of fractional extension increases the possibility of meeting of end beads. In shear flow, the inducement of the assembly mainly results from the progress of the probability distribution of fractional extension. We also calculate properties such as the radius of gyration, stretch, and so on. As the Weissenberg number increases, the radius of gyration at the formation time also increases rapidly, reaching a plateau at We = 1000.     

NMR analysis of interaction between styrene-divinylbenzene gel beads and small molecules

The interaction of small molecules such as cyclohexane, tetrahydrofuran, and acetonitrile with styrene-divinylbenzene copolymer gel beads were investigated by ¹³C-NMR spectroscopy. When the gel content was 0.1–0.3 g/mL, the ¹³C-NMR spectrum of the solvent displayed two peaks. A sharp peak at lower magnetic field was assigned to the free solvent, and a broad one to the solvent affected by the gel. This signal splitting is attributed to the upfield shift caused by aromatic rings of styrene units in the polymer chain. The nitrile carbon of acetonitrile showed the largest upfield shift. The mobility of small molecules in gel beads was also investigated using the nuclear magnetic relaxation method. In the case of good solvents for gel beads, the mobility was affected by the cross-linking density, whereas poor solvents exhibited little dependence of the cross-linking density. © 1995 John Wiley & Sons, Inc.     

Dynamics of polymer packaging

We use the stochastic rotation dynamics algorithm to investigate the packaging of flexible and semiflexible polymers into a capsid that is permeable to solvent molecules. The model takes into account hydrodynamic interactions arising due to local flow. The flexible chain maintains a random configuration as it is being fed into the capsid, in contrast to the semiflexible chain, whose configuration is initially spool-like, becoming more random at high packing. We measure the packing rate, which is found to decrease with the percentage of the chain packed and highlight the difference between the flexible and semiflexible chains. Reflecting experiments, we find pauses in the packing process for individual chains as the motor loses grip of the fluctuating beads. We also find that hydrodynamics is important, in that the packaging rate is faster when flow is included.     

Molecular dynamics simulations of polyelectrolyte-polyampholyte complexes. Effect of solvent quality and salt concentration

Complexation between polyelectrolyte and polyampholyte chains in poor solvent conditions for the polyelectrolyte backbone has been studied by molecular dynamics simulations. In a poor solvent a polyelectrolyte forms a necklace-like structure consisting of polymeric globules (beads) connected by strings of monomers. The simulation results can be explained by assuming the existence of two different mechanisms leading to the necklace formation. In the case of weak electrostatic interactions, the necklace formation is driven by optimization of short-range monomer-monomer attraction and electrostatic repulsion between charged monomers on the polymer backbone. In the case of strong electrostatic interactions, the necklace structure appears as a result of counterion condensation. While the short-range attractions between monomers are still important, the correlation-induced attraction between condensed counterions and charged monomers and electrostatic repulsion between uncompensated charges provide significant contribution to optimization of the necklace structure. Upon forming a complex with both random and diblock polyampholytes, a polyelectrolyte chain changes its necklace conformation by forming one huge bead. The collapse of the polyelectrolyte chain occurs due to the neutralization of the polyelectrolyte charge by polyampholytes. In the case of the random polyampholyte, the more positively charged sections of the chain mix with negatively charged polyelectrolyte forming the globular bead while more negatively charged chain sections form loops surrounding the collapsed core of the aggregate. In the case of diblock polyampholyte, the positively charged block, a part of the negatively charged block, and a polyelectrolyte chain form a core of the aggregate with a substantial section of the negatively charged block sticking out from the collapsed core of the aggregate. In both cases the core of the aggregate has a layered structure that is characterized by the variations in the excess of concentration of monomers belonging to polyampholyte and polyelectrolyte chains throughout the core radius. These structures appear as a result of optimization of the net electrostatic energy of the complex and short-range attractive interactions between monomers of the polyelectrolyte chain.     

Using Markov models to simulate electron spin resonance spectra from molecular dynamics trajectories

Simulating electron spin resonance (ESR) spectra directly from molecular dynamics simulations of a spin-labeled protein necessitates a large number (hundreds or thousands) of relatively long (hundreds of nanoseconds) trajectories. To meet this challenge, we explore the possibility of constructing accurate stochastic models of the spin label dynamics from atomistic trajectories. A systematic, two-step procedure, based on the probabilistic framework of hidden Markov models, is developed to build a discrete-time Markov chain process that faithfully captures the internal spin label dynamics on time scales longer than about 150 ps. The constructed Markov model is used both to gain insight into the long-lived conformations of the spin label and to generate the stochastic trajectories required for the simulation of ESR spectra. The methodology is illustrated with an application to the case of a spin-labeled poly alanine alpha helix in explicit solvent.     

高分子动力学的单链模型

高分子单链模型是高分子稀溶液理论研究的基本模型.对其进行深入地分析,不仅有助于解决高分子稀溶液体系中溶液黏度和分子链扩散等基本问题,而且能够增进人们对高分子链结构与溶液性质间关联性的理解.虽然基于经典连续性介质力学的流体动力学理论可以定性,甚至半定量地获得稀溶液的一些重要性质,但是,随着科学技术的发展,人们从分子水平上建立了许多描述高分子稀溶液性质的模型和理论,期望能够定量地描述高分子稀溶液的性质.本文以高分子稀溶液中3个典型的单链模型为例(包括:不含流体力学相互作用的Rouse模型、含二体流体力学相互作用的Zimm模型和含多体流体力学相互作用的部分穿透球模型),综述高分子稀溶液的重要性质,并详细地给出其动力学方程的推导过程及其重要的研究进展.特别是,对于Rouse模型,本文还将其预言结果拓展到了短链高分子流体体系;此外,还介绍了这一领域的关键科学问题、发展前景和研究方向.     

Foldamer dynamics expressed via Markov state models. I. Explicit solvent molecular-dynamics simulations in acetonitrile, chloroform, methanol, and water

In this article, we analyze the folding dynamics of an all-atom model of a polyphenylacetylene (pPA) 12-mer in explicit solvent for four common organic and aqueous solvents: acetonitrile, chloroform, methanol, and water. The solvent quality has a dramatic effect on the time scales in which pPA 12-mers fold. Acetonitrile was found to manifest ideal folding conditions as suggested by optimal folding times on the order of approximately 100-200 ns, depending on temperature. In contrast, chloroform and water were observed to hinder the folding of the pPA 12-mer due to extreme solvation conditions relative to acetonitrile; chloroform denatures the oligomer, whereas water promotes aggregation and traps. The pPA 12-mer in a pure methanol solution folded in approximately 400 ns at 300 K, compared relative to the experimental 12-mer folding time of approximately 160 ns measured in a 1:1 v/v THF/methanol solution. Requisite in drawing the aforementioned conclusions, analysis techniques based on Markov state models are applied to multiple short independent trajectories to extrapolate the long-time scale dynamics of the 12-mer in each respective solvent. We review the theory of Markov chains and derive a method to impose detailed balance on a transition-probability matrix computed from simulation data.