Bhupesh Taneja

Open source drug discovery– A new paradigm of collaborative research in tuberculosis drug development

It is being realized that the traditional closed-door and market driven approaches for drug discovery may not be the best suited model for the diseases of the developing world such as tuberculosis and malaria, because most patients suffering from these diseases have poor paying capacity. To ensure that new drugs are created for patients suffering from these diseases, it is necessary to formulate an alternate paradigm of drug discovery process. The current model constrained by limitations for collaboration and for sharing of resources with confidentiality hampers the opportunities for bringing expertise from diverse fields. These limitations hinder the possibilities of lowering the cost of drug discovery. The Open Source Drug Discovery project initiated by Council of Scientific and Industrial Research, India has adopted an open source model to power wide participation across geographical borders. Open Source Drug Discovery emphasizes integrative science through collaboration, open-sharing, taking up multi-faceted approaches and accruing benefits from advances on different fronts of new drug discovery. Because the open source model is based on community participation, it has the potential to self-sustain continuous development by generating a storehouse of alternatives towards continued pursuit for new drug discovery. Since the inventions are community generated, the new chemical entities developed by Open Source Drug Discovery will be taken up for clinical trial in a non-exclusive manner by participation of multiple companies with majority funding from Open Source Drug Discovery. This will ensure availability of drugs through a lower cost community driven drug discovery process for diseases afflicting people with poor paying capacity. Hopefully what LINUX the World Wide Web have done for the information technology, Open Source Drug Discovery will do for drug discovery.

Structure of Mycobacterium tuberculosis chaperonin-10 at 3.5 ˚ A resolution

Chaperonin-60 (cpn60) and chaperonin-10 (cpn10) are essential proteins involved in ATP-dependent folding of several intracellular proteins in the bacterial cell. Folding of the nascent substrate polypeptide takes place in the large central cavity formed by each ring of the tetradecameric cpn60. This large cavity is closed upon capping by the heptameric cpn10. Cpn10s interact with cpn60s primarily through a 17-residue mobile loop and regulate the release and binding of the substrate polypeptide from the cpn60 surface. Here, the structure of M. tuberculosis cpn10 is reported at 3.5 A resolution. The overall structure of the cpn10 monomer is formed of a four-stranded beta-barrel and two long stretches of highly flexible segments: the dome loop and the mobile loop. The seven subunits in the heptamer show very little conformational difference and exhibit nearly perfect sevenfold geometry. The binding sites for metal ions in the dome loop of cpn10 have been identified, suggesting the role of metal ions in the stabilization of the protein. Comparisons with the available cpn10 structures indicate several interesting features.

Structural and Functional Characterization of Rv2966c Protein Reveals an RsmD-like Methyltransferase from Mycobacterium tuberculosis and the Role of Its N-terminal Domain in Target Recognition

Nine of ten methylated nucleotides of Escherichia coli 16 S rRNA are conserved in Mycobacterium tuberculosis. All the 10 different methyltransferases are known in E. coli, whereas only TlyA and GidB have been identified in mycobacteria. Here we have identified Rv2966c of M. tuberculosis as an ortholog of RsmD protein of E. coli. We have shown that rv2966c can complement rsmD-deleted E. coli cells. Recombinant Rv2966c can use 30 S ribosomes purified from rsmD-deleted E. coli as substrate and methylate G966 of 16 S rRNA in vitro. Structure determination of the protein shows the protein to be a two-domain structure with a short hairpin domain at the N terminus and a C-terminal domain with the S-adenosylmethionine-MT-fold. We show that the N-terminal hairpin is a minimalist functional domain that helps Rv2966c in target recognition. Deletion of the N-terminal domain prevents binding to nucleic acid substrates, and the truncated protein fails to carry out the m2G966 methylation on 16 S rRNA. The N-terminal domain also binds DNA efficiently, a property that may be utilized under specific conditions of cellular growth.

Unique subunit packing in mycobacterial nanoRNase leads to alternate substrate recognitions in DHH phosphodiesterases

