Binding mode prediction for a flexible ligand in a flexible pocket using multi-conformation simulated annealing pseudo crystallographic refinement

We describe multi-conformation simulated annealing-pseudo-crystallographic refinement (MCSA-PCR), a technique developed for predicting the binding mode of a flexible ligand in a flexible binding pocket. To circumvent the local-minimum problem efficiently, this method performs multiple independent cycles of simulated annealing with explicit solvent, "growing" the ligand in the binding pocket each time. From the ensemble of structures, a pseudo-crystallographic electron density map is calculated, and then conventional crystallographic refinement methods are used to best fit a single, optimal structure into the density map. The advantage of the MCSA-PCR method is that it provides a direct means to evaluate the accuracy and uniqueness of the calculated solution, provides a measure of ligand and protein dynamics from the refined B-factors, and facilitates comparison with X-ray crystallographic data. Here, we show that our MCSA-PCR method succeeds in predicting the correct binding mode of the VSV8 peptide to the major histocompatibility complex (MHC) receptor. Importantly, there is a significant correlation between the experimentally determined crystallographic water molecules and water density observed in the pseudo map by MCSA-PCR. Furthermore, comparison of different approaches for extracting a single, most probable structure from the calculated ensemble reveals the power of the PCR method and provides insights into the nature of the energetic landscape.     

Structure of porcine heart cytoplasmic malate dehydrogenase: combining X-ray diffraction and chemical sequence data in structural studies

The amino acid sequence of cytoplasmic malate dehydrogenase (sMDH) has been determined by a combination of X-ray crystallographic and chemical sequencing methods. The initial molecular model incorporated an "X-ray amino acid sequence" that was derived primarily from an evaluation of a multiple isomorphous replacement phased electron density map calculated at 2.5-A resolution. Following restrained least-squares crystallographic refinement, difference electron density maps were calculated from model phases, and attempts were made to upgrade the X-ray amino acid sequence. The method used to find the positions of peptides in the X-ray structure was similar to those used for studying protein homology and was shown to be successful for large fragments. For sMDH, X-ray methods by themselves were insufficient to derive a complete amino acid sequence, even with partial chemical sequence data. However, for this relatively large molecule at medium resolution, the electron density maps were of considerable help in determining the linear position of peptide fragments. The N-acetylated polypeptide chain of sMDH has 331 amino acids and has been crystallographically refined to an R factor of 19% for 2.5-A resolution diffraction data.     

Crystal structure and electron transfer properties of cytochrome c3

The crystal structure of cytochrome c3 from the sulfate-reducing bacteria Desulfovibrio desulfuricans, Norway strain, has been determined through the fitting of the recently completed primary structure to a 2.5 A resolution electron density map. The phase calculations were based on three mercurial derivatives; anomalous scattering data were used to refine the four heme iron positions. A preliminary refinement of the molecular model has led to a conventional crystallographic R factor of 34%. Cytochrome c3 is folded in two structural domains with one heme in each, the two other heme moieties lying in a large groove dividing the molecule. The core of the protein is the compact four-heme cluster which presents a relatively high degree of solvent exposure. The structural pattern of redox centers suggests that electron transfer might occur through direct contacts between some of the heme groups, via the overlapping system of pi oribitals or via intervening amino acid side chains or both.     

Cryo‐EM map interpretation and protein model‐building using iterative map segmentation

A procedure for building protein chains into maps produced by single‐particle electron cryo‐microscopy (cryo‐EM) is described. The procedure is similar to the way an experienced structural biologist might analyze a map, focusing first on secondary structure elements such as helices and sheets, then varying the contour level to identify connections between these elements. Since the high density in a map typically follows the main‐chain of the protein, the main‐chain connection between secondary structure elements can often be identified as the unbranched path between them with the highest minimum value along the path. This chain‐tracing procedure is then combined with finding side‐chain positions based on the presence of density extending away from the main path of the chain, allowing generation of a C_α model. The C_α model is converted to an all‐atom model and is refined against the map. We show that this procedure is as effective as other existing methods for interpretation of cryo‐EM maps and that it is considerably faster and produces models with fewer chain breaks than our previous methods that were based on approaches developed for crystallographic maps.     

