De novo backbone and sequence design of an idealized alpha/beta-barrel protein: evidence of stable tertiary structure

We have designed, synthesized, and characterized a 216 amino acid residue sequence encoding a putative idealized alpha/beta-barrel protein. The design was elaborated in two steps. First, the idealized backbone was defined with geometric parameters representing our target fold: a central eight parallel-stranded beta-sheet surrounded by eight parallel alpha-helices, connected together with short structural turns on both sides of the barrel. An automated sequence selection algorithm, based on the dead-end elimination theorem, was used to find the optimal amino acid sequence fitting the target structure. A synthetic gene coding for the designed sequence was constructed and the recombinant artificial protein was expressed in bacteria, purified and characterized. Far-UV CD spectra with prominent bands at 222nm and 208nm revealed the presence of alpha-helix secondary structures (50%) in fairly good agreement with the model. A pronounced absorption band in the near-UV CD region, arising from immobilized aromatic side-chains, showed that the artificial protein is folded in solution. Chemical unfolding monitored by tryptophan fluorescence revealed a conformational stability (DeltaG(H2O)) of 35kJ/mol. Thermal unfolding monitored by near-UV CD revealed a cooperative transition with an apparent T(m) of 65 degrees C. Moreover, the artificial protein did not exhibit any affinity for the hydrophobic fluorescent probe 1-anilinonaphthalene-8-sulfonic acid (ANS), providing additional evidence that the artificial barrel is not in the molten globule state, contrary to previously designed artificial alpha/beta-barrels. Finally, 1H NMR spectra of the folded and unfolded proteins provided evidence for specific interactions in the folded protein. Taken together, the results indicate that the de novo designed alpha/beta-barrel protein adopts a stable three-dimensional structure in solution. These encouraging results show that de novo design of an idealized protein structure of more than 200 amino acid residues is now possible, from construction of a particular backbone conformation to determination of an amino acid sequence with an automated sequence selection algorithm.     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.     

Hydrogen exchange demonstrates three domains in aldolase unfold sequentially.

Rabbit muscle aldolase is a homotetramer in which the subunits have a classical alpha/beta-barrel structure and Mr 39,212 Da. We have previously reported that aldolase incubated in 3 M urea has three unfolding domains distinguished by their different unfolding rates. The unfolding rates of these domains were determined from isotope patterns in the mass spectra of peptic fragments derived from aldolase incubated in 3 M urea and pulse labeled in (2)H2O. The present study extends this investigation to more thoroughly characterize the structures of these unfolding intermediates. Mass spectra of intact monomers had four envelopes of isotope peaks corresponding to four structural forms of aldolase. Analysis of the present results suggests that these structural forms consist of native aldolase and three forms that have one to three domains unfolded. The molecular masses of these four structural forms indicate that there are 107 residues in each of the three unfolding domains of aldolase. Present results also show that aldolase remains a tetramer in 4 M urea, even though hydrogen exchange and circular dichroism indicate that it has lost most of its secondary and tertiary structure. The abundances of unfolded domains, which were determined from mass spectra of either intact aldolase or its peptic fragments, were used to determine the abundances of specific, partially unfolded forms of aldolase. Kinetic modeling of the abundances of these structures suggests that these structures are formed sequentially as aldolase unfolds in urea.     

Estimation of protein secondary structure from the laser Raman amide I spectrum

A new method for estimating protein secondary structure from the laser Raman spectrum has been developed whereby the amide I Raman band of a protein is analyzed directly as a linear combination of amide I bands of proteins whose secondary structures are known. For 14 proteins, analyzed by removing each one from the reference set and analyzing its structure in terms of the remaining proteins, the average correlation coefficients between the Raman and X-ray diffraction estimates of helix, beta-strand, turn, and undefined were 0.98, 0.98, 0.82 and 0.35, respectively. Significant correlations were also observed for distinctions between alpha-helix (0.98) and disordered helix (0.82), and between parallel (0.82) and antiparallel (0.97) beta-sheets. The average standard deviation of these Raman estimates from the X-ray values is less than 4%. In addition, a singular value analysis of 20 Raman amide I spectra indicates that there may be as many as nine significant independent pieces of information present in the amide I region.     

