DNA Methyltransferases Directed Anti-Cancerous Plant Medicine ( Xanthomicrol and Galloyl ) Based Molecular Docking and Dynamics Simulation

DNA methyltransferases 1 (DNMT1) has been looked as crucial targets against various types of cancers. MD simulations have advanced to a point where the atomic level information of biological macromolecule (protein or DNA-protein or protein-protein) can easily be advantageous to predict the functionality. In this study we utilize xanthomicrol and galloyl compounds to investigate potential compounds for the inhibition of DNMT1, and the results of these two compounds are compared with drug decitabine. Xanthomicrol and galloyl are found to dock successfully within the active site of DNMT1. A comparison of the inhibitory potential of screened xanthomicrol inhibited DNMT1 approximately is identical with those of their corresponding drugs, decitabine. The stability of the DNMT1 with the best docked xanthomicrol, were further analysed in molecular dynamics (MD) simulation and compared with those of the respective drugs namely decitabine which revealed stabilization of these complexes within 300 ns of simulation with better stability of DNMT1.


Introduction
DNA methylation refers to the addition of a methyl group to one of the four bases that constitute the coding sequence of DNA.In humans, methylation is normally added only on the 5 position of the cytosine base in a post-DNA synthesis reaction catalyzed by one of several DNA methyltransferases.DNA methylation plays a key role in chromatin structure, suppression of the activity of endogenous parasitic sequences, and stable suppression of gene expression (epigenetic silencing), and a process normally is reserved for special situations such as the inactive X-chromosome and imprinted genes [1].Epigenetic aberrations are common in the development and progression of cancer cells and are mediated by DNA methyltransferases (DNMTs), histone deacetylases (HDACs) and microRNA (miRNA).These modifications alter gene function and malignant cellular transformation.DNA methylation is a reversible reaction, and catalyzed by three major DNMTs.One is DNMT1, which preserves the methylation patterns throughout each cell division [2] [3].The others are DNMT3a and DNMT3b, which transfer a methyl group to previously unmethylated genomic regions [4].Hypermethylation of DNA is a key epigenetic mechanism for the silencing of many genes, including those for cell cycle regulations, receptors, DNA repair, and apoptosis [5].In this aberrant methylation, the cytosine of the CpG Island in or near the promoter region of the newly synthesized DNA strand is methylated by DNMT.The methylated CpG Island has higher binding affinity to methyl-CpG binding domain proteins that recruit transcriptional corepressors, such as histone deacetylases, resulting in chromosome condensation and transcription repression [6] [7].The inhibition of DNMT, especially DNMT1, would block the hypermethylation of the newlysynthesized DNA strand, and results in the reversal of the hypermethylation and there-expression of the silenced genes [8]- [10].The therapeutic activity of three cytosine nucleoside analogues, 5-azacytidine (azacitidine, 5 AC) and 5-aza-20-deoxycytidine (Decitabine, DAC) in acute myeloid leukemia (AML) and myelodysplastic syndromes (MDS) has been the focus of numerous recent studies.When applied at low doses, the compounds reduce genomic DNA methylation as a consequence of their irreversible binding to DNMT after incorporation into newly synthesized DNA [11].Interference with aberrant DNA methylation, leading to the reactivation of silenced tumor suppressor genes, accounts, at least in part, for the DNMT inhibitors' antitumoral effects [12].Agents that preferentially inhibit proliferation of neoplastic cells may offer means of developing drugs effectively against cancer in human.One of the strategies for finding such agents is screening plant material, especially those which have been reported to have anti-cancer, anti-inflammatory, anti-fungal or anti-bacterial effects.Spinal-Z is a traditional Iranian anticancer remedy.Spinal-Z is an extract consisting of two plants: Peganum harmala Linn.Seeds and Dracocephalum kotschyii Boiss leaves.The active ingredient in the leaf extract of D. kotschyii appears to be a flavone identified as xanthomicrol.Xanthomicrol is able to inhibit proliferation of a number of malignant cells [13].Four kinds of galloyl flavonol glycosides are found in the leaf extract of Pemphis acidula, a plant growing on the subtropical seashore.The galloyl flavonoids show more efficient antioxidant activity than that of the corresponding flavonol glycosides [14].Docking is the process by which two molecules fit together in three-dimensional space and molecular docking is a useful tool in structural molecular biology and computer-assisted drug design, and bioinformatics is seen as an emerging field with the potential to significantly improve how drugs are found and brought to the clinical trials and eventually released to the marketplace.This method is widely used to predict the predominant binding mode(s) oligand with a protein.Besides, reliability in docking of ligand molecules to protein or other targets is an important challenge for molecular modeling.One applications of the docking technique includes the prediction of the binding mode of novel drugs.Docking permits the scientist to virtually screen a database of compounds and bode the strongest binders based on variant scoring functions.The docking analyses are carried by AutoDock (4.2) docking software.In this research, we utilize xanthomicrol and galloyl compounds to investigate potential lead drugs for DNMT1.The database virtual screen provides a rapid search to identify potent compounds for DNMT1 interactions.Molecular dynamics (MD) simulation is used to simulate the dynamic changes between DNMT1 protein and docked ligand.
In this study we utilized computer-aided drug design (CADD) in this research in order to design potential lead drugs for AD therapy.CADD is an efficient approach for the rapid identification of potential lead compounds in target therapy and has been widely used in drug design including virus therapy, cancer therapy, treatment of sleeplessness, neuropathic therapy, weight loss therapy, diabetic treatment, inflammation treatment, erectile dysfunction treatment, and hair loss therapy [15].

