Small RNAs direct attack and defense mechanisms in a quorum sensing phage and its host.

Many, if not all, bacteria use quorum sensing (QS) to control collective behaviors, and more recently, QS has also been discovered in bacteriophages (phages). Phages can produce communication molecules of their own, or "listen in" on the host's communication processes, to switch between lytic and lysogenic modes of infection. Here, we study the interaction of Vibrio cholerae with the lysogenic phage VP882, which is activated by the QS molecule DPO. We discover that induction of VP882 results in the binding of phage transcripts to the major RNA chaperone Hfq, which in turn outcompetes and downregulates host-encoded small RNAs (sRNAs). VP882 itself also encodes Hfq-binding sRNAs, and we demonstrate that one of these sRNAs, named VpdS, promotes phage replication by regulating host and phage mRNA levels. We further show that host-encoded sRNAs can antagonize phage replication by downregulating phage mRNA expression and thus might be part of the host's phage defense arsenal.

ChimericFragments: computation, analysis and visualization of global RNA networks.

RNA-RNA interactions are a key feature of post-transcriptional gene regulation in all domains of life. While ever more experimental protocols are being developed to study RNA duplex formation on a genome-wide scale, computational methods for the analysis and interpretation of the underlying data are lagging behind. Here, we present ChimericFragments, an analysis framework for RNA-seq experiments that produce chimeric RNA molecules. ChimericFragments implements a novel statistical method based on the complementarity of the base-pairing RNAs around their ligation site and provides an interactive graph-based visualization for data exploration and interpretation. ChimericFragments detects true RNA-RNA interactions with high precision and is compatible with several widely used experimental procedures such as RIL-seq, LIGR-seq or CLASH. We further demonstrate that ChimericFragments enables the systematic detection of novel RNA regulators and RNA-target pairs with crucial roles in microbial physiology and virulence. ChimericFragments is written in Julia and available at: https://github.com/maltesie/ChimericFragments.

The global RNA-RNA interactome of Klebsiella pneumoniae unveils a small RNA regulator of cell division.

The ubiquitous RNA chaperone Hfq is involved in the regulation of key biological processes in many species across the bacterial kingdom. In the opportunistic human pathogen Klebsiella pneumoniae, deletion of the hfq gene affects the global transcriptome, virulence, and stress resistance; however, the ligands of the major RNA-binding protein in this species have remained elusive. In this study, we have combined transcriptomic, co-immunoprecipitation, and global RNA interactome analyses to compile an inventory of conserved and species-specific RNAs bound by Hfq and to monitor Hfq-mediated RNA-RNA interactions. In addition to dozens of RNA-RNA pairs, our study revealed an Hfq-dependent small regulatory RNA (sRNA), DinR, which is processed from the 3' terminal portion of dinI mRNA. Transcription of dinI is controlled by the master regulator of the SOS response, LexA. As DinR accumulates in K. pneumoniae in response to DNA damage, the sRNA represses translation of the ftsZ transcript by occupation of the ribosome binding site. Ectopic overexpression of DinR causes depletion of ftsZ mRNA and inhibition of cell division, while deletion of dinR antagonizes cell elongation in the presence of DNA damage. Collectively, our work highlights the important role of RNA-based gene regulation in K. pneumoniae and uncovers the central role of DinR in LexA-controlled division inhibition during the SOS response.

An RNA sponge controls quorum sensing dynamics and biofilm formation in Vibrio cholerae.

Small regulatory RNAs (sRNAs) acting in concert with the RNA chaperone Hfq are prevalent in many bacteria and typically act by base-pairing with multiple target transcripts. In the human pathogen Vibrio cholerae, sRNAs play roles in various processes including antibiotic tolerance, competence, and quorum sensing (QS). Here, we use RIL-seq (RNA-interaction-by-ligation-and-sequencing) to identify Hfq-interacting sRNAs and their targets in V. cholerae. We find hundreds of sRNA-mRNA interactions, as well as RNA duplexes formed between two sRNA regulators. Further analysis of these duplexes identifies an RNA sponge, termed QrrX, that base-pairs with and inactivates the Qrr1-4 sRNAs, which are known to modulate the QS pathway. Transcription of qrrX is activated by QrrT, a previously uncharacterized LysR-type transcriptional regulator. Our results indicate that QrrX and QrrT are required for rapid conversion from individual to community behaviours in V. cholerae.

Conserved hydrogen-bond motifs of membrane transporters and receptors.

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

Interactive Interface for Graph-Based Analyses of Dynamic H-Bond Networks: Application to Spike Protein S.

Dynamic hydrogen-bond networks are key determinants of protein conformational dynamics. In the case of macromolecular protein complexes, which can have a large number of hydrogen bonds giving rise to extensive hydrogen-bond networks, efficient algorithms are required to analyze interactions that could be important for the dynamics and biological function of the complex. We present here a highly efficient, standalone interface designed for analyses of dynamical hydrogen-bond networks of biomolecules and macromolecular complexes. To facilitate a comprehensive description of protein dynamics, the interface includes analyses of hydrophobic interactions. We illustrate the usefulness and workflow of the interface by dissecting the dynamics of the ectodomain of SARS-CoV-2 protein S in its closed conformation. We find that protein S contains numerous local clusters of dynamic hydrogen bonds and identify hydrogen bonds that are sampled persistently. The receptor binding domain of the spike protein hosts only a handful of persistent hydrogen-bond clusters, suggesting structural plasticity. Our data analysis interface is released here for open use.

