Fluctuations of primary ubiquitin folding intermediates in a force clamp

Insights into Fluctuations of Structure of Proteins: Significance of Intermediary States in Regulating Biological Functions

Proteins are indispensable to cellular communication and metabolism. The structure on which cells and tissues are developed is deciphered from proteins. To perform functions, proteins fold into a three-dimensional structural design, which is specific and fundamentally determined by their characteristic sequence of amino acids. Few of them have structural versatility, allowing them to adapt their shape to the task at hand. The intermediate states appear momentarily, while protein folds from denatured (D) ⇔ native (N), which plays significant roles in cellular functions. Prolific effort needs to be taken in characterizing these intermediate species if detected during the folding process. Protein folds into its native structure through definite pathways, which involve a limited number of transitory intermediates. Intermediates may be essential in protein folding pathways and assembly in some cases, as well as misfolding and aggregation folding pathways. These intermediate states help to understand the machinery of proper folding in proteins. In this review article, we highlight the various intermediate states observed and characterized so far under in vitro conditions. Moreover, the role and significance of intermediates in regulating the biological function of cells are discussed clearly.

Coarse-Grained Simulation of Mechanical Properties of Single Microtubules With Micrometer Length

Microtubules are one of the most important components in the cytoskeleton and play a vital role in maintaining the shape and function of cells. Because single microtubules are some micrometers long, it is difficult to simulate such a large system using an all-atom model. In this work, we use the newly developed convolutional and K-means coarse-graining (CK-CG) method to establish an ultra-coarse-grained (UCG) model of a single microtubule, on the basis of the low electron microscopy density data of microtubules. We discuss the rationale of the micro-coarse-grained microtubule models of different resolutions and explore microtubule models up to 12-micron length. We use the devised microtubule model to quantify mechanical properties of microtubules of different lengths. Our model allows mesoscopic simulations of micrometer-level biomaterials and can be further used to study important biological processes related to microtubule function.

A new method for the construction of coarse-grained models of large biomolecules from low-resolution cryo-electron microscopy data

The rapid development of cryo-electron microscopy (cryo-EM) has led to the generation of significant low-resolution electron density data of biomolecules. However, the atomistic details of huge biomolecules usually cannot be obtained because it is very difficult to construct all-atom models for MD simulations. Thus, it is still a challenge to make use of the rich low-resolution cryo-EM data for computer simulation and functional study. In this study, we proposed a new method called Convolutional and K-means Coarse-Graining (CK-CG) for the efficient coarse-graining of large biological systems. Using the CK-CG method, we could directly map the cryo-EM data into coarse-grained (CG) beads. Furthermore, the CG beads were parameterized with an empirical harmonic potential to construct a new CG model. We subjected the CK-CG models of the fibrillar protein assemblies F-actin and collagen to external forces in pulling dynamic simulations to assess their mechanical response. The agreement between the estimated tensile stiffness between CG models and experiments demonstrates the validity of the CK-CG method. Thus, our method provides a practical strategy for the direct construction of a structural model from low-resolution data for biological function studies.

Force generation by titin folding

Titin is a giant protein that provides elasticity to muscle. As the sarcomere is stretched, titin extends hierarchically according to the mechanics of its segments. Whether titin's globular domains unfold during this process and how such unfolded domains might contribute to muscle contractility are strongly debated. To explore the force-dependent folding mechanisms, here we manipulated skeletal-muscle titin molecules with high-resolution optical tweezers. In force-clamp mode, after quenching the force (<10 pN), extension fluctuated without resolvable discrete events. In position-clamp experiments, the time-dependent force trace contained rapid fluctuations and a gradual increase of average force, indicating that titin can develop force via dynamic transitions between its structural states en route to the native conformation. In 4 M urea, which destabilizes H-bonds hence the consolidated native domain structure, the net force increase disappeared but the fluctuations persisted. Thus, whereas net force generation is caused by the ensemble folding of the elastically-coupled domains, force fluctuations arise due to a dynamic equilibrium between unfolded and molten-globule states. Monte-Carlo simulations incorporating a compact molten-globule intermediate in the folding landscape recovered all features of our nanomechanics results. The ensemble molten-globule dynamics delivers significant added contractility that may assist sarcomere mechanics, and it may reduce the dissipative energy loss associated with titin unfolding/refolding during muscle contraction/relaxation cycles.

