The Structure and Function of αI Domains in Collagen Receptor and Leukocyte Integrins

Integrins are heterodimeric, signaling transmembrane adhesion receptors that connect the intracellular actin microfilaments to the extracellular matrix composed of collagens and other matrix molecules. Bidirectional signaling is mediated via drastic conformational changes in integrins. These changes also occur in the integrin αI domains, which are responsible for ligand binding by collagen receptor and leukocyte specific integrins. Like intact integrins, soluble αI domains exist in the closed, low affinity form and in the open, high affinity form, and so it is possible to use isolated αI domains to study the factors and mechanisms involved in integrin activation/deactivation. Integrins are found in all mammalian tissues and cells, where they play crucial roles in growth, migration, defense mechanisms and apoptosis. Integrins are involved in many human diseases, such as inflammatory, cardiovascular and metastatic diseases, and so plenty of effort has been invested into developing integrin specific drugs. Humans have 24 different integrins, four of which are collagen receptor (α1β1, α2β1, α10β1, α11β1) and five leukocyte specific integrins (αLβ2, αMβ2, αXβ2, αDβ2, αEβ7). These two integrin groups are quite unselective having both primary and secondary ligands. This work presents the first systematic studies performed on these integrin groups to find out how integrin activation affects ligand binding and selectivity. These kinds of studies are important not only for understanding the partially overlapping functions of integrins, but also for drug development. In general, our results indicated that selectivity in ligand recognition is greatly reduced upon integrin activation. Interestingly, in some cases the ligand binding properties of integrins have been shown to be cell type specific. The reason for this is not known, but our observations suggest that cell types with a higher integrin activation state have lower ligand selectivity, and vice versa. Furthermore, we solved the three-dimensional structure for the activated form of the collagen receptor α1I domain. This structure revealed a novel intermediate conformation not previously seen with any other integrin αI domain. This is the first 3D structure for an activated collagen receptor αI domain without ligand. Based on the differences between the open and closed conformation of the αI domain we set structural criteria for a search for effective collagen receptor drugs. By docking a large number of molecules into the closed conformation of the α2I domain we discovered two polyketides, which best fulfilled the set structural criteria, and by cell adhesion studies we showed them to be specific inhibitors of the collagen receptor integrins.

NMR structure shows that the SARS-unique domain contains a macrodomain fold

The NMR structure of a central segment of the previously annotated "SARS-unique domain" (SUD-M; "middle of the SARS-unique domain") in the SARS coronavirus (SARS-CoV) non-structural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3-residues 528-648, and there is a flexibly extended N-terminal tail with the residues 513-527 and a C-terminal flexible tail of residues 649-651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527-651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly-A and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows 3D structure homology with several helicases and NTP-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

Combining 3D and 2D for less constrained periocular recognition

Periocular recognition has recently become an active topic in biometrics. Typically it uses 2D image data of the periocular region. This paper is the first description of combining 3D shape structure with 2D texture. A simple and effective technique using iterative closest point (ICP) was applied for 3D periocular region matching. It proved its strength for relatively unconstrained eye region capture, and does not require any training. Local binary patterns (LBP) were applied for 2D image based periocular matching. The two modalities were combined at the score-level. This approach was evaluated using the Bosphorus 3D face database, which contains large variations in facial expressions, head poses and occlusions. The rank-1 accuracy achieved from the 3D data (80%) was better than that for 2D (58%), and the best accuracy (83%) was achieved by fusing the two types of data. This suggests that significant improvements to periocular recognition systems could be achieved using the 3D structure information that is now available from small and inexpensive sensors.

Cupin: A candidate molecular structure for the Nep1-like protein family

