DOT2: Macromolecular docking with improved biophysical models.

Computational docking is a useful tool for predicting macromolecular complexes, which are often difficult to determine experimentally. Here, we present the DOT2 software suite, an updated version of the DOT intermolecular docking program. DOT2 provides straightforward, automated construction of improved biophysical models based on molecular coordinates, offering checkpoints that guide the user to include critical features. DOT has been updated to run more quickly, allow flexibility in grid size and spacing, and generate an infinitive complete list of favorable candidate configurations. Output can be filtered by experimental data and rescored by the sum of electrostatic and atomic desolvation energies. We show that this rescoring method improves the ranking of correct complexes for a wide range of macromolecular interactions and demonstrate that biologically relevant models are essential for biologically relevant results. The flexibility and versatility of DOT2 accommodate realistic models of complex biological systems, improving the likelihood of a successful docking outcome.

Constitutive activation of extracellular signal-regulated kinase 2 by synergistic point mutations.

Constitutively active mutant forms of signaling enzymes provide insight into mechanisms of activation as well as useful molecular tools for probing downstream targets. In this study, point mutations in ERK2 at conserved residues L73P and S151D were identified that individually led to 8-12-fold increased specific activity and in combination reached 50-fold, indicating synergistic interactions between these residues. Examination by mass spectrometry, phosphatase sensitivity, and Western blotting revealed that the mutations enhanced ERK2 activity by facilitating intramolecular autophosphorylation predominantly at Tyr-185 and to a lesser extent at Thr-183 and that phosphorylation at both sites is required for activation. A set of short molecular dynamics simulations were carried out using different random seeds to sample locally accessible configurations. Simulations of the active mutant showed potential hydrogen bonding interactions between the phosphoryl acceptor and catalytic nucleophile, which could account for enhanced intramolecular autophosphorylation. In intact cells, the ERK2 mutants were functionally active in phosphorylating Elk-1 and RSK1 and activating the c-fos promoter. This activity was only partially reduced upon treatment of cells with the MKK1/2 inhibitor, U0126, indicating that in vivo the mechanism of ERK2 activation occurs substantially through autophosphorylation and partially through phosphorylation by MKK1/2.

Rapid atomic density methods for molecular shape characterization.

Two methods for rapid characterization of molecular shape are presented. Both techniques are based on the density of atoms near the molecular surface. The Fast Atomic Density Evaluation (FADE) algorithm uses fast Fourier transforms to quickly estimate densities. The Pairwise Atomic Density Reverse Engineering (PADRE) method derives modified density measures from the relationship between atomic density and total potentials. While many shape-characterization techniques define shape relative to a surface, the descriptors returned by FADE and PADRE can measure local geometry from points within the three-dimensional space surrounding a molecule. The methods can be used to find crevices and protrusions near the surface of a molecule and to test shape complementarity at the interface between docking molecules.

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.

Eigensystem analysis of the refinement of a small metalloprotein.

The eigenvalues and eigenvectors of the least-squares normal matrix for the full-matrix refinement problem contain a great deal of information about the quality of a model; in particular the precision of the model parameters and correlations between those parameters. They also allow the isolation of those parameters or combinations of parameters which are not determined by the available data. Since a protein refinement is usually under-determined without the application of geometric restraints, such indicators of the reliability of a model offer an important contribution to structural knowledge. Eigensystem analysis is applied to the normal matrices for the refinement of a small metalloprotein using two data sets and models determined at different resolutions. The eigenvalue spectra reveal considerable information about the conditioning of the problem as the resolution varies. In the case of a restrained refinement, it also provides information about the impact of various restraints on the refinement. Initial results support conclusions drawn from the free R factor. Examination of the eigenvectors provides information about which regions of the model are poorly determined. In the case of a restrained refinement, it is also possible to isolate places where X-ray and geometric restraints are in disagreement, usually indicating a problem in the model.

600 ps molecular dynamics reveals stable substructures and flexible hinge points in cAMP dependent protein kinase.

