A workflow for deriving chemical entities from crystallographic data and its application to the Crystallography Open Database.

Knowledge about the 3-dimensional structure, orientation and interaction of chemical compounds is important in many areas of science and technology. X-ray crystallography is one of the experimental techniques capable of providing a large amount of structural information for a given compound, and it is widely used for characterisation of organic and metal-organic molecules. The method provides precise 3D coordinates of atoms inside crystals, however, it does not directly deliver information about certain chemical characteristics such as bond orders, delocalization, charges, lone electron pairs or lone electrons. These aspects of a molecular model have to be derived from crystallographic data using refined information about interatomic distances and atom types as well as employing general chemical knowledge. This publication describes a curated automatic pipeline for the derivation of chemical attributes of molecules from crystallographic models. The method is applied to build a catalogue of chemical entities in an open-access crystallographic database, the Crystallography Open Database (COD). The catalogue of such chemical entities is provided openly as a derived database. The content of this catalogue and the problems arising in the fully automated pipeline are discussed, along with the possibilities to introduce manual data curation into the process.

Building protein structure-specific rotamer libraries.

MOTIVATION: Identifying the probable positions of the protein side-chains is one of the protein modelling steps that can improve the prediction of protein-ligand and protein-protein interactions. Most of the strategies predicting the side-chain conformations use predetermined dihedral angle lists, also called rotamer libraries, that are usually generated from a subset of high-quality protein structures. Although these methods are fast to apply, they tend to average out geometries instead of taking into account the surrounding atoms and molecules and ignore structures not included in the selected subset. Such simplifications can result in inaccuracies when predicting possible side-chain atom positions. RESULTS: We propose an approach that takes into account both of these circumstances by scanning through sterically accessible side-chain conformations and generating dihedral angle libraries specific to the target proteins. The method avoids the drawbacks of lacking conformations due to unusual or rare protein structures and successfully suggests potential rotamers with average RMSD closer to the experimentally determined side-chain atom positions than other widely used rotamer libraries. AVAILABILITY AND IMPLEMENTATION: The technique is implemented in open-source software package rotag and available at GitHub: https://www.github.com/agrybauskas/rotag, under GNU Lesser General Public License.

PLBD: protein-ligand binding database of thermodynamic and kinetic intrinsic parameters.

We introduce a protein-ligand binding database (PLBD) that presents thermodynamic and kinetic data of reversible protein interactions with small molecule compounds. The manually curated binding data are linked to protein-ligand crystal structures, enabling structure-thermodynamics correlations to be determined. The database contains over 5500 binding datasets of 556 sulfonamide compound interactions with the 12 catalytically active human carbonic anhydrase isozymes defined by fluorescent thermal shift assay, isothermal titration calorimetry, inhibition of enzymatic activity and surface plasmon resonance. In the PLBD, the intrinsic thermodynamic parameters of interactions are provided, which account for the binding-linked protonation reactions. In addition to the protein-ligand binding affinities, the database provides calorimetrically measured binding enthalpies, providing additional mechanistic understanding. The PLBD can be applied to investigations of protein-ligand recognition and could be integrated into small molecule drug design. Database URL https://plbd.org/.

Structural basis for sequence-specific recognition of guide and target strands by the Archaeoglobus fulgidus Argonaute protein.

Argonaute (Ago) proteins are found in all three domains of life. The best-characterized group is eukaryotic Argonautes (eAgos). Being the structural core of RNA interference machinery, they use guide RNA molecules for RNA targeting. Prokaryotic Argonautes (pAgos) are more diverse, both in terms of structure (there are eAgo-like 'long' and truncated 'short' pAgos) and mechanism, as many pAgos are specific for DNA, not RNA guide and/or target strands. Some long pAgos act as antiviral defence systems. Their defensive role was recently demonstrated for short pAgo-encoding systems SPARTA and GsSir2/Ago, but the function and action mechanisms of all other short pAgos remain unknown. In this work, we focus on the guide and target strand preferences of AfAgo, a truncated long-B Argonaute protein encoded by an archaeon Archaeoglobus fulgidus. We demonstrate that AfAgo associates with small RNA molecules carrying 5'-terminal AUU nucleotides in vivo, and characterize its affinity to various RNA and DNA guide/target strands in vitro. We also present X-ray structures of AfAgo bound to oligoduplex DNAs that provide atomic details for base-specific AfAgo interactions with both guide and target strands. Our findings broaden the range of currently known Argonaute-nucleic acid recognition mechanisms.

