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.