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

Obtaining anomalous and ensemble information from protein crystals from 220 K up to physiological temperatures

X-ray crystallography has been invaluable in delivering structural information about proteins. Previously, an approach has been developed that allows high-quality X-ray diffraction data to be obtained from protein crystals at and above room temperature. Here, this previous work is built on and extended by showing that high-quality anomalous signal can be obtained from single protein crystals using diffraction data collected at 220 K up to physiological temperatures. The anomalous signal can be used to directly determine the structure of a protein, i.e. to phase the data, as is routinely performed under cryoconditions. This ability is demonstrated by obtaining diffraction data from model lysozyme, thaumatin and proteinase K crystals, the anomalous signal from which allowed their structures to be solved experimentally at 7.1 keV X-ray energy and at room temperature with relatively low data redundancy. It is also demonstrated that the anomalous signal from diffraction data obtained at 310 K (37°C) can be used to solve the structure of proteinase K and to identify ordered ions. The method provides useful anomalous signal at temperatures down to 220 K, resulting in an extended crystal lifetime and increased data redundancy. Finally, we show that useful anomalous signal can be obtained at room temperature using X-rays of 12 keV energy as typically used for routine data collection, allowing this type of experiment to be carried out at widely accessible synchrotron beamline energies and enabling the simultaneous extraction of high-resolution data and anomalous signal. With the recent emphasis on obtaining conformational ensemble information for proteins, the high resolution of the data allows such ensembles to be built, while the anomalous signal allows the structure to be experimentally solved, ions to be identified, and water molecules and ions to be differentiated. Because bound metal-, phosphorus- and sulfur-containing ions all have anomalous signal, obtaining anomalous signal across temperatures and up to physiological temperatures will provide a more complete description of protein conformational ensembles, function and energetics.

CMM-An enhanced platform for interactive validation of metal binding sites

Metal ions bound to macromolecules play an integral role in many cellular processes. They can directly participate in catalytic mechanisms or be essential for the structural integrity of proteins and nucleic acids. However, their unique nature in macromolecules can make them difficult to model and refine, and a substantial portion of metal ions in the PDB are misidentified or poorly refined. CheckMyMetal (CMM) is a validation tool that has gained widespread acceptance as an essential tool for researchers working on metal-macromolecule complexes. CMM can be used during structure determination or to validate metal binding sites in structural models within the PDB. The functionalities of CMM have recently been greatly enhanced and provide researchers with additional information that can guide modeling decisions. The new version of CMM shows metals in the context of electron density maps and allows for on-the-fly refinement of metal binding sites. The improvements should increase the reproducibility of biomedical research. The web server is available at https://cmm.minorlab.org.

Graph-Convolutional Neural Net Model of the Statistical Torsion Profiles for Small Organic Molecules

We present a graph-convolutional neural network (GCNN)-based method for learning and prediction of statistical torsional profiles (STP) in small organic molecules based on the experimental X-ray structure data. A specialized GCNN torsion profile model is trained using the structures in the Crystallography Open Database (COD). The GCNN-STP model captures torsional preferences over a wide range of torsion rotor chemotypes and correctly predicts a variety of effects from the vicinal atoms and moieties. GCNN-STP statistical profiles also show good agreement with quantum chemically (DFT) calculated torsion energy profiles. Furthermore, we demonstrate the application of the GCNN-STP statistical profiles for conformer generation. A web server that allows interactive profile prediction and viewing is made freely available at https://www.molsoft.com/tortool.html.

Structure-aided optimization of non-nucleoside M. tuberculosis thymidylate kinase inhibitors

Mycobacterium tuberculosis thymidylate kinase (MtTMPK) has emerged as an attractive target for rational drug design. We recently investigated new families of non-nucleoside MtTMPK inhibitors in an effort to diversify MtTMPK inhibitor chemical space. We here report a new series of MtTMPK inhibitors by combining the Topliss scheme with rational drug design approaches, fueled by two co-crystal structures of MtTMPK in complex with developed inhibitors. These efforts furnished the most potent MtTMPK inhibitors in our assay, with two analogues displaying low micromolar MIC values against H37Rv Mtb. Prepared inhibitors address new sub-sites in the MtTMPK nucleotide binding pocket, thereby offering new insights into its druggability. We studied the role of efflux pumps as well as the impact of cell wall permeabilizers for selected compounds to potentially provide an explanation for the lack of correlation between potent enzyme inhibition and whole-cell activity.

Temperature artifacts in protein structures bias ligand-binding predictions

X-ray crystallography is the gold standard to resolve conformational ensembles that are significant for protein function, ligand discovery, and computational methods development. However, relevant conformational states may be missed at common cryogenic (cryo) data-collection temperatures but can be populated at room temperature. To assess the impact of temperature on making structural and computational discoveries, we systematically investigated protein conformational changes in response to temperature and ligand binding in a structural and computational workhorse, the T4 lysozyme L99A cavity. Despite decades of work on this protein, shifting to RT reveals new global and local structural changes. These include uncovering an apo helix conformation that is hidden at cryo but relevant for ligand binding, and altered side chain and ligand conformations. To evaluate the impact of temperature-induced protein and ligand changes on the utility of structural information in computation, we evaluated how temperature can mislead computational methods that employ cryo structures for validation. We find that when comparing simulated structures just to experimental cryo structures, hidden successes and failures often go unnoticed. When using structural information in ligand binding predictions, both coarse docking and rigorous binding free energy calculations are influenced by temperature effects. The trend that cryo artifacts limit the utility of structures for computation holds across five distinct protein classes. Our results suggest caution when consulting cryogenic structural data alone, as temperature artifacts can conceal errors and prevent successful computational predictions, which can mislead the development and application of computational methods in discovering bioactive molecules.

Structural Analysis of Strigolactone-Related Gene Products

Structural knowledge of biological macromolecules is essential for understanding their function and for modifying that function by engineering. Protein crystallography is a powerful method for elucidating molecular structures of proteins, but it is essential that the investigator has a basic knowledge of good practices and of the major pitfalls in the technique. Here we describe issues specific for the case of structural studies of strigolactone (SL) receptor structure and function, and in particular the difficulties associated with capturing complexes of SL receptors with the SL hormone ligand in the crystal.

Recognizing and validating ligands with CheckMyBlob

