Algorithms for prediction of alpha-helical and beta-structural regions in globular proteins

Algorithms are suggested for identifying α-helical and β-structural regions in native globular proteins. α-Helical and β-structural regions are predicted, with accuracy of ~80 and ~85% respectively, for 25 proteins, the three-dimensional structures of which have been determined by X-ray diffraction crystallography. Secondary structure is predicted in 25 proteins with unknown three-dimensional structure.     

Automatic identification of secondary structure in globular proteins

A computer program is used to analyse automatically and objectively the atomic co-ordinates of a large number of globular proteins in order to identify the regions of α-helix, β-sheet and reverse-turn secondary structure. Several different criteria for the assignment of secondary structure are tested for accuracy, reproducibility and efficiency. The most successful criterion, which is based on patterns of peptide hydrogen bonds, inter-C^α distances and inter-C^α torsion angles, is used to find the secondary structure of all the proteins studied. The accuracy of the derived assignments is assessed by comparing them with the secondary structure reported in the literature for each protein. The reliability of the methods is assessed by comparing the secondary structures derived from the independently determined sets of co-ordinates available for some proteins.We provide the first objective and consistent compilation of α-helix, β-sheet and reverse-turn secondary structure in almost all globular proteins of known tertiary structure. These data will be invaluable for analysing the relative tendencies of different amino acids to occur in different types of secondary structure, for analysing the regularity of the secondary structure itself, and for analysing how the pieces of secondary structure fit together to form the globular tertiary structure of each protein.     

Origins of globular structure in proteins

Thermodynamic incompatibility of polymers in a common solvent is possibly a driving force for formation and evolution of globular protein structures. Folding of polypeptide chains leads to a decrease in both excluded volume of molecules and chemical differences between surfaces of globular molecules with chemical information hidden in the hydrophobic interior. Folding of polypeptide chains results in 'molecular or thermodynamic mimicry' of globular proteins and in at least more than 10-fold higher phase separation threshold values of mixed protein solutions compared to those of classical polymers. Unusually high co-solubility might be necessary for efficient biological functioning of proteins, e.g. enzymes, enzyme inhibitors, etc.     

The role of local hydrogen bonds in formation of irregular regions of globular protein polypeptide chains]

An assumption is made on the substantial role of local hydrogen bonds in formation of irregular regions of globular protein polypeptide chains. The statistics of the amino acid composition of irregular regions is examined from this point of view. A statistical analysis of side group-backbone hydrogen bonds is carried out for three proteins: alpha-chy-motrypsin, lysozyme and myoglobin. It is shown that short side groups participate in formation of local hydrogen bonds more often than long ones. Conformations of amino acid residues in the first and the last positions are studied in beta-bends of 9 proteins. It is shown that over 70% of these residues are in conformations corresponding to the formation of local hydrogen bonds of three types: backbone-backbone, side groupbackbone, backbone-water molecule-backbone. Thus, the participation of the cooperative hydrogen-bonding network in stabilization of beta-bends is demonstrated.     

Modular evolution and the origins of symmetry: reconstruction of a three-fold symmetric globular protein

The high frequency of internal structural symmetry in common protein folds is presumed to reflect their evolutionary origins from the repetition and fusion of ancient peptide modules, but little is known about the primary sequence and physical determinants of this process. Unexpectedly, a sequence and structural analysis of symmetric subdomain modules within an abundant and ancient globular fold, the β-trefoil, reveals that modular evolution is not simply a relic of the ancient past, but is an ongoing and recurring mechanism for regenerating symmetry, having occurred independently in numerous existing β-trefoil proteins. We performed a computational reconstruction of a β-trefoil subdomain module and repeated it to form a newly three-fold symmetric globular protein, ThreeFoil. In addition to its near perfect structural identity between symmetric modules, ThreeFoil is highly soluble, performs multivalent carbohydrate binding, and has remarkably high thermal stability. These findings have far-reaching implications for understanding the evolution and design of proteins via subdomain modules.     

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.     

Selective observation of the disordered import signal of a globular protein by in-cell NMR: The example of frataxins

