Multi-ligand molecular docking, simulation, free energy calculations and wavelet analysis of the synergistic effects between natural compounds baicalein and cubebin for the inhibition of the main protease of SARS-CoV-2.

[Display omitted] • Baicalein (HQA004) has the most favorable inhibitory ligand for Mpro of SARS-CoV-2. • Cubebin (CHA008) acted to support the proposed HQA004 inhibitor at the active site. • CHA008 served to enhance HQA004 effects via allosteric effects at two sites. • CHA008 enhanced the stability of HQA004 despite diffusion to a distant site. • Wavelet analysis showed that because of CHA008 binding, the dynamics and structure of Mpro have more subtle changes. Combination drugs have been used for several diseases for many years since they produce better therapeutic effects. However, it is still a challenge to discover candidates to form a combination drug. This study aimed to investigate whether using a comprehensive in silico approach to identify novel combination drugs from a Chinese herbal formula is an appropriate and creative strategy. We, therefore, used Toujie Quwen Granules for the main protease (Mpro) of SARS-CoV-2 as an example. We first used molecular docking to identify molecular components of the formula which may inhibit Mpro. Baicalein (HQA004) is the most favorable inhibitory ligand. We also identified a ligand from the other component, cubebin (CHA008), which may act to support the proposed HQA004 inhibitor. Molecular dynamics simulations were then performed to further elucidate the possible mechanism of inhibition by HQA004 and synergistic bioactivity conferred by CHA008. HQA004 bound strongly at the active site and that CHA008 enhanced the contacts between HQA004 and Mpro. However, CHA008 also dynamically interacted at multiple sites, and continued to enhance the stability of HQA004 despite diffusion to a distant site. We proposed that HQA004 acted as a possible inhibitor, and CHA008 served to enhance its effects via allosteric effects at two sites. Additionally, our novel wavelet analysis showed that as a result of CHA008 binding, the dynamics and structure of Mpro were observed to have more subtle changes, demonstrating that the inter-residue contacts within Mpro were disrupted by the synergistic ligand. This work highlighted the molecular mechanism of synergistic effects between different herbs as a result of allosteric crosstalk between two ligands at a protein target, as well as revealed that using the multi-ligand molecular docking, simulation, free energy calculations and wavelet analysis to discover novel combination drugs from a Chinese herbal remedy is an innovative pathway. [ABSTRACT FROM AUTHOR]