SuperStar: comparison of CSD and PDB-based interaction fields as a basis for the prediction of protein-ligand interactions.

SuperStar is an empirical method for identifying interaction sites in proteins, based entirely on the experimental information about non-bonded interactions, present in the IsoStar database. The interaction information in IsoStar is contained in scatterplots, which show the distribution of a chosen probe around structure fragments. SuperStar breaks a template molecule (e.g. a protein binding site) into structural fragments which correspond to those in the scatterplots. The scatterplots are then superimposed on the corresponding parts of the template and converted into a composite propensity map.The original version of SuperStar was based entirely on scatterplots from the CSD. Here, scatterplots based on protein-ligand interactions are implemented in SuperStar, and validated on a test set of 122 X-ray structures of protein-ligand complexes. In this validation, propensity maps are compared with the experimentally observed positions of ligand atoms of comparable types. Although non-bonded interaction geometries in small molecule structures are similar to those found in protein-ligand complexes, their relative frequencies of occurrence are different. Polar interactions are more common in the first class of structures, while interactions between hydrophobic groups are more common in protein crystals. In general, PDB and CSD-based SuperStar maps appear equally successful in the prediction of protein-ligand interactions. PDB-based maps are more suitable to identify hydrophobic pockets, and inherently take into account the experimental uncertainties of protein atomic positions. If the protonation state of a histidine, aspartate or glutamate protein side-chain is known, specific CSD-based maps for that protonation state are preferred over PDB-based maps which represent an ensemble of protonation states.     

SuperStar: a knowledge-based approach for identifying interaction sites in proteins.

An empirical method for identifying interaction sites in proteins is described and validated. The method is based entirely on experimental information about non-bonded interactions occurring in small-molecule crystal structures. These data are used in the form of scatterplots that show the experimentally observed distribution of one functional group (the "contact group" or "probe") around another. A template molecule (e.g. a protein binding site) is broken down into structure fragments and the scatterplots, showing the distribution of a chosen probe around these structure fragments, are superimposed on the corresponding parts of the template. The scatterplots are then translated into a three-dimensional map that shows the propensity of the probe at different positions around the template molecule. The method is illustrated for l -arabinose-binding protein, complexed with l -arabinose and with d -fucose, and for dihydrofolate reductase complexed with methotrexate. The method is validated on 122 X-ray structures of protein-ligand complexes. For all the binding sites of these proteins, propensity maps are generated for four different probes: a charged NH+3nitrogen, a carbonyl oxygen, a hydroxyl oxygen and a methyl carbon atom. Next, the maps are compared with the experimentally observed positions of ligand atoms of these types. For 74% of these ligand atoms (84% of the solvent-inaccessible ones) the calculated propensity of the matching probe at the experimental positions is higher than expected by chance. For 68% of the atoms (82% of the solvent-inaccessible ones) the propensity of the matching probe is higher than that of the other three probes. These results indicate that the approach generally gives good predictions for protein-ligand interactions. The potential applications of the propensity maps range from an aid in manual docking and structure-based drug design to their use in pharmacophore development.     

A new test set for validating predictions of protein-ligand interaction

