Parmodel: a web server for automated comparative modeling of proteins

Parmodel is a web server for automated comparative modeling and evaluation of protein structures. The aim of this tool is to help inexperienced users to perform modeling, assessment, visualization, and optimization of protein models as well as crystallographers to evaluate structures solved experimentally. It is subdivided in four modules: Parmodel Modeling, Parmodel Assessment, Parmodel Visualization, and Parmodel Optimization. The main module is the Parmodel Modeling that allows the building of several models ford a same protein in a reduced time, through the distribution of modeling processes on a Beowulf cluster. Parmodel automates and integrates the main softwares used in comparative modeling as MODELLER, Whatcheck, Procheck, Raster3D, Molscript, and Gromacs. This web server is freely accessible at http://www.biocristalografia.df.ibilce.unesp.br/tools/parmodel. (C) 2004 Elsevier B.V. All rights reserved.

Biomechanics of the rostrum in crocodilians: A comparative analysis using finite-element modeling

This article reports the use of simple beam and finite-element models to investigate the relationship between rostral shape and biomechanical performance in living crocodilians under a range of loading conditions. Load cases corresponded to simple biting, lateral head shaking, and twist feeding behaviors. The six specimens were chosen to reflect, as far as possible, the full range of rostral shape in living crocodilians: a juvenile Caiman crocodilus, subadult Alligator mississippiensis and Crocodylus johnstoni, and adult Caiman crocodilus, Melanosuchus niger, and Paleosuchus palpebrosus. The simple beam models were generated using morphometric landmarks from each specimen. Three of the finite-element models, the A. mississippiensis, juvenile Caiman crocodilus, and the Crocodylus johnstoni, were based on CT scan data from respective specimens, but these data were not available for the other models and so these-the adult Caiman crocodilus, M. niger, and P. palpebrosus-were generated by morphing the juvenile Caiman crocodilus mesh with reference to three-dimensional linear distance measured from specimens. Comparison of the mechanical performance of the six finite-element models essentially matched results of the simple beam models: relatively tall skulls performed best under vertical loading and tall and wide skulls performed best under torsional loading. The widely held assumption that the platyrostral (dorsoventrally flattened) crocodilian skull is optimized for torsional loading was not supported by either simple beam theory models or finite-element modeling. Rather than being purely optimized against loads encountered while subduing and processing food, the shape of the crocodilian rostrum may be significantly affected by the hydrodynamic constraints of catching agile aquatic prey. This observation has important implications for our understanding of biomechanics in crocodilians and other aquatic reptiles.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

Characterization of the Entamoeba histolytica Ornithine Decarboxylase-Like Enzyme

Background: The polyamines putrescine, spermidine, and spermine are organic cations that are required for cell growth and differentiation. Ornithine decarboxylase (ODC), the first and rate-limiting enzyme in the polyamine biosynthetic pathway, is a highly regulated enzyme. Methodology and Results: To use this enzyme as a potential drug target, the gene encoding putative ornithine decarboxylase (ODC)-like sequence was cloned from Entamoeba histolytica, a protozoan parasite causing amoebiasis. DNA sequence analysis revealed an open reading frame (ORF) of similar to 1,242 bp encoding a putative protein of 413 amino acids with a calculated molecular mass of 46 kDa and a predicted isoelectric point of 5.61. The E. histolytica putative ODC-like sequence has 33% sequence identity with human ODC and 36% identity with the Datura stramonium ODC. The ORF is a single-copy gene located on a 1.9-Mb chromosome. The recombinant putative ODC protein (48 kDa) from E. histolytica was heterologously expressed in Escherichia coli. Antiserum against recombinant putative ODC protein detected a band of anticipated size similar to 46 kDa in E. histolytica whole-cell lysate. Difluoromethylornithine (DFMO), an enzyme-activated irreversible inhibitor of ODC, had no effect on the recombinant putative ODC from E. histolytica. Comparative modeling of the three-dimensional structure of E. histolytica putative ODC shows that the putative binding site for DFMO is disrupted by the substitution of three amino acids-aspartate-332, aspartate-361, and tyrosine-323-by histidine-296, phenylalanine-305, and asparagine-334, through which this inhibitor interacts with the protein. Amino acid changes in the pocket of the E. histolytica enzyme resulted in low substrate specificity for ornithine. It is possible that the enzyme has evolved a novel substrate specificity. Conclusion: To our knowledge this is the first report on the molecular characterization of putative ODC-like sequence from E. histolytica. Computer modeling revealed that three of the critical residues required for binding of DFMO to the ODC enzyme are substituted in E. histolytica, resulting in the likely loss of interactions between the enzyme and DFMO.

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.

