Molecular docking and molecular dynamics study on the effect of ERCC1 deleterious polymorphisms in ERCC1-XPF heterodimer

Enhancement of non-covalent interaction between soy protein isolate and quercetin by sodium alginate

Effects of sodium alginate (SA) on the non-covalent interaction between soybean protein isolate (SPI) and quercetin (Que) were investigated by multispectral technology, molecular docking and dynamics simulation technology. Structural alterations of the binary complexes were observed after SA addition, characterized by a red shift of maximum fluorescence emission wavelength. The introduction of 0.1% (w/v) SA led to a reduction of 12.3% in the α-helix and β-sheet structures, accompanied by 12.6% increase in the β-turn and random coil conformations. The binding of SA to SPI provided electrostatic interactions and facilitated the subsequent binding of SPI to Que. Molecular docking confirmed that hydrophobic interactions and electrostatic interactions were also the main driving force. Molecular dynamics simulation emphasized that the ternary complexes with SA exhibited greater stability compared to the binary ones. The foaming and emulsifying properties of SPI-Que complexes were enhanced by 33.76% and 68.28%, respectively, due to the addition of SA.

In silico design of novel CDK2 inhibitors through QSAR, ADMET, molecular docking and molecular dynamics simulation studies

The present study aims to investigate about the quantitative structure-activity relationship (QSAR) of a series of Thiazole derivatives reported as anticancer agents (hepatocellular carcinoma), using principally the electronic descriptors calculated by the DFT method and by applying the multiple linear regression method. The developed model showed good statistical parameters (R2 = 0.725, R2adj = 0.653, MSE = 0.060, R2test = 0.827, Q2cv = 0.536). The energy EHOMO orbital, electronic energy (TE), shape coefficient (I), number of rotatable bonds (NROT), and index of refraction (n) were revealed to be the main descriptors influencing the anti-cancer activity. Further, new Thiazole derivatives have been designed and their activities and pharmacokinetic properties have been predicted using the validated QSAR model. The designed molecules were then assessed to molecular docking (MD), and molecular dynamic (MDs) simulation accompanied by the calculation of the binding affinity using MMPBSA script according to 100 ns a simulation trajectory, to study both their affinity and their stability towards CDK2 as a target protein for the cancer disease treatment. This research concluded with the identification of four new CDK2 inhibitors which are A1, A3, A5, and A6 showing good pharmacokinetic properties. The MDs results revealed that the newly designed compound A5 remained stable in the active center of the discovered CDK2 protein, indicating its potential as a novel inhibitor for the treatment of hepatocellular carcinoma. The current findings may eventually contribute to the development of robust CDK2 inhibitors in the future.Communicated by Ramaswamy H. Sarma.

The nucleotide excision repair proteins through the lens of molecular dynamics simulations

Mutations that affect the proteins responsible for the nucleotide excision repair (NER) pathway can lead to diseases such as xeroderma pigmentosum, trichothiodystrophy, Cockayne syndrome, and Cerebro-oculo-facio-skeletal syndrome. Hence, understanding their molecular behavior is needed to elucidate these diseases' phenotypes and how the NER pathway is organized and coordinated. Molecular dynamics techniques enable the study of different protein conformations, adaptable to any research question, shedding light on the dynamics of biomolecules. However, as important as they are, molecular dynamics studies focused on DNA repair pathways are still becoming more widespread. Currently, there are no review articles compiling the advancements made in molecular dynamics approaches applied to NER and discussing: (i) how this technique is currently employed in the field of DNA repair, focusing on NER proteins; (ii) which technical setups are being employed, their strengths and limitations; (iii) which insights or information are they providing to understand the NER pathway or NER-associated proteins; (iv) which open questions would be suited for this technique to answer; and (v) where can we go from here. These questions become even more crucial considering the numerous 3D structures published regarding the NER pathway's proteins in recent years. In this work, we tackle each one of these questions, revising and critically discussing the results published in the context of the NER pathway.

A computational study of structural analysis of Class I human glucose-6-phosphate dehydrogenase (G6PD) variants: Elaborating the correlation to chronic non-spherocytic hemolytic anemia (CNSHA)

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common human enzyme defect that affects more than 500 million people worldwide. Individuals affected with G6PD deficiency may occasionally suffer mild-to-severe chronic hemolytic anemia. Chronic non-spherocytic hemolytic anemia (CNSHA) is a potential result of the Class I G6PD variants. This comparative computational study attempted to correct the defect in variants structure by docking the AG1 molecule to selected Class I G6PD variants [G6PD^Nashville (Arg393His), G6PD^Alhambra (Val394Leu), and G6PD^Durham (Lys238Arg)] at the dimer interface and structural NADP⁺ binding site. It was followed by an analysis of the enzyme conformations before and after binding to the AG1 molecule using the molecular dynamics simulation (MDS) approach, while the severity of CNSHA was determined via root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), hydrogen bonds, salt bridges, radius of gyration (Rg), solvent accessible surface area analysis (SASA), and principal component analysis (PCA). The results revealed that G6PD^Nashville (Arg393His) and G6PD^Durham (Lys238Arg) had lost the direct contact with structural NADP⁺ and salt bridges at Glu419 - Arg427 and Glu206 - Lys407 were disrupted in all selected variants. Furthermore, the AG1 molecule re-stabilized the enzyme structure by restoring the missing interactions. Bioinformatics approaches were also used to conduct a detailed structural analysis of the G6PD enzyme at a molecular level to understand the implications of these variants toward enzyme function. Our findings suggest that despite the lack of treatment for G6PDD to date, AG1 remains a novel molecule that promotes activation in a variety of G6PD variants.

