Mouse-human comparative map resources on the Web

Comparative maps have been a valuable resource for extrapolating biological information among organisms. The relationship between mouse and human maps provides a framework for integrating information from each species and thereby increasing the utility of all available data such as gene location, structure and function. This review describes the various public resources, both databases and web sites, containing genome-wide mouse-human comparative map information available through the World-Wide Web. We will focus on the use and applicability of these resources in their current form and consider future potential directions.     

Efficient molecular docking of NMR structures: application to HIV-1 protease

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.     

Nomencurator: a nomenclatural history model to handle multiple taxonomic views

Evolutionary studies are generating increasing numbers of phylogenies which, in turn, sometimes result in changes to hierarchical organization and therefore changes in taxonomic nomenclature. A three-layered data model for a nomenclature database has been developed in order to elucidate the information structure in nomenclature and as a means to organize and manage a large, dynamic knowledge-base. In contrast to most other taxonomic databases, the model is publication-oriented rather than taxon-oriented and dynamic rather than static, in order to mimic the processes that taxonomists use naturally. The three-layered structure requires data integrity localized to each publication, instead of global data integrity, which relaxes constraints common to taxonomic databases and permits multiple taxonomic opinions: taxon names are made available as metadata within the model. Its prototype implementation, written in C ⁺⁺, has an autonomous self-identification mechanism to avoid spurious data-inflation in a publication-oriented data model. Self-identification is also desirable for distributed implementations of the nomenclature database. Publication-oriented design also will make maintenance easier than for taxon-oriented databases, much of the maintenance workload being amenable to automation. The three-layered data model was designed for use by taxonomists, but is also able to provide concise, reduced expression for non-experts required in biodiversity research, for example.     

Prediction of 5-HT3 receptor agonist-binding residues using homology modeling

5-HT(3) receptors demonstrate significant structural and functional homology to other members of the Cys-loop ligand-gated ion channel superfamily. The extracellular domains of these receptors share similar sequence homology (approximately 20%) with Limnaea acetylcholine binding protein, for which an x-ray crystal structure is available. We used this structure as a template for computer-based homology modeling of the 5-HT(3) receptor extracellular domain. AutoDock software was used to dock 5-HT into the putative 5-HT(3) receptor ligand-binding site, resulting in seven alternative energetically favorable models. Residues located no more than 5 A from the docked 5-HT were identified for each model; of these, 12 were found to be common to all seven models with five others present in only certain models. Some docking models reflected the cation-pi interaction previously demonstrated for W183, and data from these and other studies were used to define our preferred models.     

Extended X-ray absorption fine structure studies on the iron-containing subunit of ribonucleotide reductase from Escherichia coli

Iron K-edge X-ray absorption spectra were obtained on the protein B2, the small subunit of ribonucleotide reductase from Escherichia coli. Protein B2 contains a binuclear iron center with many properties in common with the iron center of oxidized hemerythrins. The extended X-ray absorption fine structure (EXAFS) measurements on protein B2 were analyzed and compared with published data for oxyhemerythrin. In protein B2 there are, in the first coordination shell around each Fe atom, five or six oxygen or nitrogen atoms that are directly coordinated ligands. In oxyhemerythrin there are six ligands to each iron. As in oxyhemerythrin, one of the ligands in the first shell of protein B2 is at a short distance, about 1.78 A, confirming the existence of a mu-oxo bridge. The other atoms of the first shell are at an average distance of 2.04 A, which is about 0.1 A shorter than in oxyhemerythrin. In protein B2 the Fe-Fe distance is in the range 3.26-3.48 A, and the bridging angle falls between 130 and 150 degrees. On the basis of these data, there is no direct evidence for any histidine ligands in protein B2, but the noise level leaves way for the possibility of a maximum of about three histidines for each Fe pair. The X-ray absorption spectrum of a hydroxyurea-treated sample was not significantly different from that of the native protein B2, which implies that no significant alteration in the structure of the iron site occurs upon destruction of the tyrosine radical.     