Graph isomorphism-based algorithm for cross-checking chemical and crystallographic descriptions.

Published reports of chemical compounds often contain multiple machine-readable descriptions which may supplement each other in order to yield coherent and complete chemical representations. This publication presents a method to cross-check such descriptions using a canonical representation and isomorphism of molecular graphs. If immediate agreement between compound descriptions is not found, the algorithm derives the minimal set of simplifications required for both descriptions to arrive to a matching form (if any). The proposed algorithm is used to cross-check chemical descriptions from the Crystallography Open Database to identify coherently described entries as well as those requiring further curation.

Light-Responsive Oligothiophenes Incorporating Photochromic Torsional Switches.

We present a quaterthiophene and sexithiophene that can reversibly change their effective π-conjugation length through photoexcitation. The reported compounds make use of light-responsive molecular actuators consisting of an azobenzene attached to a bithiophene unit by both direct and linker-assisted bonding. Upon exposure to 350 nm light, the azobenzene undergoes trans-to-cis isomerization, thus mechanically inducing the oligothiophene to assume a planar conformation (extended π-conjugation). Exposure to 254 nm wavelength promotes azobenzene cis-to-trans isomerization, forcing the thiophenic backbones to twist out of planarity (confined π-conjugation). Twisted conformations are also reached by cis-to-trans thermal relaxation at a rate that increases proportionally with the conjugation length of the oligothiophene moiety. The molecular conformations of quaterthiophene and sexithiophene were characterized by using steady-state UV-vis spectroscopy, X-ray crystallography and quantum-chemical modeling. Finally, we tested the proposed light-responsive oligothiophenes in field-effect transistors to probe the photo-induced tuning of their electronic properties.

Beta and Gamma Amino Acid-Substituted Benzenesulfonamides as Inhibitors of Human Carbonic Anhydrases.

A series of novel benzenesulfonamide derivatives were synthesized bearing para-N ß,γ-amino acid or para-N ß-amino acid and thiazole moieties and their binding to the human carbonic anhydrase (CA) isozymes determined. These enzymes are involved in various illnesses, such as glaucoma, altitude sickness, epilepsy, obesity, and even cancer. There are numerous compounds that are inhibitors of CA and used as pharmaceuticals. However, most of them bind to most CA isozymes with little selectivity. The design of high affinity and selectivity towards one CA isozyme remains a significant challenge. The beta and gamma amino acid-substituted compound affinities were determined by the fluorescent thermal shift assay and isothermal titration calorimetry for all 12 catalytically active human carbonic anhydrase isozymes, showing the full affinity and selectivity profile. The structures of several compounds were determined by X-ray crystallography, and the binding mode in the active site of CA enzyme was shown.

OPTIMADE, an API for exchanging materials data.

The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API through worked examples on each of the public materials databases that support the full API specification.

Structure and mechanism of secondary sulfonamide binding to carbonic anhydrases.

Zinc-containing metalloenzyme carbonic anhydrase (CA) binds primary sulfonamides with extremely high, up to picomolar, affinity by forming a coordination bond between the negatively charged amino group and the zinc ion and making hydrogen bonds and hydrophobic contacts with other parts of the inhibitor molecule. However, N-methyl-substituted, secondary or tertiary sulfonamides bind CA with much lower affinity. In search for an explanation for this diminished affinity, a series of secondary sulfonamides were synthesized and, together with analogous primary sulfonamides, the affinities for 12 recombinant catalytically active human CA isoforms were determined by the fluorescent thermal shift assay, stopped-flow assay of the inhibition of enzymatic activity and isothermal titration calorimetry. The binding profile of secondary sulfonamides as a function of pH showed the same U-shape dependence seen for primary sulfonamides. This dependence demonstrated that there were protein binding-linked protonation reactions that should be dissected for the estimation of the intrinsic binding constants to perform structure-thermodynamics analysis. X-ray crystallographic structures of secondary sulfonamides and computational modeling dissected the atomic contributions to the binding energetics. Secondary sulfonamides bind to carbonic anhydrases via coordination bond between the negatively charged nitrogen of alkylated amino group and Zn(II) in the active site of CA. The binding reaction is linked to deprotonation of the amino group and protonation of the Zn(II)-bound hydroxide. To perform the structure-thermodynamics analysis, contributions of these linked reactions must be subtracted to determine the intrinsic energetics. In this aspect, the secondary sulfonamides are similar to primary sulfonamides as CA inhibitors.

