BSP-SLIM: a blind low-resolution ligand-protein docking approach using predicted protein structures

We developed BSP‐SLIM, a new method for ligand–protein blind docking using low‐resolution protein structures. For a given sequence, protein structures are first predicted by I‐TASSER; putative ligand binding sites are transferred from holo‐template structures which are analogous to the I‐TASSER models; ligand–protein docking conformations are then constructed by shape and chemical match of ligand with the negative image of binding pockets. BSP‐SLIM was tested on 71 ligand–protein complexes from the Astex diverse set where the protein structures were predicted by I‐TASSER with an average RMSD 2.92 Å on the binding residues. Using I‐TASSER models, the median ligand RMSD of BSP‐SLIM docking is 3.99 Å which is 5.94 Å lower than that by AutoDock; the median binding‐site error by BSP‐SLIM is 1.77 Å which is 6.23 Å lower than that by AutoDock and 3.43 Å lower than that by LIGSITE^CSC. Compared to the models using crystal protein structures, the median ligand RMSD by BSP‐SLIM using I‐TASSER models increases by 0.87 Å, while that by AutoDock increases by 8.41 Å; the median binding‐site error by BSP‐SLIM increase by 0.69Å while that by AutoDock and LIGSITE^CSC increases by 7.31 Å and 1.41 Å, respectively. As case studies, BSP‐SLIM was used in virtual screening for six target proteins, which prioritized actives of 25% and 50% in the top 9.2% and 17% of the library on average, respectively. These results demonstrate the usefulness of the template‐based coarse‐grained algorithms in the low‐resolution ligand–protein docking and drug‐screening. An on‐line BSP‐SLIM server is freely available at http://zhanglab.ccmb.med.umich.edu/BSP‐SLIM . Proteins 2012. © 2011 Wiley Periodicals, Inc.     

EXAFS and Mössbauer characterization of the Diiron(III) site in stearoyl-acyl carrier protein Δ9– desaturase

Stearoyl-acyl carrier protein (ACP) Δ⁹-desaturase (Δ9D) from the castor plant is the best characterized soluble acyl-ACP desaturase. This enzyme utilizes a diiron center to catalyze the O₂- and NADPH-dependent introduction of a cis double bond between carbons 9 and 10 of stearoyl-ACP, yielding oleoyl-ACP. In the present study, we have used X-ray absorption spectroscopy to provide the first metrical information for the diferric oxidation state. These studies reveal distinct diiron clusters that have Fe-Fe distances of either 3.12 or 3.41?Å. The species having the 3.12?Å Fe-Fe distance also exhibits a 1.8?Å Fe-O bond and is thus proposed to represent Δ9D molecules containing a (μ-oxo)bis(μ-carboxylato)diiron(III) cluster. The species having the 3.41?Å Fe-Fe distance exhibits no short Fe-O bond, and thus likely represents Δ9D molecules containing a (μ-hydroxo)diiron(III) cluster. Mössbauer studies of the extended X-ray absorption fine structure (EXAFS) samples revealed three quadrupole doublets (ΔE _Q(1)=1.53?mm/s, 72%;ΔE _Q(2)=0.72?mm/s, 21%;ΔE _Q(3)=2.20?mm/s, 7%) that originate from three distinct dinuclear clusters. From analysis of spectral intensities and by comparison with previous studies of (μ-oxo)- and (μ-hydroxo)diiron(III) clusters in both model complexes and proteins, doublet 1, the Mössbauer majority species, is likely associated with the EXAFS majority species having a 3.12?Å Fe-Fe separation and a 1.8?Å Fe-μ-oxo bond, while doublet 2 likely results from one iron site (or both) of a cluster associated with the EXAFS species having a 3.41?Å Fe-Fe separation. The presence of multiple diiron center conformations in diferric Δ9D may reflect the necessity for the active site to allow access of the substrate stearoyl-ACP (～9?kDa) during desaturation catalysis.     

Preliminary X-ray crystallographic analysis of Tritrichomonas foetus inosine-5′-monophosphate dehydrogenase

