The effects of macromolecular crowding on the mechanical stability of protein molecules

Macromolecular crowding, a common phenomenon in the cellular environments, can significantly affect the thermodynamic and kinetic properties of proteins. A single-molecule method based on atomic force microscopy (AFM) was used to investigate the effects of macromolecular crowding on the forces required to unfold individual protein molecules. It was found that the mechanical stability of ubiquitin molecules was enhanced by macromolecular crowding from added dextran molecules. The average unfolding force increased from 210 pN in the absence of dextran to 234 pN in the presence of 300 g/L dextran at a pulling speed of 0.25 microm/sec. A theoretical model, accounting for the effects of macromolecular crowding on the native and transition states of the protein molecule by applying the scaled-particle theory, was used to quantitatively explain the crowding-induced increase in the unfolding force. The experimental results and interpretation presented could have wide implications for the many proteins that experience mechanical stresses and perform mechanical functions in the crowded environment of the cell.     

Mechanical stability of multidomain proteins and novel mechanical clamps

We estimate the size of mechanostability for 318 multidomain proteins which are single-chain and contain up to 1021 amino acids. We predict existence of novel types of mechanical clamps in which interdomain contacts play an essential role. Mechanical clamps are structural regions which are the primary source of a protein's resistance to pulling. Among these clamps there is one that opposes tensile stress due to two domains swinging apart. This movement strains and then ruptures the contacts that hold the two domains together. Another clamp also involves tensile stress but it originates from an immobilization of a structural region by a surrounding knot-loop (without involving any disulfide bonds). Still another mechanism involves shear between helical regions belonging to two domains. We also consider the amyloid-prone cystatin C which provides an example of a two-chain 3D domain-swapped protein. We predict that this protein should withstand remarkably large stress, perhaps of order 800 pN, when inducing a shearing strain. The survey is generated through molecular dynamics simulations performed within a structure-based coarse grained model.     

Mechanical unfolding of covalently linked GroES: Evidence of structural subunit intermediates

It is difficult to determine the structural stability of the individual subunits or protomers of many proteins in the cell that exist in an oligomeric or complexed state. In this study, we used single‐molecule force spectroscopy on seven subunits of covalently linked cochaperonin GroES (ESC7) to evaluate the structural stability of the subunit. A modified form of ESC7 was immobilized on a mica surface. The force‐extension profile obtained from the mechanical unfolding of this ESC7 showed a distinctive sawtooth pattern that is typical for multimodular proteins. When analyzed according to the worm‐like chain model, the contour lengths calculated from the peaks in the profile suggested that linked‐GroES subunits unfold in distinct steps after the oligomeric ring structure of ESC7 is disrupted. The evidence that structured subunits of ESC7 withstand external force to some extent even after the perturbation of the oligomeric ring structure suggests that a stable monomeric intermediate is an important component of the equilibrium unfolding reaction of GroES.     

Mechanical measurement of single-molecule binding rates: kinetics of DNA helix-destabilization by T4 gene 32 protein

Bacteriophage T4 gene 32 protein (gp32) is a single-stranded DNA (ssDNA) binding protein, and is essential for DNA replication, recombination and repair. While gp32 binds preferentially and cooperatively to ssDNA, it has not been observed to lower the thermal melting temperature of natural double-stranded DNA (dsDNA). However, in single-molecule stretching experiments, gp32 significantly destabilizes lambda DNA. In this study, we develop a theory of the effect of the protein on single dsDNA stretching curves, and apply it to the measured dependence of the DNA overstretching force on pulling rate in the presence of the full-length and two truncated forms of the protein. This allows us to calculate the rate of cooperative growth of single clusters of protein along ssDNA that are formed as the dsDNA molecule is stretched, as well as determine the site size of the protein binding to ssDNA. The rate of cooperative binding (ka) of both gp32 and of its proteolytic fragment *I (which lacks 48 residues from the C terminus) varies non-linearly with protein concentration, and appears to exceed the diffusion limit. We develop a model of protein association with the ends of growing clusters of cooperatively bound protein enhanced by 1-D diffusion along dsDNA, under the condition of protein excess. Upon globally fitting ka versus protein concentration, we determine the binding site size and the non-cooperative binding constants to dsDNA for gp32 and I. Our experiment mimics the growth of clusters of gp32 that likely exist at the DNA replication fork in vivo, and explains the origin of the "kinetic block" to dsDNA melting by gene 32 protein observed in thermal melting experiments.     

