JUCHMME: a Java Utility for Class Hidden Markov Models and Extensions for biological sequence analysis.

SUMMARY: JUCHMME is an open-source software package designed to fit arbitrary custom Hidden Markov Models (HMMs) with a discrete alphabet of symbols. We incorporate a large collection of standard algorithms for HMMs as well as a number of extensions and evaluate the software on various biological problems. Importantly, the JUCHMME toolkit includes several additional features that allow for easy building and evaluation of custom HMMs, which could be a useful resource for the research community. AVAILABILITY AND IMPLEMENTATION: http://www.compgen.org/tools/juchmme, https://github.com/pbagos/juchmme. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization of the c.793-1G > A splicing variant in CHEK2 gene as pathogenic: a case report.

BACKGROUND: CHEK2 is involved in the DNA damage repair response Fanconi anemia (FA)-BRCA pathway. An increased risk for breast and other cancers has been documented in individuals who carry a single pathogenic CHEK2 variant. As for other genes involved in cancer predisposition, different types of pathogenic variants have been observed, including single nucleotide variations, short insertions/deletions, large genomic rearrangements and splicing variants. Splicing variants occurring in the splicing acceptor or donor site result in alternative mature mRNA produced and can cause intron retention, exon skipping, or creation of alternative 3' and 5' splice site. Thus, the pathogenicity of this type of alterations should always be explored experimentally and their effect in the mRNA and consequently the protein produced, should be defined. The aim of this study was the delineation of the effect of a splicing variant in the CHEK2 gene. CASE PRESENTATION: A healthy 28-year-old woman with a family history of breast and ovarian cancer was referred for genetic testing. The variant c.793-1G > A (rs730881687) was identified by Next Generation Sequencing (NGS) using a solution-based capture method, targeting 33 cancer predisposition genes (SeqCap EZ Probe library, Roche NimbleGen). Experimental analysis in patient-derived leukocytes using RT-PCR of mRNA followed by cDNA sequencing revealed the deletion of one base from the alternative transcript created (r.793del). This resulted in a frameshift leading to premature termination codon within exon 7 (p.(Asp265Thrfs*10)). CONCLUSIONS: This finding suggests that the CHEK2 splicing variant c.793-1G > A is a deleterious variant. Our case shows that RNA analysis is a valuable tool for uncharacterized splice site variants in individuals referred for testing and facilitates their personalized management.

Analysis of hereditary cancer syndromes by using a panel of genes: novel and multiple pathogenic mutations.

BACKGROUND: Hereditary cancer predisposition syndromes are responsible for approximately 5-10% of all diagnosed cancer cases. In the past, single-gene analysis of specific high risk genes was used for the determination of the genetic cause of cancer heritability in certain families. The application of Next Generation Sequencing (NGS) technology has facilitated multigene panel analysis and is widely used in clinical practice, for the identification of individuals with cancer predisposing gene variants. The purpose of this study was to investigate the extent and nature of variants in genes implicated in hereditary cancer predisposition in individuals referred for testing in our laboratory. METHODS: In total, 1197 individuals from Greece, Romania and Turkey were referred to our laboratory for genetic testing in the past 4 years. The majority of referrals included individuals with personal of family history of breast and/or ovarian cancer. The analysis of genes involved in hereditary cancer predisposition was performed using a NGS approach. Genomic DNA was enriched for targeted regions of 36 genes and sequencing was carried out using the Illumina NGS technology. The presence of large genomic rearrangements (LGRs) was investigated by computational analysis and Multiplex Ligation-dependent Probe Amplification (MLPA). RESULTS: A pathogenic variant was identified in 264 of 1197 individuals (22.1%) analyzed while a variant of uncertain significance (VUS) was identified in 34.8% of cases. Clinically significant variants were identified in 29 of the 36 genes analyzed. Concerning the mutation distribution among individuals with positive findings, 43.6% were located in the BRCA1/2 genes whereas 21.6, 19.9, and 15.0% in other high, moderate and low risk genes respectively. Notably, 25 of the 264 positive individuals (9.5%) carried clinically significant variants in two different genes and 6.1% had a LGR. CONCLUSIONS: In our cohort, analysis of all the genes in the panel allowed the identification of 4.3 and 8.1% additional pathogenic variants in other high or moderate/low risk genes, respectively, enabling personalized management decisions for these individuals and supporting the clinical significance of multigene panel analysis in hereditary cancer predisposition.