Secondary structure of substance P bound to liposomes in organic solvents and in solution from Raman and CD spectroscopy.

The membrane lipid phase may be an important mediator of the peptide-receptor interaction. In order to understand the mechanism of this interaction, it is important to know the peptide structure, not only in the hydrophobic lipid bilayer environment, but also at the bilayer surface and in solution. To investigate this problem we have measured the secondary structure of the 11-residue neuropeptide substance P (SP) and its fragments in aqueous solutions, in membrane mimetic solvents, and associated with lipid bilayers using Raman and CD spectroscopy. Raman and CD spectra of SP bound to liposomes indicate a less than 20% helix content. We interpret these results to indicate that SP contains virtually no helix when bound to negatively charged liposomes. These spectra are similar to spectra of peptides in type I and III beta-turns. SP forms between 10 and 30% (1-3 residues) helical structure in sodium dodecyl sulfate micelles and less than 10% helix in methanol and trifluoroethanol. The binding of SP to negatively charged liposomes significantly changes the structure of the lipid acyl chains, decreasing order in some cases and increasing it in others. Raman spectra of SP in water indicates that SP near 30 mM forms an ensemble of structures in water that is distinct from completely unfolded peptide and from the aggregated beta-sheet form observed in saline solutions. We conclude from our CD results that methods used to quantitate secondary structure from CD spectra of short peptides cannot be used to distinguish between very short helical segments and beta-turns.     

Probing the folding capacity and residual structures in 1-79 residues fragment of staphylococcal nuclease by biophysical and NMR methods

1-79 residues SNase fragment (SNase79) has chain length containing a sequence for helix alpha(1), omega-loop, beta(I)-sheet, and partial beta(II)-sheet of native SNase. The incomplete "beta-barrel" structural region of SNase79 makes this fragment to be interested in investigation of its conformation. For this study, we use CD, fluorescence, and NMR spectroscopy to probe the folding capacity and the residual structures in SNase79. The optical spectra obtained for SNase79 and its mutants reveal the presence of retained capacity for folding of the fragment. The NMR derived (13)C(alpha) secondary chemical shifts, (3)J(NH-Halpha) coupling constants, amide-proton temperature coefficients, interresidue NOEs, and (15)N relaxation data determine the intrinsic propensities for helix- and turn- or beta-sheet-like conformations of SNase79, which is not the result of stabilizing inter-molecular interactions by oligomerization effects. The residual turn- and helix-like structures may serve as potential local nucleation sites, whereas the residual beta(I)-sheet-like structure can be regarded as a potential non-local nucleation site in the folding of SNase79. The intrinsic local and non-local interactions in these potential initiation sites are insufficient to stabilize the folding of SNase79 due to the shortage of relevant long-range interactions from other part of the fragment. The conformational ensemble of SNase79 is a highly heterogeneous collection of interconverting conformations having transiently populated helix- and beta-sheet- or turn-like structures.     

Secondary structure of human leukocyte interferon from Raman spectroscopy

Raman spectra were taken of human alpha (leukocyte) interferon subtype A (HuIFN-alpha A) purified from extracts of transformed Escherichia coli. Quantitative analysis of the conformationally sensitive amide I band indicates that IFN (interferon)-alpha A is 75 +/- 5% helical and 7 +/- 4% beta-strand. An independent analysis of the amide III spectrum indicates 71 +/- 5% helix and 10 +/- 6% beta-strand. These results differ with a recently proposed three-dimensional model based on secondary structure predictions derived from sequence and with circular dichroism measurements. The Raman spectrum of IFN-alpha A is compared with the spectra of several other helical proteins: hemerythrin, intestinal calcium-binding protein, melittin, and insulin.     

