The amounts of thermal vibrations and static disorder in protein X-ray crystallographic B-factors.

Crystallographic B-factors provide direct dynamical information on the internal mobility of proteins that is closely linked to function, and are also widely used as a benchmark in assessing elastic network models. A significant question in the field is: what is the exact amount of thermal vibrations in protein crystallographic B-factors? This work sets out to answer this question. First, we carry out a thorough, statistically sound analysis of crystallographic B-factors of over 10 000 structures. Second, by employing a highly accurate all-atom model based on the well-known CHARMM force field, we obtain computationally the magnitudes of thermal vibrations of nearly 1000 structures. Our key findings are: (i) the magnitude of thermal vibrations, surprisingly, is nearly protein-independent, as a corollary to the universality for the vibrational spectra of globular proteins established earlier; (ii) the magnitude of thermal vibrations is small, less than 0.1 Å2 at 100 K; (iii) the percentage of thermal vibrations in B-factors is the lowest at low resolution and low temperature (<10%) but increases to as high as 60% for structures determined at high resolution and at room temperature. The significance of this work is that it provides for the first time, using an extremely large dataset, a thorough analysis of B-factors and their thermal and static disorder components. The results clearly demonstrate that structures determined at high resolution and at room temperature have the richest dynamics information. Since such structures are relatively rare in the PDB database, the work naturally calls for more such structures to be determined experimentally.

An integrated pipeline for sample preparation and characterization at the EMBL@PETRA3 synchrotron facilities.

The characterization of macromolecular samples at synchrotrons has traditionally been restricted to direct exposure to X-rays, but beamline automation and diversification of the user community has led to the establishment of complementary characterization facilities off-line. The Sample Preparation and Characterization (SPC) facility at the EMBL@PETRA3 synchrotron provides synchrotron users access to a range of biophysical techniques for preliminary or parallel sample characterization, to optimize sample usage at the beamlines. Here we describe a sample pipeline from bench to beamline, to assist successful structural characterization using small angle X-ray scattering (SAXS) or macromolecular X-ray crystallography (MX). The SPC has developed a range of quality control protocols to assess incoming samples and to suggest optimization protocols. A high-throughput crystallization platform has been adapted to reach a broader user community, to include chemists and biologists that are not experts in structural biology. The SPC in combination with the beamline and computational facilities at EMBL Hamburg provide a full package of integrated facilities for structural biology and can serve as model for implementation of such resources for other infrastructures.

How precise are reported protein coordinate data?

Atomic coordinates in the Worldwide Protein Data Bank (wwPDB) are generally reported to greater precision than the experimental structure determinations have actually achieved. By using information theory and data compression to study the compressibility of protein atomic coordinates, it is possible to quantify the amount of randomness in the coordinate data and thereby to determine the realistic precision of the reported coordinates. On average, the value of each C(α) coordinate in a set of selected protein structures solved at a variety of resolutions is good to about 0.1 Å.

Molecular replacement: tricks and treats.

