Computational Analysis of Gly482Ser Single-Nucleotide Polymorphism in PPARGC1A Gene Associated with CAD, NAFLD, T2DM, Obesity, Hypertension, and Metabolic Diseases

The rs17782313 polymorphism near MC4R gene confers a high risk of obesity and hyperglycemia, while PGC1α rs8192678 polymorphism is weakly correlated with glucometabolic disorder: a systematic review and meta-analysis

Background

The relationships of the rs17782313 polymorphism near melanocortin 4 receptor gene (MC4R) and the rs8192678 polymorphism in peroxisome proliferator-activated receptor gamma coactivator 1 alpha gene (PGC1α) with metabolic abnormalities have been explored in many populations around the world, but the findings were not all consistent and sometimes even a bit contradictory.

Methods

Electronic databases including Medline, Scopus, Embase, Web of Science, CNKI and Google Scholar were checked for studies that met the inclusion criteria. Data were carefully extracted from eligible studies. Standardized mean differences (SMDs) were calculated by using a random-effects model to examine the differences in the indexes of obesity, glucometabolic disorder and dyslipidemia between the genotypes of the rs17782313 and rs8192678 polymorphisms. Cochran's Q-statistic test and Begg's test were employed to identify heterogeneity among studies and publication bias, respectively.

Results

Fifty studies (58,716 subjects) and 51 studies (18,660 subjects) were respectively included in the pooled meta-analyses for the rs17782313 and rs8192678 polymorphisms. The C-allele carriers of the rs17782313 polymorphism had a higher average level of body mass index (SMD = 0.21 kg/m2, 95% confidence interval [95% CI] = 0.12 to 0.29 kg/m2, p < 0.001), waist circumference (SMD = 0.14 cm, 95% CI = 0.06 to 0.23 cm, p < 0.001) and blood glucose (SMD = 0.09 mg/dL, 95% CI = 0.02 to 0.16 mg/dL, p = 0.01) than the TT homozygotes. Regarding the rs8192678 polymorphism, no significant associations with the indexes of obesity, glucometabolic disorder and dyslipidemia were detected. However, significant correlations between the rs8192678 polymorphism and multiple glucometabolic indexes were observed in subgroup analyses stratified by sex, age, ethnicity and health status.

Conclusion

The meta-analysis demonstrates that the C allele of the MC4R rs17782313 polymorphism confers a higher risk of obesity and hyperglycemia, and the PGC1α rs8192678 polymorphism is weakly correlated with glucometabolic disorder. These findings may partly explain the relationships between these variants and diabetes as well as cardiovascular disease.

Systematic review registration

https://www.crd.york.ac.uk/prospero/, identifier CRD42022373543.

Computational analysis of spike protein of SARS-CoV-2 (Omicron variant) for development of peptide-based therapeutics and diagnostics

In the last few years, the worldwide population has suffered from the SARS-CoV-2 pandemic. The WHO dashboard indicated that around 504,079,039 people were infected and 6,204,155 died from COVID-19 caused by different variants of SARS-CoV-2. Recently, a new variant of SARS-CoV-2 (B.1.1.529) was reported by South Africa known as Omicron. The high transmissibility rate and resistance towards available anti-SARS-CoV-2 drugs/vaccines/monoclonal antibodies, make Omicron a variant of concern. Because of various mutations in spike protein, available diagnostic and therapeutic treatments are not reliable. Therefore, the present study explored the development of some therapeutic peptides that can inhibit the SARS-CoV-2 virus interaction with host ACE2 receptors and can also be used for diagnostic purposes. The screened linear B cell epitopes derived from receptor-binding domain of spike protein of Omicron variant were evaluated as peptide inhibitor/vaccine candidates through different bioinformatics tools including molecular docking and simulation to analyze the interaction between Omicron peptide and human ACE2 receptor. Overall, in-silico studies revealed that Omicron peptides OP1-P12, OP14, OP20, OP23, OP24, OP25, OP26, OP27, OP28, OP29, and OP30 have the potential to inhibit Omicron interaction with ACE2 receptor. Moreover, Omicron peptides OP20, OP22, OP23, OP24, OP25, OP26, OP27, and OP30 have shown potential antigenic and immunogenic properties that can be used in design and development vaccines against Omicron. Although the in-silico validation was performed by comparative analysis with the control peptide inhibitor, further validation through wet lab experimentation is required before its use as therapeutic peptides.Communicated by Ramaswamy H. Sarma.