Weak conservation of structural features in the interfaces of homologous transient protein-protein complexes

Residue types at the interface of protein-protein complexes (PPCs) are known to be reasonably well conserved. However, we show, using a dataset of known 3-D structures of homologous transient PPCs, that the 3-D location of interfacial residues and their interaction patterns are only moderately and poorly conserved, respectively. Another surprising observation is that a residue at the interface that is conserved is not necessarily in the interface in the homolog. Such differences in homologous complexes are manifested by substitution of the residues that are spatially proximal to the conserved residue and structural differences at the interfaces as well as differences in spatial orientations of the interacting proteins. Conservation of interface location and the interaction pattern at the core of the interfaces is higher than at the periphery of the interface patch. Extents of variability of various structural features reported here for homologous transient PPCs are higher than the variation in homologous permanent homomers. Our findings suggest that straightforward extrapolation of interfacial nature and inter-residue interaction patterns from template to target could lead to serious errors in the modeled complex structure. Understanding the evolution of interfaces provides insights to improve comparative modeling of PPC structures.

Characterization of the six glycosyltransferases involved in the biosynthesis of Yersinia enterocolitica serotype O:3 lipopolysaccharide outer core

Yersinia enterocolitica (Ye) is a gram-negative bacterium; Ye serotype O:3 expresses lipopolysaccharide (LPS) with a hexasaccharide branch known as the outer core (OC). The OC is important for the resistance of the bacterium to cationic antimicrobial peptides and also functions as a receptor for bacteriophage phiR1-37 and enterocoliticin. The biosynthesis of the OC hexasaccharide is directed by the OC gene cluster that contains nine genes (wzx, wbcKLMNOPQ, and gne). In this study, we inactivated the six OC genes predicted to encode glycosyltransferases (GTase) one by one by nonpolar mutations to assign functions to their gene products. The mutants expressed no OC or truncated OC oligosaccharides of different lengths. The truncated OC oligosaccharides revealed that the minimum structural requirements for the interactions of OC with bacteriophage phiR1-37, enterocoliticin, and OC-specific monoclonal antibody 2B5 were different. Furthermore, using chemical and structural analyses of the mutant LPSs, we could assign specific functions to all six GTases and also revealed the exact order in which the transferases build the hexasaccharide. Comparative modeling of the catalytic sites of glucosyltransferases WbcK and WbcL followed by site-directed mutagenesis allowed us to identify Asp-182 and Glu-181, respectively, as catalytic base residues of these two GTases. In general, conclusive evidence for specific GTase functions have been rare due to difficulties in accessibility of the appropriate donors and acceptors; however, in this work we were able to utilize the structural analysis of LPS to get direct experimental evidence for five different GTase specificities.

Determinação de motivos de ligação à quitina em vicilinas de Canavalia ensiformis e Vigna unguiculata através de métodos in silico e relação com suas toxicidades para o bruquídeo Callosobruchus maculatus (Coleoptera:Bruchidae)

