Second codon positions of genes and the secondary structures of proteins. Relationships and implications for the origin of the genetic code

The nucleotide frequencies in the second codon positions of genes are remarkably different for the coding regions that correspond to different secondary structures in the encoded proteins, namely, helix, beta-strand and aperiodic structures. Indeed, hydrophobic and hydrophilic amino acids are encoded by codons having U or A, respectively, in their second position. Moreover, the beta-strand structure is strongly hydrophobic, while aperiodic structures contain more hydrophilic amino acids. The relationship between nucleotide frequencies and protein secondary structures is associated not only with the physico-chemical properties of these structures but also with the organisation of the genetic code. In fact, this organisation seems to have evolved so as to preserve the secondary structures of proteins by preventing deleterious amino acid substitutions that could modify the physico-chemical properties required for an optimal structure.     

Role of mutational bias and natural selection on genome-wide nucleotide bias in prokaryotic organisms

Correlations between genomic GC contents and amino acid frequencies were studied in the homologous sequences of 12 eubacterial genomes. Results show that amino acids encoded by GC-rich codons increases significantly with genomic GC contents, whereas opposite trend was observed in case of amino acids encoded by GC-poor codons. Further studies show all the amino acids do not change in the predicted direction according to their genomic GC pressure, suggesting that protein evolution is not entirely dictated by their nucleotide frequencies. Amino acid substitution matrix calculated among hydrophobic, amphipathic and hydrophilic amino acid groups' shows that amphipathic and hydrophilic amino acids are more frequently substituted by hydrophobic amino acids than from hydrophobic to hydrophilic or amphipathic amino acids. This indicates that nucleotide bias induces a directional changes in proteome composition in such a way that underwent strong changes in hydropathy values. In fact, significant increases in hydrophobicity values have also been observed with the increase of genomic GC contents. Correlations between GC contents and amino acid compositions in three different predicted protein secondary structures show that hydropathy values increases significantly with GC contents in aperiodic and helix structures whereas strand structure remains insensitive with the genomic GC levels. The relative importance of mutation and selection on the evolution of proteins have been discussed on the basis of these results.     

The formation of protein secondary structure. Its connection with amino acid sequence

Statistical analysis of the occurrence of tetrapeptides in 35 globular proteins was performed. It was found that the amino acids along the polypeptide chain are close to being randomly distributed and that the same tetrapeptide segments exist in different types of secondary structure. Therefore, a new method was proposed for locating 'microdomains' in protein interiors. Amino acid replacements in the hydrophobic core of six proteins were analyzed. The results show that the locations of amino acids belonging to defined microdomains are extremely conserved. It is suggested that the structures found may play a role as nucleation centers in protein folding.     

Conformational analysis of apolipoprotein A-I and E-3 based on primary sequence and circular dichroism.

The primary and secondary structure of human plasma apolipoprotein A-I and apolipoprotein E-3 have been analyzed to further our understanding of the secondary and tertiary conformation of these proteins and the structure and function of plasma lipoprotein particles. The methods used to analyze the primary sequence of these proteins used computer programs: (a) to identify repeated patterns within these proteins on the basis of conservative substitutions and similarities within the physicochemical properties of each residue; (b) for local averaging, hydrophobic moment, and Fourier analysis of the physicochemical properties; and (c) for secondary structure prediction of each protein carried out using homology, statistical, and information theory based methods. Circular dichroism was used to study purified lipid-protein complexes of each protein and quantitate the secondary structure in a lipid environment. The data from these analyses were integrated into a single secondary structure prediction to derive a model of each protein. The sequence homology within apolipoproteins A-I, E-3, and A-IV is used to derive a consensus sequence for two 11 amino acid repeating sequences in this family of proteins.     

A two-component nucleation model of protein hydrophobicity

A fundamental characteristic of soluble globular protein structure is a hydrophobic core and protein exterior comprised predominantly of hydrophilic residues. This distribution of amino acid residue hydrophobicity, from protein interior to exterior, has recently been profiled with the use of hydrophobic moments. The calculations enable comparison of the radial hydrophobicity distribution of different proteins and had revealed two features common to 30 proteins of diverse size and structure. One, a global feature, is the overall shape of the second-order ellipsoidal hydrophobic moment. The second, a specific feature, is a quasi-invariant hydrophobic-ratio of distances. Both features are dependent upon the rates of increase, from protein interior to exterior, of the accumulated numbers of hydrophobic and hydrophilic amino acid residues. These rates can be simulated simply with a two-component nucleation model of protein hydrophobicity. The model provides insight into the origin of the shape of the observed hydrophobic moment profiles and of the observed range of hydrophobic ratios. Consistent with observation, it is shown that a relatively wide range of hydrophobic and hydrophilic rates of increase yield a relatively narrow range of hydrophobic ratios. Furthermore, the model identifies one factor, the decrease in residue density with increasing distance from the protein interior, that is critical in providing the range of values that is comparable with the observed range.     

