Protein unfolding under force: crack propagation in a network

The relative stability of trpzip1 and its mutants determined by computation and experiment

The single-point mutations of tprzip1 are indicated at left, and their relative energetics are compared at right.

Enhancement of protein mechanical stability: Correlated deformations are handcuffed by ligand binding

As revealed by previous experiments, protein mechanical stability can be effectively regulated by ligand binding with the binding site distant from the force-bearing region. However, the mechanism for such long-range allosteric control of protein mechanics is still largely unknown. In this work, we use protein topology-based elastic network model (ENM) and all-atomic steered molecular dynamics (SMD) simulations to study the impact of ligand binding on protein mechanical stability in two systems, i.e., GB1 and CheY-binding P2-domain of CheA (CBDCheA). Both ENM and SMD results show that the ligand binding has considerable and negligible effects on the mechanical stability of these two proteins, respectively. These results are consistent with the experimental observations. A physical mechanism for the enhancement of protein mechanical stability was then proposed: the correlated deformations of the force-bearing region and the binding site are handcuffed by the binding of ligand. The handcuff effect suppresses the propagation of internal force in the force-bearing region, thus improving the resistance to the loading force. Our study indicates that ENM method can effectively identify the structure motifs allosterically related to the deformation in the force bearing region, as well as the force propagation pathway within the structure of the studied proteins. Hence, it should be helpful to understand the molecular origin of the different mechanical properties in response to ligand binding for GB1 and CBDCheA.

The physics of pulling polyproteins: a review of single molecule force spectroscopy using the AFM to study protein unfolding

One of the most exciting developments in the field of biological physics in recent years is the ability to manipulate single molecules and probe their properties and function. Since its emergence over two decades ago, single molecule force spectroscopy has become a powerful tool to explore the response of biological molecules, including proteins, DNA, RNA and their complexes, to the application of an applied force. The force versus extension response of molecules can provide valuable insight into its mechanical stability, as well as details of the underlying energy landscape. In this review we will introduce the technique of single molecule force spectroscopy using the atomic force microscope (AFM), with particular focus on its application to study proteins. We will review the models which have been developed and employed to extract information from single molecule force spectroscopy experiments. Finally, we will end with a discussion of future directions in this field.

Coarse-Grained Simulations of Topology-Dependent Mechanisms of Protein Unfolding and Translocation Mediated by ClpY ATPase Nanomachines

Clp ATPases are powerful ring shaped nanomachines which participate in the degradation pathway of the protein quality control system, coupling the energy from ATP hydrolysis to threading substrate proteins (SP) through their narrow central pore. Repetitive cycles of sequential intra-ring ATP hydrolysis events induce axial excursions of diaphragm-forming central pore loops that effect the application of mechanical forces onto SPs to promote unfolding and translocation. We perform Langevin dynamics simulations of a coarse-grained model of the ClpY ATPase-SP system to elucidate the molecular details of unfolding and translocation of an α/β model protein. We contrast this mechanism with our previous studies which used an all-α SP. We find conserved aspects of unfolding and translocation mechanisms by allosteric ClpY, including unfolding initiated at the tagged C-terminus and translocation via a power stroke mechanism. Topology-specific aspects include the time scales, the rate limiting steps in the degradation pathway, the effect of force directionality, and the translocase efficacy. Mechanisms of ClpY-assisted unfolding and translocation are distinct from those resulting from non-allosteric mechanical pulling. Bulk unfolding simulations, which mimic Atomic Force Microscopy-type pulling, reveal multiple unfolding pathways initiated at the C-terminus, N-terminus, or simultaneously from both termini. In a non-allosteric ClpY ATPase pore, mechanical pulling with constant velocity yields larger effective forces for SP unfolding, while pulling with constant force results in simultaneous unfolding and translocation.

