Quaternion-based discriminant analysis method for color face recognition

Pattern recognition techniques have been used to automatically recognize the objects, personal identities, predict the function of protein, the category of the cancer, identify lesion, perform product inspection, and so on. In this paper we propose a novel quaternion-based discriminant method. This method represents and classifies color images in a simple and mathematically tractable way. The proposed method is suitable for a large variety of real-world applications such as color face recognition and classification of the ground target shown in multispectrum remote images. This method first uses the quaternion number to denote the pixel in the color image and exploits a quaternion vector to represent the color image. This method then uses the linear discriminant analysis algorithm to transform the quaternion vector into a lower-dimensional quaternion vector and classifies it in this space. The experimental results show that the proposed method can obtain a very high accuracy for color face recognition.     

Characterization of the backbone geometry of protein native state structures

We present a simple method for analyzing the geometry of noncovalent residue-residue interactions stabilizing the protein structure, which takes into account the constraints on the local backbone geometry. We find that the principal geometrical constraints are amino acid aspecific and are associated with hydrogen bond formation in helices and sheets. In contrast, amino acid residues in nonhelical and nonextended conformations, which make noncovalent interactions stabilizing the protein tertiary structure, display greater flexibility. We apply the method to an analysis of the packing of helices in helical bundle proteins requiring an efficient packing of amino acid side-chains of the interacting helices.     

Computation of 3D structure and distribution of helical parameters for D-period segments of human fibrillar collagen III

The structures of D-period segments of collagen (234 amino residues or ～1/4 of whole length) are established by methods of molecular mechanics and geometry analysis. Each D-period segment proves to have a unique spatial structure. The distributions of local helical parameters along the molecule are calculated. It is found that a second hydrogen bond is formed in every case when the second residue in the tripeptide G-X-Y is an amino acid. With such a combined H-bond network, all the peptide CO groups of glycines and of the third residues in tripeptides have quasi-equivalent positions on the surface of the collagen molecule. The local deformations of the polyproline II helix in the triple complex give rise to the observed modulation of structure at the macromolecular level, which may be important for the mutual orientation of collagen molecules during fibrillogenesis.     

Exploring the sequence patterns in the alpha-helices of proteins

This paper reports an extensive sequence analysis of the alpha-helices of proteins. alpha-Helices were extracted from the Protein Data Bank (PDB) and were divided into groups according to their sizes. It was found that some amino acids had differential propensity values for adopting helical conformation in short, medium and long alpha-helices. Pro and Trp had a significantly higher propensity for helical conformation in short helices than in medium and long helices. Trp was the strongest helix conformer in short helices. Sequence patterns favoring helical conformation were derived from a neighbor-dependent sequence analysis of proteins, which calculated the effect of neighboring amino acid type on the propensity of residues for adopting a particular secondary structure in proteins. This method produced an enhanced statistical significance scale that allowed us to explore the positional preference of amino acids for alpha-helical conformations. It was shown that the amino acid pair preference for alpha-helix had a unique pattern and this pattern was not always predictable by assuming proportional contributions from the individual propensity values of the amino acids. Our analysis also yielded a series of amino acid dyads that showed preference for alpha-helix conformation. The data presented in this study, along with our previous study on loop sequences of proteins, should prove useful for developing potential 'codes' for recognizing sequence patterns that are favorable for specific secondary structural elements in proteins.     

Insights from modeling protein evolution with context-dependent mutation and asymmetric amino acid selection

We develop an approximate maximum likelihood method to estimate flanking nucleotide context-dependent mutation rates and amino acid exchange-dependent selection in orthologous protein-coding sequences and use it to analyze genome-wide coding sequence alignments from mammals and yeast. Allowing context-dependent mutation provides a better fit to coding sequence data than simpler (context-independent or CpG "hotspot") models and significantly affects selection parameter estimates. Allowing asymmetric (nonreciprocal) selection on amino acid exchanges gives a better fit than simple dN/dS or symmetric selection models. Relative selection strength estimates from our models show good agreement with independent estimates derived from human disease-causing and engineered mutations. Selection strengths depend on local protein structure, showing expected biophysical trends in helical versus nonhelical regions and increased asymmetry on polar-hydrophobic exchanges with increased burial. The more stringent selection that has previously been observed for highly expressed proteins is primarily concentrated in buried regions, supporting the notion that such proteins are under stronger than average selection for stability. Our analyses indicate that a highly parameterized model of mutation and selection is computationally tractable and is a useful tool for exploring a variety of biological questions concerning protein and coding sequence evolution.     

