The mechanical stability of ubiquitin is linkage dependent

Ubiquitin chains are formed through the action of a set of enzymes that covalently link ubiquitin either through peptide bonds or through isopeptide bonds between their C terminus and any of four lysine residues. These naturally occurring polyproteins allow one to study the mechanical stability of a protein, when force is applied through different linkages. Here we used single-molecule force spectroscopy techniques to examine the mechanical stability of N-C-linked and Lys48-C-linked ubiquitin chains. We combined these experiments with steered molecular dynamics (SMD) simulations and found that the mechanical stability and unfolding pathway of ubiquitin strongly depend on the linkage through which the mechanical force is applied to the protein. Hence, a protein that is otherwise very stable may be easily unfolded by a relatively weak mechanical force applied through the right linkage. This may be a widespread mechanism in biological systems.     

Mutagenesis data in the automated prediction of transmembrane helix dimers

We present a molecular modeling protocol that selects modeled protein structures based on experimental mutagenesis results. The computed effect of a point mutation should be consistent with its experimental effect for correct models; mutations that do not affect protein stability and function should not affect the computed energy of a correct model while destabilizing mutations should have unfavorable computed energies. On the other hand, an incorrect model will likely display computed energies that are inconsistent with experimental results. We added terms to our energy function which penalize models that are inconsistent with experimental results. This creates a selective advantage for models that are consistent with experimental results in the Monte Carlo simulated annealing protocol we use to search conformational space. We calibrated our protocol to predict the structure of transmembrane helix dimers using glycophorin A as a model system. Inclusion of mutational data in this protocol compensates for the limitations of our force field and the limitations of our conformational search. We demonstrate an application of this structure prediction protocol by modeling the transmembrane region of the BNIP3 apoptosis factor.     

Force-clamp spectroscopy of single-protein monomers reveals the individual unfolding and folding pathways of I27 and ubiquitin

Single-protein force experiments have relied on a molecular fingerprint based on tethering multiple single-protein domains in a polyprotein chain. However, correlations between these domains remain an issue in interpreting force spectroscopy data, particularly during protein folding. Here we first show that force-clamp spectroscopy is a sensitive technique that provides a molecular fingerprint based on the unfolding step size of four single-monomer proteins. We then measure the force-dependent unfolding rate kinetics of ubiquitin and I27 monomers and find a good agreement with the data obtained for the respective polyproteins over a wide range of forces, in support of the Markovian hypothesis. Moreover, with a large statistical ensemble at a single force, we show that ubiquitin monomers also exhibit a broad distribution of unfolding times as a signature of disorder in the folded protein landscape. Furthermore, we readily capture the folding trajectories of monomers that exhibit the same stages in folding observed for polyproteins, thus eliminating the possibility of entropic masking by other unfolded modules in the chain or domain-domain interactions. On average, the time to reach the I27 folded length increases with increasing quenching force at a rate similar to that of the polyproteins. Force-clamp spectroscopy at the single-monomer level reproduces the kinetics of unfolding and refolding measured using polyproteins, which proves that there is no mechanical effect of tethering proteins to one another in the case of ubiquitin and I27.     

A novel culture morphology resulting from applied mechanical strain

Summary To demonstrate that cells both perceive and respond to external force, a strain/relaxation regimen was applied to normal human fetal and aged dermal fibroblasts cultured as monolayers on flexible membranes. The precisely controlled protocol of stretch (20% elongation of the culture membrane) at 6.67 cycles/min caused a progressive change in the monolayers, such that the original randomly distributed pattern of cells became a symmetric, radial distribution as the cell bodies aligned parallel to the applied force. High cell density interfered with the success of re-alignment in the fetal cell cultures observed, which may reflect a preference in this cell strain for cell-cell over cell-matrix contacts. The chronologically aged cells observed did not demonstrate this feature, aligning efficiently at all seeding densities examined. The role of microfilaments in force perception and transmission was investigated through the addition of cytochalasin D in graded doses. Both intercellular interactions and cytoskeletal integrity mediate the morphological response to mechanical strain.     