DHH superfamily includes RecJ, nanoRNases (NrnA), cyclic nucleotide phosphodiesterases and pyrophosphatases. In this study, we have carried out in vitro and in vivo investigations on the bifunctional NrnA-homolog from Mycobacterium smegmatis, MSMEG_2630. The crystal structure of MSMEG_2630 was determined to 2.2-Å resolution and reveals a dimer consisting of two identical subunits with each subunit folding into an N-terminal DHH domain and a C-terminal DHHA1 domain. The overall structure and fold of the individual domains is similar to other members of DHH superfamily. However, MSMEG_2630 exhibits a distinct quaternary structure in contrast to other DHH phosphodiesterases. This novel mode of subunit packing and variations in the linker region that enlarge the domain interface are responsible for alternate recognitions of substrates in the bifunctional nanoRNases. MSMEG_2630 exhibits bifunctional 3′-5′ exonuclease [on both deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) substrates] as well as CysQ-like phosphatase activity (on pAp) in vitro with a preference for nanoRNA substrates over single-stranded DNA of equivalent lengths. A transposon disruption of MSMEG_2630 in M. smegmatis causes growth impairment in the presence of various DNA-damaging agents. Further phylogenetic analysis and genome organization reveals clustering of bacterial nanoRNases into two distinct subfamilies with possible role in transcriptional and translational events during stress.

The structure of Rv3717 reveals a novel amidase from Mycobacterium tuberculosis

The structure of Rv3717 determined to 1.7 Å resolution by Pt-SAD phasing reveals a unique autolysin that lacks a cell-wall-binding domain. Rv3717 utilizes its net positive charge for substrate binding and exhibits activity towards a broad spectrum of substrate cell walls. Structural analysis reveals that Rv3717 utilizes a β-hairpin turn at its N-terminus to autoregulate its enzymatic activity.

Integration Host Factor of Mycobacterium tuberculosis, mIHF, Compacts DNA by a Bending Mechanism

The bacterial chromosomal DNA is folded into a compact structure called as ‘nucleoid’ so that the bacterial genome can be accommodated inside the cell. The shape and size of the nucleoid are determined by several factors including DNA supercoiling, macromolecular crowding and nucleoid associated proteins (NAPs). NAPs bind to different sites of the genome in sequence specific or non-sequence specific manner and play an important role in DNA compaction as well as regulation. Until recently, few NAPs have been discovered in mycobacteria owing to poor sequence similarities with other histone-like proteins of eubacteria. Several putative NAPs have now been identified in Mycobacteria on the basis of enriched basic residues or histone-like “PAKK” motifs. Here, we investigate mycobacterial Integration Host Factor (mIHF) for its architectural roles as a NAP using atomic force microscopy and DNA compaction experiments. We demonstrate that mIHF binds DNA in a non-sequence specific manner and compacts it by a DNA bending mechanism. AFM experiments also indicate a dual architectural role for mIHF in DNA compaction as well as relaxation. These results suggest a convergent evolution in the mechanism of E. coli and mycobacterial IHF in DNA compaction.

Identification of L84F mutation with a novel nucleotide change c.255G > T in the superoxide dismutase gene in a North Indian family with amyotrophic lateral sclerosis

Abstract Mutations in the superoxide dismutase (SOD1) gene account for ∼15% and in the transactive response DNA binding protein (TARDBP) gene for ∼5% of familial amyotrophic lateral sclerosis (FALS) cases. These two genes were analysed in two siblings from North India with ALS and a positive family history. The coding region of SOD1 and TARDBP genes was sequenced in both siblings. Genetic variation identified in SOD1 was typed in unaffected family members (n = 11), sporadic ALS patients (n = 48) and healthy controls (n = 35). Molecular dynamic (MD) simulations were performed on wild-type (WT) and mutant monomers of SOD1 to determine structural changes due to the identified mutation. A novel heterozygous nucleotide variation (c.255G > T) was identified in exon 4 of SOD1 in the two siblings and two asymptomatic family members but not in SALS patients and healthy controls. This variation results in a known non-synonymous substitution from leucine to phenylalanine at position 84 (L84F), making it a triallelic variation. Large conformational changes were observed in the zinc loop and electrostatic loop in an L84F mutant compared to WT SOD1 in MD simulations. In conclusion, this is the first report of mutation in SOD1 associated with FALS in India. Structural perturbations in L84F SOD1 may cause dimer destabilization, with decreased metal affinity leading to oligomerization.

