Two monoclonal antibodies to precisely the same epitope of type ii collagen select non-crossreactive phage clones by phage display: Implications for autoimmunity and molecular mimicry

Two monoclonal antibodies (mAb) CB268 and CII-C1 to type II collagen (CII) react with precisely the same conformational epitope constituted by the residues ARGLT on the three chains of the CII triple helix. The antibodies share structural similarity, with most differences in the complementarity determining region 3 of the heavy chain (HCDR3). The fine reactivity of these mAbs was investigated by screening two nonameric phage-displayed random peptide libraries. For each mAb, there were phage clones (phagotopes) that reacted strongly by ELISA only with the selecting mAb, and inhibited binding to CII only for that mAb, not the alternate mAb. Nonetheless, a synthetic peptide RRLPFGSQM corresponding to an insert from a highly reactive CII-C1-selected phagotope, which was unreactive (and non-inhibitory) with CB268, inhibited the reactivity of CB268 with CII. Most phage-displayed peptides contained a motif in the first part of the molecule that consisted of two basic residues adjacent to at least one hydrophobic residue (e.g. RRL or LRR), but the second portion of the peptides differed for the two mAbs. We predict that conserved CDR sequences interact with the basic-basic-hydrophobic motif, whereas non-conserved amino acids in the binding sites (especially HCDR3) interact with unique peptide sequences and limit cross-reactivity. The observation that two mAbs can react identically with a single epitope on one antigen (CII), but show no cross-reactivity when tested against a second (phagotope) indicates that microorganisms could exhibit mimics capable of initiating autoimmunity without this being evident from conventional assays.

Conformational epitopes on the diabetes autoantigen GAD65 identified by peptide phage display and molecular modeling

The major diabetes autoantigen, glutamic acid decarboxylase (GAD65), contains a region of sequence similarity, including six identical residues PEVKEK, to the P2C protein of coxsackie B virus, suggesting that cross-reactivity between coxsackie B virus and GAD65 can initiate autoimmune diabetes. We used the human islet cell mAbs MICA3 and MICA4 to identify the Ab epitopes of GAD65 by screening phage-displayed random peptide libraries. The identified peptide sequences could be mapped to a homology model of the pyridoxal phosphate (PLP) binding domain of GAD65. For MICA3, a surface loop containing the sequence PEVKEK and two adjacent exposed helixes were identified in the PLP binding domain as well as a region of the C terminus of GAD65 that has previously been identified as critical for MICA3 binding. To confirm that the loop containing tile PEVKEK sequence contributes to the MICA3 epitope, this loop was deleted by mutagenesis. This reduced binding of MICA3 by 70%. Peptide sequences selected using MICA4 were rich in basic or hydroxyl-containing amino acids, and the surface of the GAD65 PLP-binding domain surrounding Lys358, which is known to be critical for MICA4 binding, was likewise rich in these amino acids. Also, the two phage most reactive width MICA4 encoded the motif VALxG, and the reverse of this sequence, LAV, was located in this same region. Thus, we have defined the MICA3 and MICA4 epitopes on GAD65 using the combination of phage display, molecular modeling, and mutagenesis and have provided compelling evidence for the involvement of the PEVKEK loop in the MICA3 epitope.

Crystal structures of the endoplasmic reticulum aminopeptidase-1 (ERAP1) reveal the molecular basis for N-terminal peptide trimming

Endoplasmatic reticulum aminopeptidase 1 (ERAP1) is a multifunctional enzyme involved in trimming of peptides to an optimal length for presentation by major histocompatibility complex (MHC) class I molecules. Polymorphisms in ERAP1 have been associated with chronic inflammatory diseases, including ankylosing spondylitis (AS) and psoriasis, and subsequent in vitro enzyme studies suggest distinct catalytic properties of ERAP1 variants. To understand structure-activity relationships of this enzyme we determined crystal structures in open and closed states of human ERAP1, which provide the first snapshots along a catalytic path. ERAP1 is a zinc-metallopeptidase with typical H-E-X-X-H-(X)18-E zinc binding and G-A-M-E-N motifs characteristic for members of the gluzincin protease family. The structures reveal extensive domain movements, including an active site closure as well as three different open conformations, thus providing insights into the catalytic cycle. A K 528R mutant strongly associated with AS in GWAS studies shows significantly altered peptide processing characteristics, which are possibly related to impaired interdomain interactions.

