Accurate, conformation-dependent predictions of solvent effects on protein ionization constants

Predicting how aqueous solvent modulates the conformational transitions and influences the pKa values that regulate the biological functions of biomolecules remains an unsolved challenge. To address this problem, we developed FDPB_MF, a rotamer repacking method that exhaustively samples side chain conformational space and rigorously calculates multibody protein-solvent interactions. FDPB_MF predicts the effects on pKa values of various solvent exposures, large ionic strength variations, strong energetic couplings, structural reorganizations and sequence mutations. The method achieves high accuracy, with root mean square deviations within 0.3 pH unit of the experimental values measured for turkey ovomucoid third domain, hen lysozyme, Bacillus circulans xylanase, and human and Escherichia coli thioredoxins. FDPB_MF provides a faithful, quantitative assessment of electrostatic interactions in biological macromolecules.

Reverse engineering the (beta/alpha )8 barrel fold

The (beta/alpha)(8) barrel is the most commonly occurring fold among protein catalysts. To lay a groundwork for engineering novel barrel proteins, we investigated the amino acid sequence restrictions at 182 structural positions of the prototypical (beta/alpha)(8) barrel enzyme triosephosphate isomerase. Using combinatorial mutagenesis and functional selection, we find that turn sequences, alpha-helix capping and stop motifs, and residues that pack the interface between beta-strands and alpha-helices are highly mutable. Conversely, any mutation of residues in the central core of the beta-barrel, beta-strand stop motifs, and a single buried salt bridge between amino acids R189 and D227 substantially reduces catalytic activity. Four positions are effectively immutable: conservative single substitutions at these four positions prevent the mutant protein from complementing a triosephosphate isomerase knockout in Escherichia coli. At 142 of the 182 positions, mutation to at least one amino acid of a seven-letter amino acid alphabet produces a triosephosphate isomerase with wild-type activity. Consequently, it seems likely that (beta/alpha)(8) barrel structures can be encoded with a subset of the 20 amino acids. Such simplification would greatly decrease the computational burden of (beta/alpha)(8) barrel design.

Modular enzymes

Although modular macromolecular devices are encountered frequently in a variety of biological situations, their occurrence in biocatalysis has not been widely appreciated. Three general classes of modular biocatalysts can be identified: enzymes in which catalysis and substrate specificity are separable, multisubstrate enzymes in which binding sites for individual substrates are modular, and multienzyme systems that can catalyse programmable metabolic pathways. In the postgenomic era, the discovery of such systems can be expected to have a significant impact on the role of enzymes in synthetic and process chemistry.

Tanford-Kirkwood electrostatics for protein modeling

Solvent plays a significant role in determining the electrostatic potential energy of proteins, most notably through its favorable interactions with charged residues and its screening of electrostatic interactions. These energetic contributions are frequently ignored in computational protein design and protein modeling methodologies because they are difficult to evaluate rapidly and accurately. To address this deficiency, we report a revised form of the original Tanford-Kirkwood continuum electrostatic model [Tanford, C. & Kirkwood, J. G. (1957) J. Am. Chem. Soc. 79, 5333-5339], which accounts for the effects of solvent polarization on charged atoms in proteins. The Tanford-Kirkwood model was modified to increase its speed and to improve its sensitivity to the details of protein structure. For the 37 electrostatic self-energies of the polar side-chains in bovine pancreatic trypsin inhibitor, and their 666 interaction energies, the modified Tanford-Kirkwood potential of mean force differs from a computationally intensive numerical potential (DelPhi) by root-mean-square errors of 0.6 kcal/mol and 0.08 kcal/mol, respectively. The Tanford-Kirkwood approach makes possible a realistic treatment of electrostatics in computationally demanding protein modeling calculations. For example, pH titration calculations for ovomucoid third domain that model polar side-chain relaxation (including >2 x 10(23) rotamer conformations of the protein) provide pKa values of unprecedented accuracy.

High-resolution protein design with backbone freedom

Recent advances in computational techniques have allowed the design of precise side-chain packing in proteins with predetermined, naturally occurring backbone structures. Because these methods do not model protein main-chain flexibility, they lack the breadth to explore novel backbone conformations. Here the de novo design of a family of alpha-helical bundle proteins with a right-handed superhelical twist is described. In the design, the overall protein fold was specified by hydrophobic-polar residue patterning, whereas the bundle oligomerization state, detailed main-chain conformation, and interior side-chain rotamers were engineered by computational enumerations of packing in alternate backbone structures. Main-chain flexibility was incorporated through an algebraic parameterization of the backbone. The designed peptides form alpha-helical dimers, trimers, and tetramers in accord with the design goals. The crystal structure of the tetramer matches the designed structure in atomic detail.

