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.

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.

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.

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.

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.

Thermodynamic, kinetic, and structural parameterization of human carbonic anhydrase interactions toward enhanced inhibitor design.

The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure-thermodynamics and structure-kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure-energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure-thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design.

Restriction endonuclease AgeI is a monomer which dimerizes to cleave DNA.

Although all Type II restriction endonucleases catalyze phosphodiester bond hydrolysis within or close to their DNA target sites, they form different oligomeric assemblies ranging from monomers, dimers, tetramers to higher order oligomers to generate a double strand break in DNA. Type IIP restriction endonuclease AgeI recognizes a palindromic sequence 5΄-A/CCGGT-3΄ and cuts it ('/' denotes the cleavage site) producing staggered DNA ends. Here, we present crystal structures of AgeI in apo and DNA-bound forms. The structure of AgeI is similar to the restriction enzymes that share in their target sites a conserved CCGG tetranucleotide and a cleavage pattern. Structure analysis and biochemical data indicate, that AgeI is a monomer in the apo-form both in the crystal and in solution, however, it binds and cleaves the palindromic target site as a dimer. DNA cleavage mechanism of AgeI is novel among Type IIP restriction endonucleases.

Intrinsic thermodynamics of high affinity inhibitor binding to recombinant human carbonic anhydrase IV.

Membrane-associated carbonic anhydrase (CA) isoform IV participates in carbon metabolism and pH homeostasis and is implicated in the development of eye diseases such as retinitis pigmentosa and glaucoma. A series of substituted benzenesulfonamides were designed and their binding affinity to CA IV was determined by fluorescent thermal shift assay and isothermal titration calorimetry (ITC). Compound [(4-chloro-2-phenylsulfanyl-5-sulfamoyl-benzoyl)amino]propyl acetate (19) bound CA IV with the K d of 1.0 nM and exhibited significant selectivity over the remaining 11 human CA isoforms. The compound could be developed as a drug targeting CA IV. Various forms of recombinant CA IV were produced in Escherichia coli and mammalian cell cultures. Comparison of their temperature stability in various buffers and salt solutions demonstrated that CA IV is most stable at slightly alkaline conditions and at elevated sodium sulfate concentrations. High-resolution X-ray crystallographic structures of ortho-Cl and meta-thiazole-substituted benzene sulfonamide in complex with CA IV revealed the position of and interactions between the ligand and the protein. Sulfonamide inhibitor binding to CA IV is linked to several reactions-the deprotonation of the sulfonamide amino group, the protonation of CA-Zn(II)-bound hydroxide at the active site of CA IV, and the compensating reactions of the buffer. The dissection of binding-linked reactions yielded the intrinsic thermodynamic parameters, characterizing the interaction between CA IV and the sulfonamides in the binding-able protonation forms, including Gibbs energy, enthalpy, and entropy, that could be used for the characterization of binding to any CA in the process of drug design.

Functional significance of protein assemblies predicted by the crystal structure of the restriction endonuclease BsaWI.

Type II restriction endonuclease BsaWI recognizes a degenerated sequence 5'-W/CCGGW-3' (W stands for A or T, '/' denotes the cleavage site). It belongs to a large family of restriction enzymes that contain a conserved CCGG tetranucleotide in their target sites. These enzymes are arranged as dimers or tetramers, and require binding of one, two or three DNA targets for their optimal catalytic activity. Here, we present a crystal structure and biochemical characterization of the restriction endonuclease BsaWI. BsaWI is arranged as an 'open' configuration dimer and binds a single DNA copy through a minor groove contacts. In the crystal primary BsaWI dimers form an indefinite linear chain via the C-terminal domain contacts implying possible higher order aggregates. We show that in solution BsaWI protein exists in a dimer-tetramer-oligomer equilibrium, but in the presence of specific DNA forms a tetramer bound to two target sites. Site-directed mutagenesis and kinetic experiments show that BsaWI is active as a tetramer and requires two target sites for optimal activity. We propose BsaWI mechanism that shares common features both with dimeric Ecl18kI/SgrAI and bona fide tetrameric NgoMIV/SfiI enzymes.

Crystal structure of the R-protein of the multisubunit ATP-dependent restriction endonuclease NgoAVII.

