Cross-validated maximum likelihood enhances crystallographic simulated annealing refinement

Recently, the target function for crystallographic refinement has been improved through a maximum likelihood analysis, which makes proper allowance for the effects of data quality, model errors, and incompleteness. The maximum likelihood target reduces the significance of false local minima during the refinement process, but it does not completely eliminate them, necessitating the use of stochastic optimization methods such as simulated annealing for poor initial models. It is shown that the combination of maximum likelihood with cross-validation, which reduces overfitting, and simulated annealing by torsion angle molecular dynamics, which simplifies the conformational search problem, results in a major improvement of the radius of convergence of refinement and the accuracy of the refined structure. Torsion angle molecular dynamics and the maximum likelihood target function interact synergistically, the combination of both methods being significantly more powerful than each method individually. This is demonstrated in realistic test cases at two typical minimum Bragg spacings (dmin = 2.0 and 2.8 Å, respectively), illustrating the broad applicability of the combined method. In an application to the refinement of a new crystal structure, the combined method automatically corrected a mistraced loop in a poor initial model, moving the backbone by 4 Å.

Crystallographic comparison of the estrogen and progesterone receptor’s ligand binding domains

The 2.8-Å crystal structure of the complex formed by estradiol and the human estrogen receptor-α ligand binding domain (hERαLBD) is described and compared with the recently reported structure of the progesterone complex of the human progesterone receptor ligand binding domain, as well as with similar structures of steroid/nuclear receptor LBDs solved elsewhere. The hormone-bound hERαLBD forms a distinctly different and probably more physiologically important dimer interface than its progesterone counterpart. A comparison of the specificity determinants of hormone binding reveals a common structural theme of mutually supported van der Waals and hydrogen-bonded interactions involving highly conserved residues. The previously suggested mechanism by which the estrogen receptor distinguishes estradiol’s unique 3-hydroxy group from the 3-keto function of most other steroids is now described in atomic detail. Mapping of mutagenesis results points to a coactivator-binding surface that includes the region around the “signature sequence” as well as helix 12, where the ligand-dependent conformation of the activation function 2 core is similar in all previously solved steroid/nuclear receptor LBDs. A peculiar crystal packing event displaces helix 12 in the hERαLBD reported here, suggesting a higher degree of dynamic variability than expected for this critical substructure.

Peptide hybrids containing α- and β-amino acids: Structure of a decapeptide β-hairpin with two facing β-phenylalanine residues

A β-hairpin conformation has been characterized in crystals of the decapeptide t-butoxycarbonyl-Leu-Val-βPhe-Val-DPro-Gly-Leu-βPhe-Val-Val-methyl ester [βPhe; (S)-β3 homophenylalanine] by x-ray diffraction. The polypeptide chain reversal is nucleated by the centrally positioned DPro-Gly segment, which adopts a type-I′ β-turn conformation. Four intramolecular cross-strand hydrogen bonds stabilize the peptide fold. The βPhe(3) and βPhe(8) residues occupy facing positions on the hairpin, with the side chains projecting on opposite faces of the β-sheet. At the site of insertion of β-residues, the polarity of the peptide units along each strand reverses, as compared with the α-peptide segments. In this analog, a small segment of a polar sheet is observed, where adjacent CO and NH groups line up in opposite directions in each strand. In the crystal, an extended β-sheet is formed by hydrogen bonding between strands of antiparallel pairs of β-hairpins. The crystallographic parameters for C65H102N10O13⋅ 3H2O are: space group P212121; a = 19.059(8) Å, b = 19.470(2) Å, c = 21.077(2) Å; Z = 4; agreement factor R1 = 9.12% for 3,984 data observed >4σ(F) and a resolution of 0.90 Å.

The binding conformation of Taxol in β-tubulin: A model based on electron crystallographic density

The chemotherapeutic drug Taxol is known to interact within a specific site on β-tubulin. Although the general location of the site has been defined by photoaffinity labeling and electron crystallography, the original data were insufficient to make an absolute determination of the bound conformation. We have now correlated the crystallographic density with analysis of Taxol conformations and have found the unique solution to be a T-shaped Taxol structure. This T-shaped or butterfly structure is optimized within the β-tubulin site and exhibits functional similarity to a portion of the B9-B10 loop in the α-tubulin subunit. The model provides structural rationalization for a sizeable body of Taxol structure–activity relationship data, including binding affinity, photoaffinity labeling, and acquired mutation in human cancer cells.