Investigation of the interactions between aptamer and misfolded proteins: From monomer and oligomer to fibril by single‐molecule force spectroscopy

Increasing knowledge on the understanding interactions of aptamer with misfolded proteins (including monomer, oligomer, and amyloid fibril) is crucial for development of aggregation inhibitors and diagnosis of amyloid diseases. Herein, the interactions of lysozyme monomer–, oligomer‐, and amyloid fibril–aptamer were investigated using single‐molecule force spectroscopy. The results revealed that the aptamer screened against lysozyme monomer could also bind to oligomer and amyloid fibril, in spite of the recognition at a lower binding probability. It may be attributed to the inherent structural differences of misfolded proteins and the flexible conformation of aptamer. In addition, dynamic force spectra showed that there were similar dissociation paths in the dissociation process of lysozyme monomer–, oligomer‐, and amyloid fibril–aptamer complexes. It showed that the dissociation only passed 1 energy barrier from the binding state to the detachment. However, the dynamic parameters suggested that the oligomer‐ and amyloid fibril–aptamer were more stable than lysozyme monomer–aptamer. The phenomena may result from the exposure of aptamer‐recognized sequences on the surface and the electrostatic interactions. This work demonstrated that single‐molecule force spectroscopy could be a powerful tool to study the binding behavior of the aptamer with misfolded proteins at single‐molecule level, providing abundant information for researches and comprehensive applications of aptamer probes in diagnosis of amyloid diseases.     

Salt dependent binding of T4 gene 32 protein to single and double-stranded DNA: single molecule force spectroscopy measurements

Bacteriophage T4 gene 32 protein (gp32) is a well-studied representative of the large family of single-stranded DNA (ssDNA) binding proteins, which are essential for DNA replication, recombination and repair. Surprisingly, gp32 has not previously been observed to melt natural dsDNA. At the same time, *I, a truncated version of gp32 lacking its C-terminal domain (CTD), was shown to decrease the melting temperature of natural DNA by about 50 deg. C. This profound difference in the duplex destabilizing ability of gp32 and *I is especially puzzling given that the previously measured binding of both proteins to ssDNA was similar. Here, we resolve this apparent contradiction by studying the effect of gp32 and *I on the thermodynamics and kinetics of duplex DNA melting. We use a previously developed single molecule technique for measuring the non-cooperative association constants (K(ds)) to double-stranded DNA to determine K(ds) as a function of salt concentration for gp32 and *I. We then develop a new single molecule method for measuring K(ss), the association constant of these proteins to ssDNA. Comparing our measured binding constants to ssDNA for gp32 and *I we see that while they are very similar in high salt, they strongly diverge at [Na+] < 0.2 M. These results suggest that intact protein must undergo a conformational rearrangement involving the CTD that is in pre-equilibrium to its non-cooperative binding to both dsDNA and ssDNA. This lowers the effective concentration of protein available for binding, which in turn lowers the rate at which it can destabilize dsDNA. For the first time, we quantify the free energy of this CTD unfolding, and show it to be strongly salt dependent and associated with sodium counter-ion condensation on the CTD.     

Mechanical stability and differentially conserved physical–chemical properties of titin Ig‐domains

