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

A database overview of metal-coordination distances in metalloproteins

Metalloproteins are ubiquitous in all living organisms and take part in a very wide range of biological processes. For this reason, their experimental characterization is crucial to obtain improved knowledge of their structure and biological functions. The three-dimensional structure represents highly relevant information since it provides insight into the interaction between the metal ion(s) and the protein fold. Such interactions determine the chemical reactivity of the bound metal. The available PDB structures can contain errors due to experimental factors such as poor resolution and radiation damage. A lack of use of distance restraints during the refinement and validation process also impacts the structure quality. Here, the aim was to obtain a thorough overview of the distribution of the distances between metal ions and their donor atoms through the statistical analysis of a data set based on more than 115 000 metal-binding sites in proteins. This analysis not only produced reference data that can be used by experimentalists to support the structure-determination process, for example as refinement restraints, but also resulted in an improved insight into how protein coordination occurs for different metals and the nature of their binding interactions. In particular, the features of carboxylate coordination were inspected, which is the only type of interaction that is commonly present for nearly all metals.

Protein purification, crystallization, and structure determination of transcription factor YhaJ in complex with DNT metabolites

The microbial transcription factor YhaJ responds to 2,4-dinitrotoluene (DNT) derivatives. Here, we describe steps for overexpression and purification of the protein, characterization for the binding of a DNT derivative methylhydroquinone, and crystallization by using a random seeding technique. We then detail procedures for structure determination by employing the crystal-twin resolving processes. This protocol can also be performed using other DNT derivatives. For complete details on the use and execution of this protocol, please refer to Kim et al.1.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Automated benchmarking of combined protein structure and ligand conformation prediction

The prediction of protein-ligand complexes (PLC), using both experimental and predicted structures, is an active and important area of research, underscored by the inclusion of the Protein-Ligand Interaction category in the latest round of the Critical Assessment of Protein Structure Prediction experiment CASP15. The prediction task in CASP15 consisted of predicting both the three-dimensional structure of the receptor protein as well as the position and conformation of the ligand. This paper addresses the challenges and proposed solutions for devising automated benchmarking techniques for PLC prediction. The reliability of experimentally solved PLC as ground truth reference structures is assessed using various validation criteria. Similarity of PLC to previously released complexes are employed to judge PLC diversity and the difficulty of a PLC as a prediction target. We show that the commonly used PDBBind time-split test-set is inappropriate for comprehensive PLC evaluation, with state-of-the-art tools showing conflicting results on a more representative and high quality dataset constructed for benchmarking purposes. We also show that redocking on crystal structures is a much simpler task than docking into predicted protein models, demonstrated by the two PLC-prediction-specific scoring metrics created. Finally, we introduce a fully automated pipeline that predicts PLC and evaluates the accuracy of the protein structure, ligand pose, and protein-ligand interactions.

Three structures of PSI-LHCI from Chlamydomonas reinhardtii suggest a resting state re-activated by ferredoxin

Photosystem I (PSI) from the green alga Chlamydomonas reinhardtii, with various numbers of membrane bound antenna complexes (LHCI), has been described in great detail. In contrast, structural characterization of soluble binding partners is less advanced. Here, we used X-ray crystallography and single particle cryo-EM to investigate three structures of the PSI-LHCI supercomplex from Chlamydomonas reinhardtii. An X-ray structure demonstrates the absence of six chlorophylls from the luminal side of the LHCI belts, suggesting these pigments were either physically absent or less stably associated with the complex, potentially influencing excitation transfer significantly. CryoEM revealed extra densities on luminal and stromal sides of the supercomplex, situated in the vicinity of the electron transfer sites. These densities disappeared after the binding of oxidized ferredoxin to PSI-LHCI. Based on these structures, we propose the existence of a PSI-LHCI resting state with a reduced active chlorophyll content, electron donors docked in waiting positions and regulatory binding partners positioned at the electron acceptor site. The resting state PSI-LHCI supercomplex would be recruited to its active form by the availability of oxidized ferredoxin.

Development of Novel Small-Molecule Activators of Pyruvate Kinase Muscle Isozyme 2, PKM2, to Reduce Photoreceptor Apoptosis

