m6ACali: machine learning-powered calibration for accurate m6A detection in MeRIP-Seq.

We present m6ACali, a novel machine-learning framework aimed at enhancing the accuracy of N6-methyladenosine (m6A) epitranscriptome profiling by reducing the impact of non-specific antibody enrichment in MeRIP-Seq. The calibration model serves as a genomic feature-based classifier that refines the identification of m6A sites, distinguishing those genuinely present from those that can be detected in in-vitro transcribed (IVT) control experiments. We find that m6ACali effectively identifies non-specific binding peaks reported by exomePeak2 and MACS2 in novel MeRIP-Seq datasets without the need for paired IVT controls. The model interpretation revealed that off-target antibody binding sites commonly occur at short exons and short mRNAs, originating from high read coverage regions that share the motif sequence with true m6A sites. We also reveal that the ML strategy can efficiently adjust differentially methylated peaks and other antibody-dependent, base-resolution m6A detection techniques. As a result, m6ACali offers a promising method for the universal enhancement of m6A profiles generated by MeRIP-Seq experiments, elevating the benchmark for omics-level m6A data integration.

KH-like Domains in PARP9/DTX3L and PARP14 Coordinate Protein-Protein Interactions to Promote Cancer Cell Survival.

Certain members of the ADP-ribosyltransferase superfamily (ARTD or PARP enzymes) catalyse ADP-ribosylation in response to cellular stress, DNA damage and viral infection and are upregulated in various tumours. PARP9, its binding partner DTX3L and PARP14 protein levels are significantly correlated in head and neck squamous cell carcinoma (HNSCC) and other tumour types though a mechanism where PARP9/DTX3L regulates PARP14 post-transcriptionally. Depleting PARP9, DTX3L or PARP14 expression in HNSCC or HeLa cell lines decreases cell survival through a reduction of proliferation and an increase in apoptosis. A partial rescue of survival was achieved by expressing a PARP14 truncation containing a predicted eukaryotic type I KH domain. KH-like domains were also found in PARP9 and in DTX3L and contributed to protein-protein interactions between PARP9-DTX3L and PARP14-DTX3L. Homodimerization of DTX3L was also coordinated by a KH-like domain and was disrupted by site-specific mutation. Although, cell survival promoted by PARP14 did not require ADP-ribosyltransferase activity, interaction of DTX3L in vitro suppressed PARP14 auto-ADP-ribosylation and promoted trans-ADP-ribosylation of PARP9 and DTX3L. In summary, we characterised PARP9-DTX3L-PARP14 interactions important to pro-survival signalling in HNSCC cells, albeit in PARP14 catalytically independent fashion.

The 2023 Nucleic Acids Research Database Issue and the online molecular biology database collection.

The 2023 Nucleic Acids Research Database Issue contains 178 papers ranging across biology and related fields. There are 90 papers reporting on new databases and 82 updates from resources previously published in the Issue. Six more papers are updates from databases most recently published elsewhere. Major nucleic acid databases reporting updates include Genbank, ENA, ChIPBase, JASPAR, mirDIP and the Issue's first Breakthrough Article, NACDDB for Circular Dichroism data. Updates from BMRB and RCSB cover experimental protein structural data while AlphaFold 2 computational structure predictions feature widely. STRING and REBASE are stand-out updates in the signalling and enzymes section. Immunology-related databases include CEDAR, the second Breakthrough Article, for cancer epitopes and receptors alongside returning IPD-IMGT/HLA and the new PGG.MHC. Genomics-related resources include Ensembl, GWAS Central and UCSC Genome Browser. Major returning databases for drugs and their targets include Open Targets, DrugCentral, CTD and Pubchem. The EMPIAR image archive appears in the Issue for the first time. The entire database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). The NAR online Molecular Biology Database Collection has been updated, revisiting 463 entries, adding 92 new resources and eliminating 96 discontinued URLs so bringing the current total to 1764 databases. It is available at http://www.oxfordjournals.org/nar/database/c/.

