DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

The destabilization of human GCAP1 by a proline to leucine mutation might cause cone-rod dystrophy

Guanylate cyclase activating protein-1 (GCAP1) is required for activation of retinal guanylate cyclase-1 (RetGC1), which is essential for recovery of photoreceptor cells to the dark state. In this paper, experimentally derived observations are reported that help in explaining why a proline→leucine mutation at position 50 of human GCAP1 results in cone–rod dystrophy in a family carrying this mutation. The primary amino acid sequence of wild-type GCAP1 was mutated using site-directed mutagenesis to give a leucine at position 50. In addition, serine replaced a glutamic acid residue at position 6 to promote N‐terminal myristoylation, yielding the construct GCAP1 E6S/P50L. The enzyme was over-expressed in Escherichia coli cells, isolated and purified before being used in assays with RetGC1, characterized by circular dichroism (CD) spectroscopy, and investigated for protease resistance and thermal stability. Assays of cyclic guanosine monophosphate (cGMP) synthesis from guanosine triphosphate by RetGC1 in the presence of E6S/P50L showed that E6S/P50L could activate RetGC1 and displayed similar calcium sensitivity to wild-type GCAP1. In addition, E6S/P50L and wild-type GCAP1 possess similar CD spectra. However, there was a marked increase in the susceptibility to protease degradation and also a reduction in the thermal stability of E6S/P50L as observed by both the cGMP assay and CD spectroscopy. It is therefore suggested that although GCAP1 E6S/P50L has a similar activity and calcium dependency profile to the wild-type GCAP1, its lower stability could reduce its cellular concentration, which would in turn alter [Ca2+] and result in death of cells.

Stability of domain structures in multi-domain proteins

Multi-domain proteins have many advantages with respect to stability and folding inside cells. Here we attempt to understand the intricate relationship between the domain-domain interactions and the stability of domains in isolation. We provide quantitative treatment and proof for prevailing intuitive ideas on the strategies employed by nature to stabilize otherwise unstable domains. We find that domains incapable of independent stability are stabilized by favourable interactions with tethered domains in the multi-domain context. Stability of such folds to exist independently is optimized by evolution. Specific residue mutations in the sites equivalent to inter-domain interface enhance the overall solvation, thereby stabilizing these domain folds independently. A few naturally occurring variants at these sites alter communication between domains and affect stability leading to disease manifestation. Our analysis provides safe guidelines for mutagenesis which have attractive applications in obtaining stable fragments and domain constructs essential for structural studies by crystallography and NMR.

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.

Mechanism of Fragility at BCL2 Gene Minor Breakpoint Cluster Region during t(14;18) Chromosomal Translocation

The t(14;18) translocation in follicular lymphoma is one of the most common chromosomal translocations. Breaks in chromosome 18 are localized at the 3'-UTR of BCL2 gene or downstream and are mainly clustered in either the major breakpoint region or the minor breakpoint cluster region (mcr). The recombination activating gene (RAG) complex induces breaks at IgH locus of chromosome 14, whereas the mechanism of fragility at BCL2 mcr remains unclear. Here, for the first time, we show that RAGs can nick mcr; however, the mechanism is unique. Three independent nicks of equal efficiency are generated, when both Mg2+ and Mn2+ are present, unlike a single nick during V(D)J recombination. Further, we demonstrate that RAG binding and nicking at the mcr are independent of nonamer, whereas a CCACCTCT motif plays a critical role in its fragility, as shown by sequential mutagenesis. More importantly, we recapitulate the BCL2 mcr translocation and find that mcr can undergo synapsis with a standard recombination signal sequence within the cells, in a RAG-dependent manner. Further, mutation to the CCACCTCT motif abolishes recombination within the cells, indicating its vital role. Hence, our data suggest a novel, physiologically relevant, nonamer-independent mechanism of RAG nicking at mcr, which may be important for generation of chromosomal translocations in humans.

