PubMed Health⌕ Search

Biomedical subjects

K C Chou

Publications and source records attributed to K C Chou.

14 recordsLinked to original sources
PubMed1992

Role of loop-helix interactions in stabilizing four-helix bundle proteins.

One of the critical issues regarding proteins with a four-helix bundle motif is which interactions play the major role in stabilizing this type of folded structure: the interaction among the four alpha-helices or the interaction between the loop and helix segments. To answer this question, an energetic analysis has been carried out for three proteins with a four-helix bundle--namely, methemerythrin, cytochrome b-562, and cytochrome c'. The structures on which the analysis has been made were derived from their respective crystallographic coordinates. All three proteins have long helices (16-26 residues) and most of their loops are short (3-5 residues). However, it was found in all three proteins that loop-helix interactions were stronger than helix-helix interactions. Moreover, not only the nonbonded component but also the electrostatic component of the interaction energy were dominated by loop-helix interactions rather than by interhelix interactions, although the latter involve favorable helix-dipole interactions due to the antiparallel arrangement of neighboring helices. The results of the energetic analysis indicate that the loop segments, whether they are in a theoretical model or in real proteins, play a significant role in stabilizing proteins with four-helix bundles.

Cytochrome b Group↗
PubMed1992

A correlation-coefficient method to predicting protein-structural classes from amino acid compositions.

A protein is usually classified into one of the following four structural classes: all alpha, all beta, (alpha + beta) and alpha/beta. In this paper, based on the maximum correlation-coefficient principle, a new formulation is proposed for predicting the structural class of a protein according to its amino acid composition. Calculations have been made for a development set of proteins from which the amino acid compositions for the standard structural classes were derived, and an independent set of proteins which are outside the development set. The former can test the self consistency of a method and the latter can test its extrapolating effectiveness. In both cases, the results showed that the new method gave a considerably higher rate of correct prediction than any of the previous methods, implying that a significant improvement has been achieved by implementing the maximum-correlation-coefficient principle in the new method.

Algorithms↗
PubMed1992

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.

Amino Acid Sequence↗
PubMed1992

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.

Amino Acids↗
PubMed1992

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)

Amino Acid Sequence↗
PubMed1992

Strong electrostatic loop-helix interactions in bundle motif protein structures.

Based on CHARMM potential (Brooks et al., 1983) an energetic analysis has been carried out for four typical 4-alpha-helix bundle proteins, i.e., methemerythrin, cytochrome b-562, cytochrome c', and bovine somatotropin. The bovine somatotropin possesses long loops, but all the other three proteins have short loops. It was found that in all these four 4-alpha-helix bundle motif structures the interaction between loops and helices was much stronger than the interaction among the four helices themselves. Particularly for the electrostatic interaction energy, the loop-helix interaction is overwhelmingly stronger than the interhelix interaction although the latter involves the favorable helix dipole interaction due to the antiparallel arrangement of neighboring alpha-helices. The present study indicates that such a conclusion holds true regardless of what loops, long or short, are in the 4-alpha-helix bundle protein, and also regardless of which empirical potential, ECEPP or CHARMM, is used for calculations although in CHARMM the electrostatic energy is much more heavily emphasized than in ECEPP. Therefore, no appropriate conclusion can be drawn in arguing whether the dipole interaction among the four alpha-helices play a stabilizing role or destabilizing role for a 4-alpha-helix bundle protein without taking into consideration the effect of interaction between helices and loops. The calculated results reported here provide, from a different point of view, insights that might be useful for revealing the essence of the driving forces during the folding of proteins.

Animals↗
PubMed1992

Monte Carlo simulation studies on the prediction of protein folding types from amino acid composition.

In the methodology development for statistical prediction of protein structures, the founders of different methods usually selected different sets of proteins to test their predicted results. Therefore, it is hard to make a fair comparison according to the results they reported. Even if the predictions by different methods are performed for the same set of proteins, there is still such a problem: a method better that the other for one set of proteins would not necessarily remain so when applied to another set of proteins. To tackle this problem, a Monte Carlo simulation method is proposed to establish an objective criterion to measure the accuracy of prediction for the protein folding type. Such an objective accuracy is actually corresponding to the asymptotical limit genereated during the Monte Carlo simulation process. Based on that, it has been found that the average objective accuracy for predicting the all-alpha, all-beta, alpha + beta, and alpha/beta proteins by the least Euclid's distance method (Nakashima, H., K. Nishikawa, and T. Ooi. 1986. J. Biochem. 99:152-162) is 73.0% and that by the least Minkowski's distance method (Chou, P.Y. 1989. Prediction in Protein Structure and the Principles of Protein Conformation. Plenum Press. New York. 549-586) is 70.9%, indicating that the former is better than the latter. However, according to the original reports, the latter claimed a rate of correct prediction with 79.7% but the former with only 70.2%, leading to a completely opposite conclusion. This indicates the necessity of establishing an objective criterion, and a comparison is meaningful only when it is based on the objective criterion. The simulation method and the idea developed here also can be applied to examine any other statistical prediction methods.