MBPpred: Proteome-wide detection of membrane lipid-binding proteins using profile Hidden Markov Models.

A large number of modular domains that exhibit specific lipid binding properties are present in many membrane proteins involved in trafficking and signal transduction. These domains are present in either eukaryotic peripheral membrane or transmembrane proteins and are responsible for the non-covalent interactions of these proteins with membrane lipids. Here we report a profile Hidden Markov Model based method capable of detecting Membrane Binding Proteins (MBPs) from information encoded in their amino acid sequence, called MBPpred. The method identifies MBPs that contain one or more of the Membrane Binding Domains (MBDs) that have been described to date, and further classifies these proteins based on their position in respect to the membrane, either as peripheral or transmembrane. MBPpred is available online at http://bioinformatics.biol.uoa.gr/MBPpred. This method was applied in selected eukaryotic proteomes, in order to examine the characteristics they exhibit in various eukaryotic kingdoms and phyla.

HMMpTM: improving transmembrane protein topology prediction using phosphorylation and glycosylation site prediction.

During the last two decades a large number of computational methods have been developed for predicting transmembrane protein topology. Current predictors rely on topogenic signals in the protein sequence, such as the distribution of positively charged residues in extra-membrane loops and the existence of N-terminal signals. However, phosphorylation and glycosylation are post-translational modifications (PTMs) that occur in a compartment-specific manner and therefore the presence of a phosphorylation or glycosylation site in a transmembrane protein provides topological information. We examine the combination of phosphorylation and glycosylation site prediction with transmembrane protein topology prediction. We report the development of a Hidden Markov Model based method, capable of predicting the topology of transmembrane proteins and the existence of kinase specific phosphorylation and N/O-linked glycosylation sites along the protein sequence. Our method integrates a novel feature in transmembrane protein topology prediction, which results in improved performance for topology prediction and reliable prediction of phosphorylation and glycosylation sites. The method is freely available at http://bioinformatics.biol.uoa.gr/HMMpTM.

Analysis of Molecular Recognition Features (MoRFs) in membrane proteins.

Molecular Recognition Features (MoRFs) are defined as short, intrinsically disordered regions in proteins that undergo disorder-to-order transition upon binding to their partners. As their name suggests, they are implicated in molecular recognition, which serves as the initial step for protein-protein interactions. Membrane proteins constitute approximately 30% of fully sequenced proteomes and are responsible for a wide variety of cellular functions. The aim of the current study was to identify and analyze MoRFs in membrane proteins. Two datasets of MoRFs, transmembrane and peripheral membrane protein MoRFs, were constructed from the Protein Data Bank, and sequence, structural and functional analysis was performed. Characterization of our datasets revealed their unique compositional biases and membrane protein MoRFs were categorized depending on their secondary structure after the interaction with their partners. Moreover, the position of transmembrane protein MoRFs in relation with the protein's topology was determined. Further studies were focused on functional analyses of MoRF-containing proteins and MoRFs' partners, associating them with protein binding, regulation and cell signaling, indicating half of them as putative hubs in protein-protein interaction networks. In conclusion, we provide insights into the disorder-based protein-protein interactions involving membrane proteins.

mpMoRFsDB: a database of molecular recognition features in membrane proteins.

SUMMARY: Molecular recognition features (MoRFs) are small, intrinsically disordered regions in proteins that undergo a disorder-to-order transition on binding to their partners. MoRFs are involved in protein-protein interactions and may function as the initial step in molecular recognition. The aim of this work was to collect, organize and store all membrane proteins that contain MoRFs. Membrane proteins constitute ∼30% of fully sequenced proteomes and are responsible for a wide variety of cellular functions. MoRFs were classified according to their secondary structure, after interacting with their partners. We identified MoRFs in transmembrane and peripheral membrane proteins. The position of transmembrane protein MoRFs was determined in relation to a protein's topology. All information was stored in a publicly available mySQL database with a user-friendly web interface. A Jmol applet is integrated for visualization of the structures. mpMoRFsDB provides valuable information related to disorder-based protein-protein interactions in membrane proteins. AVAILABILITY: http://bioinformatics.biol.uoa.gr/mpMoRFsDB

Germline Co-deletion of CDKN2A and CDKN2B Genes in Pleomorphic Xanthoastrocytoma: Case Report.

