Use of genetic algorithms and solvation potential to study peptides structure

In this work, genetic algorithms concepts along with a rotamer library for proteins side chains and implicit solvation potential are used to optimize the tertiary structure of peptides. We starting from the known PDB structure of its backbone which is kept fixed while the side chains allowed adopting the conformations present in the rotamer library. It was used rotamer library independent of backbone and a implicit solvation potential. The structure of Mastoporan-X was predicted using several force fields with a growing complexity; we started it with a field where the only present interaction was Lennard-Jones. We added the Coulombian term and we considered the solvation effects through a term proportional to the solvent accessible area. This paper present good and interesting results obtained using the potential with solvation term and rotamer library. Hence, the algorithm (called YODA) presented here can be a good tool to the prediction problem. (c) 2007 Elsevier B.V. All rights reserved.

Predicting peptides structure with solvation potential and rotamer library dependent of the backbone

In this work, genetic algorithms concepts along with a rotamer library dependent of backbone and implicit solvation potential are used to study the tertiary structure of peptides. We starting from known primary sequence and optimize the structure of the backbone while the side chains allowed adopting the conformations present in a rotamer library. The GA, implemented with two force fields with a growing complexity, was used predict the structure of a polyalanine and a poly-isolueucine. This paper presents good and interesting results about the study of peptides structures and about the development of computational tools to study peptides structures. (C) 2007 Elsevier B.V. All rights reserved.

DBMODELING - A database applied to the study of protein targets from genome projects

Genome sequencing efforts are providing us with complete genetic blueprints for hundreds of organisms. We are now faced with assigning, understanding, and modifying the functions of proteins encoded by these genomes. DBMODELING is a relational database of annotated comparative protein structure models and their metabolic pathway characterization, when identified. This procedure was applied to complete genomes such as Mycobacteritum tuberculosis and Xylella fastidiosa. The main interest in the study of metabolic pathways is that some of these pathways are not present in humans, which makes them selective targets for drug design, decreasing the impact of drugs in humans. In the database, there are currently 1116 proteins from two genomes. It can be accessed by any researcher at http://www.biocristalografia.df.ibilce.unesp.br/tools/. This project confirms that homology modeling is a useful tool in structural bioinformatics and that it can be very valuable in annotating genome sequence information, contributing to structural and functional genomics, and analyzing protein-ligand docking.

Kaposi's sarcoma-associated herpesvirus latency-associated nuclear antigen 1 mimics Epstein-Barr virus EBNA1 immune evasion through central repeat domain effects on protein processing

Kaposi's sarcoma-associated herpesvirus (KSHV/human herpesvirus 8 [HHV8]) and Epstein-Barr virus (EBV/HHV4) are distantly related gammaherpesviruses causing tumors in humans. KSHV latency-associated nuclear antigen 1 (LANA1) is functionally similar to the EBV nuclear antigen-1 (EBNA1) protein expressed during viral latency, although they have no amino acid similarities. EBNA1 escapes cytotoxic lymphocyte (CTL) antigen processing by inhibiting its own proteosomal degradation and retarding its own synthesis to reduce defective ribosomal product processing. We show here that the LANA1 QED-rich central repeat (CR) region, particularly the CR2CR3 subdomain, also retards LANA1 synthesis and markedly enhances LANA1 stability in vitro and in vivo. LANA1 isoforms have half-lives greater than 24 h, and fusion of the LANA1 CR2CR3 domain to a destabilized heterologous protein markedly decreases protein turnover. Unlike EBNA1, the LANA1 CR2CR3 subdomain retards translation regardless of whether it is fused to the 5′ or 3′ end of a heterologous gene construct. Manipulation of sequence order, orientation, and composition of the CR2 and CR3 subdomains suggests that specific peptide sequences rather than RNA structures are responsible for synthesis retardation. Although mechanistic differences exist between LANA1 and EBNA1, the primary structures of both proteins have evolved to minimize provoking CTL immune responses. Simple strategies to eliminate these viral inhibitory regions may markedly improve vaccine effectiveness by maximizing CTL responses. Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Pyrearinus termitilluminans larval click beetle luciferase: Active site properties, structure and function relationships and comparison with other beetle luciferases