Validation of the Crystallography Open Database using the Crystallographic Information Framework.

Data curation practices of the Crystallography Open Database (COD) are described with additional focus being placed on the formal validation using the Crystallographic Information Framework (CIF). The cif_validate program, capable of validating CIF files against both the DDL1 and the DDLm dictionaries, is presented and used to process the entirety of the COD. Validation results collected from over 450 000 CIF files are demonstrated to be a useful resource in the data maintenance process as well as the development of the underlying ontologies. A set of programs intended to aid in the dictionary migration from DDL1 to DDLm is also presented.

Crystal structure of restriction endonuclease Kpn2I of CCGG-family.

BACKGROUND: Restriction endonucleases belong to prokaryotic restriction-modification systems, that protect host cells from invading DNA. Type II restriction endonucleases recognize short 4-8 bp sequences in the target DNA and cut both DNA strands producing double strand breaks. Type II restriction endonuclease Kpn2I cleaves 5'-T/CCGGA DNA sequence ("/" marks the cleavage position). Analysis of protein sequences suggested that Kpn2I belongs to the CCGG-family, which contains ten enzymes that recognize diverse nucleotides outside the conserved 5'-CCGG core and share similar motifs for the 5'-CCGG recognition and cleavage. METHODS: We solved a crystal structure of Kpn2I in a DNA-free form at 2.88 Å resolution. From the crystal structure we predicted active center and DNA recognition residues and tested them by mutational analysis. We estimated oligomeric state of Kpn2I by SEC-MALS and performed plasmid DNA cleavage assay to elucidate DNA cleavage mechanism. RESULTS: Structure comparison confirmed that Kpn2I shares a conserved active site and structural determinants for the 5'-CCGG tetranucleotide recognition with other restriction endonucleases of the CCGG-family. Guided by structural similarity between Kpn2I and the CCGG-family restriction endonucleases PfoI and AgeI, Kpn2I residues involved in the outer base pair recognition were proposed. CONCLUSIONS: Kpn2I is an orthodox Type IIP restriction endonuclease, which acts as a dimer. Kpn2I shares structural similarity to the CCGG-family restriction endonucleases PfoI, AgeI and PspGI. GENERAL SIGNIFICANCE: The Kpn2I structure concluded the studies of the CCGG-family, covering detailed structural and biochemical characterization of eleven restriction enzymes and their complexes with DNA.

Switching the Inhibitor-Enzyme Recognition Profile via Chimeric Carbonic Anhydrase XII.

A key part of the optimization of small molecules in pharmaceutical inhibitor development is to vary the molecular design to enhance complementarity of chemical features of the compound with the positioning of amino acids in the active site of a target enzyme. Typically this involves iterations of synthesis, to modify the compound, and biophysical assay, to assess the outcomes. Selective targeting of the anti-cancer carbonic anhydrase isoform XII (CA XII), this process is challenging because the overall fold is very similar across the twelve CA isoforms. To enhance drug development for CA XII we used a reverse engineering approach where mutation of the key six amino acids in the active site of human CA XII into the CA II isoform was performed to provide a protein chimera (chCA XII) which is amenable to structure-based compound optimization. Through determination of structural detail and affinity measurement of the interaction with over 60 compounds we observed that the compounds that bound CA XII more strongly than CA II, switched their preference and bound more strongly to the engineered chimera, chCA XII, based on CA II, but containing the 6 key amino acids from CA XII, behaved as CA XII in its compound recognition profile. The structures of the compounds in the chimeric active site also resembled those determined for complexes with CA XII, hence validating this protein engineering approach in the development of new inhibitors.

Methyl 2-Halo-4-Substituted-5-Sulfamoyl-Benzoates as High Affinity and Selective Inhibitors of Carbonic Anhydrase IX.

