Identification of potential molecular targets and repurposed drugs for tuberculosis using network-based screening approach, molecular docking, and simulation.

The spread of drug-resistant strains of tuberculosis has hampered efforts to control the disease worldwide. The Mycobacterium tuberculosis cell wall envelope is dynamic, with complex features that protect it from the host immunological response. As a result, the bacterial cell wall components represent a potential target for drug discovery. Protein-protein interaction networks (PPIN) are critical for understanding disease conditions and identifying precise therapeutic targets. We used a rational theoretical approach by constructing a PPIN with the proteins involved in cell wall biosynthesis. The PPIN was constructed through the STRING database and embB was identified as a key protein by using four topological measures, betweenness, closeness, degree, and eigenvector, in the CytoNCA tool in Cytoscape. The 'Drug repurposing' approach was employed to find suitable inhibitors against embB. We used the Schrödinger suites for molecular docking, molecular dynamics simulation, and binding free energy calculations to validate the binding of protein with the ligand. FDA-approved drugs from the ZINC database and DrugBank were screened against embB (PDB ID: 7BVF) using high-throughput virtual screening, standard precision, and extra precision docking. The drugs were screened based on the XP docking score of the standard drug ethambutol. Accordingly, from the top five hits, azilsartan and dihydroergotamine were selected based on the binding free energy values and were further subjected to Molecular Dynamics Simulation studies for 100 ns. Our study confirms that Azilsartan and Dihydroergotamine form stable complexes with embB and can be used as potential lead molecules based on further in vitro and in vivo experimental validation.Communicated by Ramaswamy H. Sarma.

Conserved acidic second shell residue modulates the structure, stability and activity of non-seleno human peroxiredoxin 6.

1-Cys peroxiredoxin6 (Prdx6) is unique and inducible bifunctional enzyme in the mammalian lungs and plays a role in the progression and inhibition of cancerous cells at different stages. The enzyme possesses two distinct active sites for phospholipase A2 and peroxidase activity. The conserved residues surrounding the peroxidase active site, also called as second shell residues are Glu50, Leu71, Ser72, His79 and Arg155. Since there is no study done about the active site stabilization of the transition state of Prdx6, there are a lot of questions unanswered regarding the Prdx6 peroxidase activity. In order to evaluate the role of second shell conserved residue Glu50, present in close vicinity to peroxidatic active site, we substituted this negatively charged residue with Alanine and Lysine. To explore the effect of mutation on the biophysical parameters, the mutant proteins were compared with Wild-Type by using biochemical, biophysical, and in silico methods. Comparative spectroscopic methods and enzyme activity demonstrate that the Glu50 plays a significant role in maintaining the structure, stability, and function of protein. From the results we conclude that Glu50 significantly controls the structure; stability and may be involved in the active site stabilization of transition state for proper position of diverse peroxides.

Investigating the binding mechanism of topiramate with bovine serum albumin using spectroscopic and computational methods.

Various spectroscopic techniques involving fluorescence spectroscopy, circular dichroism (CD), and computational approaches were used to elucidate the molecular aspects of interaction between the antiepileptic drug topiramate and the multifunctional transport protein bovine serum albumin (BSA) under physiological conditions. Topiramate quenched BSA fluorescence in a static quenching mode, according to the Stern-Volmer quenching constant (Ksv ) data derived from fluorescence spectroscopy for the topiramate-BSA complex. The binding constant was also used to calculate the binding affinity for the topiramate-BSA interaction. Fluorescence and circular dichroism experiments demonstrate that the protein's tertiary structure is affected by the microenvironmental alterations generated by topiramate binding to BSA. To establish the exact binding site, interacting residues, and interaction forces involved in the binding of topiramate to BSA, molecular modeling and simulation approaches were used. According to the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) calculations, the average binding energy between topiramate and BSA is -421.05 kJ/mol. Topiramate was discovered to have substantial interactions with BSA, changing the structural dynamic and Gibbs free energy landscape patterns.

Mechanistic insights into the urea-induced denaturation of a non-seleno thiol specific antioxidant human peroxiredoxin 6.

