Ligand Binding Mechanics of Maltose Binding Protein

Protein engineering using circular permutation - structure, function, stability, and applications

Protein engineering is important for creating novel variants from natural proteins, enabling a wide range of applications. Approaches such as rational design and directed evolution are routinely used to make new protein variants. Computational tools like de novo design can introduce new protein folds. Expanding the amino acid repertoire to include unnatural amino acids with non-canonical side chains in vitro by native chemical ligation and in vivo via codon expansion methods broadens sequence and structural possibilities. Circular permutation (CP) is an invaluable approach to redesigning a protein by rearranging the amino acid sequence, where the connectivity of the secondary structural elements is altered without changing the overall structure of the protein. Artificial CP proteins (CPs) are employed in various applications such as biocatalysis, sensing of small molecules by fluorescence, genome editing, ligand-binding protein switches, and optogenetic engineering. Many studies have shown that CP can lead to either reduced or enhanced stability or catalytic efficiency. The effects of CP on a protein's energy landscape cannot be predicted a priori. Thus, it is important to understand how CP can affect the thermodynamic and kinetic stability of a protein. In this review, we discuss the discovery and advancement of techniques to create protein CP, and existing reviews on CP. We delve into the plethora of biological applications for designed CP proteins. We subsequently discuss the experimental and computational reports on the effects of CP on the thermodynamic and kinetic stabilities of proteins of various topologies. An understanding of the various aspects of CP will allow the reader to design robust CP proteins for their specific purposes.

Probing the stability and interdomain interactions in the ABC transporter OpuA using single-molecule optical tweezers

Transmembrane transporter proteins are essential for maintaining cellular homeostasis and, as such, are key drug targets. Many transmembrane transporter proteins are known to undergo large structural rearrangements during their functional cycles. Despite the wealth of detailed structural and functional data available for these systems, our understanding of their dynamics and, consequently, how they function is generally limited. We introduce an innovative approach that enables us to directly measure the dynamics and stability of interdomain interactions of transmembrane proteins using optical tweezers. Focusing on the osmoregulatory ATP-binding cassette transporter OpuA from Lactococcus lactis, we examine the mechanical properties and potential interactions of its substrate-binding domains. Our measurements are performed in lipid nanodiscs, providing a native-mimicking environment for the transmembrane protein. The technique provides high spatial and temporal resolution and allows us to study the functionally relevant motions and interdomain interactions of individual transmembrane transporter proteins in real time in a lipid bilayer.

Inhibitor screening for volume-sensitive LRRC8A chloride channel

Leucine-rich repeat-containing protein 8 A (LRRC8A) protein is a critical member of volume-regulated anion channels. It plays a critical roles in the regulation of cellular volume and involves in the development of diseases like osteoarthritis. Screening of lead compounds to modulate its function may provide potential therapeutics of related diseases. Here, we employ virtual screening techniques and molecular dynamics (MD) simulation to screen potential inhibitors against LRRC8A. LRRC8A was regarded as the drug target to investigate potential compounds from the ZINC15 database via molecular docking. The final compound was selected among the top 10 Autodock Vina score (-8.8 Kcal/mol) with the ZINC ID ZINC000018195627 after druggability prediction. The docked complex from the virtual screening was subjected to MD simulation to analyze the stability of the LRRC8A protein-ligand complex, with parameters including root mean square deviation, root mean square fluctuation and radius of gyration. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method was further employed to predict the binding free energies from MD simulation trajectory. Our study provides insightful analysis for the potential compound to modulate LRRC8A and lay the foundation of therapeutics development against osteoarthritis.Communicated by Ramaswamy H. Sarma.

Tuning Protein Hydrogel Mechanics through Modulation of Nanoscale Unfolding and Entanglement in Postgelation Relaxation

Globular folded proteins are versatile nanoscale building blocks to create biomaterials with mechanical robustness and inherent biological functionality due to their specific and well-defined folded structures. Modulating the nanoscale unfolding of protein building blocks during network formation (in situ protein unfolding) provides potent opportunities to control the protein network structure and mechanics. Here, we control protein unfolding during the formation of hydrogels constructed from chemically cross-linked maltose binding protein using ligand binding and the addition of cosolutes to modulate protein kinetic and thermodynamic stability. Bulk shear rheology characterizes the storage moduli of the bound and unbound protein hydrogels and reveals a correlation between network rigidity, characterized as an increase in the storage modulus, and protein thermodynamic stability. Furthermore, analysis of the network relaxation behavior identifies a crossover from an unfolding dominated regime to an entanglement dominated regime. Control of in situ protein unfolding and entanglement provides an important route to finely tune the architecture, mechanics, and dynamic relaxation of protein hydrogels. Such predictive control will be advantageous for future smart biomaterials for applications which require responsive and dynamic modulation of mechanical properties and biological function.

Control of Nanoscale In Situ Protein Unfolding Defines Network Architecture and Mechanics of Protein Hydrogels

