Genome-based approaches to develop epitope-driven subunit vaccines against pathogens of infective endocarditis

Infective endocarditis (IE) has emerged as a public health problem due to changes in the etiologic spectrum and due to involvement of resistant bacterial strains with increased virulence. Developing potent vaccine is an important strategy to tackle IE. Complete genome sequences of eight selected pathogens of IE paved the way to design common T-cell driven subunit vaccines. Comparative genomics and subtractive genomic analysis were applied to identify adinosine tri phosphate (ATP)-binding cassette (ABC) transporter ATP-binding protein from Streptococcus mitis (reference organism) as common vaccine target. Reverse vaccinology technique was implemented using computational tools such as ProPred, SYFPEITHI, and Immune epitope database. Twenty-one T-cell epitopes were predicted from ABC transporter ATP-binding protein. Multiple sequence alignment of ABC transporter ATP-binding protein from eight selected IE pathogens was performed to identify six conserved T-cell epitopes. The six selected T-cell epitopes were further evaluated at structure level for HLA-DRB binding through homology modeling and molecular docking analysis using Maestro v9.2. The proposed six T-cell epitopes showed better binding affinity with the selected HLA-DRB alleles. Subsequently, the docking complexes of T-cell epitope and HLA-DRBs were ranked based on XP Gscore. The T-cell epitope (208-LNYITPDVV-216)–HLA-DRB1^?0101 (1T5?W) complex having the best XP Gscore (?13.25?kcal/mol) was assessed for conformational stability and interaction stability through molecular dynamic simulation for 10?ns using Desmond v3.2. The simulation results revealed that the HLA-DRB–epitope complex was stable throughout the simulation time. Thus, the epitope would be ideal candidate for T-cell driven subunit vaccine design against infective endocarditis.     

184 Small molecule identification against novel MDRA protein of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is an obstinate pathogen causing tuberculosis (TB) in Homo sapiens. One third of the World population is affected by Mtb (James et al., 2008). The multidrug-resistant protein-A (MDRA) belongs to ABC transporter family. The protein MDRA and the membrane integral protein MDRB together form the efflux pump (MDRA₂B₂ complex) that confers resistance by transport of the drugs out of the cell. The MDRB protein expression depends on the expression of MDRA (Baisakhee et al., 2002). In the present study, MDRA 3-D model (Figure) was generated with the help of comparative homology modeling techniques using pair-wise sequence alignment. The predicted 3-D model was subjected to refinement and validated. The active site of the protein was predicted. The virtual screening (VS) studies were performed at MDRB binding site with an in-house library of small molecules to identify a lead molecule that can inhibits the MDRA protein. The results of VS project competitive inhibitors of MDRB, for its binding with MDRA, and its drug-resistant activity. Hence, the MDRA protein may be treated as a novel target for the development of new chemical entities for tuberculosis therapy (Bhargavi et al., 2010; Malkhed et al., 2011).     

187 Plucking the high hanging fruit: a systematic approach for targeting protein interfaces

Development of specific ligands for protein targets that help decode the complexities of protein–protein interaction networks is a key goal for the field of chemical biology. Despite the emergence of powerful in silico and experimental high-throughput screening strategies, the discovery of synthetic ligands that selectively modulate protein–protein interactions remains a challenge for the chemical biologists. Proteins often utilize small folded domains for recognition of other biomolecules. The basic hypothesis guiding our research is that by mimicking these domains, we can modulate the function of a particular protein with metabolically-stable synthetic molecules (Raj et al., 2013). This presentation will discuss computational approaches (Bullock et al., 2011; Jochim & Arora, 2010) to identify targetable interfaces along with synthetic methods (Patgiri et al., 2008; Tosovska & Arora, 2010) to develop protein domain mimics (PDMs) as modulators of intracellular protein–protein interactions (Henchey et al., 2010; Patgiri et al., 2011).     

34 Cyo-EM visualization of Mycobacterium 70S ribosome reveals unique structural components at the function sites