Background: NEP1-like proteins (NLPs) are a novel family of microbial elicitors of plant necrosis. Some NLPs induce a hypersensitive-like response in dicot plants though the basis for this response remains unclear. In addition, the spatial structure and the role of these highly conserved proteins are not known.Results: We predict a 3d-structure for the beta-rich section of the NLPs based on alignments, prediction tools and molecular dynamics. We calculated a consensus sequence from 42 NLPs proteins, predicted its secondary structure and obtained a high quality alignment of this structure and conserved residues with the two Cupin superfamily motifs. The conserved sequence GHRHDWE and several common residues, especially some conserved histidines, in NLPs match closely the two cupin motifs. Besides other common residues shared by dicot Auxin-Binding Proteins (ABPs) and NLPs, an additional conserved histidine found in all dicot ABPs was also found in all NLPs at the same position.Conclusion: We propose that the necrosis inducing protein class belongs to the Cupin superfamily. Based on the 3d-structure, we are proposing some possible functions for the NLPs.

Investigation of the relationship between the structure and function of Ts2, a neurotoxin from Tityus serrulatus venom

Scorpion toxins targeting voltage-gated sodium (NaV) channels are peptides that comprise 6076 amino acid residues cross-linked by four disulfide bridges. These toxins can be divided in two groups (a and beta toxins), according to their binding properties and mode of action. The scorpion a-toxin Ts2, previously described as a beta-toxin, was purified from the venom of Tityus serrulatus, the most dangerous Brazilian scorpion. In this study, seven mammalian NaV channel isoforms (rNaV1.2, rNaV1.3, rNaV1.4, hNaV1.5, mNaV1.6, rNaV1.7 and rNaV1.8) and one insect NaV channel isoform (DmNaV1) were used to investigate the subtype specificity and selectivity of Ts2. The electrophysiology assays showed that Ts2 inhibits rapid inactivation of NaV1.2, NaV1.3, NaV1.5, NaV1.6 and NaV1.7, but does not affect NaV1.4, NaV1.8 or DmNaV1. Interestingly, Ts2 significantly shifts the voltage dependence of activation of NaV1.3 channels. The 3D structure of this toxin was modeled based on the high sequence identity (72%) shared with Ts1, another T. serrulatus toxin. The overall fold of the Ts2 model consists of three beta-strands and one a-helix, and is arranged in a triangular shape forming a cysteine-stabilized a-helix/beta-sheet (CSa beta) motif.

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.

Spectral-Element and Adjoint 3D Full-Wave Tomography for the Lithosphere of Central Italy

The primary objective of this thesis is to obtain a better understanding of the 3D velocity structure of the lithosphere in central Italy. To this end, I adopted the Spectral-Element Method to perform accurate numerical simulations of the complex wavefields generated by the 2009 Mw 6.3 L’Aquila event and by its foreshocks and aftershocks together with some additional events within our target region. For the mainshock, the source was represented by a finite fault and different models for central Italy, both 1D and 3D, were tested. Surface topography, attenuation and Moho discontinuity were also accounted for. Three-component synthetic waveforms were compared to the corresponding recorded data. The results of these analyses show that 3D models, including all the known structural heterogeneities in the region, are essential to accurately reproduce waveform propagation. They allow to capture features of the seismograms, mainly related to topography or to low wavespeed areas, and, combined with a finite fault model, result into a favorable match between data and synthetics for frequencies up to ~0.5 Hz. We also obtained peak ground velocity maps, that provide valuable information for seismic hazard assessment. The remaining differences between data and synthetics led us to take advantage of SEM combined with an adjoint method to iteratively improve the available 3D structure model for central Italy. A total of 63 events and 52 stations in the region were considered. We performed five iterations of the tomographic inversion, by calculating the misfit function gradient - necessary for the model update - from adjoint sensitivity kernels, constructed using only two simulations for each event. Our last updated model features a reduced traveltime misfit function and improved agreement between data and synthetics, although further iterations, as well as refined source solutions, are necessary to obtain a new reference 3D model for central Italy tomography.

Toward a 3D in vitro model based on decellularized thymus to maintain adult thymic ephitelial cells functionality

