The molecular mechanism underlying mechanical anisotropy of the protein GB1

Engineering an artificial catch bond using mechanical anisotropy

Catch bonds are a rare class of protein-protein interactions where the bond lifetime increases under an external pulling force. Here, we report how modification of anchor geometry generates catch bonding behavior for the mechanostable Dockerin G:Cohesin E (DocG:CohE) adhesion complex found on human gut bacteria. Using AFM single-molecule force spectroscopy in combination with bioorthogonal click chemistry, we mechanically dissociate the complex using five precisely controlled anchor geometries. When tension is applied between residue #13 on CohE and the N-terminus of DocG, the complex behaves as a two-state catch bond, while in all other tested pulling geometries, including the native configuration, it behaves as a slip bond. We use a kinetic Monte Carlo model with experimentally derived parameters to simulate rupture force and lifetime distributions, achieving strong agreement with experiments. Single-molecule FRET measurements further demonstrate that the complex does not exhibit dual binding mode behavior at equilibrium but unbinds along multiple pathways under force. Together, these results show how mechanical anisotropy and anchor point selection can be used to engineer artificial catch bonds.

Uncovering and engineering the mechanical properties of the adhesion GPCR ADGRG1 GAIN domain

Key cellular functions depend on the transduction of extracellular mechanical signals by specialized membrane receptors including adhesion G-protein coupled receptors (aGPCRs). While recently solved structures support aGPCR activation through shedding of the extracellular GAIN domain, the molecular mechanisms underpinning receptor mechanosensing remain poorly understood. When probed using single-molecule atomic force spectroscopy and molecular simulations, ADGRG1 GAIN dissociated from its tethered agonist at forces significantly higher than other reported signaling mechanoreceptors. Strong mechanical resistance was achieved through specific structural deformations and force propagation pathways under mechanical load. ADGRG1 GAIN variants computationally designed to lock the alpha and beta subdomains and rewire mechanically-induced structural deformations were found to modulate the GPS-Stachel rupture forces. Our study provides unprecedented insights into the molecular underpinnings of GAIN mechanical stability and paves the way for engineering mechanosensors, better understanding aGPCR function, and informing drug-discovery efforts targeting this important receptor class.

The transient expression of recombinant proteins in plant cell packs facilitates stable isotope labelling for NMR spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy can be used to determine the structure, dynamics and interactions of proteins. However, protein NMR requires stable isotope labelling for signal detection. The cells used for the production of recombinant proteins must therefore be grown in medium containing isotopically labelled substrates. Stable isotope labelling is well established in Escherichia coli, but bacteria are only suitable for the production of simple proteins without post-translational modifications. More complex proteins require eukaryotic production hosts, but their growth can be impaired by labelled media, thus reducing product yields and increasing costs. To address this limitation, we used media supplemented with isotope-labelled substrates to cultivate the tobacco-derived cell line BY-2, which was then cast into plant cell packs (PCPs) for the transient expression of a labelled version of the model protein GB1. Mass spectrometry confirmed the feasibility of isotope labelling with ¹⁵ N and ² H using this approach. The resulting NMR spectrum featured a signal dispersion comparable to recombinant GB1 produced in E. coli. PCPs therefore offer a rapid and cost-efficient alternative for the production of isotope-labelled proteins for NMR analysis, especially suitable for complex proteins that cannot be produced in microbial systems.

Contribution of Hydrophobic Interactions to Protein Mechanical Stability

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

Tug-of-War between Internal and External Frictions and Viscosity Dependence of Rate in Biological Reactions

The role of water in biological processes is studied in three reactions, namely, the Fe-CO bond rupture in myoglobin, GB1 unfolding, and insulin dimer dissociation. We compute both internal and external components of friction on relevant reaction coordinates. In all of the three cases, the cross-correlation between forces from protein and water is found to be large and negative that serves to reduce the total friction significantly, increase the calculated reaction rate, and weaken solvent viscosity dependence. The computed force spectrum reveals bimodal 1/f noise, suggesting the use of a non-Markovian rate theory.

