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.

Structure of a peptide inhibitor bound to the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinase.

The structure of a 20-amino acid peptide inhibitor bound to the catalytic subunit of cyclic AMP-dependent protein kinase, and its interactions with the enzyme, are described. The x-ray crystal structure of the complex is the basis of the analysis. The peptide inhibitor, derived from a naturally occurring heat-stable protein kinase inhibitor, contains an amphipathic helix that is followed by a turn and an extended conformation. The extended region occupies the cleft between the two lobes of the enzyme and contains a five-residue consensus recognition sequence common to all substrates and peptide inhibitors of the catalytic subunit. The helical portion of the peptide binds to a hydrophobic groove and conveys high affinity binding. Loops from both domains converge at the active site and contribute to a network of conserved residues at the sites of magnesium adenosine triphosphate binding and catalysis. Amino acids associated with peptide recognition, nonconserved, extend over a large surface area.

Crystal structure of the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinase.

The crystal structure of the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinase complexed with a 20-amino acid substrate analog inhibitor has been solved and partially refined at 2.7 A resolution to an R factor of 0.212. The magnesium adenosine triphosphate (MgATP) binding site was located by difference Fourier synthesis. The enzyme structure is bilobal with a deep cleft between the lobes. The cleft is filled by MgATP and a portion of the inhibitor peptide. The smaller lobe, consisting mostly of amino-terminal sequence, is associated with nucleotide binding, and its largely antiparallel beta sheet architecture constitutes an unusual nucleotide binding motif. The larger lobe is dominated by helical structure with a single beta sheet at the domain interface. This lobe is primarily involved in peptide binding and catalysis. Residues 40 through 280 constitute a conserved catalytic core that is shared by more than 100 protein kinases. Most of the invariant amino acids in this conserved catalytic core are clustered at the sites of nucleotide binding and catalysis.

Structural basis of the intrasteric regulation of myosin light chain kinases.

The smooth muscle myosin light chain kinase (smMLCK) catalytic core was modeled by using the crystallographic coordinates of the cyclic AMP-dependent protein kinase catalytic subunit (cAPK) and a bound pseudosubstrate inhibitor peptide, PKI(5-24). Despite only 30% identity in amino acid sequence, the MLCK sequence can be readily accommodated in this structure. With the exception of the short B-helix, all major elements of secondary structure in the core are very likely conserved. The active site of the modeled MLCK complements the known requirements for peptide substrate recognition. MLCK contains a pseudosubstrate sequence that overlaps the calmodulin binding domain and has been proposed to act as an intrasteric inhibitor and occupy the substrate binding site in the absence of Ca(2+)-calmodulin. The pseudosubstrate sequence can be modeled easily into the entire backbone of PKI(5-24). The results demonstrate that the intrasteric model for regulation of MLCK by intramolecular competitive inhibition is structurally plausible.

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.

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.

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.

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.

Fast Fourier transform calculation of electron density maps.

This chapter has described the mathematical basis of the fast Fourier transform as applied to the calculation of crystallographic Fourier syntheses. The relationship between real space and reciprocal space symmetry operators has been described. Finally, program organizations have been presented for performing general crystallographic Fourier transforms on computer systems ranging from the very largest systems down to minicomputers. Programs are available from the author, written in FORTRAN IV and in Ratfor, which are suitable for building blocks in these program designs.

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.

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.

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.

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.

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.

cAMP-dependent protein kinase and the protein kinase family.

The structure of the catalytic subunit of cAMP-dependent protein kinase, the first protein kinase structure to be solved, is reviewed. The general architecture of the enzyme is described as well as the active site regions associated with substrate binding and catalysis. In particular, the unique features of the protein kinase nucleotide fold are outlined. While the catalytic subunit is one of the simplest of the protein kinases, it nevertheless serves as a structural framework for the catalytic core of the entire protein kinase family which now includes over 200 important regulatory enzymes. The essential and conserved features of this core are summarized, and a preliminary model of myosin light-chain kinase, based on the structure of the catalytic subunit, is also discussed.

Structural framework for the protein kinase family.

In this review, we have summarized the general structural features of the catalytic subunit of cAMP-dependent protein kinase, emphasizing those features that will very likely be conserved in all members of the protein kinase family. The overall secondary structure of the catalytic core will probably be conserved throughout the catalytic core, as will the active site regions associated with MgATP binding and catalysis. The mechanisms for activation and the role of protein phosphorylation are unique for each kinase. The structure of the catalytic subunit now provides a general framework for modeling other protein kinases. Although this is no substitute for a crystal structure for each protein kinase, this one structure, nevertheless, does provide major insights to the molecular organization of each of these enzymes.

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.

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.