STING Millennium Suite: integrated software for extensive analyses of 3d structures of proteins and their complexes

DistAA: Database of amino acid distances in proteins and web application for statistical review of distances

Three-dimensional structure of a protein chain is determined by its amino acid interactions. One approach to the analysis of amino acid interactions refers to geometric distances of amino acid pairs in polypeptide chains. For a detailed analysis of the amino acid distances, the database with three types of amino acid distances in a set of chains was created. Web application Distances of Amino Acids has also been developed to enable scientists to explore interactions of amino acids with different properties based on distances stored in the database. Web application calculates and displays descriptive statistics and graphs of amino acid pair distances with selected properties, such as geometric distance threshold, corresponding SCOP class of proteins and secondary structure types. In addition to the analysis of pre-calculated distances stored in the database, the amino acid distances of a single protein with the specified PDB identifier can also be analyzed. The web application is available at http://andromeda.matf.bg.ac.rs/aadis_dynamic/.

De novo assembly and characterization of the Trichuris trichiura adult worm transcriptome using Ion Torrent sequencing

Infection with helminthic parasites, including the soil-transmitted helminth Trichuris trichiura (human whipworm), has been shown to modulate host immune responses and, consequently, to have an impact on the development and manifestation of chronic human inflammatory diseases. De novo derivation of helminth proteomes from sequencing of transcriptomes will provide valuable data to aid identification of parasite proteins that could be evaluated as potential immunotherapeutic molecules in near future. Herein, we characterized the transcriptome of the adult stage of the human whipworm T. trichiura, using next-generation sequencing technology and a de novo assembly strategy. Nearly 17.6 million high-quality clean reads were assembled into 6414 contiguous sequences, with an N50 of 1606bp. In total, 5673 protein-encoding sequences were confidentially identified in the T. trichiura adult worm transcriptome; of these, 1013 sequences represent potential newly discovered proteins for the species, most of which presenting orthologs already annotated in the related species T. suis. A number of transcripts representing probable novel non-coding transcripts for the species T. trichiura were also identified. Among the most abundant transcripts, we found sequences that code for proteins involved in lipid transport, such as vitellogenins, and several chitin-binding proteins. Through a cross-species expression analysis of gene orthologs shared by T. trichiura and the closely related parasites T. suis and T. muris it was possible to find twenty-six protein-encoding genes that are consistently highly expressed in the adult stages of the three helminth species. Additionally, twenty transcripts could be identified that code for proteins previously detected by mass spectrometry analysis of protein fractions of the whipworm somatic extract that present immunomodulatory activities. Five of these transcripts were amongst the most highly expressed protein-encoding sequences in the T. trichiura adult worm. Besides, orthologs of proteins demonstrated to have potent immunomodulatory properties in related parasitic helminths were also predicted from the T. trichiura de novo assembled transcriptome.

Bioinformatic analysis of pathogenic missense mutations of activin receptor like kinase 1 ectodomain

Activin A receptor, type II-like kinase 1 (also called ALK1), is a serine-threonine kinase predominantly expressed on endothelial cells surface. Mutations in its ACVRL1 encoding gene (12q11-14) cause type 2 Hereditary Haemorrhagic Telangiectasia (HHT2), an autosomal dominant multisystem vascular dysplasia. The study of the structural effects of mutations is crucial to understand their pathogenic mechanism. However, while an X-ray structure of ALK1 intracellular domain has recently become available (PDB ID: 3MY0), structure determination of ALK1 ectodomain (ALK1_EC) has been elusive so far. We here describe the building of a homology model for ALK1_EC, followed by an extensive bioinformatic analysis, based on a set of 38 methods, of the effect of missense mutations at the sequence and structural level. ALK1_EC potential interaction mode with its ligand BMP9 was then predicted combining modelling and docking data. The calculated model of the ALK1_EC allowed mapping and a preliminary characterization of HHT2 associated mutations. Major structural changes and loss of stability of the protein were predicted for several mutations, while others were found to interfere mainly with binding to BMP9 or other interactors, like Endoglin (CD105), whose encoding ENG gene (9q34) mutations are known to cause type 1 HHT. This study gives a preliminary insight into the potential structure of ALK1_EC and into the structural effects of HHT2 associated mutations, which can be useful to predict the potential effect of each single mutation, to devise new biological experiments and to interpret the biological significance of new mutations, private mutations, or non-synonymous polymorphisms.

Straight-chain alkyl isocyanides open the distal histidine gate in crystal structures of myoglobin

