Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Crystal structure of a Staphylococcus aureus protein A domain complexed with the Fab fragment of a human IgM antibody: Structural basis for recognition of B-cell receptors and superantigen activity

Staphylococcus aureus produces a virulence factor, protein A (SpA), that contains five homologous Ig-binding domains. The interactions of SpA with the Fab region of membrane-anchored Igs can stimulate a large fraction of B cells, contributing to lymphocyte clonal selection. To understand the molecular basis for this activity, we have solved the crystal structure of the complex between domain D of SpA and the Fab fragment of a human IgM antibody to 2.7-Å resolution. In the complex, helices II and III of domain D interact with the variable region of the Fab heavy chain (VH) through framework residues, without the involvement of the hypervariable regions implicated in antigen recognition. The contact residues are highly conserved in human VH3 antibodies but not in other families. The contact residues from domain D also are conserved among all SpA Ig-binding domains, suggesting that each could bind in a similar manner. Features of this interaction parallel those reported for staphylococcal enterotoxins that are superantigens for many T cells. The structural homology between Ig VH regions and the T-cell receptor Vβ regions facilitates their comparison, and both types of interactions involve lymphocyte receptor surface remote from the antigen binding site. However, T-cell superantigens reportedly interact through hydrogen bonds with T-cell receptor Vβ backbone atoms in a primary sequence-independent manner, whereas SpA relies on a sequence-restricted conformational binding with residue side chains, suggesting that this common bacterial pathogen has adopted distinct molecular recognition strategies for affecting large sets of B and T lymphocytes.

Solution structure of a 142-residue recombinant prion protein corresponding to the infectious fragment of the scrapie isoform

The scrapie prion protein (PrPSc) is the major, and possibly the only, component of the infectious prion; it is generated from the cellular isoform (PrPC) by a conformational change. N-terminal truncation of PrPSc by limited proteolysis produces a protein of ≈142 residues designated PrP 27–30, which retains infectivity. A recombinant protein (rPrP) corresponding to Syrian hamster PrP 27–30 was expressed in Escherichia coli and purified. After refolding rPrP into an α-helical form resembling PrPC, the structure was solved by multidimensional heteronuclear NMR, revealing many structural features of rPrP that were not found in two shorter PrP fragments studied previously. Extensive side-chain interactions for residues 113–125 characterize a hydrophobic cluster, which packs against an irregular β-sheet, whereas residues 90–112 exhibit little defined structure. Although identifiable secondary structure is largely lacking in the N terminus of rPrP, paradoxically this N terminus increases the amount of secondary structure in the remainder of rPrP. The surface of a long helix (residues 200–227) and a structured loop (residues 165–171) form a discontinuous epitope for binding of a protein that facilitates PrPSc formation. Polymorphic residues within this epitope seem to modulate susceptibility of sheep and humans to prion disease. Conformational heterogeneity of rPrP at the N terminus may be key to the transformation of PrPC into PrPSc, whereas the discontinuous epitope near the C terminus controls this transition.

A CDKN2-like polymorphism in Xiphophorus LG V is associated with UV-B-induced melanoma formation in platyfish–swordtail hybrids

The genetic basis of spontaneous melanoma formation in spotted dorsal (Sd) Xiphophorus platyfish–swordtail hybrids has been studied for decades, and is adequately explained by a two-gene inheritance model involving a sex-linked oncogene, Xmrk, and an autosomal tumor suppressor, DIFF. The Xmrk oncogene encodes a receptor tyrosine kinase related to EGFR; the nature of the DIFF tumor suppressor gene is unknown. We analyzed the genetic basis of UV-B-induced melanoma formation in closely related, spotted side platyfish–swordtail hybrids, which carry a different sex-linked pigment pattern locus, Sp. We UV-irradiated spotted side Xiphophorus platyfish–swordtail backcross hybrids to induce melanomas at frequencies 6-fold higher than occur spontaneously in unirradiated control animals. To identify genetic determinants of melanoma susceptibility in this UV-inducible Xiphophorus model, we genotyped individual animals from control and UV-irradiated experimental regimes using allozyme and DNA restriction fragment length polymorphisms and tested for joint segregation of genetic markers with pigmentation phenotype and UV-induced melanoma formation. Joint segregation results show linkage of a CDKN2-like DNA polymorphism with UV-B-induced melanoma formation in these hybrids. The CDKN2-like polymorphism maps to Xiphophorus linkage group V and exhibits recombination fractions with ES1 and MDH2 allozyme markers consistent with previous localization of the DIFF tumor suppressor locus. Our results indicate that the CDKN2-like sequence we have cloned and mapped is a candidate for the DIFF tumor suppressor gene.

Major recent and independent changes in levels and patterns of expression have occurred at the b gene, a regulatory locus in maize

