MAbTope: A Method for Improved Epitope Mapping.

Abs are very efficient drugs, ∼70 of them are already approved for medical use, over 500 are in clinical development, and many more are in preclinical development. One important step in the characterization and protection of a therapeutic Ab is the determination of its cognate epitope. The gold standard is the three-dimensional structure of the Ab/Ag complex by crystallography or nuclear magnetic resonance spectroscopy. However, it remains a tedious task, and its outcome is uncertain. We have developed MAbTope, a docking-based prediction method of the epitope associated with straightforward experimental validation procedures. We show that MAbTope predicts the correct epitope for each of 129 tested examples of Ab/Ag complexes of known structure. We further validated this method through the successful determination, and experimental validation (using human embryonic kidney cells 293), of the epitopes recognized by two therapeutic Abs targeting TNF-α: certolizumab and golimumab.

Natriuretic peptides appeared after their receptors in vertebrates.

BACKGROUND: In mammals, the natriuretic system contains three natriuretic peptides, NPPA, NPPB and NPPC, that bind to three transmembrane receptors, NPR1, NPR2 and NPR3. The natriuretic peptides are known only in vertebrates. In contrast, the receptors have orthologs in all the animal taxa and in plants. However, in non-vertebrates, these receptors do not have natriuretic properties, and most of their ligands are unknown. How was the interaction of the NP receptors and the NP established in vertebrates? Do natriuretic peptides have orthologs in non-vertebrates? If so, what was the function of the interaction? How did that function change? If not, are the NP homologous to ancestral NPR ligands? Or did the receptor's binding pocket completely change during evolution? METHODS: In the present study, we tried to determine if the pairs of natriuretic receptors and their ligands come from an ancestral pair, or if the interaction only appeared in vertebrates. Alignments, modeling, docking, research of positive selection, and motif research were performed in order to answer this question. RESULTS: We discovered that the binding pocket of the natriuretic peptide receptors was completely remodeled in mammals. We found several peptides in non vertebrates that could be related to human natriuretic peptides, but a set of clues, as well as modeling and docking analysis, suggest that the natriuretic peptides undoubtedly appeared later than their receptors during animal evolution. We suggest here that natriuretic peptide receptors in non vertebrates bind to other ligands. CONCLUSIONS: The present study further support that vertebrate natriuretic peptides appeared after their receptors in the tree of life. We suggest the existence of peptides that resemble natriuretic peptides in non-vertebrate species, that might be the result of convergent evolution.

Genomic analysis of codon usage shows influence of mutation pressure, natural selection, and host features on Marburg virus evolution.

BACKGROUND: The Marburg virus (MARV) has a negative-sense single-stranded RNA genome, belongs to the family Filoviridae, and is responsible for several outbreaks of highly fatal hemorrhagic fever. Codon usage patterns of viruses reflect a series of evolutionary changes that enable viruses to shape their survival rates and fitness toward the external environment and, most importantly, their hosts. To understand the evolution of MARV at the codon level, we report a comprehensive analysis of synonymous codon usage patterns in MARV genomes. Multiple codon analysis approaches and statistical methods were performed to determine overall codon usage patterns, biases in codon usage, and influence of various factors, including mutation pressure, natural selection, and its two hosts, Homo sapiens and Rousettus aegyptiacus. RESULTS: Nucleotide composition and relative synonymous codon usage (RSCU) analysis revealed that MARV shows mutation bias and prefers U- and A-ended codons to code amino acids. Effective number of codons analysis indicated that overall codon usage among MARV genomes is slightly biased. The Parity Rule 2 plot analysis showed that GC and AU nucleotides were not used proportionally which accounts for the presence of natural selection. Codon usage patterns of MARV were also found to be influenced by its hosts. This indicates that MARV have evolved codon usage patterns that are specific to both of its hosts. Moreover, selection pressure from R. aegyptiacus on the MARV RSCU patterns was found to be dominant compared with that from H. sapiens. Overall, mutation pressure was found to be the most important and dominant force that shapes codon usage patterns in MARV. CONCLUSIONS: To our knowledge, this is the first detailed codon usage analysis of MARV and extends our understanding of the mechanisms that contribute to codon usage and evolution of MARV.

Accurate determination of epitope for antibodies with unknown 3D structures.

