Dynamics of amino acid residues in globular proteins

The dynamic differential equation model developed and tested for bovine pancreatic trypsin inhibitor and tuna ferrocytochrome c in Ponnuswamy, P.K. & Bhaskaran, R. (Int. J. Peptide Protein Res. 24, 168-179, 1984) is extended for 17 more protein crystals in this work. Average displacements are computed for 20 amino acid residues observed in 19 proteins. Detailed information on the dynamic behaviour of the individual proteins and individual residues is presented. The effect of atomic packing on the fluctuations of the amino acid residues in alpha-chymotrypsin is illustrated. A number of general points on the dynamic characteristics of globular protein molecules are presented.     

Amino acid sequence of bovine neurophysin-II: a reinvestigation.

The entire amino acid sequence of bovine neurophysin-II has been redetermined by manual Edman degradation of tryptic peptides obtained from performic acid-oxidized neurophysin. Electrophoretically homogeneous bovine neurophysin-II was found to be a mixture of two species of protein molecules both containing 95 amino acid residues. The only difference between the two species of the neurophysin molecules is a single amino acid substitution at residue 89. Of the bovine neurophysin-II used in this work 70% of the protein material contained valine and 30% contained isoleucine at residue 89 in their sequences. The redetermined sequences of bovine neurophysin-II shown in Fig. 2 differ substantially from the reported sequence of bovine neurophysin-II but resemble closely that of porcine neurophysin-I and ovine neurophysin-III (Fig. 3).     

Protein packing: dependence on protein size, secondary structure and amino acid composition

We have used the occluded surface algorithm to estimate the packing of both buried and exposed amino acid residues in protein structures. This method works equally well for buried residues and solvent-exposed residues in contrast to the commonly used Voronoi method that works directly only on buried residues. The atomic packing of individual globular proteins may vary significantly from the average packing of a large data set of globular proteins. Here, we demonstrate that these variations in protein packing are due to a complex combination of protein size, secondary structure composition and amino acid composition. Differences in protein packing are conserved in protein families of similar structure despite significant sequence differences. This conclusion indicates that quality assessments of packing in protein structures should include a consideration of various parameters including the packing of known homologous proteins. Also, modeling of protein structures based on homologous templates should take into account the packing of the template protein structure.     

N-Glycans protect proteins from protease digestion through their binding affinities for aromatic amino acid residues

It was previously revealed [Yamaguchi, H., Nishiyama, T., and Uchida, M. (1999) J. Biochem. 126, 261-265] that N-glycans of both the high-mannose and complex types have binding affinity for aromatic amino acid residues. This study shows that free N-glycans protect proteins from protease digestion through their binding affinities for the aromatic amino acid residues exposed on protein molecules. Protease digestion of bovine pancreatic RNase A and bovine a-lactalbumin was depressed in solutions (1 mM or so) of free N-glycans of both the high-mannose and complex types. The increasing order of the protective effects of the N-glycans paralleled that of their affinities for aromatic amino acid residues; and the presence of aromatic amino acids practically abolished the protective effects of the N-glycans. The N-glycans also depressed the protease digestion of metallothionein, an aromatic amino acid-free protein, in agreement with the observation that the N-glycans also interact with the solvent-exposed aromatic amino acid residues of the proteases. Thus it seems probable that the N-glycans protect proteins from protease digestion by steric hindrance attributable to their binding affinity for the solvent-exposed aromatic amino acid residues of both substrate proteins and proteases.     

Protein Structures and Neutral Theory of Evolution

Abstract

The neutral theory of evolution is extended to the origin of protein molecules. Arguments are presented which suggest that the amino acid sequences of many globular proteins mainly represent “memorized” random sequences while biological evolution reduces to the “editing” these random sequences. Physical requirements for a functional globular protein are formulated and it is shown that many of these requirements do not involve strategical selection of amino acid sequences during biological evolution but are inherent also for typical random sequences. In particular, it is shown that random sequences of polar and unpolar amino acid residues can form α-helices and β-strands with lengths and arrangement along the chain similar to those in real globular proteins. These α- and β-regions in random sequences can form three-dimensional folding patterns also similar to those in real proteins. The arguments are presented suggesting that even the tight packing of side groups inside protein core do not require very strong biological selection of amino acid sequences either. Thus many structural features of real proteins can exist also in random sequences and the biological selection is needed mainly for the creation of active sites of proteins and for their stability under physiological conditions.     