Anatomy of the herpes simplex virus 1 strain F glycoprotein B gene: primary sequence and predicted protein structure of the wild type and of monoclonal antibody-resistant mutants.

In this paper we report the nucleotide sequence and predicted amino acid sequence of glycoprotein B of herpes simplex virus 1 strain F and the amino acid substitutions in the domains of the glycoprotein B gene of three mutants selected for resistance to monoclonal antibody H126-5 or H233 but not to both. Analyses of the amino acid sequence with respect to hydropathicity and secondary structure yielded a two-dimensional model of the protein. The model predicts an N-terminal, 29-amino-acid cleavable signal sequence, a 696-amino-acid hydrophilic surface domain containing six potential sites for N-linked glycosylation, a 69-amino-acid hydrophobic domain containing three segments traversing the membrane, and a charged 109-amino-acid domain projecting into the cytoplasm and previously shown to marker rescue glycoprotein B syn mutations. The nucleotide sequence of the mutant glycoprotein B DNA fragments previously shown to marker transfer or rescue the mutations revealed that the amino acid substitutions cluster in the hydrophilic surface domain between amino acids 273 and 305. Analyses of the secondary structure of these regions, coupled with the experimentally derived observation that the H126-5- and H233-antibody cognitive sites do not overlap, indicate the approximate locations of the epitopes of these neutralizing, surface-reacting, and immune-precipitating monoclonal antibodies. The predicted perturbations in the secondary structure introduced by the amino acid substitutions correlate with the extent of loss of reactivity with monoclonal antibodies in various immunoassays.     

Sequence periodicity and secondary structure propensity in model proteins

We explore the question of whether local effects (originating from the amino acids intrinsic secondary structure propensities) or nonlocal effects (reflecting the sequence of amino acids as a whole) play a larger role in determining the fold of globular proteins. Earlier circular dichroism studies have shown that the pattern of polar, non polar amino acids (nonlocal effect) dominates over the amino acid intrinsic propensity (local effect) in determining the secondary structure of oligomeric peptides. In this article, we present a coarse grained computational model that allows us to quantitatively estimate the role of local and nonlocal factors in determining both the secondary and tertiary structure of small, globular proteins. The amino acid intrinsic secondary structure propensity is modeled by a dihedral potential term. This dihedral potential is parametrized to match with experimental measurements of secondary structure propensity. Similarly, the magnitude of the attraction between hydrophobic residues is parametrized to match the experimental transfer free energies of hydrophobic amino acids. Under these parametrization conditions, we systematically explore the degree of frustration a given polar, non polar pattern can tolerate when the secondary structure intrinsic propensities are in opposition to it. When the parameters are in the biophysically relevant range, we observe that the fold of small, globular proteins is determined by the pattern of polar, non polar amino acids regardless of their instrinsic secondary structure propensities. Our simulations shed new light on previous observations that tertiary interactions are more influential in determining protein structure than secondary structure propensity. The fact that this can be inferred using a simple polymer model that lacks most of the biochemical details points to the fundamental importance of binary patterning in governing folding.     

Multi-State Proteins: Approach Allowing Experimental Determination of the Formation Order of Structure Elements in the Green Fluorescent Protein

The most complex problem in studying multi-state protein folding is the determination of the sequence of formation of protein intermediate states. A far more complex issue is to determine at what stages of protein folding its various parts (secondary structure elements) develop. The structure and properties of different intermediate states depend in particular on these parts. An experimental approach, named μ-analysis, which allows understanding the order of formation of structural elements upon folding of a multi-state protein was used in this study. In this approach the same elements of the protein secondary structure are “tested” by substitutions of single hydrophobic amino acids and by incorporation of cysteine bridges. Single substitutions of hydrophobic amino acids contribute to yielding information on the late stages of protein folding while incorporation of ss-bridges allows obtaining data on the initial stages of folding. As a result of such an μ-analysis, we have determined the order of formation of beta-hairpins upon folding of the green fluorescent protein.     

Sequence shuffle controls morphological consequences in a self-assembling tetrapeptide.

