Entropy and barrier-controlled fluctuations determine conformational viscoelasticity of single biomolecules

Effects of Proline on Internal Friction in Simulated Folding Dynamics of Several Alanine-Based α-Helical Peptides

We have studied in silico the effect of proline, a model cosolvent, on local and global friction coefficients in (un)folding of several typical alanine-based α-helical peptides. Local friction is related to dwell times of a single, ensemble-averaged hydrogen bond (HB) within each peptide. Global friction is related to energy dissipated in a series of configurational changes of each peptide experienced by increasing the number of HBs during folding. Both of these approaches are important in relation to future atomic force microscopic-based measurements of internal friction via force-clamp single-molecule force spectroscopy. Molecular dynamics (MD) simulations for six peptides, namely, ALA5, ALA8, ALA15, ALA21, (AAQAA)3, and H2N-GN(AAQAA)2G-COONH2, have been conducted at 2 and 5 M proline solutions in water. Using previously obtained MD data for these peptides in pure water as well as upgraded theoretical models, we obtained variations of local and global internal friction coefficients as a function of solution viscosity. The results showed the substantial role of proline in stabilizing the folded state and slowing the overall folding dynamics. Consequently, larger friction coefficients were obtained at larger viscosities. The local and global internal friction, i.e., respective, friction coefficients approximated to zero viscosity, was also obtained. The evolution of friction coefficients with viscosity was weakly dependent on the number of concurrent folding pathways but was rather dominated by a stabilizing effect of proline on the folded states. Obtained values of local and global internal friction showed qualitatively similar results and a clear dependency on the structure of the studied peptide.

Analytical Approaches for Deriving Friction Coefficients for Selected α-Helical Peptides Based Entirely on Molecular Dynamics Simulations

In this paper we derive analytically from molecular dynamics (MD) simulations the friction coefficients related to conformational transitions within several model peptides with α-helical structures. We study a series of alanine peptides with various length from ALA₅ to ALA₂₁ as well as their two derivatives, the (AAQAA)₃ peptide and a 13-residue KR1 peptide that is a derivative of the (AAQAA)₂ peptide with the formula GN(AAQAA)₂G. We use two kinds of approaches to derive their friction coefficients. In the local approach, friction associated with fluctuations of single hydrogen bonds are studied. In the second approach, friction coefficients associated with a folding transitions within the studied peptides are obtained. In both cases, the respective friction coefficients differentiated very well the subtle structural changes between studied peptides and compared favorably to experimentally available data.

How important are fluctuations in the treatment of internal friction in polymers?

The Rouse model with internal friction (RIF), a widely used theoretical framework to interpret the effects of internal friction on conformational transitions in biomolecules, is shown to be an approximate treatment that is based on preaveraging internal friction. By comparison with Brownian dynamics simulations of an exact coarse-grained model that incorporates fluctuations in internal friction, the accuracy of the preaveraged model predictions is examined both at and away from equilibrium. While the two models predict intrachain autocorrelations that approach each other for long enough chain segments, they differ in their predictions for shorter segments. Furthermore, the two models differ qualitatively in their predictions for the chain extension and viscosity in shear flow, which is taken to represent a prototypical out-of-equilibrium condition.

Validity of point-mass model in off-resonance dynamic atomic force microscopy

The quantitative measurement of viscoelasticity of nano-scale entities is an important goal of nanotechnology research and there is considerable progress with advent of dynamic atomic force microscopy. The hydrodynamics of cantilever, the force sensor in AFM measurements, plays a pivotal role in quantitative estimates of nano-scale viscoelasticity. The point-mass (PM) model, wherein the AFM cantilever is approximated as a point-mass with mass-less spring is widely used in dynamic AFM analysis and its validity, particularly in liquid environments, is debated. It is suggested that the cantilever must be treated as a continuous rectangular beam to obtain accurate estimates of nano-scale viscoelasticity of materials it is probing. Here, we derived equations, which relate stiffness and damping coefficient of the material under investigation to measured parameters, by approximating cantilever as a point-mass and also considering the full geometric details. These equations are derived for both tip-excited as well as base-excited cantilevers. We have performed off-resonance dynamic atomic force spectroscopy on a single protein molecule to investigate the validity of widely used PM model. We performed measurements with AFMs equipped with different cantilever excitation methods as well as detection schemes to measure cantilever response. The data was analyzed using both, continuous beam model and the PM model. We found that both models yield same results when the experiments are performed in truly off-resonance regime with small amplitudes and the cantilever stiffness is much higher than the interaction stiffness. Our findings suggest that a simple PM approximation based model is adequate to describe the dynamics, provided care is taken while performing experiments so that the approximations used in these models are valid.