Energy dependence of methyl-radical adsorption on diamond (001)-(2×1) surface

The deposition of hyperthermal CH3 on diamond (001)-(2×1) surface at room temperature has been studied by means of molecular dynamics simulation using the many-body hydrocarbon potential. The energy threshold effect has been observed. That is, with fixed collision geometry, chemisorption can occur only when the incident energy of CH3 is above a critical value (Eth). Increasing the incident energy, dissociation of hydrogen atoms from the incident molecule was observed. The chemisorption probability of CH3 as a function of its incident energy was calculated and compared with that of C2H2. We found that below 10 eV, the chemisorption probability of C2H2 is much lower than that of CH3 on the same surface. The interesting thing is that it is even lower than that of CH3 on a hydrogen covered surface at the same impact energy. It indicates that the reactive CH3 molecule is the more important species than C2H2 in diamond synthesis at low energy, which is in good agreement with the experimental observation.

Nucleotide sequences of an Australian and a Canadian isolate of potato leafroll luteovirus and their relationships with two European isolates

The genomes of an Australian and a Canadian isolate of potato leafroll virus have been cloned and sequenced. The sequences of both isolates are similar (about 93%), but the Canadian isolate (PLRV-C) is more closely related (about 98% identity) to a Scottish (PLRV-S) and a Dutch isolate (PLRV-N) than to the Australian isolate (PLRV-A). The 5'-terminal 18 nucleotide residues of PLRV-C, PLRV-A, PLRV-N and beet western yellows virus have 17 residues in common. In contrast, PLRV-S shows no obvious similarity in this region. PLRV-A and PLRV-C genomic sequences have localized regions of marked diversity, in particular a 600 nucleotide residue sequence in the polymerase gene. These data provide a world-wide perspective on the molecular biology of PLRV strains and their comparison with other luteoviruses and related RNA plant viruses suggests that there are two major subgroups in the plant luteoviruses.

Predicting structural disruption of proteins caused by crossover

We present a machine learning model that predicts a structural disruption score from a protein s primary structure. SCHEMA was introduced by Frances Arnold and colleagues as a method for determining putative recombination sites of a protein on the basis of the full (PDB) description of its structure. The present method provides an alternative to SCHEMA that is able to determine the same score from sequence data only. Circumventing the need for resolving the full structure enables the exploration of yet unresolved and even hypothetical sequences for protein design efforts. Deriving the SCHEMA score from a primary structure is achieved using a two step approach: first predicting a secondary structure from the sequence and then predicting the SCHEMA score from the predicted secondary structure. The correlation coefficient for the prediction is 0.88 and indicates the feasibility of replacing SCHEMA with little loss of precision.

Monoclonal antibody screening of a phage-displayed random peptide library reveals mimotopes of chemokine receptor CCR5: Implications for the tertiary structure of the receptor and for an n-terminal binding site for HIV-1 gp120

