Characterization, docking and molecular dynamics simulation of gonadotropin inhibitory hormone receptor (GnIHR2) in Labeo catla

Abstract

GnIH receptors (GnIHRs) belong to the family of G-protein coupled receptors (GPCRs) and play a key role in the regulation of reproduction from fishes to mammals, either by inhibiting or stimulating the expression of gonadotropins. The aim of this study was to characterize GnIH receptor (GnIHR2) from Indian Major Carp, Labeo catla and its docking and simulation with GnIH antagonist RF313. Methods: The full length sequence of GnIHR2 was obtained with RACE PCR. The docking analysis of RF313 with GnIHR2 receptor was performed with AutoDock v. 4.2.6 and molecular dynamics (MD) simulation with GROMACS 5.0. Results: In the present study, we cloned full-length cDNA (1733 bp) of GnIHR2 from the brain of L. catla. The phylogenetic analysis showed clustering of catla GnIHR2 with goldfish and zebrafish in the GPR147 group. L. catla GnIHR2 receptor comprised seven transmembrane domains and the 3D-structure was predicted by I-TASSER tool. The docking analysis revealed high binding affinity (-11.6 kcal/mol) of GnIH antagonist, RF313 towards GnIHR2 receptor. The primary bonds involved were alkyl and hydrogen bonds while the amino acids participated were proline 43, 210, 339, cysteine 214, leucine 211, serine 213 and phenylalanine 338. The MD simulation analysis of docked complex for 100 nano-seconds (ns) in the lipid membrane environment showed the stability of the complex with time. Conclusion: Our study showed that GnIH antagonist, RF313 interact tightly with the GnIH receptor, GnIHR2 of L. catla. To the best of our knowledge, this is the first report on computational modelling and MD simulation of GnIH receptor in fishes. This will help in functional characterization studies of GnIH/GnIHR system in vertebrates.