A probabilistic approach to protein backbone tracing in electron density maps

One particularly time-consuming step in protein crystallography is interpreting the electron density map; that is, fitting a complete molecular model of the protein into a 3D image of the protein produced by the crystallographic process. In poor-quality electron density maps, the interpretation may require a significant amount of a crystallographer's time. Our work investigates automating the time-consuming initial backbone trace in poor-quality density maps. We describe ACMI (Automatic Crystallographic Map Interpreter), which uses a probabilistic model known as a Markov field to represent the protein. Residues of the protein are modeled as nodes in a graph, while edges model pairwise structural interactions. Modeling the protein in this manner allows the model to be flexible, considering an almost infinite number of possible conformations, while rejecting any that are physically impossible. Using an efficient algorithm for approximate inference--belief propagation--allows the most probable trace of the protein's backbone through the density map to be determined. We test ACMI on a set of ten protein density maps (at 2.5 to 4.0 A resolution), and compare our results to alternative approaches. At these resolutions, ACMI offers a more accurate backbone trace than current approaches.     

Cluster analysis of time-dependent crystallographic data: Direct identification of time-independent structural intermediates

The initial output of a time-resolved macromolecular crystallography experiment is a time-dependent series of difference electron density maps that displays the time-dependent changes in underlying structure as a reaction progresses. The goal is to interpret such data in terms of a small number of crystallographically refinable, time-independent structures, each associated with a reaction intermediate; to establish the pathways and rate coefficients by which these intermediates interconvert; and thereby to elucidate a chemical kinetic mechanism. One strategy toward achieving this goal is to use cluster analysis, a statistical method that groups objects based on their similarity. If the difference electron density at a particular voxel in the time-dependent difference electron density (TDED) maps is sensitive to the presence of one and only one intermediate, then its temporal evolution will exactly parallel the concentration profile of that intermediate with time. The rationale is therefore to cluster voxels with respect to the shapes of their TDEDs, so that each group or cluster of voxels corresponds to one structural intermediate. Clusters of voxels whose TDEDs reflect the presence of two or more specific intermediates can also be identified. From such groupings one can then infer the number of intermediates, obtain their time-independent difference density characteristics, and refine the structure of each intermediate. We review the principles of cluster analysis and clustering algorithms in a crystallographic context, and describe the application of the method to simulated and experimental time-resolved crystallographic data for the photocycle of photoactive yellow protein.     

Three-dimensional structure of ferricytochrome c' from Rhodospirillum rubrum at 2.8 A resolution.

The structure of ferricytochrome c' extracted from Rhodospirillum rubrum has been determined by the X-ray crystallographic method. Crystals in hexagonal space group P6(1), with unit-cell dimensions a = b = 51.72 A and c = 155.49 A, contain one dimer molecule composed of chemically identical polypeptide chains (monomer I and monomer II) per asymmetric unit. An electron density map has been calculated at a resolution of 2.8 A by the multiple isomorphous replacement method using four-circle diffractometer data from native crystals and two heavy-atom derivatives. The quality of the map was improved by averaging the electron density about the non-crystallographic 2-fold axis relating the two monomers. The initial three-dimensional model of monomer I was built on a computer graphics system and that of monomer II was derived from monomer I using the non-crystallographic symmetry matrices. The dimer structure has been refined using a combination of simulated annealing and conventional restrained least-squares crystallographic refinement. The current model includes 244 amino acid residues (122 x 2) and 2 hemes, with a root-mean-square deviation in bond lengths from ideal values of 0.022 A. The current crystallographic R-factor is 23.3% for 4,481 independent reflections [magnitude of Fo greater than or equal to sigma (F)] between 5.0 and 2.8 A resolution. The monomer molecule is structurally organized as an array of four nearly parallel alpha-helices which construct a left-twisted bundle. One end of the bundle, in which a covalently bound protoheme IX prosthetic group is incorporated, is more divergent than the other.(ABSTRACT TRUNCATED AT 250 WORDS)     

XtalView, protein structure solution and protein graphics, a short history