Computational study of SENP1 in cancer by novel natural compounds and ZINC database screening

Introduction: Sentrin-specific protease 1 (SENP1) is a protein whose main function is deSUMOylation. SENP1 inhibits apoptosis, and increases angiogenesis, estrogen and androgen receptor transcription and c-Jun transcription factor, proliferation, growth, cell migration, and invasion of cancer. The in vivo and in vitro studies also demonstrated which natural compounds, especially phytochemicals, minerals, and vitamins, prevent cancer. More than 3,000 plant species have been reported in modern medicine. Natural compounds have many anti-cancerous andanti-turmeric properties such as antioxidative, antiangiogenic, antiproliferative, and pro-apoptotic properties. Methods: In this study, we investigated the interaction of some natural compounds with SENP1 to inhibit its activity. We also screened the ZINC database including natural compounds. Molecular docking was performed, and toxicity of compounds was determined; then, molecular dynamics simulation (MDS) and essential dynamics (ED) were performed on natural compounds with higher free binding energies and minimal side effects. By searching in a large library, virtual screening of the ZINC database was performed using LibDock and CDOCKER, and the final top 20 compounds were allowed for docking against SENP1. According to the docking study, the top three leading molecules were selected and further analyzed by MDS and ED. Results: The results suggest that resveratrol (from the selected compounds) and ZINC33916875 (from the ZINC database) could be more promising SENP1 inhibitory ligands. Discussion: Because these compounds can inhibit SENP1 activity, then they can be novel candidates for cancer treatment. However, wet laboratory experiments are needed to validate their efficacy as SENP1 inhibitors.

Predicting the effects of rare genetic variants on oncogenic signaling pathways: A computational analysis of HRAS protein function

The HRAS gene plays a crucial role in regulating essential cellular processes for life, and this gene's misregulation is linked to the development of various types of cancers. Nonsynonymous single nucleotide polymorphisms (nsSNPs) within the coding region of HRAS can cause detrimental mutations that disrupt wild-type protein function. In the current investigation, we have employed in-silico methodologies to anticipate the consequences of infrequent genetic variations on the functional properties of the HRAS protein. We have discovered a total of 50 nsSNPs, of which 23 were located in the exon region of the HRAS gene and denoting that they were expected to cause harm or be deleterious. Out of these 23, 10 nsSNPs ([G60V], [G60D], [R123P], [D38H], [I46T], [G115R], [R123G], [P11OL], [A59L], and [G13R]) were identified as having the most delterious effect based on results of SIFT analysis and PolyPhen2 scores ranging from 0.53 to 69. The DDG values -3.21 kcal/mol to 0.87 kcal/mol represent the free energy change associated with protein stability upon mutation. Interestingly, we identified that the three mutations (Y4C, T58I, and Y12E) were found to improve the structural stability of the protein. We performed molecular dynamics (MD) simulations to investigate the structural and dynamic effects of HRAS mutations. Our results showed that the stable model of HRAS had a significantly lower energy value of -18756 kj/mol compared to the initial model of -108915 kj/mol. The RMSD value for the wild-type complex was 4.40 Å, and the binding energies for the G60V, G60D, and D38H mutants were -107.09 kcal/mol, -109.42 kcal/mol, and -107.18 kcal/mol, respectively as compared to wild-type HRAS protein had -105.85 kcal/mol. The result of our investigation presents convincing corroboration for the potential functional significance of nsSNPs in augmenting HRAS expression and adding to the activation of malignant oncogenic signalling pathways.

Effects of PPARG and PPARGC1A gene polymorphisms on obesity markers

