Quasirandom distribution of α and β regions in globular proteins

Distributions of α and β regions in globular proteins among clusters containing different numbers of adjacent α-helices, adjacent β regions, and “overlapping” βαβ units are considered. It is shown that these distributions do not differ greatly from what can be expected for random distributions of α and β regions along protein chains. In particular, it is shown that the amounts of relatively long α, β and βαβ clusters (which provide the basis for the conventional classification of globular proteins or domains into α, β, or α/β types) in random sequences also do not differ very much from those in real globular proteins. It follows that the possibility of structural classification of globular proteins (domains) does not imply the existence of a correlation in protein primary structure. This possibility exists even in random sequences of amino acid residues and therefore may not be the result of biological evolution.     

Importance of long-range interactions in protein folding

Long-range interactions play an active role in the stability of protein molecules. In this work, we have analyzed the importance of long-range interactions in different structural classes of globular proteins in terms of residue distances. We found that 85% of residues are involved in long-range contacts. The residues occurring in the range of 4-10 residues apart contribute more towards long-range contacts in all-alpha proteins while the range is 11-20 in all-beta proteins. The hydrophobic residues Cys, Ile and Val prefer the 11-20 range and all other residues prefer the 4-10 range. The residues in all-beta proteins have an average of 3-8 long-range contacts whereas the residues in other classes have 1-4 long-range contracts. Furthermore, the preference of residue pairs to the folding and stability will be discussed.     

Influence of medium- and long-range interactions in different folding types of globular proteins

Recognition of protein fold from amino acid sequence is a challenging task. The structure and stability of proteins from different fold are mainly dictated by inter-residue interactions. In our earlier work, we have successfully used the medium- and long-range contacts for predicting the protein folding rates, discriminating globular and membrane proteins and for distinguishing protein structural classes. In this work, we analyze the role of inter-residue interactions in commonly occurring folds of globular proteins in order to understand their folding mechanisms. In the medium-range contacts, the globin fold and four-helical bundle proteins have more contacts than that of DNA-RNA fold although they all belong to all-alpha class. In long-range contacts, only the ribonuclease fold prefers 4-10 range and the other folding types prefer the range 21-30 in alpha/beta class proteins. Further, the preferred residues and residue pairs influenced by these different folds are discussed. The information about the preference of medium- and long-range contacts exhibited by the 20 amino acid residues can be effectively used to predict the folding type of each protein.     

Sequence-specific 1H NMR assignments and determination of the secondary structure for the activation domain isolated from pancreatic procarboxypeptidase B

Nearly complete sequence-specific 1H NMR assignments are presented for amino acid residues 3-81 in the 81-residue globular activation domain of porcine pancreatic procarboxypeptidase B isolated after limited tryptic proteolysis of the zymogen. These resonance assignments are consistent with the chemically determined amino acid sequence. Regular secondary structure elements were identified from nuclear Overhauser effects and the sequence locations of slowly exchanging backbone amide protons. The molecule contains two alpha-helices, including residues 20-30 and approximately residues 58-72, and a three-stranded antiparallel beta-sheet with the individual strands extending approximately from 12 to 17, 50 to 55, and 75 to 77. The identification of these secondary structures and a preliminary analysis of additional long-range NOE distance constraints show that isolated activation domain B forms a stable structure with the typical traits of a globular protein. The data presented here are the basis for the determination of the complete three-dimensional structure of activation domain B, which is currently in progress.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

Distribution of amino acid residues in the primary protein structure

The amino acid composition and sequence in primary structure of 180 proteins have been studied. It is shown that the distribution of amino acid residues is near to a random one, i.e. it is determined by the amino acid composition. The ratio between statistical and unique character of protein primary structures has been discussed. The amino acid sequence is suggested to be unique in fibrous proteins. In contrast the amino acid sequence in globular proteins is a statistical one. The statistical character of amino acids distribution in globular proteins explains the possibility of sensible text generation under the frame shift mutations, deletions and insertions.     

Protein as an edited random copolymer

