Henry van den Bedem

Accessing protein conformational ensembles using room-temperature X-ray crystallography

Modern protein crystal structures are based nearly exclusively on X-ray data collected at cryogenic temperatures (generally 100 K). The cooling process is thought to introduce little bias in the functional interpretation of structural results, because cryogenic temperatures minimally perturb the overall protein backbone fold. In contrast, here we show that flash cooling biases previously hidden structural ensembles in protein crystals. By analyzing available data for 30 different proteins using new computational tools for electron-density sampling, model refinement, and molecular packing analysis, we found that crystal cryocooling remodels the conformational distributions of more than 35% of side chains and eliminates packing defects necessary for functional motions. In the signaling switch protein, H-Ras, an allosteric network consistent with fluctuations detected in solution by NMR was uncovered in the room-temperature, but not the cryogenic, electron-density maps. These results expose a bias in structural databases toward smaller, overpacked, and unrealistically unique models. Monitoring room-temperature conformational ensembles by X-ray crystallography can reveal motions crucial for catalysis, ligand binding, and allosteric regulation.

Integrative, dynamic structural biology at atomic resolution[mdash]it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing—in atomic detail—how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

Automated diffraction image analysis and spot searching for high-throughput crystal screening

A new software package, DISTL (Diffraction Image Screening Tool and Library), for the rapid analysis of X-ray diffraction patterns collected from macromolecular crystals is presented. Within seconds, the program characterizes the strength and quality of the Bragg spots, determines the limiting resolution of the image, and identifies deleterious features such as ice-rings and intense salt reflections. The procedure also generates a reliable set of intense peaks for auto-indexing. The ability to classify a large number of crystals quickly will be especially useful at synchrotron and home-laboratory X-ray sources where automated crystal screening and data collection systems have been implemented.

Automated identification of functional dynamic contact networks from X-ray crystallography.

Protein function often depends on the exchange between conformational substates. Allosteric ligand binding or distal mutations can stabilize specific active-site conformations and consequently alter protein function. Observing alternative conformations at low levels of electron density, in addition to comparison of independently determined X-ray crystal structures, can provide mechanistic insights into conformational dynamics. Here we report a new algorithm, CONTACT, that identifies contact networks of conformationally heterogeneous residues directly from high-resolution X-ray crystallography data. Contact networks determined for Escherichia coli dihydrofolate reductase (ecDHFR) predict the observed long-range pattern of NMR chemical shift perturbations of an allosteric mutation. A comparison of contact networks in wild-type and mutant ecDHFR suggests that mutations that alter optimized contact networks of coordinated motions can impair catalytic function. CONTACT-guided mutagenesis can exploit the structure-dynamics-function relationship in protein engineering and design.

Integrated description of protein dynamics from room-temperature X-ray crystallography and NMR.

Significance Most proteins are inherently flexible and their dynamics play a central role in their biological functions. A molecular level understanding of protein function and mechanism requires an accurate description of the atomic coordinates in both time and space. Here we show, through studies of the enzyme dihydrofolate reductase, that multiconformer models derived from room-temperature X-ray crystallographic data can be used synergistically with nuclear magnetic resonance relaxation measurements to provide a detailed description of both the amplitude and timescale of fluctuations in atomic coordinates. This hybrid approach provides a more complete description of protein dynamics than can be obtained from either method alone. The room-temperature crystallographic ensemble accurately reflects the picosecond–nanosecond motions of the protein backbone and side chains. Detailed descriptions of atomic coordinates and motions are required for an understanding of protein dynamics and their relation to molecular recognition, catalytic function, and allostery. Historically, NMR relaxation measurements have played a dominant role in the determination of the amplitudes and timescales (picosecond–nanosecond) of bond vector fluctuations, whereas high-resolution X-ray diffraction experiments can reveal the presence of and provide atomic coordinates for multiple, weakly populated substates in the protein conformational ensemble. Here we report a hybrid NMR and X-ray crystallography analysis that provides a more complete dynamic picture and a more quantitative description of the timescale and amplitude of fluctuations in atomic coordinates than is obtainable from the individual methods alone. Order parameters (S2) were calculated from single-conformer and multiconformer models fitted to room temperature and cryogenic X-ray diffraction data for dihydrofolate reductase. Backbone and side-chain order parameters derived from NMR relaxation experiments are in excellent agreement with those calculated from the room-temperature single-conformer and multiconformer models, showing that the picosecond timescale motions observed in solution occur also in the crystalline state. These motions are quenched in the crystal at cryogenic temperatures. The combination of NMR and X-ray crystallography in iterative refinement promises to provide an atomic resolution description of the alternate conformational substates that are sampled through picosecond to nanosecond timescale fluctuations of the protein structure. The method also provides insights into the structural heterogeneity of nonmethyl side chains, aromatic residues, and ligands, which are less commonly analyzed by NMR relaxation measurements.

Crystal structures of two novel dye-decolorizing peroxidases reveal a beta-barrel fold with a conserved heme-binding motif.