The mechanisms that determine mechanical stabilities of protein folds remain elusive. Our understanding of these mechanisms is vital to both bioengineering efforts and to the better understanding and eventual treatment of pathogenic mutations affecting mechanically important proteins such as titin. We present a new approach to analyze data from single‐molecule force spectroscopy for different domains of the giant muscle protein titin. The region of titin found in the I‐band of a sarcomere is composed of about 40 Ig‐domains and is exposed to force under normal physiological conditions and connects the free‐hanging ends of the myosin filaments to the Z‐disc. Recent single‐molecule force spectroscopy data show a mechanical hierarchy in the I‐band domains. Domains near the C‐terminus in this region unfold at forces two to three times greater than domains near the beginning of the I‐band. Though all of these Ig‐domains are thought to share a fold and topology common to members of the Ig‐like fold family, the sequences of neighboring domains vary greatly with an average sequence identity of only 25%. We examine in this study the relation of these unique mechanical stabilities of each I‐band Ig domain to specific, conserved physical–chemical properties of amino acid sequences in related Ig domains. We find that the sequences of each individual titin Ig domain are very highly conserved, with an average sequence identity of 79% across species that are divergent as humans, chickens, and zebra fish. This indicates that the mechanical properties of each domain are well conserved and tailored to its unique position in the titin molecule. We used the PCPMer software to determine the conservation of amino acid properties in titin Ig domains grouped by unfolding forces into “strong” and “weak” families. We found two motifs unique to each family that may have some role in determining the mechanical properties of these Ig domains. A detailed statistical analysis of properties of individual residues revealed several positions that displayed differentially conserved properties in strong and weak families. In contrast to previous studies, we find evidence that suggests that the mechanical stability of Ig domains is determined by several residues scattered across the β‐sandwich fold, and force sensitive residues are not only confined to the A′‐G region. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Determination of the elastic properties of short ssDNA molecules by mechanically folding and unfolding DNA hairpins

The characterization of elastic properties of biopolymers is crucial to understand many molecular reactions determined by conformational bending fluctuations of the polymer. Direct measurement of such elastic properties using single‐molecule methods is usually hindered by the intrinsic tendency of such biopolymers to form high‐order molecular structures. For example, single‐stranded deoxyribonucleic acids (ssDNA) tend to form secondary structures such as local double helices that prevent the direct measurement of the ideal elastic response of the ssDNA. In this work, we show how to extract the ideal elastic response in the entropic regime of short ssDNA molecules by mechanically pulling two‐state DNA hairpins of different contour lengths. This is achieved by measuring the force dependence of the molecular extension and stiffness on mechanically folding and unfolding the DNA hairpin. Both quantities are fit to the worm‐like chain elastic model giving values for the persistence length and the interphosphate distance. This method can be used to unravel the elastic properties of short ssDNA and RNA sequences and, more generally, any biopolymer that can exhibit a cooperative two‐state transition between mechanically folded and unfolded states (such as proteins). © 2014 Wiley Periodicals, Inc. Biopolymers 101: 1193–1199, 2014.     

Kinetic regulation of single DNA molecule denaturation by T4 gene 32 protein structural domains

Bacteriophage T4 gene 32 protein (gp32) specifically binds single-stranded DNA, a property essential for its role in DNA replication, recombination, and repair. Although on a thermodynamic basis, single-stranded DNA binding proteins should lower the thermal melting temperature of double-stranded DNA (dsDNA), gp32 does not. Using single molecule force spectroscopy, we show for the first time that gp32 is capable of slowly destabilizing natural dsDNA. Direct measurements of single DNA molecule denaturation and renaturation kinetics in the presence of gp32 and its proteolytic fragments reveal three types of kinetic behavior, attributable to specific protein structural domains, which regulate gp32's helix-destabilizing capabilities. Whereas the full-length protein exhibits very slow denaturation kinetics, a truncate lacking the acidic C-domain exhibits much faster kinetics. This may reflect a steric blockage of the DNA binding site and/or a conformational change associated with this domain. Additional removal of the N-domain, which is needed for binding cooperativity, further increases the DNA denaturation rate, suggesting that both of these domains are critical to the regulation of gp32's helix-destabilization capabilities. This regulation is potentially biologically significant because uncontrolled helix-destabilization would be lethal to the cell. We also obtain equilibrium measurements of the helix-coil transition free energy in the presence of these proteins for the first time.     

