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.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Stem Cells as a therapy for myocardial infarction in animal models

Advances in stem cell biology have challenged the notion that infarcted myocardium is irreparable. The pluripotent ability of stem cells to differentiate into specialized cell lines began to garner intense interest within cardiology when it was shown in animal models that intramyocardial injection of bone marrow stem cells (MSCs), or the mobilization of bone marrow stem cells with spontaneous homing to myocardium, could improve cardiac function and survival after induced myocardial infarction (MI) [1, 2]. Furthermore, the existence of stem cells in myocardium has been identified in animal heart [3, 4], and intense research is under way in an attempt to clarify their potential clinical application for patients with myocardial infarction. To date, in order to identify the best one, different kinds of stem cells have been studied; these have been derived from embryo or adult tissues (i.e. bone marrow, heart, peripheral blood etc.). Currently, three different biologic therapies for cardiovascular diseases are under investigation: cell therapy, gene therapy and the more recent “tissue-engineering” therapy . During my Ph.D. course, first I focalised my study on the isolation and characterization of Cardiac Stem Cells (CSCs) in wild-type and transgenic mice and for this purpose I attended, for more than one year, the Cardiovascular Research Institute of the New York Medical College, in Valhalla (NY, USA) under the direction of Doctor Piero Anversa. During this period I learnt different Immunohistochemical and Biomolecular techniques, useful for investigating the regenerative potential of stem cells. Then, during the next two years, I studied the new approach of cardiac regenerative medicine based on “tissue-engineering” in order to investigate a new strategy to regenerate the infracted myocardium. Tissue-engineering is a promising approach that makes possible the creation of new functional tissue to replace lost or failing tissue. This new discipline combines isolated functioning cells and biodegradable 3-dimensional (3D) polymeric scaffolds. The scaffold temporarily provides the biomechanical support for the cells until they produce their own extracellular matrix. Because tissue-engineering constructs contain living cells, they may have the potential for growth and cellular self-repair and remodeling. In the present study, I examined whether the tissue-engineering strategy within hyaluron-based scaffolds would result in the formation of alternative cardiac tissue that could replace the scar and improve cardiac function after MI in syngeneic heterotopic rat hearts. Rat hearts were explanted, subjected to left coronary descending artery occlusion, and then grafted into the abdomen (aorta-aorta anastomosis) of receiving syngeneic rat. After 2 weeks, a pouch of 3 mm2 was made in the thickness of the ventricular wall at the level of the post-infarction scar. The hyaluronic scaffold, previously engineered for 3 weeks with rat MSCs, was introduced into the pouch and the myocardial edges sutured with few stitches. Two weeks later we evaluated the cardiac function by M-Mode echocardiography and the myocardial morphology by microscope analysis. We chose bone marrow-derived mensenchymal stem cells (MSCs) because they have shown great signaling and regenerative properties when delivered to heart tissue following a myocardial infarction (MI). However, while the object of cell transplantation is to improve ventricular function, cardiac cell transplantation has had limited success because of poor graft viability and low cell retention, that’s why we decided to combine MSCs with a biopolimeric scaffold. At the end of the experiments we observed that the hyaluronan fibres had not been substantially degraded 2 weeks after heart-transplantation. Most MSCs had migrated to the surrounding infarcted area where they were especially found close to small-sized vessels. Scar tissue was moderated in the engrafted region and the thickness of the corresponding ventricular wall was comparable to that of the non-infarcted remote area. Also, the left ventricular shortening fraction, evaluated by M-Mode echocardiography, was found a little bit increased when compared to that measured just before construct transplantation. Therefore, this study suggests that post-infarction myocardial remodelling can be favourably affected by the grafting of MSCs delivered through a hyaluron-based scaffold

Coastal Forests of Kenya-Ecology, Biodiversity & Conservation

Aims: the broad objective of this study is to investigate the ecological, biodiversity and conservation status of the coastal forests of Kenya fragments. The specific aims of the study are: (1) to investigate current quantitative trends in plant diversity; (2) develop a spatial and standardised vegetation database for the coastal forests Kenya; (3) investigate forest structure, species diversity and composition across the forests; (4) investigate the effect of forest fragment area on plant species diversity; (5) investigate phylogenetic diversity across these coastal remnants (6) assess vulnerability and provide conservation perspectives to concrete policy issues; (7) investigate plant and butterfly diversity correlation. Methods: I performed various analytical methods including species diversity metrics; multiple regression models for species-area relationship and small island effect; non-metric multidimensional scaling; ANOSIM; PERMANOVA; multiplicative beta diversity partitioning; species accumulation curve and species indicator analysis; statistical tests, rarefaction of species richness; phylogenetic diversity metrics of Phylogenetic diversity index, mean pairwise distance, mean nearest taxon distance, and their null-models: and Co-correspondence analysis. Results: developed the first large standardised, spatial and geo-referenced vegetation database for coastal forests of Kenya consisting of 600 plant species, across 25 forest fragments using 158 plots subdivided into 3160 subplots, 18 sacred forests and seven forest reserves; species diversity, composition and forest structure was significantly different across forest sites and between forest reserves and sacred forests, higher beta diversity, species-area relationship explained significant variability of plant diversity, small Island effect was not evident; sacred forests exhibited higher phylogenetic diversity compared to forest reserves; the threatened Red List species contributed higher evolutionary history; a strong correlation between plants and butterfly diversity. Conclusions: This study provides for the first time a standardized and large vegetation data. Results emphasizes need to improve sacred forests protection status and enhance forest connectivity across forest reserves and sacred forests.

Trajectory planning of single and dual-arm robots for time-optimal handling of liquids and objects

This Thesis studies the optimal control problem of single-arm and dual-arm serial robots to achieve the time-optimal handling of liquids and objects. The first topic deals with the planning of time-optimal anti-sloshing trajectories of an industrial robot carrying a cylindrical container filled with a liquid, considering 1-dimensional and 2-dimensional planar motions. A technique for the estimation of the sloshing height is presented, together with its extension to 3-dimensional motions. An experimental validation campaign is provided and discussed to assess the thoroughness of such a technique. As far as anti-sloshing trajectories are concerned, 2-dimensional paths are considered and, for each one of them, three constrained optimizations with different values of the sloshing-height thresholds are solved. Experimental results are presented to compare optimized and non-optimized motions. The second part focuses on the time-optimal trajectory planning for dual-arm object handling, employing two collaborative robots (cobots) and adopting an admittance-control strategy. The chosen manipulation approach, known as cooperative grasping, is based on unilateral contact between the cobots and the object, and it may lead to slipping during motion if an internal prestress along the contact-normal direction is not prescribed. Thus, a virtual penetration is considered, aimed at generating the necessary internal prestress. The stability of cooperative grasping is ensured as long as the exerted forces on the object remain inside the static-friction cone. Constrained-optimization problems are solved for 3-dimensional paths: the virtual penetration is chosen among the control inputs of the problem and friction-cone conditions are treated as inequality constraints. Also in this case experiments are presented in order to prove evidence of the firm handling of the object, even for fast motions.