Molecular Docking
At first 3D atomic coordinates of DNMT1 (ID: 2QRV) were obtained from the Protein Data bank [16] and prepared for docking.AutoDock (4.2), a similar suite of programs involved genetic algorithm, was employed to gain an insight into the compounds (Figure 1) bindings with DNMT1.All nonpolar hydrogens were merged and partial atomic charges were assigned using the Gasteiger-Marsili method.Different grid boxes with different grid points in size with a grid-point spacing of 0.375 Å were considered for docking.Each map was centered such that it covered the entire protein including all possible binding sites.

Docking Analysis
The generated ligands were docked into the defined binding site on the DNMT1protein structure.Ligplot plus was used to analysis docking poses for H-bond and hydrophobic interactions.

Molecular Dynamics Simulation
Then molecular dynamics simulations were performed using the Gromacs 4.5.4software, under G43a1 force field.Topology files and parameters of small compounds in protein-ligand complexes were generated for GROMACS simulation by SwissParam web server.The calculations were carried out using the crystallographic structure of domain DNMT1.The complex with the best position obtained from docking, was embedded into a box containing SPC216 model water.At first the entire system was minimized using the steepest descent followed by conjugate gradients algorithms, then the solvent and Na + , CL − were allowed to evolve using minimization and molecular dynamics at NVT ensemble for 500 ps and at NPT ensemble for 1000 ps at 100 K in which initial configuration of the structures was kept fixed.Finally in order to obtain equilibrium geometry at 300 K and 1 atm, 1 ns MD simulation run was performed at NPT ensemble.The time step for the simulation was 2 fs.Van der Waals and electrostatic forces cutoff were 14 Å and 12 Å respectively.

Docking Analysis
We used the Autodock

Molecular Dynamics Simulation Analyses
The stability of the DNMT1 during the simulation in the presence of Decitabine, Xanthomicrol and Galloyl was evaluated by the root mean square deviation (RMSD) and by monitoring its secondary structure elements during the simulations.As documented in Figure 3, the backbone RMSD was steady over time, exhibiting fluctuations that reached maximal values of 4.0 Å in the simulation.This consistency dictates that the DNMT1was stable under the simulation conditions.For total energy analysis, no significantly increased values were observed among all simulation times (Figure 4).The total energy of all systems remained in −9695000 kJ/mol.These results suggest that all structures of the complexes tend to become constant after the initial simulation time.In the present study, it may be concluded that that the Xanthomicrol qualifies for its further testing, using in-vivo/in-vitro studies, to develop as a potential DNMT1 directed anti-cancerous drug.
As DNMT1 have one hydrogen bond donor (OH) and one hydrogen bond acceptor (O), the probability to

Conclusion
DNMT1 plays the key role in cell proliferation phase of all types of developing cancer.Conceptually, targeting DNMT1 will disrupt both DNA replication and transcription, thereby leading to inhibition of cell division and consequently stopping the growth of actively dividing cancerous cells.Therefore, any therapeutic approach targeted against DNA replication in cancerous cells will certainly be a key to the development of a novel anti-cancerous agent.Furthermore, long term use of synthetic drugs such as Decitabine has been reported to develop severe side effects.Therefore, for therapeutic applications, there is a need to identify more specific inhibitors of DNMT1, coming from natural sources, which are relatively safer in use with little or no side effects.In view of these, in the present study, plant derived one hundred anti-cancerous has been selected and their pharmacokinetic properties are predicted through insilico analyses.Thus Xanthomicrol and Galloyl are selected.Furthermore, docking analyses reveal that Xanthomicrol inhibites DNMT1 more strongly, while Galloyl inhibites preferentially DNMT1 weaker than those of respective drugs, Decitabine outstanding.The results of best docked, Xanthomicrol, are further validated and compared with those of respective drugs namely Decitabine using molecular dynamics (MD) simulation.Thus, the present study makes a foundation for further investigations based on the experimental data (wet lab data) for therapeutic application of screened Plant Medicine in general and Xanthomicrol in particular.