Forecasting most deleterious nsSNPs in human TLR9 gene and their cumulative impact on biophysical features of the protein using in silico approaches.

In women, the uterine cervix and corpus uteri are two main suspects, playing a major role in cancer-associated-mortality. Immunologically, Toll-like receptors (TLRs) associated with the innate immune system, can recognize pathogens and induce immune responses against pathogens. Cellularly, TLR9 expression occurs in immune system cells including macrophages, natural killer cells, dendritic cells, and other antigen-presenting cells. TLR9 recognizes and interacts with viral and bacterial DNA comprising cytosine-phosphate-guanine (CpG) dideoxynucleotide motif. The current study is designed to identify the most deleterious nonsynonymous single nucleotide polymorphisms (nsSNPs) in the TLR9 gene and to delineate their deleterious effect on the structural and functional features of proteins at the molecular level. Based on the implementation of various computational tools and algorithms eight most deleterious nsSNPs (P139H, R257C, C265Y, L283P, G514D, L544Q, H566Y, and W670R) have been identified in the human TLR9 gene as potentially damaging SNPs. Further, our study suggests highly conserved patterns at deleterious nsSNPs sites could influence protein stability and its functional features. Additionally, this study identifies two nsSNPs (G514D and W670R) associated with the severity of Uterine corpus endometrial carcinoma. In support of our computational findings, the validation of key results using polymerase chain reaction and other experimental methods is warranted in the Indian population. In general, this study might be able to delineate the guideline for identifying the most damaging SNPs and enhances the understating of the risk factors for cancer and disease susceptibilities.

Structure-based study to identify alkaloids as promising cytochrome P450 (CYP1A1) inhibitors: An in silico approach using virtual screening, molecular dynamic simulations, and binding free energy calculation.

Carcinogens present in smokeless tobacco (SLT) like tobacco-specific nitrosamines can be metabolized by the cytochrome P450 (CYP450) enzyme. Functionally, the CYP450 enzyme resides in a heme pigment to perform the catalytic activity. The CYP1A1 is one of the main extrahepatic CYP450 enzymes known to detoxify toxic substances and activate carcinogens. The CYP1A1 inhibition by potential inhibitors reduce the chance of oral cancer. The current study aimed to explore more about the inhibitor binding site and identification of lead alkaloids, that could work as putative inhibitors against target CYP1A1. In respect, we have performed docking studies, virtual screening of alkaloids, and natural product libraries against CYP1A1 followed by molecular dynamic simulations and binding free energy calculations. Docking studies of tobacco-specific nitrosamine (TSNA) products and their similar carcinogen analogs revealed that the heme group is bound to the floor of the bowl-shaped cavity whereas carcinogens are bound to the roof of the rounded shape cavity. Furthermore, virtual screening and binding free energy calculations revealed Tomatidine as a putative inhibitor against CYP1A1. On the basis of altogether outcomes of the current study, we have concluded that the addition of lead-hit alkaloid Tomatidine and others in SLT products may be working as a supplement that could be able to reduce the expression of human CYP1A1 and suppresses carcinogenic by-products formations.

In silico prediction and interaction of resveratrol on methyl-CpG binding proteins by molecular docking and MD simulations study.

Resveratrol enhances the BRCA1 gene expression and the MBD family of proteins bind to the promoter region of the BRCA1 gene. However, the molecular interaction is not yet reported. Here we have analyzed the binding affinity of resveratrol with MBD proteins. Our results suggest that resveratrol binds to the MBD proteins with higher binding affinity toward MeCP2 protein (ΔG = -6.5) by sharing four hydrogen bonds as predicted by molecular docking studies. Further, the molecular dynamics simulations outcomes showed that the backbones of all three protein-ligand complexes are stabilized after the period of 75 ns, constantly fluctuating around the deviations of 0.4 Å, 0.5 Å and 0.7 Å for MBD1, MBD2 and MeCP2, respectively. The inter-molecular hydrogen bonding trajectory analysis for protein-ligand complexes also support the strong binding between MeCP2-resveratrol complex. Further, binding free energy calculations showed binding energy of -94.764 kJ mol-1, -53.826 kJ mol-1 and -36.735 kJ mol-1 for MeCP2-resveratrol, MBD2-resveratrol and MBD1-resveratrol complexes, respectively, which also supported our docking results. Our study also highlighted that the MBD family of proteins forms a binding interaction with other signaling proteins that are involved in various cancer initiation pathways.

