Tutorial on quantifying and sampling biomolecular ensembles with ShapeGMM

Here, we present a detailed workflow for clustering and enhanced sampling of biomolecular conformations using the ShapeGMM methodology. This approach fits a probabilistic model of biomolecular conformations rooted in the idea that the free energy can be expressed in terms of local fluctuations in atomic positions around metastable states. We demonstrate using a single model system how to generate and fit equilibrium molecular dynamics simulation data. We then demonstrate how to use the resulting model to generate a reaction coordinate between two states, how to sample along that coordinate using metadynamics using our size-and-shape PLUMED module, and how to cluster those biased conformations to obtain a refined equilibrium ShapeGMM model.