m6A-TSHub: Unveiling the Context-specific m6A Methylation and m6A-affecting Mutations in 23 Human Tissues.

As the most pervasive epigenetic marker present on mRNA and lncRNA, N6-methyladenosine (m6A) RNA methylation has been shown to participate in essential biological processes. Recent studies have revealed the distinct patterns of m6A methylome across human tissues, and a major challenge remains in elucidating the tissue-specific presence and circuitry of m6A methylation. We present here a comprehensive online platform m6A-TSHub for unveiling the context-specific m6A methylation and genetic mutations that potentially regulate m6A epigenetic mark. m6A-TSHub consists of four core components, including: (1) m6A-TSDB, a comprehensive database of 184,554 functionally annotated m6A sites derived from 23 human tissues and 499,369 m6A sites from 25 tumor conditions, respectively; (2) m6A-TSFinder, a web server for high-accuracy prediction of m6A methylation sites within a specific tissue from RNA sequences, which was constructed using multi-instance deep neural networks with gated attention; (3) m6A-TSVar, a web server for assessing the impact of genetic variants on tissue-specific m6A RNA modifications; and (4) m6A-CAVar, a database of 587,983 The Cancer Genome Atlas (TCGA) cancer mutations (derived from 27 cancer types) that were predicted to affect m6A modifications in the primary tissue of cancers. The database should make a useful resource for studying the m6A methylome and the genetic factors of epitranscriptome disturbance in a specific tissue (or cancer type). m6A-TSHub is accessible at www.xjtlu.edu.cn/biologicalsciences/m6ats.

Bisphenol A Inhibits the Transporter Function of the Blood-Brain Barrier by Directly Interacting with the ABC Transporter Breast Cancer Resistance Protein (BCRP).

The breast cancer resistance protein (BCRP) is an important efflux transporter in the blood-brain barrier (BBB), protecting the brain from a wide range of substances. In this study, we investigated if BCRP function is affected by bisphenol A (BPA), a high production volume chemical used in common consumer products, as well as by bisphenol F (BPF) and bisphenol S (BPS), which are used to substitute BPA. We employed a transwell-based in vitro cell model of iPSC-derived brain microvascular endothelial cells, where BCRP function was assessed by measuring the intracellular accumulation of its substrate Hoechst 33342. Additionally, we used in silico modelling to predict if the bisphenols could directly interact with BCRP. Our results showed that BPA significantly inhibits the transport function of BCRP. Additionally, BPA was predicted to bind to the cavity that is targeted by known BCRP inhibitors. Taken together, our findings demonstrate that BPA inhibits BCRP function in vitro, probably by direct interaction with the transporter. This effect might contribute to BPA's known impact on neurodevelopment.

Evolution, structure and emerging roles of C1ORF112 in DNA replication, DNA damage responses, and cancer.

The C1ORF112 gene initially drew attention when it was found to be strongly co-expressed with several genes previously associated with cancer and implicated in DNA repair and cell cycle regulation, such as RAD51 and the BRCA genes. The molecular functions of C1ORF112 remain poorly understood, yet several studies have uncovered clues as to its potential functions. Here, we review the current knowledge on C1ORF112 biology, its evolutionary history, possible functions, and its potential relevance to cancer. C1ORF112 is conserved throughout eukaryotes, from plants to humans, and is very highly conserved in primates. Protein models suggest that C1ORF112 is an alpha-helical protein. Interestingly, homozygous knockout mice are not viable, suggesting an essential role for C1ORF112 in mammalian development. Gene expression data show that, among human tissues, C1ORF112 is highly expressed in the testes and overexpressed in various cancers when compared to healthy tissues. C1ORF112 has also been shown to have altered levels of expression in some tumours with mutant TP53. Recent screens associate C1ORF112 with DNA replication and reveal possible links to DNA damage repair pathways, including the Fanconi anaemia pathway and homologous recombination. These insights provide important avenues for future research in our efforts to understand the functions and potential disease relevance of C1ORF112.

