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

2021 ◽  
Vol 22 (13) ◽  
pp. 6830
Author(s):  
Sohini Chakraborti ◽  
Kaushik Hatti ◽  
Narayanaswamy Srinivasan

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.

Blood ◽  
2018 ◽  
Vol 132 (Supplement 1) ◽  
pp. 5802-5802
Author(s):  
Derek Essegian ◽  
Tyler A. Cunningham ◽  
Jonathan H. Schatz ◽  
Stephan Schurer

Abstract Targeted signaling inhibitors for hematologic malignancies often lead to limited clinical efficacy due to the outgrowth of subpopulations with alternative pathways independent of the drug target. The eIF4F complex responsible for translation initiation is a convergence point for many cancer-promoting signaling pathways and its inhibition leads to decreased expression of key oncoproteins and apoptosis. Lymphomas and leukemias show particular dependence on constitutive eIF4F activation. Indeed, natural compounds targeting the eIF4F enzymatic component, eIF4A1, demonstrate activities both in vitro and in vivo against lymphoma and leukemia model systems, among other tumor types. eIF4A1 is a noteworthy target for hematologic malignancies based on the finding that BCR stimulation leads to increased mRNA translation primary CLL patient samples. Additionally, eIF4F components eIF4A1 and eIF4G1 had increased expression upon IgM-induced BCR activation. The natural compound silvestrol is a potent inhibitor of eIF4A1, results in cancer cell cytotoxicity, and has an established therapeutic window in vivo. Silvestrol shows potent antitumor activity against 924 pan-cancer tumor cell lines with 830/924 (90%) sensitive at IC50 <100nM with lymphoma and leukemia cell lines being particularly sensitive. Silvestrol and other natural compounds, however, lack core drug-like properties and synthetic tractability. To discover new, specific and tractable inhibitors of eIF4A1 that are more drug-like, we have constructed several molecular models that we used to virtually screen more than 20 million compounds. eIF4A1 is the founding member of the DEAD-box RNA helicases, which include its paralogs eIF4A2 (91% amino-acid identity with eIF4A1) and eIF4A3 (60% identity). All DEAD-box helicases contain two RecA-like domains separated by a flexible linker. The cleft between these domains is lined with helicase motifs that mediate nucleotide binding and hydrolysis. In an absence of RNA or nucleotide, eIF4A proteins adopt diffuse open conformations; binding of RNA and ATP triggers transition to a more stable closed state. Modeling small-molecule interactions in the nucleotide cleft of eIF4A1 therefore assesses ability of molecules to lock eIF4A1 in a conformation unable to cycle through ATPase and helicase activities. Although no experimentally derived structure of human eIF4A1 co-crystalized with ATP exists, crystal structures of other DEAD-box family members with similar motifs permit detailed studies of nucleotide and ligand-binding and the development of homology models. We have used four available high-resolution crystal structures to build models predicting interactions of small molecules in the interdomain nucleotide-binding cleft. We identified nucleotide binding-site residues and accurately reproduced ATP interactions for all four models (derived from PDB: 2J0S, 1FUU, 2VSO, 2DB3). We then performed all-atom explicit-water molecular dynamics (MD) simulations for 500-700 ns to study conformational dynamics and atomic interactions of ATP-bound and ATP-unbound states. Yeast eif4A crystal structure (PDB:1FUU) in the open state, for example, illustrated closure of the two RecA domains upon ATP binding. As expected, ATP makes strong interactions with the N-terminal, while phosphate groups extend to the C-terminal interacting with arginines, bringing the two RecA domains together. In contrast, we did not observe domain closure in the same simulation with 1FUU without ATP bound. We also assessed 2J0S, a crystal structure of human closed eIF4A3 bound to ANP. ATP docked to this active site followed by 500 ns MD held the protein in the closed state with several interdomain interactions. Upon nucleotide removal, marked RecA separation occurred. We observed similar domain closure and opening for PDB: 2VSO and 2DB3. Pooling these results, we constructed two homology models of human eIF4A1 with both open (2VSO, 1FUU) and closed conformations (2J0S) as structural templates. We therefore have developed accurate and novel in silico models of eIF4A1 highly useful in assessing interactions of small-molecule ATPase inhibitors, with focus on the ATP-binding cleft. Disclosures No relevant conflicts of interest to declare.