Treatment options are lacking to prevent photoreceptor death and subsequent vision loss. Previously, we demonstrated that reprogramming metabolism via the pharmacologic activation of PKM2 is a novel photoreceptor neuroprotective strategy. However, the features of the tool compound used in those studies, ML-265, preclude its advancement as an intraocular, clinical candidate. This study sought to develop the next generation of small-molecule PKM2 activators, aimed specifically for delivery into the eye. Compounds were developed that replaced the thienopyrrolopyridazinone core of ML-265 and modified the aniline and methyl sulfoxide functional groups. Compound 2 demonstrated that structural changes to the ML-265 scaffold are tolerated from a potency and efficacy standpoint, allow for a similar binding mode to the target, and circumvent apoptosis in models of outer retinal stress. To overcome the low solubility and problematic functional groups of ML-265, compound 2's efficacious and versatile core structure for the incorporation of diverse functional groups was then utilized to develop novel PKM2 activators with improved solubility, lack of structural alerts, and retained potency. No other molecules are in the pharmaceutical pipeline for the metabolic reprogramming of photoreceptors. Thus, this study is the first to cultivate the next generation of novel, structurally diverse, small-molecule PKM2 activators for delivery into the eye.

Interactions of hydrolyzed β-lactams with the L1 metallo-β-lactamase: Crystallography supports stereoselective binding of cephem/carbapenem products

L1 is a dizinc subclass B3 metallo-β-lactamase (MBL) that hydrolyzes most β-lactam antibiotics and is a key resistance determinant in the Gram-negative pathogen Stenotrophomonas maltophilia, an important cause of nosocomial infections in immunocompromised patients. L1 is not usefully inhibited by MBL inhibitors in clinical trials, underlying the need for further studies on L1 structure and mechanism. We describe kinetic studies and crystal structures of L1 in complex with hydrolyzed β-lactams from the penam (mecillinam), cephem (cefoxitin/cefmetazole) and carbapenem (tebipenem, doripenem and panipenem) classes. Despite differences in their structures, all the β-lactam-derived products hydrogen bond to Tyr33, Ser221 and Ser225 and are stabilized by interactions with a conserved hydrophobic pocket. The carbapenem products were modelled as Δ¹-imines, with (2S)-stereochemistry. Their binding mode is determined by the presence of a 1β-methyl substituent: the Zn-bridging hydroxide either interacts with the C-6 hydroxyethyl group (1β-hydrogen-containing carbapenems), or is displaced by the C-6 carboxylate (1β-methyl-containing carbapenems). Unexpectedly, the mecillinam product is a rearranged N-formyl amide rather than penicilloic acid, with the N-formyl oxygen interacting with the Zn-bridging hydroxide. NMR studies imply mecillinam rearrangement can occur non-enzymatically in solution. Cephem-derived imine products are bound with (3R)-stereochemistry and retain their 3' leaving groups, likely representing stable endpoints, rather than intermediates, in MBL-catalyzed hydrolysis. Our structures show preferential complex formation by carbapenem- and cephem-derived species protonated on the equivalent (β) faces, and so identify interactions that stabilize diverse hydrolyzed antibiotics. These results may be exploited in developing antibiotics, and β-lactamase inhibitors, that form long-lasting complexes with dizinc MBLs.

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.

Structure and mechanism of the bacterial lipid ABC transporter, MlaFEDB

The cell envelope of Gram-negative bacteria is composed of an inner membrane, outer membane, and an intervening periplasmic space. How the outer membrane lipids are trafficked and assembled there, and how the asymmetry of the outer membrane is maintained is an area of intense research. The Mla system has been implicated in the maintenance of lipid asymmetry in the outer membrane, and is generally thought to drive the removal of mislocalized phospholipids from the outer membrane and their retrograde transport to the inner membrane. At the heart of the Mla pathway is a structurally unique ABC transporter complex in the inner membrane, called MlaFEDB. Recently, an explosion of cryo-EM studies has begun to shed light on the structure and lipid translocation mechanism of MlaFEDB, with many parallels to other ABC transporter families, including human ABCA and ABCG, as well as bacterial lipopolysaccharide and O-antigen transporters. Here we synthesize information from all available structures, and propose a model for lipid trafficking across the cell envelope by MlaFEDB.

Mechanistic Insights into the Functioning of a Two-Subunit GMP Synthetase, an Allosterically Regulated, Ammonia Channeling Enzyme

