Paulo Sergio Lopes de Oliveira

KVFinder: steered identification of protein cavities as a PyMOL plugin

BackgroundThe characterization of protein binding sites is a major challenge in computational biology. Proteins interact with a wide variety of molecules and understanding of such complex interactions is essential to gain deeper knowledge of protein function. Shape complementarity is known to be important in determining protein-ligand interactions. Furthermore, these protein structural features have been shown to be useful in assisting medicinal chemists during lead discovery and optimization.ResultsWe developed KVFinder, a highly versatile and easy-to-use tool for cavity prospection and spatial characterization. KVFinder is a geometry-based method that has an innovative customization of the search space. This feature provides the possibility of cavity segmentation, which alongside with the large set of customizable parameters, allows detailed cavity analyses. Although the main focus of KVFinder is the steered prospection of cavities, we tested it against a benchmark dataset of 198 known drug targets in order to validate our software and compare it with some of the largely accepted methods. Using the one click mode, we performed better than most of the other methods, staying behind only of hybrid prospection methods. When using just one of KVFinder’s customizable features, we were able to outperform all other compared methods. KVFinder is also user friendly, as it is available as a PyMOL plugin, or command-line version.ConclusionKVFinder presents novel usability features, granting full customizable and highly detailed cavity prospection on proteins, alongside with a friendly graphical interface. KVFinder is freely available on http://lnbio.cnpem.br/bioinformatics/main/software/.

Improved Proteomic Approach for the Discovery of Potential Vaccine Targets in Trypanosoma cruzi

Chagas disease, caused by Trypanosoma cruzi, is a devastating parasitic infection affecting millions of people. Although many efforts have been made for the development of immunotherapies, there is no available vaccine against this deadly infection. One major hurdle for the rational approach to develop a T. cruzi vaccine is the limited information about the proteins produced by different phylogenetic lineages, strains, and stages of the parasite. Here, we have adapted a 1D nanoHPLC system to perform online 2D LC-MS/MS, using the autosampler to inject the eluting salt solutions in the first dimension separation. The application of this methodology for the proteomic analysis of the infective trypomastigote stage of T. cruzi led to the identification of 1448 nonredundant proteins. Furthermore, about 14% of the identified sequences comprise surface proteins, most of them glycosylphosphatidylinositol (GPI)-anchored and related to parasite pathogenesis. Immunoinformatic analysis revealed thousands of potential peptides with predicted high-binding affinity for major histocompatibility complex (MHC) class I and II molecules. The high diversity of proteins expressed on the trypomastigote surface may have many implications for host-cell invasion and immunoevasion mechanisms triggered by the parasite. Finally, we performed a rational approach to filter potential T-cell epitopes that could be further tested and validated for development of a Chagas disease vaccine.

Structural shifts of aldehyde dehydrogenase enzymes were instrumental for the early evolution of retinoid-dependent axial patterning in metazoans

Aldehyde dehydrogenases (ALDHs) catabolize toxic aldehydes and process the vitamin A-derived retinaldehyde into retinoic acid (RA), a small diffusible molecule and a pivotal chordate morphogen. In this study, we combine phylogenetic, structural, genomic, and developmental gene expression analyses to examine the evolutionary origins of ALDH substrate preference. Structural modeling reveals that processing of small aldehydes, such as acetaldehyde, by ALDH2, versus large aldehydes, including retinaldehyde, by ALDH1A is associated with small versus large substrate entry channels (SECs), respectively. Moreover, we show that metazoan ALDH1s and ALDH2s are members of a single ALDH1/2 clade and that during evolution, eukaryote ALDH1/2s often switched between large and small SECs after gene duplication, transforming constricted channels into wide opened ones and vice versa. Ancestral sequence reconstructions suggest that during the evolutionary emergence of RA signaling, the ancestral, narrow-channeled metazoan ALDH1/2 gave rise to large ALDH1 channels capable of accommodating bulky aldehydes, such as retinaldehyde, supporting the view that retinoid-dependent signaling arose from ancestral cellular detoxification mechanisms. Our analyses also indicate that, on a more restricted evolutionary scale, ALDH1 duplicates from invertebrate chordates (amphioxus and ascidian tunicates) underwent switches to smaller and narrower SECs. When combined with alterations in gene expression, these switches led to neofunctionalization from ALDH1-like roles in embryonic patterning to systemic, ALDH2-like roles, suggesting functional shifts from signaling to detoxification.