We present a large test set of protein-ligand complexes for the purpose of validating algorithms that rely on the prediction of protein-ligand interactions. The set consists of 305 complexes with protonation states assigned by manual inspection. The following checks have been carried out to identify unsuitable entries in this set: (1) assessing the involvement of crystallographically related protein units in ligand binding; (2) identification of bad clashes between protein side chains and ligand; and (3) assessment of structural errors, and/or inconsistency of ligand placement with crystal structure electron density. In addition, the set has been pruned to assure diversity in terms of protein-ligand structures, and subsets are supplied for different protein-structure resolution ranges. A classification of the set by protein type is available. As an illustration, validation results are shown for GOLD and SuperStar. GOLD is a program that performs flexible protein-ligand docking, and SuperStar is used for the prediction of favorable interaction sites in proteins. The new CCDC/Astex test set is freely available to the scientific community (http://www.ccdc.cam.ac.uk).     

A fragment library based on Gaussian mixtures predicting favorable molecular interactions

Here, a protein atom-ligand fragment interaction library is described. The library is based on experimentally solved structures of protein-ligand and protein-protein complexes deposited in the Protein Data Bank (PDB) and it is able to characterize binding sites given a ligand structure suitable for a protein. A set of 30 ligand fragment types were defined to include three or more atoms in order to unambiguously define a frame of reference for interactions of ligand atoms with their receptor proteins. Interactions between ligand fragments and 24 classes of protein target atoms plus a water oxygen atom were collected and segregated according to type. The spatial distributions of individual fragment - target atom pairs were visually inspected in order to obtain rough-grained constraints on the interaction volumes. Data fulfilling these constraints were given as input to an iterative expectation-maximization algorithm that produces as output maximum likelihood estimates of the parameters of the finite Gaussian mixture models. Concepts of statistical pattern recognition and the resulting mixture model densities are used (i) to predict the detailed interactions between Chlorella virus DNA ligase and the adenine ring of its ligand and (ii) to evaluate the "error" in prediction for both the training and validation sets of protein-ligand interaction found in the PDB. These analyses demonstrate that this approach can successfully narrow down the possibilities for both the interacting protein atom type and its location relative to a ligand fragment.     

Automated prediction of ligand-binding sites in proteins

We present a method, termed AutoLigand, for the prediction of ligand-binding sites in proteins of known structure. The method searches the space surrounding the protein and finds the contiguous envelope with the specified volume of atoms, which has the largest possible interaction energy with the protein. It uses a full atomic representation, with atom types for carbon, hydrogen, oxygen, nitrogen and sulfur (and others, if desired), and is designed to minimize the need for artificial geometry. Testing on a set of 187 diverse protein-ligand complexes has shown that the method is successful in predicting the location and approximate volume of the binding site in 73% of cases. Additional testing was performed on a set of 96 protein-ligand complexes with crystallographic structures of apo and holo forms, and AutoLigand was able to predict the binding site in 80% of the apo structures.     

The environment of amide groups in protein–ligand complexes: H-bonds and beyond

A comprehensive structural analysis of interactions involving amide NH and C=O groups in protein-ligand complexes has been performed based on 3,275 published crystal structures (resolution < or =2.5 A). Most of the amide C=O and NH groups at the protein-ligand interface are highly buried within the binding site and involved in H-bonds with corresponding counter-groups. Small percentages of C=O and NH groups are solvated or embedded in hydrophobic environments. In particular, C=O groups show a higher propensity to be solvated or embedded in a hydrophobic environment than NH groups do. A small percentage of carbonyl groups is involved in weak hydrogen bonds with CH. Cases of dipolar interactions, involving carbonyl oxygen and electrophilic carbon atoms, such as amide, amidinium, guanidium groups, are also identified. A higher percentage of NH are in contact with aromatic carbons, interacting either through hydrogen bonds (preferably with the NH group pointing towards a ring carbon atom) or through stacking between amide plane and ring plane. Comprehensive studies such as the present one are thought to be important for future improvements in the molecular design area, in particular for the development of new scoring functions. [Figure: see text].     

Synthesis, characterization and molecular structure of 1,1′ bis(diphenylphosphino ferrocene)dioxide complex of uranyl nitrate

A 1,1′ bis(diphenylphosphino ferrocene)dioxide complex of uranyl nitrate was synthesized and characterized by IR, ¹H and ³¹P{¹H} NMR spectroscopic and X-ray diffraction methods. The structure of the compound shows that the uranium atom is surrounded by eight oxygen atoms in a hexagonal bi-pyramidal geometry. Two oxygen atoms from 1,1′ bis(diphenylphosphino ferrocene)dioxide ligand and four oxygen atoms from the nitrate groups form a planar hexagon. The two uranyl oxygen atoms occupy the axial position. The 1,1′ bis(diphenylphosphino ferrocene)dioxide ligand acts as a bidentate chelating ligand with a bite angle of 71.56(8)° around the uranium(VI) atom, which is much smaller in value compare to any of the previously reported values (90.1°-154.0°) for this ligand.     

Synthesis and Characterization of a Ditriflate-Bridged, Diiron(II) Complex with Syn-N-Donor Ligands: [Fe(2)(μ-OTf)(2)(PIC(2)DET)(2)](BARF)(2)

The synthesis and characterization of the diiron(II) complex [Fe(2)(μ-OTf)(2)-(PIC(2)DET)(2)](BARF)(2) (2), where PIC(2)DET is a 2,3-diethynyltriptycene-linked dipicolinic methyl ester ligand, are described. The dication in 2, contains, [Fe(2)(μ-OTf)(2)(PIC(2)DET)(2)](2+) two symmetry-equivalent iron atoms with octahedral coordination geometries. Each metal ion has a N(2)O(4) atom donor set that includes four atoms from two picolinic ester N,O chelate rings, as well as two oxygen atoms from the bridging trifluoromethanesulfonate groups. The Fe(2)(μ-OTf)(2) core of 2 is stabilized by two PIC(2)DET ligands that bind the two metal ions in a head-to-head fashion, leading to an Fe···Fe distance of 5.173(1)?. Molar conductivity data for 2 are consistent with Fe(2)(μ-OTf)(2)(PIC(2)DET)(2)](2+) retaining its identity in acetone solutions, where it behaves as a 2:1 electrolyte. (1)H NMR spectroscopic, solution (d(6)-acetone) and solid-state magnetic susceptibility data all indicate that the iron atoms of 2 are high-spin (S = 2). A fit of the magnetic data (2 - 300K) to a spin-only isotropic exchange Hamiltonian H = -2JS(1)·S(2) are consistent with weak antiferromagnetic coupling between the two iron atoms with J ~ -0.99(2) cm(-1) and g = 2.10(1).     

