Derivation and testing of pair potentials for protein folding. When is the quasichemical approximation correct?

Many existing derivations of knowledge-based statistical pair potentials invoke the quasichemical approximation to estimate the expected side-chain contact frequency if there were no amino acid pair-specific interactions. At first glance, the quasichemical approximation that treats the residues in a protein as being disconnected and expresses the side-chain contact probability as being proportional to the product of the mole fractions of the pair of residues would appear to be rather severe. To investigate the validity of this approximation, we introduce two new reference states in which no specific pair interactions between amino acids are allowed, but in which the connectivity of the protein chain is retained. The first estimates the expected number of side-chain contracts by treating the protein as a Gaussian random coil polymer. The second, more realistic reference state includes the effects of chain connectivity, secondary structure, and chain compactness by estimating the expected side-chain contrast probability by placing the sequence of interest in each member of a library of structures of comparable compactness to the native conformation. The side-chain contact maps are not allowed to readjust to the sequence of interest, i.e., the side chains cannot repack. This situation would hold rigorously if all amino acids were the same size. Both reference states effectively permit the factorization of the side-chain contact probability into sequence-dependent and structure-dependent terms. Then, because the sequence distribution of amino acids in proteins is random, the quasichemical approximation to each of these reference states is shown to be excellent. Thus, the range of validity of the quasichemical approximation is determined by the magnitude of the side-chain repacking term, which is, at present, unknown. Finally, the performance of these two sets of pair interaction potentials as well as side-chain contact fraction-based interaction scales is assessed by inverse folding tests both without and with allowing for gaps.     

Structure characterization of the central repetitive domain of high molecular weight gluten proteins. I. Model studies using cyclic and linear peptides.

The high molecular weight (HMW) proteins from wheat contain a repetitive domain that forms 60-80% of their sequence. The consensus peptides PGQGQQ and GYYPTSPQQ form more than 90% of the domain; both are predicted to adopt beta-turn structure. This paper describes the structural characterization of these consensus peptides and forms the basis for the structural characterization of the repetitive HMW domain, described in the companion paper. The cyclic peptides cyclo-[PGQGQQPGQGQQ] (peptide 1), cyclo-[GYYPTSPQQGA] (peptide 2), and cyclo-[PGQGQQGYYPTSPQQ] (peptide 3) were prepared using a novel synthesis route. In addition, the linear peptides (PGQGQQ)n (n = 1, 3, 5) were prepared. CD, FTIR, and NMR data demonstrated a type II beta-turn structure at QPGQ in the cyclic peptide 1 that was also observed in the linear peptides 9PGQGQQ)n. A type I beta-turn was observed at YPTS and SPQQ in peptides 2 and 3, with additional beta-turns of either type I or II at GAGY (peptide 2) and QQGY (peptide 3). The proline in YPTS showed considerable cis/trans isomerization, with up to 50% of the population in the cis-conformation; the other prolines were more than 90% in the trans conformation. The conversion from trans to cis destroys the type I beta-turn at YPTS, but leads to an increase in turn character at SPQQ and GAGY (peptide 2) or QQGY (peptide 3).     

Crystal structure of the hydroxylase component of methane monooxygenase from Methylosinus trichosporium OB3b.

Methane monooxygenase (MMO), found in aerobic methanotrophic bacteria, catalyzes the O2-dependent conversion of methane to methanol. The soluble form of the enzyme (sMMO) consists of three components: a reductase, a regulatory "B" component (MMOB), and a hydroxylase component (MMOH), which contains a hydroxo-bridged dinuclear iron cluster. Two genera of methanotrophs, termed Type X and Type II, which differ markedly in cellular and metabolic characteristics, are known to produce the sMMO. The structure of MMOH from the Type X methanotroph Methylococcus capsulatus Bath (MMO Bath) has been reported recently. Two different structures were found for the essential diiron cluster, depending upon the temperature at which the diffraction data were collected. In order to extend the structural studies to the Type II methanotrophs and to determine whether one of the two known MMOH structures is generally applicable to the MMOH family, we have determined the crystal structure of the MMOH from Type II Methylosinus trichosporium OB3b (MMO OB3b) in two crystal forms to 2.0 A resolution, respectively, both determined at 18 degrees C. The crystal forms differ in that MMOB was present during crystallization of the second form. Both crystal forms, however, yielded very similar results for the structure of the MMOH. Most of the major structural features of the MMOH Bath were also maintained with high fidelity. The two irons of the active site cluster of MMOH OB3b are bridged by two OH (or one OH and one H2O), as well as both carboxylate oxygens of Glu alpha 144. This bis-mu-hydroxo-bridged "diamond core" structure, with a short Fe-Fe distance of 2.99 A, is unique for the resting state of proteins containing analogous diiron clusters, and is very similar to the structure reported for the cluster from flash frozen (-160 degrees C) crystals of MMOH Bath, suggesting a common active site structure for the soluble MMOHs. The high-resolution structure of MMOH OB3b indicates 26 consecutive amino acid sequence differences in the beta chain when compared to the previously reported sequence inferred from the cloned gene. Fifteen additional sequence differences distributed randomly over the three chains were also observed, including D alpha 209E, a ligand of one of the irons.     

