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.

Orchestration of multiscale model for computational oncology

Cancer is a challenging disease that involves multiple types of biological interactions in different time and space scales. Often computational modelling has been facing problems that, in the current technology level, is impracticable to represent in a single space-time continuum. To handle this sort of problems, complex orchestrations of multiscale models is frequently done. PRIMAGE is a large EU project that aims to support personalized childhood cancer diagnosis and prognosis. The goal is to do so predicting the growth of the solid tumour using multiscale in-silico technologies. The project proposes an open cloud-based platform to support decision making in the clinical management of paediatric cancers. The orchestration of predictive models is in general complex and would require a software framework that support and facilitate such task. The present work, proposes the development of an updated framework, referred herein as the VPH-HFv3, as a part of the PRIMAGE project. This framework, a complete re-writing with respect to the previous versions, aims to orchestrate several models, which are in concurrent development, using an architecture as simple as possible, easy to maintain and with high reusability. This sort of problem generally requires unfeasible execution times. To overcome this problem was developed a strategy of particularisation, which maps the upper-scale model results into a smaller number and homogenisation which does the inverse way and analysed the accuracy of this approach.

Electrochemical sensing strategies for the detection of interactions between biological macromolecules

Electrochemical biosensors provide an attractive means to analyze the content of a biological sample due to the direct conversion of a biological event to an electronic signal, enabling the development of cheap, small, portable and simple devices, that allow multiplex and real-time detection. At the same time nanobiotechnology is drastically revolutionizing the biosensors development and different transduction strategies exploit concepts developed in these field to simplify the analysis operations for operators and end users, offering higher specificity, higher sensitivity, higher operational stability, integrated sample treatments and shorter analysis time. The aim of this PhD work has been the application of nanobiotechnological strategies to electrochemical biosensors for the detection of biological macromolecules. Specifically, one project was focused on the application of a DNA nanotechnology called hybridization chain reaction (HCR), to amplify the hybridization signal in an electrochemical DNA biosensor. Another project on which the research activity was focused concerns the development of an electrochemical biosensor based on a biological model membrane anchored to a solid surface (tBLM), for the recognition of interactions between the lipid membrane and different types of target molecules.

Ocean-Atmosphere Interactions: Linking Oceanic Biological Activity to Marine Aerosol Physico-Chemical and Cloud Properties

In this work, in-situ measurements of aerosol chemical composition, particle number size distribution, cloud-relevant properties and ground-based cloud observations were combined with high-resolution satellite sea surface chlorophyll-a concentration and air mass back-trajectory data to investigate the impact of the marine biota on aerosol physico-chemical and cloud properties. Studies were performed over the North-Eastern Atlantic Ocean, the central Mediterranean Sea, and the Arctic Ocean, by deploying both multi-year datasets and short-time scale observations. All the data were chosen to be representative of the marine atmosphere, reducing to a minimum any anthropogenic input. A relationship between the patterns of marine biological activity and the time evolution of marine aerosol properties was observed, under a variety of aspects, from chemical composition to number concentration and size distribution, up to the most cloud‐relevant properties. At short-time scales (1-2 months), the aerosol properties tend to respond to biological activity variations with a delay of about one to three weeks. This delay should be considered in model applications that make use of Chlorophyll-a to predict marine aerosol properties at high temporal resolution. The impact of oceanic biological activity on the microphysical properties of marine stratiform clouds is also evidenced by our analysis, over the Eastern North Atlantic Ocean. Such clouds tend to have a higher number of smaller cloud droplets in periods of high biological activity with respect to quiescent periods. This confirms the possibility of feedback interactions within the biota-aerosol-cloud climate system. Achieving a better characterization of the time and space relationships linking oceanic biological activity to marine aerosol composition and properties may significantly impact our future capability of predicting the chemical composition of the marine atmosphere, potentially contributing to reducing the uncertainty of future climate predictions, through a better understanding of the natural climate system.

Molecular interactions: metal ions and protein chaperones in the urease system from Helicobacter pylori

