Crystal structure of a protein-toxin alpha 1-purothionin at 2.5A and a comparison with predicted models

Alpha 1-Purothionin (alpha 1-P), a wheatgerm protein and lytic toxin, has a secondary and tertiary structure similar to that of crambin as revealed by CD and NMR studies. alpha 1-P crystallizes in the tetragonal space group 1422 with unit cell dimensions: a = b = 53.59 and c = 69.79 A. X-ray diffraction data have been measured to 2.5 A Bragg spacing. The crystal structure has been determined by molecular replacement methods, using an energy-minimized alpha 1-P model structure derived from crambin (Whitlow and Teeter: Journal of Biomolecular Structure and Dynamics 2:831-848, 1985, Journal of the American Chemical Society 108:7163-7172, 1986). The energy-minimized model gives a slightly cleaner rotation solution and better refinement against the x-ray data than do the crambin or unminimized alpha 1-P structures. The final crystallographic residual with the data in the 10-2.5 A resolution range is 0.216. The refined alpha 1-P structure has a backbone rms difference of 0.74 A from crambin and 0.55 A from the energy-minimized alpha 1-P model. A low resolution NMR model of alpha 1-P calculated from metric matrix distance geometry and restrained molecular dynamics differs from crambin's backbone by 2.3 A rms deviation (Clore et al.: EMBO Journal 5:2729-2735, 1986). Backbone dihedral angles for our predicted model differ from the refined alpha 1-P structure in only one region (at a turn where there is a deletion relative to crambin). The NMR model had differences in four regions.     

Nuclear magnetic resonance study of the solution structure of alpha 1-purothionin. Sequential resonance assignment, secondary structure and low resolution tertiary structure

The solution structure of the 45-residue plant protein, alpha 1-purothionin, is investigated by nuclear magnetic resonance (n.m.r.) spectroscopy. Using a combination of two-dimensional n.m.r. techniques to demonstrate through-bond and through-space (less than 5 A) connectivities, the 1H n.m.r. spectrum of alpha 1-purothionin is assigned in a sequential manner. The secondary structure elements are then delineated on the basis of a qualitative interpretation of short-range nuclear Overhauser effects (NOE) involving the NH, C alpha H and C beta H protons. There are two helices extending from residues 10 to 19 and 23 to 28, two short beta-strands from residues 3 to 5 and 31 to 34 which form a mini anti-parallel beta-sheet, and five turns. In addition, a number of long-range NOE connectivities are assigned and a low resolution tertiary structure is proposed.     

Test of circular dichroism (CD) methods for crambin and CD-assisted secondary structure prediction of its homologous toxins

Methods that analyze protein circular dichroism (CD) spectra for fractions of secondary structure are evaluated for the plant protein crambin, which has a known high-resolution crystal structure. In addition, a two-step secondary structure prediction scheme is presented and used for the toxins homologous to crambin, shown by others to have secondary structures similar to crambin. The test of CD spectral analysis methods with the protein crambin employed two computer programs and several CD basis sets. Crambin's crystal structure, known to 0.945A resolution (Hendrickson, W.A., Teeter, M.M. Nature 290:107-113, 1981), allows accurate evaluation of results. Analysis with the protein spectra basis sets (Provencher, S.W., Gl?ckner, J. Biochemistry 20:33-37, 1981) as modified (Manavalan, P., Johnson, W.C., Jr. Anal. Biochem. 167:76-85, 1987) agreed most closely with crambin's crystal structure. This method was then applied to the CD spectra of the membrane-active toxins homologous to crambin (alpha 1- and beta-purothionin, phoratoxin A and B, and viscotoxin A3 and B). The new program SEQ (pronounced "seek") was developed to assign the secondary structure along the protein chain in a hierarchical fashion and applied to the plant toxins. The method constrained the secondary structure fractions to those from CD analysis and combined standard statistical methods with amphipathic helix location. Both CD-arrived secondary structure percentages and sequence assignment indicate that the viscotoxins are structurally most similar to crambin. Purothionin's secondary structure was predicted to be fundamentally similar to crambin's with a difference at the start of the first helix. This assignment agreed with Raman and NMR analyses of purothionin and lends validity to the method presented here. Differences from the NMR in the CD secondary structure fraction analysis for phoratoxin suggest interference in the CD from tryptophan residues.     

Tertiary structure prediction of mixed α/β proteins via energy minimization