Molecular dynamics simulations of the catalytic subunit of cAMP dependent protein kinase (cAPK) have been performed in an aqueous environment. The relations among the protein hydrogen-bonding network, secondary structural elements, and the internal motions of rigid domains were examined. The values of fluctuations of protein dihedral angles during dynamics show quite distinct maxima in the regions of loops and minima in the regions of alpha-helices and beta-strands. Analyses of conformation snapshots throughout the run show stable subdomains and indicate that these rigid domains are constrained during the dynamics by a stable network of hydrogen bonds. The most stable subdomain during the dynamics was in the small lobe including part of the carboxy-terminal tail. The most significant flexible region was the highly conserved glycine-rich loop between beta strands 1 and 2 in the small lobe. Many of the main chain dihedral angle changes measured in a comparison of the crystallographic structures of "open" and "closed" conformations of cAPK correspond to the highly flexible residues found during dynamics.

Crystal structure of the potent natural product inhibitor balanol in complex with the catalytic subunit of cAMP-dependent protein kinase.

Endogenous protein kinase inhibitors are essential for a wide range of physiological functions. These endogenous inhibitors may mimic peptide substrates as in the case of the heat-stable protein kinase inhibitor (PKI), or they may mimic nucleotide triphosphates. Natural product inhibitors, endogenous to the unique organisms producing them, can be potent exogenous inhibitors against foreign protein kinases. Balanol is a natural product inhibitor exhibiting low nanomolar Ki values against serine and threonine specific kinases, while being ineffective against protein tyrosine kinases. To elucidate balanol's specific inhibitory effects and provide a basis for understanding inhibition-regulated biological processes, a 2.1 A resolution crystal structure of balanol in complex with cAMP-dependent protein kinase (cAPK) was determined. The structure reveals conserved binding regions and displays extensive complementary interactions between balanol and conserved cAPK residues. This report describes the structure of a protein kinase crystallized with a natural ATP mimetic in the absence of metal ions and peptide inhibitor.

Conformation-invariant structures of the alpha1beta1 human hemoglobin dimer.

Analysis of the conformational differences between the oxy and deoxy forms of hemoglobin is complicated by shifting coordinate systems and correlated motions between different parts of the molecule. Methods independent of any frame of reference were used to study the differences in structure between the oxy and deoxy forms of the human hemoglobin alphabeta dimer. Differences between the deoxy and oxy dimer structures can be characterized as rearrangements of 15 substructures persisting between the two conformations. Such substructures are of two kinds, either rigid domains or tertiary substructures. Rigid domains are groups of residues for which all inter-residue distances are conformationally invariant. Residues belonging to a rigid domain do not have to be spatially contiguous nor must they have consecutive sequence numbers. The largest such substructure is a rigid core that spans both the alpha and beta monomers and includes 44% of the dimer. Other rigid domains exist within the heme pockets. An alternative but closely related view of the molecule is based on tertiary substructures. Unlike a rigid domain, a tertiary substructure must have consecutively numbered residues and the residue that ends one tertiary substructure begins the next. The decomposition of the dimer into tertiary substructures represents the dimer as a framework of connected stiff structural elements. Viewed as a set of tertiary substructures, the hemoglobin dimer has the same three principal functional elements: the dimer core and the alpha and beta heme pockets, with the heme pockets held to the dimer core by CD and FG corners. The tertiary substructures that comprise the dimer core include 51% of the molecule. When ligands bind at the hemes, the FG corners communicate structural changes in the hemes to the dimer cores, which may mediate heme-heme cooperativity.

A binary complex of the catalytic subunit of cAMP-dependent protein kinase and adenosine further defines conformational flexibility.

