Search

One of the simplest ways to predict ligand binding sites is to identify pocket-shaped regions on the protein surface. Many programs have already been proposed to identify these pocket regions. Examination of their algorithms revealed that a pocket intrinsically has two arbitrary properties, "size" and "depth". We proposed a new definition for pockets using two explicit adjustable parameters that correspond to these two arbitrary properties. A pocket region is defined as a space into which a small probe can enter, but a large probe cannot. The radii of small and large probe spheres are the two parameters that correspond to the "size" and "depth" of the pockets, respectively. These values can be adjusted individual putative ligand molecule. To determine the optimal value of the large probe spheres radius, we generated pockets for thousands of protein structures in the database, using several size of large probe spheres, examined the correspondence of these pockets with known binding site positions. A new measure of shallowness, a minimum inaccessible radius, R(inaccess), indicated that binding sites of coenzymes are very deep, while those for adenine/guanine mononucleotide have only medium shallowness and those for short peptides and oligosaccharides are shallow. The optimal radius of large probe spheres was 3-4 A for the coenzymes, 4 A for adenine/guanine mononucleotides, and 5 A or more for peptides/oligosaccharides. Comparison of our program with two other popular pocket-finding programs showed that our program had a higher performance of detecting binding pockets, although it required more computational time., ((c) 2007 Wiley-Liss, Inc.)

Despite sharing many common features, adenine-binding and guanine-binding sites in proteins often show a clear preference for the cognate over the non-cognate ligand. We have analyzed electrostatic potential (ESP) patterns at adenine and guanine-binding sites of a large number of non-redundant proteins where each binding site was first annotated as adenine/guanine-specific or non-specific from a survey of primary literature. We show that more than 90% of ESP variance at the binding sites is accounted for by only two principal component ESP vectors, each aligned to molecular dipoles of adenine and guanine. Projected on these principal component vectors, the adenine/guanine-specific and non-specific binding sites, including adenine-containing dinucleotides, show non-overlapping distributions. Adenine or guanine specificities of the binding sites also show high correlation with the corresponding electrostatic replacement (cognate by non-cognate ligand) energies. High correlation coefficients (0.94 for 35 adenine-binding sites and 1.0 for 20 guanine-binding sites) were obtained when adenine/guanine specificities were predicted using the replacement energies. Our results demonstrate that ligand-free protein ESP is an excellent indicator for discrimination between adenine and guanine-specific binding sites and that ESP of ligand-free protein can be used as a tool to annotate known and putative purine-binding sites in proteins as adenine or guanine-specific.

Superfamily classifications are based variably on similarity of sequences, global folds, local structures, or functions. We have examined the possibility of defining superfamilies purely from the viewpoint of the global fold/function relationship. For this purpose, we first classified protein domains according to the beta-sheet topology. We then introduced the concept of kinship relations among the classified beta-sheet topology by assuming that the major elementary event leading to creation of a new beta-sheet topology is either an addition or deletion of one beta-strand at the edge of an existing beta-sheet during the molecular evolution. Based on this kinship relation, a network of protein domains was constructed so that the distance between a pair of domains represents the number of evolutionary events that lead one from the other domain. We then mapped on it all known domains with a specific core chemical function (here taken, as an example, that involving ATP or its analogs). Careful analyses revealed that the domains are found distributed on the network as >20 mutually disjointed clusters. The proteins in each cluster are defined to form a fold-based superfamily. The results indicate that >20 ATP-binding protein superfamilies have been invented independently in the process of molecular evolution, and the conservative evolutionary diffusion of global folds and functions is the origin of the relationship between them.

In this paper we present a new approach to study the volume fluctuations of proteins. From a 1 ns molecular dynamics simulation, the volume fluctuation of human lysozyme has been calculated. We used two different ways for the calculation. In the first one, the volume fluctuation is extracted directly from the trajectory. For the second one, a newly developed formalism based on principal component analysis is used. The r.m.s. volume fluctuations obtained from the two analyses agree well with each other. The isothermal intrinsic compressibility was found to be larger than the one reported by experiment. The difference is discussed and suggested to exist in the assumed uncertainty of the compressibility of hydrated water to deduce the isothermal intrinsic compressibility from the experimental value. Spectral analysis shows that low-frequency dynamics dominate the total volume fluctuation. The same aspect is found in the study using principal component analysis. This low-frequency region is related to large and slow motions of proteins. Therefore a long time dynamics simulation is necessary to describe the volume fluctuations of proteins.

