Receptor based virtual screening of potential novel inhibitors of tigar [TP53 (tumour protein 53)-induced glycolysis and apoptosis regulator.

TP53 (tumor protein 53)-induced glycolysis and apoptosis regulator (TIGAR) belongs to the phosphatases family of proteins that modulates the level of reactive oxygen species in tumor cells. This protein plays a vital role as a negative regulator of glycolysis, thus lowering ROS levels in the cells, which helps the cancerous cells to resist programmed cell death. Besides, TIGAR also mediates the DNA damage repair in cancer cells by increasing tumor cell survival. In the current study, we have screened natural products that compete with the substrate to bind to the active site of TIGAR. Extra precision and MMGBSA scoring function were used to screen the lead molecules. Five compounds were considered as lead molecules with 2-(2-(3,4-dihydroxy phenyl)-3,5-dihydroxy-8-(4-hydroxyphenyl)-4-oxo-4H-furo[2,3-h]chromen-9-yl) acetic acid(DDFA) as a top lead with a docking score of -9.428, and -53.16 MMGBSA, bind to the positively charged amino acids present in the active site. Further, the molecular dynamics simulation studies indicated the structural stability attained by TIGAR protein upon the binding of DDFA, suggesting it to be a potent inhibitor of TIGAR, and could be employed as an anticancer drug during combinational therapy.

Comparison of anticancer activity of biocompatible ZnO nanoparticles prepared by solution combustion synthesis using aqueous leaf extracts of Abutilon indicum, Melia azedarach and Indigofera tinctoria as biofuels.

Recently, there has been an upsurge in the use of naturally available fuels for solution combustion synthesis (SCS) of nanoparticles. Although many reports suggest that these biofuels pose less harm to the environment, their strategic advantages and reliability for making NPs has not been discussed. In the present work, we try to address this issue using plant extracts as biofuels for the SCS of zinc oxide nanoparticles as a model system. In the present work, combustion synthesis of ZnO NPs using lactose and aqueous leaf extracts of Abutilon indicum, Melia azedarach, Indigofera tinctoria as biofuels has been carried out. A comparative analysis of the obtained powders has been conducted to understand the strategic advantages of using plant extracts over a chemical as combustion fuel for the synthesis of zinc oxide nanoparticles. The X-ray diffractograms of the samples revealed the presence of Wurtzite hexagonal structure with varying crystallite sizes. Morphological studies indicated that samples prepared using biofuels had smaller diameter than those prepared using lactose as fuel. Surface characteristics of the samples were measured by X-ray photoelectron spectroscopy. Qualitative phytochemical screening of aqueous leaf extracts revealed the presence of many phytochemicals in them, which might be responsible for combustion. Gas chromatography mass spectrum was carried out to detect the phytochemicals present in the aqueous extracts of the leaves. Further, anticancer evaluation carried out against DU-145 and Calu-6 cancer cells indicated higher anticancer activity of zinc oxide nanoparticles prepared using biofuels. The results of blood haemolysis revealed the biocompatibility of zinc oxide nanoparticles at lower concentrations. In conclusion, we propose that multiple other studies would be required in order to vindicate the potential advantages of using naturally available fuels in SCS.

Sequence Analysis and Phylogenetic Studies of Hypoxia-Inducible Factor-1α.

Hypoxia-inducible factors (HIF) belong to the basic helix loop helix-PER ARNT SIM (bHLH-PAS) family of transcription factors that induce metabolic reprogramming under hypoxic condition. The phylogenetic studies of hypoxia-inducible factor-1α (HIF-1α) sequences across different organisms/species may leave a clue on the evolutionary relationships and its probable correlation to tumorigenesis and adaptation to low oxygen environments. In this study, we have aimed at the evolutionary investigation of the protein HIF-1α across different species to decipher their sequence variations/mutations and look into the probable causes and abnormal behaviour of this molecule under exotic conditions. In total, 16 homologous sequences for HIF-1α were retrieved from the National Center for Biotechnology Information. Sequence identity was performed using the Needle program. Multiple aligned sequences were used to construct the phylogeny using the neighbour-joining method. Most of the changes were observed in oxygen-dependent degradation domain and inhibitory domain. Sixteen sequences were clustered into 5 groups. The phylogenetic analysis clearly highlighted the variations that were observed at the sequence level. Comparisons of the HIF-1α sequence among cancer-prone and cancer-resistant animals enable us to find out the probable clues towards potential risk factors in the development of cancer.

Experimental conformational energy maps of proteins and peptides.

We have presented an extensive analysis of the peptide backbone dihedral angles in the PDB structures and computed experimental Ramachandran plots for their distributions seen under a various constraints on X-ray resolution, representativeness at different sequence identity percentages, and hydrogen bonding distances. These experimental distributions have been converted into isoenergy contour plots using the approach employed previously by F. M. Pohl. This has led to the identification of energetically favored minima in the Ramachandran (ϕ, ψ) plots in which global minima are predominantly observed either in the right-handed α-helical or the polyproline II regions. Further, we have identified low energy pathways for transitions between various minima in the (ϕ,ψ) plots. We have compared and presented the experimental plots with published theoretical plots obtained from both molecular mechanics and quantum mechanical approaches. In addition, we have developed and employed a root mean square deviation (RMSD) metric for isoenergy contours in various ranges, as a measure (in kcal.mol-1 ) to compare any two plots and determine the extent of correlation and similarity between their isoenergy contours. In general, we observe a greater degree of compatibility with experimental plots for energy maps obtained from molecular mechanics methods compared to most quantum mechanical methods. The experimental energy plots we have investigated could be helpful in refining protein structures obtained from X-ray, NMR, and electron microscopy and in refining force field parameters to enable simulations of peptide and protein structures that have higher degree of consistency with experiments. Proteins 2017; 85:979-1001. © 2017 Wiley Periodicals, Inc.

Characterization of a hypothetical protein YVRE from Bacillus subtilis indicates its key role as glucono-lactonase in pentose phosphate pathway and glucose metabolism.

Hypothetical proteins are functionally uncharacterized proteins with assigned function using sequence annotation tools. Almost half of the coding regions of several genomes are hypothetical proteins. Therefore, it is of our interest to characterize a hypothetical protein YVRE from the model system Bacillus subtilis using known data. YVRE is assigned the function as a glucono-lactonase using prediction and phylogenetic analysis. A molecular dynamics simulated homology model of YVRE (with calcium) using human senescence marker protein 30 /SMP30 (PDB ID: 3G4E) as template is reported for functional inference. It is observed that the protein possesses bivalent metal binding domain. Molecular docking studies with the substrate glucono-δ-lactone show YVRE binding with the substrate. This data was further validated using cloning and sub-cloning in pUC57 and pET22b+ respectively, followed by expression and purification using nickel affinity chromatography. The activity of YVRE using the substrate glucono-δ-lactone was calculated. The results show the function of YVRE as a gluconolactonase, with higher preference to zinc than calcium or magnesium. Thus, YVRE is shown to play key role in three metabolic pathways namely, pentose phosphate pathway, ascorbate and aldarate metabolism, and caprolactam degradation.

Characterization of hypothetical protein VNG0128C from Halobacterium NRC-1 reveals GALE like activity and its involvement in Leloir pathway of galactose metabolism.

VNG0128C, a hypothetical protein from Halobacterium NRC-1, was chosen for detailed insilico and experimental investigations. Computational exercises revealed that VNG0128C functions as NAD(+) binding protein. The phylogenetic analysis with the homolog sequences of VNG0128C suggested that it could act as UDP-galactose 4-epimerase. Hence, the VNG0128C sequence was modeled using a suitable template and docking studies were performed with NAD and UDP-galactose as ligands. The binding interactions strongly indicate that VNG0128C could plausibly act as UDP-galactose 4-epimerase. In order to validate these insilico results, VNG0128C was cloned in pUC57, subcloned in pET22b(+), expressed in BL21 cells and purified using nickel affinity chromatography. An assay using blue dextran was performed to confirm the presence of NAD binding domain. To corroborate the epimerase like enzymatic role of the hypothetical protein, i.e. the ability of the enzyme to convert UDP-galactose to UDP-glucose, the conversion of NAD to NADH was measured. The experimental assay significantly correlated with the insilico predictions, indicating that VNG0128C has a NAD(+) binding domain with epimerase activity. Consequently, its key role in nucleotide-sugar metabolism was thus established. Additionally, the work highlights the need for a methodical characterization of hypothetical proteins (less studied class of biopolymers) to exploit them for relevant applications in the field of biology.

Preferential selectivity of inhibitors with human tau protein kinase gsk3ß elucidates their potential roles for off-target Alzheimer's therapy.

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by the accumulation of amyloid beta peptides (Aß) and neurofibrillary tangles (NFTs). The abnormal phosphorylation of tau leads to the formation of NFTs produced by the action of tau kinases, resulting in the loss of neurons and synapse, leading to dementia. Hence, tau kinases have become potential drug target candidates for small molecule inhibitors. With an aim to explore the identification of a common inhibitor, this investigation was undertaken towards analyzing all 10 tau kinases which are implicated in phosphorylation of AD. A set of 7 inhibitors with varied scaffolds were collected from the Protein Data Bank (PDB). The analysis, involving multiple sequence alignment, 3D structural alignment, catalytic active site overlap, and docking studies, has enabled elucidation of the pharmacophoric patterns for the class of 7 inhibitors. Our results divulge that tau protein kinases share a specific set of conserved structural elements for the binding of inhibitors and ATP, respectively. The scaffold of 3-aminopyrrolidine (inhibitor 6) exhibits high preferential affinity with GSK3ß. Surprisingly, the PDB does not contain the structural details of GSK3ß with this specific inhibitor. Thus, our investigations provide vital clues towards design of novel off-target drugs for Alzheimer's.

Identification of potential drug targets implicated in Parkinson's disease from human genome: insights of using fused domains in hypothetical proteins as probes.

High-throughput genome sequencing has led to data explosion in sequence databanks, with an imbalance of sequence-structure-function relationships, resulting in a substantial fraction of proteins known as hypothetical proteins. Functions of such proteins can be assigned based on the analysis and characterization of the domains that they are made up of. Domains are basic evolutionary units of proteins and most proteins contain multiple domains. A subset of multidomain proteins is fused domains (overlapping domains), wherein sequence overlaps between two or more domains occur. These fused domains are a result of gene fusion events and their implication in diseases is well established. Hence, an attempt has been made in this paper to identify the fused domain containing hypothetical proteins from human genome homologous to parkinsonian targets present in KEGG database. The results of this research identified 18 hypothetical proteins, with domains fused with ubiquitin domains and having homology with targets present in parkinsonian pathway.

Sequence analyses and comparative modeling of fly and worm fibroblast growth factor receptors indicate that the determinants for FGF and heparin binding are retained in evolution.

The presence of a large number of fibroblast growth factors (FGFs) and multiple splice forms of their receptors (FGFRs) in higher vertebrates makes the three-dimensional (3D) analysis of FGF interactions with their receptors a formidable task. The situation differs in Caenorhabditis elegans (worm) and Drosophila melanogaster (fruit fly), where only one or two FGF and FGFR sequences have been identified. Structural studies of the FGF-FGFR complexes in such primitive organisms should reveal the basic features of the ligand-receptor interactions as they first emerged through evolution. We have analysed the sequences of worm and fly FGFs and FGFRs and used the recently determined crystal structure of the human FGF1-FGFR2-heparin ternary complex [Pellegrini, L., Burke, D.F., von Delft, F., Mulloy, B. and Blundell, T.L. (2000) Nature 407, 1029-34] to construct 3D models of the homologous complexes. In spite of a low sequence similarity with their human counterparts, key structural features required for ligand-receptor and protein-heparin binding in humans are conserved in the fly and worm FGF-FGFR-heparin complexes. Analyses of the models show that tertiary interactions that are not conserved in sequence are maintained through novel interactions or complementary mutations in the fly and worm sequences. The overall charge distributions observed in the human FGF-FGFR-heparin complex are retained in the fly and worm models. The arginine residue at position 253 in the linker region between the Ig-like domains D2 and D3 in the wild type fly and worm sequences is particularly striking, as the Pro253Arg mutation in humans is responsible for Apert syndrome. This change may enhance the affinity of receptors for their FGF molecules as observed in Apert mutants.

Sequence-structure homology recognition by iterative alignment refinement and comparative modeling.

Our approach to fold recognition for the fourth critical assessment of techniques for protein structure prediction (CASP4) experiment involved the use of the FUGUE sequence-structure homology recognition program (http://www-cryst.bioc.cam.ac.uk/fugue), followed by model building. We treat models as hypotheses and examine these to determine whether they explain the available data. Our method depends heavily on environment-specific substitution tables derived from our database of structural alignments of homologous proteins (HOMSTRAD, http://www-cryst.bioc.cam.ac.uk/homstrad/). FUGUE uses these tables to incorporate structural information into profiles created from HOMSTRAD alignments that are matched against a profile created for the target from multiple sequence alignment. In addition, environment-specific substitution tables are used throughout the modeling procedure and as part of the model evaluation. Annotation of sequence alignments with JOY, to reflect local structural features, proved valuable, both for modifying hypotheses, and for rejecting predictions when the expected pattern of conservation is not observed. Our stringency in rejecting incorrect predictions led us to submit a relatively small number of models, including only a low number of false positives, resulting in a high average score.

Role of water in plasticity, stability, and action of proteins: the crystal structures of lysozyme at very low levels of hydration.

Earlier studies involving water-mediated transformations in lysozyme and ribonuclease A have shown that the overall movements in the protein molecule consequent to the reduction in the amount of surrounding water are similar to those that occur during enzyme action, thus highlighting the relationship among hydration, plasticity, and action of these enzymes. Monoclinic lysozyme retains its crystallinity even when the level of hydration is reduced further below that necessary for activity (about 0.2 gram of water per gram of protein). In order to gain insights into the role of water in the stability and the plasticity of the protein molecule and the geometrical basis for the loss of activity that accompanies dehydration, the crystal structures of monoclinic lysozyme with solvent contents of 17.6%, 16.9%, and 9.4% were determined and refined. A detailed comparison of these forms with the normally hydrated forms show that the C-terminal segment (residues 88-129) of domain I and the main loop (residues 65-73) in domain II exhibit large deviations in atomic positions when the solvent content is reduced, although the three-dimensional structure is essentially preserved. Many crucial water bridges between different regions of the molecule are conserved in spite of differences in detail, even when the level of hydration is reduced well below that required for activity. The loss of activity that accompany dehydration appears to be caused by the removal of functionally important water molecules from the active-site region and the reduction in the size of the substrate binding cleft.

Plasticity, hydration and accessibility in ribonuclease A. The structure of a new crystal form and its low-humidity variant.

The structures of a new crystal form of ribonuclease A and its low-humidity variant, each containing two crystallographically independent molecules, have been determined and refined. A detailed comparison of these structures with those of the other known crystal forms of the enzyme, which have different packing arrangements and solvent composition, leads to a meaningful delineation of the rigid and flexible regions of the protein molecule and the nature of its plasticity. Many of the water molecules which are common to all the structures are involved in bridging different regions of the protein molecule, thus emphasizing the role of water in stabilizing the tertiary structure. The analysis of the structures shows that for a given N or O atom, the level of hydration increases with accessible surface area, but levels off at an area of about 10 A2. Generally, the hydration level tends to drop when the area increases beyond about 20 A2. This drop correlates with an increase in the displacement parameter. The analysis also suggests that the van der Waals radii and probe radius normally used in accessible surface area calculations are not appropriate for dealing with all situations.

X-ray studies on crystalline complexes involving amino acids and peptides. XXXI. Effect of chirality on ionization state, stoichiometry and aggregation in the complexes of oxalic acid with L- and DL-histidine.

Crystals of the oxalic acid complex of L-histidine (orthorhombic P2(1)2(1)2(1); a = 5.535(4), b = 6.809(4), c = 26.878(3) A; R = 3.6% for 1188 observed reflections) contain histidine molecules and semi-oxalate ions in the 1:1 ratio, while the ratio is 1:2 in the crystals of the DL-histidine complex (monoclinic P2(1)lc; a = 6.750(7), b = 10.139(2), c = 19.352(2) A, beta = 90.8 degrees; R = 3.7% for 3176 observed reflections). The histidine molecule in the latter has an unusual ionization state with positively charged amino and imidazole groups and a neutral carboxyl group. The molecule has the sterically least favourable allowed conformation with the side chain imidazole ring staggered between the alpha-amino and the alpha- carboxyl (carboxylate) groups, in both the structures. The unlike molecules aggregate into separate alternating layers in both of them. There are elements of similarity in the aggregation patterns in the semi-oxalate layers in the two complexes, but the patterns in the amino acid layers are entirely different. Interestingly, the crystal structure of L-histidine semi-oxalate has broad similarities with that of DL-histidine = glycolate, demonstrating how broad features of aggregation could be retained inspite of changes in chirality and composition. The unusual ionization state of the amino acid molecule in the DL-histidine complex is reflected in a hitherto unobserved aggregation pattern in its crystal structure.

An X-ray analysis of native monoclinic lysozyme. A case study on the reliability of refined protein structures and a comparison with the low-humidity form in relation to mobility and enzyme action.

The atomic models of native monoclinic lysozyme obtained by refinement at Bangalore and elsewhere [Young, Dewan, Nave & Tilton (1993). J. Appl. Cryst. 26, 309-319] differed significantly in the flexible regions of the protein molecule. The two models were reconciled starting from regions where they were in reasonable agreement to produce an improved model which yielded an R value of 0.169 for 12 816 observed reflections in the 10-2 A resolution range. The reconciled model was compared with the structure of the 88% relative humidity form obtained through a water-mediated transformation [Madhusudan, Kodandapani & Vijayan (1993). Acta Cryst. D49, 234-245]. Parts of the flexible regions of the molecule register significant movements during the transformation. The changes resulting from the transformation from the native to the low-humidity forms are pronounced in many of the side chains in the active-site region, thus indicating the relationship between hydration, mobility and enzyme action. The fact that the overall changes in molecular geometry resulting from water-mediated transformation are similar to those which occur during enzyme action, further emphasizes this relationship.

Characterization of lysozyme crystals with unusually low solvent content.

Studies on the low-humidity (88%) forms of tetragonal and monoclinic lysozyme, resulting from water-mediated transformations, have provided a wealth of information on the variability in protein hydration, its structural consequences and the water structure associated with proteins, in addition to facilitating the delineation of the rigid and the flexible regions in the protein molecule and the invariant features in its hydration shell. Surprisingly, monoclinic lysozyme continues to diffract even when the environmental humidity is drastically reduced, thus permitting the structural study of the enzyme at different levels of hydration. As part of a study in this direction, three very low humidity forms, two of them occurring at a nominal relative humidity of 38% and the other at 5% relative humidity, have been characterized. These have unprecedented low solvent contents of 16.9, 17.6 and 9.4%, respectively, as determined by the Matthews method.