Peptide and protein self-assembly is a well-studied phenomenon in chemistry and biology, where nanoscopic building blocks exhibit rapid self-association to reveal supramolecular aggregates of defined structural features. These superstructures are stabilized by hydrophobic interactions, hydrogen bonding and a host of other noncovalent interactions. Thus, amino acid side chains in the primary structure hold importance in dictating secondary structures and preference for particular conformational signatures in peptide aggregates. This report describes contrasting nanoscale morphologies in antamanide-derived synthetic tetrapeptide mutants, which are composed by shuffling only two amino acids: phenylalanine and proline. Remarkable differences in ultrastructures in primary sequence-shuffled tetrapeptides suggest dissimilar aggregational pathways due to context-dependent location of proline and phenylalanine residues with respect to one another.     

A hydrophobic heptad repeat of the core protein of woodchuck hepatitis virus is required for capsid assembly.

The capsid particle of hepadnaviruses is assembled from its dimer precursors. However, the mechanism of the protein-protein interaction is still poorly understood. A small region in the capsid protein of woodchuck hepatitis virus (WHV) contains four hydrophobic residues, including leucine 101, leucine 108, valine 115, and phenylalanine 122, that are conserved and spaced every seventh residue in the primary sequence to form a hydrophobic heptad repeat (hhr). A hydrophobic force often plays an important role in the interaction of proteins. Therefore, to investigate the role of this region in capsid assembly, we individually changed the codons specifying these four hydrophobic amino acids to codons specifying alanine or proline. In addition, we examined the in vivo infectivity of a WHV genome bearing a naturally occurring single amino acid change (histidine 104-->proline) in the hhr region. The phenotype of each altered genome was determined in both eukaryotic and prokaryotic systems by a capsid protein assay and electron microscopic examination. We show that replacement of any one of the four hydrophobic residues with alanine did not prevent capsid assembly. However, assembled capsid particles were not detected if combinations of any two of the four residues were substituted with alanines or if the spacing of these four hydrophobic residues was changed. An individual introduction of a proline (which dramatically changes the secondary structure of proteins) into different positions of this small region also abolished capsid assembly in vitro or viral replication in vivo. These results suggested that the hhr region of the core protein of WHV was critical for capsid assembly.     

Nearest-neighbor classifier as a tool for classification of protein families

Knowledge about protein function is essential in understanding the biological processes. A specific class or family of protein shares common structural and chemical properties amongst its member sequences. The set of properties that display its unique characteristics for clearly classifying a protein sequence into its corresponding protein family needs to be studied. Our study of these important properties conducted on four major classes of proteins namely Globins, Homeoboxes, Heat Shock proteins (HSP) and Kinase have shown that frequency of twenty naturally occurring amino acids, hydrophobic content of protein, molecular weight of protein, isoelectric point of protein, secondary structure composition of amino acid residues as helices, coils and sheets and the composition of helices, coils and sheets in the secondary structure topology plays a significant role in correctly classifying the protein into its corresponding class or family as indicated by the overall efficiency of Nearest Neighbor Classifier as 84.92%.     

苏太猪ELOVL7基因的生物信息学分析

不同于人、鼠等物种ELOVL7基因的高相似度,不同品种的猪ELOVL7基因相似度较低。为了探究该基因的特性,本研究运用生物信息学的方法对苏太猪ELOVL7基因及其氨基酸序列的同源性、理化性质、保守结构域、亚细胞定位、信号肽、跨膜结构域、亲水性/疏水性、二级结构、功能预测以及磷酸化位点等进行预测分析。结果表明：在苏太猪中,ELOVL7全长2 324 bp,编码区为846 bp,共编码281个氨基酸。其结构稳定,分子量为33 387.4 Da,带正电荷,偏碱性。该基因所编码的蛋白质最可能位于细胞膜上,主要的功能是运输和结合,为跨膜、非分泌型疏水蛋白质,含有1个GNS1/SUR4家族的保守结构域,并有15个丝氨酸激酶、15个苏氨酸激酶和20个酪氨酸激酶潜在磷酸化位点。α螺旋是ELOVL7二级结构和三级结构中最主要的结构元件。另外ELOVL7与大部分物种的氨基酸序列相似性达90%以上,且亲缘关系较近。分析ELOVL7基因及其氨基酸序列的特征,能够为进一步挖掘该基因内的突变对长链脂肪酸表型的影响以及合成、代谢机理提供分子依据。     

