In silico identification of drug targets and vaccine candidates against Bartonella quintana: a subtractive proteomics approach.

BACKGROUND: The availability of genes and protein sequences for parasites has provided valuable information for drug target identification and vaccine development. One such parasite is Bartonella quintana, a Gram-negative, intracellular pathogen that causes bartonellosis in mammalian hosts. OBJECTIVE: Despite progress in understanding its pathogenesis, limited knowledge exists about the virulence factors and regulatory mechanisms specific to B. quintana. METHODS AND FINDINGS: To explore these aspects, we have adopted a subtractive proteomics approach to analyse the proteome of B. quintana. By subtractive proteins between the host and parasite proteome, a set of proteins that are likely unique to the parasite but absent in the host were identified. This analysis revealed that out of the 1197 protein sequences of the parasite, 660 proteins are non-homologous to the human host. Further analysis using the Database of Essential Genes predicted 159 essential proteins, with 28 of these being unique to the pathogen and predicted as potential putative targets. Subcellular localisation of the predicted targets revealed 13 cytoplasmic, eight membranes, one periplasmic, and multiple location proteins. The three-dimensional structure and B cell epitopes of the six membrane antigenic protein were predicted. Four B cell epitopes in KdtA and mraY proteins, three in lpxB and BQ09550, whereas the ftsl and yidC proteins were located with eleven and six B cell epitopes, respectively. MAINS CONCLUSIONS: This insight prioritises such proteins as novel putative targets for further investigations on their potential as drug and vaccine candidates.

Structural Characteristics of Glycocins: Unraveling the Role of S-Linked Carbohydrates and Other Structural Elements.

Glycocins are antimicrobial peptides with glycosylations, often an S-linked monosaccharide. Their recent structure elucidation has brought forth questions about their mechanisms of action as well as the impact of S-glycosylation on their structural behavior. Here, we investigated structural characteristics of glycocins using a computational approach. Depending on the peptide's class (sublancin- or glycocin F-like), the sugar changes the peptide's flexibility. Also, the presence of glycosylation is necessary for the lack of structure of Asm1. The C-terminal tail in glycocin F-like peptides influenced their structured regions, acting like a regulator. These findings corroborate the versatility of these post-translational modifications, pointing toward their potential use in molecular engineering.

ConfID: an analytical method for conformational characterization of small molecules using molecular dynamics trajectories.

MOTIVATION: The conformational space of small molecules can be vast and difficult to assess. Molecular dynamics (MD) simulations of free ligands in solution have been applied to predict conformational populations, but their characterization is often based on clustering algorithms or manual efforts. RESULTS: Here, we introduce ConfID, an analytical tool for conformational characterization of small molecules using MD trajectories. The evolution of conformational sampling and population frequencies throughout trajectories is calculated to check for sampling convergence while allowing to map relevant conformational transitions. The tool is designed to track conformational transition events and calculate time-dependent properties for each conformational population detected. AVAILABILITY AND IMPLEMENTATION: Toolkit and documentation are freely available at http://sbcb.inf.ufrgs.br/confid. CONTACT: marcelo.poleto@ufv.br or bigrisci@inf.ufrgs.br. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A cell surface arabinogalactan-peptide influences root hair cell fate.

Root hairs (RHs) develop from specialized epidermal trichoblast cells, whereas epidermal cells that lack RHs are known as atrichoblasts. The mechanism controlling RH cell fate is only partially understood. RH cell fate is regulated by a transcription factor complex that promotes the expression of the homeodomain protein GLABRA 2 (GL2), which blocks RH development by inhibiting ROOT HAIR DEFECTIVE 6 (RHD6). Suppression of GL2 expression activates RHD6, a series of downstream TFs including ROOT HAIR DEFECTIVE 6 LIKE-4 (RSL4) and their target genes, and causes epidermal cells to develop into RHs. Brassinosteroids (BRs) influence RH cell fate. In the absence of BRs, phosphorylated BIN2 (a Type-II GSK3-like kinase) inhibits a protein complex that regulates GL2 expression. Perturbation of the arabinogalactan peptide (AGP21) in Arabidopsis thaliana triggers aberrant RH development, similar to that observed in plants with defective BR signaling. We reveal that an O-glycosylated AGP21 peptide, which is positively regulated by BZR1, a transcription factor activated by BR signaling, affects RH cell fate by altering GL2 expression in a BIN2-dependent manner. Changes in cell surface AGP disrupts BR responses and inhibits the downstream effect of BIN2 on the RH repressor GL2 in root epidermis.

