Synergistic Impacts of Organic Acids and pH on Growth of Pseudomonas aeruginosa: A Comparison of Parametric and Bayesian Non-parametric Methods to Model Growth

Different weak organic acids have significant potential as topical treatments for wounds infected by opportunistic pathogens that are recalcitrant to standard treatments. These acids have long been used as bacteriostatic compounds in the food industry, and in some cases are already being used in the clinic. The effects of different organic acids vary with pH, concentration, and the specific organic acid used, but no studies to date on any opportunistic pathogens have examined the detailed interactions between these key variables in a controlled and systematic way. We have therefore comprehensively evaluated the effects of several different weak organic acids on growth of the opportunistic pathogen Pseudomonas aeruginosa. We used a semi-automated plate reader to generate growth profiles for two different strains (model laboratory strain PAO1 and clinical isolate PA1054 from a hospital burns unit) in a range of organic acids at different concentrations and pH, with a high level of replication for a total of 162,960 data points. We then compared two different modeling approaches for the interpretation of this time-resolved dataset: parametric logistic regression (with or without a component to include lag phase) vs. non-parametric Gaussian process (GP) regression. Because GP makes no prior assumptions about the nature of the growth, this method proved to be superior in cases where growth did not follow a standard sigmoid functional form, as is common when bacteria grow under stress. Acetic, propionic and butyric acids were all more detrimental to growth than the other acids tested, and although PA1054 grew better than PAO1 under non-stress conditions, this difference largely disappeared as the levels of stress increased. As expected from knowledge of how organic acids behave, their effect was significantly enhanced in combination with low pH, with this interaction being greatest in the case of propionic acid. Our approach lends itself to the characterization of combinatorial interactions between stressors, especially in cases where their impacts on growth render logistic growth models unsuitable

Using microbiological data to improve the use of antibiotics for respiratory tract infections: A protocol for an individual patient data meta-analysis

Background: Resistance to antibiotics is rising and threatens future antibiotic effectiveness. ‘Antibiotic targeting’ ensures patients who may benefit from antibiotics receive them, while being safely withheld from those who may not. Point-of-care tests may assist with antibiotic targeting by allowing primary care clinicians to establish if symptomatic patients have a viral, bacterial, combined, or no infection. However, because organisms can be harmlessly carried, it is important to know if the presence of the virus/bacteria is related to the illness for which the patient is being assessed. One way to do this is to look for associations with more severe/prolonged symptoms and test results. Previous research to answer this question for acute respiratory tract infections has given conflicting results with studies has not having enough participants to provide statistical confidence. Aim: To undertake a synthesis of IPD from both randomised controlled trials (RCTs) and observational cohort studies of respiratory tract infections (RTI) in order to investigate the prognostic value of microbiological data in addition to, or instead of, clinical symptoms and signs. Methods: A systematic search of Cochrane Central Register of Controlled Trials, Ovid Medline and Ovid Embase will be carried out for studies of acute respiratory infection in primary care settings. The outcomes of interest are duration of disease, severity of disease, repeated consultation with new/worsening illness and complications requiring hospitalisation. Authors of eligible studies will be contacted to provide anonymised individual participant data. The data will be harmonised and aggregated. Multilevel regression analysis will be conducted to determine key outcome measures for different potential pathogens and whether these offer any additional information on prognosis beyond clinical symptoms and signs. Trial registration: PROSPERO Registration number: CRD42023376769

Using microbiological data to improve the use of antibiotics for respiratory tract infections: A protocol for an individual patient data meta-analysis

Background Resistance to antibiotics is rising and threatens future antibiotic effectiveness. ‘Antibiotic targeting’ ensures patients who may benefit from antibiotics receive them, while being safely withheld from those who may not. Point-of-care tests may assist with antibiotic targeting by allowing primary care clinicians to establish if symptomatic patients have a viral, bacterial, combined, or no infection. However, because organisms can be harmlessly carried, it is important to know if the presence of the virus/bacteria is related to the illness for which the patient is being assessed. One way to do this is to look for associations with more severe/prolonged symptoms and test results. Previous research to answer this question for acute respiratory tract infections has given conflicting results with studies has not having enough participants to provide statistical confidence. Aim To undertake a synthesis of IPD from both randomised controlled trials (RCTs) and observational cohort studies of respiratory tract infections (RTI) in order to investigate the prognostic value of microbiological data in addition to, or instead of, clinical symptoms and signs. Methods A systematic search of Cochrane Central Register of Controlled Trials, Ovid Medline and Ovid Embase will be carried out for studies of acute respiratory infection in primary care settings. The outcomes of interest are duration of disease, severity of disease, repeated consultation with new/worsening illness and complications requiring hospitalisation. Authors of eligible studies will be contacted to provide anonymised individual participant data. The data will be harmonised and aggregated. Multilevel regression analysis will be conducted to determine key outcome measures for different potential pathogens and whether these offer any additional information on prognosis beyond clinical symptoms and signs. Trial registration PROSPERO Registration number: CRD42023376769

Detecting differential growth of microbial populations with Gaussian process regression