The nano-scale viscoelasticity using atomic force microscopy in liquid environment

We measured viscoelasticity of two nanoscale systems, single protein molecules and molecular layers of water confined between solid walls. In order to quantify the viscoelastic response of these nanoscale systems in liquid environment, the measurements are performed using two types of atomic force microscopes (AFMs), which employ different detection schemes to measure the cantilever response. We used a deflection detection scheme, available in commercial AFMs, that measures cantilever bending and a fibre-interferometer based detection which measures cantilever displacement. The hydrodynamics of the cantilever is modelled using Euler-Bernoulli equation with appropriate boundary conditions which accommodate both detection schemes. In a direct contradiction with many reports in the literature, the dissipation coefficient of a single octomer of titin I27₈ is found to be immeasurably low. The upper bound on the dissipation coefficient is 5 × 10^-7 kg s^-1, which is much lower than the reported values. The entropic stiffness of single unfolded domains of protein measured using both methods is in the range of 10 mN m^-1. We show that in a conventional deflection detection measurement, the phase of the bending signal can be a primary source of artefacts in the dissipation estimates. It is recognized that the measurement of cantilever displacement, which has negligibly small phase lag due to hydrodynamics of the cantilever at low excitation frequencies, is better suited for ensuring artefact-free measurement of viscoelasticity compared to the measurement of the cantilever bending. Further, it was possible to measure dissipation in molecular layers of water confined between the tip and the substrate using fibre interferometer based AFM with similar experimental parameters. It confirms that the dissipation coefficient of a single I27₈ is below the detection limit of AFM. The results shed light on the discrepancy observed in the measured diffusional dynamics of protein collapse measured using Force spectroscopic techniques and single-molecule optical techniques.

The non-classical kinetics and the mutual information of polymer loop formation

The loop formation of a single polymer chain has served as a model system for various biological and chemical processes. Theories based on the Smoluchowski equation proposed that the rate constant (k_loop) of the loop formation would be inversely proportional to viscosity (η), i.e., k_loop ∼ η^-1. Experiments and simulations showed, however, that k_loop showed the fractional viscosity dependence of k_loop ∼ η^-β with β < 1 either in glasses or in low-viscosity solutions. The origin of the fractional viscosity dependence remains elusive and has been attributed to phenomenological aspects. In this paper, we illustrate that the well-known failure of classical kinetics of the loop formation results from the breakdown of the local thermal equilibrium (LTE) approximation and that the mutual information can quantify the breakdown of the LTE successfully.

Pectin biopolymer mechanics and microstructure associated with polysaccharide phase transitions

Polysaccharide polymers like pectin can demonstrate striking and reversible changes in their physical properties depending upon relatively small changes in water content. Recent interest in using pectin polysaccharides as mesothelial sealants suggests that water content, rather than nonphysiologic changes in temperature, may be a practical approach to optimize the physical properties of the pectin biopolymers. Here, we used humidified environments to manipulate the water content of dispersed solution of pectins with a high degree of methyl esterification (high-methoxyl pectin; HMP). The gel phase transition was identified by a nonlinear increase in compression resistance at a water content of 50% (w/w). The gel phase was associated with a punched-out fracture pattern and scanning electron microscopy (SEM) images that revealed a cribiform (Swiss cheese-like) pectin microstructure. The glass phase transition was identified by a marked increase in resilience and stiffness. The glass phase was associated with a star-burst fracture pattern and SEM images that demonstrated a homogeneous pectin microstructure. In contrast, the burst strength of the pectin films was largely independent of water content over a range from 5 to 30% (w/w). These observations indicate the potential to use water content in the selective regulation of the physical properties of HMP biopolymers.

Conformational entropy of a single peptide controlled under force governs protease recognition and catalysis

An immense repertoire of protein chemical modifications catalyzed by enzymes is available as proteomics data. Quantifying the impact of the conformational dynamics of the modified peptide remains challenging to understand the decisive kinetics and amino acid sequence specificity of these enzymatic reactions in vivo, because the target peptide must be disordered to accommodate the specific enzyme-binding site. Here, we were able to control the conformation of a single-molecule peptide chain by applying mechanical force to activate and monitor its specific cleavage by a model protease. We found that the conformational entropy impacts the reaction in two distinct ways. First, the flexibility and accessibility of the substrate peptide greatly increase upon mechanical unfolding. Second, the conformational sampling of the disordered peptide drives the specific recognition, revealing force-dependent reaction kinetics. These results support a mechanism of peptide recognition based on conformational selection from an ensemble that we were able to quantify with a torsional free-energy model. Our approach can be used to predict how entropy affects site-specific modifications of proteins and prompts conformational and mechanical selectivity.