The 3D structures of prokaryotic and eukaryotic ribosomes by crystallography and electron microscopy have revealed that they share an evolutionarily conserved core (Schmeing & Ramakrishnan, 2009), but each of the ribosomes contains its own set of specific proteins (or extensions of conserved proteins) and expansion segments of rRNAs (Melnikov et al., 2012). How these differences correlate to function still remains largely unknown. A 3D cryo-EM map of the 70S ribosome from Mycobacterium smegmatis (Msm70S) unveiled striking new structural features (Shasmal & Sengupta, 2012). The core of the Msm70S shows overall similarity with the core of the Escherichia coli 70S ribosome while containing additional mass in the periphery and solvent exposed sides. Some of the Mycobacterium ribosomal proteins are significantly bigger as compared to the E. coli counter parts. The rRNAs also contain extra helices, also revealed by their secondary structures. Most of the additional density of the Msm70S can be largely attributed to the extra helices present in the rRNAs, and extra domains of homologous proteins. One of the most notable features appears in the large subunit near L1 stalk as a structure forming a long helix with its upper end located in the vicinity of the mRNA exit channel (which we term the ‘steeple’). We propose that the prominent helical structure in mycobacterium 23S rRNA participates in modulating different steps of translation, especially the E site tRNA exit mechanism and propagation of mRNA 5′ end.     

160 A dynamic model for the linker histone

Linker histones play an important role in the packing of chromatin. This family of proteins generally consists of a short, unstructured N-terminal domain, a central globular domain, and a C-terminal domain (CTD). The CTD, which makes up roughly half of the protein, is intrinsically disordered in solution but adopts a specific fold upon interaction with DNA (Fang et al., 2012). While the globular domain structure is well characterized, the structure of the CTD remains unknown. Sequence alignment alone does not reveal any significant homologs for this region of the protein. Construction of a model thus requires additional information. For example, the atomic model for the rat histone H1d CTD, proposed over a decade ago, used novel bioinformatics tools and biochemical data (Bharath et al., 2002). New fluorescence resonance energy transfer (FRET) studies of the folding of the CTD in the presence of linear DNA, single nucleosomes, and oligonucleosomal arrays (Caterino et al., 2011; Fang et al., 2012) have stimulated our interest in constructing a dynamic model of the protein. We have obtained preliminary information about the structure and dynamics of the linker histone CTD through ab initio folding simulations using the Rosetta modeling package (Rohl et al., 2004). By analyzing a large number of conformations sampled through a Monte Carlo procedure, we get a clearer picture of the preferred states of the protein and its dynamics. Our results show that the CTD may frequently adopt a structure with 3–5 helices and helix-turn-helix motifs in specific regions. Some of the best scoring structures show high similarity with the HMG-box-containing proteins previously used as templates by Bharath et al. Further clustering analysis of our results hints of a preferred set of conformations for the CTD of the linker histone. Comparison of these models with distances measured by FRET may help account for the distinct structures of the CTD observed upon binding to different macromolecular partners.     

20 Supramolecular polymerization of nucleobase-like monomers in water

Elucidating the physiochemical principles that govern molecular self-assembly is of great importance for understanding biological systems and may provide insight into the emergence of the earliest macromolecules of life, an important challenge facing the RNA World hypothesis. Self-assembly results from a delicate balance between multiple noncovalent interactions and solvent effects, but achieving efficient self-assembly in aqueous solution with synthetic molecules has proven particularly challenging. Here, we demonstrate how two physical properties – monomer solubility and large hydrophobic surfaces of intermediate structures – are key elements to achieving supramolecular polymers in aqueous solution (Cafferty et al., 2013). Applying these two principles, we report the highly cooperative self-assembly of two weakly interacting, low molecular weight monomers [cyanuric acid and a modified triaminopyrimidine] into a water-soluble supramolecular assembly (see scheme below). The observed equilibrium between only two appreciably populated states – free monomers and supramolecular assemblies – is in excellent agreement with the values previously determined for the free energy of hydrogen bonding (Klostermeier & Millar, 2002), π???π stacking (Frier et al., 1985), and the calculated free energy penalty for the solvation of hydrophobic structures in water (Chandler, 2005). The similarity of the molecules used in this study for the nucleobases found in contemporary nucleic acids and the demonstration that these monomers assemble while the natural nucleobases do not, suggests that the first informational polymers may have emerged from a similar self-assembly process, if the nucleobases were different then they are today (Hud et al., 2013).     

