Estimating residue evolutionary conservation by introducing von Neumann entropy and a novel gap-treating approach.

Evolutionary conservation derived from a multiple sequence alignment is a powerful indicator of the functional significance of a residue, and it can help to predict active sites, ligand-binding sites, and protein interaction interfaces. The results of the existing algorithms in identifying the residue's conservation strongly depend on the sequence alignment, making the results highly variable. Here, by introducing the amino acid similarity matrix, we propose a novel gap-treating approach by combining the evolutionary information and von Neumann entropies to compute the residue conservation scores. It is indicated through a series of tested results that the new approach is quite encouraging and promising and may become a useful tool in complementing the existing methods.

Cleavage mechanism of the H5N1 hemagglutinin by trypsin and furin.

The cleavage property of hemagglutinin (HA) by different proteases was the prime determinant for influenza A virus pathogenicity. In order to understand the cleavage mechanism, molecular modeling tools were utilized to study the coupled model systems of the proteases, i.e., trypsin and furin and peptides of the cleavage sites specific to H5N1 and H1 HAs, which constitute models of HA precursor in complex with cleavage proteases. The peptide segments 'RERRRKKR downward arrow G' and 'SIQSR downward arrow G' from the high pathogenic H5N1 H5 and the low pathogenic H1N1 H1 cleavage sites were docking to the trypsin and furin active pockets, respectively. It was observed through the docking studies that trypsin was able to recognize and cleave both the high pathogenic and low pathogenic hemagglutinin, while furin could only cleave the high pathogenic hemagglutinin. An analysis of binding energies indicated that furin got most of its selectivity due to the interactions with P(1), P(4), and P(6), while having less interaction with P(2) and little interactions with P(3), P(5), P(7), and P(8). Some mutations of H5N1 H5 cleavage sequence fitted less well into furin and would reduce high pathogenicity of the virus. These findings hint that we should focus at the subsites P(1), P(4), and P(6) for developing drugs against H5N1 viruses.

Euk-PLoc: an ensemble classifier for large-scale eukaryotic protein subcellular location prediction.

With the avalanche of newly-found protein sequences emerging in the post genomic era, it is highly desirable to develop an automated method for fast and reliably identifying their subcellular locations because knowledge thus obtained can provide key clues for revealing their functions and understanding how they interact with each other in cellular networking. However, predicting subcellular location of eukaryotic proteins is a challenging problem, particularly when unknown query proteins do not have significant homology to proteins of known subcellular locations and when more locations need to be covered. To cope with the challenge, protein samples are formulated by hybridizing the information derived from the gene ontology database and amphiphilic pseudo amino acid composition. Based on such a representation, a novel ensemble hybridization classifier was developed by fusing many basic individual classifiers through a voting system. Each of these basic classifiers was engineered by the KNN (K-Nearest Neighbor) principle. As a demonstration, a new benchmark dataset was constructed that covers the following 18 localizations: (1) cell wall, (2) centriole, (3) chloroplast, (4) cyanelle, (5) cytoplasm, (6) cytoskeleton, (7) endoplasmic reticulum, (8) extracell, (9) Golgi apparatus, (10) hydrogenosome, (11) lysosome, (12) mitochondria, (13) nucleus, (14) peroxisome, (15) plasma membrane, (16) plastid, (17) spindle pole body, and (18) vacuole. To avoid the homology bias, none of the proteins included has > or =25% sequence identity to any other in a same subcellular location. The overall success rates thus obtained via the 5-fold and jackknife cross-validation tests were 81.6 and 80.3%, respectively, which were 40-50% higher than those performed by the other existing methods on the same strict dataset. The powerful predictor, named "Euk-PLoc", is available as a web-server at http://202.120.37.186/bioinf/euk . Furthermore, to support the need of people working in the relevant areas, a downloadable file will be provided at the same website to list the results predicted by Euk-PLoc for all eukaryotic protein entries (excluding fragments) in Swiss-Prot database that do not have subcellular location annotations or are annotated as being uncertain. The large-scale results will be updated twice a year to include the new entries of eukaryotic proteins and reflect the continuous development of Euk-PLoc.