We have exploited the capability of in-cell NMR to selectively observe flexible regions within folded proteins to carry out a comparative study of two members of the highly conserved frataxin family which are found both in prokaryotes and in eukaryotes. They all contain a globular domain which shares more than 50% identity, which in eukaryotes is preceded by an N-terminal tail containing the mitochondrial import signal. We demonstrate that the NMR spectrum of the bacterial ortholog CyaY cannot be observed in the homologous E. coli system, although it becomes fully observable as soon as the cells are lysed. This behavior has been observed for several other compact globular proteins as seems to be the rule rather than the exception. The NMR spectrum of the yeast ortholog Yfh1 contains instead visible signals from the protein. We demonstrate that they correspond to the flexible N-terminal tail indicating that this is flexible and unfolded. This flexibility of the N-terminus agrees with previous studies of human frataxin, despite the extensive sequence diversity of this region in the two proteins. Interestingly, the residues that we observe in in-cell experiments are not visible in the crystal structure of a Yfh1 mutant designed to destabilize the first helix. More importantly, our results show that, in cell, the protein is predominantly present not as an aggregate but as a monomeric species.     

The relationship between molecular weight and quaternary structure of globular proteins

The relationship between the molecular weight and the number of subunits in oligomeric globular proteins consisting of identical subunits has been analyzed. It has been shown that the molecular weights of the subunits are distributed about a mean value of 48,000 and consequently that the molecular weights of the native oligomeric proteins are distributed in clearly distinguishable molecular weight regions. This observation allows the probability of a particular oligomeric structure to be predicted from a measurement of the oligomer molecular weight alone, which is useful in a number of types of study of protein structure, particularly comparative studies. Calculations have been performed which suggest that there is no thermodynamic limitation, in terms of the subunit interactions themselves, to the size of an oligomeric protein with a given number of subunits. Rather, an individual polypeptide chain itself has inherent size limitations, which consequently limits the molecular weight of the corresponding oligomer.     

Optical activity of polypeptides and proteins

Using methods described in a previous publication the optical activity of a number of polypeptides and proteins has been calculated. The systems included the α-helix, the two β-structures, polyproline I, polyproline II, collagen and collagen models, and poly-N-methylalanine. In addition to these orderded structures, calculations were also performed on the α, β and nonperiodic regions of myoglobin, lysozyme, ribonuclease-S and β-chymotrypsin. The α and β structures in prteins differ from the polypeptide models by being very short and partially disordered. It is concluded that the 222-nm band of the α-helix is a good method for detecting helices in proteins but that the 207-and 191-nm bands of the helix will not fit a linear superposition model. The circular dichroism of the so-called β regions of proteins differs markedly from that for ideal β structure because of breakdownin symmetry. As a result estimates of β-structure in proteins based on polypeptide models are not likely to be quantitative. The theoretical methods give an adequate account of the optical activity of all the ordered polypeptides except polyproline II and collagen and (by inference) the nonperiodic chains in the various proteins. This difficulty is the remaining barrier to a complete theory of the optical activity of the polypeptide backbone in globular proteins.     

Noncooperative temperature melting of a globular protein without specific tertiary structure: acid form of bovine carbonic anhydrase B

Heat denaturation of native globular proteins is a cooperative process usually connected with the melting of the main part of their regular secondary structure. In this paper, a noncooperative temperature-induced melting of the regular secondary structure in the carbonic anhydrase B at pH 2.6 in heavy water is observed by ir spectroscopy. The molecules of carbonic anhydrase B in an acid medium, unlike the native ones, do not have a specific tertiary structure. Nevertheless, the β-structure content is about the same in both of these states. A temperature-induced noncooperative melting process takes place from 10 to 67°C with a decrease of the antiparallel β-form content by about one third. The remaining part of the β-form melts with a more intensive heat absorption, with a maximum at 87°C. The whole melting process is practically reversible. We assume that the observed noncooperative process displays a general property of a new type of structural state of the globular protein—the “molten globule state.”     

Detection of disordered regions in globular proteins using 13C‐detected NMR

Characterization of disordered regions in globular proteins constitutes a significant challenge. Here, we report an approach based on ¹³C‐detected nuclear magnetic resonance experiments for the identification and assignment of disordered regions in large proteins. Using this method, we demonstrate that disordered fragments can be accurately identified in two homologs of menin, a globular protein with a molecular weight over 50 kDa. Our work provides an efficient way to characterize disordered fragments in globular proteins for structural biology applications.     

