Peptide design by artificial neural networks and computer-based evolutionary search

A technique for systematic peptide variation by a combination of rational and evolutionary approaches is presented. The design scheme consists of five consecutive steps: (i) identification of a “seed peptide” with a desired activity, (ii) generation of variants selected from a physicochemical space around the seed peptide, (iii) synthesis and testing of this biased library, (iv) modeling of a quantitative sequence-activity relationship by an artificial neural network, and (v) de novo design by a computer-based evolutionary search in sequence space using the trained neural network as the fitness function. This strategy was successfully applied to the identification of novel peptides that fully prevent the positive chronotropic effect of anti-β1-adrenoreceptor autoantibodies from the serum of patients with dilated cardiomyopathy. The seed peptide, comprising 10 residues, was derived by epitope mapping from an extracellular loop of human β1-adrenoreceptor. A set of 90 peptides was synthesized and tested to provide training data for neural network development. De novo design revealed peptides with desired activities that do not match the seed peptide sequence. These results demonstrate that computer-based evolutionary searches can generate novel peptides with substantial biological activity.

A peptide export–import control circuit modulating bacterial development regulates protein phosphatases of the phosphorelay

The phosphorelay signal transduction system activates developmental transcription in sporulation of Bacillus subtilis by phosphorylation of aspartyl residues of the Spo0F and Spo0A response regulators. The phosphorylation level of these response regulators is determined by the opposing activities of protein kinases and protein aspartate phosphatases that interpret positive and negative signals for development in a signal integration circuit. The RapA protein aspartate phosphatase of the phosphorelay is regulated by a peptide that directly inhibits its activity. This peptide is proteolytically processed from an inactive pre-inhibitor protein encoded in the phrA gene. The pre-inhibitor is cleaved by the protein export apparatus to a putative pro-inhibitor that is further processed to the active inhibitor peptide and internalized by the oligopeptide permease. This export–import circuit is postulated to be a mechanism for timing phosphatase activity where the processing enzymes regulate the rate of formation of the active inhibitor. The processing events may, in turn, be controlled by a regulatory hierarchy. Chromosome sequencing has revealed several other phosphatase–prepeptide gene pairs in B. subtilis, suggesting that the use of this mechanism may be widespread in signal transduction.

The identification of CD4+ T cell epitopes with dedicated synthetic peptide libraries

For a large number of T cell-mediated immunopathologies, the disease-related antigens are not yet identified. Identification of T cell epitopes is of crucial importance for the development of immune-intervention strategies. We show that CD4+ T cell epitopes can be defined by using a new system for synthesis and screening of synthetic peptide libraries. These libraries are designed to bind to the HLA class II restriction molecule of the CD4+ T cell clone of interest. The screening is based on three selection rounds using partial release of 14-mer peptides from synthesis beads and subsequent sequencing of the remaining peptide attached to the bead. With this approach, two peptides were identified that stimulate the β cell-reactive CD4+ T cell clone 1c10, which was isolated from a newly diagnosed insulin-dependent diabetes mellitus patient. After performing amino acid-substitution studies and protein database searches, a Haemophilus influenzae TonB-derived peptide was identified that stimulates clone 1c10. The relevance of this finding for the pathogenesis of insulin-dependent diabetes mellitus is currently under investigation. We conclude that this system is capable of determining epitopes for (autoreactive) CD4+ T cell clones with previously unknown peptide specificity. This offers the possibility to define (auto)antigens by searching protein databases and/or to induce tolerance by using the peptide sequences identified. In addition the peptides might be used as leads to develop T cell receptor antagonists or anergy-inducing compounds.

Interaction of pigeon cytochrome c-(43–58) peptide analogs with either T cell antigen receptor or I-Ab molecule

We determined that a pigeon cytochrome c-derived peptide, p43–58, possesses two anchor residues, 46 and 54, for binding with the I-Ab molecule that are compatible to the position 1 (P1) and position 9 (P9) of the core region in the major histocompatibility complex (MHC) class II binding peptides, respectively. In the present study to analyze each binding site between P1 and P9 of p43–58 to either I-Ab or T cell antigen receptor (TCR), we investigated T cell responses to a series of peptides (P2K, P3K, P4K, P5K, P6K, P7K, and P8E) that sequentially substituted charged amino acid residues for the residues at P2 to P8 of p43–58. T cells from C57BL/10 (I-Ab) mice immunized with P4K or P6K did not mount appreciable proliferative responses to the immunogens, but those primed with other peptides (P2K, P3K, P5K, P7K, and P8E) showed substantial responses in an immunogen-specific manner. It was demonstrated by binding studies that P1 and P9 functioned as main anchors and P4 and P6 functioned as secondary anchors to I-Ab. Analyses of Vβ usage of T cell lines specific for these analogs suggested that P8 interacts with the complementarity-determining region 1 (CDR1)/CDR2 of the TCR β chain. Furthermore, sequencing of the TCR on T cell hybridomas specific for these analogs indicated that P5 interacts with the CDR3 of the TCR β chain. The present findings are consistent with the three-dimensional structure of the trimolecular complex that has been reported for TCR/peptide/MHC class I molecules.

