Structure and dynamics of γ-secretase with presenilin 2 compared to presenilin 1

Severe early-onset familial Alzheimer's disease (FAD) is caused by more than 200 different mutations in the genes coding for presenilin, the catalytic subunit of the 4-subunit protease complex γ-secretase, which cleaves the C99 fragment of the amyloid precursor protein (APP) to produce Aβ peptides. γ-Secretase exists with either of two homologues, PS1 and PS2. All cryo-electron microscopic structures and computational work has so far focused on γ-secretase with PS1, yet PS2 mutations also cause FAD. A central question is thus whether there are structural and dynamic differences between PS1 and PS2. To address this question, we use the cryo-electron microscopic data for PS1 to develop the first structural and dynamic model of PS2-γ-secretase in the catalytically relevant mature membrane-bound state at ambient temperature, equilibrated by three independent 500 ns molecular dynamics simulations. We find that the characteristic nicastrin extra-cellular domain breathing mode and major movements in the cytosolic loop between TM6 and TM7 occur in both PS2- and PS1-γ-secretase. The overall structures and conformational states are similar, suggesting similar catalytic activities. However, at the sequence level, charge-controlled membrane-anchoring is extracellular for PS1 and intracellular for PS2, which suggests different subcellular locations. The tilt angles of the TM2, TM6, TM7 and TM9 helices differ in the two forms of γ-secretase, suggesting that the two proteins have somewhat different substrate processing and channel sizes. Our MD simulations consistently indicated that PS2 retains several water molecules near the catalytic site at the bilayer, as required for catalysis. The possible reasons for the differences of PS1 and PS2 are discussed in relation to their location and function.

Computational investigation on effect of mutations in PCNA resulting in structural perturbations and inhibition of mismatch repair pathway

From bacteria to mammals, DNA mismatch repair (MMR) pathway plays an essential role in eliminating mismatched nucleotides and insertion-deletion mismatches during the process of DNA replication. Among many of the proteins which participate in the mismatch repair process, proliferating cell nuclear antigen (PCNA) remains the principal conductor at the replication fork. The pol30-201 and pol30-204 are the two mutated alleles which encode for C22Y and C81R mutant forms of PCNA proteins. We performed long term molecular dynamics (MD) simulations analysis (0.8 μs) to understand the dynamic behavior and alterations in the structure of wild type and mutated forms of PCNA at the atomic level. We observed changes in the structural characteristics like length, radius, rise per residue of alpha helices in both the mutated forms of PCNA. Apart from it, disfigurement of the charge distribution which effects binding with the dsDNA due to mutant C22Y and other structural perturbations were also seen in regions significant for the formation of a biologically active trimeric form of PCNA due to mutant C81R. Our analysis of native and mutated forms of PCNA provides an insight into the essential structural and functional features required for proper and well-coordinated DNA mismatch repair process and consequences of the mutation leading to an impaired process of MMR. These structural characteristics are fundamental for the MMR process and hence our analysis likely contributes to or presents the novel mechanism involved in the process of MMR.Communicated by Ramaswamy H. Sarma.

Effects of Pro1266Leu mutation on structure and function of glycoprotein Ib binding domain of von Willebrand factor

Glycoprotein Ibα (GpIbα) binding ability of A1 domain of von Willebrand factor (vWF) facilitates platelet adhesion that plays a crucial role in maintaining hemostasis and thrombosis at the site of vascular damage. There are both "loss as well as gain of function" mutations observed in this domain. Naturally occurring "gain of function" mutations leave self-activating impacts on the A1 domain which turns the normal binding to characteristic constitutive binding with GPIbα. These "gain of function" mutations are associated with the von Willebrand disease type 2B. In recent years, studies focused on understanding the mechanism and conformational patterns attached to these phenomena have been conducted, but the conformational pathways leading to such binding patterns are poorly understood as of now. To obtain a microscopic picture of such events for the better understanding of pathways, we used molecular dynamics (MD) simulations along with principal component analysis and normal mode analysis to study the effects of Pro1266Leu (Pro503Leu in structural context) mutation on the structure and function of A1 domain of vWF. MD simulations have provided atomic-level details of intermolecular motions as a function of time to understand the dynamic behavior of A1 domain of vWF. Comparative analysis of the trajectories obtained from MD simulations of both the wild type and Pro503Leu mutant suggesting appreciable conformational changes in the structure of mutant which might provide a basis for assuming the "gain of function" effects of these mutations on the A1 domain of vWF, resulting in the constitutive binding with GpIbα.

Insights into the mode of flavin mononucleotide binding and catalytic mechanism of bacterial chromate reductases: A molecular dynamics simulation study