Two crystal forms of the restriction enzyme MspI-DNA complex show the same novel structure

The crystal structure of the Type IIP restriction endonuclease MspI bound to DNA containing its cognate recognition sequence has been determined in both monoclinic and orthorhombic space groups. Significantly, these two independent crystal forms present an identical structure of a novel monomer-DNA complex, suggesting a functional role for this novel enzyme-DNA complex. In both crystals, MspI interacts with the CCGG DNA recognition sequence as a monomer, using an asymmetric mode of recognition by two different structural motifs in a single polypeptide. In the crystallographic asymmetric unit, the two DNA molecules in the two MspI-DNA complexes appear to stack with each other forming an end-to-end pseudo-continuous 19-mer duplex. They are primarily B-form and no major bends or kinks are observed. For DNA recognition, most of the specific contacts between the enzyme and the DNA are preserved in the orthorhombic structure compared with the monoclinic structure. A cation is observed near the catalytic center in the monoclinic structure at a position homologous to the 74/45 metal site of EcoRV, and the orthorhombic structure also shows signs of this same cation. However, the coordination ligands of the metal are somewhat different from those of the 74/45 metal site of EcoRV. Combined with structural information from other solved structures of Type II restriction enzymes, the possible relationship between the structures of the enzymes and their cleavage behaviors is discussed.     

Secondary and tertiary structures of the transmembrane domains of the translocator protein TSPO determined by NMR. Stabilization of the TSPO tertiary fold upon ligand binding

Numerous biological functions are attributed to the peripheral-type benzodiazepine receptor (PBR) recently renamed translocator protein (TSPO). The best characterized function is the translocation of cholesterol from the outer to inner mitochondrial membrane, which is a rate-determining step in steroid biosynthesis. TSPO drug ligands have been shown to stimulate pregnenolone formation by inducing TSPO-mediated translocation of cholesterol. Until recently, no direct structural data on this membrane protein was available. In a previous paper, we showed that a part of the mouse TSPO (mTSPO) C-terminal region adopts a helical conformation, the side-chain distribution of which provides a groove able to fit a cholesterol molecule. We report here on the overall structural properties of mTSPO. This study was first undertaken by dissecting the protein sequence and studying the conformation of five peptides encompassing the five putative transmembrane domains from ¹H-NMR data. The secondary structure of the recombinant protein in micelles was then studied using CD spectroscopy. In parallel, the stability of its tertiary fold was probed using ¹H-¹⁵N NMR. This study provides the first experimental evidence for a five-helix fold of mTSPO and shows that the helical conformation of each transmembrane domain is mainly formed through local short-range interactions. Our data show that, in micelles, mTSPO exhibits helix content close to what is expected but an unstable tertiary fold. They reveal that the binding of a drug ligand that stimulates cholesterol translocation is able to stabilize the mTSPO tertiary structure.     

Molecular Dynamics Simulations and Structure-Guided Mutagenesis Provide Insight into the Architecture of the Catalytic Core of the Ectoine Hydroxylase

Many bacteria amass compatible solutes to fend-off the detrimental effects of high osmolarity on cellular physiology and water content. These solutes also function as stabilizers of macromolecules, a property for which they are referred to as chemical chaperones. The tetrahydropyrimidine ectoine is such a compatible solute and is widely synthesized by members of the Bacteria. Many ectoine producers also synthesize the stress protectant 5-hydroxyectoine from the precursor ectoine, a process that is catalyzed by the ectoine hydroxylase (EctD). The EctD enzyme is a member of the non-heme-containing iron(II) and 2-oxoglutarate-dependent dioxygenase superfamily. A crystal structure of the EctD protein from the moderate halophile Virgibacillus salexigens has previously been reported and revealed the coordination of the iron catalyst, but it lacked the substrate ectoine and the co-substrate 2-oxoglutarate. Here we used this crystal structure as a template to assess the likely positioning of the ectoine and 2-oxoglutarate ligands within the active site by structural comparison, molecular dynamics simulations, and site-directed mutagenesis. Collectively, these approaches suggest the positioning of the iron, ectoine, and 2-oxoglutarate ligands in close proximity to each other and with a spatial orientation that will allow the region-selective and stereo-specific hydroxylation of (4S)-ectoine to (4S,5S)-5-hydroxyectoine. Our study thus provides a view into the catalytic core of the ectoine hydroxylase and suggests an intricate network of interactions between the three ligands and evolutionarily highly conserved residues in members of the EctD protein family.     

