Internal friction of single polypeptide chains at high stretch

Molecular Origin of Internal Friction in Intrinsically Disordered Proteins

Protein folding and dynamics are controlled by an interplay of thermal and viscosity effects. The effect of viscous drag through the solvent molecules is described by the classic Kramers theory in the high friction limit, which considers the dampening of the reactant molecules in the solution and quantifies the dependence of the reaction rate on the frictional drag. In addition to the external energy dissipation originating from the surrounding solvent molecules, there is an additional mode of internal energy dissipative force operative within the polypeptide chain reflecting the internal resistance of the chain to its conformational alterations. This dry, solvent-independent intrinsic frictional drag is termed internal friction. In the case of natively folded proteins, the physical origin of internal friction is primarily attributed to the intrachain interactions or other nonnative interactions in their compact states. However, the molecular origin of internal friction in intrinsically disordered proteins (IDPs) remains elusive.In this Account, we address this fundamental issue: what are the principal drivers of viscosity-independent (dry) friction in highly solvated, expanded, conformationally flexible, rapidly fluctuating IDPs that do not possess persistent intrachain interactions? IDPs exhibit diffusive conformational dynamics that is predominantly dominated by the segmental motion of the backbone arising due to the dihedral rotations in the Ramachandran Φ-Ψ space. The physical origin of friction in a complex biopolymeric system such as IDPs can be described by classic polymer models, namely, Rouse/Zimm models with internal friction. These one-dimensional models do not invoke torsional fluctuation components. Kuhn's classic description includes the connection between internal friction and microscopic dihedral hopping. Based on our time-resolved fluorescence anisotropy results, we describe that the sequence-dependent, collective, short-range backbone dihedral rotations govern localized internal friction in an archetypal IDP, namely, α-synuclein. The highly sensitive, residue-specific fluorescence depolarization kinetics offers a potent methodology to characterize and quantify the directional decorrelation engendered due to the short-range dihedral relaxation of the polypeptide backbone in the dihedral space. We utilized this characteristic relaxation time scale as our dynamic readout to quantify the site-specific frictional component. Our linear viscosity-dependent model of torsional relaxation time scale furnished a finite nonzero time constant at the zero solvent viscosity representing the solvent-independent internal friction. These results unveil the effect of the degree of dihedral restraining parameter on the internal friction component by showing that a restrained proline residue imparts higher torsional stiffness in the chain segments and, therefore, exhibits higher internal friction. This Account sheds light on the molecular underpinning of the sequence-specific internal friction in IDPs and will be of interest to unmask the role of internal friction in a diverse range of biomolecular processes involving binding-induced folding, allosteric interaction, protein misfolding and aggregation, and biomolecular condensation via phase separation.

Direct and Simultaneous Measurement of the Stiffness and Internal Friction of a Single Folded Protein

The nanomechanical response of a folded single protein, the natural nanomachine responsible for myriad biological processes, provides insight into its function. The conformational flexibility of a folded state, characterized by its viscoelasticity, allows proteins to adopt different shapes to perform their function. Despite efforts, its direct measurement has not been possible so far. We present a direct and simultaneous measurement of the stiffness and internal friction of the folded domains of the protein titin using a special interferometer based atomic force microscope. We analyzed the data by carefully separating different contributions affecting the response of the experimental probe to obtain the folded state's viscoelasticity. Above ∼95 pN of force, the individual immunoglobulins of titin transition from an elastic solid-like native state to a soft viscoelastic intermediate.

Molecular interpretation of single-molecule force spectroscopy experiments with computational approaches

Single molecule force-spectroscopy techniques have granted access to unprecedented molecular-scale details about biochemical and biological mechanisms. However, the interpretation of the experimental data is often challenging. Computational and simulation approaches (all-atom steered MD simulations in particular) are key to provide molecular details about the associated mechanisms, to help test different hypotheses and to predict experimental results. In this review, particular recent efforts directed towards the molecular interpretation of single-molecule force spectroscopy experiments on proteins and protein-related systems (often in close collaboration with experimental groups) will be presented. These results will be discussed in the broader context of the field, highlighting the recent achievements and the ongoing challenges for computational biophysicists and biochemists. In particular, I will focus on the input gained from molecular simulations approaches to rationalize the origin of the unfolded protein elasticity and the protein conformational behavior under force, to understand how force denaturation differs from chemical, thermal or shear unfolding, and to unravel the molecular details of unfolding events for a variety of systems. I will also discuss the use of models based on Langevin dynamics on a 1-D free-energy surface to understand the effect of protein segmentation on the work exerted by a force, or, at the other end of the spectrum of computational techniques, how quantum calculations can help to understand the reactivity of disulfide bridges exposed to force.

