Role of a PAS sensor domain in the Mycobacterium tuberculosis transcription regulator Rv1364c

The Mycobacterium tuberculosis transcriptional regulator Rv1364c regulates the activity of the stress response sigma factor sigma(F). This multi-domain protein has several components: a signaling PAS domain and an effector segment comprising of a phosphatase, a kinase and an anti-anti-sigma factor domain. Based on Small Angle X-ray Scattering (SAXS) data, Rv1364c was recently shown to be a homo-dimer and adopt an elongated conformation in solution. The PAS domain could not be modeled into the structural envelope due to poor sequence similarity with known PAS proteins. The crystal structure of the PAS domain described here provides a structural basis for the dimerization of Rv1364c. It thus appears likely that the PAS domain regulates the anti-sigma activity of Rv1364c by oligomerization. A structural comparison with other characterized PAS domains reveal several sequence and conformational features that could facilitate ligand binding - a feature which suggests that the function of Rv1364c could potentially be governed by specific cellular signals or metabolic cues. (C) 2010 Elsevier Inc. All rights reserved.

Crystal structures, dynamics and functional implications of molybdenum-cofactor biosynthesis protein MogA from two thermophilic organisms

Molybdenum-cofactor (Moco) biosynthesis is an evolutionarily conserved pathway in almost all kingdoms of life, including humans. Two proteins, MogA and MoeA, catalyze the last step of this pathway in bacteria, whereas a single two-domain protein carries out catalysis in eukaryotes. Here, three crystal structures of the Moco-biosynthesis protein MogA from the two thermophilic organisms Thermus thermophilus (TtMogA; 1.64 angstrom resolution, space group P2(1)) and Aquifex aeolicus (AaMogA; 1.70 angstrom resolution, space group P2(1) and 1.90 angstrom resolution, space group P1) have been determined. The functional roles and the residues involved in oligomerization of the protein molecules have been identified based on a comparative analysis of these structures with those of homologous proteins. Furthermore, functional roles have been proposed for the N- and C-terminal residues. In addition, a possible protein-protein complex of MogA and MoeA has been proposed and the residues involved in protein-protein interactions are discussed. Several invariant water molecules and those present at the subunit interfaces have been identified and their possible structural and/or functional roles are described in brief. In addition, molecular-dynamics and docking studies with several small molecules (including the substrate and the product) have been carried out in order to estimate their binding affinities towards AaMogA and TtMogA. The results obtained are further compared with those obtained for homologous eukaryotic proteins.

Denaturant mediated unfolding of both native and molten globule states of maltose binding protein are accompanied by large [Delta]Cp's

Maltose binding protein (MBP) is a large, monomeric two domain protein containing 370 amino acids. In the absence of denaturant at neutral pH, the protein is in the native state, while at pH 3.0 it forms a molten globule. The molten globule lacks a tertiary circular dichroism signal but has secondary structure similar to that of the native state. The molten globule binds 8-anilino-1-naphthalene sulfonate (ANS). The unfolding thermodynamics of MBP at both pHs were measured by carrying out a series of isothermal urea melts at temperatures ranging from 274–329 K. At 298 K, values of [Delta]G°, [Delta]Cp, and Cm were 3.1 ± 0.2 kcal mol−1, 5.9 ± 0.8 kcal mol−1 K−1 (15.9 cal (mol-residue)−1 K−1), and 0.8 M, respectively, at pH 3.0 and 14.5 ± 0.4 kcal mol−1, 8.3 ± 0.7 kcal mol−1 K−1 (22.4 kcal (mol-residue)−1 K−1), and 3.3 M, respectively, at pH 7.1. Guanidine hydrochloride denaturation at pH 7.1 gave values of [Delta]G° and [Delta]Cp similar to those obtained with urea. The m values for denaturation are strongly temperature dependent, in contrast to what has been previously observed for small globular proteins. The value of [Delta]Cp per mol-residue for the molten globule is comparable to corresponding values of [Delta]Cp for the unfolding of typical globular proteins and suggests that it is a highly ordered structure, unlike molten globules of many small proteins. The value of [Delta]Cp per mol-residue for the unfolding of the native state is among the highest currently known for any protein.