We describe an improved algorithm for protein structure prediction, assuming that the location of secondary structural elements is known, with particular focus on prediction for proteins containing β-strands. Hydrogen bonding terms are incorporated into the potential function, supplementing our previously developed residue-residue potential which is based on a combination of database statistics and an excluded volume term. Two small mixed α/β proteins, 1-CTF and BPTI, are studied. In order to obtain native-like structures, it is necessary to allow the β-strands in BPTI to distort substantially from an ideal geometry, and an automated algorithm to carry this out efficiently is presented. Simulated annealing Monte Carlo methods, which contain a genetic algorithm component as well, are used to produce an ensemble of low-energy structures. For both proteins, a cluster of structures with low RMS deviation from the native structure is generated and the energetic ranking of this cluster is in the top 2 or 3 clusters obtained from simulations. These results are encouraging with regard to the possibility of constructing a robust procedure for tertiary folding which is applicable to β-strand containing proteins. Proteins 33:240–252, 1998. © 1998 Wiley-Liss, Inc.     

TMBpro: secondary structure, beta-contact and tertiary structure prediction of transmembrane beta-barrel proteins

MOTIVATION: Transmembrane beta-barrel (TMB) proteins are embedded in the outer membranes of mitochondria, Gram-negative bacteria and chloroplasts. These proteins perform critical functions, including active ion-transport and passive nutrient intake. Therefore, there is a need for accurate prediction of secondary and tertiary structure of TMB proteins. Traditional homology modeling methods, however, fail on most TMB proteins since very few non-homologous TMB structures have been determined. Yet, because TMB structures conform to specific construction rules that restrict the conformational space drastically, it should be possible for methods that do not depend on target-template homology to be applied successfully. RESULTS: We develop a suite (TMBpro) of specialized predictors for predicting secondary structure (TMBpro-SS), beta-contacts (TMBpro-CON) and tertiary structure (TMBpro-3D) of transmembrane beta-barrel proteins. We compare our results to the recent state-of-the-art predictors transFold and PRED-TMBB using their respective benchmark datasets, and leave-one-out cross-validation. Using the transFold dataset TMBpro predicts secondary structure with per-residue accuracy (Q(2)) of 77.8%, a correlation coefficient of 0.54, and TMBpro predicts beta-contacts with precision of 0.65 and recall of 0.67. Using the PRED-TMBB dataset, TMBpro predicts secondary structure with Q(2) of 88.3% and a correlation coefficient of 0.75. All of these performance results exceed previously published results by 4% or more. Working with the PRED-TMBB dataset, TMBpro predicts the tertiary structure of transmembrane segments with RMSD <6.0 A for 9 of 14 proteins. For 6 of 14 predictions, the RMSD is <5.0 A, with a GDT_TS score greater than 60.0. AVAILABILITY: http://www.igb.uci.edu/servers/psss.html.     

An all-atom force field for tertiary structure prediction of helical proteins

We have developed an all-atom free-energy force field (PFF01) for protein tertiary structure prediction. PFF01 is based on physical interactions and was parameterized using experimental structures of a family of proteins believed to span a wide variety of possible folds. It contains empirical, although sequence-independent terms for hydrogen bonding. Its solvent-accessible surface area solvent model was first fit to transfer energies of small peptides. The parameters of the solvent model were then further optimized to stabilize the native structure of a single protein, the autonomously folding villin headpiece, against competing low-energy decoys. Here we validate the force field for five nonhomologous helical proteins with 20-60 amino acids. For each protein, decoys with 2-3 A backbone root mean-square deviation and correct experimental Cbeta-Cbeta distance constraints emerge as those with the lowest energy.     

Secondary structure prediction and unrefined tertiary structure prediction for cyclin A,B, and D

We present heuristic-based predictions of the secondary and tertiary structures of the cyclins A, B, and D, representatives of the cyclin superfamily. The list of suggested constraints for tertiary structure assembly was left unrefined in order to submit this report before an announced crystal structure for cyclin A becomes available. To predict these constraints, a master sequence alignment over 270 positions of cyclin types A, B, and D was adjusted based on individual secondary structure predictions for each type. We used new heuristics for predicting aromatic residues at protein-protein interfaces and to identify sequentially distinct regions in the protein chain that cluster in the folded structure. The boundaries of two conjectured domains in the cyclin fold were predicted based on experimental data in the literature. The domain that is important for interaction of the cyclins with cyclin-dependent kinases (CDKs) is predicted to contain six helices; the second domain in the consensus model contains both helices and a β-sheet that is formed by sequentially distant regions in the protein chain. A plausible phosphorylation site is identified. This work represents a blinded test of the method for prediction of secondary and, to a lesser extent, tertiary structure from a set of homologous protein sequences. Evaluation of our predictions will become possible with the publication of the announced crystal structure.     

