OBJECTIVE: This study aimed to identify by in silico methods tropomyosin consensus B and T epitopes of shrimp species, house dust mites, insects, and nematodes associated with allergic diseases in tropical countries.
METHODS: In silico analysis included tropomyosin from mites (Der p 10, Der f 10, Blo t 10), insects (Aed a 10, Per a 7, Bla g 7), shrimp (Lit v 1, Pen m 1, Pen a 1), and nematode (Asc l 3) all sequences were taken from the UniProt database. Linear IgE epitopes were predicted with AlgPred 2.0 and validated with BepiPred 3.0. MHC-II binding T cell epitopes were predicted using the IEDB server, which implements nine predictive methods (consensus method, combinatorial library, NN-align-2.3, NN- align-2.2, SMM-align, Sturniolo, NetMHCIIpan 3.1, and NetMHCIIpan 3.2) these predictions focused on 10 HLA-DR and 2 HLA-DQ alleles associated with allergic diseases. Subsequently, consensus B and T epitopes present in all species were identified.
RESULTS: We identified 12 sequences that behaved as IgE-epitopes and B-cell epitopes, three of them: 160RKYDEVARKLAMVEA174, 192ELEEELRVVGNNLKSLEVSEEKAN215, 251KEVDRLEDELV261 were consensus in all species. Eleven peptides (T-epitopes) showed strong binding (percentile rank ≤ 2.0) to HLA-DRB1*0301, *0402, *0411, *0701, *1101, *1401, HLA-DQA1*03:01/DQB1*03:02, and HLA- DQA1*05:01/DQB1*02:01. Only two T-epitopes were consensus in all species: 167RKLAMVEADLERAEERAEt GEsKIVELEEELRV199, and 218EEeY KQQIKT LTaKLKEAEARAEFAERSV246. Subsequently, we identified 2 B and T epitope sequences and reached a consensus between species 167RKLAMVEA174 and 192ELEEELRV199.
CONCLUSIONS: These data describe three sequences that may explain the IgE cross-reactivity between the analyzed species. In addition, the consensus B and T epitopes can be used for further in vitro investigations and may help to design multiple-epitope protein-based immunotherapy for tropomyosin-related allergic diseases.
OBJECTIVE: Este estudio tuvo como objetivo identificar mediante métodos in silico epítopes B y T consenso de tropomiosina de especies de camarón, ácaros del polvo doméstico, insectos y nematodos asociados a enfermedades alérgicas en países tropicales.
UNASSIGNED: El análisis in silico incluyó tropomiosina de ácaros (Der p 10, Der f 10, Blo t 10), insectos (Aed a 10, Per a 7, Bla g 7), camarones (Lit v 1, Pen m 1, Pen a 1), y nematodo (Asc l 3). Todas las secuencias se tomaron de la base de datos UniProt. Los epítopes IgE lineales se predijeron con AlgPred 2.0 y se validaron con BepiPred 3.0. Los epítopes de células T de unión a MHC-II se predijeron utilizando el servidor IEDB, que implementa nueve métodos predictivos (método de consenso, biblioteca combinatoria, NN-align-2.3, NN-align-2.2, SMM-align, Sturniolo, NetMHCIIpan 3.1 y NetMHCIIpan 3.2). Estas predicciones se centraron en diez alelos HLA-DR y 2 HLA-DQ asociados con enfermedades alérgicas. Posteriormente, se identificaron epítopes consenso B y T presentes en todas las especies.
RESULTS: Se identificaron 12 secuencias que se comportaron como epítopes de IgE y, también, como epítopes de células B. Tres de ellas: 160RKYDEVARKLAMVEA174, 192ELEEELRVVGNNLKSLEVSEEKAN213 y 251KEVDRLEDELV261, fueron consenso en todas las especies. Once péptidos mostraron una fuerte unión (rango percentil ≤ 2,0) a HLA-DRB1*0301, *0402, *0411, *0701, *1101, *1401 y a HLA HLA-DQA1*03:01/DQB1*03:02, o HLA-DQA1*05:01/DQB1*02:01. Solo se encontraron dos secuencias: 167RKLAMVEADLERAEERAEtGEsKIVELEEELRV199 con fuerte afinidad por HLA-DQA1*03:01/DQB1*03:02, y HLA-DQA1*05:01/DQB1*02:01. Se identificaron dos secuencias que son epítopos B y T, y son consenso entre especies: 167RKLAMVEA174 y 192ELEEELRV199.
CONCLUSIONS: Estos datos describen tres secuencias que pueden explicar la reactividad cruzada de IgE entre las especies analizadas. Además, los epítopos B y T consenso se pueden usar para investigaciones in vitro adicionales, y pueden ayudar a diseñar inmunoterapia basada en proteínas de múltiepítopes para enfermedades alérgicas relacionadas con la tropomiosina.