Mapping Mechanostable Pulling Geometries of a Therapeutic Anticalin/CTLA-4 Protein Complex

We used single-molecule AFM force spectroscopy (AFM-SMFS) in combination with click chemistry to mechanically dissociate anticalin, a non-antibody protein binding scaffold, from its target (CTLA-4), by pulling from eight different anchor residues. We found that pulling on the anticalin from residue 60 or 87 resulted in significantly higher rupture forces and a decrease in koff by 2-3 orders of magnitude over a force range of 50-200 pN. Five of the six internal anchor points gave rise to complexes significantly more stable than N- or C-terminal anchor points, rupturing at up to 250 pN at loading rates of 0.1-10 nN s-1. Anisotropic network modeling and molecular dynamics simulations helped to explain the geometric dependency of mechanostability. These results demonstrate that optimization of attachment residue position on therapeutic binding scaffolds can provide large improvements in binding strength, allowing for mechanical affinity maturation under shear stress without mutation of binding interface residues.

Anisotropy in mechanical unfolding of protein upon partner-assisted pulling and handle-assisted pulling

Proteins as force-sensors respond to mechanical cues and regulate signaling in physiology. Proteins commonly connect the source and response points of mechanical cues in two conformations, independent proteins in end-to-end geometry and protein complexes in handshake geometry. The force-responsive property of independent proteins in end-to-end geometry is studied extensively using single-molecule force spectroscopy (SMFS). The physiological significance of the complex conformations in force-sensing is often disregarded as mere surge protectors. However, with the potential of force-steering, protein complexes possess a distinct mechano-responsive property over individual force-sensors. To decipher, we choose a force-sensing protein, cadherin-23, from tip-link complex and perform SMFS using end-to-end geometry and handshake complex geometry. We measure higher force-resilience of cadherin-23 with preferential shorter extensions in handshake mode of pulling over the direct mode. The handshake geometry drives the force-response of cadherin-23 through different potential-energy landscapes than direct pulling. Analysis of the dynamic network structure of cadherin-23 under tension indicates narrow force-distributions among residues in cadherin-23 in direct pulling, resulting in low force-dissipation paths and low resilience to force. Overall, the distinct and superior mechanical responses of cadherin-23 in handshake geometry than single protein geometry highlight a probable evolutionary drive of protein-protein complexes as force-conveyors over independent ones.

Estimating the Directional Flexibility of Proteins from Equilibrium Thermal Fluctuations

The directional flexibility of proteins is an equilibrium molecular property which is accessible to both experiment and computation. Single molecule force spectroscopy (SMFS) experiments report effective directional spring constants to describe the collective anisotropic response of a protein structure to mechanical pulling forces applied along selected axes. On the other hand, computational methods have thus far employed either indirect force based nonequilibrium simulations or coarse-grained elastic network models (ENM) to predict protein directional spring constants. Here, we examine the ability of equilibrium atomistic Molecular Dynamics (MD) simulations to estimate the directional flexibility and mechanical anisotropy of proteins. MD-derived effective directional spring constants are found to correlate well with SMFS spring constants (ρ² = 0.97-0.99; Adj R² = 0.92-0.99) and unfolding forces (ρ² = 0.85-0.97; Adj R² = 0.63-0.91) for five different globular proteins. Specifically, the computed spring constants reproduce the mechanical anisotropy reported by SMFS along five different directions of green fluorescence protein (GFP) and six directions of the immunoglobulin-binding B1 domain of streptococcal protein G (GB1). Further, protein dynamics as captured in MD can be translated into spring constants which can distinguish the N-C directional flexibility of ubiquitin (Ub) from two structurally homologous small ubiquitin-like modifier (SUMO1 and SUMO2) isoforms. We apply our computational framework to study the mechanical anisotropy of Ub along the seven lysine-C-term directions which are functionally relevant. We show that Ub possesses two distinct flexibility scales along these directions which roughly differ by an order of magnitude. Further, our studies reveal that the mechanical anisotropy of Ub is modified in contrasting ways by the binding of two partner proteins (UBCH5A and UEV) which attach and recognize these biomolecular tag proteins. On the basis of equilibrium MD benchmarks for flexibility along 2485 bond vectors in Ub, we propose and validate a new covariance-propagation scheme to extract spring constants from ENM normal modes. We also critically examine the ability of ENM to predict directional flexibility of proteins and suggest modifications to improve these intuitive and scalable descriptions.