Molecular Plasticity of the Human Voltage-Dependent Anion Channel Embedded Into a Membrane

The voltage-dependent anion channel (VDAC) regulates the flux of metabolites and ions across the outer mitochondrial membrane. Regulation of ion flow involves conformational transitions in VDAC, but the nature of these changes has not been resolved to date. By combining single-molecule force spectroscopy with nuclear magnetic resonance spectroscopy we show that the β barrel of human VDAC embedded into a membrane is highly flexible. Its mechanical flexibility exceeds by up to one order of magnitude that determined for β strands of other membrane proteins and is largest in the N-terminal part of the β barrel. Interaction with Ca(2+), a key regulator of metabolism and apoptosis, considerably decreases the barrel's conformational variability and kinetic free energy in the membrane. The combined data suggest that physiological VDAC function depends on the molecular plasticity of its channel.

Elastic properties and heterogeneous stiffness of the phi29 motor connector channel

The DNA packaging motor of the bacteriophage ϕ29, comprising head-tail connector, ATPase, and pRNA, transports the viral DNA inside the procapsid against pressure differences of up to ∼60 atm during replication. Several models for the DNA packaging mechanism have been proposed, which attribute different roles to the connector, and require specific mechanical properties of the connector. To characterize these properties at the atomic level, and to understand how the connector withstands this large pressure, we have carried out molecular dynamics simulations of the whole connector both in equilibrium and under mechanical stress. The simulations revealed a quite heterogeneous distribution of stiff and soft regions, resembling that of typical composite materials that are also optimized to resist mechanical stress. In particular, the conserved middle α-helical region is found to be remarkably stiff, similar only to structural proteins forming viral shell, silk, or collagen. In contrast, large parts of the peripheral interface to the ϕ29 procapsid turned out to be rather soft. Force probe and umbrella sampling simulations showed that large connector deformations are remarkably reversible, and served to calculate the free energies required for these deformations. In particular, for an untwisting deformation by 12°, as postulated by the untwist-twist model, more than four times' larger energy is required than is available from hydrolysis of one ATP molecule. Combined with previous experiments, this result is incompatible with the untwist-twist model. In contrast, our simulations support the recently proposed one-way revolution model and suggest in structural terms how the connector blocks DNA leakage. In particular, conserved loops at the rim of the central channel, which are in direct contact with the DNA, are found to be rather flexible and tightly anchored to the rigid central region. These findings suggest a check-valve mechanism, with the flexible loops obstructing the channel by interacting with the viral DNA.

Directly investigating the interaction between aptamers and thrombin by atomic force microscopy

Aptamers are single-stranded nucleic acid molecules that can be used for protein recognition, detection, and inhibition. Over the past decades, two thrombin-binding aptamers (15apt and 27apt) were reported by systemic evolution of ligands by exponential enrichment technique. Though many studies have been reported about the interactions between the aptamers and thrombin by atomic force microscopy, the thrombins in those studies were all immobilized by chemical agents. Recently, we developed a new method using atomic force microscopy to directly investigate the specific interactions between thrombin and its two aptamers without immobilizing the thrombin. Furthermore, the unbinding dynamics and dissociation energy landscapes of aptamer/thrombin were discussed. The results indicate that the underlying interaction mechanisms of thrombin with its two aptamers will be similar despite that the structures of 15apt and 27apt are different in buffer solution.

Tracking unfolding and refolding reactions of single proteins using atomic force microscopy methods

During the last two decades single-molecule manipulation techniques such as atomic force microscopy (AFM) has risen to prominence through their unique capacity to provide fundamental information on the structure and function of biomolecules. Here we describe the use of single-molecule AFM to track protein unfolding and refolding pathways, enzymatic catalysis and the effects of osmolytes and chaperones on protein stability and folding. We will outline the principles of operation for two different AFM pulling techniques: length clamp and force-clamp and discuss prominent applications. We provide protocols for the construction of polyproteins which are amenable for AFM experiments, the preparation of different coverslips, choice and calibration of AFM cantilevers. We also discuss the selection criteria for AFM recordings, the calibration of AFM cantilevers, protein sample preparations and analysis of the obtained data.