The restriction endonuclease (REase) NgoAVII is composed of two proteins, R.NgoAVII and N.NgoAVII, and shares features of both Type II restriction enzymes and Type I/III ATP-dependent restriction enzymes (see accompanying paper Zaremba et al., 2014). Here we present crystal structures of the R.NgoAVII apo-protein and the R.NgoAVII C-terminal domain bound to a specific DNA. R.NgoAVII is composed of two domains: an N-terminal nucleolytic PLD domain; and a C-terminal B3-like DNA-binding domain identified previously in BfiI and EcoRII REases, and in plant transcription factors. Structural comparison of the B3-like domains of R.NgoAVII, EcoRII, BfiI and the plant transcription factors revealed a conserved DNA-binding surface comprised of N- and C-arms that together grip the DNA. The C-arms of R.NgoAVII, EcoRII, BfiI and plant B3 domains are similar in size, but the R.NgoAVII N-arm which makes the majority of the contacts to the target site is much longer. The overall structures of R.NgoAVII and BfiI are similar; however, whilst BfiI has stand-alone catalytic activity, R.NgoAVII requires an auxiliary cognate N.NgoAVII protein and ATP hydrolysis in order to cleave DNA at the target site. The structures we present will help formulate future experiments to explore the molecular mechanisms of intersubunit crosstalk that control DNA cleavage by R.NgoAVII and related endonucleases.

Structural insight into the specificity of the B3 DNA-binding domains provided by the co-crystal structure of the C-terminal fragment of BfiI restriction enzyme.

The B3 DNA-binding domains (DBDs) of plant transcription factors (TF) and DBDs of EcoRII and BfiI restriction endonucleases (EcoRII-N and BfiI-C) share a common structural fold, classified as the DNA-binding pseudobarrel. The B3 DBDs in the plant TFs recognize a diverse set of target sequences. The only available co-crystal structure of the B3-like DBD is that of EcoRII-N (recognition sequence 5'-CCTGG-3'). In order to understand the structural and molecular mechanisms of specificity of B3 DBDs, we have solved the crystal structure of BfiI-C (recognition sequence 5'-ACTGGG-3') complexed with 12-bp cognate oligoduplex. Structural comparison of BfiI-C-DNA and EcoRII-N-DNA complexes reveals a conserved DNA-binding mode and a conserved pattern of interactions with the phosphodiester backbone. The determinants of the target specificity are located in the loops that emanate from the conserved structural core. The BfiI-C-DNA structure presented here expands a range of templates for modeling of the DNA-bound complexes of the B3 family of plant TFs.

The recognition domain of the methyl-specific endonuclease McrBC flips out 5-methylcytosine.

DNA cytosine methylation is a widespread epigenetic mark. Biological effects of DNA methylation are mediated by the proteins that preferentially bind to 5-methylcytosine (5mC) in different sequence contexts. Until now two different structural mechanisms have been established for 5mC recognition in eukaryotes; however, it is still unknown how discrimination of the 5mC modification is achieved in prokaryotes. Here we report the crystal structure of the N-terminal DNA-binding domain (McrB-N) of the methyl-specific endonuclease McrBC from Escherichia coli. The McrB-N protein shows a novel DNA-binding fold adapted for 5mC-recognition. In the McrB-N structure in complex with methylated DNA, the 5mC base is flipped out from the DNA duplex and positioned within a binding pocket. Base flipping elegantly explains why McrBC system restricts only T4-even phages impaired in glycosylation [Luria, S.E. and Human, M.L. (1952) A nonhereditary, host-induced variation of bacterial viruses. J. Bacteriol., 64, 557-569]: flipped out 5-hydroxymethylcytosine is accommodated in the binding pocket but there is no room for the glycosylated base. The mechanism for 5mC recognition employed by McrB-N is highly reminiscent of that for eukaryotic SRA domains, despite the differences in their protein folds.

Structural mechanisms of the degenerate sequence recognition by Bse634I restriction endonuclease.

Restriction endonuclease Bse634I recognizes and cleaves the degenerate DNA sequence 5'-R/CCGGY-3' (R stands for A or G; Y for T or C, '/' indicates a cleavage position). Here, we report the crystal structures of the Bse634I R226A mutant complexed with cognate oligoduplexes containing ACCGGT and GCCGGC sites, respectively. In the crystal, all potential H-bond donor and acceptor atoms on the base edges of the conserved CCGG core are engaged in the interactions with Bse634I amino acid residues located on the α6 helix. In contrast, direct contacts between the protein and outer base pairs are limited to van der Waals contact between the purine nucleobase and Pro203 residue in the major groove and a single H-bond between the O2 atom of the outer pyrimidine and the side chain of the Asn73 residue in the minor groove. Structural data coupled with biochemical experiments suggest that both van der Waals interactions and indirect readout contribute to the discrimination of the degenerate base pair by Bse634I. Structure comparison between related enzymes Bse634I (R/CCGGY), NgoMIV (G/CCGGC) and SgrAI (CR/CCGGYG) reveals how different specificities are achieved within a conserved structural core.

Crystallography Open Database (COD): an open-access collection of crystal structures and platform for world-wide collaboration.

