High-throughput virtual search of small molecules for controlling the mechanical stability of human CD4.

Protein mechanical stability determines the function of a myriad of proteins, especially proteins from the extracellular matrix. Failure to maintain protein mechanical stability may result in diseases and disorders such as cancer, cardiomyopathies, or muscular dystrophy. Thus, developing mutation-free approaches to enhance and control the mechanical stability of proteins using pharmacology-based methods may have important implications in drug development and discovery. Here, we present the first approach that employs computational high-throughput virtual screening and molecular docking to search for small molecules in chemical libraries that function as mechano-regulators of the stability of human cluster of differentiation 4, receptor of HIV-1. Using single-molecule force spectroscopy, we prove that these small molecules can increase the mechanical stability of CD4D1D2 domains over 4-fold in addition to modifying the mechanical unfolding pathways. Our experiments demonstrate that chemical libraries are a source of mechanoactive molecules and that drug discovery approaches provide the foundation of a new type of molecular function, that is, mechano-regulation, paving the way toward mechanopharmacology.

Unraveling the molecular determinants of the anti-phagocytic protein cloak of plague bacteria.

The pathogenic bacterium Yersina pestis is protected from macrophage engulfment by a capsule like antigen, F1, formed of long polymers of the monomer protein, Caf1. However, despite the importance of this pathogen, the mechanism of protection was not understood. Here we demonstrate how F1 protects the bacteria from phagocytosis. First, we show that Escherichia coli expressing F1 showed greatly reduced adherence to macrophages. Furthermore, the few cells that did adhere remained on the macrophage surface and were not engulfed. We then inserted, by mutation, an "RGDS" integrin binding motif into Caf1. This did not change the number of cells adhering to macrophages but increased the fraction of adherent cells that were engulfed. Therefore, F1 protects in two separate ways, reducing cell adhesion, possibly by acting as a polymer brush, and hiding innate receptor binding sites needed for engulfment. F1 is very robust and we show that E. coli expressing weakened mutant polymers are engulfed like the RGDS mutant. This suggests that innate attachment sites on the native cell surface are exposed if F1 is weakened. Single-molecule force spectroscopy (SMFS) experiments revealed that wild-type F1 displays a very high mechanical stability of 400 pN. However, the mechanical resistance of the destabilised mutants, that were fully engulfed, was only 20% weaker. By only marginally exceeding the mechanical force applied to the Caf1 polymer during phagocytosis it may be that the exceptional tensile strength evolved to resist the forces applied at this stage of engulfment.

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.

The influence of disulfide bonds on the mechanical stability of proteins is context dependent.

Disulfide bonds play a crucial role in proteins, modulating their stability and constraining their conformational dynamics. A particularly important case is that of proteins that need to withstand forces arising from their normal biological function and that are often disulfide bonded. However, the influence of disulfides on the overall mechanical stability of proteins is poorly understood. Here, we used single-molecule force spectroscopy (smFS) to study the role of disulfide bonds in different mechanical proteins in terms of their unfolding forces. For this purpose, we chose the pilus protein FimG from Gram-negative bacteria and a disulfide-bonded variant of the I91 human cardiac titin polyprotein. Our results show that disulfide bonds can alter the mechanical stability of proteins in different ways depending on the properties of the system. Specifically, disulfide-bonded FimG undergoes a 30% increase in its mechanical stability compared with its reduced counterpart, whereas the unfolding force of I91 domains experiences a decrease of 15% relative to the WT form. Using a coarse-grained simulation model, we rationalized that the increase in mechanical stability of FimG is due to a shift in the mechanical unfolding pathway. The simple topology-based explanation suggests a neutral effect in the case of titin. In summary, our results indicate that disulfide bonds in proteins act in a context-dependent manner rather than simply as mechanical lockers, underscoring the importance of considering disulfide bonds both computationally and experimentally when studying the mechanical properties of proteins.

Monitoring Oxidative Folding of a Single Protein Catalyzed by the Disulfide Oxidoreductase DsbA.