It is widely believed that the unique primary structure of a given protein is quite necessary for its folding into a certain three-dimensional structure as well as for its functioning and is a result of a directed selection in the course of biological evolution. The present paper provides arguments in favour of an alternative point of view according to which typical three-dimensional structures of globular proteins are characteristic even for random sequences of amino acid residues. Therefore it may be possible that primary structures of proteins are mainly examples of random amino acid sequences slightly edited in the course of biological evolution to impart them some additional (functional) meaning.     

Sequence of contactin, a 130-kD glycoprotein concentrated in areas of interneuronal contact, defines a new member of the immunoglobulin supergene family in the nervous system

The primary amino acid sequence of contactin, a neuronal cell surface glycoprotein of 130 kD that is isolated in association with components of the cytoskeleton (Ranscht, B., D. J. Moss, and C. Thomas. 1984. J. Cell Biol. 99:1803-1813), was deduced from the nucleotide sequence of cDNA clones and is reported here. The cDNA sequence contains an open reading frame for a 1,071-amino acid transmembrane protein with 962 extracellular and 89 cytoplasmic amino acids. In its extracellular portion, the polypeptide features six type 1 and two type 2 repeats. The six amino-terminal type 1 repeats (I-VI) each consist of 81-99 amino acids and contain two cysteine residues that are in the right context to form globular domains as described for molecules with immunoglobulin structure. Within the proposed globular region, contactin shares 31% identical amino acids with the neural cell adhesion molecule NCAM. The two type 2 repeats (I-II) are each composed of 100 amino acids and lack cysteine residues. They are 20-31% identical to fibronectin type III repeats. Both the structural similarity of contactin to molecules of the immunoglobulin supergene family, in particular the amino acid sequence resemblance to NCAM, and its relationship to fibronectin indicate that contactin could be involved in some aspect of cellular adhesion. This suggestion is further strengthened by its localization in neuropil containing axon fascicles and synapses.     

Sequential resonance assignments in protein 1H nuclear magnetic resonance spectra: Computation of sterically allowed proton-proton distances and statistical analysis of proton-proton distances in single crystal protein conformations

Two different, theoretical studies of intramolecular proton-proton distances in polypeptide chains are described. Firstly, the distances between amide, C^α and C^β protons of neighbouring residues in the amino acid sequence, which correspond to the sterically allowed values for the dihedral angles φ_i, ψ_i and χ_i¹, were computed. Secondly, the frequency with which short distances occur between amide, C^α and C^β protons of neighbouring and distant residues in the amino acid sequence were statistically evaluated in a representative sample of globular protein crystal structures. Both approaches imply that semi-quantitative measurements of short, non-bonding proton-proton distances, e.g. by nuclear Overhauser experiments, should present a reliable and generally applicable method for sequential, individual resonance assignments in protein ¹H nuclear magnetic resonance spectra. Similar calculations imply that corresponding distance measurements can be used for resonance assignments in the side-chains of the aromatic amino acid residues, asparagine and glutamine, where the complete spin systems cannot usually be identified from through-bond spin-spin coupling connectivities.     

Sequence analysis of alpha 1(VI) and alpha 2(VI) chains of human type VI collagen reveals internal triplication of globular domains similar to the A domains of von Willebrand factor and two alpha 2(VI) chain variants that differ in the carboxy terminus.

Amino acid sequences of human collagen alpha 1(VI) and alpha 2(VI) chains were completed by cDNA sequencing and Edman degradation demonstrating that the mature polypeptides contain 1009 and 998 amino acid residues respectively. In addition, they contain small signal peptide sequences. Both chains show 31% identity in the N-terminal (approximately 235 residues) and C-terminal (approximately 430 residues) globular domains which are connected by a triple helical segment (335-336 residues). Internal alignment of the globular sequences indicates a repetitive 200-residue structure (15-23% identity) occurring three times (N1, C1, C2) in each chain. These repeating subdomains are connected to each other and to the triple helix by short (15-30 residues) cysteine-rich segments. The globular domains possess several N-glycosylation sites but no cell-binding RGD sequences, which are exclusively found in the triple helical segment. Sequencing of alpha 2(VI) cDNA clones revealed two variant chains with a distinct C2 subdomain and 3' non-coding region. The repetitive segments C1, C2 and, to a lesser extent, N1 show significant identity (15-18%) to the collagen-binding A domains of von Willebrand factor (vWF) and they are also similar to some integrin receptors, complement components and a cartilage matrix protein. Since the globular domains of collagen VI come into close contact with triple helical segments during the formation of tissue microfibrils it suggests that the globular domains bind to collagenous structures in a manner similar to the binding of vWF to collagen I.     

