Microsimulation based quantitative analysis of COVID-19 management strategies.

Pandemic management requires reliable and efficient dynamical simulation to predict and control disease spreading. The COVID-19 (SARS-CoV-2) pandemic is mitigated by several non-pharmaceutical interventions, but it is hard to predict which of these are the most effective for a given population. We developed the computationally effective and scalable, agent-based microsimulation framework PanSim, allowing us to test control measures in multiple infection waves caused by the spread of a new virus variant in a city-sized societal environment using a unified framework fitted to realistic data. We show that vaccination strategies prioritising occupational risk groups minimise the number of infections but allow higher mortality while prioritising vulnerable groups minimises mortality but implies an increased infection rate. We also found that intensive vaccination along with non-pharmaceutical interventions can substantially suppress the spread of the virus, while low levels of vaccination, premature reopening may easily revert the epidemic to an uncontrolled state. Our analysis highlights that while vaccination protects the elderly from COVID-19, a large percentage of children will contract the virus, and we also show the benefits and limitations of various quarantine and testing scenarios. The uniquely detailed spatio-temporal resolution of PanSim allows the design and testing of complex, specifically targeted interventions with a large number of agents under dynamically changing conditions.

Nuclear architecture dictates HIV-1 integration site selection.

Long-standing evidence indicates that human immunodeficiency virus type 1 (HIV-1) preferentially integrates into a subset of transcriptionally active genes of the host cell genome. However, the reason why the virus selects only certain genes among all transcriptionally active regions in a target cell remains largely unknown. Here we show that HIV-1 integration occurs in the outer shell of the nucleus in close correspondence with the nuclear pore. This region contains a series of cellular genes, which are preferentially targeted by the virus, and characterized by the presence of active transcription chromatin marks before viral infection. In contrast, the virus strongly disfavours the heterochromatic regions in the nuclear lamin-associated domains and other transcriptionally active regions located centrally in the nucleus. Functional viral integrase and the presence of the cellular Nup153 and LEDGF/p75 integration cofactors are indispensable for the peripheral integration of the virus. Once integrated at the nuclear pore, the HIV-1 DNA makes contact with various nucleoporins; this association takes part in the transcriptional regulation of the viral genome. These results indicate that nuclear topography is an essential determinant of the HIV-1 life cycle.

Motif oriented high-resolution analysis of ChIP-seq data reveals the topological order of CTCF and cohesin proteins on DNA.

BACKGROUND: ChIP-seq provides a wealth of information on the approximate location of DNA-binding proteins genome-wide. It is known that the targeted motifs in most cases can be found at the peak centers. A high resolution mapping of ChIP-seq peaks could in principle allow the fine mapping of the protein constituents within protein complexes, but the current ChIP-seq analysis pipelines do not target the basepair resolution strand specific mapping of peak summits. RESULTS: The approach proposed here is based on i) locating regions that are bound by a sufficient number of proteins constituting a complex; ii) determining the position of the underlying motif using either a direct or a de novo motif search approach; and iii) determining the exact location of the peak summits with respect to the binding motif in a strand specific manner. We applied this method for analyzing the CTCF/cohesin complex, which holds together DNA loops. The relative positions of the constituents of the complex were determined with one-basepair estimated accuracy. Mapping the positions on a 3D model of DNA made it possible to deduce the approximate local topology of the complex that allowed us to predict how the CTCF/cohesin complex locks the DNA loops. As the positioning of the proteins was not compatible with previous models of loop closure, we proposed a plausible "double embrace" model in which the DNA loop is held together by two adjacent cohesin rings in such a way that the ring anchored by CTCF to one DNA duplex encircles the other DNA double helix and vice versa. CONCLUSIONS: A motif-centered, strand specific analysis of ChIP-seq data improves the accuracy of determining peak positions. If a genome contains a large number of binding sites for a given protein complex, such as transcription factor heterodimers or transcription factor/cofactor complexes, the relative position of the constituent proteins on the DNA can be established with an accuracy that allow one to deduce the local topology of the protein complex. The proposed high resolution mapping approach of ChIP-seq data is applicable for detecting the contact topology of DNA-binding protein complexes.

PTMTreeSearch: a novel two-stage tree-search algorithm with pruning rules for the identification of post-translational modification of proteins in MS/MS spectra.

