Progressive structure-based alignment of homologous proteins: Adopting sequence comparison strategies

Comparison of multiple protein structures has a broad range of applications in the analysis of protein structure, function and evolution. Multiple structure alignment tools (MSTAs) are necessary to obtain a simultaneous comparison of a family of related folds. In this study, we have developed a method for multiple structure comparison largely based on sequence alignment techniques. A widely used Structural Alphabet named Protein Blocks (PBs) was used to transform the information on 3D protein backbone conformation as a ID sequence string. A progressive alignment strategy similar to CLUSTALW was adopted for multiple PB sequence alignment (mulPBA). Highly similar stretches identified by the pairwise alignments are given higher weights during the alignment. The residue equivalences from PB based alignments are used to obtain a three dimensional fit of the structures followed by an iterative refinement of the structural superposition. Systematic comparisons using benchmark datasets of MSTAs underlines that the alignment quality is better than MULTIPROT, MUSTANG and the alignments in HOMSTRAD, in more than 85% of the cases. Comparison with other rigid-body and flexible MSTAs also indicate that mulPBA alignments are superior to most of the rigid-body MSTAs and highly comparable to the flexible alignment methods. (C) 2012 Elsevier Masson SAS. All rights reserved.

Structure-Based Design of DevR Inhibitor Active against Nonreplicating Mycobacterium tuberculosis

