Molecular dynamic (MD) studies on Gln233Arg (rs1137101) polymorphism of leptin receptor gene and associated variations in the anthropometric and metabolic profiles of Saudi women

The Gln233Arg (A>G; rs1137101) polymorphism of the leptin receptor gene (LEPR) has been investigated extensively and is reported to be associated with different metabolic states. In this investigation, we aimed to study the frequency of Gln233Arg genotypes and alleles in a group of Saudi women stratified by their body mass index (BMI), to correlate the LEPR genotypes with variations in anthropometric, lipid and hormonal parameters and to investigate conformational and structural variations in the mutant LEPR using molecular dynamic (MD) investigations. The study group included 122 Saudi women (normal weight = 60; obese = 62) attending the clinics for a routine checkup. Anthropometric data: height, weight, waist and hip circumference were recorded and fasting serum sample was used to estimate glucose, lipids, ghrelin, leptin and insulin. BMI, W/H ratio, and HOMA-IR values were calculated. Whole blood sample was used to extract DNA; exon 6 of the LEPR gene was amplified by PCR and sequencing was conducted on an ABI 3100 Avant Genetic Analyser. Molecular Dynamic Simulation studies were carried out using different softwares. The results showed the presence of all three genotypes of Gln233Arg in Saudi women, but the frequencies were significantly different when compared to reports from some populations. No differences were seen in the genotype and allele frequencies between the normal weight and obese women. Stratification by the genotypes showed significantly higher BMI, waist and hip circumference, leptin, insulin, fasting glucose and HOMA-IR and lower ghrelin levels in obese women carrying the GG genotype. Even in the normal weight group, individuals with GG genotype had higher BMI, waist and hip circumference and significantly lower ghrelin levels. The MD studies showed a significant effect of the Gln/Arg substitution on the conformation, flexibility, root-mean-square fluctuation (RMSF), radius of gyration (Rg) values, solvent-accessible surface area (SASA) and number of inter- and intra-molecular H-bonds. The results suggest that the structural changes brought about by the mutation, influence the signaling pathways by some unknown mechanism, which may be contributing to the abnormalities seen in the individuals carrying the G allele of rs1137101.