Guanosine 5'-monophosphate (GMP) synthetases, enzymes that catalyze the conversion of xanthosine 5'-monophosphate (XMP) to GMP, are composed of two different catalytic units, which are either two domains of a polypeptide chain or two subunits that associate to form a complex. The glutamine amidotransferase (GATase) unit hydrolyzes glutamine generating ammonia, and the ATP pyrophosphatase (ATPPase) unit catalyzes the formation of an AMP-XMP intermediate. The substrate-bound ATPPase allosterically activates GATase, and the ammonia thus generated is tunneled to the ATPPase active site where it reacts with AMP-XMP generating GMP. In ammonia channeling enzymes reported thus far, a tight complex of the two subunits is observed, while the interaction of the two subunits of Methanocaldococcus jannaschii GMP synthetase (MjGMPS) is transient with the underlying mechanism of allostery and substrate channeling largely unclear. Here, we present a mechanistic model encompassing the various steps in the catalytic cycle of MjGMPS based on biochemical experiments, crystal structure, and cross-linking mass spectrometry guided integrative modeling. pH dependence of enzyme kinetics establishes that ammonia is tunneled across the subunits with the lifetime of the complex being ≤0.5 s. The crystal structure of the XMP-bound ATPPase subunit reported herein highlights the role of conformationally dynamic loops in enabling catalysis. The structure of MjGMPS derived using restraints obtained from cross-linking mass spectrometry has enabled the visualization of subunit interactions that enable allostery under catalytic conditions. We integrate the results and propose a functional mechanism for MjGMPS detailing the various steps involved in catalysis.

Structures of Multisubunit Membrane Complexes With the CRYO ARM 200

Progress in structural membrane biology has been significantly accelerated by the ongoing ‘Resolution Revolution’ in cryo-electron microscopy (cryo-EM). In particular, structure determination by single-particle analysis has evolved into the most powerful method for atomic model building of multisubunit membrane protein complexes. This has created an ever-increasing demand in cryo-EM machine time, which to satisfy is in need of new and affordable cryo-electron microscopes. Here, we review our experience in using the JEOL CRYO ARM 200 prototype for the structure determination by single-particle analysis of three different multisubunit membrane complexes: the Thermus thermophilus V-type ATPase V_O complex, the Thermosynechococcus elongatus photosystem I monomer and the flagellar motor lipopolysaccharide peptidoglycan ring (LP ring) from Salmonella enterica.

Structural Insight into Effective Inhibitors' Binding to Toxoplasma gondii Dihydrofolate Reductase Thymidylate Synthase

Pyrimethamine (Pyr), a known dihydrofolate reductase (DHFR) inhibitor, has long been used to treat toxoplasmosis caused by Toxoplasma gondii (Tg) infection. However, Pyr is effective only at high doses with associated toxicity to patients, calling for safer alternative treatments. In this study, we investigated a series of Pyr analogues, previously developed as DHFR inhibitors of Plasmodium falciparum bifunctional DHFR-thymidylate synthase (PfDHFR-TS), for their activity against T. gondii DHFR-TS (TgDHFR-TS). Of these, a set of compounds with a substitution at the C⁶ position of the pyrimidine ring exhibited high binding affinities (in a low nanomolar range) against TgDHFR-TS and in vitro T. gondii inhibitory activity. Three-dimensional structures of TgDHFR-TS reported here include the ternary complexes with Pyr, P39, or P40. A comparison of these structures showed the minor steric strain between the p-chlorophenyl group of Pyr and Thr83 of TgDHFR-TS. Such a conflict was relieved in the complexes with the two analogues, P39 and P40, explaining their highest binding affinities described herein. Moreover, these structures suggested that the hydrophobic environment in the active-site pocket could be used for drug design to increase the potency and selectivity of antifolate inhibitors. These findings would accelerate the development of new antifolate drugs to treat toxoplasmosis.

Simplified quality assessment for small-molecule ligands in the Protein Data Bank

More than 70% of the experimentally determined macromolecular structures in the Protein Data Bank (PDB) contain small-molecule ligands. Quality indicators of ∼643,000 ligands present in ∼106,000 PDB X-ray crystal structures have been analyzed. Ligand quality varies greatly with regard to goodness of fit between ligand structure and experimental data, deviations in bond lengths and angles from known chemical structures, and inappropriate interatomic clashes between the ligand and its surroundings. Based on principal component analysis, correlated quality indicators of ligand structure have been aggregated into two largely orthogonal composite indicators measuring goodness of fit to experimental data and deviation from ideal chemical structure. Ranking of the composite quality indicators across the PDB archive enabled construction of uniformly distributed composite ranking score. This score is implemented at RCSB.org to compare chemically identical ligands in distinct PDB structures with easy-to-interpret two-dimensional ligand quality plots, allowing PDB users to quickly assess ligand structure quality and select the best exemplars.