Prediction of linear B-cell epitopes using amino acid pair antigenicity scale.

Identification of antigenic sites on proteins is of vital importance for developing synthetic peptide vaccines, immunodiagnostic tests and antibody production. Currently, most of the prediction algorithms rely on amino acid propensity scales using a sliding window approach. These methods are oversimplified and yield poor predicted results in practice. In this paper, a novel scale, called the amino acid pair (AAP) antigenicity scale, is proposed that is based on the finding that B-cell epitopes favor particular AAPs. It is demonstrated that, using SVM (support vector machine) classifier, the AAP antigenicity scale approach has much better performance than the existing scales based on the single amino acid propensity. The AAP antigenicity scale can reflect some special sequence-coupled feature in the B-cell epitopes, which is the essence why the new approach is superior to the existing ones. It is anticipated that with the continuous increase of the known epitope data, the power of the AAP antigenicity scale approach will be further enhanced.

Virtual screening for finding natural inhibitor against cathepsin-L for SARS therapy.

Recently Simmons et al. reported a new mechanism for SARS virus entry into target cells, where MDL28170 was identified as an efficient inhibitor of CTSL-meditated substrate cleavage with IC(50) of 2.5 nmol/l. Based on the molecule fingerprint searching method, 11 natural molecules were found in the Traditional Chinese Medicines Database (TCMD). Molecular simulation indicates that the MOL376 (a compound derived from a Chinese medicine herb with the therapeutic efficacy on the human body such as relieving cough, removing the phlegm, and relieving asthma) has not only the highest binding energy with the receptor but also the good match in geometric conformation. It was observed through docking studies that the van der Waals interactions made substantial contributions to the affinity, and that the receptor active pocket was too large for MDL21870 but more suitable for MOL736. Accordingly, MOL736 might possibly become a promising lead compound for CTSL inhibition for SARS therapy.

Inverse Association Between Serum Osteoprotegerin and Bone Mineral Density in Renal Transplant Recipients.

BACKGROUND: Osteoprotegerin (OPG) has pleiotropic effects on bone metabolism as well as endocrine function. Our aim was to evaluate the relationship between bone mineral density (BMD) and serum OPG concentration in renal transplant recipients. METHODS: Fasting blood samples were obtained from 69 renal transplant recipients. BMD was measured in lumbar vertebrae (L2-L4) by dual-energy X-ray absorptiometry. Eight patients (11.6%) had BMD values indicative of osteoporosis, 28 patients (40.6%) had BMD values indicative of osteopenia, and 33 patients had normal BMD values. Increased serum OPG levels (P < .001), decreased body mass index (BMI) (P = .033), and decreased body weight (P = .010) were significantly correlated with low lumbar T-score cut-off points between groups (normal, osteopenia, and osteoporosis). RESULTS: Women had significantly lower lumbar BMD values than men (P = .013). Menopause (P = .005), use of tacrolimus (P = .020), and use of cyclosporine (P = .046) were associated with lower lumbar BMD in renal transplant recipients. Univariate linear regression analysis revealed that lumbar BMD was positively correlated with height (P = .016), body weight (P = .001), and BMI (P = .015) and negatively correlated with age (P = .039) and log-OPG (P = .001). Multivariate linear regression analysis revealed that log-OPG (ß: -0.275, R(2) change = 0.154, P = .014), body weight (ß: 0.334, R(2) change = 0.073, P = .004), and age (ß: -0.285, R(2) change = 0.079, P = .008) were independent predictors of lumbar BMD values in renal transplant recipients. CONCLUSIONS: Serum OPG concentration correlated negatively with lumbar BMD values in renal transplant recipients.

A semi-analytical expression for the concentration distribution of substrate molecules in fast, enzyme-catalysed reaction systems.

