The molecular basis for beta-lactamase gene expression in B. anthracis, B. cereus and B. thuringiensis

The susceptibility of most Bacillus anthracis strains to β-lactam antibiotics is intriguing considering that the B. anthracis genome harbors two β-lactamase genes, bla1 and bla2, and closely-related species, Bacillus cereus and Bacillus thuringiensis, typically produce β-lactamases. This work demonstrates that B. anthracis bla expression is affected by two genes, sigP and rsp, predicted to encode an extracytoplasmic function sigma factor and an antisigma factor, respectively. Deletion of the sigP/rsp locus abolished bla expression in a penicillin-resistant clinical isolate and had no effect on bla expression in a prototypical penicillin-susceptible strain. Complementation with sigP/rsp from the penicillin-resistant strain, but not the penicillin-susceptible strain, conferred β-lactamase activity upon both mutants. These results are attributed to a nucleotide deletion near the 5' end of rsp in the penicillin-resistant strain that is predicted to result in a nonfunctional protein. B. cereus and B. thuringiensis sigP and rsp homologues are required for inducible penicillin resistance in those species. Expression of the B. cereus or B. thuringiensis sigP and rsp genes in a B. anthracis sigP/rsp-null mutant confers resistance to β-lactam antibiotics, suggesting that while B. anthracis contains the genes necessary for sensing β-lactam antibiotics, the B. anthracis sigP/rsp gene products are insufficient for bla induction. ^ Because alternative sigma factors recognize unique promoter sequence, direct targets can be elucidated by comparing transcriptional profiling results with an in silico search using the sigma factor binding sequence. Potential σP -10 and -35 promoter elements were identified upstream from bla1 bla2 and sigP. Results obtained from searching the B. anthracis genome with the conserved sequences were evaluated against transcriptional profiling results comparing B. anthracis 32 and an isogenic sigP/rsp -null strain. Results from these analyses indicate that while the absence of the sigP gene significantly affects the transcript levels of 16 genes, only bla1, bla2 and sigP are directly regulated by σP. The genomes of B. cereus and B. thuringiensis strains were also analyzed for the potential σP binding elements. The sequence was located upstream from the sigP and bla genes, and previously unidentified genes predicted to encode a penicillin-binding protein (PBP) and a D-alanyl-D-alanine carboxypeptidase, indicating that the σ P regulon in these species responds to cell-wall stress caused by β-lactam antibiotics. ^ β-lactam antibiotics prevent attachment of new peptidoglycan to the cell wall by blocking the active site of PBPs. A B. cereus and B. thuringiensis pbp-encoding gene located near bla1 contains a potential σP recognition sequence upstream from the annotated translational start. Deletion of this gene abolished β-lactam resistance in both strains. Mutations in the active site of the PBP were detrimental to β-lactam resistance in B. cereus, but not B. thuringiensis, indicating that the transpeptidase activity is only important in B. cereus. I also found that transcript levels of the PBP-encoding gene are not significantly affected by the presence of β-lactam antibiotic. Based on these data I hypothesize that the gene product acts a sensor of β-lactam antibiotic. ^

Inhibition of gene expression in human cells through small molecule-RNA interactions

Small molecules that bind their biological receptors with high affinity and selectivity can be isolated from randomized pools of combinatorial libraries. RNA-protein interactions are important in many cellular functions, including transcription, RNA splicing, and translation. One example of such interactions is the mechanism of trans-activation of HIV-1 gene expression that requires the interaction of Tat protein with the trans-activation responsive region (TAR) RNA, a 59-base stem-loop structure located at the 5′ end of all nascent HIV-1 transcripts. Here we demonstrate the isolation of small TAR RNA-binding molecules from an encoded combinatorial library. We have made an encoded combinatorial tripeptide library of 24,389 possible members from d-and l-alpha amino acids on TentaGel resin. Using on-bead screening we have identified a small family of mostly heterochiral tripeptides capable of structure-specific binding to the bulge loop of TAR RNA. In vitro binding studies reveal stereospecific discrimination when the best tripeptide ligand is compared with diastereomeric peptide sequences. In addition, the most strongly binding tripeptide was shown to suppress transcriptional activation by Tat protein in human cells with an IC50 of ≈50 nM. Our results indicate that tripeptide RNA ligands are cell permeable, nontoxic to cells, and capable of inhibiting expression of specific genes by interfering with RNA-protein interactions.

An artificial cell-cycle inhibitor isolated from a combinatorial library