Functional divergence and comparative in-silico study of Cas4 proteins of DUF83 class.

Clustered Regularly Interspaced Short Palindromic Repeats-CRISPR associated (CRISPR-Cas) systems present in genomes of bacteria and archaea have been the focus of many research studies recently. The Cas4 proteins of these systems are thought to be responsible for the adaptation step in the CRISPR mechanism. Cas4 proteins exhibit low sequence similarity among themselves and are currently classified into 2 main classes: DUF83 and DUF911. The characteristic features of Cas4 proteins belonging to DUF83 class have been elucidated by determining the structures of Cas4 protein from Sulfolobus solfataricus and Pyrobaculum calidifontis. Although, both Cas4 proteins characterized structurally are of same DUF83 class, these 2 proteins do exhibit significant biochemical and functional differences. The aim of the present study was to explore the structural and evolutionary features responsible for these differences. Our study predicts residues which might be responsible for such differences. Functional divergence analysis was used to predict sites exhibiting type I divergence, where certain amino acids are conserved in 1 clade whereas the same site is highly variable in the other clade. Our intra-molecular interaction analysis reinforces the influence of such divergence sites on the other functionally important amino acids. In general, this study identifies some of the divergence hotspots that could be the focus of future experimental studies for better understanding of Cas4 enzymatic activity in CRISPR mechanism.

Targeted mutations and MD simulations of a methanol-stable lipase YLIP9 from Yarrowia lipolytica MSR80 to develop a biodiesel enzyme.

Biodiesel, an environment friendly alternative for fuels, contains methyl esters of long-chain fatty acids. Our group has reported a methanol-stable YLIP9 from Yarrowia lipolytica MSR80 that shows poor catalysis of long-chain fatty acids. To shift its substrate specificity, residues within lid and binding pocket were identified for sequential mutations using YLIP2 as the template. Of the two point mutations (Glu116Leu and Ser119Val) introduced in the lid, the former mutation (YLIP9L1) increased the catalytic rate by ∼2-fold without any change in substrate specificity. In this mutant, six binding pocket residues (Bp2-Bp7) were further mutated to obtain six double mutants. YLIP9L1Bp3 showed significant shift in substrate specificity towards long-chain pNPesters with 11-fold increase in catalytic efficiency than YLIP9. Double mutations also led to increased thermostability and lowered activation energy of YLIP9L1Bp3 thereby shifting its optimum temperature from 60°C to 50°C. In silico molecular dynamics simulations revealed improved lid flexibility and increased catalytic triad volume in YLIP9L1Bp3. The enzyme YLIP9L1Bp3 was methanol-stable having selectivity for long-chain fatty acids with improved catalytic efficiency. Its application as a biodiesel enzyme was validated by transesterification of palm oil in presence of methanol, where it showed 8-fold increase in conversion of oil to methyl esters.

Dual targeting of MDM2 with a novel small-molecule inhibitor overcomes TRAIL resistance in cancer.