Molecular evolution of serpins: homologous structure of the human alpha 1-antichymotrypsin and alpha 1-antitrypsin genes

alpha 1-Antichymotrypsin belongs to a supergene family that includes alpha 1-antitrypsin, antithrombin III, ovalbumin, and angiotensinogen. The human chromosomal alpha 1-antichymotrypsin gene has been cloned and its molecular structure established. The gene is approximately 12 kb in length and contains five exons and four introns. The locations of the introns within the alpha 1-antichymotrypsin gene are identical with those of the human alpha 1-antitrypsin and angiotensinogen genes. Other members of this supergene family contain introns located at nonhomologous positions of the genes. The homologous organization of the alpha 1-antichymotrypsin and alpha 1-antitrypsin genes corresponds with the high degree of homology between their protein sequences and suggests that these loci arose by recent gene duplication. A model is presented for the evolution of both the genomic structure and the protein sequences of the serine protease inhibitor superfamily.     

P-cats: prediction of catalytic residues in proteins from their tertiary structures

SUMMARY: P-cats is a web server that predicts the catalytic residues in proteins from the atomic coordinates. P-cats receives a coordinate file of the tertiary structure and sends out analytical results via e-mail. The reply contains a summary and two URLs to allow the user to examine the conserved residues: one for interactive images of the prediction results and the other for a graphical view of the multiple sequence alignment. AVAILABILITY: P-cats is freely available at http://p-cats.hgc.jp/p-cats CONTACT: kino@ims.u-tokyo.ac.jp     

Conservation of polyproline II helices in homologous proteins: implications for structure prediction by model building.

Left-handed polyproline II (PPII) helices commonly occur in globular proteins in segments of 4-8 residues. This paper analyzes the structural conservation of PPII-helices in 3 protein families: serine proteinases, aspartic proteinases, and immunoglobulin constant domains. Calculations of the number of conserved segments based on structural alignment of homologous molecules yielded similar results for the PPII-helices, the alpha-helices, and the beta-strands. The PPII-helices are consistently conserved at the level of 100-80% in the proteins with sequence identity above 20% and RMS deviation of structure alignments below 3.0 A. The most structurally important PPII segments are conserved below this level of sequence identity. These results suggest that the PPII-helices, in addition to the other 2 secondary structure classes, should be identified as part of structurally conserved regions in proteins. This is supported by similar values for the local RMS deviations of the aligned segments for the structural classes of PPII-helices, alpha-helices, and beta-strands. The PPII-helices are shown to participate in supersecondary elements such as PPII-helix/alpha-helix. The conservation of PPII-helices depends on the conservation of a supersecondary element as a whole. PPII-helices also form links, possibly flexible, in the interdomain regions. The role of the PPII-helices in model building by homology is 2-fold; they serve as additional conserved elements in the structure allowing improvement of the accuracy of a model and provide correct chain geometry for modeling of the segments equivalenced to them in a target sequence. The improvement in model building is demonstrated in 2 test studies.     

Calmodulin binding to alpha 1-purothionin: solution binding and modeling of the complex.

CD and fluorescence spectroscopic measurements show that calmodulin (CaM) binds to purothionins (alpha 1-purothionin: alpha 1-PT; beta-purothionin: beta-PT) in 1:1 stoichiometry with an affinity similar to that exhibited with the tightest binding CaM-binding peptides. Using the available crystal structures of CaM and alpha 1-PT, a model has been built for the interaction of CaM and alpha 1-PT and subjected to potential energy minimization. In the model, there is a bend in the central helix of CaM similar to that suggested by Persechini and Kretsinger (J. Card. Pharm. 12:501-512, 1988). alpha 1-PT fits snugly into the cavity formed by the bent CaM molecule with each of its two helices making apolar interactions with each of the two hydrophobic clefts situated at the terminal domains of CaM. The complex is further stabilized by numerous polar and electrostatic interactions on the rims of the clefts. Our model is compared with two other similar models previously reported for the CaM complexes with other helical peptides and generalizations about the mode of CaM binding to target proteins are made, which have wide relevance to the function of CaM. By analogy, a similar model is predicted for a CaM-beta-PT complex.     

Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry

Globular proteins have high packing densities as a result of residue side chains in the core achieving a tight, complementary packing. The internal packing is considered the main determinant of native protein structure. From that point of view, we present here a method of energy minimization using an automata network to predict a set of amino acid sequences and their side-chain conformations from a desired backbone geometry for de novo design of proteins. Using discrete side-chain conformations, that is, rotamers, the sequence generation problem from a given backbone geometry becomes one of combinatorial problems. We focused on the residues composing the interior core region and predicted a set of amino acid Sequences and their side-chain conformations only from a given backbone geometry. The kinds of residues were restricted to six hydrophobic amino acids (Ala, Ile, Met, Leu, Phe, and Val) because the core regions are almost always composed of hydrophobic residues. The obtained sequences were well packed as was the native sequence. The method can be used for automated sequence generation in the de novo design of proteins. © 1994 Wiley-Liss, Inc.     