Hidden Intermediate State and Second Pathway Determining Folding and Unfolding Dynamics of GB1 Protein at Low Forces

Atomic force microscopy experiments found that GB1, a typical two-state model protein used for study of folding and unfolding dynamics, can sustain forces of more than 100 pN, but its response to low forces still remains unclear. Using ultrastable magnetic tweezers, we discovered that GB1 has an unexpected nonmonotonic force-dependent unfolding rate at 5-160 pN, from which a free energy landscape with two main barriers and a hidden intermediate state was constructed. A model combining two separate models by Dudko et al. with two pathways between the native state and this intermediate state is proposed to rebuild the unfolding dynamics over the full experimental force range. One candidate of this transient intermediate state is the theoretically proposed molten globule state with a loosely collapsed conformation, which might exist universally in the folding and unfolding processes of two-state proteins.

Binding of a Telomestatin Derivative Changes the Mechanical Anisotropy of a Human Telomeric G-Quadruplex

Mechanical anisotropy is an essential property for biomolecules to assume structural and functional roles in mechanobiology. However, there is insufficient information on the mechanical anisotropy of ligand-biomolecule complexes. Herein, we investigated the mechanical property of individual human telomeric G-quadruplexes bound to telomestatin, using optical tweezers. Stacking of the ligand to the G-tetrad planes changes the conformation of the G-quadruplex, which resembles a balloon squeezed in certain directions. Such a squeezed balloon effect strengthens the G-tetrad planes, but dislocates and weakens the loops in the G-quadruplex upon ligand binding. These dynamic interactions indicate that the binding between the ligand and G-quadruplex follows the induced-fit model. We anticipate that the altered mechanical anisotropy of the ligand-G-quadruplex complex can add additional level of regulations on the motor enzymes that process DNA or RNA molecules.

Mechanical Anisotropy in GNNQQNY Amyloid Crystals

Mapping the nanomechanical properties of amyloids can provide valuable insights into structure and assembly mechanisms of protein aggregates that underlie the development of various human diseases. Although it is well-known that amyloids exhibit an intrinsic stiffness comparable to that of silk (1-10 GPa), a detailed understanding of the directional dependence (anisotropy) of the stiffness of amyloids and how it relates to structural features in these protein aggregates is missing. Here we used steered molecular dynamics (SMD) simulations and amplitude modulation-frequency modulation (AM-FM) atomic force microscopy to measure the directional variation in stiffness of GNNQQNY amyloid crystals. We reveal that individual crystals display significant mechanical anisotropy and relate this anisotropy to subtle but mechanically important differences in interactions between interfaces that define the crystal architecture. Our results provide detailed insights into the structure-mechanics relationship of amyloid that may help in designing amyloid-based nanomaterials with tailored mechanical properties.

Mg2+ -Dependent High Mechanical Anisotropy of Three-Way-Junction pRNA as Revealed by Single-Molecule Force Spectroscopy

Mechanical anisotropy is ubiquitous in biological tissues but is hard to reproduce in synthetic biomaterials. Developing molecular building blocks with anisotropic mechanical response is the key towards engineering anisotropic biomaterials. The three-way-junction (3WJ) pRNA, derived from ϕ29 DNA packaging motor, shows strong mechanical anisotropy upon Mg²⁺ binding. In the absence of Mg²⁺ , 3WJ-pRNA is mechanically weak without noticeable mechanical anisotropy. In the presence of Mg²⁺ , the unfolding forces can differ by more than 4-fold along different pulling directions, ranging from about 47 pN to about 219 pN. Mechanical anisotropy of 3WJ-pRNA stems from pulling direction dependent cooperativity for the rupture of two Mg²⁺ binding sites, which is a novel mechanism for the mechanical anisotropy of biomacromolecules. It is anticipated that 3WJ-pRNA can be used as a key element for the construction of biomaterials with controllable mechanical anisotropy.