Variation of amino acid properties in all-beta globular and outer membrane protein structures

Discriminating outer membrane (OM) proteins from globular proteins is an important task. The structural analysis of β-strands dominating globular (all-β) proteins and OM proteins provides useful insight to distinguish between them. In this work, we analyze the characteristic features of the 20 amino acid residues in all-β and OM proteins. We set up numerical indices for several properties of amino acid residues, such as, conformational parameters, surrounding hydrophobicity, accessible surface area and reduction in accessibility, and inter-residue contacts. We found that all the aromatic residues prefer to be in β-strands of both globular and OM proteins. The surrounding hydrophobicity of aromatic and non-polar amino acid residues in globular proteins is significantly higher than that of OM proteins. The residues Trp, Arg, Phe and Gln show a remarkable difference of reduction in accessibility between all-β globular (βG) and OM proteins. The positively charged residues, Lys and Arg in the membrane part of OM proteins have more number of contacts than globular proteins. Further, the behavior of the 20 amino acid residues in β-strand segments of globular and OM proteins have been discussed. The parameters developed in this work can be used for identifying transmembrane β-strands in OM proteins and for discriminating βG proteins from OM proteins.     

Letter: Recognition of structural domains in globular proteins

The polypeptide chain of many globular proteins is folded into two or more structural domains which are spatially separated from one another and are frequently based on different architectural principles. The conformation of such structural domains is, in general, conserved more faithfully (assuming divergent evolution) than the amino acid sequence. The preservation of structure may originate in the three-dimensional requirements for the conservation of essential functions. A method is discussed for the recognition and comparison of structural domains.     

A comparative study of the relationship between protein structure and beta-aggregation in globular and intrinsically disordered proteins

A growing number of proteins are being identified that are biologically active though intrinsically disordered, in sharp contrast with the classic notion that proteins require a well-defined globular structure in order to be functional. At the same time recent work showed that aggregation and amyloidosis are initiated in amino acid sequences that have specific physico-chemical properties in terms of secondary structure propensities, hydrophobicity and charge. In intrinsically disordered proteins (IDPs) such sequences would be almost exclusively solvent-exposed and therefore cause serious solubility problems. Further, some IDPs such as the human prion protein, synuclein and Tau protein are related to major protein conformational diseases. However, this scenario contrasts with the large number of unstructured proteins identified, especially in higher eukaryotes, and the fact that the solubility of these proteins is often particularly good. We have used the algorithm TANGO to compare the beta aggregation tendency of a set of globular proteins derived from SCOP and a set of 296 experimentally verified, non-redundant IDPs but also with a set of IDPs predicted by the algorithms DisEMBL and GlobPlot. Our analysis shows that the beta-aggregation propensity of all-alpha, all-beta and mixed alpha/beta globular proteins as well as membrane-associated proteins is fairly similar. This illustrates firstly that globular structures possess an appreciable amount of structural frustration and secondly that beta-aggregation is not determined by hydrophobicity and beta-sheet propensity alone. We also show that globular proteins contain almost three times as much aggregation nucleating regions as IDPs and that the formation of highly structured globular proteins comes at the cost of a higher beta-aggregation propensity because both structure and aggregation obey very similar physico-chemical constraints. Finally, we discuss the fact that although IDPs have a much lower aggregation propensity than globular proteins, this does not necessarily mean that they have a lower potential for amyloidosis.     

The ESR determination of protein alpha-helicity.

Experimental data are summarized on decay kinetics of uv-induced paramagnetic centers in globular proteins. The PC decay rate in the absence of oxygen in α helical regions is shown to differ always by about two orders from that in nonhelical regions of the polypeptide chain. This effect can be utilized to estimate the degree of α helicity in the dry state, in solution, and in the frozen solution of proteins.     

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.     

Monte Carlo studies on equilibrium globular protein folding. I. Homopolymeric lattice models of beta-barrel proteins