BACKGROUND/AIM: Gliomas are highly heterogeneous malignancies originating from diverse cell types within the brain. Although their precise etiology is frequently unknown, risk factors, such as chemical exposure, radiation, and specific uncommon genetic disorders have been identified. Diagnosis typically entails imaging tests, such as magnetic resonance imaging and computed tomography, complemented by a biopsy for confirmation, which may be further validated through genetic testing. CASE REPORT: Next-generation sequencing technology revealed germline co-deletion deletion of cyclin-dependent kinase inhibitor 2 A and B genes (CDKN2A and CDKN2B) in a patient diagnosed with pleomorphic xanthoastrocytoma based on the tumor's molecular characteristics. Following this result, we performed focused genetic analysis with use of multiplex ligation-dependent probe amplification technology for the mother that revealed the same co-deletion. Moreover, due to the father's neuroendocrine pancreatic cancer, application of the NGS technology detected a pathogenic variant in the BRCA1-interacting helicase 1 (BRIP1) gene. Comprehensive multi-gene testing conducted within the familial context, marked by a varied spectrum of cancer type, revealed a constellation of genetic predispositions. CONCLUSION: This case study underscores the critical importance of molecular testing for tumor characterization and highlights the pivotal role of genetic testing in facilitating early intervention and screening for at-risk family members. Furthermore, the identification of germline co-deletions in cancer lays the foundation for the development of targeted therapeutic strategies aimed at restoring normal cellular regulation and improving patient management.

Reclassification of Splicing Gene Variants in Hereditary Cancer: Cases Report and Literature Review.

Alternative splicing (AS), a crucial cellular process, is a source of transcriptomic expansion and protein variability. Its contribution to cancer development and progression among a vast repertoire of human diseases, is highlighted lately and is under extensive investigation. In this review, the relative recent aspects of AS as a hallmark of cancer are described. In parallel, the importance of the identification of splicing-related variants through next-generation sequencing technologies is discussed. Cancer therapy and the management of patients and their families can highly benefit by the classification of these variants.

A novel deletion of exon 4 in the Ectodysplasin A gene associated with X-linked hypohidrotic ectodermal dysplasia.

OBJECTIVE: Identify the disease-causing mutation in a patient with features of X-linked hypohidrotic ectodermal dysplasia, which is a genetic disorder characterized by hypodontia, hypohidrosis and hypotrichosis. It is caused by mutations in Ectodysplasin A gene, which encodes ectodysplasin A, a member of the tumor necrosis factor superfamily. DESIGN: Genetic analysis, was performed using chromosomal microarray analysis, whole exome sequencing and multiplex ligation-dependent probe amplification analysis in a 4-year-old boy with hypohidrotic ectodermal dysplasia features. Moreover, the boy's parents were tested for clinically significant findings identified in order to elucidate the pattern of inheritance of the finding detected in the proband. RESULTS: A novel deletion of entire exon 4 in Ectodysplasin A gene identified in the 4-year-old patient. This deletion was found in heterozygous state in the mother of the proband and was not detected in his father. RNA analysis revealed an in-frame deletion r.527_706del, p.(176_236del) in exon 4 of the Ectodysplasin A gene. CONCLUSION: We identified a novel gross deletion in the Ectodysplasin A gene in a male patient with X-linked hypohidrotic ectodermal dysplasia. Clinical and molecular genetic analysis are crucial to set an accurate diagnosis in patients with hypohidrotic ectodermal dysplasia. These results highlight the importance of the collagen domain of Ectodysplasin A, encoded by exon 4, for its function in vivo.

Digging into the NGS Information from a Large-Scale South European Population with Metastatic/Unresectable Pancreatic Ductal Adenocarcinoma: A Real-World Genomic Depiction.