Simulated Force Quench Dynamics Shows GB1 Protein Is Not a Two State Folder

Single molecule force spectroscopy is a useful technique for investigating mechanically induced protein unfolding and refolding under reduced forces by monitoring the end-to-end distance of the protein. The data is often interpreted via a "two-state" model based on the assumption that the end-to-end distance alone is a good reaction coordinate and the thermodynamic behavior is then ascribed to the free energy as a function of this one reaction coordinate. In this paper, we determined the free energy surface (PMF) of GB1 protein from atomistic simulations in explicit solvent under different applied forces as a function of two collective variables (the end-to-end-distance, and the fraction of native contacts ρ). The calculated 2-d free energy surfaces exhibited several distinct states, or basins, mostly visible along the ρ coordinate. Brownian dynamics (BD) simulations on the smoothed free energy surface show that the protein visits a metastable molten globule state and is thus a three state folder, not the two state folder inferred using the end-to-end distance as the sole reaction coordinate. This study lends support to recent experiments that suggest that GB1 is not a two-state folder.

The physics of pulling polyproteins: a review of single molecule force spectroscopy using the AFM to study protein unfolding

One of the most exciting developments in the field of biological physics in recent years is the ability to manipulate single molecules and probe their properties and function. Since its emergence over two decades ago, single molecule force spectroscopy has become a powerful tool to explore the response of biological molecules, including proteins, DNA, RNA and their complexes, to the application of an applied force. The force versus extension response of molecules can provide valuable insight into its mechanical stability, as well as details of the underlying energy landscape. In this review we will introduce the technique of single molecule force spectroscopy using the atomic force microscope (AFM), with particular focus on its application to study proteins. We will review the models which have been developed and employed to extract information from single molecule force spectroscopy experiments. Finally, we will end with a discussion of future directions in this field.

Coarse-Grained Simulations of Topology-Dependent Mechanisms of Protein Unfolding and Translocation Mediated by ClpY ATPase Nanomachines

Clp ATPases are powerful ring shaped nanomachines which participate in the degradation pathway of the protein quality control system, coupling the energy from ATP hydrolysis to threading substrate proteins (SP) through their narrow central pore. Repetitive cycles of sequential intra-ring ATP hydrolysis events induce axial excursions of diaphragm-forming central pore loops that effect the application of mechanical forces onto SPs to promote unfolding and translocation. We perform Langevin dynamics simulations of a coarse-grained model of the ClpY ATPase-SP system to elucidate the molecular details of unfolding and translocation of an α/β model protein. We contrast this mechanism with our previous studies which used an all-α SP. We find conserved aspects of unfolding and translocation mechanisms by allosteric ClpY, including unfolding initiated at the tagged C-terminus and translocation via a power stroke mechanism. Topology-specific aspects include the time scales, the rate limiting steps in the degradation pathway, the effect of force directionality, and the translocase efficacy. Mechanisms of ClpY-assisted unfolding and translocation are distinct from those resulting from non-allosteric mechanical pulling. Bulk unfolding simulations, which mimic Atomic Force Microscopy-type pulling, reveal multiple unfolding pathways initiated at the C-terminus, N-terminus, or simultaneously from both termini. In a non-allosteric ClpY ATPase pore, mechanical pulling with constant velocity yields larger effective forces for SP unfolding, while pulling with constant force results in simultaneous unfolding and translocation.