Crystal Structure of Escherichia coli Diaminopropionate Ammonia-lyase Reveals Mechanism of Enzyme Activation and Catalysis

Pyridoxal 5'-phosphate (PLP)-dependent enzymes utilize the unique chemistry of a pyridine ring to carry out diverse reactions involving amino acids. Diaminopropionate (DAP) ammonia-lyase (DAPAL) is a prokaryotic PLP-dependent enzyme that catalyzes the degradation of D-and L-forms of DAP to pyruvate and ammonia. Here, we report the first crystal structure of DAPAL from Escherichia coli (EcDAPAL) in tetragonal and monoclinic forms at 2.0 and 2.2 angstrom resolutions, respectively. Structures of EcDAPAL soaked with substrates were also determined. EcDAPAL has a typical fold type II PLP-dependent enzyme topology consisting of a large and a small domain with the active site at the interface of the two domains. The enzyme is a homodimer with a unique biological interface not observed earlier. Structure of the enzyme in the tetragonal form had PLP bound at the active site, whereas the monoclinic structure was in the apo-form. Analysis of the apo and holo structures revealed that the region around the active site undergoes transition from a disordered to ordered state and assumes a conformation suitable for catalysis only upon PLP binding. A novel disulfide was found to occur near a channel that is likely to regulate entry of ligands to the active site. EcDAPAL soaked with DL-DAP revealed density at the active site appropriate for the reaction intermediate aminoacrylate, which is consistent with the observation that EcDAPAL has low activity under crystallization conditions. Based on the analysis of the structure and results of site-directed mutagenesis, a two-base mechanism of catalysis involving Asp(120) and Lys(77) is suggested.

Cloning, expression, purification, and biochemical characterisation of the FIC motif containing protein of Mycobacterium tuberculosis

The role of FIC (Filamentation induced by cAMP)(2) domain containing proteins in the regulation of many vital pathways, mostly through the transfer of NMPs from NTPs to specific target proteins (NMPylation), in microorganisms, higher eukaryotes, and plants is emerging. The identity and function of FIC domain containing protein of the human pathogen, Mycobacterium tuberculosis, remains unknown. In this regard, M. tuberculosis fic gene (Mtfic) was cloned, overexpressed, and purified to homogeneity for its biochemical characterisation. It has the characteristic FIC motif, HPFREGNGRSTR (HPFxxGNGRxxR), spanning 144th to 155th residue. Neither the His-tagged nor the GST-tagged MtFic protein, overexpressed in Escherichia coil, nor expression of Mtfic in Mycobacterium smegmatis, yielded the protein in the soluble fraction. However, the maltose binding protein (MBP) tagged MtFic (MBP-MtFic) could be obtained partly in the soluble fraction. The cloned, overexpressed, and purified recombinant MBP-MtFic showed conversion of ATP, GTP, CTP, and UTP into AMP. GMP, CMP, and UMP, respectively. Sequence alignment with several FIC motif containing proteins, complemented with homology modeling on the FIC motif containing protein, VbhT of Bartonella schoenbuchensis as the template, showed conservation and interaction of residues constituting the FIC domain. Site-specific mutagenesis of the His144, or Glu148, or Asn150 of the FIC motif, or of Arg87 residue that constitutes the FIC domain, or complete deletion of the FIC motif, abolished the NTP to NMP conversion activity. The design of NMP formation assay using the recombinant, soluble MtFic would enable identification of its target substrate for NMPylation. (C) 2012 Elsevier Inc. All rights reserved.

Site-specific N-Linked Glycosylation of Receptor Guanylyl Cyclase C Regulates Ligand Binding, Ligand-mediated Activation and Interaction with Vesicular Integral Membrane Protein 36, VIP36