Virtual screening pipeline and ligand modelling for H5N1 neuraminidase

The H5N1 virus neuraminidase structure was solved in two different conformations depending on the inhibitor concentration. In the absence of oseltamivir or at a low concentration, the neuraminidase structure assumes an open form that closes at a high oseltamivir concentration due to the shift of the so-called 150-loop near the active site. Although the close conformation is similar to all the other structurally known neuraminidase types, it doesn’t appear to be the most likely physiological condition for N1.To investigate the specific ligand binding properties of the open form, we screened by docking simulation, a large dataset of ligands and compared the results with closed form. The virtual screening procedure was implemented in a docking pipeline that also performs a step-by-step, target specific, filtering approach for data reduction. The selected ligands display binding ability involving multiple sites of interaction including the active site and an adjacent cavity made available by the 150-loop shift. Two ligands are especially interesting and are proposed as substituents to design oseltamivir derivatives specifically suited for the open conformation.     

Phosphate binding sites identification in protein structures

Nearly half of known protein structures interact with phosphate-containing ligands, such as nucleotides and other cofactors. Many methods have been developed for the identification of metal ions-binding sites and some for bigger ligands such as carbohydrates, but none is yet available for the prediction of phosphate-binding sites. Here we describe Pfinder, a method that predicts binding sites for phosphate groups, both in the form of ions or as parts of other non-peptide ligands, in proteins of known structure. Pfinder uses the Query3D local structural comparison algorithm to scan a protein structure for the presence of a number of structural motifs identified for their ability to bind the phosphate chemical group. Pfinder has been tested on a data set of 52 proteins for which both the apo and holo forms were available. We obtained at least one correct prediction in 63% of the holo structures and in 62% of the apo. The ability of Pfinder to recognize a phosphate-binding site in unbound protein structures makes it an ideal tool for functional annotation and for complementing docking and drug design methods. The Pfinder program is available at http://pdbfun.uniroma2.it/pfinder.     

Crystal structural analysis and metal-dependent stability and activity studies of the ColE7 endonuclease domain in complex with DNA/Zn2+ or inhibitor/Ni2+

The nuclease domain of ColE7 (N-ColE7) contains an H-N-H motif that folds in a beta beta alpha-metal topology. Here we report the crystal structures of a Zn2+-bound N-ColE7 (H545E mutant) in complex with a 12-bp duplex DNA and a Ni2+-bound N-ColE7 in complex with the inhibitor Im7 at a resolution of 2.5 A and 2.0 A, respectively. Metal-dependent cleavage assays showed that N-ColE7 cleaves double-stranded DNA with a single metal ion cofactor, Ni2+, Mg2+, Mn2+, and Zn2+. ColE7 purified from Escherichia coli contains an endogenous zinc ion that was not replaced by Mg2+ at concentrations of <25 mM, indicating that zinc is the physiologically relevant metal ion in N-ColE7 in host E. coli. In the crystal structure of N-ColE7/DNA complex, the zinc ion is directly coordinated to three histidines and the DNA scissile phosphate in a tetrahedral geometry. In contrast, Ni2+ is bound in N-ColE7 in two different modes, to four ligands (three histidines and one phosphate ion), or to five ligands with an additional water molecule. These data suggest that the divalent metal ion in the His-metal finger motif can be coordinated to six ligands, such as Mg2+ in I-PpoI, Serratia nuclease and Vvn, five ligands or four ligands, such as Ni2+ or Zn2+ in ColE7. Universally, the metal ion in the His-metal finger motif is bound to the DNA scissile phosphate and serves three roles during hydrolysis: polarization of the P-O bond for nucleophilic attack, stabilization of the phosphoanion transition state and stabilization of the cleaved product.     