Cell survival is critically dependent on tightly regulated protein quality control, which includes chaperone-mediated folding and degradation. In the degradation pathway, AAA+ nanomachines, such as bacterial Clp proteases, use ATP-driven mechanisms to mechanically unfold, translocate, and destroy excess or defective proteins. Understanding these remodeling mechanisms is of central importance for deciphering the details of essential cellular processes. We perform coarse-grained computer simulations to extensively probe the effect of substrate protein topology on unfolding and translocation actions of the ClpY ATPase nanomachine. We find that, independent of SP topology, unfolding proceeds from the tagged C-terminus, which is engaged by the ATPase, and translocation involves coordinated steps. Topology-specific aspects include more complex unfolding and translocation pathways of the α/β SP compared with the all-α SP due to high stability of β-hairpins and interplay of tertiary contacts. In addition, directionality of the mechanical force applied by the Clp ATPase gives rise to distinct unfolding pathways.

Monovalent Strep-Tactin for strong and site-specific tethering in nanospectroscopy

Strep-Tactin, an engineered form of streptavidin, binds avidly to the genetically encoded peptide Strep-tag II in a manner comparable to streptavidin binding to biotin. These interactions have been used in protein purification and detection applications. However, in single-molecule studies, for example using atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS), the tetravalency of these systems impedes the measurement of monodispersed data. Here, we introduce a monovalent form of Strep-Tactin that harbours a unique binding site for Strep-tag II and a single cysteine that allows Strep-Tactin to specifically attach to the atomic force microscope cantilever and form a consistent pulling geometry to obtain homogeneous rupture data. Using AFM-SMFS, the mechanical properties of the interaction between Strep-tag II and monovalent Strep-Tactin were characterized. Rupture forces comparable to biotin:streptavidin unbinding were observed. Using titin kinase and green fluorescent protein, we show that monovalent Strep-Tactin is generally applicable to protein unfolding experiments. We expect monovalent Strep-Tactin to be a reliable anchoring tool for a range of single-molecule studies.

Mapping Mechanical Force Propagation through Biomolecular Complexes

Here we employ single-molecule force spectroscopy with an atomic force microscope (AFM) and steered molecular dynamics (SMD) simulations to reveal force propagation pathways through a mechanically ultrastable multidomain cellulosome protein complex. We demonstrate a new combination of network-based correlation analysis supported by AFM directional pulling experiments, which allowed us to visualize stiff paths through the protein complex along which force is transmitted. The results implicate specific force-propagation routes nonparallel to the pulling axis that are advantageous for achieving high dissociation forces.

Nucleotides regulate the mechanical hierarchy between subdomains of the nucleotide binding domain of the Hsp70 chaperone DnaK

The regulation of protein function through ligand-induced conformational changes is crucial for many signal transduction processes. The binding of a ligand alters the delicate energy balance within the protein structure, eventually leading to such conformational changes. In this study, we elucidate the energetic and mechanical changes within the subdomains of the nucleotide binding domain (NBD) of the heat shock protein of 70 kDa (Hsp70) chaperone DnaK upon nucleotide binding. In an integrated approach using single molecule optical tweezer experiments, loop insertions, and steered coarse-grained molecular simulations, we find that the C-terminal helix of the NBD is the major determinant of mechanical stability, acting as a glue between the two lobes. After helix unraveling, the relative stability of the two separated lobes is regulated by ATP/ADP binding. We find that the nucleotide stays strongly bound to lobe II, thus reversing the mechanical hierarchy between the two lobes. Our results offer general insights into the nucleotide-induced signal transduction within members of the actin/sugar kinase superfamily.

Mechanical dissociation of the M-band titin/obscurin complex is directionally dependent

Titin and obscurin, two giant muscle proteins, bind to each other in an antiparallel Ig-Ig fashion at the M-band. This interaction must be able to withstand the mechanical strain that the M-band typically experiences and remain intact. The mechanical force on these domains is likely exerted along one of two axes: a longitudinal axis, resulting in a 'shearing' force, or a lateral axis, resulting in a 'peeling' force. Here we present molecular dynamics data suggesting that these forces result in distinct unraveling pathways of the titin/obscurin complex and that peeling the domains apart requires less work and force.

Polyproteins in structural biology

•

Structures have been determined for natural and recombinant polyproteins.

•