Unbinding and unfolding of adhesion protein complexes through stretching: Interplay between shear and tensile mechanical clamps

Using coarse‐grained molecular dynamics simulations, we analyze mechanically induced dissociation and unfolding of the protein complex CD48‐2B4. This heterodimer is an indispensable component of the immunological system: 2B4 is a receptor on natural killer cells whereas CD48 is expressed on surfaces of various immune cells. So far, its mechanostability has not been assessed either experimentally or theoretically. We find that the dissociation processes strongly depend on the direction of pulling and may take place in several pathways. Interestingly, the CD48‐2B4 interface can be divided into three distinct patches that act as units when resisting the pulling forces. At experimentally accessible pulling speeds, the characteristic mechanostability forces are in the range between 100 and 200 pN, depending on the pulling direction. These characteristic forces need not be associated with tensile forces involved in the act of separation of the complex because prior shear‐involving unraveling within individual proteins may give rise to a higher force peak. Proteins 2014; 82:3144–3153. © 2014 Wiley Periodicals, Inc.     

A divide-and-conquer approach to determine the Pareto frontier for optimization of protein engineering experiments

In developing improved protein variants by site-directed mutagenesis or recombination, there are often competing objectives that must be considered in designing an experiment (selecting mutations or breakpoints): stability versus novelty, affinity versus specificity, activity versus immunogenicity, and so forth. Pareto optimal experimental designs make the best trade-offs between competing objectives. Such designs are not "dominated"; that is, no other design is better than a Pareto optimal design for one objective without being worse for another objective. Our goal is to produce all the Pareto optimal designs (the Pareto frontier), to characterize the trade-offs and suggest designs most worth considering, but to avoid explicitly considering the large number of dominated designs. To do so, we develop a divide-and-conquer algorithm, Protein Engineering Pareto FRontier (PEPFR), that hierarchically subdivides the objective space, using appropriate dynamic programming or integer programming methods to optimize designs in different regions. This divide-and-conquer approach is efficient in that the number of divisions (and thus calls to the optimizer) is directly proportional to the number of Pareto optimal designs. We demonstrate PEPFR with three protein engineering case studies: site-directed recombination for stability and diversity via dynamic programming, site-directed mutagenesis of interacting proteins for affinity and specificity via integer programming, and site-directed mutagenesis of a therapeutic protein for activity and immunogenicity via integer programming. We show that PEPFR is able to effectively produce all the Pareto optimal designs, discovering many more designs than previous methods. The characterization of the Pareto frontier provides additional insights into the local stability of design choices as well as global trends leading to trade-offs between competing criteria.     

Structural impact of cations on lipid bilayer models: Nanomechanical properties by AFM-force spectroscopy

Atomic Force Microscopy (AFM) has become an invaluable tool for studying the micro- and nanoworlds. As a stand-alone, high-resolution imaging technique and force transducer, it defies most other surface instrumentation in ease of use, sensitivity and versatility. The main strength of AFM relies on the possibility to operate in an aqueous environment on a wide variety of biological samples, from single molecules – DNA or proteins – to macromolecular assemblies like biological membranes. Understanding the effect of mechanical stress on membranes is of primary importance in biophysics, since cells are known to perform their function under a complex combination of forces. In the later years, AFM-based Force-Spectroscopy (AFM-FS) has provided a new vista on membrane mechanics in a confined area within the nanometer realm, where most of the specific molecular interactions take place. Lipid membranes are electrostatically charged entities that physiologically coexist with electrolyte solutions. Thus, specific interactions with ions are a matter of considerable interest. The distribution of ions in the solution and their interaction with the membranes are factors that substantially modify the structure and dynamics of the cell membranes. Furthermore, signaling processes are modified by the membrane capability of retaining ions. Supported Lipid Bilayers (SLBs) are a versatile tool to investigate phospholipid membranes mimicking biological surfaces. In the present contribution, we review selected experiments on the mechanical stability of SLBs as models of lipid membranes by means of AFM-FS, with special focus on the effect of cations and ionic strength in the overall nanomechanical stability.     

