Bacterial cytochrome P450s: a bioinformatics odyssey of substrate discovery

Bacterial P450 cytochromes (BacCYPs) are versatile heme-containing proteins responsible for oxidation reactions on a wide range of substrates, contributing to the production of valuable natural products with limitless biotechnological potential. While the sequencing of microbial genomes has provided a wealth of BacCYP sequences, functional characterization lags behind, hindering our understanding of their roles. This study employs a comprehensive approach to predict BacCYP substrate specificity, bridging the gap between sequence and function. We employed an integrated approach combining sequence and functional data analysis, genomic context exploration, 3D structural modeling with molecular docking, and phylogenetic clustering. The research begins with an in-depth analysis of BacCYP sequence diversity and structural characteristics, revealing conserved motifs and recurrent residues in the active site. Phylogenetic analysis identifies distinct groups within the BacCYP family based on sequence similarity. However, our study reveals that sequence alone does not consistently predict substrate specificity, necessitating additional perspectives. The study delves into the genetic context of BacCYPs, utilizing neighboring gene information to infer potential substrates, a method proven very effective in many cases. Molecular docking is employed to assess BacCYP-substrate interactions, confirming potential substrates and providing insights into selectivity. Finally, a comprehensive strategy is proposed for predicting BacCYP substrates, involving all the evaluated approaches. The effectiveness of this strategy is demonstrated with two case studies, highlighting its potential for substrate discovery.

Spatial and Seasonal Variations in the Bacterial Community of an Anthropogenic Impacted Urban Stream

Environmental changes and human activities can alter the structure and diversity of aquatic microbial communities. In this work, we analyzed the bacterial community dynamics of an urban stream to understand how these factors affect the composition of river microbial communities. Samples were taken from a stream situated in Buenos Aires, Argentina, which flows through residential, peri-urban horticultural, and industrial areas. For sampling, two stations were selected: one influenced by a series of industrial waste treatment plants and horticultural farms (PL), and the other influenced by residential areas (R). Microbial communities were analyzed by sequence analysis of 16S rRNA gene amplicons along an annual cycle. PL samples showed high nutrient content compared with R samples. The diversity and richness of the R site were more affected by seasonality than those of the PL site. At the amplicon sequence variants level, beta diversity analysis showed a differentiation between cool-season (fall and winter) and warm-season (spring and summer) samples, as well as between PL and R sites. This demonstrated that there is spatial and temporal heterogeneity in the composition of the bacterial community, which should be considered if a bioremediation strategy is applied. The taxonomic composition analysis also revealed a differential seasonal cycle of phototrophs and chemoheterotrophs between the sampling sites, as well as different taxa associated with each sampling site. This analysis, combined with a comparative analysis of global rivers, allowed us to determine the genera Arcobacter, Simplicispira, Vogesella, and Sphingomonas as potential bioindicators of anthropogenic disturbance.

Assessing the hidden diversity underlying consensus sequences of SARS-CoV-2 using VICOS, a novel bioinformatic pipeline for identification of mixed viral populations

INTRODUCTION

Coinfection with two SARS-CoV-2 viruses is still a very understudied phenomenon. Although next generation sequencing methods are very sensitive to detect heterogeneous viral populations in a sample, there is no standardized method for their characterization, so their clinical and epidemiological importance is unknown.

MATERIAL AND METHODS

We developed VICOS (Viral COinfection Surveillance), a new bioinformatic algorithm for variant calling, filtering and statistical analysis to identify samples suspected of being mixed SARS-CoV-2 populations from a large dataset in the framework of a community genomic surveillance. VICOS was used to detect SARS-CoV-2 coinfections in a dataset of 1,097 complete genomes collected between March 2020 and August 2021 in Argentina.

RESULTS

We detected 23 cases (2%) of SARS-CoV-2 coinfections. Detailed study of VICOS's results together with additional phylogenetic analysis revealed 3 cases of coinfections by two viruses of the same lineage, 2 cases by viruses of different genetic lineages, 13 were compatible with both coinfection and intra-host evolution, and 5 cases were likely a product of laboratory contamination.

DISCUSSION

Intra-sample viral diversity provides important information to understand the transmission dynamics of SARS-CoV-2. Advanced bioinformatics tools, such as VICOS, are a necessary resource to help unveil the hidden diversity of SARS-CoV-2.

Selecting putative drought-tolerance markers in two contrasting soybeans