The b locus encodes a transcription factor that regulates the expression of genes that produce purple anthocyanin pigment. Different b alleles are expressed in distinct tissues, causing tissue-specific anthocyanin production. Understanding how phenotypic diversity is produced and maintained at the b locus should provide models for how other regulatory genes, including those that influence morphological traits and development, evolve. We have investigated how different levels and patterns of pigmentation have evolved by determining the phenotypic and evolutionary relationships between 18 alleles that represent the diversity of b alleles in Zea mays. Although most of these alleles have few phenotypic differences, five alleles have very distinct tissue-specific patterns of pigmentation. Superimposing the phenotypes on the molecular phylogeny reveals that the alleles with strong and distinctive patterns of expression are closely related to alleles with weak expression, implying that the distinctive patterns have arisen recently. We have identified apparent insertions in three of the five phenotypically distinct alleles, and the fourth has unique upstream restriction fragment length polymorphisms relative to closely related alleles. The insertion in B-Peru has been shown to be responsible for its unique expression and, in the other two alleles, the presence of the insertion correlates with the phenotype. These results suggest that major changes in gene expression are probably the result of large-scale changes in DNA sequence and/or structure most likely mediated by transposable elements.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

The crystal structure of the nitrogen regulation fragment of the yeast prion protein Ure2p

The yeast nonchromosomal gene [URE3] is due to a prion form of the nitrogen regulatory protein Ure2p. It is a negative regulator of nitrogen catabolism and acts by inhibiting the transcription factor Gln3p. Ure2p residues 1–80 are necessary for prion generation and propagation. The C-terminal fragment retains nitrogen regulatory activity, albeit somewhat less efficiently than the full-length protein, and it also lowers the frequency of prion generation. The crystal structure of this C-terminal fragment, Ure2p(97–354), at 2.3 Å resolution is described here. It adopts the same fold as the glutathione S-transferase superfamily, consistent with their sequence similarity. However, Ure2p(97–354) lacks a properly positioned catalytic residue that is required for S-transferase activity. Residues within this regulatory fragment that have been indicated by mutational studies to influence prion generation have been mapped onto the three-dimensional structure, and possible implications for prion activity are discussed.

Half tetrad analysis in alfalfa using multiple restriction fragment length polymorphism markers.

A maximum likelihood approach of half tetrad analysis (HTA) based on multiple restriction fragment length polymorphism (RFLP) markers was developed. This procedure estimates the relative frequencies of 2n gametes produced by mechanisms genetically equivalent to first division restitution (FDR) or second division restitution and simultaneously locates the centromere within a linkage group of RFLP marker loci. The method was applied to the diploid alfalfa clone PG-F9 (2n = 2x = 16) previously selected because of its high frequency of 2n egg production. HTA was based on four RFLP loci for which PG-F9 was heterozygous with codominant alleles that were absent in the tetraploid tester. Models including three linked and one unlinked RFLP loci were developed and tested. Results of the HTA showed that PG-F9 produced 6% FDR and 94% second division restitution 2n eggs. Information from a marker locus belonging to one linkage group was used to more precisely locate the centromere on a different linkage group. HTA, together with previous cytological analysis, indicated that in PG-F9, FDR 2n eggs are likely produced by diplospory, a mechanism common among apomictic species. The occurrence of FDR 2n eggs in plant species and their importance for crop evolution and breeding is discussed together with the potential applicability of multilocus HTA in the study of reproductive mutants.

Chromogranin B (secretogranin I) promotes sorting to the regulated secretory pathway of processing intermediates derived from a peptide hormone precursor.

Chromogranin B (CgB, secretogranin I) is a widespread constituent of neuroendocrine secretory granules whose function is unknown. To determine whether CgB affects the sorting of peptide hormone and neuropeptide precursors to secretory granules, we overexpressed CgB in AtT-20 cells, which exhibit an only moderate capacity to sort proopiomelanocortin and proteolytic fragments derived therefrom. In mock-transfected AtT-20 cells, a substantial proportion of newly synthesized proopiomelanocortin and its two primary proteolytic products generated in the trans-Golgi network, the N-terminal 23-kDa fragment containing adrenocorticotropin and the C-terminal beta-lipotropin fragment, was secreted via the constitutive pathway. Two- to three-fold overexpression of CgB markedly reduced the constitutive secretion of the 23-kDa fragment, but not beta-lipotropin and tripled the amount of adrenocorticotropin generated and stored in secretory granules. Our results indicate the existence of neuroendocrine-specific helper proteins which promote the sorting from the trans-Golgi network to secretory granules of certain processing intermediates derived from peptide hormone and neuropeptide precursors and demonstrate that CgB functions as such.

Tumor dormancy and cell signaling: anti-mu-induced apoptosis in human B-lymphoma cells is not caused by an APO-1-APO-1 ligand interaction.

