The PSSH database of alignments between protein sequences and tertiary structures

We introduce the PSSH ('Protein Sequence-to-Structure Homologies') database derived from HSSP2, an improved version of the HSSP ('Homology-derived Secondary Structure of Proteins') database [Dodge et al. (1998) Nucleic Acids Res., 26, 313-315]. Whereas each HSSP entry lists all protein sequences related to a given 3D structure, PSSH is the 'inverse', with each entry listing all structures related to a given sequence. In addition, we introduce two other derived databases: HSSPchain, in which each entry lists all sequences related to a given PDB chain, and HSSPalign, in which each entry gives details of one sequence aligned onto one PDB chain. This re-organization makes it easier to navigate from sequence to structure, and to map sequence features onto 3D structures. Currently (September 2002), PSSH provides structural information for over 400 000 protein sequences, covering 48% of SWALL and 61% of SWISS-PROT sequences; HSSPchain provides sequence information for over 25 000 PDB chains, and HSSPalign gives over 14 million sequence-to-structure alignments. The databases can be accessed via SRS 3D, an extension to the SRS system, at http://srs3d.ebi.ac.uk/.     

Common active site architecture and binding strategy of four phenylpropanoid P450s from Arabidopsis thaliana as revealed by molecular modeling

Despite extensive primary sequence diversity, crystal structures of several bacterial cytochrome P450 monooxygenases (P450s) and a single eukaryotic P450 indicate that these enzymes share a structural core of alpha-helices and beta-sheets and vary in the loop regions contacting individual substrates. To determine the extent to which individual structural features are conserved among divergent P450s existing in a single biosynthetic pathway, we have modeled the structures of four highly divergent P450s (CYP73A5, CYP84A1, CYP75B1, CYP98A3) in the Arabidopsis phenylpropanoid pathway synthesizing lignins, flavonoids and anthocyanins. Analysis of these models has indicated that, despite primary sequence identities as low as 13%, the structural cores and several loop regions of these P450s are highly conserved. Substrate docking indicated that all four enzymes employ a common strategy to identify their substrates in that their cinnamate-derived substrates align along helix I with their aromatic ring positioned towards the C-terminus of this helix and their aliphatic tails positioned towards the N-terminus. Further similarity was observed in the way the substrates contact the consensus P450 substrate recognition sites (SRS). Residues predicted to contact the aromatic ring region exist in SRS5, SRS6 and the C-terminal portion of SRS4 and residues contacting the distal end of each substrate exist in SRS1, SRS2 and the N-terminal portion of SRS4. Alignments of the regions contacting the aromatic ring region indicate that SRS4, SRS5 and SRS6 share higher degrees of sequence conservation than found in SRS1, SRS2 or the full-length protein.     

IMGT-ONTOLOGY for immunogenetics and immunoinformatics