MOTIVATION: Tandem mass spectrometry has become a standard tool for identifying post-translational modifications (PTMs) of proteins. Algorithmic searches for PTMs from tandem mass spectrum data (MS/MS) tend to be hampered by noisy data as well as by a combinatorial explosion of search space. This leads to high uncertainty and long search-execution times. RESULTS: To address this issue, we present PTMTreeSearch, a new algorithm that uses a large database of known PTMs to identify PTMs from MS/MS data. For a given peptide sequence, PTMTreeSearch builds a computational tree wherein each path from the root to the leaves is labeled with the amino acids of a peptide sequence. Branches then represent PTMs. Various empirical tree pruning rules have been designed to decrease the search-execution time by eliminating biologically unlikely solutions. PTMTreeSearch first identifies a relatively small set of high confidence PTM types, and in a second stage, performs a more exhaustive search on this restricted set using relaxed search parameter settings. An analysis of experimental data shows that using the same criteria for false discovery, PTMTreeSearch annotates more peptides than the current state-of-the-art methods and PTM identification algorithms, and achieves this at roughly the same execution time. PTMTreeSearch is implemented as a plugable scoring function in the X!Tandem search engine. AVAILABILITY: The source code of PTMTreeSearch and a demo server application can be found at http://net.icgeb.org/ptmtreesearch

The basic helix-loop-helix region of the transcriptional repressor hairy and enhancer of split 1 is preorganized to bind DNA.

Hairy and enhancer of split 1, one of the main downstream effectors in Notch signaling, is a transcriptional repressor of the basic helix-loop-helix (bHLH) family. Using nuclear magnetic resonance methods, we have determined the structure and dynamics of a recombinant protein, H1H, which includes an N-terminal segment, b1, containing functionally important phosphorylation sites, the basic region b2, required for binding to DNA, and the HLH domain. We show that a proline residue in the sequence divides the protein in two parts, a flexible and disordered N-terminal region including b1 and a structured, mainly helical region comprising b2 and the HLH domain. Binding of H1H to a double strand DNA oligonucleotide was monitored through the chemical shift perturbation of backbone amide resonances, and showed that the interaction surface involves not only the b2 segment but also several residues in the b1 and HLH regions.

Chemical rule-based filtering of MS/MS spectra.

MOTIVATION: Identification of proteins by mass spectrometry-based proteomics requires automated interpretation of peptide tandem mass spectrometry spectra. The effectiveness of peptide identification can be greatly improved by filtering out extraneous noise peaks before the subsequent database searching steps. RESULTS: Here we present a novel chemical rule-based filtering algorithm, termed CRF, which makes use of the predictable patterns (rules) of collision-induced peptide fragmentation. The algorithm selects peak pairs that obey the common fragmentation rules within plausible limits of mass tolerance as well as peak intensity and produces spectra that can be subsequently submitted to any search engine. CRF increases the positive predictive value and decreases the number of random matches and thus improves performance by 15-20% in terms of peptide annotation using search engines, such as X!Tandem. Importantly, the algorithm also achieves data compression rates of ∼75%. AVAILABILITY: The MATLAB source code and a web server are available at http://hydrax.icgeb.trieste.it/CRFilter/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Flexibility of the PDZ-binding motif in the micelle-bound form of Jagged-1 cytoplasmic tail.

Human Jagged-1, one of the ligands of Notch receptors, is a transmembrane protein composed of a large extracellular region and a 125-residue cytoplasmic tail which bears a C-terminal PDZ recognition motif. To investigate the interaction between Jagged-1 cytoplasmic tail and the inner leaflet of the plasma membrane we determined, by solution NMR, the secondary structure and dynamics of the recombinant protein corresponding to the intracellular region of Jagged-1, J1_tmic, bound to negatively charged lysophospholipid micelles. NMR showed that the PDZ binding motif is preceded by four alpha-helical segments and that, despite the extensive interaction between J1_tmic and the micelle, the PDZ binding motif remains highly flexible. Binding of J1_tmic to negatively charged, but not to zwitterionic vesicles, was confirmed by surface plasmon resonance. To study the PDZ binding region in more detail, we prepared a peptide corresponding to the last 24 residues of Jagged-1, J1C24, and different phosphorylated variants of it. J1C24 displays a marked helical propensity and undergoes a coil-helix transition in the presence of negatively charged, but not zwitterionic, lysophospholipid micelles. Phosphorylation at different positions drastically decreases the helical propensity of the peptides and abolishes the coil-helix transition triggered by lysophospholipid micelles. We propose that phosphorylation of residues upstream of the PDZ binding motif may shift the equilibrium from an ordered, membrane-bound, interfacial form of Jagged-1 C-terminal region to a more disordered form with an increased accessibility of the PDZ recognition motif, thus playing an indirect role in the interaction between Jagged-1 and the PDZ-containing target protein.