Oxidative folding, the process by which proteins fold and acquire disulfide bonds concurrently, is of critical importance for a wide range of biological processes. Generally, this process is catalyzed by oxidoreductase enzymes that facilitate oxidation and also bear chaperone functionality. Although this process has been well described qualitatively, fine yet important details remain obscured by a limited quantitative perspective, arising from the limitations in the application of bulk biochemical methods to the study of oxidative folding. In this work, we have applied single molecule force spectroscopy techniques to monitor in real time the process of oxidative folding as catalyzed by DsbA, the enzyme solely responsible for the catalysis of oxidative folding in the bacterial periplasm. We provide a quantitative and detailed description of the catalytic mechanism utilized by DsbA that offers insight into the entire sequence of events that occurs in the periplasm from the unfolded-reduced state to the folded-oxidized protein. We have compared our results with those of protein disulfide-isomerase, the eukaryotic counterpart of DsbA, allowing us to devise a general mechanism for oxidative folding that also reflects upon the physiological functions and demands of these enzymes in vivo.

Protein folding under mechanical forces: a physiological view.

Mechanical forces regulate the function of numerous proteins relevant to physiology. The functions and folding of proteins have been under scrutiny for decades, but it was not until recently that mechanical forces have been considered. Here, we review different techniques for studying protein folding, highlighting their physiological significance.

The strongest protein binder is surprisingly labile.

Bacterial adhesins are cell-surface proteins that anchor to the cell wall of the host. The first stage of infection involves the specific attachment to fibrinogen (Fg), a protein found in human blood. This attachment allows bacteria to colonize tissues causing diseases such as endocarditis. The study of this family of proteins is hence essential to develop new strategies to fight bacterial infections. In the case of the Gram-positive bacterium Staphylococcus aureus, there exists a class of adhesins known as microbial surface components recognizing adhesive matrix molecules (MSCRAMMs). Here, we focus on one of them, the clumping factor A (ClfA), which has been found to bind Fg through the dock-lock-latch mechanism. Interestingly, it has recently been discovered that MSCRAMM proteins employ a catch-bond to withstand forces exceeding 2 nN, making this type of interaction as mechanically strong as a covalent bond. However, it is not known whether this strength is an evolved feature characteristic of the bacterial protein or is typical only of the interaction with its partner. Here, we combine single-molecule force spectroscopy, biophysical binding assays, and molecular simulations to study the intrinsic mechanical strength of ClfA. We find that despite the extremely high forces required to break its interactions with Fg, ClfA is not by itself particularly strong. Integrating the results from both theory and experiments we dissect contributions to the mechanical stability of this protein.

Probing the chemistry of thioredoxin catalysis with force.

Thioredoxins are enzymes that catalyse disulphide bond reduction in all living organisms. Although catalysis is thought to proceed through a substitution nucleophilic bimolecular (S(N)2) reaction, the role of the enzyme in modulating this chemical reaction is unknown. Here, using single-molecule force-clamp spectroscopy, we investigate the catalytic mechanism of Escherichia coli thioredoxin (Trx). We applied mechanical force in the range of 25-600 pN to a disulphide bond substrate and monitored the reduction of these bonds by individual enzymes. We detected two alternative forms of the catalytic reaction, the first requiring a reorientation of the substrate disulphide bond, causing a shortening of the substrate polypeptide by 0.79 +/- 0.09 A (+/- s.e.m.), and the second elongating the substrate disulphide bond by 0.17 +/- 0.02 A (+/- s.e.m.). These results support the view that the Trx active site regulates the geometry of the participating sulphur atoms with sub-ångström precision to achieve efficient catalysis. Our results indicate that substrate conformational changes may be important in the regulation of Trx activity under conditions of oxidative stress and mechanical injury, such as those experienced in cardiovascular disease. Furthermore, single-molecule atomic force microscopy techniques, as shown here, can probe dynamic rearrangements within an enzyme's active site during catalysis that cannot be resolved with any other current structural biological technique.

Evolution of CRISPR-associated endonucleases as inferred from resurrected proteins.

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated Cas9 is an effector protein that targets invading DNA and plays a major role in the prokaryotic adaptive immune system. Although Streptococcus pyogenes CRISPR-Cas9 has been widely studied and repurposed for applications including genome editing, its origin and evolution are poorly understood. Here, we investigate the evolution of Cas9 from resurrected ancient nucleases (anCas) in extinct firmicutes species that last lived 2.6 billion years before the present. We demonstrate that these ancient forms were much more flexible in their guide RNA and protospacer-adjacent motif requirements compared with modern-day Cas9 enzymes. Furthermore, anCas portrays a gradual palaeoenzymatic adaptation from nickase to double-strand break activity, exhibits high levels of activity with both single-stranded DNA and single-stranded RNA targets and is capable of editing activity in human cells. Prediction and characterization of anCas with a resurrected protein approach uncovers an evolutionary trajectory leading to functionally flexible ancient enzymes.