Dynamic study of small toxic hydrophobic proteins PepA1 and PepG1 of Staphylococcus aureus

Toxin-antitoxin (TA) systems are small genetic elements which encode toxin proteins that interfere with vital cellular functions. PepA1 and PepG1 toxin proteins, known also as SprA1 and SprG1, are type I TA. In Staphylococcus aureus (S. aureus), their expression without the antitoxin counterparts (SprA1_AS and SprF1), is lethal to the pathogen. Molecular Dynamics (MD) simulation was performed for PepA1 and PepG1 to understand their dynamic state, conformational changes, and their toxicity. The protein structures were constructed and used for MD simulation and the conformational changes, stability, flexibility, fluctuations, hydrophobicity, and role of their dynamic state on function prediction were studied extensively by GROMACS MD simulation analysis tools. In silico study indicated that the PepA1 and PepG1 proteins change their structural conformation from an open to closed state where PepA1 conformational changes were faster (10 ns) than PepG1 (20 ns) while PepG1 exerted more stability and flexibility than PepA1. According to SASA values, PepG1 is more hydrophobic than the PepA1 and forms fewer hydrogen bonds than PepA1. The in vivo study with PepA1 and PepG1 proteins provided evidence that both the conformation changes between the open and closed states and the amino acid sequence are crucial for peptide toxicity.

Impact of Reactive Oxygen and Nitrogen Species Produced by Plasma on Mdm2-p53 Complex

The study of protein–protein interactions is of great interest. Several early studies focused on the murine double minute 2 (Mdm2)–tumor suppressor protein p53 interactions. However, the effect of plasma treatment on Mdm2 and p53 is still absent from the literature. This study investigated the structural changes in Mdm2, p53, and the Mdm2–p53 complex before and after possible plasma oxidation through molecular dynamic (MD) simulations. MD calculation revealed that the oxidized Mdm2 bounded or unbounded showed high flexibility that might increase the availability of tumor suppressor protein p53 in plasma-treated cells. This study provides insight into Mdm2 and p53 for a better understanding of plasma oncology.

XPF -673C>T variation is associated with the susceptibility to breast cancer

PURPOSE

XPF variations might decrease the DNA repair capacity and further contribute to cancer development. This study aimed to investigate the association of XPF polymorphisms with risk of developing breast cancer.

METHODS

TCGA, the Human Protein Atlas and Kaplan-Meier plotter were used to analyze the expression of XPF in breast cancer tissues and its effect on the survival of breast cancer patients. The expression of XPF in breast cancer tissues was detected by qRT-PCR. This case-control study included 467 breast cancer patients and 467 healthy controls. The genotype of genetic variation was detected by polymerase chain reaction restriction fragment length polymorphism. Odds ratios and 95 % confidence intervals were calculated. Correlations between XPF variation and clinicopathological parameters were assessed through Kendall's Tau-b test. The relationship between XPF gene function variation and XPF gene expression was analyzed by GTEx.

RESULTS

The expression of XPF in breast cancer tissues is higher than that in normal tissues. Breast cancer patients with high XPF expression have a higher relapse free survival rate (HR = 0.88, 95 % CI = 0.80-0.97), but have no effect on the overall survival rate (logrank P = 0.28). XPF -673C > T variant can reduce the risk of breast cancer patients (OR = 0.35, 95 %CI = 0.20-0.63 for codominant mode; OR = 0.66, 95 %CI = 0.51-0.85 for dominant model; OR = 0.40, 95 %CI = 0.23-0.70 for recessive model). The XPF 11985 GG genotype reduced the risk of early breast cancer (OR = 0.49, 95 %CI = 0.24-0.97), but not the risk of advanced breast cancer (OR = 1.20, 95 % CI = 0.58-2.48). XPF 11985A > G variant can also reduce the risk of ERBB2 expression in patients (OR = 0.50, 95 %CI = 0.27-0.94). There is no correlation between XPF -673C > T/XPF11985A > G variants and ER and PR. XPF -673C > T variant can reduce XPF expression (P < 0.05).

CONCLUSIONS

Genetic variations of XPF gene may affect its expression and the risk of breast cancer in the Chinese population.

Identification and structural characterization of deleterious non-synonymous single nucleotide polymorphisms in the human SKP2 gene

