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.

Characterization of an Allosteric Pocket in Zika Virus NS2B-NS3 Protease.

The NS2B-NS3 protease from Zika virus (ZIKV NS2B-NS3pro) cleaves the viral polyprotein, being essential for its replication and a therapeutic target. Inhibitors that target the active site of ZIKV NS2B-NS3pro have been developed, but they tend to have unfavorable pharmacokinetic properties due to their highly positive charge. Thus, the characterization of allosteric sites in this protease provides new strategies for inhibitor development. Here, we characterized a new allosteric pocket in ZIKV NS2B-NS3pro, analogous to the one previously described for the dengue virus protease. Molecular dynamics simulations indicate the presence of cavities around the residue Ala125, sampling protein conformations in which they are connected to the active site. This link between the residue Ala125 and the active site residues was reinforced by correlation network analysis. To experimentally verify the existence of this allosteric mechanism, we expressed and purified the Ala125Cys mutant of ZIKV NS2B-NS3pro and demonstrated that this variant is inhibited by the thiol-containing chemical probes 5,5'-dithiobis-(2-nitrobenzoic acid) and aldrithiol, which do not affect the activity of the wild-type protein. Inhibition of the mutant protein is reversed by the addition of strong reducing agents, supporting the involvement of Cys125 in covalent bond formation and enzyme inhibition. Together, our results provide experimental evidence for an allosteric pocket in ZIKV NS2B-NS3pro, in the region around Ala125, and computational insights on the structural connection between this region and the enzyme active site.

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.

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.

A new method bridging graph theory and residue co-evolutionary networks for specificity determinant positions detection.

MOTIVATION: Computational studies of molecular evolution are usually performed from a multiple alignment of homologous sequences, on which sequences resulting from a common ancestor are aligned so that equivalent residues are placed in the same position. Residues frequency patterns of a full alignment or from a subset of its sequences can be highly useful for suggesting positions under selection. Most methods mapping co-evolving or specificity determinant sites are focused on positions, however, they do not consider the case for residues that are specificity determinants in one subclass, but variable in others. In addition, many methods are impractical for very large alignments, such as those obtained from Pfam, or require a priori information of the subclasses to be analyzed. RESULTS: In this paper we apply the complex networks theory, widely used to analyze co-affiliation systems in the social and ecological contexts, to map groups of functional related residues. This methodology was initially evaluated in simulated environments and then applied to four different protein families datasets, in which several specificity determinant sets and functional motifs were successfully detected. AVAILABILITY AND IMPLEMENTATION: The algorithms and datasets used in the development of this project are available on http://www.biocomp.icb.ufmg.br/biocomp/software-and-databases/networkstats/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Coevolved Positions Represent Key Functional Properties in the Trypsin-Like Serine Proteases Protein Family.

Trypsin-like serine proteases are a group of homologous enzymes which exert multiple roles in both vertebrate and invertebrate organisms. Key properties of these enzymes include their activation from an inactive zymogen form to their active form by cleavage of residues in their N-terminus, the presence of a conserved catalytic triad of residues, and the existence of different patterns of substrate selectivity for residue cleavage between the various members of this protein family. In this article, we apply the decomposition of residue coevolution networks computational method to find sets of residues related to some of these key properties, especially to zymogen activation. Positive selection detection, normal modes analysis, and the calculation of thermal couplings between the bovine trypsinogen and bovine trypsin structures residues yielded further information for understanding the zymogen activation process and highlighted the importance of some of the coevolved set residues during these transitions.

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.

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.

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.

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.

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.

Her2p molecular modeling, mutant analysis and intramitochondrial localization.

Bacterial GatCAB amidotransferases are responsible for the transamidation of mischarged glutamyl-tRNA(Gln) into glutaminyl-tRNA(Gln). Mitochondria matrix also has a multienzymatic complex necessary for the transamidation of glutamyl-tRNA(Gln). Gtf1p, Her2p and Pet112p are the constituents of mitochondrial GatFAB amidotransferase complex. Her2p is subunit A of GatFAB complex, while Gtf1p is subunit F, a connector protein between Pet112p (subunit B) and Her2p. Here we evaluate through molecular modeling and amino acid correlation analysis the HER2 protein family. Localization studies indicated that Her2p is predominantly localized in the mitochondrial outer membrane, but it is also located in the mitochondrial matrix where together with Pet112p and Gtf1p constitutes the GatFAB complex. Finally, HER2 random mutagenesis unveiled important residues that provide thermo stability for the complex and are differently suppressed by overexpression of GTF1 or PET112. For instance, her2/ts11 mutant showed its fermentative growth impaired, and poorly rescued by GTF1 indicating that Her2p unknown function in the mitochondria outer membrane affects cell viability.

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.

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.

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.

Resurrecting Golgi proteins to grasp Golgi ribbon formation and self-association under stress.

The Golgi complex is an essential organelle of the eukaryotic exocytic pathway. A subfamily of Golgi matrix proteins, called GRASPs, is central in stress-induced unconventional secretion, Golgi dynamics during mitosis/apoptosis, and Golgi ribbon formation. The Golgi ribbon is vertebrate-specific and correlates with the appearance of two GRASP paralogues and two Golgins (GM130/Golgin45), which form specific GRASP-Golgin pairs. The molecular details of their appearance only in Metazoans are unknown. Moreover, despite new functionalities supported by GRASP paralogy, little is known about their structural and evolutionary differences. Here, we used ancestor sequence reconstruction and biophysical/biochemical approaches to assess the evolution of GRASPs structure/dynamics, fibrillation, and how they started anchoring their Golgin partners. Our data showed that a GRASP ancestor anchored Golgins before gorasp gene duplication in Metazoans. After gene duplication, variations within the GRASP binding pocket determined which paralogue would recruit which Golgin. These interactions are responsible for their specific Golgi location and Golgi ribbon appearance. We also suggest that GRASPs have a long-standing capacity to form supramolecular structures, affecting their participation in stress-induced processes.

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.

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.

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.