An automated method for modeling proteins on known templates using distance geometry.

We present an automated method incorporated into a software package, FOLDER, to fold a protein sequence on a given three-dimensional (3D) template. Starting with the sequence alignment of a family of homologous proteins, tertiary structures are modeled using the known 3D structure of one member of the family as a template. Homologous interatomic distances from the template are used as constraints. For nonhomologous regions in the model protein, the lower and the upper bounds for the interatomic distances are imposed by steric constraints and the globular dimensions of the template, respectively. Distance geometry is used to embed an ensemble of structures consistent with these distance bounds. Structures are selected from this ensemble based on minimal distance error criteria, after a penalty function optimization step. These structures are then refined using energy optimization methods. The method is tested by simulating the alpha-chain of horse hemoglobin using the alpha-chain of human hemoglobin as the template and by comparing the generated models with the crystal structure of the alpha-chain of horse hemoglobin. We also test the packing efficiency of this method by reconstructing the atomic positions of the interior side chains beyond C beta atoms of a protein domain from a known 3D structure. In both test cases, models retain the template constraints and any additionally imposed constraints while the packing of the interior residues is optimized with no short contacts or bond deformations. To demonstrate the use of this method in simulating structures of proteins with nonhomologous disulfides, we construct a model of murine interleukin (IL)-4 using the NMR structure of human IL-4 as the template. The resulting geometry of the nonhomologous disulfide in the model structure for murine IL-4 is consistent with standard disulfide geometry.     

Subcellular Distribution and Chemical Forms of Cadmium in the Medicine Food Homology Plant <i>Platycodon grandiflorum</i> (Jacq.) A.DC.

Although Platycodon grandiflorum (Jacq.) A.DC. is a renowned medicine food homology plant, reports of excessive cadmium (Cd) levels are common, which affects its safety for clinical use and food consumption. To enable its Cd levels to be regulated or reduced, it is necessary to first elucidate the mechanism of Cd uptake and accumulation in the plant, in addition to its detoxification mechanisms. This present study used inductively couple plasma-mass-spectrometry to analyze the subcellular distribution and chemical forms of Cd in different tissues of P. grandiflorum. The experimental results showed that Cd was mainly accumulated in the roots [predominantly in the cell wall (50.96%–61.42%)], and it was found primarily in hypomobile and hypotoxic forms. The proportion of Cd in the soluble fraction increased after Cd exposure, and the proportion of insoluble phosphate Cd and oxalate Cd increased in roots and leaves, with a higher increase in oxalate Cd. Therefore, it is likely that root retention mechanisms, cell wall deposition, vacuole sequestration, and the formation of low mobility and low toxicity forms are tolerance strategies for Cd detoxification used by P. grandiflorum. The results of this study provide a theoretical grounding for the study of Cd accumulation and detoxification mechanisms in P. grandiflorum, and they can be used as a reference for developing Cd limits and standards for other medicine food homology plants.     

Homology among RAPD fragments in interspecific comparisons

The use of RAPDs for comparative purposes relies on the assumption that similarity of fragment size is a dependable indicator of homology. To test the validity of this assumption, homology among 220 pairs of comigrating fragments from three wild sunflower species was determined. Ninety-one per cent cross-hybridized and/or displayed congruent restriction fragment profiles suggestive of homology. However, comparative linkage mapping data indicated that 13% of the homologous loci mapped to genomic locations that were incongruent with the majority of loci, suggestive of paralogous rather than orthologous relationships. Thus, of the 220 pairwise comparisons, only 174 (79.1%) identified loci that are useful for comparative genetic studies. These problems, as well as several other factors discussed in the text, will introduce noise into RAPD data sets and thereby reduce the probability of generating accurate estimates of genetic relationships. Recommended methods for reducing noise in RAPD data sets include increasing gel resolution and/or testing fragment homology. However, implementation of these approaches will not eliminate all uncertainties, and it is also recommended that RAPD data sets be tested for structure and reliability.     

Structure-function relationships of the alternative oxidase of plant mitochondria: A model of the active site