Highly Stretchable Fully Biomass Autonomic Self-Healing Polyamide Elastomers and Their Foam for Selective Oil Absorption

Renewable polymers with self-healing ability, excellent elongation, hydrophobicity, and selective oil absorption attributes are of interest for an extensive range of applications, such as e-skin, soft robots, wearable devices, and cleaning up oil spills. Herein, two fully renewable eco-friendly polyamide (PA)-based self-healing elastomers (namely, PA36,IA, and PA36,36) were prepared by a facile and green one-pot melt polycondensation of itaconic acid (IA), PripolTM 1009, and PriamineTM 1075 monomers. The molecular structures of these PAs were analyzed by FITR, 1H NMR, and 13C NMR. The distinct structure of these PAs shows superior strain values (above 2300%) and high ambient temperature autonomous self-healing ability. Interestingly, the synthesized renewable PA36,36 showed zero water absorption values and hydrophobic properties with a contact angle of θ = 91° compared to the synthesized PA36,IA and other previously reported PAs. These excellent attributes are due to the low concentration of amide groups, the highly entangled main chains, the intermolecular diffusion, the manifold dangling chains, and the numerous reversible physical bonds within the renewable PAs. Furthermore, the hydrophobic properties may aid in the selective oil absorption of the PA36,36-based foam, for which PA36,36 foam is produced by the green supercritical carbon dioxide (scCO2) batch foaming process. The PA36,36 foam with a microporous cellular structure showed better absorption capacity and high stability in repeated use. Due to these advantages, these bio-based PAs have potential for the production of eco-friendly self-healing materials, superabsorbent foams, and other polymeric materials.

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.

Interplay between Viscoelasticity and Force Rate Affects Sequential Unfolding in Polyproteins Pulled at Constant Velocity

Polyproteins are unique constructs, comprised of folded protein domains in tandem and polymeric linkers. These macromolecules perform under biological stresses by modulating their response through partial unfolding and extending. Although these unfolding events are considered independent, a history dependence of forced unfolding within polyproteins was reported. Here we measure the unfolding of single poly(I91) octamers, complemented with Brownian dynamics simulations, displaying increasing hierarchy in unfolding-foces, accompanied by a decrease in the effective stiffness. This counters the existing understanding that relates stiffness with variations in domain size and probe stiffness, which is expected to reduce the unfolding forces with every consecutive unfolding event. We utilize a simple mechanistic viscoelastic model to show that two effects are combined within a sequential forced unfolding process: the viscoelastic properties of the growing linker chain lead to a hierarchy of the unfolding events, and force-rate application governs the unfolding kinetics.

The role of near-wall drag effects in the dynamics of tethered DNA under shear flow

We utilized single-molecule tethered particle motion (TPM) tracking, optimized for studying the behavior of short (0.922 μm) dsDNA molecules under shear flow conditions, in the proximity of a wall (surface). These experiments track the individual trajectories through a gold nanobead (40 nm in radius), attached to the loose end of the DNA molecules. Under such circumstances, local interactions with the wall become more pronounced, manifested through hydrodynamic interactions. To elucidate the mechanical mechanism that affects the statistics of the molecular trajectories of the tethered molecules, we estimate the resting diffusion coefficient of our system. Using this value and our measured data, we calculate the orthogonal distance of the extended DNA molecules from the surface. This calculation considers the hydrodynamic drag effect that emerges from the proximity of the molecule to the surface, using the Faxén correction factors. Our finding enables the construction of a scenario according to which the tension along the chain builds up with the applied shear force, driving the loose end of the DNA molecule away from the wall. With the extension from the wall, the characteristic times of the system decrease by three orders of magnitude, while the drag coefficients decay to a plateau value that indicates that the molecule still experiences hydrodynamic effects due to its proximity to the wall.

Can Dissipative Properties of Single Molecules Be Extracted from a Force Spectroscopy Experiment?

We performed dynamic force spectroscopy of single dextran and titin I27 molecules using small-amplitude and low-frequency (40-240 Hz) dithering of an atomic force microscope tip excited by a sine wave voltage fed onto the tip-carrying piezo. We show that for such low-frequency dithering experiments, recorded phase information can be unambiguously interpreted within the framework of a transparent theoretical model that starts from a well-known partial differential equation to describe the dithering of an atomic force microscope cantilever and a single molecule attached to its end system, uses an appropriate set of initial and boundary conditions, and does not exploit any implicit suggestions. We conclude that the observed phase (dissipation) signal is due completely to the dissipation related to the dithering of the cantilever itself (i.e., to the change of boundary conditions in the course of stretching). For both cases, only the upper bound of the dissipation of a single molecule has been established as not exceeding 3⋅10(-7)kg/s. We compare our results with previously reported measurements of the viscoelastic properties of single molecules, and we emphasize that extreme caution must be taken in distinguishing between the dissipation related to the stretched molecule and the dissipation that originates from the viscous damping of the dithered cantilever. We also present the results of an amplitude channel data analysis, which reveal that the typical values of the spring constant of a I27 molecule at the moment of module unfolding are equal to 4±1.5mN/m, and the typical values of the spring constant of dextran at the moment of chair-boat transition are equal to 30-50mN/m.