Enzymes from natural sources protect the environment via complex biological mechanisms, which aid in reductive immobilization of toxic metals including chromium. Nevertheless, progress was being made in elucidating high-resolution crystal structures of reductases and their binding with flavin mononucleotide (FMN) to understand the underlying mechanism of chromate reduction. Therefore, herein, we employed molecular dynamics (MD) simulations, principal component analysis (PCA), and binding free energy calculations to understand the dynamics behavior of these enzymes with FMN. Six representative chromate reductases in monomeric and dimeric forms were selected to study the mode, dynamics, and energetic component that drive the FMN binding process. As evidenced by MD simulation, FMN prefers to bind the cervix formed between the catalytic domain surrounded by strong conserved hydrogen bonding, electrostatic, and hydrophobic contacts. The slight movement and reorientation of FMN resulted in breakage of some crucial H-bonds and other nonbonded contacts, which were well compensated with newly formed H-bonds, electrostatic, and hydrophobic interactions. The critical residues aiding in tight anchoring of FMN within dimer were found to be strongly conserved in the bacterial system. The molecular mechanics combined with the Poisson-Boltzmann surface area binding free energy of the monomer portrayed that the van der Waals and electrostatic energy contribute significantly to the total free energy, where, the polar solvation energy opposes the binding of FMN. The proposed proximity relationships between enzyme and FMN binding site presented in this study will open up better avenues to engineer enzymes with optimized chromate reductase activity for sustainable bioremediation of heavy metals.

Identification of prognostic biomarkers for breast cancer based on miRNA and mRNA co-expression network

PURPOSE

Breast cancer (BC) remains a serious health threat for women due to its high incidence and the trend of rejuvenation. Accumulating evidence has highlighted that microRNAs (miRNAs) and messenger RNAs (mRNAs) could play important roles in various biological processes involved in the pathogenesis of BC. The present study aimed to identify potential prognostic biomarkers associated with BC.

METHODS

Here, original gene expression profiles of patients with BC was downloaded from The Cancer Genome Atlas (TCGA) database. TargetScan, miRDB, and miRTarBase databases were used to predict the target genes of prognostic-related differentially expressed miRNAs (DEMs). Subsequently, functional enrichment analysis and topological analysis were performed on the overlaps of target genes and differentially expressed mRNAs (DEGs), and Kaplan-Meier analysis was used to predict prognosis-related target genes to identify prognostic biomarkers.

RESULTS

A total of 218 DEMs and 2222 DEGs were extracted in which eight miRNAs were associated with prognosis, and 278 target DEGs were screened out incorporated into functional enrichment analysis and protein-protein interaction network visualization studies. Additionally, five hub genes (CXCL12, IGF1, LEF1, MMP1, and RACGAP1) were observed as potential biomarkers for BC prognosis through survival analysis.

CONCLUSION

We performed a distinctive correlation analysis of miRNA-mRNA in BC patients, and identified eight miRNAs and five hub genes may be effective biomarkers for the prognosis of BC patients.

Identification of natural compound inhibitors against PfDXR: A hybrid structure-based molecular modeling approach and molecular dynamics simulation studies

In the present contribution, multicomplex-based pharmacophore studies were carried out on the structural proteome of Plasmodium falciparum 1-deoxy-D-xylulose-5-phosphate reductoisomerase. Among the constructed models, a representative model with complementary features, accountable for the inhibition was used as a primary filter for the screening of database molecules. Auxiliary evaluations of the screened molecules were performed via drug-likeness and molecular docking studies. Subsequently, the stability of the docked inhibitors was envisioned by molecular dynamics simulations, principle component analysis, and molecular mechanics-Poisson-Boltzmann surface area-based free binding energy calculations. The stability assessment of the hits was done by comparing with the reference (beta-substituted fosmidomycin analog, LC5) to prioritize more potent candidates. All the complexes showed stable dynamic behavior while three of them displayed higher binding free energy compared with the reference. The work resulted in the identification of the compounds with diverse scaffolds, which could be used as initial leads for the design of novel PfDXR inhibitors.

Molecular interaction studies of Deguelin and its derivatives with Cyclin D1 and Cyclin E in cancer cell signaling pathway: The computational approach

Deguelin is a major active ingredient and principal component in several plants and it is a potential molecule to target proteins of cancer cell signaling pathway. As a complex natural extract, deguelin interacts with various molecular targets to exert its anti-tumor properties at nanomolar level. It induces cell apoptosis by blocking anti-apoptotic pathways, while inhibiting tumor cell multiplication and malignant transformation through p27-cyclin-E-pRb-E2F1- cell cycle control and HIF-1alphaVEGF antiangiogenic pathways. In silico studies of deguelin and its derivatives is performed to explore interactions with Cyclin D1 and Cyclin E to understand the molecular insights of derivatives with the receptors. Deguelin and its derivatives are minimized by Avogadro to achieve stable conformation. All docking simulation are performed with AutoDockVina and virtual screening of docked ligands are carried out based on binding energy and number of hydrogen bonds. Molecular dynamics (MD) and Simulation of Cyclin D1 and Cyclin E1 is performed for 100 ns and stable conformation is obtained at 78 ns and 19 ns respectively. Ligands thus obtained from docking studies may be probable target to inhibit cancer cell signaling pathways.

Styryl-cinnamate hybrid inhibits glioma by alleviating translation, bioenergetics and other key cellular responses leading to apoptosis