Hierarchical assemblies of proteins exhibit a wide-range of material properties that are exploited both in nature and by artificially by humankind. However, little is understood about the importance of protein unfolding on the network assembly, severely limiting opportunities to utilize this nanoscale transition in the development of biomimetic and bioinspired materials. Here we control the force lability of a single protein building block, bovine serum albumin (BSA), and demonstrate that protein unfolding plays a critical role in defining the architecture and mechanics of a photochemically cross-linked native protein network. The internal nanoscale structure of BSA contains "molecular reinforcement" in the form of 17 covalent disulphide "nanostaples", preventing force-induced unfolding. Upon addition of reducing agents, these nanostaples are broken rendering the protein force labile. Employing a combination of circular dichroism (CD) spectroscopy, small-angle scattering (SAS), rheology, and modeling, we show that stapled protein forms reasonably homogeneous networks of cross-linked fractal-like clusters connected by an intercluster region of folded protein. Conversely, in situ protein unfolding results in more heterogeneous networks of denser fractal-like clusters connected by an intercluster region populated by unfolded protein. In addition, gelation-induced protein unfolding and cross-linking in the intercluster region changes the hydrogel mechanics, as measured by a 3-fold enhancement of the storage modulus, an increase in both the loss ratio and energy dissipation, and markedly different relaxation behavior. By controlling the protein's ability to unfold through nanoscale (un)stapling, we demonstrate the importance of in situ unfolding in defining both network architecture and mechanics, providing insight into fundamental hierarchical mechanics and a route to tune biomaterials for future applications.

From folding to function: complex macromolecular reactions unraveled one-by-one with optical tweezers

Single-molecule manipulation with optical tweezers has uncovered macromolecular behaviour hidden to other experimental techniques. Recent instrumental improvements have made it possible to expand the range of systems accessible to optical tweezers. Beyond focusing on the folding and structural changes of isolated single molecules, optical tweezers studies have evolved into unraveling the basic principles of complex molecular processes such as co-translational folding on the ribosome, kinase activation dynamics, ligand-receptor binding, chaperone-assisted protein folding, and even dynamics of intrinsically disordered proteins (IDPs). In this mini-review, we illustrate the methodological principles of optical tweezers before highlighting recent advances in studying complex protein conformational dynamics - from protein synthesis to physiological function - as well as emerging future issues that are beginning to be addressed with novel approaches.

Role of Ligand Binding Site in Modulating the Mechanical Stability of Proteins with β-Grasp Fold

Despite many studies on ligand-modulated protein mechanics, a comparative analysis of the role of ligand binding site on any specific protein fold is yet to be made. In this study, we explore the role of ligand binding site on the mechanical properties of β-grasp fold proteins, namely, ubiquitin and small ubiquitin related modifier 1 (SUMO1). The terminal segments directly connected through hydrogen bonds constitute the β-clamp geometry (or mechanical clamp), which confers high mechanical resilience to the β-grasp fold. Here, we study ubiquitin complexed with CUE2-1, a ubiquitin-binding domain (UBD) from yeast endonuclease protein Cue2, using a combination of single-molecule force spectroscopy (SMFS) and steered molecular dynamics (SMD) simulations. Our study reveals that CUE2-1 does not alter the mechanical properties of ubiquitin, despite directly interacting with its β-clamp. To explore the role of ligand binding site, we compare the mechanical properties of the ubiquitin/CUE2-1 complex with that of previously studied SUMO1/S12, another β-grasp protein complex, using SMD simulations. Simulations on the SUMO1/S12 complex corroborate previous experimentally observed enhancement in the mechanical stability of SUMO1, even though S12 binds away from the β-clamp. Differences in ligand binding-induced structural impact at the transition state of the two complexes explain the differences in ligand modulated protein mechanics. Contrary to previous reports, our study demonstrates that direct binding of ligands to the mechanical clamp does not necessarily alter the mechanical stability of β-grasp fold proteins. Rather, binding interactions away from the clamp can reinforce protein stability provided by the β-grasp fold. Our study highlights the importance of binding site and binding modes of ligands in modulating the mechanical stability of β-grasp fold proteins.

Effects of Ligand Binding on the Energy Landscape of Acyl-CoA-Binding Protein

Binding of ligands is often crucial for function yet the effects of ligand binding on the mechanical stability and energy landscape of proteins are incompletely understood. Here, we use a combination of single-molecule optical tweezers and MD simulations to investigate the effect of ligand binding on the energy landscape of acyl-coenzyme A (CoA)-binding protein (ACBP). ACBP is a topologically simple and highly conserved four-α-helix bundle protein that acts as an intracellular transporter and buffer for fatty-acyl-CoA and is active in membrane assembly. We have previously described the behavior of ACBP under tension, revealing a highly extended transition state (TS) located almost halfway between the unfolded and native states. Here, we performed force-ramp and force-jump experiments, in combination with advanced statistical analysis, to show that octanoyl-CoA binding increases the activation free energy for the unfolding reaction of ACBP without affecting the position of the transition state along the reaction coordinate. It follows that ligand binding enhances the mechanical resistance and thermodynamic stability of the protein, without changing its mechanical compliance. Steered molecular dynamics simulations allowed us to rationalize the results in terms of key interactions that octanoyl-CoA establishes with the four α-helices of ACBP and showed that the unfolding pathway is marginally affected by the ligand. The results show that ligand-induced mechanical stabilization effects can be complex and may prove useful for the rational design of stabilizing ligands.

Single molecule protein stabilisation translates to macromolecular mechanics of a protein network