The chemokine receptor CCR5 contains seven transmembrane-spanning domains. It binds chemokines and acts as co-receptor for macrophage (m)-tropic (or R5) strains of HIV-1. Monoclonal antibodies (mAb) to CCR5, 3A9 and 5C7, were used for biopanning a nonapeptide cysteine (C)-constrained phage-displayed random peptide library to ascertain contact residues and define tertiary structures of possible epitopes on CCR5. Reactivity of antibodies with phagotopes was established by enzyme-linked immunosorbent assay (ELISA). mAb 3A9 identified a phagotope C-HASIYDFGS-C (3A9/1), and 5C7 most frequently identified C-PHWLRDLRV-C (5C7/1). Corresponding peptides were synthesized. Phagotopes and synthetic peptides reacted in ELISA with corresponding antibodies and synthetic peptides inhibited antibody binding to the phagotopes. Reactivity by immunofluorescence of 3A9 with CCR5 was strongly inhibited by the corresponding peptide. Both mAb 3A9 and 5C7 reacted similarly with phagotopes and the corresponding peptide selected by the alternative mAb. The sequences of peptide inserts of phagotopes could be aligned as mimotopes of the sequence of CCR5. For phage 3A9/1, the motif SIYD aligned to residues at the N terminus and FG to residues on the first extracellular loop; for 5C7/1, residues at the N terminus, first extracellular loop, and possibly the third extracellular loop could be aligned and so would contribute to the mimotope. The synthetic peptides corresponding to the isolated phagotopes showed a CD4-dependent reactivity with gp120 of a primary, m-tropic HIV-1 isolate. Thus reactivity of antibodies raised to CCR5 against phage-displayed peptides defined mimotopes that reflect binding sites for these antibodies and reveal a part of the gp120 binding sites on CCR5.

Phagotopes derived by antibody screening of phage-displayed random peptide libraries vary in immunoreactivity: Studies using an exemplary monoclonal antibody, CII-C1, to type II collagen

Antibody screening of phage-displayed random peptide libraries to identify mimotopes of conformational epitopes is promising. However, because interpretations can be difficult, an exemplary system has been used in the present study to investigate whether variation in the peptide sequences of selected phagotopes corresponded with variation in immunoreactivity. The phagotopes, derived using a well-characterized monoclonal antibody, CII-C1, to a known conformational epitope on type II collagen, C1, were tested by direct and inhibition ELISA for reactivity with CII-C1. A multiple sequence alignment algorithm, PILEUP, was used to sort the peptides expressed by the phagotopes into clusters. A model was prepared of the C1 epitope on type II collagen. The 12 selected phagotopes reacted with CII-C1 by both direct ELISA (titres from < 100-11 200) and inhibition ELISA (20-100% inhibition); the reactivity varied according to the peptide sequence and assay format. The differences in reactivity between the phagotopes were mostly in accord with the alignment, by PILEUP, of the peptide sequences. The finding that the phagotopes functionally mimicked the C1 epitope on collagen was validated in that amino acids RRL at the amino terminal of many of the peptides were topographically demonstrable on the model of the C1 epitope. Notably, one phagotope that expressed the widely divergent peptide C-IAPKRHNSA-C also mimicked the C1 epitope, as judged by reactivity in each of the assays used: these included cross-inhibition of CII-C1 reactivity with each of the other phagotopes and inhibition by a synthetic peptide corresponding to that expressed by the most frequently selected phagotope, RRLPFGSQM. Thus, it has been demonstrated that multiple phage-displayed peptides can mimic the same epitope and that observed immunoreactivity of selected phagotopes with the selecting mAb can depend on the primary sequence of the expressed peptide and also on the assay format used.

Investigating the genetic association between ERAP1 and ankylosing spondylitis

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 × 103). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 × 10-9). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case-control study, the most significant of which was with rs27434 (P = 4.7 × 10-7). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 × 10-9) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression. © 2009 The Author(s).

Genetic, functional, and histopathological evaluation of two C-terminal BRCA1 missense variants