Crystallization and preliminary X-ray diffraction studies of aSFP, a bovine seminal plasma protein with a single CUB domain architecture.

Bovine acidic seminal fluid protein (aSFP) is a 1.29 kDa polypeptide of the spermadhesin family built by a single CUB domain architecture. The CUB domain is an extracellular module present in 16 functionally diverse proteins. To determine the three-dimensional structure of aSFP, the protein was crystallized at 21 degrees C by vapor diffusion in hanging drops, using ammonium sulfate, pH 4.7, and polyethyleneglycol 4,000 as precipitants, containing 10% dioxane to avoid the formation of clustered crystals. Elongated prismatic crystals with maximal size of 0.6 x 0.3 x 0.2 mm3 diffract to beyond 1.9 A resolution and belong to space group P2(1)2(1)2(1), with cell parameters a = 52.4 A, b = 41.5 A, c = 48.2 A. There is one aSFP molecule per asymmetric unit, which corresponds to a crystal volume per unit molecular mass of 2.04 A3/Da, and analytical ultracentrifugation analysis show that aSFP is a monomeric protein.     

A flavodoxin that is required for enzyme activation: the structure of oxidized flavodoxin from Escherichia coli at 1.8 A resolution.

In Escherichia coli, flavodoxin is the physiological electron donor for the reductive activation of the enzymes pyruvate formate-lyase, anaerobic ribonucleotide reductase, and B12-dependent methionine synthase. As a basis for studies of the interactions of flavodoxin with methionine synthase, crystal structures of orthorhombic and trigonal forms of oxidized recombinant flavodoxin from E. coli have been determined. The orthorhombic form (space group P2(1)2(1)2(1), a = 126.4, b = 41.10, c = 69.15 A, with two molecules per asymmetric unit) was solved initially by molecular replacement at a resolution of 3.0 A, using coordinates from the structure of the flavodoxin from Synechococcus PCC 7942 (Anacystis nidulans). Data extending to 1.8-A resolution were collected at 140 K and the structure was refined to an Rwork of 0.196 and an Rfree of 0.250 for reflections with I > 0. The final model contains 3,224 non-hydrogen atoms per asymmetric unit, including 62 flavin mononucleotide (FMN) atoms, 354 water molecules, four calcium ions, four sodium ions, two chloride ions, and two Bis-Tris buffer molecules. The structure of the protein in the trigonal form (space group P312, a = 78.83, c = 52.07 A) was solved by molecular replacement using the coordinates from the orthorhombic structure, and was refined with all data from 10.0 to 2.6 A (R = 0.191; Rfree = 0.249). The sequence Tyr 58-Tyr 59, in a bend near the FMN, has so far been found only in the flavodoxins from E. coli and Haemophilus influenzae, and may be important in interactions of flavodoxin with its partners in activation reactions. The tyrosine residues in this bend are influenced by intermolecular contacts and adopt different orientations in the two crystal forms. Structural comparisons with flavodoxins from Synechococcus PCC 7942 and Anaebaena PCC 7120 suggest other residues that may also be critical for recognition by methionine synthase.     

Structure characterization of the central repetitive domain of high molecular weight gluten proteins. II. Characterization in solution and in the dry state.

The structure of the central repetitive domain of high molecular weight HMW) wheat gluten proteins was characterized in solution and in the dry state using HMW proteins Bx6 and Bx7 and a subcloned, bacterially expressed part of the repetitive domain of HMW Dx5. Model studies of the HMW consensus peptides PGQGQQ and GYYPTSPQQ formed the basis for the data analysis (van Dijk AA et al., 1997, Protein Sci 6:637-648). In solution, the repetitive domain contained a continuous nonoverlapping series of both type I and type II II beta-turns at positions predicted from the model studies; type II beta-turns occurred at QPGQ and QQGY sequences and type I beta-turns at YPTS and SPQQ. The subcloned part of the HMW Dx5 repetitive domain sometimes migrated as two bands on SDS-PAGE; we present evidence that this may be caused by a single amino acid insertion that disturbs the regular structure of beta-turns. The type I beta-turns are lost when the protein is dried on a solid surface, probably by conversion to type II beta-turns. The homogeneous type II beta-turn distribution is compatible with the formation of a beta-spiral structure, which provides the protein with elastic properties. The beta-turns and thus the beta-spiral are stabilized by hydrogen bonds within and between turns. Reformation of this hydrogen bonding network after, e.g., mechanical disruption may be important for the elastic properties of gluten proteins.     