Folded globular proteins are attractive building blocks for biopolymer-based materials, as their mechanically resistant structures carry out diverse biological functionality. While much is now understood about the mechanical response of single folded proteins, a major challenge is to understand and predictably control how single protein mechanics translates to the collective response of a network of connected folded proteins. Here, by utilising the binding of maltose to hydrogels constructed from photo-chemically cross-linked maltose binding protein (MBP), we investigate the effects of protein stabilisation at the molecular level on the macroscopic mechanical and structural properties of a protein-based hydrogel. Rheological measurements show an enhancement in the mechanical strength and energy dissipation of MBP hydrogels in the presence of maltose. Circular dichroism spectroscopy and differential scanning calorimetry measurements show that MBP remains both folded and functional in situ. By coupling these mechanical measurements with mesoscopic structural information obtained by small angle scattering, we propose an occupation model in which higher proportions of stabilised, ligand occupied, protein building blocks translate their increased stability to the macroscopic properties of the hydrogel network. This provides powerful opportunities to exploit environmentally responsive folded protein-based biomaterials for many broad applications.

Unraveling the Mechanical Unfolding Pathways of a Multidomain Protein: Phosphoglycerate Kinase

Phosphoglycerate kinase (PGK) is a highly conserved enzyme that is crucial for glycolysis. PGK is a monomeric protein composed of two similar domains and has been the focus of many studies for investigating interdomain interactions within the native state and during folding. Previous studies used traditional biophysical methods (such as circular dichroism, tryptophan fluorescence, and NMR) to measure signals over a large ensemble of molecules, which made it difficult to observe transient changes in stability or structure during unfolding and refolding of single molecules. Here, we unfold single molecules of PGK using atomic force spectroscopy and steered molecular dynamic computer simulations to examine the conformational dynamics of PGK during its unfolding process. Our results show that after the initial forced separation of its domains, yeast PGK (yPGK) does not follow a single mechanical unfolding pathway; instead, it stochastically follows two distinct pathways: unfolding from the N-terminal domain or unfolding from the C-terminal domain. The truncated yPGK N-terminal domain unfolds via a transient intermediate, whereas the structurally similar isolated C-terminal domain has no detectable intermediates throughout its mechanical unfolding process. The N-terminal domain in the full-length yPGK displays a strong unfolding intermediate 13% of the time, whereas the truncated domain (yPGKNT) transitions through the intermediate 81% of the time. This effect indicates that the mechanical properties of yPGK cannot be simply deduced from the mechanical properties of its constituents. We also find that Escherichia coli PGK is significantly less mechanically stable as compared to yPGK, contrary to bulk unfolding measurements. Our results support the growing body of observations that the folding behavior of multidomain proteins is difficult to predict based solely on the studies of isolated domains.

An Active, Ligand-Responsive Pulling Geometry Reports on Internal Signaling between Subdomains of the DnaK Nucleotide-Binding Domain in Single-Molecule Mechanical Experiments

Single-molecule mechanical experiments have proven to be ideal tools for probing the energetics and mechanics of large proteins and domains. In this paper, we investigate the nucleotide-dependent unfolding mechanics of the nucleotide-binding domain (NBD) of the Hsp70 chaperone DnaK. The NBD binds ADP or ATP in the binding cleft formed by lobe I and lobe II, which consists of two subdomains each. When force is applied to the termini of the NBD, the observed unfolding forces are independent of the nucleotide state. In contrast, when force is applied across the nucleotide-binding pocket, the unfolding forces report specifically on the nucleotide-phosphate state. In this active, ligand-responsive pulling geometry, we observed a bifurcation of the unfolding pathway; the pathway proceeds either through a cooperative "coupled pathway" or through a noncooperative "uncoupled pathway". The partitioning between individual unfolding pathways can be effectively tuned by mutation or by the nucleotide exchange factor GrpE, i.e., by the factors affecting the strength of the lobe I-lobe II interactions within the native NBD. These experiments provide important insight into the molecular origin of the internal signaling between the subdomains of the nucleotide-binding domain of Hsp70 proteins and how signals are efficiently transferred inside the protein molecule.

DNA Binding Induces a Nanomechanical Switch in the RRM1 Domain of TDP-43

Understanding the molecular mechanisms governing protein-nucleic acid interactions is fundamental to many nuclear processes. However, how nucleic acid binding affects the conformation and dynamics of the substrate protein remains poorly understood. Here we use a combination of single molecule force spectroscopy AFM and biochemical assays to show that the binding of TG-rich ssDNA triggers a mechanical switch in the RRM1 domain of TDP-43, toggling between an entropic spring devoid of mechanical stability and a shock absorber bound-form that resists unfolding forces of ∼40 pN. The fraction of mechanically resistant proteins correlates with an increasing length of the TG _n oligonucleotide, demonstrating that protein mechanical stability is a direct reporter of nucleic acid binding. Steered molecular dynamics simulations on related RNA oligonucleotides reveal that the increased mechanical stability fingerprinting the holo-form is likely to stem from a unique scenario whereby the nucleic acid acts as a "mechanical staple" that protects RRM1 from mechanical unfolding. Our approach highlights nucleic acid binding as an effective strategy to control protein nanomechanics.

The life of proteins under mechanical force

