A clustering method for robust and reliable large scale functional and structural protein sequence annotation

Bioinformatics, in the last few decades, has played a fundamental role to give sense to the huge amount of data produced. Obtained the complete sequence of a genome, the major problem of knowing as much as possible of its coding regions, is crucial. Protein sequence annotation is challenging and, due to the size of the problem, only computational approaches can provide a feasible solution. As it has been recently pointed out by the Critical Assessment of Function Annotations (CAFA), most accurate methods are those based on the transfer-by-homology approach and the most incisive contribution is given by cross-genome comparisons. In the present thesis it is described a non-hierarchical sequence clustering method for protein automatic large-scale annotation, called “The Bologna Annotation Resource Plus” (BAR+). The method is based on an all-against-all alignment of more than 13 millions protein sequences characterized by a very stringent metric. BAR+ can safely transfer functional features (Gene Ontology and Pfam terms) inside clusters by means of a statistical validation, even in the case of multi-domain proteins. Within BAR+ clusters it is also possible to transfer the three dimensional structure (when a template is available). This is possible by the way of cluster-specific HMM profiles that can be used to calculate reliable template-to-target alignments even in the case of distantly related proteins (sequence identity < 30%). Other BAR+ based applications have been developed during my doctorate including the prediction of Magnesium binding sites in human proteins, the ABC transporters superfamily classification and the functional prediction (GO terms) of the CAFA targets. Remarkably, in the CAFA assessment, BAR+ placed among the ten most accurate methods. At present, as a web server for the functional and structural protein sequence annotation, BAR+ is freely available at http://bar.biocomp.unibo.it/bar2.0.

Identification of a fibronectin interaction site in the extracellular matrix protein ameloblastin

Mammalian teeth are composed of hydroxyapatite crystals that are embedded in a rich extracellular matrix. This matrix is produced by only two cell types, the mesenchymal odontoblasts and the ectodermal ameloblasts. Ameloblasts secrete the enamel proteins amelogenin, ameloblastin, enamelin and amelotin. Odontoblasts secrete collagen type I and several calcium-binding phosphoproteins including dentin sialophosphoprotein, dentin matrix protein, bone sialoprotein and osteopontin. The latter four proteins have recently been grouped in the family of the SIBLINGs (small integrin-binding ligand, N-linked glycoproteins) because they display similar gene structures and because they contain an RGD tripeptide sequence that binds to integrin receptors and thus mediates cell adhesion. We have prepared all the other tooth-specific proteins in recombinant form and examined whether they might also promote cell adhesion similar to the SIBLINGs. We found that only ameloblastin consistently mediated adhesion of osteoblastic and fibroblastic cells to plastic or titanium surfaces. The activity was dependent on the intact three-dimensional structure of ameloblastin and required de novo protein synthesis of the adhering cells. By deletion analysis and in vitro mutagenesis, the active site could be narrowed down to a sequence of 13 amino acid residues (VPIMDFADPQFPT) derived from exon 7 of the rat ameloblastin gene or exons 7-9 of the human gene. Kinetic studies and RNA interference experiments further demonstrated that this sequence does not directly bind to a cell surface receptor but that it interacts with cellular fibronectin, which in turn binds to integrin receptors. The identification of a fibronectin-binding domain in ameloblastin might permit interesting applications for dental implantology. Implants could be coated with peptides containing the active sequence, which in turn would recruit fibronectin from the patient's blood. The recruited fibronectin should then promote cell adhesion on the implant surface, thereby accelerating osseointegration of the implant.

Neuroelectromagnetic correlates of perceptual closure processes

Perceptual closure refers to the coherent perception of an object under circumstances when the visual information is incomplete. Although the perceptual closure index observed in electroencephalography reflects that an object has been recognized, the full spatiotemporal dynamics of cortical source activity underlying perceptual closure processing remain unknown so far. To address this question, we recorded magnetoencephalographic activity in 15 subjects (11 females) during a visual closure task and performed beamforming over a sequence of successive short time windows to localize high-frequency gamma-band activity (60–100 Hz). Two-tone images of human faces (Mooney faces) were used to examine perceptual closure. Event-related fields exhibited a magnetic closure index between 250 and 325 ms. Time-frequency analyses revealed sustained high-frequency gamma-band activity associated with the processing of Mooney stimuli; closure-related gamma-band activity was observed between 200 and 300 ms over occipitotemporal channels. Time-resolved source reconstruction revealed an early (0–200 ms) coactivation of caudal inferior temporal gyrus (cITG) and regions in posterior parietal cortex (PPC). At the time of perceptual closure (200–400 ms), the activation in cITG extended to the fusiform gyrus, if a face was perceived. Our data provide the first electrophysiological evidence that perceptual closure for Mooney faces starts with an interaction between areas related to processing of three-dimensional structure from shading cues (cITG) and areas associated with the activation of long-term memory templates (PPC). Later, at the moment of perceptual closure, inferior temporal cortex areas specialized for the perceived object are activated, i.e., the fusiform gyrus related to face processing for Mooney stimuli.