From a user's point-of-view we are in the Golden Age of protein crystallographic software. In the past few decades, solving protein structures has gone from a task requiring man-months of effort to a process requiring minutes on an ordinary laptop with no human intervention required. The birth of XtalView coincided with the mainstream use of synchrotron radiation, seleno-Met phasing and it continues to be used in this age of robotic crystallization, Fed-Ex data collection and fully automated structure solution "pipelines". This article is a retrospective history of protein crystallographic computing and a discussion of the current state of the art.     

The PDZ2 domain of syntenin at ultra-high resolution: bridging the gap between macromolecular and small molecule crystallography

The crystal structure of the second PDZ domain of the scaffolding protein syntenin was solved using data extending to 0.73 A resolution. The crystallographic model, including the hydrogen atoms and the anisotropic displacement parameters, was refined to a conventional R-factor of 7.5% and Rfree of 8.7%, making it the most precise crystallographic model of a protein molecule to date. The model reveals discrete disorder in several places in the molecule, and significant plasticity of the peptide bond, with some omega angles deviating by nearly 20 degrees from planarity. Most hydrogen atoms are easily identifiable in the electron density and weak hydrogen bonds of the C-H...O type are clearly visible between the beta-strands. The study sets a new standard for high-resolution protein crystallography.     

Using cryoEM Reconstruction and Phase Extension to Determine Crystal Structure of Bacteriophage ϕ6 Major Capsid Protein

Bacteriophage ?6 is a double-stranded RNA virus that has been extensively studied as a model organism. Here we describe structure determination of ?6 major capsid protein P1. The protein crystallized in base centered orthorhombic space group C222₁. Matthews’s coefficient indicated that the crystals contain from four to seven P1 subunits in the crystallographic asymmetric unit. The self-rotation function had shown presence of fivefold axes of non-crystallographic symmetry in the crystals. Thus, electron density map corresponding to a P1 pentamer was excised from a previously determined cryoEM reconstruction of the ?6 procapsid at 7 Å resolution and used as a model for molecular replacement. The phases for reflections at higher than 7 Å resolution were obtained by phase extension employing the fivefold non-crystallographic symmetry present in the crystal. The averaged 3.6 Å-resolution electron density map was of sufficient quality to allow model building.     

Automated electron-density sampling reveals widespread conformational polymorphism in proteins

Although proteins populate large structural ensembles, X-ray diffraction data are traditionally interpreted using a single model. To search for evidence of alternate conformers, we developed a program, Ringer, which systematically samples electron density around the dihedral angles of protein side chains. In a diverse set of 402 structures, Ringer identified weak, nonrandom electron-density features that suggest of the presence of hidden, lowly populated conformations for >18% of uniquely modeled residues. Although these peaks occur at electron-density levels traditionally regarded as noise, statistically significant (P < 10⁻⁵) enrichment of peaks at successive rotameric χ angles validates the assignment of these features as unmodeled conformations. Weak electron density corresponding to alternate rotamers also was detected in an accurate electron density map free of model bias. Ringer analysis of the high-resolution structures of free and peptide-bound calmodulin identified shifts in ensembles and connected the alternate conformations to ligand recognition. These results show that the signal in high-resolution electron density maps extends below the traditional 1 σ cutoff, and crystalline proteins are more polymorphic than current crystallographic models. Ringer provides an objective, systematic method to identify previously undiscovered alternate conformations that can mediate protein folding and function.     

YUP.SCX: coaxing atomic models into medium resolution electron density maps

The structures of large macromolecular complexes in different functional states can be determined by cryo-electron microscopy, which yields electron density maps of low to intermediate resolutions. The maps can be combined with high-resolution atomic structures of components of the complex, to produce a model for the complex that is more accurate than the formal resolution of the map. To this end, methods have been developed to dock atomic models into density maps rigidly or flexibly, and to refine a docked model so as to optimize the fit of the atomic model into the map. We have developed a new refinement method called YUP.SCX. The electron density map is converted into a component of the potential energy function to which terms for stereochemical restraints and volume exclusion are added. The potential energy function is then minimized (using simulated annealing) to yield a stereochemically-restrained atomic structure that fits into the electron density map optimally. We used this procedure to construct an atomic model of the 70S ribosome in the pre-accommodation state. Although some atoms are displaced by as much as 33 Å, they divide themselves into nearly rigid fragments along natural boundaries with smooth transitions between the fragments.     