Despite ongoing oncological advances, pancreatic ductal adenocarcinoma (PDAC) continues to have an extremely poor prognosis with limited targeted and immunotherapeutic options. Its genomic background has not been fully characterized yet in large-scale populations all over the world. Methods: Replicating a recent study from China, we collected tissue samples from consecutive Greek patients with pathologically-confirmed metastatic/unresectable PDAC and retrospectively investigated their genomic landscape using next generation sequencing (NGS). Findings: From a cohort of 409 patients, NGS analysis was successfully achieved in 400 cases (56.50% males, median age: 61.8 years). Consistent with a previous study, KRAS was the most frequently mutated gene in 81.50% of tested samples, followed by TP53 (50.75%), CDKN2 (8%), and SMAD4 (7.50%). BRCA1/2 variants with on-label indications were detected in 2%, and 87.50% carried a variant associated with off-label treatment (KRAS, ERBB2, STK11, or HRR-genes), while 3.5% of the alterations had unknown/preliminary-studied actionability (TP53/CDKN2A). Most of HRR-alterations were in intermediate- and low-risk genes (CHEK2, RAD50, RAD51, ATM, FANCA, FANCL, FANCC, BAP1), with controversial actionability: 8% harbored a somatic non-BRCA1/2 alteration, 6 cases had a high-risk alteration (PALB2, RAD51C), and one co-presented a PALB2/BRCA2 alteration. Elevated LOH was associated with HRR-mutated status and TP53 mutations while lowered LOH was associated with KRAS alterations. Including TMB/MSI data, the potential benefit from an NGS-oriented treatment was increased from 1.91% to 13.74% (high-MSI: 0.3%, TMB > 10 muts/MB: 12.78%). TMB was slightly increased in females (4.75 vs. 4.46 muts/MB) and in individuals with age > 60 (4.77 vs. 4.40 muts/MB). About 28.41% showed PD-L1 > 1% either in tumor or immune cells, 15.75% expressed PD-L1 ≥ 10%, and only 1.18% had PD-L1 ≥ 50%. This is the largest depiction of real-world genomic characteristics of European patients with PDAC, which offers some useful clinical and research insights.

Only 32.3% of Breast Cancer Families with Pathogenic Variants in Cancer Genes Utilized Cascade Genetic Testing.

BACKGROUND: Hereditary cancer predisposition syndromes are responsible for approximately 5-10% of all diagnosed cancer cases. In order to identify individuals at risk in a cost-efficient manner, family members of individuals carrying pathogenic alterations are tested only for the specific variant that was identified in their carrier relative. The purpose of this study was to investigate the clinical use and implementation of cascade family testing (CFT) in families of breast cancer patients with pathogenic/likely pathogenic variants (PVs/LPVs) in cancer-related predisposition genes. METHODS: Germline sequencing was carried out with NGS technology using a 52-gene panel, and cascade testing was performed by Sanger sequencing or MLPA. RESULTS: In a cohort of 1785 breast cancer patients (families), 20.3% were found to have PVs/LPVs. Specifically, 52.2%, 25.1%, and 22.7% of patients had positive findings in high-, intermediate-, and low-penetrance breast cancer susceptibility genes, respectively. Although CFT was recommended to all families, only 117 families (32.3%) agreed to proceed with genetic testing. Among the first-degree relatives who underwent CFT, 70.3% were female, and 108 of 121 (89.3%) were cancer free. Additionally, 42.7%, 36.7%, and 20.6% were offspring, siblings, and parents of the subject, respectively. Our data suggest that CFT was mostly undertaken (104/117, 88.8%) in families with positive findings in high-risk genes. CONCLUSIONS: Cascade family testing can be a powerful tool for primary cancer prevention by identifying at-risk family members. It is of utmost importance to implement genetic counseling approaches leading to increased awareness and communication of genetic testing results.

Copy Number Variations (CNVs) Account for 10.8% of Pathogenic Variants in Patients Referred for Hereditary Cancer Testing.

BACKGROUND/AIM: Germline copy number variation (CNV) is a type of genetic variant that predisposes significantly to inherited cancers. Today, next-generation sequencing (NGS) technologies have contributed to multi gene panel analysis in clinical practice. MATERIALS AND METHODS: A total of 2,163 patients were screened for cancer susceptibility, using a solution-based capture method. A panel of 52 genes was used for targeted NGS. The capture-based approach enables computational analysis of CNVs from NGS data. We studied the performance of the CNV module of the commercial software suite SeqPilot (JSI Medical Systems) and of the non-commercial tool panelcn.MOPS. Additionally, we tested the performance of digital multiplex ligation-dependent probe amplification (digitalMLPA). RESULTS: Pathogenic/likely pathogenic variants (P/LP) were identified in 464 samples (21.5%). CNV accounts for 10.8% (50/464) of pathogenic variants, referring to deletion/duplication of one or more exons of a gene. In patients with breast and ovarian cancer, CNVs accounted for 10.2% and 6.8% of pathogenic variants, respectively. In colorectal cancer patients, CNV accounted for 28.6% of pathogenic/likely pathogenic variants. CONCLUSION: In silico CNV detection tools provide a viable and cost-effective method to identify CNVs from NGS experiments. CNVs constitute a substantial percentage of P/LP variants, since they represent up to one of every ten P/LP findings identified by NGS multigene analysis; therefore, their evaluation is highly recommended to improve the diagnostic yield of hereditary cancer analysis.