Background: The vast majority of BRCA1 missense sequence variants remain uncharacterised for their possible effect on protein expression and function, and therefore are unclassified in terms of their pathogenicity. BRCA1 plays diverse cellular roles and it is unlikely that any single functional assay will accurately reflect the total cellular implications of missense mutations in this gene. Objective: To elucidate the effect of two BRCA1 variants, 5236G>C (G1706A) and 5242C>A (A1708E) on BRCA1 function, and to survey the relative usefulness of several assays to direct the characterisation of other unclassified variants in BRCA genes. Methods and Results: Data from a range of bioinformatic, genetic, and histopathological analyses, and in vitro functional assays indicated that the 1708E variant was associated with the disruption of different cellular functions of BRCA1. In transient transfection experiments in T47D and 293T cells, the 1708E product was mislocalised to the cytoplasm and induced centrosome amplification in 293T cells. The 1708E variant also failed to transactivate transcription of reporter constructs in mammalian transcriptional transactivation assays. In contrast, the 1706A variant displayed a phenotype comparable to wildtype BRCA1 in these assays. Consistent with functional data, tumours from 1708E carriers showed typical BRCA1 pathology, while tumour material from 1706A carriers displayed few histopathological features associated with BRCA1 related tumours. Conclusions: A comprehensive range of genetic, bioinformatic, and functional analyses have been combined for the characterisation of BRCA1 unclassified sequence variants. Consistent with the functional analyses, the combined odds of causality calculated for the 1706A variant after multifactorial likelihood analysis (1:142) indicates a definitive classification of this variant as "benign". In contrast, functional assays of the 1708E variant indicate that it is pathogenic, possibly through subcellular mislocalisation. However, the combined odds of 262:1 in favour of causality of this variant does not meet the minimal ratio of 1000:1 for classification as pathogenic, and A1708E remains formally designated as unclassified. Our findings highlight the importance of comprehensive genetic information, together with detailed functional analysis for the definitive categorisation of unclassified sequence variants. This combination of analyses may have direct application to the characterisation of other unclassified variants in BRCA1 and BRCA2.

Models of protein linear molecular motors for dynamic nanodevices

Protein molecular motors are natural nano-machines that convert the chemical energy from the hydrolysis of adenosine triphosphate into mechanical work. These efficient machines are central to many biological processes, including cellular motion, muscle contraction and cell division. The remarkable energetic efficiency of the protein molecular motors coupled with their nano-scale has prompted an increasing number of studies focusing on their integration in hybrid micro- and nanodevices, in particular using linear molecular motors. The translation of these tentative devices into technologically and economically feasible ones requires an engineering, design-orientated approach based on a structured formalism, preferably mathematical. This contribution reviews the present state of the art in the modelling of protein linear molecular motors, as relevant to the future design-orientated development of hybrid dynamic nanodevices. © 2009 The Royal Society of Chemistry.

The molecular mechanisms utilised by mesenchymal stem cells in migration to acute myocardial infarction