BACKGROUND: cAMP-dependent protein kinase (cAPK), a ubiquitous protein in eukaryotic cells, is one of the simplest members of the protein kinase family. It was the first protein kinase to be crystallized and continues to serve as a biochemical and structural prototype for this family of enzymes. To further understand the conformational changes that occur in different liganded and unliganded states of cAPK, the catalytic subunit of cAPK was crystallized in the absence of peptide inhibitor. RESULTS: The crystal structure of the catalytic subunit of mouse recombinant cAPK (rC) complexed with adenosine was solved at 2.6 A resolution and refined to a crystallographic R factor of 21.9% with good stereochemical parameters. This is the first structure of the rC subunit that lacks a bound inhibitor or substrate peptide. The structure was solved by molecular replacement and comprises two lobes (large and small) which contain a number of conserved loops. CONCLUSIONS: The binary complex of rC and adenosine adopts an 'intermediate' conformation relative to the previously described 'closed' and 'open' conformations of other rC complexes. Based on a comparison of these structures, the induced fit that is necessary for catalysis and closing of the active-site cleft appears to be confined to the small lobe, as in the absence of the peptide the conformation of the large lobe, including the peptide-docking surface, does not change. Three specific components contribute to the closing of the cleft: rotation of the small lobe; movement of the C-terminal tail; and closing of the so-called glycine-rich loop. There is no induced fit in the large lobe to accommodate the peptide and the closing of the cleft. A portion of the C-terminal tail, residues 315-334, serves as a gate for the entry or exit of the nucleotide into the hydrophobic active-site cleft.

Catalytic subunit of cAMP-dependent protein kinase: electrostatic features and peptide recognition.

The electrostatic field was calculated for the mammalian cAMP-dependent protein kinase (PKA) catalytic subunit (C-subunit) complexed with a 20-residue peptide from a heat stable protein kinase inhibitor (PKI: 5-24). The electrostatic field was also calculated for the C-subunit complexed with a modeled heptapeptide substrate that has been used extensively in structure/function studies for the C-subunit. Perturbations in the electrostatic free energy were calculated when single ionizable active site residues were mutated to alanine. These perturbations in electrostatic free energy were correlated to changes in the binding energy measured in a charge-to-alanine scan of the homologous yeast C-subunit by M. J. Zoller and C. S. Gibbs [(1991) Journal of Biological Chemistry, Vol. 266, pp. 8923-8931; C. S. Gibbs and M. J. Zoller (1991) Biochemistry, Vol. 30, p. 22]. This analysis indicated that the substrate binding parameters primarily depend on electrostatic interactions between a substrate or inhibitor and the C-subunit. Amino acid replacements that led to large perturbations in the electrostatic field are listed in the text. pKa shifts were also calculated for the substrate's phosphate accepting atom, the serine hydroxyl oxygen, when the active site ionizable residues were changed to structurally similar uncharged amino acids. The theoretical mutation of three active site residues caused large shifts in this parameter: E91Q, D166N, and D184N. The calculated pKa shifts for these mutants indicate that the rate of phosphotransfer should be markedly reduced in these cases. This prediction has been experimentally confirmed for the D166N mutant. The correlation between calculated electrostatic free energy changes and measured binding energy, and pKa shifts with phosphotransfer for C-subunit mutants were within experimental error of the measurements. The calculations of electrostatic energy and delta pKa have identified previously unconsidered active site residues in the mammalian C-subunit that contribute to binding energy and phosphotransfer.

Rigid domains in proteins: an algorithmic approach to their identification.

A rigid domain, defined here as a tertiary structure common to two or more different protein conformations, can be identified numerically from atomic coordinates by finding sets of residues, one in each conformation, such that the distance between any two residues within the set belonging to one conformation is the same as the distance between the two structurally equivalent residues within the set belonging to any other conformation. The distance between two residues is taken to be the distance between their respective alpha carbon atoms. With the methods of this paper we have found in the deoxy and oxy conformations of the human hemoglobin alpha 1 beta 1 dimer a rigid domain closely related to that previously identified by Baldwin and Chothia (J. Mol. Biol. 129: 175-220, 1979). We provide two algorithms, both using the difference-distance matrix, with which to search for rigid domains directly from atomic coordinates. The first finds all rigid domains in a protein but has storage and processing demands that become prohibitively large with increasing protein size. The second, although not necessarily finding every rigid domain, is computationally tractable for proteins of any size. Because of its efficiency we are able to search protein conformations recursively for groups of non-intersecting domains. Different protein conformations, when aligned by superimposing their respective domain structures, can be examined for structural differences in regions complementing a rigid domain.

cAMP-dependent protein kinase: crystallographic insights into substrate recognition and phosphotransfer.