Understanding the genetic networks that operate inside cells will require the dissection of interactions among network members. Here we describe a peptide aptamer isolated from a combinatorial library that distinguishes among such interactions. This aptamer binds to cyclin-dependent kinase 2 (Cdk2) and inhibits its kinase activity. In contrast to naturally occurring inhibitors, such as p21Cip1, which inhibit the activity of Cdk2 on all its substrates, inhibition by pep8 has distinct substrate specificity. We show that the aptamer binds to Cdk2 at or near its active site and that its mode of inhibition is competitive. Expression of pep8 in human cells retards their progression through the G1 phase of the cell cycle. Our results suggest that the aptamer inhibits cell-cycle progression by blocking the activity of Cdk2 on substrates needed for the G1-to-S transition. This work demonstrates the feasibility of selection of artificial proteins to perform functions not developed during evolution. The ability to select proteins that block interactions between a gene product and some partners but not others should make sophisticated genetic manipulations possible in human cells and other currently intractable systems.

The DNMT3B DNA methyltransferase gene is mutated in the ICF immunodeficiency syndrome

DNA methylation is an important regulator of genetic information in species ranging from bacteria to humans. DNA methylation appears to be critical for mammalian development because mice nullizygous for a targeted disruption of the DNMT1 DNA methyltransferase die at an early embryonic stage. No DNA methyltransferase mutations have been reported in humans until now. We describe here the first example of naturally occurring mutations in a mammalian DNA methyltransferase gene. These mutations occur in patients with a rare autosomal recessive disorder, which is termed the ICF syndrome, for immunodeficiency, centromeric instability, and facial anomalies. Centromeric instability of chromosomes 1, 9, and 16 is associated with abnormal hypomethylation of CpG sites in their pericentromeric satellite regions. We are able to complement this hypomethylation defect by somatic cell fusion to Chinese hamster ovary cells, suggesting that the ICF gene is conserved in the hamster and promotes de novo methylation. ICF has been localized to a 9-centimorgan region of chromosome 20 by homozygosity mapping. By searching for homologies to known DNA methyltransferases, we identified a genomic sequence in the ICF region that contains the homologue of the mouse Dnmt3b methyltransferase gene. Using the human sequence to screen ICF kindreds, we discovered mutations in four patients from three families. Mutations include two missense substitutions and a 3-aa insertion resulting from the creation of a novel 3′ splice acceptor. None of the mutations were found in over 200 normal chromosomes. We conclude that mutations in the DNMT3B are responsible for the ICF syndrome.

Making artificial antibodies: A format for phage display of combinatorial heterodimeric arrays

The gene VII protein (pVII) and gene IX protein (pIX) are associated closely on the surface of filamentous bacteriophage that is opposite of the end harboring the widely exploited pIII protein. We developed a phagemid format wherein antibody heavy- and light-chain variable regions were fused to the amino termini of pVII and pIX, respectively. Significantly, the fusion proteins interacted to form a functional Fv-binding domain on the phage surface. Our approach will be applicable to the display of generic peptide and protein libraries that can form combinatorial heterodimeric arrays. Consequently, it represents a first step toward artificial antibodies and the selection of novel biological activities.

Combinatorial roles for pRB, p107, and p130 in E2F-mediated cell cycle control

Numerous studies have implicated the pRB family of nuclear proteins in the control of cell cycle progression. Although over-expression experiments have revealed that each of these proteins, pRB, p107, and p130, can induce a G1 cell cycle arrest, mouse knockouts demonstrated distinct developmental requirements for these proteins, as well as partial functional redundancy between family members. To study the mechanism by which the closely related pRB family proteins contribute to cell cycle progression, we generated 3T3 fibroblasts derived from embryos that lack one or more of these proteins (pRB−/−, p107−/−, p130−/−, pRB−/−/p107−/−, pRB−/−/p130−/−, and p107−/−/p130−/−). By comparing the growth and cell cycle characteristics of these cells, we have observed clear differences in the manner in which they transit through the G1 and S phases as well as exit from the cell cycle. Deletion of Rb, or more than one of the family members, results in a shortening of G1 and a lengthening of S phase, as well as a reduction in growth factor requirements. In addition, the individual cell lines showed differential regulation of a subset of E2F-dependent gene promoters, as well as differences in cell cycle-dependent kinase activity. Taken together, these observations suggest that the closely related pRB family proteins affect cell cycle progression through distinct biochemical mechanisms and that their coordinated action may contribute to their diverse functions in various physiological settings.

