Kinetics of the polymer collapse transition: the role of hydrodynamics

Collapse and expansion kinetics of a single polyelectrolyte chain with hydrodynamic interactions

We investigate the collapse and expansion dynamics of a linear polyelectrolyte (PE) with hydrodynamic interactions. Using dissipative particle dynamics with a bead-spring PE model, long-range electrostatics, and explicit ions, we examine how the timescales of collapse tcol and expansion texp depend on the chain length N and obtain scaling relationships tcol ∼ Nα and texp ∼ Nβ. For neutral polymers, we derive values of α = 0.94 ± 0.01 and β = 1.97 ± 0.10. Interestingly, the introduction of electrostatic interaction markedly shifts α to α ≈ 1.4 ± 0.1 for salt concentrations within c = 10-4 to 10-2 M. A reduction in the ion-to-monomer size ratio noticeably reduces α. On the other hand, the expansion scaling remains approximately constant, β ≈ 2, regardless of the salt concentration or ion size considered. We find β > α for all conditions considered, implying that expansion is always slower than collapse in the limit of long polymers. This asymmetry is explained by distinct kinetic pathways of collapse and expansion processes.

Impact of Hydrodynamic Interactions on the Kinetic Pathway of Protein Folding

Protein folding is a fundamental process critical to cellular function and human health, but it remains a grand challenge in biophysics. Hydrodynamic interaction (HI) plays a vital role in the self-organization of soft and biological materials, yet its role in protein folding is not fully understood despite folding occurring in a fluid environment. Here, we use the fluid particle dynamics method to investigate many-body hydrodynamic couplings between amino acid residues and fluid motion in the folding kinetics of a coarse-grained four-α-helices bundle protein. Our results reveal that HI helps select fast folding pathways to the native state without being kinetically trapped, significantly speeding up the folding kinetics compared to its absence. First, the directional flow along the protein backbone expedites protein collapse. Then, the incompressibility-induced squeezing flow effects retard the accumulation of non-native hydrophobic contacts, thus preventing the protein from being trapped in local energy minima during the conformational search of the native structure. We also find that the significance of HI in folding kinetics depends on temperature, with a pronounced effect under biologically relevant conditions. Our findings suggest that HI, particularly the short-range squeezing effect, may be crucial in avoiding protein misfolding.

The role of flow in the self-assembly of dragline spider silk proteins

Hydrodynamic flow in the spider duct induces conformational changes in dragline spider silk proteins (spidroins) and drives their assembly, but the underlying physical mechanisms are still elusive. Here we address this challenging multiscale problem with a complementary strategy of atomistic and coarse-grained molecular dynamics simulations with uniform flow. The conformational changes at the molecular level were analyzed for single-tethered spider silk peptides. Uniform flow leads to coiled-to-stretch transitions and pushes alanine residues into β sheet and poly-proline II conformations. Coarse-grained simulations of the assembly process of multiple semi-flexible block copolymers using multi-particle collision dynamics reveal that the spidroins aggregate faster but into low-order assemblies when they are less extended. At medium-to-large peptide extensions (50%-80%), assembly slows down and becomes reversible with frequent association and dissociation events, whereas spidroin alignment increases and alanine repeats form ordered regions. Our work highlights the role of flow in guiding silk self-assembly into tough fibers by enhancing alignment and kinetic reversibility, a mechanism likely relevant also for other proteins whose function depends on hydrodynamic flow.

Activity mediated globule to coil transition of a flexible polymer in a poor solvent