Creating protein models from electron-density maps using particle-filtering methods

MOTIVATION: One bottleneck in high-throughput protein crystallography is interpreting an electron-density map, that is, fitting a molecular model to the 3D picture crystallography produces. Previously, we developed ACMI (Automatic Crystallographic Map Interpreter), an algorithm that uses a probabilistic model to infer an accurate protein backbone layout. Here, we use a sampling method known as particle filtering to produce a set of all-atom protein models. We use the output of ACMI to guide the particle filter's sampling, producing an accurate, physically feasible set of structures. RESULTS: We test our algorithm on 10 poor-quality experimental density maps. We show that particle filtering produces accurate all-atom models, resulting in fewer chains, lower sidechain RMS error and reduced R factor, compared to simply placing the best-matching sidechains on ACMI's trace. We show that our approach produces a more accurate model than three leading methods--Textal, Resolve and ARP/WARP--in terms of main chain completeness, sidechain identification and crystallographic R factor. AVAILABILITY: Source code and experimental density maps available at http://ftp.cs.wisc.edu/machine-learning/shavlik-group/programs/acmi/     

Structure of wheat serine carboxypeptidase II at 3.5-A resolution. A new class of serine proteinase

The structure of serine carboxypeptidase II from wheat bran has been determined to 3.5-A resolution by multiple isomorphous replacement, solvent flattening, and crystallographic refinement. The amino acid sequence has been fit to the electron density map and the model refined to a conventional crystallographic R factor of 20.9%. The molecule is an alpha + beta protein and contains a "catalytic triad" (Asp338, His397, and Ser146) similar in arrangement to those in chymotrypsin and subtilisin. The -fold of the polypeptide backbone is, however, completely different from those enzymes. This suggests that this is a third example of convergent evolution to a common enzymatic mechanism. The -fold is, on the other hand, surprisingly similar to that of the zinc proteinase carboxypeptidase A.     

Structural characterization of components of protein assemblies by comparative modeling and electron cryo-microscopy

We explore structural characterization of protein assemblies by a combination of electron cryo-microscopy (cryoEM) and comparative protein structure modeling. Specifically, our method finds an optimal atomic model of a given assembly subunit and its position within an assembly by fitting alternative comparative models into a cryoEM map. The alternative models are calculated by MODELLER [J. Mol. Biol. 234 (1993) 313] from different sequence alignments between the modeled protein and its template structures. The fitting of these models into a cryoEM density map is performed either by FOLDHUNTER [J. Mol. Biol. 308 (2001) 1033] or by a new density fitting module of MODELLER (Mod-EM). Identification of the most accurate model is based on the correlation between the model accuracy and the quality of fit into the cryoEM density map. To quantify this correlation, we created a benchmark consisting of eight proteins of different structural folds with corresponding density maps simulated at five resolutions from 5 to 15 angstroms, with three noise levels each. Each of the proteins in the set was modeled based on 300 different alignments to their remotely related templates (12-32% sequence identity), spanning the range from entirely inaccurate to essentially accurate alignments. The benchmark revealed that one of the most accurate models can usually be identified by the quality of its fit into the cryoEM density map, even for noisy maps at 15 angstroms resolution. Therefore, a cryoEM density map can be helpful in improving the accuracy of a comparative model. Moreover, a pseudo-atomic model of a component in an assembly may be built better with comparative models of the native subunit sequences than with experimentally determined structures of their homologs.     

Low-resolution structure refinement in electron microscopy

