Purification and partial physical-chemical characterization of a new bovine trypsin proteoform (zeta-trypsin)

Recent advancements in enzyme research have unveiled a new proteoform of bovine trypsin, expanding our understanding of this well-characterized enzyme. While generally similar to other trypsins, this novel proteoform comprises three polypeptide chains, marking a significant difference in activity, kinetic properties, and conformational stability. Compared with the already known bovine trypsin proteoforms, the results showed a lower: activity, kcat and kcat.KM-1 and protein 'foldedness' ratio for the new proteoform. Molecular autolysis, a common feature in trypsin and chymotrypsin, has been explored through comparative physical chemistry properties with other proteoforms. This new proteoform of trypsin not only enriches the existing enzyme repertoire but also promises to shed light on the intricate physiological pathway for enzyme inactivation. Our results suggest that the new trypsin proteoform is one of the likely final pathways for enzyme inactivation in a physiological environment. This discovery opens up new avenues for further research into the functional implications of this new trypsin proteoform.

Design of a multi-epitope vaccine (vme-VAC/MST-1) against cholera and vibriosis based on reverse vaccinology and immunoinformatics approaches

Vibriosis and cholera are serious diseases distributed worldwide and caused by six marine bacteria of the Vibrio genus. Thousands of deaths occur each year due to these illnesses, necessitating the development of new preventive measures. Presently, the existing cholera vaccine demonstrates an effectiveness of approximately 60%. Here we describe a new multi-epitope vaccine, 'vme-VAC/MST-1' based on vaccine targets identified by reverse vaccinology and epitopes predicted by immunoinformatics, two currently effective tools for predicting new vaccines for bacterial pathogens. The vaccine was designed to combat vibriosis and cholera by incorporating epitopes predicted for CTL, HTL, and B cells. These epitopes were identified from six vaccine targets revealed through subtractive genomics, combined with reverse vaccinology, and were further filtered using immunoinformatics approaches based on their predicted immunogenicity. To construct the vaccine, 28 epitopes (24 CTL/B and 4 HTL/B) were linked to the sequence of the cholera toxin B subunit adjuvant. In silico analyses indicate that the resulting immunogen is stable, soluble, non-toxic, and non-allergenic. Furthermore, it exhibits no homology to the host and demonstrates a strong capacity to elicit innate, B-cell, and T-cell immune responses. Our analysis suggests that it is likely to elicit immune reactions mediated through the TLR5 pathway, as evidenced by the molecular docking of the vaccine with the receptor, which revealed high affinity and a favorable reaction. Thus, vme-VAC/MST-1 is predicted to be a safe and effective solution against pathogenic Vibrio spp. However, further experimental analyses are required to measure the vaccine's effects In vivo.Communicated by Ramaswamy H. Sarma.

Prediction of enzymatic function with high efficiency and a reduced number of features using genetic algorithm

The post-genomic era has raised a growing demand for efficient procedures to identify protein functions, which can be accomplished by applying machine learning to the characteristics set extracted from the protein. This approach is feature-based and has been the focus of several works in bioinformatics. In this work, we investigated the characteristics of proteins, representing the primary, secondary, tertiary, and quaternary structures of the protein, that improve the model's quality by applying dimensionality reduction techniques and using the Support Vector Machine classifier for predicting the enzymes' classes. During the investigation, two approaches were evaluated: feature extraction/transformation, which was performed using the statistical technique Factor Analysis, and feature selection methods. For feature selection, we proposed an approach based on a genetic algorithm to face the optimization conflict between the simplicity and reliability of an ideal representation of the characteristics of the enzymes and also compared and employed other methods for this purpose. The best result was accomplished using a feature subset generated by our implementation of a multi-objective genetic algorithm enriched with features that this work identified as relevant to represent the enzymes. This subset representation reduced the dataset by about 87% and reached 85.78% of F-measure performance, improving the overall quality of the model classification. In addition, we verified in this work a subset addressed with only 28 features out of a total of 424 that reached a performance above 80% of F-measure for four of the six evaluated classes, showing that satisfactory classification performance can be achieved with a reduced number of enzymes's characteristics. The datasets and implementations are openly available.

Structural and biochemical characterization of Leptospira interrogans Lsa45 reveals a penicillin-binding protein with esterase activity

Leptospirosis is a bacterial disease that affects humans and animals and is caused by Leptospira. The recommended treatment for leptospirosis is antibiotic therapy, which should be given early in the course of the disease. Despite the use of these antibiotics, their role during the course of the disease is still not completely clear because of the lack of effective clinical trials, particularly for severe cases of the disease. Here, we present the characterization of L. interrogans Lsa45 protein by gel filtration, protein crystallography, SAXS, fluorescence and enzymatic assays. The oligomeric studies revealed that Lsa45 is monomeric in solution. The crystal structure of Lsa45 revealed the presence of two subdomains: a large α/β subdomain and a small α-helical subdomain. The large subdomain contains the amino acids Ser122, Lys125, and Tyr217, which correspond to the catalytic triad that is essential for β-lactamase or serine hydrolase activity in similar enzymes. Additionally, we also confirmed the bifunctional promiscuity of Lsa45, in hydrolyzing both the 4-nitrophenyl acetate (p-NPA) and nitrocefin β-lactam antibiotic. Therefore, this study provides novel insights into the structure and function of enzymes from L. interrogans, which furthers our understanding of this bacterium and the development of new therapies for the prevention and treatment of leptospirosis.

UNEXPECTED PLASTICITY OF THE QUATERNARY STRUCTURE OF IRON-MANGANESE SUPEROXIDE DISMUTASES

