Integrative Bioinformatics Analysis: Unraveling Variant Signatures and Single-Nucleotide Polymorphism Markers Associated with 5-FU-Based Chemotherapy Resistance in Colorectal Cancer Patients.

BACKGROUND: Drug resistance in colorectal cancer (CRC) is modulated by multiple molecular factors, which can be ascertained through genetic investigation. Single nucleotide polymorphisms (SNPs) within key genes have the potential to impair the efficacy of chemotherapeutic agents such as 5-fluorouracil (5-FU). Therefore, the identification of SNPs linked to drug resistance can significantly contribute to the advancement of tailored therapeutic approaches and the enhancement of treatment outcomes in patients with CRC. MATERIAL AND METHOD: To identify dysregulated genes in 5-FU-based chemotherapy responder or non-responder CRC patients, a meta-analysis was performed. Next, the protein-protein interaction (PPI) network of the identified genes was analyzed using the STRING database. The most significant module was chosen for further analysis. In addition, a literature review was conducted to identify drug resistance-related genes. Enrichment analysis was conducted to validate the main module genes and the genes identified from the literature review. The associations between SNPs and drug resistance were investigated, and the consequences of missense variants were assessed using in silico tools. RESULT: The meta-analysis identified 796 dysregulated genes. Then, to conduct PPI analysis and enrichment analysis, we were able to discover 23 genes that are intricately involved in the cell cycle pathway. Consequently, these 23 genes were chosen for SNP analysis. By using the dbSNP database and ANNOVAR, we successfully detected and labeled SNPs in these specific genes. Additionally, after careful exclusion of SNPs with allele frequencies below 0.01, we evaluated 6 SNPs from the HDAC1, MCM2, CDK1, BUB1B, CDC14B, and CCNE1 genes using 8 bioinformatics tools. Therefore, these SNPs were identified as potentially harmful by multiple computational tools. Specifically, rs199958833 in CDK1 (Val124Gly) was predicted to be damaging by all tools used. Our analysis strongly indicates that this specific SNP could negatively affect the stability and functionality of the CDK1 protein. CONCLUSION: Based on our current understanding, the evaluation of CDK1 polymorphisms in the context of drug resistance in CRC has yet to be undertaken. In this investigation, we showed that rs199958833 variant in the CDK1 gene may favor resistance to 5-FU-based chemotherapy. However, these findings need validation in an independent cohort of patients.

Identification of high-risk non-synonymous SNPs (nsSNPs) in DNAH1 and DNAH17 genes associated with male infertility: a bioinformatics analysis.

Male infertility is a significant reproductive issue affecting a considerable number of couples worldwide. While there are various causes of male infertility, genetic factors play a crucial role in its development. We focused on identifying and analyzing the high-risk nsSNPs in DNAH1 and DNAH17 genes, which encode proteins involved in sperm motility. A total of 20 nsSNPs for DNAH1 and 10 nsSNPs for DNAH17 were analyzed using various bioinformatics tools including SIFT, PolyPhen-2, CADD, PhD-SNPg, VEST-4, and MutPred2. As a result, V1287G, L2071R, R2356W, R3169C, R3229C, E3284K, R4096L, R4133C, and A4174T in DNAH1 gene and C1803Y, C1829Y, R1903C, and L3595P in DNAH17 gene were identified as high-risk nsSNPs. These nsSNPs were predicted to decrease protein stability, and almost all were found in highly conserved amino acid positions. Additionally, 4 nsSNPs were observed to alter post-translational modification status. Furthermore, the interaction network analysis revealed that DNAH1 and DNAH17 interact with DNAH2, DNAH3, DNAH5, DNAH7, DNAH8, DNAI2, DNAL1, CFAP70, DNAI3, DNAI4, ODAD1, and DNAI7, demonstrating the importance of DNAH1 and DNAH17 proteins in the overall functioning of the sperm motility machinery. Taken together, these findings revealed the detrimental effects of identified high-risk nsSNPs on protein structure and function and highlighted their potential relevance to male infertility. Further studies are warranted to validate these findings and to elucidate the underlying mechanisms.

The effects of photobiomodulation on the improvement of sperm parameters: A review study.

The prevalence of male infertility has become a significant clinical concern worldwide, with a noticeable upward trend in recent times. The rates of fertilization and subsequent development of embryos are dependent on many parameters associated with the quality and viability of sperm. Photobiomodulation (PBM) is a promising approach with a great potential for translational applications in the treatment of spermatozoa exhibiting low quality and motility. In this study, a comprehensive analysis of the existing literature, specifically examining the mechanisms of action of PBM has been presented. Our objective was to enhance knowledge in the field of laser light therapy in order to promote the usage of irradiation in clinical settings in a more effective way. Within the realm of reproductive science, the utilization of PBM has been employed to enhance the metabolic processes, motility, and viability of spermatozoa. This is attributed to its advantageous effects on mitochondria, resulting in the activation of the mitochondrial respiratory chain and subsequent synthesis of ATP. This therapeutic approach can be highly advantageous in circumventing the reliance on chemical substances within the culture medium for spermatozoa while also facilitating the viability and motility of spermatozoa, particularly in circumstances involving thawing or samples with significant immotility.