Mycobacterium leprae RecA is structurally analogous but functionally distinct from Mycobacterium tuberculosis RecA protein

Mycobacterium leprae is closely related to Mycobacterium tuberculosis, yet causes a very different illness. Detailed genomic comparison between these two species of mycobacteria reveals that the decaying M. leprae genome contains less than half of the M. tuberculosis functional genes. The reduction of genome size and accumulation of pseudogenes in the M. leprae genome is thought to result from multiple recombination events between related repetitive sequences, which provided the impetus to investigate the recombination-like activities of RecA protein. In this study, we have cloned, over-expressed and purified M. leprae RecA and compared its activities with that of M. tuberculosis RecA. Both proteins, despite being 91% identical at the amino acid level, exhibit strikingly different binding profiles for single-stranded DNA with varying GC contents, in the ability to catalyze the formation of D-loops and to promote DNA strand exchange. The kinetics and the extent of single-stranded DNA-dependent ATPase and coprotease activities were nearly equivalent between these two recombinases. However, the degree of inhibition exerted by a range of ATP:ADP ratios was greater on strand exchange promoted by M. leprae RecA compared to its M. tuberculosis counterpart. Taken together, our results provide insights into the mechanistic aspects of homologous recombination and coprotease activity promoted by M. lepare RecA, and further suggests that it differs from the M. tuberculosis counterpart. These results are consistent with an emerging concept of DNA-sequence influenced structural differences in RecA nucleoprotein filaments and how these differences reflect on the multiple activities associated with RecA protein. (C) 2011 Elsevier B.V. All rights reserved.

Tethering preferences of domain families co-occurring in multi-domain proteins

Genomic data of several organisms have revealed the presence of a vast repertoire of multi-domain proteins. The role played by individual domains in a multi-domain protein has a profound influence on the overall function of the protein. In the present analysis an attempt has been made to better understand the tethering preferences of domain families that occur in multi-domain proteins. The analysis has been carried out on an exhaustive dataset of 2 961 898 sequences of proteins from 930 organisms, where 741 274 proteins are comprised of at least two domain families. For every domain family, the number of other domain families with which it co-occurs within a protein in this dataset has been enumerated and is referred to as the tethering number of the domain family. It was found that, in the general dataset, the AAA ATPase family and the family of Ser/Thr kinases have the highest tethering numbers of 450 and 444 respectively. Further analysis reveals significant correlation between the number of members in a family and its tethering number. Positive correlation was also observed for the extent of a sequence and functional diversity within a family and the tethering numbers of domain families. Domain families that are present ubiquitously in diverse organisms tend to have large tethering numbers, while organism/kingdom-specific families have low tethering numbers. Thus, the analysis uncovers how domain families recombine and evolve to give rise to multi-domain proteins.

The relationship between classification of multi-domain proteins using an alignment-free approach and their functions: a case study with immunoglobulins

Establishing functional relationships between multi-domain protein sequences is a non-trivial task. Traditionally, delineating functional assignment and relationships of proteins requires domain assignments as a prerequisite. This process is sensitive to alignment quality and domain definitions. In multi-domain proteins due to multiple reasons, the quality of alignments is poor. We report the correspondence between the classification of proteins represented as full-length gene products and their functions. Our approach differs fundamentally from traditional methods in not performing the classification at the level of domains. Our method is based on an alignment free local matching scores (LMS) computation at the amino-acid sequence level followed by hierarchical clustering. As there are no gold standards for full-length protein sequence classification, we resorted to Gene Ontology and domain-architecture based similarity measures to assess our classification. The final clusters obtained using LMS show high functional and domain architectural similarities. Comparison of the current method with alignment based approaches at both domain and full-length protein showed superiority of the LMS scores. Using this method we have recreated objective relationships among different protein kinase sub-families and also classified immunoglobulin containing proteins where sub-family definitions do not exist currently. This method can be applied to any set of protein sequences and hence will be instrumental in analysis of large numbers of full-length protein sequences.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