Understanding the role of self-propulsion on the conformational properties of active filamentous objects has relevance in biology. In this work, we consider a flexible bead-spring model for active polymers with both attractive and repulsive interactions among the non-bonded monomers. The activity for each monomer works along its intrinsic direction of self-propulsion which changes diffusively with time. We study its kinetics in the overdamped limit, following quenching from good to poor solvent conditions. We observe that with low activities, though the kinetic pathways remain similar, the scaling exponent for the relaxation time of globule formation becomes smaller than that for the case with no activity. Interestingly, for higher activities when self-propulsion dominates over interaction energy, the polymer conformation becomes extended coil-like. There, in the steady state, the variation of the spatial extension of the polymer, measured via its gyration radius, shows two completely different scaling regimes: the corresponding Flory exponent ν changes from 1/3 to 3/5 similar to a transition of the polymer from a globular state to a self-avoiding walk. This can be explained by an interplay among the three energy scales present in the system, viz., the "ballistic", thermal, and interaction energy.

Delayed collapse transitions in a pinned polymer system

Employing Langevin dynamics simulations, we investigated the kinetics of the collapse transition for a polymer of length N when a particular monomer at a position 1≤X≤N is pinned. The results are compared with the kinetics of a free polymer. The equilibrium θ-point separating the coil from the globule phase is located by a crossover in 〈R_{g}^{2}〉/N plots of different chain lengths. Our simulation supports a three-stage mechanism for free and pinned polymer collapse: the formation of pearls, the coarsening of pearls, and the formation of a compact globule. Pinning the central monomer has negligible effects on the kinetics as it does not break the symmetry. However, pinning a monomer elsewhere causes the process to be delayed by a constant factor ϕ_{X} depending linearly upon X. The total collapse time scales with N as τ_{c}∼ϕ_{X}N^{1.60±0.03}, which implies τ_{c} is maximum when an end monomer is pinned (X=1 or N), while when pinning the central monomer (X=N/2) it is minimum and identical to that of a free polymer. The average cluster size N_{c}(t) grows in time as t^{z}, where z=1.00±0.04 for a free particle, whereas we identify two time regimes separated by a plateau for pinned polymers. At longer times, z=1.00±0.04, while it deviates in early time regimes significantly, depending on the value of X.

Kinetics of charged polymer collapse in poor solvents

Extensive molecular dynamics simulations, using simple charged polymer models, have been employed to probe the collapse kinetics of a single flexible polyelectrolyte (PE) chain under implicit poor solvent conditions. We investigate the role of the charged nature of PE chain (A), valency of counterions (Z) on the kinetics of such PE collapse. Our study shows that the collapse kinetics of charged polymers are significantly different from those of the neutral polymer and that the finite-size scaling behavior of PE collapse times does not follow the Rouse scaling as observed in the case of neutral polymers. The critical exponent for charged PE chains is found to be less than that of neutral polymers and also exhibits dependence on counterion valency. The coarsening of clusters along the PE chain suggests a multi-stage collapse and exhibits opposite behavior of exponents compared to neutral polymers: faster in the early stages and slower in the later stages of collapse.

Mesoscale modelling of polymer aggregate digestion

We use mesoscale simulations to gain insight into the digestion of biopolymers by studying the break-up dynamics of polymer aggregates (boluses) bound by physical cross-links. We investigate aggregate evolution, establishing that the linking bead fraction and the interaction energy are the main parameters controlling stability with respect to diffusion. We show via a simplified model that chemical breakdown of the constituent molecules causes aggregates that would otherwise be stable to disperse. We further investigate breakdown of biopolymer aggregates in the presence of fluid flow. Shear flow in the absence of chemical breakdown induces three different regimes depending on the flow Weissenberg number ( W i ). i) At W i ≪ 1 , shear flow has a negligible effect on the aggregates. ii) At W i ∼ 1 , the aggregates behave approximately as solid bodies and move and rotate with the flow. iii) At W i ≫ 1 , the energy input due to shear overcomes the attractive cross-linking interactions and the boluses are broken up. Finally, we study bolus evolution under the combined action of shear flow and chemical breakdown, demonstrating a synergistic effect between the two at high reaction rates.

Kinetic mechanisms of crumpled globule formation

