Localization of hydrogen-bonds within modules in barnase

Proteins in eukaryotes are composed of structural units, each encoded by discrete exons. The protein module is one such structural unit; it has been defined as the least extended or the most compact contiguous segment in a globular domain. To elucidate roles of modules in protein evolution and folding, we examined roles of hydrogen bonds and hydrophobic cores, as related to the stability of these modules. For this purpose we studied barnase, a bacterial Rnase from Bacillus amylolique-faciens. Barnase is decomposed into at least six modules, M1–M6; the module boundaries are identified at amino acid residues 24, 52, 73, 88, and 98. Hydrogen bonds are localized mainly within each of the modules, with only a few between them, thereby indicating that their locations are designed to primarily stabilize each individual module. To obtain support for this notion, an analysis was made of hypothetical modules defined as segments starting at a center of one module and ending at the center of the following one. We found that the hydrogen bonds did not localize in each hypothetical module and that many formed between the hypothetical modules. The native conformations of modules of barnase may be specified predominantly by interactions within the modules. © 1993 Wiley-Liss, Inc.     

Foldability of barnase mutants obtained by permutation of modules or secondary structure units

Modules, defined as stable, compact structure units in a globular protein, are good candidates for the construction of novel foldable proteins by permutation. Here we decomposed barnase into six modules (M1-M6) and constructed 23 barnase mutants containing permutations of the internal four (M2-M5) out of six modules. Globular proteins can also be subdivided into secondary structure units based on the extended structures that control the mutual relationships of the modules. We also decomposed barnase into six secondary structure units (S1-S6) and constructed 21 barnase mutants containing permutations of the internal four (S2-S5) out of six secondary structure units. Foldability of these two types of mutants was assessed by means of circular dichroism, fluorescence, and 1H-NMR measurements. A total of 15 of 23 module mutants and 15 of 21 secondary structure unit mutants formed definite secondary structures, such as alpha-helix and beta-sheet, at 20 microM owing to intermolecular interactions, but most of them converted to random coil structures at a lower concentration (1 microM). Of the 44 mutants, only two, M3245 and S2543, gave distinct near-UV CD spectra. S2543 especially showed definite signal dispersion in the amide and methyl regions of the 1H-NMR spectrum, though M3245 did not. Furthermore, urea-induced unfolding of S2543 monitored by far-UV CD and fluorescence measurements showed a distinct cooperative transition. These results strongly suggest that S2543 takes partially folded conformations in aqueous solution. Our results also suggest that building blocks such as secondary structure units capable of taking different stable conformations by adapting themselves to the surrounding environment, rather than building blocks such as modules having a specified stable conformation, are required for the formation of foldable proteins. Therefore, the use of secondary structure units for the construction of novel globular proteins is likely to be an effective approach.     

Dynamic structures of globular proteins with respect to correlative movements of residues calculated in the normal mode analysis

Dynamic structures of globular proteins are studied on the basis of correlative movements of residues around their native conformations, which are computed by means of the normal mode analysis. To describe the dynamic structures of a protein, the core regions moving with strong positive or negative correlations to other regions of the polypeptide chain are detected from the correlation maps of the movements of residues. Such core regions are different, according to the definition, from the regions defined from a geometrical point of view, such as secondary structures, domains, modules, and so on. The core regions are actually detected for four proteins, myoglobin, Bence-Jones protein, flavodoxin, and hen egg-white lysozyme, with different folding types from each other. The results show that some of them coincide with the secondary structures, domains, or modules, but others do not. Then, the dynamic structure of each protein is discussed in terms of the dynamic cores detected, as compared with the secondary structures, domains, and modules.     

Folding and Stability of Chimeric Immunofusion VL-Barstar

A chimeric protein VL-barstar that comprises the VL domain of anti-human ferritin monoclonal antibody F11 and barstar, the naturally occurring inhibitor of bacterial RNase barnase, has been constructed for study of structure-function characteristics of chimeric immunoglobulin fused proteins. Such chimeric constructs may be potentially employed for development of bivalent/bispecific antibodies on the basis of the high affinity interaction between barstar and barnase (the association constant is about 10(14) M(-1)). We have developed a protocol for VL-barstar expression in E. coli and purification and refolding from inclusion bodies that yields a homogeneous and soluble form of this protein. Differential scanning calorimetry in combination with fluorescence and CD spectroscopy revealed that the VL-barstar formed well-resolved ordered secondary and compact tertiary structures. However, partial loss of tertiary interactions resulted in low stability of the recombinant protein and the lack of functional activity of the two constituent modules. These conformational features suggest that the protein might be referred to the class of native molten globules, which comprises partially unfolded conformations stabilized under physiological conditions. Since individually expressed VL domain and barstar retain completely folded conformation and stable spatial structure, the incomplete folding of the chimeric protein may be attributed to interaction between heterologous domains, which appears at the folding stage preceding formation of a system of tertiary interactions in both structural modules. The results provide evidence for non-native interactions between heterologous modules that may occur in chimeric proteins composed of taxonomically distinct fusion partners.     