The Lazy Life of Lipid-Linked Oligosaccharides in All Life Domains.

Lipid-linked oligosaccharides (LLOs) play an important role in the N-glycosylation pathway as the donor substrate of oligosaccharyltransferases (OSTs), which are responsible for the en bloc transfer of glycan chains onto a nascent polypeptide. The lipid component of LLO in both eukarya and archaea consists of a dolichol, and an undecaprenol in prokarya, whereas the number of isoprene units may change between species. Given the potential relevance of LLOs and their related enzymes to diverse biotechnological applications, obtaining reliable LLO models from distinct domains of life could support further studies on complex formation and their processing by OSTs, as well as protein engineering on such systems. In this work, molecular modeling techniques, such as quantum mechanics calculations, molecular dynamics simulations, and metadynamics were employed to study eukaryotic (Glc3-Man9-GlcNAc2-PP-Dolichol), bacterial (Glc1-GalNAc5-Bac1-PP-Undecaprenol), and archaeal (Glc1-Man1-Gal1-Man1-Glc1-Gal1-Glc1-P-Dolichol) LLOs in membrane bilayers. Microsecond molecular dynamics simulations and metadynamics calculations of LLOs revealed that glycan chains are more prone to interact with the membrane lipid head groups, while the PP linkages are positioned at the lipid phosphate head groups level. The dynamics of isoprenoid chains embedded within the bilayer are described, and membrane dynamics and related properties are also investigated. Overall, there are similarities regarding the structure and dynamics of the eukaryotic, the bacterial, and the archaeal LLOs in bilayers, which can support the comprehension of their association with OSTs. These data may support future studies on the transferring mechanism of the oligosaccharide chain to an acceptor protein.

Current Status of Carbohydrates Information in the Protein Data Bank.

Carbohydrates are well known for their physicochemical, biological, functional, and therapeutic characteristics. Unfortunately, their chemical nature imposes severe challenges for the structural elucidation of these phenomena, impairing not only the depth of our understanding of carbohydrates but also the development of new biotechnological and therapeutic applications based on these molecules. In the recent past, the amount of structural information, obtained mainly from X-ray crystallography, has increased progressively, as well as its quality. In this context, the current work presents a global analysis of the carbohydrate information available in the Protein Data Bank (PDB). From high quality structures, it is clear that most of the data are highly concentrated on a few sets of residue types, on their monosaccharidic forms, and connected by a small diversity of glycosidic linkages. The geometries of these linkages can be mostly associated with the types of linkages instead of residues, while the level of puckering distortion was characterized, quantified, and located in a pseudorotational equilibrium landscape, not only to local minima but also to transitional states. These qualitative and quantitative analyses offer a global picture of the carbohydrate structural content in the PDB, potentially supporting the building of new models for carbohydrate-related biological phenomena at the atomistic level, including new developments on force field parameters.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations.

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

Design, synthesis and pharmacological assessment of new pyrazole compounds.

