CORA is a suite of programs for multiply aligning and analyzing protein structural families to identify the consensus positions and capture their most conserved structural characteristics (e.g., residue accessibility, torsional angles, and global geometry as described by inter-residue vectors/contacts). Knowledge of these structurally conserved positions, which are mostly in the core of the fold and of their properties, significantly improves the identification and classification of newly-determined relatives. Information is encoded in a consensus three-dimensional (3D) template and relatives found by a sensitive alignment method, which employs a new scoring scheme based on conserved residue contacts. By encapsulating these critical “core” features, templates perform more reliably in recognizing distant structural relatives than searches with representative structures.

Parameters for 3D-template generation and alignment were optimized for each structural class (mainly-α, mainly-β, α-β), using representative superfold families. For all families selected, the templates gave significant improvements in sensitivity and selectivity in recognizing distant structural relatives. Furthermore, since templates contain less than 70% of fold positions and compare fewer positions when aligning structures, scans are at least an order of magnitude faster than scans using selected structures. CORA was subsequently tested on eight other broad structural families from the CATH database.

Diagnostics plots are generated automatically and provide qualitative assistance for classifying newly determined relatives. They are demonstrated here by application to the large globin-like fold family. CORA templates for both homologous superfamilies and fold families will be stored in CATH and used to improve the classification and analysis of newly determined structures.

We describe an extensive test of Geocore, an ab initio peptide folding algorithm. We studied 18 short molecules for which there are structures in the Protein Data Bank; chains are up to 31 monomers long. Except for the very shortest peptides, an extremely simple energy function is sufficient to discriminate the true native state from more than 108 lowest energy conformations that are searched explicitly for each peptide. A high incidence of native-like structures is found within the best few hundred conformations generated by Geocore for each amino acid sequence. Predictions improve when the number of discrete φ/ψ choices is increased.

The reoccurrence of water molecules in crystal structures of RNase T1 was investigated. Five waters were found to be invariant in RNase T1 as well as in six other related fungal RNases. The structural, dynamical, and functional characteristics of one of these conserved hydration sites (WAT1) were analyzed by protein engineering, X-ray crystallography, and 17O and 2H nuclear magnetic relaxation dispersion (NMRD). The position of WAT1 and its surrounding hydrogen bond network are unaffected by deletions of two neighboring side chains. In the mutant Thr93Gln, the Gln93Nε2 nitrogen replaces WAT1 and participates in a similar hydrogen bond network involving Cys6, Asn9, Asp76, and Thr91. The ability of WAT1 to form four hydrogen bonds may explain why evolution has preserved a water molecule, rather than a side-chain atom, at the center of this intricate hydrogen bond network. Comparison of the 17O NMRD profiles from wild-type and Thr93Gln RNase T1 yield a mean residence time of 7 ns at 27 °C and an orientational order parameter of 0.45. The effects of mutations around WAT1 on the kinetic parameters of RNase T1 are small but significant and probably relate to the dynamics of the active site.

The tumor suppressor p53 is conformationally unstable at physiological temperature. Even the activated p53Δ30 variant, which lacks the self-inhibiting carboxy terminal domain, has a half-life of only 8 min at 37 °C in vitro. We have developed a genetic approach to identify p53 variants that stabilize the active conformation. The human p53Δ30 gene was randomly mutated, and the resulting library was expressed in Escherichia coli under conditions that apparently denatured the parental protein. Stable p53 variants were identified based on their ability to specifically bind a p53 consensus site. The initial thermostable variants were randomly recombined by DNA shuffling, and substitutions that were functionally additive or synergistic were identified in a second more stringent round of screening. The DNA binding activity of N239Y/N268D/E336V p53Δ30 variant has a half-life of 100 min at 37 °C, 12 times longer than that of the parental protein. The thermostable variants should be more amenable to crystallographic studies and more effective in gene therapies than the wild-type protein.

The effects of a number of cryoprotectants on the kinetic and structural properties of glycogen phosphorylase b have been investigated. Kinetic studies showed that glycerol, one of the most commonly used cryoprotectants in X-ray crystallographic studies, is a competitive inhibitor with respect to substrate glucose-1-P with an apparent Ki value of 3.8% (v/v). Cryogenic experiments, with the enzyme, have shown that glycerol binds at the catalytic site and competes with glucose analogues that bind at the catalytic site, thus preventing the formation of complexes. This necessitated a change in the conditions for cryoprotection in crystallographic binding experiments with glycogen phosphorylase. It was found that 2-methyl-2,4-pentanediol (MPD), polyethylene glycols (PEGs) of various molecular weights, and dimethyl sulfoxide (DMSO) activated glycogen phosphorylase b to different extents, by stabilizing its most active conformation, while sucrose acted as a noncompetitive inhibitor and elethylene glycol as an uncompetitve inhibitor with respect to glucose-1-P. A parallel experimental investigation by X-ray crystallography showed that, at 100 K, both MPD and DMSO do not bind at the catalytic site, do not induce any significant conformational change on the enzyme molecule, and hence, are more suitable cryoprotectants than glycerol for binding studies with glycogen phosphorylase.