Inosine-5′-monophosphate dehydrogenase (IMPDH) from the protozoan parasite Tritrichomonas foetus has been expressed in E. coli and crystallized. Crystals were grown to 0.1 mm in each dimension in 18 to 72 h using ammonium sulfate and low-molecular-weight polyethylene glycols. The crystals belong to the cubic space group P432 with unit cell edge = 157.25 Å. The enzyme is a homotetramer with each monomer having a molecular weight of 55,534 Da. There is one monomer per asymmetric unit, based on a volume/mass ratio of 2.7 ų/Da and self-rotation analysis. The crystals are adequately stable to allow a complete data set to be collected from a single crystal. Complete native data sets have been collected to 2.3 Å resolution at 4°C using synchrotron radiation. High-quality complete data extending to 3.0 Å resolution have been collected from crystals of four putative derivatives, and the data appear to be isomorphous with that of the native crystals in each case. Efforts to solve the derivatives for use in MIR phasing are underway. © 1995 Wiley-Liss, Inc.     

An X-ray diffraction study of poly-L-arginine hydrochloride

An x-ray study has been made of polyarginine hydrochloride to investigate whether, like polylysine hydrochloride, it can undergo conformational changes merely from variations in the degree of hydration. X-ray powder and fiber photographs of specimens containing up to about five molecules of water per arginine residue show features characteristic of α-helical structures including a 5.4-Å layer line and a meridional 1.5-Å reflection. Increasing the water content from ¹/₂ to 6¹/₂ molecules per residue causes the a axis of the hexagonal unit cell to increase from 14.4 Å to 15.8 Å, with no appreciable change in the 27.0 Å c axis. Removal of the last half molecule of water results in a very diffuse α pattern, but on rehydration the sharp pattern reappears. Specimens containing five to twenty water molecules per residue show quite a different pattern, the spacing of which do not vary appreciably with hydration. This pattern includes a meridional 3.4-Å reflection, a feature commonly shown by β structures, and indeed all the reflections can be satisfactorily indexed in terms of a monoclinic unit cell with a = 9.26 Å, b = 22.05 Å, c = 6.76 Å, and γ = 108.9°. These dimensions are shown by models to be compatible with a β pleated-sheet structure.     

Crystallization,molecular replacement solution,and refinement of tetrameric β-amylase from sweet potato

Sweet potato β-amylase is a tetramer of identical subunits, which are arranged to exhibit 222 molecular symmetry. Its subunit consists of 498 amino acid residues (Mr 55,880). It has been crystallized at room temperature using polyethylene glycol 1500 as precipitant. The crystals, growing to dimensions of 0.4 mm × 0.4 mm × 1.0 mm within 2 weeks, belong to the tetragonal space group P4₂2₁2 with unit cell dimensions of a = b = 129.63 Å and c = 68.42 Å. The asymmetric unit contains 1 subunit of β-amylase, with a crystal volume per protein mass (V_M) of 2.57 ų/Da and a solvent content of 52% by volume. The three-dimensional structure of the tetrameric β-amylase from sweet potato has been determined by molecular replacement methods using the monomeric structure of soybean enzyme as the starting model. The refined subunit model contains 3,863 nonhydrogen protein atoms (488 amino acid residues) and 319 water oxygen atoms. The current R-value is 20.3% for data in the resolution range of 8–2.3 Å (with 2 σ cut-off) with good stereochemistry. The subunit structure of sweet potato β-amylase (crystallized in the absence of α-cyclodextrin) is very similar to that of soybean β-amylase (complexed with α-cyclodextrin). The root-mean-square (RMS) difference for 487 equivalent Cα atoms of the two β-amylases is 0.96 Å. Each subunit of sweet potato β-amylase is composed of a large (α/β)₈ core domain, a small one made up of three long loops [L3 (residues 91–150), LA (residues 183–258), and L5 (residues 300–327)], and a long C-terminal loop formed by residues 445–493. Conserved Glu 187, believed to play an important role in catalysis, is located at the cleft between the (α/β)₈ barrel core and a small domain made up of three long loops (L3, L4, and L5). Conserved Cys 96, important in the inactivation of enzyme activity by sulfhydryl reagents, is located at the entrance of the (α/β)₈ barrel. © 1995 Wiley-Liss, Inc.     