Gliomas are lethal and aggressive form of brain tumors with resistance to conventional radiation and cytotoxic chemotherapies; inviting continuous efforts for drug discovery and drug delivery. Interestingly, small molecule hybrids are one such pharmacophore that continues to capture interest owing to their pluripotent medicinal effects. Accordingly, we earlier reported synthesis of potent Styryl-cinnamate hybrids (analogues of Salvianolic acid F) along with its plausible mode of action (MOA). We explored iTRAQ-LC/MS-MS technique to deduce differentially expressed landscape of native & phospho-proteins in treated glioma cells. Based on this, Protein-Protein Interactome (PPI) was looked into by employing computational tools and further validated in vitro. We hereby report that the Styryl-cinnamate hybrid, an analogue of natural Salvianolic acid F, alters key regulatory proteins involved in translation, cytoskeleton development, bioenergetics, DNA repair, angiogenesis and ubiquitination. Cell cycle analysis dictates arrest at G0/G1 stage along with reduced levels of cyclin D; involved in G1 progression. We discovered that Styryl-cinnamate hybrid targets glioma by intrinsically triggering metabolite-mediated stress. Various oncological circuits alleviated by the potential drug candidate strongly supports the role of such pharmacophores as anticancer drugs. Although, further analysis of SC hybrid in treating xenografts or solid tumors is yet to be explored but their candidature has gained huge impetus through this study. This study equips us better in understanding the shift in proteomic landscape after treating glioma cells with SC hybrid. It also allows us to elicit molecular targets of this potential drug before progressing to preclinical studies.

Impact of missense mutations in survival motor neuron protein (SMN1) leading to Spinal Muscular Atrophy (SMA): A computational approach

Spinal muscular atrophy (SMA) is a neuromuscular disorder caused by the mutations in survival motor neuron 1 gene (SMN1). The molecular pathology of missense mutations in SMN1 is not thoroughly investigated so far. Therefore, we collected all missense mutations in the SMN1 protein, using all possible search terms, from three databases (PubMed, PMC and Google Scholar). All missense mutations were subjected to in silico pathogenicity, conservation, and stability analysis tools. We used statistical analysis as a QC measure for validating the specificity and accuracy of these tools. PolyPhen-2 demonstrated the highest specificity and accuracy. While PolyPhen-1 showed the highest sensitivity; overall, PolyPhen2 showed better measures in comparison to other in silico tools. Three mutations (D44V, Y272C, and Y277C) were identified as the most pathogenic and destabilizing. Further, we compared the physiochemical properties of the native and the mutant amino acids and observed loss of H-bonds and aromatic stacking upon the cysteine to tyrosine substitution, which led to the loss of aromatic rings and may reduce protein stability. The three mutations were further subjected to Molecular Dynamics Simulation (MDS) analysis using GROMACS to understand the structural changes. The Y272C and Y277C mutants exhibited maximum deviation pattern from the native protein as compared to D44V mutant. Further MDS analysis predicted changes in the stability that may have been contributed due to the loss of hydrogen bonds as observed in intramolecular hydrogen bond analysis and physiochemical analysis. A loss of function/structural impact was found to be severe in the case of Y272C and Y277C mutants in comparison to D44V mutation. Correlating the results from in silico predictions, physiochemical analysis, and MDS, we were able to observe a loss of stability in all the three mutants. This combinatorial approach could serve as a platform for variant interpretation and drug design for spinal muscular dystrophy resulting from missense mutations.

Identification of key candidate genes and miRNA‑mRNA target pairs in chronic lymphocytic leukemia by integrated bioinformatics analysis

Chronic lymphocytic leukemia (CLL) is a malignant clonal proliferative disorder of B cells. Inhibition of cell apoptosis and cell cycle arrest are the main pathological causes of this disease, but its molecular mechanism requires further investigation. The purpose of the present study was to identify biomarkers for the early diagnosis and treatment of CLL, and to explore the molecular mechanisms of CLL progression. A total of 488 differentially expressed genes (DEGs) and 32 differentially expressed microRNAs (miRNAs; DEMs) for CLL were identified by analyzing the gene chips GSE22529, GSE39411 and GSE62137. Functional and pathway enrichment analyses of DEGs demonstrated that DEGs were mainly involved in transcriptional dysregulation and multiple signaling pathways, such as the nuclear factor‑κB and mitogen‑activated protein kinase signaling pathways. In addition, Cytoscape software was used to visualize the protein‑protein interactions of these DEGs in order to identify hub genes, which could be used as biomarkers for the early diagnosis and treatment of CLL. Cytoscape software was also used to analyze the association between the predicted target mRNAs of DEMs and DEGs and increase knowledge about the miRNA‑mRNA regulatory network associated with the progression of CLL. Taken together, the present study provided a bioinformatics basis for advancing our understanding of the pathogenesis of CLL by identifying differentially expressed hub genes, miRNA‑mRNA target pairs and molecular pathways. In addition, hub genes may be used as novel biomarkers for the diagnosis of CLL and to guide the selection of CLL drug combinations.

The variation of mitochondrial NADH dehydrogenase subunit 4 (mtND4) and molecular dynamics simulation of SNPs among Iranian women with breast cancer