Identifying high-yield genotypes under low water availability is essential for soybean climate-smart breeding. However, a major bottleneck lies in phenotyping, particularly in selecting cost-efficient markers associated with stress tolerance and yield stabilization. Here, we conducted in-depth phenotyping experiments in two soybean genotypes with contrasting drought tolerance, MUNASQA (tolerant) and TJ2049 (susceptible), to better understand soybean stress physiology and identify/statistically validate drought-tolerance and yield-stabilization traits as potential breeding markers. Firstly, at the critical reproductive stage (R5), the molecular differences between the genotype's responses to mild water deficit were explored through massive analysis of cDNA ends (MACE)-transcriptomic and gene ontology. MUNASQA transcriptional profile, compared to TJ2049, revealed significant differences when responding to drought. Next, both genotypes were phenotyped under mild water deficit, imposed in vegetative (V3) and R5 stages, by evaluating 22 stress-response, growth, and water-use markers, which were subsequently correlated between phenological stages and with yield. Several markers showed high consistency, independent of the phenological stage, demonstrating the effectiveness of the phenotyping methodology and its possible use for early selection. Finally, these markers were classified and selected according to their cost-feasibility, statistical weight, and correlation with yield. Here, pubescence, stomatal density, and canopy temperature depression emerged as promising breeding markers for the early selection of drought-tolerant soybeans.

Pathway Driven Target Selection in Klebsiella pneumoniae: Insights Into Carbapenem Exposure

Carbapenem-resistant Klebsiella pneumoniae (CR-KP) represents an emerging threat to public health. CR-KP infections result in elevated morbidity and mortality. This fact, coupled with their global dissemination and increasingly limited number of therapeutic options, highlights the urgency of novel antimicrobials. Innovative strategies linking genome-wide interrogation with multi-layered metabolic data integration can accelerate the early steps of drug development, particularly target selection. Using the BioCyc ontology, we generated and manually refined a metabolic network for a CR-KP, K. pneumoniae Kp13. Converted into a reaction graph, we conducted topological-based analyses in this network to prioritize pathways exhibiting druggable features and fragile metabolic points likely exploitable to develop novel antimicrobials. Our results point to the aptness of previously recognized pathways, such as lipopolysaccharide and peptidoglycan synthesis, and casts light on the possibility of targeting less explored cellular functions. These functions include the production of lipoate, trehalose, glycine betaine, and flavin, as well as the salvaging of methionine. Energy metabolism pathways emerged as attractive targets in the context of carbapenem exposure, targeted either alone or in conjunction with current therapeutic options. These results prompt further experimental investigation aimed at controlling this highly relevant pathogen.

Cost-Effective Method to Perform SARS-CoV-2 Variant Surveillance: Detection of Alpha, Gamma, Lambda, Delta, Epsilon, and Zeta in Argentina

SARS-CoV-2 variants with concerning characteristics have emerged since the end of 2020. Surveillance of SARS-CoV-2 variants was performed on a total of 4,851 samples from the capital city and 10 provinces of Argentina, during 51 epidemiological weeks (EWs) that covered the end of the first wave and the ongoing second wave of the COVID-19 pandemic in the country (EW 44/2020 to EW 41/2021). The surveillance strategy was mainly based on Sanger sequencing of a Spike coding region that allows the identification of signature mutations associated with variants. In addition, whole-genome sequences were obtained from 637 samples. The main variants found were Gamma and Lambda, and to a lesser extent, Alpha, Zeta, and Epsilon, and more recently, Delta. Whereas, Gamma dominated in different regions of the country, both Gamma and Lambda prevailed in the most populated area, the metropolitan region of Buenos Aires. The lineages that circulated on the first wave were replaced by emergent variants in a term of a few weeks. At the end of the ongoing second wave, Delta began to be detected, replacing Gamma and Lambda. This scenario is consistent with the Latin American variant landscape, so far characterized by a concurrent increase in Delta circulation and a stabilization in the number of cases. The cost-effective surveillance protocol presented here allowed for a rapid response in a resource-limited setting, added information on the expansion of Lambda in South America, and contributed to the implementation of public health measures to control the disease spread in Argentina.

Identification of Potential Kinase Inhibitors within the PI3K/AKT Pathway of Leishmania Species

