IntFOLD: an integrated server for modelling protein structures and functions from amino acid sequences

IntFOLD is an independent web server that integrates our leading methods for structure and function prediction. The server provides a simple unified interface that aims to make complex protein modelling data more accessible to life scientists. The server web interface is designed to be intuitive and integrates a complex set of quantitative data, so that 3D modelling results can be viewed on a single page and interpreted by non-expert modellers at a glance. The only required input to the server is an amino acid sequence for the target protein. Here we describe major performance and user interface updates to the server, which comprises an integrated pipeline of methods for: tertiary structure prediction, global and local 3D model quality assessment, disorder prediction, structural domain prediction, function prediction and modelling of protein-ligand interactions. The server has been independently validated during numerous CASP (Critical Assessment of Techniques for Protein Structure Prediction) experiments, as well as being continuously evaluated by the CAMEO (Continuous Automated Model Evaluation) project. The IntFOLD server is available at: http://www.reading.ac.uk/bioinf/IntFOLD/

Integrating protein structural information

Dissertação apresentada para obtenção de Grau de Doutor em Bioquímica,Bioquímica Estrutural, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Three-Dimensional Modelling of Honeybee Venom Allergenic Proteases: Relation to Allergenicity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelling of the nitrogen deposition and dietary lysine requirements of Redbro broilers

The aim of this study was to determine the coefficients of the Goettingen model for Redbro birds and estimate the digestible lysine requirements. To determine the model parameters, three nitrogen balance trials were performed in Periods I (14-28 days), II (42-56 days) and III (70-84 days), using 42 birds per trial. The birds were individually housed and subjected to six diets with increasing levels of nitrogen, with lysine as the limiting amino acid (deficient by 20% in relation to other amino acids). Dietary nitrogen concentrations were 8, 16, 24, 32, 40 and 48 g/kg. A control diet was added to confirm lysine as the first limiting amino acid. Nitrogen balance trials were divided into 5 days of adaptation and two periods of excreta collection, each one of 5 days. The response of the birds to a control diet confirmed that lysine was the first limiting amino acid. The adjustment of the exponential functions between nitrogen retention or excretion and nitrogen intake allowed estimation of parameters of the Goettingen model. The maximum potential for nitrogen retention was 3276, 2585 and 2603 mg/BWkg0.67.day, nitrogen maintenance requirement was 225, 135 and 122 mg/BWkg0.67.day and efficiency of nitrogen utilisation was 313 x 10(-6), 406 x 10(-6) and 415 x 10(-6) in the phases of 14-28, 42-56 and 70-84 days. The digestible lysine intake for Periods I, II and III, based on 60% of the maximum potential for nitrogen retention, was 711, 989 and 1272 mg/day (1.225%, 1.137% and 1.09% of lysine in the diet for a daily feed intake of 58, 87 and 117 g/day), respectively.

Molecular characterization of a miraculin-like gene differentially expressed during coffee development and coffee leaf miner infestation

The characterization of a coffee gene encoding a protein similar to miraculin-like proteins, which are members of the plant Kunitz serine trypsin inhibitor (STI) family of proteinase inhibitors (PIs), is described. PIs are important proteins in plant defence against insects and in the regulation of proteolysis during plant development. This gene has high identity with the Richadella dulcifica taste-modifying protein miraculin and with the tomato protein LeMir; and was named as CoMir (Coffea miraculin). Structural protein modelling indicated that CoMir had structural similarities with the Kunitz STI proteins, but suggested specific folding structures. CoMir was up-regulated after coffee leaf miner (Leucoptera coffella) oviposition in resistant plants of a progeny derived from crosses between C. racemosa (resistant) and C. arabica (susceptible). Interestingly, this gene was down-regulated during coffee leaf miner herbivory in susceptible plants. CoMir expression was up-regulated after abscisic acid application and wounding stress and was prominent during the early stages of flower and fruit development. In situ hybridization revealed that CoMir transcripts accumulated in the anther tissues that display programmed cell death (tapetum, endothecium and stomium) and in the metaxylem vessels of the petals, stigma and leaves. In addition, the recombinant protein CoMir shows inhibitory activity against trypsin. According to the present results CoMir may act in proteolytic regulation during coffee development and in the defence against L. coffeella. The similarity of CoMir with other Kunitz STI proteins and the role of CoMir in plant development and plant stress are discussed.