The three-dimensional structure of alpha1-purothionin in solution: combined use of nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The determination of the three-dimensional solution structure of α1-purothionin using a combination of metric matrix distance geometry and restrained molecular dynamics calculations based on n.m.r. data is presented. The experimental data comprise complete sequence-specific proton resonance assignments, a set of 310 approximate interproton distance restraints derived from nuclear Overhauser effects, 27 Ø backbone torsion angle restraints derived from vicinal coupling constants, 4 distance restraints from hydrogen bonds and 12 distance restraints from disulphide bridges. The average atomic rms difference between the final nine converged structures and the mean structure obtained by averaging their coordinates is 1.5 ± 0.1 å for the backbone atoms and 2.0 ± 0.1 å for all atoms. The overall shape of α1-purothionin is that of the capital letter L, similar to that of crambin, with the longer arm comprising two approximately parallel α-helices and the shorter arm a strand and a mini anti-parallel β sheet.     

Active site studies of bovine alpha1-->3-galactosyltransferase and its secondary structure prediction

The catalytic domain of bovine alpha1-->3-galactosyltransferase (alpha3GalT), residues 80-368, have been cloned and expressed, in Escherichia coli. Using a sequential purification protocol involving a Ni(2+) affinity column followed by a UDP-hexanolamine affinity column, we have obtained a pure and active protein from the soluble fraction which catalyzes the transfer of galactose (Gal) from UDP-Gal to N-acetyllactosamine (LacNAc) with a specific activity of 0.69 pmol/min/ng. The secondary structural content of alpha3GalT protein was analyzed by Fourier transform infrared (FTIR) spectroscopy, which shows that the enzyme has about 35% beta-sheet and 22% alpha-helix. This predicted secondary structure content by FTIR spectroscopy was used in the protein sequence analysis algorithm, developed by the Biomolecular Engineering Research Center at Boston University and Tasc Inc., for the assignment of secondary structural elements to the amino acid sequence of alpha3GalT. The enzyme appears to have three major and three minor helices and five sheet-like structures. The studies on the acceptor substrate specificity of the enzyme, alpha3GalT, show that in addition to LacNAc, which is the natural substrate, the enzyme accepts various other disaccharides as substrates such as lactose and Gal derivatives, beta-O-methylgalactose and beta-D-thiogalactopyranoside, albeit with lower specific activities. There is an absolute requirement for Gal to be at the non-reducing end of the acceptor molecule which has to be beta1-->4-linked to a second residue that can be more diverse in structure. The kinetic parameters for four acceptor molecules were determined. Lactose binds and functions in a similar way as LacNAc. However, beta-O-methylgalactose and Gal do not bind as tightly as LacNAc or lactose, as their K(ia) and K(A) values indicate, suggesting that the second monosaccharide is critical for holding the acceptor molecule in place. The 2' and 4' hydroxyl groups of the receiving Gal moiety are important in binding. Even though there is large structural variability associated with the second residue of the acceptor molecule, there are constraints which do not allow certain Gal-R sugars to be good acceptors for the enzyme. The beta1-->4-linked residue at the second position of the acceptor molecule is preferred, but the interactions between the enzyme and the second residue are likely to be non-specific.     

Study and prediction of secondary structure for membrane proteins

In this paper we present a novel approach to membrane protein secondary structure prediction based on the statistical stepwise discriminant analysis method. A new aspect of our approach is the possibility to derive physical-chemical properties that may affect the formation of membrane protein secondary structure. The certain physical-chemical properties of protein chains can be used to clarify the formation of the secondary structure types under consideration. Another aspect of our approach is that the results of multiple sequence alignment, or the other kinds of sequence alignment, are not used in the frame of the method. Using our approach, we predicted the formation of three main secondary structure types (alpha-helix, beta-structure and coil) with high accuracy, that is Q(3) = 76%. Predicting the formation of alpha-helix and non-alpha-helix states we reached the accuracy which was measured as Q(2) = 86%. Also we have identified certain protein chain properties that affect the formation of membrane protein secondary structure. These protein properties include hydrophobic properties of amino acid residues, presence of Gly, Ala and Val amino acids, and the location of protein chain end.     

