Analyzing the dependence structure of microarray data: a copula–based approach

The main aim of this Ph.D. dissertation is the study of clustering dependent data by means of copula functions with particular emphasis on microarray data. Copula functions are a popular multivariate modeling tool in each field where the multivariate dependence is of great interest and their use in clustering has not been still investigated. The first part of this work contains the review of the literature of clustering methods, copula functions and microarray experiments. The attention focuses on the K–means (Hartigan, 1975; Hartigan and Wong, 1979), the hierarchical (Everitt, 1974) and the model–based (Fraley and Raftery, 1998, 1999, 2000, 2007) clustering techniques because their performance is compared. Then, the probabilistic interpretation of the Sklar’s theorem (Sklar’s, 1959), the estimation methods for copulas like the Inference for Margins (Joe and Xu, 1996) and the Archimedean and Elliptical copula families are presented. In the end, applications of clustering methods and copulas to the genetic and microarray experiments are highlighted. The second part contains the original contribution proposed. A simulation study is performed in order to evaluate the performance of the K–means and the hierarchical bottom–up clustering methods in identifying clusters according to the dependence structure of the data generating process. Different simulations are performed by varying different conditions (e.g., the kind of margins (distinct, overlapping and nested) and the value of the dependence parameter ) and the results are evaluated by means of different measures of performance. In light of the simulation results and of the limits of the two investigated clustering methods, a new clustering algorithm based on copula functions (‘CoClust’ in brief) is proposed. The basic idea, the iterative procedure of the CoClust and the description of the written R functions with their output are given. The CoClust algorithm is tested on simulated data (by varying the number of clusters, the copula models, the dependence parameter value and the degree of overlap of margins) and is compared with the performance of model–based clustering by using different measures of performance, like the percentage of well–identified number of clusters and the not rejection percentage of H0 on . It is shown that the CoClust algorithm allows to overcome all observed limits of the other investigated clustering techniques and is able to identify clusters according to the dependence structure of the data independently of the degree of overlap of margins and the strength of the dependence. The CoClust uses a criterion based on the maximized log–likelihood function of the copula and can virtually account for any possible dependence relationship between observations. Many peculiar characteristics are shown for the CoClust, e.g. its capability of identifying the true number of clusters and the fact that it does not require a starting classification. Finally, the CoClust algorithm is applied to the real microarray data of Hedenfalk et al. (2001) both to the gene expressions observed in three different cancer samples and to the columns (tumor samples) of the whole data matrix.

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.

Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

Identification of eight novel coagulation factor XIII subunit A mutations: implied consequences for structure and function

Severe hereditary coagulation factor XIII deficiency is a rare homozygous bleeding disorder affecting one person in every two million individuals. In contrast, heterozygous factor XIII deficiency is more common, but usually not associated with severe hemorrhage such as intracranial bleeding or hemarthrosis. In most cases, the disease is caused by F13A gene mutations. Causative mutations associated with the F13B gene are rarer.

Recent advances into understanding some aspects of the structure and function of mammalian and avian lungs

Recent findings are reported about certain aspects of the structure and function of the mammalian and avian lungs that include (a) the architecture of the air capillaries (ACs) and the blood capillaries (BCs); (b) the pulmonary blood capillary circulatory dynamics; (c) the adaptive molecular, cellular, biochemical, compositional, and developmental characteristics of the surfactant system; (d) the mechanisms of the translocation of fine and ultrafine particles across the airway epithelial barrier; and (e) the particle-cell interactions in the pulmonary airways. In the lung of the Muscovy duck Cairina moschata, at least, the ACs are rotund structures that are interconnected by narrow cylindrical sections, while the BCs comprise segments that are almost as long as they are wide. In contrast to the mammalian pulmonary BCs, which are highly compliant, those of birds practically behave like rigid tubes. Diving pressure has been a very powerful directional selection force that has influenced phenotypic changes in surfactant composition and function in lungs of marine mammals. After nanosized particulates are deposited on the respiratory tract of healthy human subjects, some reach organs such as the brain with potentially serious health implications. Finally, in the mammalian lung, dendritic cells of the pulmonary airways are powerful agents in engulfing deposited particles, and in birds, macrophages and erythrocytes are ardent phagocytizing cellular agents. The morphology of the lung that allows it to perform different functions-including gas exchange, ventilation of the lung by being compliant, defense, and secretion of important pharmacological factors-is reflected in its "compromise design."