The 2021 Nucleic Acids Research database issue and the online molecular biology database collection.

The 2021 Nucleic Acids Research database Issue contains 189 papers spanning a wide range of biological fields and investigation. It includes 89 papers reporting on new databases and 90 covering recent changes to resources previously published in the Issue. A further ten are updates on databases most recently published elsewhere. Seven new databases focus on COVID-19 and SARS-CoV-2 and many others offer resources for studying the virus. Major returning nucleic acid databases include NONCODE, Rfam and RNAcentral. Protein family and domain databases include COG, Pfam, SMART and Panther. Protein structures are covered by RCSB PDB and dispersed proteins by PED and MobiDB. In metabolism and signalling, STRING, KEGG and WikiPathways are featured, along with returning KLIFS and new DKK and KinaseMD, all focused on kinases. IMG/M and IMG/VR update in the microbial and viral genome resources section, while human and model organism genomics resources include Flybase, Ensembl and UCSC Genome Browser. Cancer studies are covered by updates from canSAR and PINA, as well as newcomers CNCdatabase and Oncovar for cancer drivers. Plant comparative genomics is catered for by updates from Gramene and GreenPhylDB. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). The NAR online Molecular Biology Database Collection has been substantially updated, revisiting nearly 1000 entries, adding 90 new resources and eliminating 86 obsolete databases, bringing the current total to 1641 databases. It is available at https://www.oxfordjournals.org/nar/database/c/.

m6A-Atlas: a comprehensive knowledgebase for unraveling the N6-methyladenosine (m6A) epitranscriptome.

N 6-Methyladenosine (m6A) is the most prevalent RNA modification on mRNAs and lncRNAs. It plays a pivotal role during various biological processes and disease pathogenesis. We present here a comprehensive knowledgebase, m6A-Atlas, for unraveling the m6A epitranscriptome. Compared to existing databases, m6A-Atlas features a high-confidence collection of 442 162 reliable m6A sites identified from seven base-resolution technologies and the quantitative (rather than binary) epitranscriptome profiles estimated from 1363 high-throughput sequencing samples. It also offers novel features, such as; the conservation of m6A sites among seven vertebrate species (including human, mouse and chimp), the m6A epitranscriptomes of 10 virus species (including HIV, KSHV and DENV), the putative biological functions of individual m6A sites predicted from epitranscriptome data, and the potential pathogenesis of m6A sites inferred from disease-associated genetic mutations that can directly destroy m6A directing sequence motifs. A user-friendly graphical user interface was constructed to support the query, visualization and sharing of the m6A epitranscriptomes annotated with sites specifying their interaction with post-transcriptional machinery (RBP-binding, microRNA interaction and splicing sites) and interactively display the landscape of multiple RNA modifications. These resources provide fresh opportunities for unraveling the m6A epitranscriptomes. m6A-Atlas is freely accessible at: www.xjtlu.edu.cn/biologicalsciences/atlas.

RMDisease: a database of genetic variants that affect RNA modifications, with implications for epitranscriptome pathogenesis.

Deciphering the biological impacts of millions of single nucleotide variants remains a major challenge. Recent studies suggest that RNA modifications play versatile roles in essential biological mechanisms, and are closely related to the progression of various diseases including multiple cancers. To comprehensively unveil the association between disease-associated variants and their epitranscriptome disturbance, we built RMDisease, a database of genetic variants that can affect RNA modifications. By integrating the prediction results of 18 different RNA modification prediction tools and also 303,426 experimentally-validated RNA modification sites, RMDisease identified a total of 202,307 human SNPs that may affect (add or remove) sites of eight types of RNA modifications (m6A, m5C, m1A, m5U, Ψ, m6Am, m7G and Nm). These include 4,289 disease-associated variants that may imply disease pathogenesis functioning at the epitranscriptome layer. These SNPs were further annotated with essential information such as post-transcriptional regulations (sites for miRNA binding, interaction with RNA-binding proteins and alternative splicing) revealing putative regulatory circuits. A convenient graphical user interface was constructed to support the query, exploration and download of the relevant information. RMDisease should make a useful resource for studying the epitranscriptome impact of genetic variants via multiple RNA modifications with emphasis on their potential disease relevance. RMDisease is freely accessible at: www.xjtlu.edu.cn/biologicalsciences/rmd.