A designed protein as experimental model of primordial folding.

How do proteins accomplish folding during early evolution? Theoretically the mechanism involves the selective stabilization of the native structure against all other competing compact conformations in a process that involves cumulative changes in the amino acid sequence along geological timescales. Thus, an evolved protein folds into a single structure at physiological temperature, but the conformational competition remains latent. For natural proteins such competition should emerge only near cryogenic temperatures, which places it beyond experimental testing. Here, we introduce a designed monomeric miniprotein (FSD-1ss) that within biological temperatures (330-280 K) switches between simple fast folding and highly complex conformational dynamics in a structurally degenerate compact ensemble. Our findings demonstrate the physical basis for protein folding evolution in a designed protein, which exhibits poorly evolved or primordial folding. Furthermore, these results open the door to the experimental exploration of primitive folding and the switching between alternative protein structures that takes place in evolutionary branching points and prion diseases, as well as the benchmarking of de novo design methods.

Mechanochemical Evolution of Disulfide Bonds in Proteins.

Disulfide bonds play a pivotal role in the mechanical stability of proteins. Numerous proteins that are known to be exposed to mechanical forces in vivo contain disulfide bonds. The presence of cryptic disulfide bonds in a protein structure may be related to its resistance to an applied mechanical force. Disulfide bonds in proteins tend to be highly conserved but their evolution might be directly related to the evolution of the protein mechanical stability. Hence, tracking the evolution of disulfide bonds in a protein can help to derive crucial stability/function correlations in proteins that are exposed to mechanical forces. Phylogenic analysis and ancestral sequence reconstruction (ASR) allow tracking the evolution of proteins from the past ancestors to our modern days and also establish correlations between proteins from different species. In addition, ASR can be combined with single-molecule force spectroscopy (smFS) to investigate the mechanical properties of proteins including the occurrence and function of disulfide bonds. Here we present a detailed protocol to study the mechanochemical evolution of proteins using a fragment of the giant muscle protein titin as example. The protocol can be easily adapted to AFS studies of any resurrected mechanical force bearing protein of interest.

Basal oxidation of conserved cysteines modulates cardiac titin stiffness and dynamics.

Titin, as the main protein responsible for the passive stiffness of the sarcomere, plays a key role in diastolic function and is a determinant factor in the etiology of heart disease. Titin stiffness depends on unfolding and folding transitions of immunoglobulin-like (Ig) domains of the I-band, and recent studies have shown that oxidative modifications of cryptic cysteines belonging to these Ig domains modulate their mechanical properties in vitro. However, the relevance of this mode of titin mechanical modulation in vivo remains largely unknown. Here, we describe the high evolutionary conservation of titin mechanical cysteines and show that they are remarkably oxidized in murine cardiac tissue. Mass spectrometry analyses indicate a similar landscape of basal oxidation in murine and human myocardium. Monte Carlo simulations illustrate how disulfides and S-thiolations on these cysteines increase the dynamics of the protein at physiological forces, while enabling load- and isoform-dependent regulation of titin stiffness. Our results demonstrate the role of conserved cysteines in the modulation of titin mechanical properties in vivo and point to potential redox-based pathomechanisms in heart disease.

Single-molecule force spectroscopy approach to enzyme catalysis.

Enzyme catalysis has been traditionally studied using a diverse set of techniques such as bulk biochemistry, x-ray crystallography, and NMR. Recently, single-molecule force spectroscopy by atomic force microscopy has been used as a new tool to study the catalytic properties of an enzyme. In this approach, a mechanical force ranging up to hundreds of piconewtons is applied to the substrate of an enzymatic reaction, altering the conformational energy of the substrate-enzyme interactions during catalysis. From these measurements, the force dependence of an enzymatic reaction can be determined. The force dependence provides valuable new information about the dynamics of enzyme catalysis with sub-angstrom resolution, a feat unmatched by any other current technique. To date, single-molecule force spectroscopy has been applied to gain insight into the reduction of disulfide bonds by different enzymes of the thioredoxin family. This minireview aims to present a perspective on this new approach to study enzyme catalysis and to summarize the results that have already been obtained from it. Finally, the specific requirements that must be fulfilled to apply this new methodology to any other enzyme will be discussed.

