Improved sampling for simulations of interfacial membrane proteins: application of generalized shadow hybrid Monte Carlo to a peptide toxin/bilayer system

The computational costs associated with performing molecular dynamics (MD) simulations are still somewhat prohibitive and therefore limit the time and length scales that can be currently achieved. One approach to overcoming the limited size and duration of a simulation is to reduce the amount of detail when representing a system of interest, generally termed "coarse-graining". An alternative approach is via more efficient sampling methods that offer an enhanced search of a complex multidimensional energy landscape. One could also combine enhanced sampling methods with a coarse-grained (CG) force field. Here, we apply generalized shadow hybrid Monte Carlo (GSHMC), a recently proposed simulation protocol, to a biomolecular system of moderate size and show that GSHMC offers improved sampling compared to standard MD simulation. Our test system is a CG representation of a small peptide toxin interacting with a phospholipid bilayer. Specifically, we show that GSHMC allows for a quicker localization of the toxin to its equilibrium location of interaction at the headgroup/water interface of the bilayer. GSHMC therefore potentially allows for future exploration of larger and more complex systems over longer periods, which would otherwise be impractical to perform using conventional simulation methodology.