Influence of the hydrophilic face on the folding ability and stability of alpha-helix bundles: relevance to the peptide catalytic activity.

Although not the sole feature responsible, the packing of amino acid side chains in the interior of proteins is known to contribute to protein conformational specificity. While a number of amphipathic peptide sequences with optimized hydrophobic domains has been designed to fold into a desired aggregation state, the contribution of the amino acids located on the hydrophilic side of such peptides to the final packing has not been investigated thoroughly. A set of self-aggregating 18-mer peptides designed previously to adopt a high level of alpha-helical conformation in benign buffer is used here to evaluate the effect of the nature of the amino acids located on the hydrophilic face on the packing of a four alpha-helical bundle. These peptides differ from one another by only one to four amino acid mutations on the hydrophilic face of the helix and share the same hydrophobic core. The secondary and tertiary structures in the presence or absence of denaturants were determined by circular dichroism in the far- and near-UV regions, fluorescence and nuclear magnetic resonance spectroscopy. Significant differences in folding ability, as well as chemical and thermal stabilities, were found between the peptides studied. In particular, surface salt bridges may form which would increase both the stability and extent of the tertiary structure of the peptides. The structural behavior of the peptides may be related to their ability to catalyze the decarboxylation of oxaloacetate, with peptides that have a well-defined tertiary structure acting as true catalysts.     

Structure of a 16 kDa integral membrane protein that has identity to the putative proton channel of the vacuolar H(+)-ATPase.

A 16 kDa protein has been isolated in a homogeneous form as the major component of a paracrystalline paired membrane structure closely resembling the gap junction. The primary structure of this protein from arthropod and vertebrate species has been determined by protein and cDNA sequencing. The amino acid sequences are highly conserved and virtually identical to the amino acid sequence of the proteolipid subunit of the vacuolar H(+)-ATPases. The disposition of the protein in the membrane has been studied using proteases and the N,N'-dicyclohexylcarbodiimide reactive site identified. These data, together with secondary structure predictions, suggest that the 16 kDa protein is for the most part buried in the membrane, arranged in a bundle of four hydrophobic alpha-helices. Using computer graphics, a model has been constructed based on this arrangement and on the electron microscopic images of the paracrystalline arrays.     

Proteins of the kidney microvillar membrane. The amphipathic form of dipeptidyl peptidase IV.

Dipeptidyl peptidase IV was solubilized from the microvillar membrane of pig kidney by Triton X-100. The purified enzyme was homogeneous on polyacrylamide-gel electrophoresis and ultracentrifugation, although immunoelectrophoresis indicated that amino-peptidase M was a minor contaminant. A comparison of the detergent-solubilized and proteinase (autolysis)-solubilized forms of the enzyme was undertaken to elucidate the structure and function of the hydrophobic domain that serves to anchor the protein to the membrane. No differences in catalytic properties, nor in sensitivity to inhibition by di-isopropyl phosphorofluoridate were found. On the other hand, several structural differences could be demonstrated. Both forms were about 130,000 subunit mol.wt., but the detergent form appeared to be larger by no more than about 4,000. Electron microscopy showed both forms to be dimers, and gel filtration revealed a difference in the dimeric mol.wt. of about 38 000, mainly attributable to detergent molecules bound to the hydrophobic domain. Papain converted the detergent form into a hydrophilic form that could not be distinguished in properties from the autolysis form. A hydrophobic peptide of about 3500 mol.wt. was identified as a product of papain treatment. The detergent and proteinase forms differed in primary structure. Partial N-terminal amino acid sequences were shown to be different, and the pattern of release of amino acids from the C-terminus by carboxypeptidase Y was essentially similar. The results are consistent with a model in which the protein is anchored to the microvillar membrane by a small hydrophobic domain located within the N-terminal amino acid sequence of the polypeptide chain. The significance of these results in relation to biosynthesis of the enzyme and assembly in the membrane is discussed.     

不同类型氨基酸网络参量与蛋白质折叠的关系

理论和实验研究表明,蛋白质天然拓扑结构对其折叠过程具有重要的影响．采用复杂网络的方法分析蛋白质天然结构的拓扑特征,并探索蛋白质结构特征与折叠速率之间的内在联系．分别构建了蛋白质氨基酸网络、疏水网、亲水网、亲水-疏水网以及相应的长程网络,研究了这些网络的匹配系数(assortativity coefficient)和聚集系数(clustering coefficient)的统计特性．结果表明,除了亲水-疏水网,上述各网络的匹配系数均为正值,并且氨基酸网和疏水网的匹配系数与折叠速率表现出明显的线性正相关,揭示了疏水残基间相互作用的协同性有助于蛋白质的快速折叠．同时,研究发现疏水网的聚集系数与折叠速率有明显的线性负相关关系,这表明疏水残基间三角结构(triangle construction)的形成不利于蛋白质快速折叠．还进一步构建了相应的长程网络,发现序列上间距较远的残基接触对的形成将使蛋白质折叠进程变慢．     