Although much of our understanding of protein folding comes from studies of isolated protein domains in bulk, in the cellular environment the intervention of external molecular machines is essential during the protein life cycle. During the past decade single molecule force spectroscopy techniques have been extremely useful to deepen our understanding of these interventional molecular processes, as they allow for monitoring and manipulating mechanochemical events in individual protein molecules. Here, we review some of the critical steps in the protein life cycle, starting with the biosynthesis of the nascent polypeptide chain in the ribosome, continuing with the folding supported by chaperones and the translocation into different cell compartments, and ending with proteolysis in the proteasome. Along these steps, proteins experience molecular forces often combined with chemical transformations, affecting their folding and structure, which are measured or mimicked in the laboratory by the application of force with a single molecule apparatus. These mechanochemical reactions can potentially be used as targets for fighting against diseases. Inspired by these insightful experiments, we devise an outlook on the emerging field of mechanopharmacology, which reflects an alternative paradigm for drug design.

Resolving dual binding conformations of cellulosome cohesin-dockerin complexes using single-molecule force spectroscopy

Receptor-ligand pairs are ordinarily thought to interact through a lock and key mechanism, where a unique molecular conformation is formed upon binding. Contrary to this paradigm, cellulosomal cohesin-dockerin (Coh-Doc) pairs are believed to interact through redundant dual binding modes consisting of two distinct conformations. Here, we combined site-directed mutagenesis and single-molecule force spectroscopy (SMFS) to study the unbinding of Coh:Doc complexes under force. We designed Doc mutations to knock out each binding mode, and compared their single-molecule unfolding patterns as they were dissociated from Coh using an atomic force microscope (AFM) cantilever. Although average bulk measurements were unable to resolve the differences in Doc binding modes due to the similarity of the interactions, with a single-molecule method we were able to discriminate the two modes based on distinct differences in their mechanical properties. We conclude that under native conditions wild-type Doc from Clostridium thermocellum exocellulase Cel48S populates both binding modes with similar probabilities. Given the vast number of Doc domains with predicted dual binding modes across multiple bacterial species, our approach opens up new possibilities for understanding assembly and catalytic properties of a broad range of multi-enzyme complexes.

Small peptide binding stiffens the ubiquitin-like protein SUMO1

Posttranslational modification by small ubiquitin-like modifiers (SUMOs), known as SUMOylation, is a key regulatory event in many eukaryotic cellular processes in which SUMOs interact with a large number of target proteins. SUMO binding motifs (SBMs) are small peptides derived from these target proteins that interact noncovalently with SUMOs and induce conformational changes. To determine the effect of SBMs on the mechanical properties of SUMO1 (the first member of the human SUMO family), we performed single-molecule force spectroscopy experiments on SUMO1/SBM complexes. The unfolding force of SUMO1 (at a pulling speed of 400 nm/s) increased from ∼ 130 pN to ∼ 170 pN upon binding to SBMs, indicating mechanical stabilization upon complexation. Pulling-speed-dependent experiments and Monte Carlo simulations measured a large decrease in distance to the unfolding transition state for SUMO1 upon SBM binding, which is by far the largest change measured for any ligand binding protein. The stiffness of SUMO1 (measured as a spring constant for the deformation response along the line joining the N- and C-termini) increased upon SBM binding from ∼ 1 N/m to ∼ 3.5 N/m. The relatively higher flexibility of ligand-free SUMO1 might play a role in accessing various conformations before binding to a target.

Rapid and Robust Polyprotein Production Facilitates Single-Molecule Mechanical Characterization of β-Barrel Assembly Machinery Polypeptide Transport Associated Domains

Single-molecule force spectroscopy by atomic force microscopy exploits the use of multimeric protein constructs, namely, polyproteins, to decrease the impact of nonspecific interactions, to improve data accumulation, and to allow the accommodation of benchmarking reference domains within the construct. However, methods to generate such constructs are either time- and labor-intensive or lack control over the length or the domain sequence of the obtained construct. Here, we describe an approach that addresses both of these shortcomings that uses Gibson assembly (GA) to generate a defined recombinant polyprotein rapidly using linker sequences. To demonstrate the feasibility of this approach, we used GA to make a polyprotein composed of alternating domains of I27 and TmCsp, (I27-TmCsp)3-I27)(GA), and showed the mechanical fingerprint, mechanical strength, and pulling speed dependence are the same as an analogous polyprotein constructed using the classical approach. After this benchmarking, we exploited this approach to facilitiate the mechanical characterization of POTRA domain 2 of BamA from E. coli (EcPOTRA2) by assembling the polyprotein (I27-EcPOTRA2)3-I27(GA). We show that, as predicted from the α + β topology, EcPOTRA2 domains are mechanically robust over a wide range of pulling speeds. Furthermore, we identify a clear correlation between mechanical robustness and brittleness for a range of other α + β proteins that contain the structural feature of proximal terminal β-strands in parallel geometry. We thus demonstrate that the GA approach is a powerful tool, as it circumvents the usual time- and labor-intensive polyprotein production process and allows for rapid production of new constructs for single-molecule studies. As shown for EcPOTRA2, this approach allows the exploration of the mechanical properties of a greater number of proteins and their variants. This improves our understanding of the relationship between structure and mechanical strength, increasing our ability to design proteins with tailored mechanical properties.

The Mechanochemistry of a Structural Zinc Finger