Solution structure and backbone dynamics of Mason-Pfizer monkey virus (MPMV) nucleocapsid protein.

Retroviral nucleocapsid proteins (NCPs) are CCHC-type zinc finger proteins that mediate virion RNA binding activities associated with retrovirus assembly and genomic RNA encapsidation. Mason-Pfizer monkey virus (MPMV), a type D retrovirus, encodes a 96-amino acid nucleocapsid protein, which contains two Cys-X2-Cys-X4-His-X4-Cys (CCHC) zinc fingers connected by an unusually long 15-amino acid linker. Homonuclear, two-dimensional sensitivity-enhanced 15N-1H, three-dimensional 15N-1H, and triple resonance NMR spectroscopy have been used to determine the solution structure and residue-specific backbone dynamics of the structured core domain of MPMV NCP containing residues 21-80. Structure calculations and spectral density mapping of N-H bond vector mobility reveal that MPMV NCP 21-80 is best described as two independently folded, rotationally uncorrelated globular domains connected by a seven-residue flexible linker consisting of residues 42-48. The N-terminal CCHC zinc finger domain (residues 24-37) appears to adopt a fold like that described previously for HIV-1 NCP; however, residues within this domain and the immediately adjacent linker region (residues 38-41) are characterized by extensive conformational averaging on the micros-ms time scale at 25 degrees C. In contrast to other NCPs, residues 49-77, which includes the C-terminal CCHC zinc-finger (residues 53-66), comprise a well-folded globular domain with the Val49-Pro-Gly-Leu52 sequence and C-terminal tail residues 67-77 characterized by amide proton exchange properties and 15N R1, R2, and (1H-15N) NOE values indistinguishable to residues in the core C-terminal finger. Twelve refined structural models of MPMV NCP residues 49-80 (pairwise backbone RMSD of 0.77 A) reveal that the side chains of the conserved Pro50 and Trp62 are in van der Waals contact with one another. Residues 70-73 in the C-terminal tail adopt a reverse turn-like structure. Ile77 is involved in extensive van der Waals contact with the core finger domain, while the side chains of Ser68 and Asn75 appear to form hydrogen bonds that stabilize the overall fold of this domain. These residues outside of the core finger structure are conserved in D-type and related retroviral NCPs, e.g., MMTV NCP, suggesting that the structure of MPMV NCP may be representative of this subclass of retroviral NCPs.     

NMR structure of the human doppel protein

The NMR structure of the recombinant human doppel protein, hDpl(24-152), contains a flexibly disordered "tail" comprising residues 24-51, and a globular domain extending from residues 52 to 149 for which a detailed structure was obtained. The globular domain contains four alpha-helices comprising residues 72-80 (alpha1), 101-115 (alpha2(a)), 117-121 (alpha2(b)), and 127-141 (alpha3), and a short two-stranded anti-parallel beta-sheet comprising residues 58-60 (beta1) and 88-90 (beta2). The fold of the hDpl globular domain thus coincides nearly identically with the structure of the murine Dpl protein. There are close similarities with the human prion protein (hPrP) but, similar to the situation with the corresponding murine proteins, hDpl shows marked local differences when compared to hPrP: the beta-sheet is flipped by 180 degrees with respect to the molecular scaffold formed by the four helices, and the beta1-strand is shifted by two residues toward the C terminus. A large solvent-accessible hydrophobic cleft is formed on the protein surface between beta2 and alpha3, which has no counterpart in hPrP. The helix alpha2 of hPrP is replaced by two shorter helices, alpha2(a) and alpha2(b). The helix alpha3 is shortened by more than two turns when compared with alpha3 of hPrP, which is enforced by the positioning of the second disulfide bond in hDpl. The C-terminal peptide segment 144-149 folds back onto the loop connecting beta2 and alpha2. All but four of the 20 conserved residues in the globular domains of hPrP and hDpl appear to have a structural role in maintaining a PrP-type fold. The conservation of R76, E96, N110 and R134 in hDpl, corresponding to R148, E168, N183 and R208 in hPrP suggests that these amino acid residues might have essential roles in the so far unknown functions of PrP and Dpl in healthy organisms.     