Chitin is an important structural component of the cellular wall of fungi and exoskeleton of many invertebrate plagues, such as insects and nematodes. In digestory systems of insects it forms a named matrix of peritrophic membrane. One of the most studied interaction models protein-carbohydrate is the model that involves chitin-binding proteins. Among the involved characterized domains already in this interaction if they detach the hevein domain (HD), from of Hevea brasiliensis (Rubber tree), the R&R consensus domain (R&R), found in cuticular proteins of insects, and the motif called in this study as conglicinin motif (CD), found in the cristallography structure of the β-conglicinin bounded with GlcNac. These three chitin-binding domains had been used to determine which of them could be involved in silico in the interaction of Canavalia ensiformis and Vigna unguiculata vicilins with chitin, as well as associate these results with the WD50 of these vicilins for Callosobruchus maculatus larvae. The technique of comparative modeling was used for construction of the model 3D of the vicilin of V. unguiculata, that was not found in the data bases. Using the ClustalW program it was gotten localization of these domains in the vicilins primary structure. The domains R&R and CD had been found with bigger homology in the vicilins primary sequences and had been target of interaction studies. Through program GRAMM models of interaction ( dockings ) of the vicilins with GlcNac had been gotten. The results had shown that, through analysis in silico, HD is not part of the vicilins structures, proving the result gotten with the alignment of the primary sequences; the R&R domain, although not to have structural similarity in the vicilins, probably it has a participation in the activity of interaction of these with GlcNac; whereas the CD domain participates directly in the interaction of the vicilins with GlcNac. These results in silico show that the amino acid number, the types and the amount of binding made for the CD motif with GlcNac seem to be directly associates to the deleterious power that these vicilins show for C. maculatus larvae. This can give an initial step in the briefing of as the vicilins interact with alive chitin in and exert its toxic power for insects that possess peritrophic membrane

Modelagem molecular da Chiquimato quinase de Mycobacterium leprae

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Bioinformática estrutural aplicada ao estudo de proteínas alvo do genoma do Mycobacterium tuberculosis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Function inferences from a molecular structural model of bacterial ParE toxin

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Purification, characterization and structural determination of chitinases produced by Moniliophthora perniciosa

The enzyme chitinase from Moniliophthora perniciosa the causative agent of the witches' broom disease in Theobroma cacao, was partially purified with ammonium sulfate and filtration by Sephacryl S-200 using sodium phosphate as an extraction buffer. Response surface methodology (RSM) was used to determine the optimum pH and temperature conditions. Four different isoenzymes were obtained: ChitMp I, ChitMp II, ChitMp III and ChitMp IV. ChitMp I had an optimum temperature at 44-73ºC and an optimum pH at 7.0-8.4. ChitMp II had an optimum temperature at 45-73ºC and an optimum pH at 7.0-8.4. ChitMp III had an optimum temperature at 54-67ºC and an optimum pH at 7.3-8.8. ChitMp IV had an optimum temperature at 60ºC and an optimum pH at 7.0. For the computational biology, the primary sequence was determined in silico from the database of the Genome/Proteome Project of M. perniciosa, yielding a sequence with 564 bp and 188 amino acids that was used for the three-dimensional design in a comparative modeling methodology. The generated models were submitted to validation using Procheck 3.0 and ANOLEA. The model proposed for the chitinase was subjected to a dynamic analysis over a 1 ns interval, resulting in a model with 91.7% of the residues occupying favorable places on the Ramachandran plot and an RMS of 2.68.

Structural model of the CopA copper ATPase of Enterococcus hirae based on chemical cross-linking

