Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

The photosynthetic bacterial reaction center in native and artificial envirnoments: effects on light-induced electron transfer

In this thesis we focussed on the characterization of the reaction center (RC) protein purified from the photosynthetic bacterium Rhodobacter sphaeroides. In particular, we discussed the effects of native and artificial environment on the light-induced electron transfer processes. The native environment consist of the inner antenna LH1 complex that copurifies with the RC forming the so called core complex, and the lipid phase tightly associated with it. In parallel, we analyzed the role of saccharidic glassy matrices on the interplay between electron transfer processes and internal protein dynamics. As a different artificial matrix, we incorporated the RC protein in a layer-by-layer structure with a twofold aim: to check the behaviour of the protein in such an unusual environment and to test the response of the system to herbicides. By examining the RC in its native environment, we found that the light-induced charge separated state P+QB - is markedly stabilized (by about 40 meV) in the core complex as compared to the RC-only system over a physiological pH range. We also verified that, as compared to the average composition of the membrane, the core complex copurifies with a tightly bound lipid complement of about 90 phospholipid molecules per RC, which is strongly enriched in cardiolipin. In parallel, a large ubiquinone pool was found in association with the core complex, giving rise to a quinone concentration about ten times larger than the average one in the membrane. Moreover, this quinone pool is fully functional, i.e. it is promptly available at the QB site during multiple turnover excitation of the RC. The latter two observations suggest important heterogeneities and anisotropies in the native membranes which can in principle account for the stabilization of the charge separated state in the core complex. The thermodynamic and kinetic parameters obtained in the RC-LH1 complex are very close to those measured in intact membranes, indicating that the electron transfer properties of the RC in vivo are essentially determined by its local environment. The studies performed by incorporating the RC into saccharidic matrices evidenced the relevance of solvent-protein interactions and dynamical coupling in determining the kinetics of electron transfer processes. The usual approach when studying the interplay between internal motions and protein function consists in freezing the degrees of freedom of the protein at cryogenic temperature. We proved that the “trehalose approach” offers distinct advantages with respect to this traditional methodology. We showed, in fact, that the RC conformational dynamics, coupled to specific electron transfer processes, can be modulated by varying the hydration level of the trehalose matrix at room temperature, thus allowing to disentangle solvent from temperature effects. The comparison between different saccharidic matrices has revealed that the structural and dynamical protein-matrix coupling depends strongly upon the sugar. The analyses performed in RCs embedded in polyelectrolyte multilayers (PEM) structures have shown that the electron transfer from QA - to QB, a conformationally gated process extremely sensitive to the RC environment, can be strongly modulated by the hydration level of the matrix, confirming analogous results obtained for this electron transfer reaction in sugar matrices. We found that PEM-RCs are a very stable system, particularly suitable to study the thermodynamics and kinetics of herbicide binding to the QB site. These features make PEM-RC structures quite promising in the development of herbicide biosensors. The studies discussed in the present thesis have shown that, although the effects on electron transfer induced by the native and artificial environments tested are markedly different, they can be described on the basis of a common kinetic model which takes into account the static conformational heterogeneity of the RC and the interconversion between conformational substates. Interestingly, the same distribution of rate constants (i.e. a Gamma distribution function) can describe charge recombination processes in solutions of purified RC, in RC-LH1 complexes, in wet and dry RC-PEM structures and in glassy saccharidic matrices over a wide range of hydration levels. In conclusion, the results obtained for RCs in different physico-chemical environments emphasize the relevance of the structure/dynamics solvent/protein coupling in determining the energetics and the kinetics of electron transfer processes in a membrane protein complex.

Effetti del chitosano su composti polifenolici in colture cellulari di vite: aspetti molecolari e metabolici

