Advanced graph mining methods for protein analysis

As one of the primary substances in a living organism, protein defines the character of each cell by interacting with the cellular environment to promote the cell’s growth and function [1]. Previous studies on proteomics indicate that the functions of different proteins could be assigned based upon protein structures [2,3]. The knowledge on protein structures gives us an overview of protein fold space and is helpful for the understanding of the evolutionary principles behind structure. By observing the architectures and topologies of the protein families, biological processes can be investigated more directly with much higher resolution and finer detail. For this reason, the analysis of protein, its structure and the interaction with the other materials is emerging as an important problem in bioinformatics. However, the determination of protein structures is experimentally expensive and time consuming, this makes scientists largely dependent on sequence rather than more general structure to infer the function of the protein at the present time. For this reason, data mining technology is introduced into this area to provide more efficient data processing and knowledge discovery approaches.

Unlike many data mining applications which lack available data, the protein structure determination problem and its interaction study, on the contrary, could utilize a vast amount of biologically relevant information on protein and its interaction, such as the protein data bank (PDB) [4], the structural classification of proteins (SCOP) databases [5], CATH databases [6], UniProt [7], and others. The difficulty of predicting protein structures, specially its 3D structures, and the interactions between proteins as shown in Figure 6.1, lies in the computational complexity of the data. Although a large number of approaches have been developed to determine the protein structures such as ab initio modelling [8], homology modelling [9] and threading [10], more efficient and reliable methods are still greatly needed.

In this chapter, we will introduce a state-of-the-art data mining technique, graph mining, which is good at defining and discovering interesting structural patterns in graphical data sets, and take advantage of its expressive power to study protein structures, including protein structure prediction and comparison, and protein-protein interaction (PPI). The current graph pattern mining methods will be described, and typical algorithms will be presented, together with their applications in the protein structure analysis.

The rest of the chapter is organized as follows: Section 6.2 will give a brief introduction of the fundamental knowledge of protein, the publicly accessible protein data resources and the current research status of protein analysis; in Section 6.3, we will pay attention to one of the state-of-the-art data mining methods, graph mining; then Section 6.4 surveys several existing work for protein structure analysis using advanced graph mining methods in the recent decade; finally, in Section 6.5, a conclusion with potential further work will be summarized.

Identifying foldable regions in protein sequence from the hydrophobic signal

Structural genomics initiatives aim to elucidate representative 3D structures for the majority of protein families over the next decade, but many obstacles must be overcome. The correct design of constructs is extremely important since many proteins will be too large or contain unstructured regions and will not be amenable to crystallization. It is therefore essential to identify regions in protein sequences that are likely to be suitable for structural study. Scooby-Domain is a fast and simple method to identify globular domains in protein sequences. Domains are compact units of protein structure and their correct delineation will aid structural elucidation through a divide-and-conquer approach. Scooby-Domain predictions are based on the observed lengths and hydrophobicities of domains from proteins with known tertiary structure. The prediction method employs an A*-search to identify sequence regions that form a globular structure and those that are unstructured. On a test set of 173 proteins with consensus CATH and SCOP domain definitions, Scooby-Domain has a sensitivity of 50% and an accuracy of 29%, which is better than current state-of-the-art methods. The method does not rely on homology searches and, therefore, can identify previously unknown domains.

Using genetic algorithm to design protein sequence