A Gaussian network model study suggests that structural fluctuations are higher for inactive states than active states of protein kinases

We performed Gaussian network model based normal mode analysis of 3-dimensional structures of multiple active and inactive forms of protein kinases. In 14 different kinases, a more number of residues (1095) show higher structural fluctuations in inactive states than those in active states (525), suggesting that, in general, mobility of inactive states is higher than active states. This statistically significant difference is consistent with higher crystallographic B-factors and conformational energies for inactive than active states, suggesting lower stability of inactive forms. Only a small number of inactive conformations with the DFG motif in the ``in'' state were found to have fluctuation magnitudes comparable to the active conformation. Therefore our study reports for the first time, intrinsic higher structural fluctuation for almost all inactive conformations compared to the active forms. Regions with higher fluctuations in the inactive states are often localized to the aC-helix, aG-helix and activation loop which are involved in the regulation and/or in structural transitions between active and inactive states. Further analysis of 476 kinase structures involved in interactions with another domain/protein showed that many of the regions with higher inactive-state fluctuation correspond to contact interfaces. We also performed extensive GNM analysis of (i) insulin receptor kinase bound to another protein and (ii) holo and apo forms of active and inactive conformations followed by multi-factor analysis of variance. We conclude that binding of small molecules or other domains/proteins reduce the extent of fluctuation irrespective of active or inactive forms. Finally, we show that the perceived fluctuations serve as a useful input to predict the functional state of a kinase.

Proteomic analysis of the Cdc48/Ubx network identifies a role for Ubx2 in the regulation of lipid biosynthesis

Cdc48/p97 is an essential, highly abundant hexameric member of the AAA (ATPase associated with various cellular activities) family. It has been linked to a variety of processes throughout the cell but it is best known for its role in the ubiquitin proteasome pathway. In this system it is believed that Cdc48 behaves as a segregase, transducing the chemical energy of ATP hydrolysis into mechanical force to separate ubiquitin-conjugated proteins from their tightly-bound partners.

Current models posit that Cdc48 is linked to its substrates through a variety of adaptor proteins, including a family of seven proteins (13 in humans) that contain a Cdc48-binding UBX domain. As such, due to the complexity of the network of adaptor proteins for which it serves as the hub, Cdc48/p97 has the potential to exert a profound influence on the ubiquitin proteasome pathway. However, the number of known substrates of Cdc48/p97 remains relatively small, and smaller still is the number of substrates that have been linked to a specific UBX domain protein. As such, the goal of this dissertation research has been to discover new substrates and better understand the functions of the Cdc48 network. With this objective in mind, we established a proteomic screen to assemble a catalog of candidate substrate/targets of the Ubx adaptor system.

Here we describe the implementation and optimization of a cutting-edge quantitative mass spectrometry method to measure relative changes in the Saccharomyces cerevisiae proteome. Utilizing this technology, and in order to better understand the breadth of function of Cdc48 and its adaptors, we then performed a global screen to identify accumulating ubiquitin conjugates in cdc48-3 and ubxΔ mutants. In this screen different ubx mutants exhibited reproducible patterns of conjugate accumulation that differed greatly from each other, pointing to various unexpected functional specializations of the individual Ubx proteins.

As validation of our mass spectrometry findings, we then examined in detail the endoplasmic-reticulum bound transcription factor Spt23, which we identified as a putative Ubx2 substrate. In these studies ubx2Δ cells were deficient in processing of Spt23 to its active p90 form, and in localizing p90 to the nucleus. Additionally, consistent with reduced processing of Spt23, ubx2Δ cells demonstrated a defect in expression of their target gene OLE1, a fatty acid desaturase. Overall, this work demonstrates the power of proteomics as a tool to identify new targets of various pathways and reveals Ubx2 as a key regulator lipid membrane biosynthesis.