A graph-based approach identifies dynamic H-bond communication networks in spike protein S of SARS-CoV-2.

Corona virus spike protein S is a large homo-trimeric protein anchored in the membrane of the virion particle. Protein S binds to angiotensin-converting-enzyme 2, ACE2, of the host cell, followed by proteolysis of the spike protein, drastic protein conformational change with exposure of the fusion peptide of the virus, and entry of the virion into the host cell. The structural elements that govern conformational plasticity of the spike protein are largely unknown. Here, we present a methodology that relies upon graph and centrality analyses, augmented by bioinformatics, to identify and characterize large H-bond clusters in protein structures. We apply this methodology to protein S ectodomain and find that, in the closed conformation, the three protomers of protein S bring the same contribution to an extensive central network of H-bonds, and contribute symmetrically to a relatively large H-bond cluster at the receptor binding domain, and to a cluster near a protease cleavage site. Markedly different H-bonding at these three clusters in open and pre-fusion conformations suggest dynamic H-bond clusters could facilitate structural plasticity and selection of a protein S protomer for binding to the host receptor, and proteolytic cleavage. From analyses of spike protein sequences we identify patches of histidine and carboxylate groups that could be involved in transient proton binding.

Mechanism of Inward Proton Transport in an Antarctic Microbial Rhodopsin.

Although the outward-directed proton transport across biological membranes is well studied and its importance for bioenergetics is clearly understood, inward-directed light-driven proton pumping by microbial rhodopsins has remained a mystery both physiologically and mechanistically. A new family of Antarctic rhodopsins, which is a subgroup within a novel class of schizorhodopsins reported recently, includes a member, denoted as AntR, which proved amenable to extensive characterization with experiments and computation. Phylogenetic analyses identify AntR as distinct from the well-studied microbial rhodopsins that function as outward-directed ion pumps, and bioinformatics sequence analyses reveal amino acid substitutions at conserved sites essential for outward proton pumping. Modeling and numerical simulations of AntR, combined with advanced analyses using the graph theory and centrality measures from social sciences, identify the dynamic three-dimensional network of hydrogen-bonded water molecules and amino acid residues that function as communication hubs in AntR. This network undergoes major rearrangement upon retinal isomerization, showing important changes in the connectivity of the active center, retinal Schiff base, to the opposing sides of the membrane, as required for proton transport. Numerical simulations and experimental studies of the photochemical cycle of AntR by spectroscopy and site-directed mutagenesis allowed us to identify pathways that could conduct protons in the direction opposite to that commonly known for outward-directed pumps.

Bridge: A Graph-Based Algorithm to Analyze Dynamic H-Bond Networks in Membrane Proteins.

Membrane proteins that function as transporters or receptors must communicate with both sides of the lipid bilayer in which they sit. This long distance communication enables transporters to move protons or other ions and small molecules across the bilayer and receptors to transmit an external signal to the cell. Hydrogen bonds, hydrogen-bond networks, and lipid-protein interactions are essential for the motions and functioning of the membrane protein and, consequently, of outmost interest to structural biology and numerical simulations. We present here Bridge, an algorithm tailored for efficient analyses of hydrogen-bond networks in membrane transporter and receptor proteins. For channelrhodopsin, a membrane protein whose functioning involves proton-transfer reactions, Bridge identifies extensive networks of protein-water hydrogen bonds and an unanticipated network that can bridge transiently two proton donors across a distance of ∼20 Å. Graphs of the protein hydrogen bonds reveal rapid propagation of structural changes within hydrogen-bond networks of mutant transporters and identify protein groups potentially important for the proton transfer activity. The algorithm is made available as a plugin for PyMol.

Dynamics of Long-Distance Hydrogen-Bond Networks in Photosystem II.

Photosystem II uses the energy of absorbed light to split water molecules, generating molecular oxygen, electrons, and protons. The four protons generated during each reaction cycle are released to the lumen via mechanisms that are poorly understood. Given the complexity of photosystem II, which consists of multiple protein subunits and cofactor molecules and hosts numerous waters, a fundamental issue is finding transient networks of hydrogen bonds that bridge potential proton donor and acceptor groups. Here, we address this issue by performing all-atom molecular dynamics simulations of wild-type and mutant photosystem II monomers, which we analyze using a new protocol designed to facilitate efficient analysis of hydrogen-bond networks. Our computations reveal that local protein/water hydrogen-bond networks can assemble transiently in photosystem II such that the reaction center connects to the lumen. The dynamics of the hydrogen-bond networks couple to the protonation state of specific carboxylate groups and are altered in a mutant with defective proton transfer. Simulations on photosystem II without its extrinsic PsbO subunit provide a molecular interpretation of the elusive functional role of this subunit.