Several beetle luciferases have been cloned and sequenced. However, most studies on structure and function relationships and bioanalytical applications were done with firefly luciferases, which are pH sensitive. Several years ago we cloned Pyrearinus termitilluminans larval click beetle luciferase, which displays the most blue-shifted bioluminescence among beetle luciferases and is pH insensitive. This enzyme was expressed in E. coli, purified, and its properties investigated. This luciferase shows slower luminescence kinetics, KM values comparable to other beetle luciferases and high catalytic constant. Fluorescence studies with 8-anilino-1-naphtalene-sulfonic acid (1,8-ANS) and modeling studies suggest that the luciferin binding site of this luciferase is very hydrophobic, supporting the solvent and orientation polarizability effects as determining mechanisms for bioluminescence colors. Although pH insensitive in the range between pH 6-8, at pH 10 this luciferase displays a remarkable red-shift and broadening of the bioluminescence spectrum. Modeling studies suggest that the residue C312 may play an important role in bioluminescence color modulation. Compared to other beetle luciferases, Pyrearinus termitilluminans luciferase also displays higher thermostability and sustained luminescence in a bacterial cell environment, which makes this luciferase particularly suitable for in vivo cell analysis and bioimaging. © The Royal Society of Chemistry and Owner Societies 2009.

Chemometric analysis of Hymenoptera toxins and defensins: A model for predicting the biological activity of novel peptides from venoms and hemolymph

When searching for prospective novel peptides, it is difficult to determine the biological activity of a peptide based only on its sequence. The trial and error approach is generally laborious, expensive and time consuming due to the large number of different experimental setups required to cover a reasonable number of biological assays. To simulate a virtual model for Hymenoptera insects, 166 peptides were selected from the venoms and hemolymphs of wasps, bees and ants and applied to a mathematical model of multivariate analysis, with nine different chemometric components: GRAVY, aliphaticity index, number of disulfide bonds, total residues, net charge, pI value, Boman index, percentage of alpha helix, and flexibility prediction. Principal component analysis (PCA) with non-linear iterative projections by alternating least-squares (NIPALS) algorithm was performed, without including any information about the biological activity of the peptides. This analysis permitted the grouping of peptides in a way that strongly correlated to the biological function of the peptides. Six different groupings were observed, which seemed to correspond to the following groups: chemotactic peptides, mastoparans, tachykinins, kinins, antibiotic peptides, and a group of long peptides with one or two disulfide bonds and with biological activities that are not yet clearly defined. The partial overlap between the mastoparans group and the chemotactic peptides, tachykinins, kinins and antibiotic peptides in the PCA score plot may be used to explain the frequent reports in the literature about the multifunctionality of some of these peptides. The mathematical model used in the present investigation can be used to predict the biological activities of novel peptides in this system, and it may also be easily applied to other biological systems. © 2011 Elsevier Inc.

Structural and Phylogenetic Studies with MjTX-I Reveal a Multi-Oligomeric Toxin - a Novel Feature in Lys49-PLA2s Protein Class

The mortality caused by snakebites is more damaging than many tropical diseases, such as dengue haemorrhagic fever, cholera, leishmaniasis, schistosomiasis and Chagas disease. For this reason, snakebite envenoming adversely affects health services of tropical and subtropical countries and is recognized as a neglected disease by the World Health Organization. One of the main components of snake venoms is the Lys49-phospholipases A2, which is catalytically inactive but possesses other toxic and pharmacological activities. Preliminary studies with MjTX-I from Bothrops moojeni snake venom revealed intriguing new structural and functional characteristics compared to other bothropic Lys49-PLA2s. We present in this article a comprehensive study with MjTX-I using several techniques, including crystallography, small angle X-ray scattering, analytical size-exclusion chromatography, dynamic light scattering, myographic studies, bioinformatics and molecular phylogenetic analyses.Based in all these experiments we demonstrated that MjTX-I is probably a unique Lys49-PLA2, which may adopt different oligomeric forms depending on the physical-chemical environment. Furthermore, we showed that its myotoxic activity is dramatically low compared to other Lys49-PLA2s, probably due to the novel oligomeric conformations and important mutations in the C-terminal region of the protein. The phylogenetic analysis also showed that this toxin is clearly distinct from other bothropic Lys49-PLA2s, in conformity with the peculiar oligomeric characteristics of MjTX-I and possible emergence of new functionalities inresponse to environmental changes and adaptation to new preys. © 2013 Salvador et al.