We carry out a systematic analysis of the correlation between similarity of protein three-dimensional structures and their evolutionary relationships. The structural similarity is quantitatively identified by an all-against-all comparison of the spatial arrangement of secondary structural elements in nonredundant 967 representative proteins, and the evolutionary relationship is judged according to the definition of superfamily in the SCOP database. We find the following symmetry rule: a protein pair that has similar folds but belong to different superfamilies has (with a very rare exception) certain internal symmetry in its common similar folds. Possible reasons behind the symmetry rule are discussed.

Currently, collective coordinates are commonly employed in order to examine protein dynamics. In recent studies, they have been successfully applied to finding functionally relevant motions, to investigating the physical nature of protein dynamics, to sampling of the conformational space and to the analysis of experimental data. Collective coordinates also have other possible applications.

In order to search for a common structural motif in the phosphate-binding sites of protein-mononucleotide complexes, we investigated the structural variety of phosphate-binding schemes by an all-against-all comparison of 491 binding sites found in the Protein Data Bank. We found four frequently occurring structural motifs composed of protein atoms interacting with phosphate groups, each of which appears in different protein superfamilies with different folds. The most frequently occurring motif, which we call the structural P-loop, is shared by 13 superfamilies and is characterized by a four-residue fragment, GXXX, interacting with a phosphate group through the backbone atoms. Various sequence motifs, including Walker's A motif or the P-loop, turn out to be a structural P-loop found in a few specific superfamilies. The other three motifs are found in pairs of superfamilies: protein kinase and glutathione synthetase ATPase domain like, actin-like ATPase domain and nucleotidyltransferase, and FMN-linked oxidoreductase and PRTase.

We have developed a method of searching for similar spatial arrangements of atoms around a given chemical moiety in proteins that bind a common ligand. The first step in this method is to consider a set of atoms that closely surround a given chemical moiety. Then, to compare the spatial arrangements of such surrounding atoms in different proteins, they are translated and rotated so that the chemical moieties are superposed on each other. Spatial arrangements of surrounding atoms in a pair of proteins are judged to be similar, when there are many corresponding atoms occupying similar spatial positions. Because the method focuses on the arrangements of surrounding atoms, it can detect structural similarities of binding sites in proteins that are dissimilar in their amino acid sequences or in their chain folds. We have applied this method to identify modes of nucleotide base recognition by proteins. An all-against-all comparison of the arrangements of atoms surrounding adenine moieties revealed an unexpected structural similarity between protein kinases, cAMP-dependent protein kinase (cAPK), and casein kinase-1 (CK1), and D-Ala:D-Ala ligase (DD-ligase) at their adenine-binding sites, despite a lack of similarity in their chain folds. The similar local structure consists of a four-residue segment and three sequentially separated residues. In particular the four-residue segments of these enzymes were found to have nearly identical conformations in their backbone parts, which are involved in the recognition of adenine. This common local structure was also found in substrate-free three-dimensional structures of other proteins that are similar to DD-ligase in the chain fold and of other protein kinases. As the proteins with different folds were found to share a common local structure, these proteins seem to constitute a remarkable example of convergent evolution for the same recognition mechanism.

A comparison of a normal mode analysis and principal component analysis of a 200-ps molecular dynamics trajectory of bovine pancreatic trypsin inhibitor in vacuum has been made in order to further elucidate the harmonic and anharmonic aspects in the dynamics of proteins. An anharmonicity factor is defined which measures the degree of anharmonicity in the modes, be they principal modes or normal modes, and it is shown that the principal mode system naturally divides into anharmonic modes with peak frequencies below 80 cm-1, and harmonic modes with frequencies above this value. In general the larger the mean-square fluctuation of a principal mode, the greater the degree of anharmonicity in its motion. The anharmonic modes represent only 12% of the total number of variables, but account for 98% of the total mean-square fluctuation. The transitional nature of the anharmonic motion is demonstrated. The results strongly suggest that in a large subspace, the free energy surface, as probed by the simulation, is approximated by a multi-dimensional parabola which is just a rescaled version of the parabola corresponding to the harmonic approximation to the conformational energy surface at a single minimum. After 200 ps, the rescaling factor, termed the "normal mode rescaling factor," has apparently converged to a value whereby the mean-square fluctuation within the subspace is about twice that predicted by the normal mode analysis.

What is the significance of three-dimensional structural similarity? This fundamental question still remains unanswered in spite of advances in automatic structure comparison methods that have been made in the last few years. The answer to this question will give us a much deeper insight into the principles of protein architecture.