IMGT, the international ImMunoGeneTics information system(R) (http://imgt.cines.fr), is a high quality integrated knowledge resource specializing in immunoglobulins (IG), T cell receptors (TR), major histocompatibility complex (MHC) and related proteins of the immune system (RPI) of human and other vertebrates, created in 1989, by the Laboratoire d'ImmunoGenetique Moleculaire LIGM. IMGT provides a common access to standardized data which include nucleotide and protein sequences, oligonucleotide primers, gene maps, genetic polymorphisms, specificities, 2D and 3D structures. IMGT consists of several sequence databases (IMGT/LIGM-DB, IMGT/MHC-DB, IMGT/PRIMER-DB), one genome database (IMGT/GENE-DB) and one three-dimensional structure database (IMGT/3Dstructure-DB), interactive tools for sequence analysis (IMGT/V-QUEST, IMGT/JunctionAnalysis, IMGT/PhyloGene, IMGT/Allele-Align), for genome analysis (IMGT/GeneSearch, IMGT/GeneView, IMGT/LocusView) and for 3D structure analysis (IMGT/StructuralQuery), and Web resources ("IMGT Marie-Paule page") comprising 8000 HTML pages. IMGT other accesses include SRS, FTP, search by BLAST, etc. By its high quality and its easy data distribution, IMGT has important implications in medical research (repertoire in autoimmune diseases, AIDS, leukemias, lymphomas, myelomas), veterinary research, genome diversity and genome evolution studies of the adaptive immune responses, biotechnology related to antibody engineering (scFv, phage displays, combinatorial libraries) and therapeutical approaches (grafts, immunotherapy). IMGT is freely available at http://imgt.cines.fr.     

Single mutations change CYP2F3 from a dehydrogenase of 3-methylindole to an oxygenase

Pulmonary cytochrome P450 2F3 (CYP2F3) catalyzes the dehydrogenation of the pneumotoxin 3-methylindole (3MI) to an electrophilic intermediate, 3-methyleneindolenine, which is responsible for the toxicity of the parent compound. Members of the CYP2F subfamily are the only enzymes known to exclusively dehydrogenate 3MI, without detectable formation of oxygenation products. Thus, CYP2F3 is an attractive model to study dehydrogenation mechanisms. The purpose of this study was to identify specific residues that could facilitate 3MI dehydrogenation. Both single and double mutations were constructed to study the molecular mechanisms that direct dehydrogenation. Double mutations in substrate recognition sites (SRS) 1 produced an inactive enzyme, while double mutants in SRS 4 did not alter 3MI metabolism. However, double mutations in SRS 5 and SRS 6 successfully introduced oxygenase activity to CYP2F3. Single mutations in SRS 5, SRS 6 and near SRS 2 also introduced 3MI oxygenase activity. Mutants S474H and D361T oxygenated 3MI but also increased dehydrogenation rates, while G214L, E215Q and S475I catalyzed 3MI oxygenation exclusively. A homology model of CYP2F3 was precisely consistent with specific dehydrogenation of 3MI via initial hydrogen atom abstraction from the methyl group. In addition, intramolecular kinetic deuterium isotope studies demonstrated an isotope effect ( K H/ K D) of 6.8. This relatively high intramolecular deuterium isotope effect confirmed the initial hydrogen abstraction step; a mutant (D361T) that retained the dehydrogenation reaction exhibited the same deuterium isotope effect. The results showed that a single alteration, such as a serine to isoleucine change at residue 475, dramatically switched catalytic preference from dehydrogenation to oxygenation.     

Survival in Nuclear Waste,Extreme Resistance,and Potential Applications Gleaned from the Genome Sequence of Kineococcus radiotolerans SRS30216

Kineococcus radiotolerans {"type":"entrez-protein","attrs":{"text":"SRS30216","term_id":"1414251709","term_text":"SRS30216"}}SRS30216 was isolated from a high-level radioactive environment at the Savannah River Site (SRS) and exhibits γ-radiation resistance approaching that of Deinococcus radiodurans. The genome was sequenced by the U.S. Department of Energy''s Joint Genome Institute which suggested the existence of three replicons, a 4.76 Mb linear chromosome, a 0.18 Mb linear plasmid, and a 12.92 Kb circular plasmid. Southern hybridization confirmed that the chromosome is linear. The K. radiotolerans genome sequence was examined to learn about the physiology of the organism with regard to ionizing radiation resistance, the potential for bioremediation of nuclear waste, and the dimorphic life cycle. K. radiotolerans may have a unique genetic toolbox for radiation protection as it lacks many of the genes known to confer radiation resistance in D. radiodurans. Additionally, genes involved in the detoxification of reactive oxygen species and the excision repair pathway are overrepresented. K. radiotolerans appears to lack degradation pathways for pervasive soil and groundwater pollutants. However, it can respire on two organic acids found in SRS high-level nuclear waste, formate and oxalate, which promote the survival of cells during prolonged periods of starvation. The dimorphic life cycle involves the production of motile zoospores. The flagellar biosynthesis genes are located on a motility island, though its regulation could not be fully discerned. These results highlight the remarkable ability of K radiotolerans to withstand environmental extremes and suggest that in situ bioremediation of organic complexants from high level radioactive waste may be feasible.     

BLAST2SRS,a web server for flexible retrieval of related protein sequences in the SWISS-PROT and SPTrEMBL databases

SRS (Sequence Retrieval System) is a widely used keyword search engine for querying biological databases. BLAST2 is the most widely used tool to query databases by sequence similarity search. These tools allow users to retrieve sequences by shared keyword or by shared similarity, with many public web servers available. However, with the increasingly large datasets available it is now quite common that a user is interested in some subset of homologous sequences but has no efficient way to restrict retrieval to that set. By allowing the user to control SRS from the BLAST output, BLAST2SRS (http://blast2srs.embl.de/) aims to meet this need. This server therefore combines the two ways to search sequence databases: similarity and keyword.     

Identification of unique amino acids that modulate CYP4A7 activity

A multifamily sequence alignment of the rabbit CYP4A members with the known structure of CYP102 indicates amino acid differences falling within the so-called substrate recognition site(s) (SRS). Chimeric proteins constructed between CYP4A4 and CYP4A7 indicate that laurate activity is affected by the residues within SRS1 and prostaglandin activity is influenced by SRS2-3. Site-directed mutant proteins of CYP4A7 found laurate and arachidonate activity markedly diminished in the R90W mutant (SRS1) and somewhat decreased in W93S. While PGE(1) activity was only slightly increased, the mutant proteins H206Y and S255F (SRS2-3), on the other hand, exhibited remarkable increases in laurate and arachidonate metabolism (3-fold) above wild-type substrate metabolism. Mutant proteins H206Y, S255F, and H206Y/S255F but not R90W/W93S, wild-type CYP4A4, or CYP4A7 metabolized arachidonic acid in the absence of cytochrome b(5). Stopped-flow kinetic experiments were performed in a CO-saturated environment performed to estimate interaction rates of the monooxygenase reaction components. The mutant protein H206Y, which exhibits 3-fold higher than wild-type substrate activity, interacts with CPR at a rate at least 10 times faster than that of wild-type CYP4A7. These experimental results provide insight regarding the residues responsible for modulation of substrate specificity, affinity, and kinetics, as well as possible localization within the enzyme structure based on comparisons with homologous, known cytochrome P450 structures.     

PROMALS3D: a tool for multiple protein sequence and structure alignments

Although multiple sequence alignments (MSAs) are essential for a wide range of applications from structure modeling to prediction of functional sites, construction of accurate MSAs for distantly related proteins remains a largely unsolved problem. The rapidly increasing database of spatial structures is a valuable source to improve alignment quality. We explore the use of 3D structural information to guide sequence alignments constructed by our MSA program PROMALS. The resulting tool, PROMALS3D, automatically identifies homologs with known 3D structures for the input sequences, derives structural constraints through structure-based alignments and combines them with sequence constraints to construct consistency-based multiple sequence alignments. The output is a consensus alignment that brings together sequence and structural information about input proteins and their homologs. PROMALS3D can also align sequences of multiple input structures, with the output representing a multiple structure-based alignment refined in combination with sequence constraints. The advantage of PROMALS3D is that it gives researchers an easy way to produce high-quality alignments consistent with both sequences and structures of proteins. PROMALS3D outperforms a number of existing methods for constructing multiple sequence or structural alignments using both reference-dependent and reference-independent evaluation methods.     

3dswap-pred: prediction of 3D domain swapping from protein sequence using Random Forest approach

3D domain swapping is a protein structural phenomenon that mediates the formation of the higher order oligomers in a variety of proteins with different structural and functional properties. 3D domain swapping is associated with a variety of biological functions ranging from oligomerization to pathological conformational diseases. 3D domain swapping is realised subsequent to structure determination where the protein is observed in the swapped conformation in the oligomeric state. This is a limiting step to understand this important structural phenomenon in a large scale from the growing sequence data. A new machine learning approach, 3dswap-pred, has been developed for the prediction of 3D domain swapping in protein structures from mere sequence data using the Random Forest approach. 3Dswap-pred is implemented using a positive sequence dataset derived from literature based structural curation of 297 structures. A negative sequence dataset is obtained from 462 SCOP domains using a new sequence data mining approach and a set of 126 sequencederived features. Statistical validation using an independent dataset of 68 positive sequences and 313 negative sequences revealed that 3dswap-pred achieved an accuracy of 63.8%. A webserver is also implemented using the 3dswap-pred Random Forest model. The server is available from the URL: http://caps.ncbs.res.in/3dswap-pred.     

TESS: a geometric hashing algorithm for deriving 3D coordinate templates for searching structural databases. Application to enzyme active sites.

It is well established that sequence templates such as those in the PROSITE and PRINTS databases are powerful tools for predicting the biological function and tertiary structure for newly derived protein sequences. The number of X-ray and NMR protein structures is increasing rapidly and it is apparent that a 3D equivalent of the sequence templates is needed. Here, we describe an algorithm called TESS that automatically derives 3D templates from structures deposited in the Brookhaven Protein Data Bank. While a new sequence can be searched for sequence patterns, a new structure can be scanned against these 3D templates to identify functional sites. As examples, 3D templates are derived for enzymes with an O-His-O "catalytic triad" and for the ribonucleases and lysozymes. When these 3D templates are applied to a large data set of nonidentical proteins, several interesting hits are located. This suggests that the development of a 3D template database may help to identify the function of new protein structures, if unknown, as well as to design proteins with specific functions.     

The Medium-Chain Dehydrogenase/reductase Engineering Database: a systematic analysis of a diverse protein family to understand sequence-structure-function relationship

The Medium-Chain Dehydrogenase/Reductase Engineering Database (MDRED, http://www.mdred.uni-stuttgart.de) has been established to serve as an analysis tool for a systematic investigation of sequence-structure-function relationships. It includes sequence and structure information of 2684 and 42 medium-chain dehydrogenases/reductases (MDRs), respectively. Although MDRs are very diverse in sequence, they have a conserved tertiary structure. MDRs are assigned to 199 homologous families and 29 superfamilies. For each family, annotated multiple sequence alignments are provided, and functionally relevant residues are annotated. Twenty-five superfamilies were classified as zinc-containing MDRs, four as non-zinc-containing MDRs. For the zinc-containing MDRs, three subclasses were identified by systematic analysis of a variable loop region, the quaternary structure determining loop (QSDL): the class of short, medium, and long QSDL, which include 11, 3, and 5 superfamilies, respectively. The length of the QSDL is predictive for tetramer (short QSDL) and dimer (long QSDL) formation. The class of medium QSDL includes both tetrameric and dimeric MDRs. The shape of the substrate-binding site is highly conserved in all zinc-containing MDRs with the exception of two variable regions, the substrate recognition sites (SRS): two residues located on the QSDL (SRS1) and, for the class of long QSDL, one residue located in the catalytic domain (SRS2). The MDRED is the first online-accessible resource of MDRs that integrates information on sequence, structure, and function. Annotation of functionally relevant residues assist the understanding of sequence-structure-function relationships. Thus, the MDRED serves as a valuable tool to identify potential hotspots for engineering properties such as substrate specificity.     

3D-Printed high-NA catadioptric thin lens for suppression of XPM background in Stimulated Raman Scattering microscopy

Stimulated Raman Scattering (SRS) is a fast chemical imaging technique with remarkable bioscience applications. Cross Phase Modulation (XPM) is a ubiquitous nonlinear phenomenon that can create spurious background signals that render difficult a high-contrast imaging in SRS measurements. The XPM-induced signal is usually suppressed using high numerical aperture (NA) microscope objectives or condensers to collect the transmitted excitation beam. However, these high NA optics feature short working distances, hence they are not compatible with stage-top incubators, that are necessary to perform live-cell time-lapse experiments in controlled environments. Here, we show a 3D printed high NA compact catadioptric lens that fits inside stage-top incubators and allows the collection of XPM-free SRS signals. The lens delivers SRS images and spectra with a quality comparable to a signal collection with a high-NA microscope objective. We also demonstrate the compatibility of the 3D printed lens with other nonlinear microscopies usually associated with SRS in multimodal microscopes.     

Use of a structural alphabet to find compatible folds for amino acid sequences

The structural annotation of proteins with no detectable homologs of known 3D structure identified using sequence‐search methods is a major challenge today. We propose an original method that computes the conditional probabilities for the amino‐acid sequence of a protein to fit to known protein 3D structures using a structural alphabet, known as “Protein Blocks” (PBs). PBs constitute a library of 16 local structural prototypes that approximate every part of protein backbone structures. It is used to encode 3D protein structures into 1D PB sequences and to capture sequence to structure relationships. Our method relies on amino acid occurrence matrices, one for each PB, to score global and local threading of query amino acid sequences to protein folds encoded into PB sequences. It does not use any information from residue contacts or sequence‐search methods or explicit incorporation of hydrophobic effect. The performance of the method was assessed with independent test datasets derived from SCOP 1.75A. With a Z‐score cutoff that achieved 95% specificity (i.e., less than 5% false positives), global and local threading showed sensitivity of 64.1% and 34.2%, respectively. We further tested its performance on 57 difficult CASP10 targets that had no known homologs in PDB: 38 compatible templates were identified by our approach and 66% of these hits yielded correctly predicted structures. This method scales‐up well and offers promising perspectives for structural annotations at genomic level. It has been implemented in the form of a web‐server that is freely available at http://www.bo‐protscience.fr/forsa .     

Impact of incorporating the 2C5 crystal structure into comparative models of cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) metabolizes approximately one third of the drugs in current clinical use. To gain insight into its structure and function, we have produced four different sets of comparative models of 2D6: one based on the structures of P450s from four different microorganisms (P450 terp, P450 eryF, P450 cam, and P450 BM3), another on the only mammalian P450 (2C5) structure available, and the other two based on alternative amino acid sequence alignments of 2D6 with all five of these structures. Principal component analysis suggests that inclusion of the 2C5 crystal structure has a profound effect on the modeling process, altering the general topology of the active site, and that the models produced differ significantly from all of the templates. The four models of 2D6 were also used in conjunction with molecular docking to produce complexes with the substrates codeine and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP); this identified Glu 216 [in the F-helix; substrate recognition site (SRS) 2] as a key determinant in the binding of the basic moiety of the substrate. Our studies suggest that both Asp 301 and Glu 216 are required for metabolism of basic substrates. Furthermore, they suggest that Asp 301 (I-helix, SRS-4), a residue thought from mutagenesis studies to bind directly to the basic moiety of substrates, may play a key role in positioning the B'-C loop (SRS-1) and that the loss of activity on mutating Asp 301 may therefore be the result of an indirect effect (movement of the B'-C loop) on replacing this residue.     