Structure and function of Hib pili from Haemophilus influenzae type b

Pathogenic bacteria are specifically adapted to bind to their customary host. Disease is then caused by subsequent colonization and/or invasion of the local environmental niche. Initial binding of Haemophilus influenzae type b to the human nasopharynx is facilitated by Hib pili, filaments expressed on the bacterial surface. With three-dimensional reconstruction of electron micrograph images, we show that Hib pili comprise a helix 70 A in diameter with threefold symmetry. The Hib pilus filament has 3.0 subunits per turn, with each set of three subunits translated 26.9 A along and rotated 53 degrees about the helical axis. Amino acid sequence analysis of pilins from Hib pili and from P-pili expressed on uropathogenic Escherichia coli were used to predict the physical location of the highly variable and immunogenic region of the HifA pilin in the Hib pilus structure. Structural differences between Hib pili and P-pili suggest a difference in the strategies by which bacteria remain bound to their host cells: P-pili were shown to be capable of unwinding to five times their original length (E. Bullitt and L. Makowski, Nature 373:164-167, 1995), while damage to Hib pili occurs by slight shearing of subunits with respect to those further along the helical axis. This capacity to resist unwinding may be important for continued adherence of H. influenzae type b to the nasopharynx, where the three-stranded Hib pilus filaments provide a robust tether to withstand coughs and sneezes.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. I. Myohemerythrin.

In an attempt to delineate potential folding initiation sites for different protein structural motifs, we have synthesized series of peptides that span the entire length of the polypeptide chain of two proteins, and examined their conformational preferences in aqueous solution using proton nuclear magnetic resonance and circular dichroism spectroscopy. We describe here the behavior of peptides derived from a simple four-helix bundle protein, myohemerythrin. The peptides correspond to the sequences of the four long helices (the A, B, C and D helices), the N- and C-terminal loops and the connecting sequences between the helices. The peptides corresponding to the helices of the folded protein all exhibit preferences for helix-like conformations in solution. The conformational ensembles of the A- and D-helix peptides contain ordered helical forms, as shown by extensive series of medium-range nuclear Overhauser effect connectivities, while the B- and C-helix peptides exhibit conformational preferences for nascent helix. All four peptides adopt ordered helical conformations in mixtures of trifluoroethanol and water. The terminal and interconnecting loop peptides also appear to contain appreciable populations of conformers with backbone phi and psi angles in the alpha-region and include highly populated hydrophobic cluster and/or turn conformations in some cases. Trifluoroethanol is unable to drive these peptides towards helical conformations. Overall, the peptide fragments of myohemerythrin have a marked preference towards secondary structure formation in aqueous solution. In contrast, peptide fragments derived from the beta-sandwich protein plastocyanin are relatively devoid of secondary structure in aqueous solution (see accompanying paper). These results suggest that the two different protein structural motifs may require different propensities for formation of local elements of secondary structure to initiate folding, and that there is a prepartitioning of conformational space determined by the local amino acid sequence that is different for the helical and beta-sandwich structural motifs.     

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.     

Insights into the local residual entropy of proteins provided by NMR relaxation.

A simple model is used to illustrate the relationship between the dynamics measured by NMR relaxation methods and the local residual entropy of proteins. The expected local dynamic behavior of well-packed extended amino acid side chains are described by employing a one-dimensional vibrator that encapsulates both the spatial and temporal character of the motion. This model is then related to entropy and to the generalized order parameter of the popular "model-free" treatment often used in the analysis of NMR relaxation data. Simulations indicate that order parameters observed for the methyl symmetry axes in, for example, human ubiquitin correspond to significant local entropies. These observations have obvious significance for the issue of the physical basis of protein structure, dynamics, and stability.     