A real-space structure refinement method, originally developed for macromolecular X-ray crystallography, has been applied to protein structure analysis by electron microscopy (EM). This method simultaneously optimizes the fit of an atomic model to a density map and the stereo-chemical properties of the model by minimizing an energy function. The performance of this method is characterized at different resolution and signal-to-noise ratio conditions typical for EM electron density maps. A multi-resolution scheme is devised to improve the convergence of the refinement on the global energy minimum. Applications of the method to various model systems are demonstrated here. The first case is the arrangement of FlgE molecules in the helical filament of flagellar hook, in which refinement with segmented rigid bodies improves the density correlation and reduces severe van der Waals contacts among the symmetry-related subunits. The second case is a conformational analysis of the NSF AAA ATPase in which a multi-conformer model is used in the refinement to investigate the arrangement of the two ATPase domains in the molecule. The third case is a docking simulation in which the crystal structure of actin and the NOE data from NMR experiments on the dematin headpiece are combined with a low-resolution EM density map to generate an atomic model of the F-actin-dematin headpiece structure.     

Structure of ferricytochrome c' from Rhodospirillum rubrum at 6 A resolution

The structure of a ferricytochrome c' extracted from Rhodospirillum rubrum has been determined at 6 A resolution by the X-ray crystallographic method. The crystals, obtained by dialyzing the protein solution against polyethylene glycol 4000, belong to the hexagonal space group P6(1). Two heavy atom derivatives were obtained by soaking the native crystals in K2PtCl6 and CH3HgCl solution. The phases calculated by the multiple isomorphous replacement method gave an overall figure of merit of 0.90 at 6 A resolution. The resulting electron density map showed the molecular boundary clearly, and gave molecular dimensions of 50 X 25 X 30 A for a monomer molecule. From visual examination of this map, the cytochrome c' from Rhodospirillum rubrum has a similar chain-folding pattern to the cytochrome c' from Rhodospirillum molischianum, the structure determination of which has already been carried out.     

Crystallization and initial crystallographic results for pepstatin A inhibited bovine cathepsin D.

Cathepsin D was purified from bovine liver by a method using two pepstatin A affinity columns. The eluted protein was combined with pepstatin A and the complex crystallized from 15% polyethylene glycol 8000 at pH 5.9. The crystals diffract to a resolution of 3.0 A and have space group P2(1)2(1)2(1) with unit cell dimensions a = 74.8 A, b = 76.0 A, c = 157.7 A. There are two molecules in the asymmetric unit. The structure was solved by molecular replacement using a pepsin search model and both molecules showed clearly interpretable density in the position expected for pepstatin A in a preliminary difference map. The refined model has r.m.s. deviations from ideal bond lengths and angles of 0.014 A and 3.2 degrees, respectively, and a crystallographic R factor of 17%.     

Three-dimensional structure of recombinant human muscle fatty acid-binding protein.

The three-dimensional structure of recombinant human muscle fatty acid-binding protein with a bound fatty acid has been solved and refined with x-ray diffraction data to 2.1 A resolution. The refined model has a crystallographic R factor of 19.5% for data between 9.0 and 2.1 A (7243 unique reflections) and root-mean-square deviations in bond length and bond angle of 0.013 A and 2.7 degrees. The protein contains 10 antiparallel beta-strands and two short alpha-helices which are arranged into two approximately orthogonal beta-sheets. Difference electron density maps and a multiple isomorphous derivative electron density map showed the presence of a single bound molecule of a long chain fatty acid within the interior core of the protein. The hydrocarbon tail of the fatty acid was found to be in a "U-shaped" conformation. Seven ordered water molecules were also identified within the interior of the protein in a pocket on the pseudo-si face of the fatty acid's bent hydrocarbon tail. The methylene tail of the fatty acid forms van der Waals interactions with atoms from 13 residues and three ordered waters. The carboxylate of the fatty acid is located in the interior of the protein where it forms hydrogen bonds with the side chains of Tyr128 and Arg126 and two ordered water molecules. A comparison of the three-dimensional structure of human muscle fatty acid-binding protein and rat intestinal fatty acid-binding protein shows strong similarity. Both proteins bind a single fatty acid within their interior cores, but the bound fatty acids are very different in their conformations and interactions. These findings suggest that the intestinal and muscle fatty acid-binding proteins have evolved distinct binding sites in order to satisfy different requirements within the tissues where they are expressed.