Forced Unfolding Mechanism of Bacteriorhodopsin as Revealed by Coarse-Grained Molecular Dynamics.
Biophys J 2017;
111:2086-2098. [PMID:
27851934 DOI:
10.1016/j.bpj.2016.09.051]
[Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2015] [Revised: 09/23/2016] [Accepted: 09/26/2016] [Indexed: 12/18/2022] Open
Abstract
Developments in atomic force microscopy have opened up a new path toward single-molecular phenomena; in particular, during the process of pulling a membrane protein out of a lipid bilayer. However, the characteristic features of the force-distance (F-D) curve of a bacteriorhodopsin in purple membrane, for instance, have not yet been fully elucidated in terms of physicochemical principles. To address the issue, we performed a computer simulation of bacteriorhodopsin with, to our knowledge, a novel coarse-grained (C-G) model. Peptide planes are represented as rigid spheres, while the surrounding environment consisting of water solvents and lipid bilayers is represented as an implicit continuum. Force-field parameters were determined on the basis of auxiliary simulations and experimental values of transfer free energy of each amino acid from water to membrane. According to Popot's two-stage model, we separated molecular interactions involving membrane proteins into two parts: I) affinity of each amino acid to the membrane and intrahelical hydrogen bonding between main chain peptide bonds; and II) interhelix interactions. Then, only part I was incorporated into the C-G model because we assumed that the part plays a dominant role in the forced unfolding process. As a result, the C-G simulation has successfully reproduced the key features, including peak positions, of the experimental F-D curves in the literature, indicating that the peak positions are essentially determined by the residue-lipid and intrahelix interactions. Furthermore, we investigated the relationships between the energy barrier formation on the forced unfolding pathways and the force peaks of the F-D curves.
Collapse