Crystal structures of methyl, ethyl, propyl, and butyl isocyanide bound to sperm whale myoglobin (Mb) reveal two major conformations. In the in conformer, His(E7) is in a "closed" position, forcing the ligand alkyl chain to point inward. In the out conformer, His(E7) is in an "open" position, allowing the ligand side chain to point outward. A progressive increase in the population of the out conformer is observed with increasing ligand length in P2(1) crystals of native Mb at pH 7.0. This switch from in to out with increasing ligand size also occurs in solution as measured by the decrease in the relative intensity of the low-frequency ( approximately 2075 cm(-1)) versus high-frequency ( approximately 2125 cm(-1)) isocyano bands. In contrast, all four isocyanides in P6 crystals of wild-type recombinant Mb occupy the in conformation. However, mutating either His64 to Ala, creating a "hole" to solvent, or Phe46 to Val, freeing rotation of His64, causes bound butyl isocyanide to point completely outward in P6 crystals. Thus, the unfavorable hindrance caused with crowding a large alkyl side chain into the distal pocket appears to be roughly equal to that for pushing open the His(E7) gate and is easily affected by crystal packing. This structural conclusion supports the "side path" kinetic mechanism for O(2) release, in which the dissociated ligand first moves toward the protein interior and then encounters steric resistance, which is roughly equal to that for escaping to solvent through the His(E7) channel.

Actin-interacting and flagellar proteins in Leishmania spp.: Bioinformatics predictions to functional assignments in phagosome formation

Several motile processes are responsible for the movement of proteins into and within the flagellar membrane, but little is known about the process by which specific proteins (either actin-associated or not) are targeted to protozoan flagellar membranes. Actin is a major cytoskeleton protein, while polymerization and depolymerization of parasite actin and actin-interacting proteins (AIPs) during both processes of motility and host cell entry might be key events for successful infection. For a better understanding the eukaryotic flagellar dynamics, we have surveyed genomes, transcriptomes and proteomes of pathogenic Leishmania spp. to identify pertinent genes/proteins and to build in silico models to properly address their putative roles in trypanosomatid virulence. In a search for AIPs involved in flagellar activities, we applied computational biology and proteomic tools to infer from the biological meaning of coronins and Arp2/3, two important elements in phagosome formation after parasite phagocytosis by macrophages. Results presented here provide the first report of Leishmania coronin and Arp2/3 as flagellar proteins that also might be involved in phagosome formation through actin polymerization within the flagellar environment. This is an issue worthy of further in vitro examination that remains now as a direct, positive bioinformatics-derived inference to be presented.

Characterization of the genome sequence of an oncolytic Newcastle disease virus strain Italien

We determined the complete genome sequence of strain Italien, a virulent and oncolytic strain of Newcastle disease virus (NDV) by direct nucleotide sequencing of RT-PCR products, a size of 15,186 nucleotides (nt). Comparison of six coding genes and non-coding regions of Italien with those of the other 25 sequenced strains revealed NDV Herts/33 was the most similar strain with Italien. The gene encoding the RNA dependent RNA polymerase was the most highly conserved, while the gene encoding phosphoprotein was the most highly variable. The HN and F proteins of Italien have been modeled on the crystal structure in order to study the structural characteristics. Interaction between the HN protein and the heptad repeat B (HRB) region of F protein was analyzed in silico by molecular docking predicted five critical residues I133, V142, D143, R480, and K567 on HN protein. Identification of amino acid residues that could be crucial for this interaction provides working hypotheses for subsequent studies.

Tools for integrated sequence-structure analysis with UCSF Chimera

Background

Comparing related structures and viewing the structures in the context of sequence alignments are important tasks in protein structure-function research. While many programs exist for individual aspects of such work, there is a need for interactive visualization tools that: (a) provide a deep integration of sequence and structure, far beyond mapping where a sequence region falls in the structure and vice versa; (b) facilitate changing data of one type based on the other (for example, using only sequence-conserved residues to match structures, or adjusting a sequence alignment based on spatial fit); (c) can be used with a researcher's own data, including arbitrary sequence alignments and annotations, closely or distantly related sets of proteins, etc.; and (d) interoperate with each other and with a full complement of molecular graphics features. We describe enhancements to UCSF Chimera to achieve these goals.

Results

The molecular graphics program UCSF Chimera includes a suite of tools for interactive analyses of sequences and structures. Structures automatically associate with sequences in imported alignments, allowing many kinds of crosstalk. A novel method is provided to superimpose structures in the absence of a pre-existing sequence alignment. The method uses both sequence and secondary structure, and can match even structures with very low sequence identity. Another tool constructs structure-based sequence alignments from superpositions of two or more proteins. Chimera is designed to be extensible, and mechanisms for incorporating user-specific data without Chimera code development are also provided.

Conclusion

The tools described here apply to many problems involving comparison and analysis of protein structures and their sequences. Chimera includes complete documentation and is intended for use by a wide range of scientists, not just those in the computational disciplines. UCSF Chimera is free for non-commercial use and is available for Microsoft Windows, Apple Mac OS X, Linux, and other platforms from .

Effects of 3-methylcholanthrene on gene expression profiling in the rat using cDNA microarray analyses