A homopolymer chain with beads forming pairwise reversible bonds is a well-known model in polymer physics. We studied the kinetics of homopolymer chain collapse, which was induced by pairwise reversible bond formation. We compared the kinetic mechanism of this coil-globule transition with the mechanism of collapse in a poor solvent. We discovered that if collapse of a spatially confined chain is induced by formation of long-living pairwise reversible bonds, the transient structures do not resemble classical pearl-necklace conformations formed during collapse in a poor solvent. However, both types of collapse lead to formation of a metastable state of a crumpled globule, which is one of the well-known models of interphase chromatin structure in different organisms. Moreover, we found that the stability and dynamics of this state can be controlled by the fraction of reversible bonds and bond lifetime.

Pearl-Necklace-Like Local Ordering Drives Polypeptide Collapse

The collapse of the polypeptide backbone is an integral part of protein folding. Using polyglycine as a probe, we explore the nonequilibrium pathways of protein collapse in water. We find that the collapse depends on the competition between hydration effects and intrapeptide interactions. Once intrapeptide van der Waal interactions dominate, the chain collapses along a nonequilibrium pathway characterized by formation of pearl-necklace-like local clusters as intermediates that eventually coagulate into a single globule. By describing this coarsening through the contact probability as a function of distance along the chain, we extract a time-dependent length scale that grows in a linear fashion. The collapse dynamics is characterized by a dynamical critical exponent z ≈ 0.5 that is much smaller than the values of z = 1–2 reported for nonbiological polymers. This difference in the exponents is explained by the instantaneous formation of intrachain hydrogen bonds and local ordering that may be correlated with the observed fast folding times of proteins.

A hybrid method for micro-mesoscopic stochastic simulation of reaction-diffusion systems

The present paper introduces a new micro-meso hybrid algorithm based on the Ghost Cell Method concept in which the microscopic subdomain is governed by the Reactive Multi-Particle Collision (RMPC) dynamics. The mesoscopic subdomain is modeled using the Reaction-Diffusion Master Equation (RDME). The RDME is solved by means of the Inhomogeneous Stochastic Simulation Algorithm. No hybrid algorithm has hitherto used the RMPC dynamics for modeling reactions and the trajectories of each individual particle. The RMPC is faster than other molecular based methods and has the advantage of conserving mass, energy and momentum in the collision and free streaming steps. The new algorithm is tested on three reaction-diffusion systems. In all the systems studied, very good agreement with the deterministic solutions of the corresponding differential equations is obtained. In addition, it has been shown that proper discretization of the computational domain results in significant speed-ups in comparison with the full RMPC algorithm.

Impact of hydrodynamic interactions on protein folding rates depends on temperature

We investigated the impact of hydrodynamic interactions (HI) on protein folding using a coarse-grained model. The extent of the impact of hydrodynamic interactions, whether it accelerates, retards, or has no effect on protein folding, has been controversial. Together with a theoretical framework of the energy landscape theory (ELT) for protein folding that describes the dynamics of the collective motion with a single reaction coordinate across a folding barrier, we compared the kinetic effects of HI on the folding rates of two protein models that use a chain of single beads with distinctive topologies: a 64-residue α/β chymotrypsin inhibitor 2 (CI2) protein, and a 57-residue β-barrel α-spectrin Src-homology 3 domain (SH3) protein. When comparing the protein folding kinetics simulated with Brownian dynamics in the presence of HI to that in the absence of HI, we find that the effect of HI on protein folding appears to have a "crossover" behavior about the folding temperature. This means that at a temperature greater than the folding temperature, the enhanced friction from the hydrodynamic solvents between the beads in an unfolded configuration results in lowered folding rate; conversely, at a temperature lower than the folding temperature, HI accelerates folding by the backflow of solvent toward the folded configuration of a protein. Additionally, the extent of acceleration depends on the topology of a protein: for a protein like CI2, where its folding nucleus is rather diffuse in a transition state, HI channels the formation of contacts by favoring a major folding pathway in a complex free energy landscape, thus accelerating folding. For a protein like SH3, where its folding nucleus is already specific and less diffuse, HI matters less at a temperature lower than the folding temperature. Our findings provide further theoretical insight to protein folding kinetic experiments and simulations.