Polyphenols, including flavonoids and stilbenes, are an essential part of human diet and constitute one of the most abundant and ubiquitous group of plant secondary metabolites. The level of these compounds is inducible by stress or fungal attack, so attempts are being made to identify likely biotic and abiotic elicitors and to better understand the underlying mechanism. Resveratrol (3,5,4’-trihydroxystilbene), which belongs to the stilbene family, is a naturally occurring polyphenol, found in several fruits, vegetables and beverages including red wine. It is one of the most important plant polyphenols with proved benefic activity on animal health. In the last two decades, the potential protective effects of resveratrol against cardiovascular and neurodegenerative diseases, as well as the chemopreventive properties against cancer, have been largely investigated. The most important source of polyphenols and in particular resveratrol for human diet is grape (Vitis vinifera). Since stilbenes and flavonoids play a very important role in plant defence responses and enviromental interactions, and their effects on human health seem promising, the aim of the research of this Thesis was to study at different levels the activation and the regulation of their biosynthetic pathways after chitosan treatment. Moreover, the polyphenol production in grape cells and the optimisation of cultural conditions bioreactor scale-up, were also investigated. Cell suspensions were obtained from cv. Barbera (Vitis vinifera L.) petioles and were treated with a biotic elicitor, chitosan (50 μg/mL, dissolved in acetic acid) to promote phenylpropanoid metabolism. Chitosan is a D-glucosamine polymer from fungi cell wall and therefore mimes fungal pathogen attack. Liquid cultures have been monitored for 15 days, measuring cell number, cell viability, pH and grams of fresh weight. The endogenous and released amounts of 7 stilbenes (trans and cis isomers of resveratrol, piceid and resveratroloside, and piceatannol), gallic acid, 6 hydroxycinnamic acids (trans-cinnamic, p-coumaric, caffeic, ferulic, sinapic and chlorogenic acids), 5 catechines (catechin, epicatechin, epigallocatechin-gallate (EGCG), epigallocatechin and epicatechin-gallate) and other 5 flavonoids (chalcon, naringenin, kaempferol, quercetin and rutin) in cells and cultural medium, were measured by HPLC-DAD analysis and total anthocyanins were quantified by spectrophotometric analysis. Chitosan was effective in stimulating trans-resveratrol endogenous accumulation with a sharp peak at day 4 (exceeding acetic acid and water controls by 36% and 63%, respectively), while it did not influence the production of the cis-isomer. Compared to both water and acetic acid controls, chitosan decreased the release of both trans- and cis-resveratrol respect to controls. No effect was shown on the accumulation of single resveratrol mono-glucoside isomers, but considering their total amount, normalized for the relative water control, it was possible to evidence an increase in both accumulation and release of those compounds, in chitosan-treated cells, throughout the culture period and particularly during the second week. Many of the analysed flavonoids and hydroxycinnamic acids were not present or detectable in trace amounts. Catechin, epicatechin and epigallocatechin-gallate (EGCG) were detectable both inside the cells and in the culture media, but chitosan did not affect their amounts. On the contrary, total anthocyanins have been stimulated by chitosan and their level, from day 4 to 14, was about 2-fold higher than in both controls, confirming macroscopic observations that treated suspensions showed an intense brown-red color, from day 3 onwards. These elicitation results suggest that chitosan selectively up-regulates specific biosynthetic pathways, without modifying the general accumulation pattern of other flavonoids. Proteins have been extracted from cells at day 4 of culture (corresponding to the production peak of trans-resveratrol) and separated by bidimensional electrophoresis. The 73 proteins that showed a consistently changed amount between untreated, chitosan and acetic acid (chitosan solvent) treated cells, have been identified by mass spectrometry. Chitosan induced an increase in stilbene synthase (STS, the resveratrol biosynthetic enzyme), chalcone-flavanone isomerase (CHI, that switches the pathway from chalcones to flavones and anthocyanins), pathogenesis-related proteins 10 (PRs10, a large family of defence proteins), and a decrease in many proteins belonging to primary metabolisms. A train of six distinct spots of STS encoded by the same gene and increased by chitosan, was detected on the 2-D gels, and related to the different phosphorylation degree of STS spots. Northern blot analyses have been performed on RNA extracted from cells treated with chitosan and relative controls, using probes for STS, PAL (phenylalanine ammonia lyase, the first enzyme of the biosynthetic pathway), CHS (chalcone synthase, that shares with STS the same precursors), CHI and PR-10. The up-regulation of PAL, CHS and CHI transcript expression levels correlated with the accumulation of anthocyanins. The strong increase of different molecular weight PR-10 mRNAs, correlated with the 11 PR-10 protein spots identified in proteomic analyses. The sudden decrease in trans-resveratrol endogenous accumulation after day 4 of culture, could be simply explained by the diminished resveratrol biosynthetic activity due to the lower amount of biosynthetic enzymes. This might be indirectly demonstrated by northern blot expression analyses, that showed lower levels of phenylalanine ammonia lyase (PAL) and stilbene synthase (STS) mRNAs starting from day 4. Other possible explanations could be a resveratrol oxidation process and/or the formation of other different mono-, di-glucosides and resveratrol oligomers such as viniferins. Immunolocalisation experiments performed on grape protoplasts and the subsequent analyses by confocal microscope, showed that STS, and therefore the resveratrol synthetic site, is mostly associated to intracellular membranes close to the cytosolic side of plasma membrane and in a smaller amount is localized in the cytosol. STS seemed not to be present inside vacuole and nucleus. There were no differences in the STS intracellular localisation between the different treatments. Since it was shown that stilbenes are largely released in the culture medium and that STS is a soluble protein, a possible interaction of STS with a plasma membrane transporter responsible for the extrusion of stilbenes in the culture medium, might be hypothesized. Proteomic analyses performed on subcellular fractions identified in the microsomial fraction 5 proteins taking part in channel complexes or associated with channels, that significantly changed their amount after chitosan treatment. In soluble and membrane fractions respectively 3 and 4 STS and 6 and 3 PR-10 have been identified. Proteomic results obtained from subcellular fractions substantially confirmed previous result obtained from total cell protein extracts and added more information about protein localisation and co-localisation. The interesting results obtained on Barbera cell cultures with the aim to increase polyphenol (especially stilbenes) production, have encouraged scale up tests in 1 litre bioreactors. The first trial fermentation was performed in parallel with a normal time-course in 20 mL flasks, showing that the scale-up (bigger volume and different conditions) process influenced in a very relevant way stilbenes production. In order to optimise culture parameters such as medium sucrose amount, fermentation length and inoculum cell concentration, few other fermentations were performed. Chitosan treatments were also performed. The modification of each parameter brought relevant variations in stilbenes and catechins levels, so that the production of a certain compound (or class of compounds) could be hypothetically promoted by modulating one or more culture parameters. For example the catechin yield could be improved by increasing sucrose content and the time of fermentation. The best results in stilbene yield were obtained in a 800 mL fermentation inoculated with 10.8 grams of cells and supplemented with chitosan. The culture was fed with MS medium added with 30 g/L sucrose, 25 μg/mL rifampicin and 50 μg/mL of chitosan, and was maintained at 24°C, stirred by marine impeller at 100 rpm and supplied of air at 0.16 L/min rate. Resveratroloside was the stilbene present in the larger amount, 3-5 times more than resveratrol. Because resveratrol glucosides are similarly active and more stable than free resveratrol, their production using a bioreactor could be a great advantage in an hypothetical industrial process. In my bioreactor tests, stilbenes were mainly released in the culture medium (60-80% of the total) and this fact could be another advantage for industrial applications, because it allows recovering the products directly from the culture medium without stopping the fermentation and/or killing the cells. In my best cultural conditions, it was possible to obtain 3.95 mg/L of stilbenes at day 4 (maximum resveratrol accumulation) and 5.13 mg/L at day 14 (maximum resveratroloside production). In conclusion, chitosan effect in inducing Vitis vinifera defense mechanisms can be related to its ability to increase the intracellular content of a large spectrum of antioxidants, and in particular of resveratrol, its derivates and anthocyanins. Its effect can be observed at transcriptional, proteomic (variation of soluble and membrane protein amounts) and metabolic (polyphenols production) level. The chitosan ability to elicit specific plant matabolisms can be useful to produce large quantities of antioxidant compounds from cell culture in bioreactor.