Among the twelve catalytically active carbonic anhydrase isozymes present in the human body, the CAIX is highly overexpressed in various solid tumors. The enzyme acidifies the tumor microenvironment enabling invasion and metastatic processes. Therefore, many attempts have been made to design chemical compounds that would exhibit high affinity and selective binding to CAIX over the remaining eleven catalytically active CA isozymes to limit undesired side effects. It has been postulated that such drugs may have anticancer properties and could be used in tumor treatment. Here we have designed a series of compounds, methyl 5-sulfamoyl-benzoates, which bear a primary sulfonamide group, a well-known marker of CA inhibitors, and determined their affinities for all twelve CA isozymes. Variations of substituents on the benzenesulfonamide ring led to compound 4b, which exhibited an extremely high observed binding affinity to CAIX; the Kd was 0.12 nM. The intrinsic dissociation constant, where the binding-linked protonation reactions have been subtracted, reached 0.08 pM. The compound also exhibited more than 100-fold selectivity over the remaining CA isozymes. The X-ray crystallographic structure of compound 3b bound to CAIX showed the structural position, while several structures of compounds bound to other CA isozymes showed structural reasons for compound selectivity towards CAIX. Since this series of compounds possess physicochemical properties suitable for drugs, they may be developed for anticancer therapeutic purposes.

Isoform-Selective Enzyme Inhibitors by Exploring Pocket Size According to the Lock-and-Key Principle.

In the design of high-affinity and enzyme isoform-selective inhibitors, we applied an approach of augmenting the substituents attached to the benzenesulfonamide scaffold in three ways, namely, substitutions at the 3,5- or 2,4,6-positions or expansion of the condensed ring system. The increased size of the substituents determined the spatial limitations of the active sites of the 12 catalytically active human carbonic anhydrase (CA) isoforms until no binding was observed because of the inability of the compounds to fit in the active site. This approach led to the discovery of high-affinity and high-selectivity compounds for the anticancer target CA IX and antiobesity target CA VB. The x-ray crystallographic structures of compounds bound to CA IX showed the positions of the bound compounds, whereas computational modeling confirmed that steric clashes prevent the binding of these compounds to other isoforms and thus avoid undesired side effects. Such an approach, based on the Lock-and-Key principle, could be used for the development of enzyme-specific drug candidate compounds.

Halogenated and di-substituted benzenesulfonamides as selective inhibitors of carbonic anhydrase isoforms.

By applying an approach of a "ring with two tails", a series of novel inhibitors possessing high-affinity and significant selectivity towards selected carbonic anhydrase (CA) isoforms has been designed. The "ring" consists of 2-chloro/bromo-benzenesulfonamide, where the sulfonamide group is as an anchor coordinating the Zn(II) in the active site of CAs, and halogen atom orients the ring affecting the affinity and selectivity. First "tail" is a substituent containing carbonyl, carboxyl, hydroxyl, ether groups or hydrophilic amide linkage. The second "tail" contains aryl- or alkyl-substituents attached through a sulfanyl or sulfonyl group. Both "tails" are connected to the benzene ring and play a crucial role in selectivity. Varying the substituents, we designed compounds selective for CA VII, CA IX, CA XII, or CA XIV. Since due to binding-linked protonation reactions the binding-ready fractions of the compound and protein are much lower than one, the "intrinsic" affinities were calculated that should be used to study correlations between crystal structures and the thermodynamics of binding for rational drug design. The "intrinsic" affinities together with the intrinsic enthalpies and entropies of binding together with co-crystal structures were used demonstrate structural factors determining major contributions for compound affinity and selectivity.

Raman Open Database: first interconnected Raman-X-ray diffraction open-access resource for material identification.