Although nickel is a toxic metal for living organisms in its soluble form, its importance in many biological processes recently emerged. In this view, the investigation of the nickel-dependent enzymes urease and [NiFe]-hydrogenase, especially the mechanism of nickel insertion into their active sites, represent two intriguing case studies to understand other analogous systems and therefore to lead to a comprehension of the nickel trafficking inside the cell. Moreover, these two enzymes have been demonstrated to ensure survival and colonization of the human pathogen H. pylori, the only known microorganism able to proliferate in the gastric niche. The right nickel delivering into the urease active site requires the presence of at least four accessory proteins, UreD, UreE, UreF and UreG. Similarly, analogous process is principally mediated by HypA and HypB proteins in the [NiFe]-hydrogenase system. Indeed, HpHypA and HpHypB also have been proposed to act in the activation of the urease enzyme from H. pylori, probably mobilizing nickel ions from HpHypA to the HpUreE-HpUreG complex. A complete comprehension of the interaction mechanism between the accessory proteins and the crosstalk between urease and hydrogenase accessory systems requires the determination of the role of each protein chaperone that strictly depends on their structural and biochemical properties. The availability of HpUreE, HpUreG and HpHypA proteins in a pure form is a pre-requisite to perform all the subsequent protein characterizations, thus their purification was the first aim of this work. Subsequently, the structural and biochemical properties of HpUreE were investigated using multi-angle and quasi-elastic light scattering, as well as NMR and circular dichroism spectroscopy. The thermodynamic parameters of Ni2+ and Zn2+ binding to HpUreE were principally established using isothermal titration calorimetry and the importance of key histidine residues in the process of binding metal ions was studied using site-directed mutagenesis. The molecular details of the HpUreE-HpUreG and HpUreE-HpHypA protein-protein assemblies were also elucidated. The interaction between HpUreE and HpUreG was investigated using ITC and NMR spectroscopy, and the influence of Ni2+ and Zn2+ metal ions on the stabilization of this association was established using native gel electrophoresis, light scattering and thermal denaturation scanning followed by CD spectroscopy. Preliminary HpUreE-HpHypA interaction studies were conducted using ITC. Finally, the possible structural architectures of the two protein-protein assemblies were rationalized using homology modeling and docking computational approaches. All the obtained data were interpreted in order to achieve a more exhaustive picture of the urease activation process, and the correlation with the accessory system of the hydrogenase enzyme, considering the specific role and activity of the involved protein players. A possible function for Zn2+ in the chaperone network involved in Ni2+ trafficking and urease activation is also envisaged.

Multi-level models and infrastructures for simulating biological system development

The hierarchical organisation of biological systems plays a crucial role in the pattern formation of gene expression resulting from the morphogenetic processes, where autonomous internal dynamics of cells, as well as cell-to-cell interactions through membranes, are responsible for the emergent peculiar structures of the individual phenotype. Being able to reproduce the systems dynamics at different levels of such a hierarchy might be very useful for studying such a complex phenomenon of self-organisation. The idea is to model the phenomenon in terms of a large and dynamic network of compartments, where the interplay between inter-compartment and intra-compartment events determines the emergent behaviour resulting in the formation of spatial patterns. According to these premises the thesis proposes a review of the different approaches already developed in modelling developmental biology problems, as well as the main models and infrastructures available in literature for modelling biological systems, analysing their capabilities in tackling multi-compartment / multi-level models. The thesis then introduces a practical framework, MS-BioNET, for modelling and simulating these scenarios exploiting the potential of multi-level dynamics. This is based on (i) a computational model featuring networks of compartments and an enhanced model of chemical reaction addressing molecule transfer, (ii) a logic-oriented language to flexibly specify complex simulation scenarios, and (iii) a simulation engine based on the many-species/many-channels optimised version of Gillespie’s direct method. The thesis finally proposes the adoption of the agent-based model as an approach capable of capture multi-level dynamics. To overcome the problem of parameter tuning in the model, the simulators are supplied with a module for parameter optimisation. The task is defined as an optimisation problem over the parameter space in which the objective function to be minimised is the distance between the output of the simulator and a target one. The problem is tackled with a metaheuristic algorithm. As an example of application of the MS-BioNET framework and of the agent-based model, a model of the first stages of Drosophila Melanogaster development is realised. The model goal is to generate the early spatial pattern of gap gene expression. The correctness of the models is shown comparing the simulation results with real data of gene expression with spatial and temporal resolution, acquired in free on-line sources.

Synthesis of Modified Amino Acids and Insertion in Peptides and Mimetics. Structural Aspects and Impact on Biological Activity.

I studied the effects exerted by the modifications on structures and biological activities of the compounds so obtained. I prepared peptide analogues containing unusual amino acids such as halogenated, alkylated (S)- or (R)-tryptophans, useful for the synthesis of mimetics of the endogenous opioid peptide endomorphin-1, or 2-oxo-1,3-oxazolidine-4-carboxylic acids, utilized as pseudo-prolines having a clear all-trans configuration of the preceding peptide bond. The latter gave access to a series of constrained peptidomimetics with potential interest in medicinal chemistry and in the field of the foldamers. In particular, I have dedicated much efforts to the preparation of cyclopentapeptides containing D-configured, alfa-, or beta-aminoacids, and also of cyclotetrapeptides including the retro-inverso modification. The conformational analyses confirmed that these cyclic compounds can be utilized as rigid scaffolds mimicking gamma- or beta-turns, allowing to generate new molecular and 3D diversity. Much work has been dedicated to the structural analysis in solution and in the receptor-bound state, fundamental for giving a rationale to the experimentally determined bioactivity, as well as for predicting the activity of virtual compounds (in silico pre-screen). The conformational analyses in solution has been done mostly by NMR (2D gCosy, Roesy, VT, molecular dynamics, etc.). A special section is dedicated to the prediction of plausible poses of the ligands when bound to the receptors by Molecular Docking. This computational method proved to be a powerful tool for the investigation of ligand-receptor interactions, and for the design of selective agonists and antagonists. Another practical use of cyclic peptidomimetics was the synthesis and biological evaluation of cyclic analogues of endomorphin-1 lacking in a protonable amino group. The studies revealed that a inverse type II beta-turn on D-Trp-Phe constituted the bioactive conformation.

