Nonlinear elasticity of an alpha -helical polypeptide

Steady-state running rate sets the speed and accuracy of accumulation of swimming bacteria

We study the chemotaxis of a population of genetically identical swimming bacteria undergoing run and tumble dynamics driven by stochastic switching between clockwise and counterclockwise rotation of the flagellar rotary system, where the steady-state rate of the switching changes in different environments. Understanding chemotaxis quantitatively requires that one links the measured steady-state switching rates of the rotary system, as well as the directional changes of individual swimming bacteria in a gradient of chemoattractant/repellant, to the efficiency of a population of bacteria in moving up/down the gradient. Here we achieve this by using a probabilistic model, parametrized with our experimental data, and show that the response of a population to the gradient is complex. We find the changes to the steady-state switching rate in the absence of gradients affect the average speed of the swimming bacterial population response as well as the width of the distribution. Both must be taken into account when optimizing the overall response of the population in complex environments.

Statistics and topology of fluctuating ribbons

Ribbons are a class of slender structures whose length, width, and thickness are widely separated from each other. This scale separation gives a ribbon unusual mechanical properties in athermal macroscopic settings, for example, it can bend without twisting, but cannot twist without bending. Given the ubiquity of ribbon-like biopolymers in biology and chemistry, here we study the statistical mechanics of microscopic inextensible, fluctuating ribbons loaded by forces and torques. We show that these ribbons exhibit a range of topologically and geometrically complex morphologies exemplified by three phases-a twist-dominated helical phase (HT), a writhe-dominated helical phase (HW), and an entangled phase-that arise as the applied torque and force are varied. Furthermore, the transition from HW to HT phases is characterized by the spontaneous breaking of parity symmetry and the disappearance of perversions (that correspond to chirality-reversing localized defects). This leads to a universal response curve of a topological quantity, the link, as a function of the applied torque that is similar to magnetization curves in second-order phase transitions.

Orientational Fluctuations and Bimodality in Semiflexible Nunchucks

Semiflexible nunchucks are block copolymers consisting of two long blocks with high bending rigidity jointed by a short block of lower bending stiffness. Recently, the DNA nanotube nunchuck was introduced as a simple nanoinstrument that mechanically magnifies the bending angle of short double-stranded (ds) DNA and allows its measurement in a straightforward way [Fygenson et al., Nano Lett. 2020, 20, 2, 1388-1395]. It comprises two long DNA nanotubes linked by a dsDNA segment, which acts as a hinge. The semiflexible nunchuck geometry also appears in dsDNA with a hinge defect (e.g., a quenched denaturation bubble or a nick), and in end-linked stiff filaments. In this article, we theoretically investigate various aspects of the conformations and the tensile elasticity of semiflexible nunchucks. We analytically calculate the distribution of bending fluctuations of a wormlike chain (WLC) consisting of three blocks with different bending stiffness. For a system of two weakly bending WLCs end-jointed by a rigid kink, with one end grafted, we calculate the distribution of positional fluctuations of the free end. For a system of two weakly bending WLCs end-jointed by a hinge modeled as harmonic bending spring, with one end grafted, we calculate the positional fluctuations of the free end. We show that, under certain conditions, there is a pronounced bimodality in the transverse fluctuations of the free end. For a semiflexible nunchuck under tension, under certain conditions, there is bimodality in the extension as a function of the hinge position. We also show how steric repulsion affects the bending fluctuations of a rigid-rod nunchuck.

Tensile elasticity of a freely jointed chain with reversible hinges

Many biopolymers exhibit reversible conformational transitions within the chain, which affect their bending stiffness and their response to a stretching force. For example, double stranded DNA may have denatured "bubbles" of unzipped single strands which open and close randomly. In other polymers, the transitions may be due to the reversible attachment and detachment of ligands on ligand-receptor complexes along the backbone. Semiflexible bundles under tension formed by the reversible attachment of cross-linkers, on a coarse-grained level, exhibit similar behaviour. The simplest theoretical model which captures what the above mentioned systems have in common is a freely jointed chain (FJC) with reversible hinges. Each hinge can be open, as in the usual FJC, or closed forcing the adjacent segments to align (stretch). In this article, we analyse it in the Gibbs ensemble. Remarkably, even though the usual FJC in the thermodynamic limit exhibits ensemble equivalence, the reversible FJC exhibits ensemble inequivalence. Even though a mean field treatment suggests a continuous phase transition to a fully hinged state at a certain force, the generating function method ("necklace model") shows that there is no phase transition. However, there is a crossover between the two states with clearly different responses. In the low force (linear response) regime, the reversible FJC has higher tensile compliance than its usual counterpart. In contrast, in the strong force regime, the tensile compliance of the reversible FJC is much lower than that of the usual FJC.