Advances in the microrheology of complex fluids

New developments in the microrheology of complex fluids are considered. Firstly the requirements for a simple modern particle tracking microrheology experiment are introduced, the error analysis methods associated with it and the mathematical techniques required to calculate the linear viscoelasticity. Progress in microrheology instrumentation is then described with respect to detectors, light sources, colloidal probes, magnetic tweezers, optical tweezers, diffusing wave spectroscopy, optical coherence tomography, fluorescence correlation spectroscopy, elastic- and quasi-elastic scattering techniques, 3D tracking, single molecule methods, modern microscopy methods and microfluidics. New theoretical techniques are also reviewed such as Bayesian analysis, oversampling, inversion techniques, alternative statistical tools for tracks (angular correlations, first passage probabilities, the kurtosis, motor protein step segmentation etc), issues in micro/macro rheological agreement and two particle methodologies. Applications where microrheology has begun to make some impact are also considered including semi-flexible polymers, gels, microorganism biofilms, intracellular methods, high frequency viscoelasticity, comb polymers, active motile fluids, blood clots, colloids, granular materials, polymers, liquid crystals and foods. Two large emergent areas of microrheology, non-linear microrheology and surface microrheology are also discussed.

Dynamics of single-stranded DNA tethered to a solid

Tethering is used to deliver specific biological and industrial functions. For example, single-stranded DNA (ssDNA) is tethered to polymerases and long sequences of double-stranded DNA (dsDNA) during replication, and to solids in DNA microarrays. However, tethering ssDNA to a large object limits not only the available ssDNA conformations, but also the range of time-scales over which the mechanical responses of ssDNA are important. In this work we examine the effect of tethering by measurement of the mechanical response of ssDNA that is tethered at each end to two separate atomic force microscope cantilevers in aqueous solution. Thermal motion of the cantilevers drives the ends of the ssDNA chain at frequencies near 2 kHz. The presence of a tethered molecule makes a large difference to the asymmetric cross-correlation of two cantilevers, which enables resolution of the mechanical properties in our experiments. By analysis of the correlated motion of the cantilevers we extract the friction and stiffness of the ssDNA. We find that the measured friction is much larger than the friction that is usually associated with the unencumbered motion of ssDNA. We also find that the measured relaxation time, ∼30 μs, is much greater than prior measurements of the free-molecule relaxation time. We attribute the difference to the loss of conformational possibilities as a result of constraining the ends of the ssDNA.

Assessment of local friction in protein folding dynamics using a helix cross-linker

Internal friction arising from local steric hindrance and/or the excluded volume effect plays an important role in controlling not only the dynamics of protein folding but also conformational transitions occurring within the native state potential well. However, experimental assessment of such local friction is difficult because it does not manifest itself as an independent experimental observable. Herein, we demonstrate, using the miniprotein trp-cage as a testbed, that it is possible to selectively increase the local mass density in a protein and hence the magnitude of local friction, thus making its effect directly measurable via folding kinetic studies. Specifically, we show that when a helix cross-linker, m-xylene, is placed near the most congested region of the trp-cage it leads to a significant decrease in both the folding rate (by a factor of 3.8) and unfolding rate (by a factor of 2.5 at 35 °C) but has little effect on protein stability. Thus, these results, in conjunction with those obtained with another cross-linked trp-cage and two uncross-linked variants, demonstrate the feasibility of using a nonperturbing cross-linker to help quantify the effect of internal friction. In addition, we estimate that a m-xylene cross-linker could lead to an increase in the roughness of the folding energy landscape by as much as 0.4-1.0k(B)T.

Changes of protein stiffness during folding detect protein folding intermediates

Single-molecule force-quench atomic force microscopy (FQ-AFM) is used to detect folding intermediates of a simple protein by detecting changes of molecular stiffness of the protein during its folding process. Those stiffness changes are obtained from shape and peaks of an autocorrelation of fluctuations in end-to-end length of the folding molecule. The results are supported by predictions of the equipartition theorem and agree with existing Langevin dynamics simulations of a simplified model of a protein folding. In the light of the Langevin simulations the experimental data probe an ensemble of random-coiled collapsed states of the protein, which are present both in the force-quench and thermal-quench folding pathways.

JC polyomavirus (JCV) and monoclonal antibodies: friends or potential foes?