Maturation of the functional mouse CRES amyloid from globular form.

The epididymal lumen contains a complex cystatin-rich nonpathological amyloid matrix with putative roles in sperm maturation and sperm protection. Given our growing understanding for the biological function of this and other functional amyloids, the problem still remains: how functional amyloids assemble including their initial transition to early oligomeric forms. To examine this, we developed a protocol for the purification of nondenatured mouse CRES, a component of the epididymal amyloid matrix, allowing us to examine its assembly to amyloid under conditions that may mimic those in vivo. Herein we use X-ray crystallography, solution-state NMR, and solid-state NMR to follow at the atomic level the assembly of the CRES amyloidogenic precursor as it progressed from monomeric folded protein to an advanced amyloid. We show the CRES monomer has a typical cystatin fold that assembles into highly branched amyloid matrices, comparable to those in vivo, by forming ß-sheet assemblies that our data suggest occur via two distinct mechanisms: a unique conformational switch of a highly flexible disulfide-anchored loop to a rigid ß-strand and by traditional cystatin domain swapping. Our results provide key insight into our understanding of functional amyloid assembly by revealing the earliest structural transitions from monomer to oligomer and by showing that some functional amyloid structures may be built by multiple and distinctive assembly mechanisms.

WHISTLE: a high-accuracy map of the human N6-methyladenosine (m6A) epitranscriptome predicted using a machine learning approach.

N 6-methyladenosine (m6A) is the most prevalent post-transcriptional modification in eukaryotes, and plays a pivotal role in various biological processes, such as splicing, RNA degradation and RNA-protein interaction. We report here a prediction framework WHISTLE for transcriptome-wide m6A RNA-methylation site prediction. When tested on six independent datasets, our approach, which integrated 35 additional genomic features besides the conventional sequence features, achieved a major improvement in the accuracy of m6A site prediction (average AUC: 0.948 and 0.880 under the full transcript or mature messenger RNA models, respectively) compared to the state-of-the-art computational approaches MethyRNA (AUC: 0.790 and 0.732) and SRAMP (AUC: 0.761 and 0.706). It also out-performed the existing epitranscriptome databases MeT-DB (AUC: 0.798 and 0.744) and RMBase (AUC: 0.786 and 0.736), which were built upon hundreds of epitranscriptome high-throughput sequencing samples. To probe the putative biological processes impacted by changes in an individual m6A site, a network-based approach was implemented according to the 'guilt-by-association' principle by integrating RNA methylation profiles, gene expression profiles and protein-protein interaction data. Finally, the WHISTLE web server was built to facilitate the query of our high-accuracy map of the human m6A epitranscriptome, and the server is freely available at: www.xjtlu.edu.cn/biologicalsciences/whistle and http://whistle-epitranscriptome.com.

Critical assessment of approaches for molecular docking to elucidate associations of HLA alleles with adverse drug reactions.

