Phylogenetic relationships and classification of thiolases and thiolase-like proteins of Mycobacterium tuberculosis and Mycobacterium smegmatis

Thiolases are enzymes involved in lipid metabolism. Thiolases remove the acetyl-CoA moiety from 3-ketoacyl-CoAs in the degradative reaction. They can also catalyze the reverse Claisen condensation reaction, which is the first step of biosynthetic processes such as the biosynthesis of sterols and ketone bodies. In human, six distinct thiolases have been identified. Each of these thiolases is different from the other with respect to sequence, oligomeric state, substrate specificity and subcellular localization. Four sequence fingerprints, identifying catalytic loops of thiolases, have been described. In this study genome searches of two mycobacterial species (Mycobacterium tuberculosis and Mycobacterium smegmatis), were carried out, using the six human thiolase sequences as queries. Eight and thirteen different thiolase sequences were identified in M. tuberculosis and M. smegmatis, respectively. In addition, thiolase-like proteins (one encoded in the Mtb and two in the Msm genome) were found. The purpose of this study is to classify these mostly uncharacterized thiolases and thiolase-like proteins. Several other sequences obtained by searches of genome databases of bacteria, mammals and the parasitic protist family of the Trypanosomatidae were included in the analysis. Thiolase-like proteins were also found in the trypanosomatid genomes, but not in those of mammals. In order to study the phylogenetic relationships at a high confidence level, additional thiolase sequences were included such that a total of 130 thiolases and thiolase-like protein sequences were used for the multiple sequence alignment. The resulting phylogenetic tree identifies 12 classes of sequences, each possessing a characteristic set of sequence fingerprints for the catalytic loops. From this analysis it is now possible to assign the mycobacterial thiolases to corresponding homologues in other kingdoms of life. The results of this bioinformatics analysis also show interesting differences between the distributions of M. tuberculosis and M. smegmatis thiolases over the 12 different classes. (C) 2014 Elsevier Ltd. All rights reserved.

The Evolution of Guanylyl Cyclases as Multidomain Proteins: Conserved Features of Kinase-Cyclase Domain Fusions

Guanylyl cyclases (GCs) are enzymes that generate cyclic GMP and regulate different physiologic and developmental processes in a number of organisms. GCs possess sequence similarity to class III adenylyl cyclases (ACs) and are present as either membrane-bound receptor GCs or cytosolic soluble GCs. We sought to determine the evolution of GCs using a large-scale bioinformatic analysis and found multiple lineage-specific expansions of GC genes in the genomes of many eukaryotes. Moreover, a few GC-like proteins were identified in prokaryotes, which come fused to a number of different domains, suggesting allosteric regulation of nucleotide cyclase activity Eukaryotic receptor GCs are associated with a kinase homology domain (KHD), and phylogenetic analysis of these proteins suggest coevolution of the KHD and the associated cyclase domain as well as a conservation of the sequence and the size of the linker region between the KHD and the associated cyclase domain. Finally, we also report the existence of mimiviral proteins that contain putative active kinase domains associated with a cyclase domain, which could suggest early evolution of the fusion of these two important domains involved in signa transduction.

Crystal Structure of a Monomeric Thiolase-Like Protein Type 1 (TLP1) from Mycobacterium smegmatis

An analysis of the Mycobacterium smegmatis genome suggests that it codes for several thiolases and thiolase-like proteins. Thiolases are an important family of enzymes that are involved in fatty acid metabolism. They occur as either dimers or tetramers. Thiolases catalyze the Claisen condensation of two acetyl-Coenzyme A molecules in the synthetic direction and the thiolytic cleavage of 3-ketoacyl-Coenzyme A molecules in the degradative direction. Some of the M. smegmatis genes have been annotated as thiolases of the poorly characterized SCP2-thiolase subfamily. The mammalian SCP2-thiolase consists of an N-terminal thiolase domain followed by an additional C-terminal domain called sterol carrier protein-2 or SCP2. The M. smegmatis protein selected in the present study, referred to here as the thiolase-like protein type 1 (MsTLP1), has been biochemically and structurally characterized. Unlike classical thiolases, MsTLP1 is a monomer in solution. Its structure has been determined at 2.7 angstrom resolution by the single wavelength anomalous dispersion method. The structure of the protomer confirms that the N-terminal domain has the thiolase fold. An extra C-terminal domain is indeed observed. Interestingly, it consists of six beta-strands forming an anti-parallel beta-barrel which is completely different from the expected SCP2-fold. Detailed sequence and structural comparisons with thiolases show that the residues known to be essential for catalysis are not conserved in MsTLP1. Consistent with this observation, activity measurements show that MsTLP1 does not catalyze the thiolase reaction. This is the first structural report of a monomeric thiolase-like protein from any organism. These studies show that MsTLP1 belongs to a new group of thiolase related proteins of unknown function.

