Impact of emerging mutations on the dynamic properties the SARS-CoV-2 main protease: anin silicoinvestigation

The new coronavirus (SARS-CoV-2) is a global threat to world health and its economy. Its main protease (M^pro), which functions as a dimer, cleaves viral precursor proteins in the process of viral maturation. It is a good candidate for drug development owing to its conservation and the absence of a human homolog. An improved understanding of the protein behaviour can accelerate the discovery of effective therapies in order to reduce mortality. 100 ns all-atom molecular dynamics simulations of 50 homology modelled mutant M^prodimers were performed at pH 7 from filtered sequences obtained from the GISAID database. Protease dynamics were analysed using RMSD, RMSF, R_g, the averagedbetweenness centralityand geometry calculations. Domains from each M^proprotomer were found to generally have independent motions, while the dimer-stabilising N-finger region was found to be flexible in most mutants. A mirrored interprotomer pocket was found to be correlated to the catalytic site using compaction dynamics, and can be a potential allosteric target. The high number of titratable amino acids of M^promay indicate an important role of pH on enzyme dynamics, as previously reported for SARS-CoV. Independent coarse-grained Monte Carlo simulations suggest a link between rigidity/mutability and enzymatic function.