Structure-guided drug design depends on the correct identification of ligands in crystal structures of protein complexes. However, the interpretation of the electron density maps is challenging and often burdened with confirmation bias. Ligand identification can be aided by automatic methods such as CheckMyBlob, a machine learning algorithm that learns to generalize ligand descriptions from sets of moieties deposited in the Protein Data Bank. Here, we present the CheckMyBlob web server, a platform that can identify ligands in unmodeled fragments of electron density maps or validate ligands in existing models. The server processes PDB/mmCIF and MTZ files and returns a ranking of 10 most likely ligands for each detected electron density blob along with interactive 3D visualizations. Additionally, for each prediction/validation, a plugin script is generated that enables users to conduct a detailed analysis of the server results in Coot. The CheckMyBlob web server is available at https://checkmyblob.bioreproducibility.org.

'All That Glitters Is Not Gold': High-Resolution Crystal Structures of Ligand-Protein Complexes Need Not Always Represent Confident Binding Poses

Our understanding of the structure–function relationships of biomolecules and thereby applying it to drug discovery programs are substantially dependent on the availability of the structural information of ligand–protein complexes. However, the correct interpretation of the electron density of a small molecule bound to a crystal structure of a macromolecule is not trivial. Our analysis involving quality assessment of ~0.28 million small molecule–protein binding site pairs derived from crystal structures corresponding to ~66,000 PDB entries indicates that the majority (65%) of the pairs might need little (54%) or no (11%) attention. Out of the remaining 35% of pairs that need attention, 11% of the pairs (including structures with high/moderate resolution) pose serious concerns. Unfortunately, most users of crystal structures lack the training to evaluate the quality of a crystal structure against its experimental data and, in general, rely on the resolution as a ‘gold standard’ quality metric. Our work aims to sensitize the non-crystallographers that resolution, which is a global quality metric, need not be an accurate indicator of local structural quality. In this article, we demonstrate the use of several freely available tools that quantify local structural quality and are easy to use from a non-crystallographer’s perspective. We further propose a few solutions for consideration by the scientific community to promote quality research in structural biology and applied areas.

Crystallographic models of SARS-CoV-2 3CLpro: in-depth assessment of structure quality and validation

The appearance at the end of 2019 of the new SARS-CoV-2 coronavirus led to an unprecedented response by the structural biology community, resulting in the rapid determination of many hundreds of structures of proteins encoded by the virus. As part of an effort to analyze and, if necessary, remediate these structures as deposited in the Protein Data Bank (PDB), this work presents a detailed analysis of 81 crystal structures of the main protease 3CLpro, an important target for the design of drugs against COVID-19. The structures of the unliganded enzyme and its complexes with a number of inhibitors were determined by multiple research groups using different experimental approaches and conditions; the resulting structures span 13 different polymorphs representing seven space groups. The structures of the enzyme itself, all determined by molecular replacement, are highly similar, with the exception of one polymorph with a different inter-domain orientation. However, a number of complexes with bound inhibitors were found to pose significant problems. Some of these could be traced to faulty definitions of geometrical restraints for ligands and to the general problem of a lack of such information in the PDB depositions. Several problems with ligand definition in the PDB itself were also noted. In several cases extensive corrections to the models were necessary to adhere to the evidence of the electron-density maps. Taken together, this analysis of a large number of structures of a single, medically important protein, all determined within less than a year using modern experimental tools, should be useful in future studies of other systems of high interest to the biomedical community.

Synchrotron Radiation as a Tool for Macromolecular X-Ray Crystallography: a XXI Century Perspective

Intense X-rays available at powerful synchrotron beamlines provide macromolecular crystallographers with an incomparable tool for investigating biological phenomena on an atomic scale. The resulting insights into the mechanism's underlying biological processes have played an essential role and shaped biomedical sciences during the last 30 years, considered the "golden age" of structural biology. In this review, we analyze selected aspects of the impact of synchrotron radiation on structural biology. Synchrotron beamlines have been used to determine over 70% of all macromolecular structures deposited into the Protein Data Bank (PDB). These structures were deposited by over 13,000 different research groups. Interestingly, despite the impressive advances in synchrotron technologies, the median resolution of macromolecular structures determined using synchrotrons has remained constant throughout the last 30 years, at about 2 Å. Similarly, the median times from the data collection to the deposition and release have not changed significantly. We describe challenges to reproducibility related to recording all relevant data and metadata during the synchrotron experiments, including diffraction images. Finally, we discuss some of the recent opinions suggesting a diminishing importance of X-ray crystallography due to impressive advances in Cryo-EM and theoretical modeling. We believe that synchrotrons of the future will increasingly evolve towards a life science center model, where X-ray crystallography, Cryo-EM, and other experimental and computational resources and knowledge are encompassed within a versatile research facility. The recent response of crystallographers to the COVID-19 pandemic suggests that X-ray crystallography conducted at synchrotron beamlines will continue to play an essential role in structural biology and drug discovery for years to come.

Ligand Incorporation into Protein Microcrystals for MicroED by On-Grid Soaking

A high throughout method for soaking ligands into protein microcrystals on TEM grids is presented. Every crystal on the grid is soaked simultaneously using only standard cryoelectron microscopy vitrification equipment. The method is demonstrated using proteinase K microcrystals soaked with the 5-amino-2,4,6-triodoisophthalic acid (I3C) magic triangle. A soaked microcrystal is milled to a thickness of approximately 200 nm using a focused ion beam, and MicroED data are collected. A high-resolution structure of the protein with four ligands at high occupancy is determined. Both the number of ligands bound and their occupancy is higher using on-grid soaking of microcrystals compared with much larger crystals treated similarly and investigated by X-ray crystallography. These results indicate that on-grid soaking ligands into microcrystals results in efficient uptake of ligands into protein microcrystals.

Exploring Ligand Stability in Protein Crystal Structures Using Binding Pose Metadynamics

Identification of correct protein–ligand binding poses is important in structure-based drug design and crucial for the evaluation of protein–ligand binding affinity. Protein–ligand coordinates are commonly obtained from crystallography experiments that provide a static model of an ensemble of conformations. Binding pose metadynamics (BPMD) is an enhanced sampling method that allows for an efficient assessment of ligand stability in solution. Ligand poses that are unstable under the bias of the metadynamics simulation are expected to be infrequently occupied in the energy landscape, thus making minimal contributions to the binding affinity. Here, the robustness of the method is studied using crystal structures with ligands known to be incorrectly modeled, as well as 63 structurally diverse crystal structures with ligand fit to electron density from the Twilight database. Results show that BPMD can successfully differentiate compounds whose binding pose is not supported by the electron density from those with well-defined electron density.