Dimerization of aurein 1.2: Effects in structure, antimicrobial activity and aggregation of Cândida albicans cells

Antimicrobial peptides (AMPs) are a promising solution to face the antibiotic-resistant problem because they display little or no resistance effects. Dimeric analogues of select AMPs have shown pharmacotechnical advantages, making these molecules promising candidates for the development of novel antibiotic agents. Here, we evaluate the effects of dimerization on the structure and biological activity of the AMP aurein 1.2 (AU). AU and the C- and N-terminal dimers, (AU)2K and E(AU)2, respectively, were synthesized by solid-phase peptide synthesis. Circular dichroism spectra indicated that E(AU)2 has a coiled coil structure in water while (AU)2K has an α-helix structure. In contrast, AU displayed typical spectra for disordered structures. In LPC micelles, all peptides acquired a high amount of α-helix structure. Hemolytic and vesicle permeabilization assays showed that AU has a concentration dependence activity, while this effect was less pronounced for dimeric versions, suggesting that dimerization may change the mechanism of action of AU. Notably, the antimicrobial activity against bacteria and yeast decreased with dimerization. However, dimeric peptides promoted the aggregation of C. albicans. The ability to aggregate yeast cells makes dimeric versions of AU attractive candidates to inhibit the adhesion of C. albicans to biological targets and medical devices, preventing disease caused by this fungus. © 2013 Springer-Verlag Wien.

Ineraction Model Between HRSV G-Protein and Flavonoids

Background: Acute respiratory infections (ARI) are the leading cause of infant mortality in the world, and human respiratory syncytial virus (HRSV) is one of the main agents of ARI. One of the key targets of the adaptive host immune response is the RSV G-protein, which is responsible for attachment to the host cell. There is evidence that compounds such as flavonoids can inhibit viral infection in vitro. With this in mind, the main purpose of this study was to determine, using computational tools, the potential sites for interactions between G-protein and flavonoids. Results: Our study allowed the recognition of an hRSV G-protein model, as well as a model of the interaction with flavonoids. These models were composed, mainly, of -helix and random coil proteins. The docking process showed that molecular interactions are likely to occur. The flavonoid kaempferol-3-O-α-L-arabinopyranosil-(2 → 1)-α-L-apiofuranoside-7-O-α-L-rhamnopyranoside was selected as a candidate inhibitor. The main forces of the interaction were hydrophobic, hydrogen and electrostatic. Conclusions: The model of G-protein is consistent with literature expectations, since it was mostly composed of random coils (highly glycosylated sites) and -helices (lipid regions), which are common in transmembrane proteins. The docking analysis showed that flavonoids interact with G-protein in an important ectodomain region, addressing experimental studies to these sites. The determination of the G-protein structure is of great importance to elucidate the mechanism of viral infectivity, and the results obtained in this study will allow us to propose mechanisms of cellular recognition and to coordinate further experimental studies in order to discover effective inhibitors of attachment proteins.

Estudo da Estrutura da Proteína SH do Vírus Sincicial Respiratório Humano: análise funcional da estrutura pentamérica por ferramentas de bioinformática

Pós-graduação em Biofísica Molecular - IBILCE

Recombinant hepatitis C virus-envelope protein 2 interactions with low-density lipoprotein/CD81 receptors

