Role of pulling direction in understanding the energy landscape of proteins

Contribution of Hydrophobic Interactions to Protein Mechanical Stability

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

Biaxial Dielectrophoresis Force Spectroscopy: A Stoichiometric Approach for Examining Intermolecular Weak Binding Interactions

The direct quantification of weak intermolecular binding interactions is very important for many applications in biology and medicine. Techniques that can be used to investigate such interactions under a controlled environment, while varying different parameters such as loading rate, pulling direction, rupture event measurements, and the use of different functionalized probes, are still lacking. Herein, we demonstrate a biaxial dielectrophoresis force spectroscopy (BDFS) method that can be used to investigate weak unbinding events in a high-throughput manner under controlled environments and by varying the pulling direction (i.e., transverse and/or vertical axes) as well as the loading rate. With the BDFS system, we can quantitatively analyze binding interactions related to hydrogen bonding or ionic attractions between functionalized microbeads and a surface within a microfluidic device. Our BDFS system allowed for the characterization of the number of bonds involved in an interaction, bond affinity, kinetic rates, and energy barrier heights and widths from different regimes of the energy landscape.

Nanomechanics of β-rich proteins related to neuronal disorders studied by AFM, all-atom and coarse-grained MD methods

Computer simulations of protein unfolding substantially help to interpret force-extension curves measured in single-molecule atomic force microscope (AFM) experiments. Standard all-atom (AA) molecular dynamics simulations (MD) give a good qualitative mechanical unfolding picture but predict values too large for the maximum AFM forces with the common pulling speeds adopted here. Fine tuned coarse-grain MD computations (CG MD) offer quantitative agreement with experimental forces. In this paper we address an important methodological aspect of MD modeling, namely the impact of numerical noise generated by random assignments of bead velocities on maximum forces (F_max) calculated within the CG MD approach. Distributions of CG forces from 2000 MD runs for several model proteins rich in β structures and having folds with increasing complexity are presented. It is shown that F_max have nearly Gaussian distributions and that values of F_max for each of those β-structures may vary from 93.2 ± 28.9 pN (neurexin) to 198.3 ± 25.2 pN (fibronectin). The CG unfolding spectra are compared with AA steered MD data and with results of our AFM experiments for modules present in contactin, fibronectin and neurexin. The stability of these proteins is critical for the proper functioning of neuronal synaptic clefts. Our results confirm that CG modeling of a single molecule unfolding is a good auxiliary tool in nanomechanics but large sets of data have to be collected before reliable comparisons of protein mechanical stabilities are made.

Figure

Experimental and computational characterization of biological liquid crystals: a review of single-molecule bioassays

Quantitative understanding of the mechanical behavior of biological liquid crystals such as proteins is essential for gaining insight into their biological functions, since some proteins perform notable mechanical functions. Recently, single-molecule experiments have allowed not only the quantitative characterization of the mechanical behavior of proteins such as protein unfolding mechanics, but also the exploration of the free energy landscape for protein folding. In this work, we have reviewed the current state-of-art in single-molecule bioassays that enable quantitative studies on protein unfolding mechanics and/or various molecular interactions. Specifically, single-molecule pulling experiments based on atomic force microscopy (AFM) have been overviewed. In addition, the computational simulations on single-molecule pulling experiments have been reviewed. We have also reviewed the AFM cantilever-based bioassay that provides insight into various molecular interactions. Our review highlights the AFM-based single-molecule bioassay for quantitative characterization of biological liquid crystals such as proteins.