Davide Pirolli

Modeling the ternary complex TCR-Vβ/collagenii(261-273)/HLA-DR4 associated with rheumatoid arthritis

Background It is known that genetic predisposition to rheumatoid arthritis (RA) is associated with the MHC class II allele HLA-DR4 and that residues 261–273 of type II collagen (huCollp261) represent an immunodominant T cell epitope restricted by the DR4 molecule. Despite recent advances in characterization of MHC and T cell receptor (TCR) contacts to this epitope, the atomic details of TCR/huCollp261/HLA-DR4 ternary complex are not known. Methodology/Principal Findings Here we have used computational modeling to get insight into this interaction. A three-dimensional model of the TCR Vβ domain from a DR4+ patient affected by RA has been derived by homology modeling techniques. Subsequently, the structure of the TCR Vβ domain in complex with huCollp261/HLA-DR4 was obtained from a docking approach in conjunction with a filtering procedure based on biochemical information. The best complex from the docking experiments was then refined by 20 ns of molecular dynamics simulation in explicit water. The predicted model is consistent with available experimental data. Our results indicate that residues 97–101 of CDR3β are critical for recognition of huCollp261/HLA-DR4 by TCR. We also show that TCR contacts on p/MHC surface affect the conformation of the shared epitope expressed by DR alleles associated with RA susceptibility. Conclusions/Significance This work presents a three-dimensional model for the ternary complex TCR-Vβ/collagenII(261–273)/HLA-DR4 associated with rheumatoid arthritis that can provide insights into the molecular mechanisms of self reactivity.

Chrono-Proteomics of Human Saliva: Variations of the Salivary Proteome during Human Development

An important contribution to the variability of any proteome is given by the time dimension that should be carefully considered to define physiological modifications. To this purpose, whole saliva proteome was investigated in a wide age range. Whole saliva was collected from 17 preterm newborns with a postconceptional age at birth of 178-217 days. In these subjects sample collection was performed serially starting immediately after birth and within about 1 year follow-up, gathering a total of 111 specimens. Furthermore, whole saliva was collected from 182 subjects aged between 0 and 17 years and from 23 adults aged between 27 and 57 years. The naturally occurring intact salivary proteome of the 316 samples was analyzed by low- and high-resolution HPLC-ESI-MS platforms. Proteins peculiar of the adults appeared in saliva with different time courses during human development. Acidic proline-rich proteins encoded by PRH2 locus and glycosylated basic proline-rich proteins encoded by PRB3 locus appeared following 180 days of postconceptional age, followed at 7 months (±2 weeks) by histatin 1, statherin, and P-B peptide. The other histatins and acidic proline-rich proteins encoded by PRH1 locus appeared in whole saliva of babies from 1 to 3 weeks after the normal term of delivery, S-type cystatins appeared at 1 year (±3 months), and basic proline-rich proteins appeared at 4 years (±1 year) of age. All of the proteinases involved in the maturation of salivary proteins were more active in preterm than in at-term newborns, on the basis of the truncated forms detected. The activity of the Fam20C kinase, involved in the phosphorylation of various proteins, started around 180 days of postconceptional age, slowly increased reaching values comparable to adults at about 2 years (±6 months) of age. Instead, MAPK14 involved in the phosphorylation of S100A9 was fully active since birth also in preterm newborns.

A second Ig-like domain identified in dystroglycan by molecular modelling and dynamics

Dystroglycan (DG) is a cell surface receptor which is composed of two subunits that interact noncovalently, namely α- and β-DG. In skeletal muscle, DG is the central component of the dystrophin-glycoprotein complex (DGC) that anchors the actin cytoskeleton to the extracellular matrix. To date only the three-dimensional structure of the N-terminal region of α-DG has been solved by X-ray crystallography. To expand such a structural analysis, a theoretical molecular model of the murine α-DG C-terminal region was built based on folding recognition/threading techniques. Although there is no a significant (<30%) sequence homology with the N-terminal region of α-DG, protein fold recognition methods found a significant resemblance to the α-DG N-terminal crystallographic structure. Our in silico structural prediction identified two subdomains in this region. Amino acid residues ∼ 500-600 of α-DG were predicted to adopt an immunoglobulin-like (Ig-like) β-sandwich fold. Such modeled domain includes the β-DG binding epitope of α-DG and, confirming our previous experimental results, suggests that the linear epitope (residues 550-565) assumes a β-strand conformation. The remaining segment of the α-DG C-terminal region (residues 601-653) is organized in a coil-helix-coil motif. A 20-ns molecular dynamics simulation in explicit water solvent provided support to the predicted Ig-like model structure. The identification of a second Ig-like domain in DG represents another important step towards a full structural and functional description of the α/β DG interface. Preliminary characterization of a novel recombinant peptide (505-600) encompassing this second Ig-like domain demonstrates that it is soluble and stable, further corroborating our in silico analysis.

Polymorphism, selection and tandem duplication of transferrin genes in Atlantic cod (Gadus morhua) - Conserved synteny between fish monolobal and tetrapod bilobal transferrin loci