Structure and function of CinD (YtjD) of Lactococcus lactis, a copper-induced nitroreductase involved in defense against oxidative stress

In Lactococcus lactis IL1403, 14 genes are under the control of the copper-inducible CopR repressor. This so-called CopR regulon encompasses the CopR regulator, two putative CPx-type copper ATPases, a copper chaperone, and 10 additional genes of unknown function. We addressed here the function of one of these genes, ytjD, which we renamed cinD (copper-induced nitroreductase). Copper, cadmium, and silver induced cinD in vivo, as shown by real-time quantitative PCR. A knockout mutant of cinD was more sensitive to oxidative stress exerted by 4-nitroquinoline-N-oxide and copper. Purified CinD is a flavoprotein and reduced 2,6-dichlorophenolindophenol and 4-nitroquinoline-N-oxide with k(cat) values of 27 and 11 s(-1), respectively, using NADH as a reductant. CinD also exhibited significant catalase activity in vitro. The X-ray structure of CinD was resolved at 1.35 A and resembles those of other nitroreductases. CinD is thus a nitroreductase which can protect L. lactis against oxidative stress that could be exerted by nitroaromatic compounds and copper.

Diversity of structure and function of alpha1alpha6beta3delta GABAA receptors: comparison with alpha1beta3delta and alpha6beta3delta receptors

delta subunit-containing gamma-aminobutyric acid, type A (GABA(A))receptors are expressed extrasynaptically and mediate tonic inhibition. In cerebellar granule cells, they often form receptors together with alpha(1) and/or alpha(6) subunits. We were interested in determining the architecture of receptors containing both subunits. We predefined the subunit arrangement of several different GABA(A) receptor pentamers by concatenation. These receptors composed of alpha(1), alpha(6), beta(3), and delta subunits were expressed in Xenopus oocytes. Currents elicited in response to GABA were determined in the presence and absence of 3alpha,21-dihydroxy-5alpha-pregnan-20-one (THDOC) or ethanol, or currents were elicited by 4,5,6,7-tetrahydroisoxazolo[5,4-c]-pyridin-3-ol (THIP). Several subunit configurations formed active channels. We therefore conclude that delta can assume multiple positions in a receptor pentamer made up of alpha(1), alpha(6), beta(3), and delta subunits. The different receptors differ in their functional properties. Functional expression of one receptor type was only evident in the combined presence of the neurosteroid THDOC with the channel agonist GABA. Most, but not all, receptors active with GABA/THDOC responded to THIP. None of the receptors was modulated by ethanol concentrations up to 30 mm. Several observations point to a preferred position of delta subunits between two alpha subunits in alpha(1)alpha(6)beta(3)delta receptors. This property is shared by alpha(1)beta(3)delta and alpha(6)beta(3)delta receptors, but there are differences in the additionally expressed isoforms.

Abandoned Coal Mine Drainage and Its Remediation: Impacts on Stream Ecosystem Structure and Function