Method of Determining Protein Conformations by the Two-Dimensional Nuclear Overhauser Enhancement Spectroscopy Data

Abstract

A method is suggested to determine the most probable values of the angles Φ, Ψ of the protein backbone by the data on the availability and absence of d connectivities in the two-dimensional nuclear Overhauser enhancement spectra. In view of this, the dependences of the protonproton distances in dipeptide units of L-amino acid residues on the dihedral angles Φ, Ψ, χ, are considered and the conformational states of amino acid residues of the proteins with the known spatial structure are analysed statistically. The potentialities of the method are assessed with the aid of model spectral nuclear magnetic resonance (NMR) parameters obtained from the X-ray data for the bovine pancreatic trypsin inhibitor and avian pancreatic polypeptide.

It is shown that the developed procedure of structural interpretation of the NMR data allows one to correctly reproduce the local conformation of the protein backbone. The obtained backbone conformation may serve as a starting point to build and refine molecular three- dimensional structure.     

Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry

Globular proteins have high packing densities as a result of residue side chains in the core achieving a tight, complementary packing. The internal packing is considered the main determinant of native protein structure. From that point of view, we present here a method of energy minimization using an automata network to predict a set of amino acid sequences and their side-chain conformations from a desired backbone geometry for de novo design of proteins. Using discrete side-chain conformations, that is, rotamers, the sequence generation problem from a given backbone geometry becomes one of combinatorial problems. We focused on the residues composing the interior core region and predicted a set of amino acid Sequences and their side-chain conformations only from a given backbone geometry. The kinds of residues were restricted to six hydrophobic amino acids (Ala, Ile, Met, Leu, Phe, and Val) because the core regions are almost always composed of hydrophobic residues. The obtained sequences were well packed as was the native sequence. The method can be used for automated sequence generation in the de novo design of proteins. © 1994 Wiley-Liss, Inc.     

Determination and comparative analysis of the conformation of bovine pancreatic trypsin inhibitor and trypsin inhibitors E and K from the data of two-dimensional 1H-NMR spectroscopy

On the basis of joint consideration of distance dependences between amide proton NH and protons C alpha H, NH, C beta H of the preceding in amino acid sequence residue from the torsion angles phi psi, chi 1, the correlation diagram of these proton-proton distances with the regions of sterically allowed conformational space (phi, psi) is presented and the method for the determination of the L-amino acid residues backbone conformations is proposed. The diagram was used for the determination of backbone conformations of bovine pancreatic trypsin inhibitor and trypsin inhibitors E and K from Dendroaspis polylepis using the data from two-dimensional 1H-NMR spectroscopy. The analysis of backbone conformations was carried out. The individual elements of these protein molecules secondary structure were characterized and their high conformational homology was shown. The inference about qualitative coincidence of three protein molecules conformation in solution, preservation of secondary structure basic elements and their similarity with bovine pancreatic trypsin inhibitor crystalline structure was made.     

The spatial distribution of physical, chemical, energetic and conformational properties of amino acid residues in globular proteins

An attempt is made to obtain information as to the spatial distribution of amino acid residues in globular proteins in terms of their chemical, physical, energetic and conformational properties. The crystallographic data on twenty-one protein molecules form the basis for the study. The properties of the residues, namely, hydrophobicity, polarity, acidity, molecular weight, bulkness, chromatographic index, refractivity, short/medium/long-range energetics, and powers to adopt a-helical, extended and bend structures are analysed by dividing the protein globule into six concentric shells containing equal numbers of residues. The results show that the decisive factor in determining the spatial position of a residue in a protein molecule is of composite nature involving a compromise between the various properties of the residue. The observed deviations from the general hydrophobic interior and hydrophilic exterior in globular proteins are nicely brought out. A valuable clue is obtained as to the directional properties of α-helical and extended structure segments in proteins. The relative buried, exposed and intermediate characters of the residues in protein globules are obtained in a very realistic approach.     

Distributions of water around amino acid residues in proteins