MAbTope is a docking-based method for the determination of epitopes. It has been used to successfully determine the epitopes of antibodies with known 3D structures. However, during the antibody discovery process, this structural information is rarely available. Although we already have evidence that homology models of antibodies could be used instead of their 3D structure, the choice of the template, the methodology for homology modeling and the resulting performance still have to be clarified. Here, we show that MAbTope has the same performance when working with homology models of the antibodies as compared to crystallographic structures. Moreover, we show that even low-quality models can be used. We applied MAbTope to determine the epitope of dupilumab, an anti- interleukin 4 receptor alpha subunit therapeutic antibody of unknown 3D structure, that we validated experimentally. Finally, we show how the MAbTope-determined epitopes for a series of antibodies targeting the same protein can be used to predict competitions, and demonstrate the accuracy with an experimentally validated example.3D: three-dimensionalRMSD: root mean square deviationCDR: complementary-determining regionCPU: central processing unitsVH: heavy chain variable regionVL: light chain variable regionscFv: single-chain variable fragmentsVHH: single-chain antibody variable regionIL4Rα: Interleukin 4 receptor alpha chainSPR: surface plasmon resonancePDB: protein data bankHEK293: Human embryonic kidney 293 cellsEDTA: Ethylenediaminetetraacetic acidFBS: Fetal bovine serumANOVA: Analysis of varianceEGFR: Epidermal growth factor receptorPE: PhycoerythrinAPC: AllophycocyaninFSC: forward scatterSSC: side scatterWT: wild typeKeywords: MAbTope, Epitope Mapping, Molecular docking, Antibody modeling, Antibody-antigen docking.

ß-arrestin signalling and bias in hormone-responsive GPCRs.

G protein-coupled receptors (GPCRs) play crucial roles in the ability of target organs to respond to hormonal cues. GPCRs' activation mechanisms have long been considered as a two-state process connecting the agonist-bound receptor to heterotrimeric G proteins. This view is now challenged as mounting evidence point to GPCRs being connected to large arrays of transduction mechanisms involving heterotrimeric G proteins as well as other players. Amongst the G protein-independent transduction mechanisms, those elicited by ß-arrestins upon their recruitment to the active receptors are by far the best characterized and apply to most GPCRs. These concepts, in conjunction with remarkable advances made in the field of GPCR structural biology and biophysics, have supported the notion of ligand-selective signalling also known as pharmacological bias. Interestingly, recent reports have opened intriguing prospects to the way ß-arrestins control GPCR-mediated signalling in space and time within the cells. In the present paper, we review the existing evidence linking endocrine-related GPCRs to ß-arrestin recruitement, signalling, pathophysiological implications and selective activation by biased ligands and/or receptor modifications. Emerging concepts surrounding ß-arrestin-mediated transduction are discussed in the light of the peculiarities of endocrine systems.

Investigation of novel chemical inhibitors of human lysosomal acid lipase: virtual screening and molecular docking studies.

In the current study, identification of new potent small inhibitors of human lysosomal acid lipase using structure-based methods has been reported. Virtual Screening (VS), compounds from literature and molecular docking studies were employed to find the suitable inhibitors against lysosomal acid lipase (LAL). Specifically for this study a homology model of LipA enzyme was generated based on the structure of dog gastric lipase. As a result of structurebased virtual screening 28 inhibitors were identified from ZINC database. Rest of the inhibitors were selected from literature. Among the studied 65 inhibitors, compound having zinc ID ZINC15707335 exhibiting minimum binding affinity and hydrogen bond and hydrophobic interactions with specific amino acid residues was selected as lead compound.

Comparative genomics analysis of Mycobacterium ulcerans for the identification of putative essential genes and therapeutic candidates.

Mycobacterium ulcerans, the causative agent of Buruli ulcer, is the third most common mycobacterial disease after tuberculosis and leprosy. The present treatment options are limited and emergence of treatment resistant isolates represents a serious concern and a need for better therapeutics. Conventional drug discovery methods are time consuming and labor-intensive. Unfortunately, the slow growing nature of M. ulcerans in experimental conditions is also a barrier for drug discovery and development. In contrast, recent advancements in complete genome sequencing, in combination with cheminformatics and computational biology, represent an attractive alternative approach for the identification of therapeutic candidates worthy of experimental research. A computational, comparative genomics workflow was defined for the identification of novel therapeutic candidates against M. ulcerans, with the aim that a selected target should be essential to the pathogen, and have no homology in the human host. Initially, a total of 424 genes were predicted as essential from the M. ulcerans genome, via homology searching of essential genome content from 20 different bacteria. Metabolic pathway analysis showed that the most essential genes are associated with carbohydrate and amino acid metabolism. Among these, 236 proteins were identified as non-host and essential, and could serve as potential drug and vaccine candidates. Several drug target prioritization parameters including druggability were also calculated. Enzymes from several pathways are discussed as potential drug targets, including those from cell wall synthesis, thiamine biosynthesis, protein biosynthesis, and histidine biosynthesis. It is expected that our data will facilitate selection of M. ulcerans proteins for successful entry into drug design pipelines.