Cell Host-Microbe Interactions: Turning Pathogen Mechanisms Into Cell's Advantages

Host-Pathogen Interaction is a very vast field of biological sciences, indeed every year many un- known pathogens are uncovered leading to an exponential growth of this field. The present work lyes between its boundaries, touching different aspects of host-pathogen interaction: We have evaluate the permissiveness of Mesenchimal Stem cell (FM-MSC from now on) to all known human affecting herpesvirus. Our study demonstrate that FM-MSC are full permissive to HSV1, HSV2, HCMV and VZV. On the other hand HHV6, HHV7, EBV and HHV8 are susceptible, but failed to activate a lytic infection program. FM-MSC are pluripotent stem cell and have been studied intensely in last decade. FM-MSC are employed in some clinical applications. For this reason it is important to known the degree of susceptibility to transmittable pathogens. Our atten- tion has then moved to bacterial pathogens: we have performed a proteome-wide in silico analy- sis of Chlamydiaceae family, searching for putative Nuclear localization Signal (NLS). Chlamy- diaceae are a family of obligate intracellular parasites. It’s reasonably to think that its members could delivered to nucleus effector proteins via NLS sequences: if that were the case the identifi- cation of NLS carrying proteins could open the way to therapeutic approaches. Our results strengthen this hypothesis: we have identified 72 protein bearing NLS, and verified their func- tionality with in vivo assays. Finally we have conceived a molecular scissor, creating a fusion protein between HIV-1 IN protein and FokI catalytic domain (a deoxyexonuclease domain). Our aim is to obtain chimeric enzyme (trojIN) which selectively identify IN naturally occurring target (HIV LTR sites) and cleaves subsequently LTR carrying DNA (for example integrated HIV1 DNA). Our preliminary results are promising since we have identified trojIN mutated version capable to selectively recognize LTR carrying DNA in an in vitro experiments.

Contribution of Non-Covalent Interactions and Electronic Effects on the Conformational Landscape and Tautomeric Equilibria of Molecules and Molecular Complexes: Structural and Dynamical Data from Rotational Spectroscopy

This thesis concerns the study of complex conformational surfaces and tautomeric equilibria of molecules and molecular complexes by quantum chemical methods and rotational spectroscopy techniques. In particular, the focus of this research is on the effects of substitution and noncovalent interactions in determining the energies and geometries of different conformers, tautomers or molecular complexes. The Free-Jet Absorption Millimeter Wave spectroscopy and the Pulsed-Jet Fourier Transform Microwave spectroscopy have been applied to perform these studies and the obtained results showcase the suitability of these techniques for the study of conformational surfaces and intermolecular interactions. The series of investigations of selected medium-size molecules and complexes have shown how different instrumental setups can be used to obtain a variety of results on molecular properties. The systems studied, include molecules of biological interest such as anethole and molecules of astrophysical interest such as N-methylaminoethanol. Moreover halogenation effects have been investigated on halogen substituted tautomeric systems (5-chlorohydroxypyridine and 6-chlorohydroxypyridine), where it has shown that the position of the inserted halogen atom affects the prototropic equilibrium. As for fluorination effects, interesting results have been achieved investigating some small complexes where a molecule of water is used as a probe to reveal the changes on the electrostatic potential of different fluorinated compounds: 2-fluoropyridine, 3-fluoropyridine and penta-fluoropyridine. While in the case of the molecular complex between water and 2-fluoropyridine and 3-fluoropyridine the geometry of the complex with one water molecule is analogous to that of pyridine with the water molecule linked to the pyridine nitrogen, the case of pentafluoropyridine reveals the effect of perfluorination and the water oxygen points towards the positive center of the pyridine ring. Additional molecular adducts with a molecule of water have been analyzed (benzylamine-water and acrylic acid-water) in order to reveal the stabilizing driving forces that characterize these complexes.

Photoactive nanodevices for potential biological applications