The atomic co-ordinates from 16 high-resolution (less than or equal to 1.7 A = 0.1 nm), non-homologous proteins have been used to study the distributions of water molecule sites around the 20 different amino acid residues. The proportion of residues whose main-chain atoms are in contact with water molecules was fairly constant (between 40% and 60%), irrespective of the nature of the side-chain. However, the proportion of residues whose side-chain atoms were in contact with water molecules showed a clear (inverse) correlation with the hydrophobicity of the residue, being as low as 14% for leucine and isoleucine but greater than 80% for asparagine and arginine. Despite the problems in determining accurate water molecule sites from X-ray diffraction data and the complexity of the protein surface, distinct non-random distributions of water molecules were found. These hydration patterns are consistent with the expected stereochemistry of the potential hydrogen-bonding sites on the polar side-chains. The water molecules around apolar side-chains lie predominantly at van der Waals' contact distances, but most of these have a primary, shorter contact with a neighbouring polar atom. Further analysis of these distributions, combined with energy minimization techniques, should lead to improved modelling of protein structures, including their primary shells of hydration.     

Primary structure of bovine conglutinin, a member of the C-type animal lectin family

We report the complete amino acid sequence of bovine conglutinin obtained by structural characterization of peptides derived from the protein by various chemical and enzymatic fragmentation methods. The protein consists of 351 amino acid residues including 55 apparent Gly-X-Y repeats with two interruptions. This 171-residue-long collagenous domain separates a short noncollagenous NH2-terminal region of 25 residues from the 155-residue-long globular COOH terminus revealing the structural relation of conglutinin with mannose-binding proteins, pulmonary surfactant-associated proteins, and a complement component C1q. Eight hydroxylysine residues were found in the collagenous domain. All of these hydroxylysine residues which occupy a Y position in a Gly-X-Y triplet are possible glycosylation sites since no phenylthiohydantoin amino acid was identified in automated Edman degradation cycles corresponding to these sites. The noncollagenous COOH domain of conglutinin, on the other hand, contains a carbohydrate recognition domain which shares substantial sequence homology with C-type animal lectins. Conglutinin has the greatest sequence similarity with mannose-binding proteins and pulmonary surfactant-associated proteins.     

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

Characterization of Human Brain S100 Protein Fraction: Amino Acid Sequence of S100β

Two major components of human brain S100 fraction were purified by HPLC and an amino acid sequence was elucidated for the S100 beta component. Human S100 proteins showed absorption spectra and amino acid compositions similar to S100 alpha and S100 beta from bovine brain. However, the relative amounts of the human proteins were 4% S100 alpha and 96% S100 beta by weight, while the bovine protein distribution was 47% S100 alpha and 53% S100 beta by weight. An amino acid sequence of human S100 beta was established by analysis of overlapping fragments generated by cyanogen bromide and trypsin cleavage. Three amino acid sequence differences between the human and bovine S100 beta were found at residues 7, 62, and 80. These differences were chemically conservative and compatible with minimum single base changes in the codon structures. These results document that S100 beta is a conserved protein among mammals and provide the necessary foundation for current clinical studies.     

Nuclear basic protein transition during sperm differentiation. Amino acid sequence of a spermatid-specific protein from the dog-fish Scylliorhinus caniculus

During dog-fish spermatogenesis, chromatin undergoes a continuous processing which involves two basic protein transitions: the first from somatic-type histones to spermatid-specific proteins and the second leading to protamines. Two spermatid-specific proteins S1 and S2 were isolated from nuclei of spermatid-enriched testis zone and the amino acid sequence of S1 has been determined. S1 contains 87 amino acids and has a molecular mass of 11179 Da. It is mainly characterized by a high content of basic residues (45%) and the presence of one residue of cysteine. Its primary structure shows that the N-terminal half is highly basic while the hydrophobic residues are preferentially localized in the C-terminal region. Three forms of S1 are present in testis which correspond to di-, mono- and nonphosphorylated molecules. This spermatid-specific protein shares no common structural feature with either histones and dog-fish protamines or rat spermatid-specific protein which has been previously described.     

beta-Hydroxyaspartic acid or beta-hydroxyasparagine in bovine low density lipoprotein receptor and in bovine thrombomodulin