Selection of peptides that functionally replace a zinc finger in the Sp1 transcription factor by using a yeast combinatorial library

We have developed a strategy for the identification of peptides able to functionally replace a zinc finger domain in a transcription factor. This strategy could have important ramifications for basic research on gene regulation and for the development of therapeutic agents. In this study in yeast, we expressed chimeric proteins that included a random peptide combinatorial library in association with two zinc finger domains and a transactivating domain. The library was screened for chimeric proteins capable of activating transcription from a target sequence in the upstream regulatory regions of selectable or reporter genes. In a screen of approximately 1.5 × 107 transformants we identified 30 chimeric proteins that exhibited transcriptional activation, some of which were able to discriminate between wild-type and mutant DNA targets. Chimeric library proteins expressed as glutathione S-transferase fusions bound to double-stranded oligonucleotides containing the target sequence, suggesting that the chimeras bind directly to DNA. Surprisingly, none of the peptides identified resembled a zinc finger or other well-known transcription factor DNA binding domain.

Analysis of endoplasmic reticulum trafficking signals by combinatorial screening in mammalian cells

To improve the accuracy of predicting membrane protein sorting signals, we developed a general methodology for defining trafficking signal consensus sequences in the environment of the living cell. Our approach uses retroviral gene transfer to create combinatorial expression libraries of trafficking signal variants in mammalian cells, flow cytometry to sort cells based on trafficking phenotype, and quantitative trafficking assays to measure the efficacy of individual signals. Using this strategy to analyze arginine- and lysine-based endoplasmic reticulum localization signals, we demonstrate that small changes in the local sequence context dramatically alter signal strength, generating a broad spectrum of trafficking phenotypes. Finally, using sequences from our screen, we found that the potency of di-lysine, but not di-arginine, mediated endoplasmic reticulum localization was correlated with the strength of interaction with α-COP.

Identification and dissection of an enhancer controlling epithelial gene expression in skin

Keratins 14 and 5 are the structural hallmarks of the basal keratinocytes of the epidermis and outer root sheath (ORS) of the hair follicle. Their genes are controlled in a tissue-specific manner and thus serve as useful tools to elucidate the regulatory mechanisms involved in keratinocyte-specific transcription. Previously we identified several keratinocyte-specific DNase I hypersensitive sites (HSs) in the 5′ regulatory sequences of the K14 gene and showed that a 700-bp regulatory domain encompassing HSs II and III can confer epidermal and ORS-specific gene expression in transgenic mice in vivo. Although HS II harbored much of the transactivation activity in vitro, it was not sufficient to restrict expression to keratinocytes in vivo. We now explore the HS III regulatory element. Surprisingly, this element on its own confers gene expression to the keratinocytes of the inner root sheath (IRS) of the hair follicle, whereas a 275-bp DNA fragment containing both HSs II and III shifts the expression from the IRS to the basal keratinocytes and ORS in vivo. Electrophoretic mobility-shift assays and mutational studies of HSs III reveal a role for CACCC-box binding proteins, Sp1 family members, and other factors adding to the list of previously described factors that are involved in keratinocyte-specific gene expression. These studies highlight a cooperative interaction of the two HSs domains and strengthen the importance of combinatorial play of transcription factors that govern keratinocyte-specific gene regulation.

Pseudoknots in prion protein mRNAs confirmed by comparative sequence analysis and pattern searching

The human prion gene contains five copies of a 24 nt repeat that is highly conserved among species. An analysis of folding free energies of the human prion mRNA, in particular in the repeat region, suggested biased codon selection and the presence of RNA patterns. In particular, pseudoknots, similar to the one predicted by Wills in the human prion mRNA, were identified in the repeat region of all available prion mRNAs available in GenBank, but not those of birds and the red slider turtle. An alignment of these mRNAs, which share low sequence homology, shows several co-variations that maintain the pseudoknot pattern. The presence of pseudoknots in yeast Sup35p and Rnq1 suggests acquisition in the prokaryotic era. Computer generated three-dimensional structures of the human prion pseudoknot highlight protein and RNA interaction domains, which suggest a possible effect in prion protein translation. The role of pseudoknots in prion diseases is discussed as individuals with extra copies of the 24 nt repeat develop the familial form of Creutzfeldt–Jakob disease.

Tbx19, a tissue-selective regulator of POMC gene expression