Zinc fingers are highly ubiquitous structural motifs that provide stability to proteins, thus contributing to their correct folding. Despite the high thermodynamic stability of the ZnCys4 centers, their kinetic properties display remarkable lability. Here, we use a combination of protein engineering with single molecule force spectroscopy atomic force microscopy (AFM) to uncover the surprising mechanical lability (∼90 pN) of the individual Zn-S bonds that form the two equivalent zinc finger motifs embedded in the structure of the multidomain DnaJ chaperone. Rational mutations within the zinc coordinating residues enable direct identification of the chemical determinants that regulate the interplay between zinc binding-requiring the presence of all four cysteines-and disulfide bond formation. Finally, our observations show that binding to hydrophobic short peptides drastically increases the mechanical stability of DnaJ. Altogether, our experimental approach offers a detailed, atomistic vista on the fine chemical mechanisms that govern the nanomechanics of individual, naturally occurring zinc finger.

Identifying sequential substrate binding at the single-molecule level by enzyme mechanical stabilization

Enzyme-substrate binding is a dynamic process intimately coupled to protein structural changes, which in turn changes the unfolding energy landscape. By the use of single-molecule force spectroscopy (SMFS), we characterize the open-to-closed conformational transition experienced by the hyperthermophilic adenine diphosphate (ADP)-dependent glucokinase from Thermococcus litoralis triggered by the sequential binding of substrates. In the absence of substrates, the mechanical unfolding of TlGK shows an intermediate 1, which is stabilized in the presence of Mg·ADP(-), the first substrate to bind to the enzyme. However, in the presence of this substrate, an additional unfolding event is observed, intermediate 1*. Finally, in the presence of both substrates, the unfolding force of intermediates 1 and 1* increases as a consequence of the domain closure. These results show that SMFS can be used as a powerful experimental tool to investigate binding mechanisms of different enzymes with more than one ligand, expanding the repertoire of protocols traditionally used in enzymology.

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.

Mechanical unfolding of ribose binding protein and its comparison with other periplasmic binding proteins

Folding and unfolding studies on large, multidomain proteins are still rare despite their high abundance in genomes of prokaryotes and eukaryotes. Here, we investigate the unfolding properties of a 271 residue, two-domain ribose binding protein (RBP) from the bacterial periplasm using single-molecule force spectroscopy. We observe that RBP predominately unfolds via a two-state pathway with an unfolding force of ∼80 pN and an unfolding contour length of ∼95 nm. Only a small population (∼15%) of RBP follows three-state pathways. The ligand binding neither increases the mechanical stability nor influences the unfolding flux of RBP through different pathways. The kinetic partitioning between two-state and three-state pathways, which has been reported earlier for other periplasmic proteins, is also observed in RBP, albeit to a lesser extent. These results provide important insights into the mechanical stability and unfolding processes of large two-domain proteins and highlight the contrasting features upon ligand binding. Protein structural topology diagrams are used to explain the differences in the mechanical unfolding behavior of RBP with other periplasmic binding proteins.

Force spectroscopy studies on protein-ligand interactions: a single protein mechanics perspective

Protein-ligand interactions are ubiquitous and play important roles in almost every biological process. The direct elucidation of the thermodynamic, structural and functional consequences of protein-ligand interactions is thus of critical importance to decipher the mechanism underlying these biological processes. A toolbox containing a variety of powerful techniques has been developed to quantitatively study protein-ligand interactions in vitro as well as in living systems. The development of atomic force microscopy-based single molecule force spectroscopy techniques has expanded this toolbox and made it possible to directly probe the mechanical consequence of ligand binding on proteins. Many recent experiments have revealed how ligand binding affects the mechanical stability and mechanical unfolding dynamics of proteins, and provided mechanistic understanding on these effects. The enhancement effect of mechanical stability by ligand binding has been used to help tune the mechanical stability of proteins in a rational manner and develop novel functional binding assays for protein-ligand interactions. Single molecule force spectroscopy studies have started to shed new lights on the structural and functional consequence of ligand binding on proteins that bear force under their biological settings.

Ca2+ binding enhanced mechanical stability of an archaeal crystallin

Structural topology plays an important role in protein mechanical stability. Proteins with β-sandwich topology consisting of Greek key structural motifs, for example, I27 of muscle titin and ¹⁰FNIII of fibronectin, are mechanically resistant as shown by single-molecule force spectroscopy (SMFS). In proteins with β-sandwich topology, if the terminal strands are directly connected by backbone H-bonding then this geometry can serve as a “mechanical clamp”. Proteins with this geometry are shown to have very high unfolding forces. Here, we set out to explore the mechanical properties of a protein, M-crystallin, which belongs to β-sandwich topology consisting of Greek key motifs but its overall structure lacks the “mechanical clamp” geometry at the termini. M-crystallin is a Ca²⁺ binding protein from Methanosarcina acetivorans that is evolutionarily related to the vertebrate eye lens β and γ-crystallins. We constructed an octamer of crystallin, (M-crystallin)₈, and using SMFS, we show that M-crystallin unfolds in a two-state manner with an unfolding force ∼90 pN (at a pulling speed of 1000 nm/sec), which is much lower than that of I27. Our study highlights that the β-sandwich topology proteins with a different strand-connectivity than that of I27 and ¹⁰FNIII, as well as lacking “mechanical clamp” geometry, can be mechanically resistant. Furthermore, Ca²⁺ binding not only stabilizes M-crystallin by 11.4 kcal/mol but also increases its unfolding force by ∼35 pN at the same pulling speed. The differences in the mechanical properties of apo and holo M-crystallins are further characterized using pulling speed dependent measurements and they show that Ca²⁺ binding reduces the unfolding potential width from 0.55 nm to 0.38 nm. These results are explained using a simple two-state unfolding energy landscape.