Guanylyl cyclase C (GC-C) is a multidomain, membrane-associated receptor guanylyl cyclase. GC-C is primarily expressed in the gastrointestinal tract, where it mediates fluid-ion homeostasis, intestinal inflammation, and cell proliferation in a cGMP-dependent manner, following activation by its ligands guanylin, uroguanylin, or the heat-stable enterotoxin peptide (ST). GC-C is also expressed in neurons, where it plays a role in satiation and attention deficiency/hyperactive behavior. GC-C is glycosylated in the extracellular domain, and differentially glycosylated forms that are resident in the endoplasmic reticulum (130 kDa) and the plasma membrane (145 kDa) bind the ST peptide with equal affinity. When glycosylation of human GC-C was prevented, either by pharmacological intervention or by mutation of all of the 10 predicted glycosylation sites, ST binding and surface localization was abolished. Systematic mutagenesis of each of the 10 sites of glycosylation in GC-C, either singly or in combination, identified two sites that were critical for ligand binding and two that regulated ST-mediated activation. We also show that GC-C is the first identified receptor client of the lectin chaperone vesicular integral membrane protein, VIP36. Interaction with VIP36 is dependent on glycosylation at the same sites that allow GC-C to fold and bind ligand. Because glycosylation of proteins is altered in many diseases and in a tissue-dependent manner, the activity and/or glycan-mediated interactions of GC-C may have a crucial role to play in its functions in different cell types.

Dramatic Structural Changes Resulting from the Loss of a Crucial Hydrogen Bond in the Hinge Region Involved in C-Terminal Helix Swapping in SurE: A Survival Protein from Salmonella typhimurium

Domain swapping is an interesting feature of some oligomeric proteins in which each protomer of the oligomer provides an identical surface for exclusive interaction with a segment or domain belonging to another protomer. Here we report results of mutagenesis experiments on the structure of C-terminal helix swapped dimer of a stationary phase survival protein from Salmonella typhimurium (StSurE). Wild type StSurE is a dimer in which a large helical segment at the C-terminus and a tetramerization loop comprising two beta strands are swapped between the protomers. Key residues in StSurE that might promote C-terminal helix swapping were identified by sequence and structural comparisons. Three mutants in which the helix swapping is likely to be avoided were constructed and expressed in E. coli. Three-dimensional X-ray crystal structures of the mutants H234A and D230A/H234A could be determined at 2.1 angstrom and 2.35 angstrom resolutions, respectively. Contrary to expectations, helix swapping was mostly retained in both the mutants. The loss of the crucial D230 OD2- H234 NE2 hydrogen bond (2.89 angstrom in the wild type structure) in the hinge region was compensated by new inter and intra-chain interactions. However, the two fold molecular symmetry was lost and there were large conformational changes throughout the polypeptide. In spite of these changes, the dimeric structure and an approximate tetrameric organization were retained, probably due to the interactions involving the tetramerization loop. Mutants were mostly functionally inactive, highlighting the importance of precise inter-subunit interactions for the symmetry and function of StSurE.

Mechanistic features of Salmonella typhimurium propionate kinase (TdcD): Insights from kinetic and crystallographic studies

Short-chain fatty acids (SCFAs) play a major role in carbon cycle and can be utilized as a source of carbon and energy by bacteria. Salmonella typhimurium propionate kinase (StTdcD) catalyzes reversible transfer of the gamma-phosphate of ATP to propionate during L-threonine degradation to propionate. Kinetic analysis revealed that StTdcD possesses broad ligand specificity and could be activated by various SCFAs (propionate > acetate approximate to butyrate), nucleotides (ATP approximate to GTP > CTP approximate to TTP; dATP > dGTP > dCTP) and metal ions (Mg2+ approximate to Mn2+ > Co2+). Inhibition of StTdcD by tricarboxylic acid (TCA) cycle intermediates such as citrate, succinate, alpha-ketoglutarate and malate suggests that the enzyme could be under plausible feedback regulation. Crystal structures of StTdcD bound to PO4 (phosphate), AMP, ATP, Ap4 (adenosine tetraphosphate), GMP, GDP, GTP, CMP and CTP revealed that binding of nucleotide mainly involves hydrophobic interactions with the base moiety and could account for the broad biochemical specificity observed between the enzyme and nucleotides. Modeling and site-directed mutagenesis studies suggest Ala88 to be an important residue involved in determining the rate of catalysis with SCFA substrates. Molecular dynamics simulations on monomeric and dimeric forms of StTdcD revealed plausible open and closed states, and also suggested role for dimerization in stabilizing segment 235-290 involved in interfacial interactions and ligand binding. Observation of an ethylene glycol molecule bound sufficiently close to the gamma-phosphate in StTdcD complexes with triphosphate nucleotides supports direct in-line phosphoryl transfer. (C) 2013 Elsevier B.V. All rights reserved.