Crystallization and preliminary crystallographic analysis of an amylopullulanase from the hyperthermophilic archaeon Pyrococcus woesei

The thermostable amylopullulanase from Pyrococcus woesei was crystallized. Crystals, suitable for a crystallographic analysis up to a size of 0.6 mm in their longest dimension, have been obtained by the vapor diffusion method in a solution containing polyethyleneglycol 4000 (PEG 4000), isopropanol, and Tris/Cl^? buffer pH 7.5. Crystals grown under these conditions form hexagonal rods and diffract to a maximum resolution of 3 Å. The crystals belong to the trigonal lattice type with the spacegroup P3₁21 or P3₂21, respectively, have the cell dimensions a = b = 96.8 Å, c = 196.2 Å, α = β = 90°,γ = 120°. The crystals have a theoretical packing density of 2.7 ų/Da, assuming one molecule with a molecular weight of 88.8 kDa in the asymmetric unit. Furthermore the self-rotation analysis of the dataset revealed only crystallographic symmetries. The merged native data of two crystals resulted in a 88% complete dataset. © 1995 Wiley-Liss, Inc.     

Structure of GlyGly peptide in the crystalline state as studied by X-ray diffraction and solid state 13C-NMR methods

A new type of crystal of glycylglycine (GlyGly) hydrate was crystallized from an aqueous solution, and the structure of the crystal has been determined by x-ray diffraction. The crystal is monoclinic, and the space group is C2/c, with the cell constants of a = 15.941(2) Å, b = 4.774(2) Å, c = 19.428(2) Å and β = 109.884(7)° at 296 K. There are eight GlyGly molecules and six water solvent in the cell. The GlyGly molecules are packed in a parallel β-sheet arrangement. The single crystal was obtained with a maximum size of 10 × 7 × 4 mm and is not stable under atmospheric conditions. The transparent crystal turned to turbid with the elapse of time. The isotropic ¹³C chemical shifts obtained from the ¹³C cross polarization magic angle spinning nmr experiments reveal that GlyGly hydrate was changed into GlyGly (form α) by dehydration. © 1998 John Wiley & Sons, Inc. Biopoly 45: 333–339, 1998     

Three-dimensional image reconstruction of straight flagella from a mutant Salmonella typhimurium.

The structure of the straight flagella from a mutant Salmonella typhimurium was studied by electron microscopy using digital image processing, including three-dimensional reconstruction, to an effective resolution of about 14 Å.Three-dimensional studies suggest that there are two sets of intersubunit bonds, i.e. intraprotofilament bonds along the (n = 11, l = 1) helix at a radius of about 55 Å and interprotofilament bonds along the (n = ?5, l = 7) helix at radii of about 10 to 15 Å and 50 Å, and along the (n = 6, l = 8) helix at a radius of about 45 Å and along the (n = 1, l = 15) helix at a radius of about 20 Å. There are four high density regions in a morphological subunit. These regions are situated at radii of about 15 Å, 40 Å, 70 Å and 80 Å. Variation was seen in the position of the high density regions at radii of about 15 Å and 40 Å among the ten models that were reconstituted individually. The regions at radii of 40 Å and 70 Å are the highest in density. The radial distance between these two regions is consistent with the 32 Å feature of a cylindrically averaged Patterson function calculated using equatorial data from X-ray diffraction pattern (Champness, 1968,1971).At the outer radii of the flagellum the shape of the morphological subunit roughly corresponded to that of the “chevron” described by O'Brien &; Bennett (1972), but there was no corresponding structure at the inner radii; the appearance of chevrons in that region might arise from the superposition of the two sides of the helical lattice.The biological significance of the “beaded” submolecular structure of flagellin and the presence of two sets of intersubunit bonds at the different radii is discussed with reference to the waveform and polymorphic behaviour of flagellar filaments.     

The utility of geometrical and chemical restraint information extracted from predicted ligand-binding sites in protein structure refinement