Mutations in the human phospholamban gene in patients with heart failure

BACKGROUND Phospholamban (PLN) is a crucial Ca(2+) cycling protein and a primary mediator of the β-adrenergic effects resulting in enhanced cardiac output. Mutations in the gene encoding PLN have been associated with idiopathic dilated cardiomyopathy; however, no systematic search for PLN mutations in heart failure has been conducted. METHODS We screened a cohort of 1,014 Brazilian patients with heart failure for mutations in the PLN gene. Molecular modeling studies of the mutations found were developed. Different disease etiologies were present in our sample: idiopathic, ischemic, Chagas, valvular, hypertensive, and others. RESULTS We identified 4 unrelated patients with PLN mutations (prevalence of 0.4%), 3 of them in the same amino acid residue (R9). Two patients presented a G-T missense mutation at the G26 nucleotide, which encodes an Arg-Leu substitution at codon 9 (R9L). One patient presented a G-A missense mutation at the same nucleotide, which encodes an Arg-His substitution at codon 9 (R9H). The fourth affected patient presented a T-G nonsense mutation at the nucleotide 116, substituting a termination codon for Leu-39 (L39stop). Molecular modeling studies suggested that R9L and R9H mutations might affect the region involved in protein kinase A docking and probably affect the mechanism modulating the release of phosphorylated PLN from the substrate binding site of protein kinase A. CONCLUSIONS Mutations in the PLN gene are a rare cause of heart failure, present almost exclusively in patients with dilated cardiomyopathy etiology. The Arg9 and Leu39 residues are the leading location of mutations described at this locus to date. Despite the few mutated residues described to date, the clinical spectrum of presentation appears to vary considerably.

Hydrocephalus and arthrogryposis in an immunocompetent mouse model of ZIKA teratogeny: a developmental study

The teratogenic mechanisms triggered by ZIKV are still obscure due to the lack of a suitable animal model. Here we present a mouse model of developmental disruption induced by ZIKV hematogenic infection. The model utilizes immunocompetent animals from wild-type FVB/NJ and C57BL/6J strains, providing a better analogy to the human condition than approaches involving immunodeficient, genetically modified animals, or direct ZIKV injection into the brain. When injected via the jugular vein into the blood of pregnant females harboring conceptuses from early gastrulation to organogenesis stages, akin to the human second and fifth week of pregnancy, ZIKV infects maternal tissues, placentas and embryos/fetuses. Early exposure to ZIKV at developmental day 5 (second week in humans) produced complex manifestations of anterior and posterior dysraphia and hydrocephalus, as well as severe malformations and delayed development in 10.5 days post-coitum (dpc) embryos. Exposure to the virus at 7.5–9.5 dpc induces intra-amniotic hemorrhage, widespread edema, and vascular rarefaction, often prominent in the cephalic region. At these stages, most affected embryos/fetuses displayed gross malformations and/or intrauterine growth restriction (IUGR), rather than isolated microcephaly. Disrupted conceptuses failed to achieve normal developmental landmarks and died in utero. Importantly, this is the only model so far to display dysraphia and hydrocephalus, the harbinger of microcephaly in humans, as well as arthrogryposis, a set of abnormal joint postures observed in the human setting. Late exposure to ZIKV at 12.5 dpc failed to produce noticeable malformations. We have thus characterized a developmental window of opportunity for ZIKV-induced teratogenesis encompassing early gastrulation, neurulation and early organogenesis stages. This should not, however, be interpreted as evidence for any safe developmental windows for ZIKV exposure. Late developmental abnormalities correlated with damage to the placenta, particularly to the labyrinthine layer, suggesting that circulatory changes are integral to the altered phenotypes.

Active Glutaminase C Self-assembles into a Supratetrameric Oligomer That Can Be Disrupted by an Allosteric Inhibitor

