While successfully reproducing hydrophobic and hydrophilic interactions the Martini model is insufficient to keep a protein folded as it lacks electrostatic interactions. Using split charge equilibration at each time step can yield realistic dynamic bead charges. Combining this with a Drude oscillator based polarization model for all beads will permit modeling of hydrogen bonds to maintain secondary protein structure and will enable more accurate coarse-grained simulations of protein-protein and protein-membrane interactions.


Talk Number PIRSA:13120031
Speaker Profile Yuriy Khalak