Automatic recognition of ligands in electron density by machine learning

Motivation

The correct identification of ligands in crystal structures of protein complexes is the cornerstone of structure-guided drug design. However, cognitive bias can sometimes mislead investigators into modeling fictitious compounds without solid support from the electron density maps. Ligand identification can be aided by automatic methods, but existing approaches are based on time-consuming iterative fitting.

Results

Here we report a new machine learning algorithm called CheckMyBlob that identifies ligands from experimental electron density maps. In benchmark tests on portfolios of up to 219 931 ligand binding sites containing the 200 most popular ligands found in the Protein Data Bank, CheckMyBlob markedly outperforms the existing automatic methods for ligand identification, in some cases doubling the recognition rates, while requiring significantly less time. Our work shows that machine learning can improve the automation of structure modeling and significantly accelerate the drug screening process of macromolecule-ligand complexes.

Availability and implementation

Code and data are available on GitHub at https://github.com/dabrze/CheckMyBlob.

Supplementary information

Supplementary data are available at Bioinformatics online.

Analysis and validation of macromolecular B values

This paper describes a global analysis of macromolecular B values. It is shown that the distribution of B values generally follows the shifted inverse-gamma distribution (SIGD). The parameters of the SIGD are estimated using the Fisher scoring technique with the expected Fisher information matrix. It is demonstrated that a contour plot based on the parameters of the SIGD can play a role in the validation of macromolecular structures. The dependence of the peak-height distribution on resolution and atomic B values is also analysed. It is demonstrated that the B-value distribution can have a dramatically different effect on peak heights at different resolutions. Consequently, a comparative analysis of the B values of neighbouring atoms must account for resolution. A combination of the SIGD, peak-height distribution and outlier detection was used to identify a number of entries from the PDB that require attention. It is also shown that the presence of a multimodal B-value distribution often indicates that some loops or parts of the molecule have either been mismodelled or have dramatically different mobility, depending on their environment within the crystal. These distributions can also indicate the level of sharpening/blurring used before atomic structure refinement. It is recommended that procedures such as sharpening/blurring should be avoided during refinement, although they can play important roles in map visualization and model building.

Ice formation and solvent nanoconfinement in protein crystals

Ice formation within protein crystals is a major obstacle to the cryocrystallographic study of protein structure, and has limited studies of how the structural ensemble of a protein evolves with temperature in the biophysically interesting range from ∼260 K to the protein-solvent glass transition near 200 K. Using protein crystals with solvent cavities as large as ∼70 Å, time-resolved X-ray diffraction was used to study the response of protein and internal solvent during rapid cooling. Solvent nanoconfinement suppresses freezing temperatures and ice-nucleation rates so that ice-free, low-mosaicity diffraction data can be reliably collected down to 200 K without the use of cryoprotectants. Hexagonal ice (Ih) forms in external solvent, but internal crystal solvent forms stacking-disordered ice (Isd) with a near-random stacking of cubic and hexagonal planes. Analysis of powder diffraction from internal ice and single-crystal diffraction from the host protein structure shows that the maximum crystallizable solvent fraction decreases with decreasing crystal solvent-cavity size, and that an ∼6 Å thick layer of solvent adjacent to the protein surface cannot crystallize. These results establish protein crystals as excellent model systems for the study of nanoconfined solvent. By combining fast cooling, intense X-ray beams and fast X-ray detectors, complete structural data sets for high-value targets, including membrane proteins and large complexes, may be collected at ∼220-240 K that have much lower mosaicities and comparable B factors, and that may allow more confident identification of ligand binding than in current cryocrystallographic practice.

Nucleobase carbonyl groups are poor Mg2+ inner-sphere binders but excellent monovalent ion binders-a critical PDB survey

Precise knowledge of Mg²⁺ inner-sphere binding site properties is vital for understanding the structure and function of nucleic acid systems. Unfortunately, the PDB, which represents the main source of Mg²⁺ binding sites, contains a substantial number of assignment issues that blur our understanding of the functions of these ions. Here, following a previous study devoted to Mg²⁺ binding to nucleobase nitrogens, we surveyed nucleic acid X-ray structures from the PDB with resolutions ≤2.9 Å to classify the Mg²⁺ inner-sphere binding patterns to nucleotide carbonyl, ribose hydroxyl, cyclic ether, and phosphodiester oxygen atoms. From this classification, we derived a set of "prior-knowledge" nucleobase Mg²⁺ binding sites. We report that crystallographic examples of trustworthy nucleobase Mg²⁺ binding sites are fewer than expected since many of those are associated with misidentified Na⁺ or K⁺ We also emphasize that binding of Na⁺ and K⁺ to nucleic acids is much more frequent than anticipated. Overall, we provide evidence derived from X-ray structures that nucleobases are poor inner-sphere binders for Mg²⁺ but good binders for monovalent ions. Based on strict stereochemical criteria, we propose an extended set of guidelines designed to help in the assignment and validation of ions directly contacting nucleobase and ribose atoms. These guidelines should help in the interpretation of X-ray and cryo-EM solvent density maps. When borderline Mg²⁺ stereochemistry is observed, alternative placement of Na⁺, K⁺, or Ca²⁺ must be considered. We also critically examine the use of lanthanides (Yb³⁺, Tb³⁺) as Mg²⁺ substitutes in crystallography experiments.

High-throughput quantum-mechanics/molecular-mechanics (ONIOM) macromolecular crystallographic refinement with PHENIX/DivCon: the impact of mixed Hamiltonian methods on ligand and protein structure

