Natarajan, C.; Srivastava, A.

Accurately sampling the membrane-bound conformations of a peripheral membrane proteins (PMP) using classical all-atom molecular dynamics simulations (AAMD) is a formidable enterprise due to the rugged free energy landscape of the protein-membrane system. In general, AAMD-based extraction of binding geometry requires simulations of multiple systems with different initial user-defined binding poses that may not be exhaustive. As an alternative, advanced sampling methods are also applied to elucidate the membrane-binding mechanism of PMPs but these techniques are generally computationally expensive and often depend on choice of the collective variables (CV). In this work, we showcase the utility of CV-free replica exchange with hybrid tempering (REHT) method in capturing the membrane-bound conformations of PMPs by testing it on the Osh4 amphipathic lipid-packing sensor (ALPS) motif, a 27 amino-acid membrane binding peptide. We show that REHT samples all the membrane-bound conformations of the Osh4 ALPS peptide at their correct populations and does it in a highly efficient manner with minimum computational time. We clearly show that, out of the two significant conformations, the peptide prefers horizontal conformations over vertical ones. In both the conformations, REHT captures all the vital residue-wise membrane contacts. The transition between the two configuration is not uncommon as our calculations reveal a ~ 2 kT free energy difference between the two conformations. Interestingly, from our simulations we also find that the transition from vertical to horizontal conformation involves limited unfolding the main helix\'s last turn. From our findings, we conclude that REHT samples the membrane-bound conformations of Osh4 ALPS peptide very efficiently and also provides additional insights and information that are often not available with regular piece-wise AAMD simulations. The method can be used as an efficient tool to explore the membrane-binding mechanisms of PMPs.