Search

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.

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.

The dynamic structure of a protein, human lysozyme, is determined by the normal mode refinement of X-ray crystal structure. This method uses the normal modes of both internal and external motions to distinguish the real internal dynamics from the external terms such as lattice disorder, and gives an anisotropic and concerted picture of atomic fluctuations. The refinement is carried out with diffraction data of 5.0 to 1.8 A resolution, which are collected on an imaging plate. The results of the refinement show: (1) Debye-Waller factor consists of two parts, highly anisotropic internal fluctuations and almost isotropic external terms. The former is smaller than the latter by a factor of 0.72 in the scale of B-factor. Therefore, the internal dynamics cannot be recognized directly from the apparent electron density distribution. (2) The internal fluctuations show basically similar features as those predicted by the normal mode analysis, with almost the same amplitude and a similar level of anisotropy. (3) Correlations of fluctuations are detected between two lobes forming the active site cleft, which move simultaneously in opposite directions. This corresponds to the hinge-bending motion of lysozyme.

A dynamic structure refinement method for X-ray crystallography, referred to as the normal mode refinement, is proposed. The Debye-Waller factor is expanded in terms of the low-frequency normal modes whose amplitudes and eigenvectors are experimentally optimized in the process of the crystallographic refinement. In this model, the atomic fluctuations are treated as anisotropic and concerted. The normal modes of the external motion (TLS model) are also introduced to cover the factors other than the internal fluctuations, such as the lattice disorder and diffusion. A program for the normal mode refinement (NM-REF) has been developed. The method has first been tested against simulated diffraction data for human lysozyme calculated by a Monte Carlo simulation. Applications of the method have demonstrated that the normal mode refinement has: (1) improved the fitting to the diffraction data, even with fewer adjustable parameters; (2) distinguished internal fluctuations from external ones; (3) determined anisotropic thermal factors; and (4) identified concerted fluctuations in the protein molecule.

We have developed a new method of detecting common spatial arrangements of backbone fragments in proteins. This method allows corresponding fragments to occur in a different order in respective amino acid sequences. We applied this method to detect structural similarities between an acid protease, endothiapepsin, and all other proteins in the protein data bank. Significant similarities were found not only with other acid proteases but also with virus proteases and with proteins having different functions. The possible biological meaning of these similarities is discussed.

Dynamic properties of deoxymyoglobin are studied theoretically by the analysis of conformational fluctuations. Root-mean-square atomic fluctuations and distance fluctuations between different segments reveal the mechanical construction of the molecule. Eight alpha-helices behave as relatively rigid bodies and corner regions are more flexible, showing larger fluctuations. More particularly, corner regions EF and GH are specific in that flanking alpha-helices extend their rigidity up to a point in the corner region and the two rigid segments are connected flexibly at that point. The FG corner is exceptional. A segment from the F helix to the beginning of the G helix, in which the FG corner is included, becomes relatively rigid by means of strong interactions with the heme group. The whole myoglobin molecule is divided into two large units of motion, one extending from the B to the E helix, and the other from the F to the H helix. These two units are connected covalently by the EF corner. However, dynamic interactions between these two units take place mainly through contacts between helices B and G and not through the EF corner. From correlation coefficients between fluctuational motions of residues and the heme group, 55 residues are identified as having strong dynamic interactions with the heme moiety. Among them, 18 residues in the three segments, one consisting of residues from the C helix to the CD corner, a second consisting of the E helix, and a third from the F helix to the beginning of the G helix, are in close contact with the heme group. Twenty-two of the 55 residues are within four residues of the 18 residues in their sequential residue number and are more than 3 A away from the heme group. The other 15 residues are located further in the sequential residue number and are all found in helices A and H. They are more than 6 A away from the heme group. By the use of correlation coefficients of fluctuations between residues, it is found that dynamic interaction with the heme group is transmitted to the A helix and the beginning of the H helix in the direction Leu(E15)----[Val(All) and Trp(A12)]. The transmission to the C-terminal end of the H helix is mediated by a long segment, from the end of the EF corner to the beginning of the G helix, that lies on the heme group and has close contacts over a wide range.(ABSTRACT TRUNCATED AT 400 WORDS)