Kinetics of polymer collapse: effect of temperature on cluster growth and aging

Using state of the art Monte Carlo simulations of a bead-spring model we investigate both the equilibrium and the nonequilibrium behavior of the homopolymer collapse. The equilibrium properties obtained via multicanonical sampling recover the well-known finite-size scaling behavior of collapse for our model polymer. For the nonequilibrium dynamics we study the collapse by quenching the homopolymer from an expanded coiled state into the globular phase. The sequence of events observed during the collapse is independent of the quench depth. In particular, we focus on finding out universal scaling behaviors related to the growth or coarsening of clusters of monomers, by drawing phenomenological analogies with ordering kinetics. We distinguish the cluster coarsening stage from the initial stage of primary cluster formation. By successful application of a nonequilibrium finite-size scaling analysis we show that at all quench temperatures, during the coarsening stage, the cluster growth is roughly linear and can be characterised by a universal finite-size scaling function. In addition, we provide evidence of aging by constructing a suitable autocorrelation function and its corresponding dynamical power-law scaling with respect to the growing cluster sizes. The predicted theoretical bound for the exponent governing such scaling is strictly obeyed by the numerical data irrespective of the quench temperature. The results and methods presented here in general should find application in similar phenomena such as the collapse of a protein molecule preceding its folding.

Thermoosmotic microfluidics

Microchannels with asymmetrically ratcheted walls are here shown to behave as effective and versatile microfluidic pumps if locally heated. When the boundary walls have different temperatures, the confined liquid experiences a temperature gradient along the sawtooth edges, which can induce a thermoosmotic flow. A mesoscale molecular simulation approach is here employed to investigate the flows which are contrasted using an analytical approach. Microchannels can be composed by one or two ratcheted walls which can be straight or cylindrical. Varying the channel geometry can not only change the overall fluid flux, but also vary the flow patters from shear to capillary type, or even to extensional type flows. This scheme does not require multiphase fluids or any movable channel parts, although they are possible to be implemented. The proposed principle is then very versatile to locally manipulate complex fluids, and a promising tool to recover waste heat, to facilitate cooling of microchips, and to manufacture portable lab-on-a-chip devices.

Molecular dynamics of protein A and a WW domain with a united-residue model including hydrodynamic interaction

The folding of the N-terminal part of the B-domain of staphylococcal protein A (PDB ID: 1BDD, a 46-residue three-α-helix bundle) and the formin-binding protein 28 WW domain (PDB ID: 1E0L, a 37-residue three-stranded anti-parallel β protein) was studied by means of Langevin dynamics with the coarse-grained UNRES force field to assess the influence of hydrodynamic interactions on protein-folding pathways and kinetics. The unfolded, intermediate, and native-like structures were identified by cluster analysis, and multi-exponential functions were fitted to the time dependence of the fractions of native and intermediate structures, respectively, to determine bulk kinetics. It was found that introducing hydrodynamic interactions slows down both the formation of an intermediate state and the transition from the collapsed structures to the final native-like structures by creating multiple kinetic traps. Therefore, introducing hydrodynamic interactions considerably slows the folding, as opposed to the results obtained from earlier studies with the use of Gō-like models.

Electrophoretic Mobility of Polyelectrolytes within a Confining Well

We present a numerical study of polyelectrolytes electrophoresing in free solution while squeezed by an axisymmetric confinement force transverse to their net displacement. Hybrid multiparticle collision dynamics and molecular dynamics simulations with mean-field finite Debye layers show that even though the polyelectrolyte chains remain "free-draining" their electrophoretic mobility increases with confinement in nanoconfining potential wells. The primary mechanism leading to the increase in mobility above the free-solution value, despite long-range hydrodynamic screening by counterion layers, is the orientation of polymer segments within Debye layers. The observed length dependence of the electrophoretic mobility arises due to secondary effects of counterion condensation related to confinement compactification.