The SRS superfamily of Toxoplasma surface proteins

The surface of the protozoan parasite Toxoplasma gondii is coated with developmentally expressed, glycosylphosphatidylinositol-linked proteins structurally related to the highly immunogenic surface antigen SAG1. Collectively, these surface antigens are known as the SRS (SAG1-related sequences) superfamily of proteins. SRS proteins are thought to mediate attachment to host cells and activate host immunity to regulate the parasite's virulence. To better understand the number, evolution and developmental expression of SRS genes, this study has bioinformatically identified 161 unique SRS DNA sequences present in the T. gondii type II Me49 genome. The SRS superfamily of sequences phylogenetically bifurcates into two subfamilies, the prototypic members being SAG1 and SAG2A, respectively. Paralogous SRS sequences are 24-99% identical, are tandemly arrayed throughout the genome, and are present on most, if not all, chromosomes. All 11 SRS sequences on chromosomes Ia and Ib are clustered at sub-telomeric expression sites. Messenger RNA expression in the majority of SRS sequences for which multiple Expressed Sequence Tags exist is developmentally regulated. A consensus nucleotide sequence surrounding both the splice acceptor and donor sites was identified in those SRS sequences possessing an intron. Genotypic differences among SRS sequences are present at several loci (e.g. the absence of SAG5B, the truncation of SAG2D in Me49 compared with RH) indicating that different genotypes possess distinct sets of SRS sequences. Orthologous genes are restricted to tissue-dwelling coccidia (Neospora, Sarcocystis) with no related sequences present in other more distant apicomplexa such as Eimeria, Cryptosporidia, and Plasmodium spp.     