Amino Acids↗
PubMed1992

Diagrammatization of codon usage in 339 human immunodeficiency virus proteins and its biological implication.

The occurrence frequencies of bases A (adenine), C (cytosine, G (guanine), and T (thymine) occurring in the 1st, 2nd, and 3rd codon positions in the codon usage table of viral genes for the 339 human immunodeficiency virus (HIV) proteins compiled recently have been calculated and diagrammatized. For comparison, the corresponding diagrammatic representations for the 2681 human proteins from the codon usage table for primate genes are also presented. The analyzed results based on these characteristic diagrams indicate that considerably similar features have been found between HIV and human proteins for the 1st and 2nd codon positions; i.e., they are all occupied predominantly by purine, especially base A. However, a significant difference in the 3rd codon position between HIV and human proteins has been observed; i.e., human proteins are of high C + G content and low A + G content in the 3rd codon position, whereas the case is just the opposite for HIV proteins. The biological implication of such a duality on the codon bias of HIV against human proteins is discussed. It is suggested that the 1st and 2nd codon positions can be termed as the structure-determining position, and the 3rd codon position termed as the species-determining position. The diagrammatic representation and analysis method described here possess a great potential for the study of molecular evolution from the viewpoint of the genetic code for which data have been accumulated rapidly and will continue to grow at a much faster pace.

Base Composition↗
PubMed1991

Mass-weighted molecular dynamics simulation of cyclic polypeptides.

A modified molecular dynamics (MD) method in which atomic masses are weighted was developed previously for studying the conformational flexibility of neuroregulating tetrapeptide Phe-Met-Arg-Phe-amide (FMRF-amide). The method has now been applied to longer and constrained molecules, namely a disulfide-linked cyclic hexapeptide, c[CYFQNC], and its linear and "pseudo-cyclic" analogues. The sampling of dehedral conformational space of teh linear hexapeptide in mass-weighted MD simulations was found to be improved significantly over conventional MD simulations, as in the case of the shorter FMRF-amide molecule studied previously. In the cyclic hexapeptide, the internal constraint of the molecule due to the intramolecular disulfide bond (hence the absence of free terminals in the molecule) does not adversely affect the significant improvement of conformational sampling in mass-weighted MD simulations over normal MD simulations. The pseudo-cyclic polypeptide is identical to the linear CYFQNC molecule in amino acid sequence (i.e., side chains of the cysteine residues are reduced), but the positions of its two terminal heavy atoms were held fixed in space such that the molecule has a nearly cyclic conformation. For this molecule, the mass-weighted MD simulation generated a wide range of polypeptide backbone conformations covering the internal dihedral degrees of freedom; moreover, the physical space of the pseudo-cyclic structure was also sampled in a complete revolution of the entire molecular fragment about the two fixed termini during the simulation. These characteristics suggest that mass-weighted MD can also be an extremely useful method for conformational analyses of constrained molecules and, in particular, for modeling loops on protein surfaces.

Amino Acid Sequence↗
PubMed1991

Energetic approach to the folding of alpha/beta barrels.

The folding of a polypeptide into a parallel (alpha/beta)8 barrel (which is also called a circularly permuted beta 8 alpha 8 barrel) has been investigated in terms of energy minimization. According to the arrangement of hydrogen bonds between two neighboring beta-strands of the central barrel therein, such an alpha/beta barrel structure can be folded into six different types: (1) left-tilted, left-handed crossover; (2) left-tilted, right-handed crossover; (3) nontilted, left-handed crossover; (4) nontilted, right-handed crossover; (5) right-tilted, left-handed crossover; and (6) right-tilted, right-handed crossover. Here "tilt" refers to the orientational relation of the beta-strands to the axis of the central beta-barrel, and "crossover" to the beta alpha beta folding connection feature of the parallel beta-barrel. It has been found that the right-tilted, right-handed crossover alpha/beta barrel possesses much lower energy than the other five types of alpha/beta barrels, elucidating why the observed alpha/beta barrels in proteins always assume the form of right tilt and right-handed crossover connection. As observed, the beta-strands in the energy-minimized right-tilted, right-handed crossover (alpha/beta)8-barrel are of strong right-handed twist. The value of root-mean-square fits also indicates that the central barrel contained in the lowest energy (alpha/beta)8 structure thus found coincides very well with the observed 8-stranded parallel beta-barrel in triose phosphate isomerase (TIM). Furthermore, an energetic analysis has been made demonstrating why the right-tilt, right-handed crossover barrel is the most stable structure. Our calculations and analysis support the principle that it is possible to account for the main features of frequently occurring folding patterns in proteins by means of conformational energy calculations even for very complicated structures such as (alpha/beta)8 barrels.

Hydrogen Bonding↗
PubMed1991

Simulated annealing approach to the study of protein structures.

One of the most difficult problems in predicting the three dimensional structure of proteins is how to deal with the local minimum problem. In many cases of practical interest this problem has been reduced to how to select an appropriate set of starting conformations for carrying out energy minimizations. How these starting conformations are selected, however, is often based on the physical intuition of the person doing the calculations, and hence it is hard to avoid bearing some sort of arbitrariness. To improve such a situation, we introduced the simulated annealing Monte Carlo algorithm to locate the optimal starting conformations for energy minimizations. The method developed here is valid for both single and multiple polypeptide chain systems. The annealing process can be conducted with respect to either the internal dihedral angles of a polypeptide chain or the external rotations and translations of various constituent polypeptide chains, and hence is particularly useful for studying the packing arrangements of secondary structures in proteins, such as helix/helix packing, helix/sheet packing and sheet/sheet packing. It was shown via a number of comparative calculations that the final structures obtained through the annealing process not only had lower energies than the corresponding energy-minimized structures reported previously, but also assumed the forms closer to the observations in proteins. All these results indicate that a better result can be obtained in search of low-energy structures of proteins by incorporating the simulated annealing approach.(ABSTRACT TRUNCATED AT 250 WORDS)

Algorithms↗
PubMed1991

Geometric and energy parameters in lysine-retinal chromophores.

The parameters used in the computer program ECEPP (Empirical Conformational Energy Program for Peptides) have been expanded to cover some key elements in retinal-containing proteins. These elements are 'all-trans retinal lysine with unprotonated imine', 'all-trans retinal lysine with protonated imine', '13-cis retinal lysine with unprotonated imine' and '13-cis retinal lysine with protonated imine' respectively. The geometric parameters of these four new 'amino acid residues' were derived by optimizing their molecular structures with the AM1 Hamiltonian included in MOPAC (Molecular Orbital PACkage), and their partial atomic charges were determined with a CNDO/2 (Complete Neglect of Differential Overlap) calculation. The parameters for nonbonded interactions and torsional potentials were obtained from the existing ECEPP parameters through a logical extension. The augmented ECEPP system thus obtained can be employed to investigate the conformation of bacteriorhodopsin and its proton-pumping mechanism from an energetic point of view. The computer modeling study on bacteriorhodopsin and other seven-helix membrane proteins, e.g. serotonin receptor and dopamine receptor, is under way in the Upjohn Laboratories.

Bacteriorhodopsins↗
PubMedDate not supplied

The quantitative relations between diffusion-controlled reaction rate and characteristic parameters in enzyme-substrate reaction systems. I. Neutral substrates.

With the development of techniques for studying fast reactions, we are confronted with some experimental observations which, according to the classical diffusion-controlled reaction theory based on the assumption of pherical symmetry, will lead to paradoxical consequences. Recently, some investigators have attempted to present a theory of non-spherically symmetric diffusion-controlled reaction. Owing to the difficulty in mathematics, only the formal discussion has been given, and the numerical results have not yet been reported. The combination between enzyme and substrate is not only non-spherically symmetric, but also related to the force field existing between the reacting molecules. In the previous article we introduced the spatial factor and the force field factor and derived a general equation for these reaction systems. We present in this article a general method for finding the numerical solution of the above equation. With the numerical results obtained by this method we may interpret the experimental facts which can not be explained from the viewpoint of the classical diffusion-controlled reaction theory. In order to examine more closely the kinetic behaviour of enzyme-substrate reaction systems, we defined several characteristic parameters. The quantitative relation between these parameters and the rate of diffusion-controlled reaction has been discussed.

Binding Sites↗