During my PhD,I have been develop an innovative technique to reproduce in vitro the 3D thymic microenvironment, to be used for growth and differentiation of thymocytes, and possible transplantation replacement in conditions of depressed thymic immune regulation. The work has been developed in the laboratory of Tissue Engineering at the University Hospital in Basel, Switzerland, under the tutorship of Prof.Ivan Martin. Since a number of studies have suggested that the 3D structure of the thymic microenvironment might play a key role in regulating the survival and functional competence of thymocytes, I’ve focused my effort on the isolation and purification of the extracellular matrix of the mouse thymus. Specifically, based on the assumption that TEC can favour the differentiation of pre-T lymphocytes, I’ve developed a specific decellularization protocol to obtain the intact, DNA-free extracellular matrix of the adult mouse thymus. Two different protocols satisfied the main characteristics of a decellularized matrix, according to qualitative and quantitative assays. In particular, the quantity of DNA was less than 10% in absolute value, no positive staining for cells was found and the 3D structure and composition of the ECM were maintained. In addition, I was able to prove that the decellularized matrixes were not cytotoxic for the cells themselves, and were able to increase expression of MHC II antigens compared to control cells grown in standard conditions. I was able to prove that TECs grow and proliferate up to ten days on top the decellularized matrix. After a complete characterization of the culture system, these innovative natural scaffolds could be used to improve the standard culture conditions of TEC, to study in vitro the action of different factors on their differentiation genes, and to test the ability of TECs to induce in vitro maturation of seeded T lymphocytes.

Three-dimensional consensus structure of sst2-selective somatostatin (SRIF) antagonists by NMR

The three-dimensional NMR structures of seven octapeptide analogs of somatostatin (SRIF), based on octreotide, with the basic sequence H-Cpa/Phe2-c[DCys3-Xxx7-DTrp/DAph(Cbm)8-Lys9-Thr10-Cys14]-Yyy-NH2 (the numbering refers to the position in native SRIF), with Xxx7 being Aph(Cbm)/Tyr/Agl(NMe,benzoyl) and Yyy being Nal/DTyr/Thr, are presented here. Most of these analogs exhibit potent and highly selective binding to sst2 receptors, and all of the analogs are antagonists inhibiting receptor signaling. Based on their consensus 3D structure, the pharmacophore of the sst2-selective antagonist has been defined. The pharmacophore involves the side chains of Cpa2, DTrp/DAph(Cbm)8, and Lys9, with the backbone for most of the sst2-selective antagonists comprised a Type-II' beta-turn. Hence, the sst2-selective antagonist pharmacophore is very similar to the sst2-selective agonist pharmacophore previously described.

Structural conservations and diversities of dsRNA viruses: Structural studies of the cytoplasmic polyhedrosis virus by electron cryomicroscopy and 3D reconstruction