The effects of abandoned mine drainage (AMD) on streams and responses to remediation efforts were studied using three streams (AMD-impacted, remediated, reference) in both the anthracite and the bituminous coal mining regions of Pennsylvania (USA). Response variables included ecosystem function as well as water chemistry and macroinvertebrate community composition. The bituminous AMD stream was extremely acidic with high dissolved metals concentrations, a prolific mid-summer growth of the filamentous alga, Mougeotia, and .10-fold more chlorophyll than the reference stream. The anthracite AMD stream had a higher pH, substrata coated with iron hydroxide(s), and negligible chlorophyll. Macroinvertebrate communities in the AMD streams were different from the reference streams, the remediated streams, and each other. Relative to the reference stream, the AMD stream(s) had (1) greater gross primary productivity (GPP) in the bituminous region and undetectable GPP in the anthracite region, (2) greater ecosystem respiration in both regions, (3) greatly reduced ammonium uptake and nitrification in both regions, (4) lower nitrate uptake in the bituminous (but not the anthracite) region, (5) more rapid phosphorus removal from the water column in both regions, (6) activities of phosphorus-acquiring, nitrogenacquiring, and hydrolytic-carbon-acquiring enzymes that indicated extreme phosphorus limitation in both regions, and (7) slower oak and maple leaf decomposition in the bituminous region and slower oak decomposition in the anthracite region. Remediation brought chlorophyll concentrations and GPP nearer to values for respective reference streams, depressed ecosystem respiration, restored ammonium uptake, and partially restored nitrification in the bituminous (but not the anthracite) region, reduced nitrate uptake to an undetectable level, restored phosphorus uptake to near normal rates, and brought enzyme activities more in line with the reference stream in the bituminous (but not the anthracite) region. Denitrification was not detected in any stream. Water chemistry and macroinvertebrate community structure analyses capture the impact of AMD at the local reach scale, but functional measures revealed that AMD has ramifications that can cascade to downstream reaches and perhaps to receiving estuaries.

Structure and function of the glucose PTS transporter from Escherichia coli

The glucose transporter IICB of the Escherichia coli phosphotransferase system (PTS) consists of a polytopic membrane domain (IIC) responsible for substrate transport and a hydrophilic C-terminal domain (IIB) responsible for substrate phosphorylation. We have overexpressed and purified a triple mutant of IIC (mut-IIC), which had recently been shown to be suitable for crystallization purposes. Mut-IIC was homodimeric as determined by blue native-PAGE and gel-filtration, and had an eyeglasses-like structure as shown by negative-stain transmission electron microscopy (TEM) and single particle analysis. Glucose binding and transport by mut-IIC, mut-IICB and wildtype-IICB were compared with scintillation proximity and in vivo transport assays. Binding was reduced and transport was impaired by the triple mutation. The scintillation proximity assay allowed determination of substrate binding, affinity and specificity of wildtype-IICB by a direct method. 2D crystallization of mut-IIC yielded highly-ordered tubular crystals and made possible the calculation of a projection structure at 12Å resolution by negative-stain TEM. Immunogold labeling TEM revealed the sidedness of the tubular crystals, and high-resolution atomic force microscopy the surface structure of mut-IIC. This work presents the structure of a glucose PTS transporter at the highest resolution achieved so far and sets the basis for future structural studies.

From Structure to Function: Characterization of Cu(I) Adducts in Leveler Additives by DFT Calculations

We present the first molecular model of the coordination complex formed by Cu(I) and imidazole-epichlorohydrin polymers. Our calculations show that the Cu(I) ion has linear coordination and the whole complex has neutral charge. Our model suggests salt couple pairing as the driving force for the formation of the surface-confined precipitation, which is crucial to obtain flat surfaces in industrial copper deposition processes, required for mass fabrication of state-of-the-art electronic and memory devices.

Relationship between glomerular filtration rate and the adipokines adiponectin, resistin and leptin in coronary patients with predominantly normal or mildly impaired renal function