Jagged-1 juxtamembrane region: biochemical characterization and cleavage by ADAM17 (TACE) catalytic domain.

Ectodomain shedding of membrane receptors and ligands carried out by ADAMs (A disintegrin and metalloprotease) plays a major role in several signaling pathways, including Notch. The grounds of substrate recognition, however, are poorly understood. We demonstrate that a recombinant protein corresponding to the juxtamembrane region of Jagged-1, one of the Notch ligands, behaves as a structured module and is cleaved by ADAM17 catalytic domain at E1054. A short synthetic peptide is cleaved at the same site but at a much higher rate, implying that the structure of the cleavage site in the native protein is a key determinant for substrate recognition. We also show that an Alagille syndrome-associated mutation near E1054 increases the cleavage rate, which suggests that this mutation may lead to an unbalance in Notch signaling due to a higher level of Jagged-1 shedding.

The organization of the quorum sensing luxI/R family genes in Burkholderia.

Members of the Burkholderia genus of Proteobacteria are capable of living freely in the environment and can also colonize human, animal and plant hosts. Certain members are considered to be clinically important from both medical and veterinary perspectives and furthermore may be important modulators of the rhizosphere. Quorum sensing via N-acyl homoserine lactone signals (AHL QS) is present in almost all Burkholderia species and is thought to play important roles in lifestyle changes such as colonization and niche invasion. Here we present a census of AHL QS genes retrieved from public databases and indicate that the local arrangement (topology) of QS genes, their location within chromosomes and their gene neighborhoods show characteristic patterns that differ between the known Burkholderia clades. In sequence phylogenies, AHL QS genes seem to cluster according to the local gene topology rather than according to the species, which suggests that the basic topology types were present prior to the appearance of current Burkholderia species. The data are available at http://net.icgeb.org/burkholderia/.

Wastewater-based modeling, reconstruction, and prediction for COVID-19 outbreaks in Hungary caused by highly immune evasive variants.

(MOTIVATION): Wastewater-based epidemiology (WBE) has emerged as a promising approach for monitoring the COVID-19 pandemic, since the measurement process is cost-effective and is exposed to fewer potential errors compared to other indicators like hospitalization data or the number of detected cases. Consequently, WBE was gradually becoming a key tool for epidemic surveillance and often the most reliable data source, as the intensity of clinical testing for COVID-19 drastically decreased by the third year of the pandemic. Recent results suggests that the model-based fusion of wastewater measurements with clinical data and other indicators is essential in future epidemic surveillance. (METHOD): In this work, we developed a wastewater-based compartmental epidemic model with a two-phase vaccination dynamics and immune evasion. We proposed a multi-step optimization-based data assimilation method for epidemic state reconstruction, parameter estimation, and prediction. The computations make use of the measured viral load in wastewater, the available clinical data (hospital occupancy, delivered vaccine doses, and deaths), the stringency index of the official social distancing rules, and other measures. The current state assessment and the estimation of the current transmission rate and immunity loss allow a plausible prediction of the future progression of the pandemic. (RESULTS): Qualitative and quantitative evaluations revealed that the contribution of wastewater data in our computational epidemiological framework makes predictions more reliable. Predictions suggest that at least half of the Hungarian population has lost immunity during the epidemic outbreak caused by the BA.1 and BA.2 subvariants of Omicron in the first half of 2022. We obtained a similar result for the outbreaks caused by the subvariant BA.5 in the second half of 2022. (APPLICABILITY): The proposed approach has been used to support COVID management in Hungary and could be customized for other countries as well.

The virtue of temperance: built-in negative regulators of quorum sensing in Pseudomonas.

Many bacteria are now believed to produce small signal molecules in order to communicate in a process called quorum sensing (QS), which mediates cooperative traits and a co-ordinated behaviour. Pseudomonads have been extensively studied for their QS response highlighting that it plays a major role in determining their lifestyle. The main QS signal molecules produced by Pseudomonas belong to the family of N-acyl-homoserine lactones (AHLs); these are synthesized by a LuxI-family synthase and sensed by a LuxR-family regulator. Most often in Pseudomonas, repressor genes intergenically located between luxI and luxR form an integral part of QS system. Recent studies have highlighted an important role of these repressors (called RsaL and RsaM) in containing the QS response within cost-effective levels; this is central for pseudomonads as they have very versatile genomes allowing them to live in constantly changing and highly dynamic environments. This review focuses on the role played by RsaL and RsaM repressors and discusses the important implications of this control of the QS response.

