HydroPaCe: understanding and predicting cross-inhibition in serine proteases through hydrophobic patch centroids.

MOTIVATION: Protein-protein interfaces contain important information about molecular recognition. The discovery of conserved patterns is essential for understanding how substrates and inhibitors are bound and for predicting molecular binding. When an inhibitor binds to different enzymes (e.g. dissimilar sequences, structures or mechanisms what we call cross-inhibition), identification of invariants is a difficult task for which traditional methods may fail. RESULTS: To clarify how cross-inhibition happens, we model the problem, propose and evaluate a methodology called HydroPaCe to detect conserved patterns. Interfaces are modeled as graphs of atomic apolar interactions and hydrophobic patches are computed and summarized by centroids (HP-centroids), and their conservation is detected. Despite sequence and structure dissimilarity, our method achieves an appropriate level of abstraction to obtain invariant properties in cross-inhibition. We show examples in which HP-centroids successfully predicted enzymes that could be inhibited by the studied inhibitors according to BRENDA database. AVAILABILITY: www.dcc.ufmg.br/~raquelcm/hydropace CONTACT: valdetemg@ufmg.br; raquelcm@dcc.ufmg.br; santoro@icb.ufmg.br SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

On the characterization of energy networks of proteins.

The construction of a realistic theoretical model of proteins is determinant for improving the computational simulations of their structural and functional aspects. Modeling proteins as a network of non-covalent connections between the atoms of amino acid residues has shown valuable insights into these macromolecules. The energy-related properties of protein structures are known to be very important in molecular dynamics. However, these same properties have been neglected when the protein structures are modeled as networks of atoms and amino acid residues. A new approach for the construction of protein models based on a network of atoms is presented. This method, based on interatomic interaction, takes into account the energy and geometric aspects of the protein structures that were not employed before, such as atomic occlusion inside the protein, the use of solvation, protein modeling and analysis, and the use of energy potentials to estimate the energies of interatomic non-covalent contacts. As a result, we achieved a more realistic network model of proteins. This model has the virtue of being more robust in face of different unknown variables that usually are arbitrarily estimated. We were able to determine the most connected residues of all the proteins studied, so that we are now in a better condition to study their structural role.

Finding protein-protein interaction patterns by contact map matching.

We propose a novel method for defining patterns of contacts present in protein-protein complexes. A new use of the traditional contact maps (more frequently used for representation of the intra-chain contacts) is presented for analysis of inter-chain contacts. Using an algorithm based on image processing techniques, we can compare protein-protein interaction maps and also obtain a dissimilarity score between them. The same algorithm used to compare the maps can align the contacts of all the complexes and be helpful in the determination of a pattern of conserved interactions at the interfaces. We present an example for the application of this method by analyzing the pattern of interaction of bovine pancreatic trypsin inhibitors and trypsins, chymotrypsins, a thrombin, a matriptase, and a kallikrein - all classified as serine proteases. We found 20 contacts conserved in trypsins and chymotrypsins and 3 specific ones are present in all the serine protease complexes studied. The method was able to identify important contacts for the protein family studied and the results are in agreement with the literature.

Random amplified polymorphic DNA profiles of Trypanosoma cruzi isolates from chagasic patients with different clinical forms.

The genetic variability of 61 Trypanosoma cruzi isolates from 47 chronic chagasic patients of Minas Gerais state was analyzed by random amplified polymorphic DNA polymerase chain reaction (RAPD-PCR) using M13-40, lambdagt11-F, and L15996 primers. Cluster analysis by unweighted pair group method analysis was applied to RAPD profiles, and cluster analysis used to verify a possible correlation among different clinical forms of the disease from these patients. The T. cruzi isolates showed distinct grouping on tree topology, with the isolates not being possible to establish a correlation to the clinical forms of Chagas' disease. These data showed that the T. cruzi isolates from these patients would compose a group of populations well correlated genetically.

A contact map matching approach to protein structure similarity analysis

We modeled the problem of identifying how close two proteins are structurally by measuring the dissimilarity of their contact maps. These contact maps are colored images, in which the chromatic information encodes the chemical nature of the contacts. We studied two conceptually distinct image-processing algorithms to measure the dissimilarity between these contact maps; one was a content-based image retrieval method, and the other was based on image registration. In experiments with contact maps constructed from the protein data bank, our approach was able to identify, with greater than 80% precision, instances of monomers of apolipoproteins, globins, plastocyanins, retinol binding proteins and thioredoxins, among the monomers of Protein Data Bank Select. The image registration approach was only slightly more accurate than the content-based image retrieval approach.

Purification of a cysteine proteinase from Carica candamarcensis L. and cloning of a genomic putative fragment coding for this enzyme.

We describe the purification of a cysteine proteinase from latex of Carica candamarcensis, hereby designated CC23. The enzyme has been purified by ion-exchange chromatography and behaves electrophoretically as a monomer of M(r) 23,000 and optimal pH of 8.0. It displays a basic isoelectric point, has one cysteine residue in the active site by titration with E-64, confirmed by DNA sequencing, and responds to proteinase inhibitors as a classic cysteine proteinase. The K(m) and k(cat)/K(m) for CC23 using BAPNA were respectively 14.7 +/- 1.8 x 10(-4) M and 1.3 x 10(3) M(-1) s(-1). Therefore, the catalytic efficiency of CC23 is sixfold higher than that of CC-I, another proteinase from the same plant. DNA primers were designed to amplify by PCR a genomic sequence related to this enzyme. An 895-bp DNA fragment was cloned and sequenced. It shows strong homology with chymopapain isoform IV from C. papaya. The translated sequence is similar to that of chymopapain isoform II (73%) and CC-III (77%) from C. candamarcensis.

Evaluation of cDNA libraries from different developmental stages of Schistosoma mansoni for production of expressed sequence tags (ESTs).

A comparative study of the gene expression profile in different developmental stages of Schistosoma mansoni has been initiated based on the expressed sequence tag (EST) approach. A total of 1401 ESTs were generated from seven different cDNA libraries constructed from four distinct stages of the parasite life cycle. The libraries were first evaluated for their quality for a large-scale cDNA sequencing program. Most of them were shown to have less than 20% useless clones and more than 50% new genes. The redundancy of each library was also analyzed, showing that one adult worm cDNA library was composed of a small number of highly frequent genes. When comparing ESTs from distinct libraries, we could detect that most genes were present only in a single library, but others were expressed in more than one developmental stage and may represent housekeeping genes in the parasite. When considering only once the genes present in more than one library, a total of 466 unique genes were obtained, corresponding to 427 new S. mansoni genes. From the total of unique genes, 20.2% were identified based on homology with genes from other organisms, 8.3% matched S. mansoni characterized genes and 71.5% represent unknown genes.

Update of the gene discovery program in Schistosoma mansoni with the expressed sequence tag approach.

Continuing the Schistosoma mansoni Genome Project 363 new templates were sequenced generating 205 more ESTs corresponding to 91 genes. Seventy four of these genes (81%) had not previously been described in S. mansoni. Among the newly discovered genes there are several of significant biological interest such as synaptophysin, NIFs-like and rho-GDP dissociation inhibitor.