The membrane topology of the amino-terminal domain of the red cell calcium pump.

A systematic study of the membrane-associated regions in the plasma membrane Ca2+ pump of erythrocytes has been performed by hydrophobic photolabeling. Purified Ca2+ pump was labeled with 3-(trifluoromethyl)-3-(m-[125I]iodophenyl)-diazirine ([125I]TID), a generic photoactivatable hydrophobic probe. These results were compared with the enzyme labeled with a strictly membrane-bound probe, [3H]bis-phosphatidylethanolamine (trifluoromethyl) phenyldiazirine. A significant light-dependent labeling of an M(r) 135,000-140,000 peptide, corresponding to the full Ca2+ pump, was observed with both probes. After proteolysis of the pump labeled with each probe and isolation of fragments by SDS-PAGE, a common pattern of labeled peptides was observed. Similarly, labeling of the Ca2+ pump with [125I]TID, either in isolated red blood cell membranes or after the enzyme was purified, yields a similar pattern of labeled peptides. Taken together, these results validate the use of either probe to study the lipid interface of the membrane-embedded region of this protein, and sustain the notion that the conformation of the pump is maintained throughout the procedures of solubilization, affinity purification, and reconstitution into proteoliposomes. In this work, we put special emphasis on a detailed analysis of the N-terminal domain of the Ca2+ pump. A labeled peptide of M(r) 40,000 belonging to this region was purified and further digested with V8 protease. The specific incorporation of [125I]TID to proteolytic fragments pertaining to the amino-terminal region indicates the existence of two transmembrane stretches in this domain. A theoretical analysis based on the amino acid sequence 1-322 predicts two segments with high probability of membrane insertion, in agreement with the experimental data. Each segment shows a periodicity pattern of hydrophobicity and variability compatible with alpha-helical structure. These results strongly suggest the existence of a transmembrane helical hairpin motif near the N-terminus of the Ca2+ pump.     

Three-dimensional structure prediction of the NAD binding site of proton-pumping transhydrogenase from Escherichia coli

A three-dimensional structure of the NAD site of Escerichia coli transhydrogenase has been predicted. The model is based on analysis of conserved residues among the transhydrogenases from five different sources, homologies with enzymes using NAD as cofactors or substrates, hydrophilicity profiles, and secondary structure predictions. The present model supports the hypothesis that there is one binding site, located relatively close to the N-terminus of the α-subunit. The proposed structure spans residues α145 to α287, and it includes five β-strands and five α-helices oriented in a typical open twisted α/β conformation. The amino acid sequence following the GXGXXG dinucleotide binding consensus sequence (residues α172 to α177) correlates exactly to a typical fingerprint region for ADP binding βαβ folds in dinucleotide binding enzymes. In the model, aspartic acid α195 forms hydrogen bonds to one or both hydroxyl groups on the adenosine ribose sugar moiety. Threonine α196 and alanine α256, located at the end of βB and βD, respectively, create a hydrophobic sandwich with the adenine part of NAD buried inside. The nicotinamide part is located in a hydrophobic cleft between αA and βE. Mutagenesis work has been carried out in order to test the predicted model and to determine whether residues within this domain are important for proton pumping directly. All data support the predicted structure, and no residue crucial for proton pumping Was detected. Since no three-dimensional structure of transhydrogenase has been solved, a well based tertiary structure prediction is of great value for further experimental design in trying to elucidate the mechanism of the energy-linked proton pump. © 1995 Wiley-Liss, Inc.     

Knowledge-based modeling of the D-lactate dehydrogenase three-dimensional structure

A three-dimensional structure of the NAD-dependent D -lactate dehydrogenase of Lactobacillus bulgaricus is modeled using the structure of the formate dehydrogenase of Pseudomonas sp. as template. Both sequences share only 22% of identical residues. Regions for knowledge-based modeling are defined from the structurally conserved regions predicted by multiple alignment of a set of related protein sequences with low homology. The model of the D -LDH subunit shows, as for the formate dehydrogenase, an α/β structure, with a catalytic domain and a coenzyme binding domain. It points out the catalytic histidine (His-296) and supports the hypothetical catalytic mechanism. It also suggests that the other residues involved in the active site are Arg-235, possibly involved in the binding of the carboxyl group of the pyruvate, and Phe-299, a candidate for stabilizing the methyl group of the substrate. © 1995 Wiley-Liss, Inc.     

A predicted consensus structure for the protein kinase C2 homology (C2H) domain,the repeating unit of synaptotagmin

A secondary structure has been predicted for the protein kinase C2 regulatory domain found in homologous form in synaptotagmin, some phospholipases, and some GTP activated proteins. The proposed structure is built from seven consecutive beta strands followed by a terminal alpha helix. Considerations of overall surface exposure of individual secondary structural elements suggest that these are packed into a 2-sheet beta sandwich structure, with one of only three of the many possible folds being preferred. © 1995 Wiley-Liss, Inc.     

A Continuum theory for the prediction of lateral and rotational positioning of α-helices in membrane proteins: Bacteriorhodopsin