The reaction of Octopus vulgaris hemocyanin with exogenous ligands: proposal of an allosteric model for the binding of cyanide and thiourea to the 11 S subunit

Octopus vulgaris hemocyanin in 11 S aggregation state binds oxygen following a noncooperative oxygen saturation curve with Hill coefficient n = 1. Under the same conditions the equilibrium and kinetics of the reaction with cyanide and other ligands are indicative of an anticooperative behavior displaying different characteristics for the different ligands. The data are consistent with an induced-fit type allosteric model which assumes for the 11 S subunit of O. vulgaris hemocyanin an annular structure made up by five identical domains each containing one binding site whose reactivity is near-neighbor regulated.     

The interplay of ligand binding and quaternary structure in the diverse interactions of dynein light chain LC8

Dynein light chain LC8 is a small, dimeric, and very highly conserved globular protein that is an integral part of the dynein and myosin molecular motors but appears to have a broader role in multiple protein complexes unrelated to molecular motors. LC8 binds to two families of targets: those having a KXTQT sequence fingerprint and those having a GIQVD fingerprint. All known LC8 binding partners containing these fingerprints share a common binding site on LC8 that raises the question of what determines binding specificity. Here, we present the crystal structure of apo-LC8 at 1.7-Å resolution, which, when compared with the crystal structures of several LC8 complexes, gives insight into the mechanism underlying the binding diversity of LC8. Peptide binding is associated with a shift in quaternary structure that expands the hydrophobic binding surface available to the ligand, in addition to changes in tertiary structure and ordering of LC8 around the binding groove. The observed quaternary shift suggests a mechanism by which binding at one of the two identical sites can influence binding at the other. NMR spectra of titrations with peptides from each fingerprint family show evidence of allosteric interaction between the two binding sites, to a differing degree in the two ligand families. Allosteric interaction between the binding sites may be a mechanism to promote simultaneous binding of ligands from the same family, providing a physiological role for the two fingerprints.     

Structural study of biologically significant ligands with major birch pollen allergen Betv1 by docking and molecular dynamics simulation

The major birch pollen allergen, Betv1 of Betula verrucosa is the main causative agent of birch pollen allergy in humans. Betv1 is capable of binding several physiological ligands including fatty acids, flavones, cytokinins and sterols. Until now, no structural information from crystallography or NMR is available regarding binding mode of any of these ligands into the binding pocket of Betv1. In the present study thirteen ligands have been successfully docked into the hydrophobic cavity of Betv1 and binding free energies of the complexes have been calculated using AutoDock 3.0.5. A linear relationship with correlation coefficient (R2) of 0.6 is obtained between ΔG(b)s values plotted against their corresponding IC50 values. The complex formed between Betv1 and the best docking pose for each ligand has been optimized by molecular dynamics simulation. Here, we describe the ligand binding of Betv1, which provides insight into the biological function of this protein. This knowledge is required for structural alteration or inhibition of some of these ligands in order to modify the allergenic properties of this protein.     

Enlarged FAMSBASE: protein 3D structure models of genome sequences for 41 species

Enlarged FAMSBASE is a relational database of comparative protein structure models for the whole genome of 41 species, presented in the GTOP database. The models are calculated by Full Automatic Modeling System (FAMS). Enlarged FAMSBASE provides a wide range of query keys, such as name of ORF (open reading frame), ORF keywords, Protein Data Bank (PDB) ID, PDB heterogen atoms and sequence similarity. Heterogen atoms in PDB include cofactors, ligands and other factors that interact with proteins, and are a good starting point for analyzing interactions between proteins and other molecules. The data may also work as a template for drug design. The present number of ORFs with protein 3D models in FAMSBASE is 183 805, and the database includes an average of three models for each ORF. FAMSBASE is available at http://famsbase.bio.nagoya-u.ac.jp/famsbase/.     