Modernized uniform representation of carbohydrate molecules in the Protein Data Bank

Since 1971, the Protein Data Bank (PDB) has served as the single global archive for experimentally determined 3D structures of biological macromolecules made freely available to the global community according to the FAIR principles of Findability-Accessibility-Interoperability-Reusability. During the first 50 years of continuous PDB operations, standards for data representation have evolved to better represent rich and complex biological phenomena. Carbohydrate molecules present in more than 14,000 PDB structures have recently been reviewed and remediated to conform to a new standardized format. This machine-readable data representation for carbohydrates occurring in the PDB structures and the corresponding reference data improves the findability, accessibility, interoperability and reusability of structural information pertaining to these molecules. The PDB Exchange MacroMolecular Crystallographic Information File data dictionary now supports (i) standardized atom nomenclature that conforms to International Union of Pure and Applied Chemistry-International Union of Biochemistry and Molecular Biology (IUPAC-IUBMB) recommendations for carbohydrates, (ii) uniform representation of branched entities for oligosaccharides, (iii) commonly used linear descriptors of carbohydrates developed by the glycoscience community and (iv) annotation of glycosylation sites in proteins. For the first time, carbohydrates in PDB structures are consistently represented as collections of standardized monosaccharides, which precisely describe oligosaccharide structures and enable improved carbohydrate visualization, structure validation, robust quantitative and qualitative analyses, search for dendritic structures and classification. The uniform representation of carbohydrate molecules in the PDB described herein will facilitate broader usage of the resource by the glycoscience community and researchers studying glycoproteins.

'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.

Predictable cholesterol binding sites in GPCRs lack consensus motifs

A rich diversity of transmembrane G protein-coupled receptors (GPCRs) are used by eukaryotes to sense physical and chemical signals. In humans alone, 800 GPCRs comprise the largest and most therapeutically targeted receptor class. Recent advances in GPCR structural biology have produced hundreds of GPCR structures solved by X-ray diffraction and increasingly, cryo-electron microscopy (cryo-EM). Many of these structures are stabilized by site-specific cholesterol binding, but it is unclear whether these interactions are a product of recurring cholesterol-binding motifs and if observed patterns of cholesterol binding differ by experimental technique. Here, we comprehensively analyze the location and composition of cholesterol binding sites in the current set of 473 human GPCR structural chains. Our findings establish that cholesterol binds similarly in cryo-EM and X-ray structures and show that 92% of cholesterol molecules on GPCR surfaces reside in predictable locations that lack discernable cholesterol-binding motifs.

Macromolecular refinement of X-ray and cryoelectron microscopy structures with Phenix/OPLS3e for improved structure and ligand quality

With the advent of the resolution revolution in cryoelectron microscopy (cryo-EM), low-resolution refinement is common, and likewise increases the need for a reliable force field. Here, we report on the incorporation of the OPLS3e force field with the VSGB2.1 solvation model in the structure determination package Phenix. Our results show significantly improved structure quality and reduced ligand strain at lower resolution for X-ray refinement. For refinement of cryo-EM-based structures, we find comparable quality structures, goodness-of-fit, and reduced ligand strain. We also show how structure quality and ligand strain are related to the map-model cross-correlation as a function of data weight, and how that can detect overfitting. Signs of overfitting are found in over half of our cryo-EM dataset, which can be remedied by a re-refinement at a lower data weight. Finally, a start-to-end script for refining structures with Phenix/OPLS3e is available in the Schrödinger 2020-3 distribution.

Enhanced validation of small-molecule ligands and carbohydrates in the Protein Data Bank

The Worldwide Protein Data Bank (wwPDB) has provided validation reports based on recommendations from community Validation Task Forces for structures in the PDB since 2013. To further enhance validation of small molecules as recommended from the 2016 Ligand Validation Workshop, wwPDB, Global Phasing Ltd., and the Noguchi Institute, recently formed a public/private partnership to incorporate some of their software tools into the wwPDB validation package. Augmented wwPDB validation report features include: two-dimensional (2D) diagrams of small-molecule ligands and carbohydrates, highlighting geometric validation outcomes; 2D topological diagrams of oligosaccharides present in branched entities generated using 2D Symbol Nomenclature for Glycan representation; and views of 3D electron density maps for ligands and carbohydrates, illustrating the goodness-of-fit between the atomic structure and experimental data (X-ray crystallographic structures only). These improvements will impact confidence in ligand conformation and ligand-macromolecular interactions that will aid in understanding biochemical function and contribute to small-molecule drug discovery.