Conventional macromolecular crystallographic refinement relies on often dubious stereochemical restraints, the preparation of which often requires human validation for unusual species, and on rudimentary energy functionals that are devoid of nonbonding effects owing to electrostatics, polarization, charge transfer or even hydrogen bonding. While this approach has served the crystallographic community for decades, as structure-based drug design/discovery (SBDD) has grown in prominence it has become clear that these conventional methods are less rigorous than they need to be in order to produce properly predictive protein-ligand models, and that the human intervention that is required to successfully treat ligands and other unusual chemistries found in SBDD often precludes high-throughput, automated refinement. Recently, plugins to the Python-based Hierarchical ENvironment for Integrated Xtallography (PHENIX) crystallographic platform have been developed to augment conventional methods with the in situ use of quantum mechanics (QM) applied to ligand(s) along with the surrounding active site(s) at each step of refinement [Borbulevych et al. (2014), Acta Cryst D70, 1233-1247]. This method (Region-QM) significantly increases the accuracy of the X-ray refinement process, and this approach is now used, coupled with experimental density, to accurately determine protonation states, binding modes, ring-flip states, water positions and so on. In the present work, this approach is expanded to include a more rigorous treatment of the entire structure, including the ligand(s), the associated active site(s) and the entire protein, using a fully automated, mixed quantum-mechanics/molecular-mechanics (QM/MM) Hamiltonian recently implemented in the DivCon package. This approach was validated through the automatic treatment of a population of 80 protein-ligand structures chosen from the Astex Diverse Set. Across the entire population, this method results in an average 3.5-fold reduction in ligand strain and a 4.5-fold improvement in MolProbity clashscore, as well as improvements in Ramachandran and rotamer outlier analyses. Overall, these results demonstrate that the use of a structure-wide QM/MM Hamiltonian exhibits improvements in the local structural chemistry of the ligand similar to Region-QM refinement but with significant improvements in the overall structure beyond the active site.

Role of Computational Methods in Going beyond X-ray Crystallography to Explore Protein Structure and Dynamics

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.

Crystallographic snapshots of ligand binding to hexameric purine nucleoside phosphorylase and kinetic studies give insight into the mechanism of catalysis

Purine nucleoside phosphorylase (PNP) catalyses the cleavage of the glycosidic bond of purine nucleosides using phosphate instead of water as a second substrate. PNP from Escherichia coli is a homohexamer, build as a trimer of dimers, and each subunit can be in two conformations, open or closed. This conformational change is induced by the presence of phosphate substrate, and very likely a required step for the catalysis. Closing one active site strongly affects the others, by a yet unclear mechanism and order of events. Kinetic and ligand binding studies show strong negative cooperativity between subunits. Here, for the first time, we managed to monitor the sequence of nucleoside binding to individual subunits in the crystal structures of the wild-type enzyme, showing that first the closed sites, not the open ones, are occupied by the nucleoside. However, two mutations within the active site, Asp204Ala/Arg217Ala, are enough not only to significantly reduce the effectiveness of the enzyme, but also reverse the sequence of the nucleoside binding. In the mutant the open sites, neighbours in a dimer of those in the closed conformation, are occupied as first. This demonstrates how important for the effective catalysis of Escherichia coli PNP is proper subunit cooperation.

A close look onto structural models and primary ligands of metallo-β-lactamases

β-Lactamases are hydrolytic enzymes capable of opening the β-lactam ring of antibiotics such as penicillin, thus endowing the bacteria that produce them with antibiotic resistance. Of particular medical concern are metallo-β-lactamases (MBLs), with an active site built around coordinated Zn cations. MBLs are pan-reactive enzymes that can break down almost all classes of β-lactams, including such last-resort antibiotics as carbapenems. They are not only broad-spectrum-reactive but are often plasmid-borne (e.g., the New Delhi enzyme, NDM), and can spread horizontally even among unrelated bacteria. Acquired MBLs are encoded by mobile genetic elements, which often include other resistance genes, making the microbiological situation particularly alarming. There is an urgent need to develop MBL inhibitors in order to rescue our antibiotic armory. A number of such efforts have been undertaken, most notably using the 3D structures of various MBLs as drug-design targets. Structure-guided drug discovery depends on the quality of the structures that are collected in the Protein Data Bank (PDB) and on the consistency of the information in dedicated β-lactamase databases. We conducted a careful review of the crystal structures of class B β-lactamases, concluding that the quality of these structures varies widely, especially in the regions where small molecules interact with the macromolecules. In a number of examples the interpretation of the bound ligands (e.g., inhibitors, substrate/product analogs) is doubtful or even incorrect, and it appears that in some cases the modeling of ligands was not supported by electron density. For ten MBL structures, alternative interpretations of the original diffraction data could be proposed and the new models have been deposited in the PDB. In four cases, these models, prepared jointly with the authors of the original depositions, superseded the previous deposits. This review emphasizes the importance of critical assessment of structural models describing key drug design targets at the level of the raw experimental data. Since the structures reviewed here are the basis for ongoing design of new MBL inhibitors, it is important to identify and correct the problems with ambiguous crystallographic interpretations, thus enhancing reproducibility in this highly medically relevant area.

Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

Structural biology data archiving - where we are and what lies ahead

For almost 50 years, structural biology has endeavoured to conserve and share its experimental data and their interpretations (usually, atomistic models) through global public archives such as the Protein Data Bank, Electron Microscopy Data Bank and Biological Magnetic Resonance Data Bank (BMRB). These archives are treasure troves of freely accessible data that document our quest for molecular or atomic understanding of biological function and processes in health and disease. They have prepared the field to tackle new archiving challenges as more and more (combinations of) techniques are being utilized to elucidate structure at ever increasing length scales. Furthermore, the field has made substantial efforts to develop validation methods that help users to assess the reliability of structures and to identify the most appropriate data for their needs. In this Review, we present an overview of public data archives in structural biology and discuss the importance of validation for users and producers of structural data. Finally, we sketch our efforts to integrate structural data with bioimaging data and with other sources of biological data. This will make relevant structural information available and more easily discoverable for a wide range of scientists.

Against Method: Table 1-Cui Bono?

The almost universally required "Table 1," summarizing data-collection and data-processing statistics, has in its present form outlived its usefulness in almost all publications of biomolecular crystal structure reports. Information contained in "Table 1" is insufficient to evaluate or repeat the experiment; is redundant with information extractable from deposited diffraction data; and includes data items whose meaning is under increased scrutiny in the crystallographic community. Direct and consistent extraction and analysis of data quality metrics from preferably unmerged intensity data with graphical presentation of reciprocal space features, including impact on map and model features, should replace "Table 1."

Validation of ligands in macromolecular structures determined by X-ray crystallography

Crystallographic studies of ligands bound to biological macromolecules (proteins and nucleic acids) play a crucial role in structure-guided drug discovery and design, and also provide atomic level insights into the physical chemistry of complex formation between macromolecules and ligands. The quality with which small-molecule ligands have been modelled in Protein Data Bank (PDB) entries has been, and continues to be, a matter of concern for many investigators. Correctly interpreting whether electron density found in a binding site is compatible with the soaked or co-crystallized ligand or represents water or buffer molecules is often far from trivial. The Worldwide PDB validation report (VR) provides a mechanism to highlight any major issues concerning the quality of the data and the model at the time of deposition and annotation, so the depositors can fix issues, resulting in improved data quality. The ligand-validation methods used in the generation of the current VRs are described in detail, including an examination of the metrics to assess both geometry and electron-density fit. It is found that the LLDF score currently used to identify ligand electron-density fit outliers can give misleading results and that better ligand-validation metrics are required.

