Energy Landscape of Ubiquitin Is Weakly Multidimensional

The Effect of Pulling and Twisting Forces on Chameleon Sequence Peptides

Chameleon sequences are amino acid sequences found in several distinct configurations in experiment. They challenge our understanding of the link between sequence and structure, and provide insight into structural competition in proteins. Here, we study the energy landscapes for three such sequences, and interrogate how pulling and twisting forces impact the available structural ensembles. Chameleon sequences do not necessarily exhibit multiple structural ensembles on a multifunnel energy landscape when we consider them in isolation. The application of even small forces leads to drastic changes in the energy landscapes. For pulling forces, we observe transitions from helical to extended structures in a very small span of forces. For twisting forces, the picture is much more complex, and highly dependent on the magnitude and handedness of the applied force as well as the reference angle for the twist. Depending on these parameters, more complex and more simplistic energy landscapes are observed alongside more and less diverse structural ensembles. The impact of even small forces is significant, confirming their likely role in folding events. In addition, small forces exerted by the remaining scaffold of a protein may be sufficient to lead to the adoption of a specific structural ensemble by a chameleon sequence.

Conformational Fluctuations and Induced Orientation of a Protein, Its Solvation Shell, and Bulk Water in Weak Non-Unfolding External Electric Fields

Extreme external electric fields have been reported to disrupt the tertiary structure of stably folded proteins; however, the effects of weaker electric fields on many biomolecules, especially net-uncharged proteins, and on the surrounding aqueous environment have been rarely discussed. To explore these effects at the atomic level, here, we have used molecular dynamics simulations to estimate rotational motion and induced structural fluctuations in the model protein ubiquitin and its hydration layer due to applied non-unfolding electrostatic fields. When exposed to weak electric fields of up to 0.2 V nm-1, ubiquitin displayed competition between internal structure-maintaining molecular interactions and the external orienting force, which disrupted the local structure in certain regions of the protein. Moreover, relative to hydration water, bulk water showed a greater tendency to align with the electric field, indicating that the presence of protein caused hydration water to acquire rotational mobility different from that in a pure-water system. The differential influence of the applied electric field on the hydration and bulk water surrounding ubiquitin will be common to almost all (nonmembrane) biomacromolecules. Our findings highlight the importance of local dipoles and their electric polarizability even in net-uncharged biomolecules.

Molecular interpretation of single-molecule force spectroscopy experiments with computational approaches

Single molecule force-spectroscopy techniques have granted access to unprecedented molecular-scale details about biochemical and biological mechanisms. However, the interpretation of the experimental data is often challenging. Computational and simulation approaches (all-atom steered MD simulations in particular) are key to provide molecular details about the associated mechanisms, to help test different hypotheses and to predict experimental results. In this review, particular recent efforts directed towards the molecular interpretation of single-molecule force spectroscopy experiments on proteins and protein-related systems (often in close collaboration with experimental groups) will be presented. These results will be discussed in the broader context of the field, highlighting the recent achievements and the ongoing challenges for computational biophysicists and biochemists. In particular, I will focus on the input gained from molecular simulations approaches to rationalize the origin of the unfolded protein elasticity and the protein conformational behavior under force, to understand how force denaturation differs from chemical, thermal or shear unfolding, and to unravel the molecular details of unfolding events for a variety of systems. I will also discuss the use of models based on Langevin dynamics on a 1-D free-energy surface to understand the effect of protein segmentation on the work exerted by a force, or, at the other end of the spectrum of computational techniques, how quantum calculations can help to understand the reactivity of disulfide bridges exposed to force.

Native Salt Bridges Are a Key Regulator of Ubiquitin's Mechanical Stability

Although it is known that various intramolecular interactions determine protein mechanical stability, a detailed molecular-level understanding of the key regulators of protein mechanical stability is still lacking. Here, we present evidence for salt bridges in ubiquitin as important intramolecular interactions that can affect protein mechanical stability. Ubiquitin has two salt bridges: one relatively surface-exposed (SB1:K11-E34) and the other relatively buried (SB2:K27-D52). Ubiquitin is a reversible post-translational modifier and is stable mechanically (F_avg^u = 185 pN). On breaking SB1, the mechanical stability of ubiquitin is slightly enhanced (F_avg^u = 193 pN). In contrast, the mechanical stability significantly decreased upon breaking SB2 (F_avg^u = 158 pN). These results suggest that SB1 are SB2 are regulators of the mechanical stability of ubiquitin. Interestingly, the mechanical stability decreased further (F_avg^u = 145 pN) for the double salt bridge (DB) null variant. Monte Carlo simulations elucidate that the main regulating factor is the spontaneous unfolding rate constant (k_u⁰), being the highest for the DB null variant followed by the SB2 null variant, and it remains unaltered for the SB1 null variant, while the native-to-transition-state distance (x_u) remains unchanged. Our study provides mechanistic understanding on how two native salt bridges can independently regulate the mechanical stability in a protein, which has implications in designing protein-based robust biomaterials in the future.

Reconstruction of the Free Energy Profile for SUMO1 from Nonequilibrium Single-Molecule Pulling Experiments

Free energy profiles form the cornerstone in the study of protein folding and function. In this study, the free energy profile of SUMO1 protein is directly reconstructed using an extension of the Jarzynski equality from atomic force microscope (AFM) based single-molecule force spectroscopy (SMFS) experiments. SUMO1 is a ubiquitin-like posttranslational modifier protein having a β clamp motif in its structure, imparting it with mechanical stability. We use the Jarzynski equality to obtain the equilibrium free energy profile from repeated nonequilibrium single-molecule pulling experiments. Indeed, the free energy values determined by the Jarzynski equality are lesser than the normal work average at all extensions. The free energy profiles constructed for the two velocities (100 and 400 nm/s) overlap with each other. The unfolding free energy barrier is estimated to be ∼7.5 kcal/mol. We anticipate that the Jarzynski equality can be applied in a similar manner to other ubiquitin-like proteins to extract their differences in the free energy profile, and hence, the effect of sequence diversity of structurally homologous proteins on the free energy landscape can be studied.

Recent Advances and Emerging Challenges in the Molecular Modeling of Mechanobiological Processes

Many biological processes result from the effect of mechanical forces on macromolecular structures and on their interactions. In particular, the cell shape, motion, and differentiation directly depend on mechanical stimuli from the extracellular matrix or from neighboring cells. The development of experimental techniques that can measure and characterize the tiny forces acting at the cellular scale and down to the single-molecule, biomolecular level has enabled access to unprecedented details about the involved mechanisms. However, because the experimental observables often do not provide a direct atomistic picture of the corresponding phenomena, particle-based simulations performed at various scales are instrumental in complementing these experiments and in providing a molecular interpretation. Here, we will review the recent key achievements in the field, and we will highlight and discuss the many technical challenges these simulations are facing, as well as suggest future directions for improvement.