Abrogation of HMX1 function causes rare oculoauricular syndrome associated with congenital cataract, anterior segment dysgenesis, and retinal dystrophy.

PURPOSE: To define the phenotypic manifestation, confirm the genetic basis, and delineate the pathogenic mechanisms underlying an oculoauricular syndrome (OAS). METHODS: Two individuals from a consanguineous family underwent comprehensive clinical phenotyping and electrodiagnostic testing (EDT). Genome-wide microarray analysis and Sanger sequencing of the candidate gene were used to identify the likely causal variant. Protein modelling, Western blotting, and dual luciferase assays were used to assess the pathogenic effect of the variant in vitro. RESULTS: Complex developmental ocular abnormalities of congenital cataract, anterior segment dysgenesis, iris coloboma, early-onset retinal dystrophy, and abnormal external ear cartilage presented in the affected family members. Genetic analyses identified a homozygous c.650A>C; p.(Gln217Pro) missense mutation within the highly conserved homeodomain of the H6 family homeobox 1 (HMX1) gene. Protein modelling predicts that the variant may have a detrimental effect on protein folding and/or stability. In vitro analyses were able to demonstrate that the mutation has no effect on protein expression but adversely alters function. CONCLUSIONS: Oculoauricular syndrome is an autosomal recessive condition that has a profound effect on the development of the external ear, anterior segment, and retina, leading to significant visual loss at an early age. This study has delineated the phenotype and confirmed HMX1 as the gene causative of OAS, enabling the description of only the second family with the condition. HMX1 is a key player in ocular development, possibly in both the pathway responsible for lens and retina development, and via the gene network integral to optic fissure closure.

The cytochrome P450 aromatase lacking exon 5 is associated with a phenotype of nonclassic aromatase deficiency and is also present in normal human steroidogenic tissues

OBJECTIVE: The previously described c655G>A mutation of the human cytochrome P450 aromatase gene (P450aro, CYP19) results in aberrant splicing due to disruption of a donor splice site. To explain the phenotype of partial aromatase deficiency observed in a female patient described with this mutation, molecular consequences of the c655G>A mutation were investigated. DESIGN: To investigate whether the c655G>A mutation causes an aberrant spliced mRNA lacking exon 5 (-Ex5), P450aro RNA was analysed from the patient's lymphocytes by reverse transcription polymerase chain reaction (RT-PCR) and by splicing assays performed in Y1 cells transfected with a P450aro -Ex5 expression vector. Aromatase activity of the c655G>A mutant was predicted by three dimensional (3D) protein modelling studies and analysed in transiently transfected Y1 cells. Exon 5 might be predicted as a poorly defined exon suggesting a susceptibility to both splicing mutations and physiological alternative splicing events. Therefore, expression of the -Ex5 mRNA was also assessed as a possibly naturally occurring alternative splicing transcript in normal human steroidogenic tissues. PATIENTS: An aromatase deficient girl was born with ambiguous genitalia. Elevated serum LH, FSH and androgens, as well as cystic ovaries, were found during prepuberty. At the age of 8.4 years, spontaneous breast development and a 194.6 pmol/l serum oestradiol level was observed. RESULTS: The -Ex5 mRNA was found in lymphocytes of the P450aro deficient girl and her father, who was a carrier of the mutation. Mutant minigene expression resulted in complete exon 5 skipping. As expected from 3D protein modelling, -Ex5 cDNA expression in Y1 cells resulted in loss of P450aro activity. In addition, the -Ex5 mRNA was present in placenta, prepubertal testis and adrenal tissues. CONCLUSIONS: Alternative splicing of exon 5 of the CYP19 gene occurs in the wild type (WT) as well as in the c655G>A mutant. We speculate that for the WT it might function as a regulatory mechanism for aromatization, whereas for the mutant a relative prevalence of the shorter over the full-length protein might explain the phenotype of partial aromatase deficiency.