AIMS: This study investigated the antinociceptive and anti-inflammatory effects of new pyrazole compounds LQFM011(5), LQFM043(6) and LQFM044(7) as well as the mechanisms of action and acute in vitro toxicity. MAIN METHODS: The antinociceptive activity was evaluated using the acetic acid-induced abdominal writhing test, formalin-induced pain test and the Randall-Selitto test. The anti-inflammatory activity was evaluated using models of paw oedema and pleurisy induced by carrageenan; cell migration, the levels of tumour necrosis factor α (TNF-α) and myeloperoxidase (MPO) enzyme activity were evaluated. In addition, the ability to inhibit phospholipase A2 (PLA2) in vitro and docking in PLA2 were used. Acute oral systemic toxicity in mice was evaluated through the neutral red uptake assay. KEY FINDINGS: The synthesised compounds (5-7), delivered via gavage (p.o.) at 70, 140 or 280 µmol/kg, decreased the number of writhings induced by acetic acid; the three compounds (280 µmol/kg p.o.) reduced the paw licking time in the first and second phase of the formalin test and decreased the nociceptive threshold variation in the Randall-Selitto test. Furthermore, this dose reduced oedema formation, leucocyte migration (specifically through reduction in polymorphonuclear cell movement) and increased mononuclear cells. MPO activity and the levels of pro-inflammatory cytokines TNF-α were decreased. Evaluation of PLA2 inhibition via the docking simulation revealed more interactions of LQFM043R(6) and LQFM044(7), data that corroborated the half-maximal inhibitory concentration (IC50) of PLA2 inhibition in vitro. Therefore, LQFM011(5), LQFM043(6) and LQFM044(7) were classified with the Globally Harmonized System of Classification and Labelling of Chemicals (GHS) as category 4.

All-Hydrocarbon Staples and Their Effect over Peptide Conformation under Different Force Fields.

Olefinic staples enhance α-helical content and conformational stability in peptides, maintaining a structural scaffold that allows the emulation of specific regions of protein surfaces for therapeutical purposes. The ability to anticipate the efficacy of adding a staple to a peptide through computational simulations may contribute to lowering the costs associated with rational drug design. We evaluated the capabilities of different force fields to reproduce the effect of all-hydrocarbon staples in molecular dynamics simulations. Using the AMBER99SB-ILDN, CHARMM36, and GROMOS54A7 force fields and two distinct initial conformations, we compared our results to experimentally obtained circular dichroism data. The GROMOS54A7 united-atom force field seems to be more accurate compared with all-atom force fields, despite being unable to reproduce the effect of the staple in some of the simulated systems. With further force field enhancements, MD simulations may be used to anticipate conformational effects of all-hydrocarbon staples in peptides.

Molecular docking and pharmacological/toxicological assessment of a new compound designed from celecoxib and paracetamol by molecular hybridization.

Nonsteroidal anti-inflammatory drugs are commonly used worldwide; however, they have several adverse effects, evidencing the need for the development of new, more effective and safe anti-inflammatory and analgesic drugs. This research aimed to design, synthesize and carry out a pharmacological/toxicological investigation of LQFM-102, which was designed from celecoxib and paracetamol by molecular hybridization. To evaluate the analgesic effect of this compound, we performed formalin-induced pain, hot plate and tail flick tests. The anti-inflammatory effect of LQFM-102 was evaluated in carrageenan-induced paw oedema and pleurisy tests. The biochemical markers indicative of toxicity-AST, ALT, GSH, urea and creatinine-as well as the index of gastric lesion after prolonged administration of LQFM-102 were also analyzed. In addition, the interaction of LQFM-102 with COX enzymes was evaluated by molecular docking. In all experimental protocols, celecoxib or paracetamol was used as a positive control at equimolar doses to LQFM-102. LQFM-102 reduced the pain induced by formalin in both phases of the test. However, this compound did not increase the latency to thermal stimuli in the hot plate and tail flick tests, suggesting an involvement of peripheral mechanisms in this effect. Furthermore, LQFM-102 reduced paw oedema, the number of polymorphonuclear cells, myeloperoxidase activity and TNF-α and IL-1ß levels. Another interesting finding was the absence of alterations in the markers of hepatic and renal toxicity or lesions of gastric mucosa. In molecular docking simulations, LQFM-102 interacted with the key residues for activity and potency of cyclooxygenase enzymes, suggesting an inhibition of the activity of these enzymes.

An Unusual Intramolecular Halogen Bond Guides Conformational Selection.