Characterization of β-turns in cyclic hexapeptides in solution by fourier transform IR spectroscopy

The β-turn represents a structural element frequently encountered in globular proteins. However, in spite of various theoretical and experimental studies the ir signature bands of pure β-turns are still not established beyond doubt. Although considerable information exists now on the ir spectra of β-helical and β-sheet structures, the lack of knowledge concerning turn structures in general, and that of β-turns in particular, presents a major uncertainty in the estimation of global protein secondary structures from ir spectroscopic data. To obtain more specific information about the characteristic amide bands in β-turns, we report herein an ir spectroscopic analysis of a series of five cyclic pseudo-hexapeptides known to form β-turns from previous CD and nmr studies [A. Perczel, M. Hollósi, B. M. Foxman, and G. D. Fasman (1991) Journal of the American Chemical Society, Volume 113, pp. 9772-9784 ]. We show here that in these cyclic peptides the amide groups involved in β-turns that comprise a ten-membered hydrogen-bonded ring (and represent the first H-bond pair in a β-sheet), give rise to characteristic amide I bands in the range 1638–1646 cm^?1, with the exact position depending on the solvent and the nature of the side-chain substituents. © 1993 John Wiley & Sons, Inc.     

A mini-protein designed by removing a module from barnase: molecular modeling and NMR measurements of the conformation.

A globular domain can be decomposed into compact modules consisting of contiguous 10-30 amino acid residues. The correlation between modules and exons observed in different proteins suggests that each module was encoded by an ancestral exon and that modules were combined into globular domains by exon fusion. Barnase is a single domain RNase consisting of 110 amino acid residues and was decomposed into six modules. We designed a mini-protein by removing the second module, M2, from barnase in order to gain an insight into the structural and functional roles of the module. In the molecular modeling of the mini-protein, we evaluated thermodynamic stability and aqueous solubility together with mechanical stability of the model. We chemically synthesized a mini-barnase with (15)N-labeling at 10 residues, whose corresponding residues in barnase are all found in the region around the hydrophobic core. Circular dichroism and NMR measurements revealed that mini-barnase takes a non-random specific conformation that has a similar hydrophobic core structure to that of barnase. This result, that a module could be deleted without altering the structure of core region of barnase, supports the view that modules act as the building blocks of protein design.     

Foldability, enzymatic activity, and interacting ability of barnase mutants obtained by permutation of secondary structure units

Barnase, a well-characterized ribonuclease, has been decomposed into six modules (M1-M6) or secondary structure units (S1-S6). We have studied the foldability and activity of the barnase mutants obtained by permutation of the four internal modules (M2-M5) or secondary structure units (S2-S5) to investigate whether permutation of these building blocks is a useful way to create foldable and/or functional proteins. In this study, we found that one of the secondary structure unit mutants was expressed in Escherichia coli only when His102 was substituted by alanine, which is a catalytic residue of wild-type barnase. This mutant (S2354H102A) had ordered conformations, which unfolded cooperatively during urea-induced unfolding experiments. S2354H102A interacted with other barnase mutants to show a distinct RNase activity, although its own activity was quite weak. This interaction was specific, because S2354H102A interacted with only barnase mutants having His 102 and certain orders of the secondary structure units giving a distinct RNase activity. These results suggest that secondary structure units permuted in barnase mutants maintain their intrinsic "interacting ability" that is used for the folding of wild-type barnase, and the units can form certain conformations that complement those of the appropriate counterparts. Seven of 23 secondary structure unit mutants and only 2 of 23 module mutants had RNase activity. On the basis of the results of analyses of foldability and RNase activity of the mutants performed in this and previous studies, we conclude that secondary structure units are more suitable than modules as building blocks to create novel foldable and/or functional proteins in the case of barnase.     

Conformational characterization of designed minibarnase