Toward the discovery of vaccine adjuvants:coupling in silico screening and in vitro analysis of antagonist binding to human and mouse CCR4 receptors

Background Adjuvants enhance or modify an immune response that is made to an antigen. An antagonist of the chemokine CCR4 receptor can display adjuvant-like properties by diminishing the ability of CD4+CD25+ regulatory T cells (Tregs) to down-regulate immune responses. Methodology Here, we have used protein modelling to create a plausible chemokine receptor model with the aim of using virtual screening to identify potential small molecule chemokine antagonists. A combination of homology modelling and molecular docking was used to create a model of the CCR4 receptor in order to investigate potential lead compounds that display antagonistic properties. Three-dimensional structure-based virtual screening of the CCR4 receptor identified 116 small molecules that were calculated to have a high affinity for the receptor; these were tested experimentally for CCR4 antagonism. Fifteen of these small molecules were shown to inhibit specifically CCR4-mediated cell migration, including that of CCR4+ Tregs. Significance Our CCR4 antagonists act as adjuvants augmenting human T cell proliferation in an in vitro immune response model and compound SP50 increases T cell and antibody responses in vivo when combined with vaccine antigens of Mycobacterium tuberculosis and Plasmodium yoelii in mice.

Influence of Plasma Protein Binding on Pharmacodynamics: Estimation of In Vivo Receptor Affinities of beta Blockers Using a New Mechanism-Based PK-PD Modelling Approach

The objective of this investigation was to examine in a systematic manner the influence of plasma protein binding on in vivo pharmacodynamics. Comparative pharmacokinetic-pharmacodynamic studies with four beta blockers were performed in conscious rats, using heart rate under isoprenaline-induced tachycardia as a pharmacodynamic endpoint. A recently proposed mechanism-based agonist-antagonist interaction model was used to obtain in vivo estimates of receptor affinities (K(B),(vivo)). These values were compared with in vitro affinities (K(B),(vitro)) on the basis of both total and free drug concentrations. For the total drug concentrations, the K(B),(vivo) estimates were 26, 13, 6.5 and 0.89 nM for S(-)-atenolol, S(-)-propranolol, S(-)-metoprolol and timolol. The K(B),(vivo) estimates on the basis of the free concentrations were 25, 2.0, 5.2 and 0.56 nM, respectively. The K(B),(vivo)-K(B),(vitro) correlation for total drug concentrations clearly deviated from the line of identity, especially for the most highly bound drug S(-)-propranolol (ratio K(B),(vivo)/K(B),(vitro) similar to 6.8). For the free drug, the correlation approximated the line of identity. Using this model, for beta-blockers the free plasma concentration appears to be the best predictor of in vivo pharmacodynamics. (C) 2008 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 98:3816-3828, 2009

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.

Diclofenac plasma protein binding: PK-PD modelling in cardiac patients submitted to cardiopulmonary bypass