Mycoplasma genitalium: a comparative genomics study of metabolic pathways for the identification of drug and vaccine targets.

Increasing emergence of antibiotic-resistant pathogenic microorganisms is one of the biggest challenges for biomedical research and drug development. Traditional drug discovery methods are time-consuming, expensive and often yield few drug targets. In contrast, advances in complete genome sequencing, bioinformatics and cheminformatics represent an attractive alternative approach to identify drug targets worthy of experimental follow-up. Mycoplasma genitalium is a human parasitic pathogen that is associated with several sexually transmitted diseases. Recently, emergence of treatment-resistant isolates has been reported, which raises serious concern and a need for identification of additional drug targets. In the present study, a systematic workflow consisting of comparative genomics, metabolic pathways analysis and additional drug prioritizing parameters was defined for the identification of novel drug and vaccine targets that are essential for M. genitalium, but absent in its human host. In silico analyses and manual mining identified 79 proteins of M. genitalium, which showed no similarity to human proteins. Among these, 67 proteins were identified as non-homologous essential proteins that could serve as potential drug and vaccine targets. Subcellular localization, molecular weight, and three-dimensional structural characteristics that could facilitate filtering of attractive drug targets were also calculated for the non-homologous essential proteins. Enzymes from thiamine biosynthesis, protein biosynthesis, and folate biosynthesis were identified as attractive candidates for drug development. Furthermore, druggability of each of the identified drug targets was also evaluated by the DrugBank database. Results from this study could facilitate selection of M. genitalium proteins for entry into drug design and vaccine production pipelines.

Computational identification of interplay between phosphorylation and O-ß-glycosylation of human occludin as potential mechanism to impair hepatitis C virus entry.

Hepatitis C virus (HCV) is one of the leading causes of liver diseases. Several host factors that facilitate the attachment and entry of HCV have been discovered, of which human occludin seems to be the most promising. Studies have shown that activity of occludin is dependent upon its phosphorylation status, and that during HCV infection deregulation of phosphorylated occludin collectively leads to a reduction in tight junction (TJ) integrity of hepatocytes and favors HCV entry. However, detailed information of the posttranslational modifications (PTMs) of occludin still remains largely unknown. In addition to phosphorylation, serine/threonine residues of several proteins are also regulated by a unique type of modification known as O-ß-glycosylation and this crosstalk serves as a functional switch. To identify the O-ß-glycosylation potential and how interplay between phosphorylation and O-ß-glycosylation can be exploited for the inhibition of HCV entry, here we report a computational analysis of PTMs of human occludin. Several conserved phosphorylation residues and kinases that can alter the ability of occludin to regulate the integrity of TJs were identified. In addition to previously reported Tyr residues, two additional Tyr residues (Tyr29 and Tyr287) were identified as target sites of Src kinase. To our knowledge, this is the first study to report the O-ß-GlcNAc potential of occludin and target sites of ERK (Ser8, Ser310, and Thr345), GSK-3 (Ser8, Ser341) and Cdk5 (Thr376). Furthermore, based on findings from this study, a potential novel interplay between phosphorylation and O-ß-glycosylation at the two Yin Yang sites (Ser408 and Ser490) is also proposed.

Comparative sequence, antigenic and phylogenetic analysis of avian influenza (H9N2) surface proteins isolated in Pakistan between 1999 and 2008.

INTRODUCTION: Influenza A viruses possess a unique genomic structure which leads to genetic instability, especially in products of neuraminidase and hemagglutinin genes. These surface proteins play major roles in viral entry and release, and in the activation of the host immune system. METHODOLOGY: This study involved an in silico sequence, phylogenetic and antigenic analyses of hemagglutinin and neuraminidase proteins of avian influenza A (H9N2) strains that circulated in Pakistan's poultry flocks from 1999 to 2008 and determined variations among these sequences at different levels. RESULTS: Sequence and phylogenetic analysis revealed a large number of similar substitution mutations and close evolutionary relation among sequences of both proteins. Changes were observed in the N-glycosylation sites of both surface proteins, along with the appearance of a new glycosylation site in the neuraminidase sequence isolated in 2007. Epitopes for hemagglutinin remained conserved, whereas for neuraminidase, epitopes from older strains reappeared in present sequences. CONCLUSIONS: Because of the rapid mutating nature of avian influenza subtype H9N2, constant surveillance of annual sequence variations is important. Preventive measures and vaccine products can be evaluated by keeping track of changes that may lead to reassortment among different circulating strains in Pakistan's commercial poultry flocks or in humans.