Multi-dimensional structural footprint identification for the design of potential scaffolds targeting METTL3 in cancer treatment from natural compounds

Abstract

Context: [Formula: see text]-adenosine-methyltransferase (METTL3) is the catalytic domain of the 'writer' proteins which is involved in the post modifications of [Formula: see text]-methyladinosine ([Formula: see text]). Though its activities are essential in many biological processes, it has been implicated in several types of cancer. Thus, drug developers and researchers are relentlessly in search of small molecule inhibitors that can ameliorate the oncogenic activities of METTL3. Currently, STM2457 is a potent, highly selective inhibitor of METTL3 but is yet to be approved.

Methods: In this study, we employed structure-based virtual screening through consensus docking by using AutoDock Vina in PyRx interface and Glide virtual screening workflow of Schrodinger Glide. Thermodynamics via MM-PBSA calculations was further used to rank the compounds based on their total free binding energies. All atom molecular dynamics simulations were performed using AMBER 18 package. FF14SB force fields and Antechamber were used to parameterize the protein and compounds respectively. Post analysis of generated trajectories was analyzed with CPPTRAJ and PTRAJ modules incorporated in the AMBER package while Discovery studio and UCSF Chimera were used for visualization, and origin data tool used to plot all graphs.

Results: Three compounds with total free binding energies higher than STM2457 were selected for extended molecular dynamics simulations. The compounds, SANCDB0370, SANCDB0867, and SANCDB1033, exhibited stability and deeper penetration into the hydrophobic core of the protein. They engaged in relatively stronger intermolecular interactions involving hydrogen bonds with resultant increase in stability, reduced flexibility, and decrease in the surface area of the protein available for solvent interactions suggesting an induced folding of the catalytic domain. Furthermore, in silico pharmacokinetics and physicochemical analysis of the compounds revealed good properties suggesting these compounds could serve as promising MEETL3 entry inhibitors upon modifications and optimizations as presented by natural compounds. Further biochemical testing and experimentations would aid in the discovery of effective inhibitors against the berserk activities of METTL3.

Keywords: Consensus docking; MD simulations; METTL3 inhibitors; MM-PBSA; Natural compounds.