Cryo-EM structure of a functional monomeric Photosystem I from Thermosynechococcus elongatus reveals red chlorophyll cluster

A high-resolution structure of trimeric cyanobacterial Photosystem I (PSI) from Thermosynechococcus elongatus was reported as the first atomic model of PSI almost 20 years ago. However, the monomeric PSI structure has not yet been reported despite long-standing interest in its structure and extensive spectroscopic characterization of the loss of red chlorophylls upon monomerization. Here, we describe the structure of monomeric PSI from Thermosynechococcus elongatus BP-1. Comparison with the trimer structure gave detailed insights into monomerization-induced changes in both the central trimerization domain and the peripheral regions of the complex. Monomerization-induced loss of red chlorophylls is assigned to a cluster of chlorophylls adjacent to PsaX. Based on our findings, we propose a role of PsaX in the stabilization of red chlorophylls and that lipids of the surrounding membrane present a major source of thermal energy for uphill excitation energy transfer from red chlorophylls to P700.

2-Mercaptomethyl-thiazolidines use conserved aromatic-S interactions to achieve broad-range inhibition of metallo-β-lactamases

Infections caused by multidrug resistant (MDR) bacteria are a major public health threat. Carbapenems are among the most potent antimicrobial agents that are commercially available to treat MDR bacteria. Bacterial production of carbapenem-hydrolysing metallo-β-lactamases (MBLs) challenges their safety and efficacy, with subclass B1 MBLs hydrolysing almost all β-lactam antibiotics. MBL inhibitors would fulfil an urgent clinical need by prolonging the lifetime of these life-saving drugs. Here we report the synthesis and activity of a series of 2-mercaptomethyl-thiazolidines (MMTZs), designed to replicate MBL interactions with reaction intermediates or hydrolysis products. MMTZs are potent competitive inhibitors of B1 MBLs in vitro (e.g., K i = 0.44 μM vs. NDM-1). Crystal structures of MMTZ complexes reveal similar binding patterns to the most clinically important B1 MBLs (NDM-1, VIM-2 and IMP-1), contrasting with previously studied thiol-based MBL inhibitors, such as bisthiazolidines (BTZs) or captopril stereoisomers, which exhibit lower, more variable potencies and multiple binding modes. MMTZ binding involves thiol coordination to the Zn(ii) site and extensive hydrophobic interactions, burying the inhibitor more deeply within the active site than d/l-captopril. Unexpectedly, MMTZ binding features a thioether-π interaction with a conserved active-site aromatic residue, consistent with their equipotent inhibition and similar binding to multiple MBLs. MMTZs penetrate multiple Enterobacterales, inhibit NDM-1 in situ, and restore carbapenem potency against clinical isolates expressing B1 MBLs. Based on their inhibitory profile and lack of eukaryotic cell toxicity, MMTZs represent a promising scaffold for MBL inhibitor development. These results also suggest sulphur-π interactions can be exploited for general ligand design in medicinal chemistry.

Lysozyme crystals dyed with bromophenol blue: where has the dye gone?

Protein crystals can easily be coloured by adding dyes to their mother liquor, but most structures of these protein-dye complexes remain unsolved. Here, structures of lysozyme in complex with bromophenol blue obtained by soaking orthorhombic and tetragonal crystals in a saturated solution of the dye at different pH values from 5.0 to 7.5 are reported. Two different binding sites can be found in the lysozyme-bromophenol blue crystals: binding site I is located near the amino- and carboxyl-termini, while binding site II is located adjacent to helices α1 (residues 4-15) and α3 (residues 88-100). In the orthorhombic crystals soaked at pH 7.0, binding of the dye takes place in both sites without significant changes in the unit cell. However, soaking tetragonal crystals with bromophenol blue results in two different complexes. Crystals soaked at pH 5.5 (HEWL-T1) show a single dye molecule bound to site II, and the crystals belong to space group P4₃2₁2 without significant changes in the unit cell (a = b = 78.50, c = 37.34 Å). On the other hand, crystals soaked at pH 6.5 in the presence of imidazole (HEWL-T2) show up to eight molecules of the dye bound to site II, and display changes in space group (P2₁2₁2₁) and unit cell (a = 38.00, b = 76.65, c = 84.86 Å). In all of the structures, the dye molecules are placed at the surface of the protein near to positively charged residues accessible through the main solvent channels of the crystal. Differences in the arrangement of the dye molecules at the surface of the protein suggest that the binding is not specific and is mainly driven by electrostatic interactions.

