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

Cytosolic Interactome Protects Against Protein Unfolding in a Single Molecule Experiment

Single molecule techniques are particularly well suited for investigating the processes of protein folding and chaperone assistance. However, current assays provide only a limited perspective on the various ways in which the cellular environment can influence the folding pathway of a protein. In this study, a single molecule mechanical interrogation assay is developed and used to monitor protein unfolding and refolding within a cytosolic solution. This allows to test the cumulative topological effect of the cytoplasmic interactome on the folding process. The results reveal a stabilization against forced unfolding for partial folds, which are attributed to the protective effect of the cytoplasmic environment against unfolding and aggregation. This research opens the possibility of conducting single molecule molecular folding experiments in quasi-biological environments.

Retracing the evolution of a modern periplasmic binding protein

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

Native Salt Bridges Are a Key Regulator of Ubiquitin's Mechanical Stability

Although it is known that various intramolecular interactions determine protein mechanical stability, a detailed molecular-level understanding of the key regulators of protein mechanical stability is still lacking. Here, we present evidence for salt bridges in ubiquitin as important intramolecular interactions that can affect protein mechanical stability. Ubiquitin has two salt bridges: one relatively surface-exposed (SB1:K11-E34) and the other relatively buried (SB2:K27-D52). Ubiquitin is a reversible post-translational modifier and is stable mechanically (F_avg^u = 185 pN). On breaking SB1, the mechanical stability of ubiquitin is slightly enhanced (F_avg^u = 193 pN). In contrast, the mechanical stability significantly decreased upon breaking SB2 (F_avg^u = 158 pN). These results suggest that SB1 are SB2 are regulators of the mechanical stability of ubiquitin. Interestingly, the mechanical stability decreased further (F_avg^u = 145 pN) for the double salt bridge (DB) null variant. Monte Carlo simulations elucidate that the main regulating factor is the spontaneous unfolding rate constant (k_u⁰), being the highest for the DB null variant followed by the SB2 null variant, and it remains unaltered for the SB1 null variant, while the native-to-transition-state distance (x_u) remains unchanged. Our study provides mechanistic understanding on how two native salt bridges can independently regulate the mechanical stability in a protein, which has implications in designing protein-based robust biomaterials in the future.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

Pulling the springs of a cell by single-molecule force spectroscopy

The fundamental unit of the human body comprises of the cells which remain embedded in a fibrillar network of extracellular matrix proteins which in turn provides necessary anchorage the cells. Tissue repair, regeneration and reprogramming predominantly involve a traction force mediated signalling originating in the ECM and travelling deep into the cell including the nucleus via circuitry of spring-like filamentous proteins like microfilaments or actin, intermediate filaments and microtubules to elicit a response in the form of mechanical movement as well as biochemical changes. The 'springiness' of these proteins is highlighted in their extension-contraction behaviour which is manifested as an effect of differential traction force. Atomic force microscope (AFM) provides the magic eye to visualize and quantify such force-extension/indentation events in these filamentous proteins as well as in whole cells. In this review, we have presented a summary of the current understanding and advancement of such measurements by AFM based single-molecule force spectroscopy in the context of cytoskeletal and nucleoskeletal proteins which act in tandem to facilitate mechanotransduction.

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.

Structure dictates the mechanism of ligand recognition in the histidine and maltose binding proteins

Two mechanisms, induced fit (IF) and conformational selection (CS), have been proposed to explain ligand recognition coupled conformational changes. The histidine binding protein (HisJ) adopts the CS mechanism, in which a pre-equilibrium is established between the open and the closed states with the ligand binding to the closed state. Despite being structurally similar to HisJ, the maltose binding protein (MBP) adopts the IF mechanism, in which the ligand binds the open state and induces a transition to the closed state. To understand the molecular determinants of this difference, we performed molecular dynamics (MD) simulations of coarse-grained dual structure based models. We find that intra-protein contacts unique to the closed state are sufficient to promote the conformational transition in HisJ, indicating a CS-like mechanism. In contrast, additional ligand-mimicking contacts are required to “induce” the conformational transition in MBP suggesting an IF-like mechanism. In agreement with experiments, destabilizing modifications to two structural features, the spine helix (SH) and the balancing interface (BI), present in MBP but absent in HisJ, reduce the need for ligand-mimicking contacts indicating that SH and BI act as structural restraints that keep MBP in the open state. We introduce an SH like element into HisJ and observe that this can impede the conformational transition increasing the importance of ligand-mimicking contacts. Similarly, simultaneous mutations to BI and SH in MBP reduce the barrier to conformational transitions significantly and promote a CS-like mechanism. Together, our results show that structural restraints present in the protein structure can determine the mechanism of conformational transitions and even simple models that correctly capture such structural features can predict their positions. MD simulations of such models can thus be used, in conjunction with mutational experiments, to regulate protein ligand interactions, and modulate ligand binding affinities.