We have developed a new method of protein structure comparison based on spatial arrangements of secondary structural elements (SSEs). Each SSE is represented by a single vector, and common spatial arrangements of vectors in a pair of proteins are detected. The method allows not only insertions and deletions of SSEs, but also topological permutations. It has a flexible target function that can be adjusted depending upon particular levels or definitions of structural similarity, and it is fast enough to allow structural comparisons for many pairs of proteins. The parameters for the target function are determined based on distributions of the geometrical variables for the spatial arrangements of the equivalent SSEs in well-known structural motifs. The obtained parameter set is tuned for detecting relatively strong structural similarity. We report several tests on examples including comparisons of known structural similarity and database searches for a target structure, and examine the results when this parameter set is used for the comparison of distantly related structures.

We have completed an exhaustive search for the common spatial arrangements of backbone fragments (SARFs) in nonhomologous proteins. This type of local structural similarity, incorporating short fragments of backbone atoms, arranged not necessarily in the same order along the polypeptide chain, appears to be important for protein function and stability. To estimate the statistical significance of the similarities, we have introduced a similarity score. We present several locally similar structures, with a large similarity score, which have not yet been reported. On the basis of the results of pairwise comparison, we have performed hierarchical cluster analysis of protein structures. Our analysis is not limited by comparison of single chains but also includes complex molecules consisting of several subunits. The SARFs with backbone fragments from different polypeptide chains provide a stable interaction between subunits in protein molecules. In many cases the active site of enzyme is located at the same position relative to the common SARFs, implying a function of the certain SARFs as a universal interface of the protein-substrate interaction.

Two molecular dynamics simulations on bovine pancreatic trypsin inhibitor (BPTI), have been made, one in vacuum, the other in water, in order to assess the effect of the solvent water on collective motions. Principal component analysis has been performed to determine collective modes, the principal components, which are assumed to behave as effectively independent harmonic oscillators. Projection of the protein's motion in water onto the plane defined by the first two principal components shows a clustering effect in the trajectory, absent in the vacuum trajectory. This is thought to be due to many local minima in the free energy surface caused by solute-solvent interactions. In order to assess the viscous effect of the solvent, friction coefficients for the principal components were determined by analyzing their velocity correlation functions in terms of the Langevin equation for an independent damped oscillator. Consistent with this analysis is that all modes have friction coefficients centered on the value of 47 cm-1 in a range of +/- 10 cm-1. With this friction coefficient, all modes of effective frequencies below 23.5 cm-1 display overdamped motion. By assuming the harmonic approximation for the conformational energy surface for BPTI in vacuum to be valid for BPTI in water, and treating each mode as an independent damped oscillator with a friction coefficient of 47 cm-1, the shift to higher frequencies in the water spectrum relative to the vacuum spectrum could be almost exactly reproduced, indicating this shift is due solely to the viscous effect of the solvent. By analyzing the time correlation functions of the first four principal components it is found that they can be very well described as independent damped oscillators each with a friction coefficient of 47 cm-1.

An x-ray crystallographic refinement method, referred to as the normal mode refinement, is proposed. The Debye-Waller factor is expanded in terms of the effective normal modes whose amplitudes and eigenvectors are experimentally determined by the crystallographic refinement. In contrast to the conventional method, the atomic motions are treated generally as anisotropic and concerted. This method is assessed by using the simulated x-ray data given by a Monte Carlo simulation of human lysozyme. In this article, we refine the dynamic structure by fixing the average static structure to exact coordinates. It is found that the normal mode refinement, using a smaller number of variables, gives a better R factor and more information on the dynamics (anisotropy and collectivity in the motion).

A long record of computer simulation of folding-unfolding transition in a two-dimensional lattice model of protein as monitored by one conformational order parameter was studied to see if it could be approximated by a markoffian process. For this purpose the normalized time correlation functions of the order parameter were calculated (i) directly from the record of simulation and (ii) by assuming the markoffian behavior of the record. Both of them can be well approximated by a sum of two simple exponential terms. The relaxation time of the slow relaxing term, which corresponds to the overall folding-unfolding transition, becomes very short when the markoffian assumption is made. From this observation we conclude that intermediate states defined by one more-or-less arbitrarily chosen conformational parameter are, in general, collections of very heterogeneous conformations and therefore transitions between them cannot be markoffian. This indicates the importance of multi-parameter observation of dynamic process of folding. Characteristic features of the method of trapping disulfide intermediates are discussed.