Pediatric obesity presents a multifactorial etiology, which involves genetic traits as well, including single nucleotide polymorphisms. The aim of the study is to investigate the contribution of PPARG gene polymorphisms (namely Pro12Ala rs1801282, His447His rs3856806, and Pro115Gln rs1800571) and PPARGC1A rs8192678 SNP on the anthropometric and metabolic parameters in a population of Romanian children. We conducted a cross-sectional study of 295 Caucasian children, divided according to the body mass index (BMI) z-score into the study (obese and overweight) group of 130 children and the control (normoponderal) group of 165 children. Anthropometric parameters were greater in the obese and overweight population as opposed to controls, with significant differences (p < 0.01) found for the weight (2.77 ± 1.54 SD vs. -0.04 ± 1.15 SD), body mass index (BMI) (2.28 ± 0.97 SD vs. -0.18 ± 1.19 SD), mid-upper arm circumference (MUAC) (4.59 ± 2.28 SD vs. 0.28 ± 3.45 SD), tricipital skin-fold (TSF) (3.31 ± 3.09 SD vs. 0.62 ± 7.28 SD) and waist-to-height ratio (WHtR) (0.61 ± 1.51 SD vs. -0.35 ± 1.35 SD) z-scores. Moreover, triglyceride values were higher in the study group (118.70 ± 71.99 SD vs. 77.09 ± 37.39 SD). No significant difference in the allele and genotype distribution of investigates gene polymorphisms was observed between the studied groups (p > 0.05). PPARG (rs1801282, rs3856806, and rs1800571) were not associated with demographic, anthropometric, and laboratory parameters. However, PPARGC1A rs8192678 CC genotype was associated with TSF z-score (p = 0.03), whereas total and LDL cholesterol levels were significantly higher among TT homozygotes (p < 0.01). Our data suggest that PPARG (rs1801282, rs3856806, and rs1800571) and PPARGC1A (rs8192678) gene polymorphisms were not associated with childhood and adolescence overweight and obesity. The present study identified a significant increase in fasting glucose levels, triglyceride, albumin, and ALT levels in children with excess weight, as well as expected important upward variation of anthropometric parameters (BMI, MUAC, TSF z-scores).

Profile analysis and functional modeling identify circular RNAs in nonalcoholic fatty liver disease as regulators of hepatic lipid metabolism

Nonalcoholic fatty liver disease (NAFLD) is the leading cause of chronic liver disease, associated with an outcome of hepatic fibrosis/cirrhosis and hepatocellular carcinoma. However, limited exploration of the underlying mechanisms hinders its prevention and treatment. To investigate the mechanisms of epigenetic regulation in NAFLD, the expression profile of circular RNA (circRNA) of rodents in which NAFLD was induced by a high-fat, high-cholesterol (HFHC) diet was studied. Modeling of the circRNA-microRNA (miRNA) -mRNA regulatory network revealed the functional characteristics of NAFLD-specific circRNAs. The targets and effects in the liver of such NAFLD-specific circRNAs were further assessed. Our results uncovered that the downregulation of 28 annotated circRNAs characterizes HFHC diet-induced NAFLD. Among the downregulated circRNAs, long intergenic non-protein coding RNA, P53 induced transcript (LNCPINT) -derived circRNAs (circ_0001452, circ_0001453, and circ_0001454) targeted both miR-466i-3p and miR-669c-3p. Their deficiency in NAFLD abrogated the circRNA-based inhibitory effect on both miRNAs, which further inactivated the AMPK signaling pathway via AMPK-α1 suppression. Inhibition of the AMPK signaling pathway promotes hepatic steatosis, depending on the transcriptional and translational upregulation of lipogenic genes, such as those encoding sterol regulatory element-binding protein 1 (SREBP1) and fatty acid synthase (FASN) in hepatocytes. The levels of LNCPINT-derived circRNAs displayed a negative association with hepatic triglyceride (TG) concentration. These findings suggest that loss of LNCPINT-derived circRNAs may underlie NAFLD via miR-466i-3p- and miR-669c-3p-dependent inactivation of the AMPK signaling pathway.