BtDyP from Bacteroides thetaiotaomicron (strain VPI‐5482) and TyrA from Shewanella oneidensis are dye‐decolorizing peroxidases (DyPs), members of a new family of heme‐dependent peroxidases recently identified in fungi and bacteria. Here, we report the crystal structures of BtDyP and TyrA at 1.6 and 2.7 Å, respectively. BtDyP assembles into a hexamer, while TyrA assembles into a dimer; the dimerization interface is conserved between the two proteins. Each monomer exhibits a two‐domain, α+β ferredoxin‐like fold. A site for heme binding was identified computationally, and modeling of a heme into the proposed active site allowed for identification of residues likely to be functionally important. Structural and sequence comparisons with other DyPs demonstrate a conservation of putative heme‐binding residues, including an absolutely conserved histidine. Isothermal titration calorimetry experiments confirm heme binding, but with a stoichiometry of 0.3:1 (heme:protein). Proteins 2007.

Crystal structure of thy1, a thymidylate synthase complementing protein from Thermotoga maritima at 2.25 Å resolution

Peter Kuhn, Scott A. Lesley, Irimpan I. Mathews, Jaume M. Canaves, Linda S. Brinen, Xiaoping Dai, Ashley M. Deacon, Marc A. Elsliger, Said Eshaghi, Ross Floyd, Adam Godzik, Carina Grittini, Slawomir K. Grzechnik, Chittibabu Guda, Keith O. Hodgson, Lukasz Jaroszewski, Cathy Karlak, Heath E. Klock, Eric Koesema, John M. Kovarik, Andreas T. Kreusch, Daniel McMullan, Timothy M. McPhillips, Mark A. Miller, Mitchell Miller, Andrew Morse, Kin Moy, Jie Ouyang, Alyssa Robb, Kevin Rodrigues, Thomas L. Selby, Glen Spraggon, Raymond C. Stevens, Susan S. Taylor, Henry van den Bedem, Jeff Velasquez, Juli Vincent, Xianhong Wang, Bill West, Guenter Wolf, John Wooley, and Ian A. Wilson* The Joint Center for Structural Genomics Stanford Synchrotron Radiation Laboratory, Stanford University, Menlo Park, California The Genomics Institute of Novartis Foundation, San Diego, California The San Diego Supercomputer Center, La Jolla, California The University of California, San Diego, La Jolla, California The Scripps Research Institute, La Jolla, California

Real-space protein-model completion: an inverse-kinematics approach.

Rapid protein-structure determination relies greatly on software that can automatically build a protein model into an experimental electron-density map. In favorable circumstances, various software systems are capable of building over 90% of the final model. However, completeness falls off rapidly with the resolution of the diffraction data. Manual completion of these partial models is usually feasible, but is time-consuming and prone to subjective interpretation. Except for the N- and C-termini of the chain, the end points of each missing fragment are known from the initial model. Hence, fitting fragments reduces to an inverse-kinematics problem. A method has been developed that combines fast inverse-kinematics algorithms with a real-space torsion-angle refinement procedure in a two-stage approach to fit missing main-chain fragments into the electron density between two anchor points. The first stage samples a large number of closing conformations, guided by the electron density. These candidates are ranked according to density fit. In a subsequent refinement stage, optimization steps are projected onto a carefully chosen subspace of conformation space to preserve rigid geometry and closure. Experimental results show that fitted fragments are in excellent agreement with the final refined structure for lengths of up to 12-15 residues in areas of weak or ambiguous electron density, even at medium to low resolution.

Structural Basis of Murein Peptide Specificity of a γ-D-glutamyl-L-diamino Acid Endopeptidase

The crystal structures of two homologous endopeptidases from cyanobacteria Anabaena variabilis and Nostoc punctiforme were determined at 1.05 and 1.60 A resolution, respectively, and contain a bacterial SH3-like domain (SH3b) and a ubiquitous cell-wall-associated NlpC/P60 (or CHAP) cysteine peptidase domain. The NlpC/P60 domain is a primitive, papain-like peptidase in the CA clan of cysteine peptidases with a Cys126/His176/His188 catalytic triad and a conserved catalytic core. We deduced from structure and sequence analysis, and then experimentally, that these two proteins act as gamma-D-glutamyl-L-diamino acid endopeptidases (EC 3.4.22.-). The active site is located near the interface between the SH3b and NlpC/P60 domains, where the SH3b domain may help define substrate specificity, instead of functioning as a targeting domain, so that only muropeptides with an N-terminal L-alanine can bind to the active site.

Structure of the γ-d-glutamyl-l-diamino acid endopeptidase YkfC from Bacillus cereus in complex with l-Ala-γ-d-Glu: insights into substrate recognition by NlpC/P60 cysteine peptidases

The crystal structure of the highly specific γ-d-glutamyl-l-diamino acid endopeptidase YkfC from Bacillus cereus in complex with l-Ala-γ-d-Glu reveals the structural basis for the substrate specificity of NlpC/P60-family cysteine peptidases.