169 Polymorphic assembly of computationally designed hydrophobic-rich collagen peptides

We seek to understand how the position and length of hydrophobic content within a collagen peptide sequence dictates morphology of self-assembly. We modeled collagen assembly using diffusion limited aggregation¹ (DLA) (Parkinson et al. 1995). of discretized, rigid rods composed of hydrophilic and hydrophobic spheres. Simulations predicted that the inclusion of short hydrophobic domains should direct the assembly of lamellar structures. We designed a set of collagen peptide sequences with six, five and four contiguous nonpolar residues. Electron microscopy of aggregates revealed the peptide with six nonpolar residues self-assembled into uniform fibrils and the peptide with five residues assembled into both fibrils and plates, while including four hydrophobic residues that formed only plates. This polymorphic behavior can be explained by packing models of rod vs. screw-like-particles.     

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).     

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).     

177 T-cell vaccine design for Streptococcus pneumoniae: an in silico approach

Streptococcus pneumoniae (pneumococcus) remains an important cause of meningitis, bacteremia, acute otitis media, community acquired pneumonia associated with significant morbidity, and mortality world wide. Conjugated polysaccharide, glycoconjugated, and capsular polysaccharide based vaccines were existent for pneumococcal disease but are still specific and restricted to serotypes of S. pneumoniae. Proteome of eight serotypes of S. pneumoniae was retrieved and identified in common proteins (Munikumar et al., 2012). 18 membrane proteins were distinguished from 1657 common proteins of eight serotypes of S. pneumoniae. Implementing comparative genomic approach and subtractive genomic approach, three membrane proteins were predicted as essential for bacterial survival and non-homologous to human (Munikumar et al., 2012; Umamaheswari et al., 2011). ProPred server was used to propose four promiscuous T-cell epitopes from three membrane proteins and validated through published positive control, SYFPEITHI and immune epitope database (Munikumar et al., in press). The four epitopes docked into peptide binding region of predominant HLA-DRB alleles with good binding affinity in Maestro v9.2. The T-cell epitope 89-VVYLLPILI-97 and HLA-DRB5^?0101 docking complex was with best XPG score (?13.143?kcal/mol). Further, the stability of the complex was checked through molecular dynamics simulations in Desmond v3.3. The simulation results had revealed that the complex was stable throughout 5000?ps (Munikumar et al., in press). Thus, the epitope would be the ideal candidate for T-cell driven subunit vaccine design against selected serotypes of S. pneumoniae.     

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.     

136 Dendritic Alkyl Chains on DNA Cages: A Geometry-Dependent Inter- or Intramolecular “Handshake”

The selective association of hydrophobic sidechains is a strong determinant of protein organization. We have observed a parallel mode of assembly in DNA nanotechnology. Firstly, dendritic DNA amphiphiles (D-DNA) were synthesized (Carneiro, Aldaye, & Sleiman, 2009) comprising an addressable oligonucleotide portion and a hydrophobic alkyl dendron at the 5’ terminus. DNA amphiphiles have gathered interest recently as they can self-assemble in aqueous media to form well defined micelles while also retaining the ability to hybridize to their complement (Kwak & Herrmann, 2011; Patwa, et al. 2011) Two variations of alkyl D-DNA were hybridized to the single-stranded edges of a DNA cube (McLaughlin, et al., 2012). It was found that anisotropic organization of these hydrophobic domains on the 3D scaffold results in a new set of assembly rules, dependent on spatial orientation, number, and chemical identity of the D-DNA on the cubic structure (Edwardson. et al. 2012). When four amphiphiles are organized on one cube face, the hydrophobic residues engage in an intermolecular “handshake” between two cubes, resulting in a dimer. When eight amphiphiles are organized on the top and bottom faces of the cube, they engage in a “handshake” inside the cube. Combining the highly specific recognition of the oligonucleotide sequence with the orthogonal association of hydrophobic moieties can lead to a variety of structures with such diverse applications as membrane anchoring, cell uptake, directed hydrophobic assembly, and encapsulation and release of small molecules.     

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.     

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).      

The In Vitro Self-Assembly of Supramolecular Structures after the Spontaneous Disassembly of Carotovoricins

The self-assembly of supramolecular structures (empty sheaths and polysheaths of the macromolecular Erwinia carotovora bacteriocins) was studied by electron microscopy in the course of 1- to 2-year incubation of phage particles at 4°C. This study showed that the empty sheaths and polysheaths of the bacteriocins of eight E. carotovora strains spontaneously assemble at the self-assembly centers (or crystallization centers), which have a diameter of 26–65 nm and contain a dense proteinaceous material. The self-assembly center consists of two components, a primer and the structural protein of contracted sheaths. Empty sheaths assembled in the crystallization centers are polar structures synthesized through the stepwise head-to-tail polymerization of monomeric units. The supramolecular structures of two E. carotovora 62A bacteriocins are assembled in a different way. At the early stages of their self-assembly, a reticular structure is formed, which then transforms into very long polysheaths composed of monomers. Along with polysheaths, rounded or lamplike structures 33–117 nm in size composed of the subunits of contracted sheaths are produced. Carotovoricins may serve as suitable objects for the study of the self-assembly of elementary biological structures.     

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.     

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).     

150 Getting a grip on alpha-synuclein amyloid oligomers

