Comparing Multislice Projections of MD Simulations with CryoEM Exposes Structural Prediction Errors

Cryo-electron microscopy (cryoEM) is a powerful tool for atomic- and molecular-resolution structure determination, while molecular dynamics (MD) simulations are similarly powerful tools for predicting molecular trajectories. Given the challenges in estimating biomolecule dynamics with cryoEM alone, MD simulations are employed to forecast molecular motions and to interpret cryoEM reconstructions. Few methods, however, can evaluate MD predictions directly. Here, we use multislice wave propagation to project sampled snapshots of MD trajectories, either coarse-grained (CG) or all-atom (AA), into simulated cryoEM 3D reconstructions. We compared simulated and experimental images of low- and high-curvature membranes to show that MD simulations qualitatively reflect the fluidity and thus the contrast of biological membranes observed by cryoEM. MD simulations also correctly predicted bilayer dimensions for single component flat bilayers observed in cryoEM images. However, Martini3 CG-MD simulations failed to predict changes in membrane thickness induced by high curvature and with heterogeneous lipid compositions. We pinpointed the misbehavior of polyunsaturated lipid tails and cholesterol in Martini3 simulations as the main error sources contributing to inaccurate bilayer thicknesses. Our comparisons also explain membrane structure discrepancies between cryoEM and small angle X-ray scattering (SAXS). Further testing of MD predictions by direct comparisons between simulated and experimental cryoEM images should lead to the development of more accurate MD force fields.