Forced Unfolding Mechanism of Bacteriorhodopsin as Revealed by Coarse-Grained Molecular Dynamics

Mechanistic Insight into the Mechanical Unfolding of the Integral Membrane Diacylglycerol Kinase

The essential forces stabilizing membrane proteins and governing their folding and unfolding are difficult to decipher. Single-molecule atomic force spectroscopy mechanically unfolds individual membrane proteins and quantifies their dynamics and energetics. However, it remains challenging to structurally assign unfolding intermediates precisely and to deduce dominant interactions between specific residues that facilitate either the localized stabilization of these intermediates or the global assembly of membrane proteins. Here, we performed force spectroscopy experiments and multiscale molecular dynamics simulations to study the unfolding pathway of diacylglycerol kinase (DGK), a small trimeric multispan transmembrane enzyme. The remarkable agreement between experiments and simulations allowed precise structural assignment and interaction analysis of unfolding intermediates, bypassing existing limitations on structural mapping, and thus provided mechanistic explanations for the formation of these states. DGK unfolding was found to proceed with structural segments varying in size that do not correlate with its secondary structure. We identified intermolecular side-chain packing interactions as one of the major contributions to the stability of unfolding intermediates. Mutagenesis creating packing defects induced a dramatic decrease in the mechano-stability of corresponding intermediates and also in the thermo-stability of DGK trimer, in good agreement with predictions from simulations. Hence, the molecular determinants of the mechano- and thermo-stability of a membrane protein can be identified at residue resolution. The accurate structural assignment established and microscopic mechanism revealed in this work may substantially expand the scope of single-molecule studies of membrane proteins.

Single-Molecule Force Spectroscopy of Membrane Protein Folding

Single-molecule force spectroscopy is a unique method that can probe the structural changes of single proteins at a high spatiotemporal resolution while mechanically manipulating them over a wide force range. Here, we review the current understanding of membrane protein folding learned by using the force spectroscopy approach. Membrane protein folding in lipid bilayers is one of the most complex biological processes in which diverse lipid molecules and chaperone proteins are intricately involved. The approach of single protein forced unfolding in lipid bilayers has produced important findings and insights into membrane protein folding. This review provides an overview of the forced unfolding approach, including recent achievements and technical advances. Progress in the methods can reveal more interesting cases of membrane protein folding and clarify general mechanisms and principles.

Factors That Control the Force Needed to Unfold a Membrane Protein in Silico Depend on the Mode of Denaturation

Single-molecule force spectroscopy methods, such as AFM and magnetic tweezers, have proved extremely beneficial in elucidating folding pathways for soluble and membrane proteins. To identify factors that determine the force rupture levels in force-induced membrane protein unfolding, we applied our near-atomic-level Upside molecular dynamics package to study the vertical and lateral pulling of bacteriorhodopsin (bR) and GlpG, respectively. With our algorithm, we were able to selectively alter the magnitudes of individual interaction terms and identify that, for vertical pulling, hydrogen bond strength had the strongest effect, whereas other non-bonded protein and membrane-protein interactions had only moderate influences, except for the extraction of the last helix where the membrane-protein interactions had a stronger influence. The up-down topology of the transmembrane helices caused helices to be pulled out as pairs. The rate-limiting rupture event often was the loss of H-bonds and the ejection of the first helix, which then propagated tension to the second helix, which rapidly exited the bilayer. The pulling of the charged linkers across the membrane had minimal influence, as did changing the bilayer thickness. For the lateral pulling of GlpG, the rate-limiting rupture corresponded to the separation of the helices within the membrane, with the H-bonds generally being broken only afterward. Beyond providing a detailed picture of the rupture events, our study emphasizes that the pulling mode greatly affects the factors that determine the forces needed to unfold a membrane protein.

Unfolding and identification of membrane proteins in situ

Single-molecule force spectroscopy (SMFS) uses the cantilever tip of an AFM to apply a force able to unfold a single protein. The obtained force-distance curve encodes the unfolding pathway, and from its analysis it is possible to characterize the folded domains. SMFS has been mostly used to study the unfolding of purified proteins, in solution or reconstituted in a lipid bilayer. Here, we describe a pipeline for analyzing membrane proteins based on SMFS, that involves the isolation of the plasma membrane of single cells and the harvesting of force-distance curves directly from it. We characterized and identified the embedded membrane proteins combining, within a Bayesian framework, the information of the shape of the obtained curves, with the information from Mass Spectrometry and proteomic databases. The pipeline was tested with purified/reconstituted proteins and applied to five cell types where we classified the unfolding of their most abundant membrane proteins. We validated our pipeline by overexpressing 4 constructs, and this allowed us to gather structural insights of the identified proteins, revealing variable elements in the loop regions. Our results set the basis for the investigation of the unfolding of membrane proteins in situ, and for performing proteomics from a membrane fragment.