Exhaustive exploration of molecular interactions at the level of complete proteomes requires efficient and reliable computational approaches to protein function inference. Ligand docking and ranking techniques show considerable promise in their ability to quantify the interactions between proteins and small molecules. Despite the advances in the development of docking approaches and scoring functions, the genome-wide application of many ligand docking/screening algorithms is limited by the quality of the binding sites in theoretical receptor models constructed by protein structure prediction. In this study, we describe a new template-based method for the local refinement of ligand-binding regions in protein models using remotely related templates identified by threading. We designed a Support Vector Regression (SVR) model that selects correct binding site geometries in a large ensemble of multiple receptor conformations. The SVR model employs several scoring functions that impose geometrical restraints on the Cα positions, account for the specific chemical environment within a binding site and optimize the interactions with putative ligands. The SVR score is well correlated with the RMSD from the native structure; in 47% (70%) of the cases, the Pearson’s correlation coefficient is >0.5 (>0.3). When applied to weakly homologous models, the average heavy atom, local RMSD from the native structure of the top-ranked (best of top five) binding site geometries is 3.1 Å (2.9 Å) for roughly half of the targets; this represents a 0.1 (0.3) Å average improvement over the original predicted structure. Focusing on the subset of strongly conserved residues, the average heavy atom RMSD is 2.6 Å (2.3 Å). Furthermore, we estimate the upper bound of template-based binding site refinement using only weakly related proteins to be ～2.6 Å RMSD. This value also corresponds to the plasticity of the ligand-binding regions in distant homologues. The Binding Site Refinement (BSR) approach is available to the scientific community as a web server that can be accessed at http://cssb.biology.gatech.edu/bsr/.     

Crystallization and preliminary X-ray investigation of the recombinant Trypanosoma brucei rhodesiense calmodulin

Bipyramidal crystals of the recombinant calmodulin from Trypanosoma brucei rhodesiense were obtained by vapor diffusion against 55% (v/v) 2-methyl-2,4-pentanediol in 0.05 M cacodylate buffer, pH 5.6. When few nucleation events occurred, crystals grew to 0.25 × 0.25 × 1.20 mm. The space group of the crystal is I4₁22, with unit cell dimensions a = b = 56.88 Å, c = 230.11 Å, α = β = γ = 90°, z = 16. The molecular mass and volume of the unit cell suggest that there is one molecule in the asymmetric unit. The I/σ(I) ratio for data at 3.0 Å resolution was 3.67, indicating that the final structure can be refined at higher resolution. Molecular replacement methods and the PC-refinement technique have not yet yielded the structure under a variety of search conditions. We are currently investigating the multiple isomorphous replacement approach to determine this crystal structure. © 1995 Wiley-Liss, Inc.     

Preliminary Crystallographic Analysis of a Proteolytically Modified Form of E. coli Single Stranded DNA Binding Protein

Abstract

A proteolytically modified form of the Escherichia coli single-stranded DNA-Binding protein (SSB) has been crystallized from 15% saturated sodium citrate. Crystals as large as 1.0mm × 0.3mm × 0.2mm were obtained and these diffract beyond 3Å resolution. X-ray photographic analysis demonstrated a rhombohedral unit cell of space group R3 with an equivalent triple centered hexagonal unit cell having dimensions of a = b = 62.9 A and c = 264.3A. These crystals were judged to be adequate for a three dimensional structure determination.     

Thermodynamic resolution: How do errors in modeled protein structures affect binding affinity predictions?

The present study addresses the effect of structural distortion, caused by protein modeling errors, on calculated binding affinities toward small molecules. The binding affinities to a total of 300 distorted structures based on five different protein–ligand complexes were evaluated to establish a broadly applicable relationship between errors in protein structure and errors in calculated binding affinities. Relatively accurate protein models (less than 2 Å RMSD within the binding site) demonstrate a 14.78 (±7.5)% deviation in binding affinity from that calculated by using the corresponding crystal structure. For structures of 2–3 Å, 3–4 Å, and >4 Å RMSD within the binding site, the error in calculated binding affinity increases to 20.8 (±5.98), 22.79 (±11.3), and 29.43 (±11.47)%, respectively. The results described here may be used in combination with other tools to evaluate the utility of modeled protein structures for drug development or other ligand‐binding studies. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