Spatial and functional relationships between air conduits and blood capillaries in the pulmonary gas exchange tissue of adult and developing chickens

The documented data regarding the three-dimensional structure of the air capillaries (ACs), the ultimate sites of gas exchange in the avian lung is contradictory. Further, the mode of gas exchange, described as cross-current has not been clearly elucidated. We studied the temporal and spatial arrangement of the terminal air conduits of the chicken lung and their relationship with the blood capillaries (BCs) in embryos as well as the definitive architecture in adults. Several visualization techniques that included corrosion casting, light microscopy as well as scanning and transmission electron microscopy were used. Two to six infundibulae extend from each atrium and give rise to numerous ACs that spread centrifugally. Majority of the ACs are tubular structures that give off branches, which anastomose with their neighboring cognates. Some ACs have globular shapes and a few are blind-ending tapering tubes. During inauguration, the luminal aspects of the ACs are characterized by numerous microvillus-like microplicae, which are formed during the complex processes of cell attenuation and canalization of the ACs. The parabronchial exchange BCs, initially inaugurated as disorganized meshworks, are reoriented via pillar formation to lie predominantly orthogonal to the long axes of the ACs. The remodeling of the retiform meshworks by intussusceptive angiogenesis essentially accomplishes a cross-current system at the gas exchange interface in the adults, where BCs form ring-like patterns around the ACs, thus establishing a cross-current system. Our findings clarify the mode of gas exchange in the parabronchial mantle and illuminate the basis for the functional efficiency of the avian lung.

2,3-Dichloro-1,4-hydroquinone 2,3-dichloro-1,4-benzoquinone monohydrate: a quinhydrone-type 1:1 donor-acceptor [D-A] charge-transfer complex

In the crystal structure of the title compound (systematic name: 2,3-dichlorobenzene-1,4-diol 2,3-dichlorocyclohexa-2,5-diene-1,4-dione monohydrate), C(6)H(4)Cl(2)O(2)center dot C(6)H(2)Cl(2)O(2)center dot H(2)O, the 2,3-dichloro-1,4-hydroquinone donor (D) and the 2,3-dichloro-1,4-benzoquinone acceptor (A) molecules form alternating stacks along [100]. Their molecular planes [maximum deviations for non-H atoms: 0.0133 (14) (D) and 0.0763 (14) angstrom (A)] are inclined to one another by 1.45 (3)degrees and are thus almost parallel. There are pi-pi interactions involving the D and A molecules, with centroid-centroid distances of 3.5043 (9) and 3.9548 (9) angstrom. Intermolecular O-H center dot center dot center dot O hydrogen bonds involving the water molecule and the hydroxy and ketone groups lead to the formation of two-dimensional networks lying parallel to (001). These networks are linked by C-H center dot center dot center dot O interactions, forming a three-dimensional structure.

Besnoitia besnoiti protein disulfide isomerase (BbPDI): molecular characterization, expression and in silico modelling

Besnoitia besnoiti is an apicomplexan parasite responsible for bovine besnoitiosis, a disease with a high prevalence in tropical and subtropical regions and re-emerging in Europe. Despite the great economical losses associated with besnoitiosis, this disease has been underestimated and poorly studied, and neither an effective therapy nor an efficacious vaccine is available. Protein disulfide isomerase (PDI) is an essential enzyme for the acquisition of the correct three-dimensional structure of proteins. Current evidence suggests that in Neosporacaninum and Toxoplasmagondii, which are closely related to B. besnoiti, PDI play an important role in host cell invasion, is a relevant target for the host immune response, and represents a promising drug target and/or vaccine candidate. In this work, we present the nucleotide sequence of the B. besnoiti PDI gene. BbPDI belongs to the thioredoxin-like superfamily (cluster 00388) and is included in the PDI_a family (cluster defined cd02961) and the PDI_a_PDI_a'_c subfamily (cd02995). A 3D theoretical model was built by comparative homology using Swiss-Model server, using as a template the crystallographic deduced model of Tapasin-ERp57 (PDB code 3F8U chain C). Analysis of the phylogenetic tree for PDI within the phylum apicomplexa reinforces the close relationship among B. besnoiti, N. caninum and T. gondii. When subjected to a PDI-assay based on the polymerisation of reduced insulin, recombinant BbPDI expressed in E. coli exhibited enzymatic activity, which was inhibited by bacitracin. Antiserum directed against recombinant BbPDI reacted with PDI in Western blots and by immunofluorescence with B. besnoiti tachyzoites and bradyzoites.