Unfolding mechanism and free energy landscape of single, stable, alpha helices at low pull speeds

Single alpha helices (SAHs) stable in isolated form are often found in motor proteins where they bridge functional domains. Understanding the mechanical response of SAHs is thus critical to understand their function. The quasi-static force-extension relation of a small number of SAHs is known from single-molecule experiments. Unknown, or still controversial, are the molecular scale details behind those observations. We show that the deformation mechanism of SAHs pulled from the termini at pull speeds approaching the quasi-static limit differs from that of typical helices found in proteins, which are stable only when interacting with other protein domains. Using molecular dynamics simulations with atomistic resolution at low pull speeds previously inaccessible to simulation, we show that SAHs start unfolding from the termini at all pull speeds we investigated. Unfolding proceeds residue-by-residue and hydrogen bond breaking is not the main event determining the barrier to unfolding. We use the molecular simulation data to test the cooperative sticky chain model. This model yields excellent fits of the force-extension curves and quantifies the distance, xE = 0.13 nm, to the transition state, the natural frequency of bond vibration, ν0 = 0.82 ns-1, and the height, V0 = 2.9 kcal mol-1, of the free energy barrier associated with the deformation of single residues. Our results demonstrate that the sticky chain model could advantageously be used to analyze experimental force-extension curves of SAHs and other biopolymers.

Single-molecule mechanical unfolding experiments reveal a critical length for the formation of α-helices in peptides

α-Helix is the most predominant secondary structure in proteins and supports many functions in biological machineries. The conformation of the helix is dictated by many factors such as its primary sequence, intramolecular interactions, or the effect of the close environment. Several computational studies have proposed that there is a critical maximum length for the formation of intact compact helical structures, supporting the fact that most intact α-helices in proteins are constituted of a small number of amino acids. To obtain a detailed picture on the formation of α-helices in peptides and their mechanical stability, we have synthesized a long homopolypeptide of about 90 amino acids, poly(γ-benzyl-l-glutamate), and investigated its mechanical behaviour by AFM-based single-molecule force spectroscopy. The characteristic plateaus observed in the force-extension curves reveal the unfolding of a series of small helices (from 1 to 4) of about 20 amino acid residues connected to each other, rather than a long helix of 90 residues. Our results suggest the formation of a tertiary structure made of short helices with kinks, instead of an intact compact helical structure for sequences of more than 20 amino acid residues. To our knowledge, this is the first experimental evidence supporting the concept of a helical critical length previously proposed by several computational studies.

Bubble-raft collapse and the nonequilibrium dynamics of two-state elastica

We report on the collapse of bubble rafts under compression in a closed rectangular geometry. A bubble raft is a single layer of bubbles at the air-water interface. A collapse event occurs when bubbles submerge beneath the neighboring bubbles under compression, causing the structure of the bubble raft to go from single-layer to multilayer. We studied the collapse dynamics as a function of compression velocity. At higher compression velocity we observe a more uniform distribution of collapse events, whereas at lower compression velocities the collapse events accumulate at the system boundaries. We propose that this system can be understood in terms of a linear elastic sheet coupled to a local internal (Ising) degree of freedom. The two internal states, which represent one bubble layer versus two, couple to the elasticity of the sheet by locally changing the reference state of the material. By exploring the collapse dynamics of the bubble raft, one may address the basic nonlinear mechanics of a number of complex systems in which elastic stress is coupled to local internal variables.

Buckling transition in long α-helices

The treatment of bending and buckling of stiff biopolymer filaments by the popular worm-like chain model does not provide adequate understanding of these processes at the microscopic level. Thus, we have used the atomistic molecular-dynamic simulations and the Amber03 force field to examine the compression buckling of α-helix (AH) filaments at room temperature. It was found that the buckling instability occurs in AHs at the critical force f(c) in the range of tens of pN depending on the AH length. The decrease of the force f(c) with the contour length follows the prediction of the classic thin rod theory. At the force f(c) the helical filament undergoes the swift and irreversible transition from the smoothly bent structure to the buckled one. A sharp kink in the AH contour arises at the transition, accompanied by the disruption of the hydrogen bonds in its vicinity. The kink defect brings in an effective softening of the AH molecule at buckling. Nonbonded interactions between helical branches drive the rearrangement of a kinked AH into the ultimate buckled structure of a compact helical hairpin described earlier in the literature.