•

MBP operates by induced fit, HisJ by the conformational selection mechanism.

•

Dual structure based models (dSBMs) encode two structures of a protein.

•

MD simulations of dSBMs can identify the mechanism of conformational transitions.

•

Locks, absent in HisJ, hold MBP open with ligand contacts required for closing.

•

Binding mechanisms can be modified by altering such structural locks.

Next Generation Methods for Single-Molecule Force Spectroscopy on Polyproteins and Receptor-Ligand Complexes

Single-molecule force spectroscopy with the atomic force microscope provides molecular level insights into protein function, allowing researchers to reconstruct energy landscapes and understand functional mechanisms in biology. With steadily advancing methods, this technique has greatly accelerated our understanding of force transduction, mechanical deformation, and mechanostability within single- and multi-domain polyproteins, and receptor-ligand complexes. In this focused review, we summarize the state of the art in terms of methodology and highlight recent methodological improvements for AFM-SMFS experiments, including developments in surface chemistry, considerations for protein engineering, as well as theory and algorithms for data analysis. We hope that by condensing and disseminating these methods, they can assist the community in improving data yield, reliability, and throughput and thereby enhance the information that researchers can extract from such experiments. These leading edge methods for AFM-SMFS will serve as a groundwork for researchers cognizant of its current limitations who seek to improve the technique in the future for in-depth studies of molecular biomechanics.

Thermodynamics and folding landscapes of large proteins from a statistical mechanical model

Statistical mechanical models that afford an intermediate resolution between macroscopic chemical models and all-atom simulations have been successful in capturing folding behaviors of many small single-domain proteins. However, the applicability of one such successful approach, the Wako-Saitô-Muñoz-Eaton (WSME) model, is limited by the size of the protein as the number of conformations grows exponentially with protein length. In this work, we surmount this size limitation by introducing a novel approximation that treats stretches of 3 or 4 residues as blocks, thus reducing the phase space by nearly three orders of magnitude. The performance of the 'bWSME' model is validated by comparing the predictions for a globular enzyme (RNase H) and a repeat protein (IκBα), against experimental observables and the model without block approximation. Finally, as a proof of concept, we predict the free-energy surface of the 370-residue, multi-domain maltose binding protein and identify an intermediate in good agreement with single-molecule force-spectroscopy measurements. The bWSME model can thus be employed as a quantitative predictive tool to explore the conformational landscapes of large proteins, extract the structural features of putative intermediates, identify parallel folding paths, and thus aid in the interpretation of both ensemble and single-molecule experiments.

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.

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.

Deciphering the Mechanical Properties of Type III Secretion System EspA Protein by Single Molecule Force Spectroscopy

Bacterial pathogens inject virulence factors into host cells during bacterial infections using type III secretion systems. In enteropathogenic Escherichia coli, this system contains an external filament, formed by a self-oligomerizing protein called E. coli secreted protein A (EspA). The EspA filament penetrates the thick viscous mucus layer to facilitate the attachment of the bacteria to the gut-epithelium. To do that, the EspA filament requires noteworthy mechanical endurance considering the mechanical shear stresses found within the intestinal tract. To date, the mechanical properties of the EspA filament and the structural and biophysical knowledge of monomeric EspA are very limited, mostly due to the strong tendency of the protein to self-oligomerize. To overcome this limitation, we employed a single molecule force spectroscopy (SMFS) technique and studied the mechanical properties of EspA. Force extension dynamic of (I91)₄-EspA-(I91)₄ chimera revealed two structural unfolding events occurring at low forces during EspA unfolding, thus indicating no unique mechanical stability of the monomeric protein. SMFS examination of purified monomeric EspA protein, treated by a gradually refolding protocol, exhibited similar mechanical properties as the EspA protein within the (I91)₄-EspA-(I91)₄ chimera. Overall, our results suggest that the mechanical integrity of the EspA filament likely originates from the interactions between EspA monomers and not from the strength of an individual monomer.