Breast cancer is the second cause of death among women all around the world. One out of every eight women is diagnosed with breast cancer in Iran. There are many reasons for cancer, one of which is the mutations in the mitochondrial genome observed in most breast cancer studies. However, the aim of this study is to evaluate the genetic region of NADH dehydrogenase subunit 4 in patients with breast cancer. First, the genomic DNA was extracted from a tissue. The NADH dehydrogenase subunit 4 coding region was amplified by PCR, and then the SSCP was sequenced. After that, the molecular dynamics were employed. The association between the mutations and the prognostic factors such as ER, PR, HER-2, and age were statistically examined. The sequence of the ND4 area was determined in 24 suspected patients, and 15 nucleotide changes were reported. The role of this variations was investigated by in-silico. The harmful mutations were predicted based on some servers. The molecular dynamics results showed that there is a significant relationship between the mutant protein and the changes in the structural conformation. Our results showed that the mutation in the ND4 area plays an important role in developing breast cancer. So, it can be concluded that the mitochondrial NADH dehydrogenase analysis may help to detect breast cancer in the early stages.

Identification of candidate aberrantly methylated and differentially expressed genes in thyroid cancer

Aberrant methylation of DNA sequences plays a criticle role in finding novel aberrantly methylated genes and pathways in thyroid cancer (THCA). This study aimed to integrate three cohorts profile datasets to find novel aberrantly methylated genes and pathways in THCA. Data of gene expression profiling microarrays (GSE33630 and GSE65144) and gene methylation profiling microarrays (GSE51090) were downloaded from the Gene Expression Omnibus database. Aberrantly methylated and differentially expressed genes were sorted and pathways were analyzed. Functional and enrichment analyses of selected genes were performed using the String database. A protein‐protein interaction network was constructed using the Cytoscape software, and module analysis was performed using Molecular Complex detection. In total, we identified 12 hypomethylation/high‐expression genes and 30 hypermethylation/low‐expression genes at the screening step and, finally, found 6 mostly changed hub genes including PPARGC1A, CREBBP, EP300, CD44, SPP1, and MMP9. Pathway analysis showed that aberrantly methylated differentially expressed genes were mainly associated with the thyroid hormone signaling pathway, AMP‐activated protein kinase (AMPK) signaling pathway, and cell cycle process in THCA. After validation in the Cancer Genome Atlas database, the methylation and expression status of hub genes was significantly altered and the same with our results. Taken together, we identified novel aberrantly methylated genes and pathways in THCA, which could improve our understanding of the cause and underlying molecular events, and these candidate genes could serve as aberrant methylation‐based biomarkers for precise diagnosis and treatment of THCA.

Targeting Outer Membrane Protein Component AdeC for the Discovery of Efflux Pump Inhibitor against AdeABC Efflux Pump of Multidrug Resistant Acinetobacter baumannii

The structure and functioning of multidrug efflux systems provide us with a better understanding of the transport of various antibiotics, thus giving a path for the discovery of effective compounds for combating the multidrug resistance in Acinetobacter baumannii. In the present study, a number of computational techniques have been used to search for an inhibitor for the RND efflux pump, AdeABC, of A. baumannii targeting specifically its outermost component, i.e., AdeC. We have prepared the three-dimensional structure for AdeC using MODELLER v9.16 and identified its active binding site using SiteMap. Using high-throughput virtual screening, we identified compounds from a large library of biogenic compounds on the basis of their effective interaction at the binding site of AdeC. The validation of docking step was performed by plotting ROC curve (enrichment calculations). The docked complexes were further analyzed for their binding free energies by molecular mechanics using Generalized Born model and Solvent Accessibility (MMGBSA). The molecular dynamics simulation was performed for AdeC-ZINC77257599 complex using GROMACS. The present rational drug designing, molecular mechanics and molecular dynamics data provided an inhibitor, i.e, ZINC77257599 [(3R,4Z,6E,8E)-3-hydroxy-2,2,4-trimethyl-10-oxazol-5-yl-deca-4,6,8-trienamide], for the outer membrane protein component (AdeC) of efflux pump AdeABC of A. baumannii.

Structural analysis of missense mutations in galactokinase 1 (GALK1) leading to galactosemia type-2

Galactosemia type 2 is an autosomal recessive disorder characterized by the deficiency of galactokinase (GALK) enzyme due to missense mutations in GALK1 gene, which is associated with various manifestations such as hyper galactosemia and formation of cataracts. GALK enzyme catalyzes the adenosine triphosphate (ATP)-dependent phosphorylation of α-d-galactose to galactose-1-phosphate. We searched 4 different literature databases (Google Scholar, PubMed, PubMed Central, and Science Direct) and 3 gene-variant databases (Online Mendelian Inheritance in Man, Human Gene Mutation Database, and UniProt) to collect all the reported missense mutations associated with GALK deficiency. Our search strategy yielded 32 missense mutations. We used several computational tools (pathogenicity and stability, biophysical characterization, and physiochemical analyses) to prioritize the most significant mutations for further analyses. On the basis of the pathogenicity and stability predictions, 3 mutations (P28T, A198V, and L139P) were chosen to be tested further for physicochemical characterization, molecular docking, and simulation analyses. Molecular docking analysis revealed a decrease in interaction between the protein and ATP in all the 3 mutations, and molecular dynamic simulations of 50 ns showed a loss of stability and compactness in the mutant proteins. As the next step, comparative physicochemical changes of the native and the mutant proteins were carried out using essential dynamics. Overall, P28T and A198V were predicted to alter the structure and function of GALK protein when compared to the mutant L139P. This study demonstrates the power of computational analysis in variant classification and interpretation and provides a platform for developing targeted therapeutics.