All of the vitamin K-dependent plasma proteins with domains that are homologous to the epidermal growth factor (EGF) precursor have 1 hydroxylated aspartic acid residue in the NH2-terminal EGF-homology region. In addition, protein S has 1 hydroxylated asparagine residue in each of the three COOH-terminal EGF-homology regions. All of these proteins have been found to have the amino acid sequence, CX(D or N)XXXX(F or Y)XCXC (corresponding to residues 20 to 33 in EGF), where the Asp or Asn residue is hydroxylated. This sequence also appears in two of the three EGF-homology regions of the human low density lipoprotein receptor and in two of the six EGF-homology regions of bovine thrombomodulin so far identified, suggesting that they may have the modified amino acid. We have now identified beta-hydroxyaspartic acid in acid hydrolysates of both these proteins.     

Mammalian hemodynamics: a new similarity principle

An attempt is made to obtain information as to the spatial distribution of amino acid residues in globular proteins in terms of their chemical, physical, energetic and conformational properties. The crystallographic data on twenty-one protein molecules form the basis for the study. The properties of the residues, namely, hydrophobicity, polarity, acidity, molecular weight, bulkness, chromatographic index, refractivity, short/medium/long-range energetics, and powers to adopt a-helical, extended and bend structures are analysed by dividing the protein globule into six concentric shells containing equal numbers of residues. The results show that the decisive factor in determining the spatial position of a residue in a protein molecule is of composite nature involving a compromise between the various properties of the residue. The observed deviations from the general hydrophobic interior and hydrophilic exterior in globular proteins are nicely brought out. A valuable clue is obtained as to the directional properties of α-helical and extended structure segments in proteins. The relative buried, exposed and intermediate characters of the residues in protein globules are obtained in a very realistic approach.     

A numerical measure of amino acid residues similarity based on the analysis of their surroundings in natural protein sequences

A measure of similarity between amino acid residues based on the analysis of the surroundings of each residue in primary structures of native proteins is proposed. The statistical data used for this purpose were obtained from the analysis of 168,808 protein sequences, which comprise the Protein Identification Research database (release 63). Using various threshold values of the proposed measure, amino acid residues were classified into several groups. The classification elaborated differs essentially from groupings previously used. The numerical measure of amino acid residues similarity can be used in site-directed mutagenesis studies for the prediction of probability of local spatial rearrangements in proteins.     

New amino acid residue type identification experiments valid for protonated and deuterated proteins

Two experiments are presented that yield amino acid type identification of individual residues in a protein by editing the ¹H?C¹⁵N correlations into four different 2D subspectra, each corresponding to a different amino acid type class, and that can be applied to deuterated proteins. One experiment provides information on the amino acid type of the residue preceding the detected amide ¹H?C¹⁵N correlation, while the other gives information on the type of its own residue. Versions for protonated proteins are also presented, and in this case it is possible to classify the residues into six different classes. Both sequential and intraresidue experiments provide highly complementary information, greatly facilitating the assignment of protein resonances. The experiments will also assist in transferring the assignment of a protein to the spectra obtained under different experimental conditions (e.g. temperature, pH, presence of ligands, cofactors, etc.).     

Amino acid sequence of P-57, a neurospecific calmodulin-binding protein

The amino acid sequence was determined for bovine brain P-57, a neurospecific, membrane-associated, calmodulin-binding protein. It consists of a single 239-residue polypeptide chain blocked at its amino terminus and containing an unusually hydrophilic amino acid composition. Seventy percent of the molecule is composed of Glu/Gln, Ala, Lys, Asp/Asn, and Pro; there is only one aromatic residue. A lack of favorable cleavage sites required that a particularly wide variety of digests and subdigests be performed to obtain appropriate sets of overlapping peptides. This protein is clearly homologous with the cDNA-derived sequence of mouse brain P-57, although the bovine protein is 12 amino acid residues longer; the homology is less obvious in the middle sections of the two sequences. Bovine brain P-57 lacks homology with any other protein in an updated sequence database. A segment reported to interact with calmodulin (Arg-Gly-His-Ile-Thr-Arg-Lys-Lys-Leu) is placed at residues 43-51 within the only extended segment of P-57 that carries the net positive charge that would favor that interaction. There is no hydrophobic segment characteristic of many proteins that interact with membranes.