Effect of ligand binding on a protein with a complex folding landscape

Ligand binding to a protein can stabilize it significantly against unfolding. The variation of the folding free energy, ΔΔG⁰, due to ligand binding can be derived from a simple reaction scheme involving exclusive binding to the native state. One obtains the following expression: , where K_d is the ligand dissociation constant and L is its concentration, R is the universal gas constant and T is the temperature. This expression has been shown to correctly describe experimental results on multiple proteins. In the current work we studied the effect of ligand binding on the stability of the multi-domain protein adenylate kinase from E. coli (AKE). Unfolding experiments were conducted using single-molecule FRET spectroscopy, which allowed us to directly obtain the fraction of unfolded protein in a model-free way from FRET efficiency histograms. Surprisingly, it was found that the effect of two inhibitors (Ap₅A and AMPPNP) and a substrate (AMP) on the stability of AKE was much smaller than expected based on K_d values obtained independently using microscale thermophoresis. To shed light on this issue, we measured the K_d for Ap₅A over a range of chemical denaturant concentrations where the protein is still folded. It was found that K_d increases dramatically over this range, likely due to the population of folding intermediates, whose binding to the ligand is much weaker than that of the native state. We propose that binding to folding intermediates may dominate the effect of ligands on the stability of multi-domain proteins, and could therefore have a strong impact on protein homeostasis in vivo.

Differences in the mechanical unfolding pathways of apo- and copper-bound azurins

Metalloproteins carry out diverse biological functions including metal transport, electron transfer, and catalysis. At present, the influence of metal cofactors on metalloprotein stability is not well understood. Here, we report the mechanical stability and unfolding pathway of azurin, a cupredoxin family protein with β-barrel topology and type I copper-binding centre. Single-molecule force spectroscopy (SMFS) experiments reveal 2-state and 3-state unfolding pathways for apo-azurin. The intermediate in the 3-state pathway occurs at an unfolding contour length of 7.5 nm from the native state. Steered molecular dynamics (SMD) simulations show that apo-azurin unfolds via a first transition state (TS) where β2Β–β8 and β7–β8 strand pairs rupture to form the intermediate, which subsequently unfolds by the collective rupture of remaining strands. SMFS experiments on holo-azurin exhibit an additional 4-state pathway besides the 2-state and 3-state pathways. The unfolding contour length leading to the first intermediate is 6.7 nm suggesting a sequestration of ~1 nm polypeptide chain length by the copper. SMD simulations reveal atomistic details of the copper sequestration and predict a combined β4–β7 pair and copper coordination sphere rupture to create the third TS in the 4-state pathway. Our systematic studies provide detailed mechanistic insights on modulation of protein mechanical properties by metal-cofactors.

Binding of interferon reduces the force of unfolding for interferon receptor 1

Differential signaling of the type I interferon receptor (IFNAR) has been correlated with the ability of its subunit, IFNAR1, to differentially recognize a large spectrum of different ligands, which involves intricate conformational re-arrangements of multiple interacting domains. To shed light onto the structural determinants governing ligand recognition, we compared the force-induced unfolding of the IFNAR1 ectodomain when bound to interferon and when free, using the atomic force microscope and steered molecular dynamics simulations. Unexpectedly, we find that IFNAR1 is easier to mechanically unfold when bound to interferon than when free. Analysis of the structures indicated that the origin of the reduction in unfolding forces is a conformational change in IFNAR1 induced by ligand binding.

Force Spectroscopy of the Plasmodium falciparum Vaccine Candidate Circumsporozoite Protein Suggests a Mechanically Pliable Repeat Region

The most effective vaccine candidate of malaria is based on the Plasmodium falciparum circumsporozoite protein (CSP), a major surface protein implicated in the structural strength, motility, and immune evasion properties of the infective sporozoites. It is suspected that reversible conformational changes of CSP are required for infection of the mammalian host, but the detailed structure and dynamic properties of CSP remain incompletely understood, limiting our understanding of its function in the infection. Here, we report the structural and mechanical properties of the CSP studied using single-molecule force spectroscopy on several constructs, one including the central region of CSP, which is rich in NANP amino acid repeats (CSP_rep), and a second consisting of a near full-length sequence without the signal and anchor hydrophobic domains (CSP_ΔHP). Our results show that the CSP_rep is heterogeneous, with 40% of molecules requiring virtually no mechanical force to unfold (<10 piconewtons (pN)), suggesting that these molecules are mechanically compliant and perhaps act as entropic springs, whereas the remaining 60% are partially structured with low mechanical resistance (∼70 pN). CSP_ΔHP having multiple force peaks suggests specifically folded domains, with two major populations possibly indicating the open and collapsed forms. Our findings suggest that the overall low mechanical resistance of the repeat region, exposed on the outer surface of the sporozoites, combined with the flexible full-length conformations of CSP, may provide the sporozoites not only with immune evasion properties, but also with lubricating capacity required during its navigation through the mosquito and vertebrate host tissues. We anticipate that these findings would further assist in the design and development of future malarial vaccines.