Detect, correct, retract: How to manage incorrect structural models

The massive technical and computational progress of biomolecular crystallography has generated some adverse side effects. Most crystal structure models, produced by crystallographers or well-trained structural biologists, constitute useful sources of information, but occasional extreme outliers remind us that the process of structure determination is not fail-safe. The occurrence of severe errors or gross misinterpretations raises fundamental questions: Why do such aberrations emerge in the first place? How did they evade the sophisticated validation procedures which often produce clear and dire warnings, and why were severe errors not noticed by the depositors themselves, their supervisors, referees and editors? Once detected, what can be done to either correct, improve or eliminate such models? How do incorrect models affect the underlying claims or biomedical hypotheses they were intended, but failed, to support? What is the long-range effect of the propagation of such errors? And finally, what mechanisms can be envisioned to restore the validity of the scientific record and, if necessary, retract publications that are clearly invalidated by the lack of experimental evidence? We suggest that cognitive bias and flawed epistemology are likely at the root of the problem. By using examples from the published literature and from public repositories such as the Protein Data Bank, we provide case summaries to guide correction or improvement of structural models. When strong claims are unsustainable because of a deficient crystallographic model, removal of such a model and even retraction of the affected publication are necessary to restore the integrity of the scientific record.

The science is in the data

Understanding published research results should be through one's own eyes and include the opportunity to work with raw diffraction data to check the various decisions made in the analyses by the original authors. Today, preserving raw diffraction data is technically and organizationally viable at a growing number of data archives, both centralized and distributed, which are empowered to register data sets and obtain a preservation descriptor, typically a 'digital object identifier'. This introduces an important role of preserving raw data, namely understanding where we fail in or could improve our analyses. Individual science area case studies in crystallography are provided.

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein-ligand crystal structures? This value is important biophysically, for better understanding the protein-ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh-resolution crystal structures that can be refined without restraints.

X-ray dynamical diffraction in amino acid crystals: a step towards improving structural resolution of biological molecules via physical phase measurements

X-ray phase measurements have been applied to study hydrogen bonds and radiation damage in amino acid crystals.

In this work, experimental and data analysis procedures were developed and applied for studying amino acid crystals by means of X-ray phase measurements. The results clearly demonstrated the sensitivity of invariant triplet phases to electronic charge distribution in d-alanine crystals, providing useful information for molecular dynamics studies of intermolecular forces. The feasibility of using phase measurements to investigate radiation damage mechanisms is also discussed on experimental and theoretical grounds.

Resolving the cofactor-binding site in the proline biosynthetic enzyme human pyrroline-5-carboxylate reductase 1

Pyrroline-5-carboxylate reductase (PYCR) is the final enzyme in proline biosynthesis, catalyzing the NAD(P)H-dependent reduction of Δ1-pyrroline-5-carboxylate (P5C) to proline. Mutations in the PYCR1 gene alter mitochondrial function and cause the connective tissue disorder cutis laxa. Furthermore, PYCR1 is overexpressed in multiple cancers, and the PYCR1 knock-out suppresses tumorigenic growth, suggesting that PYCR1 is a potential cancer target. However, inhibitor development has been stymied by limited mechanistic details for the enzyme, particularly in light of a previous crystallographic study that placed the cofactor-binding site in the C-terminal domain rather than the anticipated Rossmann fold of the N-terminal domain. To fill this gap, we report crystallographic, sedimentation-velocity, and kinetics data for human PYCR1. Structures of binary complexes of PYCR1 with NADPH or proline determined at 1.9 Å resolution provide insight into cofactor and substrate recognition. We see NADPH bound to the Rossmann fold, over 25 Å from the previously proposed site. The 1.85 Å resolution structure of a ternary complex containing NADPH and a P5C/proline analog provides a model of the Michaelis complex formed during hydride transfer. Sedimentation velocity shows that PYCR1 forms a concentration-dependent decamer in solution, consistent with the pentamer-of-dimers assembly seen crystallographically. Kinetic and mutational analysis confirmed several features seen in the crystal structure, including the importance of a hydrogen bond between Thr-238 and the substrate as well as limited cofactor discrimination.

A multi-crystal method for extracting obscured crystallographic states from conventionally uninterpretable electron density

In macromolecular crystallography, the rigorous detection of changed states (for example, ligand binding) is difficult unless signal is strong. Ambiguous (‘weak' or ‘noisy') density is experimentally common, since molecular states are generally only fractionally present in the crystal. Existing methodologies focus on generating maximally accurate maps whereby minor states become discernible; in practice, such map interpretation is disappointingly subjective, time-consuming and methodologically unsound. Here we report the PanDDA method, which automatically reveals clear electron density for the changed state—even from inaccurate maps—by subtracting a proportion of the confounding ‘ground state'; changed states are objectively identified from statistical analysis of density distributions. The method is completely general, implying new best practice for all changed-state studies, including the routine collection of multiple ground-state crystals. More generally, these results demonstrate: the incompleteness of atomic models; that single data sets contain insufficient information to model them fully; and that accuracy requires further map-deconvolution approaches.

Building a ligand into a weak region of an electron density map of a protein is a subjective process. Here, the authors present a new method to obtain a clear electron density for a bound ligand based on multi-crystal experiments and 3D background correction.

The XChemExplorer graphical workflow tool for routine or large-scale protein-ligand structure determination