Biological systems are complex and highly organized architectures governed by noncovalent interactions, which are responsible for molecular recognition, self-assembly, self-organization, adaptation and evolution processes. These systems provided the inspiration for the development of supramolecular chemistry, that aimed at the design of artificial multicomponent molecular assemblies, namely supramolecular systems, properly designed to perform different operations: each constituting unit performs a single act, whereas the entire supramolecular system is able to execute a more complex function, resulting from the cooperation of the constituting components. Supramolecular chemistry deals with the development of molecular systems able to mimic naturally occurring events, for example complexation and self-assembly through the establishment of noncovalent interactions. Moreover, the application of external stimuli, such as light, allows to perform these operations in a time- and space-controlled manner. These systems can interact with biological systems and, thus, can be applied for bioimaging, therapeutic and drug delivery purposes. In this work the study of biocompatible supramolecular species able to interact with light is presented. The first part deals with the photophysical, photochemical and electrochemical characterization of water-soluble blue emitting triazoloquinolinium and triazolopyridinium salts. Moreover, their interaction with DNA has been explored, in the perspective of developing water-soluble systems for bioimaging applications. In the second part, the effect exerted by the presence of azobenzene-bearing supramolecular species in liposomes, inserted both in the phospholipid bilayer and in the in the aqueous core of vesicles has been studied, in order to develop systems able to deliver small molecules and ions in a photocontrolled manner. Moreover, the versatility of azobenzene and its broad range of applications have been highlighted, since conjugated oligoazobenzene derivatives proved not to be adequate to be inserted in the phospholipid bilayer of liposomes, but their electrochemical properties made them interesting candidates as electron acceptor materials for photovoltaic applications.

The role of biological processes in affecting the dynamics and fate of microplastics in coastal marine systems

Microplastics (MP) are omnipresent contaminants in the marine environment. Ingestion of MP has been reported for a wide range of marine biota, but to what extent the uptake by organisms affects the dynamics and fate of MP in the marine system has received little attention. My thesis explored this topic by integrating laboratory tests and experiments, field quantitative surveys of MP distribution and dynamics, and the use of specialised analytical techniques such as Attenuated-Total-Reflectance- (ATR) and imaging- Fourier Transformed Infrared Spectroscopy (FTIR). I compared different methodologies to extract MP from wild invertebrate specimens, and selected the use of potassium hydroxide (KOH) as the most cost-effective approach. I used this approach to analyse the MP contamination in various invertebrate species with different ecological traits from European salt marshes. I found that 96% of the analysed specimens (330) did not contain any MP. As preliminary environmental analyses showed high levels of environmental MP contamination, I hypothesised that most MP do not accumulate into organisms but are rather fast egested. I subsequently used laboratory multi-trophic experiments and a long-term field experiment using the filter-feeding mussel Mytilus galloprovincialis and the detritus feeding polychaete Hediste diversicolor to test the aforementioned hypothesis. Overall, results showed that MP are ingested but rapidly egested by marine invertebrates, which may limit MP transfer via predator-prey interactions but at the same time enhance their transfer via detrital pathways in the sediments. These processes seem to be extremely variable over time, with potential unexplored environmental consequences. This rapid dynamics also limits the conclusions that can be derived from static observations of MP contents in marine organisms, not fully capturing the real levels of potential contaminations by marine species. This emphasises the need to consider such dynamics in future work to measure the uptake rates by organisms in natural systems.

Beyond pollinator reward: steps forward and knowledge gaps on the role of floral nectar in plant-animal interactions.

The work done within the framework of my PhD project has been carried out between November 2019 and January 2023 at the Department of Biological, Geological and Environmental Sciences of the University of Bologna, under the supervision of Prof. Marta Galloni and PhD Gherardo Bogo. A period of three months was spent at the Natural History Museum of Rijeka, under the supervision of Prof. Boštjan Surina. The main aim of the thesis was to investigate further the so-called pollinator manipulation hypothesis, which states that when a floral visitor gets in contact with a specific nectar chemistry, the latter affects its behavior of visit on flowers, with potential repercussions on the plant reproductive fitness. To the purpose, the topic was tackled by means of three main approaches: field studies, laboratory assessments, and bibliographic reviews. This research project contributes to two main aspects. First, when insects encounter nectar-like concentrations of a plethora of secondary metabolites in their food-environment, various aspects of their behavior relevant to flower visitation can be affected. In addition, the results I gained confirm that the combination of field studies and laboratory assessments allows to get more realistic pictures of a given phenomenon than the single approaches. Second, reviewing the existent literature in the field of nectar ecology has highlighted how crucial is to establish the origin of nectar biogenic amines to either confirm or reject the multiple speculations made on the role of nectar microbes in shaping plant-animal interactions.