The design of idealized alpha/beta-barrels: analysis of beta-sheet closure requirements

The 8-fold parallel alpha/beta-barrel topology is encountered in proteins that display an impressive variety of functions, suggesting that this topology may be a rather nonspecific and stable folding motif. Consequently, this motif can be considered as an interesting framework to design novel proteins. It has been shown that the shape of the beta-sheet portion of the barrel can be approximated by a hyperboloid. This geometric object may therefore be used as a scaffold to construct an idealized eight-stranded beta-barrel. To facilitate the de novo design of such structures, a collection of modeling tools has been developed allowing secondary structure elements to be mapped onto the scaffold surface and rotation and translation operations to be performed about user defined axes while evaluating their contribution to the conformational energy of the system. These tools have been applied in a systematic study assessing the phi, psi requirements to design symmetric eight stranded beta barrels with optimal hydrogen bonding between adjacent beta-strands. It is observed that: (a) the beta-sheet structure can be closed without introducing irregular stagger between beta-strands and (b) the region of phi, psi dihedral angle space compatible with the formation of regular symmetric eight stranded beta-barrels coincides with the phi, psi region corresponding to average beta-strands in known protein structures, suggesting that barrel closure does not impose gross constraints on beta-strand geometry.     

The pattern of amino acid replacements in alpha/beta-barrels

The determinants of site-to-site variability in the rate of amino acid replacement in alpha/beta-barrel enzyme structures are investigated. Of 125 available alpha/beta-barrel structures, only 25 meet a variety of phylogenetic and statistical criteria necessary to ensure sufficient data for reliable analysis. These 25 enzyme structures (from a wide variety of taxa with diverse lifestyles in diverse habitats) differ greatly in size, number, and topology of domains in addition to the alpha/beta-barrel, quaternary structure, metabolic role, reaction catalyzed, presence of prosthetic groups, regulatory mechanisms, use of cofactors, and catalytic mechanisms. Yet, with the exception of ribulose-1,5-bisphosphate carboxylase, all structures have similar frequency distributions of amino acid replacement rates. Hence, site-specific variability in rates of evolution is largely independent of differences in biology, biochemistry, and molecular structure. A correlation between site-specific rate variation and (1) distance from the active site, (2) solvent accessibility, and (3) treating glycines in unusual main-chain conformations as a separate class, explains approximately half the causal variation. Secondary structure exerts little influence on the pattern and distribution of replacements. Additional domains and subunits, side-chain hydrogen bonds, unusual side-chain rotamers, nonplanar peptide bonds, strained main-chain conformations, and buried hydrophilic-charged residues contribute little to variability among sites because they are rare. Nonlinear models do not improve the fits. In several enzymes, deviations from the typical pattern of replacements suggest the possible action of natural selection. A statistical analysis shows that, in all cases, much of the remaining unexplained variation is not attributable to chance and that other, as yet unidentified, causal relations must exist.     

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.     

Catalytic and structural role of a metal-free histidine residue in bovine Cu-Zn superoxide dismutase

Cu-Zn superoxide dismutase (SOD) contains a conserved, metal-free His residue at an opening of the backbone beta-barrel in addition to six Cu- and/or Zn-bound His residues in the active site. We examined the protonation and hydrogen bonding state of the metal-free His residue (His41) in bovine SOD by UV Raman spectroscopy. Analysis of the His Raman intensity at 1406 cm(-1) in a D2O solution has shown that His41 has a pKa of 9.4, consistent with the NMR and X-ray structures at acidic to neutral pH, in which two imidazole nitrogen atoms of cationic His41 are hydrogen bonded to the main chain C=O groups of Thr37 and His118. Upon deprotonation of His41 at pH 9.4, the Thr37-His41-His118 hydrogen bond bridge breaks on the His118 side and SOD loses 70% of its activity. Concomitantly, hydrogen-deuterium exchange is accelerated for amide groups of beta-strands, indicating an increased conformational fluctuation of the beta-barrel. Thr37 and His41 are in direct contact with Leu36, whose hydrophobic side chain closes off the opening of the beta-barrel, while His118 is indirectly connected to Arg141 that assists the docking of superoxide to Cu. These Raman findings strongly suggest that the His41-mediated hydrogen bond bridge plays a crucial role in keeping the protein structure suitable for highly efficient catalytic reactions. The catalytic and structural role of His41 is consistent with the observation that the mutation of His43 in human SOD (equivalent to His41 in bovine SOD) to Arg largely reduces the dismutase activity and the protein structural stability.     

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.     