A consensus-binding structure for adenine at the atomic level permits searching for the ligand site in a wide spectrum of adenine-containing complexes

Attempts to derive structural features of ligand-binding sites have traditionally involved seeking commonalities at the residue level. Recently, structural studies have turned to atomic interactions of small molecular fragments to extract common binding-site properties. Here, we explore the use of larger ligand elements to derive a consensus binding structure for the ligand as a whole. We superimposed multiple molecular structures from a nonredundant set of adenosine-5'-triphosphate (ATP) protein complexes, using the adenine moiety as template. Clustered binding-site atoms of compatible atomic classes forming attractive contacts with the adenine probe were extracted. A set of atomic clusters characterizing the adenine binding pocket was then derived. Among the clusters are three vertices representing the interactions of adenine atom N6 with its protein-binding niche. These vertices, together with atom C6 of the purine ring system, complete the set of four vertices for the pyramid-like structure of the N6 anchor atom. Also, the sequence relationship for the adenine-binding loop interacting with the C2-N6 end of the conjugated ring system is expanded to include a third hydrophilic cluster interacting with atom N1. A search procedure involving interatomic distances between cluster centers was formulated and applied to seek putative binding sites in test cases. The results show that a consensus network of clusters, based on an adenine probe and an ATP-complexed training set of proteins, is sufficient to recognize the experimental cavity for adenine in a wide spectrum of ligand-protein complexes.     

Active site modeling in copper azurin molecular dynamics simulations

Active site modeling in molecular dynamics simulations is investigated for the reduced state of copper azurin. Five simulation runs (5 ns each) were performed at room temperature to study the consequences of a mixed electrostatic/constrained modeling for the coordination between the metal and the polypeptide chain, using for the ligand residues a set of charges that is modified with respect to the apo form of the protein by the presence of the copper ion.The results show that the different charge values do not lead to relevant effects on the geometry of the active site of the protein, as long as bond distance constraints are used for all the five ligand atoms. The distance constraint on the O atom of Gly45 can be removed without altering the active site geometry. The coordination between Cu and the other axial ligand Met121 is outlined as being flexible. Differences are found between the bonds of the copper ion with the two apparently equivalent N¹ atoms of His46 and His117.The overall findings are discussed in connection with the issue of determining a model for the active site of azurin suitable to be used in molecular dynamics simulations under unfolding conditions. Figure Model of azurin active site. Copper ligand residues are cut off at C position except Gly45, for which the portion of backbone connecting it to His46 is shown. Only polar H atoms are shown. All atoms are in standard colors (Cu in violet), and the five ligands are labeled     