At4g24160, a Soluble Acyl-Coenzyme A-Dependent Lysophosphatidic Acid Acyltransferase

Human CGI-58 (for comparative gene identification-58) and YLR099c, encoding Ict1p in Saccharomyces cerevisiae, have recently been identified as acyl-CoA-dependent lysophosphatidic acid acyltransferases. Sequence database searches for CGI-58 like proteins in Arabidopsis (Arabidopsis thaliana) revealed 24 proteins with At4g24160, a member of the alpha/beta-hydrolase family of proteins being the closest homolog. At4g24160 contains three motifs that are conserved across the plant species: a GXSXG lipase motif, a HX4D acyltransferase motif, and V(X)(3)HGF, a probable lipid binding motif. Dendrogram analysis of yeast ICT1, CGI-58, and At4g24160 placed these three polypeptides in the same group. Here, we describe and characterize At4g24160 as, to our knowledge, the first soluble lysophosphatidic acid acyltransferase in plants. A lipidomics approach revealed that At4g24160 has additional triacylglycerol lipase and phosphatidylcholine hydrolyzing enzymatic activities. These data establish At4g24160, a protein with a previously unknown function, as an enzyme that might play a pivotal role in maintaining the lipid homeostasis in plants by regulating both phospholipid and neutral lipid levels.

Characterization of DNA Strand Transfer Promoted by Mycobacterium smegmatis RecA Reveals Functional Diversity with Mycobacterium tuberculosis RecA†

The RecA-like proteins constitute a group of DNA strand transfer proteins ubiquitous in eubacteria, eukarya, and archaea. However, the functional relationship among RecA proteins is poorly understood. For instance, Mycobacterium tuberculosis RecA is synthesized as a large precursor, which undergoes an unusual protein-splicing reaction to generate an active form. Whereas the precursor was inactive, the active form promoted DNA strand transfer less efficiently compared to EcRecA. Furthermore, gene disruption studies have indicated that the frequencies of allele exchange are relatively lower in Mycobacterium tuberculosis compared to Mycobacterium smegmatis. The mechanistic basis and the factors that contribute to differences in allele exchange remain to be understood. Here, we show that the extent of DNA strand transfer promoted by the M. smegmatis RecA in vitro differs significantly from that of M. tuberculosis RecA. Importantly, M. smegmatis RecA by itself was unable to promote strand transfer, but cognate or noncognate SSBs rendered it efficient even when added prior to RecA. In the presence of SSB, MsRecA or MtRecA catalyzed strand transfer between ssDNA and varying lengths of linear duplex DNA with distinctly different pH profiles. The factors that were able to suppress the formation of DNA networks greatly stimulated strand transfer reactions promoted by MsRecA or MtRecA. Although the rate and pH profiles of dATP hydrolysis catalyzed by MtRecA and MsRecA were similar, only MsRecA was able to couple dATP hydrolysis to DNA strand transfer. Together, these results provide insights into the functional diversity in DNA strand transfer promoted by RecA proteins of pathogenic and nonpathogenic species of mycobacteria.

Structure of the putative mutarotase YeaD from Salmonella typhimurium in two different forms: in vivo binding of a sugar phosphate

Salmonella typhimurium YeaD (stYeaD), annotated as a putative aldose 1-epimerase, has a very low sequence identity to other well characterized mutarotases. Sequence analysis suggested that the catalytic residues and a few of the substrate-binding residues of galactose mutarotases (GalMs) are conserved in stYeaD. Determination of the crystal structure of stYeaD in an orthorhombic form at 1.9 angstrom resolution and in a monoclinic form at 2.5 angstrom resolution revealed this protein to adopt the beta-sandwich fold similar to GalMs. Structural comparison of stYeaD with GalMs has permitted the identification of residues involved in catalysis and substrate binding. In spite of the similar fold and conservation of catalytic residues, minor but significant differences were observed in the substrate- binding pocket. These analyses pointed out the possible role of Arg74 and Arg99, found only in YeaD-like proteins, in ligand anchoring and suggested that the specificity of stYeaD may be distinct from those of GalMs