Antitubercular treatment is directed against actively replicating organisms. There is an urgent need to develop drugs targeting persistent subpopulations of Mycobacterium tuberculosis. The DevR response regulator is believed to play a key role in bacterial dormancy adaptation during hypoxia. We developed a homology-based model of DevR and used it for the rational design of inhibitors. A phenylcoumarin derivative (compound 10) identified by in silico pharmacophore-based screening of 2.5 million compounds employing protocols with some novel features including a water-based pharmacophore query, was characterized further. Compound 10 inhibited DevR binding to target DNA, down-regulated dormancy genes transcription, and drastically reduced survival of hypoxic but not nutrient-starved dormant bacteria or actively growing organ ` isms. Our findings suggest that compound 10 ``locks'' DevR in an inactive conformation that is unable to bind cognate DNA and induce the dormancy regulon. These results provide proof-of-concept for DevR as a novel target to develop molecules with sterilizing activity against tubercle bacilli.

Structure-Based Phylogeny as a Diagnostic for Functional Characterization of Proteins with a Cupin Fold

Background: The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved beta-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as ``hypothetical''. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function. Methodology/Principal Findings: A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer. Conclusions/Significance: Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.

AlignHUSH: Alignment of HMMs using structure and hydrophobicity information

Background: Sensitive remote homology detection and accurate alignments especially in the midnight zone of sequence similarity are needed for better function annotation and structural modeling of proteins. An algorithm, AlignHUSH for HMM-HMM alignment has been developed which is capable of recognizing distantly related domain families The method uses structural information, in the form of predicted secondary structure probabilities, and hydrophobicity of amino acids to align HMMs of two sets of aligned sequences. The effect of using adjoining column(s) information has also been investigated and is found to increase the sensitivity of HMM-HMM alignments and remote homology detection. Results: We have assessed the performance of AlignHUSH using known evolutionary relationships available in SCOP. AlignHUSH performs better than the best HMM-HMM alignment methods and is observed to be even more sensitive at higher error rates. Accuracy of the alignments obtained using AlignHUSH has been assessed using the structure-based alignments available in BaliBASE. The alignment length and the alignment quality are found to be appropriate for homology modeling and function annotation. The alignment accuracy is found to be comparable to existing methods for profile-profile alignments. Conclusions: A new method to align HMMs has been developed and is shown to have better sensitivity at error rates of 10% and above when compared to other available programs. The proposed method could effectively aid obtaining clues to functions of proteins of yet unknown function. A web-server incorporating the AlignHUSH method is available at http://crick.mbu.iisc.ernet.in/similar to alignhush/

Targeting Mycobacterium tuberculosis nucleoid-associated protein HU with structure-based inhibitors

The nucleoid-associated protein HU plays an important role in maintenance of chromosomal architecture and in global regulation of DNA transactions in bacteria. Although HU is essential for growth in Mycobacterium tuberculosis (Mtb), there have been no reported attempts to perturb HU function with small molecules. Here we report the crystal structure of the N-terminal domain of HU from Mtb. We identify a core region within the HU-DNA interface that can be targeted using stilbene derivatives. These small molecules specifically inhibit HU-DNA binding, disrupt nucleoid architecture and reduce Mtb growth. The stilbene inhibitors induce gene expression changes in Mtb that resemble those induced by HU deficiency. Our results indicate that HU is a potential target for the development of therapies against tuberculosis.

Structure-Based Phylogenetic Analysis of the Lipocalin Superfamily

Lipocalins constitute a superfamily of extracellular proteins that are found in all three kingdoms of life. Although very divergent in their sequences and functions, they show remarkable similarity in 3-D structures. Lipocalins bind and transport small hydrophobic molecules. Earlier sequence-based phylogenetic studies of lipocalins highlighted that they have a long evolutionary history. However the molecular and structural basis of their functional diversity is not completely understood. The main objective of the present study is to understand functional diversity of the lipocalins using a structure-based phylogenetic approach. The present study with 39 protein domains from the lipocalin superfamily suggests that the clusters of lipocalins obtained by structure-based phylogeny correspond well with the functional diversity. The detailed analysis on each of the clusters and sub-clusters reveals that the 39 lipocalin domains cluster based on their mode of ligand binding though the clustering was performed on the basis of gross domain structure. The outliers in the phylogenetic tree are often from single member families. Also structure-based phylogenetic approach has provided pointers to assign putative function for the domains of unknown function in lipocalin family. The approach employed in the present study can be used in the future for the functional identification of new lipocalin proteins and may be extended to other protein families where members show poor sequence similarity but high structural similarity.

Structure-function relationship in serine hydroxymethyltransferase

Serine hydroxymethyltransferase (SHMT), a pyridoxal-5V-phosphate (PLP)-dependent enzyme catalyzes thetetrahydrofolate (H4-folate)- dependent retro-aldol cleavage of serine to form 5,10-methylene H4-folate and glycine. The structure–function relationship of SHMT wasstudied in our laboratory initially by mutation of residues that are conserved in all SHMTs and later by structure-based mutagenesis of residues located in the active site. The analysis of mutants showed that K71, Y72, R80, D89, W110, S202, C203, H304, H306 and H356 residues are involved in maintenance of the oligomeric structure. The mutation of D227, a residue involved in charge relay system, led to the formation of inactive dimers, indicating that this residue has a role in maintaining the tetrameric structure and catalysis. E74, a residue appropriately positioned in the structure of the enzyme to carry out proton abstraction, was shown by characterization of E74Q and E74K mutants to be involved in conversion of the enzyme from an ‘open’ to ‘closed’ conformation rather than proton abstraction from the hydroxylgroup of serine. K256, the residue involved in the formation of Schiffs base with PLP, also plays a crucial role in the maintenance of the tetrameric structure. Mutation of R262 residue established the importance of distal interactions in facilitating catalysis and Y82 is not involved in the formaldehyde transfer via the postulated hemiacetal intermediate but plays a role in stabilizing the quinonoid intermediate.The mutational analysis of scSHMT along with the structure of recombinant Bacillus stearothermophilus SHMT and its substrate(s)complexes was used to provide evidence for a direct transfer mechanism rather than retro-aldol cleavage for the reaction catalyzed by SHMT.

PocketAnnotate: towards site-based function annotation

A computational pipeline PocketAnnotate for functional annotation of proteins at the level of binding sites has been proposed in this study. The pipeline integrates three in-house algorithms for site-based function annotation: PocketDepth, for prediction of binding sites in protein structures; PocketMatch, for rapid comparison of binding sites and PocketAlign, to obtain detailed alignment between pair of binding sites. A novel scheme has been developed to rapidly generate a database of non-redundant binding sites. For a given input protein structure, putative ligand-binding sites are identified, matched in real time against the database and the query substructure aligned with the promising hits, to obtain a set of possible ligands that the given protein could bind to. The input can be either whole protein structures or merely the substructures corresponding to possible binding sites. Structure-based function annotation at the level of binding sites thus achieved could prove very useful for cases where no obvious functional inference can be obtained based purely on sequence or fold-level analyses. An attempt has also been made to analyse proteins of no known function from Protein Data Bank. PocketAnnotate would be a valuable tool for the scientific community and contribute towards structure-based functional inference. The web server can be freely accessed at http://proline.biochem.iisc.ernet.in/pocketannotate/.

Integrated magnetics based multi-port bidirectional DC-DC converter topology for discontinuous-mode operation

A new type of multi-port isolated bidirectional DC-DC converter is proposed in this study. In the proposed converter, transfer of power takes place through addition of magnetomotive forces generated by multiple windings on a common transformer core. This eliminates the need for a centralised storage capacitor to interface all the ports. Hence, the requirement of an additional power transfer stage from the centralised capacitor can also be eliminated. The converter can be used for a multi-input, multi-output (MIMO) system. A pulse width modulation (PWM) strategy for controlling simultaneous power flow in the MIMO converter is also proposed. The proposed PWM scheme works in the discontinuous conduction mode. The leakage inductance can be chosen to aid power transfer. By using the proposed converter topology and PWM scheme, the need to compute power flow equations to determine the magnitude and direction of power flow between ports is alleviated. Instead, a simple controller structure based on average current control can be used to control the power flow. This study discusses the operating phases of the proposed multi-port converter along with its PWM scheme, the design process for each of the ports and finally experimental waveforms that validate the multi-port scheme.

Structural Basis of Drug Resistance by Genetic Variants of HIV Type 1 Clade C Protease from India

Using computer modeling of three-dimensional structures and structural information available on the crystal structures of HIV-1 protease, we investigated the structural effects of mutations, in treatment-naive and treatment-exposed individuals from India and postulated mechanisms of resistance in clade C variants. A large number of models (14) have been generated by computational mutation of the available crystal structures of drug bound proteases. Localized energy minimization was carried out in and around the sites of mutation in order to optimize the geometry of interactions present. Most of the mutations result in structural differences at the flap that favors the semiopen state of the enzyme. Some of the mutations were also found to confer resistance by affecting the geometry of the active site. The E35D mutation affects the flap structure in clade B strains and E35N and E35K mutation, seen in our modeled strains, have a more profound effect. Common polymorphisms at positions 36 and 63 in clade C also affected flap structure. Apart from a few other residues Gln-58, Asn-83, Asn-88, and Gln-92 and their interactions are important for the transition from the closed to the open state. Development of protease inhibitors by structure-based design requires investigation of mechanisms operative for clade C to improve the efficacy of therapy.

Design, Development, Synthesis, and Docking Analysis of 2-substituted Triclosan Analogs as Inhibitors for Plasmodium falciparum Enoyl-ACP Reductase

A structure-based approach has been adopted to develop 2'substituted analogs of triclosan. The Cl at position 2' in ring B of triclosan was chemically substituted with other functional groups like NH2, NO2 and their inhibitory potencies against PfENR were determined. The binding energies of the 2' substituted analogs of triclosan for enoyl-acyl carrier protein reductase (ENR) of Plasmodium falciparum were determined using Autodock. Based on the autodock results, we synthesized the potential compounds. The IC50 and inhibition constant (K-i) of 2' substituted analogs of triclosan were determined against purified PfENR. Among them, two compounds,2-(2'-Amino-4'-chloro-phenoxy)-5-chloro-phenol (compound 4) and 5-chloro-2-(4'-chloro-2'-nitro-phenoxy)-phenol) (compound 5) exhibited good potencies. Compound 4 followed uncompetitive inhibition kinetics with crotonoyl CoA and competitive with NADH. It was shown to have an IC50 of 110 nM; inhibition constant was 104 nM with the substrate and 61 nM with the cofactor. IC50 Of compound 5 was determined to be 229 nM. Compounds 4 and 5 showed significant inhibition of the parasite growth in P. falciparum culture. (C) 2009 IUBMB IUBMB Life, 61(11):1083-1091, 2009.

Toward the Discovery of Vaccine Adjuvants: Coupling In Silico Screening and In Vitro Analysis of Antagonist Binding to Human and Mouse CCR4 Receptors

Background: Adjuvants enhance or modify an immune response that is made to an antigen. An antagonist of the chemokine CCR4 receptor can display adjuvant-like properties by diminishing the ability of CD4+CD25+ regulatory T cells (Tregs) to down-regulate immune responses. Methodology: Here, we have used protein modelling to create a plausible chemokine receptor model with the aim of using virtual screening to identify potential small molecule chemokine antagonists. A combination of homology modelling and molecular docking was used to create a model of the CCR4 receptor in order to investigate potential lead compounds that display antagonistic properties. Three-dimensional structure-based virtual screening of the CCR4 receptor identified 116 small molecules that were calculated to have a high affinity for the receptor; these were tested experimentally for CCR4 antagonism. Fifteen of these small molecules were shown to inhibit specifically CCR4-mediated cellmigration, including that of CCR4(+) Tregs. Significance: Our CCR4 antagonists act as adjuvants augmenting human T cell proliferation in an in vitro immune response model and compound SP50 increases T cell and antibody responses in vivo when combined with vaccine antigens of Mycobacterium tuberculosis and Plasmodium yoelii in mice.

Parallel implementation of AutoDock

Computational docking of ligands to protein structures is a key step in structure-based drug design. Currently, the time required for each docking run is high and thus limits the use of docking in a high-throughput manner, warranting parallelization of docking algorithms. AutoDock, a widely used tool, has been chosen for parallelization. Near-linear increases in speed were observed with 96 processors, reducing the time required for docking ligands to HIV-protease from 81 min, as an example, on a single IBM Power-5 processor ( 1.65 GHz), to about 1 min on an IBM cluster, with 96 such processors. This implementation would make it feasible to perform virtual ligand screening using AutoDock.

Reversible electrochemistry of cytochrome c on recompressed, binderless exfoliated graphite electrodes

The direct electrochemistry of cytochrome c (cyt-c) has been investigated on exfoliated graphite (EG) electrodes. The as-polished and roughened (using SiC emery sheet) EG surfaces are inactive for the direct electron transfer. However, when the EG electrode was sonicated before the experiment, a pair of redox waves were obtained for freely diffusing cyt-c in the solution phase. The formal potential was found to be 0.01 V (vs. SCE) in 0.1 M phosphate buffer at a pH of 7.1. The electrochemical response for the adsorbed cyt-c on sonicated EG electrodes, which is shown to have carbonyl functional groups on its surface, shows nearly reversible voltammograms in the same electrolyte. However, the formal potential in the adsorbed state is more negative than that observed for the solution phase cyt-c. A structure based on an open heme conformation proposed by Hildebrandt and Stockburger is probably present on the EG surface. It is suggested that the electrochemistry at the EG electrode is essentially governed by favourable electrostatic interactions.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.