Shoshana Wodak | Search PubMed

PubMed2004

Progressive combinatorial algorithm for multiple structural alignments: application to distantly related proteins.

MALECON is a progressive combinatorial procedure for multiple alignments of protein structures. It searches a library of pairwise alignments for all three-protein alignments in which a specified number of residues is consistently aligned. These alignments are progressively expanded to include additional proteins and more spatially equivalent residues, subject to certain criteria. This action involves superimposing the aligned proteins by their hitherto equivalent residues and searching for additional Calpha atoms that lie close in space. The performance of MALECON is illustrated and compared with several extant multiple structure alignment methods by using as test the globin homologous superfamily, the OB and the Jellyrolls folds. MALECON gives better definitions of the common structural features in the structurally more diverse proteins of the OB and Jellyrolls folds, but it yields comparable results for the more similar globins. When no consistent multiple alignments can be derived for all members of a protein group, our procedure is still capable of automatically generating consistent alignments and common core definitions for subgroups of the members. This finding is illustrated for proteins of the OB fold and SH3 domains, believed to share common structural features, and should be very instrumental in homology modeling and investigations of protein evolution.

Algorithms↗

PubMed2002

Modeling the three-dimensional structure of H+-ATPase of Neurospora crassa.

Homology modeling in combination with transmembrane topology predictions are used to build the atomic model of Neurospora crassa plasma membrane H+-ATPase, using as template the 2.6 A crystal structure of rabbit sarcoplasmic reticulum Ca2+-ATPase [Toyoshima, C., Nakasako, M., Nomura, H. & Ogawa, H. (2000) Nature 405, 647-655]. Comparison of the two calcium-binding sites in the crystal structure of Ca2+-ATPase with the equivalent region in the H+-ATPase model shows that the latter is devoid of most of the negatively charged groups required to bind the cations, suggesting a different role for this region. Using the built model, a pathway for proton transport is then proposed from computed locations of internal polar cavities, large enough to contain at least one water molecule. As a control, the same approach is applied to the high-resolution crystal structure of halorhodopsin and the proton pump bacteriorhodopsin. This revealed a striking correspondence between the positions of internal polar cavities, those of crystallographic water molecules and, in the case of bacteriorhodopsin, the residues mediating proton translocation. In our H+-ATPase model, most of these cavities are in contact with residues previously shown to affect coupling of proton translocation to ATP hydrolysis. A string of six polar cavities identified in the cytoplasmic domain, the most accurate part of the model, suggests a proton entry path starting close to the phosphorylation site. Strikingly, members of the haloacid dehalogenase superfamily, which are close structural homologs of this domain but do not share the same function, display only one polar cavity in the vicinity of the conserved catalytic Asp residue.

Amino Acid Sequence↗