Biomimetic systems for studying actin-based motility

Actin polymerization provides a major driving force for eukaryotic cell motility. Successive intercalation of monomeric actin subunits between the plasma membrane and the filamentous actin network results in protrusions of the membrane enabling the cell to move or to change shape. One of the challenges in understanding eukaryotic cell motility is to dissect the elementary biochemical and biophysical steps that link actin polymerization to mechanical force generation. Recently, significant progress was made using biomimetic, in vitro systems that are inspired by the actin-based motility of bacterial pathogens such as Listeria monocytogenes. Polystyrene microspheres and synthetic phospholipid vesicles coated with proteins that initiate actin polymerization display motile behavior similar to Listeria, mimicking the leading edge of lamellipodia and filopodia. A major advantage of these biomimetic systems is that both biochemical and physical parameters can be controlled precisely. These systems provide a test bed for validating theoretical models on force generation and polarity establishment resulting from actin polymerization. In this review, we discuss recent experimental progress using biomimetic systems propelled by actin polymerization and discuss these results in the light of recent theoretical models on actin-based motility.     

The backbone conformational entropy of protein folding: experimental measures from atomic force microscopy

The energy dissipated during the atomic force microscopy-based mechanical unfolding and extension of proteins is typically an order of magnitude greater than their folding free energy. The vast majority of the "excess" energy dissipated is thought to arise due to backbone conformational entropy losses as the solvated, random-coil unfolded state is stretched into an extended, low-entropy conformation. We have investigated this hypothesis in light of recent measurements of the energy dissipated during the mechanical unfolding of "polyproteins" comprised of multiple, homogeneous domains. Given the assumption that backbone conformational entropy losses account for the vast majority of the energy dissipated (an assumption supported by numerous lines of experimental evidence), we estimate that approximately 19(+/-2)J/(mol K residue) of entropy is lost during the extension of three mechanically stable beta-sheet polyproteins. If, as suggested by measured peak-to-peak extension distances, pulling proceeds to near completion, this estimate corresponds to the absolute backbone conformational entropy of the unfolded state. As such, it is exceedingly close to previous theoretical and semi-empirical estimates that place this value at approximately 20J/(mol K residue). The estimated backbone conformational entropy lost during the extension of two helical polyproteins, which, in contrast to the mechanically stable beta-sheet polyproteins, rupture at very low applied forces, is three- to sixfold less. Either previous estimates of the backbone conformational entropy are significantly in error, or the reduced mechanical strength of the helical proteins leads to the rupture of a subsequent domain before full extension (and thus complete entropy loss) is achieved.     

Engineering of an elastic scaffolding polyprotein based on an SH3-binding intrinsically disordered titin PEVK module

Titin is a large elastic protein found in muscle that maintains the elasticity and structural integrity of the sarcomere. The PEVK region of titin is intrinsically disordered, highly elastic and serves as a hub to bind signaling proteins. Systematic investigation of the structure and affinity profile of the PEVK region will provide important information about the functions of titin. Since PEVK is highly heterogeneous due to extensive differential splicing from more than one hundred exons, we engineered and expressed polyproteins that consist of a defined number of identical single exon modules. These customized polyproteins reduce heterogeneity, amplify interactions of less dominant modules, and most importantly, provide tags for atomic force microscopy and allow more readily interpretable data from single-molecule techniques. Expression and purification of recombinant polyprotein with repeat regions presented many technical challenges: recombination events in tandem repeats of identical DNA sequences exacerbated by high GC content, toxicity of polymer plasmid and expressed protein to the bacteria; early truncation of proteins expressed with different numbers of modules; and extreme sensitivity to proteolysis. We have investigated a number of in vitro and in vivo bacterial and yeast expression systems, as well as baculoviral systems as potential solutions to these problems. We successfully expressed and purified in gram quantities a polyprotein derived from human titin exon 172 using Pichia pastoris yeast. This study provides valuable insights into the technical challenges regarding the engineering and purification of a tandem repeat sequence of an intrinsically disordered biopolymer.     