Importance of multigene panel test in patients with consanguineous marriage and family history of breast cancer.

Next-generation sequencing (NGS) technology is used to evaluate hereditary cancer risks of patients worldwide; however, information concerning the germline multigene mutational spectrum among patients with breast cancer (BC) with consanguineous marriage (CM) is limited. Therefore, this prospective study aimed to determine the molecular characteristics of patients with BC who were tested with multigene hereditary cancer predisposition NGS panel and to show the effect of CM on cancer-related genes. Patients with BC with or without CM and family history (FH) of BC treated in our breast center were selected according to The National Comprehensive Cancer Network (NCCN) criteria for hereditary BC. In these patients, the analysis of a panel of 33 genes involved in hereditary cancer predisposition was performed after genetic counseling by using NGS. The pathogenic variant (PV) and the variant of uncertain significance (VUS) were found to be 15.8 and 47.4%, respectively. PVs were identified in 10/33 genes in 34 patients; 38.2% in BRCA1/2 genes; 6, 24, and 14% in other high, moderate and low-risk genes, respectively. The CM rate was 17.7% among the 215 patients with BC. The PV rate was 13.2% in patients with CM and 16.4% in patients without CM (P=0.80). When PV and VUS were evaluated together, the PV+VUS ratio was significantly higher in patients with CM and FH of BC than patients without CM and FH of BC (88.2 vs. 63.3%, P=0.045). Analysis of multigene panel provided 9.76% additional PVs in moderate/low-risk genes. The PV rate was similar in patients with BC with or without CM. A high PV+VUS ratio in patients with CM and FH of BC suggests that genes whose importance are unknown are likely to be pathogenic genes later.

Revisiting the Implications of Positive Germline Testing Results Using Multi-gene Panels in Breast Cancer Patients.

BACKGROUND/AIM: The use of multi-gene panels for germline testing in breast cancer enables the estimation of cancer risk and guides risk-reducing management options. The aim of this study was to present data that demonstrate the different levels of actionability for multi-gene panels used in genetic testing of breast cancer patients and their family members. MATERIALS AND METHODS: We performed an analysis in our clinical database to identify breast cancer patients undergoing genetic testing. We reviewed positive results in respect of risk estimation and management, cascade family testing, secondary findings and information for treatment decision-making. RESULTS: A total of 415 positive test reports were identified with 57.1%, 18.1%, 10.8% and 13.5% of individuals having pathogenic/likely pathogenic variants in high, moderate, low and with insufficient evidence for breast cancer risk genes, respectively. Six point seven percent of individuals were double heterozygotes. CONCLUSION: Germline findings in 92% of individuals are linked to evidence-based treatment information and risk estimates for predisposition to breast and/or other cancer types. The use of germline findings for treatment decision making expands the indication of genetic testing to include individuals that could benefit from targeted treatments.

Genetic Predisposition to Male Breast Cancer: A Case Series.

BACKGROUND/AIM: Male breast cancer (MBC) is a very rare disorder affecting approximately 1 in 833 men. Genetic predisposition is one of the most important risk factors of MBC with BRCA2 being the most commonly mutated gene in males diagnosed with breast cancer. However, a large part of MBC heritability is still unexplained. This study sought to add to the data already available on the genetics of MBC. MATERIALS AND METHODS: Our study initially involved comprehensive analysis of BRCA1 and BRCA2, followed by analysis of 43 genes implicated in cancer predisposition in a series of 100 Greek patients diagnosed with MBC between 1995-2015. RESULTS: Pathogenic variants were identified in 13 patients, with BRCA2 being the most commonly affected gene, followed by BRCA1, RAD50, RAD51B, and MSH3. CONCLUSION: In agreement with previous reports, BRCA2 is the most important genetic factor of MBC predisposition, while the remaining known cancer predisposition genes are each very rarely involved, rendering conclusions as to their cumulative effect difficult to draw.

ExTopoDB: a database of experimentally derived topological models of transmembrane proteins.

UNLABELLED: ExTopoDB is a publicly accessible database of experimentally derived topological models of transmembrane proteins. It contains information collected from studies in the literature that report the use of biochemical methods for the determination of the topology of α-helical transmembrane proteins. Transmembrane protein topology is highly important in order to understand their function and ExTopoDB provides an up to date, complete and comprehensive dataset of experimentally determined topologies of α-helical transmembrane proteins. Topological information is combined with transmembrane topology prediction resulting in more reliable topological models. AVAILABILITY: http://bioinformatics.biol.uoa.gr/ExTopoDB.