First molluscan TNFR homologue in Zhikong scallop: Molecular characterization and expression analysis

Tumor necrosis factor receptors (TNFRs) are a superfamily of proteins characterized by the unique cysteine-rich domain (CRD) and their important roles in diverse physiological and pathological events such as inflammation, apoptosis, autoimmunity and organogenesis. The first member of the molluscan TNFR family, designated as CfTNFR, was identified from Zhikong scallop Chlamys farreri by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of CfTNFR was of 1334 bp, consisting of a 5' UTR of 17 bp, a 3'UTR of 69 by with a poly (A) tail, and an open reading frame (ORE) of 1248 by encoding a polypeptide of 415 amino acids with a theoretical isoelectric point of 8.33 and predicted molecular weight of 47.07 kDa. There were a signal peptide, a CRD, a transmembrane region and a death domain in the deduced amino acid sequence of CfTNFR, suggesting that it was a typical type 1 membrane protein. The high identities (22-40%) of CfTNFR with other TNFR superfamily members indicated that CfTNFR should be a member of TNFR superfamily, and moreover, it should be the first death domain-containing TNFR found in invertebrates. Phylogenetic analysis revealed that CfTNFR was closely related to TNFR-like proteins from Strongylocentrotus purpuratus, Drosophila melanogaster and Ciona intestinalis, and they formed a separate branch apart from vertebrate TNFRs. The spatial expression of CfTNFR transcripts in healthy and bacteria challenged scallops was examined by quantitative real-time PCR. CfTNFR transcripts could be detected in all tested tissues, including haemocytes, gonad, gill, mantle and hepatopancreas, and significantly up-regulated in the tissues of gonad, gill, mantle and hepatopancreas after Listonella anguillarum challenge, indicating that CfTNFR was constitutive and inducible acute-phase protein involved in immune defence. The present results suggested the existence of the TNFR-like molecules and TNF-TNFR system in low invertebrates, and provided new insights into the role of CfTNFR in scallop innate immune responses to invading microorganisms. (C) 2009 Elsevier Ltd. All rights reserved.

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

CONTRIBUTION OF A SPERM PROTEIN, PAWP, TO THE SIGNAL TRANSDUCT PATHWAY DURING VERTEBRATE FERTILIZATION

PAWP, postacrosomal sheath WW domain binding protein, is a novel sperm protein identified as a candidate sperm borne, oocyte-activating factor (SOAF). PAWP induces both early and later egg activation events including meiotic resumption, pronuclear formation and egg cleavage. Based on the fact that calcium increase is universally accepted as the sole requirement for egg activation, we hypothesized that PAWP is an upstream regulator of the calcium signaling pathway during fertilization. Intracellular calcium increase was detected by two-photon laser scanning fluorescence microscopy following microinjection of recombinant PAWP into Xenopus oocytes, bolstering our hypothesis and suggesting the involvement of a novel PAWP-mediated signaling pathway during fertilization. The N-terminal of PAWP shares a high homology to WW domain binding protein while the C-terminal half contains a functional PPXY motif, which allows it to interact with group I WW domain proteins. These structural considerations together with published data indicating that PPXY synthetic peptide derived from PAWP inhibits ICSI-induced fertilization led to the hypothesis that PAWP triggers egg activation by binding to a group I WW domain protein in the oocyte. By far-Western analysis of oocyte cytoplasmic fraction, PAWP was found to bind to a 52 kDa protein. The competitive inhibition studies with PPXY synthetic peptide, WW domain constructs, and their point mutants demonstrated that the interaction between PAWP and its binding partner is specifically via the PPXY-WW domain module. The 52 kDa protein band crossreacted with antibodies against group I WW domain protein YAP in Western blot assay, indicating that this 52 kDa PAWP binding partner is either YAP or a YAP-related protein. In addition, the far-Western competitive inhibition studies with recombinant GST fusion protein YAP and another WW domain-containing protein, TAZ, demonstrated that the binding of PAWP to its binding partner was significantly reduced by TAZ, providing evidence that TAZ could be the 52 kDa protein candidate. Mass spectrometry was employed to identify this PAWP binding partner candidate. However, due to the low abundance of the candidate protein and the complexity of the sample, several strategies are still needed to enrich this protein. This study correlates PAWP induced meiotic resumption and calcium efflux at fertilization and uncovers a 52 kDa candidate WW domain protein in the oocyte cytoplasm that most likely interacts with PAWP to trigger egg activation.