Crosslinking Studies To Probe Substrate Binding By Soybean Lipoxygenase-1

Soybean lipoxygenase-1 (SBLO-1) catalyzes the oxygenation of polyunsaturated fatty acids into conjugated diene hydroperoxides. The three dimensional structure of SBLO-1 is known, but it is not certain how substrates bind. One hypothesis involves the transient separation of helix-2 and helix-11 located on the exterior of the molecule in front of the active site iron. A second hypothesis involves a conformational change in the side chains of residues leucine 541 and threonine 259. To test these hypotheses, site directed mutagenesis was used to create a cysteine mutation on each helix, which could allow for the formation of a disulfide linkage. Disulfide formation between the two cysteines in the T259C,S545C mutant was found to be unfavorable, but later shown to be present at higher pH values using SDS-PAGE. Treatment of the T259C,S545C with the crosslinker 2,3-dibromomaleimide (DBM) resulted in a 50% reduction in catalytic activity. No loss of activity was observed when the single mutant, S545C, or the wild type was treated with DBM. Single mutants T259C and L541C both showed approximately 20% reduction in the rate after addition of DBM. Double mutants T259C,L541C and S263C,S545C showed approximately 30% reduction in the rate after addition of DBM. Single mutants T259C and L541C showed an increase in activity after incubation with NEM. Double mutants T259C,S545C and T259C,L541C showed an increase in activity after incubation with NEM. The S263C,S545C double mutant showed a slight decrease in activity in the presence of NEM. It is unclear how the NEM and DBM are interacting with the molecule, but this can easily be determined through mass spectrometry experiments.

Conformational analysis of a potent SSTR3-selective somatostatin analogue by NMR in water solution

The three-dimensional structure of a potent SSTR3-selective analogue of somatostatin, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-Agl(8)(N(beta) Me, 2-naphthoyl)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-Agl(8)(N(beta) Me, 2-naphthoyl)]-SRIF) (peptide 1) has been determined by (1)H NMR in water and molecular dynamics (MD) simulations. The peptide exists in two conformational isomers differing mainly by the cis/trans isomerization of the side chain in residue 8. The structure of 1 is compared with the consensus structural motifs of other somatostatin analogues that bind predominantly to SSTR1, SSTR2/SSTR5 and SSTR4 receptors, and to the 3D structure of a non-selective SRIF analogue, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-2Nal(8)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-2Nal(8)]-SRIF) (peptide 2). The structural determinant factors that could explain selectivity of peptide 1 for SSTR3 receptors are discussed.

Different amino acid substitutions at the same position in rhodopsin lead to distinct phenotypes