Molecular replacement is the method of choice for X-ray crystallographic structure determination provided that suitable structural homologues are available in the PDB. Presently, there are ~80,000 structures in the PDB (8074 were deposited in the year 2012 alone), of which ~70% have been solved by molecular replacement. For successful molecular replacement the model must cover at least 50% of the total structure and the Cα r.m.s.d. between the core model and the structure to be solved must be less than 2 Å. Here, an approach originally implemented in the CaspR server (http://www.igs.cnrs-mrs.fr/Caspr2/index.cgi) based on homology modelling to search for a molecular-replacement solution is discussed. How the use of as much information as possible from different sources can improve the model(s) is briefly described. The combination of structural information with distantly related sequences is crucial to optimize the multiple alignment that will define the boundaries of the core domains. PDB clusters (sequences with ≥30% identical residues) can also provide information on the eventual changes in conformation and will help to explore the relative orientations assumed by protein subdomains. Normal-mode analysis can also help in generating series of conformational models in the search for a molecular-replacement solution. Of course, finding a correct solution is only the first step and the accuracy of the identified solution is as important as the data quality to proceed through refinement. Here, some possible reasons for failure are discussed and solutions are proposed using a set of successful examples.

Techniques, tools and best practices for ligand electron-density analysis and results from their application to deposited crystal structures.

As a result of substantial instrumental automation and the continuing improvement of software, crystallographic studies of biomolecules are conducted by non-experts in increasing numbers. While improved validation almost ensures that major mistakes in the protein part of structure models are exceedingly rare, in ligand-protein complex structures, which in general are most interesting to the scientist, ambiguous ligand electron density is often difficult to interpret and the modelled ligands are generally more difficult to properly validate. Here, (i) the primary technical reasons and potential human factors leading to problems in ligand structure models are presented; (ii) the most common categories of building errors or overinterpretation are classified; (iii) a few instructive and specific examples are discussed in detail, including an electron-density-based analysis of ligand structures that do not contain any ligands; (iv) means of avoiding such mistakes are suggested and the implications for database validity are discussed and (v) a user-friendly software tool that allows non-expert users to conveniently inspect ligand density is provided.

Localization and orientation of heavy-atom cluster compounds in protein crystals using molecular replacement.

Heavy-atom clusters (HA clusters) containing a large number of specifically arranged electron-dense scatterers are especially useful for experimental phase determination of large complex structures, weakly diffracting crystals or structures with large unit cells. Often, the determination of the exact orientation of the HA cluster and hence of the individual heavy-atom positions proves to be the critical step in successful phasing and subsequent structure solution. Here, it is demonstrated that molecular replacement (MR) with either anomalous or isomorphous differences is a useful strategy for the correct placement of HA cluster compounds. The polyoxometallate cluster hexasodium α-metatungstate (HMT) was applied in phasing the structure of death receptor 6. Even though the HA cluster is bound in alternate partially occupied orientations and is located at a special position, its correct localization and orientation could be determined at resolutions as low as 4.9 Å. The broad applicability of this approach was demonstrated for five different derivative crystals that included the compounds tantalum tetradecabromide and trisodium phosphotungstate in addition to HMT. The correct placement of the HA cluster depends on the length of the intramolecular vectors chosen for MR, such that both a larger cluster size and the optimal choice of the wavelength used for anomalous data collection strongly affect the outcome.

A quality comparison of protein crystals grown under containerless conditions generated by diamagnetic levitation, silicone oil and agarose gel.

High-quality crystals are key to obtaining accurate three-dimensional structures of proteins using X-ray diffraction techniques. However, obtaining such protein crystals is often a challenge. Several containerless crystallization techniques have been reported to have the ability to improve crystal quality, but it is unknown which is the most favourable way to grow high-quality protein crystals. In this paper, a quality comparison of protein crystals which were grown under three containerless conditions provided by diamagnetic levitation, silicone oil and agarose gel was conducted. A control experiment on a vessel wall was also simultaneously carried out. Seven different proteins were crystallized under the four conditions, and the crystal quality was assessed in terms of the resolution limit, the mosaicity and the Rmerge. It was found that the crystals grown under the three containerless conditions demonstrated better morphology than those of the control. X-ray diffraction data indicated that the quality of the crystals grown under the three containerless conditions was better than that of the control. Of the three containerless crystallization techniques, the diamagnetic levitation technique exhibited the best performance in enhancing crystal quality. This paper is to our knowledge the first report of improvement of crystal quality using a diamagnetic levitation technique. Crystals obtained from agarose gel demonstrated the second best improvement in crystal quality. The study indicated that the diamagnetic levitation technique is indeed a favourable method for growing high-quality protein crystals, and its utilization is thus potentially useful in practical efforts to obtain well diffracting protein crystals.

On optimal placement of molecules in the unit cell.

There are currently no rules for a unified, standard way of placing macromolecular structures in the crystal lattice. An analysis of all possible symmetry-equivalent representations of molecular structures in various space groups leads to the concept of the anti-Cheshire symmetry and suggests that the center of a unique structural motif can always be placed within the selected asymmetric unit of the anti-Cheshire cell. The placement of structures according to this suggestion will ensure uniformity of presentation of all structurally equivalent Protein Data Bank models and will therefore diminish the possibility of confusing less crystallographically knowledgeable users of the PDB. The anti-Cheshire cells and their asymmetric units are defined and tabulated for all 65 space groups relevant to macromolecular crystallography that exhibit only rotational symmetry operations.

High-resolution structure prediction and the crystallographic phase problem.

The energy-based refinement of low-resolution protein structure models to atomic-level accuracy is a major challenge for computational structural biology. Here we describe a new approach to refining protein structure models that focuses sampling in regions most likely to contain errors while allowing the whole structure to relax in a physically realistic all-atom force field. In applications to models produced using nuclear magnetic resonance data and to comparative models based on distant structural homologues, the method can significantly improve the accuracy of the structures in terms of both the backbone conformations and the placement of core side chains. Furthermore, the resulting models satisfy a particularly stringent test: they provide significantly better solutions to the X-ray crystallographic phase problem in molecular replacement trials. Finally, we show that all-atom refinement can produce de novo protein structure predictions that reach the high accuracy required for molecular replacement without any experimental phase information and in the absence of templates suitable for molecular replacement from the Protein Data Bank. These results suggest that the combination of high-resolution structure prediction with state-of-the-art phasing tools may be unexpectedly powerful in phasing crystallographic data for which molecular replacement is hindered by the absence of sufficiently accurate previous models.

The structures of T6, T3R3 and R6 bovine insulin: combining X-ray diffraction and absorption spectroscopy.

The crystal structures of three conformations, T(6), T(3)R(3) and R(6), of bovine insulin were solved at 1.40, 1.30 and 1.80 Å resolution, respectively. All conformations crystallized in space group R3. In contrast to the T(6) and T(3)R(3) structures, different conformations of the N-terminal B-chain residue PheB1 were observed in the R(6) insulin structure, resulting in an eightfold doubling of the unit-cell volume upon cooling. The zinc coordination in each conformation was studied by X-ray absorption spectroscopy (XAS), including both EXAFS and XANES. Zinc adopts a tetrahedral coordination in all R(3) sites and an octahedral coordination in T(3) sites. The coordination distances were refined from XAS with a standard deviation of <0.01 Å. In contrast to the distances determined from the medium-resolution crystal structures, the XAS results were in good agreement with similar coordination geometries found in small molecules, as well as in other high-resolution insulin structures. As the radiation dose for XRD experiments is two orders of magnitude higher compared with that of XAS experiments, the single crystals were exposed to a higher degree of radiation damage that affected the zinc coordination in the T(3) sites in particular. Furthermore, XANES spectra for the zinc sites in T(6) and R(6) insulin were successfully calculated using finite difference methods and the bond distances and angles were optimized from a quantitative XANES analysis.

Physico-mechanical and stability evaluation of carbamazepine cocrystal with nicotinamide.

The focus of this investigation was to prepare the cocrystal of carbamazepine (CBZ) using nicotinamide as a coformer and to compare its preformulation properties and stability profile with CBZ. The cocrystal was prepared by solution cooling crystallization, solvent evaporation, and melting and cryomilling methods. They were characterized for solubility, intrinsic dissolution rate, chemical identification by Fourier transform infrared spectroscopy, crystallinity by differential scanning calorimetry, powder X-ray diffraction, and morphology by scanning electron microscopy. Additionally, mechanical properties were evaluated by tensile strength and Heckel analysis of compacts. The cocrystal and CBZ were stored at 40°C/94% RH, 40°C/75% RH, 25°C/60% RH, and 60°C to determine their stability behavior. The cocrystals were fluffy, with a needle-shaped crystal, and were less dense than CBZ. The solubility profiles of the cocrystals were similar to CBZ, but its intrinsic dissolution rate was lower due to the high tensile strength of its compacts. Unlike CBZ, the cocrystals were resistant to hydrate transformation, as revealed by the stability studies. Plastic deformation started at a higher compression pressure in the cocrystals than CBZ, as indicated by the high yield pressure. In conclusion, the preformulation profile of the cocrystals was similar to CBZ, except that it had an advantageous resistance to hydrate transformation.

PAIREF: paired refinement also for Phenix users.

In macromolecular crystallography, paired refinement is generally accepted to be the optimal approach for the determination of the high-resolution cutoff. The software tool PAIREF provides automation of the protocol and associated analysis. Support for phenix.refine as a refinement engine has recently been implemented in the program. This feature is presented here using previously published data for thermolysin. The results demonstrate the importance of the complete cross-validation procedure to obtain a thorough and unbiased insight into the quality of high-resolution data.

Optimizing the Growth of Endothiapepsin Crystals for Serial Crystallography Experiments.

Here, a protocol is presented to facilitate the creation of large volumes (> 100 µL) of micro-crystalline slurries suitable for serial crystallography experiments at both synchrotrons and XFELs. The method is based upon an understanding of the protein crystal phase diagram, and how that knowledge can be utilized. The method is divided into three stages: (1) optimizing crystal morphology, (2) transitioning to batch, and (3) scaling. Stage 1 involves finding well diffracting, single crystals, hopefully but not necessarily, presenting in a cube-like morphology. In Stage 2, the Stage 1 condition is optimized by crystal growth time. This strategy can transform crystals grown by vapor diffusion to batch. Once crystal growth can occur within approximately 24 h, a morphogram of the protein and precipitant mixture can be plotted and used as the basis for a scaling strategy (Stage 3). When crystals can be grown in batch, scaling can be attempted, and the crystal size and concentration optimized as the volume is increased. Endothiapepsin has been used as a demonstration protein for this protocol. Some of the decisions presented are specific to endothiapepsin. However, it is hoped that the way they have been applied will inspire a way of thinking about this procedure that others can adapt to their own projects.

Characterization of wax as a potential diffraction intensity standard for macromolecular crystallography beamlines.

A number of commercially available waxes in the form of thin disc samples have been investigated as possible diffraction intensity standards for macromolecular crystallography synchrotron beamlines. Synchrotron X-ray powder diffraction measurements show that beeswax offers the best performance of these waxes owing to its polycrystallinity. Crystallographic lattice parameters and diffraction intensities were examined between 281 and 309 K, and show stable and predictable thermal behaviour. Using an X-ray beam of known incident flux at lambda = 1 A, the diffraction power of two strong Bragg reflections for beeswax were quantified as a function of sample thickness and normalized to 10(10) photons s(-1). To demonstrate its feasibility as a diffraction intensity standard, test measurements were then performed on a new third-generation macromolecular crystallography synchrotron beamline.

Detecting the nature and solving the crystal structure of a contaminant protein from an opportunistic pathogen.

The unintentional crystallization of contaminant proteins in the place of target recombinant proteins is sporadically reported, despite the availability of stringent expression/purification protocols and of software for the detection of contaminants. Typically, the contaminant protein originates from the expression organism (for example Escherichia coli), but in rare circumstances contaminants from different sources have been reported. Here, a case of contamination from a Serratia bacterial strain that occurred while attempting to crystallize an unrelated protein from Burkholderia pseudomallei (overexpressed in E. coli) is presented. The contamination led to the unintended crystallization and structure analysis of a cyanase hydratase from a bacterial strain of the Serratia genus, an opportunistic enterobacterium that grows under conditions similar to those of E. coli and that is found in a variety of habitats, including the laboratory environment. In this context, the procedures that were adopted to identify the contaminant based on crystallographic data only are presented and the crystal structure of Serrata spp. cyanase hydratase is briefly discussed.