There is significant human exposure to polycyclic aromatic hydrocarbons (PAHs), many of which are bioactivated by the cytochrome P450 (P450) 1A family of enzymes to metabolites that are capable of covalently binding to DNA, a critical step in the initiation of carcinogenesis. We reported earlier that exposure of rats to 3-methylcholanthrene (MC) causes sustained induction of hepatic cytochrome P4501A expression for up to 45 days. Here, we tested the hypothesis that MC elicits persistent induction of other genes that are regulated by the Ah receptor (AHR). Female Sprague-Dawley rats were treated with MC (100 micromol/kg) ip once daily for 4 days, and gene expression patterns were investigated using total liver RNA isolated from animals at 1, 15, and 28 days after MC withdrawal. Gene expression was studied by cDNA microarray analyses using 4608 unique clones from liver-derived expressed sequence tag (EST) libraries fortified with clones of known liver genes representing approximately 4000 genes. Several phase I (P4501A1, -1A2) and phase II [e.g., glutathione-S-transferase (GST)-M1, UDP-glucuronosyl transferases (UGT)] genes were persistently induced (3-10-fold) by MC for 15-28 days. The persistent induction of P4501A1 gene expression was confirmed by real time reverse transcriptase polymerase chain reaction (RT-PCR) experiments. MC also elicited a 5-fold persistent augmentation of acute phase genes such as orosomucoid 1 and alpha-1-acid glycoprotein (AGP), and this was accompanied by sustained liver damage and inflammation in the MC-exposed rats. In conclusion, our results strongly suggest that sustained induction of P4501A1 by MC is accompanied by persistent expression of other genes belonging to the Ah gene battery, as well as certain other genes involved in toxic responses. Elucidating the mechanisms of persistent induction of P4501A1 and other genes by MC might lead to a better understanding of the mechanisms of toxicity mediated by PAHs.

Pathogenomics analysis of Leishmania spp.: flagellar gene families of putative virulence factors

The trypanosomatid flagellar apparatus contains conventional and unique features, whose roles in infectivity are still enigmatic. Although the flagellum and the flagellar pocket are critical organelles responsible for all vesicular trafficking between the cytoplasm and cell surface, still very little is known about their roles in pathogenesis and how molecules get to and from the flagellar pocket. The ongoing analysis of the genome sequences and proteome profiles of Leishmania major and L infantum, Trypanosoma cruzi, T. brucei, and T. gambiensi ( www.genedb.org ), coupled with our own work on L. chagasi (as part of the Brazilian Northeast Genome Program- www.progene.ufpe.br ), prompted us to scrutinize flagellar genes and proteins of Leishmania spp. promastigotes that could be virulence factors in leishmaniasis. We have identified some overlooked parasite factors such as the MNUDC-1 (a protein involved in nuclear development and genomic fusion) and SQS (an enzyme of sterol biosynthesis), among the described flagellar gene families. A database concerning the results of this work, as well as of other studies of Leishmania and its organelles, is available at http://nugen.lcc.uece.br/LPGate . It will serve as a convenient bioinformatics resource on genomics and pathology of the etiological agents of leishmaniasis.

A molecular docking model of SARS-CoV S1 protein in complex with its receptor, human ACE2

The exact residues within severe acute respiratory syndrome coronavirus (SARS-CoV) S1 protein and its receptor, human ACE2, involved in their interaction still remain largely undetermined. Identification of exact amino acid residues that are crucial for the interaction of S1 with ACE2 could provide working hypotheses for experimental studies and might be helpful for the development of antiviral inhibitor. In this paper, a molecular docking model of SARS-CoV S1 protein in complex with human ACE2 was constructed. The interacting residue pairs within this complex model and their contact types were also identified. Our model, supported by significant biochemical evidence, suggested receptor-binding residues were concentrated in two segments of S1 protein. In contrast, the interfacial residues in ACE2, though close to each other in tertiary structure, were found to be widely scattered in the primary sequence. In particular, the S1 residue ARG453 and ACE2 residue LYS341 might be the key residues in the complex formation.

The Diamond STING server

Diamond STING is a new version of the STING suite of programs for a comprehensive analysis of a relationship between protein sequence, structure, function and stability. We have added a number of new functionalities by both providing more structure parameters to the STING Database and by improving/expanding the interface for enhanced data handling. The integration among the STING components has also been improved. A new key feature is the ability of the STING server to handle local files containing protein structures (either modeled or not yet deposited to the Protein Data Bank) so that they can be used by the principal STING components: (Java)Protein Dossier ((J)PD) and STING Report. The current capabilities of the new STING version and a couple of biologically relevant applications are described here. We have provided an example where Diamond STING identifies the active site amino acids and folding essential amino acids (both previously determined by experiments) by filtering out all but those residues by selecting the numerical values/ranges for a set of corresponding parameters. This is the fundamental step toward a more interesting endeavor-the prediction of such residues. Diamond STING is freely accessible at http://sms.cbi.cnptia.embrapa.br and http://trantor.bioc.columbia.edu/SMS.