The RESID Database of protein structure modifications and the NRL-3D Sequence–Structure Database

The RESID Database is a comprehensive collection of annotations and structures for protein post-translational modifications including N-terminal, C-terminal and peptide chain cross-link modifications. The RESID Database includes systematic and frequently observed alternate names, Chemical Abstracts Service registry numbers, atomic formulas and weights, enzyme activities, taxonomic range, keywords, literature citations with database cross-references, structural diagrams and molecular models. The NRL-3D Sequence–Structure Database is derived from the three-dimensional structure of proteins deposited with the Research Collaboratory for Structural Bioinformatics Protein Data Bank. The NRL-3D Database includes standardized and frequently observed alternate names, sources, keywords, literature citations, experimental conditions and searchable sequences from model coordinates. These databases are freely accessible through the National Cancer Institute–Frederick Advanced Biomedical Computing Center at these web sites: http://www.ncifcrf.gov/RESID, http://www.ncifcrf.gov/ NRL-3D; or at these National Biomedical Research Foundation Protein Information Resource web sites: http://pir.georgetown.edu/pirwww/dbinfo/resid.html, http://pir.georgetown.edu/pirwww/dbinfo/nrl3d.html

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs.bionet.nsc.ru/systems/Activity/.

Mapping the active site in vasoactive intestinal peptide to a core of four amino acids: Neuroprotective drug design

The understanding of the molecular mechanisms leading to peptide action entails the identification of a core active site. The major 28-aa neuropeptide, vasoactive intestinal peptide (VIP), provides neuroprotection. A lipophilic derivative with a stearyl moiety at the N-terminal and norleucine residue replacing the Met-17 was 100-fold more potent than VIP in promoting neuronal survival, acting at femtomolar–picomolar concentration. To identify the active site in VIP, over 50 related fragments containing an N-terminal stearic acid attachment and an amidated C terminus were designed, synthesized, and tested for neuroprotective properties. Stearyl-Lys-Lys-Tyr-Leu-NH2 (derived from the C terminus of VIP and the related peptide, pituitary adenylate cyclase activating peptide) captured the neurotrophic effects offered by the entire 28-aa parent lipophilic derivative and protected against β-amyloid toxicity in vitro. Furthermore, the 4-aa lipophilic peptide recognized VIP-binding sites and enhanced choline acetyltransferase activity as well as cognitive functions in Alzheimer’s disease-related in vivo models. Biodistribution studies following intranasal administration of radiolabeled peptide demonstrated intact peptide in the brain 30 min after administration. Thus, lipophilic peptide fragments offer bioavailability and stability, providing lead compounds for drug design against neurodegenerative diseases.

Synthetic genes for glycoprotein design and the elucidation of hydroxyproline-O-glycosylation codes

Design of hydroxyproline (Hyp)-rich glycoproteins (HRGPs) offers an approach for the structural and functional analysis of these wall components, which are broadly implicated in plant growth and development. HRGPs consist of multiple small repetitive “glycomodules” extensively O-glycosylated through the Hyp residues. The patterns of Hyp-O-glycosylation are putatively coded by the primary sequence as described by the Hyp contiguity hypothesis, which predicts contiguous Hyp residues to be attachment sites of small arabinooligosaccharides (1–5 Ara residues/Hyp); while clustered, noncontiguous Hyp residues are sites of arabinogalactan polysaccharide attachment. As a test, we designed two simple HRGPs as fusion proteins with green fluorescent protein. The first was a repetitive Ser-Hyp motif that encoded only clustered noncontiguous Hyp residues, predicted polysaccharide addition sites. The resulting glycoprotein had arabinogalactan polysaccharide O-linked to all Hyp residues. The second construct, based on the consensus sequence of a gum arabic HRGP, contained both arabinogalactan and arabinooligosaccharide addition sites and, as predicted, gave a product that contained both saccharide types. These results identify an O-glycosylation code of plants.

Topoisomerase II-mediated site-directed alkylation of DNA by psorospermin and its use in mapping other topoisomerase II poison binding sites

Psorospermin is a plant natural product that shows significant in vivo activity against P388 mouse leukemia. The molecular basis for this selectivity is unknown, although psorospermin has been demonstrated to intercalate into DNA and alkylate N7 of guanine. Significantly, the alkylation reactivity of psorospermin at specific sites on DNA increased 25-fold in the presence of topoisomerase II. In addition, psorospermin trapped the topoisomerase II-cleaved complex formation at the same site. These results imply that the efficacy of psorospermin is related to its interaction with the topoisomerase II–DNA complex. Because thermal treatment of (N7 guanine)–DNA adducts leads to DNA strand breakage, we were able to determine the site of alkylation of psorospermin within the topoisomerase II gate site and infer that intercalation takes place at the gate site between base pairs at the +1 and +2 positions. These results provide not only additional mechanistic information on the mode of action of the anticancer agent psorospermin but also structural insights into the design of an additional class of topoisomerase II poisons. Because the alkylation site for psorospermin in the presence of topoisomerase II can be assigned unambiguously and the intercalation site inferred, this drug is a useful probe for other topoisomerase poisons where the sites for interaction are less well defined.