Protein 3D structure can be remarkably robust to the accumulation of mutations during evolution. On the other hand, sometimes a single amino acid substitution can be sufficient to generate dramatic and completely unpredictable structural consequences. In an attempt to rationally alter the preferences for the metal ion at the active site of a member of the Iron/Manganese superoxide dismutase family, two examples of the latter phenomenon were identified. Site directed mutants of SOD from Trichoderma reesei were generated and studied crystallographically together with the wild type enzyme. Despite being chosen for their potential impact on the redox potential of the metal, two of the mutations (D150G and G73A) in fact resulted in significant alterations to the protein quaternary structure. The D150G mutant presented alternative inter-subunit contacts leading to a loss of symmetry of the wild type tetramer, whereas the G73A mutation transformed the tetramer into an octamer despite not participating directly in any of the inter-subunit interfaces. We conclude that there is considerable intrinsic plasticity in the Fe/MnSOD fold that can be unpredictably affected by single amino acid substitutions. In much the same way as phenotypic defects at the organism level can reveal much about normal function, so too can such mutations teach us much about the subtleties of protein structure.

Naview: A d3.js Based JavaScript Library for Drawing and Annotating Voltage-Gated Sodium Channels Membrane Diagrams

Voltage-gated sodium channels (Nav) are membrane proteins essential to initiating and propagating action potential in neurons and other excitable cells. For a given organism there are often multiple, specialized sodium channels found in different tissues, whose mutations can cause deleterious effects observed in numerous diseases. Consequently, there is high medical and pharmacological interest in these proteins. Scientific literature often uses membrane diagrams to depict important patterns in these channels including the six transmembrane segments (S1–S6) present in four different homologous domains (D1–D4), the S4 voltage sensors, the pore-lining residue segments and the ion selectivity filter residues, glycosylation and phosphorylation residues, toxin binding sites and the inactivation loop, among others. Most of these diagrams are illustrated either digitally or by hand and programs specifically dedicated to the interactive and data-friendly generation of such visualizations are scarce or non-existing. This paper describes Naview, an open-source javascript visualization compatible with modern web browsers for the dynamic drawing and annotation of voltage-gated sodium channels membrane diagrams based on the D3.js library. By using a graphical user interface and combining user-defined annotations with optional UniProt code as inputs, Naview allows the creation and customization of membrane diagrams. In this interface, a user can also map and display important sodium channel properties, residues, regions and their relationships through symbols, colors, and edge connections. Such features can facilitate data exploration and provide fast, high-quality publication-ready graphics for this highly active area of research.

Coq3p relevant residues for protein activity and stability

Coenzyme Q (CoQ) is an essential molecule that consists of a highly substituted benzene ring attached to a polyprenyl tail anchored in the inner mitochondrial membrane. CoQ transfers electrons from NADH dehydrogenase and succinate dehydrogenase complexes toward ubiquinol-cytochrome c reductase, and that allows aerobic growth of cells. In Saccharomyces cerevisiae, the synthesis of CoQ depends on fourteen proteins Coq1p-Co11p, Yah1p, Arh1p, and Hfd1p. Some of these proteins are components of CoQ synthome. Using ab initio molecular modeling and site-directed mutagenesis, we identified the functional residues of the O-methyltransferase Coq3p, which depends on S-adenosylmethionine for catalysis and is necessary for two O-methylation steps required for CoQ maturation. Conserved residues as well as those that coevolved in the protein structure were found to have important roles in respiratory growth, CoQ biosynthesis, and also in the stability of CoQ synthome proteins. Finally, a multiple sequence alignment showed that S. cerevisiae Coq3p has a 45 amino acid residues insertion that is poorly conserved or absent in oleaginous yeast, cells that can store up to 20% of their dry weight as lipids. These results point to the Coq3p structural determinants of its biological and catalytic function and could contribute to the development of lipid-producing yeast for biotechnology.

CONAN: a web application to detect specificity determinants and functional sites by amino acids co-variation network analysis

SUMMARY

CONAN is a web application developed to detect specificity determinants and function-related sites by amino acids co-variation networks analysis, emphasizing local coevolutionary constraints. The software allows the characterization of structurally and functionally relevant groups of residues and their relationship with subsets of sequences by automatic cross-referencing with GO terms, UniprotKb annotations and INTERPRO.

AVAILABILITY AND IMPLEMENTATION

CONAN is free and open-source, being distributed in the terms of the GPLV3 license. The software is available as a web application and python script versions and can be accessed at http://bioinfo.icb.ufmg.br/conan. We also provide running instructions, the source code and a user guide.

Reenacting the Birth of a Function: Functional Divergence of HIUases and Transthyretins as Inferred by Evolutionary and Biophysical Studies

Transthyretin was discovered in the 1940s, named after its ability to bind thyroid hormones and retinol. In the genomic era, transthyretins were found to be part of a larger family with homologs of no obvious function, then called transthyretin-related proteins. Thus, it was proposed that the transthyretin gene could be the result of gene duplication of an ancestral of this newly identified homolog, later found out to be an enzyme involved in uric acid degradation, then named HIUase (5-hydroxy-isourate hydrolase). Here, we sought to re-enact the evolutionary history of this protein family by reconstructing, from a phylogeny inferred from 123 vertebrate sequences, three ancestors corresponding to key moments in their evolution-before duplication; the common transthyretin ancestor after gene duplication and the common ancestor of Eutheria transthyretins. Experimental and computational characterization showed the reconstructed ancestor before duplication was unable to bind thyroxine and likely presented the modern HIUase reaction mechanism, while the substitutions after duplication prevented that activity and were enough to provide stable thyroxine binding, as confirmed by calorimetry and x-ray diffraction. The Eutheria transthyretin ancestor was less prone to characterization, but limited data suggested thyroxine binding as expected. Sequence/structure analysis suggests an early ability to bind the Retinol Binding Protein. We solved the X-ray structures from the two first ancestors, the first at 1.46 resolution, the second at 1.55 resolution with well-defined electron density for thyroxine, providing a useful tool for the understanding of structural adaptation from enzyme to hormone distributor.