Leishmaniasis is a public health disease that requires the development of more effective treatments and the identification of novel molecular targets. Since blocking the PI3K/AKT pathway has been successfully studied as an effective anticancer strategy for decades, we examined whether the same approach would also be feasible in Leishmania due to their high amount and diverse set of annotated proteins. Here, we used a best reciprocal hits protocol to identify potential protein kinase homologues in an annotated human PI3K/AKT pathway. We calculated their ligandibility based on available bioactivity data of the reported homologues and modelled their 3D structures to estimate the druggability of their binding pockets. The models were used to run a virtual screening method with molecular docking. We found and studied five protein kinases in five different Leishmania species, which are AKT, CDK, AMPK, mTOR and GSK3 homologues from the studied pathways. The compounds found for different enzymes and species were analysed and suggested as starting point scaffolds for the design of inhibitors. We studied the kinases’ participation in protein–protein interaction networks, and the potential deleterious effects, if inhibited, were supported with the literature. In the case of Leishmania GSK3, an inhibitor of its human counterpart, prioritized by our method, was validated in vitro to test its anti-Leishmania activity and indirectly infer the presence of the enzyme in the parasite. The analysis contributes to improving the knowledge about the presence of similar signalling pathways in Leishmania, as well as the discovery of compounds acting against any of these kinases as potential molecular targets in the parasite.

From Genome to Drugs: New Approaches in Antimicrobial Discovery

Decades of successful use of antibiotics is currently challenged by the emergence of increasingly resistant bacterial strains. Novel drugs are urgently required but, in a scenario where private investment in the development of new antimicrobials is declining, efforts to combat drug-resistant infections become a worldwide public health problem. Reasons behind unsuccessful new antimicrobial development projects range from inadequate selection of the molecular targets to a lack of innovation. In this context, increasingly available omics data for multiple pathogens has created new drug discovery and development opportunities to fight infectious diseases. Identification of an appropriate molecular target is currently accepted as a critical step of the drug discovery process. Here, we review how diverse layers of multi-omics data in conjunction with structural/functional analysis and systems biology can be used to prioritize the best candidate proteins. Once the target is selected, virtual screening can be used as a robust methodology to explore molecular scaffolds that could act as inhibitors, guiding the development of new drug lead compounds. This review focuses on how the advent of omics and the development and application of bioinformatics strategies conduct a "big-data era" that improves target selection and lead compound identification in a cost-effective and shortened timeline.

Target-Pathogen: a structural bioinformatic approach to prioritize drug targets in pathogens