Unusual size-dependence of effective interactions between collapsed polymers in crowded environments

We investigate the influence of macromolecular crowding on interactions between collapsed polymers using computer simulations, to gain insights into biomacromolecular interactions in crowded biological environments. The effective attraction is induced between two collapsed polymers due to the macromolecular crowding, and it is found that the strength of the effective attraction decreases as the crowder size is reduced for a fixed crowder volume fraction, which is sharply contrasted with the conventional viewpoint based on the depletion attraction observed for hard-core spherical colloids. This unusual trend of size-dependence is interpreted by dividing the effective interaction into the polymer-mediated repulsion and crowder-mediated attraction. It is found that the ranges of repulsive and attractive contributions overlap significantly due to the flexible nature of polymer boundaries, resulting in partial cancellation over this range which leads to the observed size-dependence. Thus, this work suggests that the effective interactions between biomacromolecules in crowded environments may be qualitatively different from the depletion interactions predicted for hard-core spherical colloids.

Dynamics of solutes with hydrodynamic interactions: comparison between Brownian dynamics and stochastic rotation dynamics simulations

The dynamics of particles in solution or suspension is influenced by thermal fluctuations and hydrodynamic interactions. Several mesoscale methods exist to account for these solvent-induced effects such as Brownian dynamics with hydrodynamic interactions and hybrid molecular dynamics-stochastic rotation dynamics methods. Here we compare two ways of coupling solutes to the solvent with stochastic rotation dynamics (SRD) to Brownian dynamics with and without explicit hydrodynamic interactions. In the first SRD scheme [SRD with collisional coupling (CC)] the solutes participate in the collisional step with the solvent and in the second scheme [SRD with central force coupling (CFC)] the solutes interact through direct forces with the solvent, generating slip boundary conditions. We compare the transport coefficients of neutral and charged solutes in a model system obtained by these simulation schemes. Brownian dynamics without hydrodynamic interactions is used as a reference to quantify the influence of hydrodynamics on the transport coefficients as modeled by the different methods. We show that, in the dilute range, the SRD CFC method provides results similar to those of Brownian dynamics with hydrodynamic interactions for the diffusion coefficients and for the electrical conductivity. The SRD CC scheme predicts diffusion coefficients close to those obtained by Brownian dynamic simulations without hydrodynamic interactions, but accounts for part of the influence of hydrodynamics on the electrical conductivity.

A molecule-centered method for accelerating the calculation of hydrodynamic interactions in Brownian dynamics simulations containing many flexible biomolecules

Inclusion of hydrodynamic interactions (HIs) is essential in simulations of biological macromolecules that treat the solvent implicitly if the macromolecules are to exhibit correct translational and rotational diffusion. The present work describes the development and testing of a simple approach aimed at allowing more rapid computation of HIs in coarse-grained Brownian dynamics simulations of systems that contain large numbers of flexible macromolecules. The method combines a complete treatment of intramolecular HIs with an approximate treatment of the intermolecular HIs which assumes that the molecules are effectively spherical; all of the HIs are calculated at the Rotne-Prager-Yamakawa level of theory. When combined with Fixman's Chebyshev polynomial method for calculating correlated random displacements, the proposed method provides an approach that is simple to program but sufficiently fast that it makes it computationally viable to include HIs in large-scale simulations. Test calculations performed on very coarse-grained models of the pyruvate dehydrogenase (PDH) E2 complex and on oligomers of ParM (ranging in size from 1 to 20 monomers) indicate that the method reproduces the translational diffusion behavior seen in more complete HI simulations surprisingly well; the method performs less well at capturing rotational diffusion but its discrepancies diminish with increasing size of the simulated assembly. Simulations of residue-level models of two tetrameric protein models demonstrate that the method also works well when more structurally detailed models are used in the simulations. Finally, test simulations of systems containing up to 1024 coarse-grained PDH molecules indicate that the proposed method rapidly becomes more efficient than the conventional BD approach in which correlated random displacements are obtained via a Cholesky decomposition of the complete diffusion tensor.