In the past few years, a new generation of fold recognition methods has been developed, in which the classical sequence information is combined with information obtained from secondary structure and, sometimes, accessibility predictions. The results are promising, indicating that this approach may compete with potential-based methods (Rost B et al., 1997, J Mol Biol 270:471–480). Here we present a systematic study of the different factors contributing to the performance of these methods, in particular when applied to the problem of fold recognition of remote homologues.

Our results indicate that secondary structure and accessibility prediction methods have reached an accuracy level where they are not the major factor limiting the accuracy of fold recognition. The pattern degeneracy problem is confirmed as the major source of error of these methods. On the basis of these results, we study three different options to overcome these limitations: normalization schemes, mapping of the coil state into the different zones of the Ramachandran plot, and post-threading graphical analysis.

Phage-display peptide library analysis of an anti-F actin polyclonal antibody identified 12 amino acid residues of actin that appear, in its X-ray crystal structure, to be grouped together in a surface accessible conformational epitope. Phage epitope mapping was carried out by isolating immune complexes containing members of the J404 nonapeptide phage-display library formed in diluted antiserum and isolated on a protein A affinity matrix. Immunoreactive clones were grown as plaques, replica plated onto nitrocellulose, and labeled with anti-actin immune serum. One hundred and forty-four positively staining clones identified in this way were sequenced. Of these, 54 displayed peptides with sequence similarities. When the most abundantly selected sequence, KQTWQQLWD, was produced as a synthetic peptide and derivatized to ovalbumin, the complex was strongly recognized by the antiserum on Western blots and inhibited the binding of the antibody to immobilized F-actin by 60%. A scrambled version of this sequence WQDKWLQTQ, when coupled to ovalbumin, was not recognized by the antiserum and minimally inhibited binding of antiserum to immobilized F-actin by 10%. KQTWQQLWD contained four residues that corresponded, in frame, to a highly conserved six residue region of the chicken β-actin sequence 351TFQQMW356 (identical residues are shown in bold). Examination of the rabbit skeletal muscle X-ray crystal structure suggested that within a 15 Å radius of W356, nine additional residues were arranged on the actin surface in such a way that they could be mimicked by several of the selected phage sequences with root-mean-square deviation fits of 2.1–2.5 Å. We conclude that phage-display analysis can provide information about the relative location of amino acids on the surfaces of proteins using antibody imprints of the protein surface structure.

We describe the results of a procedure for maximizing the number of sequences that can be reliably linked to a protein of known three-dimensional structure. Unlike other methods, which try to increase sensitivity through the use of fold recognition software, we only use conventional sequence alignment tools, but apply them in a manner that significantly increases the number of relationships detected. We analyzed 11 genomes and found that, depending on the genome, between 23 and 32% of the ORFs had significant matches to proteins of known structure. In all cases, the aligned region consisted of either >100 residues or >50% of the smaller sequence. Slightly higher percentages could be attained if smaller motifs were also included. This is significantly higher than most previously reported methods, even those that have a fold-recognition component. We survey the biochemical and structural characteristics of the most frequently occurring proteins, and discuss the extent to which alignment methods can realistically assign function to gene products.

The 2.5 Å crystal structure of the full length human placental isoform of the Gly12 to Val mutant Cdc42 protein (Cdc42(G12V)) bound to both GDP/Mg2+ and GDPNH2 (guanosine-5′-diphospho-β-amidate) is reported. The crystal contains two molecules in the asymmetric unit, of which one has bound GDP/Mg2+, while the other has bound GDPNH2 without a Mg2+ ion. Crystallization of the protein was induced via hydrolysis of the Cdc42·GppNHp complex by the presence of contaminating alkaline phosphatase activity in combination with the crystallization conditions. This prompted us to compare the binding characteristics of GDPNH2 vs. GDP. The amino group of GDPNH2 drastically reduces the affinity to Cdc42 in comparison with that of GDP, causes the loss of the Mg2+ ion, and apparently also increases the conformational flexibility of the protein as seen in the crystal. Both the switch I and switch II regions are visible in the electron density of the GDP-bound molecule, but not in the molecule bound to GDPNH2. The C-terminus containing the CaaX-motif is partly ordered in both molecules due to an intramolecular disulfide bond formed between Cys105/Cys188 and Cys305/Cys388, respectively.