Progressive multifocal leukoencephalopathy (PML) is a demyelinating disease of the central nervous system (CNS), observed in immunodeficient patients and caused by JC virus ((JCV), also called JC polyomavirus (JCPyV)). After the HIV pandemic and the introduction of immunomodulatory therapy, the PML incidence significantly increased. The correlation between the use of natalizumab, a drug used in multiple sclerosis (MS), and the PML development of particular relevance. The high incidence of PML in natalizumab-treated patients has highlighted the importance of two factors: the need of PML risk stratification among natalizumab-treated patients and the need of effective therapeutic options. In this review, we discuss these two needs under the light of the major viral models of PML etiopathogenesis.

Elasticity, structure, and relaxation of extended proteins under force

Force spectroscopies have emerged as a powerful and unprecedented tool to study and manipulate biomolecules directly at a molecular level. Usually, protein and DNA behavior under force is described within the framework of the worm-like chain (WLC) model for polymer elasticity. Although it has been surprisingly successful for the interpretation of experimental data, especially at high forces, the WLC model lacks structural and dynamical molecular details associated with protein relaxation under force that are key to the understanding of how force affects protein flexibility and reactivity. We use molecular dynamics simulations of ubiquitin to provide a deeper understanding of protein relaxation under force. We find that the WLC model successfully describes the simulations of ubiquitin, especially at higher forces, and we show how protein flexibility and persistence length, probed in the force regime of the experiments, are related to how specific classes of backbone dihedral angles respond to applied force. Although the WLC model is an average, backbone model, we show how the protein side chains affect the persistence length. Finally, we find that the diffusion coefficient of the protein's end-to-end distance is on the order of 10(8) nm(2)/s, is position and side-chain dependent, but is independent of the length and independent of the applied force, in contrast with other descriptions.

Rate limit of protein elastic response is tether dependent

The elastic restoring force of tissues must be able to operate over the very wide range of loading rates experienced by living organisms. It is surprising that even the fastest events involving animal muscle tissues do not surpass a few hundred hertz. We propose that this limit is set in part by the elastic dynamics of tethered proteins extending and relaxing under a changing load. Here we study the elastic dynamics of tethered proteins using a fast force spectrometer with sub-millisecond time resolution, combined with Brownian and Molecular Dynamics simulations. We show that the act of tethering a polypeptide to an object, an inseparable part of protein elasticity in vivo and in experimental setups, greatly reduces the attempt frequency with which the protein samples its free energy. Indeed, our data shows that a tethered polypeptide can traverse its free-energy landscape with a surprisingly low effective diffusion coefficient D(eff) ~ 1,200 nm(2)/s. By contrast, our Molecular Dynamics simulations show that diffusion of an isolated protein under force occurs at D(eff) ~ 10(8) nm(2)/s. This discrepancy is attributed to the drag force caused by the tethering object. From the physiological time scales of tissue elasticity, we calculate that tethered elastic proteins equilibrate in vivo with D(eff) ~ 10(4)-10(6) nm(2)/s which is two to four orders magnitude smaller than the values measured for untethered proteins in bulk.

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.

Analyzing single-molecule manipulation experiments

Single-molecule manipulation studies can provide quantitative information about the physical properties of complex biological molecules without ensemble artifacts obscuring the measurements. We demonstrate computational techniques which aim at more fully utilizing the wealth of information contained in noisy experimental time series. The "noise" comes from multiple sources e.g., inherent thermal motion, instrument measurement error, etc. The primary focus of this paper is a methodology that uses time domain based methods to extract the effective molecular friction from single-molecule pulling data. We studied molecules composed of eight tandem repeat titin I27 domains, but the modeling approaches have applicability to other single-molecule mechanical studies. The merits and challenges associated with applying such a computational approach to existing single-molecule manipulation data are also discussed.

Extracting Kinetic and Stationary Distribution Information from Short MD Trajectories via a Collection of Surrogate Diffusion Models

Low-dimensional stochastic models can summarize dynamical information and make long time predictions associated with observables of complex atomistic systems. Maximum likelihood based techniques for estimating low-dimensional surrogate diffusion models from relatively short time series are presented. It is found that a heterogeneous population of slowly evolving conformational degrees of freedom modulates the dynamics. This underlying heterogeneity results in a collection of estimated low-dimensional diffusion models. Numerical techniques for exploiting this finding to approximate skewed histograms associated with the simulation are presented. In addition, statistical tests are also used to assess the validity of the models and determine physically relevant sampling information, e.g. the maximum sampling frequency at which one can discretely sample from an atomistic time series and have a surrogate diffusion model pass goodness-of-fit tests. The information extracted from such analyses can possibly be used to assist umbrella sampling computations as well as help in approximating effective diffusion coefficients. The techniques are demonstrated on simulations of Adenylate Kinase.