The theoretical model of proteins on the two-dimensional square lattice, introduced previously, is extended to include the hydrophobic interactions. Two proteins, whose native conformations have different folded patterns, are studied. Units in the protein chains are classified into polar units and nonpolar units. If there is a vacant lattice point next to a nonpolar unit, it is interpreted as being occupied by solvent water and the entropy of the system is assumed to decrease by a certain amount. Besides these hydrophobic free energies, the specific long-range interactions studied in previous papers are assumed to be operative in a protein chain. Equilibrium properties of the folding and unfolding transitions of the two proteins are found to be similar, even though one of them was predicted, based on the one globule model of the transitions, to unfold through a significant intermediate state (or at least to show a tendency toward such a behavior), when the hydrophobic interactions are strongly weighted. The failure of this prediction led to the development of a more refined model of transitions; a non-interacting local structure model. The hydrophobic interactions assumed here have a character of non-specific long-range interactions. Because of this character the hydrophobic interactions have the effect of decelerating the folding kinetics. The deceleration effect is less pronounced in one of the two proteins, whose native conformation is stabilized by many pairs of medium-range interactions. It is therefore inferred that the medium-range interactions have the power to cope with the decelerating effect of the non-specific hydrophobic interactions.

How specific or definite are pathways of folding and unfolding in globular proteins? In order to study this question, computer simulation of the folding-unfolding transition was carried out in a two-dimensional lattice model of proteins in which it is assumed that strongly specific intramolecular interactions contribute to the stability of the native conformation. This specificity of the interactions should tend to make the pathways of folding and unfolding more definite than reality. Yet, the analysis of the record of simulation indicated the process of transition to be stochastic rather than definite. This poses a fundamental problem of how to describe the pathways of folding and unfolding transition. It is argued that the description should consist of (i) a definition of intermediate states in terms of characteristic conformational features and (ii) stochastic rules of transitions between these intermediate states. The simplest would be the case in which the transitions occur as a markoffian process.

Globular proteins in the native state are assumed to behave as continuous elastic spheres in the low frequency breathing motions. Reasonable values of Young's modulus E = 10(11) dyne/cm2 and the radius of the sphere ro = 20 A, yield a wave number of 26 cm-1 for the fundamental vibration of the sphere. The peak at around 30 cm-1 in the laser Raman spectra of native alpha-chymotrypsin and pepsin observed by Brown et al. might be assigned to the breathing motion which the native proteins undergo as continuous elastic bodies.

It is pointed out that various energy terms contributing to stabilize the native state of globular proteins are consistent in the first approximation with each other in the native state. This means that each energy term is individually minimized at the minimum point of the total energy. I proposed (1) to call this fact "the consistency principle in protein structure." The fair success of various methods of prediction of the secondary structures in globular proteins from their amino acid sequence is often interpreted as indicating the dominance of the short-range interactions in determining the local structures of the polypeptide chains. Partly from such a point of view, the hierarchic condensation model has been popular for the process of protein folding. However the consistency principle indicates that the short-range interactions are just one type of intramolecular interaction which contributes to stabilization of the native structure together with other mutually consistent types of intramolecular interactions. Therefore the hierarchic condensation model is not necessarily a unique model of protein folding. Roles of a possible nonspecific globular state, stabilized by nonspecific long-range intramolecular interactions, in the folding process are discussed. It is expected that this nonspecific globular state is observed either as an equilibrium or a kinetic intermediate state between the unfolded and the folded native states. Observation as a kinetic intermediate state is expected to occur in experiments done under strongly refolding conditions. In this case the polypeptide chain in the unfolded state collapses into a nonspecific globule by the action of nonspecific long-range intramolecular interactions. Two possible mechanisms of the transition from the nonspecific globular state to the specific native folded state are discussed. In an experiment done under weakly refolding conditions, folding is expected to occur according to the embryo-nucleus model. This model is a refined version of the hierarchic condensation model. Refinement is done by taking into account the fact that the intermediate structures assumed in the hierarchic condensation model are unstable against both the native folded state and the unfolded state. A nucleus is an ordered structure of a certain size. Ordered structures of a size larger than a nucleus tend to fold further to become the native specific globule. Ordered structures of a size smaller than a nucleus tend to unfold. Embryos are intrinsically unstable ordered structures smaller than a nucleus. Folding occurs when embryos grow in size to become a nucleus.(ABSTRACT TRUNCATED AT 400 WORDS)

Multicopper blue proteins (MCBPs) are multidomain proteins that utilize the distinctive redox ability of copper ions. There are a variety of MCBPs that have been roughly classified into three different groups, based on their domain organization and functions: (i) nitrite reductase-type with two domains, (ii) laccase-type with three domains, and (iii) ceruloplasmin-type with six domains. Together, the second and third group are often commonly called multicopper oxidases (MCOs). The rapid accumulation of genome sequence information in recent years has revealed several new types of proteins containing MCBP domains, mainly from bacteria. In this review, the recent research on the functions and structures of MCBPs is summarized, mainly focusing on the new types. The latter half of this review focusses on the twodomain MCBPs, which we propose as the evolutionary intermediate of the MCBP family. [ABSTRACT FROM AUTHOR]