Dynamic Monte Carlo studies have been performed on various diamond lattice models of β-proteins. Unlike previous work, no bias toward the native state is introduced; instead, the protein is allowed to freely hunt through all of phase space to find the equilibrium conformation. Thus, these systems may aid in the elucidation of the rules governing protein folding from a given primary sequence; in particular, the interplay of short- vs long-range interaction can be explored. Three distinct models (A? C) were examined. In model A, in addition to the preference for trans (t) over gauche states (g⁺ and g^?) (thereby perhaps favoring β-sheet formation), attractive interactions are allowed between all nonbonded, nearest neighbor pairs of segments. If the molecules possess a relatively large fraction of t states in the denatured form, on cooling spontaneous collapse to a well-defined β-barrel is observed. Unfortunately, in model A the denatured state exhibits too much secondary structure to correctly model the globular protein collapse transition. Thus in models B and C, the local stiffness is reduced. In model B, in the absence of long-range interactions, t and g states are equally weighted, and cooperativity is introduced by favoring formation of adjacent pairs of nonbonded (but not necessarily parallel) t states. While the denatured state of these systems behaves like a random coil, their native globular structure is poorly defined. Model C retains the cooperativity of model B but allows for a slight preference of t over g states in the short-range interactions. Here, the denatured state is indistinguishable from a random coil, and the globular state is a well-defined β-barrel. Over a range of chain lengths, the collapse is well represented by an all-or-none model. Hence, model C possesses the essential qualitative features observed in real globular proteins. These studies strongly suggest that the uniqueness of the globular conformation requires some residual secondary structure to be present in the denatured state.     

Effect of the fluctuations in electron density within a globular protein on the excess small-angle X-ray scattering

Excess small angle X-ray scattering in solvents of differing electron density has been calculated from the crystal structures obtained for rubredoxin, trypsin inhibitor, myogen, ferricytochrome c₂, ribonuclease S, lysozyme, nuclease, myoglobin, α-chymotrypsin, elastase, subtilisin, carboxypeptidase A, thermolysin, methemoglobin, deoxyhemoglobin, and a single polypeptide chain of M₄ lactate dehydrogenase. The scattering curves for each protein can be reproduced by the sum of three curves, with the weighting of the three curves depending on the electron density of the solvent. The radius of gyration obtained from the small angle X-ray scattering by globular proteins in aqueous solution will usually exceed the values defined by the shape of the macromolecule. Deviations for certain of the proteins cited are calculated to be as large as 6%. These deviations arise from the tendency for the amino acid residues with low electron density to be situated closer to the center of the protein than the amino acid residues of high electron density. An upper limit of 19% is obtained for the discrepancy between the radius of gyration defined by the shape of a spherical globular protein of typical amino acid composition and the apparent radius of gyration measured for that protein in water by small angle X-ray scattering.     

Biosynthesis and cocoon-export of a recombinant globular protein in transgenic silkworms

A gene construct was made by fusing the coding sequence of the red fluorescent protein (DsRed) to the exon 2 of the fibrohexamerin gene (fhx), that encodes a subunit of fibroin, the major silk protein of the silkworm Bombyx mori. The fusion gene was inserted into a piggyBac vector to establish a series of transgenic lines. The expression of the transgene was monitored during the course of larval life and was found restricted to the posterior silk gland cells as the endogenous fhx gene, in all the selected transgenic lines. The exogenous polypeptide was secreted into the lumen of the posterior silk gland together with fibroin, and further exported with the silk proteins as a foreign constituent of the cocoon fiber. The capacity of DsRed to emit fluorescence in the air-dried silk thread led to show that the recombinant protein was distributed over the whole length of the fiber. A remarkable property of the system lies in the localization of the globular protein at the periphery of the silk thread, allowing its rapid and easy recovery in aqueous solutions, without dissolving fibroin. The procedure represents a novel and promising strategy for the production of massive recombinant proteins of biomedical and pharmaceutical interest, with reduced cost.     

Far-infrared spectra of some globular proteins

The far-infrared absorption spectrum (40–400 cm^?1) of solid pellets and films of several globular proteins (lysozyme, myoglobin, hemoglobin, serum albumin, ribonuclease, chymotrypsinogen, subtilisin) and of some representative polypeptides [nylon 66, poly (γ-benzyl L -glutamate)] have been investigated by using a Michelson interferometer. While polypeptides are known to present several peaks which can be assigned mostly to hydrogen-bond modes, all the investigated globular proteins display only one broad, intense baud in the 100–200 cm^?1 region. The origin of this band, which persists even after denaturation or partial digestion, is discussed.