Hepatitis C virus (HCV) envelope protein 2 (E2) is involved in viral binding to host cells. The aim of this work was to produce recombinant E2B and E2Y HCV proteins in Escherichia coli and Pichia pastoris, respectively, and to study their interactions with low-density lipoprotein receptor (LDLr) and CD81 in human umbilical vein endothelial cells (HUVEC) and the ECV304 bladder carcinoma cell line. To investigate the effects of human LDL and differences in protein structure (glycosylated or not) on binding efficiency, the recombinant proteins were either associated or not associated with lipoproteins before being assayed. The immunoreactivity of the recombinant proteins was analysed using pooled serum samples that were either positive or negative for hepatitis C. The cells were immunophenotyped by LDLr and CD81 using flow cytometry. Binding and binding inhibition assays were performed in the presence of LDL, foetal bovine serum (FCS) and specific antibodies. The results revealed that binding was reduced in the absence of FCS, but that the addition of human LDL rescued and increased binding capacity. In HUVEC cells, the use of antibodies to block LDLr led to a significant reduction in the binding of E2B and E2Y. CD81 antibodies did not affect E2B and E2Y binding. In ECV304 cells, blocking LDLr and CD81 produced similar effects, but they were not as marked as those that were observed in HUVEC cells. In conclusion, recombinant HCV E2 is dependent on LDL for its ability to bind to LDLr in HUVEC and ECV304 cells. These findings are relevant because E2 acts to anchor HCV to host cells; therefore, high blood levels of LDL could enhance viral infectivity in chronic hepatitis C patients.

The subdomain structure of human serum albumin in solution under different pH conditions studied by small angle X-ray scattering

Small-angle X-ray scattering (SAXS) was used to study structural characteristics of human serum albumin (HSA) in solution under different pH conditions. Guinier analysis of SAXS results yielded values of the molecular radius of gyration ranging from 26.7 Å to 34.5 Å for pH varying from 2.5 to 7.0. This suggests the existence of significant differences in the overall shape of the molecule at different pH. Molecular models based on subdomains with different spatial configurations were proposed. The distance distribution functions associated with these models were calculated and compared with those determined from the experimental SAXS intensity functions. The conclusion of this SAXS study is that the arrangement of molecular subdomains is clearly pH dependent; the molecule adopting more or less compact configuration for different pH conditions. The conclusions of this systematic study on the modification in molecular shape of HSA as a response to pH changes is consistent with those of previous investigations performed for particular pH conditions.

Extensive structural change of the envelope protein of dengue virus induced by a tuned ionic strength: conformational and energetic analyses

The Dengue has become a global public health threat, with over 100 million infections annually; to date there is no specific vaccine or any antiviral drug. The structures of the envelope (E) proteins of the four known serotype of the dengue virus (DENV) are already known, but there are insufficient molecular details of their structural behavior in solution in the distinct environmental conditions in which the DENVs are submitted, from the digestive tract of the mosquito up to its replication inside the host cell. Such detailed knowledge becomes important because of the multifunctional character of the E protein: it mediates the early events in cell entry, via receptor endocytosis and, as a class II protein, participates determinately in the process of membrane fusion. The proposed infection mechanism asserts that once in the endosome, at low pH, the E homodimers dissociate and insert into the endosomal lipid membrane, after an extensive conformational change, mainly on the relative arrangement of its three domains. In this work we employ all-atom explicit solvent Molecular Dynamics simulations to specify the thermodynamic conditions in that the E proteins are induced to experience extensive structural changes, such as during the process of reducing pH. We study the structural behavior of the E protein monomer at acid pH solution of distinct ionic strength. Extensive simulations are carried out with all the histidine residues in its full protonated form at four distinct ionic strengths. The results are analyzed in detail from structural and energetic perspectives, and the virtual protein movements are described by means of the principal component analyses. As the main result, we found that at acid pH and physiological ionic strength, the E protein suffers a major structural change; for lower or higher ionic strengths, the crystal structure is essentially maintained along of all extensive simulations. On the other hand, at basic pH, when all histidine residues are in the unprotonated form, the protein structure is very stable for ionic strengths ranging from 0 to 225 mM. Therefore, our findings support the hypothesis that the histidines constitute the hot points that induce configurational changes of E protein in acid pH, and give extra motivation to the development of new ideas for antivirus compound design.