Homodimeric Escherichia coli Toxin CcdB (Controller of Cell Division or Death B Protein) Folds via Parallel Pathways

The existence of parallel pathways in the folding of proteins seems intuitive, yet remains controversial. We explore the folding kinetics of the homodimeric Escherichia coli toxin CcdB (Controller of Cell Division or Death B protein) using multiple optical probes and approaches. Kinetic studies performed as a function of protein and denaturant concentrations demonstrate that the folding of CcdB is a four-state process. The two intermediates populated during folding are present on parallel pathways. Both form by rapid association of the monomers in a diffusion limited manner and appear to be largely unstructured, as they are silent to the optical probes employed in the current study. The existence of parallel pathways is supported by the insensitivity of the amplitudes of the refolding kinetic phases to the different probes used in the study. More importantly, interrupted refolding studies and ligand binding studies clearly demonstrate that the native state forms in a biexponential manner, implying the presence of at least two pathways. Our studies indicate that the CcdA antitoxin binds only to the folded CcdB dimer and not to any earlier folding intermediates. Thus, despite being part of the same operon, the antitoxin does not appear to modulate the folding pathway of the toxin encoded by the downstream cistron. This study highlights the utility of ligand binding in distinguishing between sequential and parallel pathways in protein folding studies, while also providing insights into molecular interactions during folding in Type II toxin-antitoxin systems.

Significance of 1B and 2B domains in modulating elastic properties of lamin A

Nuclear lamins are type V intermediate filament proteins which form an elastic meshwork underlying the inner nuclear membrane. Lamins directly contribute to maintain the nuclear shape and elasticity. More than 400 mutations have been reported in lamin A that are involved in diseases known as laminopathies. These mutations are scattered mainly in the lamin rod domain along with some in its C-terminal domain. The contribution of the rod domain towards the elasticity of lamin A molecule was hitherto unknown. Here, we have elucidated the significance of the 1B and 2B domains of the rod in modulating the elastic behavior of lamin A by single-molecule force spectroscopy. In addition, we have also studied the network forming capacity of these domains and their corresponding viscoelastic behavior. We have shown that the 1B domain has the ability to form a lamin-like network and resists larger deformation. However at the single-molecular level, both the domains have comparable mechanical properties. The self-assembly of the 1B domain contributes to the elasticity of the lamin A network.

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.

The mechanochemistry of copper reports on the directionality of unfolding in model cupredoxin proteins

Understanding the directionality and sequence of protein unfolding is crucial to elucidate the underlying folding free energy landscape. An extra layer of complexity is added in metalloproteins, where a metal cofactor participates in the correct, functional fold of the protein. However, the precise mechanisms by which organometallic interactions are dynamically broken and reformed on (un)folding are largely unknown. Here we use single molecule force spectroscopy AFM combined with protein engineering and MD simulations to study the individual unfolding pathways of the blue-copper proteins azurin and plastocyanin. Using the nanomechanical properties of the native copper centre as a structurally embedded molecular reporter, we demonstrate that both proteins unfold via two independent, competing pathways. Our results provide experimental evidence of a novel kinetic partitioning scenario whereby the protein can stochastically unfold through two distinct main transition states placed at the N and C termini that dictate the direction in which unfolding occurs.

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.

Application of merged spectroscopic data combined with chemometric analysis for resolution of hemoglobin intermediates during chemical unfolding