Analysis of the REJ Module of Polycystin-1 Using Molecular Modeling and Force-Spectroscopy Techniques

Polycystin-1 is a large transmembrane protein, which, when mutated, causes autosomal dominant polycystic kidney disease, one of the most common life-threatening genetic diseases that is a leading cause of kidney failure. The REJ (receptor for egg lelly) module is a major component of PC1 ectodomain that extends to about 1000 amino acids. Many missense disease-causing mutations map to this module; however, very little is known about the structure or function of this region. We used a combination of homology molecular modeling, protein engineering, steered molecular dynamics (SMD) simulations, and single-molecule force spectroscopy (SMFS) to analyze the conformation and mechanical stability of the first ~420 amino acids of REJ. Homology molecular modeling analysis revealed that this region may contain structural elements that have an FNIII-like structure, which we named REJd1, REJd2, REJd3, and REJd4. We found that REJd1 has a higher mechanical stability than REJd2 (~190 pN and 60 pN, resp.). Our data suggest that the putative domains REJd3 and REJd4 likely do not form mechanically stable folds. Our experimental approach opens a new way to systematically study the effects of disease-causing mutations on the structure and mechanical properties of the REJ module of PC1.

Effects of ligand binding on the mechanical properties of ankyrin repeat protein gankyrin

Ankyrin repeat proteins are elastic materials that unfold and refold sequentially, repeat by repeat, under force. Herein we use atomistic molecular dynamics to compare the mechanical properties of the 7-ankyrin-repeat oncoprotein Gankyrin in isolation and in complex with its binding partner S6-C. We show that the bound S6-C greatly increases the resistance of Gankyrin to mechanical stress. The effect is specific to those repeats of Gankyrin directly in contact with S6-C, and the mechanical ‘hot spots’ of the interaction map to the same repeats as the thermodynamic hot spots. A consequence of stepwise nature of unfolding and the localized nature of ligand binding is that it impacts on all aspects of the protein's mechanical behavior, including the order of repeat unfolding, the diversity of unfolding pathways accessed, the nature of partially unfolded intermediates, the forces required and the work transferred to the system to unfold the whole protein and its parts. Stepwise unfolding thus provides the means to buffer repeat proteins and their binding partners from mechanical stress in the cell. Our results illustrate how ligand binding can control the mechanical response of proteins. The data also point to a cellular mechano-switching mechanism whereby binding between two partner macromolecules is regulated by mechanical stress.

Here we use molecular dynamics simulation to compare the mechanical properties of the 7-ankyrin-repeat oncoprotein Gankyrin in isolation and in complex with binding partner S6-C. Tandem repeat proteins like Gankyrin comprise tandem arrays of small structural motifs that pack linearly to produce elongated architectures. They are elastic, mechanically weak molecules and they unfold and refold repeat by repeat under force. We show that S6-C binding greatly increases the resistance of Gankyrin to mechanical stress. The enhanced mechanical stability is specific to those ankyrin repeats in contact with S6-C, and the localized nature of the effect results in fundamental changes in the way the protein responds to force. Thus, the forced unfolding of isolated Gankryin involves a diverse set of pathways with a preference for a C- to N-terminus unfolding mechanism whereas this diversity is reduced upon complex formation with the central repeats, which are those most tightly bound to the ligand, tending to unfold last. Our study shows how stepwise unfolding can buffer repeat proteins and their binding partners from mechanical stress in the cell. It also points to a mechano-switching mechanism whereby binding between two partner macromolecules is regulated by mechanical stress.

Single-molecule experiments reveal the flexibility of a Per-ARNT-Sim domain and the kinetic partitioning in the unfolding pathway under force

Per-ARNT-Sim (PAS) domains serve as versatile binding motifs in many signal-transduction proteins and are able to respond to a wide spectrum of chemical or physical signals. Despite their diverse functions, PAS domains share a conserved structure. It has been suggested that the structure of PAS domains is flexible and thus adaptable to many binding partners. However, direct measurement of the flexibility of PAS domains has not yet been provided. Here, we quantitatively measure the mechanical unfolding of a PAS domain, ARNT PAS-B, using single-molecule atomic force microscopy. Our force spectroscopy results indicate that the structure of ARNT PAS-B can be unraveled under mechanical forces as low as ~30 pN due to its broad potential well for the mechanical unfolding transition of ~2 nm. This allows the PAS-B domain to extend by up to 75% of its resting end-to-end distance without unfolding. Moreover, we found that the ARNT PAS-B domain unfolds in two distinct pathways via a kinetic partitioning mechanism. Sixty-seven percent of ARNT PAS-B unfolds through a simple two-state pathway, whereas the other 33% unfolds with a well-defined intermediate state in which the C-terminal β-hairpin is detached. We propose that the structural flexibility and force-induced partial unfolding of PAS-B domains may provide a unique mechanism for them to recruit diverse binding partners and lower the free-energy barrier for the formation of the binding interface.