The Reoviridae virus family is a group of economically and pathologically important viruses that have either single-, double-, or triple-shelled protein layers enclosing a segmented double stranded RNA genome. Each virus particle in this family has its own viral RNA dependent RNA polymerase and the enzymatic activities necessary for the mature RNA synthesis. Based on the structure of the inner most cores of the viruses, the Reoviridae viruses can be divided into two major groups. One group of viruses has a smooth surfaced inner core, surrounded by complete outer shells of one or two protein layers. The other group has an inner core decorated with turrets on the five-fold vertices, and could either completely lack or have incomplete outer protein layers. The structural difference is one of the determinant factors for their biological differences during the infection. ^ Cytoplasmic polyhedrosis virus (CPV) is a single-shelled, turreted virus and the structurally simplest member in Reoviridae. It causes specific chronic infections in the insect gut epithelial cells. Due to its wide range of insect hosts, CPV has been engineered as a potential insecticide for use in fruit and vegetable farming. Its unique structural simplicity, unparalleled capsid stability and ease of purification make CPV an ideal model system for studying the structural basis of dsRNA virus assembly at the highest possible resolution by electron cryomicroscopy (cryoEM) and three-dimensional (3D) reconstruction. ^ In this thesis work, I determined the first 3D structure of CPV capsids using 100 kV cryoEM. At an effective resolution of 17 Å, the full capsid reveals a 600-Å diameter, T = 1 icosahedral shell decorated with A and B spikes at the 5-fold vertices. The internal space of the empty CPV is unoccupied except for 12 mushroom-shaped densities that are attributed to the transcriptional enzyme complexes. The inside of the full capsid is packed with icosahedrally-ordered viral genomic RNA. The interactions of viral RNA with the transcriptional enzyme complexes and other capsid proteins suggest a mechanism for RNA transcription and subsequent release. ^ Second, the interactions between the turret proteins (TPs) and the major capsid shell protein (CSPs) have been identified through 3D structural comparisons of the intact CPV capsids with the spikeless CPV capsids, which were generated by chemical treatments. The differential effects of these chemical treatment experiments also indicated that CPV has a significantly stronger structural integrity than other dsRNA viruses, such as the orthoreovirus subcores, which are normally enclosed within outer protein shells. ^ Finally, we have reconstructed the intact CPV to an unprecendented 8 Å resolution from several thousand of 400kV cryoEM images. The 8 Å structure reveals interactions among the 120 molecules of each of the capsid shell protein (CSP), the large protrusion protein (LPP), and 60 molecules of the turret protein (TP). A total of 1980 α-helices and 720 β-sheets have been identified in these capsid proteins. The CSP structure is largely conserved, with the majority of the secondary structures homologous to those observed in the x-ray structures of corresponding proteins of other reoviruses, such as orthoreovirus and bluetongue virus. The three domains of TP are well positioned to play multifunctional roles during viral transcription. The completely non-equivalent interactions between LPP and CSP and those between the anchoring domain of TP and CSP account for the unparalleled stability of this structurally simplest member of the Reoviridae. ^

Fast 2D to 3D conversion using a clustering-based hierarchical search in a machine learning framework

Automatic 2D-to-3D conversion is an important application for filling the gap between the increasing number of 3D displays and the still scant 3D content. However, existing approaches have an excessive computational cost that complicates its practical application. In this paper, a fast automatic 2D-to-3D conversion technique is proposed, which uses a machine learning framework to infer the 3D structure of a query color image from a training database with color and depth images. Assuming that photometrically similar images have analogous 3D structures, a depth map is estimated by searching the most similar color images in the database, and fusing the corresponding depth maps. Large databases are desirable to achieve better results, but the computational cost also increases. A clustering-based hierarchical search using compact SURF descriptors to characterize images is proposed to drastically reduce search times. A significant computational time improvement has been obtained regarding other state-of-the-art approaches, maintaining the quality results.

Large-scale protein structure modeling of the Saccharomyces cerevisiae genome

The function of a protein generally is determined by its three-dimensional (3D) structure. Thus, it would be useful to know the 3D structure of the thousands of protein sequences that are emerging from the many genome projects. To this end, fold assignment, comparative protein structure modeling, and model evaluation were automated completely. As an illustration, the method was applied to the proteins in the Saccharomyces cerevisiae (baker’s yeast) genome. It resulted in all-atom 3D models for substantial segments of 1,071 (17%) of the yeast proteins, only 40 of which have had their 3D structure determined experimentally. Of the 1,071 modeled yeast proteins, 236 were related clearly to a protein of known structure for the first time; 41 of these previously have not been characterized at all.

Herbicide-binding Sites Revealed in the Structure of Plant Acetohydroxyacid Synthase