MicroRNA expression in cervical cancer: Novel diagnostic and prognostic biomarkers

Growing evidence has shown that a large number of miRNAs are abnormally expressed in cervical cancer (CC) tissues and play irreplaceable roles in tumorigenesis, progression, and metastasis. This study aimed to identify new biomarkers and pivotal genes associated with CC prognosis through comprehensive bioinformatics analysis. At first, the data of gene expression microarray (GSE30656) was downloaded from GEO database and differential miRNAs were obtained. Additionally, 4 miRNAs associated with the survival time of patients with CC were screened through TCGA differential data analysis, Kaplan-Meier, and Landmark analysis. Among them, the low expression of miR-188 and high expression of miR-223 correlated with the short survival of CC patients, while the down-regulation of miR-99a and miR-125b was closely related to the 5-year survival rate of patients. Then, based on the correspondence between the differentially expressed genes (DEGs) in CC from the TCGA data and the 4 miRNAs target genes, 58 target genes were screened to perform the analysis of function enrichment and the visualization of protein-protein interaction (PPI) networks. The seven pivotal genes of the PPI network as the target genes of four miRNAs related to prognosis, they were directly or indirectly involved in the development of CC. In this study, based on high-throughput data mining, differentially expressed miRNAs and related target genes were analyzed to provide an effective bioinformatics basis for further understanding of the pathogenesis and prognosis of CC. And the results may be a promising biomarker for the early screening of high-risk populations and early diagnosis of cervical cancer.

Nearly neutral evolution in IFNL3 gene retains the immune function to detect and clear the viral infection in HCV

IFNL3 gene plays a crucial role in immune defense against viruses. It induces the interferon stimulated genes (ISGs) with antiviral properties by activating the JAK-STAT pathway. In this study, we investigated the evolutionary force involved in shaping the IFNL3 gene to perform its downstream function as a regulatory gene in HCV clearance. We have selected 25 IFNL3 coding sequences with human gene as a reference sequence and constructed a phylogeny. Furthermore, rate of variation, substitution saturation test, phylogenetic informativeness and differential selection were also analysed. The codon evolution result suggests that nearly neutral mutation is the key pattern in shaping the IFNL3 evolution. The results were validated by subjecting the human IFNL3 protein variants to that of the native through a molecular dynamics simulation study. The molecular dynamics simulation clearly depicts the negative impact on the reported variants in human IFNL3 protein. However, these detrimental mutations (R157Q and R157W) were shown to be negatively selected in the evolutionary study of the mammals. Hence, the variation revealed a mild impact on the IFNL3 function and may be removed from the population through negative selection due to its high functional constraints. In a nutshell, our study may contribute the overall evidence in phylotyping and structural transformation that takes place in the non-synonymous substitutions of IFNL3 protein. Substantially, our obtained theoretical knowledge will lay the path to extend the experimental validation in HCV clearance.

Novel Insights for Inhibiting Mutant Heterodimer IDH1wt-R132H in Cancer: An In-Silico Approach

BACKGROUND

Isocitrate dehydrogenase 1 (IDH1) is a dimeric enzyme responsible for supplying the cell's nicotinamide adenine dinucleotide phosphate (NADPH) reserves via dehydrogenation of isocitrate (ICT) and reduction of NADP+. Mutations in position R132 trigger cancer by enabling IDH1 to produce D-2-hydroxyglutarate (2-HG) and reduce inhibition by ICT. Mutant IDH1 can be found as a homodimer or a heterodimer.

OBJECTIVE

We propose a novel strategy to inhibit IDH1 R132 variants as a means not to decrease the concentration of 2-HG but to provoke a cytotoxic effect, as the cell malignancy at this point no longer depends on 2-HG. We aim to inhibit the activity of the mutant heterodimer to block the wild-type subunit. Limiting the NADPH reserves in a cancerous cell will enhance its susceptibility to the oxidative stress provoked by chemotherapy.

METHODS

We performed a virtual screening using all US FDA-approved drugs to replicate the loss of inhibition of mutant IDH1 by ICT. We characterized our results based on molecular interactions and correlated them with the described phenotypes.

RESULTS

We replicated the loss of inhibition by ICT in mutant IDH1. We identified 20 drugs with the potential to inhibit the heterodimeric isoform. Six of them are used in cancer treatment.

CONCLUSIONS