In silico predictions of 3D structures of linear and cyclic peptides with natural and non‐proteinogenic residues

We extended the use of Peplook, an in silico procedure for the prediction of three‐dimensional (3D) models of linear peptides to the prediction of 3D models of cyclic peptides and thanks to the ab initio calculation procedure, to the calculation of peptides with non‐proteinogenic amino acids. Indeed, such peptides cannot be predicted by homology or threading. We compare the calculated models with NMR and X‐ray models and for the cyclic peptides, with models predicted by other in silico procedures (Pep‐Fold and I‐Tasser). For cyclic peptides, on a set of 38 peptides, average root mean square deviation of backbone atoms (BB‐RMSD) was 3.8 and 4.1 Å for Peplook and Pep‐Fold, respectively. The best results are obtained with I‐Tasser (2.5 Å) although evaluations were biased by the fact that the resolved Protein Data Bank models could be used as template by the server. Peplook and Pep‐Fold give similar results, better for short (up to 20 residues) than for longer peptides. For peptides with non‐proteinogenic residues, performances of Peplook are sound with an average BB‐RMSD of 3.6 Å for ‘non‐natural peptides’ and 3.4 Å for peptides combining non‐proteinogenic residues and cyclic structure. These results open interesting possibilities for the design of peptidic drugs. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Crystallization and preliminary X-ray crystallographic analysis of chitinase from barley seeds

Chitinase from barley seeds has been crystallized at room temperature using polyethylene glycol as precipitant. The crystal is monoclinic, belonging to the space group P2₁, with unit cell parameters of a = 69.43 Å, b = 44.55 Å, c = 81.41 Å, and β = 111.95 Å. The asymmetric unit seems to contain two molecules of chitinase with a corresponding crystal volume per protein mass (V_M) of 2.25 ų/Da and a solvent content of 45% by volume. The crystal diffracts to at least 2.0 Å with X-rays from a rotating anode source and is very stable in the X-ray beam. X-ray data have been collected to better than 2.2 Å Bragg spacing from a native crystal. © 1993 Wiley-Liss, Inc.     

Preliminary X-ray diffraction analysis of crystals of turnip yellow mosaic virus (TYMV)

Turnip yellow mosaic virus (TYMV) was purified from Chinese cabbage and crystallized in a form that permits high resolution structural analysis using X-ray diffraction. The crystals have a hexagonal bipyramidal morphology and often achieve dimensions of 1.0 × 1.0 × 0.5 mm. The crystals appear to be of hexagonal space group P6₂22 with a = b = 525 Å, c=315 Å, but we cannot strictly rule out the possibility that the space group is P622. They appear different than any crystals of TYMV previously reported. There are three T = 3 virus particles in the unit cell, which implies that one quarter of the particle, or 45 protein subunits, comprises the asymmetric unit of the crystal. Native data have been collected using synchrotron radiation to a resolution of 3.2 Å. © 1995 Wiley-Liss, Inc.     

FINDSITE-metal: integrating evolutionary information and machine learning for structure-based metal-binding site prediction at the proteome level

The rapid accumulation of gene sequences, many of which are hypothetical proteins with unknown function, has stimulated the development of accurate computational tools for protein function prediction with evolution/structure‐based approaches showing considerable promise. In this article, we present FINDSITE‐metal, a new threading‐based method designed specifically to detect metal‐binding sites in modeled protein structures. Comprehensive benchmarks using different quality protein structures show that weakly homologous protein models provide sufficient structural information for quite accurate annotation by FINDSITE‐metal. Combining structure/evolutionary information with machine learning results in highly accurate metal‐binding annotations; for protein models constructed by TASSER, whose average Cα RMSD from the native structure is 8.9 Å, 59.5% (71.9%) of the best of top five predicted metal locations are within 4 Å (8 Å) from a bound metal in the crystal structure. For most of the targets, multiple metal‐binding sites are detected with the best predicted binding site at rank 1 and within the top two ranks in 65.6% and 83.1% of the cases, respectively. Furthermore, for iron, copper, zinc, calcium, and magnesium ions, the binding metal can be predicted with high, typically 70% to 90%, accuracy. FINDSITE‐metal also provides a set of confidence indexes that help assess the reliability of predictions. Finally, we describe the proteome‐wide application of FINDSITE‐metal that quantifies the metal‐binding complement of the human proteome. FINDSITE‐metal is freely available to the academic community at http://cssb.biology.gatech.edu/findsite‐metal/ . Proteins 2011. © 2010 Wiley‐Liss, Inc.     