Isolation and sequencing of cDNAs for splice variants of growth hormone-releasing hormone receptors from human cancers

The proliferation of various tumors is inhibited by the antagonists of growth hormone-releasing hormone (GHRH) in vitro and in vivo, but the receptors mediating the effects of GHRH antagonists have not been identified so far. Using an approach based on PCR, we detected two major splice variants (SVs) of mRNA for human GHRH receptor (GHRH-R) in human cancer cell lines, including LNCaP prostatic, MiaPaCa-2 pancreatic, MDA-MB-468 breast, OV-1063 ovarian, and H-69 small-cell lung carcinomas. In addition, high-affinity, low-capacity binding sites for GHRH antagonists were found on the membranes of cancer cell lines such as MiaPaCa-2 that are negative for the vasoactive intestinal peptide/pituitary adenylate cyclase-activating polypeptide receptor (VPAC-R) or lines such as LNCaP that are positive for VPAC-R. Sequence analysis of cDNAs revealed that the first three exons in SV1 and SV2 are replaced by a fragment of retained intron 3 having a new putative in-frame start codon. The rest of the coding region of SV1 is identical to that of human pituitary GHRH-R, whereas in SV2 exon 7 is spliced out, resulting in a 1-nt upstream frameshift, which leads to a premature stop codon in exon 8. The intronic sequence may encode a distinct 25-aa fragment of the N-terminal extracellular domain, which could serve as a proposed signal peptide. The continuation of the deduced protein sequence coded by exons 4–13 in SV1 is identical to that of pituitary GHRH-R. SV2 may encode a GHRH-R isoform truncated after the second transmembrane domain. Thus SVs of GHRH-Rs have now been identified in human extrapituitary cells. The findings support the view that distinct receptors are expressed on human cancer cells, which may mediate the antiproliferative effect of GHRH antagonists.

Covalent HLA-B27/peptide complex induced by specific recognition of an aziridine mimic of arginine.

The class I major histocompatibility complex (MHC) glycoprotein HLA-B27 binds short peptides containing arginine at peptide position 2 (P2). The HLA-B27/peptide complex is recognized by T cells both as part of the development of the repertoire of T cells in the cellular immune system and during activation of cytotoxic T cells. Based on the three-dimensional structure of HLA-B27, we have synthesized a ligand with an aziridine-containing side chain designed to mimic arginine and to bind covalently in the arginine-specific P2 pocket of HLA-B27. Using tryptic digestion followed by mass spectrometry and amino acid sequencing, the aziridine-containing ligand is shown to alkylate specifically cysteine 67 of HLA-B27. Neither free cysteine in solution nor an exposed cysteine on a class II MHC molecule can be alkylated, showing that specific recognition between the anchor side-chain pocket of an MHC class I protein and the designed ligand (propinquity) is necessary to induce the selective covalent reaction with the MHC class I molecule.

Conformational dependence of electron transfer across de novo designed metalloproteins.

Flash photolysis and pulse radiolysis measurements demonstrate a conformational dependence of electron transfer rates across a 16-mer helical bundle (three-helix metalloprotein) modified with a capping CoIII(bipyridine)3 electron acceptor at the N terminus and a 1-ethyl-1'-ethyl-4,4'- bipyridinium donor at the C terminus. For the CoIII(peptide)3-1-ethyl-1'-ethyl-4,4'-bipyridinium maquettes, the observed transfer is a first order, intramolecular process, independent of peptide concentration or laser pulse energy. In the presence of 6 M urea, the random coil bundle (approximately 0% helicity) has an observed electron transfer rate constant of kobs = 900 +/- 100 s-1. In the presence of 25% trifluoroethanol (TFE), the helicity of the peptide is 80% and the kobs increases to 2000 +/- 200 s-1. Moreover, the increase in the rate constant in TFE is consistent with the observed decrease in donor-acceptor distance in this solvent. Such bifunctional systems provide a class of molecules for testing the effects of conformation on electron transfer in proteins and peptides.