A major characteristic of plant mitochondria is the presence of a cyanide-insensitive alternative oxidase which catalyzes the reduction of oxygen to water. Current information on the properties of the oxidase is reviewed. Conserved amino acid motifs have been identified which suggest the presence of a hydroxo-bridged di-iron center in the active site of the alternative oxidase. On the basis of sequence comparison with other di-iron center proteins, a structural model for the active site of the alternative oxidase has been developed that has strong similarity to that of methane monoxygenase. Evidence is presented to suggest that the alternative oxidase of plant mitochondria is the newest member of the class II group of di-iron center proteins.     

Peptide recognition by PTB and PDZ domains

Protein tyrosine binding (PTB) and ‘post synaptic density disc-large zo-1’ (PDZ) domains bind to short peptidic ligands by augmentation of one of the domain's β sheets and other recognition mechanisms. The two domain classes have a superficial resemblance to each other, even though no sequential homology exists. The structural bases of the interactions are well understood for the domains now experimentally determined, and ligand—target pairs can probably be identified in favorable cases by analogy with the known domains. For both PTB and PDZ classes, functional activities are still not fully defined: it is possible that these domain classes, along with pleckstrin homology domains, have multiple roles.     

Subtilases: the superfamily of subtilisin-like serine proteases.

Subtilases are members of the clan (or superfamily) of subtilisin-like serine proteases. Over 200 subtilases are presently known, more than 170 of which with their complete amino acid sequence. In this update of our previous overview (Siezen RJ, de Vos WM, Leunissen JAM, Dijkstra BW, 1991, Protein Eng 4:719-731), details of more than 100 new subtilases discovered in the past five years are summarized, and amino acid sequences of their catalytic domains are compared in a multiple sequence alignment. Based on sequence homology, a subdivision into six families is proposed. Highly conserved residues of the catalytic domain are identified, as are large or unusual deletions and insertions. Predictions have been updated for Ca(2+)-binding sites, disulfide bonds, and substrate specificity, based on both sequence alignment and three-dimensional homology modeling.     

Evidence for PDZ domains in bacteria, yeast, and plants.

Several dozen signaling proteins are now known to contain 80-100 residue repeats, called PDZ (or DHR or GLGF) domains, several of which interact with the C-terminal tetrapeptide motifs X-Ser/Thr-X-Val-COO- of ion channels and/or receptors. PDZ domains have previously been noted only in mammals, flies, and worms, suggesting that the primordial PDZ domain arose relatively late in eukaryotic evolution. Here, techniques of sequence analysis-including local alignment, profile, and motif database searches-indicate that PDZ domain homologues are present in yeast, plants, and bacteria. It is suggested that two PDZ domains occur in bacterial high-temperature requirement A (htrA) and one in tail-specific protease (tsp) homologues, and that a yeast htrA homologue contains four PDZ domains. Sequence comparisons suggest that the spread of PDZ domains in these diverse organisms may have occurred via horizontal gene transfer. The known affinity of Escherichia coli tsp for C-terminal polypeptides is proposed to be mediated by its PDZ-like domain, in a similar manner to the binding of C-terminal polypeptides by animal PDZ domains.     

Characterization of a Small Family (CAIII) of Microsatellite-Containing Sequences with X-Y Homology

Four X-linked loci showing homology with a previously described Y-linked polymorphic locus (DYS413) were identified and characterized. By fluorescent in situ hybridization (FISH), somatic cell hybrids, and YAC screening, the X-linked members of this small family of sequences (CAIII) all map in Xp22, while the Y members map in Yq11. These loci contribute to the overall similarity of the two genomic regions. All of the CAIII loci contain an internal microsatellite of the (CA)_n type. The microsatellites display extensive length polymorphism in two of the X-linked members as well as in the Y members. In addition, common sequence variants are found in the portions flanking the microsatellites in two of the X-linked members. Our results indicate that, during the evolution of this family, length variation on the Y chromosome was accumulated at a rate not slower than that on the X chromosome. Finally, these sequences represent a model system with which to analyze human populations for similar X- and Y-linked polymorphisms. Received: 29 July 1996 / Accepted: 15 January 1997     

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.