XChemExplorer (XCE) is a data-management and workflow tool to support large-scale simultaneous analysis of protein-ligand complexes during structure-based ligand discovery (SBLD). The user interfaces of established crystallographic software packages such as CCP4 [Winn et al. (2011), Acta Cryst. D67, 235-242] or PHENIX [Adams et al. (2010), Acta Cryst. D66, 213-221] have entrenched the paradigm that a `project' is concerned with solving one structure. This does not hold for SBLD, where many almost identical structures need to be solved and analysed quickly in one batch of work. Functionality to track progress and annotate structures is essential. XCE provides an intuitive graphical user interface which guides the user from data processing, initial map calculation, ligand identification and refinement up until data dissemination. It provides multiple entry points depending on the need of each project, enables batch processing of multiple data sets and records metadata, progress and annotations in an SQLite database. XCE is freely available and works on any Linux and Mac OS X system, and the only dependency is to have the latest version of CCP4 installed. The design and usage of this tool are described here, and its usefulness is demonstrated in the context of fragment-screening campaigns at the Diamond Light Source. It is routinely used to analyse projects comprising 1000 data sets or more, and therefore scales well to even very large ligand-design projects.

Estimation of the protein-ligand interaction energy for model building and validation

Macromolecular X-ray crystallography is one of the main experimental techniques to visualize protein-ligand interactions. The high complexity of the ligand universe, however, has delayed the development of efficient methods for the automated identification, fitting and validation of ligands in their electron-density clusters. The identification and fitting are primarily based on the density itself and do not take into account the protein environment, which is a step that is only taken during the validation of the proposed binding mode. Here, a new approach, based on the estimation of the major energetic terms of protein-ligand interaction, is introduced for the automated identification of crystallographic ligands in the indicated binding site with ARP/wARP. The applicability of the method to the validation of protein-ligand models from the Protein Data Bank is demonstrated by the detection of models that are `questionable' and the pinpointing of unfavourable interatomic contacts.

Twilight reloaded: the peptide experience

The de facto commoditization of biomolecular crystallography as a result of almost disruptive instrumentation automation and continuing improvement of software allows any sensibly trained structural biologist to conduct crystallographic studies of biomolecules with reasonably valid outcomes: that is, models based on properly interpreted electron density. Robust validation has led to major mistakes in the protein part of structure models becoming rare, but some depositions of protein-peptide complex structure models, which generally carry significant interest to the scientific community, still contain erroneous models of the bound peptide ligand. Here, the protein small-molecule ligand validation tool Twilight is updated to include peptide ligands. (i) The primary technical reasons and potential human factors leading to problems in ligand structure models are presented; (ii) a new method used to score peptide-ligand models is presented; (iii) a few instructive and specific examples, including an electron-density-based analysis of peptide-ligand structures that do not contain any ligands, are discussed in detail; (iv) means to avoid such mistakes and the implications for database integrity are discussed and (v) some suggestions as to how journal editors could help to expunge errors from the Protein Data Bank are provided.

Ligand fitting with CCP4

Crystal structures of protein-ligand complexes are often used to infer biology and inform structure-based drug discovery. Hence, it is important to build accurate, reliable models of ligands that give confidence in the interpretation of the respective protein-ligand complex. This paper discusses key stages in the ligand-fitting process, including ligand binding-site identification, ligand description and conformer generation, ligand fitting, refinement and subsequent validation. The CCP4 suite contains a number of software tools that facilitate this task: AceDRG for the creation of ligand descriptions and conformers, Lidia and JLigand for two-dimensional and three-dimensional ligand editing and visual analysis, Coot for density interpretation, ligand fitting, analysis and validation, and REFMAC5 for macromolecular refinement. In addition to recent advancements in automatic carbohydrate building in Coot (LO/Carb) and ligand-validation tools (FLEV), the release of the CCP4i2 GUI provides an integrated solution that streamlines the ligand-fitting workflow, seamlessly passing results from one program to the next. The ligand-fitting process is illustrated using instructive practical examples, including problematic cases such as post-translational modifications, highlighting the need for careful analysis and rigorous validation.

AceDRG: a stereochemical description generator for ligands

