Protein–protein
interactions between transmembrane helices
are essential elements for membrane protein structures and functions.
To understand the effects of peptide sequences and lipid compositions
on these interactions, single-molecule experiments using model systems
comprising artificial peptides and membranes have been extensively
performed. However, their dynamic behavior at the atomic level remains
largely unclear. In this study, we applied the all-atom molecular
dynamics (MD) method to simulate the interactions of single-transmembrane
helical peptide dimers in membrane environments, which has previously
been analyzed by single-molecule experiments. The simulations were
performed with two peptides (Ala- and Leu-based artificially designed
peptides, termed “host peptide”, and the host peptide
added with the GXXXG motif, termed “GXXXG peptide”),
two membranes (pure-POPC and POPC mixed with 30% cholesterols), and
two dimer directions (parallel and antiparallel), consistent with
those in the previous experiment. As a result, the MD simulations
with parallel dimers reproduced the experimentally observed tendency
that introducing cholesterols weakened the interactions in the GXXXG
dimer and facilitated those in the host dimer. Our simulation suggested
that the host dimer formed hydrogen bonds but the GXXXG dimer did
not. However, some discrepancies were also observed between the experiments
and simulations. Limitations in the space and time scales of simulations
restrict the large-scale undulation and peristaltic motions of the
membranes, resulting in differences in lateral pressure profiles.
This effect could cause a discrepancy in the rotation angles of helices
against the membrane normal.