2014 ◽  
Vol 70 (a1) ◽  
pp. C631-C631
Author(s):  
Elena Boldyreva

Supramolecular interactions in the solid state attract much attention. Different experimental and computational approaches are used, to predict and to design crystal structures, to predict the properties based on molecular and crystal structures, to range different types of intermolecular interactions. Analysis of the crystal structures at fixed (e.g. ambient) temperature and pressure conditions is most common for experiments, whereas most DFT calculations are limited to 0 K, to minimize computational costs. At the same time, evolution of a crystal structure as a function of experimental conditions can contribute significantly to understanding the structure-forming role and relative energies of different types of intermolecular interactions in the same crystal structure and of similar interactions in a series of different but structurally or chemically related compounds. In the present invited contribution I attempt to illustrate this using several selected examples from my own practice and from the papers published by other research groups. I consider, in particular, the results of variable-temperature and variable-pressure studies of continuous lattice strain and phase transitions in small-molecule organic compounds, the results of variable-temperature and variable-pressure crystallization, the results of comparing the dissolution profiles of mono- and multi-component small-molecule organic crystals. I shall also discuss how variable-temperature and variable-pressure experimental diffraction data can assist in optimizing the calculations aimed at comparing the relative stability of polymorphs and predicting polymorph transitions. The study was supported by Russian Ministry of Science and Education and Russian Academy of Sciences.


2002 ◽  
Vol 58 (3) ◽  
pp. 380-388 ◽  
Author(s):  
Frank H. Allen

The Cambridge Structural Database (CSD) now contains data for more than a quarter of a million small-molecule crystal structures. The information content of the CSD, together with methods for data acquisition, processing and validation, are summarized, with particular emphasis on the chemical information added by CSD editors. Nearly 80% of new structural data arrives electronically, mostly in CIF format, and the CCDC acts as the official crystal structure data depository for 51 major journals. The CCDC now maintains both a CIF archive (more than 73000 CIFs dating from 1996), as well as the distributed binary CSD archive; the availability of data in both archives is discussed. A statistical survey of the CSD is also presented and projections concerning future accession rates indicate that the CSD will contain at least 500000 crystal structures by the year 2010.


2014 ◽  
Vol 70 (a1) ◽  
pp. C495-C495
Author(s):  
Colin Groom ◽  
Suzanna Ward ◽  
Neil Feeder ◽  
Elna Pidcock ◽  
Peter Wood ◽  
...  

The crystallographic community has done something remarkable and almost unique in science. It has operated in such a way that the data generated in virtually every experiment reported in a publication is available to all. This data, in the form of individual crystal structures, is valuable not just in itself, but as a collection. To fully exploit the results of a new structure determination, we never analyse a single structure, we analyse it in the context of every previous crystal structure. Our knowledge of molecular geometry and molecular interactions derived from these structural databases is put to routine use in almost every field of chemistry. This presentation will specifically highlight what we can learn from the world's database of small molecule crystal structures and demonstrate how we can apply that knowledge not just to increase our understanding in structural chemistry, but in structural biology too.


2020 ◽  
Author(s):  
Eric Boittier ◽  
Jed Burns ◽  
Neha Gandhi ◽  
Vito Ferro