161 Discovery of potent KdsA inhibitors of Leptospira interrogans through homology modeling,docking, and molecular dynamics simulations

Leptospira interrogans is the foremost cause of human leptospirosis. Discovery of novel lead molecules for common drug targets of more than 250 Leptospira serovars is of significant research interest. Lipopolysaccharide (LPS) layer prevent entry of hydrophobic agents into the cell and protect structural integrity of the bacterium. KDO-8-phosphate synthase (KdsA) catalyzes the first step of KDO biosynthesis that leads to formation of inner core of LPS. KdsA was identified as a potential drug target against Leptospira interrogans through subtractive genomic approach, metabolic pathway analysis, and comparative analysis (Amineni et al., 2010). The present study rationalizes a systematic implementation of homology modeling, docking, and molecular dynamics simulations to discover potent KdsA inhibitors (Pradhan et al., 2013; Umamaheswari et al., 2010). A reliable tertiary structure of KdsA in complex with substrate PEP was constructed based on co-crystal structure of Aquifex aeolicus KdsA synthase with PEP using Modeller9v10. Geometry-based analog search for PEP was performed from LigandInfo database to generate an in house library of 352 ligands. The ligand data-set was docked into KdsA active site through three-stage docking technique (HTVS, SP, and XP) using Glidev5.7. Thirteen lead molecules were found to have better binding affinity compared to PEP (XP Gscore?=??7.38?kcal/mol; Figure 1). The best lead molecule (KdsA- lead1 docking complex) showed XP Gscore of ?10.26?kcal/mol and the binding interactions (Figure 2) were correlated favorably with PEP–KdsA interactions (Figure 1). Molecular dynamics simulations of KdsA– lead1 docking complex for 10?ns had revealed that the complex (Figure 3) remained stable in closer to physiological environmental condition. The predicted pharmacological properties of lead1 were well within the range of a drug molecule with good ADME profile, hence, would be intriguing towards development of potent inhibitor molecule against KdsA of Leptospira.     

Vaccinomics approach for developing multi-epitope peptide pneumococcal vaccine

Streptococcus pneumoniae is a leading cause of some diseases such as pneumonia, sepsis, and meningitis mostly in children less than 5?years of age. Presently, two types of pneumococcal vaccine are available on the market: polysaccharide vaccines (PPV) that are based on capsular polysaccharides of at least 92 different serotypes, and protein-conjugated polysaccharide vaccine (PCV). The PPVs such as PPV23 do not stimulate efficient protective immunity in children under 2?years old, while the PCVs such as PCV7, PCV10, and PCV13 that cover 7, 10, and 13 serotypes, respectively, highly protect newborns, but have some disadvantages such as complications in manufacturing, costly production, and also requires refrigeration and multiple injections. Epitope-based vaccines, including varied mixtures of conserved virulence proteins, are a promising alternative to the existing capsular antigen vaccines. In this study, it has been tried to design an efficient subunit vaccine in order to elicit both CTL and HTL responses. The immunodominant epitopes from highly protective antigens of S. pneumoniae (PspA, CbpA, PiuA, and PhtD) were selected from different databanks, such as IEDB, PROPRED, RANKPEP, and MHCPRED. The PspA and CbpA were chosen as CTL epitope stimulants, and PhtD and PiuA were defined as helper epitopes. Because of low immunogenicity of epitope vaccines, PorB protein as a TLR2 agonist was employed to increase the immunogenicity of the vaccine. All the peptide segments were fused to each other by proper linkers, and the physicochemical, structural, and immunological characteristics of the construct were also evaluated. To achieve a high-quality 3?D structure of the protein, modeling, refinement, and validation of the final construct were done. Docking and molecular dynamics analyses demonstrated an appropriate and stable interaction between the vaccine and TLR2 during the simulation period. The computational studies suggested the designed vaccine as a novel construct, capable to elicit efficient humoral and cellular immunities, which are crucial for protection against S. pneumoniae.

Communicated by Ramaswamy H. Sarma     

203 Polymers through pores: single-molecule experiments with nucleic acids,polypeptides and polysaccharides

