The IDSM mass spectrometry extension: searching mass spectra using SPARQL.

SUMMARY: The Integrated Database of Small Molecules (IDSM) integrates data from small-molecule datasets, making them accessible through the SPARQL query language. Its unique feature is the ability to search for compounds through SPARQL based on their molecular structure. We extended IDSM to enable mass spectra databases to be integrated and searched for based on mass spectrum similarity. As sources of mass spectra, we employed the MassBank of North America database and the In Silico Spectral Database of natural products. AVAILABILITY AND IMPLEMENTATION: The extension is an integral part of IDSM, which is available at https://idsm.elixir-czech.cz. The manual and usage examples are available at https://idsm.elixir-czech.cz/docs/ms. The source codes of all IDSM parts are available under open-source licences at https://github.com/idsm-src.

PAPerFly: Partial Assembly-based Peak Finder for ab initio binding site reconstruction.

BACKGROUND: The specific recognition of a DNA locus by a given transcription factor is a widely studied issue. It is generally agreed that the recognition can be influenced not only by the binding motif but by the larger context of the binding site. In this work, we present a novel heuristic algorithm that can reconstruct the unique binding sites captured in a sequencing experiment without using the reference genome. RESULTS: We present PAPerFly, the Partial Assembly-based Peak Finder, a tool for the binding site and binding context reconstruction from the sequencing data without any prior knowledge. This tool operates without the need to know the reference genome of the respective organism. We employ algorithmic approaches that are used during genome assembly. The proposed algorithm constructs a de Bruijn graph from the sequencing data. Based on this graph, sequences and their enrichment are reconstructed using a novel heuristic algorithm. The reconstructed sequences are aligned and the peaks in the sequence enrichment are identified. Our approach was tested by processing several ChIP-seq experiments available in the ENCODE database and comparing the results of Paperfly and standard methods. CONCLUSIONS: We show that PAPerFly, an algorithm tailored for experiment analysis without the reference genome, yields better results than an aggregation of ChIP-seq agnostic tools. Our tool is freely available at https://github.com/Caeph/paperfly/ or on Zenodo ( https://doi.org/10.5281/zenodo.7116424 ).

Amino Acid Interactions (INTAA) web server v2.0: a single service for computation of energetics and conservation in biomolecular 3D structures.