Machine learning methods for prediction of disulphide bonding states of cysteine residues in proteins

The goal of this thesis work is to develop a computational method based on machine learning techniques for predicting disulfide-bonding states of cysteine residues in proteins, which is a sub-problem of a bigger and yet unsolved problem of protein structure prediction. Improvement in the prediction of disulfide bonding states of cysteine residues will help in putting a constraint in the three dimensional (3D) space of the respective protein structure, and thus will eventually help in the prediction of 3D structure of proteins. Results of this work will have direct implications in site-directed mutational studies of proteins, proteins engineering and the problem of protein folding. We have used a combination of Artificial Neural Network (ANN) and Hidden Markov Model (HMM), the so-called Hidden Neural Network (HNN) as a machine learning technique to develop our prediction method. By using different global and local features of proteins (specifically profiles, parity of cysteine residues, average cysteine conservation, correlated mutation, sub-cellular localization, and signal peptide) as inputs and considering Eukaryotes and Prokaryotes separately we have reached to a remarkable accuracy of 94% on cysteine basis for both Eukaryotic and Prokaryotic datasets, and an accuracy of 90% and 93% on protein basis for Eukaryotic dataset and Prokaryotic dataset respectively. These accuracies are best so far ever reached by any existing prediction methods, and thus our prediction method has outperformed all the previously developed approaches and therefore is more reliable. Most interesting part of this thesis work is the differences in the prediction performances of Eukaryotes and Prokaryotes at the basic level of input coding when ‘profile’ information was given as input to our prediction method. And one of the reasons for this we discover is the difference in the amino acid composition of the local environment of bonded and free cysteine residues in Eukaryotes and Prokaryotes. Eukaryotic bonded cysteine examples have a ‘symmetric-cysteine-rich’ environment, where as Prokaryotic bonded examples lack it.