Expression cloning and functional characterization of the kidney cortex high-affinity proton-coupled peptide transporter.

The presence of a proton-coupled electrogenic high-affinity peptide transporter in the apical membrane of tubular cells has been demonstrated by microperfusion studies and by use of brush border membrane vesicles. The transporter mediates tubular uptake of filtered di- and tripeptides and aminocephalosporin antibiotics. We have used expression cloning in Xenopus laevis oocytes for identification and characterization of the renal high-affinity peptide transporter. Injection of poly(A)+ RNA isolated from rabbit kidney cortex into oocytes resulted in expression of a pH-dependent transport activity for the aminocephalosporin antibiotic cefadroxil. After size fractionation of poly(A)+ RNA the transport activity was identified in the 3.0- to 5.0-kb fractions, which were used for construction of a cDNA library. The library was screened for expression of cefadroxil transport after injection of complementary RNA synthesized in vitro from different pools of clones. A single clone (rPepT2) was isolated that stimulated cefadroxil uptake into oocytes approximately 70-fold at a pH of 6.0. Kinetic analysis of cefadroxil uptake expressed by the transporter's complementary RNA showed a single saturable high-affinity transport system shared by dipeptides, tripeptides, and selected amino-beta-lactam antibiotics. Electrophysiological studies established that the transport activity is electrogenic and affected by membrane potential. Sequencing of the cDNA predicts a protein of 729 amino acids with 12 membrane-spanning domains. Although there is a significant amino acid sequence identity (47%) to the recently cloned peptide transporters from rabbit and human small intestine, the renal transporter shows distinct structural and functional differences.

Cloning and sequencing of CATR1.3, a human gene associated with tumorigenic conversion.

The human squamous cell carcinoma cell line SCC83-01-82 (SCC) contains mutations in both the H-ras and p53 genes, but it exhibits a nontumorigenic phenotype in nude mice. This cell line can be converted into a cell line with a tumorigenic phenotype, SCC83-01-82CA (CA), by treatment with the mutagen methyl methanesulfonate (MMS). This indicates that additional genetic events leading to expression of a cooperating tumor susceptibility gene(s) may be required for tumorigenicity. To identify the cooperating gene(s), an expression cDNA library was made from tumorigenic Ca cells. The library DNA was transfected into nontumorigenic SCC cells and the transfected SCC cells were then injected into nude mice for the selection of a tumorigenic phenotype. Tumors developed in 3 of the 18 mice after injection. Several new cell lines were established from these transfected cell-induced tumors and designated as CATR cells. Tumor histology and karyotype analysis of these cells indicated that they were of human epithelial cell origin. All the CATR cells have the library vector sequence integrated in their genome. Cell line CATR1 expressed a single message from the integrated library representing a 1.3-kb cDNA insert that was absent from untransfected SCC cells or MMS-converted CA cells. This 1.3-kb cDNA insert was cloned by PCR amplification of reverse-transcribed CATR1 total RNA and was designated CATR1.3. The nucleotide sequence of CATR1.3 encodes a peptide of 79 amino acids, has a long 3' untranslated region, and represents an unknown gene product that was associated with the tumorigenic conversion due to the transfected expression library.

Ancylostoma caninum anticoagulant peptide: a hookworm-derived inhibitor of human coagulation factor Xa.

Human hookworm infection is a major cause of gastrointestinal blood loss and iron deficiency anemia, affecting up to one billion people in the developing world. These soil-transmitted helminths cause blood loss during attachment to the intestinal mucosa by lacerating capillaries and ingesting extravasated blood. We have isolated the major anticoagulant used by adult worms to facilitate feeding and exacerbate intestinal blood loss. This 8.7-kDa peptide, named the Ancylostoma caninum anticoagulant peptide (AcAP), was purified by using a combination of ion-exchange chromatography, gel-filtration chromatography, and reverse-phase HPLC. N-terminal sequencing of AcAP reveals no homology to any previously identified anticoagulant or protease inhibitor. Single-stage chromogenic assays reveal that AcAP is a highly potent and specific inhibitor of human coagulation, with an intrinsic K*i for the inhibition of free factor Xa of 323.5 pM. In plasma-based clotting time assays, AcAP was more effective at prolonging the prothrombin time than both recombinant hirudin and tick anticoagulant peptide. These data suggest that AcAP, a specific inhibitor of factor Xa, is one of the most potent naturally occurring anticoagulants described to date.

Lethal chondrodysplasia in a family of Holstein cattle is associated with a de novo splice site variant of COL2A1