Genetic variation in MME in relation to neprilysin protein and enzyme activity, Aβ levels, and Alzheimer's disease risk

Neprilysin (NEP), also known as membrane metalloendopeptidase (MME), is considered amongst the most important ß-amyloid (Aß)-degrading enzymes with regard to prevention of Alzheimer's disease (AD) pathology. Variation in the NEP gene (MME) has been suggested as a risk factor for AD. We conducted a genetic association study of 7MME SNPs - rs1836914, rs989692, rs9827586, rs6797911, rs61760379, rs3736187, rs701109 - with respect to AD risk in a cohort of 1057 probable and confirmed AD cases and 424 age-matched non-demented controls from the United Kingdom, Italy and Sweden. We also examined the association of these MME SNPs with NEP protein level and enzyme activity, and on biochemical measures of Aß accumulation in frontal cortex - levels of total soluble Aß, oligomeric Aß(1-42), and guanidine-extractable (insoluble) Aß - in a sub-group of AD and control cases with post-mortem brain tissue. On multivariate logistic regression analysis one of the MME variants (rs6797911) was associated with AD risk (P = 0.00052, Odds Ratio (O.R. = 1.40, 95% confidence interval (1.16-1.70)). None of the SNPs had any association with Aß levels; however, rs9827586 was significantly associated with NEP protein level (p=0.014) and enzyme activity (p=0.006). Association was also found between rs701109 and NEP protein level (p=0.026) and a marginally non-significant association was found for rs989692 (p=0.055). These data suggest that MME variation may be associated with AD risk but we have not found evidence that this is mediated through modification of NEP protein level or activity.

Parsing parallel evolution:Ecological divergence and differential gene expression in the adaptive radiations of thick-lipped Midas cichlid fishes from Nicaragua

The study of parallel evolution facilitates the discovery of common rules of diversification. Here, we examine the repeated evolution of thick lips in Midas cichlid fishes (the Amphilophus citrinellus species complex) - from two Great Lakes and two crater lakes in Nicaragua - to assess whether similar changes in ecology, phenotypic trophic traits and gene expression accompany parallel trait evolution. Using next-generation sequencing technology, we characterize transcriptome-wide differential gene expression in the lips of wild-caught sympatric thick- and thin-lipped cichlids from all four instances of repeated thick-lip evolution. Six genes (apolipoprotein D, myelin-associated glycoprotein precursor, four-and-a-half LIM domain protein 2, calpain-9, GTPase IMAP family member 8-like and one hypothetical protein) are significantly underexpressed in the thick-lipped morph across all four lakes. However, other aspects of lips' gene expression in sympatric morphs differ in a lake-specific pattern, including the magnitude of differentially expressed genes (97-510). Generally, fewer genes are differentially expressed among morphs in the younger crater lakes than in those from the older Great Lakes. Body shape, lower pharyngeal jaw size and shape, and stable isotopes (dC and dN) differ between all sympatric morphs, with the greatest differentiation in the Great Lake Nicaragua. Some ecological traits evolve in parallel (those related to foraging ecology; e.g. lip size, body and head shape) but others, somewhat surprisingly, do not (those related to diet and food processing; e.g. jaw size and shape, stable isotopes). Taken together, this case of parallelism among thick- and thin-lipped cichlids shows a mosaic pattern of parallel and nonparallel evolution. © 2012 Blackwell Publishing Ltd.