Structural model for interferons

Secondary structures of leucocyte alpha 1- and alpha 2-interferons and of fibroblast beta-interferon are calculated using the molecular theory of protein secondary structures. The common secondary structure calculated for alpha- and beta-interferons is used to predict the three-dimensional structures of fragments 1-110 and 111-166 of the chains (which are supposed to be quasi-independent domains). The predicted structure of the active domain I (1-110) is an 'up-and-down' tetrahelical complex (in which the second helix is shorter than the others and can be absent in alpha 1-interferon) similar to the mirror image of myohaemoerythrin. The predicted structure of domain II (111-166) is either a three-stranded beta-sheet screened from one side by two alpha-helices or a three-helical complex (similar to that in the N-domain of papain), the first structure being more consistent with the circular dichroism data of alpha-interferon and its C-end fragment.     

Predicting the three-dimensional structure of alpha- and beta-interferons

Secondary structures of leukocyte alpha 1- and alpha 2-interferons and of fibroblast beta-interferon are calculated using the molecular theory of protein secondary structures. The common secondary structure calculated for alpha- and beta-interferons is used to predict the three-dimensional structures of fragments 1-110 and 111-166 of the chains (which are supposed to be quasi-independent domains). The predicted structure of the active domain I (1-110) is an "up-and-down" tetrahelical complex (in which the second helix is shorter than the others and can be missed in alpha 1-interferon) similar to the mirror-image of myohaemoerythrin. The predicted structure of domain II (111-166) is either a three-stranded beta-sheet screened from one side by two alpha-helices or a three-helical complex (similar to that in the N-domain of papain), the first structure being better consistent with the circular dichroism data of alpha-interferon and its C-end fragment.     

The beta-barrel finder (BBF) program,allowing identification of outer membrane beta-barrel proteins encoded within prokaryotic genomes

Many outer membrane proteins (OMPs) in Gram-negative bacteria possess known beta-barrel three-dimensional (3D) structures. These proteins, including channel-forming transmembrane porins, are diverse in sequence but exhibit common structural features. We here report computational analyses of six outer membrane proteins of known 3D structures with respect to (1) secondary structure, (2) hydropathy, and (3) amphipathicity. Using these characteristics, as well as the presence of an N-terminal targeting sequence, a program was developed allowing prediction of integral membrane beta-barrel proteins encoded within any completely sequenced prokaryotic genome. This program, termed the beta-barrel finder (BBF) program, was used to analyze the proteins encoded within the Escherichia coli genome. Out of 4290 sequences examined, 118 (2.8%) were retrieved. Of these, almost all known outer membrane proteins with established beta-barrel structures as well as many probable outer membrane proteins were identified. This program should be useful for predicting the occurrence of outer membrane proteins in bacteria with completely sequenced genomes.     

Dynamics and reactivity of HbXL99 alpha. A cross-linked hemoglobin derivative