PIK-75 is a phosphoinositide-3-kinase (PI3K) α-isoform-selective inhibitor with high potency. Although published structure-activity relationship data show the importance of the NO2 and the Br substituents in PIK-75, none of the published studies could correctly determine the underlying reason for their importance. In this publication, we report the first X-ray crystal structure of PIK-75 in complex with the kinase GSK-3ß. The structure shows an unusual U-shaped conformation of PIK-75 within the active site of GSK-3ß that is likely stabilized by an atypical intramolecular Br⋅⋅⋅NO2 halogen bond. NMR and MD simulations show that this conformation presumably also exists in solution and leads to a binding-competent preorganization of the PIK-75 molecule, thus explaining its high potency. We therefore suggest that the site-specific incorporation of halogen bonds could be generally used to design conformationally restricted bioactive substances with increased potencies.

The role of Zn2+, dimerization and N-glycosylation in the interaction of Auxin-Binding Protein 1 (ABP1) with different auxins.

Auxin is critical for plant growth and development. The main natural auxin is indole-3-acetic acid (IAA), whereas 1-naphthalene acetic acid (NAA) is a synthetic form. Auxin-Binding Protein 1 (ABP1) specifically binds auxins, presumably playing roles as receptor in nontranscriptional cell responses. ABP1 structure was previously established from maize at 1.9 Å resolution. To gain further insight on ABP1 structural biology, this study was carried out employing molecular dynamics simulations of the complete models of the oligomeric glycosylated proteins from maize and Arabidopsis thaliana with or without auxins. In maize, both Zn2+ coordination and glycosylation promoted conformational stability and most of such stabilization effect was located on the N-terminal region. The α-helix of C-terminal regions in ABP1 of both species unfolded during simulations, assuming a more extended structure in maize. In Arabidopsis, the helix appeared more stable, being preserved in most of the monomeric simulations and unfolding when the protein was in the dimeric form. In Arabidopsis ABP1 bound to IAA or NAA, glycosylation structures arranged around the protein, covering the putative site of entrance or egress of auxin. NAA bound protein folding was more similar to the crystal structure showing higher stability compared to that of IAA bound. The molecular structural differences of ABP1 found between the species and auxin types indicate that this auxin-binding protein shows functional specificities in dicots and monocots, as well as in auxin type binding.

Homology modeling and molecular dynamics provide structural insights into tospovirus nucleoprotein.

BACKGROUND: Tospovirus is a plant-infecting genus within the family Bunyaviridae, which also includes four animal-infecting genera: Hantavirus, Nairovirus, Phlebovirus and Orthobunyavirus. Compared to these members, the structures of Tospovirus proteins still are poorly understood. Despite multiple studies have attempted to identify candidate N protein regions involved in RNA binding and protein multimerization for tospovirus using yeast two-hybrid systems (Y2HS) and site-directed mutagenesis, the tospovirus ribonucleocapsids (RNPs) remains largely uncharacterized at the molecular level and the lack of structural information prevents detailed insight into these interactions. RESULTS: Here we used the nucleoprotein structure of LACV (La Crosse virus-Orthobunyavirus) and molecular dynamics simulations to access the structure and dynamics of the nucleoprotein from tospovirus GRSV (Groundnut ringspot virus). The resulting model is a monomer composed by a flexible N-terminal and C-terminal arms and a globular domain with a positively charged groove in which RNA is deeply encompassed. This model allowed identifying the candidate amino acids residues involved in RNA interaction and N-N multimerization. Moreover, most residues predicted to be involved in these interactions are highly conserved among tospoviruses. CONCLUSIONS: Crucially, the interaction model proposed here for GRSV N is further corroborated by the all available mutational studies on TSWV (Tomato spotted wilt virus) N, so far. Our data will help designing further and more accurate mutational and functional studies of tospovirus N proteins. In addition, the proposed model may shed light on the mechanisms of RNP shaping and could allow the identification of essential amino acid residues as potential targets for tospovirus control strategies.

Glycosylation is crucial for a proper catalytic site organization in human glucocerebrosidase.