In silico screening and analysis of nonsynonymous SNPs in human CYP1A2 to assess possible associations with pathogenicity and cancer susceptibility.

Cytochrome P450 1A2 (CYP1A2) is one of the main hepatic CYPs involved in metabolism of carcinogens and clinically used drugs. Nonsynonymous single nucleotide polymorphisms (nsSNPs) of this enzyme could affect cancer susceptibility and drug efficiency. Hence, identification of human CYP1A2 pathogenic nsSNPs could be of great importance in personalized medicine and pharmacogenetics. Here, 176 nsSNPs of human CYP1A2 were evaluated using a variety of computational tools, of which 18 nsSNPs were found to be associated with pathogenicity. Further analysis suggested possible association of 9 nsSNPs (G73R, G73W, R108Q, R108W, E168K, E346K, R431W, F432S and R456H) with the risk of hepatocellular carcinoma. Molecular dynamics simulations revealed higher overall flexibility, decreased intramolecular hydrogen bonds and lower content of regular secondary structures for both cancer driver variants G73W and F432S when compared to the wild-type structure. In case of F432S, loss of the conserved hydrogen bond between Arg137 and heme propionate oxygen may affect heme stability and the observed significant rise in fluctuation of the CD loop could modify CYP1A2 interactions with its redox partners. Together, these findings propose CYP1A2 as a possible candidate for hepatocellular carcinoma and provide structural insights into how cancer driver nsSNPs could affect protein structure, heme stability and interaction network.

Purification, Characterization and Mechanistic Evaluation of Angiotensin Converting Enzyme Inhibitory Peptides Derived from Zizyphus Jujuba Fruit.

The synthetic Angiotensin Converting Enzyme (ACE) inhibitors have side effects and hence demands for natural ACE inhibitors have been rising. The aim of this study is to purify and introduce natural ACE inhibitors extracted from Zizyphus jujuba fruits. Proteins from Zizyphus jujuba were lysed by trypsin, papain and their combination. Acquired peptides were purified and evaluated for ACE inhibitory activity. Peptide fractions with inhibitory activity were sequenced using tandem mass spectrometry. To elucidate the mode of peptide binding to ACE, homology modeling, molecular docking and molecular dynamics simulation were performed. Amino acid sequence of F2 and F4 peptides, which were the most active hydrolysates, were determined to be IER and IGK with the IC50 values of 0.060 and 0.072 mg/ml, respectively. Results obtained by computational analysis revealed that similar to the common ACE competitive inhibitors such as captopril, IER tripeptide binds to the enzyme active site, in vicinity of the zinc binding site, and occupies the S1 and S2' subsites. Binding occurs through hydrogen bonding with Gln293, Lys522, His524, Tyr531 and also several hydrophobic interactions. Collectively, these findings indicate that IER tripeptide inhibits the rabbit ACE enzyme through a competitive mechanism of inhibition with IC50 values in the millimolar range.

Structure and dynamics of inactive and active MARK4: conformational switching through the activation process.

Microtubule affinity-regulating kinase 4 (MARK4) is a serine/threonine protein kinase belonging to a highly-conserved group of PAR proteins that phosphorylate microtubule-associated proteins (tau, MAP2 and MAP4) and regulate cell polarity. MARK4 isoform is mainly found in brain tissue, causing microtubule destabilization in neuronal cells and also tau-protein phosphorylation seen in Alzheimer's disease. In this study, the dynamic behavior of inactive and active structures of human MARK4 was studied by modeling and molecular dynamics simulations and motions of the protein through activation process were interpolated. Structure and dynamics of the protein active state in presence of Mg-ATP were also studied. The results suggest for occlusion of ATP binding site by activation loop, as the main inactivation mechanism. Data also justify the necessity of UBA (ubiquitin-associated) domain auto-inhibitory role. Within the inactive and active state, G-loop is highly fluctuating and assumes an open conformation and away from Mg-ATP complex. This behavior raises the assumption that this loop may experience other stabilizing interactions with residues out of the kinase core which help its stabilization during the phospho-transfer reaction. Mg-ATP complex localization is well preserved within the catalytic cleft through the electrostatic interactions. However, minor rearrangements of water molecules around this complex are allowed, which may further refine the delicate position of this complex for phospho-transfer reaction.Communicated by Ramaswamy H. Sarma.