Polyprotein Processing as a Strategy for Gene Expression in RNA Viruses

RNA viruses in many families and genera express their genomes in ways which involve the synthesis and subsequent cleavage of precursor polyproteins. This strategem allows the activation of subsets of proteins with different biochemical functions from the same precursor protein. Although the virus-encoded enzymes responsible for processing the polyproteins are structurally diverse, they are all highly specific for their substrates. The resulting processing cascade is a tightly controlled process, which in several cases involves the action of protein cofactors to modulate the activity of the proteinase.     

A Brownian Dynamics Model of Kinesin in Three Dimensions Incorporating the Force-Extension Profile of the Coiled-Coil Cargo Tether

The kinesin family of motor proteins are involved in a variety of cellular processes that transport materials and generate force. With recent advances in experimental techniques, such as optical tweezers can probe individual molecules, there has been an increasing interest in understanding the mechanisms by which motor proteins convert chemical energy into mechanical work. Here we present a mathematical model for the chemistry and three dimensional mechanics of the kinesin motor protein which captures many of the force dependent features of the motor. For the elasticity of the tether that attaches cargo to the motor we develop a method for deriving the non-linear force-extension relationship from optical trap data. For the kinesin heads, cargo, and microscope stage we formulate a three dimensional Brownian Dynamics model that takes into account excluded volume interactions. To efficiently compute statistics from the model, an algorithm is proposed which uses a two step protocol that separates the simulation of the mechanical features of the model from the chemical kinetics of the model. Using this approach for a bead transported by the motor, the force dependent average velocity and randomness parameter are computed and compared with the experimental data.     

Tuning the mechanical stability of fibronectin type III modules through sequence variations

Cells can switch the functional states of extracellular matrix proteins by stretching them while exerting mechanical force. Using steered molecular dynamics, we investigated how the mechanical stability of FnIII modules from the cell adhesion protein fibronectin is affected by natural variations in their amino acid sequences. Despite remarkably similar tertiary structures, FnIII modules share low sequence homology. Conversely, the sequence homology for the same FnIII module across multiple species is notably higher, suggesting that sequence variability is functionally significant. Our studies find that the mechanical stability of FnIII modules can be tuned through substitutions of just a few key amino acids by altering access of water molecules to hydrogen bonds that break early in the unfolding pathway. Furthermore, the FnIII hierarchy of mechanical unfolding can be changed by environmental conditions, such as pH for FnIII10, or by forming complexes with other molecules, such as heparin binding to FnIII13.     

Fast and easy protocol for the purification of recombinant S-layer protein for synthetic biology applications

A goal of synthetic biology is to make biological systems easier to engineer. One of the aims is to design, with nanometer-scale precision, biomaterials with well-defined properties. The surface-layer protein SbpA forms 2D arrays naturally on the cell surface of Lysinibacillus sphaericus, but also as the purified protein in solution upon the addition of divalent cations. The high propensity of SbpA to form crystalline arrays, which can be simply controlled by divalent cations, and the possibility to genetically alter the protein, make SbpA an attractive molecule for synthetic biology. To be a useful tool, however, it is important that a simple protocol can be used to produce recombinant wild-type and modified SbpA in large quantities and in a biologically active form. The present study addresses this requirement by introducing a mild and non-denaturing purification protocol to produce milligram quantities of recombinant, active SbpA.     

The mechanical hierarchies of fibronectin observed with single-molecule AFM

