Discovering functionally important sites in proteins.

Proteins play important roles in biology, biotechnology and pharmacology, and missense variants are a common cause of disease. Discovering functionally important sites in proteins is a central but difficult problem because of the lack of large, systematic data sets. Sequence conservation can highlight residues that are functionally important but is often convoluted with a signal for preserving structural stability. We here present a machine learning method to predict functional sites by combining statistical models for protein sequences with biophysical models of stability. We train the model using multiplexed experimental data on variant effects and validate it broadly. We show how the model can be used to discover active sites, as well as regulatory and binding sites. We illustrate the utility of the model by prospective prediction and subsequent experimental validation on the functional consequences of missense variants in HPRT1 which may cause Lesch-Nyhan syndrome, and pinpoint the molecular mechanisms by which they cause disease.

Deciphering the Alphabet of Disorder-Glu and Asp Act Differently on Local but Not Global Properties.

Compared to folded proteins, the sequences of intrinsically disordered proteins (IDPs) are enriched in polar and charged amino acids. Glutamate is one of the most enriched amino acids in IDPs, while the chemically similar amino acid aspartate is less enriched. So far, the underlying functional differences between glutamates and aspartates in IDPs remain poorly understood. In this study, we examine the differential effects of aspartate and glutamates in IDPs by comparing the function and conformational ensemble of glutamate and aspartate variants of the disordered protein Dss1, using a range of assays, including interaction studies, nuclear magnetic resonance spectroscopy, small-angle X-ray scattering and molecular dynamics simulation. First, we analyze the sequences of the rapidly growing database of experimentally verified IDPs (DisProt) and show that glutamate enrichment is not caused by a taxonomy bias in IDPs. From analyses of local and global structural properties as well as cell growth and protein-protein interactions using a model acidic IDP from yeast and three Glu/Asp variants, we find that while the Glu/Asp variants support similar function and global dimensions, the variants differ in their binding affinities and population of local transient structural elements. We speculate that these local structural differences may play roles in functional diversity, where glutamates can support increased helicity, important for folding and binding, while aspartates support extended structures and form helical caps, as well as playing more relevant roles in, e.g., transactivation domains and ion-binding.

Disease-linked mutations cause exposure of a protein quality control degron.

More than half of disease-causing missense variants are thought to lead to protein degradation, but the molecular mechanism of how these variants are recognized by the cell remains enigmatic. Degrons are stretches of amino acids that help mediate recognition by E3 ligases and thus confer protein degradation via the ubiquitin-proteasome system. While degrons that mediate controlled degradation of, for example, signaling components and cell-cycle regulators are well described, so-called protein-quality-control degrons that mediate the degradation of destabilized proteins are poorly understood. Here, we show that disease-linked dihydrofolate reductase (DHFR) missense variants are structurally destabilized and chaperone-dependent proteasome targets. We find two regions in DHFR that act as degrons, and the proteasomal turnover of one of these was dependent on the molecular chaperone Hsp70. Structural analyses by nuclear magnetic resonance (NMR) and hydrogen/deuterium exchange revealed that this degron is buried in wild-type DHFR but becomes transiently exposed in the disease-linked missense variants.

UV oxidation of cyclic AMP receptor protein, a global bacterial gene regulator, decreases DNA binding and cleaves DNA at specific sites.

UV light is a widely-employed, and environmentally-sensitive bactericide but its mechanism of action is not fully defined. Proteins are major chromophores and targets for damage due to their abundance, but the role of proteins in inducing damage to bound DNA, and the effects on DNA-protein interactions is less well characterized. In E. coli (and other Gram-negative bacteria) the cyclic AMP receptor protein (CRP/CAP) regulates more than 500 genes. In this study we show that exposure of isolated dimeric CRP-cAMP to UV modifies specific Met, Trp, Tyr, and Pro side-chains, induces inter-protein Tyr63-Tyr41 cross-links, and decreases DNA binding via oxidation of Met114/Pro110 residues in close proximity at the CRP dimer interface. UV exposure also modifies DNA-bound cAMP-CRP, with this resulting in DNA cleavage at specific G/C residues within the sequence bound to CRP, but not at other G/C sites. Oxidation also increases CRP dissociation from DNA. The modifications at the CRP dimer interface, and the site-specific DNA strand cleavage are proposed to occur via oxidation of two species Met residues (Met114 and Met189, respectively) to reactive persulfoxides that damage neighbouring amino acids and DNA bases. These data suggest that modification to CRP, and bound DNA, contributes to UV sensitivity.