Available genomic data for pathogens has created new opportunities for drug discovery and development to fight them, including new resistant and multiresistant strains. In particular structural data must be integrated with both, gene information and experimental results. In this sense, there is a lack of an online resource that allows genome wide-based data consolidation from diverse sources together with thorough bioinformatic analysis that allows easy filtering and scoring for fast target selection for drug discovery. Here, we present Target-Pathogen database (http://target.sbg.qb.fcen.uba.ar/patho), designed and developed as an online resource that allows the integration and weighting of protein information such as: function, metabolic role, off-targeting, structural properties including druggability, essentiality and omic experiments, to facilitate the identification and prioritization of candidate drug targets in pathogens. We include in the database 10 genomes of some of the most relevant microorganisms for human health (Mycobacterium tuberculosis, Mycobacterium leprae, Klebsiella pneumoniae, Plasmodium vivax, Toxoplasma gondii, Leishmania major, Wolbachia bancrofti, Trypanosoma brucei, Shigella dysenteriae and Schistosoma Smanosoni) and show its applicability. New genomes can be uploaded upon request.

Core regulon of the global anaerobic regulator Anr targets central metabolism functions in Pseudomonas species

A comparative genome analysis of the global anaerobic regulator Anr regulon in five species of Pseudomonas with different life style was performed. Expression of this regulator was detected in all analyzed Pseudomonas. The predicted Anr regulon (pan-regulon) consisted of 253 genes. However, only 11 Anr-boxes located upstream of qor/hemF, hemN, cioA/PA3931, azu, rpsL, gltP, orthologous to PA2867, cspD, tyrZ, slyD and oprG, were common to all species. Whole genome in silico prediction of metabolic pathways identified genes belonging to heme biosynthesis, cytochromes and Entner-Doudoroff pathway as members of Anr regulon in all strains. Extending genome analysis to 28 Pseudomonas spp. spanning all phylogenetic groups showed Anr-boxes conservation in genes related to these functions. When present, genes related to anaerobic metabolism were predicted to hold Anr-boxes. Focused on the genomes of eight P. aeruginosa isolates of diverse origins, we observed a conserved regulon, sharing nearly 80% of the genes, indicating its key role in this opportunistic pathogen. The results suggest that the core Anr regulon comprises genes involved in central metabolism and aerobic electron transport chain, whereas those genes related to anaerobic metabolism and other functions constitute the accessory Anr-regulon, thereby differentially contributing to the ecological fitness of each Pseudomonas species.

An integrative, multi-omics approach towards the prioritization of Klebsiella pneumoniae drug targets

Klebsiella pneumoniae (Kp) is a globally disseminated opportunistic pathogen that can cause life-threatening infections. It has been found as the culprit of many infection outbreaks in hospital environments, being particularly aggressive towards newborns and adults under intensive care. Many Kp strains produce extended-spectrum β-lactamases, enzymes that promote resistance against antibiotics used to fight these infections. The presence of other resistance determinants leading to multidrug-resistance also limit therapeutic options, and the use of 'last-resort' drugs, such as polymyxins, is not uncommon. The global emergence and spread of resistant strains underline the need for novel antimicrobials against Kp and related bacterial pathogens. To tackle this great challenge, we generated multiple layers of 'omics' data related to Kp and prioritized proteins that could serve as attractive targets for antimicrobial development. Genomics, transcriptomics, structuromic and metabolic information were integrated in order to prioritize candidate targets, and this data compendium is freely available as a web server. Twenty-nine proteins with desirable characteristics from a drug development perspective were shortlisted, which participate in important processes such as lipid synthesis, cofactor production, and core metabolism. Collectively, our results point towards novel targets for the control of Kp and related bacterial pathogens.

Achieving maximal speed of solution exchange for patch clamp experiments

Background

Resolving the kinetics of agonist binding events separately from the subsequent channel gating processes requires the ability of applying and removing the agonist before channel gating occurs. No reported system has yet achieved pulses shorter than 100 µs, necessary to study nicotinic ACh receptor or AMPA receptor activation.

Methodology/Principal Findings

Solution exchange systems deliver short agonist pulses by moving a sharp interface between a control and an experimental solution across a channel preparation. We achieved shorter pulses by means of an exchange system that combines a faster flow velocity, narrower partition between the two streams, and increased velocity and bandwidth of the movement of the interface. The measured response of the entire system was fed back to optimize the voltage signal applied to the piezoelectric actuator overcoming the spurious oscillations arising from the mechanical resonances when a high bandwidth driving function was applied. Optimization was accomplished by analyzing the transfer function of the solution exchange system. When driven by optimized command pulses the enhanced system provided pulses lasting 26 ± 1 µs and exchanging 93 ± 1% of the solution, as measured in the open tip of a patch pipette.

Conclusions/Significance

Pulses of this duration open the experimental study of the molecular events that occur between the agonist binding and the opening of the channel.

Leptin in embryos from control and diabetic rats during organogenesis: a modulator of nitric oxide production and lipid homeostasis

BACKGROUND

Leptin is involved in many metabolic and reproductive events and its levels are altered by the diabetic pathology. In this study, leptin concentrations and leptin effects on both nitric oxide (NO) and lipid concentrations were investigated in embryos from control and diabetic rats.

METHODS

Diabetes was induced by neonatal streptozotocin administration (90 mg/kg). Embryos from control and diabetic rats were obtained on days 10.5 and 13.5 of gestation, corresponding to early organogenesis and post-placentation periods respectively. Leptin was analysed by enzyme immunoanalysis and immunohistochemistry. Nitrates and nitrites were assessed as an index of NO production. Lipid concentrations were analysed by thin layer chromatography.

RESULTS

Leptin concentrations were decreased in embryos obtained from diabetic rats on days 10.5 and 13.5 of gestation when compared to controls. NO concentrations, elevated in diabetic embryopathy, were diminished in the presence of leptin in the embryos obtained from control and diabetic animals both during early organogenesis and after placentation. Leptin additions reduced phospholipid, cholesterol and cholesteryl ester concentrations in embryos obtained from diabetic rats during early organogenesis, although no leptin effects on lipid concentrations were observed in control embryos at this developmental stage. In embryos obtained on day 13.5 of gestation leptin additions reduced cholesteryl ester concentrations in controls, and diminished cholesteryl ester, triglycerides and phospholipids in embryos from diabetic rats.

CONCLUSIONS

We demonstrated that leptin plays a role in the regulation of NO concentrations and lipid homeostasis during embryo organogenesis and that the diabetic environment causes a reduction of leptin concentrations in rat embryos.