Estimation of protein folding free energy barriers from calorimetric data by multi-model Bayesian analysis.

The realization that folding free energy barriers can be small enough to result in significant population of the species at the barrier top has sprouted in several methods to estimate folding barriers from equilibrium experiments. Some of these approaches are based on fitting the experimental thermogram measured by differential scanning calorimetry (DSC) to a one-dimensional representation of the folding free-energy surface (FES). Different physical models have been used to represent the FES: (1) a Landau quartic polynomial as a function of the total enthalpy, which acts as an order parameter; (2) the projection onto a structural order parameter (i.e. number of native residues or native contacts) of the free energy of all the conformations generated by Ising-like statistical mechanical models; and (3) mean-field models that define conformational entropy and stabilization energy as functions of a continuous local order parameter. The fundamental question that emerges is how can we obtain robust, model-independent estimates of the thermodynamic folding barrier from the analysis of DSC experiments. Here we address this issue by comparing the performance of various FES models in interpreting the thermogram of a protein with a marginal folding barrier. We chose the small α-helical protein PDD, which folds-unfolds in microseconds crossing a free energy barrier previously estimated as ~1 RT. The fits of the PDD thermogram to the various models and assumptions produce FES with a consistently small free energy barrier separating the folded and unfolded ensembles. However, the fits vary in quality as well as in the estimated barrier. Applying Bayesian probabilistic analysis we rank the fit performance using a statistically rigorous criterion that leads to a global estimate of the folding barrier and its precision, which for PDD is 1.3 ± 0.4 kJ mol(-1). This result confirms that PDD folds over a minor barrier consistent with the downhill folding regime. We have further validated the multi-model Bayesian approach through the analysis of two additional protein systems: gpW, a midsize single-domain with α + ß topology that also folds in microseconds and has been previously catalogued as a downhill folder, and α-spectrin SH3, a domain of similar size but with a ß-barrel fold, slow-folding kinetics and two-state-like thermodynamics. From a general viewpoint, the Bayesian analysis developed here results in a statistically robust, virtually model-independent, method to estimate the thermodynamic free-energy barriers to protein folding from DSC thermograms. Our method appears to be sufficiently accurate to consistently detect small differences in the barrier height, and thus opens up the possibility of characterizing experimentally the changes in thermodynamic folding barriers induced by single-point mutations on proteins within the downhill regime.

Role of conservative mutations in protein multi-property adaptation.

Protein physicochemical properties must undergo complex changes during evolution, as a response to modifications in the organism environment, the result of the proteins taking up new roles or because of the need to cope with the evolution of molecular interacting partners. Recent work has emphasized the role of stability and stability-function trade-offs in these protein adaptation processes. In the present study, on the other hand, we report that combinations of a few conservative, high-frequency-of-fixation mutations in the thioredoxin molecule lead to largely independent changes in both stability and the diversity of catalytic mechanisms, as revealed by single-molecule atomic force spectroscopy. Furthermore, the changes found are evolutionarily significant, as they combine typically hyperthermophilic stability enhancements with modulations in function that span the ranges defined by the quite different catalytic patterns of thioredoxins from bacterial and eukaryotic origin. These results suggest that evolutionary protein adaptation may use, in some cases at least, the potential of conservative mutations to originate a multiplicity of evolutionarily allowed mutational paths leading to a variety of protein modulation patterns. In addition the results support the feasibility of using evolutionary information to achieve protein multi-feature optimization, an important biotechnological goal.

Large-scale modulation of thermodynamic protein folding barriers linked to electrostatics.