Synthesis, characterisation, X-ray structure and biological activity of three new 5-formyluracil thiosemicarbazone complexes

Three new complexes of transition metals as copper, nickel and cobalt with 5-formyluracil thiosemicarbazone (H3ut) have been synthesised and characterised by single-crystal X-ray diffraction. In all compounds the ligand behaves as SNO terdentate. In the copper complex the coordination geometry is square pyramidal with the ligand lying on the basal plane and two water molecules that complete the metal environment, the nickel compound is surrounded by six donor atoms (three of the ligand, two water oxygen atoms and a chlorine atom) in an octahedral fashion, and cobalt also shows an octahedral geometry but determined only by two terdentate ligand molecules. These three compounds have been tested on human leukemic cell lines K562 and CEM. The nickel and cobalt complexes have demonstrated low activity in cell growth, while the copper complex that is more active has been tested also on a third leukemic human cell line (U937), but it was not able to induce apoptosis on all cell lines.     

Structural features of the Cu(2+)-vancomycin complex

The structure of vancomycin coordinated to Cu(2+) ions is presented and structural aspects upon metal coordination are discussed. The asymmetric part of unit cell comprises two independent molecules of vancomycin-Cu(2+) complex, one of them is partially disordered. The binding site involves one imino nitrogen atom, two amide nitrogen atoms delivered by peptide bonds, and carboxyl oxygen from the peptide moiety. The identical set of donor atoms is not reflected in identical coordination geometry around individual metal ions. The studied complex presents two distinct types of conformation. Additionally, leucinyl side chain in one conformer is disordered leading to another type of conformation. The complex molecules form heterodimer with antiparallel hydrogen bonding.     

Novel lipophilic amidate oxorhenium and oxotechnetium complexes as potential brain agents: synthesis, characterization and biological evaluation

Novel oxorhenium and oxotechnetium complexes based on the tetradentate 1-(2-hydroxybenzamido)-2-(pyridinecarboxamido)benzene, H3L, ligand have been synthesized and characterized herein. Thus, by reacting equimolar quantities of the triply deprotonated ligand L3- with the suitable MO3+ precursor, the following neutral MOL complexes could be easily produced following similar synthetic routes: M = Re (1), M = 99gTc (2), and M = 99mTc (3). Complexes 1 and 2, prepared in macroscopic amounts, were chemically characterized and their structure determined by single-crystal X-ray analysis. They are isostructural metal chelates, adopting a distorted square pyramidal geometry around the metal. The N3O donor atom set of the tetradentate ligand defines the basal plane and the oxygen atom of the M = O core occupies the apex of the pyramid. Complex 3 forms quantitatively at tracer level by mixing the H3L ligand with Na99mTcO4 generator eluate in aqueous alkaline media and using tin chloride as reductant in the presence of citrate. Its structure was established by chromatographic comparison with prototypic complexes 1 and 2 using high-performance liquid chromatographic techniques. When challenged with excess glutathione in vitro, complex 3 is rapidly converted to hydrophilic unidentified metal species. Tissue distribution data after administration of complex 3 in vivo revealed a significant uptake and retention of this compound in brain tissue.     

Characterization of dioxygenated cobalt(II)-carnosine complexes by Raman and IR spectroscopy

Raman and IR studies are carried out on carnosine (beta-alanyl-L-histidine, Carnos) and its complexes with cobalt(II) at different metal/ligand ratios and basic pH. Binuclear complexes that bind molecular oxygen are formed and information regarding the O-O bridge is obtained from the Raman spectra. When the Co(II)/Carnos ratio is 相似文献   