Interactions among amino acid residues are the principal contributor to the stability of the three-dimensional structure of a protein. The Amino Acid Interactions (INTAA) web server (https://bioinfo.uochb.cas.cz/INTAA/) has established itself as a unique computational resource, which enables users to calculate the contribution of individual residues in a biomolecular structure to its total energy using a molecular mechanical scoring function. In this update, we describe major additions to the web server which help solidify its position as a robust, comprehensive resource for biomolecular structure analysis. Importantly, a new continuum solvation model was introduced, allowing more accurate representation of electrostatic interactions in aqueous media. In addition, a low-overhead pipeline for the estimation of evolutionary conservation in protein chains has been added. New visualization options were introduced as well, allowing users to easily switch between and interrelate the energetic and evolutionary views of the investigated structures.

TRPM5 Channel Binds Calcium-Binding Proteins Calmodulin and S100A1.

Melastatin transient receptor potential (TRPM) channels belong to one of the most significant subgroups of the transient receptor potential (TRP) channel family. Here, we studied the TRPM5 member, the receptor exposed to calcium-mediated activation, resulting in taste transduction. It is known that most TRP channels are highly modulated through interactions with extracellular and intracellular agents. The binding sites for these ligands are usually located at the intracellular N- and C-termini of the TRP channels, and they can demonstrate the character of an intrinsically disordered protein (IDP), which allows such a region to bind various types of molecules. We explored the N-termini of TRPM5 and found the intracellular regions for calcium-binding proteins (CBPs) the calmodulin (CaM) and calcium-binding protein S1 (S100A1) by in vitro binding assays. Furthermore, molecular docking and molecular dynamics simulations (MDs) of the discovered complexes confirmed their known common binding interface patterns and the uniqueness of the basic residues present in the TRPM binding regions for CaM/S100A1.

Fusion of two unrelated protein domains in a chimera protein and its 3D prediction: Justification of the x-ray reference structures as a prediction benchmark.

Proteins are naturally formed by domains edging their functional and structural properties. A domain out of the context of an entire protein can retain its structure and to some extent also function on its own. These properties rationalize construction of artificial fusion multidomain proteins with unique combination of various functions. Information on the specific functional and structural characteristics of individual domains in the context of new artificial fusion proteins is inevitably encoded in sequential order of composing domains defining their mutual spatial positions. So the challenges in designing new proteins with new domain combinations lie dominantly in structure/function prediction and its context dependency. Despite the enormous body of publications on artificial fusion proteins, the task of their structure/function prediction is complex and nontrivial. The degree of spatial freedom facilitated by a linker between domains and their mutual orientation driven by noncovalent interactions is beyond a simple and straightforward methodology to predict their structure with reasonable accuracy. In the presented manuscript, we tested methodology using available modeling tools and computational methods. We show that the process and methodology of such prediction are not straightforward and must be done with care even when recently introduced AlphaFold II is used. We also addressed a question of benchmarking standards for prediction of multidomain protein structures-x-ray or Nuclear Magnetic Resonance experiments. On the study of six two-domain protein chimeras as well as their composing domains and their x-ray structures selected from PDB, we conclude that the major obstacle for justified prediction is inappropriate sampling of the conformational space by the explored methods. On the other hands, we can still address particular steps of the methodology and improve the process of chimera proteins prediction.

Iterative Landmark-Based Umbrella Sampling (ILBUS) Protocol for Sampling of Conformational Space of Biomolecules.

Computer simulations of biomolecules such as molecular dynamics often suffer from insufficient sampling. Due to limited computational resources, insufficient sampling prevents obtaining proper equilibrium distributions of observed properties. To deal with this problem, we proposed a simulation protocol for efficient resampling of collected off-equilibrium trajectories. These trajectories are utilized for the initial mapping of the conformational space, which is later properly resampled by the introduced Iterative Landmark-Based Umbrella Sampling (ILBUS) method. Reconstruction of static equilibrium properties is achieved by the multistate Bennett acceptance ratio (MBAR) method, which enables efficient use of simulated data. The ILBUS protocol is geometry-based and does not demand any additional collective variable or a dimensional-reduction technique. The only requirement is a set of suitably spaced reference conformations, which serve as landmarks in the mapped conformational space. Additionally, the ILBUS protocol encompasses an iterative process that optimizes the force constant used in the umbrella sampling simulation. Such tuning is an inherent feature of the protocol and does not need to be performed by the user in advance. Furthermore, even the simulations with suboptimal force constants can be used in estimates by MBAR. We demonstrate the feasibility and the performance of this approach in the study of the conformational landscape of the alanine dipeptide, met-enkephalin, and adenylate kinase.

ShinySOM: graphical SOM-based analysis of single-cell cytometry data.

SUMMARY: ShinySOM offers a user-friendly interface for reproducible, high-throughput analysis of high-dimensional flow and mass cytometry data guided by self-organizing maps. The software implements a FlowSOM-style workflow, with improvements in performance, visualizations and data dissection possibilities. The outputs of the analysis include precise statistical information about the dissected samples, and R-compatible metadata useful for the batch processing of large sample volumes. AVAILABILITY AND IMPLEMENTATION: ShinySOM is free and open-source, available online at gitlab.com/exaexa/ShinySOM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exquisite ligand stereoselectivity of a Drosophila juvenile hormone receptor contrasts with its broad agonist repertoire.

The sesquiterpenoid juvenile hormone (JH) is vital to insect development and reproduction. Intracellular JH receptors have recently been established as basic helix-loop-helix transcription factor (bHLH)/PAS proteins in Drosophila melanogaster known as germ cell-expressed (Gce) and its duplicate paralog, methoprene-tolerant (Met). Upon binding JH, Gce/Met activates its target genes. Insects possess multiple native JH homologs whose molecular activities remain unexplored, and diverse synthetic compounds including insecticides exert JH-like effects. How the JH receptor recognizes its ligands is unknown. To determine which structural features define an active JH receptor agonist, we tested several native JHs and their nonnative geometric and optical isomers for the ability to bind the Drosophila JH receptor Gce, to induce Gce-dependent transcription, and to affect the development of the fly. Our results revealed high ligand stereoselectivity of the receptor. The geometry of the JH skeleton, dictated by two stereogenic double bonds, was the most critical feature followed by the presence of an epoxide moiety at a terminal position. The optical isomerism at carbon C11 proved less important even though Gce preferentially bound a natural JH enantiomer. The results of receptor-ligand-binding and cell-based gene activation assays tightly correlated with the ability of different geometric JH isomers to induce gene expression and morphogenetic effects in the developing insects. Molecular modeling supported the requirement for the proper double-bond geometry of JH, which appears to be its major selective mechanism. The strict stereoselectivity of Gce toward the natural hormone contrasts with the high potency of synthetic Gce agonists of disparate chemistries.

3DPatch: fast 3D structure visualization with residue conservation.

Summary: Amino acid residues showing above background levels of conservation are often indicative of functionally significant regions within a protein. Understanding how the sequence conservation profile relates in space requires projection onto a protein structure, a potentially time-consuming process. 3DPatch is a web application that streamlines this task by automatically generating multiple sequence alignments (where appropriate) and finding structural homologs, presenting the user with a choice of structures matching their query, annotated with residue conservation scores in a matter of seconds. Availability and implementation: 3DPatch is written in JavaScript and is freely available at http://www.skylign.org/3DPatch/. Mozilla Firefox, Google Chrome, and Safari web browsers are supported. Source code is available under MIT license at https://github.com/davidjakubec/3DPatch. Supplementary information: Supplementary data are available at Bioinformatics online.

Hydrophobic Amino Acids as Universal Elements of Protein-Induced DNA Structure Deformation.

Interaction with the DNA minor groove is a significant contributor to specific sequence recognition in selected families of DNA-binding proteins. Based on a statistical analysis of 3D structures of protein-DNA complexes, we propose that distortion of the DNA minor groove resulting from interactions with hydrophobic amino acid residues is a universal element of protein-DNA recognition. We provide evidence to support this by associating each DNA minor groove-binding amino acid residue with the local dimensions of the DNA double helix using a novel algorithm. The widened DNA minor grooves are associated with high GC content. However, some AT-rich sequences contacted by hydrophobic amino acids (e.g., phenylalanine) display extreme values of minor groove width as well. For a number of hydrophobic amino acids, distinct secondary structure preferences could be identified for residues interacting with the widened DNA minor groove. These results hold even after discarding the most populous families of minor groove-binding proteins.

Mapping of CaM, S100A1 and PIP2-Binding Epitopes in the Intracellular N- and C-Termini of TRPM4.

Molecular determinants of the binding of various endogenous modulators to transient receptor potential (TRP) channels are crucial for the understanding of necessary cellular pathways, as well as new paths for rational drug designs. The aim of this study was to characterise interactions between the TRP cation channel subfamily melastatin member 4 (TRPM4) and endogenous intracellular modulators-calcium-binding proteins (calmodulin (CaM) and S100A1) and phosphatidylinositol 4, 5-bisphosphate (PIP2). We have found binding epitopes at the N- and C-termini of TRPM4 shared by CaM, S100A1 and PIP2. The binding affinities of short peptides representing the binding epitopes of N- and C-termini were measured by means of fluorescence anisotropy (FA). The importance of representative basic amino acids and their combinations from both peptides for the binding of endogenous TRPM4 modulators was proved using point alanine-scanning mutagenesis. In silico protein-protein docking of both peptides to CaM and S100A1 and extensive molecular dynamics (MD) simulations enabled the description of key stabilising interactions at the atomic level. Recently solved cryo-Electron Microscopy (EM) structures made it possible to put our findings into the context of the entire TRPM4 channel and to deduce how the binding of these endogenous modulators could allosterically affect the gating of TRPM4. Moreover, both identified binding epitopes seem to be ideally positioned to mediate the involvement of TRPM4 in higher-order hetero-multimeric complexes with important physiological functions.

Characterization of AMBN I and II Isoforms and Study of Their Ca2+-Binding Properties.

Ameloblastin (Ambn) as an intrinsically disordered protein (IDP) stands for an important role in the formation of enamel-the hardest biomineralized tissue commonly formed in vertebrates. The human ameloblastin (AMBN) is expressed in two isoforms: full-length isoform I (AMBN ISO I) and isoform II (AMBN ISO II), which is about 15 amino acid residues shorter than AMBN ISO I. The significant feature of AMBN-its oligomerization ability-is enabled due to a specific sequence encoded by exon 5 present at the N-terminal part in both known isoforms. In this study, we characterized AMBN ISO I and AMBN ISO II by biochemical and biophysical methods to determine their common features and differences. We confirmed that both AMBN ISO I and AMBN ISO II form oligomers in in vitro conditions. Due to an important role of AMBN in biomineralization, we further addressed the calcium (Ca2+)-binding properties of AMBN ISO I and ISO II. The binding properties of AMBN to Ca2+ may explain the role of AMBN in biomineralization and more generally in Ca2+ homeostasis processes.

Amino Acid Interaction (INTAA) web server.

Large biomolecules-proteins and nucleic acids-are composed of building blocks which define their identity, properties and binding capabilities. In order to shed light on the energetic side of interactions of amino acids between themselves and with deoxyribonucleotides, we present the Amino Acid Interaction web server (http://bioinfo.uochb.cas.cz/INTAA/). INTAA offers the calculation of the residue Interaction Energy Matrix for any protein structure (deposited in Protein Data Bank or submitted by the user) and a comprehensive analysis of the interfaces in protein-DNA complexes. The Interaction Energy Matrix web application aims to identify key residues within protein structures which contribute significantly to the stability of the protein. The application provides an interactive user interface enhanced by 3D structure viewer for efficient visualization of pairwise and net interaction energies of individual amino acids, side chains and backbones. The protein-DNA interaction analysis part of the web server allows the user to view the relative abundance of various configurations of amino acid-deoxyribonucleotide pairs found at the protein-DNA interface and the interaction energies corresponding to these configurations calculated using a molecular mechanical force field. The effects of the sugar-phosphate moiety and of the dielectric properties of the solvent on the interaction energies can be studied for the various configurations.

TRPM6 N-Terminal CaM- and S100A1-Binding Domains.

Transient receptor potential (TRPs) channels are crucial downstream targets of calcium signalling cascades. They can be modulated either by calcium itself and/or by calcium-binding proteins (CBPs). Intracellular messengers usually interact with binding domains present at the most variable TRP regions-N- and C-cytoplasmic termini. Calmodulin (CaM) is a calcium-dependent cytosolic protein serving as a modulator of most transmembrane receptors. Although CaM-binding domains are widespread within intracellular parts of TRPs, no such binding domain has been characterised at the TRP melastatin member-the transient receptor potential melastatin 6 (TRPM6) channel. Another CBP, the S100 calcium-binding protein A1 (S100A1), is also known for its modulatory activities towards receptors. S100A1 commonly shares a CaM-binding domain. Here, we present the first identified CaM and S100A1 binding sites at the N-terminal of TRPM6. We have confirmed the L520-R535 N-terminal TRPM6 domain as a shared binding site for CaM and S100A1 using biophysical and molecular modelling methods. A specific domain of basic amino acid residues (R526/R531/K532/R535) present at this TRPM6 domain has been identified as crucial to maintain non-covalent interactions with the ligands. Our data unambiguously confirm that CaM and S100A1 share the same binding domain at the TRPM6 N-terminus although the ligand-binding mechanism is different.

Sequence Versus Composition: What Prescribes IDP Biophysical Properties?

Intrinsically disordered proteins (IDPs) represent a distinct class of proteins and are distinguished from globular proteins by conformational plasticity, high evolvability and a broad functional repertoire. Some of their properties are reminiscent of early proteins, but their abundance in eukaryotes, functional properties and compositional bias suggest that IDPs appeared at later evolutionary stages. The spectrum of IDP properties and their determinants are still not well defined. This study compares rudimentary physicochemical properties of IDPs and globular proteins using bioinformatic analysis on the level of their native sequences and random sequence permutations, addressing the contributions of composition versus sequence as determinants of the properties. IDPs have, on average, lower predicted secondary structure contents and aggregation propensities and biased amino acid compositions. However, our study shows that IDPs exhibit a broad range of these properties. Induced fold IDPs exhibit very similar compositions and secondary structure/aggregation propensities to globular proteins, and can be distinguished from unfoldable IDPs based on analysis of these sequence properties. While amino acid composition seems to be a major determinant of aggregation and secondary structure propensities, sequence randomization does not result in dramatic changes to these properties, but for both IDPs and globular proteins seems to fine-tune the tradeoff between folding and aggregation.

Substituents at the C3' and C3'N positions are critical for taxanes to overcome acquired resistance of cancer cells to paclitaxel.

We tested the role of substituents at the C3' and C3'N positions of the taxane molecule to identify taxane derivatives capable of overcoming acquired resistance to paclitaxel. Paclitaxel-resistant sublines SK-BR-3/PacR and MCF-7/PacR as well as the original paclitaxel-sensitive breast cancer cell lines SK-BR-3 and MCF-7 were used for testing. Increased expression of the ABCB1 transporter was found to be involved in the acquired resistance. We tested three groups of taxane derivatives: (1) phenyl group at both C3' and C3'N positions, (2) one phenyl at one of the C3' and C3'N positions and a non-aromatic group at the second position, (3) a non-aromatic group at both C3' and C3'N positions. We found that the presence of phenyl groups at both C3' and C3'N positions is associated with low capability of overcoming acquired paclitaxel resistance compared to taxanes containing at least one non-aromatic substituent at the C3' and C3'N positions. The increase in the ATPase activity of ABCB1 transporter after the application of taxanes from the first group was found to be somewhat higher than after the application of taxanes from the third group. Molecular docking studies demonstrated that the docking score was the lowest, i.e. the highest binding affinity, for taxanes from the first group. It was intermediate for taxanes from the second group, and the highest for taxanes from the third group. We conclude that at least one non-aromatic group at the C3' and C3'N positions of the taxane structure, resulting in reduced affinity to the ABCB1 transporter, brings about high capability of taxane to overcome acquired resistance of breast cancer cells to paclitaxel, due to less efficient transport of the taxane compound out of the cancer cells.

Insights into Unfolded Proteins from the Intrinsic ϕ/ψ Propensities of the AAXAA Host-Guest Series.

Various host-guest peptide series are used by experimentalists as reference conformational states. One such use is as a baseline for random-coil NMR chemical shifts. Comparison to this random-coil baseline, through secondary chemical shifts, is used to infer protein secondary structure. The use of these random-coil data sets rests on the perception that the reference chemical shifts arise from states where there is little or no conformational bias. However, there is growing evidence that the conformational composition of natively and nonnatively unfolded proteins fail to approach anything that can be construed as random coil. Here, we use molecular dynamics simulations of an alanine-based host-guest peptide series (AAXAA) as a model of unfolded and denatured states to examine the intrinsic propensities of the amino acids. We produced ensembles that are in good agreement with the experimental NMR chemical shifts and confirm that the sampling of the 20 natural amino acids in this peptide series is be far from random. Preferences toward certain regions of conformational space were both present and dependent upon the environment when compared under conditions typically used to denature proteins, i.e., thermal and chemical denaturation. Moreover, the simulations allowed us to examine the conformational makeup of the underlying ensembles giving rise to the ensemble-averaged chemical shifts. We present these data as an intrinsic backbone propensity library that forms part of our Structural Library of Intrinsic Residue Propensities to inform model building, to aid in interpretation of experiment, and for structure prediction of natively and nonnatively unfolded states.

Artificial proteins as allosteric modulators of PDZ3 and SH3 in two-domain constructs: A computational characterization of novel chimeric proteins.

Artificial multidomain proteins with enhanced structural and functional properties can be utilized in a broad spectrum of applications. The design of chimeric fusion proteins utilizing protein domains or one-domain miniproteins as building blocks is an important advancement for the creation of new biomolecules for biotechnology and medical applications. However, computational studies to describe in detail the dynamics and geometry properties of two-domain constructs made from structurally and functionally different proteins are lacking. Here, we tested an in silico design strategy using all-atom explicit solvent molecular dynamics simulations. The well-characterized PDZ3 and SH3 domains of human zonula occludens (ZO-1) (3TSZ), along with 5 artificial domains and 2 types of molecular linkers, were selected to construct chimeric two-domain molecules. The influence of the artificial domains on the structure and dynamics of the PDZ3 and SH3 domains was determined using a range of analyses. We conclude that the artificial domains can function as allosteric modulators of the PDZ3 and SH3 domains. Proteins 2016; 84:1358-1374. © 2016 Wiley Periodicals, Inc.

Computational methods for the description of pharmacologically relevant platinum complexes--molecular structure and bond dissociation.

Cancer is after cardiovascular disease the most frequent cause of death in Europe. In 28 of 53 countries considered in this area it is already the leading cause of death and expected to gain even more importance until the year 2020. Amongst the large arsenal of different anti-cancer drugs, platinum drugs belong to the first developed anticancer drugs and still have a large impact on cancer therapy. Nevertheless therapy with platinum-anticancer drugs is accompanied by severe adverse effects caused by frequent interactions with the amino acids of different human proteins. Computational chemistry offers methods to study such interactions and even those of not yet synthesized drugs in silico. For such studies a profound knowledge of the prediction quality of various computational methods towards platinum-drug-like complexes is necessary. By this article we are aiming on delivering important accuracy information of the frequently used computational methods. Most important findings are the high performance of the double hybrid functional B2PLYP for the calculation of geometries, even in small basis sets, followed by BP86 and PBE and the still acceptable performance of the semi-empirical Method PM6-D3H4X for extremely large systems. To follow absolute energies of the dissociation process, LPNO-CEPA and B3LYP-D3 can be suggested while SCS-MP2 shows an extremely narrow standard deviation and a low maximum error, which make it an ideal candidate for relative energy calculations in the exploration of reaction mechanisms.

Regioselective Palmitoylation of 9-(2,3-Dihydroxy- propyl)adenine Catalyzed by a Glycopolymer-enzyme Conjugate.

The enzymatic regioselective monopalmitoylation of racemic 9-(2,3-dihydroxypropyl)- adenine (DHPA), an approved antiviral agent, has been performed by an immobilized form of Candida antarctica B lipase (CAL-B) using a 4:1 DMF/hexane mixture as the reaction medium. To improve the chemical yield of the desired monopalmitoylation reaction, solid-phase chemical modifications of the lipase were evaluated. The reaction yield was successfully increased obtaining 100% product after a second treatment of the product solution with fresh immobilised chemically glycosylated-CAL-B.