G-Quadruplex Structures Formed at the HOX11 Breakpoint Region Contribute to Its Fragility during t(10;14) Translocation in T-Cell Leukemia

The t(10;14) translocation involving the HOX11 gene is found in several T-cell leukemia patients. Previous efforts to determine the causes of HOX11 fragility were not successful. The role of non-B DNA structures is increasingly becoming an important cause of genomic instability. In the present study, bioinformatics analysis revealed two G-quadruplex-forming motifs at the HOX11 breakpoint cluster. Gel shift assays showed formation of both intra- and intermolecular G-quadruplexes, the latter being more predominant. The structure formation was dependent on four stretches of guanines, as revealed by mutagenesis. Circular dichroism analysis identified parallel conformations for both quadruplexes. The non-B DNA structure could block polymerization during replication on a plasmid, resulting in consistent K K+-dependent pause sites, which were abolished upon mutation of G-motifs, thereby demonstrating the role of the stretches of guanines even on double-stranded DNA. Extrachromosomal assays showed that the G-quadruplex motifs could block transcription, leading to reduced expression of green fluorescent protein (GFP) within cells. More importantly, sodium bisulfite modification assay showed the single-stranded character at regions I and II of HOX11 in the genome. Thus, our findings suggest the occurrence of G-quadruplex structures at the HOX11 breakpoint region, which could explain its fragility during the t(10;14) translocation.

A multi-level multi-scale approach to study essential genes in Mycobacterium tuberculosis

Background: The set of indispensable genes that are required by an organism to grow and sustain life are termed as essential genes. There is a strong interest in identification of the set of essential genes, particularly in pathogens, not only for a better understanding of the pathogen biology, but also for identifying drug targets and the minimal gene set for the organism. Essentiality is inherently a systems property and requires consideration of the system as a whole for their identification. The available experimental approaches capture some aspects but each method comes with its own limitations. Moreover, they do not explain the basis for essentiality in most cases. A powerful prediction method to recognize this gene pool including rationalization of the known essential genes in a given organism would be very useful. Here we describe a multi-level multi-scale approach to identify the essential gene pool in a deadly pathogen, Mycobacterium tuberculosis. Results: The multi-level workflow analyses the bacterial cell by studying (a) genome-wide gene expression profiles to identify the set of genes which show consistent and significant levels of expression in multiple samples of the same condition, (b) indispensability for growth by using gene expression integrated flux balance analysis of a genome-scale metabolic model, (c) importance for maintaining the integrity and flow in a protein-protein interaction network and (d) evolutionary conservation in a set of genomes of the same ecological niche. In the gene pool identified, the functional basis for essentiality has been addressed by studying residue level conservation and the sub-structure at the ligand binding pockets, from which essential amino acid residues in that pocket have also been identified. 283 genes were identified as essential genes with high-confidence. An agreement of about 73.5% is observed with that obtained from the experimental transposon mutagenesis technique. A large proportion of the identified genes belong to the class of intermediary metabolism and respiration. Conclusions: The multi-scale, multi-level approach described can be generally applied to other pathogens as well. The essential gene pool identified form a basis for designing experiments to probe their finer functional roles and also serve as a ready shortlist for identifying drug targets.

Specific Sequence of a Beta Turn in Human La Protein May Contribute to Species Specificity of Hepatitis C Virus