Characterization of N-cadherin unbinding properties in non-malignant (HCV29) and malignant (T24) bladder cells

The expression of N‐cadherin, characteristic of various cancers, very often leads to changes in the cells' adhesive properties. Thus, we sought to find out if N‐cadherin expressed in various, but cancer‐related cells, differs in its functional properties that could contribute to variations in cells' phenotypes. In our work, measurements of an unbinding force of a single N‐cadherin molecule, probed with the same antibody both on a surface of living non‐malignant (HCV29) and malignant cells (T24) of bladder cancer, were carried out with the use of an atomic force microscopy. The results show the enhanced N‐cadherin level in T24 malignant cells (8.7% vs. 3.6% obtained for non‐malignant one), confirmed by the Western blot and the immunohistochemical staining. The effect was accompanied by changes in unbinding properties of an individual N‐cadherin molecule. Lower unbinding force values (26.1 ± 7.1 pN) in non‐malignant cells reveal less stable N‐cadherin complexes, as compared to malignant cells (61.7 ± 14.6 pN). This suggests the cancer‐related changes in a structure of the binding site of the antibody, located at the extracellular domain of N‐cadherin. Copyright © 2011 John Wiley & Sons, Ltd.     

Gastrodianin-like mannose-binding proteins: a novel class of plant proteins with antifungal properties

The orchid Gastrodia elata depends on the fungus Armillaria mellea to complete its life cycle. In the interaction, fungal hyphae penetrate older, nutritive corms but not newly formed corms. From these corms, a protein fraction with in vitro activity against plant-pathogenic fungi has previously been purified. Here, the sequence of gastrodianin, the main constituent of the antifungal fraction, is reported. Four isoforms that encoded two different mature proteins were identified at the cDNA level. Another isoform was detected in sequenced peptides. Because the antifungal activity of gastrodianins produced in and purified from Escherichia coli and Nicotiana tabacum was comparable to that of gastrodianin purified from the orchid, gastrodianins are the active component of the antifungal fractions. Gastrodianin accumulation is probably an important part of the mechanism by which the orchid controls Armillaria penetration. Gastrodianin was found to be homologous to monomeric mannose-binding proteins of other orchids, of which at least one (Epipactis helleborine mannose-binding protein) also displayed in vitro antifungal activity. This establishes the gastrodianin-like proteins (GLIPs) as a novel class of antifungal proteins.     

The presence of non‐native helical structure in the unfolding of a beta‐sheet protein MPT63

MPT63, a major secreted protein from Mycobacterium tuberculosis, has been shown to have immunogenic properties and has been implicated in virulence. MPT63 is a β‐sandwich protein containing 11 β strands and a very short stretch of 3₁₀ helix. The detailed experimental and computational study reported here investigates the equilibrium unfolding transition of MPT63. It is shown that in spite of being a complete β‐sheet protein, MPT63 has a strong propensity toward helix structures in its early intermediates. Far UV‐CD and FTIR spectra clearly suggest that the low‐pH intermediate of MTP63 has enhanced helical content, while fluorescence correlation spectroscopy suggests a significant contraction. Molecular dynamics simulation complements the experimental results indicating that the unfolded state of MPT63 traverses through intermediate forms with increased helical characteristics. It is found that this early intermediate contains exposed hydrophobic surface, and is aggregation prone. Although MPT63 is a complete β‐sheet protein in its native form, the present findings suggest that the secondary structure preferences of the local interactions in early folding pathway may not always follow the native conformation. Furthermore, the Gly25Ala mutant supports the proposed hypothesis by increasing the non‐native helical propensity of the protein structure.     