The translocation of polymers through pores has been examined for almost two decades with an emphasis on nucleic acids. There are also interesting circumstances in biology where polypeptides and polysaccharides pass through transmembrane pores, and our laboratory has been investigating examples of them. Single-molecule nucleic acid sequencing by nanopore technology is an emerging approach for ultrarapid genomics. Strand sequencing with engineered protein nanopores is a viable technology which has required advances in four areas: nucleic acid threading, nucleobase identification, controlled strand translocation, and nanopore arrays (Bayley, 2012). The latter remain a pressing need and our attempts to improve arrays will be described. In several physiological situations, folded proteins pass through transmembrane pores. We have developed a model system comprising mutant thioredoxins as the translocated proteins, and staphylococcal alpha-hemolysin, as the pore. Our findings support a mechanism in which there is local unfolding near the terminus of the polypeptide that enters the pore. The remainder of the protein then unfolds spontaneously and diffuses through the pore into the recipient compartment (Rodriguez-Larrea & Bayley, 2013). We have also examined the pore formed by the E. coli outer membrane protein Wza, which transports capsular polysaccharide from its site of synthesis to the outside of the cell. We made mutant open forms of the pore and screened blockers for them by electrical recording in planar bilayers. The most effective blocker binds in the alpha-helix barrel of Wza, a site accessible from the external medium, and therefore, a prospective target for antibiotics (Kong et al., 2013).     

172 Role of conserved water molecular triad in the recognition of IMP,NAD+ with Asp 274, Asn 303, Arg 322, and Asp 364 in both the isoform of hIMPDH

Inosine monophosphate dehydrogenase (IMPDH) plays an important role in the Guanosine monophosphate (GMP) biosynthesis pathway. As hIMPDH-II is involved in CML-Cancer, it is thought to be an active target for leukemic drug design. The importance of conserved water molecules in the salt-bridge-mediated interdomain recognition and loop-flap recognition of hIMPDH has already been indicated in some simulation studies (Bairagya et al., 2009, 2011a, 2011b, 2012; Mishra et al., 2012). In this work, the role of conserved water molecules in the recognition of Inosine monophosphate (IMP) and NAD⁺ (co-factor) to active site residues of both the isoforms has been investigated by all atoms MD-Simulation studies. During 25-ns dynamics of the solvated hIMPDH-II and I (1B3O and 1JCN PDB structures), the involvement of conserved water molecular triad (W _M, W _L and W _C) in the recognition of active site residues (Asp 274, Asn 303, Arg 322, and Asp 364), IMP and NAD⁺ has been observed (Figure 1). The H-bonding co-ordination of all three conserved water molecular centers is within 4–7 and their occupation frequency is 1.0. The H-bonding geometry and the electronic consequences of the water molecular interaction at the different residues (and also IMP and NAD⁺) may put forward some rational clues on antileukemic agent design.     

210 Computational study of substrates and mediators features of lacasses

Laccases are enzymes of the family multicopper oxidases, being widely used for biotechnological applications (Canas & Camarero, 2010). The enzymes’ catalytic cycle consists of the oxidation of the substrate with the concomitant reduction of molecular oxygen to water. In this process, the substrate is converted to a free radical, that can oxidize larger substrates acting as a mediator or it can undergo polymerization. Substrate binding is not specific, and there is a large diversity of substrates for laccases. Moreover, the binding site shows important differences among diverse species. The goal of the present work is to characterize the laccase binding pocket of different species, in order to establish their common pharmacophoric characteristics. For this purpose, we have carried out docking studies with a subset of substrates, covering the diversity of substrates using the Glide program (Friesner et al., 2004). We have also analyzed the characteristics of the binding site using diverse probes. We further have rationalized the differential values of km found among diverse species for a specific substrate. Finally, special attention has been devoted to the binding of the mediator 2,2′-azido-di-(3-ethylbenzothiazoline)-6-sulfonic acid (ABTS), commonly used in industrial processes. Figure 1 shows, ABTS docked onto the fungal laccase, whereas Figure 2 shows ABTS docked onto the bacterial laccase. The analysis of the protein–ligand complex together with the corresponding optimized geometries of the possible substrate species carried out using DFT suggest that the bound species is the protonated form of ABTS as previously suggested (Enguita et al., 2004). Furthermore, the results of this study also suggest that its mechanism of oxidation occurs in a similar way to the rest of substrates/mediators, in contrast to previous suggestions (Fabbrini et al., 2002).      