We have investigated the solution conformation of the functionally relevant C-terminal extremes of α- and β-tubulin, employing the model recombinant peptides RL52α3 and RL33β6, which correspond to the amino acid sequences 404–451(end) and 394–445(end) of the main vertebrate isotypes of α- and β-tubulin, respectively, and synthetic peptides with the α-tubulin(430–443) and β-tubulin(412–431) internal sequences. α(404–451) and β(394–445) are monomeric in neutral aqueous solution (as indicated by sedimentation equilibrium), and have circular dichroism (CD) spectra characteristic of nearly disordered conformation, consistent with low scores in peptide helicity prediction. Limited proteolysis of β(394–445) with subtilisin, instead of giving extensive degradation, resulted in main cleavages at positions Thr409–Glu410 and Tyr422–Gln423–Gln424, defining the proteolysis resistant segment 410–422, which corresponds to the central part of the predicted β-tubulin C-terminal helix. Both recombinant peptides inhibited microtubule assembly, probably due to sequestration of the microtubule stabilizing associated proteins. Trifluoroethanol (TFE)-induced markedly helical CD spectra in α(404–451) and β(394–445). A substantial part of the helicity of β(394–445) was found to be in the CD spectrum of the shorter peptide β(412–431) with TFE. Two-dimensional 1H-NMR parameters (nonsequential nuclear Overhauser effects (NOE) and conformational CαH shifts) in 30% TFE permitted to conclude that about 25% of α(404–451) and 40% of β(394–451) form well-defined helices encompassing residues 418–432 and 408–431, respectively, flanked by disordered N- and C-segments. The side chains of β(394–451) residues Leu418, Val419, Ser420, Tyr422, Tyr425, and Gln426 are well defined in structure calculations from the NOE distance constraints. The apolar faces of the helix in both α and β chains share a characteristic sequence of conserved residues Ala,Met(+4),Leu(+7),Tyr(+11). The helical segment of α(404–451) is the same as that described in the electron crystallographic model structure of αβ-tubulin, while in β(394–451) it extends for nine residues more, supporting the possibility of a functional coil → helix transition at the C-terminus of β-tubulin. These peptides may be employed to construct model complexes with microtubule associated protein binding sites.

The relative orientations of adjacent structural elements without many well-defined NOE contacts between them are typically poorly defined in NMR structures. For apo-S100B(ββ) and the structurally homologous protein calcyclin, the solution structures determined by conventional NMR exhibited considerable differences and made it impossible to draw unambiguous conclusions regarding the Ca2+-induced conformational change required for target protein binding. The structure of rat apo-S100B(ββ) was recalculated using a large number of constraints derived from dipolar couplings that were measured in a dilute liquid crystalline phase. The dipolar couplings orient bond vectors relative to a single-axis system, and thereby remove much of the uncertainty in NOE-based structures. The structure of apo-S100B(ββ) indicates a minimal change in the first, pseudo-EF-hand Ca2+ binding site, but a large reorientation of helix 3 in the second, classical EF-hand upon Ca2+ binding.

Members of the Ly-6/uPAR protein family share one or several repeat units of the Ly-6/uPAR domain that is defined by a distinct disulfide bonding pattern between 8 or 10 cysteine residues. The Ly-6/uPAR protein family can be divided into two subfamilies. One comprises GPI-anchored glycoprotein receptors with 10 cysteine residues. The other subfamily includes the secreted single-domain snake and frog cytotoxins, and differs significantly in that its members generally possess only eight cysteines and no GPI-anchoring signal sequence. We report the purification and structural characterization of human SLURP-1 (secreted mammalian Ly-6/uPAR related protein 1) from blood and urine peptide libraries. SLURP-1 is encoded by the ARS (component B)-81/s locus, and appears to be the first mammalian member of the Ly-6/uPAR family lacking a GPI-anchoring signal sequence. A phylogenetic analysis based on the SLURP-1 primary protein structure revealed a closer relationship to the subfamily of cytotoxins. Since the SLURP-1 gene maps to the same chromosomal region as several members of the Ly-6/uPAR subfamily of glycoprotein receptors, it is suggested that both biologically distinct subfamilies might have co-evolved from local chromosomal duplication events.