The crystal structure of ternary and binary substrate complexes of the catalytic subunit of cAMP-dependent protein kinase has been refined at 2.2 and 2.25 A resolution, respectively. The ternary complex contains ADP and a 20-residue substrate peptide, whereas the binary complex contains the phosphorylated substrate peptide. These 2 structures were refined to crystallographic R-factors of 17.5 and 18.1%, respectively. In the ternary complex, the hydroxyl oxygen OG of the serine at the P-site is 2.7 A from the OD1 atom of Asp 166. This is the first crystallographic evidence showing the direct interaction of this invariant carboxylate with a peptide substrate, and supports the predicted role of Asp 166 as a catalytic base and as an agent to position the serine -OH for nucleophilic attack. A comparison of the substrate and inhibitor ternary complexes places the hydroxyl oxygen of the serine 2.7 A from the gamma-phosphate of ATP and supports a direct in-line mechanism for phosphotransfer. In the binary complex, the phosphate on the Ser interacts directly with the epsilon N of Lys 168, another conserved residue. In the ternary complex containing ATP and the inhibitor peptide, Lys 168 interacts electrostatically with the gamma-phosphate of ATP (Zheng J, Knighton DR, Ten Eyck LF, Karlsson R, Xuong NH, Taylor SS, Sowadski JM, 1993, Biochemistry 32:2154-2161). Thus, Lys 168 remains closely associated with the phosphate in both complexes. A comparison of this binary complex structure with the recently solved structure of the ternary complex containing ATP and inhibitor peptide also reveals that the phosphate atom traverses a distance of about 1.5 A following nucleophilic attack by serine and transfer to the peptide. No major conformational changes of active site residues are seen when the substrate and product complexes are compared, although the binary complex with the phosphopeptide reveals localized changes in conformation in the region corresponding to the glycine-rich loop. The high B-factors for this loop support the conclusion that this structural motif is a highly mobile segment of the protein.

A conserved helix motif complements the protein kinase core.

Residues 40-300 of the mammalian catalytic (C) subunit of cAMP-dependent protein kinase define a conserved bilobal catalytic core shared by all eukaryotic protein kinases. Contiguous to the core is an extended amphipathic alpha-helix (A helix). Trp30, a prominent feature of this helix, fills a deep hydrophobic pocket between the two lobes on the surface opposite to the active site. The C subunit in Dictyostelium discoideum shows sequence conservation of residues 40-350 with the mouse enzyme but contains an N-terminal extension of 332 residues. A sequence corresponding to the A helix contiguous to the core is absent. However, we have now identified a remote A-helix motif (residues 77-98). When the core of the Dictyostelium C subunit was modeled, based on the mouse C subunit, complementarity between this putative A helix and the surface of the core was found to be conserved. Analysis of other protein kinases reveals that the A-helix motif is not restricted to cAMP-dependent protein kinase. In the Src-related family of protein kinases, for example, an A helix is very likely contiguous to the core, thus serving as a linker between the conserved catalytic core and the Src homology 2 domain. We predict that an A-helix motif complementary to the core will be a conserved feature of most eukaryotic protein kinases.

Crystal structures of the myristylated catalytic subunit of cAMP-dependent protein kinase reveal open and closed conformations.

Three crystal structures, representing two distinct conformational states, of the mammalian catalytic subunit of cAMP-dependent protein kinase were solved using molecular replacement methods starting from the refined structure of the recombinant catalytic subunit ternary complex (Zheng, J., et al., 1993a, Biochemistry 32, 2154-2161). These structures correspond to the free apoenzyme, a binary complex with an iodinated inhibitor peptide, and a ternary complex with both ATP and the unmodified inhibitor peptide. The apoenzyme and the binary complex crystallized in an open conformation, whereas the ternary complex crystallized in a closed conformation similar to the ternary complex of the recombinant enzyme. The model of the binary complex, refined at 2.9 A resolution, shows the conformational changes associated with the open conformation. These can be described by a rotation of the small lobe and a displacement of the C-terminal 30 residues. This rotation of the small lobe alters the cleft interface in the active-site region surrounding the glycine-rich loop and Thr 197, a critical phosphorylation site. In addition to the conformational changes, the myristylation site, absent in the recombinant enzyme, was clearly defined in the binary complex. The myristic acid binds in a deep hydrophobic pocket formed by four segments of the protein that are widely dispersed in the linear sequence. The N-terminal 40 residues that lie outside the conserved catalytic core are anchored by the N-terminal myristylate plus an amphipathic helix that spans both lobes and is capped by Trp 30. Both posttranslational modifications, phosphorylation and myristylation, contribute directly to the stable structure of this enzyme.