A semi-analytical expression for the concentration distribution of the substrate molecules in fast enzyme-catalysed reaction system is presented. On this basis, the physical pictures of the reactions are discussed.

Energy-optimized structure of antifreeze protein and its binding mechanism.

A combination of Monte Carlo simulated annealing and energy minimization was utilized to determine the conformation of the antifreeze protein from the fish winter flounder. It was found from the energy-optimized structure that the hydroxyl groups of its four threonine residues, i.e. Thr2, Thr13, Thr24, Thr35, are aligned on almost the same line parallel to the helix axis and separated successively by 16.1, 16.0 and 16.2 A, respectively, very close to the 16.6 A repeat spacing along [0112] in ice. Based on such a space match, a zipper-like model is proposed to elucidate the binding mechanism of the antifreeze protein to ice crystals. According to the current model, the antifreeze protein may bind to an ice nucleation structure in a zipper-like fashion through hydrogen bonding of the hydroxyl groups of these four Thr residues to the oxygen atoms along the [0112] direction in ice lattice, subsequently stopping or retarding the growth of ice pyramidal planes so as to depress the freeze point. The calculated results and the binding mechanism thus derived accord with recent experimental observations. The mechanistic implications derived from such a special antifreeze molecule might be generally applied to elucidate the structure-function relationship of other antifreeze proteins with the following two common features: (1) recurrence of a Thr residue (or any other polar amino acid residue whose side-chain can form a hydrogen bond with water) in an 11-amino-acid period along the sequence concerned; and (2) a high percentage of Ala residue component therein. Further experiments are suggested to test the ice binding model.

A graphic approach to analyzing codon usage in 1562 Escherichia coli protein coding sequences.

The occurrence frequencies of the four bases (adenine, cytosine, guanine and thymine) at each of the three codon positions for 1562 Escherichia coli protein coding sequences have been calculated. The 1562 x 4 x 3 = 18,744 data thus obtained have been analyzed by a graphic method in which the four base occurrence frequencies at each codon position for each coding sequence are represented by a point in a three-dimensional space. Thus, the 18,744 data, which would otherwise occupy several printed pages, can be intuitively displayed by a graphy. The point distribution pattern for each of the three codon positions has been analyzed. The results of our analysis indicate that the patterns for the first two codon positions reflect the origin for producing native folding structures of proteins. We thus come to the conclusion that the distribution patterns for the first two codon positions should be basically species-independent, as confirmed by studies for a number of other species. However, the distribution pattern for the third codon position is species-dependent. Based on the point distribution of the third codon position, six collective parameters have been defined to describe the overall feature of the pattern concerned. These collective parameters can be generally used to classify different species, and hence would be a useful vehicle for studies in taxonomy. In addition to E. coli, the collective parameters for a number of other species have been calculated and analyzed.

Role of interchain interactions in the stabilization of the right-handed twist of beta-sheets.

Conformational energy computations have been carried out for parallel and antiparallel beta-sheets composed of poly-L-Val and poly-L-Ile peptide chains, each consisting of four and of six residues, respectively, with CH3CO- and-NHCH3 end groups. The beta-sheets considered contained three and five equivalent chains, respectively. All computed minimum-energy beta-sheets were found to have a large right-handed twist of a magnitude that corresponds to the mean twist of beta-sheets observed in globular proteins. The twist has the same sign but is much larger than in beta-sheets of poly-L-Ala, because of intra- and interchain interactions between the bulky beta-branched side-chains. While the right-handed twist is a result of intrachain interactions between side-chains in the case of poly-L-Val, these interactions would favor a left-handed twist in poly-L-Ile, and the right-handed twist in the latter is a result of interchain interactions. Parallel beta-sheets are more stable than antiparallel sheets for both poly-L-Val and poly-L-Ile, in contrast to poly-L-Ala. This result agrees with observations on the preferred orientation of the chains in oligopeptides that form beta-structures. It also explains the observed high relative frequencies of occurrence of Val and Ile residues in parallel beta-sheets, as compared with antiparallel sheets, in globular proteins.