Adverse drug reactions have been linked with genetic polymorphisms in HLA genes in numerous different studies. HLA proteins have an essential role in the presentation of self and non-self peptides, as part of the adaptive immune response. Amongst the associated drugs-allele combinations, anti-HIV drug Abacavir has been shown to be associated with the HLA-B*57:01 allele, and anti-epilepsy drug Carbamazepine with B*15:02, in both cases likely following the altered peptide repertoire model of interaction. Under this model, the drug binds directly to the antigen presentation region, causing different self peptides to be presented, which trigger an unwanted immune response. There is growing interest in searching for evidence supporting this model for other ADRs using bioinformatics techniques. In this study, in silico docking was used to assess the utility and reliability of well-known docking programs when addressing these challenging HLA-drug situations. The overall aim was to address the uncertainty of docking programs giving different results by completing a detailed comparative study of docking software, grounded in the MHC-ligand experimental structural data - for Abacavir and to a lesser extent Carbamazepine - in order to assess their performance. Four docking programs: SwissDock, ROSIE, AutoDock Vina and AutoDockFR, were used to investigate if each software could accurately dock the Abacavir back into the crystal structure for the protein arising from the known risk allele, and if they were able to distinguish between the HLA-associated and non-HLA-associated (control) alleles. The impact of using homology models on the docking performance and how using different parameters, such as including receptor flexibility, affected the docking performance were also investigated to simulate the approach where a crystal structure for a given HLA allele may be unavailable. The programs that were best able to predict the binding position of Abacavir were then used to recreate the docking seen for Carbamazepine with B*15:02 and controls alleles. It was found that the programs investigated were sometimes able to correctly predict the binding mode of Abacavir with B*57:01 but not always. Each of the software packages that were assessed could predict the binding of Abacavir and Carbamazepine within the correct sub-pocket and, with the exception of ROSIE, was able to correctly distinguish between risk and control alleles. We found that docking to homology models could produce poorer quality predictions, especially when sequence differences impact the architecture of predicted binding pockets. Caution must therefore be used as inaccurate structures may lead to erroneous docking predictions. Incorporating receptor flexibility was found to negatively affect the docking performance for the examples investigated. Taken together, our findings help characterise the potential but also the limitations of computational prediction of drug-HLA interactions. These docking techniques should therefore always be used with care and alongside other methods of investigation, in order to be able to draw strong conclusions from the given results.

Large scale ab initio modeling of structurally uncharacterized antimicrobial peptides reveals known and novel folds.

Antimicrobial resistance within a wide range of infectious agents is a severe and growing public health threat. Antimicrobial peptides (AMPs) are among the leading alternatives to current antibiotics, exhibiting broad spectrum activity. Their activity is determined by numerous properties such as cationic charge, amphipathicity, size, and amino acid composition. Currently, only around 10% of known AMP sequences have experimentally solved structures. To improve our understanding of the AMP structural universe we have carried out large scale ab initio 3D modeling of structurally uncharacterized AMPs that revealed similarities between predicted folds of the modeled sequences and structures of characterized AMPs. Two of the peptides whose models matched known folds are Lebocin Peptide 1A (LP1A) and Odorranain M, predicted to form ß-hairpins but, interestingly, to lack the intramolecular disulfide bonds, cation-π or aromatic interactions that generally stabilize such AMP structures. Other examples include Ponericin Q42, Latarcin 4a, Kassinatuerin 1, Ceratotoxin D, and CPF-B1 peptide, which have α-helical folds, as well as mixed αß folds of human Histatin 2 peptide and Garvicin A which are, to the best of our knowledge, the first linear αßß fold AMPs lacking intramolecular disulfide bonds. In addition to fold matches to experimentally derived structures, unique folds were also obtained, namely for Microcin M and Ipomicin. These results help in understanding the range of protein scaffolds that naturally bear antimicrobial activity and may facilitate protein design efforts towards better AMPs.

The 24th annual Nucleic Acids Research database issue: a look back and upcoming changes.