Cell survival is critically dependent on tightly regulated protein quality control, which includes chaperone-mediated folding and degradation. In the degradation pathway, AAA+ nanomachines, such as bacterial Clp proteases, use ATP-driven mechanisms to mechanically unfold, translocate, and destroy excess or defective proteins. Understanding these remodeling mechanisms is of central importance for deciphering the details of essential cellular processes. We perform coarse-grained computer simulations to extensively probe the effect of substrate protein topology on unfolding and translocation actions of the ClpY ATPase nanomachine. We find that, independent of SP topology, unfolding proceeds from the tagged C-terminus, which is engaged by the ATPase, and translocation involves coordinated steps. Topology-specific aspects include more complex unfolding and translocation pathways of the α/β SP compared with the all-α SP due to high stability of β-hairpins and interplay of tertiary contacts. In addition, directionality of the mechanical force applied by the Clp ATPase gives rise to distinct unfolding pathways.

Partial unfolding and refolding for structure refinement: A unified approach of geometric simulations and molecular dynamics

The most successful protein structure prediction methods to date have been template-based modeling (TBM) or homology modeling, which predicts protein structure based on experimental structures. These high accuracy predictions sometimes retain structural errors due to incorrect templates or a lack of accurate templates in the case of low sequence similarity, making these structures inadequate in drug-design studies or molecular dynamics simulations. We have developed a new physics based approach to the protein refinement problem by mimicking the mechanism of chaperons that rehabilitate misfolded proteins. The template structure is unfolded by selectively (targeted) pulling on different portions of the protein using the geometric based technique FRODA, and then refolded using hierarchically restrained replica exchange molecular dynamics simulations (hr-REMD). FRODA unfolding is used to create a diverse set of topologies for surveying near native-like structures from a template and to provide a set of persistent contacts to be employed during re-folding. We have tested our approach on 13 previous CASP targets and observed that this method of folding an ensemble of partially unfolded structures, through the hierarchical addition of contact restraints (that is, first local and then nonlocal interactions), leads to a refolding of the structure along with refinement in most cases (12/13). Although this approach yields refined models through advancement in sampling, the task of blind selection of the best refined models still needs to be solved. Overall, the method can be useful for improved sampling for low resolution models where certain of the portions of the structure are incorrectly modeled.

Cooperativity in thermal and force-induced protein unfolding: integration of crack propagation and network elasticity models

We investigate force-induced and temperature-induced unfolding of proteins using the combination of a gaussian network model and a crack propagation model based on "bond"-breaking independent events. We assume the existence of threshold values for the mean strain and strain fluctuations that dictate bond rupture. Surprisingly, we find that this stepwise process usually leads to a few cooperative, first-order-like, transitions in which several bonds break simultaneously, reminiscent of the "avalanches" seen in disordered networks.

Why is F19Ap53 unable to bind MDM2? Simulations suggest crack propagation modulates binding

Why doesn't the F19A mutant of p53 bind to MDM2? Binding thermodynamics have suggested that the loss of packing interactions upon mutating Phe into Ala sidechain results in destabilizing the binding free energy between p53 and MDM2. Does this mutation also modulate the initial recognition between p53 and MDM2? We look at atomistic computer simulations of the process of the initial encounter between wild type p53 peptide and its F19A mutant with the N-terminal domain of MDM2. These simulations show that binding is characterized by a complex multistep process. It starts with the capture of F19 of wild type p53 by certain residues in the MDM2 binding pocket. This initial step anchors the peptide onto the surface of MDM2, and with the consequent reduction in the search space of the peptide, the peptide docks into the partially occluded surface of MDM2. This is similar to a crack forming in an otherwise occluded hydrophobic cavity in MDM2, and the peptide, docked through F19, modulates the propagation of this crack, which subsequently results in the stepwise docking of the rest of the peptide through insertions of W23 and L26. The lack of the bulky sidechain of F in the F19A mutant results in the absence of the initial "grasp" complex, and hence the mutant peptide diffuses randomly on the surface of MDM2 without binding. This is the first such demonstration of the possibility that a "kinetic" effect may partly underlie the destabilized thermodynamics of binding of F19A and is a feature that appears to be conserved in evolution. The observations by Wallace et al. (Mol Cell 2006; 23:251-63) that despite the inability of F19A to bind at the N-terminal domain of MDM2, it gets ubiquitinated, can now be partly understood based on a mechanism whereby the occupation of the binding pocket by ligands/peptides induces, via crack propagation and the dynamics of gatekeeper Y100, the ubiquitination signal for interactions between the acidic domain of MDM2 and the DNA binding domain of p53.