In this work, genetic algorithms concepts along with a rotamer library for proteins side chains are used to optimize the tertiary structure of the hydrophobic core of Cytochrome b(562) starting from the known PDB structure of its backbone which is kept fixed while the side chains of the hydrophobic core are allowed to adopt the conformations present in the rotamer library. The atoms of the side chains forming the core interact via van der Waals energy. Besides the prediction of the native core structure, it is also suggested a set of different amino acid sequences for this core. Comparison between these new cores and the native are made in terms of their volumes, van der Waals energies values and the numbers of contacts made by the side chains forming the cores. This paper proves that genetic algorithms area efficient to design new sequence for the protein core. (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.

Molecular modeling databases: A new way in the search of protein targets for drug development

DBMODELING is a relational database of annotated comparative protein structure models and their metabolic, pathway characterization. It is focused on enzymes identified in the genomes of Mycobacterium tuberculosis and Xylella fastidiosa. The main goal of the present database is to provide structural models to be used in docking simulations and drug design. However, since the accuracy of structural models is highly dependent on sequence identity between template and target, it is necessary to make clear to the user that only models which show high structural quality should be used in such efforts. Molecular modeling of these genomes generated a database, in which all structural models were built using alignments presenting more than 30% of sequence identity, generating models with medium and high accuracy. All models in the database are publicly accessible at http://www.biocristalografia.df.ibilce.unesp.br/tools. DBMODELING user interface provides users friendly menus, so that all information can be printed in one stop from any web browser. Furthermore, DBMODELING also provides a docking interface, which allows the user to carry out geometric docking simulation, against the molecular models available in the database. There are three other important homology model databases: MODBASE, SWISSMODEL, and GTOP. The main applications of these databases are described in the present article. © 2007 Bentham Science Publishers Ltd.

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.

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.

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.

Assessing the functional structure of molecular transporters by EPR spectroscopy

In this thesis, the self-assembled functional structure of a broad range of amphiphilic molecular transporters is studied. By employing paramagnetic probe molecules and ions, continuous-wave and pulse electron paramagnetic resonance spectroscopy reveal information about the local structure of these materials from the perspective of incorporated guest molecules. First, the transport function of human serum albumin for fatty acids is in the focus. As suggested by the crystal structure, the anchor points for the fatty acids are distributed asymmetrically in the protein. In contrast to the crystallographic findings, a remarkably symmetric entry point distribution of the fatty acid binding channels is found, which may facilitate the uptake and release of the guest molecules. Further, the metal binding of 1,2,3-triazole modified star-shaped cholic acid oligomers is studied. These biomimetic molecules are able to include and transport molecules in solvents of different polarity. A pre-arrangement of the triazole groups induces a strong chelate-like binding and close contact between guest molecule and metal ion. In absence of a preordering, each triazole moiety acts as a single entity and the binding affinity for metal ions is strongly decreased. Hydrogels based on N-isopropylacrylamide phase separate from water above a certain temperature. The macroscopic thermal collapse of these hydrogels is utilized as a tool for dynamic nuclear polarization. It is shown that a radical-free hyperpolarized solution can be achieved with a spin-labeled gel as separable matrix. On the nanoscale, these hydrogels form static heterogeneities in both structure and function. Collapsed regions protect the spin probes from a chemical decay while open, water-swollen regions act as catalytic centers. Similarly, thermoresponsive dendronized polymers form structural heterogeneities, which are, however, highly dynamic. At the critical temperature, they trigger the aggregation of the polymer into mesoglobules. The dehydration of these aggregates is a molecularly controlled non-equilibrium process that is facilitated by a hydrophobic dendritic core. Further, a slow heating rate results in a kinetically entrapped non-equilibrium state due to the formation of an impermeable dense polymeric layer at the periphery of the mesoglobule.

Characterization of novel StAR (steroidogenic acute regulatory protein) mutations causing non-classic lipoid adrenal hyperplasia

Context Steroidogenic acute regulatory protein (StAR) is crucial for transport of cholesterol to mitochondria where biosynthesis of steroids is initiated. Loss of StAR function causes lipoid congenital adrenal hyperplasia (LCAH). Objective StAR gene mutations causing partial loss of function manifest atypical and may be mistaken as familial glucocorticoid deficiency. Only a few mutations have been reported. Design To report clinical, biochemical, genetic, protein structure and functional data on two novel StAR mutations, and to compare them with published literature. Setting Collaboration between the University Children's Hospital Bern, Switzerland, and the CIBERER, Hospital Vall d'Hebron, Autonomous University, Barcelona, Spain. Patients Two subjects of a non-consanguineous Caucasian family were studied. The 46,XX phenotypic normal female was diagnosed with adrenal insufficiency at the age of 10 months, had normal pubertal development and still has no signs of hypergonodatropic hypogonadism at 32 years of age. Her 46,XY brother was born with normal male external genitalia and was diagnosed with adrenal insufficiency at 14 months. Puberty was normal and no signs of hypergonadotropic hypogonadism are present at 29 years of age. Results StAR gene analysis revealed two novel compound heterozygote mutations T44HfsX3 and G221S. T44HfsX3 is a loss-of-function StAR mutation. G221S retains partial activity (~30%) and is therefore responsible for a milder, non-classic phenotype. G221S is located in the cholesterol binding pocket and seems to alter binding/release of cholesterol. Conclusions StAR mutations located in the cholesterol binding pocket (V187M, R188C, R192C, G221D/S) seem to cause non-classic lipoid CAH. Accuracy of genotype-phenotype prediction by in vitro testing may vary with the assays employed.