Molecular docking studies of human MCT8 protein with soy isoflavones in Allan-Herndon-Dudley syndrome(AHDS)

Divya Shaji

Department of Pediatrics,Aichi Medical University,Nagakute,Aichi,Japan

Keywords:Molecular docking MCT8 AHDS Soy isoflavones Daidzein Genistein

A B S T R A C T Monocarboxylate transporter-8(MCT8)is a specific thyroid hormone transporter,essential for the uptake of thyroid hormone into target tissues.Mutations in the MCT8 gene have been identified as the cause of Allan-Herndon-Dudley syndrome(AHDS).It has been reported that soy isoflavones influence thyroid hormone system and can interact with thyroid hormone transporter proteins.Therefore,the present study aimed to find out whether soy isoflavones(genistein,daidzein and glycitein)can be used as a natural inhibitor to target MCT8 in AHDS.Docking studies were performed for soy isoflavones in order to evaluate their binding affinity to MCT8 protein using AutoDock4(version 4.2.6)and AutoDock Vina.After docking,the ligands were ranked according to their binding energy and the best lead compound was selected based on the least binding energy.The docking results indicated that daidzein possesses the lowest binding energy against MCT8.Moreover,it was found that the residues PRO-338,HIS-341,and GLU-348 were involved in hydrogen bond interactions with genistein and daidzein.This study suggests that daidzein is a promising natural inhibitor to target MCT8 in AHDS.

1.Introduction

Thyroid hormone plays an important role in our body with widespread biological actions.Adequate levels of thyroid hormone are crucial for the development and proper function of multiple organs[1,2].Thyroid hormone exists in two forms:T4(3,3’,5,5’-tetraiodothyronine)and T3(3,3’,5-triodothyronine).The biological activity of thyroid hormone is related to the level of T3 within the cell.MCT8(monocarboxylate transporter 8)protein is a transporter specific for T3[1].The MCT8 gene is located in Xq13 and mutations in MCT8 are responsible for an X-linked condition,known as the Allan–Herndon–Dudley syndrome(AHDS)showing high serum T3levels in affected male patients[1,2].

Currently,there is no effective treatment available for AHDS.Hence,there is an urgent need for the identification and validation of novel drug lead compounds for treating AHDS.Therefore,this study aimed to identify the natural potent inhibitors of MCT8 from soy isoflavones.While isoflavones occur in many types of legumes,soybean contains the highest concentration of isoflavones.Genistein,daidzein and glycitein are the soy isoflavones typically accounting for 50%,40%and 10%,respectively[3].Due to their chemical structure,the isoflavones can bind to estrogen receptors.As a result of this binding,isoflavones inhibit and promote the expression of estrogensensitive genes[4].Previous researches have shown that the occurrence of breast cancer is lower in Asian individuals than in other populations because of the high soy consumption as part of their regular diet[3,4].

Soy isoflavones influence thyroid hormone system and can interact with thyroid hormone transporter proteins[5,6].In vitro and in vivo studies have indicated that genistein and daidzein are the potent inhibitors of thyroid peroxidase[5].Moreover,soy isoflavones are the potent ligands for transthyretin in serum and cerebrospinal fluid[5].However,soy isoflavones have not been studied against MCT8.So this study made an attempt to find out whether soy isoflavones(genistein,daidzein and glycitein)can be used as a natural inhibitor to target MCT8 in AHDS.

2.Materials and methods

2.1.Drug-likeness of the ligands

Genistein,daidzein and glycitein were considered as ligands.SwissADME tool[7]was used to calculate the molecular properties of the ligands.The molecular properties were screened based on the “Lipinski's rule of five”[8,9].The total polar surface area(TPSA)and the number of rotatable bonds were also calculated using SwissADME[7].

Fig.1.Chemical structures of(A)genistein,(B)daidzein and(C)glycitein.

Table 1 Molecular properties of the ligands.

Fig.2.3D structure of MCT8 protein by multi-template homology modeling[11].This 3D model of MCT8 was predicted based on multiple templates(PDB IDs:1pw4A,4u4tA,4ikxA,4j05A,and 4gbyA)using the advanced modeling feature of MODELLERv9.17[12,13].

2.2.Protein preparation for docking

The human MCT8 protein is not available in protein data bank(PDB)[10].The protein used in the docking study was obtained through multi-template based homology modeling[11].In the previous study[11],a good-quality 3D model of MCT8 was predicted based on multiple templates(PDB IDs:1pw4A,4u4tA,4ikxA,4j05A,and 4gbyA)using the advanced modeling feature of MODELLERv9.17[12,13].The structure refinement of modeled protein was done by ModRefiner[14].For docking,all water molecules were removed and polar hydrogen atoms were added to the refined model using AutoDock Tools(ADT)[15].The prepared protein was saved in PDBQT format.