Classifying the topology of AHL-driven quorum sensing circuits in proteobacterial genomes.

Virulence and adaptability of many Gram-negative bacterial species are associated with an N-acylhomoserine lactone (AHL) gene regulation mechanism called quorum sensing (QS). The arrangement of quorum sensing genes is variable throughout bacterial genomes, although there are unifying themes that are common among the various topological arrangements. A bioinformatics survey of 1,403 complete bacterial genomes revealed characteristic gene topologies in 152 genomes that could be classified into 16 topological groups. We developed a concise notation for the patterns and show that the sequences of LuxR regulators and LuxI autoinducer synthase proteins cluster according to the topological patterns. The annotated topologies are deposited online with links to sequences and genome annotations at http://bacteria.itk.ppke.hu/QStopologies/.

Different degrees of structural order in distinct regions of the transcriptional repressor HES-1.

HES-1 is a transcriptional repressor of the basic helix-loop-helix (bHLH) family and one of the main downstream effectors in Notch signaling. Its domain architecture is composed of a bHLH region, an Orange domain, and a poorly characterized C-terminal half. We show that different degrees of structural order are present in the different regions of HES-1. The isolated bHLH domain is only marginally stable in solution, and partially folds upon dimerization. Binding to DNA promotes folding, stabilization, and protection from proteolysis of the bHLH domain. The Orange domain, on the contrary, is well folded in all conditions, forms stable dimers, and greatly increases protein resistance to thermal denaturation. The isolated proline-rich C-terminal region is mainly disordered in solution, and remains unstructured also in the full length protein. Measurements of binding constants show that HES-1 recognizes dsDNA synthetic oligonucleotides corresponding to several functional DNA targets with high affinity, but with relatively little specificity. We propose that order/disorder transitions in the different domains are associated not only with binding to DNA, but also with protein homo- and hetero-dimerization.

The quest for orthologs: finding the corresponding gene across genomes.

Orthology is a key evolutionary concept in many areas of genomic research. It provides a framework for subjects as diverse as the evolution of genomes, gene functions, cellular networks and functional genome annotation. Although orthologous proteins usually perform equivalent functions in different species, establishing true orthologous relationships requires a phylogenetic approach, which combines both trees and graphs (networks) using reliable species phylogeny and available genomic data from more than two species, and an insight into the processes of molecular evolution. Here, we evaluate the available bioinformatics tools and provide a set of guidelines to aid researchers in choosing the most appropriate tool for any situation.

The plant pathogen Pseudomonas fuscovaginae contains two conserved quorum sensing systems involved in virulence and negatively regulated by RsaL and the novel regulator RsaM.

Pseudomonas fuscovaginae is a Gram-negative fluorescent pseudomonad pathogenic towards several plant species. Despite its importance as a plant pathogen, no molecular studies of virulence have thus far been reported. In this study we show that P. fuscovaginae possesses two conserved N-acyl homoserine lactone (AHL) quorum sensing (QS) systems which we designated PfsI/R and PfvI/R. The PfsI/R system is homologous to the BviI/R system of Burkholderia vietnamiensis and produces and responds to C10-HSL and C12-HSL whereas PfvI/R is homologous to the LasI/R system of Pseudomonas aeruginosa and produces several long-chain 3-oxo-HSLs and responds to 3-oxo-C10-HSL and 3-oxo-C12-HSL and at high AHL concentrations can also respond to structurally different long-chain AHLs. Both systems were found to be negatively regulated by a repressor protein which was encoded by a gene located intergenically between the AHL synthase and LuxR-family response regulator. The pfsI/R system was regulated by a novel repressor designated RsaM while the pfvI/R system was regulated by both the RsaL repressor and by RsaM. The two systems are not transcriptionally hierarchically organized but share a common AHL response and both are required for plant virulence. Pseudomonas fuscovaginae has therefore a unique complex regulatory network composed of at least two different repressors which directly regulate the AHL QS systems and pathogenicity.

Prediction of new abiotic stress genes in Arabidopsis thaliana and Oryza sativa according to enumeration-based statistical analysis.