Relation of amino acid composition and the Moffitt parameters to the secondary structure of proteins

A statistical analysis of the relation between the amino acid composition of proteins and the amount of helical secondary structure as characterized by the Moffitt b₀ parameter has shown a high degree of correlation of the b₀ parameter with those amino acids whose homopolymers can form helical structures. Using the data for 107 proteins, a linear relation was found between b₀ and the sum of the residue percentages of alanine, arginine, aspartic acid, cysteine, glutamic acid, leucine, andlysine. A statistical analysis of the Moffitt a₀ parameter, on the other hand, showed no statistically significant grouping of amino acids in relation to the amount of secondary β structure in a protein. A plot of b₀ versus a₀–a, where a represents the a⁰ parameter for a fully denatured protein, for 55 proteins showed distinct nonlinearity. This nonlinearity was postulated to be due to presence of β structure, and a nomagram was constructed which allowed a semiquantitative estimate of the amount of helical and β-type secondary structures from the b₀ versus a₀–a plot.     

The major protein from lipid bodies of maize. Characterization and structure based on cDNA cloning

In plant seeds, the storage triacylglycerol is packed in discrete particles called lipid bodies which consist of a lipid core surrounded by a phospholipid monolayer with embedded proteins. We have cloned and sequenced a nearly full-length cDNA for the major protein (L3) associated with the lipid bodies of maize. The L3-cDNA clone was identified by hybrid-selected translation analysis and contains the complete 3' noncoding region and an open reading frame of 432 nucleotides. This open reading frame encodes a polypeptide with amino acid composition, hydrophobicity, and predicted protease digestion pattern which correlate well with those of the authentic L3 protein. Analyses of predicted secondary structure and local hydropathy of the deduced amino acid sequence suggest three structural domains in the protein. An internal domain of 72 contiguous hydrophobic or neutral amino acids is bounded at the amino-terminal side by a hydrophilic alpha-helix and on the carboxyl-terminal side by an amphipathic alpha-helix. The data suggest that L3 is uniquely suited to interact with both lipid and phospholipid moieties of the lipid body. A simple model for the topology of L3 on the lipid body is proposed. The unusual structure of the lipid body protein is discussed and compared to those of the two well-studied classes of lipid-associated proteins, apolipoproteins and intrinsic membrane proteins.     

植物萜类合成酶3-羟基-3-甲基戊二酰辅酶A还原酶的生物信息学分析

采用生物信息学的方法和工具对已在GenBank上注册的橡胶、烟草、辣椒、穿心莲等植物的萜类合成酶3-羟基-3-甲基戊二酰辅酶A还原酶的核酸及氨基酸序列进行分析,并对其组成成分、信号肽、跨膜拓朴结构域、疏水性/亲水性、蛋白质二级及三级结构、分子系统进化关系等进行预测和推断。结果表明该类酶基因的全长包括5′、3′非翻译区和一个开放阅读框,无信号肽,是一个跨膜的亲水性蛋白,包括两个功能HMG-CoA结合motif及两个功能NADPH结合motif,α-螺旋和不规则盘绕是蛋白质二级结构最大量的结构元件,β-转角和延伸链散布于整个蛋白质中,蛋白质的功能域在空间布局上折叠成“V”形,“V”形的两臂由螺旋状的N结构域和S结构域构成,中间部分由L结构域构成。     

Solvent interactions in nucleic acid crystal hydrates