Synthesis and biological activity of DNA damaging agents that form decoy binding sites for the estrogen receptor

It is a goal of cancer chemotherapy to achieve the selective killing of tumor cells while minimizing toxicity to normal tissues. We describe the design of selective toxins forming DNA adducts that attract the estrogen receptor (ER), a transcription factor that is overexpressed in many human breast and ovarian tumors. The compounds consist of 4-(3-aminopropyl)-N,N-(2-chloroethyl)-aniline linked to 2-(4′-hydroxyphenyl)-3-methyl-5-hydroxy-indole. The former moiety is a DNA damaging nitrogen mustard and the latter is a ligand for the ER. The connection between these groups was refined to permit DNA adducts formed by the mustard portion of the molecule to present the ligand domain so that it was able to interact efficiently with the ER. By using 16-mers containing specific DNA adducts, it was determined that monoadducts and putative intrastrand crosslinks were preferred targets for the ER over interstrand crosslinks. A series of structurally related 2-phenylindole mustards was prepared, some of which were selectively toxic to the ER-positive breast cancer cell line MCF-7, as compared with the ER(−) negative line MDA-MB231. The ability both to bind to DNA and to interact significantly with the ER were essential to achieve selective lethality toward ER(+) cells. Compounds forming DNA adducts without the ability to bind receptor showed similar toxicities in the two cell lines. Several models could explain the selective toxicity of the mustard–phenylindole compounds toward ER(+) cells. The favored model suggests that a mustard–DNA adduct is shielded by the ER from DNA repair enzymes and hence cells possessing an abundance of the ER selectively retain the adduct and are killed.

Design of polyzinc finger peptides with structured linkers

Zinc finger domains are perhaps the most versatile of all known DNA binding domains. By fusing up to six zinc finger modules, which normally recognize up to 18 bp of DNA, designer transcription factors can be produced to target unique sequences within large genomes. However, not all continuous DNA sequences make good zinc finger binding sites. To avoid having to target unfavorable DNA sequences, we designed multizinc finger peptides with linkers capable of spanning long stretches of nonbound DNA. Two three-finger domains were fused by using either transcription factor IIIA for the Xenopus 5S RNA gene (TFIIIA) finger 4 or a non-sequence-specific zinc finger as a “structured” linker. Our gel-shift results demonstrate that these peptides are able to bind with picomolar affinities to target sequences containing 0–10 bp of nonbound DNA. Furthermore, these peptides display greater sequence selectivity and bind with higher affinity than similar six-finger peptides containing long, flexible linkers. These peptides are likely to be of use in understanding the behavior of polydactyl proteins in nature and in the targeting of human, animal, or plant genomes for numerous applications. We also suggest that in certain polydactyl peptides an individual finger can “flip” out of the major groove to allow its neighbors to bind shorter, nontarget DNA sequences.

Identification of two short internal ribosome entry sites selected from libraries of random oligonucleotides

Sequences that control translation of mRNA may play critical roles in regulating protein levels. One such element is the internal ribosome entry site (IRES). We previously showed that a 9-nt segment in the 5′ leader sequence of the mRNA encoding Gtx homeodomain protein could function as an IRES. To identify other short sequences with similar properties, we designed a selection procedure that uses a retroviral vector to express dicistronic mRNAs encoding enhanced green and cyan fluorescent proteins as the first and second cistrons, respectively. Expression of the second cistron was dependent upon the intercistronic sequences and was indicative of IRES activity. B104 cells were infected with two retroviral libraries that contained random sequences of 9 or 18 nt in the intercistronic region. Cells expressing both cistrons were sorted, and sequences recovered from selected cells were reassayed for IRES activity in a dual luciferase dicistronic mRNA. Two novel IRESes were identified by this procedure, and both contained segments with complementarity to 18S rRNA. When multiple copies of either segment were linked together, IRES activities were dramatically enhanced. Moreover, these synthetic IRESes were differentially active in various cell types. These properties are similar to those of the previously identified 9-nt IRES module from Gtx mRNA. These results provide further evidence that short nucleotide sequences can function as IRESes and support the idea that some cellular IRESes may be composed of shorter functional modules. The ability to identify IRES modules with specific expression properties may be useful in the design of vectors for biotechnology and gene therapy.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