Associations of Polymorphisms in the Peroxisome Proliferator-Activated Receptor Gamma Coactivator-1 Alpha Gene With Subsequent Coronary Heart Disease: An Individual-Level Meta-Analysis

Background: The knowledge of factors influencing disease progression in patients with established coronary heart disease (CHD) is still relatively limited. One potential pathway is related to peroxisome proliferator–activated receptor gamma coactivator-1 alpha (PPARGC1A), a transcription factor linked to energy metabolism which may play a role in the heart function. Thus, its associations with subsequent CHD events remain unclear. We aimed to investigate the effect of three different SNPs in the PPARGC1A gene on the risk of subsequent CHD in a population with established CHD.

Methods: We employed an individual-level meta-analysis using 23 studies from the GENetIcs of sUbSequent Coronary Heart Disease (GENIUS-CHD) consortium, which included participants (n = 80,900) with either acute coronary syndrome, stable CHD, or a mixture of both at baseline. Three variants in the PPARGC1A gene (rs8192678, G482S; rs7672915, intron 2; and rs3755863, T528T) were tested for their associations with subsequent events during the follow-up using a Cox proportional hazards model adjusted for age and sex. The primary outcome was subsequent CHD death or myocardial infarction (CHD death/myocardial infarction). Stratified analyses of the participant or study characteristics as well as additional analyses for secondary outcomes of specific cardiovascular disease diagnoses and all-cause death were also performed.

Results: Meta-analysis revealed no significant association between any of the three variants in the PPARGC1A gene and the primary outcome of CHD death/myocardial infarction among those with established CHD at baseline: rs8192678, hazard ratio (HR): 1.01, 95% confidence interval (CI) 0.98–1.05 and rs7672915, HR: 0.97, 95% CI 0.94–1.00; rs3755863, HR: 1.02, 95% CI 0.99–1.06. Similarly, no significant associations were observed for any of the secondary outcomes. The results from stratified analyses showed null results, except for significant inverse associations between rs7672915 (intron 2) and the primary outcome among 1) individuals aged ≥65, 2) individuals with renal impairment, and 3) antiplatelet users.

Conclusion: We found no clear associations between polymorphisms in the PPARGC1A gene and subsequent CHD events in patients with established CHD at baseline.

Molecular Dynamics Simulation and Essential Dynamics of Deleterious Proline 12 Alanine Single-Nucleotide Polymorphism in PPARγ2 Associated with Type 2 Diabetes, Cardiovascular Disease, and Nonalcoholic Fatty Liver Disease

Background. Peroxisome proliferator-activated receptor-γ (PPARγ) gene is located at 3p25 position. PPARγ functions as the master regulator of glucose homeostasis and lipoprotein metabolism, and recent studies have reported that it is involved in various metabolic diseases such as diabetes mellitus, hyperlipidemia, coronary artery disease (CAD), and nonalcoholic fatty liver disease (NAFLD). PPARγ1 and PPARγ2 are necessary for the development of adipose tissue and insulin sensitivity regulation. But PPARγ2 is the isoform that was controlled in response to nutrient intake and obesity. Methodology. In this study, we used computational techniques to show the impact of Pro12Ala polymorphism on PPARγ2. The 3-D structure of PPARγ2 was modeled using I-TASSER server. The modeled structure was validated with the ZLab server, and the mutation was created with SPDB viewer. Stability prediction tools were used. Molecular dynamics simulation (MDS) was used to understand the structural and functional behavior of the wild type and mutant. Essential dynamics was also applied. Results and Conclusions. Stability prediction tools were showed that this mutation has a destabilizing effect on the PPARγ2 structure. The RMSD, RMSF, Rg, SASA, and DSSP were in line with H-bond results that showed less flexibility in the mutant structure. Essential dynamics was used to verify MDS results. Pro12Ala polymorphism leads to rigidity of the PPARγ2 protein and might disturb the conformational changes and interactions of PPARγ2 and results in type 2 diabetes mellitus (T2DM), CAD, and NAFLD. This study can help scientists to develop a drug therapy against these diseases.