Data-guided Multi-Map variables for ensemble refinement of molecular movies

Abstract

AbstractDriving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. Here, the 3-dimensional electron density of a protein, as it would be determined by cryo-EM or X-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work uses the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along reaction coordinates derived from electron densities can serve to induce conformational transitions between known intermediates. The simulated pathways appear reversible, with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein DataBank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution (~5 Å) density information coupled to comparatively large structural transitions. Practical considerations for generating transitions with multiple intermediate atomic density distributions are also discussed.