Twenty-four surgical patients of both sexes without cardiac, hepatic, renal or endocrine dysfunctions were divided into two groups: 10 cardiac surgical patients submitted to myocardial revascularization and cardiopulmonary bypass (CPB), 3 females and 7 males aged 65 ± 11 years, 74 ± 16 kg body weight, 166 ± 9 cm height and 1.80 ± 0.21 m2 body surface area (BSA), and control, 14 surgical patients not submitted to CPB, 11 female and 3 males aged 41 ± 14 years, 66 ± 14 kg body weight, 159 ± 9 cm height and 1.65 ± 0.16 m2 BSA (mean ± SD). Sodium diclofenac (1 mg/kg, im Voltaren 75® twice a day) was administered to patients in the Recovery Unit 48 h after surgery. Venous blood samples were collected during a period of 0-12 h and analgesia was measured by the visual analogue scale (VAS) during the same period. Plasma diclofenac levels were measured by high performance liquid chromatography. A two-compartment open model was applied to obtain the plasma decay curve and to estimate kinetic parameters. Plasma diclofenac protein binding decreased whereas free plasma diclofenac levels were increased five-fold in CPB patients. Data obtained for analgesia reported as the maximum effect (EMAX) were: 25% VAS (CPB) vs 10% VAS (control), P<0.05, median measured by the visual analogue scale where 100% is equivalent to the highest level of pain. To correlate the effect versus plasma diclofenac levels, the EMAX sigmoid model was applied. A prolongation of the mean residence time for maximum effect (MRTEMAX) was observed without any change in lag-time in CPB in spite of the reduced analgesia reported for these patients, during the time-dose interval. In conclusion, the extent of plasma diclofenac protein binding was influenced by CPB with clinically relevant kinetic-dynamic consequences

Delineation and modelling of a nucleolar retention signal in the coronavirus nucleocapsid protein

Unlike nuclear localization signals, there is no obvious consensus sequence for the targeting of proteins to the nucleolus. The nucleolus is a dynamic subnuclear structure which is crucial to the normal operation of the eukaryotic cell. Studying nucleolar trafficking signals is problematic as many nucleolar retention signals (NoRSs) are part of classical nuclear localization signals (NLSs). In addition, there is no known consensus signal with which to inform a study. The avian infectious bronchitis virus (IBV), coronavirus nucleocapsid (N) protein, localizes to the cytoplasm and the nucleolus. Mutagenesis was used to delineate a novel eight amino acid motif that was necessary and sufficient for nucleolar retention of N protein and colocalize with nucleolin and fibrillarin. Additionally, a classical nuclear export signal (NES) functioned to direct N protein to the cytoplasm. Comparison of the coronavirus NoRSs with known cellular and other viral NoRSs revealed that these motifs have conserved arginine residues. Molecular modelling, using the solution structure of severe acute respiratory (SARS) coronavirus N-protein, revealed that this motif is available for interaction with cellular factors which may mediate nucleolar localization. We hypothesise that the N-protein uses these signals to traffic to and from the nucleolus and the cytoplasm.

Gastrointestinal modelling and optimization of protein and phosphorus bioaccessibility by phytase in fish diets

Programa de doctorado: Acuicultura

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.

Functional recycling of C2 domains throughout evolution: A comparative study of synaptotagmin, protein kinase C and phospholipase C by sequence, structural and modelling approaches

The C2 domain is one of the most frequent and widely distributed calcium-binding motifs. Its structure comprises an eight-stranded beta-sandwich with two structural types as if the result of a circular permutation. Combining sequence, structural and modelling information, we have explored, at different levels of granularity, the functional characteristics of several families of C2 domains. At the coarsest level,the similarity correlates with key structural determinants of the C2 domain fold and, at the finest level, with the domain architecture of the proteins containing them, highlighting the functional diversity between the various subfamilies. The functional diversity appears as different conserved surface patches throughout this common fold. In some cases, these patches are related to substrate-binding sites whereas in others they correspond to interfaces of presumably permanent interaction between other domains within the same polypeptide chain. For those related to substrate-binding sites, the predictions overlap with biochemical data in addition to providing some novel observations. For those acting as protein-protein interfaces' our modelling analysis suggests that slight variations between families are a result of not only complementary adaptations in the interfaces involved but also different domain architecture. In the light of the sequence and structural genomic projects, the work presented here shows that modelling approaches along with careful sub-typing of protein families will be a powerful combination for a broader coverage in proteomics. (C) 2003 Elsevier Ltd. All rights reserved.