The structural origins of the specificity of the neurophysin hormone-binding site for an aromatic residue in peptide position 2 were explored by analyzing the binding of a series of peptides in the context of the crystal structure of liganded neurophysin. A new modeling method for describing the van der Waals surface of binding sites assisted in the analysis. Particular attention was paid to the unusually large (5 kcal/mol) difference in binding free energy between Phe and Leu in position 2, a value representing more than three times the maximum expected based on hydrophobicity alone, and additionally remarkable since modeling indicated that the Leu side chain was readily accommodated by the binding pocket. Although evidence was obtained of a weak thermodynamic linkage between the binding interactions of the residue 2 side chain and of the peptide α-amino group, two factors are considered central. (1) The bound Leu side chain can establish only one-third of the van der Waals contacts available to a Phe side chain. (2) The bound Phe side chain appears to be additionally stabilized relative to Leu by more favorable dipole and induced dipole interactions with nonaromatic polar and sulfur ligands in the binding pocket, as evidenced by examination of its interactions in the pocket, analysis of the detailed energetics of transfer of Phe and Leu side chains from water to other phases, and comparison with thermodynamic and structural data for the binding of residue 1 side chains in this system. While such polar interactions of aromatic rings have been previously observed, the present results suggest their potential for significant thermodynamic contributions to protein structure and ligand recognition.

The effect of methylurea, N,N′-dimethylurea, ethylurea, and butylurea as well as guanidine hydrochloride (GuHCl), urea and pH on the thermal stability, structural properties, and preferential solvation changes accompanying the thermal unfolding of ribonuclease A (RNase A) has been investigated by differential scanning calorimetry (DSC), UV, and circular dichroism (CD) spectroscopy. The results show that the thermal stability of RNase A decreases with increasing concentration of denaturants and the size of the hydrophobic group substituted on the urea molecule. From CD measurements in the near- and far-UV range, it has been observed that the tertiary structure of RNase A melts at about 3 °C lower temperature than its secondary structure, which means that the hierarchy in structural building blocks exists for RNase A even at conditions at which according to DSC and UV measurements the RNase A unfolding can be interpreted in terms of a two-state approximation. The far-UV CD spectra also show that the final denatured states of RNase A at high temperatures in the presence of different denaturants including 4.5 M GuHCl are similar to each other but different from the one obtained in 4.5 M GuHCl at 25 °C. The concentration dependence of the preferential solvation change ΔΓ23, expressed as the number of cosolvent molecules entering or leaving the solvation shell of the protein upon denaturation and calculated from DSC data, shows the same relative denaturation efficiency of alkylureas as other methods.

The objective of this study was to evaluate the suitability of the WW domain as a desirable model system to understand the folding and stability of an isolated three-stranded antiparallel β-sheet structure. The WW domain was subjected to thermal and chaotropic denaturation/reconstitution utilizing a variety of biophysical methods. This three-stranded sheet folds reversibly and cooperatively utilizing both urea and GdnHCl as denaturants; however, the denatured state retains structure in the form of a hydrophobic cluster involving at least one aromatic side chain. In contrast to chaotropic denaturation, thermal denaturation appears to be more complete and may be a two state process. The suitability of the WW domain for future studies aimed at understanding the kinetics and thermodynamics of antiparallel β-sheet folding clearly emerges from this initial study. The most exciting and significant result in this manuscript is the finding that the chaotropic denatured state of WW has a hydrophobic cluster as discerned by near-UV CD evidence. The role that the denatured state plays in the folding and stability of a three-stranded β-sheets, and its capacity for preventing aggregation may be particularly important and is the subject of ongoing studies.

Here we describe the NMR conformational study of a 20-residue linear peptide designed to fold into a monomeric three-stranded antiparallel β-sheet in aqueous solution. Experimental and statistical data on amino acid β-turn and β-sheet propensities, cross-strand side-chain interactions, solubility criteria, and our previous experience with β-hairpins were considered for a rational selection of the peptide sequence. Sedimentation equilibrium measurements and NMR dilution experiments provide evidence that the peptide is monomeric. Analysis of 1H and 13C-NMR parameters of the peptide, in particular NOEs and chemical shifts, and comparison with data obtained for two 12-residue peptides encompassing the N- and C-segments of the designed sequence indicates that the 20-residue peptide folds into the expected conformation. Assuming a two-state model, the exchange kinetics between the β-sheet and the unfolded peptide molecules is in a suitable range to estimate the folding rate on the basis of the NMR linewidths of several resonances. The time constant for the coil-β-sheet transition is of the order of several microseconds in the designed peptide. Future designs based on this peptide system are expected to contribute greatly to our knowledge of the many factors involved in β-sheet formation and stability.