Peroxiredoxin 6 (Prdx6) is a unique enzyme among mammalian peroxiredoxins as it lacks resolving cysteine. It is found to be involved in number of different diseases including tumours and its expression level is highest in lungs as compared to other organs. It has been found that Prdx6 plays a significant role different metabolic diseases, ocular damage, neurodegeneration and male infertility. It is a bifunctional protein having phospholipase A2 and peroxidase (also has the ability to reduce phospholipid hydroperoxides) activities. In order to complete the peroxidise reaction cycle it requires glutathione catalyzed by glutathione S-transferase. Equilibrium unfolding and conformational stability of Prdx6 was studied by using urea as a chemical denaturant to understand the changes it goes under cellular stress conditions. Three different spectroscopic methods were employed to monitor urea-induced denaturation. From the results obtained, it was found that the urea denaturation of Prdx6 follows a variable two state process due to non-coincidence of the normalized transition curves obtained from different optical probes. The different denaturation curves were normalized and thermodynamic parameters, ΔGDo, Gibbs free energy change related to the urea-induced denaturation, midpoint of denaturation (Cm), and m = (δΔGD / [urea]) were obtained. The structural information of Prdx6 were further analysed by several parameters obtained by 100 ns MD simulation. The results of MD simulation clearly favour the outcome of spectroscopic studies.

Emerging biomarkers and clinical significance of HPV genotyping in prevention and management of cervical cancer.

Cervical cancer is a growing and serious problem world-wide in women, but more acute in developing countries especially in Indian subcontinent. The main causative agent for the disease is Human Papilloma Virus (HPV). The history of the cervical cancer goes back to eighteenth century as the HPV infection is reported since 1800s. Presently, the genetic structure of HPV is well defined. Several screening tests including cytology and visual based screening and high risk HPV testing are available. Also available are various clinical and commercial diagnostic tests. However due to the lack of awareness and population-based screening programs, the morbidity and mortality rate is alarmingly high. There are new emerging biomarkers including E6/E7 mRNA, p16ink4a, markers of aberrant S-phase induction, chromosomal abnormalities and miRNAs along with advanced genotyping methods. These markers have clinical significance and are helpful in disease prevention and management. Further, recent advancement in the field of metagenomics has increased the prospects of identifying newer microbes, viruses hitherto reported thus far in the context of HPV infection. Analysis of HPV cases using modern tools including genotyping using more powerful biomarkers is envisaged to enhance the prospects of early diagnosis, better prognosis, more reliable treatment and eventual management of the disease.

In silico based unraveling of New Delhi metallo-ß-lactamase (NDM-1) inhibitors from natural compounds: a molecular docking and molecular dynamics simulation study.

The development of pathogenic microbial resistance toward antibiotics has become a global clinical concern. New Delhi metallo-ß-lactmase-1 (NDM-1) and its variants have recently drawn immense attention for its biological ability to catalyze the hydrolysis of almost all of ß-lactam antibiotics including the Carbapenems which are generally considered as the last-resort antibiotics. Also, the horizontal gene transfer is expediting the rapid spread of NDM-1 in bacteria. In the wake of this serious antibiotic resistance problem it becomes imperative to find inhibitors which can render the present antibiotics functional and useful. In the present study, we have used Molecular docking and Molecular Dynamics (MD) simulation approach to find out suitable inhibitors against NDM-1 from an array of different natural compounds. We have screened unique natural compounds from ZINC database and also a set of standard antibiotics and inhibitors. Based upon the highest binding affinity demonstrated by docking with NDM-1, the best binding antibiotic Meropenem and the top five natural compounds, viz., Withaferin A, Beta-Sitosterol, Aristolochic acid, Diosgenin and Guggulsterone E were selected and subjected to MD simulations study. The docked NDM-1 complex with withaferin A, beta-sitosterol and diosgenin were found to be more stable as compared to the one with meropenem throughout the MD simulation process with the relative RMSD and RMSF in acceptable range. In conclusion, these compounds can be readily tested in vitro and in vivo to fully establish and confirm their inhibition potentiality and can also serve as lead molecules for the development of future functional inhibitors.Communicated by Ramaswamy H. Sarma.