Normal mode analysis and beyond

Normal mode analysis provides a powerful tool in biophysical computations. Particularly, we shed light on its application to protein properties because they directly lead to biological functions. As a result of normal mode analysis, the protein motion is represented as a linear combination of mutually independent normal mode vectors. It has been widely accepted that the large amplitude motions throughout the entire protein molecule can be well described with a few low-frequency normal modes. Furthermore, it is possible to represent the effect of external perturbations, e.g., ligand binding, hydrostatic pressure, as the shifts of normal mode variables. Making use of this advantage, we are able to explore mechanical properties of proteins such as Young's modulus and compressibility. Within thermally fluctuating protein molecules under physiological conditions, tightly packed amino acid residues interact with each other through heat and energy exchanges. Since the structure and dynamics of protein molecules are highly anisotropic, the flow of energy and heat should also be anisotropic. Based on the harmonic approximation of the heat current operator, it is possible to analyze the communication map of a protein molecule. By using this method, the energy transfer pathways of photoactive yellow protein were calculated. It turned out that these pathways are similar to those obtained via the Green-Kubo formalism with equilibrium molecular dynamics simulations, indicating that normal mode analysis captures the intrinsic nature of the transport properties of proteins.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

On the Interpretation of Force-Induced Unfolding Studies of Membrane Proteins Using Fast Simulations

Single-molecule force spectroscopy has proven extremely beneficial in elucidating folding pathways for membrane proteins. Here, we simulate these measurements, conducting hundreds of unfolding trajectories using our fast Upside algorithm for slow enough speeds to reproduce key experimental features that may be missed using all-atom methods. The speed also enables us to determine the logarithmic dependence of pulling velocities on the rupture levels to better compare to experimental values. For simulations of atomic force microscope measurements in which force is applied vertically to the C-terminus of bacteriorhodopsin, we reproduce the major experimental features including even the back-and-forth unfolding of single helical turns. When pulling laterally on GlpG to mimic the experiment, we observe quite different behavior depending on the stiffness of the spring. With a soft spring, as used in the experimental studies with magnetic tweezers, the force remains nearly constant after the initial unfolding event, and a few pathways and a high degree of cooperativity are observed in both the experiment and simulation. With a stiff spring, however, the force drops to near zero after each major unfolding event, and numerous intermediates are observed along a wide variety of pathways. Hence, the mode of force application significantly alters the perception of the folding landscape, including the number of intermediates and the degree of folding cooperativity, important issues that should be considered when designing experiments and interpreting unfolding data.

Improved free-energy landscape reconstruction of bacteriorhodopsin highlights local variations in unfolding energy

Precisely quantifying the energetics that drive the folding of membrane proteins into a lipid bilayer remains challenging. More than 15 years ago, atomic force microscopy (AFM) emerged as a powerful tool to mechanically extract individual membrane proteins from a lipid bilayer. Concurrently, fluctuation theorems, such as the Jarzynski equality, were applied to deduce equilibrium free energies (ΔG0) from non-equilibrium single-molecule force spectroscopy records. The combination of these two advances in single-molecule studies deduced the free-energy of the model membrane protein bacteriorhodopsin in its native lipid bilayer. To elucidate this free-energy landscape at a higher resolution, we applied two recent developments. First, as an input to the reconstruction, we used force-extension curves acquired with a 100-fold higher time resolution and 10-fold higher force precision than traditional AFM studies of membrane proteins. Next, by using an inverse Weierstrass transform and the Jarzynski equality, we removed the free energy associated with the force probe and determined the molecular free-energy landscape of the molecule under study, bacteriorhodopsin. The resulting landscape yielded an average unfolding free energy per amino acid (aa) of 1.0 ± 0.1 kcal/mol, in agreement with past single-molecule studies. Moreover, on a smaller spatial scale, this high-resolution landscape also agreed with an equilibrium measurement of a particular three-aa transition in bacteriorhodopsin that yielded 2.7 kcal/mol/aa, an unexpectedly high value. Hence, while average unfolding ΔG0 per aa is a useful metric, the derived high-resolution landscape details significant local variation from the mean. More generally, we demonstrated that, as anticipated, the inverse Weierstrass transform is an efficient means to reconstruct free-energy landscapes from AFM data.