The conformational change taking place in myoglobin concomitantly with the observed geometrical change at the heme-His(F8) linkage upon oxygenation is studied by normal mode analysis, which is based on the quadratic approximation of the conformational energy function. The heme-globin interaction energy increases for this change by 8.114 kcal/mol when both the heme group and the globin molecule are held rigid. When they are permitted flexibility, the interaction energy relaxes by 7.038 kcal/mol, and the difference (1.076 kcal/mol) is distributed as strain energy within the molecule. This increase is the work necessary for the heme group to move against the force exerted by the globin. If the force is assumed to be invariable during this move, the work is small, 0.276 kcal/mol, meaning that the force is strongly variable. Furthermore, this means that the heme group is located near the equilibrium point of the potential energy of the heme-globin interaction. The change in the localized strain energy stored in the force field at the linkage between the heme and the imidazole of HisF8 is estimated to be of the same order of magnitude as the distributed energy. The largest atomic displacements are observed at the region from the F helix to the beginning of the G helix, and secondary large displacements occur at several regions, i.e, the A helix, from the C helix to the CD corner, the E helix, and the C-terminal side of the H helix. All of these regions have strong dynamic interactions with the heme group, either directly or indirectly. Their secondary structures show complex deformations. In other parts, relatively rigid segments undergo rotational and/or bending changes in a way consistent with the large changes described above and close atomic packing within the molecule. The calculated conformational change is decomposed to vibrational normal modes of deoxymyoglobin. The magnitude of the conformational change, measured by the mass-weighted mean-square atomic displacement, is accounted for up to 92.0% by the 151 normal modes with frequencies lower than 40 cm-1. In descending order of contribution, the first six modes, each of which has a frequency lower than 12 cm-1, account for up to 57.4%. This means that the functionally important conformational change can well be expressed in terms of a relatively small number of collective low frequency normal modes.

A set of conformational restraints derived from nuclear magnetic resonance (n.m.r.) measurements on solutions of the basic pancreatic trypsin inhibitor (BPTI) was used as input for distance geometry calculations with the programs DISGEO and DISMAN. Five structures obtained with each of these algorithms were systematically compared among themselves and with the crystal structure of BPTI. It is clear that the protein architecture observed in single crystals of BPTI is largely preserved in aqueous solution, with local structural differences mainly confined to the protein surface. The results confirm that protein conformations determined in solution by combined use of n.m.r. and distance geometry are a consequence of the experimental data and do not depend significantly on the algorithm used for the structure determination. The data obtained further provide an illustration that long intramolecular distances in proteins, which are comparable with the radius of gyration, are defined with high precision by relatively imprecise nuclear Overhauser enhancement measurements of a large number of much shorter distances.

A powerful method of conformational energy minimization which uses both first and second derivatives of the energy function is applied both to a small globular protein, bovine pancreatic trypsin inhibitor (BPTI), consisting of 58 amino acid residues and to its chemical derivative obtained by carboxamidomethylation of cysteinyl residues of the 14-38 disulphide bond. Conformational fluctuations are also calculated from the second derivative matrix obtained at the respective minimum energy conformations. Appreciable conformational changes upon chemical modification are observed only in the vicinity of the site of the modification. The nuclear magnetic resonance data on both BPTI and the modified BPTI are analyzed to compare with the calculated conformational changes upon chemical modification. Good correlations are found between the theoretically predicted and experimentally deduced conformational changes. The theoretical method employed here has a general application for the calculation of small conformational changes of globular proteins upon their chemical modification or an amino acid substitution.

We have developed a method to determine the three-dimensional structure of a protein molecule from such a set of distance constraints as can be determined by nuclear magnetic resonance studies. The currently popular methods for distance geometry based on the use of the metric matrix are applicable only to small systems. The method developed here is applicable to large molecules, such as proteins, with all atoms treated explicitly. This method works in the space of variable dihedral angles and determines a three-dimensional structure by minimization of a target function. We avoid difficulties hitherto inherent in this type of approach by two new devices: the use of variable target functions; and a method of rapid calculation of the gradient of the target functions. The method is applied to the determination of the structures of a small globular protein, bovine pancreatic trypsin inhibitor, from several artificial sets of distance constraints extracted from the X-ray crystal structure of this molecule. When a good set of constraints was available for both short- and long-range distances, the crystal structure was regenerated nearly exactly. When some ambiguities, such as those expected in experimental information, are allowed, the protein conformation can be determined up to a few local deformations. These ambiguities are mainly associated with the low resolving power of the short-range information.