To be folded or to be unfolded?

The lack of ordered structure in “natively unfolded” proteins raises a general question: Are there intrinsic properties of amino acid residues that are responsible for the absence of fixed structure at physiological conditions? In this article, we demonstrate that the competence of a protein to be folded or to be unfolded may be determined by the property of amino acid residues to form a sufficient number of contacts in a globular state. The expected average number of contacts per residue calculated from the amino acid sequence alone (using the average number of contacts for 20 amino acid residues in globular proteins) can be used as one of the simple indicators of natively unfolded proteins. The prediction accuracy for the sets of 80 folded and 90 natively unfolded proteins reaches 89% if the expected average number of contacts is used as a parameter and 83% in the case of hydrophobicity. An optimal set of artificial parameters for 20 amino acid residues obtained by Monte Carlo algorithm to maximally separate the sets of 90 natively unfolded and 80 folded proteins demonstrates the upper limit for prediction accuracy, which is 95%.     

Spatial sign-alternating charge clusters in globular proteins

Large sign-alternating charge clusters formed by the charged side groups of amino acid residues and N- and C-terminal groups were found in the majority of considered globular proteins, namely 235 in a total of 274 protein structures, i.e. 85.8%. The clusters were determined by the criteria proposed earlier: charged groups were included in the cluster if their charged N and O atoms were located at distances between 2.4 and 7.0 A. The set of selected proteins consisted of known non-homologous protein structures from the Protein Data Bank with a resolution less than or equal to 2.5 A and pair sequence similarity less than 25%. Molecular masses of the proteins were from 5.5 to 91.5 kDa and protein chain length from 50 to 830 residues. The distribution of charged groups on the protein surface between isolated charged groups, small clusters with two and three groups, and large clusters with four or more groups were found to be approximately similar making 33, 35 and 32% of the total amount of protein charged groups, respectively. The large sign-alternating charge clusters with four or more charged groups were studied in greater detail. The amount of such clusters depends on the protein chain length. The small proteins contain 1-3 clusters while the large proteins display 4-6 or more clusters. On average, 1.5 clusters per each 100 residues were observed. In contrast with this, the size of a cluster, i.e. the number of charged groups inside a cluster, does not depend on the protein molecular mass, and large clusters are observed for proteins from a range of molecular masses. Clusters consisting of four to six charged groups occur most frequently, although extra large clusters are also often revealed. We can conclude that sign-alternating charge clusters are a common feature of the protein surface of globular protein. They are suggested to play a general functional role as a local polar factor of protein surface.     

Novel approach for alpha-helical topology prediction in globular proteins: generation of interhelical restraints

The protein folding problem represents one of the most challenging problems in computational biology. Distance constraints and topology predictions can be highly useful for the folding problem in reducing the conformational space that must be searched by deterministic algorithms to find a protein structure of minimum conformational energy. We present a novel optimization framework for predicting topological contacts and generating interhelical distance restraints between hydrophobic residues in alpha-helical globular proteins. It should be emphasized that since the model does not make assumptions about the form of the helices, it is applicable to all alpha-helical proteins, including helices with kinks and irregular helices. This model aims at enhancing the ASTRO-FOLD protein folding approach of Klepeis and Floudas (Journal of Computational Chemistry 2003;24:191-208), which finds the structure of global minimum conformational energy via a constrained nonlinear optimization problem. The proposed topology prediction model was evaluated on 26 alpha-helical proteins ranging from 2 to 8 helices and 35 to 159 residues, and the best identified average interhelical distances corresponding to the predicted contacts fell below 11 A in all 26 of these systems. Given the positive results of applying the model to several protein systems, the importance of interhelical hydrophobic-to-hydrophobic contacts in determining the folding of alpha-helical globular proteins is highlighted.     