Predicted Folding of β-Structure in Myelin Basic Protein

Predictions of myelin basic protein secondary structure have not previously considered a major role for beta-structure in the organization of the native molecule because optical rotatory dispersion and circular dichroism studies have provided little, if any, evidence for beta-structure, and because a polycationic protein is generally considered to resist folding into a compact structure. However, the Chou-Fasman, Lim, and Robson algorithms identify a total of five beta-strands in the amino acid sequence. Four of these hydrophobic amino acid sequences (37-45, 87-95, 110-118, and 150-158) could form a hairpin intermediate that initiates folding of a Greek-key-type beta-structure. A second fold on the more hydrophobic side, with the addition of a strand from the N-terminus (residues 13-21), would complete the five-stranded antiparallel beta-sheet. A unique strand alignment can be predicted by phasing the hydrophobic residues. The unusual triproline sequence of myelin basic protein (100-102) is enclosed in the 14-residue hairpin loop. If these prolines are in the trans conformation, models show that a reverse turn could occur at residues 102-105 (Pro-Ser-Gln-Gly). Algorithms do not agree on the prediction of alpha-helices, but each of the two large loops could accommodate an alpha-helix. Myelin basic protein is known to be phosphorylated in vivo on as many as five Ser/Thr residues. Phosphorylation might alter the dynamics of folding if the nascent polypeptide were phosphorylated in the cytoplasm. In particular, phosphorylation of Thr-99 could neutralize cationic residues Lys-106 and Arg-108 within the hairpin loop. In addition, the methylation of Arg-108 might stabilize the hairpin loop structure through hydrophobic interaction with the side chain of Pro-97. The cationic side chains of arginine and lysine residues located on the faces of the beta-sheet (Arg-43, Arg-114, Lys-13, Lys-92, Lys-153, and Lys-156) could provide sites for interaction with phospholipids and other anionic structures on the surface of the myelin lipid bilayer.     

Hydrophobicity and structural classes in proteins.

The bulk hydrophobic character for the 20 natural amino acid residues, has been obtained from a database of 60 protein structures, grouped in the four structural classes alpha alpha, beta beta, alpha + beta and alpha/beta. The hydrophobicity coefficients thus obtained are compared with Ponnuswamy's original values using scales normalized to average = 0.0 and standard deviation = 1.0. Even though most of the amino acid residues do not change their hydropathic character in the different structural classes, their behaviour suggests the convenience that averaging methods should only consider proteins of the same structural class and that this information should be included in the secondary structure methods.     

Secondary structures of a new class of lipid body proteins from oilseeds.

The three main isoforms of the 19-kDa lipid body proteins (oleosin) have been purified to homogeneity from embryos of rapeseed. The secondary structures of these proteins as derived from circular dichroism (CD) and Fourier transform infrared (FTIR) spectroscopy were compared with the secondary structures predicted from the primary sequences. The salient feature of the primary sequence of all oleosins is its division into three defined structural domains: a central hydrophobic domain flanked on either side by relatively hydrophilic domains, respectively. Using a variety of predictive methods based on primary amino acid sequence data, the oleosins exhibited a high probability of beta-strand structure in the 70-residue central hydrophobic domain, with relatively little alpha-helical content. Secondary structure data derived from CD and FTIR were consistent with the predictions from primary sequence, showing that the oleosins contained about 45% beta-strand and 13% alpha-helical structure. Under high salt conditions, a 40-kDa polypeptide was obtained from purified preparations of the 19-kDa oleosins. The 40-kDa polypeptide has a very similar secondary structure, as analyzed by CD and FTIR, to that of the 19-kDa oleosins. This polypeptide is therefore probably a dimer of the 19-kDa oleosins that is formed in high salt environments. A model of the general structure of oleosins is proposed whereby the central hydrophobic domain of the protein with a predominantly beta-strand structure is embedded into the non-aqueous phase of lipid-bodies. This hydrophobic region is flanked by putative alpha-helical structures in the polar N- and C-terminal domains which are probably oriented at the lipid-water interface.