Molecular architecture of fur binding to iron-induced and - repressed genes in Helicobacter pylori

The ferric uptake regulator protein Fur regulates iron-dependent gene expression in bacteria. In the human pathogen Helicobacter pylori, Fur has been shown to regulate iron-induced and iron-repressed genes. Herein we investigate the molecular mechanisms that control this differential iron-responsive Fur regulation. Hydroxyl radical footprinting showed that Fur has different binding architectures, which characterize distinct operator typologies. On operators recognized with higher affinity by holo-Fur, the protein binds to a continuous AT-rich stretch of about 20 bp, displaying an extended protection pattern. This is indicative of protein wrapping around the DNA helix. DNA binding interference assays with the minor groove binding drug distamycin A, point out that the recognition of the holo-operators occurs through the minor groove of the DNA. By contrast, on the apo-operators, Fur binds primarily to thymine dimers within a newly identified TCATTn10TT consensus element, indicative of Fur binding to one side of the DNA, in the major groove of the double helix. Reconstitution of the TCATTn10TT motif within a holo-operator results in a feature binding swap from an holo-Fur- to an apo-Fur-recognized operator, affecting both affinity and binding architecture of Fur, and conferring apo-Fur repression features in vivo. Size exclusion chromatography indicated that Fur is a dimer in solution. However, in the presence of divalent metal ions the protein is able to multimerize. Accordingly, apo-Fur binds DNA as a dimer in gel shift assays, while in presence of iron, higher order complexes are formed. Stoichiometric Ferguson analysis indicates that these complexes correspond to one or two Fur tetramers, each bound to an operator element. Together these data suggest that the apo- and holo-Fur repression mechanisms apparently rely on two distinctive modes of operator-recognition, involving respectively the readout of a specific nucleotide consensus motif in the major groove for apo-operators, and the recognition of AT-rich stretches in the minor groove for holo-operators, whereas the iron-responsive binding affinity is controlled through metal-dependent shaping of the protein structure in order to match preferentially the major or the minor groove.

Investigating the role of single point mutations in the human proteome: a computational study