Background: GAC supplies for increased metabolic needs of tumors because of exclusive localization and kinetic properties. Results: Higher than tetramer oligomers are the active form in in vitro and in cellular assays. Bis-2-(5-phenylacetamido-1,3,4-thiadiazol-2-yl)ethyl sulfide disrupts oligomers. Conclusion: A novel molecular mechanism for GAC activation is proposed. Significance: The data affect the development of therapies targeting GAC in tumors, with emphasis on allosteric inhibitors. The phosphate-dependent transition between enzymatically inert dimers into catalytically capable tetramers has long been the accepted mechanism for the glutaminase activation. Here, we demonstrate that activated glutaminase C (GAC) self-assembles into a helical, fiber-like double-stranded oligomer and propose a molecular model consisting of seven tetramer copies per turn per strand interacting via the N-terminal domains. The loop 321LRFNKL326 is projected as the major regulating element for self-assembly and enzyme activation. Furthermore, the previously identified in vivo lysine acetylation (Lys311 in humans, Lys316 in mouse) is here proposed as an important down-regulator of superoligomer assembly and protein activation. Bis-2-(5-phenylacetamido-1,3,4-thiadiazol-2-yl)ethyl sulfide, a known glutaminase inhibitor, completely disrupted the higher order oligomer, explaining its allosteric mechanism of inhibition via tetramer stabilization. A direct correlation between the tendency to self-assemble and the activity levels of the three mammalian glutaminase isozymes was established, with GAC being the most active enzyme while forming the longest structures. Lastly, the ectopic expression of a fiber-prone superactive GAC mutant in MDA-MB 231 cancer cells provided considerable proliferative advantages to transformed cells. These findings yield unique implications for the development of GAC-oriented therapeutics targeting tumor metabolism.

Structural Insights into Human Peroxisome Proliferator Activated Receptor Delta (PPAR-Delta) Selective Ligand Binding

Peroxisome proliferator activated receptors (PPARs δ, α and γ) are closely related transcription factors that exert distinct effects on fatty acid and glucose metabolism, cardiac disease, inflammatory response and other processes. Several groups developed PPAR subtype specific modulators to trigger desirable effects of particular PPARs without harmful side effects associated with activation of other subtypes. Presently, however, many compounds that bind to one of the PPARs cross-react with others and rational strategies to obtain highly selective PPAR modulators are far from clear. GW0742 is a synthetic ligand that binds PPARδ more than 300-fold more tightly than PPARα or PPARγ but the structural basis of PPARδ:GW0742 interactions and reasons for strong selectivity are not clear. Here we report the crystal structure of the PPARδ:GW0742 complex. Comparisons of the PPARδ:GW0742 complex with published structures of PPARs in complex with α and γ selective agonists and pan agonists suggests that two residues (Val312 and Ile328) in the buried hormone binding pocket play special roles in PPARδ selective binding and experimental and computational analysis of effects of mutations in these residues confirms this and suggests that bulky substituents that line the PPARα and γ ligand binding pockets as structural barriers for GW0742 binding. This analysis suggests general strategies for selective PPARδ ligand design.

Identification of Novel Interaction between ADAM17 (a Disintegrin and Metalloprotease 17) and Thioredoxin-1

Background: The identification of potential interaction partners for TACE could be instrumental in understanding the regulation of TACE activity. Results: Trx-1 interacts with the cytoplasmic domain of ADAM17. Conclusion: Trx-1 regulates ADAM17 activity. Significance: The data suggest a negative ADAM17 regulation in the HB-EGF shedding model. ADAM17, which is also known as TNFα-converting enzyme, is the major sheddase for the EGF receptor ligands and is considered to be one of the main proteases responsible for the ectodomain shedding of surface proteins. How a membrane-anchored proteinase with an extracellular catalytic domain can be activated by inside-out regulation is not completely understood. We characterized thioredoxin-1 (Trx-1) as a partner of the ADAM17 cytoplasmic domain that could be involved in the regulation of ADAM17 activity. We induced the overexpression of the ADAM17 cytoplasmic domain in HEK293 cells, and ligands able to bind this domain were identified by MS after protein immunoprecipitation. Trx-1 was also validated as a ligand of the ADAM17 cytoplasmic domain and full-length ADAM17 recombinant proteins by immunoblotting, immunolocalization, and solid phase binding assay. In addition, using nuclear magnetic resonance, it was shown in vitro that the titration of the ADAM17 cytoplasmic domain promotes changes in the conformation of Trx-1. The MS analysis of the cross-linked complexes showed cross-linking between the two proteins by lysine residues. To further evaluate the functional role of Trx-1, we used a heparin-binding EGF shedding cell model and observed that the overexpression of Trx-1 in HEK293 cells could decrease the activity of ADAM17, activated by either phorbol 12-myristate 13-acetate or EGF. This study identifies Trx-1 as a novel interaction partner of the ADAM17 cytoplasmic domain and suggests that Trx-1 is a potential candidate that could be involved in ADAM17 activity regulation.