We present 20 FDA-approved drugs with the potential to inhibit IDH1 wild-type activity in mutated cells. We believe this work may provide important insights into current and new approaches to dealing with IDH1 mutations. In addition, it may be used as a basis for additional studies centered on drugs presenting differential sensitivities to different IDH1 isoforms.

A molecular dynamics approach to explore the structural characterization of cataract causing mutation R58H on human γD crystallin

The crystallins are a family of monomeric proteins present in the mammalian lens and mutations in these proteins cause various forms of cataracts. The aim of our current study is to emphasize the structural characterization of aggregation propensity of mutation R58H on γD crystallin using molecular dynamics (MD) approach. MD result revealed that difference in the sequence level display a wide variation in the backbone atomic position, and thus exhibits rigid conformational dynamics. Changes in the flexibility of residues favoured to increase the number of intra-molecular hydrogen bonds in mutant R58H. Moreover, notable changes in the hydrogen bonding interaction resulted to cause the misfolding of mutant R58H by introducing α-helix. Principal component analysis (PCA) result suggested that mutant R58H showed unusual conformational dynamics along the two principal components when compared to the wild-type (WT)-γD crystallin. In a nutshell, the increased surface hydrophobicity could be the cause of self-aggregation of mutant R58H leading to aculeiform cataract.

Structure-based design of selective phosphodiesterase 4B inhibitors based on ginger phenolic compounds

Phosphodiesterase 4 (PDE4) has been established as a drug target for inflammatory diseases of respiratory tract like asthma and chronic obstructive pulmonary disease. The selective inhibitors of PDE4B, a subtype of PDE4, are devoid of adverse effects like nausea and vomiting commonly associated with non-selective PDE4B inhibitors. This makes the development of PDE4B subtype selective inhibitors a desirable research goal. Thus, in the present study, molecular docking, molecular dynamic simulations and binding free energy were performed to explore potential selective PDE4B inhibitors based on ginger phenolic compounds. The results of docking studies indicate that some of the ginger phenolic compounds demonstrate higher selective PDE4B inhibition than existing selective PDE4B inhibitors. Additionally, 6-gingerol showed the highest PDE4B inhibitory activity as well as selectivity. The comparison of binding mode of PDE4B/6-gingerol and PDE4D/6-gingerol complexes revealed that 6-gingerol formed additional hydrogen bond and hydrophobic interactions with active site and control region 3 (CR3) residues in PDE4B, which were primarily responsible for its PDE4B selectivity. The results of binding free energy demonstrated that electrostatic energy is the primary factor in elucidating the mechanism of PDE4B inhibition by 6-gingerol. Dynamic cross-correlation studies also supported the results of docking and molecular dynamics simulation. Finally, a small library of molecules were designed based on the identified structural features, majority of designed molecules showed higher PDE4B selectivity than 6-gingerol. These results provide important structural features for designing new selective PDE4B inhibitors as anti-inflammatory drugs and promising candidates for synthesis and pre-clinical pharmacological investigations.

Comparative study of the efficacy of pulsed electromagnetic field and low level laser therapy on mitogen-activated protein kinases

Mitogen-Activated Protein Kinases (MAPKs) consist of three major signaling members: extracellular signal-regulated kinase (ERK), p38 and C-JUN N-terminal kinase (JNK). We investigated physiological effects of Pulsed Electromagnetic Field Therapy (PEMFT) and Low Level Laser Therapy (LLLT) on human body, adopting the expression level of mitogen-activated protein kinases as an indicator via assessment of the activation levels of three major families of MAPKS, ERK, p38 and JNK in the peripheral lymphocytes of patients before and after the therapies. Assessment for the expression levels of MAPKs families' were done, in the peripheral lymphocytes of patients recently have appendectomy, using flow cytometric analysis of multiple signaling pathways, pre and post LLLT and PEMFT application (twice daily for 6 successive days) on the appendectomy wound. There were non-significant differences in the expression levels of MAPKs families' pre- therapies application. But there were significant increase in the ERK expression levels post application of LLLT compared to its pre application (p<0.01). Also, there was significant increase in the ERK, p38 and C-Jun N terminal expression level values post application of PEMFT compared to its pre application expression levels (p<0.01 for each). The present study demonstrates that PEMFT has a powerful healing effect more than LLLT as it increase the activation of ERK, P38 and C-Jun-N Terminal while LLLT only increase the activation of ERK. LLLT has more potent pain decreasing effect than PEMFT as it does not activate P38 pathway like PEMFT.

•

PEMFT stimulates all pathways of MAPKs including ERK, P38 and C-Jun-terminal.

•

LLLT stimulates only the ERK pathway in MAPKS activation pathways.

•

PEMFT has a powerful healing effect more than LLLT as it increase the activation of ERK, P38 and C-Jun-N Terminal.

•

LLLT has more potent pain decreasing effect than PEMFT as it does not activate p38 pathway.

•

PEMFT should not be used as pain killer modality in acute pain but it's very useful to be used in chronic pain management with hypothalamic –ve Ingram due to its stimulatory effect on P38.

Clinical significance and prospective molecular mechanism of MALAT1 in pancreatic cancer exploration: a comprehensive study based on the GeneChip, GEO, Oncomine, and TCGA databases