A crystallographic map of the transition from B-DNA to A-DNA

The transition between B- and A-DNA was first observed nearly 50 years ago. We have now mapped this transformation through a set of single-crystal structures of the sequence d(GGCGCC)2, with various intermediates being trapped by methylating or brominating the cytosine bases. The resulting pathway progresses through 13 conformational steps, with a composite structure that pairs A-nucleotides with complementary B-nucleotides serving as a distinct transition intermediate. The details of each step in the conversion of B- to A-DNA are thus revealed at the atomic level, placing intermediates for this and other sequences in the context of a common pathway.

The Three-Dimensional Structure of Aspergillus niger Pectin Lyase B at 1.7-Å Resolution1

The three-dimensional structure of Aspergillus niger pectin lyase B (PLB) has been determined by crystallographic techniques at a resolution of 1.7 Å. The model, with all 359 amino acids and 339 water molecules, refines to a final crystallographic R factor of 16.5%. The polypeptide backbone folds into a large right-handed cylinder, termed a parallel β helix. Loops of various sizes and conformations protrude from the central helix and probably confer function. The largest loop of 53 residues folds into a small domain consisting of three antiparallel β strands, one turn of an α helix, and one turn of a 310 helix. By comparison with the structure of Erwinia chrysanthemi pectate lyase C (PelC), the primary sequence alignment between the pectate and pectin lyase subfamilies has been corrected and the active site region for the pectin lyases deduced. The substrate-binding site in PLB is considerably less hydrophilic than the comparable PelC region and consists of an extensive network of highly conserved Trp and His residues. The PLB structure provides an atomic explanation for the lack of a catalytic requirement for Ca2+ in the pectin lyase family, in contrast to that found in the pectate lyase enzymes. Surprisingly, however, the PLB site analogous to the Ca2+ site in PelC is filled with a positive charge provided by a conserved Arg in the pectin lyases. The significance of the finding with regard to the enzymatic mechanism is discussed.

Crystallographic analysis of endogenous peptides associated with HLA-DR1 suggests a common, polyproline II-like conformation for bound peptides.

The structure of the human major histocompatibility complex (MHC) class II molecule HLA-DR1 derived from the human lymphoblastoid cell line LG-2 has been determined in a complex with the Staphylococcus aureus enterotoxin B superantigen. The HLA-DR1 molecule contains a mixture of endogenous peptides derived from cellular or serum proteins bound in the antigen-binding site, which copurify with the class II molecule. Continuous electron density for 13 amino acid residues is observed in the MHC peptide-binding site, suggesting that this is the core length of peptide that forms common interactions with the MHC molecule. Electron density is also observed for side chains of the endogenous peptides. The electron density corresponding to peptide side chains that interact with the DR1-binding site is more clearly defined than the electron density that extends out of the binding site. The regions of the endogenous peptides that interact with DRI are therefore either more restricted in conformation or sequence than the peptide side chains or amino acids that project out of the peptide-binding site. The hydrogen-bond interactions and conformation of a peptide model built into the electron density are similar to other HLA-DR-peptide structures. The bound peptides assume a regular conformation that is similar to a polyproline type II helix. The side-chain pockets and conserved asparagine residues of the DR1 molecule are well-positioned to interact with peptides in the polyproline type II conformation and may restrict the range of acceptable peptide conformations.

Structure determination of murine mitochondrial carbonic anhydrase V at 2.45-A resolution: implications for catalytic proton transfer and inhibitor design.

The three-dimensional structure of murine mitochondrial carbonic anhydrase V has been determined and refined at 2.45-A resolution (crystallographic R factor = 0.187). Significant structural differences unique to the active site of carbonic anhydrase V are responsible for differences in the mechanism of catalytic proton transfer as compared with other carbonic anhydrase isozymes. In the prototypical isozyme, carbonic anhydrase II, catalytic proton transfer occurs via the shuttle group His-64; carbonic anhydrase V has Tyr-64, which is not an efficient proton shuttle due in part to the bulky adjacent side chain of Phe-65. Based on analysis of the structure of carbonic anhydrase V, we speculate that Tyr-131 may participate in proton transfer due to its proximity to zinc-bound solvent, its solvent accessibility, and its electrostatic environment in the protein structure. Finally, the design of isozyme-specific inhibitors is discussed in view of the complex between carbonic anhydrase V and acetazolamide, a transition-state analogue. Such inhibitors may be physiologically important in the regulation of blood glucose levels.