PURPOSE: Identification of a novel rhodopsin mutation in a family with retinitis pigmentosa and comparison of the clinical phenotype to a known mutation at the same amino acid position. METHODS: Screening for mutations in rhodopsin was performed in 78 patients with retinitis pigmentosa. All exons and flanking intronic regions were amplified by PCR, sequenced, and compared to the reference sequence derived from the National Center for Biotechnology Information (NCBI, Bethesda, MD) database. Patients were characterized clinically according to the results of best corrected visual acuity testing (BCVA), slit lamp examination (SLE), funduscopy, Goldmann perimetry (GP), dark adaptometry (DA), and electroretinography (ERG). Structural analyses of the rhodopsin protein were performed with the Swiss-Pdb Viewer program available on-line (http://www.expasy.org.spdvbv/ provided in the public domain by Swiss Institute of Bioinformatics, Geneva, Switzerland). RESULTS: A novel rhodopsin mutation (Gly90Val) was identified in a Swiss family of three generations. The pedigree indicated autosomal dominant inheritance. No additional mutation was found in this family in other autosomal dominant genes. The BCVA of affected family members ranged from 20/25 to 20/20. Fundus examination showed fine pigment mottling in patients of the third generation and well-defined bone spicules in patients of the second generation. GP showed concentric constriction. DA demonstrated monophasic cone adaptation only. ERG revealed severely reduced rod and cone signals. The clinical picture is compatible with retinitis pigmentosa. A previously reported amino acid substitution at the same position in rhodopsin leads to a phenotype resembling night blindness in mutation carriers, whereas patients reported in the current study showed the classic retinitis pigmentosa phenotype. The effect of different amino acid substitutions on the three-dimensional structure of rhodopsin was analyzed by homology modeling. Distinct distortions of position 90 (shifts in amino acids 112 and 113) and additional hydrogen bonds were found. CONCLUSIONS: Different amino acid substitutions at position 90 of rhodopsin can lead to night blindness or retinitis pigmentosa. The data suggest that the property of the substituted amino acid distinguishes between the phenotypes.

Pseudodominant Inheritance of Goitrous Congenital Hypothyroidism Caused by TPO Mutations: Molecular and in Silico Studies

Context and Objective: Most cases of goitrous congenital hypothyroidism (CH) from thyroid dyshormonogenesis 1) follow a recessive mode of inheritance and 2) are due to mutations in the thyroid peroxidase gene (TPO). We report the genetic mechanism underlying the apparently dominant inheritance of goitrous CH in a nonconsanguineous family of French Canadian origin. Design, Setting, and Participants: Two brothers identified by newborn TSH screening had severe hypothyroidism and a goiter with increased (99m)Tc uptake. The mother was euthyroid, but the father and two paternal uncles had also been diagnosed with goitrous CH. After having excluded PAX8 gene mutations, we hypothesized that the underlying defect could be TPO mutations. Results: Both compound heterozygous siblings had inherited a mutant TPO allele carried by their mother (c.1496delC; p.Pro499Argfs2X), and from their father, one brother had inherited a missense mutation (c.1978C-->G; p.Gln660Glu) and the other an insertion (c.1955insT; p.Phe653Valfs15X). The thyroid gland of one uncle who is a compound heterozygote for TPO mutations (p.Phe653Valfs15X/p.Gln660Glu) was removed because of concurrent multiple endocrine neoplasia type 2A. Immunohistochemistry revealed normal TPO staining, implying that Gln660Glu TPO is expressed properly. Modeling of this mutant in silico suggests that its three-dimensional structure is conserved, whereas the electrostatic binding energy between the Gln660Glu TPO and its heme group becomes repulsive. Conclusion: We report a pedigree presenting with pseudodominant goitrous CH due to segregation of three different TPO mutations. Although goitrous CH generally follows a recessive mode of inheritance, the high frequency of TPO mutations carriers may lead to pseudodominant inheritance.

Clinical, structural and functional implications of mutations and polymorphisms in human NADPH P450 oxidoreductase

Cytochrome P450 proteins are involved in metabolism of drugs and xenobiotics. In the endoplasmic reticulum a single nicotinamide adenine dinucleotide phosphate (NADPH) P450 oxidoreductase (POR) supplies electrons to all microsomal P450s for catalytic activity. POR is a flavoprotein that contains both flavin mononucleotide and flavin adenine dinucleotide as cofactors and uses NADPH as the source of electrons. We have recently reported a number of POR mutations in the patients with disordered steroidogenesis. In the first report we had described missense mutations (A287P, R457H, V492E, C569Y, and V608F) identified in four patients with defects in steroid production. Each POR variant was produced as recombinant N-27 form of the enzyme in bacteria and as full-length form in yeast. Membranes from bacteria or yeast expressing normal or variant POR were purified and their activities were characterized in cytochrome c and CYP17A1 assays. Later we have published a larger study that described a whole range of POR mutations and characterized the mutants/polymorphisms A115V, T142A, M263V, Y459H, A503V, G539R, L565P, R616X, V631I, and F646del from the sequencing of patient DNA. We also studied POR variants Y181D, P228L, R316W, G413S, and G504R that were available in public databases or published literature. Three-dimensional structure of rat POR is known and we have used this structure to deduce the structure-function correlation of POR mutations in human. The missense mutations found in patients with disordered steroidogenesis are generally in the co-factor binding and functionally important domains of POR and the apparent polymorphisms are found in regions with lesser structural importance. A variation in POR can alter the activity of all microsomal P450s, and therefore, can affect the metabolism of drugs and xenobiotics even when the P450s involved are otherwise normal. It is important to study the genetic and biochemical basis of POR variants in human population to gain information about possible differences in P450 mediated reactions among the individuals carrying a variant or polymorphic form of POR that could impact their metabolism.

Simulating in vivo-like Synaptic Input Patterns in Multicompartmental Models