Intramolecular Interactions Dominate the Autoregulation of Escherichia coli Stringent Factor RelA.

Amino acid starvation in Escherichia coli activates the enzymatic activity of the stringent factor RelA, leading to accumulation of the alarmone nucleotide (p)ppGpp. The alarmone acts as an intercellular messenger to regulate transcription, translation and metabolism to mediate bacterial stress adaptation. The enzymatic activity of RelA is subject to multi-layered allosteric control executed both by ligands - such as "starved" ribosomal complexes, deacylated tRNA and pppGpp - and by individual RelA domains. The auto-regulation of RelA is proposed to act either in cis (inhibition of the enzymatic activity of the N-terminal region, NTD, by regulatory C-terminal region, CTD) or in trans (CTD-mediated dimerization leading to enzyme inhibition). In this report, we probed the regulatory roles of the individual domains of E. coli RelA and our results are not indicative of RelA dimerization being the key regulatory mechanism. First, at growth-permitting levels, ectopic expression of RelA CTD does not interfere with activation of native RelA, indicating lack of regulation via inhibitory complex formation in the cell. Second, in our biochemical assays, increasing RelA concentration does not decrease the enzyme activity, as would be expected in the case of efficient auto-inhibition via dimerization. Third, while high-level CTD expression efficiently inhibits the growth, the effect is independent of native RelA and is mediated by direct inhibition of protein synthesis, likely via direct interaction with the ribosomal A-site. Finally, deletion of the RRM domain of the CTD region leads to growth inhibition mediated by accumulation of (p)ppGpp, suggesting de-regulation of the synthetic activity in this mutant.

Quorum Sensing-Regulated Phenol-Soluble Modulins Limit Persister Cell Populations in Staphylococcus aureus.

Incomplete killing of bacterial pathogens by antibiotics is an underlying cause of treatment failure and accompanying complications. Among those avoiding chemotherapy are persisters being individual cells in a population that for extended periods of time survive high antibiotic concentrations proposedly by being in a quiescent state refractory to antibiotic killing. While investigating the human pathogen Staphylococcus aureus and the influence of growth phase on persister formation, we noted that spent supernatants of stationary phase cultures of S. aureus or S. epidermidis, but not of distantly related bacteria, significantly reduced the persister cell frequency upon ciprofloxacin challenge when added to exponentially growing and stationary phase S. aureus cells. Curiously, the persister reducing activity of S. aureus supernatants was also effective against persisters formed by either S. carnosus or Listeria monocytogenes. The persister reducing component, which resisted heat but not proteases and was produced in the late growth phase in an agr quorum-sensing dependent manner, was identified to be the phenol-soluble modulin (PSM) toxins. S. aureus express several PSMs, each with distinct cytolytic and antimicrobial properties; however, the persister reducing activity was specifically linked to synthesis of the PSMα family. Correspondingly, a high-persister phenotype of a PSMα mutant was observed upon fluoroquinolone or aminoglycoside challenge, demonstrating that the persister reducing activity of PSMs can be endogenously synthesized or extrinsically added. Given that PSMs have been associated with lytic activity against bacterial membranes we propose that PSM toxins increase the susceptibility of persister cells to killing by intracellularly acting antibiotics and that chronic and re-occurring infections with quorum sensing, agr negative mutants may be difficult to treat with antibiotics because of persister cells formed in the absence of PSM toxins.

Barley plants over-expressing the NAC transcription factor gene HvNAC005 show stunting and delay in development combined with early senescence.

