Protein docking using continuum electrostatics and geometric fit.

The computer program DOT quickly finds low-energy docked structures for two proteins by performing a systematic search over six degrees of freedom. A novel feature of DOT is its energy function, which is the sum of both a Poisson-Boltzmann electrostatic energy and a van der Waals energy, each represented as a grid-based correlation function. DOT evaluates the energy of interaction for many orientations of the moving molecule and maintains separate lists scored by either the electrostatic energy, the van der Waals energy or the composite sum of both. The free energy is obtained by summing the Boltzmann factor over all rotations at each grid point. Three important findings are presented. First, for a wide variety of protein-protein interactions, the composite-energy function is shown to produce larger clusters of correct answers than found by scoring with either van der Waals energy (geometric fit) or electrostatic energy alone. Second, free-energy clusters are demonstrated to be indicators of binding sites. Third, the contributions of electrostatic and attractive van der Waals energies to the total energy term appropriately reflect the nature of the various types of protein-protein interactions studied.

Definition of the interaction domain for cytochrome c on cytochrome c oxidase. III. Prediction of the docked complex by a complete, systematic search.

The electron transfer complex between bovine cytochrome c oxidase and horse cytochrome c has been predicted with the docking program DOT, which performs a complete, systematic search over all six rotational and translational degrees of freedom. Energies for over 36 billion configurations were calculated, providing a free-energy landscape showing guidance of positively charged cytochrome c to the negative region on the cytochrome c oxidase surface formed by subunit II. In a representative configuration, the solvent-exposed cytochrome c heme edge is within 4 A of the indole ring of subunit II residue Trp(104), indicating a likely electron transfer path. These two groups are surrounded by a small, hydrophobic contact region, which is surrounded by electrostatically complementary hydrophilic interactions. Cytochrome c/cytochrome c oxidase interactions of Lys(13) with Asp(119) and Lys(72) with Gln(103) and Asp(158) are the most critical polar interactions due to their proximity to the hydrophobic region and exclusion from bulk solvent. The predicted complex matches previous mutagenesis, binding, and time-resolved kinetics studies that implicate Trp(104) in electron transfer and show the importance of specific charged residues to protein affinity. Electrostatic forces not only enhance long range protein/protein association; they also predominate in short range alignment, creating the transient interaction needed for rapid turnover.

Visualizing the future of molecular graphics.

The field of computer graphics has played an important role in the advancement of structural molecular biology and in the development of structure-based drug design. This article will provide a brief background on the development of this technology, and then focus on the current trends and future directions in molecular graphics and how they will impact the practice of molecular modeling and design. Specific areas that will be covered include: 1) the development of surface and volume based representations of molecular properties and interactions; 2) new approaches to modeling flexible and multi-component structures, and 3) the impact of object-oriented graphics-based programming and the rapidly growing use of network based computing.

AVS software for visualization in molecular microscopy.

AVS (Application Visualization System) is commercially available software for analyzing and viewing data. AVS is primarily used in the physical sciences and engineering, and here we describe the application of AVS for examining three-dimensional density maps generated by electron microscopy and image processing. For this purpose, AVS can be applied with relative ease, even though the software is indeed quite sophisticated. The primary advantage is that visualization applications can be generated by combining software components, called modules, into executable flow networks. Simple networks are described for generating ribbon diagrams of macromolecules, surface-shaded views, and contour maps. Easy to use dials, bar sliders, and buttons provide tremendous versatility for real-time manipulation of isosurface values, depth cueing, view orientation, size, and animation. In addition, AVS supplies a framework for building new modules in C or FORTRAN. Modules for excavation and cropping provide tools that are particularly useful for extracting segments of a map and for examining maps of supramolecular complexes such as viruses. We describe a number of modules we have designed for analysis of three-dimensional data sets, as well as modules for importing image data from other software packages into AVS. We also describe xformat, a stand-alone file conversion utility designed to allow import of a variety of image and map file formats into AVS.

Atomic and residue hydrophilicity in the context of folded protein structures.

Water-protein interactions drive protein folding, stabilize the folded structure, and influence molecular recognition and catalysis. We analyzed the closest protein contacts of 10,837 water molecules in crystallographic structures to define a specific hydrophilicity scale reflecting specific rather than bulk solvent interactions. The tendencies of different atom and residue types to be the nearest protein neighbors of bound water molecules correlated with other hydrophobicity scales, verified the relevance of crystallographically determined water positions, and provided a direct experimental measure of water affinity in the context of the folded protein. This specific hydrophilicity was highly correlated with hydrogen-bonding capacity, and correlated better with experimental than computationally derived measures of partitioning between aqueous and organic phases. Atoms with related chemistry clustered with respect to the number of bound water molecules. Neutral and negatively charged oxygen atoms were the most hydrophilic, followed by positively-charged then neutral nitrogen atoms, followed by carbon and sulfur atoms. Agreement between observed side-chain specific hydrophilicity values and values derived from the atomic hydrophilicity scale showed that hydrophilicity values can be synthesized for different functional groups, such as unusual side or main chains, discontinuous epitopes, and drug molecules. Two methods of atomic hydrophilicity analysis provided a measure of complementarity in the interfaces of trypsin:pancreatic trypsin inhibitor and HIV protease:U-75875 inhibitor complexes.

The interdependence of protein surface topography and bound water molecules revealed by surface accessibility and fractal density measures.

To characterize water binding to proteins, which is fundamental to protein folding, stability and activity, the relationships of 10,837 bound water positions to protein surface shape and residue type were analyzed in 56 high-resolution crystallographic structures. Fractal atomic density and accessibility algorithms provided an objective characterization of deep grooves in solvent-accessible protein surfaces. These deep grooves consistently had approximately the diameter of one water molecule, suggesting that deep grooves are formed by the interactions between protein atoms and bound water molecules. Protein surface topography dominates the chemistry and extent of water binding. Protein surface area within grooves bound three times as many water molecules as non-groove surface; grooves accounted for one-quarter of the total surface area yet bound half the water molecules. Moreover, only within grooves did bound water molecules discriminate between different side-chains. In grooves, main-chain surface was as hydrated as that of the most hydrophilic side-chains, Asp and Glu, whereas outside grooves all main and side-chains bound water to a similar, and much decreased, extent. This identification of the interdependence of protein surface shape and hydration has general implications for modelling and prediction of protein surface shape, recognition, local folding and solvent binding.

Visualization of molecular flexibility and its effects on electrostatic recognition.

To study the effect of protein flexibility on electrostatic recognition, we have devised two novel computer graphic representations of the changes in the electrostatic field of a protein resulting from its internal motions. The atomic structure of Cu, Zn superoxide dismutase was minimized, and the 200 lowest frequency normal modes of the enzyme were determined. Individual and combined normal-mode vibrations were visualized interactively with the program Flex. Normal-mode motions are fast enough (approximately 10(-11) s cycle-1) to evade solvent damping, thus allowing long-range electrostatic interactions to dominate. The changing electrostatic environment of the protein was examined by animating precalculated frames of electrostatic field vectors with GRAMPS. With Vu, changes in electrostatic potential were displayed as variations in the color-coding of dots lying on a consensus surface that maintains the protein's shape. The consensus surface was calculated with the program Sphinx, and was derived from spherical harmonic approximations of expanded molecular surfaces. The ability to view the effects of molecular motions interactively should be useful in understanding the relationships of protein structure to function.