Variation in the mechanical unfolding pathway of p53DBD induced by interaction with p53 N-terminal region or DNA

The tumor suppressor p53 plays a crucial role in the cell cycle checkpoints, DNA repair, and apoptosis. p53 consists of a natively unfolded N-terminal region (NTR), central DNA binding domain (DBD), C-terminal tetramerization domain, and regulatory region. In this paper, the interactions between the DBD and the NTR, and between the DBD and DNA were investigated by measuring changes in the mechanical unfolding trajectory of the DBD using atomic force microscopy (AFM)-based single molecule force spectroscopy. In the absence of DNA, the DBD (94–293, 200 amino acids (AA)) showed two different mechanical unfolding patterns. One indicated the existence of an unfolding intermediate consisting of approximately 60 AA, and the other showed a 100 AA intermediate. The DBD with the NTR did not show such unfolding patterns, but heterogeneous unfolding force peaks were observed. Of the heterogeneous patterns, we observed a high frequency of force peaks indicating the unfolding of a domain consisting of 220 AA, which is apparently larger than that of a sole DBD. This observation implies that a part of NTR binds to the DBD, and the mechanical unfolding happens not solely on the DBD but also accompanying a part of NTR. When DNA is bound, the mechanical unfolding trajectory of p53NTR+DBD showed a different pattern from that without DNA. The pattern was similar to that of the DBD alone, but two consecutive unfolding force peaks corresponding to 60 and 100 AA sub-domains were observed. These results indicate that interactions with the NTR or DNA alter the mechanical stability of DBD and result in drastic changes in the mechanical unfolding trajectory of the DBD.

Prying open single GroES ring complexes by force reveals cooperativity across domains

Understanding how the mechanical properties of a protein complex emerge from the interplay of intra- and interchain interactions is vital at both fundamental and applied levels. To investigate whether interdomain cooperativity affects protein mechanical strength, we employed single-molecule force spectroscopy to probe the mechanical stability of GroES, a homoheptamer with a domelike quaternary stucture stabilized by intersubunit interactions between the first and last β-strands of adjacent domains. A GroES variant was constructed in which each subunit of the GroES heptamer is covalently linked to adjacent subunits by tripeptide linkers and folded domains of protein L are introduced to the heptamer's termini as handle molecules. The force-distance profiles for GroES unfolding showed, for the first time that we know of, a mechanical phenotype whereby seven distinct force peaks, with alternating behavior of unfolding force and contour length (ΔL(c)), were observed with increasing unfolding-event number. Unfolding of (GroES)(7) is initiated by breakage of the interface between domains 1 and 7 at low force, which imparts a polarity to (GroES)(7) that results in two distinct mechanical phenotypes of these otherwise identical protein domains. Unfolding then proceeds by peeling domains off the domelike native structure by sequential repetition of the denaturation of mechanically weak (unfoldon 1) and strong (unfoldon 2) units. These results indicate that domain-domain interactions help to determine the overall mechanical strength and unfolding pathway of the oligomeric structure. These data reveal an unexpected richness in the mechanical behavior of this homopolyprotein, yielding a complex with greater mechanical strength and properties distinct from those that would be apparent for GroES domains in isolation.

Mechanical anisotropy of ankyrin repeats

Red blood cells are frequently deformed and their cytoskeletal proteins such as spectrin and ankyrin-R are repeatedly subjected to mechanical forces. While the mechanics of spectrin was thoroughly investigated in vitro and in vivo, little is known about the mechanical behavior of ankyrin-R. In this study, we combine coarse-grained steered molecular dynamics simulations and atomic force spectroscopy to examine the mechanical response of ankyrin repeats (ARs) in a model synthetic AR protein NI6C, and in the D34 fragment of native ankyrin-R when these proteins are subjected to various stretching geometry conditions. Our steered molecular dynamics results, supported by AFM measurements, reveal an unusual mechanical anisotropy of ARs: their mechanical stability is greater when their unfolding is forced to propagate from the N-terminus toward the C-terminus (repeats unfold at ~60 pN), as compared to the unfolding in the opposite direction (unfolding force ∼ 30 pN). This anisotropy is also reflected in the complex refolding behavior of ARs. The origin of this unfolding and refolding anisotropy is in the various numbers of native contacts that are broken and formed at the interfaces between neighboring repeats depending on the unfolding/refolding propagation directions. Finally, we discuss how these complex mechanical properties of ARs in D34 may affect its behavior in vivo.

Dynamics of protein folding and cofactor binding monitored by single-molecule force spectroscopy