In SCF (Skp, Cullin, F-box) ubiquitin-protein ligase complexes, S-phase kinase 2 (SKP2) is one of the major players of F-box family, that is responsible for the degradation of several important cell regulators and tumor suppressor proteins. Despite of having significant evidence for the role of SKP2 on tumorgenesis, there is a lack of available data regarding the effect of non-synonymous polymorphisms. In this communication, the structural and functional consequences of non-synonymous single nucleotide polymorphisms (nsSNPs) of SKP2 have been reported by employing various computational approaches and molecular dynamics simulation. Initially, several computational tools like SIFT, PolyPhen-2, PredictSNP, I-Mutant 2.0 and ConSurf have been implicated in this study to explore the damaging SNPs. In total of 172 nsSNPs, 5 nsSNPs were identified as deleterious and 3 of them were predicted to be decreased the stability of protein. Guided from ConSurf analysis, P101L (rs761253702) and Y346C (rs755010517) were categorized as the highly conserved and functional disrupting mutations. Therefore, these mutations were subjected to three dimensional model building and molecular dynamics simulation study for the detailed structural consequences upon the mutations. The study revealed that P101L and Y346C mutations increased the flexibility and changed the structural dynamics. As both these mutations are located in the most functional regions of SKP2 protein, these computational insights might be helpful to consider these nsSNPs for wet-lab confirmatory analysis as well as in rationalizing future population based studies and structure based drug design against SKP2.

An integrated in silico approach for functional and structural impact of non- synonymous SNPs in the MYH1 gene in Jeju Native Pigs

BACKGROUND

This study was performed to identify the non- synonymous polymorphisms in the myosin heavy chain 1 gene (MYH1) association with skeletal muscle development in economically important Jeju Native Pig (JNP) and Berkshire breeds. Herein, we present an in silico analysis, with a focus on (a) in silico approaches to predict the functional effect of non-synonymous SNP (nsSNP) in MYH1 on growth, and (b) molecular docking and dynamic simulation of MYH1 to predict the effects of those nsSNP on protein-protein association.

RESULTS

The NextGENe (V 2.3.4.) tool was used to identify the variants in MYH1 from JNP and Berkshire using RNA seq. Gene ontology analysis of MYH1 revealed significant association with muscle contraction and muscle organ development. The 95 % confidence intervals clearly indicate that the mRNA expression of MYH1 is significantly higher in the Berkshire longissimus dorsi muscle samples than JNP breed. Concordant in silico analysis of MYH1, the open-source software tools identified 4 potential nsSNP (L884T, K972C, N981G, and Q1285C) in JNP and 1 nsSNP (H973G) in Berkshire pigs. Moreover, protein-protein interactions were studied to investigate the effect of MYH1 mutations on association with hub proteins, and MYH1 was found to be closely associated with the protein myosin light chain, phosphorylatable, fast skeletal muscle MYLPF. The results of molecular docking studies on MYH1 (native and 4 mutants) and MYLFP demonstrated that the native complex showed higher electrostatic energy (-466.5 Kcal mol(-1)), van der Walls energy (-87.3 Kcal mol(-1)), and interaction energy (-835.7 Kcal mol(-1)) than the mutant complexes. Furthermore, the molecular dynamic simulation revealed that the native complex yielded a higher root-mean-square deviation (0.2-0.55 nm) and lower root-mean-square fluctuation (approximately 0.08-0.3 nm) as compared to the mutant complexes.

CONCLUSIONS

The results suggest that the variants at L884T, K972C, N981G, and Q1285C in MYH1 in JNP might represent a cause for the poor growth performance for this breed. This study is a pioneering in-depth in silico analysis of polymorphic MYH1 and will serve as a valuable resource for further targeted molecular diagnosis and population-based studies conducted for improving the growth performance of JNP.

XPF-673C>T polymorphism effect on the susceptibility to esophageal cancer in Chinese population

Purpose

Xeroderma pigmentsum group F (XPF) plays a pivotal role in DNA nucleotide excision repair and has been linked to the development of various cancers. This study aims to assess the association of XPF genetic variants with the susceptibility to esophageal squamous cell carcinoma (ESCC) in Chinese population.

Methods

This two-stage case-control study was conducted in a total of 1524 patients with ESCC and 1524 controls. Genotype of XPF -673C>T and 11985A>G variants were determined by polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP). Logistic regression analysis was performed to estimate odd ratios (ORs) and 95% confidence intervals (95% CI).

Results

Our case-control study showed that XPF -673TT genotype was associated with a decreased risk of ESCC compared with CC genotype in both case-control sets (Tangshan set: OR = 0.58; 95%CI = 0.34–0.99, P = 0.040; Beijing set: OR = 0.66; 95%CI = 0.46–0.95, P = 0.027). Stratified analyses revealed that a multiplicative interaction between -673C>T variant and age, sex or smoking status was evident (Gene-age: P_interaction = 0.002; Gene-sex: P_interaction = 0.002; Gene-smoking: P_interaction = 0.002). For XPF 11985A>G polymorphism, there was no significant difference of genotype distribution between ESCC cases and controls.

Conclusion

These findings indicated that genetic variants in XPF might contribute to the susceptibility to ESCC.