Using an open-access distribution model, the Crystallography Open Database (COD, http://www.crystallography.net) collects all known 'small molecule / small to medium sized unit cell' crystal structures and makes them available freely on the Internet. As of today, the COD has aggregated ~150,000 structures, offering basic search capabilities and the possibility to download the whole database, or parts thereof using a variety of standard open communication protocols. A newly developed website provides capabilities for all registered users to deposit published and so far unpublished structures as personal communications or pre-publication depositions. Such a setup enables extension of the COD database by many users simultaneously. This increases the possibilities for growth of the COD database, and is the first step towards establishing a world wide Internet-based collaborative platform dedicated to the collection and curation of structural knowledge.

4-amino-substituted benzenesulfonamides as inhibitors of human carbonic anhydrases.

A series of N-aryl-ß-alanine derivatives and diazobenzenesulfonamides containing aliphatic rings were designed, synthesized, and their binding to carbonic anhydrases (CA) I, II, VI, VII, XII, and XIII was studied by the fluorescent thermal shift assay and isothermal titration calorimetry. The results showed that 4-substituted diazobenzenesulfonamides were more potent CA binders than N-aryl-ß-alanine derivatives. Most of the N-aryl-ß-alanine derivatives showed better affinity for CA II while diazobenzenesulfonamides possessed nanomolar affinities towards CA I isozyme. X-ray crystallographic structures showed the modes of binding of both compound groups.

From X-ray crystallographic structure to intrinsic thermodynamics of protein-ligand binding using carbonic anhydrase isozymes as a model system.

Carbonic anhydrase (CA) was among the first proteins whose X-ray crystal structure was solved to atomic resolution. CA proteins have essentially the same fold and similar active centers that differ in only several amino acids. Primary sulfonamides are well defined, strong and specific binders of CA. However, minor variations in chemical structure can significantly alter their binding properties. Over 1000 sulfonamides have been designed, synthesized and evaluated to understand the correlations between the structure and thermodynamics of their binding to the human CA isozyme family. Compound binding was determined by several binding assays: fluorescence-based thermal shift assay, stopped-flow enzyme activity inhibition assay, isothermal titration calorimetry and competition assay for enzyme expressed on cancer cell surfaces. All assays have advantages and limitations but are necessary for deeper characterization of these protein-ligand interactions. Here, the concept and importance of intrinsic binding thermodynamics is emphasized and the role of structure-thermodynamics correlations for the novel inhibitors of CA IX is discussed - an isozyme that is overexpressed in solid hypoxic tumors, and thus these inhibitors may serve as anticancer drugs. The abundant structural and thermodynamic data are assembled into the Protein-Ligand Binding Database to understand general protein-ligand recognition principles that could be used in drug discovery.

4-Substituted-2,3,5,6-tetrafluorobenzenesulfonamides as inhibitors of carbonic anhydrases I, II, VII, XII, and XIII.

A series of 4-substituted-2,3,5,6-tetrafluorobenezenesulfonamides were synthesized and their binding potencies as inhibitors of recombinant human carbonic anhydrase isozymes I, II, VII, XII, and XIII were determined by the thermal shift assay, isothermal titration calorimetry, and stop-flow CO2 hydration assay. All fluorinated benzenesulfonamides exhibited nanomolar binding potency toward tested CAs and fluorinated benzenesulfonamides posessed higher binding potency than non-fluorinated compounds. The crystal structures of 4-[(4,6-dimethylpyrimidin-2-yl)thio]-2,3,5,6-tetrafluorobenzenesulfonamide in complex with CA II and CA XII, and 2,3,5,6-tetrafluoro-4-[(2-hydroxyethyl)sulfonyl]benzenesulfonamide in complex with CA XIII were determined. The observed dissociation constants for several fluorinated compounds reached subnanomolar range for CA I isozyme. The affinity and the selectivity of the compounds towards tested isozymes are presented.

Benzenesulfonamides with pyrimidine moiety as inhibitors of human carbonic anhydrases I, II, VI, VII, XII, and XIII.

Two groups of benzenesulfonamide derivatives, bearing pyrimidine moieties, were designed and synthesized as inhibitors of carbonic anhydrases (CA). Their binding affinities to six recombinant human CA isoforms I, II, VI, VII, XII, and XIII were determined by the thermal shift assay (TSA). The binding of several inhibitors was measured by isothermal titration calorimetry (ITC). Direct demonstration of compound inhibition was achieved by determining the inhibition constant by stopped-flow CO2 hydration assay. The most potent compounds demonstrated selectivity towards isoform I and affinities of 0.5 nM. The crystal structures of selected compounds in complex with CA II, XII, and XIII were determined to atomic resolution. Compounds described here were compared with previously published pyrimidinebenzenesulfonamides.(1) Systematic structure-activity analysis of 40 compound interactions with six isoforms yields clues for the design of compounds with greater affinities and selectivities towards target CA isoforms.

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.