The plant-specific NAC transcription factors have attracted particular attention because of their involvement in stress responses, senescence, and nutrient remobilization. The HvNAC005 gene of barley encodes a protein belonging to subgroup NAC-a6 of the NAC family. This study shows that HvNAC005 is associated with developmental senescence. It was significantly up-regulated following ABA treatment, supported by ABA-responsive elements in its promoter, but it was not up-regulated during dark-induced senescence. The C-termini of proteins closely related to HvNAC005 showed overall high divergence but also contained conserved short motifs. A serine- and leucine-containing central motif was essential for transcriptional activity of the HvNAC005 C-terminus in yeast. Over-expression of HvNAC005 in barley resulted in a strong phenotype with delayed development combined with precocious senescence. The over-expressing plants showed up-regulation of genes involved with secondary metabolism, hormone metabolism, stress, signalling, development, and transport. Up-regulation of senescence markers and hormone metabolism and signalling genes supports a role of HvNAC005 in the cross field of different hormone and signalling pathways. Binding of HvNAC005 to promoter sequences of putative target genes containing the T[G/A]CGT core motif was shown by direct protein-DNA interactions of HvNAC005 with promoters for two of the up-regulated genes. In conclusion, HvNAC005 was shown to be a strong positive regulator of senescence and so is an obvious target for the fine-tuning of gene expression in future attempts to improve nutrient remobilization related to the senescence process in barley.

A DNA-binding-site landscape and regulatory network analysis for NAC transcription factors in Arabidopsis thaliana.

Target gene identification for transcription factors is a prerequisite for the systems wide understanding of organismal behaviour. NAM-ATAF1/2-CUC2 (NAC) transcription factors are amongst the largest transcription factor families in plants, yet limited data exist from unbiased approaches to resolve the DNA-binding preferences of individual members. Here, we present a TF-target gene identification workflow based on the integration of novel protein binding microarray data with gene expression and multi-species promoter sequence conservation to identify the DNA-binding specificities and the gene regulatory networks of 12 NAC transcription factors. Our data offer specific single-base resolution fingerprints for most TFs studied and indicate that NAC DNA-binding specificities might be predicted from their DNA-binding domain's sequence. The developed methodology, including the application of complementary functional genomics filters, makes it possible to translate, for each TF, protein binding microarray data into a set of high-quality target genes. With this approach, we confirm NAC target genes reported from independent in vivo analyses. We emphasize that candidate target gene sets together with the workflow associated with functional modules offer a strong resource to unravel the regulatory potential of NAC genes and that this workflow could be used to study other families of transcription factors.

A novel indirect sequence readout component in the E. coli cyclic AMP receptor protein operator.

The cyclic AMP receptor protein (CRP) from Escherichia coli has been extensively studied for several decades. In particular, a detailed characterization of CRP interaction with DNA has been obtained. The CRP dimer recognizes a consensus sequence AANTGTGANNNNNNTCACANTT through direct amino acid nucleobase interactions in the major groove of the two operator half-sites. Crystal structure analyses have revealed that the interaction results in two strong kinks at the TG/CA steps closest to the 6-base-pair spacer (N6). This spacer exhibits high sequence variability among the more than 100 natural binding sites in the E. coli genome, but the exact role of the N6 region in CRP interaction has not previously been systematic examined. Here we employ an in vitro selection system based on a randomized N6 spacer region to demonstrate that CRP binding to the lacP1 site may be enhanced up to 14-fold or abolished by varying the N6 spacer sequences. Furthermore, on the basis of sequence analysis and uranyl (UO2(2+)) probing data, we propose that the underlying mechanism relies on N6 deformability.

ATAF1 transcription factor directly regulates abscisic acid biosynthetic gene NCED3 in Arabidopsis thaliana.

ATAF1, an Arabidopsis thaliana NAC transcription factor, plays important roles in plant adaptation to environmental stress and development. To search for ATAF1 target genes, we used protein binding microarrays and chromatin-immunoprecipitation (ChIP). This identified T[A,C,G]CGT[A,G] and TT[A,C,G]CGT as ATAF1 consensus binding sequences. Co-expression analysis across publicly available microarray experiments identified 25 genes co-expressed with ATAF1. The promoter regions of ATAF1 co-expressors were significantly enriched for ATAF1 binding sites, and TTGCGTA was identified in the promoter of the key abscisic acid (ABA) phytohormone biosynthetic gene NCED3. ChIP-qPCR and expression analysis showed that ATAF1 binding to the NCED3 promoter correlated with increased NCED3 expression and ABA hormone levels. These results indicate that ATAF1 regulates ABA biosynthesis.

Structure, function and networks of transcription factors involved in abiotic stress responses.