Combining native and 'omics' mass spectrometry to identify endogenous ligands bound to membrane proteins

Ligands bound to protein assemblies provide critical information for function, yet are often difficult to capture and define. Here we develop a top-down method, 'nativeomics', unifying 'omics' (lipidomics, proteomics, metabolomics) analysis with native mass spectrometry to identify ligands bound to membrane protein assemblies. By maintaining the link between proteins and ligands, we define the lipidome/metabolome in contact with membrane porins and a mitochondrial translocator to discover potential regulators of protein function.

Probing the Hydrogen Bond Involving Acridone Trapped in a Hydrophobic Biological Nanocavity: Integrated Spectroscopic and Docking Analyses

Spectroscopic analyses reveal that acridone (AD) penetrates through the structure and enters the hydrophobic cavity of the protein β-lactoglobulin (βLG). Although the protein contains two tryptophan (Trp) residues, AD interacts with only one (Trp-19), which is authenticated by the appearance of a single isoemissive point in TRANES. Alteration in the secondary structure of the protein while AD pierces through βLG is evident from the circular dichroism spectroscopic study. The ground-state interaction between AD and βLG is proven from the UV-vis spectroscopic study and the static nature of quenching of intrinsic fluorescence of the protein by the ligand. The steady-state fluorescence study in varied temperatures indicates the involvement of hydrogen bonding in the ligand-protein interaction. Further, the time-resolved fluorescence anisotropy study gives a hint of the presence of a hydrogen bond in AD-βLG interaction, which possibly involves the rotamers of Trp-19. In fact, the idea of involvement of rotamers of Trp-19 is obtained from the increase in fluorescence lifetime of βLG in the presence of AD. The docking study agrees to the involvement of hydrogen bonding in AD-βLG interaction. The direct evidence of hydrogen bonding between Trp and AD is obtained from the laser flash photolysis studies where the signature of formation of ADH^• and Trp^• through hydrogen abstraction between Trp and AD, loosely bound through hydrogen bonding, gets prominence. Thus, binding of AD to βLG involves hydrogen bonding in a hydrophobic pocket of the protein.

Ectopic suicide inhibition of thioredoxin glutathione reductase

Selective suicide inhibitors represent a seductively attractive approach for inactivation of therapeutically relevant enzymes since they are generally devoid of off-target toxicity in vivo. While most suicide inhibitors are converted to reactive species at enzyme active sites, theoretically bioactivation can also occur in ectopic (secondary) sites that have no known function. Here, we report an example of such an "ectopic suicide inhibition", an unprecedented bioactivation mechanism of a suicide inhibitor carried out by a non-catalytic site of thioredoxin glutathione reductase (TGR). TGR is a promising drug target to treat schistosomiasis, a devastating human parasitic disease. Utilizing hits selected from a high throughput screening campaign, time-resolved X-ray crystallography, molecular dynamics, mass spectrometry, molecular modeling, protein mutagenesis and functional studies, we find that 2-naphtholmethylamino derivatives bound to this novel ectopic site of Schistosoma mansoni (Sm)TGR are transformed to covalent modifiers and react with its mobile selenocysteine-containing C-terminal arm. In particular, one 2-naphtholmethylamino compound is able to specifically induce the pro-oxidant activity in the inhibited enzyme. Since some 2-naphtholmethylamino analogues show worm killing activity and the ectopic site is not conserved in human orthologues, a general approach to development of novel and selective anti-parasitic therapeutics against schistosoma is proposed.

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.

Discovery of New and Potent InhA Inhibitors as Antituberculosis Agents: Structure-Based Virtual Screening Validated by Biological Assays and X-ray Crystallography