2.3.Ligand preparation for docking

The ligands were downloaded from Pubchem Database[16,17]and converted to PDB[10] file format by using Openbabel software[18].The ligands were prepared using ADT[15].Gasteiger charge was assigned to the ligands.The prepared ligands were saved in PDBQT format.

2.4.Molecular docking

AutoDock4(version 4.2.6)[15]and AutoDock Vina[19]were used for molecular docking studies.AutoGrid program supplied with AutoDock4[15]was used for the preparation of grid maps.The grid box size was set at 76,70,and 76 Å for x,y,and z,respectively.The spacing between the grid points was 1.0 Å.The grid centre was set at 30.375,17.112,and-37.003 Å for x,y,and z,respectively.The Lamarckian Genetic Algorithm(LGA)was chosen to search for the best conformers.During the docking process,a maximum of 10 conformers was considered for each ligand.All the docking processes were performed with the default parameters of AutoDock 4[15].Population size was set to 150,maximum number of evaluations 2,500,000,maximum number of generations 27,000,maximum number of top individual that automatically survived 1,gene mutation rate 0.02 and crossover rate 0.8.Auto-Dock4[15]and AutoDock Vina[19]were compiled and run under Windows 10 Operating System.All figures with structure representations were produced using Discovery Studio Visualizer[20].

3.Results and discussion

3.1.Drug-likeness of the ligands

“Lipinski’s rule of five”[8,9]is used to evaluate the drug-likeness of a chemical compound.The molecular properties of a chemical compound consist of molecular weight,hydrogen bond donor,hydrogen bond acceptor,and logP. “Lipinski’s rule of five”states that an orally active drug has no more than one violation of the following criteria:(1)less than 5 hydrogen-bond donors,(2)less than 10 hydrogen-bond acceptors,(3)a molecular mass less than 500 Da,and(4)logPnot greater than 5.All ligands of the present study met the requirements of“Lipinski’s rule of five”.The other significant properties such as total polar surface area(TPSA)and the number of rotatable bonds were also calculated.TPSA of a compound should be less than 140 Å2and the number of rotatable bonds should be less than 10[21].All the ligands had the above properties.The chemical structures of the ligands are shown in Fig.1 and their molecular properties are shown in Table 1.

3.2.Molecular docking

Fig.3.(A)The output of AutoDock showing the binding site residues of MCT8 protein with the ligand genistein.The residues in the binding site are shown in red color.Genistein is shown in purple stick format.(B)2D diagram showing the types of contacts formed between MCT8 and genistein.The green dotted lines indicate H‑bond interactions between MCT8 and genistein.The values adjacent to the green dotted lines indicate their distance.

Fig.4.(A)The output of AutoDock showing the binding site residues of MCT8 protein with the ligand daidzein.The residues in the binding site are shown in red color.Daidzein is shown in green stick format.(B)2D diagram showing the types of contacts formed between MCT8 and daidzein.The green dotted lines indicate H‑bond interactions between MCT8 and daidzein.The values adjacent to the green dotted lines indicate their distance.

Fig.5.(A)The output of AutoDock showing the binding site residues of MCT8 protein with the ligand glycitein.The residues in the binding site are shown in blue color.Glycitein is shown in red stick format.(B)2D diagram showing the types of contacts formed between MCT8 and glycitein.

Fig.6.(A)The output of AutoDock Vina showing the binding site residues of MCT8 protein with the ligand genistein.The residues in the binding site are shown in green color.Genistein is shown in red stick format.(B)2D diagram showing the types of contacts formed between MCT8 and genistein.The green dotted lines indicate H‑bond interactions between MCT8 and genistein.The values adjacent to the green dotted lines indicate their distance.

Fig.7.(A)The output of AutoDock Vina showing the binding site residues of MCT8 protein with the ligand daidzein.The residues in the binding site are shown in green color.Daidzein is shown in green stick format.(B)2D diagram showing the types of contacts formed between MCT8 and daidzein.

Molecular docking is an important tool in pharmaceutical research[22].The molecular docking approach can be used to model the interaction between a ligand and a protein at the atomic level[22,23].The docking process involves two basic steps:prediction of the ligand conformation and assessment of the binding affinity[22].These two steps are related to sampling methods and scoring schemes,respectively[22].

Docking studies were performed for 3 compounds in order to evaluate their binding affinity to MCT8 protein using AutoDock4(version 4.2.6)[15]and AutoDock Vina[19].The 3D structure of the MCT8 protein is shown in Fig.2.The results were analyzed based on the binding energies of the docked complexes.Auto-Dock4[15]and AutoDock Vina[19]generated 10 poses for each ligand.The selection of the best pose was done on the least binding energy between the ligand and the protein.After docking,the ligands were ranked according to their binding energy.