Protein folding barriers, which range from zero to the tens of RT that result in classical two-state kinetics, are primarily determined by protein size and structural topology [Plaxco KW, Simons KT, Baker D (1998) J Mol Biol 277:985-994]. Here, we investigate the thermodynamic folding barriers of two relatively large proteins of the same size and topology: bovine alpha-lactalbumin (BLA) and hen-egg-white lysozyme (HEWL). From the analysis of differential scanning calorimetry experiments with the variable-barrier model [Muñoz V, Sanchez-Ruiz JM (2004) Proc Natl Acad Sci USA 101:17646-17651] we obtain a high barrier for HEWL and a marginal folding barrier for BLA. These results demonstrate a remarkable tuning range of at least 30 kJ/mol (i.e., five to six orders of magnitude in population) within a unique protein scaffold. Experimental and theoretical analyses on these proteins indicate that the surprisingly small thermodynamic folding barrier of BLA arises from the stabilization of partially unfolded conformations by electrostatic interactions. Interestingly, there is clear reciprocity between the barrier height and the biological function of the two proteins, suggesting that the marginal barrier of BLA is a product of natural selection. Electrostatic surface interactions thus emerge as a mechanism for the modulation of folding barriers in response to special functional requirements within a given structural fold.

Enzymatically produced cellulose nanocrystals as reinforcement for waterborne polyurethane and its applications.

Waterborne polyurethanes (WBPUs) have been proposed as ecofriendly elastomers with several applications in coatings and adhesives. WBPU's physicochemical properties can be enhanced by the addition of cellulose nanocrystals (CNCs). The way CNCs are isolated has a strong effect on their properties and can determine their role as reinforcement. In this work, CNCs produced using ancestral endoglucanase (EnCNCs) were used as reinforcement for WBPU and compared with CNC produced by sulfuric acid hydrolysis (AcCNC). The enzymatic method produced highly thermostable and crystalline CNCs. The addition of small contents of EnCNCs improved the thermomechanical stability and mechanical properties of WBPUs, even better than commercial AcCNCs. Besides, WBPU reinforced by adding EnCNCs was studied as a coating for paper materials, increasing its abrasion resistance and as electrospun nanocomposite mats where EnCNCs helped maintaining the morphology of the fibers.

Navigating the downhill protein folding regime via structural homologues.

Proteins that fold over free-energy barriers

Entropic bonding of the type 1 pilus from experiment and simulation.

The type 1 pilus is a bacterial filament consisting of a long coiled proteic chain of subunits joined together by non-covalent bonding between complementing ß -strands. Its strength and structural stability are critical for its anchoring function in uropathogenic Escherichia coli bacteria. The pulling and unravelling of the FimG subunit of the pilus was recently studied by atomic force microscopy experiments and steered molecular dynamics simulations (Alonso-Caballero et al. 2018 Nat. Commun. 9, 2758. (doi:10.1038/s41467-018-05107-6)). In this work, we perform a quantitative comparison between experiment and simulation, showing a good agreement in the underlying work values for the unfolding. The simulation results are then used to estimate the free energy difference for the detachment of FimG from the complementing strand of the neighbouring subunit in the chain, FimF. Finally, we show that the large free energy difference for the unravelling and detachment of the subunits which leads to the high stability of the chain is entirely entropic in nature.

Instrumental Effects in the Dynamics of an Ultrafast Folding Protein under Mechanical Force.

The analysis and interpretation of single molecule force spectroscopy (smFS) experiments is often complicated by hidden effects from the measuring device. Here we investigate these effects in our recent smFS experiments on the ultrafast folding protein gpW, which has been previously shown to fold without crossing a free energy barrier in the absence of force (i.e., downhill folding). Using atomic force microscopy (AFM) smFS experiments, we found that a very small force of ∼5 pN brings gpW near its unfolding midpoint and results in two-state (un)folding patterns that indicate the emergence of a force-induced free energy barrier. The change in the folding regime is concomitant with a 30,000-fold slowdown of the folding and unfolding times, from a few microseconds that it takes gpW to (un)fold at the midpoint temperature to seconds in the AFM. These results are puzzling because the barrier induced by force in the folding free energy landscape of gpW is far too small to account for such a difference in time scales. Here we use recently developed theoretical methods to resolve the origin of the strikingly slow dynamics of gpW under mechanical force. We find that, while the AFM experiments correctly capture the equilibrium distance distribution, the measured dynamics are entirely controlled by the response of the cantilever and polyprotein linker, which is much slower than the protein conformational dynamics. This interpretation is likely applicable to the folding of other small biomolecules in smFS experiments, and becomes particularly important in the case of systems with fast folding dynamics and small free energy barriers, and for instruments with slow response times.