Glutantβase: a database for improving the rational design of glucose-tolerant β-glucosidases

Β-glucosidases are key enzymes used in second-generation biofuel production. They act in the last step of the lignocellulose saccharification, converting cellobiose in glucose. However, most of the β-glucosidases are inhibited by high glucose concentrations, which turns it a limiting step for industrial production. Thus, β-glucosidases have been targeted by several studies aiming to understand the mechanism of glucose tolerance, pH and thermal resistance for constructing more efficient enzymes. In this paper, we present a database of β-glucosidase structures, called Glutantβase. Our database includes 3842 GH1 β-glucosidase sequences collected from UniProt. We modeled the sequences by comparison and predicted important features in the 3D-structure of each enzyme. Glutantβase provides information about catalytic and conserved amino acids, residues of the coevolution network, protein secondary structure, and residues located in the channel that guides to the active site. We also analyzed the impact of beneficial mutations reported in the literature, predicted in analogous positions, for similar enzymes. We suggested these mutations based on six previously described mutants that showed high catalytic activity, glucose tolerance, or thermostability (A404V, E96K, H184F, H228T, L441F, and V174C). Then, we used molecular docking to verify the impact of the suggested mutations in the affinity of protein and ligands (substrate and product). Our results suggest that only mutations based on the H228T mutant can reduce the affinity for glucose (product) and increase affinity for cellobiose (substrate), which indicates an increment in the resistance to product inhibition and agrees with computational and experimental results previously reported in the literature. More resistant β-glucosidases are essential to saccharification in industrial applications. However, thermostable and glucose-tolerant β-glucosidases are rare, and their glucose tolerance mechanisms appear to be related to multiple and complex factors. We gather here, a set of information, and made predictions aiming to provide a tool for supporting the rational design of more efficient β-glucosidases. We hope that Glutantβase can help improve second-generation biofuel production. Glutantβase is available at http://bioinfo.dcc.ufmg.br/glutantbase .

Structural and immunological characterization of a new nucleotidyltransferase-like antigen from Paracoccidioides brasiliensis

Pb27 antigen is an interesting alternative to immunological diagnosis of Paracoccidioidomycosis (PCM) and has demonstrated to be protective in experimental PCM. Its tertiary structure and possible function remained unknown till now. To study Pb27 at the atomic level, the recombinant protein was expressed in Escherichia coli BL21(DE3), purified, and its three-dimensional structure was solved by X-ray crystallography. Based on this structure, we performed a residue correlation analysis and in silico ligand search assays to address a possible biological function to Pb27. We identified Pb27 as a member of the extensive nucleotidyltransferase superfamily. The protein has an αβαβαβ topology with two domains (N- and C-terminal domains) and adopts a monomeric form as its biological unit in solution. Structural comparisons with similar members of the superfamily clearly indicate Pb27 C-terminal domain is singular and may play an important role in its biological function. Bioinformatics analysis suggested that Pb27 might bind to ATP and CTP. This suggestion is corroborated by the fact that a magnesium cation is coordinated by two aspartic acid residues present at the active site (between N- and C-terminal domains), as evidenced by X-ray diffraction data. Besides, NMR assays (¹H-¹⁵N HSQC spectra) confirmed the binding of CTP to Pb27, demonstrating for the first time an interaction between a nucleotide and this protein. Moreover, we evaluated the reactivity of sera from patients with Paracoccidioides brasiliensis infection against the recombinant form of Pb27 and showed that it was recognized by sera from infected and treated patients. Predicted B and T cell epitopes were synthesized and further evaluated against sera of PCM patients, providing information of the most reactive peptides in Pb27 primary structure which interact with specific Pb27 antibodies.

Closed- and open-state models of human skeletal muscle sodium channel

Voltage-gated sodium channels play important roles in human physiology. However, their complexity hinders the understanding of their physiology and pathology at atomic level. We took advantage of the structural reports of similar channels obtained by cryo-EM (EeNav1.4, and NavPaS), and constructed models of human Nav1.4 channels at closed and open states. The open-state model is very similar to the recently published cryo-EM structure of hNav1.4. The comparison of both models shows shifts of the voltage sensors (VS) of DIII and DIV. The activated position of VS-DII in the closed model was demonstrated by Ts1 docking, thereby confirming the requirement that VS-DI, VS-DII and VS-DIII must be activated for the channel to open. The interactions observed with VS-DIII suggest a stepwise, yet fast, transition from resting to activated state. These models provide structural insights on the closed-open transition of the channel.

Thyroxine-binding globulin deficiency due to a novel SERPINA7 mutation: Clinical characterization, analysis of X-chromosome inactivation pattern and protein structural modeling

OBJECTIVE

Thyroxine-binding globulin (TBG) is the major human thyroid hormone transport protein, encoded by the SERPINA7 gene (Xq22.2). We aim to investigate the molecular basis of partial TBG deficiency (TBG-PD) in a female, by evaluating the X-chromosome inactivation pattern as well as the mutant protein structural modeling.

DESIGN AND METHODS

Sequencing of the coding region of the SERPINA7 gene was performed in a female with a TBG-PD phenotype and her first-degree relatives. The proband presented with low serum levels of total T3 (TT3) and total T4 (TT4), serum TSH level of 5.4 μUI/mL (normal range, 0.35-5.5), and serum TBG level of 5.5 mg/L (normal range, 13.6-27.2). X-chromosome inactivation pattern was evaluated by methylation analysis of the androgen receptor gene (Xq11.2). Structural analysis of the SERPIN family was performed using Pymol and Areaimol, and PFSTATS for conservation analysis and family-wide investigation of equivalent positions in human homologs. Modeller was used for point mutation structural modeling.