Native HIV Gag polyprotein architecture was revealed by cryo-EM of immature capsids.

•

Recombinant polyprotein technology has resolved sample preparation bottlenecks.

•

The high-resolution structure of influenza polymerase has been solved.

•

Single-molecule analysis of polyproteins revealed their folding characteristics.

Polyproteins are chains of covalently conjoined smaller proteins that occur in nature as versatile means to organize the proteome of viruses including HIV. During maturation, viral polyproteins are typically cleaved into the constituent proteins with different biological functions by highly specific proteases, and structural analyses at defined stages of this maturation process can provide clues for antiviral intervention strategies. Recombinant polyproteins that use similar mechanisms are emerging as powerful tools for producing hitherto inaccessible protein targets such as the influenza polymerase, for high-resolution structure determination by X-ray crystallography. Conversely, covalent linking of individual protein subunits into single polypeptide chains are exploited to overcome sample preparation bottlenecks. Moreover, synthetic polyproteins provide a promising tool to dissect dynamic folding of polypeptide chains into three-dimensional architectures in single-molecule structure analysis by atomic force microscopy (AFM). The recent use of natural and synthetic polyproteins in structural biology and major achievements are highlighted in this contribution.

Extraction of accurate biomolecular parameters from single-molecule force spectroscopy experiments

The atomic force microscope (AFM) is able to manipulate biomolecules and their complexes with exquisite force sensitivity and distance resolution. This capability, complemented by theoretical models, has greatly improved our understanding of the determinants of mechanical strength in proteins and revealed the diverse effects of directional forces on the energy landscape of biomolecules. In unbinding experiments, the interacting partners are usually immobilized on their respective substrates via extensible linkers. These linkers affect both the force and contour length (Lc) of the complex at rupture. Surprisingly, while the former effect is well understood, the latter is largely neglected, leading to incorrect estimations of Lc, a parameter that is often used as evidence for the detection of specific interactions and remodeling events and for the inference of interaction regions. To address this problem, a model that predicts contour length measurements from single-molecule forced-dissociation experiments is presented that considers attachment position on the AFM tip, geometric effects, and polymer dynamics of the linkers. Modeled data are compared with measured contour length distributions from several different experimental systems, revealing that current methods underestimate contour lengths. The model enables nonspecific interactions to be identified unequivocally, allows accurate determination of Lc, and, by comparing experimental and modeled distributions, enables partial unfolding events before rupture to be identified unequivocally.

Force-induced remodelling of proteins and their complexes

Force can drive conformational changes in proteins, as well as modulate their stability and the affinity of their complexes, allowing a mechanical input to be converted into a biochemical output. These properties have been utilised by nature and force is now recognised to be widely used at the cellular level. The effects of force on the biophysical properties of biological systems can be large and varied. As these effects are only apparent in the presence of force, studies on the same proteins using traditional ensemble biophysical methods can yield apparently conflicting results. Where appropriate, therefore, force measurements should be integrated with other experimental approaches to understand the physiological context of the system under study.

Effective harmonic potentials: insights into the internal cooperativity and sequence-specificity of protein dynamics

The proper biological functioning of proteins often relies on the occurrence of coordinated fluctuations around their native structure, or on their ability to perform wider and sometimes highly elaborated motions. Hence, there is considerable interest in the definition of accurate coarse-grained descriptions of protein dynamics, as an alternative to more computationally expensive approaches. In particular, the elastic network model, in which residue motions are subjected to pairwise harmonic potentials, is known to capture essential aspects of conformational dynamics in proteins, but has so far remained mostly phenomenological, and unable to account for the chemical specificities of amino acids. We propose, for the first time, a method to derive residue- and distance-specific effective harmonic potentials from the statistical analysis of an extensive dataset of NMR conformational ensembles. These potentials constitute dynamical counterparts to the mean-force statistical potentials commonly used for static analyses of protein structures. In the context of the elastic network model, they yield a strongly improved description of the cooperative aspects of residue motions, and give the opportunity to systematically explore the influence of sequence details on protein dynamics.