Cholesterol increases kinetic, energetic, and mechanical stability of the human β2-adrenergic receptor

The steroid cholesterol is an essential component of eukaryotic membranes, and it functionally modulates membrane proteins, including G protein-coupled receptors. To reveal insight into how cholesterol modulates G protein-coupled receptors, we have used dynamic single-molecule force spectroscopy to quantify the mechanical strength and flexibility, conformational variability, and kinetic and energetic stability of structural segments stabilizing the human β(2)-adrenergic receptor (β(2)AR) in the absence and presence of the cholesterol analog cholesteryl hemisuccinate (CHS). CHS considerably increased the kinetic, energetic, and mechanical stability of almost every structural segment at sufficient magnitude to alter the structure and functional relationship of β(2)AR. One exception was the structural core segment of β(2)AR, which establishes multiple ligand binding sites, and its properties were not significantly influenced by CHS.

Structural, energetic, and mechanical perturbations in rhodopsin mutant that causes congenital stationary night blindness

Several point mutations in rhodopsin cause retinal diseases including congenital stationary night blindness and retinitis pigmentosa. The mechanism by which a single amino acid residue substitution leads to dysfunction is poorly understood at the molecular level. A G90D point mutation in rhodopsin causes constitutive activity and leads to congenital stationary night blindness. It is unclear which perturbations the mutation introduces and how they can cause the receptor to be constitutively active. To reveal insight into these mechanisms, we characterized the perturbations introduced into dark state G90D rhodopsin from a transgenic mouse model expressing exclusively the mutant rhodopsin in rod photoreceptor cells. UV-visible absorbance spectroscopy revealed hydroxylamine accessibility to the chromophore-binding pocket of dark state G90D rhodopsin, which is not detected in dark state wild-type rhodopsin but is detected in light-activated wild-type rhodopsin. Single-molecule force spectroscopy suggested that the structural changes introduced by the mutation are small. Dynamic single-molecule force spectroscopy revealed that, compared with dark state wild-type rhodopsin, the G90D mutation decreased energetic stability and increased mechanical rigidity of most structural regions in the dark state mutant receptor. The observed structural, energetic, and mechanical changes in dark state G90D rhodopsin provide insights into the nature of perturbations caused by a pathological point mutation. Moreover, these changed properties observed for dark state G90D rhodopsin are consistent with properties expected for an active state.

Investigation of the interaction between a bivalent aptamer and thrombin by AFM

Aptamers are a new class of molecular probes for protein recognition, detection, and inhibition. Multivalent aptamer-protein binding through aptamer assembly has been currently developed as an effective way to achieve higher protein affinity and selectivity. In this study, the specific interaction between bivalent aptamer Bi-8S and thrombin has been measured directly and quantitatively by atomic force microscopy to investigate the unbinding dynamics and dissociation energy landscape of the multivalent interaction. Bivalent aptamer Bi-8S contains thrombin's two aptamers, 15apt and 27apt, which are linked by eight spacer phosphoramidites. The results revealed the sequential dissociation of the two aptamers. Moreover, the dynamic force spectroscopy data revealed that the 27apt's binding to the thrombin remains largely unaffected by the eight-spacer phosphoramidites within Bi-8S. In contrast, the eight-spacer phosphoramidites stabilized the 15apt-thrombin binding.

Long-lived states to monitor protein unfolding by proton NMR

The relaxation of long-lived states (LLS) corresponds to the slow return to statistical thermal equilibrium between symmetric and antisymmetric proton spin states. This process is remarkably sensitive to the presence of external spins and can be used to obtain information about partial unfolding of proteins. We detected the appearance of a destabilized conformer of ubiquitin when urea is added to the protein in its native state. This conformer shows increased mobility in the C-terminus, which significantly extends the lifetimes of proton LLS magnetisation in Ser-65. These changes could not be detected by conventional measurements of T(1) and T(2) relaxation times of protons, and would hardly be sensed by carbon-13 or nitrogen-15 relaxation measurements. Conformers with similar dynamic and structural features, as revealed by LLS relaxation times, could be observed, in the absence of urea, in two ubiquitin mutants, L67S and L69S.