Mouse double minute 2 (MDM2) protein functionally inactivates the tumor suppressor p53 in human cancer. Conventional MDM2 inhibitors provide limited clinical application as they interfere only with the MDM2-p53 interaction to release p53 from MDM2 sequestration but do not prevent activated p53 from transcriptionally inducing MDM2 expression. Here, we report a rationally synthesized chalcone-based pyrido[ b ]indole, CPI-7c, as a unique small-molecule inhibitor of MDM2, which not only inhibited MDM2-p53 interaction but also promoted MDM2 degradation. CPI-7c bound to both RING and N-terminal domains of MDM2 to promote its ubiquitin-mediated degradation and p53 stabilization. CPI-7c-induced p53 directly recruited to the promoters of DR4 and DR5 genes and enhanced their expression, resulting in sensitization of TNF-related apoptosis-inducing ligand (TRAIL)-resistant cancer cells toward TRAIL-induced apoptosis. Collectively, we identified CPI-7c as a novel small-molecule inhibitor of MDM2 with a unique two-prong mechanism of action that sensitized TRAIL-resistant cancer cells to apoptosis by modulating the MDM2-p53-DR4/DR5 pathway.

Selective disruption of disulphide bonds lowered activation energy and improved catalytic efficiency in TALipB from Trichosporon asahii MSR54: MD simulations revealed flexible lid and extended substrate binding area in the mutant.

TALipB (33 kDa) is a solvent stable, enantioselective lipase from Trichosporon asahii MSR54. It is cysteine-rich and shows activation in presence of thiol reducing agents. DIANNA server predicted three disulphide bridges C53-C195 (S1), C89-C228 (S2) and C164-C254 (S3) in the enzyme. Selective disruption of disulphide bonds by cysteine to alanine mutations at Cys53 and Cys89 of S1 and S2 bonds resulted in enzyme activation. Mutant mTALipB (S1+S2) showed increase in specific activity by ∼4-fold (834 mM/mg) and improved Vmax of 6.27 µmol/ml/min at 40 °Con pNP caprate. Temperature optima of mTALipB shifted from 50 to 40 °C and activation energy decreased by 0.7 kcal mol(-1). However, the mutant was less thermostable with a t1/2 of 18 min at 60 °C as compared to t1/2 of 38 min for the native enzyme. Mutant also displayed an improved activity on all pNP esters and higher enantiomeric excess (61%) during esterification of (±) 1-phenylethanol. Far-UV CD analysis showed significant changes in secondary structure after S-S bridge disruption with 7.16% decrease in α-helices and 1.31% increase in ß-sheets. In silico analysis predicted two lids (α5 and α9) in TALipB. Molecular dynamic simulations at 40 °C and 50 °C revealed that in the mTALipB, both the lids opened at 40 °C with clockwise and anticlockwise rotations in Lid1 and Lid2, respectively. In the native protein, however, the lid was only partially open even at 50 °C. Concomitant to lid flexibility, there was an extension of accessible catalytic triad surface area resulting in improved catalytic efficiency of the mutant enzyme.

"Phylogenetic and evolutionary analysis of functional divergence among Gamma glutamyl transpeptidase (GGT) subfamilies".

BACKGROUND: γ-glutamyltranspeptidase (GGT) is a bi-substrate enzyme conserved in all three domains of life. It catalyzes the cleavage and transfer of γ-glutamyl moiety of glutathione to either water (hydrolysis) or substrates like peptides (transpeptidation). GGTs exhibit great variability in their enzyme kinetics although the mechanism of catalysis is conserved. Recently, GGT has been shown to be a virulence factor in microbes like Helicobacter pylori and Bacillus anthracis. In mammalian cells also, GGT inhibition prior to chemotherapy has been shown to sensitize tumors to the therapy. Therefore, lately both bacterial and eukaryotic GGTs have emerged as potential drug targets, but the efforts directed towards finding suitable inhibitors have not yielded any significant results yet. We propose that delineating the residues responsible for the functional diversity associated with these proteins could help in design of species/clade specific inhibitors. RESULTS: In the present study, we have carried out phylogenetic analysis on a set of 47 GGT-like proteins to address the functional diversity. These proteins segregate into various subfamilies, forming separate clades on the tree. Sequence conservation and motif prediction studies show that even though most of the highly conserved residues have been characterized biochemically in previous studies, a significant number of novel putative sites and motifs are discovered that vary in a clade specific manner. Many of the putative sites predicted during the functional divergence type I and type II analysis, lie close to the known catalytic residues and line the walls of the substrate binding cavity, reinforcing their role in modulating the substrate specificity, catalytic rates and stability of this protein. CONCLUSION: The study offers interesting insights into the evolution of GGT-like proteins in pathogenic vs. non-pathogenic bacteria, archaea and eukaryotes. Our analysis delineates residues that are highly specific to each GGT subfamily. We propose that these sites not only explain the differences in stability and catalytic variability of various GGTs but can also aid in design of specific inhibitors against particular GGTs. Thus, apart from the commonly used in-silico inhibitor screening approaches, evolutionary analysis identifying the functional divergence hotspots in GGT proteins could augment the structure based drug design approaches.