Crystal structure of Rv2118c: an AdoMet-dependent methyltransferase from Mycobacterium tuberculosis H37Rv

Rv2118c belongs to the class of conserved hypothetical proteins from Mycobacterium tuberculosis H37Rv. The crystal structure of Rv2118c in complex with S-adenosyl-Image -methionine (AdoMet) has been determined at 1.98 Å resolution. The crystallographic asymmetric unit consists of a monomer, but symmetry-related subunits interact extensively, leading to a tetrameric structure. The structure of the monomer can be divided functionally into two domains: the larger catalytic C-terminal domain that binds the cofactor AdoMet and is involved in the transfer of methyl group from AdoMet to the substrate and a smaller N-terminal domain. The structure of the catalytic domain is very similar to that of other AdoMet-dependent methyltransferases. The N-terminal domain is primarily a β-structure with a fold not found in other methyltransferases of known structure. Database searches reveal a conserved family of Rv2118c-like proteins from various organisms. Multiple sequence alignments show several regions of high sequence similarity (motifs) in this family of proteins. Structure analysis and homology to yeast Gcd14p suggest that Rv2118c could be an RNA methyltransferase, but further studies are required to establish its functional role conclusively.

Macromolecular Crystallography in India in the Global Context

The most spectacular applications of crystallography are currently concerned with biological macromolecules like proteins and their assemblies. Macromolecular crystallography originated in England in the thirties of the last century, but definitive results began to appear only around 1960. Since then macromolecular crystallography has grown to become central to modern biology. India has a long tradition in crystallography starting with the work of K. Banerjee in the thirties. In addition to their contributions to crystallography, G.N. Ramachandran and his colleagues gave a head start to India in computational biology, molecular modeling and what we now call bioinformatics. However, attempts to initiate macromolecular crystallography in India started only in the seventies. The work took off the ground after the Department of Science and Technology handsomely supported the group at Indian Institute of Science, Bangalore in 1983. The Bangalore group was also recognized as a national nucleus for the development of the area in the country. Since then macromolecular crystallography, practiced in more than 30 institutions in the country, has grown to become an important component of scientific research in India. The articles in this issue provide a flavor of activities in the area in the country. The area is still in an expanding phase and is poised to scale greater heights.

Spatio-temporal correlations in aqueous systems: computational studies of static and dynamic heterogeneity by 2D-IR spectroscopy