RESULTS

A novel missense SERPINA7 mutation (p.R35W; c.163C > T) was found in heterozygosity in the proband, and in hemizygosity in her affected siblings. The proband X-chromosome inactivation ratio was 20:80. The substitution of an arginine by a tryptophan is predicted to disrupt the protein surface and main electrostatic interactions. Tryptophans are extremely rare (0.1%) in this position.

CONCLUSIONS

We report a new SERPINA7 variant associated with TBG-PD in three siblings. We named this variant TBG-Brasilia. The X-chromosome inactivation pattern may have accounted for the rare phenotypic expression in a female. The hydrophobic nature of the mutant is predicted to create an apolar patch at the surface, which results in protein aggregation and/or misfolding, potentially responsible for thyroid hormone transport defect.

Mitochondrial ribosome bL34 mutants present diminished translation of cytochrome c oxidase subunits

Saccharomyces cerevisiae mitoribosomes are specialized in the translation of a few number of highly hydrophobic membrane proteins, components of the oxidative phosphorylation system. Mitochondrial characteristics, such as the membrane system and its redox state driven mitoribosomes evolution through great diversion from their bacterial and cytosolic counterparts. Therefore, mitoribosome presents a considerable number of mitochondrial-specific proteins, as well as new protein extensions. In this work we characterize temperature sensitive mutants of the subunit bL34 present in the 54S large subunit. Although bL34 has bacterial homologs, in yeast it has a long 65 aminoacids mitochondrial N-terminal addressing sequence, here we demonstrate that it can be replaced by the mitochondrial addressing sequence of Neurospora crassa ATP9 gene. The bL34 temperature sensitive mutants present lowered translation of mitochondrial COX1 and COX3, which resulted in reduced cytochrome c oxidase activity and respiratory growth deficiency. The sedimentation properties of bL34 in sucrose gradients suggest that similarly to its bacterial homolog, bL34 is also a later participant in the process of mitoribosome biogenesis.

Novel arsenic-transforming bacteria and the diversity of their arsenic-related genes and enzymes arising from arsenic-polluted freshwater sediment

Bacteria are essential in arsenic cycling. However, few studies have addressed 16S rRNA and arsenic-related functional gene diversity in long-term arsenic-contaminated tropical sediment. Here, using culture-based, metagenomic and computational approaches, we describe the diversity of bacteria, genes and enzymes involved in AsIII and AsV transformation in freshwater sediment and in anaerobic AsIII- and AsV-enrichment cultures (ECs). The taxonomic profile reveals significant differences among the communities. Arcobacter, Dechloromonas, Sedimentibacter and Clostridium thermopalmarium were exclusively found in ECs, whereas Anaerobacillus was restricted to AsV-EC. Novel taxa that are both AsV-reducers and AsIII-oxidizers were identified: Dechloromonas, Acidovorax facilis, A. delafieldii, Aquabacterium, Shewanella, C. thermopalmarium and Macellibacteroides fermentans. Phylogenic discrepancies were revealed among the aioA, arsC and arrA genes and those of other species, indicating horizontal gene transfer. ArsC and AioA have sets of amino acids that can be used to assess their functional and structural integrity and familial subgroups. The positions required for AsV reduction are conserved, suggesting strong selective pressure for maintaining the functionality of ArsC. Altogether, these findings highlight the role of freshwater sediment bacteria in arsenic mobility, and the untapped diversity of dissimilatory arsenate-reducing and arsenate-resistant bacteria, which might contribute to arsenic toxicity in aquatic environments.

Structural and functional mapping of Rtg2p determinants involved in retrograde signaling and aging of Saccharomyces cerevisiae

In Saccharomyces cerevisiae mitochondrial dysfunction induces retrograde signaling, a pathway of communication from mitochondria to the nucleus that promotes a metabolic remodeling to ensure sufficient biosynthetic precursors for replication. Rtg2p is a positive modulator of this pathway that is also required for cellular longevity. This protein belongs to the ASKHA superfamily, and contains a putative N-terminal ATP-binding domain, but there is no detailed structural and functional map of the residues in this domain that accounts for their contribution to retrograde signaling and aging. Here we use Decomposition of Residue Correlation Networks and site-directed mutagenesis to identify Rtg2p structural determinants of retrograde signaling and longevity. We found that most of the residues involved in retrograde signaling surround the ATP-binding loops, and that Rtg2p N-terminus is divided in three regions whose mutants have different aging phenotypes. We also identified E137, D158 and S163 as possible residues involved in stabilization of ATP at the active site. The mutants shown here may be used to map other Rtg2p activities that crosstalk to other pathways of the cell related to genomic stability and aging.