We have designed a minibarnase by removing one module from barnase, a bacterial RNase from Bacillus amyloliquefaciens. Barnase, consisting of 110 amino acid residues, is decomposed into six modules, M1-M6. Module is defined as a peptide segment consisting of contiguous amino acid residues that makes a small compact conformation within a globular domain. To understand the role of module in protein architecture, we analyzed NMR and CD spectra of a minibarnase, which lacked 26 amino acid residues corresponding to module M2. We demonstrated the formation of hydrophobic cores in the minibarnase similar to those of barnase. Although its conformational stability against acids and heat was reduced in comparison with barnase, the minibarnase retained cooperative folding character (two-state folding). Therefore, the folding of the minibarnase consisting of modules M1 and M3-M6 is independent to some extent of module M2. This finding may be useful for future module-based protein design.     

Refolding simulations of an isolated fragment of barnase into a native-like β hairpin: Evidence for compactness and hydrogen bonding as concurrent stabilizing factors

Experimental evidence and theoretical models both suggest that protein folding is initiated within specific fragments intermittently adopting conformations close to that found in the protein native structure. These folding initiation sites encompassing short portions of the protein are ideally suited for study in isolation by computational methods aimed at peering into the very early events of folding. We have used Molecular Dynamics (MD) technique to investigate the behavior of an isolated protein fragment formed by residues 85 to 102 of barnase that folds into a β hairpin in the protein native structure. Three independent MD simulations of 1.3 to 1.8 ns starting from unfolded conformations of the peptide portrayed with an all-atom model in water were carried out at gradually decreasing temperature. A detailed analysis of the conformational preferences adopted by this peptide in the course of the simulations is presented. Two of the unfolded peptide conformations fold into a hairpin characterized by native and a larger bulk of nonnative interactions. Both refolding simulations substantiate the close relationship between interstrand compactness and hydrogen bonding network involving backbone atoms. Persistent compactness witnessed by side-chain interactions always occurs concomitantly with the formation of backbone hydrogen bonds. No highly populated conformations generated in a third simulation starting from the remotest unfolded conformer relative to the native structure are observed. However, nonnative long-range and medium-range contacts with the aromatic moiety of Trp94 are spotted, which are in fair agreement with a former nuclear magnetic resonance study of a denaturing solution of an isolated barnase fragment encompassing the β hairpin. All this lends reason to believe that the 85–102 barnase fragment is a strong initiation site for folding. Proteins 29:212–227, 1997. © 1997 Wiley-Liss, Inc.     

A molecular dynamics study of the pentacyclo-undecane (PCU) cage polypeptides of the type Ac-3Ala-Cage-3Ala-NHMe

An extensive exploration of the conformational space of the seven-residue peptide sequences, Ac-Ala-Ala-Ala-Cage-Ala-Ala-Ala-NHMe and the model peptide Ac-Ala-Ala-Ala-Ala-Ala-Ala-Ala-NHMe, was carried out using single trajectories of molecular dynamics (MD) in the solution phase using the periodic boundary conditions. Our MD studies revealed that the majority of the motifs of the PCU cage peptide exist as type I–III β-turns along with their mirror conformations, viz. type I′–III′ β-turns. This peptide sequence adopted a U-shaped backbone, with alpha-helical characteristics. The results reported here provide further evidence that the PCU cage amino acid exhibits C_7eq, C_7aq, α_R and α_L conformations in aqueous solution.     

Stability Order of Basic Peptide Conformations Reflected by Density Functional Theory

Ab initio density functional methods at the B3LYP/6-311+G(d, p) and 6-31G(d) levels were performed on several basic peptide conformations representing typical elements of secondary structure (β-sheets, β- and γ-turns). The results are compared with those from Hartree-Fock and MP2 correlation energy calculations. Whereas the geometries of the structures are well described at all approximation levels, there are considerable discrepancies of the stability orders. Contrary to the Hartree-Fock calculations, the correlation energy methods provide the more compact structures with intramolecular hydrogen bonds distinctly favoured over extended conformations when comparing the energy differences. However, due to considerable compensation of correlation energy and entropy contributions, the stability order at the Gibbs free energy level closely corresponds to that at the Hartree-Fock level.     

The role of protein sequence and amino acid composition in amyloid formation: scrambling and backward reading of IAPP amyloid fibrils