Crystal structure of Pseudomonas aeruginosa catabolic ornithine transcarbamoylase at 3.0-A resolution: a different oligomeric organization in the transcarbamoylase family.

The crystal structure of the Glu-105-->Gly mutant of catabolic ornithine transcarbamoylase (OTCase; carbamoyl phosphate + L-ornithine = orthophosphate + L-citrulline, EC 2.1.3.3) from Pseudomonas aeruginosa has been determined at 3.0-A resolution. This mutant is blocked in the active R (relaxed) state. The structure was solved by the molecular replacement method, starting from a crude molecular model built from a trimer of the catalytic subunit of another transcarbamoylase, the extensively studied aspartate transcarbamoylase (ATCase) from Escherichia coli. This model was used to generate initial low-resolution phases at 8-A resolution, which were extended to 3-A by noncrystallographic symmetry averaging. Four phase extensions were required to obtain an electron density map of very high quality from which the final model was built. The structure, including 4020 residues, has been refined to 3-A, and the current crystallographic R value is 0.216. No solvent molecules have been added to the model. The catabolic OTCase is a dodecamer composed of four trimers organized in a tetrahedral manner. Each monomer is composed of two domains. The carbamoyl phosphate binding domain shows a strong structural homology with the equivalent ATCase part. In contrast, the other domain, mainly implicated in the binding of the second substrate (ornithine for OTCase and aspartate for ATCase) is poorly conserved. The quaternary structures of these two allosteric transcarbamoylases are quite divergent: the E. coli ATCase has pseudo-32 point-group symmetry, with six catalytic and six regulatory chains; the catabolic OTCase has 23 point-group symmetry and only catalytic chains. However, both enzymes display homotropic and heterotropic cooperativity.

A cooperative model for receptor recognition and cell adhesion: evidence from the molecular packing in the 1.6-A crystal structure of the pheromone Er-1 from the ciliated protozoan Euplotes raikovi.

The crystal structure of the pheromone Er-1 from the unicellular eukaryotic organism Euplotes raikovi was determined at 1.6 A resolution and refined to a crystallographic R factor of 19.9%. In the tightly packed crystal, two extensive intermolecular helix-helix interactions arrange the Er-1 molecules into layers. Since the putative receptor of the pheromone is a membrane-bound protein, whose extracellular C-terminal domain is identical in amino acid sequence to the soluble pheromone, the interactions found in the crystal may mimic the pheromone-receptor interactions as they occur on a cell surface. Based on this, we propose a model for the interaction between soluble pheromone molecules and their receptors. In this model, strong pheromone-receptor binding emerges as a consequence of the cooperative utilization of several weak interactions. The model offers an explanation for the results of binding studies and may also explain the adhesion between cells that occurs during mating.

The three-dimensional structure of NAD(P)H:quinone reductase, a flavoprotein involved in cancer chemoprotection and chemotherapy: mechanism of the two-electron reduction.

Quinone reductase [NAD(P)H:(quinone acceptor) oxidoreductase, EC 1.6.99.2], also called DT diaphorase, is a homodimeric FAD-containing enzyme that catalyzes obligatory NAD(P)H-dependent two-electron reductions of quinones and protects cells against the toxic and neoplastic effects of free radicals and reactive oxygen species arising from one-electron reductions. These two-electron reductions participate in the reductive bioactivation of cancer chemotherapeutic agents such as mitomycin C in tumor cells. Thus, surprisingly, the same enzymatic reaction that protects normal cells activates cytotoxic drugs used in cancer chemotherapy. The 2.1-A crystal structure of rat liver quinone reductase reveals that the folding of a portion of each monomer is similar to that of flavodoxin, a bacterial FMN-containing protein. Two additional portions of the polypeptide chains are involved in dimerization and in formation of the two identical catalytic sites to which both monomers contribute. The crystallographic structures of two FAD-containing enzyme complexes (one containing NADP+, the other containing duroquinone) suggest that direct hydride transfers from NAD(P)H to FAD and from FADH2 to the quinone [which occupies the site vacated by NAD(P)H] provide a simple rationale for the obligatory two-electron reductions involving a ping-pong mechanism.