Detailed crystallographic information provided by X-ray diffraction (XRD) is complementary to molecular information provided by Raman spectroscopy. Accordingly, the combined use of these techniques allows the identification of an unknown compound without ambiguity. However, a full combination of Raman and XRD results requires an appropriate and reliable reference database with complete information. This is already available for XRD. The main objective of this paper is to introduce and describe the recently developed Raman Open Database (ROD, http://solsa.crystallography.net/rod). It comprises a collection of high-quality uncorrected Raman spectra. The novelty of this database is its interconnectedness with other open databases like the Crystallography Open Database (http://www.crystallography.net/cod and Theoretical Crystallography Open Database (http://www.crystallography.net/tcod/). The syntax adopted to format entries in the ROD is based on the worldwide recognized and used CIF format, which offers a simple way for data exchange, writing and description. ROD also uses JCAMP-DX files as an alternative format for submitted spectra. JCAMP-DX files are compatible to varying degrees with most commercial Raman software and can be read and edited using standard text editors.

Unique mechanism of target recognition by PfoI restriction endonuclease of the CCGG-family.

Restriction endonucleases (REs) of the CCGG-family recognize a set of 4-8 bp target sequences that share a common CCGG or CCNGG core and possess PD…D/ExK nuclease fold. REs that interact with 5 bp sequence 5'-CCNGG flip the central N nucleotides and 'compress' the bound DNA to stack the inner base pairs to mimic the CCGG sequence. PfoI belongs to the CCGG-family and cleaves the 7 bp sequence 5'-T|CCNGGA ("|" designates cleavage position). We present here crystal structures of PfoI in free and DNA-bound forms that show unique active site arrangement and mechanism of sequence recognition. Structures and mutagenesis indicate that PfoI features a permuted E…ExD…K active site that differs from the consensus motif characteristic to other family members. Although PfoI also flips the central N nucleotides of the target sequence it does not 'compress' the bound DNA. Instead, PfoI induces a drastic change in DNA backbone conformation that shortens the distance between scissile phosphates to match that in the unperturbed CCGG sequence. Our data demonstrate the diversity and versatility of structural mechanisms employed by restriction enzymes for recognition of related DNA sequences.

Design of two-tail compounds with rotationally fixed benzenesulfonamide ring as inhibitors of carbonic anhydrases.

Rational design of compounds that would bind specific pockets of the target proteins is a difficult task in drug design. The 12 isoforms of catalytically active human carbonic anhydrases (CAs) have highly similar active sites that make it difficult to design inhibitors selective for one or several CA isoforms. A series of CA inhibitors based on 2-chloro/bromo-benzenesulfonamide that is largely fixed in the CA active site together with one or two tails yielded compounds that were synthesized and evaluated as inhibitors of CA isoforms. Introduction of a second tail had significant influence on the binding affinity and two-tailed compounds in most cases provided high affinity and selectivity for CA IX and CA XIV. The contacts between several compounds and CA amino acids were determined by X-ray crystallography. Together with the intrinsic enthalpy and entropy of binding they provided the structure-thermodynamics correlations for this series of compounds with the insight how to rationally build compounds with desired CA isoform as a target.

Using SMILES strings for the description of chemical connectivity in the Crystallography Open Database.

Computer descriptions of chemical molecular connectivity are necessary for searching chemical databases and for predicting chemical properties from molecular structure. In this article, the ongoing work to describe the chemical connectivity of entries contained in the Crystallography Open Database (COD) in SMILES format is reported. This collection of SMILES is publicly available for chemical (substructure) search or for any other purpose on an open-access basis, as is the COD itself. The conventions that have been followed for the representation of compounds that do not fit into the valence bond theory are outlined for the most frequently found cases. The procedure for getting the SMILES out of the CIF files starts with checking whether the atoms in the asymmetric unit are a chemically acceptable image of the compound. When they are not (molecule in a symmetry element, disorder, polymeric species,etc.), the previously published cif_molecule program is used to get such image in many cases. The program package Open Babel is then applied to get SMILES strings from the CIF files (either those directly taken from the COD or those produced by cif_molecule when applicable). The results are then checked and/or fixed by a human editor, in a computer-aided task that at present still consumes a great deal of human time. Even if the procedure still needs to be improved to make it more automatic (and hence faster), it has already yielded more than 160,000 curated chemical structures and the purpose of this article is to announce the existence of this work to the chemical community as well as to spread the use of its results.

Crystal structure correlations with the intrinsic thermodynamics of human carbonic anhydrase inhibitor binding.

The structure-thermodynamics correlation analysis was performed for a series of fluorine- and chlorine-substituted benzenesulfonamide inhibitors binding to several human carbonic anhydrase (CA) isoforms. The total of 24 crystal structures of 16 inhibitors bound to isoforms CA I, CA II, CA XII, and CA XIII provided the structural information of selective recognition between a compound and CA isoform. The binding thermodynamics of all structures was determined by the analysis of binding-linked protonation events, yielding the intrinsic parameters, i.e., the enthalpy, entropy, and Gibbs energy of binding. Inhibitor binding was compared within structurally similar pairs that differ by para- or meta-substituents enabling to obtain the contributing energies of ligand fragments. The pairs were divided into two groups. First, similar binders-the pairs that keep the same orientation of the benzene ring exhibited classical hydrophobic effect, a less exothermic enthalpy and a more favorable entropy upon addition of the hydrophobic fragments. Second, dissimilar binders-the pairs of binders that demonstrated altered positions of the benzene rings exhibited the non-classical hydrophobic effect, a more favorable enthalpy and variable entropy contribution. A deeper understanding of the energies contributing to the protein-ligand recognition should lead toward the eventual goal of rational drug design where chemical structures of ligands could be designed based on the target protein structure.