miRNA polymorphisms and risk of premature coronary artery disease.

OBJECTIVE: Several microRNA (miRNA) polymorphisms have been associated with susceptibility to specific health disorders, including cardiovascular diseases. The aim of the present study was to investigate whether four well-studied miRNA polymorphisms in non-Caucasian populations, namely miR146a G>C (rs2910164), miR149 C>T (rs2292832), miR196a2 C>T (rs11614913) and miR499 A>G (rs3746444), contribute to the risk for the development of premature Coronary Artery Disease (CAD) in the Greek population. METHODS: We used a case-control study to examine these associations in 400 individuals: 200 CAD patients [including a subgroup of myocardial infraction (MI) patients] and 200 healthy controls, all of Greek origin. MiRNA polymorphisms were genotyped using three different assays: Polymerase chain reaction - restriction fragment length polymorphism (PCR-RFLP), High resolution Melting (HRM) and Sanger sequencing. RESULTS: Two of these polymorphisms, miR196a2 C>T (rs11614913) and miR499 A>G (rs3746444) were found to be strongly associated with increased risk for CAD (p=0.0388 and p=0.0013, respectively) and for MI (p=0.0281 and p=0.0273, respectively). Furthermore, miR146C-miR149C-miR196T-miR499G allele combination appeared to be significantly related to CAD (p=0.0185) and MI (p=0.0337) prevalence. CONCLUSIONS: Our results suggest that at least two of the studied polymorphisms, miR196a2 C>T (rs11614913) and miR499 A>G (rs3746444), as well as the miR146C-miR149C-miR196T-miR499G allele combination could represent useful biomarkers of CAD and/or MI susceptibility in the Greek population. These special genetic characteristics, in combination with environmental factors and personal habits, might contribute to CAD and/or MI prevalence.

Clinical Utility of Functional RNA Analysis for the Reclassification of Splicing Gene Variants in Hereditary Cancer.

BACKGROUND: Classification of splicing variants (SVs) in genes associated with hereditary cancer is often challenging. The aim of this study was to investigate the occurrence of SVs in hereditary cancer genes and the clinical utility of RNA analysis. MATERIAL AND METHODS: 1518 individuals were tested for cancer predisposition, using a Next Generation Sequencing (NGS) panel of 36 genes. Splicing variant analysis was performed using RT-PCR and Sanger Sequencing. RESULTS: In total, 34 different SVs were identified, 53% of which were classified as pathogenic or likely pathogenic. The remaining 16 variants were initially classified as Variant of Uncertain Significance (VUS). RNA analysis was performed for 3 novel variants. CONCLUSION: The RNA analysis assisted in the reclassification of 20% of splicing variants from VUS to pathogenic. RNA analysis is essential in the case of uncharacterized splicing variants, for proper classification and personalized management of these patients.

PerMemDB: A database for eukaryotic peripheral membrane proteins.

The majority of all proteins in cells interact with membranes either permanently or temporarily. Peripheral membrane proteins form transient complexes with membrane proteins and/or lipids, via non-covalent interactions and are of outmost importance, due to numerous cellular functions in which they participate. In an effort to collect data regarding this heterogeneous group of proteins we designed and constructed a database, called PerMemDB. PerMemDB is currently the most complete and comprehensive repository of data for eukaryotic peripheral membrane proteins deposited in UniProt or predicted with the use of MBPpred - a computational method that specializes in the detection of proteins that interact non-covalently with membrane lipids, via membrane binding domains. The first version of the database contains 231,770 peripheral membrane proteins from 1009 organisms. All entries have cross-references to other databases, literature references and annotation regarding their interactions with other proteins. Moreover, additional sequence annotation of the characteristic domains that allow these proteins to interact with membranes is available, due to the application of MBPpred. Through the web interface of PerMemDB, users can browse the contents of the database, submit advanced text searches and BLAST queries against the protein sequences deposited in PerMemDB. We expect this repository to serve as a source of information that will allow the scientific community to gain a deeper understanding of the evolution and function of peripheral membrane proteins via the enhancement of proteome-wide analyses. The database is available at: http://bioinformatics.biol.uoa.gr/db=permemdb.