Crystallization and preliminary X-ray crystallographic analysis of squalene-hopene cyclase from Alicyclobacillus acidocaldarius.

The membrane-associated protein squalene-hopene cyclase from Alicyclobacillus acidocaldarius was overexposed in Escherichia coli and purified by ion exchange and gel permeation chromatography. Crystals of three interrelated forms were grown by vapor diffusion under identical conditions. The crystals diffract to about 2.3 A resolution, but they are unstable in the X-ray beam. An interpretable heavy-atom derivative was obtained.     

Characterization and crystallization of human uroporphyrinogen decarboxylase.

The cytosolic enzyme uroporphyrinogen decarboxylase (URO-D) catalyzes the fifth step in the heme biosynthetic pathway, converting uroporphyrinogen to coproporphyrinogen by decarboxylating the four acetate side chains of the substrate. Recombinant human URO-D has been expressed in Escherichia coli with a histidine tag and has been purified to homogeneity. Purified protein was determined to be a monodisperse dimer by dynamic light scattering. Equilibrium sedimentation analysis confirmed that the protein is dimeric, with a dissociation constant of 0.1 microM. URO-D containing an amino-terminal histidine tag was crystallized in space group P3(1)21 or its enantiomer P3(2)21 with unit cell dimensions a = b = 103.6 A, c = 75.2 A. There is one molecule in the asymmetric unit, consistent with generation of the dimer by the twofold axis of this crystallographic operator. Native data have been collected to 3.0 a resolution.     

Crystallization of the first three domains of the human insulin-like growth factor-1 receptor.

The insulin-like growth factor-1 receptor (IGF-1R) is a tyrosine kinase receptor of central importance in cell proliferation. A fragment (residues 1-462) comprising the L1-cysteine rich-L2 domains of the human IGF-1R ectodomain has been overexpressed in glycosylation-deficient Lec8 cells and has been affinity-purified via a c-myc tag followed by gel filtration. The fragment was recognized by two anti-IGF-1R monoclonal antibodies, 24-31 and 24-60, but showed no detectable binding of IGF-1 or IGF-2. Isocratic elution of IGF-1R/462 on anion-exchange chromatography reduced sample heterogeneity, permitting the production of crystals that diffracted to 2.6 A resolution with cell dimensions a = 77.0 A, b = 99.5 A, c = 120.1 A, and space group P2(1)2(1)2(1).     

Characterization of the allosteric binding pocket of human liver fructose-1,6-bisphosphatase by protein crystallography and inhibitor activity studies.

The structures of three complexes of human fructose-1,6-bisphosphatase (FB) with the allosteric inhibitor AMP and two AMP analogues have been determined and all fully refined. The data used for structure determination were collected at cryogenic temperature (110 K), and with the use of synchrotron radiation. The structures reveal a common mode of binding for AMP and formycine monophosphate (FMP). 5-Amino-4-carboxamido-1 beta-D-5-phosphate-ribofuranosyl-1H-imidazole (AICAR-P) shows an unexpected mode of binding to FB, different from that of the other two ligands. The imidazole ring of AICAR-P is rotated 180 degrees compared to the AMP and FMP bases. This rotation results in a slightly different hydrogen bonding pattern and minor changes in the water structure in the binding pocket. Common features of binding are seen for the ribose and phosphate moieties of all three compounds. Although binding in a different mode, AICAR-P is still capable of making all the important interactions with the residues building the allosteric binding pocket. The IC50 values of AMP, FMP, and AICAR-P were determined to be 1.7, 1.4, and 20.9 microM, respectively. Thus, the approximately 10 times lower potency of AICAR-P is difficult to explain solely from the variations observed in the binding pocket. Only one water molecule in the allosteric binding pocket was found to be conserved in all four subunits in all three structures. This water molecule coordinates to a phosphate oxygen atom and the N7 atom of the AMP molecule, and to similarly situated atoms in the FMP and AICAR-P complexes. This implies an important role of the conserved water molecule in binding of the ligand.