Using tetradecyltrimethylammonium bromide (TTAB) as a surfactant denaturant, and augmentation of different spectroscopic data, helped to detect the intermediates of hemoglobin (Hb) during unfolding process. UV-vis, fluorescence, and circular dichroism spectroscopy were used simultaneously to monitor different aspects of hemoglobin species from the tertiary or secondary structure points of view. Application of the multivariate curve resolution-alternating least square (MCR-ALS), using the initial estimates of spectral profiles and appropriate constraints on different parts of augmented spectroscopic data, showed good efficiency for characterization of intermediates during Hb unfolding. These results indicated the existence of five protein species, including three intermediate-like compounds in this process. The unfolding pathway in the presence of TTAB included conversion of oxyhemoglobin into deoxyhemoglobin, and then ferrylhemoglobin, ferrihemoglobin or aquamethemoglobin, which finally transformed into hemichrome. This is the first application of chemometric analysis on the merged spectroscopic data related to chemical denaturation of a protein. These types of analysis in multisubunit proteins not only increase the domain of information, but also can reduce the ambiguities of the obtained results.

Optimizing the calculation of energy landscape parameters from single-molecule protein unfolding experiments

Single-molecule force spectroscopy using an atomic force microscope (AFM) can be used to measure the average unfolding force of proteins in a constant velocity experiment. In combination with Monte Carlo simulations and through the application of the Zhurkov-Bell model, information about the parameters describing the underlying unfolding energy landscape of the protein can be obtained. Using this approach, we have completed protein unfolding experiments on the polyprotein (I27)(5) over a range of pulling velocities. In agreement with previous work, we find that the observed number of protein unfolding events observed in each approach-retract cycle varies between one and five, due to the nature of the interactions between the polyprotein, the AFM tip, and the substrate, and there is an unequal unfolding probability distribution. We have developed a Monte Carlo simulation that incorporates the impact of this unequal unfolding probability distribution on the median unfolding force and the calculation of the protein unfolding energy landscape parameters. These results show that while there is a significant, unequal unfolding probability distribution, the unfolding energy landscape parameters obtained from use of the Zhurkov-Bell model are not greatly affected. This result is important because it demonstrates that the minimum acceptance criteria typically used in force extension experiments are justified and do not skew the calculation of the unfolding energy landscape parameters. We further validate this approach by determining the error in the energy landscape parameters for two extreme cases, and we provide suggestions for methods that can be employed to increase the level of accuracy in single-molecule experiments using polyproteins.

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.

Preventing disulfide bond formation weakens non-covalent forces among lysozyme aggregates

Nonnative disulfide bonds have been observed among protein aggregates in several diseases like amyotrophic lateral sclerosis, cataract and so on. The molecular mechanism by which formation of such bonds promotes protein aggregation is poorly understood. Here in this work we employ previously well characterized aggregation of hen eggwhite lysozyme (HEWL) at alkaline pH to dissect the molecular role of nonnative disulfide bonds on growth of HEWL aggregates. We employed time-resolved fluorescence anisotropy, atomic force microscopy and single-molecule force spectroscopy to quantify the size, morphology and non-covalent interaction forces among the aggregates, respectively. These measurements were performed under conditions when disulfide bond formation was allowed (control) and alternatively when it was prevented by alkylation of free thiols using iodoacetamide. Blocking disulfide bond formation affected growth but not growth kinetics of aggregates which were ∼50% reduced in volume, flatter in vertical dimension and non-fibrillar in comparison to control. Interestingly, single-molecule force spectroscopy data revealed that preventing disulfide bond formation weakened the non-covalent interaction forces among monomers in the aggregate by at least ten fold, thereby stalling their growth and yielding smaller aggregates in comparison to control. We conclude that while constrained protein chain dynamics in correctly disulfide bonded amyloidogenic proteins may protect them from venturing into partial folded conformations that can trigger entry into aggregation pathways, aberrant disulfide bonds in non-amyloidogenic proteins (like HEWL) on the other hand, may strengthen non-covalent intermolecular forces among monomers and promote their aggregation.

Experimental and theoretical studies of mechanical unfolding of different proteins

Mechanical properties of proteins are important for a wide range of biological processes including cell adhesion, muscle contraction, and protein translocation across biological membranes. It is necessary to reveal how proteins achieve their required mechanical stability under natural conditions in order to understand the biological processes and also to use the knowledge for constructing new biomaterials for medical and industrial purposes. In this connection, it is important to know how a protein will behave in response to various impacts. Theoretical and experimental works on mechanical unfolding of globular proteins will be considered in detail in this review.

Single-molecule studies on PolySUMO proteins reveal their mechanical flexibility