Mechanically induced conformational changes in proteins such as fibronectin are thought to regulate the assembly of the extracellular matrix and underlie its elasticity and extensibility. Fibronectin contains a region of tandem repeats of up to 15 type III domains that play critical roles in cell binding and self-assembly. Here, we use single-molecule force spectroscopy to examine the mechanical properties of fibronectin (FN) and its individual FNIII domains. We found that fibronectin is highly extensible due to the unfolding of its FNIII domains. We found that the native FNIII region displays strong mechanical unfolding hierarchies requiring 80 pN of force to unfold the weakest domain and 200 pN for the most stable domain. In an effort to determine the identity of the weakest/strongest domain, we engineered polyproteins composed of an individual domain and measured their mechanical stability by single-protein atomic force microscopy (AFM) techniques. In contrast to chemical and thermal measurements of stability, we found that the tenth FNIII domain is mechanically the weakest and that the first and second FNIII domains are the strongest. Moreover, we found that the first FNIII domain can acquire multiple, partially folded conformations, and that their incidence is modulated strongly by its neighbor FNIII domain. The mechanical hierarchies of fibronectin demonstrated here may be important for the activation of fibrillogenesis and matrix assembly.     

A screening strategy for heterologous protein expression in Escherichia coli with the highest return of investment

Heterologous protein expression in Escherichia coli is commonly used to obtain recombinant proteins for a variety of downstream applications. However, many proteins are not, or are only poorly, expressed in soluble form. High level expression often leads to the formation of inclusion bodies and an inactive product that needs to be refolded. By screening the solubility pattern for a set of 71 target proteins in different host-strains and varying parameters such as location of purification tag, promoter and induction temperature we propose a protocol with a success rate of 77% of clones returning a soluble protein. This protocol is particularly suitable for high-throughput screening with the goal to obtain soluble protein product for e.g. structure determination.     

On the interpretation of biochemical data by molecular dynamics computer simulation.

The application of computer simulation to molecular systems of biochemical interest is reviewed. It is shown that computer simulation is a tool complementary to experimental methods, which can be used to access atomic details inaccessible to experimental probes. Examples are given in which computer simulation augments the experimental information by providing an atomic picture of high resolution with respect to space, energy or time. The usefulness of a computer simulation largely depends on its quality. The most important factors that limit the accuracy of simulated results are discussed. The accuracy of different simulation studies can differ by orders of magnitude. The accuracy will depend on the type of biomolecular system and process studied. It will also depend on the choice of force field, the simulation set-up and the protocol that is used. A list of quality-determining factors is given, which may be useful when interpreting simulation studies appearing in the literature.     

Probing protein mechanics: residue-level properties and their use in defining domains

It is becoming clear that, in addition to structural properties, the mechanical properties of proteins can play an important role in their biological activity. It nevertheless remains difficult to probe these properties experimentally. Whereas single-molecule experiments give access to overall mechanical behavior, notably the impact of end-to-end stretching, it is currently impossible to directly obtain data on more local properties. We propose a theoretical method for probing the mechanical properties of protein structures at the single-amino acid level. This approach can be applied to both all-atom and simplified protein representations. The probing leads to force constants for local deformations and to deformation vectors indicating the paths of least mechanical resistance. It also reveals the mechanical coupling that exists between residues. Results obtained for a variety of proteins show that the calculated force constants vary over a wide range. An analysis of the induced deformations provides information that is distinct from that obtained with measures of atomic fluctuations and is more easily linked to residue-level properties than normal mode analyses or dynamic trajectories. It is also shown that the mechanical information obtained by residue-level probing opens a new route for defining so-called dynamical domains within protein structures.     

Recognition Imaging of Acetylated Chromatin Using a DNA Aptamer

Histone acetylation plays an important role in the regulation of gene expression. A DNA aptamer generated by in vitro selection to be highly specific for histone H4 protein acetylated at lysine 16 was used as a recognition element for atomic force microscopy-based recognition imaging of synthetic nucleosomal arrays with precisely controlled acetylation. The aptamer proved to be reasonably specific at recognizing acetylated histones, with recognition efficiencies of 60% on-target and 12% off-target. Though this selectivity is much poorer than the >2000:1 equilibrium specificity of the aptamer, it is a large improvement on the performance of a ChIP-quality antibody, which is not selective at all in this application, and it should permit high-fidelity recognition with repeated imaging. The ability to image the precise location of posttranslational modifications may permit nanometer-scale investigation of their effect on chromatin structure.     