Locating an extracellular K+-dependent interaction site that modulates betaine-binding of the Na+-coupled betaine symporter BetP

BetP, a trimeric Na(+)-coupled betaine symporter, senses hyperosmotic stress via its cytoplasmic C-terminal domain and regulates transport activity in dependence of the cytoplasmic K(+)-concentration. This transport regulation of BetP depends on a sophisticated interaction network. Using single-molecule force spectroscopy we structurally localize and quantify these interactions changing on K(+)-dependent transport activation and substrate-binding. K(+) significantly strengthened all interactions, modulated lifetimes of functionally important structural regions, and increased the mechanical rigidity of the symporter. Substrate-binding could modulate, but not establish most of these K(+)-dependent interactions. A pronounced effect triggered by K(+) was observed at the periplasmic helical loop EH2. Tryptophan quenching experiments revealed that elevated K(+)-concentrations akin to those BetP encounters during hyperosmotic stress trigger the formation of a periplasmic second betaine-binding (S2) site, which was found to be at a similar position reported previously for the BetP homologue CaiT. In BetP, the presence of the S2 site strengthened the interaction between EH2, transmembrane α-helix 12 and the K(+)-sensing C-terminal domain resulting in a K(+)-dependent cooperative betaine-binding.

Velocity-dependent mechanical unfolding of bacteriorhodopsin is governed by a dynamic interaction network

Bacteriorhodopsin is a model system for membrane proteins. This seven transmembrane helical protein is embedded within a membrane structure called purple membrane. Its structural stability against mechanical stress was recently investigated by atomic force microscopy experiments, in which single proteins were extracted from the purple membrane. Here, we study this process by all-atom molecular dynamics simulations, in which single bacteriorhodopsin molecules were extracted and unfolded from an atomistic purple membrane model. In our simulations, key features from the experiments like force profiles and location of key residues that resist mechanical unfolding were reproduced. These key residues were seen to be stabilized by a dynamic network of intramolecular interactions. Further, the unfolding pathway was found to be velocity-dependent. Simulations in which the mechanical stress was released during unfolding revealed relaxation motions that allowed characterization of the nonequilibrium processes during fast extraction.

Minimum energy compact structures in force-quench polyubiquitin folding are domain swapped

Single molecule experiments that initiate folding using mechanical force are uniquely suited to reveal the nature of populated states in the folding process. Using a strategy proposed on theoretical grounds, which calls for repeated cycling of force from high to low values using force pulses, it was demonstrated in atomic force spectroscopy (AFM) experiments that an ensemble of minimum energy compact structures (MECS) are sampled during the folding of polyubiquitin. The structures in the ensemble are mechanically resistant to a lesser extent than the native state. Remarkably, forced unfolding of the populated intermediates reveals a broad distribution of extensions including steps up to 30 nm and beyond. We show using molecular simulations that favorable interdomain interactions leading to domain swapping between adjacent ubiquitin modules results in the formation of the ensemble of MECS, whose unfolding leads to an unusually broad distribution of steps. We obtained the domain-swapped structures using coarse-grained ubiquitin dimer models by exchanging native interactions between two monomeric ubiquitin molecules. Brownian dynamics force unfolding of the proposed domain-swapped structures, with mechanical stability that is approximately 100-fold lower than the native state, gives rise to a distribution of extensions from 2 to 30 nm. Our results, which are in quantitative agreement with AFM experiments, suggest that domain swapping may be a general mechanism in the assembly of multi-sub-unit proteins.

Dual energy landscape: the functional state of the β-barrel outer membrane protein G molds its unfolding energy landscape