The enoyl-acyl carrier protein reductase InhA of Mycobacterium tuberculosis is an attractive, validated target for antituberculosis drug development. Moreover, direct inhibitors of InhA remain effective against InhA variants with mutations associated with isoniazid resistance, offering the potential for activity against MDR isolates. Here, structure-based virtual screening supported by biological assays was applied to identify novel InhA inhibitors as potential antituberculosis agents. High-speed Glide SP docking was initially performed against two conformations of InhA differing in the orientation of the active site Tyr158. The resulting hits were filtered for drug-likeness based on Lipinski's rule and avoidance of PAINS-like properties and finally subjected to Glide XP docking to improve accuracy. Sixteen compounds were identified and selected for in vitro biological assays, of which two (compounds 1 and 7) showed MIC of 12.5 and 25 μg/mL against M. tuberculosis H37Rv, respectively. Inhibition assays against purified recombinant InhA determined IC₅₀ values for these compounds of 0.38 and 0.22 μM, respectively. A crystal structure of the most potent compound, compound 7, bound to InhA revealed the inhibitor to occupy a hydrophobic pocket implicated in binding the aliphatic portions of InhA substrates but distant from the NADH cofactor, i.e., in a site distinct from those occupied by the great majority of known InhA inhibitors. This compound provides an attractive starting template for ligand optimization aimed at discovery of new and effective compounds against M. tuberculosis that act by targeting InhA.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

High-throughput structures of protein-ligand complexes at room temperature using serial femtosecond crystallography

High-throughput X-ray crystal structures of protein-ligand complexes are critical to pharmaceutical drug development. However, cryocooling of crystals and X-ray radiation damage may distort the observed ligand binding. Serial femtosecond crystallography (SFX) using X-ray free-electron lasers (XFELs) can produce radiation-damage-free room-temperature structures. Ligand-binding studies using SFX have received only modest attention, partly owing to limited beamtime availability and the large quantity of sample that is required per structure determination. Here, a high-throughput approach to determine room-temperature damage-free structures with excellent sample and time efficiency is demonstrated, allowing complexes to be characterized rapidly and without prohibitive sample requirements. This yields high-quality difference density maps allowing unambiguous ligand placement. Crucially, it is demonstrated that ligands similar in size or smaller than those used in fragment-based drug design may be clearly identified in data sets obtained from <1000 diffraction images. This efficiency in both sample and XFEL beamtime opens the door to true high-throughput screening of protein-ligand complexes using SFX.

Finding High-Quality Metal Ion-Centric Regions Across the Worldwide Protein Data Bank

As the number of macromolecular structures in the worldwide Protein Data Bank (wwPDB) continues to grow rapidly, more attention is being paid to the quality of its data, especially for use in aggregated structural and dynamics analyses. In this study, we systematically analyzed 3.5 Å regions around all metal ions across all PDB entries with supporting electron density maps available from the PDB in Europe. All resulting metal ion-centric regions were evaluated with respect to four quality-control criteria involving electron density resolution, atom occupancy, symmetry atom exclusion, and regional electron density discrepancy. The resulting list of metal binding sites passing all four criteria possess high regional structural quality and should be beneficial to a wide variety of downstream analyses. This study demonstrates an approach for the pan-PDB evaluation of metal binding site structural quality with respect to underlying X-ray crystallographic experimental data represented in the available electron density maps of proteins. For non-crystallographers in particular, we hope to change the focus and discussion of structural quality from a global evaluation to a regional evaluation, since all structural entries in the wwPDB appear to have both regions of high and low structural quality.

Hydrogen Donation but not Abstraction by a Tyrosine (Y68) during Endoperoxide Installation by Verruculogen Synthase (FtmOx1)

Hydrogen-atom transfer (HAT) from a substrate carbon to an iron(IV)-oxo (ferryl) intermediate initiates a diverse array of enzymatic transformations. For outcomes other than hydroxylation, coupling of the resultant carbon radical and hydroxo ligand (oxygen rebound) must generally be averted. A recent study of FtmOx1, a fungal iron(II)- and 2-(oxo)glutarate-dependent oxygenase that installs the endoperoxide of verruculogen by adding O2 between carbons 21 and 27 of fumitremorgin B, posited that tyrosine (Tyr or Y) 224 serves as HAT intermediary to separate the C21 radical (C21•) and Fe(III)-OH HAT products and prevent rebound. Our reinvestigation of the FtmOx1 mechanism revealed, instead, direct HAT from C21 to the ferryl complex and surprisingly competitive rebound. The C21-hydroxylated (rebound) product, which undergoes deprenylation, predominates when low [O2] slows C21•-O2 coupling in the next step of the endoperoxidation pathway. This pathway culminates with addition of the C21-O-O• peroxyl adduct to olefinic C27 followed by HAT to the C26• from a Tyr. The last step results in sequential accumulation of Tyr radicals, which are suppressed without detriment to turnover by inclusion of the reductant, ascorbate. Replacement of each of four candidates for the proximal C26 H• donor (including Y224) with phenylalanine (F) revealed that only the Y68F variant (i) fails to accumulate the first Tyr• and (ii) makes an altered major product, identifying Y68 as the donor. The implied proximities of C21 to the iron cofactor and C26 to Y68 support a new docking model of the enzyme-substrate complex that is consistent with all available data.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