On the importance of hydrodynamic interactions in lipid membrane formation

Hydrodynamic interactions (HI) give rise to collective motions between molecules, which are known to be important in the dynamics of random coil polymers and colloids. However, their role in the biological self-assembly of many molecule systems has not been investigated. Here, using Brownian dynamics simulations, we evaluate the importance of HI on the kinetics of self-assembly of lipid membranes. One-thousand coarse-grained lipid molecules in periodic simulation boxes were allowed to assemble into stable bilayers in the presence and absence of intermolecular HI. Hydrodynamic interactions reduce the monomer-monomer association rate by 50%. In contrast, the rate of association of lipid clusters is much faster in the presence of intermolecular HI. In fact, with intermolecular HI, the membrane self-assembly rate is 3-10 times faster than that without intermolecular HI. We introduce an analytical model to describe the size dependence of the diffusive encounter rate of particle clusters, which can qualitatively explain our simulation results for the early stage of the membrane self-assembly process. These results clearly suggest that HI greatly affects the kinetics of self-assembly and that simulations without HI will significantly underestimate the kinetic parameters of such processes.

High-resolution imaging in a deep turbid medium based on an ultrasound-switchable fluorescence technique

The spatial resolution of fluorescence imaging techniques in deep optically turbid media such as tissues is limited by photon diffusion. To break the diffusion limit and achieve high-resolution and deep-tissue fluorescence imaging, a fundamentally different method was demonstrated based on a concept of ultrasound-switchable fluorescence. The results showed that a small fluorescent tube with a diameter of ∼180 μm at a depth of ∼20 mm in an optical scattering medium ([Formula: see text] and [Formula: see text] cm(-1)) can be clearly imaged with a size of ∼260 μm. The depth-to-resolution ratio is shown to be about one order of magnitude better than other deep-tissue fluorescence imaging techniques.

Collapse kinetics of vibrated granular chains

The kinetics of the collapse of the coil state into condensed states is studied with vibrated granular chain composed of N metal beads partially immersed in water. The radius of gyration of the chain, R(g) is measured. For short chains (N < 140), disk-like condensed state is formed and R(g) decreases with time such that the function ΔR(g)(2) (≡ R(g)(2) - R(g)(2)(∞)) = A e(-t/τ), where the relaxation time τ follows a power-law dependence on the chain length N with an exponent γ = 1.9 ± 0.2. For the chains with length N ≥ 300, rod-like clusters are observed during the initial stage of collapse and R(g)(2) = R(g)(2)(0) - Bt(β), with β = 0.6 ± 0.1. In the coarsening stage, the exponential dependence of ΔR(g)(2) on time still holds, however, the relaxation time τ fluctuates and has no simple dependence on N. Furthermore, the time dependence of the averaged radius of gyration of the individual clusters, R(g,cl) can be described by the theory of Lifshitz and Slyozov. A peak in the structure function of long chains is observed in the initial stage of the collapse transition. The collapse transition in the bead chains is a first order phase transition. However, features of the spinodal decomposition are also observed.

Hydrodynamic fluctuations in thermostatted multiparticle collision dynamics

In this work we study the behavior of mesoscopic fluctuations of a fluid simulated by Multiparticle Collision Dynamics when this is applied together with a local thermostatting procedure that constrains the strength of temperature fluctuations. We consider procedures in which the thermostat interacts with the fluid at every simulation step as well as cases in which the thermostat is applied only at regular time intervals. Due to the application of the thermostat temperature fluctuations are forced to relax to equilibrium faster than they do in the nonthermostatted, constant-energy case. Depending on the interval of application of the thermostat, it is demonstrated that the thermodynamic state changes gradually from isothermal to adiabatic conditions. In order to exhibit this effect we compute from simulations diverse correlation functions of the hydrodynamic fluctuating fields. These correlation functions are compared with those predicted by a linearized hydrodynamic theory of a simple fluid in which a thermostat is applied locally. We find a good agreement between the model and the numerical results, which confirms that hydrodynamic fluctuations in Multiparticle Collision Dynamics in the presence of the thermostat have the properties expected for spontaneous fluctuations in fluids in contact with a heat reservoir.