In the post genomic era with the massive production of biological data the understanding of factors affecting protein stability is one of the most important and challenging tasks for highlighting the role of mutations in relation to human maladies. The problem is at the basis of what is referred to as molecular medicine with the underlying idea that pathologies can be detailed at a molecular level. To this purpose scientific efforts focus on characterising mutations that hamper protein functions and by these affect biological processes at the basis of cell physiology. New techniques have been developed with the aim of detailing single nucleotide polymorphisms (SNPs) at large in all the human chromosomes and by this information in specific databases are exponentially increasing. Eventually mutations that can be found at the DNA level, when occurring in transcribed regions may then lead to mutated proteins and this can be a serious medical problem, largely affecting the phenotype. Bioinformatics tools are urgently needed to cope with the flood of genomic data stored in database and in order to analyse the role of SNPs at the protein level. In principle several experimental and theoretical observations are suggesting that protein stability in the solvent-protein space is responsible of the correct protein functioning. Then mutations that are found disease related during DNA analysis are often assumed to perturb protein stability as well. However so far no extensive analysis at the proteome level has investigated whether this is the case. Also computationally methods have been developed to infer whether a mutation is disease related and independently whether it affects protein stability. Therefore whether the perturbation of protein stability is related to what it is routinely referred to as a disease is still a big question mark. In this work we have tried for the first time to explore the relation among mutations at the protein level and their relevance to diseases with a large-scale computational study of the data from different databases. To this aim in the first part of the thesis for each mutation type we have derived two probabilistic indices (for 141 out of 150 possible SNPs): the perturbing index (Pp), which indicates the probability that a given mutation effects protein stability considering all the “in vitro” thermodynamic data available and the disease index (Pd), which indicates the probability of a mutation to be disease related, given all the mutations that have been clinically associated so far. We find with a robust statistics that the two indexes correlate with the exception of all the mutations that are somatic cancer related. By this each mutation of the 150 can be coded by two values that allow a direct comparison with data base information. Furthermore we also implement computational methods that starting from the protein structure is suited to predict the effect of a mutation on protein stability and find that overpasses a set of other predictors performing the same task. The predictor is based on support vector machines and takes as input protein tertiary structures. We show that the predicted data well correlate with the data from the databases. All our efforts therefore add to the SNP annotation process and more importantly found the relationship among protein stability perturbation and the human variome leading to the diseasome.

Identificazione di fattori trascrizionali associati al differenziamento neurale in cellule di sarcoma di Ewing: ruolo di NF-kB

Il sarcoma di Ewing (ES) è un tumore maligno pediatrico dell’apparato scheletrico; è associato a una traslocazione specifica codificante la proteina di fusione EWS-FLI1 e all’alta espressione di CD99, una glicoproteina di membrana fisiologicamente coinvolta in diversi processi biologici. EWS-FLI1 e CD99, sono riportati avere ruoli divergenti nella modulazione della malignità e del differenziamento di ES. CD99 inoltre è riportato modulare il pathway di MAPK, il quale interagendo con molteplici fattori di trascrizione partecipa a processi di proliferazione e differenziamento. In questo studio abbiamo investigato in due linee cellulari di ES silenziate per CD99 (TC-71shCD99 e IOR/CARshCD99) l’attività basale di diversi fattori trascrizionali quali: NF-kBp65, AP1, Elk-1, E2F e CREB. L’unico fattore trascrizionale statisticamente significativo è risultato essere NF-kBp65 e abbiamo valutato il suo ruolo nel differenziamento neurale di cellule di ES e la relazione con EWS-FLI1 e CD99. L’attività trascrizionale di NF-kB è stata valutata attraverso gene reporter assay in linee cellulari di ES a diversa espressione di CD99, EWS-FLI1 e NF-kB stesso. Il differenziamento neurale è stato valutato come espressione di βIII-Tubulin in immunofluorescenza. Il silenziamento di CD99 induce una down-modulazione dell’attività trascrizionale di NF-kB, mentre il knockdown di EWS-FLI1 ne induce un’aumento. Inoltre, il silenziamento di EWS-FLI1 non è in grado di contrastare la riduzione dell’attività di NF-kB osservata dopo silenziamento di CD99, suggerendo un ruolo dominante del CD99 nel signaling di NF-kB. Cellule deprivate di CD99 ma non di EWS-FLI1, mostrano un fenotipo differenziato in senso neurale, fenotipo che viene perso quando le cellule sono indotte a sovraesprimere NF-kB. Inoltre, in cellule CD99 positive, il silenziamento di NF-kB induce un leggero differenziamento neurale. In conclusione, questi dati hanno evidenziato il ruolo di NF-kB nel differenziamento di cellule di ES e che potrebbe essere un potenziale target nel ridurre la progressione di questo tumore.