Gaucher disease, an autosomal recessive disorder, is caused by a deficiency of glucocerebrosidase (GCase) enzyme, a peripheral membrane-associated glycoprotein that hydrolyses glucosylceramide in lysosomes. Glycosylation is essential for the development of a catalytically active enzyme, specifically in the first site, located at Asn19. However, both the molecular basis of the relevance of N-glycosylation over GCase activity and the effects of glycosylation over its structure and dynamics are still not fully understood. Thus, the present work evaluated GCase enzyme in increasing glycosylation content using triplicate unbiased molecular dynamics simulations. Accordingly, the N-linked glycan chains caused local conformational stabilization effects over the protein, as well as in regions flanking the enzyme catalytic dyad. In the case of the Asn19-linked glycan, it also occurred around region 438-444, where one of the most prevalent GCase mutations is found. Markedly, an increasing catalytic dyad organization was related to increasing glycosylation contents, offering the first atomic-level explanation for the experimental observation that GCase activity is controlled by glycosylation, especially at Asn19.

Structure-function studies on jaburetox, a recombinant insecticidal peptide derived from jack bean (Canavalia ensiformis) urease.

BACKGROUND: Ureases are metalloenzymes involved in defense mechanisms in plants. The insecticidal activity of Canavalia ensiformis (jack bean) ureases relies partially on an internal 10kDa peptide generated by enzymatic hydrolysis of the protein within susceptible insects. A recombinant version of this peptide, jaburetox, exhibits insecticidal, antifungal and membrane-disruptive properties. Molecular modeling of jaburetox revealed a prominent ß-hairpin motif consistent with either neurotoxicity or pore formation. METHODS: Aiming to identify structural motifs involved in its effects, mutated versions of jaburetox were built: 1) a peptide lacking the ß-hairpin motif (residues 61-74), JbtxΔ-ß; 2) a peptide corresponding the N-terminal half (residues 1-44), Jbtx N-ter, and 3) a peptide corresponding the C-terminal half (residues 45-93), Jbtx C-ter. RESULTS: 1) JbtxΔ-ß disrupts liposomes, and exhibited entomotoxic effects similar to the whole peptide, suggesting that the ß-hairpin motif is not a determinant of these biological activities; 2) both Jbtx C-ter and Jbtx N-ter disrupted liposomes, the C-terminal peptide being the most active; and 3) while Jbtx N-ter persisted to be biologically active, Jbtx C-ter was less active when tested on different insect preparations. Molecular modeling and dynamics were applied to the urease-derived peptides to complement the structure-function analysis. MAJOR CONCLUSIONS: The N-terminal portion of the Jbtx carries the most important entomotoxic domain which is fully active in the absence of the ß-hairpin motif. Although the ß-hairpin contributes to some extent, probably by interaction with insect membranes, it is not essential for the entomotoxic properties of Jbtx. GENERAL SIGNIFICANCE: Jbtx represents a new type of insecticidal and membrane-active peptide.

Dynamics on human Toll-like receptor 4 complexation to MD-2: the coreceptor stabilizing function.

The interaction between human Toll-like receptor 4 (hTLR4) and its coreceptor, myeloid differentiation factor 2 (MD-2), is important in Gram-negative bacteria lipopolysaccharide (LPS) recognition. In this process, MD-2 recognizes LPS and promotes the dimerization of the complex hTLR4-MD-2-LPS, triggering an intracellular immune signaling. In this study, we employed distinct computational methods to explore the dynamical properties of the hTLR4-MD-2 complex and investigated the implications of the coreceptor complexation to the structural biology of hTLR4. We characterized both global and local dynamics of free and MD-2 complexed hTLR4, in both (hTLR4-MD-2)1 and (hTLR4-MD-2)2 states. Both molecular dynamics and normal mode analysis reveled a stabilization of the terminal regions of hTLR4 upon complexation to MD-2. We are able to identify conserved important residues involved on the hTLR4-MD-2 interaction dynamics and disclose C-terminal motions that may be associated to the signaling process upon oligomerization.

In silico Investigation of the PglB Active Site Reveals Transient Catalytic States and Octahedral Metal Ion Coordination.