Pituitary cell types arise in a temporally and spatially specific fashion, in response to combinatorial actions of transcription factors induced by transient signaling gradients. The critical transcriptional determinants of the two pituitary cell types that express the pro-opiomelanocortin (POMC) gene, the anterior lobe corticotropes, producing adrenocorticotropin, and the intermediate lobe melanotropes, producing melanocyte-stimulating hormone (MSHα), have remained unknown. Here, we report that a member of the T-box gene family, Tbx19, which is expressed only in the rostral ventral diencephalon and pituitary gland, commencing on e11.5, marks pituitary cells that will subsequently express the POMC gene and is capable of altering progression of ventral cell types and inducing adrenocorticotropin in rostral tip cells. It is suggested that Tbx19, depending on the presence of synergizing transcription factors, can activate POMC gene expression and repress the α glycoprotein subunit and thyroid-stimulating hormone β promoters.

Combinatorial control of muscle development by basic helix-loop-helix and MADS-box transcription factors.

Members of the MyoD family of muscle-specific basic helix-loop-helix (bHLH) proteins function within a genetic pathway to control skeletal muscle development. Mutational analyses of these factors suggested that their DNA binding domains mediated interaction with a coregulator required for activation of muscle-specific transcription. Members of the myocyte enhancer binding factor 2 (MEF2) family of MADS-box proteins are expressed at high levels in muscle and neural cells and at lower levels in several other cell types. MEF2 factors are unable to activate muscle gene expression alone, but they potentiate the transcriptional activity of myogenic bHLH proteins. This potentiation appears to be mediated by direct interactions between the DNA binding domains of these different types of transcription factors. Biochemical and genetic evidence suggests that MEF2 factors are the coregulators for myogenic bHLH proteins. The presence of MEF2 and cell-specific bHLH proteins in other cell types raises the possibility that these proteins may also cooperate to regulate other programs of cell-specific gene expression. We present a model to account for such cooperative interactions.

Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.

Controlling Hox gene expression and activity to build the vertebrate axial skeleton

It has long been known that Hox genes are central players in patterning the vertebrate axial skeleton. Extensive genetic studies in the mouse have revealed that the combinatorial activity of Hox genes along the anterior-posterior body axis specifies different vertebral identities. In addition, Hox genes were instrumental for the evolutionary diversification of the vertebrate body plan. In this review, we focus on fundamental questions regarding the intricate mechanisms controlling Hox gene activity. In particular, we discuss the functional relevance of the precise timing of Hox gene activation in the embryo. Moreover, we provide insight into the epigenetic regulatory mechanisms that are likely to control this process and are responsible for the maintenance of spatially restricted Hox expression domains throughout embryonic development. We also analyze how specific features of each Hox protein may contribute to the functional diversity of Hox family. Altogether, the work reviewed here further supports the notion that the Hox program is far more complex than initially assumed. Exciting new findings will surely emerge in the years ahead.

Searching for convergence in phylogenetic Markov chain Monte Carlo

Markov chain Monte Carlo (MCMC) is a methodology that is gaining widespread use in the phylogenetics community and is central to phylogenetic software packages such as MrBayes. An important issue for users of MCMC methods is how to select appropriate values for adjustable parameters such as the length of the Markov chain or chains, the sampling density, the proposal mechanism, and, if Metropolis-coupled MCMC is being used, the number of heated chains and their temperatures. Although some parameter settings have been examined in detail in the literature, others are frequently chosen with more regard to computational time or personal experience with other data sets. Such choices may lead to inadequate sampling of tree space or an inefficient use of computational resources. We performed a detailed study of convergence and mixing for 70 randomly selected, putatively orthologous protein sets with different sizes and taxonomic compositions. Replicated runs from multiple random starting points permit a more rigorous assessment of convergence, and we developed two novel statistics, delta and epsilon, for this purpose. Although likelihood values invariably stabilized quickly, adequate sampling of the posterior distribution of tree topologies took considerably longer. Our results suggest that multimodality is common for data sets with 30 or more taxa and that this results in slow convergence and mixing. However, we also found that the pragmatic approach of combining data from several short, replicated runs into a metachain to estimate bipartition posterior probabilities provided good approximations, and that such estimates were no worse in approximating a reference posterior distribution than those obtained using a single long run of the same length as the metachain. Precision appears to be best when heated Markov chains have low temperatures, whereas chains with high temperatures appear to sample trees with high posterior probabilities only rarely. [Bayesian phylogenetic inference; heating parameter; Markov chain Monte Carlo; replicated chains.]