The EBI SRS server-new features

MOTIVATION: Here we report on recent developments at the EBI SRS server (http://srs.ebi.ac.uk). SRS has become an integration system for both data retrieval and sequence analysis applications. The EBI SRS server is a primary gateway to major databases in the field of molecular biology produced and supported at EBI as well as European public access point to the MEDLINE database provided by US National Library of Medicine (NLM). It is a reference server for latest developments in data and application integration. The new additions include: concept of virtual databases, integration of XML databases like the Integrated Resource of Protein Domains and Functional Sites (InterPro), Gene Ontology (GO), MEDLINE, Metabolic pathways, etc., user friendly data representation in 'Nice views', SRSQuickSearch bookmarklets. AVAILABILITY: SRS6 is a licensed product of LION Bioscience AG freely available for academics. The EBI SRS server (http://srs.ebi.ac.uk) is a free central resource for molecular biology data as well as a reference server for the latest developments in data integration.     

Molecular evolution of ion channels: Amino acid sequences and 3D structures

Comparative analysis of structure and function of macromolecules, such as proteins, is an integral part of modern evolutionary biology. The first and critical step in understanding evolution of homologous proteins is their amino acid sequence alignment. However, standard algorithms fail to provide unambiguous sequence alignment for proteins of poor homology. More reliable results can be provided by comparing experimental 3D structures obtained at atomic resolution with the aid of X-ray structural analysis. If such structures are lacking, homology modeling is used which considers indirect experimental data on functional roles of individual amino acid residues. An important problem is that sequence alignment, which reflects genetic modifications, not necessarily corresponds to functional homology, which depends on 3D structures critical for natural selection. Since the alignment techniques relying only on the analysis of primary structures carry no information on the functional properties of proteins, the inclusion of 3D structures into consideration is of utmost importance. Here we consider several ion channels as examples to demonstrate that alignment of their 3D structures can significantly improve sequence alignment obtained by traditional methods.     

苔藓动物主要类群28S rDNA多变区的分子形态分析及其系统学意义

测定了苔藓动物主要代表类群3种苔虫的28S rDNA全序列(3 459～3 613 bp),结合已公布的相关类群的28S rDNA全序列数据,分析了其一级结构特征,并分别绘制了D2、D3和D8区段的二级结构.比较分析结果显示:各类群的全序列总长均超过3 300 bp,GC含量为41.8%～57.7%;D3和D8区段的二级结构保守性较高,D2区差异显著;各区段的长度变化较大,但是,大多都具有相似的主干结构.基于二级结构多变区分子形态学特征的比较基本和基于一级结构的分子系统学结论相吻合,提示它们同样包含着重要的系统发育信息.     

SpoIIIE mechanism of directional translocation involves target search coupled to sequence‐dependent motor stimulation

SpoIIIE/FtsK are membrane‐anchored, ATP‐fuelled, directional motors responsible for chromosomal segregation in bacteria. Directionality in these motors is governed by interactions between specialized sequence‐recognition modules (SpoIIIE‐γ/FtsK‐γ) and highly skewed chromosomal sequences (SRS/KOPS). Using a new combination of ensemble and single‐molecule methods, we dissect the series of steps required for SRS localization and motor activation. First, we demonstrate that SpoIIIE/DNA association kinetics are sequence independent, with binding specificity being uniquely determined by dissociation. Next, we show by single‐molecule and modelling methods that hexameric SpoIIIE binds DNA non‐specifically and finds SRS by an ATP‐independent target search mechanism, with ensuing oligomerization and binding of SpoIIIE‐γ to SRS triggering motor stimulation. Finally, we propose a new model that provides an entirely new interpretation of previous observations for the origin of SRS/KOPS‐directed translocation by SpoIIIE/FtsK.