This year's Database Issue of Nucleic Acids Research contains 152 papers that include descriptions of 54 new databases and update papers on 98 databases, of which 16 have not been previously featured in NAR As always, these databases cover a broad range of molecular biology subjects, including genome structure, gene expression and its regulation, proteins, protein domains, and protein-protein interactions. Following the recent trend, an increasing number of new and established databases deal with the issues of human health, from cancer-causing mutations to drugs and drug targets. In accordance with this trend, three recently compiled databases that have been selected by NAR reviewers and editors as 'breakthrough' contributions, denovo-db, the Monarch Initiative, and Open Targets, cover human de novo gene variants, disease-related phenotypes in model organisms, and a bioinformatics platform for therapeutic target identification and validation, respectively. We expect these databases to attract the attention of numerous researchers working in various areas of genetics and genomics. Looking back at the past 12 years, we present here the 'golden set' of databases that have consistently served as authoritative, comprehensive, and convenient data resources widely used by the entire community and offer some lessons on what makes a successful database. The Database Issue is freely available online at the https://academic.oup.com/nar web site. An updated version of the NAR Molecular Biology Database Collection is available at http://www.oxfordjournals.org/nar/database/a/.

The centrosomal deubiquitylase USP21 regulates Gli1 transcriptional activity and stability.

USP21 is a centrosome-associated deubiquitylase (DUB) that has been implicated in the formation of primary cilia - crucial organelles for the regulation of the Hedgehog (Hh) signaling pathway in vertebrates. Here, we identify KCTD6 - a cullin-3 E3-ligase substrate adapter that has been previously linked to Hh signaling - as well as Gli1, the key transcription factor responsible for Hh signal amplification, as new interacting partners of USP21. We identify a cryptic structured protein interaction domain in KCTD6, which is predicted to have a similar fold to Smr domains. Importantly, we show that both depletion and overexpression of catalytically active USP21 suppress Gli1-dependent transcription. Gli proteins are negatively regulated through protein kinase A (PKA)-dependent phosphorylation. We provide evidence that USP21 recruits and stabilises Gli1 at the centrosome where it promotes its phosphorylation by PKA. By revealing an intriguing functional pairing between a spatially restricted deubiquitylase and a kinase, our study highlights the centrosome as an important hub for signal coordination.

The 2016 database issue of Nucleic Acids Research and an updated molecular biology database collection.

The 2016 Database Issue of Nucleic Acids Research starts with overviews of the resources provided by three major bioinformatics centers, the U.S. National Center for Biotechnology Information (NCBI), the European Bioinformatics Institute (EMBL-EBI) and Swiss Institute for Bioinformatics (SIB). Also included are descriptions of 62 new databases and updates on 95 databases that have been previously featured in NAR plus 17 previously described elsewhere. A number of papers in this issue deal with resources on nucleic acids, including various kinds of non-coding RNAs and their interactions, molecular dynamics simulations of nucleic acid structure, and two databases of super-enhancers. The protein database section features important updates on the EBI's Pfam, PDBe and PRIDE databases, as well as a variety of resources on pathways, metabolomics and metabolic modeling. This issue also includes updates on popular metagenomics resources, such as MG-RAST, EBI Metagenomics, and probeBASE, as well as a newly compiled Human Pan-Microbe Communities database. A significant fraction of the new and updated databases are dedicated to the genetic basis of disease, primarily cancer, and various aspects of drug research, including resources for patented drugs, their side effects, withdrawn drugs, and potential drug targets. A further six papers present updated databases of various antimicrobial and anticancer peptides. The entire Database Issue is freely available online on the Nucleic Acids Research website (http://nar.oxfordjournals.org/). The NAR online Molecular Biology Database Collection, http://www.oxfordjournals.org/nar/database/c/, has been updated with the addition of 88 new resources and removal of 23 obsolete websites, which brought the current listing to 1685 databases.

The 2015 Nucleic Acids Research Database Issue and molecular biology database collection.