Transcription factors (TFs) are master regulators of abiotic stress responses in plants. This review focuses on TFs from seven major TF families, known to play functional roles in response to abiotic stresses, including drought, high salinity, high osmolarity, temperature extremes and the phytohormone ABA. Although ectopic expression of several TFs has improved abiotic stress tolerance in plants, fine-tuning of TF expression and protein levels remains a challenge to avoid crop yield loss. To further our understanding of TFs in abiotic stress responses, emerging gene regulatory networks based on TFs and their direct targets genes are presented. These revealed components shared between ABA-dependent and independent signaling as well as abiotic and biotic stress signaling. Protein structure analysis suggested that TFs hubs of large interactomes have extended regions with protein intrinsic disorder (ID), referring to their lack of fixed tertiary structures. ID is now an emerging topic in plant science. Furthermore, the importance of the ubiquitin-proteasome protein degradation systems and modification by sumoylation is also apparent from the interactomes. Therefore; TF interaction partners such as E3 ubiquitin ligases and TF regions with ID represent future targets for engineering improved abiotic stress tolerance in crops.

DNA binding by the plant-specific NAC transcription factors in crystal and solution: a firm link to WRKY and GCM transcription factors.

NAC (NAM/ATAF/CUC) plant transcription factors regulate essential processes in development, stress responses and nutrient distribution in important crop and model plants (rice, Populus, Arabidopsis), which makes them highly relevant in the context of crop optimization and bioenergy production. The structure of the DNA-binding NAC domain of ANAC019 has previously been determined by X-ray crystallography, revealing a dimeric and predominantly ß-fold structure, but the mode of binding to cognate DNA has remained elusive. In the present study, information from low resolution X-ray structures and small angle X-ray scattering on complexes with oligonucleotides, mutagenesis and (DNase I and uranyl photo-) footprinting, is combined to form a structural view of DNA-binding, and for the first time provide experimental evidence for the speculated relationship between plant-specific NAC proteins, WRKY transcription factors and the mammalian GCM (Glial cell missing) transcription factors, which all use a ß-strand motif for DNA-binding. The structure shows that the NAC domain inserts the edge of its core ß-sheet into the major groove, while leaving the DNA largely undistorted. The structure of the NAC-DNA complex and a new crystal form of the unbound NAC also indicate limited flexibility of the NAC dimer arrangement, which could be important in recognizing suboptimal binding sites.

A DNA minor groove electronegative potential genome map based on photo-chemical probing.

The double-stranded DNA of the genome contains both sequence information directly relating to the protein and RNA coding as well as functional and structural information relating to protein recognition. Only recently is the importance of DNA shape in this recognition process being fully appreciated, and it also appears that minor groove electronegative potential may contribute significantly in guiding proteins to their cognate binding sites in the genome. Based on the photo-chemical probing results, we have derived an algorithm that predicts the minor groove electronegative potential in a DNA helix of any given sequence. We have validated this model on a series of protein-DNA binding sites known to involve minor groove electrostatic recognition as well as on stable nucleosome core complexes. The algorithm allows for the first time a full minor groove electrostatic description at the nucleotide resolution of any genome, and it is illustrated how such detailed studies of this sequence dependent, inherent property of the DNA may reflect on genome organization, gene expression and chromosomal condensation.

One-step isolation of plasma membrane proteins using magnetic beads with immobilized concanavalin A.

We have developed a simple method for isolating and purifying plasma membrane proteins from various cell types. This one-step affinity-chromatography method uses the property of the lectin concanavalin A (ConA) and the technique of magnetic bead separation to obtain highly purified plasma membrane proteins from crude membrane preparations or cell lines. ConA is immobilized onto magnetic beads by binding biotinylated ConA to streptavidin magnetic beads. When these ConA magnetic beads were used to enrich plasma membranes from a crude membrane preparation, this procedure resulted in 3.7-fold enrichment of plasma membrane marker 5'-nucleotidase activity with 70% recovery of the activity in the crude membrane fraction of rat liver. In agreement with the results of 5'-nucleotidase activity, immunoblotting with antibodies specific for a rat liver plasma membrane protein, CEACAM1, indicated that CEACAM1 was enriched about threefold relative to that of the original membranes. In similar experiments, this method produced 13-fold enrichment of 5'-nucleotidase activity with 45% recovery of the activity from a total cell lysate of PC-3 cells and 7.1-fold enrichment of 5'-nucleotidase activity with 33% recovery of the activity from a total cell lysate of HeLa cells. These results suggest that this one-step purification method can be used to isolate total plasma membrane proteins from tissue or cells for the identification of membrane biomarkers.