A hydrolytic γ-glutamyl transpeptidase from thermo-acidophilic archaeon Picrophilus torridus: binding pocket mutagenesis and transpeptidation.

γ-Glutamyl transpeptidase of a thermo-acidophilic archaeon Picrophilus torridus was cloned and expressed using E. coli Rosetta-pET 51b(+) expression system. The enzyme was expressed at 37 °C/200 rpm with γ-GT production of 1.99 U/mg protein after 3 h of IPTG induction. It was improved nearby 10-fold corresponding to 18.92 U/mg protein in the presence of 2 % hexadecane. The enzyme was purified by Ni(2+)-NTA with a purification fold of 3.6 and recovery of 61 %. It was synthesized as a precursor heterodimeric protein of 47 kDa with two subunits of 30 kDa and 17 kDa, respectively, as revealed by SDS-PAGE and western blot. The enzyme possesses hydrolase activity with optima at pH 7.0 and 55 °C. It was thermostable with a t (1/2) of 1 h at 50 °C and 30 min at 60 °C, and retained 100 % activity at 45 °C even after 24 h. It was inhibited by azaserine and DON and PMSF. Ptγ-GT shared 37 % sequence identity and 53 % homology with an extremophile γ-GT from Thermoplasma acidophilum. Functional residues identified by in silico approaches were further validated by site-directed mutagenesis where Tyr327 mutated by Asn327 introduced significant transpeptidase activity.

Comparative biochemical characterization and in silico analysis of novel lipases Lip11 and Lip12 with Lip2 from Yarrowia lipolytica.

Novel lipases lip11 and lip12 from Yarrowia lipolytica MSR80 were cloned and expressed in E. coli HB101 pEZZ18 system along with lip2. These enzymes were constitutively expressed as extracellular proteins with IgG tag. The enzymes were purified by affinity chromatography and analyzed by SDS-PAGE with specific activity of 314, 352 and 198 U/mg for Lip2, Lip11 and Lip12, respectively on olive oil. Biochemical characterization showed that all were active over broad range of pH 4.0-9.0 and temperature 20-80 °C with optima at pH 7 and 40 °C. All the three lipases were thermostable up to 80 °C with varying t(½). Activity on various substrates revealed that they were most active on oils > triacylglycerides > p-np-esters. Relatively Lip2 and Lip11 showed specificity for mid to long chain fatty acids, while Lip12 was mid chain specific. GC analysis of triolein hydrolysis by these lipases revealed that Lip2 and Lip11 are regioselective, while Lip12 is not. Effect of metal ions showed that Lip2 and Lip12 were activated by Ca²âº whereas Lip11 by Mg²âº. All were thiol activated and inhibited by PMSF and N-bromosuccinimide. All were activated by non polar solvents and inhibited by polar solvents. Detailed sequence analysis and structural predictions revealed Lip11 and Lip12 shared 61 and 62 % homology with Lip2 (3O0D) and three dimensional superimposition revealed Lip2 was closer to Lip11 than to Lip12 as was observed during biochemical characterization. Finally, thermostability and substrate specificity has been explained on the basis of detailed amino acid analysis.