Protein folding creates structure-based, noncontiguous consensus phosphorylation motifs recognized by kinases

Kinases can recognize consensus motifs that are noncontiguous in folded proteins. Origami-Like Substrate Recognition Proteins fold into complex three-dimensional structures, yet most sites in proteins that are modified posttranslationally have been identified within short linear consensus motifs in the primary amino acid sequence. Duarte et al. found that kinases can recognize a consensus site that is formed by distinct noncontiguous parts of the folded substrate protein. They characterized in α-tubulin an example of what they termed a “structurally formed” consensus site, a threonine residue phosphorylated by a specific member of the protein kinase C (PKC) family, and then identified structurally formed consensus sites in other substrates of PKC and PKA (protein kinase A). Thus, researchers need to look beyond the linear sequence of the protein to its three-dimensional structure to identify all of the potential consensus phosphorylation sites in a protein. Linear consensus motifs are short contiguous sequences of residues within a protein that can form recognition modules for protein interaction or catalytic modification. Protein kinase specificity and the matching of kinases to substrates have been mostly defined by phosphorylation sites that occur in linear consensus motifs. However, phosphorylation can also occur within sequences that do not match known linear consensus motifs recognized by kinases and within flexible loops. We report the identification of Thr253 in α-tubulin as a site that is phosphorylated by protein kinase C βI (PKCβI). Thr253 is not part of a linear PKC consensus motif. Instead, Thr253 occurs within a region on the surface of α-tubulin that resembles a PKC phosphorylation site consensus motif formed by basic residues in different parts of the protein, which come together in the folded protein to form the recognition motif for PKCβI. Mutations of these basic residues decreased substrate phosphorylation, confirming the presence of this “structurally formed” consensus motif and its importance for the protein kinase–substrate interaction. Analysis of previously reported protein kinase A (PKA) and PKC substrates identified sites within structurally formed consensus motifs in many substrates of these two kinase families. Thus, the concept of consensus phosphorylation site motif needs to be expanded to include sites within these structurally formed consensus motifs.

The overexpression of a single oncogene (ERBB2/HER2) alters the proteomic landscape of extracellular vesicles

ERBB2/HER2 amplification activates signaling cascades that lead to a tumor cell phenotype. However, despite its remarkable importance in oncology, the consequences of HER2 amplification over the extracellular vesicles (EVs) content have not yet been investigated. Here, we isolated EVs secreted by HB4a, a mammary luminal epithelial cell line and C5.2, its HER2‐overexpressing clone. We isolated two EV sets (20 and 100 K) by ultracentrifugation and used electron microscopy and nanoparticle tracking analysis for their morphological characterization. We employed GeLC‐MS/MS combined with isotope‐coded protein labeling to evaluate cell‐derived proteins and LC‐MS/MS label free spectral counting to quantify the EVs proteome. We found higher HER2 levels in both C5.2‐derived EVs when compared with C5.2 cells, suggesting its preferential shuttling. Proteins capable of inducing malignant transformation are enriched in both C5.2 EV subsets, including two HER2‐related proteins involved in cell motility and invasion, cofilin and CD44. MetaCore™ analysis indicated an enrichment of cell adhesion and cytoskeleton‐remodeling pathways in C5.2 EVs, as well as proteins related to HER2 signaling, such as sphingosine‐1‐phosphate pathway. Together, our data indicate that in terms of protein content, distinct vesicle sets reinforce and complement each other. Our results also suggest that HER2‐upregulated proteins from EVs may be relevant for cellular malignancy and can be potential biomarkers for HER2+ cancer patients.