Modeling of metal interaction geometries for protein-ligand docking

The accurate modeling of metal coordination geometries plays an important role for structure-based drug design applied to metalloenzymes. For the development of a new metal interaction model, we perform a statistical analysis of metal interaction geometries that are relevant to protein-ligand complexes. A total of 43,061 metal sites of the Protein Data Bank (PDB), containing amongst others magnesium, calcium, zinc, iron, manganese, copper, cadmium, cobalt, and nickel, were evaluated according to their metal coordination geometry. Based on statistical analysis, we derived a model for the automatic calculation and definition of metal interaction geometries for the purpose of molecular docking analyses. It includes the identification of the metal-coordinating ligands, the calculation of the coordination geometry and the superposition of ideal polyhedra to identify the optimal positions for free coordination sites. The new interaction model was integrated in the docking software FlexX and evaluated on a data set of 103 metalloprotein-ligand complexes, which were extracted from the PDB. In a first step, the quality of the automatic calculation of the metal coordination geometry was analyzed. In 74% of the cases, the correct prediction of the coordination geometry could be determined on the basis of the protein structure alone. Secondly, the new metal interaction model was tested in terms of predicting protein-ligand complexes. In the majority of test cases, the new interaction model resulted in an improved docking accuracy of the top ranking placements.     

A 1H nuclear relaxation study of the Mn2+ . bleomycin complex. Proximity of the metal to the DNA-binding site

The antineoplastic action of bleomycin is thought to involve the aerobic degradation of DNA by the Fe2+ . bleomycin complex. Different parts of the bleomycin molecule have been implicated in metal binding and DNA binding. To probe the structure of a metal-containing bleomycin, we studied the effects of the high spin Mn2+ ion in the Mn2+ . bleomycin complex on the longitudinal nuclear relaxation rates of various protons in the molecule. Complexation of Mn2+ to bleomycin was also studied by EPR, and a Scatchard plot of the EPR data revealed a single tight divalent cation-binding site per molecule. From the magnitudes of the paramagnetic effects of Mn2+ on the nuclear relaxation rates of several assigned resonances, we calculate the relative distances of the corresponding protons from the metal. Using a pyrimidine methyl to metal distance of 6.5 A, consistent with the metal coordination of this aromatic group of bleomycin established on the basis of other studies, we find from our data that the bithiazole and COOH-terminal portions of the molecule are located spatially very close to the metal. These groups have previously been implicated in DNA binding. Our metal to bithiazole proton distances (approximately 5.4 A) are consistent with bithiazole as a metal ligand, although possible involvement of interactions other than direct coordination in maintaining close proximity cannot be excluded. Our distance data also argue against the imidazole ring of beta-hydroxyhistidine as a ligand. The short distance between the metal- and DNA-binding sites indicated by our studies would help ensure that the reactive reduced oxygen radicals produced at the metal site during Fe2+ oxidation in the aerobic Fe2+ . bleomycin complex reach the substrate DNA before the destruction of these radicals can occur in other ways.     

Estimating protein-ligand binding free energy: atomic solvation parameters for partition coefficient and solvation free energy calculation

Solvation energy calculation is one of the main difficulties for the estimation of protein-ligand binding free energy and the correct scoring in docking studies. We have developed a new solvation energy estimation method for protein-ligand binding based on atomic solvation parameter (ASP), which has been shown to improve the power of protein-ligand binding free energy predictions. The ASP set, designed to handle both proteins and organic compounds and derived from experimental n-octanol/water partition coefficient (log P) data, contains 100 atom types (united model that treats hydrogen atoms implicitly) or 119 atom types (all-atom model that treats hydrogen atoms explicitly). By using this unified ASP set, an algorithm was developed for solvation energy calculation and was further integrated into a score function for predicting protein-ligand binding affinity. The score function reproduced the absolute binding free energies of a test set of 50 protein-ligand complexes with a standard error of 8.31 kJ/mol. As a byproduct, a conformation-dependent log P calculation algorithm named ASPLOGP was also implemented. The predictive results of ASPLOGP for a test set of 138 compounds were r = 0.968, s = 0.344 for the all-atom model and r = 0.962, s = 0.367 for the united model, which were better than previous conformation-dependent approaches and comparable to fragmental and atom-based methods. ASPLOGP also gave good predictive results for small peptides. The score function based on the ASP model can be applied widely in protein-ligand interaction studies and structure-based drug design.     