BackgroundThe two homologous iron-binding lobes of transferrins are thought to have evolved by gene duplication of an ancestral monolobal form, but any conserved synteny between bilobal and monolobal transferrin loci remains unexplored. The important role played by transferrin in the resistance to invading pathogens makes this polymorphic gene a highly valuable candidate for studying adaptive divergence among local populations.ResultsThe Atlantic cod genome was shown to harbour two tandem duplicated serum transferrin genes (Tf1, Tf2), a melanotransferrin gene (MTf), and a monolobal transferrin gene (Omp). Whereas Tf1 and Tf2 were differentially expressed in liver and brain, the Omp transcript was restricted to the otoliths. Fish, chicken and mammals showed highly conserved syntenic regions in which monolobal and bilobal transferrins reside, but contrasting with tetrapods, the fish transferrin genes are positioned on three different linkage groups. Sequence alignment of cod Tf1 cDNAs from Northeast (NE) and Northwest (NW) Atlantic populations revealed 22 single nucleotide polymorphisms (SNP) causing the replacement of 16 amino acids, including eight surface residues revealed by the modelled 3D-structures, that might influence the binding of pathogens for removal of iron. SNP analysis of a total of 375 individuals from 14 trans-Atlantic populations showed that the Tf1-NE variant was almost fixed in the Baltic cod and predominated in the other NE Atlantic populations, whereas the NW Atlantic populations were more heterozygous and showed high frequencies of the Tf- NW SNP alleles.ConclusionsThe highly conserved synteny between fish and tetrapod transferrin loci infers that the fusion of tandem duplicated Omp-like genes gave rise to the modern transferrins. The multiple nonsynonymous substitutions in cod Tf1 with putative structural effects, together with highly divergent allele frequencies among different cod populations, strongly suggest evidence for positive selection and local adaptation in trans-Atlantic cod populations.

Insights from molecular dynamics simulations: structural basis for the V567D mutation-induced instability of zebrafish alpha-dystroglycan and comparison with the murine model.

A missense amino acid mutation of valine to aspartic acid in 567 position of alpha-dystroglycan (DG), identified in dag1-mutated zebrafish, results in a reduced transcription and a complete absence of the protein. Lacking experimental structural data for zebrafish DG domains, the detailed mechanism for the observed mutation-induced destabilization of the DG complex and membrane damage, remained unclear. With the aim to contribute to a better clarification of the structure-function relationships featuring the DG complex, three-dimensional structural models of wild-type and mutant (V567D) C-terminal domain of alpha-DG from zebrafish were constructed by a template-based modelling approach. We then ran extensive molecular dynamics (MD) simulations to reveal the structural and dynamic properties of the C-terminal domain and to evaluate the effect of the single mutation on alpha-DG stability. A comparative study has been also carried out on our previously generated model of murine alpha-DG C-terminal domain including the I591D mutation, which is topologically equivalent to the V567D mutation found in zebrafish. Trajectories from MD simulations were analyzed in detail, revealing extensive structural disorder involving multiple beta-strands in the mutated variant of the zebrafish protein whereas local effects have been detected in the murine protein. A biochemical analysis of the murine alpha-DG mutant I591D confirmed a pronounced instability of the protein. Taken together, the computational and biochemical analysis suggest that the V567D/I591D mutation, belonging to the G beta-strand, plays a key role in inducing a destabilization of the alpha-DG C-terminal Ig-like domain that could possibly affect and propagate to the entire DG complex. The structural features herein identified may be of crucial help to understand the molecular basis of primary dystroglycanopathies.

Insight into a Novel p53 Single Point Mutation (G389E) by Molecular Dynamics Simulations

The majority of inactivating mutations of p53 reside in the central core DNA binding domain of the protein. In this computational study, we investigated the structural effects of a novel p53 mutation (G389E), identified in a patient with congenital adrenal hyperplasia, which is located within the extreme C-terminal domain (CTD) of p53, an unstructured, flexible region (residues 367–393) of major importance for the regulation of the protein. Based on the three-dimensional structure of a carboxyl-terminal peptide of p53 in complex with the S100B protein, which is involved in regulation of the tumor suppressor activity, a model of wild type (WT) and mutant extreme CTD was developed by molecular modeling and molecular dynamics simulation. It was found that the G389E amino acid replacement has negligible effects on free p53 in solution whereas it significantly affects the interactions of p53 with the S100B protein. The results suggest that the observed mutation may interfere with p53 transcription activation and provide useful information for site-directed mutagenesis experiments.

Characterization of the cell penetrating properties of a human salivary proline-rich peptide.