A detailed knowledge of structural and energetic aspects of water-nucleic acid interactions is essential for understanding the role of solvent in stabilizing the various helical forms of nucleic acids. In this study, computer simulation techniques have been used to predict structural properties of solvent networks in small nucleic acid crystal hydrates. A detailed comparison of predicted and experimental results on the structure of the solvent networks is presented and includes an analysis of both the local environment and hydrogen bond pattern of each water molecule. A correlation between the environment of each unique water molecule and its energetic properties (such a dipole moment and binding energy) is seen. As in the previous studies on small amino acid hydrate crystals, non-pair additive (cooperative) effects are found to be non-negligible. It is concluded that the potential functions used in this initial study lead to simulated solvent networks in reasonable agreement with experimental data. Thus, it is now feasible to use them in studies of hydration of larger helical fragments of nucleic acids of more direct biological interest.     

Three-dimensional motion of avian spermatozoa

Observations have been made on spermatozoa from the domestic fowl, quail and pigeon (non-passerine birds) and also from the starling and zebra finch (passerine birds). In free motion, all these spermatozoa roll (spin) continuously about the progression axis, whether or not they are close to a plane surface. Furthermore, the direction of roll is consistently clockwise (as seen from ahead). The flagellar wave has been shown to be helical and dextral (as predicted) for domestic fowl sperm when they swim rapidly in low viscosity salines. Calculations have shown that their forward velocity is consistent with their induced angular velocity but that the size of the sperm head is suboptimal for progression speed under these conditions. Dextrally helical waves also occur on the distal flagellum of fowl, quail and pigeon sperm in high viscosity solutions. But in other cases, the mechanism of torque-generation is more problematical. The problem is most profound for passerine sperm, in that typically these cells spin rapidly while seeming to remain virtually straight. Because there is no evidence for a helical wave on these flagella, we have considered other possible means whereby rotation about the local flagellar axis (self-spin) might be achieved. Sometimes, passerine sperm, while maintaining their spinning motion, adopt a fixed curvature; this must be an instance of bend-transfer circumferentially around the axonemal cylinder-though the mechanism is obscure. It is suggested that the self-spin phenomenon may be occurring in non-passerine sperm that in some circumstances spin persistently, yet without expressing regular helical waves. More complex waves are apparent in non-passerine sperm swimming in high viscosity solutions: added to the small scale bends is a large scale, sinistrally helical curvature of the flagellum. It is argued that the flagellum follows this sinistrally helical path (i.e. "screws" though the fluid) because of the shape of the sperm head and the angle at which the flagellum is inserted into it. These conclusions concerning avian sperm motility are thought to have relevance to other animal groups. Also reported are relevant aspects of flagellar ultrastructure for pigeon and starling sperm.     

The complete amino acid sequence of a human folate binding protein from KB cells determined from the cDNA

Two species of folate binding protein (FBP), an integral membrane-associated form and a soluble secreted form, have been previously purified from cultured human KB cells. The complete nucleotide sequence of the complementary DNA (cDNA) clone for the coding region of the mature membrane-associated FBP has now been determined, and the deduced amino acid sequence has been computer-analyzed for a prediction of the secondary structure of the protein. The clone has 857 nucleotides of which 678 comprise the coding region for 226 amino acids. The deduced amino sequence contains the identical sequence of the published 18 NH2-terminal amino acids of the purified FBP from KB cells and the published partial amino acid sequence of the human milk FBP except for 1 residue. There was also over 90% homology with the published amino acid sequence of the bovine milk FBP. A total of 16 cysteine residues has been conserved in the three proteins indicating that this amino acid may provide a tertiary structure which is required for its ligand binding function. Northern blot analysis using the cDNA probe identified a single band of 1.28-kilobase pair mRNA in KB cells which was 4.7-fold more intense in folate-depleted cells than in normal cells. These results indicate that the membrane FBP and the soluble FBP in the medium are translation products of the same gene. Computer analysis of the deduced amino acid sequence indicates that there is only one stretch of amino acids of sufficient hydrophobicity and length to span the lipid bilayer of the plasma membrane, but it lacked a predictable helical structure. Those regions of the sequence which did have a predictable helical structure lacked sufficient hydrophobicity required for a membrane anchor. Thus, it is likely that the fatty acids previously reported to be present in the membrane-associated FBP from these cells rather than a peptide sequence provide an important membrane anchoring function.