The CopA copper ATPase of Enterococcus hirae belongs to the family of heavy metal pumping CPx-type ATPases and shares 43% sequence similarity with the human Menkes and Wilson copper ATPases. Due to a lack of suitable protein crystals, only partial three-dimensional structures have so far been obtained for this family of ion pumps. We present a structural model of CopA derived by combining topological information obtained by intramolecular cross-linking with molecular modeling. Purified CopA was cross-linked with different bivalent reagents, followed by tryptic digestion and identification of cross-linked peptides by mass spectrometry. The structural proximity of tryptic fragments provided information about the structural arrangement of the hydrophilic protein domains, which was integrated into a three-dimensional model of CopA. Comparative modeling of CopA was guided by the sequence similarity to the calcium ATPase of the sarcoplasmic reticulum, Serca1, for which detailed structures are available. In addition, known partial structures of CPx-ATPase homologous to CopA were used as modeling templates. A docking approach was used to predict the orientation of the heavy metal binding domain of CopA relative to the core structure, which was verified by distance constraints derived from cross-links. The overall structural model of CopA resembles the Serca1 structure, but reveals distinctive features of CPx-type ATPases. A prominent feature is the positioning of the heavy metal binding domain. It features an orientation of the Cu binding ligands which is appropriate for the interaction with Cu-loaded metallochaperones in solution. Moreover, a novel model of the architecture of the intramembranous Cu binding sites could be derived.

Determinação de motivos de ligação à quitina em vicilinas de Canavalia ensiformis e Vigna unguiculata através de métodos in silico e relação com suas toxicidades para o bruquídeo Callosobruchus maculatus (Coleoptera:Bruchidae)

Chitin is an important structural component of the cellular wall of fungi and exoskeleton of many invertebrate plagues, such as insects and nematodes. In digestory systems of insects it forms a named matrix of peritrophic membrane. One of the most studied interaction models protein-carbohydrate is the model that involves chitin-binding proteins. Among the involved characterized domains already in this interaction if they detach the hevein domain (HD), from of Hevea brasiliensis (Rubber tree), the R&R consensus domain (R&R), found in cuticular proteins of insects, and the motif called in this study as conglicinin motif (CD), found in the cristallography structure of the β-conglicinin bounded with GlcNac. These three chitin-binding domains had been used to determine which of them could be involved in silico in the interaction of Canavalia ensiformis and Vigna unguiculata vicilins with chitin, as well as associate these results with the WD50 of these vicilins for Callosobruchus maculatus larvae. The technique of comparative modeling was used for construction of the model 3D of the vicilin of V. unguiculata, that was not found in the data bases. Using the ClustalW program it was gotten localization of these domains in the vicilins primary structure. The domains R&R and CD had been found with bigger homology in the vicilins primary sequences and had been target of interaction studies. Through program GRAMM models of interaction ( dockings ) of the vicilins with GlcNac had been gotten. The results had shown that, through analysis in silico, HD is not part of the vicilins structures, proving the result gotten with the alignment of the primary sequences; the R&R domain, although not to have structural similarity in the vicilins, probably it has a participation in the activity of interaction of these with GlcNac; whereas the CD domain participates directly in the interaction of the vicilins with GlcNac. These results in silico show that the amino acid number, the types and the amount of binding made for the CD motif with GlcNac seem to be directly associates to the deleterious power that these vicilins show for C. maculatus larvae. This can give an initial step in the briefing of as the vicilins interact with alive chitin in and exert its toxic power for insects that possess peritrophic membrane

Modeling and reasoning with qualitative comparative clinical knowledge

The number of clinical trials reports is increasing rapidly due to a large number of clinical trials being conducted; it, therefore, raises an urgent need to utilize the clinical knowledge contained in the clinical trials reports. In this paper, we focus on the qualitative knowledge instead of quantitative knowledge. More precisely, we aim to model and reason with the qualitative comparison (QC for short) relations which consider qualitatively how strongly one drug/therapy is preferred to another in a clinical point of view. To this end, first, we formalize the QC relations, introduce the notions of QC language, QC base, and QC profile; second, we propose a set of induction rules for the QC relations and provide grading interpretations for the QC bases and show how to determine whether a QC base is consistent. Furthermore, when a QC base is inconsistent, we analyze how to measure inconsistencies among QC bases, and we propose different approaches to merging multiple QC bases. Finally, a case study on lowering intraocular pressure is conducted to illustrate our approaches.