Purpose

Long noncoding RNAs (lncRNAs) are known to function as regulators in the development and occurrence of various tumors. MALAT1 is a highly conserved lncRNA and has vital functions in diverse tumors, including pancreatic cancer (PC). However, the underlying molecular regulatory mechanism involved in the occurrence and development of PC remains largely unknown. Thus, it is important to explore MALAT1 in PC and elucidate its function, which might offer a new perspective for clinical diagnosis and therapy.

Methods

First, we used the Gene Expression Omnibus, Oncomine, and The Cancer Genome Atlas databases to determine the clinical diagnostic and prognostic values of MALAT1. We next used our own GeneChip and The Cancer Genome Atlas database to collect the possible target genes of MALAT1 and further utilized a bioinformatics analysis to explore the underlying significant pathways that might be crucial in PC. Finally, we identified several key target genes of MALAT1 and hope to offer references for future research.

Results

We found that the expression of MALAT1 was significantly elevated in patients with PC. A receiver operating characteristics curve analysis showed a moderate diagnostic value (area under the curve =0.75, sensitivity =0.66, specificity =0.72). A total of 224 important overlapping genes were collected, and six hub genes (CCND1, MAPK8, VEGFA, FOS, CDH1, and HSP90AA1) were identified, of which CCND1, MAPK8, and VEGFA, are important genes in PC. Several pathways, including the mTOR signaling pathway, pathways in cancer, and the MAPK signaling pathway, were suggested to be the vital MALAT1 pathways in PC.

Conclusion

MALAT1 is suggested to be a promising diagnostic biomarker in PC. Six hub genes (CCND1, MAPK8, VEGFA, FOS, CDH1, and HSP90AA1), and specifically CCND1, MAPK8, and VEGFA, might be key MALAT1 target genes in PC. Due to their possible clinical significance in PC, several pathways, such as the mTOR signaling pathway, pathways in cancer, and the MAPK signaling pathway, are worthy of further study.

Non-active site mutations disturb the loop dynamics, dimerization, viral budding and egress of VP40 of the Ebola virus

The first account of the dynamic features of the loop region of VP40 of the Ebola virus (EboV) using accelerated molecular dynamics (aMD) simulations is reported herein. Due to its major role in the Ebola life cycle, VP40 is considered a promising therapeutic target. The available experimental data on the N-terminal domain (NTD) loop indicates that mutations K127A, T129A and N130A demonstrate an unrecognized role for NTD-plasma membrane (PM) interaction for efficient VP40-PM localization, oligomerization, matrix assembly and egress. Despite experimental results, the molecular description of VP40 and the information it can provide still remain vague. Therefore, to gain further molecular insight into the effect of mutations on the loop region of VP40 and its effects on the overall protein conformation and VP40 dimerization, aMD simulations and post-dynamic analyses were employed for wildtype (WT) and mutant systems. The results showed significant variations in the presence of mutations as per RMSF, RMSD, R_g, PCA and distance calculations in comparison to the WT. These results could provide researchers with insight with regards to the conformational aspects concerning VP40 and its close relation to the experimental data. We believe that the results presented in this study will ultimately provide a useful understanding of the structural landscape of the loop region of VP40, which would contribute towards the discovery of novel EboV inhibitors.

Simulation Based Investigation of Deleterious nsSNPs in ATXN2 Gene and Its Structural Consequence Toward Spinocerebellar Ataxia

Spinocerebellar degeneration, termed as ataxia is a neurological disorder of central nervous system, characterized by limb in-coordination and a progressive gait. The patient also demonstrates specific symptoms of muscle weakness, slurring of speech, and decreased vibration senses. Expansion of polyglutamine trinucleotide (CAG) within ATXN2 gene with 35 or more repeats, results in spinocerebellar ataxia type-2. Protein ataxin-2 coded by ATXN2 gene has been reported to have a crucial role in translation of the genetic information through sequestering the histone acetyl transferases (HAT) resulting in a state of hypo-acetylation. In the present study, we have evaluated the outcome for 122 non synonymous single nucleotide polymorphisms (nsSNPs) reported within ATXN2 gene through computational tools such as SIFT, PolyPhen 2.0, PANTHER, I-mutant 2.0, Phd-SNP, Pmut, MutPred. The apo and mutant (L305V and Q339L) form of structures for the ataxin-2 protein were modeled for gaining insights toward 3D spatial arrangement. Further, molecular dynamics simulations and structural analysis were performed to observe the brunt of disease associated nsSNPs toward the strength and secondary properties of ataxin-2 protein structure. Our results showed that, L305V is a highly deleterious and disease causing point substitution. Analysis based on RMSD, RMSF, Rg, SASA, number of hydrogen bonds (NH bonds), covariance matrix trace, projection analysis for eigen vector demonstrated a significant instability and conformation along with rise in mutant flexibility values in comparison to the apo form of ataxin-2 protein. The study provides a blue print of computational methodologies to examine the ataxin-blend SNPs. J. Cell. Biochem. 119: 499-510, 2018. © 2017 Wiley Periodicals, Inc.