The program AceDRG is designed for the derivation of stereochemical information about small molecules. It uses local chemical and topological environment-based atom typing to derive and organize bond lengths and angles from a small-molecule database: the Crystallography Open Database (COD). Information about the hybridization states of atoms, whether they belong to small rings (up to seven-membered rings), ring aromaticity and nearest-neighbour information is encoded in the atom types. All atoms from the COD have been classified according to the generated atom types. All bonds and angles have also been classified according to the atom types and, in a certain sense, bond types. Derived data are tabulated in a machine-readable form that is freely available from CCP4. AceDRG can also generate stereochemical information, provided that the basic bonding pattern of a ligand is known. The basic bonding pattern is perceived from one of the computational chemistry file formats, including SMILES, mmCIF, SDF MOL and SYBYL MOL2 files. Using the bonding chemistry, atom types, and bond and angle tables generated from the COD, AceDRG derives the `ideal' bond lengths, angles, plane groups, aromatic rings and chirality information, and writes them to an mmCIF file that can be used by the refinement program REFMAC5 and the model-building program Coot. Other refinement and model-building programs such as PHENIX and BUSTER can also use these files. AceDRG also generates one or more coordinate sets corresponding to the most favourable conformation(s) of a given ligand. AceDRG employs RDKit for chemistry perception and for initial conformation generation, as well as for the interpretation of SMILES strings, SDF MOL and SYBYL MOL2 files.

Validation and extraction of molecular-geometry information from small-molecule databases

A freely available small-molecule structure database, the Crystallography Open Database (COD), is used for the extraction of molecular-geometry information on small-molecule compounds. The results are used for the generation of new ligand descriptions, which are subsequently used by macromolecular model-building and structure-refinement software. To increase the reliability of the derived data, and therefore the new ligand descriptions, the entries from this database were subjected to very strict validation. The selection criteria made sure that the crystal structures used to derive atom types, bond and angle classes are of sufficiently high quality. Any suspicious entries at a crystal or molecular level were removed from further consideration. The selection criteria included (i) the resolution of the data used for refinement (entries solved at 0.84 Å resolution or higher) and (ii) the structure-solution method (structures must be from a single-crystal experiment and all atoms of generated molecules must have full occupancies), as well as basic sanity checks such as (iii) consistency between the valences and the number of connections between atoms, (iv) acceptable bond-length deviations from the expected values and (v) detection of atomic collisions. The derived atom types and bond classes were then validated using high-order moment-based statistical techniques. The results of the statistical analyses were fed back to fine-tune the atom typing. The developed procedure was repeated four times, resulting in fine-grained atom typing, bond and angle classes. The procedure will be repeated in the future as and when new entries are deposited in the COD. The whole procedure can also be applied to any source of small-molecule structures, including the Cambridge Structural Database and the ZINC database.

Mg2+ ions: do they bind to nucleobase nitrogens?

Given the many roles proposed for Mg2+ in nucleic acids, it is essential to accurately determine their binding modes. Here, we surveyed the PDB to classify Mg2+ inner-sphere binding patterns to nucleobase imine N1/N3/N7 atoms. Among those, purine N7 atoms are considered to be the best nucleobase binding sites for divalent metals. Further, Mg2+ coordination to N7 has been implied in several ribozyme catalytic mechanisms. We report that Mg2+ assigned near imine nitrogens derive mostly from poor interpretations of electron density patterns and are most often misidentified Na+, K+, NH4+ ions, water molecules or spurious density peaks. Consequently, apart from few documented exceptions, Mg2+ ions do not bind to N7 atoms. Without much of a surprise, Mn2+, Zn2+ and Cd2+, which have a higher affinity for nitrogens, may contact N7 atoms when present in crystallization buffers. In this respect, we describe for the first time a potential Zn2+ ribosomal binding site involving two purine N7 atoms. Further, we provide a set of guidelines to help in the assignment of Mg2+ in crystallographic, cryo-EM, NMR and model building practices and discuss implications of our findings related to ion substitution experiments.

Validation of Protein-Ligand Crystal Structure Models: Small Molecule and Peptide Ligands

Models of target proteins in complex with small molecule ligands or peptide ligands are of significant interest to the biomedical research community. Structure-guided lead discovery and structure-based drug design make extensive use of such models. The bound ligands comprise only a small fraction of the total X-ray scattering mass, and therefore particular care must be taken to properly validate the atomic model of the ligand as experimental data can often be scarce. The ligand model must be validated against both the primary experimental data and the local environment, specifically: (1) the primary evidence in the form of the electron density, (2) examined for reasonable stereochemistry, and (3) the chemical plausibility of the binding interactions must be inspected. Tools that assist the researcher in the validation process are presented.

Outcome of the First wwPDB/CCDC/D3R Ligand Validation Workshop

Crystallographic studies of ligands bound to biological macromolecules (proteins and nucleic acids) represent an important source of information concerning drug-target interactions, providing atomic level insights into the physical chemistry of complex formation between macromolecules and ligands. Of the more than 115,000 entries extant in the Protein Data Bank (PDB) archive, ∼75% include at least one non-polymeric ligand. Ligand geometrical and stereochemical quality, the suitability of ligand models for in silico drug discovery and design, and the goodness-of-fit of ligand models to electron-density maps vary widely across the archive. We describe the proceedings and conclusions from the first Worldwide PDB/Cambridge Crystallographic Data Center/Drug Design Data Resource (wwPDB/CCDC/D3R) Ligand Validation Workshop held at the Research Collaboratory for Structural Bioinformatics at Rutgers University on July 30-31, 2015. Experts in protein crystallography from academe and industry came together with non-profit and for-profit software providers for crystallography and with experts in computational chemistry and data archiving to discuss and make recommendations on best practices, as framed by a series of questions central to structural studies of macromolecule-ligand complexes. What data concerning bound ligands should be archived in the PDB? How should the ligands be best represented? How should structural models of macromolecule-ligand complexes be validated? What supplementary information should accompany publications of structural studies of biological macromolecules? Consensus recommendations on best practices developed in response to each of these questions are provided, together with some details regarding implementation. Important issues addressed but not resolved at the workshop are also enumerated.

The Intersection of Structural and Chemical Biology - An Essential Synergy

The continual improvement in our ability to generate high resolution structural models of biological molecules has stimulated and supported innovative chemical biology projects that target increasingly challenging ligand interaction sites. In this review we outline some of the recent developments in chemical biology and rational ligand design and show selected examples that illustrate the synergy between these research areas.

Improved ligand geometries in crystallographic refinement using AFITT in PHENIX

Modern crystal structure refinement programs rely on geometry restraints to overcome the challenge of a low data-to-parameter ratio. While the classical Engh and Huber restraints work well for standard amino-acid residues, the chemical complexity of small-molecule ligands presents a particular challenge. Most current approaches either limit ligand restraints to those that can be readily described in the Crystallographic Information File (CIF) format, thus sacrificing chemical flexibility and energetic accuracy, or they employ protocols that substantially lengthen the refinement time, potentially hindering rapid automated refinement workflows. PHENIX-AFITT refinement uses a full molecular-mechanics force field for user-selected small-molecule ligands during refinement, eliminating the potentially difficult problem of finding or generating high-quality geometry restraints. It is fully integrated with a standard refinement protocol and requires practically no additional steps from the user, making it ideal for high-throughput workflows. PHENIX-AFITT refinements also handle multiple ligands in a single model, alternate conformations and covalently bound ligands. Here, the results of combining AFITT and the PHENIX software suite on a data set of 189 protein-ligand PDB structures are presented. Refinements using PHENIX-AFITT significantly reduce ligand conformational energy and lead to improved geometries without detriment to the fit to the experimental data. For the data presented, PHENIX-AFITT refinements result in more chemically accurate models for small-molecule ligands.

New Biological Insights from Better Structure Models

Structure validation is a key component of all steps in the structure determination process, from structure building, refinement, deposition, and evaluation all the way to post-deposition optimisation of structures in the Protein Data Bank (PDB) by re-refinement and re-building. Today, many aspects of protein structures are understood better than 10years ago, and combined with improved software and more computing power, the automated PDB_REDO procedure can significantly improve about 85% of all X-ray structures ever deposited in the PDB. We review structure validation, structure improvement, and a series of validation resources and facilities that give access to improved PDB files and to reports on the quality of the original and the improved structures. Post-deposition optimisation generally leads to improved protein structures and a series of examples will illustrate how that, in turn, leads to improved or even novel biological insights.

Fitmunk: improving protein structures by accurate, automatic modeling of side-chain conformations

Improvements in crystallographic hardware and software have allowed automated structure-solution pipelines to approach a near-`one-click' experience for the initial determination of macromolecular structures. However, in many cases the resulting initial model requires a laborious, iterative process of refinement and validation. A new method has been developed for the automatic modeling of side-chain conformations that takes advantage of rotamer-prediction methods in a crystallographic context. The algorithm, which is based on deterministic dead-end elimination (DEE) theory, uses new dense conformer libraries and a hybrid energy function derived from experimental data and prior information about rotamer frequencies to find the optimal conformation of each side chain. In contrast to existing methods, which incorporate the electron-density term into protein-modeling frameworks, the proposed algorithm is designed to take advantage of the highly discriminatory nature of electron-density maps. This method has been implemented in the program Fitmunk, which uses extensive conformational sampling. This improves the accuracy of the modeling and makes it a versatile tool for crystallographic model building, refinement and validation. Fitmunk was extensively tested on over 115 new structures, as well as a subset of 1100 structures from the PDB. It is demonstrated that the ability of Fitmunk to model more than 95% of side chains accurately is beneficial for improving the quality of crystallographic protein models, especially at medium and low resolutions. Fitmunk can be used for model validation of existing structures and as a tool to assess whether side chains are modeled optimally or could be better fitted into electron density. Fitmunk is available as a web service at http://kniahini.med.virginia.edu/fitmunk/server/ or at http://fitmunk.bitbucket.org/.