Saliva contains hundreds of small proline-rich peptides most of which derive from the post-translational and post-secretory processing of the acidic and basic salivary proline-rich proteins. Among these peptides we found that a 20 residue proline-rich peptide (p1932), commonly present in human saliva and patented for its antiviral activity, was internalized within cells of the oral mucosa. The cell-penetrating properties of p1932 have been studied in a primary gingival fibroblast cell line and in a squamous cancer cell line, and compared to its retro-inverso form. We observed by mass-spectrometry, flow cytometry and confocal microscopy that both peptides were internalized in the two cell lines on a time scale of minutes, being the natural form more efficient than the retro-inverso one. The cytosolic localization was dependent on the cell type: both peptide forms were able to localize within nuclei of tumoral cells, but not in the nuclei of gingival fibroblasts. The uptake was shown to be dependent on the culture conditions used: peptide internalization was indeed effective in a complete medium than in a serum-free one allowing the hypothesis that the internalization could be dependent on the cell cycle. Both peptides were internalized likely by a lipid raft-mediated endocytosis mechanism as suggested by the reduced uptake in the presence of methyl-ß-cyclodextrin. These results suggest that the natural peptide may play a role within the cells of the oral mucosa after its secretion and subsequent internalization. Furthermore, lack of cytotoxicity of both peptide forms highlights their possible application as novel drug delivery agents.

Understanding the binding of daunorubicin and doxorubicin to NADPH-dependent cytosolic reductases by computational methods

The anthracycline anticancer agents daunorubicin (DAUN) and doxorubicin (DOX) are reduced by different NADPH-dependent cytosolic reductases into their corresponding alcohol metabolites daunorubicinol (DAUNol) and doxorubicinol (DOXol), which have been implicated in the development of chronic cardiomyopathy. To better understand the individual importance of each enzyme in the reduction and to provide deeper insight into the binding at atomic level we performed molecular docking and dynamics simulations of DAUN and DOX into the active sites of human carbonyl reductase 1 (CBR1) and human aldehyde reductase (AKR1A1). Such simulations evidenced a different behavior between the reductases with respect to DAUN and DOX suggesting major contribution of CBR1 in the reduction. The results are in agreement with available experimental data and for each enzyme and anthracycline pair provided the identification of key residues involved in the interactions. The structural models that we have derived could serve as a useful tool for structure-guided drug design studies.

Evolutionary history and adaptive significance of the polymorphic Pan I in migratory and stationary populations of Atlantic cod (Gadus morhua).

The synaptophysin (SYP) family comprises integral membrane proteins involved in vesicle-trafficking events, but the physiological function of several members has been enigmatic for decades. The presynaptic SYP protein controls neurotransmitter release, while SYP-like 2 (SYPL2) contributes to maintain normal Ca(2+)-signaling in the skeletal muscles. The polymorphic pantophysin (Pan I) of Atlantic cod shows strong genetic divergence between stationary and migratory populations, which seem to be adapted to local environmental conditions. We have investigated the functional involvement of Pan I in the different ecotypes by analyzing the 1) phylogeny, 2) spatio-temporal gene expression, 3) structure-function relationship of the Pan I(A) and I(B) protein variants, and 4) linkage to rhodopsin (rho) recently proposed to be associated with different light sensitivities in Icelandic populations of Atlantic cod. We searched for SYP family genes in phylogenetic key species and identified a single syp-related gene in three invertebrate chordates, while four members, Syp, Sypl1, Sypl2 and synaptoporin (Synpr), were found in tetrapods, Comoran coelacanth and spotted gar. Teleost fish were shown to possess duplicated syp, sypl2 and synpr genes of which the sypl2b paralog is identical to Pan I. The ubiquitously expressed cod Pan I codes for a tetra-spanning membrane protein possessing five amino acid substitutions in the first intravesicular loop, but only minor structural differences were shown between the allelic variants. Despite sizable genomic distance (>2.5 Mb) between Pan I and rho, highly significant linkage disequilibrium was found by genotyping shallow and deep water juvenile settlers predominated by the Pan I(A)-rho(A) and Pan I(B)-rho(B) haplotypes, respectively. However, the predicted rhodopsin protein showed no amino acid changes, while multiple polymorphic sites in the upstream region might affect the gene expression and pigment levels in stationary and migratory cod. Alternatively, other strongly linked genes might be responsible for the sharp settling stratification of juveniles and the different vertical behavior patterns of adult Atlantic cod.

Protein kinase C mediates caspase 3 activation: A role for erythrocyte morphology changes

We have previously showed that morphological alterations in Red Blood Cells (RBCs) are correlated to an impaired eNOS enzymatic activity and a concomitant reduced NO derived metabolites formation. Here we extend our previous observations, reporting that RBC morphology is regulated by a series of specific cell signaling events linked to Protein Kinase C (PKC)-mediated activation of caspase 3. Pretreatment of RBCs with the PKC inhibitor chelerythrine, prior to the addition of phorbol-12-myristate-13-acetate (PMA), an activator of PKC, blocks the appearance of the morphology alterations and the sustained decrease in nitrates and nitrites levels induced by PMA. Inhibition of PKC also completely inhibits PMA mediated caspase-3 activation. On the other hand, caspase 3 inhibition, lessens the PMA induced-effects on the appearance of RBC morphology alterations, although it enhances PMA-mediated effects on nitric oxide (NO) derived metabolites levels. These data demonstrate that PKC-mediated activation of caspase 3 is an integral and essential part of signaling pathway in RBCs, that may be a regulatory factor of RBC mechanical properties, through regulation of NO metabolism.