Human La protein is known to be an essential host factor for translation and replication of hepatitis C virus (HCV) RNA. Previously, we have demonstrated that residues responsible for interaction of human La protein with the HCV internal ribosomal entry site (IRES) around the initiator AUG within stem-loop IV form a beta-turn in the RNA recognition motif (RRM) structure. In this study, sequence alignment and mutagenesis suggest that the HCV RNA-interacting beta-turn is conserved only in humans and chimpanzees, the species primarily known to be infected by HCV. A 7-mer peptide corresponding to the HCV RNA-interacting region of human La inhibits HCV translation, whereas another peptide corresponding to the mouse La sequence was unable to do so. Furthermore, IRES-mediated translation was found to be significantly high in the presence of recombinant human La protein in vitro in rabbit reticulocyte lysate. We observed enhanced replication with HCV subgenomic and full-length replicons upon overexpression of either human La protein or a chimeric mouse La protein harboring a human La beta-turn sequence in mouse cells. Taken together, our results raise the possibility of creating an immunocompetent HCV mouse model using human-specific cell entry factors and a humanized form of La protein.

A Histidine Aspartate Ionic Lock Gates the Iron Passage in Miniferritins from Mycobacterium smegmatis

Background: DNA-binding protein from starved cells (Dps) are nano-compartments that can oxidize and store iron rendering protection from free radicals. Results: A histidine-aspartate ionic cluster in mycobaterial Dps2 modulates the rate of iron entry and exit in these proteins. Conclusion: Substitutions that disrupt the cluster interface alter the iron uptake/release properties with localized structural changes. Significance: Identifying important gating residues can help in designing nano-delivery vehicles. Dps (DNA-binding protein from starved cells) are dodecameric assemblies belonging to the ferritin family that can bind DNA, carry out ferroxidation, and store iron in their shells. The ferritin-like trimeric pore harbors the channel for the entry and exit of iron. By representing the structure of Dps as a network we have identified a charge-driven interface formed by a histidine aspartate cluster at the pore interface unique to Mycobacterium smegmatis Dps protein, MsDps2. Site-directed mutagenesis was employed to generate mutants to disrupt the charged interactions. Kinetics of iron uptake/release of the wild type and mutants were compared. Crystal structures were solved at a resolution of 1.8-2.2 for the various mutants to compare structural alterations vis a vis the wild type protein. The substitutions at the pore interface resulted in alterations in the side chain conformations leading to an overall weakening of the interface network, especially in cases of substitutions that alter the charge at the pore interface. Contrary to earlier findings where conserved aspartate residues were found crucial for iron release, we propose here that in the case of MsDps2, it is the interplay of negative-positive potentials at the pore that enables proper functioning of the protein. In similar studies in ferritins, negative and positive patches near the iron exit pore were found to be important in iron uptake/release kinetics. The unique ionic cluster in MsDps2 makes it a suitable candidate to act as nano-delivery vehicle, as these gated pores can be manipulated to exhibit conformations allowing for slow or fast rates of iron release.

TSpred: a web server for the rational design of temperature-sensitive mutants

Temperature sensitive (Ts) mutants of proteins provide experimentalists with a powerful and reversible way of conditionally expressing genes. The technique has been widely used in determining the role of gene and gene products in several cellular processes. Traditionally, Ts mutants are generated by random mutagenesis and then selected though laborious large-scale screening. Our web server, TSpred (http://mspc.bii.a-star.edu.sg/TSpred/), now enables users to rationally design Ts mutants for their proteins of interest. TSpred uses hydrophobicity and hydrophobic moment, deduced from primary sequence and residue depth, inferred from 3D structures to predict/identify buried hydrophobic residues. Mutating these residues leads to the creation of Ts mutants. Our method has been experimentally validated in 36 positions in six different proteins. It is an attractive proposition for Ts mutant engineering as it proposes a small number of mutations and with high precision. The accompanying web server is simple and intuitive to use and can handle proteins and protein complexes of different sizes.