Structural and thermodynamic characterisation of L94F mutant of horse cytochrome c.

Mammalian mitochondrial cytochromes c (cyts c) has a conserved Leu94 which is replaced by valine/isoleucine in some lower eukaryotes and prokaryotes. It is expected that nature substituted Leu94 with Val/Ile, for they have similar van der Waals volume and hydrophobicity with a difference in side chain branching only. Reports also suggested the presence of phenylalanine at position 94, which leads to questions: (i) How bulky aromatic amino acid residue fitted at position 94 in cyt c family proteins? (ii) What is the effect of L94F mutation on protein stability and folding? Here, we selected horse cyt-c as a model to answer the second question. We generated L94F mutant of horse cytochrome c and subsequently characterised using far-UV, near-UV and Soret circular dichroism, absorbance, intrinsic and extrinsic ANS (8-anilino-1-napthalenesulfonic acid) fluorescence and dynamic light scattering measurements. We observed that this mutation affects the native state and arrests the protein folding at the molten globule state. Thermal stability of L94F mutant is also measured by spectroscopic techniques and differential scanning calorimetry. The midpoint of thermal denaturation of L94F mutant is 17°C less than wild type. Molecular dynamics simulation study also supports our in vitro observation that this mutant has stable backbone conformation.

Computational screening for new inhibitors of M. tuberculosis mycolyltransferases antigen 85 group of proteins as potential drug targets.

The group of antigen 85 proteins of Mycobacterium tuberculosis is responsible for converting trehalose monomycolate to trehalose dimycolate, which contributes to cell wall stability. Here, we have used a serial enrichment approach to identify new potential inhibitors by searching the libraries of compounds using both 2D atom pair descriptors and binary fingerprints followed by molecular docking. Three different docking softwares AutoDock, GOLD, and LigandFit were used for docking calculations. In addition, we applied the criteria of selecting compounds with binding efficiency close to the starting known inhibitor and showing potential to form hydrogen bonds with the active site amino acid residues. The starting inhibitor was ethyl-3-phenoxybenzyl-butylphosphonate, which had IC50 value of 2.0 μM in mycolyltransferase inhibition assay. Our search from more than 34 million compounds from public libraries yielded 49 compounds. Subsequently, selection was restricted to compounds conforming to the Lipinski rule of five and exhibiting hydrogen bonding to any of the amino acid residues in the active site pocket of all three proteins of antigen 85A, 85B, and 85C. Finally, we selected those ligands which were ranked top in the table with other known decoys in all the docking results. The compound NIH415032 from tuberculosis antimicrobial acquisition and coordinating facility was further examined using molecular dynamics simulations for 10 ns. These results showed that the binding is stable, although some of the hydrogen bond atom pairs varied through the course of simulation. The NIH415032 has antitubercular properties with IC90 at 20 μg/ml (53.023 μM). These results will be helpful to the medicinal chemists for developing new antitubercular molecules for testing.

Physical Properties of Intact Proteins May Predict Allergenicity or Lack Thereof

Background Predicting the allergenicity of proteins is challenging. We considered the possibility that the properties of the intact protein that may alter the likelihood of being taken up by antigen presenting cells, may be useful adjuncts in predicting allergens and non-allergens in silico. It has been shown that negatively charged acidic proteins are preferentially processed by dendritic cells. Methodology Datasets (aeroallergen, food-allergen and non-allergen) for in-silico study were obtained from public databases. Isoelectric point (pI), net charge, and electrostatic potential (EP) were calculated from the protein sequence (for pI and net charge) or predicted structure (for EP). Result Allergens and non allergens differed significantly in pI, net charge, and EP (p<0.0001). Cluster analysis based on these parameters resulted in well defined clusters. Non-allergens were characterized by neutral to basic pI (mean±SE, 7.6±0.16) and positive charge. In contrast allergens were acidic (5.7±0.15) and negatively charged. Surface electrostatic potentials calculated from predicted structures were mostly negative for allergens and mostly positive for non-allergens. The classification accuracy for non-allergens was superior to that for allergens. Thus neutral to basic pI, positive charge, and positive electrostatic potentials characterize non-allergens, and seem rare in allergens (p<0.0001). It may be possible to predict reduced likelihood of allergenicity in such proteins, but this needs to be prospectively validated.