Data Mining of Macromolecular Structures

The use of macromolecular structures is widespread for a variety of applications, from teaching protein structure principles all the way to ligand optimization in drug development. Applying data mining techniques on these experimentally determined structures requires a highly uniform, standardized structural data source. The Protein Data Bank (PDB) has evolved over the years toward becoming the standard resource for macromolecular structures. However, the process selecting the data most suitable for specific applications is still very much based on personal preferences and understanding of the experimental techniques used to obtain these models. In this chapter, we will first explain the challenges with data standardization, annotation, and uniformity in the PDB entries determined by X-ray crystallography. We then discuss the specific effect that crystallographic data quality and model optimization methods have on structural models and how validation tools can be used to make informed choices. We also discuss specific advantages of using the PDB_REDO databank as a resource for structural data. Finally, we will provide guidelines on how to select the most suitable protein structure models for detailed analysis and how to select a set of structure models suitable for data mining.

Atomic resolution experimental phase information reveals extensive disorder and bound 2-methyl-2,4-pentanediol in Ca(2+)-calmodulin

Calmodulin (CaM) is the primary calcium signaling protein in eukaryotes and has been extensively studied using various biophysical techniques. Prior crystal structures have noted the presence of ambiguous electron density in both hydrophobic binding pockets of Ca(2+)-CaM, but no assignment of these features has been made. In addition, Ca(2+)-CaM samples many conformational substates in the crystal and accurately modeling the full range of this functionally important disorder is challenging. In order to characterize these features in a minimally biased manner, a 1.0 Å resolution single-wavelength anomalous diffraction data set was measured for selenomethionine-substituted Ca(2+)-CaM. Density-modified electron-density maps enabled the accurate assignment of Ca(2+)-CaM main-chain and side-chain disorder. These experimental maps also substantiate complex disorder models that were automatically built using low-contour features of model-phased electron density. Furthermore, experimental electron-density maps reveal that 2-methyl-2,4-pentanediol (MPD) is present in the C-terminal domain, mediates a lattice contact between N-terminal domains and may occupy the N-terminal binding pocket. The majority of the crystal structures of target-free Ca(2+)-CaM have been derived from crystals grown using MPD as a precipitant, and thus MPD is likely to be bound in functionally critical regions of Ca(2+)-CaM in most of these structures. The adventitious binding of MPD helps to explain differences between the Ca(2+)-CaM crystal and solution structures and is likely to favor more open conformations of the EF-hands in the crystal.

Progress in protein crystallography

Macromolecular crystallography evolved enormously from the pioneering days, when structures were solved by "wizards" performing all complicated procedures almost by hand. In the current situation crystal structures of large systems can be often solved very effectively by various powerful automatic programs in days or hours, or even minutes. Such progress is to a large extent coupled to the advances in many other fields, such as genetic engineering, computer technology, availability of synchrotron beam lines and many other techniques, creating the highly interdisciplinary science of macromolecular crystallography. Due to this unprecedented success crystallography is often treated as one of the analytical methods and practiced by researchers interested in structures of macromolecules, but not highly competent in the procedures involved in the process of structure determination. One should therefore take into account that the contemporary, highly automatic systems can produce results almost without human intervention, but the resulting structures must be carefully checked and validated before their release into the public domain.

Models of protein-ligand crystal structures: trust, but verify

X-ray crystallography provides the most accurate models of protein-ligand structures. These models serve as the foundation of many computational methods including structure prediction, molecular modelling, and structure-based drug design. The success of these computational methods ultimately depends on the quality of the underlying protein-ligand models. X-ray crystallography offers the unparalleled advantage of a clear mathematical formalism relating the experimental data to the protein-ligand model. In the case of X-ray crystallography, the primary experimental evidence is the electron density of the molecules forming the crystal. The first step in the generation of an accurate and precise crystallographic model is the interpretation of the electron density of the crystal, typically carried out by construction of an atomic model. The atomic model must then be validated for fit to the experimental electron density and also for agreement with prior expectations of stereochemistry. Stringent validation of protein-ligand models has become possible as a result of the mandatory deposition of primary diffraction data, and many computational tools are now available to aid in the validation process. Validation of protein-ligand complexes has revealed some instances of overenthusiastic interpretation of ligand density. Fundamental concepts and metrics of protein-ligand quality validation are discussed and we highlight software tools to assist in this process. It is essential that end users select high quality protein-ligand models for their computational and biological studies, and we provide an overview of how this can be achieved.

X-ray crystallography over the past decade for novel drug discovery - where are we heading next?

INTRODUCTION

Macromolecular X-ray crystallography has been the primary methodology for determining the three-dimensional structures of proteins, nucleic acids and viruses. Structural information has paved the way for structure-guided drug discovery and laid the foundations for structural bioinformatics. However, X-ray crystallography still has a few fundamental limitations, some of which may be overcome and complemented using emerging methods and technologies in other areas of structural biology.

AREAS COVERED

This review describes how structural knowledge gained from X-ray crystallography has been used to advance other biophysical methods for structure determination (and vice versa). This article also covers current practices for integrating data generated by other biochemical and biophysical methods with those obtained from X-ray crystallography. Finally, the authors articulate their vision about how a combination of structural and biochemical/biophysical methods may improve our understanding of biological processes and interactions.

EXPERT OPINION

X-ray crystallography has been, and will continue to serve as, the central source of experimental structural biology data used in the discovery of new drugs. However, other structural biology techniques are useful not only to overcome the major limitation of X-ray crystallography, but also to provide complementary structural data that is useful in drug discovery. The use of recent advancements in biochemical, spectroscopy and bioinformatics methods may revolutionize drug discovery, albeit only when these data are combined and analyzed with effective data management systems. Accurate and complete data management is crucial for developing experimental procedures that are robust and reproducible.