Crystallization and preliminary X-ray diffraction study of the green flavoenzyme 5-Hydroxyvaleryl-CoA dehydratase/dehydrogenase from Clostridium aminovalericum

The bifunctional flavoenzyme 5-hydroxyvaleryl-CoA dehydratase/ dehydrogenase has been crystallized from solutions containing ammonium sulfate (form I) or polyethylene glycol (form II) as precipitant. In both cases, the crystals grew in the monoclinic space group C2. The unit cell dimensions for form I crystals were determined as a = 162.8 Å, b = 71.8 Å, c = 83.5 Å, β = 109.1°. Corresponding values for form II crystals were a = 161.2 Å, b = 71.6 Å, c = 82.2 Å, β = 109.3°. In both cases most probably there are two monomers per asymmetric unit. The crystals diffract to about 2 Å resolution and are rather stable in the X-ray beam. © 1994 Wiley-Liss, Inc.     

Water structure in a protein crystal: rubredoxin at 1.2 A resolution

The model for rubredoxin based on X-ray diffraction data has been extensively refined with a 1.2 Å resolution data set. Water oxygen atoms were deleted from the model if B exceeded 50 Ų and occupancy was less than 0.3 eÅ^?3. The final water model consists of 127 sites with B values ranging from 15 to $\text{6}\text{?}$0 Ų and occupancies from unity down to 0.3, the most tightly bound water oxygen atoms being hydrogen bonded to two or more main-chain nitrogen or oxygen atoms. The water forms extensive hydrogen bond networks bridging the crevices on the molecular surfaces and between adjacent molecules. The minimum distances of the water sites from the protein surface are distributed about two distinct maxima, the major one at 2.5 to 3 Å and a minor one at 4 to 4.5 Å. Beyond $\text{5}\text{?}$ to 6 Å from the protein surface, the discrete water merges into the aqueous continuum.     

An evaluation of implicit and explicit solvent model systems for the molecular dynamics simulation of bacteriophage T4 lysozyme

In this report we examine several solvent models for use in molecular dynamics simulations of protein molecules with the Discover program from Biosym Technologies. Our goal was to find a solvent system which strikes a reasonable balance among theoretical rigor, computational efficiency, and experimental reality. We chose phage T4 lysozyme as our model protein and analyzed 14 simulations using different solvent models. We tested both implicit and explicit solvent models using either a linear distance-dependent dielectric or a constant dielectric. Use of a linear distance-dependent dielectric with implicit solvent significantly diminished atomic fluctuations in the protein and kept the protein close to the starting crystal structure. In systems using a constant dielectric and explicit solvent, atomic fluctuations were much greater and the protein was able to sample a larger portion of conformational space. A series of nonbonded cutoff distances (9.0, 11.5, 15.0, 20.0 Å) using both abrupt and smooth truncation of the nonbonded cutoff distances were tested. The method of dual cutoffs was also tested. We found that a minimum nonbonded cutoff distance of 15.0 Å was needed in order to properly couple solvent and solute. Distances shorter than 15.0 Å resulted in a significant temperature gradient between the solvent and solute. In all trajectories using the proprietary Discover switching function, we found significant denaturation in the protein backbone; we were able to run successful trajectories only in those simulations that used no switching function. We were able to significantly reduce the computational burden by using dual cutoffs and still calculate a quality trajectory. In this method, we found that an outer cutoff distance of 15.0 Å and an inner cutoff distance of 11.5 worked well. While a 10 Å shell of explicit water yielded the best results, a 6 A shell of water yielded satisfactory results with nearly a 40% reduction in computational cost. © 1994 John Wiley & Sons, Inc.