The specific functional structure of natural proteins is determined by the way in which amino acids are sequentially connected in the polypeptide. The tight sequence/structure relationship governing protein folding does not seem to apply to amyloid fibril formation because many proteins without any sequence relationship have been shown to assemble into very similar β-sheet-enriched structures. Here, we have characterized the aggregation kinetics, seeding ability, morphology, conformation, stability, and toxicity of amyloid fibrils formed by a 20-residue domain of the islet amyloid polypeptide (IAPP), as well as of a backward and scrambled version of this peptide. The three IAPP peptides readily aggregate into ordered, β-sheet-enriched, amyloid-like fibrils. However, the mechanism of formation and the structural and functional properties of aggregates formed from these three peptides are different in such a way that they do not cross-seed each other despite sharing a common amino acid composition. The results confirm that, as for globular proteins, highly specific polypeptide sequential traits govern the assembly pathway, final fine structure, and cytotoxic properties of amyloid conformations.     

A new computational model for protein folding based on atomic solvation.

A new model for calculating the solvation energy of proteins is developed and tested for its ability to identify the native conformation as the global energy minimum among a group of thousands of computationally generated compact non-native conformations for a series of globular proteins. In the model (called the WZS model), solvation preferences for a set of 17 chemically derived molecular fragments of the 20 amino acids are learned by a training algorithm based on maximizing the solvation energy difference between native and non-native conformations for a training set of proteins. The performance of the WZS model confirms the success of this learning approach; the WZS model misrecognizes (as more stable than native) only 7 of 8,200 non-native structures. Possible applications of this model to the prediction of protein structure from sequence are discussed.     

Critical amino acids responsible for converting specificities of proteins and for enhancing enzyme evolution are located around β-turn potentials: Data-based prediction

Various reports have described that amino acid substitutions can alter substrate, positional, inhibitory, and target gene specificities of proteins. By using the method of Chou and Fasman, the present work predicts that critical amino acids for converting these specificities are located around β-turns. Residues responsible for the alterations of substrate specificities of trypsin,l-lactate dehydrogenase, aspartate aminotransferase, β-lactamase, and cytochrome P-450 are found to exist within regions predicted as β-turns. The ratios of hydroxylation and oxygenation positions of substrates by cytochrome P-450 and lipoxygenase, respectively, are varied by changes of the protein structures, probably around turn conformations. Inhibitory specificities of bovine pancreatic trypsin inhibitor and α₁-antitrypsin and target gene specificity of glucocorticoid receptor are converted by changing turn structures. Occurrence of β-turn probabilities can be predicted around the amino acid alteration positions of an evolutionally antecedent protein of a nylon degradation enzyme. These findings will have relevance to work on protein engineering and enzyme evolution.     

Module-intron correlation and intron sliding in family F/10 xylanase genes

Xylanases are classified into two families, numbered F/10 and G/11 according to the similarity of amino acid sequences of their catalytic domain (Henrissat, B., Bairoch, A., 1993. New families in the classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem. J. 293, 781-788). Three-dimensional structure of the catalytic domain of the family F/10 xylanase was reported (White, A., Withers, S.G., Gilkes, N.R., Rose, D.R., 1994. Crystal structure of the catalytic domain of the beta-1,4-glycanase Cex from Cellulomonas fimi. Biochemistry 33, 12546-12552). The domain was decomposed into 22 modules by centripetal profiles (Go, M., Nosaka, M., 1987. Protein architecture and the origin of introns. Cold Spring Harbor Symp. Quant. Biol. 52, 915-924; Noguti, T., Sakakibara, H., Go, M., 1993. Localization of hydrogen-bonds within modules in barnase. Proteins 16, 357-363). A module is a contiguous polypeptide segment of amino acid residues having a compact conformation within a globular domain. Collected 31 intron sites of the family F/10 xylanase genes from fungus were found to be correlated to module boundaries with considerable statistical force (p values <0.001). The relationship between the intron locations and protein structures provides supporting evidence for the ancient origin of introns, because such a relationship cannot be expected by random insertion of introns into eukaryotic genes, but it rather suggests pre-existence of introns in the ancestral genes of prokaryotes and eukaryotes. A phylogenetic tree of the fungal and bacterial xylanase sequences made two clusters; one includes both the bacterial and fungal genes, but the other consists of only fungal genes. The mixed cluster of bacterial genes without introns and the fungal genes with introns further supports the ancient origin of introns. Comparison of the conserved base sequences of introns indicates that sliding of a splice site occurred in Aspergillus kawachii gene by one base from the ancestral position. Substrate-binding sites of xylanase are localized on eight modules, and introns are found at both termini of six out of these functional modules. This result suggests that introns might play a functional role in shuffling the exons encoding the substrate-binding modules.     