Unwinding and rewinding the nucleosome inner turn: force dependence of the kinetic rate constants

A simple model for the force-dependent unwinding and rewinding rates of the nucleosome inner turn is constructed and quantitatively compared to the results of recent measurements [A. H. Mack et al., J. Mol. Biol. 423, 687 (2012)]. First, a coarse-grained model for the histone-DNA free-energy landscape that incorporates both an elastic free-energy barrier and specific histone-DNA bonds is developed. Next, a theoretical expression for the rate of transitions across a piecewise linear free-energy landscape with multiple minima and maxima is presented. Then, the model free-energy landscape, approximated as a piecewise linear function, and the theoretical expression for the transition rates are combined to construct a model for the force-dependent unwinding and rewinding rates of the nucleosome inner turn. Least-mean-squares fitting of the model rates to the rates observed in recent experiments rates demonstrates that this model is able to well describe the force-dependent unwinding and rewinding rates of the nucleosome inner turn, observed in the recent experiments, except at the highest forces studied, where an additional ad hoc term is required to describe the data, which may be interpreted as an indication of an alternate high-force nucleosome disassembly pathway, that bypasses simple unwinding. The good agreement between the measurements and the model at lower forces demonstrates that both specific histone-DNA contacts and an elastic free-energy barrier play essential roles for nucleosome winding and unwinding, and quantifies their relative contributions.

Peptide chain dynamics in light and heavy water: zooming in on internal friction

Frictional effects due to the chain itself, rather than the solvent, may have a significant effect on protein dynamics. Experimentally, such "internal friction" has been investigated by studying folding or binding kinetics at varying solvent viscosity; however, the molecular origin of these effects is hard to pinpoint. We consider the kinetics of disordered glycine-serine and α-helix forming alanine peptides and a coarse-grained protein folding model in explicit-solvent molecular dynamics simulations. By varying the solvent mass over more than two orders of magnitude, we alter only the solvent viscosity and not the folding free energy. Folding dynamics at the near-vanishing solvent viscosities accessible by this approach suggests that solvent and internal friction effects are intrinsically entangled. This finding is rationalized by calculation of the polymer end-to-end distance dynamics from a Rouse model that includes internal friction. An analysis of the friction profile along different reaction coordinates, extracted from the simulation data, demonstrates that internal as well as solvent friction varies substantially along the folding pathways and furthermore suggests a connection between friction and the formation of hydrogen bonds upon folding.

Marine polysaccharides in microencapsulation and application to aquaculture: "from sea to sea"

This review's main objective is to discuss some physico-chemical features of polysaccharides as intrinsic determinants for the supramolecular structures that can efficiently provide encapsulation of drugs and other biological entities. Thus, the general characteristics of some basic polysaccharides are outlined in terms of their conformational, dynamic and thermodynamic properties. The analysis of some polysaccharide gelling properties is also provided, including the peculiarity of the charged polysaccharides. Then, the way the basic physical chemistry of polymer self-assembly is made in practice through the laboratory methods is highlighted. A description of the several literature procedures used to influence molecular interactions into the macroscopic goal of the encapsulation is given with an attempt at classification. Finally, a practical case study of specific interest, the use of marine polysaccharide matrices for encapsulation of vaccines in aquaculture, is reported.

Conformational dynamics and internal friction in homopolymer globules: equilibrium vs. non-equilibrium simulations