We have developed a new method for the prediction of the lateral and the rotational positioning of transmembrane helices, based upon the present status of knowledge about the dominant interaction of the tertiary structure formation. The basic assumption about the interaction is that the interhelix binding is due to the polar interactions and that very short extramembrane loop segments restrict the relative position of the helices. Another assumption is made for the simplification of the prediction that a helix may be regarded as a continuum rod having polar interaction fields around it. The polar interaction field is calculated by a probe helix method, using a copolymer of serine and alanine as probe helices. The lateral position of helices is determined by the strength of the interhelix binding estimated from the polar interaction field together with the length of linking loop segments. The rotational positioning is determined by the polar interaction field, assuming the optimum lateral configuration. The structural change due to the binding of a prosthetic group is calculated, fixing the rotational freedom of a helix that is connected to the prosthetic group. Applying this method to bacteriorhodopsin, the optimum lateral and rotational positioning of transmembrane helices that are very similar to the experimental configuration was obtained. This method was implemented by a software system, which was developed for this work, and automatic calculation became possible for membrane proteins comprised of several transmembrane helices. © 1995 Wiley-Liss, Inc.     

Homology modeling of histidine-containing phosphocarrier protein and eosinophil-derived neurotoxin: Construction of models and comparison with experiment

Homology modeling methods have been used to construct models of two proteins—the histidine-containing phosphocarrier protein (HPr) from Mycoplasma capricolum and human eosinophil-derived neurotoxin (EDN). Comparison of the models with the subsequently determined X-ray crystal structures indicates that the core regions of both proteins are reasonably well reproduced, although the template structures are closer to the X-ray structures in these regions—possible enhancements are discussed. The conformations of most of the side chains in the core of HPr are well reproduced in the modeled structure. As expected, the conformations of surface side chains in this protein differ significantly from the X-ray structure. The loop regions of EDN were incorrectly modeled—reasons for this and possible enhancements are discussed. © 1995 Wiley-Liss, Inc.     

The cytidylyltransferase superfamily: Identification of the nucleotide-binding site and fold prediction

The crystal structure of glycerol-3-phosphate cytidylyltransferase from B. subtilis (TagD) is about to be solved. Here, we report a testable structure prediction based on the identification by sequence analysis of a superfamily of functionally diverse but structurally similar nucleotide-binding enzymes. We predict that TagD is a member of this family. The most conserved region in this superfamily resembles the ATP-binding HiGH motif of class I aminoacyI-tRNA synthetases. The predicted secondary structure of cytidylyltransferase and its homologues is compatible with the α/β topography of the class I aminoacyl-tRNA synthetases. The hypothesis of similarity of fold is strengthened by sequence-structure alignment and 3D model building using the known structure of tyrosyl tRNA synthetase as template. The proposed 3D model of TagD is plausible both structurally, with a well packed hydrophobic core, and functionally, as the most conserved residues cluster around the putative nucleotide binding site. If correct, the model would imply a very ancient evolutionary link between class I tRNA synthetases and the novel cytidylyltransferase superfamily. © 1995 Wiley-Liss, Inc.     

Progress of 1D protein structure prediction at last

Accuracy of predicting protein secondary structure and solvent accessibility from sequence information has been improved significantly by using information contained in multiple sequence alignments as input to a neural 'network system. For the Asilomar meeting, predictions for 13 proteins were generated automatically using the publicly available prediction method PHD. The results confirm the estimate of 72% three-state prediction accuracy. The fairly accurate predictions of secondary structure segments made the tool useful as a starting point for modeling of higher dimensional aspects of protein structure. © 1995 Wiley-Liss, Inc.     

Progress in fold recognition

The prediction experiment reveals that fold recognition has become a powerful tool in structural biology. We applied our fold recognition technique to 13 target sequences. In two cases, replication terminating protein and prosequence of subtilisin, the predicted structures are very similar to the experimentally determined folds. For the first time, in a public blind test, the unknown structures of proteins have been predicted ahead of experiment to an accuracy approaching molecular detail. In two other cases the approximate folds have been predicted correctly. According to the assessors there were 12 recognizable folds among the target proteins. In our postprediction analysis we find that in 7 cases our fold recognition technique is successful. In several of the remaining cases the predicted folds have interesting features in common with the experimental results. We present our procedure, discuss the results, and comment on several fundamental and technical problems encountered in fold recognition. © 1995 Wiley-Liss, Inc.     

Homology modeling by the ICM method

Five models have been built by the ICM method for the Comparative Modeling section of the Meeting on the Critical Assessment of Techniques for Protein Structure Prediction. The targets have homologous proteins with known three-dimensional structure with sequence identity ranging from 25 to 77%. After alignment of the target sequence with the related three-dimensional structure, the modeling procedure consists of two subproblems: side-chain prediction and loop prediction. The ICM method approaches these problems with the following steps: (1) a starting model is created based on the homologous structure with the conserved portion fixed and the noncon-served portion having standard covalent geometry and free torsion angles; (2) the Biased Probability Monte Carlo (BPMC) procedure is applied to search the subspaces of either all the nonconservative side-chain torsion angles or torsion angles in a loop backbone and surrounding side chains. A special algorithm was designed to generate low-energy loop deformations. The BPMC procedure globally optimizes the energy function consisting of ECEPP/3 and solvation energy terms. Comparison of the predictions with the NMR or crystallographic solutions reveals a high proportion of correctly predicted side chains. The loops were not correctly predicted because imprinted distortions of the backbone increased the energy of the near-native conformation and thus made the solution unrecognizable. Interestingly, the energy terms were found to be reliable and the sampling of conformational space sufficient. The implications of this finding for the strategies of future comparative modeling are discussed. © 1995 Wiley-Liss, Inc.