In the laboratory of Dr. Dieter Jaeger at Emory University, we use computer simulations to study how the biophysical properties of neurons—including their three-dimensional structure, passive membrane resistance and capacitance, and active membrane conductances generated by ion channels—affect the way that the neurons transfer synaptic inputs into the action potential streams that represent their output. Because our ultimate goal is to understand how neurons process and relay information in a living animal, we try to make our computer simulations as realistic as possible. As such, the computer models reflect the detailed morphology and all of the ion channels known to exist in the particular neuron types being simulated, and the model neurons are tested with synaptic input patterns that are intended to approximate the inputs that real neurons receive in vivo. The purpose of this workshop tutorial was to explain what we mean by ‘in vivo-like’ synaptic input patterns, and how we introduce these input patterns into our computer simulations using the freely available GENESIS software package (http://www.genesis-sim.org/GENESIS). The presentation was divided into four sections: first, an explanation of what we are talking about when we refer to in vivo-like synaptic input patterns

Autosomal dominant neurohypophyseal diabetes insipidus in a Swiss family, caused by a novel mutation (C59Delta/A60W) in the neurophysin moiety of prepro-vasopressin-neurophysin II (AVP-NP II).

OBJECTIVE To study clinical, morphological and molecular characteristics in a Swiss family with autosomal dominant familial neurohypophyseal diabetes insipidus (adFNDI). PARTICIPANTS AND METHODS A 15-month-old girl presenting with symptoms of polydipsia and polyuria was investigated by water deprivation test. Evaluation of the family revealed three further family members with symptomatic vasopressin-deficient diabetes insipidus. T1-weighted magnetic resonance images of the posterior pituitary were taken in two affected adult family members and molecular genetic analysis was performed in all affected individuals. RESULTS The water deprivation test in the 15-month-old child confirmed the diagnosis of vasopressin-deficient diabetes insipidus and the pedigree was consistent with autosomal dominant inheritance. The characteristic bright spot of the normal vasopressin-containing neurophypophysis was absent in both adults with adFNDI. Direct sequence analysis revealed a new deletion (177-179DeltaCGC) in exon 2 of the AVP-NP II gene in all affected individuals. At the amino acid level, this deletion eliminates cysteine 59 (C59Delta) and substitutes alanine 60 by tryptophan (A60W) in the AVP-NP II precursor; interestingly, the remainder of the reading frame remains unchanged. According to the three-dimensional structure of neurophysin, C59 is involved in a disulphide bond with C65. CONCLUSIONS Deletion of C59 and substitution of A60W in the AVP-NP II precursor is predicted to disrupt one of the seven disulphide bridges required for correct folding of the neurophysin moiety and thus disturb the function of neurophysin as the vasopressin transport protein. These data are in line with the clinical and morphological findings in the reported family with adFNDI.

Interaction of the Trinuclear Triangular Secondary Building Unit [Cu 3 (μ 3 -OH)(μ-pz) 3 ] 2+ with 4,4′-Bipyridine. Structural Characterizations of New Coordination Polymers and Hexanuclear Cu II Clusters. 2°

By reacting 4,4′-bipyridine (bpy) with selected trinuclear triangular CuII complexes, [Cu3(μ3-OH)(μ-pz)3(RCOO)2(LL′)] [pz = pyrazolate anion; R = CH3, CH3CH2, CH2═CH, CH2═C(CH3); L, L′ = Hpz, H2O, MeOH] in MeOH, the substitution of monotopic ligands by ditopic bpy was observed. Depending on the stoichiometric reaction ratios, different compounds were isolated and structurally characterized. One- and two-dimensional coordination polymers (CPs), as well as two hexanuclear CuII clusters were identified. One of the hexanuclear clusters self-assembles into a supramolecular three-dimensional structure, and its crystal packing shows the presence of two intersecting channels, one of which is almost completely occupied by guest bpy, while in the second one guest water molecules are present. This compound also shows a reversible, thermally induced, single-crystal-to-single-crystal transition.

Luminal particles within cellular microtubules.

The regulation of microtubule dynamics is attributed to microtubule-associated proteins that bind to the microtubule outer surface, but little is known about cellular components that may associate with the internal side of microtubules. We used cryoelectron tomography to investigate in a quantitative manner the three dimensional structure of microtubules in intact mammalian cells. We show that the lumen of microtubules in this native state is filled with discrete, globular particles with a diameter of 7 nm and spacings between 8 and 20 nm in neuronal cells. Cross-sectional views of microtubules confirm the presence of luminal material in vitreous sections of brain tissue. Most of the luminal particles had connections to the microtubule wall, as revealed in tomograms. A higher accumulation of particles was seen near the retracting plus ends of microtubules. The luminal particles were abundant in neurons, but were also observed in other cells, such as astrocytes and stem cells.