The 2015 Nucleic Acids Research Database Issue contains 172 papers that include descriptions of 56 new molecular biology databases, and updates on 115 databases whose descriptions have been previously published in NAR or other journals. Following the classification that has been introduced last year in order to simplify navigation of the entire issue, these articles are divided into eight subject categories. This year's highlights include RNAcentral, an international community portal to various databases on noncoding RNA; ValidatorDB, a validation database for protein structures and their ligands; SASBDB, a primary repository for small-angle scattering data of various macromolecular complexes; MoonProt, a database of 'moonlighting' proteins, and two new databases of protein-protein and other macromolecular complexes, ComPPI and the Complex Portal. This issue also includes an unusually high number of cancer-related databases and other databases dedicated to genomic basics of disease and potential drugs and drug targets. The size of NAR online Molecular Biology Database Collection, http://www.oxfordjournals.org/nar/database/a/, remained approximately the same, following the addition of 74 new resources and removal of 77 obsolete web sites. The entire Database Issue is freely available online on the Nucleic Acids Research web site (http://nar.oxfordjournals.org/).

Conversion of a disulfide bond into a thioacetal group during echinomycin biosynthesis.

Echinomycin is a nonribosomal depsipeptide natural product with a range of interesting bioactivities that make it an important target for drug discovery and development. It contains a thioacetal bridge, a unique chemical motif derived from the disulfide bond of its precursor antibiotic triostin A by the action of an S-adenosyl-L-methionine-dependent methyltransferase, Ecm18. The crystal structure of Ecm18 in complex with its reaction products S-adenosyl-L-homocysteine and echinomycin was determined at 1.50 Å resolution. Phasing was achieved using a new molecular replacement package called AMPLE, which automatically derives search models from structure predictions based on ab initio protein modelling. Structural analysis indicates that a combination of proximity effects, medium effects, and catalysis by strain drives the unique transformation of the disulfide bond into the thioacetal linkage.

The first structure in a family of peptidase inhibitors reveals an unusual Ig-like fold.

We report the crystal structure solution of the Intracellular Protease Inhibitor (IPI) protein from Bacillus subtilis, which has been reported to be an inhibitor of the intracellular subtilisin Isp1 from the same organism. The structure of IPI is a variant of the all-beta, immunoglobulin (Ig) fold. It is possible that IPI is important for protein-protein interactions, of which inhibition of Isp1 is one. The intracellular nature of ISP is questioned, because an alternative ATG codon in the ipi gene would produce a protein with an N-terminal extension containing a signal peptide. It is possible that alternative initiation exists, producing either an intracellular inhibitor or a secreted form that may be associated with the cell surface.  Homologues of the IPI protein from other species are multi-domain proteins, containing signal peptides and domains also associated with the bacterial cell-surface. The cysteine peptidase inhibitors chagasin and amoebiasin also have Ig-like folds, but their topology differs significantly from that of IPI, and they share no recent common ancestor. A model of IPI docked to Isp1 shows similarities to other subtilisin:inhibitor complexes, particularly where the inhibitor interacts with the peptidase active site.

Prediction of antimicrobial peptides based on the adaptive neuro-fuzzy inference system application.

Antimicrobial peptides (AMPs) are widely distributed defense molecules and represent a promising alternative for solving the problem of antibiotic resistance. Nevertheless, the experimental time required to screen putative AMPs makes computational simulations based on peptide sequence analysis and/or molecular modeling extremely attractive. Artificial intelligence methods acting as simulation and prediction tools are of great importance in helping to efficiently discover and design novel AMPs. In the present study, state-of-the-art published outcomes using different prediction methods and databases were compared to an adaptive neuro-fuzzy inference system (ANFIS) model. Data from our study showed that ANFIS obtained an accuracy of 96.7% and a Matthew's Correlation Coefficient (MCC) of0.936, which proved it to be an efficient model for pattern recognition in antimicrobial peptide prediction. Furthermore, a lower number of input parameters were needed for the ANFIS model, improving the speed and ease of prediction. In summary, due to the fuzzy nature ofAMP physicochemical properties, the ANFIS approach presented here can provide an efficient solution for screening putative AMP sequences and for exploration of properties characteristic of AMPs.