Glycosaminoglycans (GAGs) are a family of anionic carbohydrates that play an essential role in the physiology and pathology of all eukaryotic life. Experimental determination of GAG-protein complexes remains difficult due to the considerable diversity in both carbohydrate linkage, and sulfation patterns. To complement existing methods of structural determination, we present our molecular docking tool, GlycoTorchVina (GTV), which demonstrates a substantial improvement at reproducing low energy conformations of GAGs compared with traditional docking programs. Based on the carbohydrate specific docking program VinaCarb (VC), GTV utilizes rotational energy functions, calculated using density functional theory (DFT), specifically designed for glycosidic linkages found in GAGs. The redocking accuracy of four programs (GTV, VC, AutoDock Vina and Glide) was tested over a set of 10 high-quality crystal structures containing co-crystallized GAGs (tetrasaccharides or longer). GTV outperformed other programs and was able to reproduce the native pose of eight structures and produced top-scoring docked poses that were on average only 1.8 Å RMSD away from the crystal structure. Although imitation of crystal structures is a standard test used for assessing the accuracy of docking programs, we illustrate how the underlying quality of the crystal structure, which is often overlooked during benchmarking, affects conclusions drawn from this approach. Statistical and theoretical investigations into charge-charge (“salt-bridge”) interactions are also presented. Again, DFT calculations were used to derive non-bonded potentials describing salt-bridges, and solvent-mediated charge-charge (“water-bridge”) interactions. These data suggest that water-bridges play an important, yet poorly understood, role in the structures of GAG-protein complexes.


CrystEngComm ◽  
2021 ◽  
Author(s):  
Leen N. Kalash ◽  
Jason C. Cole ◽  
Royston C. B. Copley ◽  
Colin M. Edge ◽  
Alexandru A. Moldovan ◽  
...  

Analysis of the molecular and structural features of the GSK crystal structure database and Cambridge Structural Database leads to improved reliability in hydrogen bond propensity models for pharmaceutical polymorphs.


Blood ◽  
2005 ◽  
Vol 106 (11) ◽  
pp. 69-69 ◽  
Author(s):  
Titus J. Boggon

Abstract Janus (Jak) family non-receptor tyrosine kinases are critical for appropriate signaling of many growth factors and cytokines. The four vertebrate Jak kinase family members demonstrate differential receptor cytoplasmic tail binding associations and transduce discrete signals on extracellular binding of ligand to the transmembrane cytokine or growth factor receptor. On ligand binding, a rapid tyrosine phosphorylation mediated signaling cascade is initiated, culminating in translocation of cytosolic latent transcription factors, signal transducer and activator of transcription (Stat) proteins, to the nucleus and targeted activation of transcription. Dysregulation of this Jak-mediated signaling pathway is documented in a number of hematological diseases: Improper upregulation of Jak activity is seen in certain hematological malignancies [1] and inability to appropriately transduce signals from the common gamma chain, γc, through Jak3 is responsible for approximately 60% of human severe combined immunodeficiency cases [2]. In addition, Jak2 point mutation V617F is frequently documented in myeloproliferative disorders including polycythemia vera; this activating mutation may disrupt an autoinhibited conformation [3]. Therapies targeting restraint of Jak family tyrosine kinase activity may be useful for treating inappropriate activation of Jak signaling cascades or for suppressing the immune response. Advances towards structure-directed drug design of Jak-specific inhibitors were made recently with solution of the Jak3 kinase domain X-ray crystal structure [4], representing the first three-dimensional structural data for any portion of the Jak family of tyrosine kinases. Here three further crystal structures of the kinase domain of Jak3 are presented: an improved resolution co-crystal structure with staurosporine analog AFN-941 and two crystal forms of Jak3 kinase domain in complex with the kinase inhibitor small molecule compound QAD-409. Comparisons between these three solved Jak3 kinase domain crystal structures illustrate conformational flexibility between the kinase domain lobes and in the area of the catalytic cleft. Further structure analysis is also presented documenting in silico modeling of the binding of small molecule CP-690,550 [5] to the different Jak3 kinase domain crystal forms. Potential binding conformations of this inhibitor to the Jak3 kinase domain are suggested with one highly scored binding conformation predicted for all crystal forms. The crystal structures and modeling studies presented further define the extent of the Jak kinase catalytic cleft, demonstrate conformational plasticity in the active conformation Jak3 kinase domain and will aid the design of higher specificity Jak inhibitors.