A comparative molecular force spectroscopy study of homophilic JAM-A interactions and JAM-A interactions with reovirus attachment protein sigma1

JAM-A belongs to a family of immunoglobulin-like proteins called junctional adhesion molecules (JAMs) that localize at epithelial and endothelial intercellular tight junctions. JAM-A is also expressed on dendritic cells, neutrophils, and platelets. Homophilic JAM-A interactions play an important role in regulating paracellular permeability and leukocyte transmigration across epithelial monolayers and endothelial cell junctions, respectively. In addition, JAM-A is a receptor for the reovirus attachment protein, sigma1. In this study, we used single molecular force spectroscopy to compare the kinetics of JAM-A interactions with itself and sigma1. A chimeric murine JAM-A/Fc fusion protein and the purified sigma1 head domain were used to probe murine L929 cells, which express JAM-A and are susceptible to reovirus infection. The bond half-life (t(1/2)) of homophilic JAM-A interactions was found to be shorter (k(off)(o) = 0.688 +/- 0.349 s(-1)) than that of sigma1/JAM-A interactions (k(off)(o) = 0.067 +/- 0.041 s(-1)). These results are in accordance with the physiological functions of JAM-A and sigma1. A short bond lifetime imparts a highly dynamic nature to homophilic JAM-A interactions for regulating tight junction permeability while stable interactions between sigma1 and JAM-A likely anchor the virus to the cell surface and facilitate viral entry.     

Computational diagnosis of protein conformational diseases: short molecular dynamics simulations reveal a fast unfolding of r-LDL mutants that cause familial hypercholesterolemia

The molecular basis of conformational diseases frequently resides in mutant proteins constituting a subset of the vast mutational space. While the subtleties of protein structure point to molecular dynamics (MD) techniques as promising tools for an efficient exploration of such a space, the average size of proteins and the time scale of unfolding events make this goal difficult with present computational capabilities. We show here, nevertheless, that an efficient approach is already feasible for modular proteins. Familial hypercholesterolemia (FH) is a conformational disease linked to mutations in the gene encoding the low density lipoprotein receptor. A high percentage of these mutations has been found in the seven small modular binding repeats of the receptor. Taking advantage of its small size, we have performed an in depth MD study of the fifth binding repeat. Fast unfolding dynamics have been observed in the absence of a structural bound calcium ion, which agrees with its reported essential role in the stability of the module. In addition, several mutations detected in FH patients have been analyzed, starting from the native conformation. Our results indicate that in contrast with the wild type protein and an innocuous control mutant, disease-related mutants experience, in short simulation times (2-8 ns), gross departures from the native state that lead to unfolded conformations and, in some cases, to binding site desorganization deriving in calcium release. Computational diagnosis of mutations leading to conformational diseases seems thus feasible, at least for small or modular pathogenic proteins.     

Direct and tunable modulation of protein levels in rice and wheat with a synthetic small molecule

Direct control of protein level enables rapid and efficient analyses of gene functions in crops. Previously, we developed the RDDK‐Shield1 (Shld1) system in the model plant Arabidopsis thaliana for direct modulation of protein stabilization using a synthetic small molecule. However, it was unclear whether this system is applicable to economically important crops. In this study, we show that the RDDK‐Shld1 system enables rapid and tunable control of protein levels in rice and wheat. Accumulation of RDDK fusion proteins can be reversibly and spatio‐temporally controlled by the synthetic small‐molecule Shld1. Moreover, RDDK‐Bar and RDDK‐Pid3 fusions confer herbicide and rice blast resistance, respectively, in a Shld1‐dependent manner. Therefore, the RDDK‐Shld1 system provides a reversible and tunable technique for controlling protein functions and conditional expression of transgenes in crops.     

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.