NAD+ facilitates high-yield reactivation of clostridial glutamate dehydrogenase (GDH) after unfolding in urea. The specificity of this effect has been explored by using analogues and fragments of NAD+. The adenine portion, unlike the nicotinamide portion, is important for reactivation. Alteration in the nicotinamide portion, in acetylpyridine adenine dinucleotide, has little effect, whereas loss of the 6-NH2 substitution on the adenine ring, in 6-deamino NAD, diminishes the effectiveness of the nucleotide in promoting refolding. Also ADP-ribose, lacking nicotinamide, promotes reactivation whereas NMN-phosphoribose, lacking the adenine, does not.

Of the smaller fragments, those containing an adenosine moiety, and especially those with one or more phosphate groups, impede the refolding ability of NAD+, and are able to bind to the folding intermediate though unable to facilitate refolding.

These results are interpreted in terms of the known 3D structure for clostridial glutamate dehydrogenase. It is assumed that the refolding intermediate has a more or less fully formed NAD+-binding domain but a partially disordered substrate-binding domain and linking region. Binding of NAD+ or ADP-ribose appears to impose new structural constraints that result in completion of the correct folding of the second domain, allowing association of enzyme molecules to form the native hexamer.

We have performed a detailed study of methanol-induced conformational transitions of horse heart apomyoglobin (apoMb) to investigate the existence of the compact and expanded denatured states. A combination of far- and near-ultraviolet circular dichroism, NMR spectroscopy, and small-angle X-ray scattering (SAXS) was used, allowing a phase diagram to be constructed as a function of pH and the methanol concentration. The phase diagram contains four conformational states, the native (N), acid-denatured (UA), compact denatured (IM), and expanded helical denatured (H) states, and indicates that the compact denatured state (IM) is stable under relatively mild denaturing conditions, whereas the expanded denatured states (UA and H) are realized under extreme conditions of pH (strong electric repulsion) or alcohol concentration (weak hydrophobic interaction). The results of this study, together with many previous studies in the literature, indicate the general existence of the compact denatured states not only in the salt-pH plane but also in the alcohol-pH plane. Furthermore, to determine the general feature of the H conformation we used several proteins including ubiquitin, ribonuclease A, α-lactalbumin, β-lactoglobulin, and Streptomyces subtilisin inhibitor (SSI) in addition to apoMb. SAXS studies of these proteins in 60% methanol showed that the H states of these all proteins have expanded and nonglobular conformations. The qualitative agreement of the experimental data with computer-simulated Kratky profiles also supports this structural feature of the H state.

We propose a new approach that permits correlation of specific domains defined by their primary sequence with their location in the structure of complex macromolecular aggregates. It is based on the combination of well-established structural analysis methods that incorporate the use of overlapping peptides on cellulose membranes for the isolation and purification of specific antibodies from a polyclonal antiserum. Monospecific antibodies to the connector protein of bacteriophage φ29 were isolated from polyclonal antisera using a new development of the spotscan method. These antibodies can be purified in quantities that allow antigenicity testing in enzyme-linked immunosorbent assays, Western blotting and immunoprecipitations, demonstrating the specificity of this isolation procedure. This approach has allowed us to generate direct antibody probes for immunoelectron microscopy mapping of different connector protein domains in a low resolution three-dimensional epitope map.

An apparent conservative mutation, Leu to Val, at the second residue of the rat liver mitochondrial aldehyde dehydrogenase (ALDH) presequence resulted in a precursor protein that was not imported into mitochondria. Additional mutants were made to substitute various amino acids with nonpolar side chains for Leu2. The Ile, Phe, and Trp mutants were imported to an extent similar to that of the native precursor, but the Ala mutant was imported only about one-fourth as well. It was shown that the N-terminal methionine was removed from the L2V mutant in a reaction catalyzed by methionine aminopeptidase. The N-terminal methionine of native pALDH and the other mutant presequences was blocked, presumably by acetylation. Because of the difference in co-translational modification, the L2V mutant sustained a significant loss in the available hydrophobic surface of the presequence. Import competence was restored to the L2V mutant when it was translated using a system that did not remove Met1. The removal of an Arg-Gly-Pro helix linker segment (residues 11–14) from the L2V mutant, which shifted three leucine residues toward the N-terminus, also restored import competence. These results lead to the conclusion that a minimum amount of hydrophobic surface area near the N-termini of mitochondrial presequences is an essential property to determine their ability to be imported. As a result, both electrostatic and hydrophobic components must be considered when trying to understand the interactions between precursor proteins and proteins of the mitochondrial import apparatus.