Local heterogeneity is ubiquitous in natural aqueous systems. It can be caused locally by external biomolecular subsystems like proteins, DNA, micelles and reverse micelles, nanoscopic materials etc., but can also be intrinsic to the thermodynamic nature of the aqueous solution itself (like binary mixtures or at the gas-liquid interface). The altered dynamics of water in the presence of such diverse surfaces has attracted considerable attention in recent years. As these interfaces are quite narrow, only a few molecular layers thick, they are hard to study by conventional methods. The recent development of two dimensional infra-red (2D-IR) spectroscopy allows us to estimate length and time scales of such dynamics fairly accurately. In this work, we present a series of interesting studies employing two dimensional infra-red spectroscopy (2D-IR) to investigate (i) the heterogeneous dynamics of water inside reverse micelles of varying sizes, (ii) supercritical water near the Widom line that is known to exhibit pronounced density fluctuations and also study (iii) the collective and local polarization fluctuation of water molecules in the presence of several different proteins. The spatio-temporal correlation of confined water molecules inside reverse micelles of varying sizes is well captured through the spectral diffusion of corresponding 2D-IR spectra. In the case of supercritical water also, we observe a strong signature of dynamic heterogeneity from the elongated nature of the 2D-IR spectra. In this case the relaxation is ultrafast. We find remarkable agreement between the different tools employed to study the relaxation of density heterogeneity. For aqueous protein solutions, we find that the calculated dielectric constant of the respective systems unanimously shows a noticeable increment compared to that of neat water. However, the `effective' dielectric constant for successive layers shows significant variation, with the layer adjacent to the protein having a much lower value. Relaxation is also slowest at the surface. We find that the dielectric constant achieves the bulk value at distances more than 3 nm from the surface of the protein.

A novel supersecondary structure in globular proteins comprising the collagen-like helix and ?-turn

A structure consisting of the polyproline-II or collagen-like helix immediately succeeded by a ?-turn is seen in several synthetic peptides and has been suggested to be the conformational requirement for proline hydroxylation in nascent procollagen. Using a simple algorithm for detecting secondary structures, we have analysed crystal structure data on 40 globular proteins and have found eight examples of the collagen-helix + ?-turn supersecondary structure in 15 proteins that contain the collagen-like helical segments.

SecB binds only to a late native-like intermediate in the folding pathway of barstar and not to the unfolded state

SecB is a cytosolic, tetrameric chaperone of Escherichia coli which maintains precursor proteins in a translocation competent state. We have investigated the effect of SecB on the refolding kinetics of the small protein barstar in I M guanidine hydrochloride at pH 7.0 and 25 degrees C using fluorescence spectroscopy. We show that SecB does not bind either the native or the unfolded states of barstar but binds to a late near-native intermediate along the folding pathway. For barstar, polypeptide collapse and formation of a hydrophobic surface are required for binding to SecB. SecB does not change the apparent rate constant of barstar refolding. The kinetic data for SecB binding to barstar are not consistent with simple kinetic partitioning models.

Nucleotide sequence of the 3' terminal region of belladonna mottle virus-Iowa (renamed Physalis mottle virus) RNA and an analysis of the relationships of tymoviral coat proteins

The 3prime terminal 1255nt sequence of Physalis mottle virus (PhMV) genomic RNA has been determined from a set of overlapping cDNA clones. The open reading frame (ORF) at the 3prime terminus corresponds to the amino acid sequence of the coat protein (CP) determined earlier except for the absence of the dipeptide, Lys-Leu, at position 110-111. In addition, the sequence upstream of the CP gene contains the message coding for 178 amino acid residues of the C-terminus of the putative replicase protein (RP). The sequence downstream of the CP gene contains an untranslated region whose terminal 80 nucleotides can be folded into a characteristic tRNA-like structure. A phylogenetic tree constructed after aligning separately the sequence of the CP, the replicase protein (RP) and the tRNA-like structure determined in this study with the corresponding sequences of other tymoviruses shows that PhMV wrongly named belladonna mottle virus [BDMV(I)] is a separate tymovirus and not another strain of BDMV(E) as originally envisaged. The phylogenetic tree in all the three cases is identical showing that any subset of genomic sequence of sufficient length can be used for establishing evolutionary relationships among tymoviruses.

Artifactual staining of proteins on polyacrylamide gels by nitrobluetetrazolium chloride and phenazine methosulfate

Different purified proteins were shown to give purple formazan bands corresponding to the protein stain following electrophoresis on polyacrylamide gels, in the presence of nitrobluetetrazolium (NBT) and phenazine methosulfate (PMS). Both PMS and NBT are needed for formazan production which has a favorable pH at 8.5. Sulfhydryl blockers in the incubation medium inhibited this color development to different extents. While proteins with free SH groups like bovine serum albumin, ovalbumin, and urease showed this pyridine nucleotide independent artifact, nonthiol proteins, viz., bovine pancreatic ribonuclease A, and riboflavin-binding protein from chicken egg white failed to do so. The nonenzymatic formazan formation observed with different proteins could also be shown in an in vitro assay system. It is clear that the “nothing dehydrogenase” phenomenon observed in several cases may be due to the thiol group-mediated artifactual staining of proteins.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of Delta C-p of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Molten globule-like state of peanut lectin monomer retains its carbohydrate specificity - Implications in protein folding and legume lectin oligomerization

A central question in biological chemistry is the minimal structural requirement of a protein that would determine its specificity and activity, the underlying basis being the importance of the entire structural element of a protein with regards to its activity vis a vis the overall integrity and stability of the protein. Although there are many reports on the characterization of protein folding/ unfolding intermediates, with considerable secondary structural elements but substantial loss of tertiary structure, none of them have been reported to show any activity toward their respective ligands. This may be a result of the conditions under which such intermediates have been isolated or due to the importance of specific structural elements for the activity. In this paper we report such an intermediate in the unfolding of peanut agglutinin that seems to retain, to a considerable degree, its carbohydrate binding specificity and activity. This result has significant implications on the molten globule state during the folding pathway(s) of proteins in general and the quaternary association in legume lectins in particular, where precise subunit topology is required for their biologic activities.