Production of recombinant protein in Pichia pastoris by fermentation

This protocol is applicable to recombinant protein expression by small-scale fermentation using the Pichia pastoris expression system. P. pastoris has the capacity to produce large quantities of protein with eukaryotic processing. Expression is controlled by a methanol-inducible promoter, which allows a biomass-generation phase before protein production is initiated. The target protein is secreted directly into a protein-free mineral salt medium, and is relatively easy to purify. The protocol is readily interfaced with expanded bed adsorption for immediate capture and purification of recombinant protein. The setting up of the bioreactor plus the fermentation itself takes 1 wk. Making the master and user seed lots takes approximately 2 wk for each individual clone.     

A Web Interface for Easy Flexible Protein-Protein Docking with ATTRACT

Protein-protein docking programs can give valuable insights into the structure of protein complexes in the absence of an experimental complex structure. Web interfaces can facilitate the use of docking programs by structural biologists. Here, we present an easy web interface for protein-protein docking with the ATTRACT program. While aimed at nonexpert users, the web interface still covers a considerable range of docking applications. The web interface supports systematic rigid-body protein docking with the ATTRACT coarse-grained force field, as well as various kinds of protein flexibility. The execution of a docking protocol takes up to a few hours on a standard desktop computer.     

Computational investigation of the effect of thermal perturbation on the mechanical unfolding of titin I27

The emergence of single-molecule force measurement experiments has facilitated a better understanding of protein folding pathways and the thermodynamics involved. Computational methods such as steered molecular dynamics (SMD) simulations are helpful in providing atomistic level information on the unfolding pathways. Recent experimental studies have showed that combinations of single-molecule experiments with traditional methods such as chemical and/or thermal denaturation yield additional insights into the folding phenomenon. In this study, we report results from extensive computations (a total of about 60 SMD simulations with a total length of about 0.4 μs) that address the effect of thermal perturbation on the mechanical stability of the I27 domain of the protein titin. A wide range of temperatures (280-340 K) were considered for the pulling, which was done at both constant velocity and constant force using SMD simulations. Good agreement with experimental data, such as for the trends in changes in average force and the maximum force with respect to the temperature, was obtained. This study identifies two competing pathways for the mechanical unfolding of I27, and illustrates the significance of combining various techniques to examine protein folding.     

Naturally occurring osmolytes modulate the nanomechanical properties of polycystic kidney disease domains

Polycystin-1 (PC1) is a large membrane protein that is expressed along the renal tubule and exposed to a wide range of concentrations of urea. Urea is known as a common denaturing osmolyte that affects protein function by destabilizing their structure. However, it is known that the native conformation of proteins can be stabilized by protecting osmolytes that are found in the mammalian kidney. PC1 has an unusually long ectodomain with a multimodular structure including 16 Ig-like polycystic kidney disease (PKD) domains. Here, we used single-molecule force spectroscopy to study directly the effects of several naturally occurring osmolytes on the mechanical properties of PKD domains. This experimental approach more closely mimics the conditions found in vivo. We show that upon increasing the concentration of urea there is a remarkable decrease in the mechanical stability of human PKD domains. We found that protecting osmolytes such as sorbitol and trimethylamine N-oxide can counteract the denaturing effect of urea. Moreover, we found that the refolding rate of a structurally homologous archaeal PKD domain is significantly slowed down in urea, and this effect was counteracted by sorbitol. Our results demonstrate that naturally occurring osmolytes can have profound effects on the mechanical unfolding and refolding pathways of PKD domains. Based on these findings, we hypothesize that osmolytes such as urea or sorbitol may modulate PC1 mechanical properties and may lead to changes in the activation of the associated polycystin-2 channel or other intracellular events mediated by PC1.