RegulonDB (version 3.2): transcriptional regulation and operon organization in Escherichia coli K-12

RegulonDB is a database on mechanisms of transcription regulation and operon organization in Escherichia coli K-12. The current version has considerably increased numbers of regulatory elements such as promoters, binding sites and terminators. The complete repertoire of known and predicted DNA-binding transcriptional regulators can be considered to be included in this version. The database now distinguishes different allosteric conformations of regulatory proteins indicating the one active in binding and regulating the different promoters. A new set of operon predictions has been incorporated. The relational design has been modified accordingly. Furthermore, a major improvement is a graphic display enabling browsing of the database with a Java-based graphic user interface with three zoom-levels connected to properties of each chromo­somal element. The purpose of these modifications is to make RegulonDB a useful tool and control set for tran­scriptome experiments. RegulonDB can be accessed on the web at the URL: http://www.cifn.unam.mx/Computational_Biology/regulondb/

SpliceDB: database of canonical and non-canonical mammalian splice sites

A database (SpliceDB) of known mammalian splice site sequences has been developed. We extracted 43 337 splice pairs from mammalian divisions of the gene-centered Infogene database, including sites from incomplete or alternatively spliced genes. Known EST sequences supported 22 815 of them. After discarding sequences with putative errors and ambiguous location of splice junctions the verified dataset includes 22 489 entries. Of these, 98.71% contain canonical GT–AG junctions (22 199 entries) and 0.56% have non-canonical GC–AG splice site pairs. The remainder (0.73%) occurs in a lot of small groups (with a maximum size of 0.05%). We especially studied non-canonical splice sites, which comprise 3.73% of GenBank annotated splice pairs. EST alignments allowed us to verify only the exonic part of splice sites. To check the conservative dinucleotides we compared sequences of human non-canonical splice sites with sequences from the high throughput genome sequencing project (HTG). Out of 171 human non-canonical and EST-supported splice pairs, 156 (91.23%) had a clear match in the human HTG. They can be classified after sequence analysis as: 79 GC–AG pairs (of which one was an error that corrected to GC–AG), 61 errors corrected to GT–AG canonical pairs, six AT–AC pairs (of which two were errors corrected to AT–AC), one case was produced from a non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two other cases left of supported non-canonical splice pairs. The information about verified splice site sequences for canonical and non-canonical sites is presented in SpliceDB with the supporting evidence. We also built weight matrices for the major splice groups, which can be incorporated into gene prediction programs. SpliceDB is available at the computational genomic Web server of the Sanger Centre: http://genomic.sanger.ac.uk/spldb/SpliceDB.html and at http://www.softberry.com/spldb/SpliceDB.html.

The Homeodomain Resource: sequences, structures, DNA binding sites and genomic information

The Homeodomain Resource is an annotated collection of non-redundant protein sequences, three-dimensional structures and genomic information for the homeodomain protein family. Release 3.0 contains 795 full-length homeodomain-containing sequences, 32 experimentally-derived structures and 143 homeo­box loci implicated in human genetic disorders. Entries are fully hyperlinked to facilitate easy retrieval of the original records from source databases. A simple search engine with a graphical user interface is provided to query the component databases and assemble customized data sets. A new feature for this release is the addition of DNA recognition sites for all human homeodomain proteins described in the literature. The Homeodomain Resource is freely available through the World Wide Web at http://genome.nhgri.nih.gov/homeodomain.

Creating a Web-accessible, point-of-care, team-based information system (PoinTIS): the librarian as publisher*

The Internet has created new opportunities for librarians to develop information systems that are readily accessible at the point of care. This paper describes the multiyear process used to justify, fund, design, develop, promote, and evaluate a rehabilitation prototype of a point-of-care, team-based information system (PoinTIS) and train health care providers to use this prototype for their spinal cord injury and traumatic brain injury patient care and education activities. PoinTIS is a successful model for librarians in the twenty-first century to serve as publishers of information created or used by their parent organizations and to respond to the opportunities for information dissemination provided by recent technological advances.

Conversion of a maltose receptor into a zinc biosensor by computational design

We have demonstrated that it is possible to radically change the specificity of maltose binding protein by converting it into a zinc sensor using a rational design approach. In this new molecular sensor, zinc binding is transduced into a readily detected fluorescence signal by use of an engineered conformational coupling mechanism linking ligand binding to reporter group response. An iterative progressive design strategy led to the construction of variants with increased zinc affinity by combining binding sites, optimizing the primary coordination sphere, and exploiting conformational equilibria. Intermediates in the design series show that the adaptive process involves both introduction and optimization of new functions and removal of adverse vestigial interactions. The latter demonstrates the importance of the rational design approach in uncovering cryptic phenomena in protein function, which cannot be revealed by the study of naturally evolved systems.