BACKGROUND: Due to their molecular weight, it is possible that the adipokines adiponectin, resistin and leptin accumulate when glomerular filtration rate (GFR) is decreased. In reduced renal clearance, altered serum concentrations of these proteins might affect cardiovascular risk. The objective of the study was to investigate the relationship between adipokine concentrations and GFR. METHODS: The association between GFR, as determined by the abbreviated MDRD equation, and the concentrations of the adipokines adiponectin, resistin and leptin was assessed in a cohort of coronary patients (n=538; 363 male, 165 female). After calculation of correlations between GFR and adipokine concentrations, the association was further assessed by analysis of covariance following adjustment for age, gender, BMI, presence of type 2 diabetes, presence of hypertension, history of smoking as well as for serum lipid concentrations. RESULTS: Mean GFR in our study population was 68.74+/-15.27 ml/min/1.73 m(2). 74.3% of the patients had a GFR >60 ml/min/1.73 m(2), 24% of the patients had a GFR between 30 and 60 ml/min/1.73 m(2), and 1.7% of the patients had a GFR <30 ml/min/1.73 m(2). There were significant inverse correlations between adiponectin (r=-0.372; p<0.001), resistin (r=-0.227; p<0.001) and leptin (r=-0.151; p=0.009) concentrations and GFR. After multivariate adjustment, the associations remained significant for adiponectin and resistin. Subgroup analysis in patients with GFR >60 ml/min/1.73 m(2) showed a significant correlation between GFR and adiponectin as well as leptin concentrations. However, after adjustment, these associations no longer were significant. CONCLUSIONS: There is an independent association between GFR and the serum concentrations of adiponectin and resistin. However, this association is not present at GFR >60 ml/min/1.73 m(2). This finding suggests that adipokine concentrations in mildly impaired and normal renal function are influenced by factors other than GFR.

Relationship between renal resistance index and renal function in liver transplant recipients after cessation of calcineurin inhibitor

End stage renal disease is a major complication after orthotopic liver transplantation (OLT). Vasoconstriction of renal arterial vessels because of calcineurin inhibitor (CNI) treatment plays a pivotal role in the development of renal insufficiency following OLT. Renal resistance can be measured non-invasively by determining the resistance index (RI) of segmental arteries by color-coded duplex ultrasonography, a measure with predictive value for future renal failure. Sixteen OLT patients on long-term CNI therapy were recruited prospectively and randomly assigned either to receive the m-TOR inhibitor sirolimus (SRL) or to continue on CNI treatment, and were followed for one yr. Serum creatinine (crea) declined after conversion to SRL, whereas it tended to increase in patients remaining on CNI (meanDelta crea SRL: -27, -18, -18, -15 micromol/L; meanDelta crea CNI: 4, 5, 8, 11 micromol/L at 1, 3, 6, 12 months, p = 0.02). RI improved after switching to SRL and was lower on SRL than on CNI (meanDeltaRI SRL: -0.04, -0.04, -0.03, -0.03; meanDeltaRI CNI: -0.006, 0.004, -0.007, -0.01 after 1, 3, 6, 12 months, p = 0.016). Individual changes of RI correlated significantly with individual changes of crea (r = 0.54, p < 0.001). Conversion from CNI to SRL can ameliorate renal function accompanied by a reduction of intrarenal RI after OLT.

Relationship between inflammation and cognitive function in obstructive sleep apnea

OBJECTIVES: Obstructive sleep apnea (OSA) can have adverse effects on cognitive functioning, mood, and cardiovascular functioning. OSA brings with it disturbances in sleep architecture, oxygenation, sympathetic nervous system function, and inflammatory processes. It is not clear which of these mechanisms is linked to the decrease in cognitive functioning. This study examined the effect of inflammatory parameters on cognitive dysfunction. MATERIALS AND METHODS: Thirty-nine patients with untreated sleep apnea were evaluated by polysomnography and completed a battery of neuropsychological tests. After the first night of evaluation in the sleep laboratory, blood samples were taken for analysis of interleukin 6, tumor necrosis factor-alpha (TNF-alpha), and soluble TNF receptor 1 (sTNF-R1). RESULTS: sTNF-R1 significantly correlated with cognitive dysfunction. In hierarchical linear regression analysis, measures of obstructive sleep apnea severity explained 5.5% of the variance in cognitive dysfunction (n.s.). After including sTNF-R1, percentage of variance explained by the full model increased more than threefold to 19.6% (F = 2.84, df = 3, 36, p = 0.05). Only sTNF-R1 had a significant individual relationship with cognitive dysfunction (beta = 0.376 t = 2.48, p = 0.02). CONCLUSIONS: sTNF-R1 as a marker of chronic inflammation may be associated with diminished neuropsychological functioning in patients with OSA.