Rapid Sampling of Folding Equilibria of β-Peptides in Methanol Using a Supramolecular Solvent Model.

Molecular dynamics simulation of biomolecules in solvent using an atomic model for both the biomolecules and the solvent molecules is still computationally rather demanding considering the time scale of the biomolecular motions. The use of a supramolecular coarse-grained (CG) model can speed up the simulation considerably, but it also reduces the accuracy inevitably. Combining an atomic fine-grained (FG) level of modeling for the biomolecules and a supramolecular CG level for the solvent into a hybrid system, the increased computational efficiency may outweigh the loss of accuracy with respect to the biomolecular properties in the hybrid FG/CG simulation. Here, a previously published CG methanol model is reparametrized, and then a 1:1 mixture of FG and CG methanol is used to calibrate the FG-CG interactions using thermodynamic and dielectric screening data for liquid methanol. The FG-CG interaction parameter set is applied in hybrid FG/CG solute/solvent simulations of the folding equilibria of three β-peptides that adopt different folds. The properties of the peptides are compared with those obtained in FG solvent simulations and with experimental NMR data. The comparison shows that the folding equilibria in the pure CG solvent simulations are different from those in the FG solvent simulations because of the lack of hydrogen-bonding partners in the supramolecular CG solvent. Next, we introduced an FG methanol layer around the peptides in CG solvent to recover the hydrogen-bonding pattern of the FG solvent simulations. The result shows that with the FG methanol layer, the folding equilibria of the three β-peptides are very similar to those in the FG solvent simulations, while the computational efficiency is at least 3 times higher and the cutoff radius for nonbonded interactions could be increased from 1.4 to 2.0 nm.