Misfolding and aggregation of proteins into nanometer-scale fibrillar assemblies is a hallmark of many neurodegenerative diseases. Aggregation of the human alpha-synuclein protein is implicated in the etiology of Parkinson’s disease. A particularly relevant question is the role of early oligomeric aggregates of alpha-synuclein in modulating the dynamics of protein aggregation, and in the interactions with essential cellular components. However, very little is known about the molecular details of these aggregate species. For large protein aggregates, such as alpha-synuclein oligomers, it is very difficult to determine the number of monomers that form an oligomer using conventional techniques. We have developed a method that uses sub-stoichiometric labeling, that is, only a fraction of the monomers contains a fluorescent label, in combination with single-molecule photobleaching to determine the number of monomers per oligomer (Zijlstra et al., 2012). The number of bleaching steps gives the number of fluorescent labels per oligomer. Knowing the exact label density, that is, the fraction of labeled monomers at the start of the aggregation, we can correlate the number of fluorescent labels per oligomer to the total number of monomers. Using this method, we can determine the composition, probe the distribution in the number of monomers per oligomer, and investigate the influence of the fluorescent label on the aggregation process. For wild-type alpha-synuclein, we find no distribution in the number of monomers per oligomer and find a single, well-defined oligomeric species consisting of ～30 monomers per oligomer. On the other hand, for oligomers formed in the presence of dopamine, we find a distinctly bimodal distribution suggesting the existence of two populations of oligomeric species.     

39 Simulation of RNA tandem GA base pairs provides insights about the forces behind conformational preference

Conformational changes are important for RNA function. We used molecular mechanics with all-atom models to understand conformational preference in RNA tandem guanine–adenine (GA) base pairs. These tandem GA base pairs play important roles in determining the stability and structural dynamics of RNA tertiary structures. Previous solution structures showed that these tandem GA base pairs adopt either imino (cis-Watson-Crick/cis-Watson-Crick interaction) or sheared (trans-Hoogsteen/trans-Hoogsteen interaction) pairing depending upon the sequence and orientation of the adjacent base pairs. In our simulations, we modeled (GCGGACGC)₂ (Wu and Turner 1996) and (GCGGAUGC)₂ (Tolbert et al., 2007), experimentally preferred as imino and sheared, respectively. Besides the experimentally preferred conformation, we constructed models of the nonnative conformations by changing cytosine to uracil or uracil to cytosine. We used explicit solvent molecular dynamics and free energy calculation with umbrella sampling to measure the free energy deference of the experimentally preferred conformation and the nonnative conformations. A modification to ff10 was required, which allowed the guanine bases’ amino group to leave the base plane (Yildirim et al., 2009). With this modification, the RMSD of unrestrained simulations and the free energy surfaces are improved, suggesting the importance of electrostatic interactions by G amino groups in stabilizing the native structures.     

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.     

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.     

43 Why G? Aggregation and gelation of GMP with XMP: the plot thickens

GMP alone, among the individual ribonucleotides, exhibits a reversible self-aggregation through hydrogen bonding to form tetrads that are the building blocks of higher order structures. These “G-tetrads” can further associate through π–π stacking to form chiral, columnar aggregates and, at higher monomer concentrations, lyotropic liquid crystalline phases. This alternate pathway for GMP should compete with its incorporation into oligonucleotides, which is why it is difficult to synthesize or amplify highly G-rich RNA or DNA with good efficiency in the absence of natural proteins, such as helicases, that function to unwind the strands. Given this competing pathway for GMP, we can ask if it came to be one of the four ribonucleotides in modern RNA in spite of, or because of, its unique properties. Our hypothesis is that the competition between reversible aggregation and covalent polymerization directed RNA toward sequences that were best suited to life on early earth. We find support in the observation that the same interactions that promote self-assembly of monomeric GMP also promote folding of G-rich RNA and DNA sequences to form inter- and intramolecular G-quadruplex structures. Such sequences are prevalent throughout the biological world and are thought to serve important functions related to genomic stability and gene regulation. G-quadruplex structures are also common motifs in aptamers, which are combinatorially derived DNA or RNA sequences that exhibit highly selective, high-affinity binding to molecular and macromolecular targets. An important consideration for GMP aggregation in a prebiotic RNA World scenario is the effect of other XMP on GMP self-assembly. In this talk, we will focus on the properties of solutions containing mixtures of GMP with AMP, CMP, and UMP. The results show that each nucleotide exerts a different influence on the self-assembly of GMP, raising interesting questions about scenarios on prebiotic Earth that would be consistent with abiotic RNA polymerization.     

128 Some derivatives of 1,3-diazaadamantane strongly stimulate strand exchange reaction between short oligonucleotides

Earlier, we characterized the behavior of HG122 (compound III) in the DNA Strand Exchange Reaction (SER) and found that 1,3-diazaadamantane derivative facilitates SER in the system of short oligonucleotides (Gabrielian et al., 2011). In the present study, a series of new derivatives of 1,3-diazaadamantane have been synthesized with the purpose to discern how small variations in the compound structure can influence its activity in SER and try to get more effective substances for stimulation of SER. We revealed that most of the small variations in the structure significantly influenced the compounds’ efficacy in accelerating SER. For example, an increase in the compounds’ aliphatic chain lengths considerably enhanced its efficiency in SER stimulation and in the series of compounds presented in the Figure HG188 (compound IV) was eminently the most potent agent in the stimulation of SER. Small modifications in other parts of the 1,3-diazaadamantane molecule also influenced the SER stimulation and several derivatives more efficient in facilitating SER than HG122 were revealed. Some of the compounds exhibited virtually negligible capability to stimulate SER but, interestingly, out of 12 derivatives characterized, agents that retard spontaneous SER were not found. Earlier, the stimulation of the DNA strand exchange was documented for different ligands of the policationic nature such as Cationic Comb Copolymers (Kim et al., 2003), linker histones (Bocharova et al., 2012), cobalt hexamine (unpublished observation) etc. The present results provide us with a novel class of SER facilitating compounds – the cationic amphiphiles that can serve as interesting objects for further understanding of different aspects of DNA SER. Biomedical implications and prospects in biomedical applications of these and similar compounds’ activity in SER remain to be investigated.