R, C. S. Expediting dynamics approach to understand the influence of 14-3-3ζ causing metastatic cancer through the interaction of YAP1 and β-TRCP

The 14-3-3ζ protein acts as a molecular switch in regulating the TGF-β pathway, which alters from a tumor suppressor in the early stage of breast cancer to a promoter of metastasis in the late stage.

Influence of MDM2 polymorphisms on squamous cell carcinoma susceptibility: a meta-analysis

PURPOSE

Controversial associations between single-nucleotide polymorphisms (rs2279744, rs937283, rs3730485) of the MDM2 gene and the etiology of squamous cell carcinomas (SCCs) have been reported. This merits further comprehensive assessment.

MATERIALS AND METHODS

We systematically reviewed the available data and conducted an updated meta-analysis to evaluate the genetic effect of MDM2 polymorphisms in SCC susceptibility, using Stata/SE 12.0 software.

RESULTS

After screening, 7,987 SCC cases and 12,954 controls from 26 eligible case-control studies were enrolled. Overall, compared with the control group, a significantly increased SCC risk was observed for the MDM2 rs2279744 polymorphism in the Asian population (test of association: odds ratio [OR] 1.12, P=0.027 for G vs T; OR 1.26, P=0.016 for GG vs TT; OR 1.25, P<0.001 for GG vs TT + TG; and OR 1.08, P=0.023 for carrier G vs T). In subgroup analysis by SCC type, a similarly increased esophageal SCC risk was detected (OR 1.19, P<0.001 for G vs T; OR 1.46, P<0.001 for GG vs TT; and OR 1.48, P=0.005 for GG vs TT + TG). Furthermore, MDM2-TP53 double mutation was statistically associated with increased SCC susceptibility overall (OR 1.52, P=0.001), especially in the Asian population (OR 1.49, P=0.022). However, no significant difference between the control and case groups was obtained for MDM2 rs937283 or rs3730485 under any genetic model (all P>0.05).

CONCLUSION

Our results highlight a positive association between the GG genotype of MDM2 rs2279744 polymorphism and an increased risk of esophageal SCC in the Asian population, which needs to be clarified by more large-scale studies.

In Silico Template Selection of Short Antimicrobial Peptide Viscotoxin for Improving Its Antimicrobial Efficiency in Development of Potential Therapeutic Drugs

Rapid increase in antibiotic resistance has posed a worldwide threat, due to increased mortality, morbidity, and expenditure caused by antibiotic-resistant microbes. Recent development of the antimicrobial peptides like viscotoxin (Vt) has been successfully comprehended as a substitute for classical antibiotics. A structurally stable peptide, Vt can enhance antimicrobial property and can be used for various developmental purposes. Thus, structural stability among the antimicrobial peptides, Vt A1 (3C8P), A2 (1JMN), A3 (1ED0), B (1JMP), and C (1ORL) of Viscus album was computationally analyzed. In specific, the static confirmation of VtA3 showed high number of intramolecular interactions, along with an increase in hydrophobicity than others comparatively. Further, conformational sampling was used to analyze various geometrical parameters such as root mean square deviation, root mean square fluctuation, radius of gyration, and ovality which also revealed the structural stability of VtA3. Moreover, the statistically validated contours of surface area, lipophilicity, and distance constraints of disulfide bonds also supported the priority of VtA3 with respect to stability. Finally, the functional activity of peptides was accessed by computing their free energy of membrane association and membrane interactions, which defined VtA3 as functionally stable. Currently, peptide-based antibiotics and nanoparticles have attracted the pharmaceutical industries for their potential therapeutic applications. Thereby, it is proposed that viscotoxin A3 (1ED0) could be used as a preeminent template for scaffolding potentially efficient antimicrobial peptide-based drugs and nanomaterials in future.

Candidate SNP Markers of Chronopathologies Are Predicted by a Significant Change in the Affinity of TATA-Binding Protein for Human Gene Promoters

Variations in human genome (e.g., single nucleotide polymorphisms, SNPs) may be associated with hereditary diseases, their complications, comorbidities, and drug responses. Using Web service SNP_TATA_Comparator presented in our previous paper, here we analyzed immediate surroundings of known SNP markers of diseases and identified several candidate SNP markers that can significantly change the affinity of TATA-binding protein for human gene promoters, with circadian consequences. For example, rs572527200 may be related to asthma, where symptoms are circadian (worse at night), and rs367732974 may be associated with heart attacks that are characterized by a circadian preference (early morning). By the same method, we analyzed the 90 bp proximal promoter region of each protein-coding transcript of each human gene of the circadian clock core. This analysis yielded 53 candidate SNP markers, such as rs181985043 (susceptibility to acute Q fever in male patients), rs192518038 (higher risk of a heart attack in patients with diabetes), and rs374778785 (emphysema and lung cancer in smokers). If they are properly validated according to clinical standards, these candidate SNP markers may turn out to be useful for physicians (to select optimal treatment for each patient) and for the general population (to choose a lifestyle preventing possible circadian complications of diseases).