Interconversion of inactive to active conformation of MARK2: Insights from molecular modeling and molecular dynamics simulation.

The Ser/Thr protein kinase MARK2, also known as Par1b, belongs to the highly-conserved family of PAR proteins which regulate cell polarity and partitioning through the animal kingdom. In the current study, inactive and active structures of human MARK2 were constructed by modeling and molecular dynamics simulation, based on available incomplete crystal structures in Protein Data Bank, to investigate local structural changes through which MARK2 switches from inactive to active state. None of the MARK2 wild type inactive crystal structures represent the position of activation segment. So, the contribution of this loop to the formation of inactive state is not clear. In the modeled structure of inactive MARK2, activation segment occludes the enzyme active site and assumes a relatively stable position. We also presented a detailed description of the major structural changes occur through the activation process and proposed a framework on how these deviations might be affected by the phosphorylation of Thr208 or existence of the UBA domain. Inspection of protein active state in the presence of Mg-ATP, demonstrated the precise arrangement of the various parts of enzyme around Mg-ATP and the importance of their stability in localization of the resulting complex. The results also confirmed the alleged mild auto-inhibitory role of the UBA domain and suggested a reason for the necessity of this role, based on structural similarities to other related kinases.

Structural insights into the effects of charge-reversal substitutions at the surface of horseradish peroxidase.

Horseradish peroxidase (HRP), has gained significant interests in biotechnology, especially in biosensor field and diagnostic test kits. Hence, its solvent-exposed lysine residues 174, 232, and 241 have been frequently modified with the aim of improving its stability and catalytic efficiency. In this computational study, we investigated the effects of Lys-to-Glu substitutions on HRP structure to model charge-reversal manipulations at the enzyme surface. Simulation results implied that upon these substitutions, the number of stable hydrogen bonds and α-helical content of HRP are increased and the proximal Ca2+ binding pocket becomes more integrated. The results revealed that although Glu174-heme hydrogen bond is lost after mutation, formation of a new hydrogen bonding network contributes to the stability of heme-protein linkage. Together, it may be concluded that these substitutions enhance the stability of the protein moiety as well as the heme-protein non-covalent interactions. In the enzyme active site, we observed increased accessibility of peroxide binding site and heme prosthetic group to the peroxide and aromatic substrates, respectively. Results also demonstrated that the bottleneck entry of the peroxide-binding site has become wider and more flexible upon substitutions. Moreover, the hydrophobic patch functioning as a binding site or trap for reducing aromatic substrates is more extended in mutated enzyme. These observations suggest that the reactivity of the enzyme to its substrates has increased. Together, the results of this simulation study could provide possible structural clues to explain those experimental observations in which the protein stability achieved upon manipulation of charge distribution on protein surface.

How modification of accessible lysines to phenylalanine modulates the structural and functional properties of horseradish peroxidase: a simulation study.

Horseradish Peroxidase (HRP) is one of the most studied peroxidases and a great number of chemical modifications and genetic manipulations have been carried out on its surface accessible residues to improve its stability and catalytic efficiency necessary for biotechnological applications. Most of the stabilized derivatives of HRP reported to date have involved chemical or genetic modifications of three surface-exposed lysines (K174, K232 and K241). In this computational study, we altered these lysines to phenylalanine residues to model those chemical modifications or genetic manipulations in which these positively charged lysines are converted to aromatic hydrophobic residues. Simulation results implied that upon these substitutions, the protein structure becomes less flexible. Stability gains are likely to be achieved due to the increased number of stable hydrogen bonds, improved heme-protein interactions and more integrated proximal Ca2+ binding pocket. We also found a new persistent hydrogen bond between the protein moiety (F174) and the heme prosthetic group as well as two stitching hydrogen bonds between the connecting loops GH and F'F″ in mutated HRP. However, detailed analysis of functionally related structural properties and dynamical features suggests reduced reactivity of the enzyme toward its substrates. Molecular dynamics simulations showed that substitutions narrow the bottle neck entry of peroxide substrate access channel and reduce the surface accessibility of the distal histidine (H42) and heme prosthetic group to the peroxide and aromatic substrates, respectively. Results also demonstrated that the area and volume of the aromatic-substrate binding pocket are significantly decreased upon modifications. Moreover, the hydrophobic patch functioning as a binding site or trap for reducing aromatic substrates is shrunk in mutated enzyme. Together, the results of this simulation study could provide possible structural clues to explain those experimental observations in which the protein stability achieved concurrent with a decrease in enzyme activity, upon manipulation of charge/hydrophobicity balance at the protein surface.