Amino acid sequence of the non-collagenous globular domain (NC1) of the alpha 1(IV) chain of basement membrane collagen as derived from complementary DNA

NC1, the C-terminal non-collagenous globular domain of collagen IV, represents one of the two end regions responsible for the assembly and cross-linking of the extracellular network of basement membrane collagen. Several cDNA clones for the NC1 domain of the alpha 1(IV) collagen chain of mouse have been isolated by using synthetic oligonucleotides as screening probes for mouse libraries. The oligonucleotides were synthesized according to known stretches of the corresponding protein sequence. Sequencing of the overlapping cDNA clones allowed the complete amino acid sequence of the NC1 domain to be deduced as well as the C-terminal 165 amino acid residues of the triple helix. It consists of 229 amino acid residues which comprise two homologous regions with a high content of cysteine. These DNA and protein sequences are compared to the corresponding sequences of other collagens and discussed with respect to their structural and biological significance.     

Two-dimensional structure of beta-amyloid(10-35) fibrils

Beta-amyloid (Abeta) peptides are the main protein component of the pathognomonic plaques found in the brains of patients with Alzheimer's disease. These heterogeneous peptides adopt a highly organized fibril structure both in vivo and in vitro. Here we use solid-state NMR on stable, homogeneous fibrils of Abeta(10-35). Specific interpeptide distance constraints are determined with dipolar recoupling NMR on fibrils prepared from a series of singly labeled peptides containing (13)C-carbonyl-enriched amino acids, and skipping no more that three residues in the sequence. From these studies, we demonstrate that the peptide adopts the structure of an extended parallel beta-sheet in-register at pH 7.4. Analysis of DRAWS data indicates interstrand distances of 5.3 +/- 0.3 A (mean +/- standard deviation) throughout the entire length of the peptide, which is compatible only with a parallel beta-strand in-register. Intrastrand NMR constraints, obtained from peptides containing labels at two adjacent amino acids, confirm the secondary structural findings obtained using DRAWS. Using peptides with (13)C incorporated at the carbonyl position of adjacent amino acids, structural transitions from alpha-helix to beta-sheet were observed at residues 19 and 20, but using similar techniques, no evidence for a turn could be found in the putative turn region comprising residues 25-29. Implications of this extended parallel organization for Abeta(10-35) for overall fibril formation, stability, and morphology based upon specific amino acid contacts are discussed.     

The globular domains of type VI collagen are related to the collagen-binding domains of cartilage matrix protein and von Willebrand factor.

Type VI collagen is a transformation-sensitive glycoprotein of the extracellular matrix of fibroblasts. We have isolated and sequenced several overlapping cDNA clones (4153 bp) which encode the entire alpha 2 subunit of chicken type VI collagen. The deduced amino acid sequence predicts that the alpha 2(VI) polypeptide consists of 1015 amino acid residues that are arranged in four domains: a hydrophobic signal peptide of 20 residues, an amino-terminal globular domain of 228 residues, a collagenous segment of 335 residues and a carboxy-terminal globular domain of 432 residues. The collagenous domain contains seven Arg-Gly-Asp tripeptide units, some of which are likely to be used as cell-binding sites. The globular domains contain three homologous repeats with an average length of 180 amino acid residues. These repeats show a striking similarity to the collagen-binding motifs found in von Willebrand factor and cartilage matrix protein. We therefore speculate that the globular domains of the alpha 2(VI) polypeptide may interact with collagenous structures.     

Structural characterization of the trypsin-resistant core in the nuclear sperm-specific protein from Spisula solidissima

Trypsin digestion of the protamine-like protein from Spisula solidissima has revealed the existence of an internal resistant core. The peptide contains 75 amino acid residues, and its primary structure shows some conserved sequences that are common to those found in the core of the somatic histone H5 from chicken erythrocytes. The secondary structure of this core exhibits 33% antiparallel beta-sheet, 18% beta-turns, 37% random coil, and only 10% alpha-helix, in contrast to histone H5. Hydrodynamic measurements indicate a compact globular assembly for the tertiary structure of this peptide, when compared to the more extended shape observed for the whole protein. The possible relatedness of this protein to the histone H1 family is discussed.