Signal transduction initiated by crosslinking of antigen-specific receptors on T- and B-lymphoma cells induces apoptosis. In T-lymphoma cells, such crosslinking results in upregulation of the APO-1 ligand, which then interacts with induced or constitutively expressed APO-1, thereby triggering apoptosis. Here we show that crosslinking the membrane immunoglobulin on human lymphoma cells (Daudi) (that constitutively express APO-1) does not induce synthesis of APO-1 ligand. Further, a noncytotoxic fragment of anti-APO-1 antibody that blocks T-cell-receptor-mediated apoptosis in T-lymphoma cells does not block anti-mu-induced apoptosis. Hence, in B-lymphoma cells, apoptosis induced by signaling via membrane IgM is not mediated by the APO-1 ligand.

Null mutation of endothelin receptor type B gene in spotting lethal rats causes aganglionic megacolon and white coat color.

Mutations in the gene encoding the endothelin receptor type B (EDNRB) produce congenital aganglionic megacolon and pigment abnormalities in mice and humans. Here we report a naturally occurring null mutation of the EDNRB gene in spotting lethal (sl) rats, which exhibit aganglionic megacolon associated with white coat color. We found a 301-bp deletion spanning the exon 1-intron 1 junction of the EDNRB gene in sl rats. A restriction fragment length polymorphism caused by this deletion perfectly cosegregates with the sl phenotype. The deletion leads to production of an aberrantly spliced EDNRB mRNA that lacks the coding sequence for the first and second putative transmembrane domains of the G-protein-coupled receptor. Radioligand binding assays revealed undetectable levels of functional EDNRB in tissues from homozygous sl/sl rats. We conclude that EDNRB plays an essential role in the normal development of two neural crest-derived cell lineages, epidermal melanocytes and enteric neurons, in three mammalian species--humans, mice, and rats. The EDNRB-deficient rat may also prove valuable in defining the postnatal physiologic role of this receptor.

Crystal structure of the complex of a catalytic antibody Fab fragment with a transition state analog: structural similarities in esterase-like catalytic antibodies.

The x-ray structure of the complex of a catalytic antibody Fab fragment with a phosphonate transition-state analog has been determined. The antibody (CNJ206) catalyzes the hydrolysis of p-nitrophenyl esters with significant rate enhancement and substrate specificity. Comparison of this structure with that of the uncomplexed Fab fragment suggests hapten-induced conformational changes: the shape of the combining site changes from a shallow groove in the uncomplexed Fab to a deep pocket where the hapten is buried. Three hydrogen-bond donors appear to stabilize the charged phosphonate group of the hapten: two NH groups of the heavy (H) chain complementarity-determining region 3 (H3 CDR) polypeptide chain and the side-chain of histidine-H35 in the H chain (His-H35) in the H1 CDR. The combining site shows striking structural similarities to that of antibody 17E8, which also has esterase activity. Both catalytic antibody ("abzyme") structures suggest that oxyanion stabilization plays a significant role in their rate acceleration. Additional catalytic groups that improve efficiency are not necessarily induced by the eliciting hapten; these groups may occur because of the variability in the combining sites of different monoclonal antibodies that bind to the same hapten.

X-ray structure of clotting factor IXa: active site and module structure related to Xase activity and hemophilia B.

Hereditary deficiency of factor IXa (fIXa), a key enzyme in blood coagulation, causes hemophilia B, a severe X chromosome-linked bleeding disorder afflicting 1 in 30,000 males; clinical studies have identified nearly 500 deleterious variants. The x-ray structure of porcine fIXa described here shows the atomic origins of the disease, while the spatial distribution of mutation sites suggests a structural model for factor X activation by phospholipid-bound fIXa and cofactor VIIIa. The 3.0-A-resolution diffraction data clearly show the structures of the serine proteinase module and the two preceding epidermal growth factor (EGF)-like modules; the N-terminal Gla module is partially disordered. The catalytic module, with covalent inhibitor D-Phe-1I-Pro-2I-Arg-3I chloromethyl ketone, most closely resembles fXa but differs significantly at several positions. Particularly noteworthy is the strained conformation of Glu-388, a residue strictly conserved in known fIXa sequences but conserved as Gly among other trypsin-like serine proteinases. Flexibility apparent in electron density together with modeling studies suggests that this may cause incomplete active site formation, even after zymogen, and hence the low catalytic activity of fIXa. The principal axes of the oblong EGF-like domains define an angle of 110 degrees, stabilized by a strictly conserved and fIX-specific interdomain salt bridge. The disorder of the Gla module, whose hydrophobic helix is apparent in electron density, can be attributed to the absence of calcium in the crystals; we have modeled the Gla module in its calcium form by using prothrombin fragment 1. The arched module arrangement agrees with fluorescence energy transfer experiments. Most hemophilic mutation sites of surface fIX residues occur on the concave surface of the bent molecule and suggest a plausible model for the membrane-bound ternary fIXa-FVIIIa-fX complex structure: fIXa and an equivalently arranged fX arch across an underlying fVIIIa subdomain from opposite sides; the stabilizing fVIIIa interactions force the catalytic modules together, completing fIXa active site formation and catalytic enhancement.