Prediction of Potential Epitopes for Peptide Vaccine Formulation Against Teschovirus A Using Immunoinformatics

Teschovirus A belongs to the family Picornaviridae and is a causal agent of the disease Teschovirus encephalomyelitis and other infections that remain asymptomatic. The present study was performed to design epitope-based peptide vaccine against Teschovirus A by identifying the potential T cell and B-cell epitopes from capsid proteins (VP1, VP3 and VP2) of the virus using reverse vaccinology and immunoinformatics approaches. In the current study, hexapeptide T-cell and octapeptide B-cell epitopes were analyzed for immunogenicity, antigenicity and hydrophilicity scores of each epitope. Each potential epitope was further characterized using ExPASy-ProtParam and Antimicrobial Peptide Database (APD3) tools for determining various physical and chemical parameters of the epitope. One linear hexapeptide T-cell epitope, i.e., RPVNDE (epitope position 77–82) and one linear octapeptide B-cell epitope, i.e., AYSRSHPQ (236–243) were identified from the viral capsid protein as they possess the capability to raise effective immunogenic reaction in the host organism against the virus. Pharmaceutical industries could harness the results of this investigation to develop epitope-based peptide vaccines by loading the identified epitopes in combination with targeting signal peptides of T-cells and B-cells and then inserting the combination into virus like particle (vlp) or constructing subunit vaccines for further trial.

     

24 Evolutionary dynamics of inteins in bacteria

Inteins are protein sequences that autocatalytically splice themselves out of protein precursors – analogous to introns – and ligate the flanking regions into a functional protein. Inteins are present in all three kingdoms of life, but have a sporadic distribution. They are found predominantly in proteins involved in DNA replication and repair such as helicases. The distribution of inteins suggests an adaptive function. The evolutionary forces which shaped the observed distribution of inteins are generally unknown. Some authors view inteins only as the selfish elements and argue that frequent horizontal transfer is behind inteins sporadic dissemination (Gogarten et al., 2002). On the other hand, the ancient nature of the inteins and the process of gain/loss could lead to the scattered distribution of inteins among species (Pietrokovski, 2001). It is necessary to note that the exclusively selfish nature of inteins is questionable; recent findings support the hypothesis of possible functional roles of inteins in protein regulation (Callahan et al., 2011). Moreover, both hypotheses were built on a limited number of the intein representatives. The amount of genomic data available for bacteria is enormous and in silico analysis for diverse inteins is warranted. We decided to take advantage of these microbial genomic data and performed comprehensive mining for the inteins using a bioinformatic pipeline. Altogether, 1757 species were analysed from 19 major phyla yielding more than 4500 intein-like sequences. The majority of these bacterial inteins were not described previously. Approximately 55% of the inteins were found in polymerases, helicases, or recombinases (Figure 1). Phylogenetic analysis indicated the complex evolutionary dynamics of inteins which includes horizontal transfers, high evolutionary rates coupled with recurrent gains, and losses. The preponderance of inteins in helicases and reductases is being investigated in terms of functional relevance.     

3 Origins and evolution of the translation machinery

The protein synthesis machinery largely evolved prior to the last common ancestor and hence its study can provide insight to early events in the origin of life, including the transition from the hypothetical RNA world to living systems as we know them. By utilizing information from primary sequences, atomic resolution structures, and functional properties of the various components, it is possible to identify timing relationships (Hsiao et al., 2009; Fox, 2010). Taken together, these timing events are used to develop a preliminary time line for major evolutionary events leading to the modern protein synthesis machinery. It has been argued that a key initial event was the hybridization of two or more RNAs that created the peptidyl transferase center, (PTC), of the ribosome (Agmon et al. 2005). The PTC, left side of figure, contains a characteristic cavity/pore that serves as the entrance to the exit tunnel and is thought to be essential to the catalysis (Fox et al., 2012). This cavity is distinct from typical RNA pores (right side of figure) in that the nitrogenous bases face towards the lumen of the pore and thus are available for hydrogen bonding interactions. In typical RNA pores, the bases carefully avoid the lumen region. In support of Agmon et al. 2005), it is argued that this key difference reflects the fact the pore was created by an early hybridization event rather than normal RNA folding.