Resonance Raman spectroscopy, transient absorption, and fluroescence techniques have been employed to investigate the structure and dynamics of the alpha-cross-linked hemoglobin derivative, HbXL99 alpha. The resonance Raman spectra of the deoxy form of HbXL99 alpha are identical to those of native NbA (VFe-His approximately 222 cm-1), which exhibit a T-state (low affinity) structure regardless of solvent conditions. The resonance Raman spectra of the transient heme photoproduct resulting from CO photolysis from HbXL99 alpha appear to have structures intermediate between deoxy-T and ligand-bound R structures (VFe-His approximately 222 cm-1). Time-resolved resonance Raman data of HbXL99 alpha-CO show that complete CO recombination occurs after approximately 5 ms, with only a small amount of the CO-bound species reforming within approximately 200 ns (geminate recombination). Transient absorption spectra of HbXL99 alpha-O2 indicate that the extent of sub-nanosecond geminate recombination of O2 is also reduced in the cross-linked derivative relative to native HbA. The decrease in tryptophan fluorescence of HbXL99 alpha upon oxygenation further indicates that tertiary structural changes at the alpha 1-beta 2 interface upon ligation are apparently reduced, but not eliminated in the cross-linked derivative relative to HbA.     

INVESTIGATION OF SECONDARY STRUCTURES AND MACROMOLECULAR INTERACTIONS IN BACTERIOPHAGE P22 BY LASER RAMAN SPECTROSCOPY

Laser Raman spectra of the DNA bacteriophage P22 and of its precursor particles and related structures have been obtained using 514.5-nm excitation. The spectra show that P22 DNA exists in the B form both inside of the phage head and after extraction from the phage. The major coat protein (gp5) contains a secondary structure composed of 18% α-helix, 20% β-sheet and 62% irregular conformations. The scaffolding protein (gp8) in the phage prohead is substantially richer than gp5 in α-helical content. Among the amino acid residues which give prominent Raman lines, the spectra show that tryptophans are exposed to solvent and most tyrosines are hydrogen bonded to positive donor groups. The above features of phage DNA and protein structures are nearly invariant to changes in temperature up to 80°C, indicating a remarkable thermal stability of the phage head and its encapsulated DNA.     

New insight into the pH-dependent conformational changes in bovine beta-lactoglobulin from Raman optical activity.

We have studied the conformation of beta-lactoglobulin in aqueous solution at room temperature over the pH range approximately 2.0-9.0 using vibrational Raman optical activity (ROA). The ROA spectra clearly show that the basic up and down beta-barrel core is preserved over the entire pH range, in agreement with other studies. However, from the shift of a sharp positive ROA band at approximately 1268 to approximately 1294 cm(-1) on going from pH values below that of the Tanford transition, which is centered at pH approximately 7.5, to values above, the Tanford transition appears to be associated with changes in the local conformations of residues in loop sequences possibly corresponding to a migration into the alpha-helical region of the Ramachandran surface from a nearby region. These changes may be related to those detected in X-ray crystal structures which revealed that the Tanford transition is associated with conformational changes in loops which form a doorway to the interior of the protein. The results illustrate how the ability of ROA to detect loop and turn structure separately from secondary structure is useful for studying conformational plasticity in proteins.     

Predicting transmembrane beta-barrels and interstrand residue interactions from sequence

Transmembrane beta-barrel (TMB) proteins are embedded in the outer membrane of Gram-negative bacteria, mitochondria, and chloroplasts. The cellular location and functional diversity of beta-barrel outer membrane proteins (omps) makes them an important protein class. At the present time, very few nonhomologous TMB structures have been determined by X-ray diffraction because of the experimental difficulty encountered in crystallizing transmembrane proteins. A novel method using pairwise interstrand residue statistical potentials derived from globular (nonouter membrane) proteins is introduced to predict the supersecondary structure of transmembrane beta-barrel proteins. The algorithm transFold employs a generalized hidden Markov model (i.e., multitape S-attribute grammar) to describe potential beta-barrel supersecondary structures and then computes by dynamic programming the minimum free energy beta-barrel structure. Hence, the approach can be viewed as a "wrapping" component that may capture folding processes with an initiation stage followed by progressive interaction of the sequence with the already-formed motifs. This approach differs significantly from others, which use traditional machine learning to solve this problem, because it does not require a training phase on known TMB structures and is the first to explicitly capture and predict long-range interactions. TransFold outperforms previous programs for predicting TMBs on smaller (相似文献