Evidence for a novel circulating alpha 1 protein in tumour-bearing rats. Differentiation from acute phase alpha 1 proteins and from alpha 1 fetoprotein

It has been found in this study that the serum from rats bearing a transplanted dibenzanthracene-induced tumour ($\text{RD}\text{3}$), has a high concentration of alpha₁ proteins compared with normal rat serum. These alpha₁ proteins have been isolated by an immunoabsorption method and have been compared by immunological methods with the acute phase alpha₁ proteins isolated by the same method from the serum of rats presenting an inflammatory reaction. It has been found that the isolated $\text{RD}\text{3}$ alpha₁ proteins were composed of two major proteins: one of these corresponded to an inflammatory protein, the alpha₁-AP-globulin. The other may be a new protein, as it is absent from the serum of rats with an acute phase inflammatory reaction and nor does it correspond to alpha₁ feto-protein, a carcino-embryonic protein presenting the same electrophoretic mobility.     

MUFOLD: A new solution for protein 3D structure prediction

There have been steady improvements in protein structure prediction during the past 2 decades. However, current methods are still far from consistently predicting structural models accurately with computing power accessible to common users. Toward achieving more accurate and efficient structure prediction, we developed a number of novel methods and integrated them into a software package, MUFOLD. First, a systematic protocol was developed to identify useful templates and fragments from Protein Data Bank for a given target protein. Then, an efficient process was applied for iterative coarse‐grain model generation and evaluation at the Cα or backbone level. In this process, we construct models using interresidue spatial restraints derived from alignments by multidimensional scaling, evaluate and select models through clustering and static scoring functions, and iteratively improve the selected models by integrating spatial restraints and previous models. Finally, the full‐atom models were evaluated using molecular dynamics simulations based on structural changes under simulated heating. We have continuously improved the performance of MUFOLD by using a benchmark of 200 proteins from the Astral database, where no template with >25% sequence identity to any target protein is included. The average root‐mean‐square deviation of the best models from the native structures is 4.28 Å, which shows significant and systematic improvement over our previous methods. The computing time of MUFOLD is much shorter than many other tools, such as Rosetta. MUFOLD demonstrated some success in the 2008 community‐wide experiment for protein structure prediction CASP8. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Kynurenine formamidase: determination of primary structure and modeling-based prediction of tertiary structure and catalytic triad

Kynurenine formamidase (KFase) (EC 3.5.1.9) hydrolyzes N-formyl-L-kynurenine, an obligatory step in the conversion of tryptophan to nicotinic acid. Low KFase activity in chicken embryos, from inhibition by organophosphorus insecticides and their metabolites such as diazoxon, leads to marked developmental abnormalities. While KFase was purportedly isolated previously, the structure and residues important for catalysis and inhibition were not established. KFase was isolated here from mouse liver cytosol by (NH4)2SO4 precipitation and three FPLC steps (resulting in 221-fold increase in specific activity for N-formyl-L-kynurenine hydrolysis) followed by conversion to [3H]diethylphosphoryl-KFase and finally isolation by C4 reverse-phase high-performance liquid chromatography. Determination of tryptic fragment amino acid sequences and cDNA cloning produced a new 305-amino-acid protein sequence. Although an amidase by function, the primary structure of KFase lacks the amidase signature sequence and is more similar to esterases and lipases. Sequence profile analysis indicates KFase is related to the esterase/lipase/thioesterase family containing the conserved active-site serine sequence GXSXG. The alpha/beta-hydrolase fold is suggested for KFase by its primary sequence and predicted secondary conformation. A three-dimensional model based on the structures of homologous carboxylesterase EST2 and brefeldin A esterase implicates Ser162, Asp247 and His279 as the active site triad.     

NMR solution structure of viscotoxin C1 from Viscum album species Coloratum ohwi: toward a structure-function analysis of viscotoxins

The high resolution three-dimensional structure of the newly discovered plant viscotoxin C1, from the Asiatic Viscum album ssp. Coloratum ohwi, has been determined in solution by (1)H NMR spectroscopy at pH 3.6 and 285 K. The viscotoxin C1-fold, consisting of a helix-turn-helix motif and a short stretch of an antiparralel beta-sheet is very similar to that found for the highly similar viscotoxins A2 and A3 and for other related thionins. Different functional properties of members of the thionin family are discussed here in light of the structural and electrostatic properties. Among the very homologous family of alpha- and beta-thionins, known for their antimicrobial activity, the viscotoxin subfamily differs from the other members because of its high toxicity against tumoral cells. Key residues for the modulation of viscotoxin cytotoxicity have been identified on the basis of sequence and structural alignment.