Plants undergo an extensive change in gene regulation during abiotic stress. It is of great agricultural importance to know which genes are affected during stress response. The genome sequence of a number of plant species has been determined, among them Arabidopsis and Oryza sativa, whose genome has been annotated most completely as of yet, and are well-known organisms widely used as experimental systems. This paper applies a statistical algorithm for predicting new stress-induced motifs and genes by analyzing promoter sets co-regulated by abiotic stress in the previously mentioned two species. After identifying characteristic putative regulatory motif sequence pairs (dyads) in the promoters of 125 stress-regulated Arabidopsis genes and 87 O. sativa genes, these dyads were used to screen the entire Arabidopsis and O. sativa promoteromes to find related stress-induced genes whose promoters contained a large number of these dyads found by our algorithm. We were able to predict a number of putative dyads, characteristic of a large number of stress-regulated genes, some of them newly discovered by our algorithm and serve as putative transcription factor binding sites. Our new motif prediction algorithm comes complete with a stand-alone program. This algorithm may be used in motif discovery in the future in other species. The more than 1,200 Arabidopsis and 1,700 Orzya sativa genes found by our algorithm are good candidates for further experimental studies in abiotic stress.

The interaction of Jagged-1 cytoplasmic tail with afadin PDZ domain is local, folding-independent, and tuned by phosphorylation.

Jagged-1, one of the five Notch ligands in man, is a membrane-spanning protein made of a large extracellular region and a 125-residue cytoplasmic tail bearing a C-terminal PDZ recognition motif ((1213) RMEYIV(1218) ). Binding of Jagged-1 intracellular region to the PDZ domain of afadin, a protein located at cell-cell adherens junctions, couples Notch signaling with the adhesion system and the cytoskeleton. Using NMR chemical shift perturbation and surface plasmon resonance, we studied the interaction between the PDZ domain of afadin (AF6_PDZ) and a series of polypeptides comprising the PDZ-binding motif. Chemical shift mapping of AF6_PDZ upon binding of ligands of different length (6, 24, and 133 residues) showed that the interaction is strictly local and involves only the binding groove in the PDZ. The recombinant protein corresponding to the entire intracellular region of Jagged-1, J1_ic, is mainly disordered in solution, and chemical shift mapping of J1_ic in the presence of AF6_PDZ showed that binding is not coupled to folding. Binding studies on a series of 24-residue peptides phosphorylated at different positions showed that phosphorylation of the tyrosine at position -2 of the PDZ-binding motif decreases its affinity for AF6_PDZ, and may play a role in the modulation of this interaction. Finally, we show that the R1213Q mutation located in the PDZ-binding motif and associated with extrahepatic biliary atresia increases the affinity for AF6_PDZ, suggesting that this syndrome may arise from an imbalance in the coupling of Notch signaling to the cytoskeleton.

Multi-netclust: an efficient tool for finding connected clusters in multi-parametric networks.

UNLABELLED: Multi-netclust is a simple tool that allows users to extract connected clusters of data represented by different networks given in the form of matrices. The tool uses user-defined threshold values to combine the matrices, and uses a straightforward, memory-efficient graph algorithm to find clusters that are connected in all or in either of the networks. The tool is written in C/C++ and is available either as a form-based or as a command-line-based program running on Linux platforms. The algorithm is fast, processing a network of > 10(6) nodes and 10(8) edges takes only a few minutes on an ordinary computer. AVAILABILITY: http://www.bioinformatics.nl/netclust/.

ProGMap: an integrated annotation resource for protein orthology.

Current protein sequence databases employ different classification schemes that often provide conflicting annotations, especially for poorly characterized proteins. ProGMap (Protein Group Mappings, http://www.bioinformatics.nl/progmap) is a web-tool designed to help researchers and database annotators to assess the coherence of protein groups defined in various databases and thereby facilitate the annotation of newly sequenced proteins. ProGMap is based on a non-redundant dataset of over 6.6 million protein sequences which is mapped to 240,000 protein group descriptions collected from UniProt, RefSeq, Ensembl, COG, KOG, OrthoMCL-DB, HomoloGene, TRIBES and PIRSF. ProGMap combines the underlying classification schemes via a network of links constructed by a fast and fully automated mapping approach originally developed for document classification. The web interface enables queries to be made using sequence identifiers, gene symbols, protein functions or amino acid and nucleotide sequences. For the latter query type BLAST similarity search and QuickMatch identity search services have been incorporated, for finding sequences similar (or identical) to a query sequence. ProGMap is meant to help users of high throughput methodologies who deal with partially annotated genomic data.