Effect of confinement on the collapsing mechanism of a flexible polymer chain

In this paper, Brownian dynamics simulation (BDS) studies are executed to demonstrate the distinctive influences of the extent of confinement on the collapsing mechanism and kinetics of a flexible hydrophobic polymer chain in a poor solvent. The collapsing behavior is quantified by the time of collapse, which below a critical dimension of the confinement (h(c)), encounters a drastic reduction with a further strengthening in the degree of confinement. For dimensions greater than this critical one, the collapse occurs through the well-known hydrodynamic interaction (HI) controlled multiple-globule-mediated mechanisms. However, for channel dimensions less than this critical one, the collapse mechanism is drastically altered. Under such circumstances, the collapse gets predominantly controlled by the confinement effects (with negligible contribution of the HIs) and occurs via the formation of a single central globule. This central globule rapidly engulfs the noncondensed polymer segments, and in the process largely hastens up the collapsing event. Under such circumstances, the collapse time is found to decrease linearly with decrements in the channel height. On the contrary, for channel heights greater than h(c), the multiple-globule-mediated collapse is characterized by a collapse time that shows an exponential dependence on the channel height, rapidly attaining a state in which the confinement effect becomes inconsequential and HIs dictate the entire collapsing behavior. We further propose detailed arguments based on physical reasoning as well as free energy estimations to conclusively support the qualitative and quantitative nature of influences of the confinement on the polymer collapse.

Hydrodynamic effects in proteins

Experimental and numerical results pertaining to flow-induced effects in proteins are reviewed. Special emphasis is placed on shear-induced unfolding and on the role of solvent mediated hydrodynamic interactions in the conformational transitions in proteins.

Probing solvation decay length in order to characterize hydrophobicity-induced bead-bead attractive interactions in polymer chains

In this paper, we quantitatively demonstrate that exponentially decaying attractive potentials can effectively mimic strong hydrophobic interactions between monomer units of a polymer chain dissolved in aqueous solvent. Classical approaches to modeling hydrophobic solvation interactions are based on invariant attractive length scales. However, we demonstrate here that the solvation interaction decay length may need to be posed as a function of the relative separation distances and the sizes of the interacting species (or beads or monomers) to replicate the necessary physical interactions. As an illustrative example, we derive a universal scaling relationship for a given solute-solvent combination between the solvation decay length, the bead radius, and the distance between the interacting beads. With our formalism, the hydrophobic component of the net attractive interaction between monomer units can be synergistically accounted for within the unified framework of a simple exponentially decaying potential law, where the characteristic decay length incorporates the distinctive and critical physical features of the underlying interaction. The present formalism, even in a mesoscopic computational framework, is capable of incorporating the essential physics of the appropriate solute-size dependence and solvent-interaction dependence in the hydrophobic force estimation, without explicitly resolving the underlying molecular level details.

Hydrodynamic selection of the kinetic pathway of a polymer coil-globule transition

Recently, the role of hydrodynamic interactions in the selection of a kinetic pathway for phase transitions has attracted considerable attention. Here we study this problem numerically by taking as an example a coil-globule transition of a single polymer, which is a prototype model of protein folding. When a swollen polymer collapses into a globule state, hydrodynamic interactions accelerate the transition. We find, on the other hand, that when a rather compact polymer collapses into the same final state, hydrodynamic interactions decelerate the transition due to a slow squeezing process of the solvent. We reveal that the degree of the initial enhancement of anisotropy of the polymer configuration determines whether hydrodynamic interactions accelerate or decelerate the collapsing dynamics. We also discuss the possible relevance of squeezing flow effects in protein folding.