Variable intron/exon structure in the oligochaete lombricine kinase gene

Lombricine kinase is an annelid enzyme that belongs to the phosphagen kinase family of which creatine kinase and arginine kinase are the typical representatives. The enzymes play important roles in the cellular energy metabolism of animals. Biochemical, physiological and molecular information with respect to lombricine kinase is limited compared to other phosphagen kinases. This study presents data on the cDNA sequences of lombricine kinase from two smaller oligochaetes, Enchytraeus sp. and Stylaria sp. The deduced amino acid sequences are analyzed and compared with other selected phosphagen kinases. The intron/exon structure of the lombricine kinase gene was determined for these two species as well as two additional oligochaetes, Lumbriculus variegatus and Tubifex tubifex, and compared with available data for annelid phosphagen kinases. The data indicate the existence of a variable organization of the proposed 8-intron/9-exon gene structure. The results provide further insights in the evolution and position of these enzymes within the phosphagen kinase family.     

Novel information on the epitope of an inverse agonist monoclonal antibody provides insight into the structure of the TSH receptor

The TSH receptor (TSHR) comprises an extracellular leucine-rich domain (LRD) linked by a hinge region to the transmembrane domain (TMD). Insight into the orientation of these components to each other is required for understanding how ligands activate the receptor. We previously identified residue E251 at the LRD-hinge junction as contributing to coupling TSH binding with receptor activation. However, a single residue cannot stabilize the LRD-hinge unit. Therefore, based on the LRD crystal structure we selected for study four other potential LRD-hinge interface charged residues. Alanine substitutions of individual residues K244, E247, K250 and R255 (as well as previously known E251A) did not affect TSH binding or function. However, the cumulative mutation of these residues in varying permutations, primarily K250A and R255A when associated with E251A, partially uncoupled TSH binding and function. These data suggest that these three residues, spatially very close to each other at the LRD base, interact with the hinge region. Unexpectedly and most important, monoclonal antibody CS-17, a TSHR inverse agonist whose epitope straddles the LRD-hinge, was found to interact with residues K244 and E247 at the base of the convex LRD surface. These observations, together with the functional data, exclude residues K244 and E247 from the TSHR LRD-hinge interface. Further, for CS-17 accessibility to K244 and E247, the concave surface of the TSHR LRD must be tilted forwards towards the hinge region and plasma membrane. Overall, these data provide insight into the mechanism by which ligands either activate the TSHR or suppress its constitutive activity.     

Crystal structure of D-erythronate-4-phosphate dehydrogenase complexed with NAD

Pyridoxal-5′-phosphate (the active form of vitamin B6) is an essential cofactor in many enzymatic reactions. While animals lack any of the pathways for de novo synthesis and salvage of vitamin B6, it is synthesized by two distinct biosynthetic routes in bacteria, fungi, parasites, and plants. One of them is the PdxA/PdxJ pathway found in the γ subdivision of proteobacteria. It depends on the pdxB gene, which encodes erythronate-4-phosphate dehydrogenase (PdxB), a member of the d-isomer specific 2-hydroxyacid dehydrogenase superfamily. Although three-dimensional structures of other functionally related dehydrogenases are available, no structure of PdxB has been reported. To provide the missing structural information and to gain insights into the catalytic mechanism, we have determined the first crystal structure of erythronate-4-phosphate dehydrogenase from Pseudomonas aeruginosa in the ligand-bound state. It is a homodimeric enzyme consisting of 380-residue subunits. Each subunit consists of three structural domains: the lid domain, the nucleotide-binding domain, and the C-terminal dimerization domain. The latter domain has a unique fold and is largely responsible for dimerization. Interestingly, two subunits of the dimeric enzyme are bound with different combinations of ligands in the crystal and they display significantly different conformations. Subunit A is bound with NAD and a phosphate ion, while subunit B, with a more open active site cleft, is bound with NAD and l(+)-tartrate. Our structural data allow a detailed understanding of cofactor and substrate recognition, thus providing substantial insights into PdxB catalysis.