Structural features that specify tyrosine kinase activity deduced from homology modeling of the epidermal growth factor receptor.

To identify structural features that distinguish protein-tyrosine kinases from protein-serine kinases, a molecular model of the kinase domain of epidermal growth factor receptor was constructed by substituting its amino acid sequence for the amino acid sequence of the catalytic subunit of cAMP-dependent protein kinase in a 2.7-A refined crystallographic model. General folding was conserved as was the configuration of invariant residues at the active site. Two sequence motifs that distinguish the two families correspond to loops that converge at the active site of the enzyme. A conserved arginine in the catalytic loop is proposed to interact with the gamma phosphate of ATP. The second loop provides a binding surface that positions the tyrosine of the substrate. A positively charged surface provides additional sites for substrate recognition.

cAMP-dependent protein kinase defines a family of enzymes.

The structure of the recombinant mouse catalytic subunit of cAMP-dependent protein kinase is reviewed with particular emphasis on the overall features and specific amino acids that are shared by all members of the eukaryotic protein kinase family. The crystal structure of a ternary complex containing both MgATP and a twenty-residue inhibitor peptide defines the precise role of the conserved residues that are clustered at the active site. In addition to catalysing the post-translational modification of other proteins, the catalytic subunit is itself subject to covalent modifications. It is a phosphoprotein and is also myristylated at its amino terminus. The enzyme when crystallized in the presence of detergent shows a detergent molecule bound to an acyl pocket that is presumably occupied by the myristyl moiety in the mammalian enzyme. When expressed in E. coli, the catalytic subunit is autophosphorylated at four sites. Two stable phosphates at Ser338 and Thr197 interact with multiple protein side chains thus explaining why they are inaccessible to phosphatases. Although all substrates and inhibitors of the catalytic subunit share a general minimum consensus sequence, the high affinity binding of protein inhibitors such as the regulatory subunits and the heat stable protein kinase inhibitors require additional determinants that lie beyond the consensus site. These two physiological inhibitors of the catalytic subunit appear to use different sites to achieve high-affinity binding.

2.0 A refined crystal structure of the catalytic subunit of cAMP-dependent protein kinase complexed with a peptide inhibitor and detergent.

. A mutant (Serl39Ala) of the mouse recombinant catalytic (C) subunit of cAMP-dependent protein kinase was co-crystallized with a peptide inhibitor, PKI(5-24), and MEGA-8 (octanoyl-N-methylglucamide) detergent. This structure was refined using all observed data (30 248 reflections) between 30 and 1.95 A resolution to an R factor of 0.186. R.m.s. deviations of bond lengths and bond angles are 0.013 A and 2.3 degrees, respectively. The final model has 3075 atoms (207 solvent) with a mean B factor of 31.9 A(2). The placement of invariant protein-kinase residues and most C:PKI(5-24) interactions were confirmed, but register errors affecting residues 55-64 and 309-339 were corrected during refinement by shifting the affected sequences toward the C terminus along the previously determined backbone path. New details of C:PKI(5-24) interactions and the Ser338 autophosphorylation site are described, and the acyl group binding site near the catalytic subunit NH(2) terminus is identified.

2.2 A refined crystal structure of the catalytic subunit of cAMP-dependent protein kinase complexed with MnATP and a peptide inhibitor.

. The crystal structure of a ternary complex containing the catalytic subunit of cAMP-dependent protein kinase, ATP and a 20-residue inhibitor peptide was refined at a resolution of 2.2 A to an R value of 0.177. In order to identify the metal binding sites, the crystals, originally grown in the presence of low concentrations of Mg(2+), were soaked in Mn(2+). Two Mn(2+) ions were identified using an anomalous Fourier map. One Mn(2+) ion bridges the gamma- and beta-phosphates and interacts with Asp184 and two water molecules. The second Mn(2+) ion interacts with the side chains of Asn171 and Asp l84 as well as with a water molecule. Modeling a serine into the P site of the inhibitor peptide suggests a mechanism for phosphotransfer.