The last step of the bacterial N-glycosylation pathway involves PglB, an oligosaccharyltransferase, which is responsible for the en bloc transfer of a fully assembled oligosaccharide chain to a protein possessing the extended motif D/E-X-N-X-S/T. Recently, this molecule had its full structure elucidated, enabling the description of its domains and the proposition of a catalytic mechanism. By employing molecular dynamics simulations, we were able to evaluate structural aspects of PglB, suggesting prevalent motions that may bring insights into the mechanism of the glycosylated peptide detachment. Additionally, we identified transient states at the catalytic site, in which the previously described carboxamide twisting mechanism was observed. Aided by quantum mechanics calculations for each different conformational states of the catalytic site, we determined the presence of an octahedral metal coordination, along with the presence of one water molecule at the catalytic site.

Structural glycobiology of human α1-acid glycoprotein and its implications for pharmacokinetics and inflammation.

Human α1-acid glycoprotein (AGP) is an abundant human plasma glycoprotein that may be N-glycosylated at five positions. AGP plays important roles on pharmacokinetics and can rise up to 5-fold in inflammatory events. In such events, the glycan chains attached to Asn54, Asn75 and Asn85 may become fucosylated, originating a sialyl-Lewis X epitope. This epitope, in turn, can bind selectin proteins. Such interplay is important for immunomodulation. While the X-ray structure of unglycosylated AGP has been reported, the absence of the glycan chains hampered the further insights into its structural biology and, ultimately, into its biological function. Thus, the current work intends to contribute in the characterization of the structural glycobiology and function of AGP by building a structural model of its fully glycosylated form, taking into account the different glycoforms that are found in vivo. The obtained data points to the absence of a major influence of glycosylation on AGP's secondary structure, in agreement with crystallography observations. However, the glycan chains seem able to interfere with the protein dynamics, mainly at the AGP-ligand-binding site, indicating a possible role in its complexation to drugs and other bioactive compounds. By examining the influence of fucosylation on AGP structure and binding to selectins, it is proposed that the latter may bind to glycan chains linked to Asn54 and Asn75, and that this binding may involve other glycans, such as the one attached to Asn15. These results point to an increased participation of carbohydrates on the observed AGP roles in pharmacokinetics and inflammation.

Ribifolin, an orbitide from Jatropha ribifolia, and its potential antimalarial activity.

A new orbitide named ribifolin was isolated and characterized from Jatropha ribifolia using mass spectrometry, NMR spectroscopy, quantitative amino acid analysis, molecular dynamics/simulated annealing, and Raman optical activity measurements and calculations. Ribifolin (1) and its linear form (1a) were synthesized by solid-phase peptide synthesis, followed by evaluation of its antiplasmodial and cytotoxicity activities. Compound 1 was moderately effective (IC50 = 42 µM) against the Plasmodium falciparum strain 3D7.

Natural Plant Alkaloid (Emetine) Inhibits HIV-1 Replication by Interfering with Reverse Transcriptase Activity.

Ipecac alkaloids are secondary metabolites produced in the medicinal plant Psychotria ipecacuanha. Emetine is the main alkaloid of ipecac and one of the active compounds in syrup of Ipecac with emetic property. Here we evaluated emetine's potential as an antiviral agent against Human Immunodeficiency Virus. We performed in vitro Reverse Transcriptase (RT) Assay and Natural Endogenous Reverse Transcriptase Activity Assay (NERT) to evaluate HIV RT inhibition. Emetine molecular docking on HIV-1 RT was also analyzed. Phenotypic assays were performed in non-lymphocytic and in Peripheral Blood Mononuclear Cells (PBMC) with HIV-1 wild-type and HIV-harboring RT-resistant mutation to Nucleoside Reverse Transcriptase Inhibitors (M184V). Our results showed that HIV-1 RT was blocked in the presence of emetine in both models: in vitro reactions with isolated HIV-1 RT and intravirion, measured by NERT. Emetine revealed a strong potential of inhibiting HIV-1 replication in both cellular models, reaching 80% of reduction in HIV-1 infection, with low cytotoxic effect. Emetine also blocked HIV-1 infection of RT M184V mutant. These results suggest that emetine is able to penetrate in intact HIV particles, and bind and block reverse transcription reaction, suggesting that it can be used as anti-HIV microbicide. Taken together, our findings provide additional pharmacological information on the potential therapeutic effects of emetine.