Coil-helix transition of polypeptide at water-lipid interface

We present the exact solution of a microscopic statistical mechanical model for the transformation of a long polypeptide between an unstructured coil conformation and an α-helix conformation. The polypeptide is assumed to be adsorbed to the interface between a polar and a non-polar environment such as realized by water and the lipid bilayer of a membrane. The interfacial coil-helix transformation is the first stage in the folding process of helical membrane proteins. Depending on the values of model parameters, the conformation changes as a crossover, a discontinuous transition, or a continuous transition with helicity in the role of order parameter. Our model is constructed as a system of statistically interacting quasiparticles that are activated from the helix pseudo-vacuum. The particles represent links between adjacent residues in coil conformation that form a self-avoiding random walk in two dimensions. Explicit results are presented for helicity, entropy, heat capacity, and the average numbers and sizes of sboth coil and helix segments.

Cavity approach for modeling and fitting polymer stretching

The mechanical properties of molecules are today captured by single molecule manipulation experiments, so that polymer features are tested at a nanometric scale. Yet devising mathematical models to get further insight beyond the commonly studied force-elongation relation is typically hard. Here we draw from techniques developed in the context of disordered systems to solve models for single and double-stranded DNA stretching in the limit of a long polymeric chain. Since we directly derive the marginals for the molecule local orientation, our approach allows us to readily calculate the experimental elongation as well as other observables at wish. As an example, we evaluate the correlation length as a function of the stretching force. Furthermore, we are able to fit successfully our solution to real experimental data. Although the model is admittedly phenomenological, our findings are very sound. For single-stranded DNA our solution yields the correct (monomer) scale and yet, more importantly, the right persistence length of the molecule. In the double-stranded case, our model reproduces the well-known overstretching transition and correctly captures the ratio between native DNA and overstretched DNA. Also in this case the model yields a persistence length in good agreement with consensus, and it gives interesting insights into the bending stiffness of the native and overstretched molecule, respectively.

High-energy deformation of filaments with internal structure and localized torque-induced melting of DNA

We develop a continuum elastic approach to examining the bending mechanics of semiflexible filaments with a local internal degree of freedom that couples to the bending modulus. We apply this model to study the nonlinear mechanics of a double-stranded DNA oligomer (shorter than its thermal persistence length) whose free ends are linked by a single-stranded DNA chain. This construct, studied by H. Qu and G. Zocchi [Europhys. Lett. 94, 18003 (2011)], displays nonlinear strain softening associated with the local melting of the double-stranded DNA under applied torque and serves as a model system with which to study the nonlinear elasticity of DNA under large energy deformations. We show that one can account quantitatively for the observed bending mechanics using an augmented wormlike chain model, the helix-coil wormlike chain. We also predict that the highly bent and partially molten dsDNA should exhibit particularly large end-to-end fluctuations associated with the fluctuation of the length of the molten region, and propose appropriate experimental tests. We suggest that the augmented wormlike chain model discussed here is a useful analytic approach to the nonlinear mechanics of DNA or other biopolymer systems.

Single molecule measurements of mechanical interactions within ternary SNARE complexes and dynamics of their disassembly: SNAP25 vs. SNAP23

Regulated exocytosis is a crucial event for intercellular communication between neurons and astrocytes within the CNS. The soluble N-ethylmaleimide-sensitive fusion protein attachment protein receptor (SNARE) complex, composed of synaptobrevin 2, syntaxin and synaptosome-associated protein of 25 kDa or 23 kDa (SNAP25 or SNAP23), is essential in this process. It was reported that SNAP25 and SNAP23 have distinct roles in exocytotic release, where SNAP25, but not SNAP23, supports an exocytotic burst. It is not clear, however, whether this is due to the intrinsic properties of the ternary SNARE complex, containing either SNAP25 or SNAP23, or perhaps due to the differential association of these proteins with ancillary proteins to the complex. Here, using force spectroscopy, we show from single molecule investigations of the SNARE complex, that SNAP23A created a local interaction at the ionic layer by cuffing syntaxin 1A and synaptobrevin 2, similar to the action of SNAP25B; thus either of the ternary complexes would allow positioning of vesicles at a maximal distance of approximately 13 nm from the plasma membrane. However, the stability of the ternary SNARE complex containing SNAP23A is less than half of that for the complex containing SNAP25B. Thus, differences in the stability of the two different ternary complexes could underlie some of the SNAP25/23 differential ability to control the exocytotic burst.

Free energies and forces in helix-coil transition of homopolypeptides under stretching