We applied dynamic single-molecule force spectroscopy to quantify the parameters (free energy of activation and distance of the transition state from the folded state) characterizing the energy barriers in the unfolding energy landscape of the outer membrane protein G (OmpG) from Escherichia coli. The pH-dependent functional switching of OmpG directs the protein along different regions on the unfolding energy landscape. The two functional states of OmpG take the same unfolding pathway during the sequential unfolding of β-hairpins I-IV. After the initial unfolding events, the unfolding pathways diverge. In the open state, the unfolding of β-hairpin V in one step precedes the unfolding of β-hairpin VI. In the closed state, β-hairpin V and β-strand S11 with a part of extracellular loop L6 unfold cooperatively, and subsequently β-strand S12 unfolds with the remaining loop L6. These two unfolding pathways in the open and closed states join again in the last unfolding step of β-hairpin VII. Also, the conformational change from the open to the closed state witnesses a rigidified extracellular gating loop L6. Thus, a change in the conformational state of OmpG not only bifurcates its unfolding pathways but also tunes its mechanical properties for optimum function.

The effect of different force applications on the protein-protein complex Barnase-Barstar

Steered molecular dynamics simulations are a tool to examine the energy landscape of protein-protein complexes by applying external forces. Here, we analyze the influence of the velocity and geometry of the probing forces on a protein complex using this tool. With steered molecular dynamics, we probe the stability of the protein-protein complex Barnase-Barstar. The individual proteins are mechanically labile. The Barnase-Barstar binding site is more stable than the folds of the individual proteins. By using different force protocols, we observe a variety of responses of the system to the applied tension.

Periodic forces trigger a complex mechanical response in ubiquitin

Mechanical forces govern physiological processes in all living organisms. Many cellular forces, for example, those generated in cyclic conformational changes of biological machines, have repetitive components. In apparent contrast, little is known about how dynamic protein structures respond to periodic mechanical information. Ubiquitin is a small protein found in all eukaryotes. We developed molecular dynamics simulations to unfold single and multimeric ubiquitins with periodic forces. By using a coarse-grained representation, we were able to model forces with periods about 2 orders of magnitude longer than the protein's relaxation time. We found that even a moderate periodic force weakened the protein and shifted its unfolding pathways in a frequency- and amplitude-dependent manner. A complex dynamic response with secondary structure refolding and an increasing importance of local interactions was revealed. Importantly, repetitive forces with broadly distributed frequencies elicited very similar molecular responses compared to fixed-frequency forces. When testing the influence of pulling geometry on ubiquitin's mechanical stability, it was found that the linkage involved in the mechanical degradation of cellular proteins renders the protein remarkably insensitive to periodic forces. We also devised a complementary kinetic energy landscape model that traces these observations and explains periodic-force, single-molecule measurements. In turn, this analytical model is capable of predicting dynamic protein responses. These results provide new insights into ubiquitin mechanics and a potential mechanical role during protein degradation, as well as first frameworks for dynamic protein stability and the modeling of repetitive mechanical processes.

Substrate binding tunes conformational flexibility and kinetic stability of an amino acid antiporter

We used single molecule dynamic force spectroscopy to unfold individual serine/threonine antiporters SteT from Bacillus subtilis. The unfolding force patterns revealed interactions and energy barriers that stabilized structural segments of SteT. Substrate binding did not establish strong localized interactions but appeared to be facilitated by the formation of weak interactions with several structural segments. Upon substrate binding, all energy barriers of the antiporter changed thereby describing the transition from brittle mechanical properties of SteT in the unbound state to structurally flexible conformations in the substrate-bound state. The lifetime of the unbound state was much shorter than that of the substrate-bound state. This leads to the conclusion that the unbound state of SteT shows a reduced conformational flexibility to facilitate specific substrate binding and a reduced kinetic stability to enable rapid switching to the bound state. In contrast, the bound state of SteT showed an increased conformational flexibility and kinetic stability such as required to enable transport of substrate across the cell membrane. This result supports the working model of antiporters in which alternate substrate access from one to the other membrane surface occurs in the substrate-bound state.

The protein folding problem

The "protein folding problem" consists of three closely related puzzles: (a) What is the folding code? (b) What is the folding mechanism? (c) Can we predict the native structure of a protein from its amino acid sequence? Once regarded as a grand challenge, protein folding has seen great progress in recent years. Now, foldable proteins and nonbiological polymers are being designed routinely and moving toward successful applications. The structures of small proteins are now often well predicted by computer methods. And, there is now a testable explanation for how a protein can fold so quickly: A protein solves its large global optimization problem as a series of smaller local optimization problems, growing and assembling the native structure from peptide fragments, local structures first.

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.