Dissecting direct and indirect readout of cAMP receptor protein DNA binding using an inosine and 2,6-diaminopurine in vitro selection system.

The DNA interaction of the Escherichia coli cyclic AMP receptor protein (CRP) represents a typical example of a dual recognition mechanism exhibiting both direct and indirect readout. We have dissected the direct and indirect components of DNA recognition by CRP employing in vitro selection of a random library of DNA-binding sites containing inosine (I) and 2,6-diaminopurine (D) instead of guanine and adenine, respectively. Accordingly, the DNA helix minor groove is structurally altered due to the 'transfer' of the 2-amino group of guanine (now I) to adenine (now D), whereas the major groove is functionally intact. The majority of the selected sites contain the natural consensus sequence TGTGAN(6)TCACA (i.e. TITIDN(6)TCDCD). Thus, direct readout of the consensus sequence is independent of minor groove conformation. Consequently, the indirect readout known to occur in the TG/CA base pair step (primary kink site) in the consensus sequence is not affected by I-D substitutions. In contrast, the flanking regions are selected as I/C rich sequences (mostly I-tracts) instead of A/T rich sequences which are known to strongly increase CRP binding, thereby demonstrating almost exclusive indirect readout of helix structure/flexibility in this region through (anisotropic) flexibility of I-tracts.

Chemokinesis by tetrahymena in response to bacterial oligopeptides.

The ciliate Tetrahymena responds very efficiently by chemoattraction to a group of trichloroacetic acid-soluble oligopeptides isolated from a commercial bioprotein from Methanococcus. When fractionated by reversed phase C18-high-pressure liquid chromatography, this group of very efficient chemoattractants turned out to consist of a heterogeneous group of oligopeptides with molecular weight ranging from 0.2 to 1.5 kDa. The peptides were very rich in the following amino acids: aspartic acid, alanine, glutamic acid, proline, glycine, lysine, and arginine. The term chemokinesis is used throughout to emphasise that chemoattraction does not necessarily include an element of orientation of cells.

Uranyl photoprobing of nonbent A/T- and bent A-tracts. A difference of flexibility?

In this study, we have systematically compared the uranyl photocleavage of a range of bent A-tracts and nonbent TA-tracts as well as interrupted A-tracts. We demonstrate that uranyl photocleavage of A-tracts and TA-tracts is almost identical, indicating a very similar minor groove conformation. Furthermore, a 10 base pair A-tract is divided into two independent tracts by an intervening TA or GC step. Uranyl probing also clearly distinguishes the bent A4T4 and the nonbent T4A4 sequences as adopting different structures, and our interpretation of the data is consistent with a structure for the bent A4T4 sequence that resembles a continuous A-tract, whereas the nonbent T4A4 sequences are closer to two independent and opposite A-tracts that cancel each other in terms of macroscopic bending. Finally, we also note that even single TA and TAT steps are highly sensitive to uranyl photocleavage and propose that in addition to average minor groove width, uranyl also senses DNA helix flexibility/deformability. Thus, the structural difference of TA-tracts and A-tracts may to a large extent reflect a difference in flexibility, and DNA curvature may consequently require a rigid narrow minor groove conformation that creates distinct A-tract-B-DNA junctions as the predominant cause of the bending.

Polyamines preferentially interact with bent adenine tracts in double-stranded DNA.

Polyamines, such as putrescine, spermidine and spermine, have indirectly been linked with the regulation of gene expression, and their concentrations are typically increased in cancer cells. Although effects on transcription factor binding to cognate DNA targets have been demonstrated, the mechanisms of the biological action of polyamines is poorly understood. Employing uranyl photo-probing we now demonstrate that polyamines at submillimolar concentrations bind preferentially to bent adenine tracts in double-stranded DNA. These results provide the first clear evidence for the sequence-specific binding of polyamines to DNA, and thereby suggest a mechanism by which the cellular effects of polyamines in terms of differential gene transcriptional activity could, at least partly, be a direct consequence of sequence-specific interactions of polyamines with promoters at the DNA sequence level.