Conformational and geometrical properties of idealized beta-barrels in proteins.

An equation for calculating the distances between the atoms involved in forming an idealized hydrogen bond in a parallel or antiparallel beta-barrel has been derived by adjusting the corresponding data given by Pauling and Corey for a beta-sheet. Based on these distances, a geometrical optimization method was developed, by which one can generate various idealized beta-barrels: parallel or antiparallel, tilted or non-tilted, right-tilted or left-tilted. For each type of idealized beta-barrel thus obtained, the corresponding conformation and characteristic geometric parameters as well as their relationship are analyzed and discussed. Since the strand in a tilted beta-barrel traces a curve rather than a straight line on a cylinder-like surface, a regular chain in which the dihedral angles of each residue are the same cannot form a tilted beta-barrel but only a non-tilted beta-barrel. As observed, the strands of a right-tilted beta-barrel possess a very strong right-handed twist. The radii of the idealized tilted parallel and antiparallel beta-barrels are greater than those of the corresponding non-tilted ones by approximately 1 A and approximately 1.5 A, respectively. Consequently, there is relatively more room for a tilted beta-barrel to accommodate the internal side-chains, suggesting that a conformational change from a non-tilted beta-barrel to a tilted one would ease the repulsion among the crowded internal side-chains so as to make the structure more stable. The values of root-mean-square fits indicate that the idealized right-tilted beta-barrels coincide quite well with the observed beta-barrels in both parallel and antiparallel cases.

Energy of stabilization of the right-handed beta alpha beta crossover in proteins.

An explanation in terms of conformational energies is provided for the observed nearly exclusive preference of the beta alpha beta structure for forming a right-handed, rather than a left-handed, crossover connection. Conformational energy computations have been carried out on a model beta alpha beta structure, consisting of two six-residue Val beta-strands and of a 12-residue Ala alpha-helix, connected by two flexible four-residue Ala links to the strands. The energy of the most favorable right-handed crossover is 15.51 kcal/mol lower than that of the corresponding left-handed cross-over. The right-handed crossover is a strain-free structure. Its energy of stabilization arises largely from the interactions of the two beta-strands with one another and with the alpha-helix. On the other hand, the left-handed crossover is either disrupted after energy minimization or it remains conformationally strained, as indicated by an energetically unfavorable left twisting of the beta-sheet and by the presence of high-energy local residue conformations. In the energetically most favorable right-handed crossover, the right twisting of the beta-sheet and its manner of interacting with the alpha-helix are identical with those computed earlier for isolated beta-sheets and for packed alpha/beta structures. This result supports a proposed principle that it is possible to account for the main features of frequently occurring structural arrangements in globular proteins in terms of the properties of their component structural elements.

Interactions between two beta-sheets. Energetics of beta/beta packing in proteins.

The analysis of the interactions between regularly folded segments of the polypeptide chain contributes to an understanding of the energetics of protein folding. Conformational energy-minimization calculations have been carried out to determine the favorable ways of packing two right-twisted beta-sheets. The packing of two five-stranded beta-sheets was investigated, with the strands having the composition CH3CO-(L-Ile)6-NHCH3 in one beta-sheet and CH3CO-(L-Val)6-NHCH3 in the other. Two distinct classes of low-energy packing arrangements were found. In the class with lowest energies, the strands of the two beta-sheets are aligned nearly parallel (or antiparallel) with each other, with a preference for a negative orientation angle, because this arrangement corresponds to the best complementary packing of the two twisted saddle-shaped beta-sheets. In the second class, with higher interaction energies, the strands of the two beta-sheets are oriented nearly perpendicular to each other. While the surfaces of the two beta-sheets are not complementary in this arrangement, there is good packing between the corner of one beta-sheet and the interior part of the surface of the other, resulting in a favorable energy of packing. Both classes correspond to frequently observed orientations of beta-sheets in proteins. In proteins, the second class of packing is usually observed when the two beta-sheets are covalently linked, i.e. when a polypeptide strand passes from one beta-sheet to the other, but we have shown here that a large contribution to the stabilization of this packing arrangement arises from noncovalent interactions.