Periodic variation in side-chain polarities of T-cell antigenic peptides correlates with their structure and activity.

We present an analysis that synthesizes information on the sequence, structure, and motifs of antigenic peptides, which previously appeared to be in conflict. Fourier analysis of T-cell antigenic peptides indicates a periodic variation in amino acid polarities of 3-3.6 residues per period, suggesting an amphipathic alpha-helical structure. However, the diffraction patterns of major histocompatibility complex (MHC) molecules indicate that their ligands are in an extended non-alpha-helical conformation. We present two mutually consistent structural explanations for the source of the alpha-helical periodicity, based on an observation that the side chains of MHC-bound peptides generally partition with hydrophobic (hydrophilic) side chains pointing into (out of) the cleft. First, an analysis of haplotype-dependent peptide motifs indicates that the locations of their defining residues tend to force a period 3-4 variation in hydrophobicity along the peptide sequence, in a manner consistent with the spacing of pockets in the MHC. Second, recent crystallographic determination of the structure of a peptide bound to a class II MHC molecule reveals an extended but regularly twisted peptide with a rotation angle of about 130 degrees. We show that similar structures with rotation angles of 100-130 degrees are energetically acceptable and also span the length of the MHC cleft. These results provide a sound physical chemical and structural basis for the existence of a haplotype-independent antigenic motif which can be particularly important in limiting the search time for antigenic peptides.

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.

The crystal structure of an N-terminal two-domain fragment of vascular cell adhesion molecule 1 (VCAM-1): a cyclic peptide based on the domain 1 C-D loop can inhibit VCAM-1-alpha 4 integrin interaction.

Vascular cell adhesion molecule 1 (VCAM-1) represents a structurally and functionally distinct class of immunoglobulin superfamily molecules that bind leukocyte integrins and are involved in inflammatory and immune functions. X-ray crystallography defines the three-dimensional structure of the N-terminal two-domain fragment that participates in ligand binding. Residues in domain 1 important for ligand binding reside in the C-D loop, which projects markedly from one face of the molecule near the contact between domains 1 and 2. A cyclic peptide that mimics this loop inhibits binding of alpha 4 beta 1 integrin-bearing cells to VCAM-1. These data demonstrate how crystallographic structural information can be used to design a small molecule inhibitor of biological function.

Structure-assisted redesign of a protein-zinc-binding site with femtomolar affinity.

We have inserted a fourth protein ligand into the zinc coordination polyhedron of carbonic anhydrase II (CAII) that increases metal affinity 200-fold (Kd = 20 fM). The three-dimensional structures of threonine-199-->aspartate (T199D) and threonine-199-->glutamate (T199E) CAIIs, determined by x-ray crystallographic methods to resolutions of 2.35 Angstrum and 2.2 Angstrum, respectively, reveal a tetrahedral metal-binding site consisting of H94, H96, H119, and the engineered carboxylate side chain, which displaces zinc-bound hydroxide. Although the stereochemistry of neither engineered carboxylate-zinc interaction is comparable to that found in naturally occurring protein zinc-binding sites, protein-zinc affinity is enhanced in T199E CAII demonstrating that ligand-metal separation is a significant determinant of carboxylate-zinc affinity. In contrast, the three-dimensional structure of threonine-199-->histidine (T199H) CAII, determined to 2.25-Angstrum resolution, indicates that the engineered imidazole side chain rotates away from the metal and does not coordinate to zinc; this results in a weaker zinc-binding site. All three of these substitutions nearly obliterate CO2 hydrase activity, consistent with the role of zinc-bound hydroxide as catalytic nucleophile. The engineering of an additional protein ligand represents a general approach for increasing protein-metal affinity if the side chain can adopt a reasonable conformation and achieve inner-sphere zinc coordination. Moreover, this structure-assisted design approach may be effective in the development of high-sensitivity metal ion biosensors.