Heart damage caused by acute myocardial infarction (AMI) is a leading cause of death and disability in Australia. Novel therapies are still required for the treatment of this condition due to the poor reparative ability of the heart. As such, cellular therapies that assist in the recovery of heart muscle are of great current interest. Culture expanded mesenchymal stem cells (MSC) represent a stem and progenitor cell population that has been shown to promote tissue recovery in pre-clinical studies of AMI. For MSC-based therapies in the clinic, an intravenous route of administration would ideally be used due to the low cost, ease of delivery and relative safety. The study of MSC migration is therefore clinically relevant for a minimally invasive cell therapy to promote regeneration of damaged tissue. C57BL/6, UBI-GFP-BL/6 and CD44-/-/GFP+/+ mice were utilised to investigate mMSC migration. To assist in murine models of MSC migration, a novel method was used for the isolation of murine MSC (mMSC). These mMSC were then expanded in culture and putative mMSC were positive for Sca-1, CD90.2, and CD44 and were negative for CD45 and CD11b. Furthermore, mMSC from C57BL/6 and UBI-GFP-BL/6 mice were shown to differentiate into cells of the mesodermal lineage. Cells from CD44-/-/GFP+/+ mice were positive for Sca-1 and CD90.2, and negative for CD44, CD45 and CD11b however, these cells were unable to differentiate into adipocytes and chondrocytes and express lineage specific genes, PLIN and ACAN. Analysis of mMSC chemokine receptor (CR) expression showed that although mMSC do express chemokine receptors, (including those specific for chemokines released after AMI), these were low or undetectable by mRNA. However, protein expression could be detected, which was predominantly cytoplasmic. It was further shown that in both healthy (unperturbed) and inflamed tissues, mMSC had very little specific migration and engraftment after intravenous injection. To determine if poor mMSC migration was due to the inability of mMSC to respond to chemotactic stimuli, chemokine expression in bone marrow, skin injury and hearts (healthy and after AMI) was analysed at various time points by quantitative real-time PCR (qRT PCR). Many chemokines were up-regulated after skin biopsy and AMI, but the highest acute levels were found for CXCL12 and CCL7. Due to their high expression in infarcted hearts, the chemokines CXCL12 and CCL7 were tested for their effect on mMSC migration. Despite CR expression at both protein and mRNA levels, migration in response to CXCL12 and CCL7 was low in mMSC cultured on Nunclon plastic. A novel tissue culture plastic technology (UpCellTM) was then used that allowed gentle non-enzymatic dissociation of mMSC, thus preserving surface expression of the CRs. Despite this the in vitro data indicated that CXCL12 fails to induce significant migration ability of mMSC, while CCL7 induces significant, but low-level migration. We speculated this may be because of low levels of surface expression of chemokine receptors. In a strategy to increase cell surface expression of mMSC chemokine receptors and enhance their in vitro and in vivo migration capacity, mMSC were pre-treated with pro-inflammatory cytokines. Increased levels of both mRNA and surface protein expression were found for CRs by pre-treating mMSC with pro-inflammatory cytokines including TNF-á, IFN-ã, IL-1á and IL-6. Furthermore, the chemotactic response of mMSC to CXCL12 and CCL7 was significantly higher with these pretreated cells. Finally, the effectiveness of this type of cell manipulation was demonstrated in vivo, where mMSC pre-treated with TNF-á and IFN-ã showed significantly increased migration in skin injury and AMI models. Therefore this thesis has demonstrated, using in vitro and in vivo models, the potential for prior manipulation of MSC as a possible means for increasing the utility of intravenously delivery for MSC-based cellular therapies.

Probability distributed time delays: integrating spatial effects into temporal models

Background In order to provide insights into the complex biochemical processes inside a cell, modelling approaches must find a balance between achieving an adequate representation of the physical phenomena and keeping the associated computational cost within reasonable limits. This issue is particularly stressed when spatial inhomogeneities have a significant effect on system's behaviour. In such cases, a spatially-resolved stochastic method can better portray the biological reality, but the corresponding computer simulations can in turn be prohibitively expensive. Results We present a method that incorporates spatial information by means of tailored, probability distributed time-delays. These distributions can be directly obtained by single in silico or a suitable set of in vitro experiments and are subsequently fed into a delay stochastic simulation algorithm (DSSA), achieving a good compromise between computational costs and a much more accurate representation of spatial processes such as molecular diffusion and translocation between cell compartments. Additionally, we present a novel alternative approach based on delay differential equations (DDE) that can be used in scenarios of high molecular concentrations and low noise propagation. Conclusions Our proposed methodologies accurately capture and incorporate certain spatial processes into temporal stochastic and deterministic simulations, increasing their accuracy at low computational costs. This is of particular importance given that time spans of cellular processes are generally larger (possibly by several orders of magnitude) than those achievable by current spatially-resolved stochastic simulators. Hence, our methodology allows users to explore cellular scenarios under the effects of diffusion and stochasticity in time spans that were, until now, simply unfeasible. Our methodologies are supported by theoretical considerations on the different modelling regimes, i.e. spatial vs. delay-temporal, as indicated by the corresponding Master Equations and presented elsewhere.