Interactions between an alpha-helix and a beta-sheet. Energetics of alpha/beta packing in proteins.

Conformational energy computations have been carried out to determine the favorable ways of packing a right-handed alpha-helix on a right-twisted antiparallel or parallel beta-sheet. Co-ordinate transformations have been developed to relate the position and orientation of the alpha-helix to the beta-sheet. The packing was investigated for a CH3CO-(L-Ala)16-NHCH3 alpha-helix interacting with five-stranded beta-sheets composed of CH3CO-(L-Val)6-NHCH3 chains. All internal and external variables for both the alpha-helix and the beta-sheet were allowed to change during energy minimization. Four distinct classes of low-energy packing arrangements were found for the alpha-helix interacting with both the parallel and the anti-parallel beta-sheet. The classes differ in the orientation of the axis of the alpha-helix relative to the direction of the strands of the right-twisted beta-sheet. In the class with the most favorable arrangement, the alpha-helix is oriented along the strands of the beta-sheet, as a result of attractive non-bonded side-chain-side-chain interactions along the entire length of the alpha-helix. A class with nearly perpendicular orientation of the helix axis to the strands is also of low energy, because it allows similarly extensive attractive interactions. In the other two classes, the helix is oriented diagonally relative to the strands of the beta-sheet. In one of them, it interacts with the convex surface near the middle of the saddle-shaped twisted beta-sheet. In the other, it is oriented along the concave diagonal of the beta-sheet and, therefore, it interacts only with the corner regions of the sheet, so that this packing is energetically less favorable. The packing arrangements involving an antiparallel and a parallel beta-sheet are generally similar, although the antiparallel beta-sheet has been found to be more flexible. The major features of 163 observed alpha/beta packing arrangements in 37 proteins are accounted for in terms of the computed structural preferences. The energetically most favored packing arrangement is similar to the right-handed beta alpha beta crossover structure that is observed in proteins; thus, the preference for this connectivity arises in large measure from this energetically favorable interaction.

Vibrational properties of CO at the Pt(111)-solution interface: the anomalous Stark-tuning slope.

Vibrational properties of CO have been studied on Pt(111) in acid and alkaline electrolytes by synchronous measurements of CO oxidation current (0.5 mV/s) and IRAS spectra (one spectrum for every 1 mV). We found that in acid solutions the frequency-tuning rate (dnu(CO)/dE) as well as the potential-dependent bandwidth (dDeltanu1/2/dE) deviates from expected linear relationships. This unusual potential-dependent behavior is interpreted in terms of compression/dissipation of CO islands during the CO oxidation, engendered by competitive adsorption between inactive anions from a supporting electrolyte and the reactive OH species.

A sequence-coupled vector-projection model for predicting the specificity of GalNAc-transferase.

The specificity of GalNAc-transferase is consistent with the existence of an extended site composed of nine subsites, denoted by R4, R3, R2, R1, R0, R1', R2', R3', and R4', where the acceptor at R0 is either Ser or Thr to which the reducing monosaccharide is being anchored. To predict whether a peptide will react with the enzyme to form a Ser- or Thr-conjugated glycopeptide, a new method has been proposed based on the vector-projection approach as well as the sequence-coupled principle. By incorporating the sequence-coupled effect among the subsites, the interaction mechanism among subsites during glycosylation can be reflected and, by using the vector projection approach, arbitrary assignment for insufficient experimental data can be avoided. The very high ratio of correct predictions versus total predictions for the data in both the training and the testing sets indicates that the method is self-consistent and efficient. It provides a rapid means for predicting O-glycosylation and designing effective inhibitors of GalNAc-transferase, which might be useful for targeting drugs to specific sites in the body and for enzyme replacement therapy for the treatment of genetic disorders.

An optimization approach to predicting protein structural class from amino acid composition.