Microbial growth curves are used to study differential effects of media, genetics, and stress on microbial population growth. Consequently, many modeling frameworks exist to capture microbial population growth measurements. However, current models are designed to quantify growth under conditions for which growth has a specific functional form. Extensions to these models are required to quantify the effects of perturbations, which often exhibit nonstandard growth curves. Rather than assume specific functional forms for experimental perturbations, we developed a general and robust model of microbial population growth curves using Gaussian process (GP) regression. GP regression modeling of high-resolution time-series growth data enables accurate quantification of population growth and allows explicit control of effects from other covariates such as genetic background. This framework substantially outperforms commonly used microbial population growth models, particularly when modeling growth data from environmentally stressed populations. We apply the GP growth model and develop statistical tests to quantify the differential effects of environmental perturbations on microbial growth across a large compendium of genotypes in archaea and yeast. This method accurately identifies known transcriptional regulators and implicates novel regulators of growth under standard and stress conditions in the model archaeal organism Halobacterium salinarum. For yeast, our method correctly identifies known phenotypes for a diversity of genetic backgrounds under cyclohexamide stress and also detects previously unidentified oxidative stress sensitivity across a subset of strains. Together, these results demonstrate that the GP models are interpretable, recapitulating biological knowledge of growth response while providing new insights into the relevant parameters affecting microbial population growth

Predicted coronavirus Nsp5 protease cleavage sites in the human proteome

BACKGROUND: The coronavirus nonstructural protein 5 (Nsp5) is a cysteine protease required for processing the viral polyprotein and is therefore crucial for viral replication. Nsp5 from several coronaviruses have also been found to cleave host proteins, disrupting molecular pathways involved in innate immunity. Nsp5 from the recently emerged SARS-CoV-2 virus interacts with and can cleave human proteins, which may be relevant to the pathogenesis of COVID-19. Based on the continuing global pandemic, and emerging understanding of coronavirus Nsp5-human protein interactions, we set out to predict what human proteins are cleaved by the coronavirus Nsp5 protease using a bioinformatics approach. RESULTS: Using a previously developed neural network trained on coronavirus Nsp5 cleavage sites (NetCorona), we made predictions of Nsp5 cleavage sites in all human proteins. Structures of human proteins in the Protein Data Bank containing a predicted Nsp5 cleavage site were then examined, generating a list of 92 human proteins with a highly predicted and accessible cleavage site. Of those, 48 are expected to be found in the same cellular compartment as Nsp5. Analysis of this targeted list of proteins revealed molecular pathways susceptible to Nsp5 cleavage and therefore relevant to coronavirus infection, including pathways involved in mRNA processing, cytokine response, cytoskeleton organization, and apoptosis. CONCLUSIONS: This study combines predictions of Nsp5 cleavage sites in human proteins with protein structure information and protein network analysis. We predicted cleavage sites in proteins recently shown to be cleaved in vitro by SARS-CoV-2 Nsp5, and we discuss how other potentially cleaved proteins may be relevant to coronavirus mediated immune dysregulation. The data presented here will assist in the design of more targeted experiments, to determine the role of coronavirus Nsp5 cleavage of host proteins, which is relevant to understanding the molecular pathology of coronavirus infection. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12863-022-01044-y

A Bayesian non-parametric mixed-effects model of microbial growth curves

Substantive changes in gene expression, metabolism, and the proteome are manifested in overall changes in microbial population growth. Quantifying how microbes grow is therefore fundamental to areas such as genetics, bioengineering, and food safety. Traditional parametric growth curve models capture the population growth behavior through a set of summarizing parameters. However, estimation of these parameters from data is confounded by random effects such as experimental variability, batch effects or differences in experimental material. A systematic statistical method to identify and correct for such confounding effects in population growth data is not currently available. Further, our previous work has demonstrated that parametric models are insufficient to explain and predict microbial response under non-standard growth conditions. Here we develop a hierarchical Bayesian non-parametric model of population growth that identifies the latent growth behavior and response to perturbation, while simultaneously correcting for random effects in the data. This model enables more accurate estimates of the biological effect of interest, while better accounting for the uncertainty due to technical variation. Additionally, modeling hierarchical variation provides estimates of the relative impact of various confounding effects on measured population growth

A regulatory hierarchy controls the dynamic transcriptional response to extreme oxidative stress in archaea.

Networks of interacting transcription factors are central to the regulation of cellular responses to abiotic stress. Although the architecture of many such networks has been mapped, their dynamic function remains unclear. Here we address this challenge in archaea, microorganisms possessing transcription factors that resemble those of both eukaryotes and bacteria. Using genome-wide DNA binding location analysis integrated with gene expression and cell physiological data, we demonstrate that a bacterial-type transcription factor (TF), called RosR, and five TFIIB proteins, homologs of eukaryotic TFs, combinatorially regulate over 100 target genes important for the response to extremely high levels of peroxide. These genes include 20 other transcription factors and oxidative damage repair genes. RosR promoter occupancy is surprisingly dynamic, with the pattern of target gene expression during the transition from rapid growth to stress correlating strongly with the pattern of dynamic binding. We conclude that a hierarchical regulatory network orchestrated by TFs of hybrid lineage enables dynamic response and survival under extreme stress in archaea. This raises questions regarding the evolutionary trajectory of gene networks in response to stress

The genotype‐phenotype landscape of an allosteric protein

Abstract Allostery is a fundamental biophysical mechanism that underlies cellular sensing, signaling, and metabolism. Yet a quantitative understanding of allosteric genotype‐phenotype relationships remains elusive. Here, we report the large‐scale measurement of the genotype‐phenotype landscape for an allosteric protein: the lac repressor from Escherichia coli, LacI. Using a method that combines long‐read and short‐read DNA sequencing, we quantitatively measure the dose‐response curves for nearly 105 variants of the LacI genetic sensor. The resulting data provide a quantitative map of the effect of amino acid substitutions on LacI allostery and reveal systematic sequence‐structure‐function relationships. We find that in many cases, allosteric phenotypes can be quantitatively predicted with additive or neural‐network models, but unpredictable changes also occur. For example, we were surprised to discover a new band‐stop phenotype that challenges conventional models of allostery and that emerges from combinations of nearly silent amino acid substitutions