Background Lethal chondrodysplasia (bulldog syndrome) is a well-known congenital syndrome in cattle and occurs sporadically in many breeds. In 2015, it was noticed that about 12 % of the offspring of the phenotypically normal Danish Holstein sire VH Cadiz Captivo showed chondrodysplasia resembling previously reported bulldog calves. Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor, suggesting the causative allele was not a rare recessive. The normal phenotype of the sire suggested a dominant inheritance with incomplete penetrance or a mosaic mutation. Results Three malformed calves were examined by necropsy, histopathology, radiology, and computed tomography scanning. These calves were morphologically similar and displayed severe disproportionate dwarfism and reduced body weight. The syndrome was characterized by shortening and compression of the body due to reduced length of the spine and the long bones of the limbs. The vicerocranium had severe dysplasia and palatoschisis. The bones had small irregular diaphyses and enlarged epiphyses consisting only of chondroid tissue. The sire and a total of four affected half-sib offspring and their dams were genotyped with the BovineHD SNP array to map the defect in the genome. Significant genetic linkage was obtained for several regions of the bovine genome including chromosome 5 where whole genome sequencing of an affected calf revealed a COL2A1 point mutation (g.32473300 G > A). This private sequence variant was predicted to affect splicing as it altered the conserved splice donor sequence GT at the 5’-end of COL2A1 intron 36, which was changed to AT. All five available cases carried the mutant allele in heterozygous state and all five dams were homozygous wild type. The sire VH Cadiz Captivo was shown to be a gonadal and somatic mosaic as assessed by the presence of the mutant allele at levels of about 5 % in peripheral blood and 15 % in semen. Conclusions The phenotypic and genetic findings are comparable to a previously reported COL2A1 missense mutation underlying lethal chondrodysplasia in the offspring of a mosaic French Holstein sire (Igale Masc). The identified independent spontaneous splice site variant in COL2A1 most likely caused chondrodysplasia and must have occurred during the early foetal development of the sire. This study provides a first example of a dominant COL2A1 splice site variant as candidate causal mutation of a severe lethal chondrodysplasia phenotype. Germline mosaicism is a relatively frequent mechanism in the origin of genetic disorders and explains the prevalence of a certain fraction of affected offspring. Paternal dominant de novo mutations are a risk in cattle breeding, especially because the ratio of defective offspring may be very high and be associated with significant animal welfare problems.

Baltikinin: A New Myotropic Tryptophyllin-3 Peptide Isolated from the Skin Secretion of the Purple-Sided Leaf Frog, Phyllomedusa baltea

Here we report the identification of a novel tryptophyllin-3 peptide with arterial smooth muscle relaxation activity from the skin secretion of the purple-sided leaf frog, Phyllomedusa baltea. This new peptide was named baltikinin and had the following primary structure, pGluDKPFGPPPIYPV, as determined by tandem mass spectrometry (MS/MS) fragmentation sequencing and from cloned skin precursor-encoding cDNA. A synthetic replicate of baltikinin was found to have a similar potency to bradykinin in relaxing arterial smooth muscle (half maximal effective concentration (EC50) is 7.2 nM). These data illustrate how amphibian skin secretions can continue to provide novel potent peptides that act through functional targets in mammalian tissues.

Coarse Graining to Invesitigate Peptide-Lipid Interactions