We show here that constant velocity steered molecular dynamics (SMD) simulations of alpha-helices in a vacuum present a well defined plateau in the force-extension relationship for homopolypeptides having more than (approximately) twenty residues. With the processes being far away from equilibrium, the energies strongly depend on the stretching velocity. Importantly, for a given velocity variation, the energy variation depends also on the helix sequence. Additionally, our observations show that homopolypeptides made of ten different amino acids (Ala, Cys, Gln, Ile, Leu, Met, Phe, Ser, Thr and Val) present a linear helix-coil transition.

Comparative energy measurements in single molecule interactions

Single molecule experiments have opened promising new avenues of investigations in biology, but the quantitative interpretation of results remains challenging. In particular, there is a need for a comparison of such experiments with theoretical methods. We experimentally determine the activation free energy for single molecule interactions between two synaptic proteins syntaxin 1A and synaptobrevin 2, using an atomic force microscope and the Jarzynski equality of nonequilibrium thermodynamics. The value obtained is shown to be reasonably consistent with that from single molecule reaction rate theory. The temperature dependence of the spontaneous dissociation lifetime along with different pulling speeds is used to confirm the approach to the adiabatic limit. This comparison of the Jarzynski equality for intermolecular interactions extends the procedure for calculation of activation energies in nonequilibrium processes.

Nonlinear elasticity of an alpha-helical polypeptide: Monte Carlo studies

We report on Monte Carlo studies of the elastic properties of the helix-coil wormlike chain model of alpha-helical polypeptides. In this model the secondary structure enters as a scalar (Ising-like) variable that controls the local chain bending modulus. We characterize the nonlinear elastic properties of these molecules including their response to applied tensile forces and bending torques both individually and in combination. We find a pronounced effect of applied torque on the extensional compliance of the molecule and a similar effect of tension on the bending compliance. Finally we speculate that the strongly nonlinear response of alpha-helical polypeptides to combinations of torque and force plays a role in allosteric transitions in proteins.

Loops in DNA: an overview of experimental and theoretical approaches

DNA loop formation plays a central role in many cellular processes. The aim of this paper is to present the state of the art and open problems regarding the experimental and theoretical approaches to DNA looping. A particular attention is devoted to the effects of the protein bridge size and of protein induced sharp DNA bending on DNA loop formation enhancement.

Ribosome exit tunnel can entropically stabilize alpha-helices

Several experiments have suggested that newly synthesized polypeptide chains can adopt helical structures deep within the ribosome exit tunnel. We hypothesize that confinement in the roughly cylindrical tunnel can entropically stabilize alpha-helices. The hypothesis is validated by using theory and simulations of coarse-grained off-lattice models. The model helix, which is unstable in the bulk, is stabilized in a cylindrical cavity provided the diameter (D) of the cylinder exceeds a critical value D*. When D < D* both the helical content and the helix-coil transition temperature (T(f)) decrease abruptly. Surprisingly, we find that the stability of the alpha-helix depends on the number (N) of amino acid residues. Entropic stabilization, as measured by changes in T(f), increases nonlinearly as N increases. The simulation results are in quantitative agreement with a standard helix-coil theory that takes into account entropy cost of confining a polypeptide chain in a cylinder. The results of this work are in qualitative accord with most of the findings of a recent experiment in which N-dependent ribosome-induced helix stabilization of transmembrane sequences was measured by fluorescence resonance energy transfer.

Protein-mediated DNA loops: effects of protein bridge size and kinks

This paper focuses on the probability that a portion of DNA closes on itself through thermal fluctuations. We investigate the dependence of this probability upon the size of a protein bridge and/or the presence of a kink at half DNA length. The DNA is modeled by the wormlike chain model, and the probability of loop formation is calculated in two ways: exact numerical evaluation of the constrained path integral and the extension of the Shimada and Yamakawa saddle point approximation. For example, we find that the looping free energy of a 100-base-pairs DNA decreases from 24 kBT to 13 kBT when the loop is closed by a protein of r=10 length. It further decreases 5 kBT to when the loop has a kink of 120 degrees at half-length.

Elasticity of semiflexible polymers in two dimensions

We study theoretically the entropic elasticity of a semiflexible polymer, such as DNA, confined to two dimensions. Using the worm-like-chain model we obtain an exact analytical expression for the partition function of the polymer pulled at one end with a constant force. The force-extension relation for the polymer is computed in the long chain limit in terms of Mathieu characteristic functions. We also present applications to the interaction between a semiflexible polymer and a nematic field, and derive the nematic order parameter and average extension of the polymer in a strong field.