We study the conformational dynamics within homopolymer globules by solvent-implicit Brownian dynamics simulations. A strong dependence of the internal chain dynamics on the Lennard-Jones cohesion strength ε and the globule size N (G) is observed. We find two distinct dynamical regimes: a liquid-like regime (for ε < ε(s) with fast internal dynamics and a solid-like regime (for ε > ε(s) with slow internal dynamics. The cohesion strength ε(s) of this freezing transition depends on N (G) . Equilibrium simulations, where we investigate the diffusional chain dynamics within the globule, are compared with non-equilibrium simulations, where we unfold the globule by pulling the chain ends with prescribed velocity (encompassing low enough velocities so that the linear-response, viscous regime is reached). From both simulation protocols we derive the internal viscosity within the globule. In the liquid-like regime the internal friction increases continuously with ε and scales extensive in N (G) . This suggests an internal friction scenario where the entire chain (or an extensive fraction thereof) takes part in conformational reorganization of the globular structure.

Deconvolution of dynamic mechanical networks

Time-resolved single-molecule biophysical experiments yield data that contain a wealth of dynamic information, in addition to the equilibrium distributions derived from histograms of the time series. In typical force spectroscopic setups the molecule is connected via linkers to a readout device, forming a mechanically coupled dynamic network. Deconvolution of equilibrium distributions, filtering out the influence of the linkers, is a straightforward and common practice. We have developed an analogous dynamic deconvolution theory for the more challenging task of extracting kinetic properties of individual components in networks of arbitrary complexity and topology. Our method determines the intrinsic linear response functions of a given object in the network, describing the power spectrum of conformational fluctuations. The practicality of our approach is demonstrated for the particular case of a protein linked via DNA handles to two optically trapped beads at constant stretching force, which we mimic through Brownian dynamics simulations. Each well in the protein free energy landscape (corresponding to folded, unfolded, or possibly intermediate states) will have its own characteristic equilibrium fluctuations. The associated linear response function is rich in physical content, because it depends both on the shape of the well and its diffusivity-a measure of the internal friction arising from such processes as the transient breaking and reformation of bonds in the protein structure. Starting from the autocorrelation functions of the equilibrium bead fluctuations measured in this force clamp setup, we show how an experimentalist can accurately extract the state-dependent protein diffusivity using a straightforward two-step procedure.

Dynamics of the coiled-coil unfolding transition of myosin rod probed by dissipation force spectrum

The motor protein myosin II plays a crucial role in muscle contraction. The mechanical properties of its coiled-coil region, the myosin rod, are important for effective force transduction during muscle function. Previous studies have investigated the static elastic response of the myosin rod. However, analogous to the study of macroscopic complex fluids, how myosin will respond to physiological time-dependent loads can only be understood from its viscoelastic response. Here, we apply atomic force microscopy using a magnetically driven oscillating cantilever to measure the dissipative properties of single myosin rods that provide unique dynamical information about the coiled-coil structure as a function of force. We find that the friction constant of the single myosin rod has a highly nontrivial variation with force; in particular, the single-molecule friction constant is reduced dramatically and increases again as it passes through the coiled-uncoiled transition. This is a direct indication of a large free-energy barrier to uncoiling, which may be related to a fine-tuned dynamic mechanosignaling response to large and unexpected physiological loads. Further, from the critical force at which the minimum in friction occurs we determine the asymmetry of the bistable landscape that controls uncoiling of the coiled coil. This work highlights the sensitivity of the dissipative signal in force unfolding to dynamic molecular structure that is hidden to the elastic signal.

Linear viscoelasticity of a single semiflexible polymer with internal friction

The linear viscoelastic behaviors of single semiflexible chains with internal friction are studied based on the wormlike-chain model. It is shown that the frequency dependence of the complex compliance in the high frequency limit is the same as that of the Voigt model. This asymptotic behavior appears also for the Rouse model with internal friction. We derive the characteristic times for both the high frequency limit and the low frequency limit and compare the results with those obtained by Khatri et al.

The elastic properties of the Cryptococcus neoformans capsule

Microbial capsules are important for virulence, but their architecture and physical properties are poorly understood. The human pathogenic fungus Cryptococcus neoformans has a large polysaccharide capsule that is necessary for virulence and is the target of protective antibody responses. To study the C. neoformans capsule we developed what we believe is a new approach whereby we probed the capsular elastic properties by applying forces using polystyrene beads manipulated with optical tweezers. This method allowed us to determine the Young's modulus for the capsule in various conditions that affect capsule growth. The results indicate that the Young's modulus of the capsule decreases with its size and increases with the Ca(2+) concentration in solution. Also, capsular polysaccharide manifests an unexpected affinity for polystyrene beads, a property that may function in attachment to host cells and environmental structures. Bead probing with optical tweezers provides a new, nondestructive method that may have wide applicability for studying the effects of growth conditions, immune components, and drugs on capsular properties.

Internal friction and nonequilibrium unfolding of polymeric globules

The stretching response of a single collapsed homopolymer is studied using Brownian dynamic simulations. The irreversibly dissipated work is found to be dominated by internal friction effects below the collapse temperature, and the internal viscosity grows exponentially with the effective cohesive strength between monomers. These results explain friction effects of globular DNA and are relevant for dissipation at intermediate stages of protein folding.

Internal friction of single polypeptide chains at high stretch

Experiments that measure the viscoelasticity of single molecules from the Brownian fluctuations of an atomic force microscope (AFM) have provided a new window onto their internal dynamics in an underlying conformational landscape. Here we develop and apply these methods to examine the internal friction of unfolded polypeptide chains at high stretch. The results reveal a power law dependence of internal friction with tension (exponent 1.3 +/- 0.5) and a relaxation time approximately independent of force. To explain these results we develop a frictional worm-like chain (FWLC) model based on the Rayleigh dissipation function of a stiff chain with dynamical resistance to local bending. We analyse the dissipation rate integrated over the chain length by its Fourier components to calculate an effective tension-dependent friction constant for the end-to-end vector of the chain. The result is an internal friction that increases as a power law with tension with an exponent 3/2, consistent with experiment. Extracting the intrinsic bending friction constant of the chain it is found to be approximately 7 orders of magnitude greater than expected from solvent friction alone; a possible explanation we offer is that the underlying energy landscape for bending amino acids and/or peptide bond is rough, consistent with recent results on both proteins and polysaccharides.

The effect of temperature on mechanical resistance of the native and intermediate states of I27

We investigated the effect of temperature on the mechanical unfolding of I27 from human cardiac titin, employing a custom-built temperature control device for single-molecule atomic force microscopy measurement. A sawtooth pattern was observed in the force curves where each force peak reports on the unfolding of an I27 domain. In early unfolding events, we observed a "hump-like" deviation due to the detachment of beta-strand A from each I27 domain. The force at which the humps appear was approximately 130 pN and showed no temperature dependence, at least in the temperature range of 2 degrees C-30 degrees C. The hump structure was successfully analyzed with a two-state worm-like chain model, and the Gibbs free energy difference of the detachment reaction was estimated to be 11.6 +/- 0.58 kcal/mol and found to be temperature independent. By contrast, upon lowering the temperature, the mean unfolding force from the partly unfolded intermediate state was found to markedly increase and the unfolding force distribution to broaden significantly, suggesting that the distance (x(u)) between the folded and transition states in the energy landscape along the pulling direction is decreased. These results suggest that the local structure of beta-strand A are stabilized by topologically simple local hydrogen-bond network and that the temperature does not affect the detachment reaction thermodynamically and kinetically, whereas the interaction between the beta-strands A' and G, which is a critical region for its mechanical stability, is strongly dependent on the temperature.

Dynamics of single-molecule force-ramp experiments: The role of fluctuations

In the force-ramp mode of the atomic force microscope, the force with which a macromolecule is stretched is increased linearly in time by properly controlling the motion of the cantilever through a feedback loop. Using a master equation approach for the coupled cantilever-macromolecule system, we mimic such a feedback loop, to study nonequilibrium effects in the measurements of force-extension curves and fluctuations. In particular, it is shown that the fluctuations are the same for force-ramp experiments and for the more commonly used constant velocity experiments. Thus the exact same statistics suffice for the explanation of either experiment. Specific results are presented for the stretching of Dextran.

Investigation of the effects of fine structure on the nanomechanical properties of pectin

Pectin is an important structural polysaccharide found in the cell walls of all land plants. While in detail its composition and its organization in muro are complex, it is predominantly a copolymer of galacturonic acid and its methylesterified counterpart. Previous single-molecule stretching studies carried out on a sparsely methylesterified pectin sample indicated the importance of force-induced conformational transitions of the pyranose ring during extension, and the possible biological role of such transitions was discussed. More heavily methylesterified samples are better biomimetic models of the polymeric components as found in the plant cell wall, in particular being less restricted by the shackles of the significant intermolecular interactions expected to constrain the behavior of bare galacturonic acid sequences. Density functional theory calculations revealed that upon extending galacturonic acid monomers, whether methylesterified or not, the initial ((4)C1) chair structure is transformed to a ((3)S5) skew boat and that subsequently upon further elongation, via an intermediate inverted skew boat ((5)S3), the inverted chair ((1)C4) is reached. Experimentally, the force-extension curve of highly methylesterified pectin was found to be solvent dependent in the same manner as the un-esterified sample, indicating that minimal changes in the strength of interring hydrogen bonding result from such a substitution, and finally, as only subtle changes in the force-extension behavior of pectin resulted from changes in the degree of methylesterification, previous speculations about the role of force-induced transformations in vivo are supported.