Experimental characterization of molecular details is challenging, and although single molecule experiments have gained prominence, oligomer characterization remains largely unexplored. The ability to monitor the time evolution of individual molecules while they self assemble is essential in providing mechanistic insights about biological events. Molecular dynamics (MD) simulations can fill the gap in knowledge between single molecule experiments and ensemble studies like NMR, and are increasingly used to gain a better understanding of microscopic properties. Coarse-grained (CG) models aid in both exploring longer length and time scale molecular phenomena, and narrowing down the key interactions responsible for significant system characteristics. Over the past decade, CG techniques have made a significant impact in understanding physicochemical processes. However, the realm of peptide-lipid interfacial interactions, primarily binding, partitioning and folding of amphipathic peptides, remains largely unexplored compared to peptide folding in solution. The main drawback of existing CG models is the inability to capture environmentally sensitive changes in dipolar interactions, which are indigenous to protein folding, and lipid dynamics. We have used the Drude oscillator approach to incorporate structural polarization and dipolar interactions in CG beads to develop a minimalistic peptide model, WEPPROM (Water Explicit Polarizable PROtein Model), and a lipid model WEPMEM (Water Explicit Polarizable MEmbrane Model). The addition of backbone dipolar interactions in a CG model for peptides enabled us to achieve alpha-beta secondary structure content de novo, without any added bias. As a prelude to studying amphipathic peptide-lipid membrane interactions, the balance between hydrophobicity and backbone dipolar interactions in driving ordered peptide aggregation in water and at a hydrophobic-hydrophilic interface, was explored. We found that backbone dipole interactions play a crucial role in driving ordered peptide aggregation, both in water and at hydrophobic-hydrophilic interfaces; while hydrophobicity is more relevant for aggregation in water. A zwitterionic (POPC: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and an anionic lipid (POPS: 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine) are used as model lipids for WEPMEM. The addition of head group dipolar interactions in lipids significantly improved structural, dynamic and dielectric properties of the model bilayer. Using WEPMEM and WEPPROM, we studied membrane-induced peptide folding of a cationic antimicrobial peptide with anticancer activity, SVS-1. We found that membrane-induced peptide folding is driven by both (a) cooperativity in peptide self interaction and (b) cooperativity in membrane-peptide interactions. The dipolar interactions between the peptide and the lipid head-groups contribute to stabilizing folded conformations. The role of monovalent ion size and peptide concentration in driving lipid domain formation in anionic/zwitterionic lipid mixtures was also investigated. Our study suggest monovalent ion size to be a crucial determinant of interaction with lipid head groups, and hence domain formation in lipid mixtures. This study reinforces the role of dipole interactions in protein folding, lipid membrane properties, membrane induced peptide folding and lipid domain formation. Therefore, the models developed in this thesis can be used to explore a multitude of biomolecular processes, both at longer time-scales and larger system sizes.

Transcriptional dynamics of the developing sweet cherry (Prunus avium L.) fruit: sequencing, annotation and expression profiling of exocarp-associated genes

The exocarp, or skin, of fleshy fruit is a specialized tissue that protects the fruit, attracts seed dispersing fruit eaters, and has large economical relevance for fruit quality. Development of the exocarp involves regulated activities of many genes. This research analyzed global gene expression in the exocarp of developing sweet cherry (Prunus avium L., 'Regina'), a fruit crop species with little public genomic resources. A catalog of transcript models (contigs) representing expressed genes was constructed from de novo assembled short complementary DNA (cDNA) sequences generated from developing fruit between flowering and maturity at 14 time points. Expression levels in each sample were estimated for 34 695 contigs from numbers of reads mapping to each contig. Contigs were annotated functionally based on BLAST, gene ontology and InterProScan analyses. Coregulated genes were detected using partitional clustering of expression patterns. The results are discussed with emphasis on genes putatively involved in cuticle deposition, cell wall metabolism and sugar transport. The high temporal resolution of the expression patterns presented here reveals finely tuned developmental specialization of individual members of gene families. Moreover, the de novo assembled sweet cherry fruit transcriptome with 7760 full-length protein coding sequences and over 20 000 other, annotated cDNA sequences together with their developmental expression patterns is expected to accelerate molecular research on this important tree fruit crop.

Efeito sacietogênico de um novo inibidor de tripsina da paçoca do amendoim com aumento plasmático de colecistocinina (cck)

Obesity is increasing, reaching epidemic levels in many regions of the world. Studies have shown that consumption of peanuts influences on weight control and this influence may be due to the action of trypsin inhibitors sacietogênica that condition increased plasma colescistocinina (CCK). Moreover, the peanut has other health benefits, and these assignments are guaranteed to increase their production and consumption of several of its products, including the paçoca peanut. The aim of this study was to identify the presence of a trypsin inhibitor in paçoca peanut and evaluate its effect on food intake, weight gain and histomorphological changes in swiss mice (n = 8) and Wistar rats (n = 6). Experimental diets were prepared based on the AIN-93G and supplemented with tack or peanut trypsin inhibitor partially purified paçoca peanut (AHTI). After each treatment, the animals were anesthetized and euthanized, their bloods were collected by cardiac puncture for the determination of CCK and other biochemical parameters (glucose, triglycerides, total cholesterol, high density lipoprotein, low density lipoprotein, glutamic-pyruvic transaminase, glutamic oxaloacetic transaminase and albumin) and their pancreas removed for histologic and morphometric analysis. The supplementation with paçoca peanut and the AHTI showed a decrease of body weight gain and food intake in both mice and rats, due to the satiety, since the animals showed no evidence of impairment of nutritional status conditioned by consumption the AHTI. There were also observed biochemical or morphological important when compared with controls. However, AHTI led to increased secretion of CCK, a peptide sacietogênico. Thus, these results indicate that AHTI present in paçoca peanut, is able to enhance the secretion of plasma CCK and thereby reduce the weight gain associated with lower food intake of experimenta animals