Many proteins in living cells require cofactors to carry out their biological functions. To reach their functional states, these proteins need to fold into their unique three-dimensional structures in the presence of their cofactors. Two processes, folding of the protein and binding of cofactors, intermingle with each other, making the direct elucidation of the folding mechanism of proteins in the presence of cofactors challenging. Here we use single-molecule atomic force microscopy to directly monitor the folding and cofactor binding dynamics of an engineered metal-binding protein G6-53 at the single-molecule level. Using the mechanical stability of different conformers of G6-53 as sensitive probes, we directly identified different G6-53 conformers (unfolded, apo- and Ni(2+)-bound) populated along the folding pathway of G6-53 in the presence of its cofactor Ni(2+). By carrying out single-molecule atomic force microscopy refolding experiments, we monitored kinetic evolution processes of these different conformers. Our results suggested that the majority of G6-53 folds through a binding-after-folding mechanism, whereas a small fraction follows a binding-before-folding pathway. Our study opens an avenue to utilizing force spectroscopy techniques to probe the folding dynamics of proteins in the presence of cofactors at the single-molecule level, and we anticipated that this method can be used to study a wide variety of proteins requiring cofactors for their function.

Ligand-modulated parallel mechanical unfolding pathways of maltose-binding proteins

Protein folding and unfolding are complex phenomena, and it is accepted that multidomain proteins generally follow multiple pathways. Maltose-binding protein (MBP) is a large (a two-domain, 370-amino acid residue) bacterial periplasmic protein involved in maltose uptake. Despite the large size, it has been shown to exhibit an apparent two-state equilibrium unfolding in bulk experiments. Single-molecule studies can uncover rare events that are masked by averaging in bulk studies. Here, we use single-molecule force spectroscopy to study the mechanical unfolding pathways of MBP and its precursor protein (preMBP) in the presence and absence of ligands. Our results show that MBP exhibits kinetic partitioning on mechanical stretching and unfolds via two parallel pathways: one of them involves a mechanically stable intermediate (path I) whereas the other is devoid of it (path II). The apoMBP unfolds via path I in 62% of the mechanical unfolding events, and the remaining 38% follow path II. In the case of maltose-bound MBP, the protein unfolds via the intermediate in 79% of the cases, the remaining 21% via path II. Similarly, on binding to maltotriose, a ligand whose binding strength with the polyprotein is similar to that of maltose, the occurrence of the intermediate is comparable (82% via path I) with that of maltose. The precursor protein preMBP also shows a similar behavior upon mechanical unfolding. The percentages of molecules unfolding via path I are 53% in the apo form and 68% and 72% upon binding to maltose and maltotriose, respectively, for preMBP. These observations demonstrate that ligand binding can modulate the mechanical unfolding pathways of proteins by a kinetic partitioning mechanism. This could be a general mechanism in the unfolding of other large two-domain ligand-binding proteins of the bacterial periplasmic space.

Inhibitor binding increases the mechanical stability of staphylococcal nuclease

Staphylococcal nuclease (SNase) catalyzes the hydrolysis of DNA and RNA in a calcium-dependent fashion. We used AFM-based single-molecule force spectroscopy to investigate the mechanical stability of SNase alone and in its complex with an SNase inhibitor, deoxythymidine 3',5'-bisphosphate. We found that the enzyme unfolds in an all-or-none fashion at ∼26 pN. Upon binding to the inhibitor, the mechanical unfolding forces of the enzyme-inhibitor complex increase to ∼50 pN. This inhibitor-induced increase in the mechanical stability of the enzyme is consistent with the increased thermodynamical stability of the complex over that of SNase. Because of its strong mechanical response to inhibitor binding, SNase, a model protein folding system, offers a unique opportunity for studying the relationship between enzyme mechanics and catalysis.

Single-molecule force-clamp spectroscopy: dwell time analysis and practical considerations

Single-molecule force-clamp spectroscopy has become a powerful tool for studying protein folding/unfolding, bond rupture, and enzymatic reactions. Different methods have been developed to analyze force-clamp spectroscopy data on polyproteins to obtain kinetic parameters characterizing the mechanical unfolding of proteins, which are often modeled as a two-state process (a Poisson process). However, because of the finite number of domains in polyproteins, the statistical analysis of the force-clamp spectroscopy data is different from that of a classical Poisson process, and the equivalency of different analysis methods remains to be proven. In this article, we show that these methods are equivalent and lead to accurate measurements of the unfolding rate constant. We also demonstrate that distinct from the constant-pulling-velocity experiments, in which the unfolding rate extracted from the data is dependent on the number of protein domains in the polyproteins (the N effect), force-clamp experiments do not show any N effect. Using a simulated data set, we also highlighted important practical considerations that one needs to take into account when using the single-molecule force-clamp spectroscopy technique to characterize the unfolding energy landscape of proteins.

Unravelling the design principles for single protein mechanical strength

In recent years single molecule manipulation techniques have improved to the extent that measurements of the mechanical strength of single proteins can now be undertaken routinely. This powerful new tool, coupled with theoretical frameworks to characterise the unfolding process, has enabled significant progress to be made in understanding the physical mechanisms that underlie protein mechanical strength. These design concepts have allowed the search for proteins with novel, mechanically strong folds to be automated and for previously mechanically characterised proteins to be engineered rationally. Methods to achieve the latter are diverse and include re-engineering of specific hydrophobic core residues, changing solvent conditions and the 'cross-linking' of side-chains that are separated in the rate-limiting unfolding transition. Predicting the mechanical behaviour of larger proteins and those with more complex structures remains a significant challenge while on-going instrument development is beginning to allow the examination of mechanical strength of protein across a wide range of force loading rates. The integral role of force in biology and the potential for exploitation of catalytic and structural proteins as functional bio-materials makes this a particularly important area of research.