1999 ◽  
Vol 32 (5) ◽  
pp. 963-967 ◽  
Author(s):  
Angela Altomare ◽  
Carmelo Giacovazzo ◽  
Antonietta Guagliardi ◽  
Anna Grazia Giuseppine Moliterni ◽  
Rosanna Rizzi

In direct procedures for crystal structure solution from powder data, information on the location and orientation of a molecular fragment may readily become available. Such information may be used retrospectively to improve the powder-pattern decomposition, with favourable effects on the phasing process. A method is described by which accurate estimation of a large number of structure-factor moduli is possible by exploiting the prior partial structural information.


2007 ◽  
Vol 22 (1) ◽  
pp. 27-34 ◽  
Author(s):  
A. Gómez ◽  
J. Rodríguez-Hernández ◽  
E. Reguera

A new structural model is proposed for cubic nitroprussides and the crystal structure for the complex salts of Fe(2+), Co(2+), and Ni(2+) refined in that model. In cubic nitroprussides the building unit, [Fe(CN)5NO]2−, and the assembling metal (M=Fe2+, Co2+, Ni2+), have ¾ occupancy with three formula units per cell (Z=3). This leads to certain structural disorder and to different local environments for the outer metal. The crystallographic results are supported by the Mössbauer and infrared data. The XRD powder patterns, index in a cubic cell (Fm3m space group), show a sinuous background because of diffuse scattering from positional disorder of the metal centers. Because of this, the crystal structures were refined allowing the metal centers to move from the (0,0,0) and (0,0,1/2) positions (away from positional symmetry restrictions). The refinement under these conditions leads to excellent agreement factors (Rwp, Rp, S), good pattern background fitting, and produced a refined structural model consistent with the crystal chemistry of nitroprussides. The studied materials are obtained as hydrates. On heating, the crystal water evolves, and below 100°C an anhydrous phase is obtained, preserving the framework of the original hydrates. The loss of the crystal water leads to cell contraction that represents around 2% of cell volume reduction. On cooling down from room temperature to 77 and 12 K, a slight expansion for the -M-N≡C-Fe-C≡N-M- chain length is observed, suggesting that at low temperature and reduction in the metals charge delocalization on the CN bridges takes place. For M=Fe and Co the crystal structure was also refined for the anhydrous phase at 12, 77, and 300 K.


2018 ◽  
Vol 293 (47) ◽  
pp. 18337-18352 ◽  
Author(s):  
Zongyang Lv ◽  
Katelyn M. Williams ◽  
Lingmin Yuan ◽  
James H. Atkison ◽  
Shaun K. Olsen

Ubiquitin (Ub) signaling plays a key regulatory role in nearly every aspect of eukaryotic biology and is initiated by E1 enzymes that activate and transfer Ub to E2 Ub-conjugating enzymes. Despite Ub E1's fundamental importance to the cell and its attractiveness as a target for therapeutic intervention in cancer and other diseases, its only available structural information is derived from yeast orthologs of human ubiquitin-like modifier–activating enzyme 1 (hUBA1). To illuminate structural differences between yeast and hUBA1 structures that might be exploited for the development of small-molecule therapeutics, we determined the first crystal structure of a hUBA1–Ub complex. Using structural analysis, molecular modeling, and biochemical analysis, we demonstrate that hUBA1 shares a conserved overall structure and mechanism with previously characterized yeast orthologs, but displays subtle structural differences, particularly within the active site. Computational analysis revealed four potential ligand-binding hot spots on the surface of hUBA1 that might serve as targets to inhibit hUBA1 at the level of Ub activation or E2 recruitment or that might potentially be used in approaches such as protein-targeting chimeric molecules. Taken together, our work enhances our understanding of the hUBA1 mechanism, provides an improved framework for the development of small-molecule inhibitors of UBA1, and serves as a stepping stone for structural studies that involve the enzymes of the human Ub system at the level of both E1 and E2.


Sign in / Sign up

Export Citation Format

Share Document