Sub-Micellar Concentration of Sodium Dodecyl Sulphate Prevents Thermal Denaturation Induced Aggregation of Plant Lectin, Jacalin.

The irreversible thermal unfolding of jacalin, the lectin purified from jackfruit seeds was accompanied by aggregation, where intermolecular interactions among the subunits are favoured over intramolecular interactions. The extent of aggregation increased as a function of temperature, time and protein concentration. The anionic surfactant, sodium dodecyl sulphate (SDS) significantly suppressed the formation of aggregates as observed by turbidity measurements and Rayleigh scattering assay. Moreover, far UV-CD spectra indicate that the protein ß sheet transforms into α helical structure, when denatured in the presence of 3 mM SDS. Further, jacalin when heated in the presence of SDS partially retained the hemagglutination activity when jacalin-SDS mixture was diluted to 1:8 factor since 3 mM SDS was found to lyse the red blood cells. Thus, SDS only altered the aggregation behaviour of jacalin by preventing intermolecular hydrogen bonding among the exposed residues but did not completely stabilize the native conformation.

Sustained mitogenic effect on K562 human chronic myelogenous leukemia cells by dietary lectin, jacalin.

Dietary lectins have been shown to affect the proliferation of human cancer cell lines. The anti-proliferative effects of lectins from varied sources have been extensively studied and in some cases, the underlying mechanism has been explored. Except for peanut agglutinin (PNA), the mitogenic effects of no other lectins have been studied in detail. In the present study, we have shown that jacalin, lectin purified from jackfruit (Artocarpus integrifolia) seeds act as a mitogen for K562, the Bcr-Abl expressing erythroleukemia cell line (K562) and the effect was found to be dose dependent. K562 cells remained in the proliferative state for a longer period even after the withdrawal of jacalin stimulation, thus jacalin was found to induce sustained mitogenic effect on K562 cells. Further, conditioned media from K562 cells treated with jacalin were observed to have the similar mitogenic effect even in the presence of galactose. Importantly, galactose which is a known ligand for jacalin will interact with functionally active jacalin present in the conditioned media and neutralise its effect. In addition, jacalin treatment also resulted in increased mRNA expression levels of pro-inflammatory cytokines including IL-1ß, IL-6 and IFN-γ. Our results indicate that jacalin induces secretion of soluble molecules, which maybe responsible for this observed increased proliferation of K562 cells.

The response of Ω-loop D dynamics to truncation of trimethyllysine 72 of yeast iso-1-cytochrome c depends on the nature of loop deformation.

Trimethyllysine 72 (tmK72) has been suggested to play a role in sterically constraining the heme crevice dynamics of yeast iso-1-cytochrome c mediated by the Ω-loop D cooperative substructure (residues 70-85). A tmK72A mutation causes a gain in peroxidase activity, a function of cytochrome c that is important early in apoptosis. More than one higher energy state is accessible for the Ω-loop D substructure via tier 0 dynamics. Two of these are alkaline conformers mediated by Lys73 and Lys79. In the current work, the effect of the tmK72A mutation on the thermodynamic and kinetic properties of wild-type iso-1-cytochrome c (yWT versus WT*) and on variants carrying a K73H mutation (yWT/K73H versus WT*/K73H) is studied. Whereas the tmK72A mutation confers increased peroxidase activity in wild-type yeast iso-1-cytochrome c and increased dynamics for formation of a previously studied His79-heme alkaline conformer, the tmK72A mutation speeds return of the His73-heme alkaline conformer to the native state through destabilization of the His73-heme alkaline conformer relative to the native conformer. These opposing behaviors demonstrate that the response of the dynamics of a protein substructure to mutation depends on the nature of the perturbation to the substructure. For a protein substructure which mediates more than one function of a protein through multiple non-native structures, a mutation could change the partitioning between these functions. The current results suggest that the tier 0 dynamics of Ω-loop D that mediates peroxidase activity has similarities to the tier 0 dynamics required to form the His79-heme alkaline conformer.