3.2.1.Docked results with AutoDock4

All the ligand molecules were docked against MCT8 using AutoDock4[15].The best selected pose of MCT8-genistein docked complex(binding energy-5.76 kcal/mol)with binding site residues is shown in Fig.3A.The hydrogen bonds and the types of contacts involved in MCT8-genistein complex are shown in Fig.3B.It was observed that PRO-338,HIS-341,LEU-342,MET-343 and GLU-348 were involved in hydrogen bond interactions.

The best selected pose of MCT8-daidzein docked complex(binding energy-6.22kcal/mol)with binding site residues is shown in Fig.4A.The hydrogen bonds and the types of contacts involved in MCT8-daidzein complex are shown in Fig.4B.It was observed that PRO-338,HIS-341,and GLU-348 were involved in hydrogen bond interactions.

The best selected pose of MCT8-glycitein docked complex(binding energy-5.54kcal/mol)with binding site residues is shown in Fig.5A.The types of contacts involved in MCT8-glycitein complex are shown in Fig.5B.The results showed that there were no hydrogen bonds formed between MCT8 and glycitein.

The docking with AutoDock4 showed that daidzein was the best scored compound against MCT8 with the lowest binding energy.

Fig.8.(A)The output of AutoDock Vina showing the binding site residues of MCT8 protein with the ligand glycitein.The residues in the binding site are shown in yellow color.Glycitein is shown in yellow stick format.(B)2D diagram showing the types of contacts formed between MCT8 and glycitein.The green dotted lines indicate H‑bond interactions between MCT8 and glycitein.The values adjacent to the green dotted lines indicate their distance.

Table 2 Binding energies of the ligands.

3.2.2.Docked results with AutoDock Vina

The ligands were docked against MCT8 using AutoDock Vina[19].The best scored conformation was selected by considering the lowest binding energy between the protein and the ligand.The best selected pose of MCT8-genistein docked complex(binding energy-8.6 kcal/mol)predicted by AutoDock Vina is shown in Fig.6A.The 2D diagram showing the hydrogen bonds and the types of contacts involved in MCT8-genistein complex is shown in Fig.6B.The results showed that genistein interacted with MCT8 by forming hydrogen bonds with LYS-133.

The best docking pose of MCT8-daidzein docked complex(binding energy-8.6 kcal/mol)is shown in Fig.7A.The types of contacts involved in MCT8-daidzein complex is shown in Fig.7B.There were no hydrogen bonds formed between MCT8 and daidzein.

The best docking pose of MCT8-glycitein complex(binding energy-8.5kcal/mol)is shown in Fig.8A.The 2D diagram showing the hydrogen bonds and the types of contacts involved in MCT8-glycitein complex is shown in Fig.8B.It was observed that CYS-110 and ASN-111 were involved in hydrogen bond interactions.

3.2.3.Analysis of the docked results

The docking results predicted by AutoDock4[15]were compared to those of AutoDock Vina[19].Docking analysis showed that there were hydrogen bonds formed between MCT8 and the inhibitors used.Opposite to the results of AutoDock4[15],Auto-Dock Vina[19]generated hydrogen bonds in MCT8-glycitein interaction.Opposite to the results of AutoDock Vina,AutoDock generated hydrogen bonds in MCT8-daidzein interaction.Moreover,van der Waals interactions were also involved in addition to hydrogen bonds.It is important to point out that both compounds,genistein and daidzein,have hydrogen bonds with residues PRO-338,HIS-341 and GLU-348.

The best lead compound was selected in terms of binding energy.The binding energies of the ligands calculated by AutoDock4[15]and AutoDock Vina[19]are shown in Table 2.Based on the analysis with AutoDock Vina[19],it was observed that binding energies of the three compounds were almost the same.Docking studies with AutoDock4[15]and AutoDock Vina[19]showed that the natural compound Daidzein showed the lowest binding energy value of-6.22 and-8.6 kcal/mol,respectively(Table 2).Based on these findings,gaidzein can be used as a natural inhibitor to target MCT8 in AHDS.However,daidzein should be subjected to further investigation using in vitro studies.

4.Conclusion

Recently,many researches have focused on the identification of inhibitors from natural sources.This study concludes that daidzein is an effective lead compound which will be useful for the design of novel less toxic and highly efficient drugs for the treatment of AHDS.Daidzein should be subjected to further experimental study in order to confirm this finding.This study also identified that PRO-338,HIS-341 and GLU-348 of MCT8 play an important role in hydrogen bonding with genistein and daidzein.

Conflicts of interest

The authors declare that there are no conflicts of interest.