Proteins with β-sandwich and β-grasp topologies are resistant to mechanical unfolding as shown by single-molecule force spectroscopy studies. Their high mechanical stability has generally been associated with the mechanical clamp geometry present at the termini. However, there is also evidence for the importance of interactions other than the mechanical clamp in providing mechanical stability, which needs to be tested thoroughly. Here, we report the mechanical unfolding properties of ubiquitin-like proteins (SUMO1 and SUMO2) and their comparison with those of ubiquitin. Although ubiquitin and SUMOs have similar size and structural topology, they differ in their sequences and structural contacts, making them ideal candidates to understand the variations in the mechanical stability of a given protein topology. We observe a two-state unfolding pathway for SUMO1 and SUMO2, similar to that of ubiquitin. Nevertheless, the unfolding forces of SUMO1 (∼130 pN) and SUMO2 (∼120 pN) are lower than that of ubiquitin (∼190 pN) at a pulling speed of 400 nm/s, indicating their lower mechanical stability. The mechanical stabilities of SUMO proteins and ubiquitin are well correlated with the number of interresidue contacts present in their structures. From pulling speed-dependent mechanical unfolding experiments and Monte Carlo simulations, we find that the unfolding potential widths of SUMO1 (∼0.51 nm) and SUMO2 (∼0.33 nm) are much larger than that of ubiquitin (∼0.19 nm), indicating that SUMO1 is six times and SUMO2 is three times mechanically more flexible than ubiquitin. These findings might also be important in understanding the functional differences between ubiquitin and SUMOs.

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.

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.

Mapping the conformational stability of maltose binding protein at the residue scale using nuclear magnetic resonance hydrogen exchange experiments

Being able to differentiate local fluctuations from global folding-unfolding dynamics of a protein is of major interest for improving our understanding of structure-function determinants. The maltose binding protein (MBP), a protein that belongs to the maltose transport system, has a structure composed of two globular domains separated by a rigid-body "hinge bending". Here we determined, by using hydrogen exchange (HX) nuclear magnetic resonance experiments, the apparent stabilization free energies of 101 residues of MBP bound to β-cyclodextrin (MBP-βCD) under native conditions. We observed that the last helix of MBP (helix α14) has a lower protection factor than the rest of the protein. Further, HX experiments were performed using guanidine hydrochloride under subdenaturing conditions to discriminate between local fluctuations and global unfolding events and to determine the MBP-βCD energy landscape. The results show that helix α4 and a part of helices α5 and α6 are clearly grouped into a subdenaturing folding unit and represent a partially folded intermediate under native conditions. In addition, we observed that amide protons located in the hinge between the two globular domains share similar ΔG(gu)(app) and m values and should unfold simultaneously. These observations provide new points of view for improving our understanding of the thermodynamic stability and the mechanisms that drive folding-unfolding dynamics of proteins.

Iterative cloning, overexpression, purification and isotopic labeling of an engineered dimer of a Ca(2+)-binding protein of the βγ-crystallin superfamily from Methanosarcina acetivorans

βγ-Crystallins are a large superfamily of proteins found in vertebrate eye lens. They are hetero-dimers (linked in tandem by a specific peptide) and are shown to bind calcium. The monomers possess two β-strand rich greek-key motifs. Recently, a structurally closest member to the family of lens βγ-crystallins has been described, for the first time, from the archaea Methanosarcina acetivorans, which is named as M-crystallin. Unlike lens βγ-crystallins, M-crystallin exits as a monomer. Here, we synthesized a dimeric gene of M-crystallin in which two monomers are linked by a 10-amino acid residue coding sequence. The linker sequence in the target protein is long and flexible enough to reduce the proximity between the individual crystallins in the dimer. This methodology would be highly beneficial in designing polyproteins (two or more proteins linked in tandem to aid mechanical stretching studies) that are regularly used in single-molecule force spectroscopy. The dimer of M-crystallin was overexpressed in Escherichia coli BLR(DE3) strain. The overexpressed protein containing an N-terminal hexa-histidine tag was purified using nickel affinity chromatography and then by size-exclusion chromatography. Further, a method to purify isotopically ((15)N) labeled protein with high yield for NMR studies is reported. The uniformly (15)N-labeled M-crystallin dimer thus produced has been characterized by recording sensitivity enhanced 2D [(15)N-(1)H] HSQC and other optical spectroscopy techniques. Observation of only one set of peaks in the HSQC, along with the structural characterization using optical spectroscopy, suggests that the domains in the dimer possess similar structure as that of the monomer.