Understanding the actual work (i.e., \"work-as-done\") rather than theorized work (i.e., \"work-as-imagined\") during complex medical processes is critical for developing approaches that improve patient outcomes. Although process mining has been used to discover process models from medical activity logs, it often omits critical steps or produces cluttered and unreadable models. In this paper, we introduce a TraceAlignment-based ProcessDiscovery method called TAD Miner to build interpretable process models for complex medical processes. TAD Miner creates simple linear process models using a threshold metric that optimizes the consensus sequence to represent the backbone process, and then identifies both concurrent activities and uncommon-but-critical activities to represent the side branches. TAD Miner also identifies the locations of repeated activities, an essential feature for representing medical treatment steps. We conducted a study using activity logs of 308 pediatric trauma resuscitations to develop and evaluate TAD Miner. TAD Miner was used to discover process models for five resuscitation goals, including establishing intravenous (IV) access, administering non-invasive oxygenation, performing back assessment, administering blood transfusion, and performing intubation. We quantitively evaluated the process models with several complexity and accuracy metrics, and performed qualitative evaluation with four medical experts to assess the accuracy and interpretability of the discovered models. Through these evaluations, we compared the performance of our method to that of two state-of-the-art process discovery algorithms: Inductive Miner and Split Miner. The process models discovered by TAD Miner had lower complexity and better interpretability than the state-of-the-art methods, and the fitness and precision of the models were comparable. We used the TAD process models to identify (1) the errors and (2)the best locations for the tentative steps in knowledge-driven expert models. The knowledge-driven models were revised based on the modifications suggested by the discovered models. The improved modeling using TAD Miner may enhance understanding of complex medical processes.

The freshwater snail Physella acuta was selected to expand the perspective of comparative snail immunology. Analysis of Physella acuta, belonging to the Physidae, taxonomic sister family to Planorbidae, affords family-level comparison of immune features characterized from Biomphalaria glabrata, the model snail often used to interpret general gastropod immunity. To capture constitutive and induced immune sequences, transcriptomes of an individual Physella acuta snail, 12 h post injection with bacteria (Gram -/+) and one sham-exposed snail were recorded with 454 pyrosequencing. Assembly yielded a combined reference transcriptome containing 24,288 transcripts. Additionally, genomic Illumina reads were obtained (∼15-fold coverage). Recovery of transcripts for two macin-like antimicrobial peptides (AMPs), 12 aplysianins, four LBP/BPIs and three physalysins indicated that Physella acuta shares a similar organization of antimicrobial defenses with Biomphalaria glabrata, contrasting a modest AMP arsenal with a diverse set of antimicrobial proteins. The lack of predicted transmembrane domains in all seven Physella acuta PGRP transcripts supports the notion that gastropods do not employ cell-bound PGRP receptors, different from ecdysozoan invertebrates yet similar to mammals (vertebrate deuterostomes). The well-documented sequence diversification by Biomphalaria glabrata FREPs (immune lectins comprising immunoglobulin superfamily domains and fibrinogen domains), resulting from somatic mutations of a large FREP gene family is hypothesized to be unique to Planorbidae; Physella acuta revealed just two bonafide FREP genes and these were not diversified. Furthermore, the flatworm parasite Echinostoma paraensei, confirmed here to infect both snail species, did not evoke from Physella acuta the abundant expression of FREP proteins at 2, 4 and 8 days post exposure that was previously observed from Biomphalaria glabrata. The Physella acuta reference transcriptome also revealed 24 unique transcripts encoding proteins consisting of a single fibrinogen-related domain (FReDs), with a short N-terminal sequence encoding either a signal peptide, transmembrane domain or no predicted features. The Physella acuta FReDs are candidate immune genes based on implication of similar sequences in immunity of bivalve molluscs. Overall, comparative analysis of snails of sister families elucidated the potential for taxon-specific immune features and investigation of strategically selected species will provide a more comprehensive view of gastropod immunity.

BACKGROUND: The current standard care therapy for hepatitis C virus (HCV) infection consists of two regimes, namely interferon-based and interferon-free treatments. The treatment through the combination of ribavirin and pegylated interferon is expensive, only mildly effective, and is associated with severe side effects. In 2011, two direct-acting antiviral (DAA) drugs, boceprevir and telaprevir, were licensed that have shown enhanced sustained virologic response (SVR) in phase III clinical trial, however, these interferon-free treatments are more sensitive to HCV genotype 1 infection. The variable nature of HCV, and the limited number of inhibitors developed thus aim in expanding the repertoire of available drug targets, resulting in targeting the virus assembly therapeutically.
OBJECTIVE: We conducted this study to predict the 3D structure of the p7 protein from the HCV genotypes 3 and 4. Approximately 63 amino acid residues encoded in HCV render this channel sensitive to inhibitors, making p7 a promising target for novel therapies. HCV p7 protein forms a small membrane known as viroporin, and is essential for effective self-assembly of large channels that conduct cation assembly and discharge infectious virion particles.
METHODS: In this study, we screened drugs and flavonoids known to disrupt translation and production of HCV proteins, targeted against the active site of p7 residues of HCV genotype 3 (GT3) (isolatek3a) and HCV genotype 4a (GT4) (isolateED43). Furthermore, we conducted a quantitative structure-activity relationship and docking interaction study.
RESULTS: The drug NB-DNJ formed the highest number of hydrogen bond interactions with both modeled p7 proteins with high interaction energy, followed by BIT225. A flavonoid screen demonstrated that Epigallocatechin gallate (EGCG), nobiletin, and quercetin, have more binding modes in GT3 than in GT4. Thus, the predicted p7 protein molecule of HCV from GT3 and GT4 provides a general avenue to target structure-based antiviral compounds.
CONCLUSIONS: We hypothesize that the inhibitors of viral p7 identified in this screen may be a new class of potent agents, but further confirmation in vitro and in vivo is essential. This structure-guided drug design for both GT3 and GT4 can lead to the identification of drug-like natural compounds, confirming p7 as a new target in the rapidly increasing era of HCV.

Phage-encoded cell wall peptidoglycan hydrolyzing enzymes, called endolysins, are essential for efficient release of virions from bacteria, and show species-specific killing of the host. We have demonstrated previously that the interaction between N-terminal catalytic and C-terminal cell wall binding domains of mycobacteriophage D29 endolysin makes the enzyme inactive in Escherichiacoli. Here, we demonstrate that such interaction occurs intramolecularly and is facilitated by a charged linker that connects the two domains. We also show that linker composition is crucial for the inactivation of PG hydrolase in E. coli. Such knowledge will immensely help in bioengineering of endolysins with narrow or broad spectrum antimicrobial activity.

Human short-chain dehydrogenases/reductases (SDRs) protein family has been the subject of recent studies for its critical role in human metabolism. Studies also found that single nucleotide polymorphisms of the SDR protein family were responsible for a variety of genetic diseases, including type II diabetes. This study reports the effect of sequence variation on the structural and functional integrities of human SDR protein family using phylogenetics and correlated mutation analysis tools. Our results indicated that (i) tyrosine, serine, and lysine are signature protein residues that have direct contribution to the structural and functional stabilities of the SDR protein family, (ii) subgroups of SDR protein family have their own signature protein combination that represent their unique functionality, and (iii) mutations of the human SDR protein family showed high correlation in terms of evolutionary history. In combination, the results inferred that over evolutionary history, the SDR protein family was able to diverge itself in order to adapt with the changes in human nutritional demands. Our study reveals understanding of structural and functional scaffolds of specific SDR subgroups that may facilitate the design of specific inhibitor.

The proteolytic activation of protein kinase Cδ (PKCδ) generates a catalytic fragment called PKCδ-CF, which induces cell death. However, the mechanisms underlying PKCδ-CF-mediated cell death are largely unknown. On the basis of an engineering leukemic cell line with inducible expression of PKCδ-CF, here we employ SILAC-based quantitative phosphoproteomics to systematically and dynamically investigate the overall phosphorylation events during cell death triggered by PKCδ-CF expression. Totally, 3000 phosphorylation sites were analyzed. Considering the fact that early responses to PKCδ-CF expression initiate cell death, we sought to identify pathways possibly related directly with PKCδ by further analyzing the data set of phosphorylation events that occur in the initiation stage of cell death. Interacting analysis of this data set indicates that PKCδ-CF triggers complicated networks to initiate cell death, and motif analysis and biochemistry verification reveal that several kinases in the downstream of PKCδ conduct these networks. By analysis of the specific sequence motif of kinase-substrate, we also find 59 candidate substrates of PKCδ from the up-regulated phosphopeptides, of which 12 were randomly selected for in vitro kinase assay and 9 were consequently verified as substrates of PKCδ. To our greatest understanding, this study provides the most systematic analysis of phosphorylation events initiated by the cleaved activated PKCδ, which would vastly extend the profound understanding of PKCδ-directed signal pathways in cell death. The MS data have been deposited to the ProteomeXchange with identifier PXD000225.

Nowadays, prokaryotic genomes are sequenced faster than the capacity to manually curate gene annotations. Automated genome annotation engines provide users a straight-forward and complete solution for predicting ORF coordinates and function. For many labs, the use of AGEs is therefore essential to decrease the time necessary for annotating a given prokaryotic genome. However, it is not uncommon for AGEs to provide different and sometimes conflicting predictions. Combining multiple AGEs might allow for more accurate predictions. Here we analyzed the ab initio open reading frame (ORF) calling performance of different AGEs based on curated genome annotations of eight strains from different bacterial species with GC% ranging from 35-52%. We present a case study which demonstrates a novel way of comparative genome annotation, using combinations of AGEs in a pre-defined order (or path) to predict ORF start codons. The order of AGE combinations is from high to low specificity, where the specificity is based on the eight genome annotations. For each AGE combination we are able to derive a so-called projected confidence value, which is the average specificity of ORF start codon prediction based on the eight genomes. The projected confidence enables estimating likeliness of a correct prediction for a particular ORF start codon by a particular AGE combination, pinpointing ORFs notoriously difficult to predict start codons. We correctly predict start codons for 90.5±4.8% of the genes in a genome (based on the eight genomes) with an accuracy of 81.1±7.6%. Our consensus-path methodology allows a marked improvement over majority voting (9.7±4.4%) and with an optimal path ORF start prediction sensitivity is gained while maintaining a high specificity.

Neph molecules are highly conserved immunoglobulin superfamily proteins (IgSF) which are essential for multiple morphogenetic processes, including glomerular development in mammals and neuronal as well as nephrocyte development in D. melanogaster. While D. melanogaster expresses two Neph-like proteins (Kirre and IrreC/Rst), three Neph proteins (Neph1-3) are expressed in the mammalian system. However, although these molecules are highly abundant, their molecular functions are still poorly understood. Here we report on a fly system in which we overexpress and replace endogenous Neph homologs with mammalian Neph1-3 proteins to identify functional Neph protein networks required for neuronal and nephrocyte development. Misexpression of Neph1, but neither Neph2 nor Neph3, phenocopies the overexpression of endogenous Neph molecules suggesting a functional diversity of mammalian Neph family proteins. Moreover, structure-function analysis identified a conserved and specific Neph1 protein motif that appears to be required for the functional replacement of Kirre. Hereby, we establish D. melanogaster as a genetic system to specifically model molecular Neph1 functions in vivo and identify a conserved amino acid motif linking Neph1 to Drosophila Kirre function.

BACKGROUND: The characterization of HIV-1 transmission strains may inform the design of an effective vaccine. Shorter variable loops with fewer predicted glycosites have been suggested as signatures enriched in envelope sequences derived during acute HIV-1 infection. Specifically, a transmission-linked lack of glycosites within the V1 and V2 loops of gp120 provides greater access to an α4β7 binding motif, which promotes the establishment of infection. Also, a histidine at position 12 in the leader sequence of Env has been described as a transmission signature that is selected against during chronic infection. The purpose of this study is to measure the association of the presence of an α4β7 binding motif, the number of N-linked glycosites, the length of the variable loops, and the prevalence of histidine at position 12 with HIV-1 transmission. A case-control study design was used to measure the prevalence of these variables between subtype B and C transmission sequences and frequency-matched randomly-selected sequences derived from chronically infected controls.
RESULTS: Subtype B transmission strains had shorter V3 regions than chronic strains (p = 0.031); subtype C transmission strains had shorter V1 loops than chronic strains (p = 0.047); subtype B transmission strains had more V3 loop glycosites (p = 0.024) than chronic strains. Further investigation showed that these statistically significant results were unlikely to be biologically meaningful. Also, there was no difference observed in the prevalence of a histidine at position 12 among transmission strains and controls of either subtype.
CONCLUSIONS: Although a genetic bottleneck is observed after HIV-1 transmission, our results indicate that summary characteristics of Env hypothesised to be important in transmission are not divergent between transmission and chronic strains of either subtype. The success of a transmission strain to initiate infection may be a random event from the divergent pool of donor viral sequences. The characteristics explored through this study are important, but may not function as genotypic signatures of transmission as previously described.