Structural characteristic of the initial unfolded state on refolding determines catalytic efficiency of the folded protein in presence of osmolytes.

Osmolytes are low molecular weight organic molecules accumulated by organisms to assist proper protein folding, and to provide protection to the structural integrity of proteins under denaturing stress conditions. It is known that osmolyte-induced protein folding is brought by unfavorable interaction of osmolytes with the denatured/unfolded states. The interaction of osmolyte with the native state does not significantly contribute to the osmolyte-induced protein folding. We have therefore investigated if different denatured states of a protein (generated by different denaturing agents) interact differently with the osmolytes to induce protein folding. We observed that osmolyte-assisted refolding of protein obtained from heat-induced denatured state produces native molecules with higher enzyme activity than those initiated from GdmCl- or urea-induced denatured state indicating that the structural property of the initial denatured state during refolding by osmolytes determines the catalytic efficiency of the folded protein molecule. These conclusions have been reached from the systematic measurements of enzymatic kinetic parameters (Km and kcat), thermodynamic stability (Tm and ΔHm) and secondary and tertiary structures of the folded native proteins obtained from refolding of various denatured states (due to heat-, urea- and GdmCl-induced denaturation) of RNase-A in the presence of various osmolytes.

Conformational properties of polyglutamine sequences in guanidine hydrochloride solutions.

Two sets of iso-1-cytochrome c variants have been prepared with N-terminal insertions of pure polyglutamine, i.e., PolyQ variants, or polyglutamine interrupted with lysine every sixth residue, i.e., Gln-rich variants. The polymer properties of these pure polyGln or Gln-rich sequences have been evaluated using equilibrium and kinetic His-heme loop formation methods for loop sizes ranging from 22 to 46 in 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). In 6.0 M GdnHCl, the scaling exponent, ν(3), for the pure polyGln sequences, is ~1.7--significantly less than ν(3) ≈ 2.15 for the Gln-rich sequences. The stability of the His-heme loops becomes progressively greater for the pure polyGln sequences relative to the Gln-rich sequences as GdnHCl concentration decreases from 6.0 to 1.5 M. Thus, the context of the sequence effects the polymer properties of Gln repeats even in denaturing concentrations of GdnHCl. Comparison of data for the Gln-rich variants with previous results for Gly-rich and Ala-rich variants shows that ν(3) ~ 2.2 for the Gln-rich, Gly-rich, and Ala-rich sequences in 6.0 M GdnHCl, whereas ν(3) remains unchanged at 3.0 M GdnHCl concentration for the Gln-rich and Ala-rich sequences but decreases to ~1.7 for the Gly-rich sequences. Thus, the polymer properties of Gln-rich and Ala-rich sequences are less sensitive to solvent quality in denaturing solutions of GdnHCl than Gly-rich sequences. Evaluation of Flory's characteristic ratio, C(n), for the Gln-rich and Ala-rich sequences relative to the Gly-rich sequences shows that Gln-rich sequences are stiffer than Ala-rich sequences at both 3.0 and 6.0 M GdnHCl.

Conformation and thermodynamic stability of pre-molten and molten globule states of mammalian cytochromes-c.

Proteins in cells fold via a number of intermediates. These intermediates are quite important as they guide the protein to attain its unique native conformation. To solve the immensely difficult problem of protein folding, it is necessary to characterize intermediates which will unravel the mystery of the steps involved in the proper folding of proteins. Cytochromes-c (cyts-c) have played an important role in studies of the earliest events and intermediates in protein folding. They have always been considered as model proteins for protein folding studies due to their intrinsic properties that can be measured by multiple probes. A large number of different solvent conditions have been employed to obtain equilibrium intermediates of cyts-c. These intermediates show structural heterogeneity which is mainly due to the different solvent conditions used to induce them. In this review we present results of conformational and thermodynamic characterization of equilibrium intermediates (molten globules and pre-molten globules) of the mammalian cyts-c under different solvent conditions.