Protein Data Bank: the single global archive for 3D macromolecular structure data

The Protein Data Bank (PDB) is the single global archive of experimentally determined three-dimensional (3D) structure data of biological macromolecules. Since 2003, the PDB has been managed by the Worldwide Protein Data Bank (wwPDB; wwpdb.org), an international consortium that collaboratively oversees deposition, validation, biocuration, and open access dissemination of 3D macromolecular structure data. The PDB Core Archive houses 3D atomic coordinates of more than 144 000 structural models of proteins, DNA/RNA, and their complexes with metals and small molecules and related experimental data and metadata. Structure and experimental data/metadata are also stored in the PDB Core Archive using the readily extensible wwPDB PDBx/mmCIF master data format, which will continue to evolve as data/metadata from new experimental techniques and structure determination methods are incorporated by the wwPDB. Impacts of the recently developed universal wwPDB OneDep deposition/validation/biocuration system and various methods-specific wwPDB Validation Task Forces on improving the quality of structures and data housed in the PDB Core Archive are described together with current challenges and future plans.

Crystal structures of VIM-1 complexes explain active site heterogeneity in VIM-class metallo-β-lactamases

Metallo-β-Lactamases (MBLs) protect bacteria from almost all β-lactam antibiotics. Verona integron-encoded MBL (VIM) enzymes are among the most clinically important MBLs, with VIM-1 increasing in carbapenem-resistant Enterobacteriaceae (Escherichia coli, Klebsiella pneumoniae) that are among the hardest bacterial pathogens to treat. VIM enzymes display sequence variation at residues (224 and 228) that in related MBLs are conserved and participate in substrate binding. How they accommodate this variability, while retaining catalytic efficiency against a broad substrate range, has remained unclear. Here, we present crystal structures of VIM-1 and its complexes with a substrate-mimicking thioenolate inhibitor, ML302F, that restores meropenem activity against a range of VIM-1 producing clinical strains, and the hydrolysed product of the carbapenem meropenem. Comparison of these two structures identifies a water-mediated hydrogen bond, between the carboxylate group of substrate/inhibitor and the backbone carbonyl of the active site zinc ligand Cys221, that is common to both complexes. Structural comparisons show that the responsible Cys221-bound water is observed in all known VIM structures, participates in carboxylate binding with other inhibitor classes, and thus effectively replicates the role of the conserved Lys224 in analogous complexes with other MBLs. These results provide a mechanism for substrate binding that permits the variation at positions 224 and 228 that is a hallmark of VIM MBLs. ENZYMES: EC 3.5.2.6 DATABASES: Co-ordinates and structure factors for protein structures described in this manuscript have been deposited in the Protein Data Bank (www.rcsb.org/pdb) with accession codes 5N5G (VIM-1), 5N5H (VIM-1:ML302F complex) and 5N5I (VIM-1-hydrolysed meropenem complex).

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.

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."

Worldwide Protein Data Bank validation information: usage and trends

Realising the importance of assessing the quality of the biomolecular structures deposited in the Protein Data Bank (PDB), the Worldwide Protein Data Bank (wwPDB) partners established Validation Task Forces to obtain advice on the methods and standards to be used to validate structures determined by X-ray crystallography, nuclear magnetic resonance spectroscopy and three-dimensional electron cryo-microscopy. The resulting wwPDB validation pipeline is an integral part of the wwPDB OneDep deposition, biocuration and validation system. The wwPDB Validation Service webserver (https://validate.wwpdb.org) can be used to perform checks prior to deposition. Here, it is shown how validation metrics can be combined to produce an overall score that allows the ranking of macromolecular structures and domains in search results. The ValTrendsDB database provides users with a convenient way to access and analyse validation information and other properties of X-ray crystal structures in the PDB, including investigating trends in and correlations between different structure properties and validation metrics.