Repacking protein cores with backbone freedom: structure prediction for coiled coils

Progress in homology modeling and protein design has generated considerable interest in methods for predicting side-chain packing in the hydrophobic cores of proteins. Present techniques are not practically useful, however, because they are unable to model protein main-chain flexibility. Parameterization of backbone motions may represent a general and efficient method to incorporate backbone relaxation into such fixed main-chain models. To test this notion, we introduce a method for treating explicitly the backbone motions of alpha-helical bundles based on an algebraic parameterization proposed by Francis Crick in 1953 [Crick, F. H. C. (1953) Acta Crystallogr. 6, 685-689]. Given only the core amino acid sequence, a simple calculation can rapidly reproduce the crystallographic main-chain and core side-chain structures of three coiled coils (one dimer, one trimer, and one tetramer) to within 0.6-A root-mean-square deviations. The speed of the predictive method [approximately 3 min per rotamer choice on a Silicon Graphics (Mountain View, CA) 4D/35 computer] permits it to be used as a design tool.

Crystal structure of an isoleucine-zipper trimer

Subunit oligomerization in many proteins is mediated by short coiled-coil motifs. These motifs share a characteristic seven-amino-acid repeat containing hydrophobic residues at the first (a) and fourth (d) positions. Despite this common pattern, different sequences form two-, three- and four-stranded helical ropes. We have investigated the basis for oligomer choice by characterizing variants of the GCN4 leucine-zipper dimerization domain that adopt trimeric or tetrameric structures in response to mutations at the a and d positions. We now report the high-resolution X-ray crystal structure of an isoleucine-containing mutant that folds into a parallel three-stranded, alpha-helical coiled coil. In contrast to the dimer and tetramer structures, the interior packing of the trimer can accommodate beta-branched residues in the most preferred rotamer at both hydrophobic positions. Compatibility of the shape of the core amino acids with the distinct packing spaces in the two-, three- and four-stranded conformations appears to determine the oligomerization state of the GCN4 leucine-zipper variants.

A switch between two-, three-, and four-stranded coiled coils in GCN4 leucine zipper mutants

Coiled-coil sequences in proteins consist of heptad repeats containing two characteristic hydrophobic positions. The role of these buried hydrophobic residues in determining the structures of coiled coils was investigated by studying mutants of the GCN4 leucine zipper. When sets of buried residues were altered, two-, three-, and four-helix structures were formed. The x-ray crystal structure of the tetramer revealed a parallel, four-stranded coiled coil. In the tetramer conformation, the local packing geometry of the two hydrophobic positions in the heptad repeat is reversed relative to that in the dimer. These studies demonstrate that conserved, buried residues in the GCN4 leucine zipper direct dimer formation. In contrast to proposals that the pattern of hydrophobic and polar amino acids in a protein sequence is sufficient to determine three-dimensional structure, the shapes of buried side chains in coiled coils are essential determinants of the global fold.

Extraction of protein-bound ATP and ADP from human platelets in plasma

Actin is the major ATP and ADP binding protein in platelets, 0.9-1.3 nmol/10(8) cells, 50-70% in the unpolymerized state. The goal of these experiments was to develop a method for extracting all protein-bound ATP and ADP from undisturbed platelets in plasma. Extraction of actin-bound ADP is routine while extraction of actin-bound ATP from platelets in buffer has been unsuccessful. Prior to extraction the platelets were exposed to 14-C adenine, to label the metabolic and actin pools of ATP and ADP. The specific activity was determined from the actin-bound ADP in the 43% ethanol precipitate. Sequential ethanol and perchlorate extractions of platelet rich plasma, and the derived supernatants and precipitates were performed. ATP concentrations were determined with the luciferase assay, and radioactive nucleotides separated by TLC. A total of 1.18 nmol/10(8) cells of protein-bound ATP and ADP was recovered, 52% ATP (0.61 nmol). The recovery of protein-bound ADP was increased from 0.3 to 0.57 nmol/10(8) cells. This approach for the first time successfully recovered protein bound ATP and ADP from platelets in a concentration expected for actin.