Proteins are generally classified into four structural classes: all-alpha proteins, all-beta proteins, alpha + beta proteins, and alpha/beta proteins. In this article, a protein is expressed as a vector of 20-dimensional space, in which its 20 components are defined by the composition of its 20 amino acids. Based on this, a new method, the so-called maximum component coefficient method, is proposed for predicting the structural class of a protein according to its amino acid composition. In comparison with the existing methods, the new method yields a higher general accuracy of prediction. Especially for the all-alpha proteins, the rate of correct prediction obtained by the new method is much higher than that by any of the existing methods. For instance, for the 19 all-alpha proteins investigated previously by P.Y. Chou, the rate of correct prediction by means of his method was 84.2%, but the correct rate when predicted with the new method would be 100%! Furthermore, the new method is characterized by an explicable physical picture. This is reflected by the process in which the vector representing a protein to be predicted is decomposed into four component vectors, each of which corresponds to one of the norms of the four protein structural classes.

Singular points of protein beta-sheets.

Protein beta-sheets can be regarded as surfaces. Two surfaces can be connected along a common edge to form a larger surface, or two edges of a surface can coalesce to form a closed sheet such as a beta-barrel. Singular points are locations where these connections are not perfect. In protein beta-sheets, a singular point is characterized by a residue separating two beta-ladders. In this paper, we study the singular points of protein beta-sheets from the surface topologic viewpoint, summarize our search results from the protein structural data in the Protein Data Bank, and present examples where singular points are near the active sites and may contribute to forming the proper relative positions of catalytic residues.

An energy-based approach to packing the 7-helix bundle of bacteriorhodopsin.

Based on the heavy-atom coordinates determined by the electron microscopy for the seven main helical regions of bacteriorhodopsin with the all-trans retinal isomer, energy optimizations were carried out for helix bundles containing the all-trans retinal and 13-cis retinal chromophores, respectively. A combination of simulated annealing and energy minimization was utilized during the process of energy optimization. It was found that the 7-helix bundle containing the all-trans isomer is about 10 kcal/mol lower in conformational energy than that containing the 13-cis isomer. An energetic analysis indicates that such a difference in energy is consistent with the observation that absorption of a 570-nm proton is required for the conversion of a bacteriorhodopsin from its all-trans to 13-cis form. It was also found that the above conversion process is accompanied by a significant conformational perturbation around the chromophore, as reflected by the fact that the beta-ionone ring of retinal moves about 5.6 A along the direction perpendicular to the membrane plane. This is consistent with the observation by Fodor et al. (Fodor, S.P.A., Ames, J.B., Gebhard, R., van der Berg, E.M.M., Stoeckenius, W., Lugtenburg, J., & Mathies, R.A., 1988, Biochemistry 27, 7097-7101). Furthermore, it is interesting to observe that although the retinal chromophore undergoes a significant change in its spatial position, the orientation of its transition dipole changes only slightly, in accord with experimental observations. In other words, even though orientation of the retinal transition dipole is very restricted, there is sufficient room, and degrees of freedom, for the retinal chromophore to readjust its position considerably. This finding provides new insight into the subtle change of the retinal microenvironment, which may be important for revealing the proton-pumping mechanism of bacteriorhodopsin. The importance of electrostatic and nonbonded interactions in stabilizing the 7-helix bundle structure has also been analyzed. Electrostatic interactions favor an antiparallel arrangement among adjacent helices. Nonbonded interactions, however, drive most of the closely packed helices into an arrangement in which the packing angles lie around -160 degrees, a value very near the -154 degrees value computed earlier as the most favorable packing arrangement of two poly(Ala) alpha-helices (Chou, K.-C., Némethy, G., & Scheraga, H.A., 1983, J. Phys. Chem. 87, 2869-2881). The structural features of the 7-helix bundle and their relationship to those found in typical 4-helix bundle proteins are also discussed.(ABSTRACT TRUNCATED AT 400 WORDS)