The sulfonylureas and imidazolinones are potent commercial herbicide families. They are among the most popular choices for farmers worldwide, because they are nontoxic to animals and highly selective. These herbicides inhibit branched-chain amino acid biosynthesis in plants by targeting acetohydroxyacid synthase (AHAS, EC 2.2.1.6). This report describes the 3D structure of Arabidopsis thaliana AHAS in complex with five sulfonylureas (to 2.5 angstrom resolution) and with the imidazolinone, imazaquin (IQ; 2.8 angstrom). Neither class of molecule has a structure that mimics the substrates for the enzyme, but both inhibit by blocking a channel through which access to the active site is gained. The sulfonylureas approach within 5 angstrom of the catalytic center, which is the C2 atom of the cofactor thiamin diphosphate, whereas IQ is at least 7 angstrom from this atom. Ten of the amino acid residues that bind the sulfonylureas also bind IQ. Six additional residues interact only with the sulfonylureas, whereas there are two residues that bind IQ but not the sulfonylureas. Thus, the two classes of inhibitor occupy partially overlapping sites but adopt different modes of binding. The increasing emergence of resistant weeds due to the appearance of mutations that interfere with the inhibition of AHAS is now a worldwide problem. The structures described here provide a rational molecular basis for understanding these mutations, thus allowing more sophisticated AHAS inhibitors to be developed. There is no previously described structure for any plant protein in complex with a commercial herbicide.

Hsc70-induced changes in clathrin-auxilin cage structure suggest a role for clathrin light chains in cage disassembly

The molecular chaperone, Hsc70, together with its cofactor, auxilin, facilitates the ATP-dependent removal of clathrin during clathrin-mediated endocytosis in cells. We have used cryo-electron microscopy to determine the 3D structure of a complex of clathrin, auxilin401-910 and Hsc70 at pH 6 in the presence of ATP, frozen within 20 seconds of adding Hsc70 in order to visualize events that follow the binding of Hsc70 to clathrin and auxilin before clathrin disassembly. In this map,we observe density beneath the vertex of the cage that we attribute to bound Hsc70. This density emerges asymmetrically from the clathrin vertex, suggesting preferential binding by Hsc70 for one of the three possible sites at the vertex. Statistical comparison with a map of whole auxilin and clathrin previously published by us reveals the location of statistically significant differences which implicate involvement of clathrin light chains in structural rearrangements which occur after Hsc70 is recruited. Clathrin disassembly assays using light scattering suggest that loss of clathrin light chains reduces the efficiency with which auxilin facilitates this reaction. These data support a regulatory role for clathrin light chains in clathrin disassembly in addition to their established role in regulating clathrin assembly. © 2013 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.

Confirmação do modelo da estrutura 3D do vale inferior do Tejo a partir de dados de ruído sísmico ambiente

Ao longo de sua história a região do Vale Inferior do Tejo VIT foi abalada por vários sismos consideravelmente destrutivas, tendo alguns deles produzido significativas deformações nas estruturas marítimas localizadas no litoral a sudoeste do território Português; outros, moderados, foram produzidos por fontes locais, como os de 1344, 1531 e 1909. Nos últimos anos, devido à melhoria dos modelos de estrutura 3D e o desenvolvimento dos métodos numéricos, foram elaborados vários estudos de síntese de movimento forte do solo para a região do Baixo Tejo utilizando o método de diferenças finitas. Para confirmar o modelo de velocidades desta bacia usámos medidas de ruído sísmico, aplicámos um método baseado na razão espectral H/V, e, a partir destas curvas, por inversão, obtivemos um modelo de velocidades para a região estudada. Os resultados revelam uma boa concordância entre o modelo obtido e os dados geofísicos e geológicos recolhidos na mesma área._ ABSTRACT: Along his history the Lower Tagus Valley (LTV) area was shaken by several earthquakes. The largest reported had their origin in the southwestern part of Iberia. Other moderate earthquakes were produced by local sources such as the 1344, 1531 and the 1909. ln the last years, due to 3D structural model improvement and development in numerical methods, several studies have successful obtained strong-ground motion synthesis for the LVT region using finite difference method. To confirm the velocity model of the LTV sedimentary basin obtained by geophysical and geological data, we use broad-band microtremor measurements and application of the horizontal to vertical (H/V) spectral ratio method. We have obtained a velocity model for the studied region by inversion of the H/V curve. The results have good agreement geological and geophysical data.