Identification of new bacteria harboring qnrS and aac(6')-Ib/cr and mutations possibly involved in fluoroquinolone resistance in raw sewage and activated sludge samples from a full-scale WWTP

Wastewater treatment plants (WWTPs) harbor bacteria and antimicrobial resistance genes, favoring gene exchange events and resistance dissemination. Here, a culture-based and metagenomic survey of qnrA, qnrB, qnrS, and aac(6')-Ib genes from raw sewage (RS) and activated sludge (AS) of a full-scale municipal WWTP was performed. A total of 96 bacterial isolates were recovered from nalidixic acid-enrichment cultures. Bacteria harboring the aac(6')-Ib gene predominated in RS, whereas qnrS-positive isolates were specific to AS. Novel qnrS- and aac(6')-Ib-cr positive species were identified: Morganella morganii, Providencia rettgeri, and Pseudomonas guangdongensis (qnrS), and Alcaligenes faecalis and P. rettgeri (aac(6')-Ib-cr). Analysis of qnrS and aac(6')-Ib sequences from isolates and clone libraries suggested that the diversity of qnrS is wider than that of aac(6')-Ib. A large number of amino acid mutations were observed in the QnrS and AAC(6')-Ib proteins at previously undetected positions, whose structural implications are not clear. An accumulation of mutations at the C72, Q73, L74, A75 and M76 positions of QnrS, and D181 of AAC(6')-Ib might be important for resistance. These findings add significant information on bacteria harboring qnrS and aac(6')-Ib genes, and the presence of novel mutations that may eventually emerge in clinical isolates.

Sequence, structure and function relationships in flaviviruses as assessed by evolutive aspects of its conserved non-structural protein domains

Flaviviruses are responsible for serious diseases such as dengue, yellow fever, and zika fever. Their genomes encode a polyprotein which, after cleavage, results in three structural and seven non-structural proteins. Homologous proteins can be studied by conservation and coevolution analysis as detected in multiple sequence alignments, usually reporting positions which are strictly necessary for the structure and/or function of all members in a protein family or which are involved in a specific sub-class feature requiring the coevolution of residue sets. This study provides a complete conservation and coevolution analysis on all flaviviruses non-structural proteins, with results mapped on all well-annotated available sequences. A literature review on the residues found in the analysis enabled us to compile available information on their roles and distribution among different flaviviruses. Also, we provide the mapping of conserved and coevolved residues for all sequences currently in SwissProt as a supplementary material, so that particularities in different viruses can be easily analyzed.

Residue correlation networks in nuclear receptors reflect functional specialization and the formation of the nematode-specific P-box

Background

Nuclear receptors (NRs) are transcription factors which bind small hormones, whose evolutionary history and the presence of different functional surfaces makes them an interesting target for a correlation based analysis.

Results

Correlation analysis of ligand binding domains shows that correlated residue subsets arise from the differences between functional sites in different nuclear receptor subfamilies. For the DNA binding domain, particularly, the analysis shows that the main source of correlation comes from residues that regulate hormone response element specificity, and one of the conserved residue sub-sets arises due to the presence of an unusual sequence for the DNA binding motif known as P-box in nematodes, suggesting the existence of different DBD-DNA specificities in nuclear receptors.

Conclusions

We conclude that DNA specificity and functional surface specialization has independently driven nuclear receptor evolution, and suggest possible binding modes for the class of divergent nuclear receptors in nematodes.

Using amino acid correlation and community detection algorithms to identify functional determinants in protein families

Correlated mutation analysis has a long history of interesting applications, mostly in the detection of contact pairs in protein structures. Based on previous observations that, if properly assessed, amino acid correlation data can also provide insights about functional sub-classes in a protein family, we provide a complete framework devoted to this purpose. An amino acid specific correlation measure is proposed, which can be used to build networks summarizing all correlation and anti-correlation patterns in a protein family. These networks can be submitted to community structure detection algorithms, resulting in subsets of correlated amino acids which can be further assessed by specific parameters and procedures that provide insight into the relationship between different communities, the individual importance of community members and the adherence of a given amino acid sequence to a given community. By applying this framework to three protein families with contrasting characteristics (the Fe/Mn-superoxide dismutases, the peroxidase-catalase family and the C-type lysozyme/α-lactalbumin family), we show how our method and the proposed parameters and procedures are related to biological characteristics observed in these protein families, highlighting their potential use in protein characterization and gene annotation.

Crystallization and preliminary structural analysis of the giant haemoglobin from Glossoscolex paulistus at 3.2 Å

Glossoscolex paulistus is a free-living earthworm encountered in south-east Brazil. Its oxygen transport requirements are undertaken by a giant extracellular haemoglobin, or erythrocruorin (HbGp), which has an approximate molecular mass of 3.6 MDa and, by analogy with its homologue from Lumbricus terrestris (HbLt), is believed to be composed of a total of 180 polypeptide chains. In the present work the full 3.6 MDa particle in its cyanomet state was purified and crystallized using sodium citrate or PEG8000 as precipitant. The crystals contain one-quarter of the full particle in the asymmetric unit of the I222 cell and have parameters of a = 270.8 Å, b = 320.3 Å and c = 332.4 Å. Diffraction data were collected to 3.15 Å using synchrotron radiation on beamline X29A at the Brookhaven National Laboratory and represent the highest resolution data described to date for similar erythrocruorins. The structure was solved by molecular replacement using a search model corresponding to one-twelfth of its homologue from HbLt. This revealed that HbGp belongs to the type I class of erythrocruorins and provided an interpretable initial electron density map in which many features including the haem groups and disulfide bonds could be identified.

Novel Zn2+-binding sites in human transthyretin: implications for amyloidogenesis and retinol-binding protein recognition

Human transthyretin (TTR) is a homotetrameric protein involved in several amyloidoses. Zn(2+) enhances TTR aggregation in vitro, and is a component of ex vivo TTR amyloid fibrils. We report the first crystal structure of human TTR in complex with Zn(2+) at pH 4.6-7.5. All four structures reveal three tetra-coordinated Zn(2+)-binding sites (ZBS 1-3) per monomer, plus a fourth site (ZBS 4) involving amino acid residues from a symmetry-related tetramer that is not visible in solution by NMR. Zn(2+) binding perturbs loop E-α-helix-loop F, the region involved in holo-retinol-binding protein (holo-RBP) recognition, mainly at acidic pH; TTR affinity for holo-RBP decreases ∼5-fold in the presence of Zn(2+). Interestingly, this same region is disrupted in the crystal structure of the amyloidogenic intermediate of TTR formed at acidic pH in the absence of Zn(2+). HNCO and HNCA experiments performed in solution at pH 7.5 revealed that upon Zn(2+) binding, although the α-helix persists, there are perturbations in the resonances of the residues that flank this region, suggesting an increase in structural flexibility. While stability of the monomer of TTR decreases in the presence of Zn(2+), which is consistent with the tertiary structural perturbation provoked by Zn(2+) binding, tetramer stability is only marginally affected by Zn(2+). These data highlight structural and functional roles of Zn(2+) in TTR-related amyloidoses, as well as in holo-RBP recognition and vitamin A homeostasis.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: implications to tetramer stability and ligand-binding

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis.

Crystal structure and statistical coupling analysis of highly glycosylated peroxidase from royal palm tree (Roystonea regia)

Royal palm tree peroxidase (RPTP) is a very stable enzyme in regards to acidity, temperature, H(2)O(2), and organic solvents. Thus, RPTP is a promising candidate for developing H(2)O(2)-sensitive biosensors for diverse applications in industry and analytical chemistry. RPTP belongs to the family of class III secretory plant peroxidases, which include horseradish peroxidase isozyme C, soybean and peanut peroxidases. Here we report the X-ray structure of native RPTP isolated from royal palm tree (Roystonea regia) refined to a resolution of 1.85A. RPTP has the same overall folding pattern of the plant peroxidase superfamily, and it contains one heme group and two calcium-binding sites in similar locations. The three-dimensional structure of RPTP was solved for a hydroperoxide complex state, and it revealed a bound 2-(N-morpholino) ethanesulfonic acid molecule (MES) positioned at a putative substrate-binding secondary site. Nine N-glycosylation sites are clearly defined in the RPTP electron-density maps, revealing for the first time conformations of the glycan chains of this highly glycosylated enzyme. Furthermore, statistical coupling analysis (SCA) of the plant peroxidase superfamily was performed. This sequence-based method identified a set of evolutionarily conserved sites that mapped to regions surrounding the heme prosthetic group. The SCA matrix also predicted a set of energetically coupled residues that are involved in the maintenance of the structural folding of plant peroxidases. The combination of crystallographic data and SCA analysis provides information about the key structural elements that could contribute to explaining the unique stability of RPTP.

Crystal structure of a soluble decoy receptor IL-22BP bound to interleukin-22

Interleukin-22 (IL-22) plays an important role in the regulation of immune and inflammatory responses in mammals. The IL-22 binding protein (IL-22BP), a soluble receptor that specifically binds IL-22, prevents the IL-22/interleukin-22 receptor 1 (IL-22R1)/interleukin-10 receptor 2 (IL-10R2) complex assembly and blocks IL-22 biological activity. Here we present the crystal structure of the IL-22/IL-22BP complex at 2.75 A resolution. The structure reveals IL-22BP residues critical for IL-22 binding, which were confirmed by site-directed mutagenesis and functional studies. Comparison of IL-22/IL-22BP and IL-22/IL-22R1 crystal structures shows that both receptors display an overlapping IL-22 binding surface, which is consistent with the inhibitory role played by IL-22 binding protein.

Crystal structure of the IL-22/IL-22R1 complex and its implications for the IL-22 signaling mechanism

Interleukin-22 (IL-22) is a member of the interleukin-10 cytokine family, which is involved in anti-microbial defenses, tissue damage protection and repair, and acute phase responses. Its signaling mechanism involves the sequential binding of IL-22 to interleukin-22 receptor 1 (IL-22R1), and of this dimer to interleukin-10 receptor 2 (IL-10R2) extracellular domain. We report a 1.9A crystal structure of the IL-22/IL-22R1 complex, revealing crucial interacting residues at the IL-22/IL-22R1 interface. Functional importance of key residues was confirmed by site-directed mutagenesis and functional studies. Based on the X-ray structure of the binary complex, we discuss a molecular basis of the IL-22/IL-22R1 recognition by IL-10R2.

Crystal structure of yeast hexokinase PI in complex with glucose: A classical "induced fit" example revised

Hexokinase is the first enzyme in the glycolytic pathway that catalyzes the transfer of a phosphoryl group from ATP to glucose to form glucose-6-phosphate and ADP. Two yeast hexokinase isozymes are known, namely PI and PII. Here we redetermined the crystal structure of yeast hexokinase PI from Saccharomyces cerevisiae as a complex with its substrate, glucose, and refined it at 2.95 A resolution. Comparison of the holo-PI yeast hexokinase and apo-hexokinase structures shows in detail the rigid body domain closure and specific loop movements as glucose binds and sheds more light on structural basis of the "induced fit" mechanism of reaction in the HK enzymatic action. We also performed statistical coupling analysis of the hexokinase family, which reveals two co-evolved continuous clusters of amino acid residues and shows that the evolutionary coupled amino acid residues are mostly confined to the active site and the hinge region, further supporting the importance of these parts of the protein for the enzymatic catalysis.

Structural basis of GC-1 selectivity for thyroid hormone receptor isoforms

Background

Thyroid receptors, TRα and TRβ, are involved in important physiological functions such as metabolism, cholesterol level and heart activities. Whereas metabolism increase and cholesterol level lowering could be achieved by TRβ isoform activation, TRα activation affects heart rates. Therefore, β-selective thyromimetics have been developed as promising drug-candidates for treatment of obesity and elevated cholesterol level. GC-1 [3,5-dimethyl-4-(4'-hydroxy-3'-isopropylbenzyl)-phenoxy acetic acid] has ability to lower LDL cholesterol with 600- to 1400-fold more potency and approximately two- to threefold more efficacy than atorvastatin (Lipitor^©) in studies in rats, mice and monkeys.

Results

To investigate GC-1 specificity, we solved crystal structures and performed molecular dynamics simulations of both isoforms complexed with GC-1. Crystal structures reveal that, in TRα Arg228 is observed in multiple conformations, an effect triggered by the differences in the interactions between GC-1 and Ser277 or the corresponding asparagine (Asn331) of TRβ. The corresponding Arg282 of TRβ is observed in only one single stable conformation, interacting effectively with the ligand. Molecular dynamics support this model: our simulations show that the multiple conformations can be observed for the Arg228 in TRα, in which the ligand interacts either strongly with the ligand or with the Ser277 residue. In contrast, a single stable Arg282 conformation is observed for TRβ, in which it strongly interacts with both GC-1 and the Asn331.

Conclusion

Our analysis suggests that the key factors for GC-1 selectivity are the presence of an oxyacetic acid ester oxygen and the absence of the amino group relative to T₃. These results shed light into the β-selectivity of GC-1 and may assist the development of new compounds with potential as drug candidates to the treatment of hypercholesterolemia and obesity.

Orphan nuclear receptor NGFI-B forms dimers with nonclassical interface

The orphan receptor nerve growth factor-induced B (NGFI-B) is a member of the nuclear receptor's subfamily 4A (Nr4a). NGFI-B was shown to be capable of binding both as a monomer to an extended half-site containing a single AAAGGTCA motif and also as a homodimer to a widely separated everted repeat, as opposed to a large number of nuclear receptors that recognize and bind specific DNA sequences predominantly as homo- and/or heterodimers. To unveil the structural organization of NGFI-B in solution, we determined the quaternary structure of the NGFI-B LBD by a combination of ab initio procedures from small-angle X-ray scattering (SAXS) data and hydrogen-deuterium exchange followed by mass spectrometry. Here we report that the protein forms dimers in solution with a radius of gyration of 2.9 nm and maximum dimension of 9.0 nm. We also show that the NGFI-B LBD dimer is V-shaped, with the opening angle significantly larger than that of classical dimer's exemplified by estrogen receptor (ER) or retinoid X receptor (RXR). Surprisingly, NGFI-B dimers formation does not occur via the classical nuclear receptor dimerization interface exemplified by ER and RXR, but instead, involves an extended surface area composed of the loop between helices 3 and 4 and C-terminal fraction of the helix 3. Remarkably, the NGFI-B dimer interface is similar to the dimerization interface earlier revealed for glucocorticoid nuclear receptor (GR), which might be relevant to the recognition of cognate DNA response elements by NGFI-B and to antagonism of NGFI-B-dependent transcription exercised by GR in cells.

Low-resolution structure and fluorescence anisotropy analysis of protein tyrosine phosphatase eta catalytic domain

The rat protein tyrosine phosphatase eta, rPTPeta, is a class I "classical" transmembrane RPTP, with an intracellular portion composed of a unique catalytic region. The rPTPeta and the human homolog DEP-1 are downregulated in rat and human neoplastic cells, respectively. However, the malignant phenotype is reverted after exogenous reconstitution of rPTPeta, suggesting that its function restoration could be an important tool for gene therapy of human cancers. Using small-angle x-ray scattering (SAXS) and biophysical techniques, we characterized the intracellular catalytic domain of rat protein tyrosine phosphatase eta (rPTPetaCD) in solution. The protein forms dimers in solution as confirmed by SAXS data analysis. The SAXS data also indicated that rPTPetaCD dimers are elongated and have an average radius of gyration of 2.65 nm and a D(max) of 8.5 nm. To further study the rPTPetaCD conformation in solution, we built rPTPetaCD homology models using as scaffolds the crystallographic structures of RPTPalpha-D1 and RPTPmicro-D1 dimers. These models were, then, superimposed onto ab initio low-resolution SAXS structures. The structural comparisons and sequence alignment analysis of the putative dimerization interfaces provide support to the notion that the rPTPetaCD dimer architecture is more closely related to the crystal structure of autoinhibitory RPTPalpha-D1 dimer than to the dimeric arrangement exemplified by RPTPmicro-D1. Finally, the characterization of rPTPetaCD by fluorescence anisotropy measurements demonstrates that the dimer dissociation is concentration dependent with a dissociation constant of 21.6 +/- 2.0 microM.

Structural Rearrangements in the Thyroid Hormone Receptor Hinge Domain and Their Putative Role in the Receptor Function

The thyroid hormone receptor (TR) D-domain links the ligand-binding domain (LBD, EF-domain) to the DNA-binding domain (DBD, C-domain), but its structure, and even its existence as a functional unit, are controversial. The D domain is poorly conserved throughout the nuclear receptor family and was originally proposed to comprise an unfolded hinge that facilitates rotation between the LBD and the DBD. Previous TR LBD structures, however, have indicated that the true unstructured region is three to six amino acid residues long and that the D-domain N terminus folds into a short amphipathic alpha-helix (H0) contiguous with the DBD and that the C terminus of the D-domain comprises H1 and H2 of the LBD. Here, we solve structures of TR-LBDs in different crystal forms and show that the N terminus of the TRalpha D-domain can adopt two structures; it can either fold into an amphipathic helix that resembles TRbeta H0 or form an unstructured loop. H0 formation requires contacts with the AF-2 coactivator-binding groove of the neighboring TR LBD, which binds H0 sequences that resemble coactivator LXXLL motifs. Structural analysis of a liganded TR LBD with small angle X-ray scattering (SAXS) suggests that AF-2/H0 interactions mediate dimerization of this protein in solution. We propose that the TR D-domain has the potential to form functionally important extensions of the DBD and LBD or unfold to permit TRs to adapt to different DNA response elements. We also show that mutations of the D domain LXXLL-like motif indeed selectively inhibit TR interactions with an inverted palindromic response element (F2) in vitro and TR activity at this response element in cell-based transfection experiments.

Recombinant human receptors and functional assays in the discovery of altinicline (SIB-1508Y), a novel acetylcholine-gated ion channel (nAChR) agonist

Neuronal nicotinic acetylcholine receptors (nAChRs) are a class of ion channels with significant potential as molecular targets for the design of drugs to treat a variety of CNS disorders. The discovery that neuronal nAChRs are further subdivided into multiple subtypes suggests that drugs which act selectively at specific nAChR subtypes might effectively treat Parkinson's disease (PD), Alzheimer's disease (AD), schizophrenia, ADHD, depression, anxiety or pain without the accompanying adverse side effects associated with non-selective agents such as nicotine (1) and epibatidine. Altinicline (SIB-1508Y) is a novel, small molecule designed to selectively activate neuronal nAChRs and is undergoing clinical evaluation for the treatment of PD. It was selected from a series of compounds primarily on the basis of results from functional assays, including (a) measurement of Ca2+ flux in stable cell lines expressing specific recombinant human neuronal nAChR subtypes; (b) determination of in vitro and in vivo neurotransmitter release; (c) in vivo models of PD. Biological data on both altinicline and the series of compounds from which it was selected are reported.

SIB-1757 and SIB-1893: selective, noncompetitive antagonists of metabotropic glutamate receptor type 5

Cell lines expressing the human metabotropic glutamate receptor subtype 5a (hmGluR5a) and hmGluR1b were used as targets in an automated high-throughput screening (HTS) system that measures changes in intracellular Ca2+ ([Ca2+]i) using fluorescence detection. This functional screen was used to identify the mGluR5-selective antagonist, SIB-1757 [6-methyl-2-(phenylazo)-3-pyridinol], which inhibited the glutamate-induced [Ca2+]i responses at hmGluR5 with an IC50 of 0.37 microM compared with an IC50 of >100 microM at hmGluR1. Schild analysis demonstrated a noncompetitive mechanism of inhibition. Pharmacophore mapping was used to identify an additional compound, SIB-1893 [(E)-2-methyl-6-(2-phenylethenyl)pyridine], which was also shown to block glutamate-induced increases in [Ca2+]i at hmGluR5 with an IC50 of 0.29 microM compared with an IC50 of >100 microM at hmGluR1. SIB-1757 and SIB-1893 showed little or no activity when tested for agonist and antagonist activity at the other recombinant human mGluR subtypes, alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid, kainate, and N-methyl-D-aspartate receptors. In rat neonatal brain slices, SIB-1757 and SIB-1893 inhibited (S)-3,5-dihydroxyphenylglycine (DHPG)-evoked inositol phosphate accumulation in hippocampus and striatum by 60% to 80%, with a potency similar to that observed on recombinant mGluR5. However, in the cerebellum, a brain region with low mGluR5 expression, SIB-1757 failed to inhibit DHPG-evoked inositol phosphate accumulation. In cultured rat cortical neurons, SIB-1757 and SIB-1893 largely inhibited DHPG-evoked [Ca2+]i signals, revealing a population of neurons that were less sensitive to SIB-1757 and SIB-1893. This is the first description of highly selective, noncompetitive mGluR5 antagonists. These compounds will be useful tools in evaluating the role of mGluR5 in normal physiology and in animal models of disease.

Pharmacological characterization of SIB-1765F: a novel cholinergic ion channel agonist

Nicotine, the prototypical agonist for neuronal nicotinic acetylcholine receptors (NAChR), nonselectively activates NAChR limiting its use in elucidating the function of NAChR subtypes. SIB-1765F is a subtype selective NAChR agonist that displaces [3H]-nicotine binding with an IC50 of 4.6 nM and [3H]-cytisine binding with an IC50 of 12.2 nM which is 2000- to 6000-fold lower than its displacement of [3H]-QNB or [125I]-alpha-bungarotoxin. SIB-1765F did not inhibit human or rat cholinesterases or the uptake of [3H]-DA in synaptosomal preparations. SIB-1765F mimicked (-)-nicotine in stimulating [3H]-DA release from rat striatal and olfactory tubercle slices, with EC50 values of 99.6 and 39.6 microM, respectively. Such stimulation was sensitive to mecamylamine and DH beta E. SIB-1765F also released endogenous DA in the striatum and the nucleus accumbens as measured by in vivo microdialysis. SIB-1765F was less efficacious than (-)-nicotine at stimulating [3H]-NE release from rat hippocampal slices; in contrast, SIB-1765F increased [3H]-NE release from rat thalamic and cortical slices with efficacies approaching those of (-)-nicotine. Similar to (-)-nicotine and (+/-)-epibatidine, subcutaneous administration of SIB-1765F increased the turnover rate of dopamine ex vivo both in the striatum and olfactory tubercles in a mecamylamine-sensitive manner. Because the release of striatal DA and hippocampal NE appears to be regulated by distinct NAChR, differential effects of SIB-1765F on striatal DA and hippocampal NE release supports the NAChR subtype selectivity of SIB-1765F compared to (-)-nicotine. This is further demonstrated by observations showing that SIB-1765F has a higher affinity for h alpha 4 beta 2 NAChR relative to h alpha 4 beta 4 NAChRs in displacing [3H]-epibatidine binding and increasing cytosolic CA+2 concentration in cell lines stably expressing h alpha 4 beta 2 or h alpha 4 beta 4.