The structure of the zinc sites of Escherichia coli DNA-dependent RNA polymerase.

X-ray absorption spectroscopy is ideally suited for the investigation of the electronic structure and the local environment (approximately 5 A) of specific atoms in biomolecules. While the edge region provides information about the valence state of the absorbing atom, the chemical identity of neighboring atoms, and the coordination geometry, the extended x-ray absorption fine structure region contains information about the number and average distance of neighboring atoms and their relative disorder. The development of sensitive detection methods has allowed studies using near physiological concentrations (as low as approximately 100 microM). RNA polymerase from Escherichia coli contains two zinc atoms: one tightly bound in the beta' subunit, the subunit that participates in template binding, and the other loosely bound in the beta subunit, the subunit that participates in substrate binding. X-ray absorption studies of these zinc sites in the native protein and of the zinc site in the beta' subunit after removal of the zinc in the beta subunit site by p-(hydroxymercuri)benzenesulfonate (Giedroc, D. P., and Coleman, J. E. (1986) Biochemistry 25, 4969-4978) indicate that both zinc sites have octahedral coordination. The zinc in the beta' subunit site has four sulfur ligands at an average distance of 2.36 +/- 0.02 A and two oxygen (or nitrogen) ligands at an average distance of 2.23 +/- 0.02 A. The beta subunit zinc site has five sulfur ligands at an average distance of 2.38 +/- 0.01 A and one histidine nitrogen ligand at 2.14 +/- 0.02 A. These results are in general agreement with earlier biochemical and spectroscopic studies.     

Synthesis, crystal structure and electroluminescent properties of a Schiff base zinc complex

A new complex of zinc with a Schiff base, zinc(N,N′-bis(salicylidene)-3, 6-dioxa-1, 8-diaminooctane monohydrate) (ZnBSO · H₂O), was synthesized and characterized by means of elemental analyses, IR spectra and DTA-TG. Its structure was determined by X-ray single crystal analysis. It was demonstrated that the zinc atom is coordinated by the two oxygen atoms in phenolate and two nitrogen atoms in imine of the ligand in a slightly distorted tetrahedral geometry, while the two oxygen atoms from the oxa-alkyl chain are not coordinated to Zn(II) atom. The energy levels of the HOMO, LUMO and the electrochemical band gap were determined by cyclic voltammeter. The electroluminescent devices with the complex as the emitter showed bright blue emission with a peak at 450 nm, which is same as the fluorescence of the complex in both solution and solid states.     

Direct determination of the hydrogen bonding arrangement in anhydrous β-chitin by neutron fiber diffraction

The hydrogen bonding arrangement in anhydrous β-chitin, a homopolymer of N-acetylglucosamine, was directly determined by neutron fiber diffraction. Data were collected from a sample prepared from the bathophilous tubeworm Lamellibrachia satsuma in which all labile hydrogen atoms had been replaced by deuterium. Initial positions of deuterium atoms on hydroxyl and acetamide groups were directly located in Fourier maps synthesized using phases calculated from the X-ray structure and amplitudes measured from the neutron data. The hydrogen bond arrangement in the refined structure is in general agreement with predictions based on the X-ray structure: O3 donates a hydrogen bond to the O5 ring oxygen atom of a neighboring residue in the same chain; N2 and O6 donate hydrogen bonds to the same carbonyl oxygen O7 of an adjacent chain. The intramolecular O3···O5 hydrogen bond has the most energetically favorable geometry with a hydrogen to acceptor distance of 1.77 ? and a hydrogen bond angle of 171°.