Size-exclusion high performance liquid chromatography of native trypsinogen, the denatured protein, and partially refolded molecules. Further evidence that non-native disulfide bonds are dominant in refolding the completely reduced protein

Size-exclusion high performance liquid chromatography was used to compare the Stokes radius of the mixed disulfide of trypsinogen refolded for 10 min with the Stokes radius of denatured trypsinogen in high concentrations of urea. After folding for 10 min, rechromatography of a collection of sequential fractions of an initial separation showed that the fractions display microheterogeneity as seen in the value of the Stokes radius of each fraction. These intermediate species differed in their Stokes radius, and each had a globular structure cross-linked by disulfide bonds. In contrast, when trypsinogen with the native disulfides intact was equilibrated at different concentrations of urea (0-8 M), a progressive increase in Stokes radius was observed with extent of unfolding. Rechromatography of a series of fractions collected at a specific urea concentration showed that each had the same Stokes radius as the fraction in the initial separation. Urea-denatured trypsinogen and partially refolded trypsinogen must therefore differ in the disulfide pairing that links regions of the polypeptide chain. These observations support the suggestion that non-native disulfide bonds are responsible for the many stable conformations that form early in the folding of the mixed disulfide of trypsinogen (Light, A., and Higaki, J.N. (1987) Biochemistry 26, 5556-5564). These intermediates initially are loose structures (large Stokes radius) that become more compact with time (decreasing Stokes radius). The intermediates must therefore undergo a continuing disulfide interchange until native disulfides form late in the process when the stable conformation of the native molecule is reached.     

Protein anatomy: structure and function of peptide fragments corresponding to the secondary structure units of barnase

Globular proteins can be decomposed into several modules or secondary structure units. It is useful to investigate the functions of such structural units in order to understand the folding units of proteins. In our previous work, barnase was divided into six peptide fragments corresponding to modules, and some of them were shown to have RNA-binding and RNase activity [Yanagawa, et al. (1993) J. Biol. Chem. 268, 5861-5865]. Barnase mutant proteins obtained by permutation of the structural units also had RNase activity [Tsuji, T. et al. (1999) J. Mol. Biol. 286, 1581-1596]. Here we investigated the structure and function of peptide fragments corresponding to secondary structure units of barnase. The results of circular dichroism spectroscopy indicated that some of the peptide fragments form helical structures in aqueous solutions containing over 30% 2,2,2-trifluoroethanol, and the S6 (94-110) peptide fragment is induced to form a beta-sheet structure in the presence of RNA. The S6 peptide fragment forms aggregate complexes with RNA. Electron microscopic analysis showed that the aggregate complexes were comprised of filaments. These results indicate that not only modules but also secondary structure units dissected from a globular protein have functional and structure-forming capabilities.     

Descriptive statistics of disallowed regions and various protein secondary structures in the context of studying twisted β-hairpins

A detailed analysis of polypeptide-chain backbone conformations was carried out for polypeptide-chain segments adjacent to β-turn regions, including the sites of disallowed conformations. A cross comparison of conformations was performed for disallowed regions of the Ramachandran plot and main types of β-turns and adjacent secondary structures. Based on the results, disallowed region 2 (II, II') in the Ramachandran plot was shown to coincide mainly with β-hairpins and, more exactly, twisted β-hairpins. The frequency of residues with angles ?_i, ψ_i that fall in region 2 (II, II') in the latter is 140 times higher than in common β-hairpins.     

Theory of protein molecule self-organization. II. A comparison of calculated thermodynamic parameters of local secondary structures with experiments

Constants of the helix–coil transition for all natural amino acid residues are evaluated on the basis of thermodynamic parameters obtained in paper I of this series. The specific effects at the termini of the helices are also considered as well as the parameters controlling the formation of β-bends in the unfolded protein chain. Evaluated s constants of the helix–coil transition agree with the experimental data on helix–coil transitions of synthetic polypeptides in water. Only a very qualitative correlation exists between s constants (both experimental and theoretical) and the occurrence of corresponding residues in internal turns of α-helices in globular proteins: residues with s > 1 occur in helices as a rule more often than residues with s < 1. At the same time a direct correlation is demonstrated between theoretical parameters of residue incorporation into α-helical termini and β-bends in an unfolded polypeptide chain and the occurrence of residues in corresponding positions of the globular protein secondary structures.