Exploration of Streptococcus core genome to reveal druggable targets and novel therapeutics against S. pneumoniae

Abstract

Streptococcus pneumoniae (S. pneumoniae), the major etiological agent of community-acquired pneumonia (CAP) contributes significantly to the global burden of infectious diseases which is getting resistant day by day. Nearly 30% of the S. pneumoniae genomes encode hypothetical proteins (HPs), and better understandings of these HPs in virulence and pathogenicity plausibly decipher new treatments. Some of the HPs are present across many Streptococcus species, systematic assessment of these unexplored HPs will disclose prospective drug targets. In this study, through a stringent bioinformatics analysis of the core genome and proteome of S. pneumoniae PCS8235, we identified and analyzed 28 HPs that are common in many Streptococcus species and might have a potential role in the virulence or pathogenesis of the bacteria. Functional annotations of the proteins were conducted based on the physicochemical properties, subcellular localization, virulence prediction, protein-protein interactions, and identification of essential genes, to find potentially druggable proteins among 28 HPs. The majority of the HPs are involved in bacterial transcription and translation. Besides, some of them were homologs of enzymes, binding proteins, transporters, and regulators. Protein-protein interactions revealed HP PCS8235_RS05845 made the highest interactions with other HPs and also has TRP structural motif along with virulent and pathogenic properties indicating it has critical cellular functions and might go under unconventional protein secretions. The second highest interacting protein HP PCS8235_RS02595 interacts with the Regulator of chromosomal segregation (RocS) which participates in chromosome segregation and nucleoid protection in S. pneumoniae. In this interacting network, 54% of protein members have virulent properties and 40% contain pathogenic properties. Among them, most of these proteins circulate in the cytoplasmic area and have hydrophilic properties. Finally, molecular docking and dynamics simulation demonstrated that the antimalarial drug Artenimol can act as a drug repurposing candidate against HP PCS8235_RS 04650 of S. pneumoniae. Hence, the present study could aid in drugs against S. pneumoniae.

Introduction

Community-acquired pneumonia (CAP) is one of the prime causes of death from infectious diseases worldwide [1]. It is an acute lung infection caused by a variety of microorganisms that cause symptoms such as shortness of breath, coughing, heavy sputum, fever, chills, and chest pain [2]. Currently, lower respiratory tract infection (LRTI) is considered to be the fourth-largest cause of death worldwide (https://www.who.int/news-room/fact-sheets/detail/the-top-10-causes-of-death), and CAP is still the leading cause of death among all infectious diseases in the United States of America [1]. Earlier studies have shown that about 900 000 children <5 years died from pneumococcal diseases (https://www.who.int/en/news-room/fact-sheets/detail/pneumonia) and 2.2–50.9% of cases belonged to pediatric CAP [3].

S. pneumoniae also known as pneumococcus is identified as the most recognizable causative agent of CAP in individuals with a compromised immune system [4]. In 2016, LRTI by S. pneumoniae engendered 1, 517, 388 deaths globally [5]. Between 1918 and 1919, pneumococcus played a dominant role in the global influenza pandemic [6]. Misuse of pneumonia-related antibiotics during the COVID-19 pandemic has the potential to generate more multi-drug resistant S. pneumoniae [7]. The fatality of CAP depends on the age and the existence of the comorbidities. The infection is more common among adults aged over 65 years and children under 2 years, or individuals who smoke, have asthma, or have Chronic obstructive pulmonary disease (COPD) [8, 9]. Apart from CAP, S. pneumoniae also causes bacterial meningitis, bacteremia, sinusitis, otitis media, septic arthritis, aortitis, gingival lesions, phlegmonous gastritis, inguinal adenitis, testicular and tubo-ovarian abscesses, and necrotizing fasciitis [10].

Currently, the therapeutics offered for patients with pneumonia are antibiotic therapies and vaccines. The application of these treatments varies from place to place and also on the severity of the disease. Usually, if the disease is not treated at the right time or is caused by virulent/ resistant strains; it leads to septic shock, empyema, parapneumonic effusion, lung abscess, necrotizing pneumonia, and often death [11–13]. There are two types of vaccines in the global market: the pneumococcal polysaccharide vaccine (23-valent pneumococcal polysaccharide vaccine (PPSV23), and the pneumococcal conjugate vaccine (10-valent PCV10 and 13-valent PCV13) [14, 15]. Owing to its poor immunogenicity in infants, PPSV23 is only recommended for adults and immunocompromised patients while PCV 13 is advised for infants <2 years old. Although PPSV23 has been widely used for over 40 years, there are some limitations. The antibody responses produced by this vaccine are T-cell independent and it was observed that the level of serotype-specific IgG and Opsonophagocytic killing Assay (OPA) lessens over time, hence revaccination is required for protection at an interval of 5 years [16]. Moreover, PPSV23 is less effective in men compared to women [17]. Thus, the efficacy of this vaccine is still unclear. Also, the existing vaccines had a major effect on CAP etiology giving rise to new S. pneumoniae serotypes. Currently, 98 serotypes with different polysaccharide capsules have been acknowledged [18]. Notably, capsule polysaccharide-based vaccines are limited to a certain number of serotypes.

Over the last few decades, a significant rise in antimicrobial resistance was observed in S. pneumoniae [19, 20]. This rise is associated with the emergence of multiple serotypes of the bacteria in various parts of the world [21]. Pneumococcus can undergo recombination-mediated genetic plasticity which allows them to acquire antibiotic resistance genes from several closely related species like Streptococcus mitis and Streptococcus oralis [22]. This in turn helps them evade vaccines and antimicrobial agents and evolve into new vaccine-escape mutants and high antimicrobial-resistant strains [23, 24]. The standard antimicrobial agent approved for CAP is a combination of β-lactam plus a macrolide or fluoroquinolone. β-lactam works by binding to the penicillin-binding proteins (PBP) of the bacteria and inhibits cell wall synthesis. Nevertheless, mutations in the mosaic genes encoding Penicillin-Binding Protein (PBP) results in resistant isolates [25]. Macrolide resistance and resistance to fluoroquinolones originate from the overuse of these broad-spectrum antibiotics. The level of macrolide resistance varies from region to region. North America and the United Kingdom are prone to low macrolide resistance by drug efflux whereas Asian countries are liable to high macrolide resistance as a result of ribosomal methylation [26]. Generally, fluoroquinolone is the only antibiotic used to target the DNA gyrase of S. pneumoniae directly bringing a halt to its protein synthesis. However, repeated use of this antimicrobial agent resulted in a spontaneous number of mutations in the chromosomal genes that encode this enzyme [27]. As a result, the growing resistance of S. pneumoniae to commonly used antibiotics and non-vaccine serotypes underlines the urgent need for a new therapeutic target.

Examination of the whole genome of the organism is an important approach in combating the regulation of antibiotic resistance, non-vaccine serotypes, and developments of therapeutics. When varieties of genus Streptococcus strains are aligned together, several core genes can be observed. Around 33% of the S. pneumoniae genome constitutes uncharacterized proteins documented as Hypothetical proteins (HPs) [28]. These proteins are encoded by computationally predicted open reading frames but lack biochemical and chemical evidence. Although they lack functional characterization, they play a crucial role in biochemical and physiological pathways [22]. Bioinformatics tools and algorithms are very efficient to explore proteins [29, 30]. Previous studies have shown some of the Streptococcus mutans’ HPs are critical for antibiotic resistance and biofilm formation [22, 31]. Additionally, this infectious agent produces several virulence factors involved in the survival of the pathogen and the progression of the disease [28]. Therefore, in the present study, we aimed to characterize the HPs encoded by the core genome of S. pneumoniae using a computational approach to uncover novel targets for drug development. Since the HPs are mutual in many Streptococcus species they might be interesting targets.

Methodology

Retrieval of the genome sequences

The core genome and proteome of Streptococcus strains from different Streptococci species were extracted using the Efficient Database framework for comparative Genome Analyses using BLAST score Ratios (EDGAR 3.0) and visualized by BioCircos [32]. EDGAR 3.0 is software designed to perform genome comparisons using a high throughput approach [33]. S. pneumoniae PCS8235 (NCBI Reference Sequence: NZ_CM001835.1) was taken as the reference strain. HPs were mined manually from the proteome datasets of the strains.

Functional enrichment and determination of physicochemical properties of the hypothetical proteins

The functions of the HPs were unveiled using the Gene Ontology Functional Enrichment Annotation Tool (GO FEAT) webserver which works through sequence homology search. The proteins were classified based on the conservation of domains, motifs, families, and superfamilies and categorized via the InterPro, UniProt, European Molecular Biology Laboratory (EMBL), Kyoto Encyclopedia of Genes and Genomes (KEGG), and the National Center for Biotechnology Information (NCBI) databases respectively [34]. The physicochemical properties of the annotated HPs were further documented using the ProtParam tool of Expasy (https://web.expasy.org/protparam/). The parameters of the proteins included the molecular weight, theoretical pI point, amino acid composition, extinction coefficient, instability index, aliphatic index, and grand average of hydropathicity (GRAVY) of the protein [35].

Analysis of subcellular localization and unconventional protein secretion

CELLO (http://cello.life.nctu.edu.tw/) and PSORTb 3.0 (https://www.psort.org/psortb/) were utilized to predict the subcellular localization of the HPs [36, 37]. Understanding subcellular localization is very important to characterize a protein as a target for a drug or vaccine [38]. The OutCyte 1.0 (http://www.outcyte.com/) and SecretomeP 2.0 (http://www.cbs.dtu.dk/services/SecretomeP/) were also used to reveal the unconventional protein secretion as well as the HPs taking non-classical secretory pathways respectively [39]. OutCyte 1.0 is an online bioinformatics tool that mediates two steps to finally generate the proteins without N-terminal signals [40].

Identification of essential proteins

All the HPs of S. pneumoniae were queried against the Database of Essential Genes (DEG) database (http://origin.tubic.org/deg/public/index.php) in search of homologous genes within the inquired sequence [41]. Only the proteins encoded by the query sequences which shared similarities with the essential genes in the DEG database and an E-value of <0.0001 and bit-score >100 as cutoff were designated as essential proteins.

Virulence properties and pathogenicity

VirulentPred (http://bioinfo.icgeb.res.in/virulent/) server was used to evaluate the virulence activity of the HPs. This online tool is an SVM-based method that calculates a virulence potential score for a given protein [42]. Simultaneously, the tool can distinguish virulent and non-virulent proteins. Moreover, the pathogenic proteins were identified from the functionally annotated proteome using the MP3 tool [43]. This software exploits a combined SVM-HMM approach while identifying proteins from both genomic and metagenomic databases with high accuracy, efficiency, and sensitivity.

Protein-protein interaction network analysis

The protein-protein interaction (PPI) (https://string-db.org/) network between the functionally annotated HPs was visualized using the STRING database [44]. In this study, the available S. pneumoniae D39 in the STRING database was selected as the reference genome and an interconnected PPI network was constructed. The basis for the network lay in high-throughput lab experiments, gene expression data, and computational data. The network was built with the default confidence parameters.

Excavation of druggable proteins

The druggability of the annotated proteins was assessed through Drugbank BlastP [45]. Drugbank (https://go.drugbank.com/) is a comprehensive bioinformatics and cheminformatics resource containing relevant information about drugs and their corresponding targets.

Assessment of Structural proteins

Models for the selected proteins were obtained from Robetta (https://robetta.bakerlab.org/). Robetta analyzes the putative domains of the submitted sequences and generates 3-dimensional structural models [46]. The models were further evaluated using SWISS Structure assessment (https://swissmodel.expasy.org/). SWISS Structure assessment includes Ramachandran Plots and MolProbity scores [47]. Ramachandran Plots visualize energetically favored regions for backbone dihedral angles against the amino acids present in the protein structure while MolProbity examines the quality of protein models for both nucleic acids and proteins at global and local levels [48].

Molecular docking simulation

The Canonical Simplified molecular-input line-entry system (SMILES) of the interacted drugs were collected from DrugBank BLASTp. These canonical smiles were converted into PDB files via the CACTUS online smile translator (https://cactus.nci.nih.gov/translate/). The druggable proteins were then docked using AutoDock within PyRx software, which is a combination of several tools necessary for Molecular Docking [49]. The protein-ligand interactions were then visualized using Discovery Studio Visualizer (https://discover.3ds.com/).

Molecular Dynamics (MD) simulation

To evaluate the stability of the complexes under physiological conditions, 100 ns Molecular Dynamics (MD) simulation was carried out using GROningen MAchine for Chemical Simulations (GROMACS version 5.1.1). The GROMOS96 43a1 force-field was applied to the protein-ligand complexes. The physiological condition of the system was defined as (300 K, pH 7.4, 0.9% NaCl). The structures were solvated in a dodecahedral box of the SPC (simple point charge) water model with its edges at a 1nm distance from the protein surface. The overall charge of the system was neutralized through the addition of 2 sodium ions using the genion module. Energy minimization of the neutralized system was carried out using the steepest descent minimization algorithm with a maximum number of minimization steps to perform was set at 50000. The ligand was restrained before carrying out the isothermal-isochoric (NVT) equilibration of the system for 100 ps with a short-range electrostatic cutoff value of 1.2 nm. Isobaric (NPT) equilibration of the system was carried out for 100 ps following the NVT with a short-range van der Waals cutoff fixed at 1.2 nm. Later, a 10 ns molecular dynamic simulation was run using periodic boundary conditions and a time integration step of 2 fs. The energy of the system was saved every 100 ps. For calculating the long-range electrostatic potential, the Particle Mesh Ewald (PME) method was applied. The short-range van der Waals cutoff was kept at 1.2 A modified Berendsen thermostat was used to control simulation temperature while the pressure was kept constant using the Parrinello-Rahman algorithm. The simulation time step was selected as 2.0 fs. The snapshot interval was set to 100 ps for analyzing the trajectory data. Finally, all of the trajectories were concatenated to calculate and plot root mean square deviation (RMSD), root mean square fluctuation (RMSF), the radius of gyration (Rg), and solvent accessible surface area (SASA) data. Root Mean Square Deviation (RMSD) calculation was performed to evaluate when a system attains equilibrium. The “rms” module built into the GROMACS software was utilized to extract RMSD information throughout the simulation. The results were plotted graphically using the ggplot2 package of R (https://ggplot2.tidyverse.org/).

Room Mean Square Fluctuation (RMSF) is used to determine the flexibility of a certain region of the protein. The radius of gyration of our proteins was measured to determine their degree of compactness. A relatively steady value of the radius of gyration means stable folding of a protein. Fluctuation of the radius of gyration implies the unfolding of the protein. The “gyrate” module was used to generate the radius of gyration graphs for our proteins.

Hydrophobic interactions composed of non-polar amino acids are crucial for maintaining the stability of the hydrophobic core of proteins. They do so by covering the non-polar amino acids within the hydrophobic cores and keeping them at a distance from the solvent. Solvent Accessible Surface Area (SASA) is used in MD simulations to predict the hydrophobic core stability of proteins.

Results

In the study, the core genome of S. pneumoniae that encodes hypothetical proteins (HPs) was computationally annotated and analyzed to identify potential drug targets. The workflow has been illustrated in Fig 1.

[Figure omitted. See PDF.]

Identification of proteins from core and pan genome

Out of 22 different Streptococci species from EDGAR 3.0, 9084 genes from pan genome were identified and several were visualized by the BioCircos tool. Subsequently, 498 genes were extracted from the core genome (S1 Data). The core genome dataset of S. pneumoniae encoded 28 HPs (Fig 2) (S2 Data).

[Figure omitted. See PDF.]

The brown bar of the circular plot represents the core genes.

Functional analysis of the hypothetical proteins

The gene ontology (GO) analysis of the 28 HPs revealed that 21 annotated proteins had been involved in biological processes (BPs), molecular functions (MFs) and cellular components (CCs) and the remaining belonged to uncharacterized domains/ protein families (S3 Data). The HPs corresponded to 20 GO terms in total while some of the HPs possessed more than one GO term (Fig 3 and Table 1). HPs were homologs of enzymes, binding proteins, transporters, and regulators (Fig 3). A high proportion of HPs exhibited MFs such as aminoacyl-tRNA ligase activity, RNA binding, DNA binding, hydrolase activity, transferase activity, tRNA binding, methylation, nucleic acid binding, ATP binding, metal ions binding, Flavin adenine dinucleotide binding, tRNA modification, and metallopeptidase activity, indicating they may take part in translation and transcription of the bacteria. The study also predicted 6 proteins as members of the integral membrane protein family.

[Figure omitted. See PDF.]

Proteins involved in molecular functions were relatively higher compared to others.

[Figure omitted. See PDF.]

Physicochemical characterization

The physicochemical properties of the 28 HPs are provided in Table 2. The length of the analyzed sequence ranged from 71 to 560 AA with a molecular weight (MW) ranging from 8462.34 and 64392.6 Da. The isoelectric point (pI) of the proteins ranged from 3.92 to 10.07 with an average value of 6.20. The isoelectric point determines the protein load. At this pH, the protein does not have any charge and therefore does not move in the electric field when a direct current is passed through it [50]. MW and pI parameters are important for purification and crystallization during in vivo experimentation. The extinction coefficient measures the amount of light that a protein absorbs at a certain wavelength. In silico identification of the extinction coefficient is necessary to evaluate the study of protein-protein interaction in solution [51]. The extinction coefficient of the provided hypothetical proteins ranged from 1615 to 57190 at 280 nm. Here, HPs possessed high extinction coefficient values that indicate the proteins constitute a high concentration of cysteine, tyrosine, and tryptophan. The instability index determines the stability of the protein in a test tube. An instability index value of less than 40 predicts the protein to be stable and over 40 adopts the protein to be unstable [52]. With variable ranges of instability index, more than 50% of the proteins were noticed to be stable with a mean value of 39.21 in this study. The aliphatic index estimates the thermostability of the protein depending on the area occupied by the aliphatic chain. Proteins with high aliphatic index show stability in high temperatures, whereas proteins with low aliphatic value are less flexible and thermally unstable [53]. Here, the aliphatic index of the proteins ranged from 32.66 to 1135.05 with a mean value of 92.16. Most of the proteins were observed with an aliphatic index value above 40 which means most of them are thermally stable. The GRAVY of the protein indicates the interaction of proteins with water. It is calculated by adding all the hydropathy values of the amino acids followed by dividing the number of residues in the sequence. The outcome for GRAVY predicts the presence of the protein-water interactions between all the proteins except two. Therefore, the majority of them are water-soluble.

[Figure omitted. See PDF.]

Prediction of essential proteins

Identification of essential genes contributes a principal role in the development of new drugs or vaccines that inhibit the dissemination of resistant strains [54]. Through a similarity search against the essential proteins of bacteria present in the DEG database, 16 HPs were considered as essential proteins with an E-value<0.0001 and bit score>100 as shown in Table 3.

[Figure omitted. See PDF.]

More than half of the proteins were found to be essential.

Virulence factors and pathogenicity

Virulence factors are the basis of bacterial infections. The virulent and pathogenic proteins help microbes invade the host and manipulate the host immune system for its survival [55]. The comprehensive predictive data for virulence is provided. The VirulentPred server identified 54% of the HPs as having virulent properties while 46% owned non virulent properties. Moreover, the MP3 webserver recognized 8 proteins out of 28 as pathogenic. Among these 8 proteins, 4 portrayed both virulence and pathogenic properties (Fig 4) (S4 Data).

[Figure omitted. See PDF.]

Analysis of the components of the PPI network

Protein-Protein interactions are crucial for all cellular functions. The association between the hypothetical proteins of S. pneumoniae and other proteins is visualized in Fig 5. Among 28 proteins, 20 proteins were involved within the network with a default confidence value. The network contained 28 nodes (proteins) and 32 edges (interactions). However, two of the proteins (PCS8235_RS05845 and PCS8235_RS02595) were observed to have the highest number of interactions with other proteins. The network also revealed some proteins that interact with Pneumococcal vaccine antigen A, Ribosomal Silencing Factor, Cof family domain, CBS domain, and cell division protein Zap A, putative Anti-CRISPR, Regulator of chromosomal segregation (RocS) proteins, and so on. Further analysis of the functional enrichments in the network exhibited proteins that work as Ras family proteins, NusA-like KH domain, Putative fluoride ion transporter CrcB and P-loop containing nucleoside triphosphate hydrolase (S5 Data).

[Figure omitted. See PDF.]

Examination of the subcellular localization and unconventional protein secretion

During drug discovery, knowledge of subcellular localization can be leveraged to significantly improve the identification of druggable targets [56]. Proteins are categorized as drugs or vaccine targets based on subcellular location. Considering the results of the CELLO and PSORTb 3.0 tools, the majority of the proteins were characterized as cytoplasmic and 6 as membrane-bound with a sole exception being located in the extracellular matrix (Table 4). This implies that most of the proteins prefer the cytoplasmic region for their activity. Furthermore, OutCyte 1.0 data predicted that 64% of all the HPs undergo potential unconventional protein secretion (UPS) and the rest possessed signal peptide, the transmembrane domain, and intracellular secretion. Additionally, SecretomeP 2.0 identified 20 non classical secreted proteins. These proteins lacked a signal peptide and had a threshold SecP value above 0.5. All the predicted results for CELLO, PSORTb 3.0, OutCyte 1.0, and SecretomeP 2.0 are presented in Table 4.

[Figure omitted. See PDF.]

Druggability of the proteins

The druggability prediction of HPs classified three proteins with homologous genes (PCS8235_RS02820, PCS8235_RS04650, and PCS8235_RS10805) available in the Drugbank database (Table 5). DrugBank contains arrays of gene families that have been successfully targeted by drugs with required affinity. One of the HPs (PCS8235_RS04650) had two homologous compounds. Three of the predictive druggable proteins were filtered into two (PCS8235_RS02820, PCS8235_RS04650) based on a cut-off E-value of <0.0001 and bit-score >100.

[Figure omitted. See PDF.]

Protein model building

Robetta predicted 5 models for each druggable protein. Robetta modeling program was implemented to predict the best model among all. The Ramachandran plot percentage and the Molprobity score of the proteins are tabulated below (S6 Data). MolProbity score should be as low as possible and Ramachandran plot >98% is favored for an ideal protein structure. Model 2 for both of the proteins was selected using these criteria. Model 2 of both of the proteins (PCS8235_RS02820 and PCS8235_RS04650) were selected based on Ramachandran plot percentages of 98.77% and 99.24% with a low morbidity score of 2.64 and 2.69 respectively.

Elucidation of molecular docking

Blind docking was performed to evaluate the binding energy between each ligand and receptors. Two of the expected potential HPs, PCS8235_RS02820 and PCS8235_RS04650 were docked against the interacted drugs from DrugBank (Fig 6). The binding affinity between druggable targets and commercially available drugs was satisfactory. Protein PCS8235_RS02820 docked against 2-(N-morpholino) ethanesulfonic acid with docking energy of -4.4 kcal/mol. Protein PCS8235_RS 04650 bonded with Artenimol and Aspartate beryllium trifluoride with the binding affinity of -9.5 kcal/mol and -4.7 kcal/mol respectively.

[Figure omitted. See PDF.]

Interaction between the best drug for the HPs in S. pneumoniae and the active site amino acid residues of each of the HPs (a) PCS8235_RS04650- Artenimol (b) PCS8235_RS04650- Aspartate beryllium trifluoride (c) PCS8235_RS02820-2-(N-morpholino) ethanesulfonic acid.

Examination of molecular dynamics simulation

The RMSD graph of PCS8235_RS02820 and PCS8235_RS02820-drug complex demonstrated that the backbone of the proteins started to deviate around 12 ns and overlapped around 75 ns (Fig 7). The RMSF and SASA graphs of this protein showed that the binding of the drug altered mobility and solvent accessibility in the protein (Figs 8 and 9). Rg analysis unveiled that the compactness of the protein and drug-protein complex was most of the time similar during the simulations (Fig 10). For PCS8235_RS04650 and PCS8235_RS04650-drug bound complexes also gave similar results. The RMSD backbones did not overlap most of the time through 100ns. The RMSF values or mobility with protein compactness (Rg) of the structures altered due to drugs. Slight solvent accessibility was observed in SASA analysis.

[Figure omitted. See PDF.]

RMSD analysis of (A) PCS8235_RS02820 (green) and PCS8235_RS02820-2-(N-morpholino) ethanesulfonic acid complex (purple) (B) PCS8235_RS04650 (red) and PCS8235_RS04650 –Artenimol (purple) and PCS8235_RS04650-Aspartate beryllium trifluoride (gray) at 100ns.

[Figure omitted. See PDF.]

RMSF analysis of (A) PCS8235_RS02820 (brown) and PCS8235_RS02820-2-(N-morpholino) ethanesulfonic acid complex (blue) (B) PCS8235_RS04650 (blue) and PCS8235_RS04650 –Artenimol (golden) and PCS8235_RS04650-Aspartate beryllium trifluoride (green) at 100 ns.

[Figure omitted. See PDF.]

SASA calculation of (A) PCS8235_RS02820 (golden) and PCS8235_RS02820-2-(N-morpholino) ethanesulfonic acid complex (green) (B) PCS8235_RS04650 (blue) and PCS8235_RS04650 –Artenimol (green) and PCS8235_RS04650-Aspartate beryllium trifluoride (red) at 100ns.

[Figure omitted. See PDF.]

Rg measurement of (A) PCS8235_RS02820 (green) and PCS8235_RS02820-2-(N-morpholino) ethanesulfonic acid complex (red) (B) PCS8235_RS04650 (green) and PCS8235_RS04650 –Artenimol (red) and PCS8235_RS04650-Aspartate beryllium trifluoride (yellow) at 100 ns.

Discussion

S. pneumoniae has been the prime cause of death via community-acquired pneumonia (CAP) worldwide [1]. Combating the mortality and morbidity rate of CAP is becoming more challenging due to the increasing prevalence of antimicrobial resistance [57]. Following the drastic rise in antimicrobial resistance, current research on producing new antibiotics seems time-consuming, and costly. For several years, genomics and evolutionary studies on this organism were endeavored to determine new druggable targets [58, 59]. Yet, few studies have been published on the unexplored realm of proteins known as HPs. About nearly 30% of S. pneumoniae encodes HPs. Some homologs of HPs are conserved within various streptococci species [60]. In the present study, while analyzing the core genome of the Streptococcus species, we have found 498 core genes among 22 Streptococci species where 28 belonged to HPs (around 5.6%). Since they are widely distributed in the Streptococcus genus through evolution, they might be critical for the survival of the bacteria. Proper characterization of these HPs could give plausible targets against S. pneumoniae pathogen. Previous studies demonstrated that HPs play a vital role in biofilm formation, pathogenesis, and virulence of tooth decaying S. mutans [61] whereas, in this study, we have also identified the protein PCS8235_RS05215 that plays a key role in biofilm formation. Understanding the functional annotation of these HPs is essential to comprehending the molecular and biological mechanisms of all species. Bioinformatics approaches aiding functional annotation and computational drug designing can be a better way to identify the potential drug targets from these HPs.

When the core HPs went under gene ontology-based characterization, 75% of the proteins revealed their characteristic. Most of them were enzymes and binding proteins. Enzymes are involved in several biochemical reactions and self-defense mechanisms playing a key role in bacterial survival inside the host. Here, protein PCS8235_RS3915 and PCS8235_RS4650 were involved in hydrolase activity. Protein PCS8235_RS3915 is a DNA repair RadC homolog while PCS8235_RS4650 is a haloacid dehalogenase homolog. Haloacid dehalogenase homologs are a superfamily of enzymes involved in cellular functions starting from amino acid biosynthesis to detoxification [62, 63]. In prokaryotes, DNA repair RadC homolog helps in DNA damage repair after UV and X-ray radiation [64]. In the same way, protein PCS8235_RS2800 is related to branched-chain amino acid aminotransferase activity and protein PCS8235_RS6025 is associated with methyltransferase activity. Hence, enzymes can be assumed as a druggable target for different therapeutics. The preceding proteins (PCS8235_RS01335, PCS8235_RS01340, PCS8235_RS02405, PCS8235_RS03345, PCS8235_RS06025 and PCS8235_RS06945) also take part in nucleic acid-binding. Other proteins (PCS8235_RS02405, PCS8235_RS03345, PCS8235_RS06025, and PCS8235_RS06945) were somehow involved in the central dogma of the bacterial pathogen. Protein PCS8235_RS05215 denoted transcriptional regulator YlbF, YmcA that formed a complex with YaaT regulated competence and biofilm formation in Bacillus subtilis [65]. Furthermore, targeted proteins must be essential to the concerned pathogen since essential genes strongly support the cellular life of living microorganisms [66]. For the study, more than half of the proteins were identified as essential proteins of S. pneumoniae. Later, the virulence and pathogenicity of the essential HPs were noted. Out of the 19 proteins analyzed for virulence and pathogenic characteristics, 21% were found to be both virulent and pathogenic.

Moreover, the integrated roles of the bacterial proteins were elucidated via protein-protein interactions. Protein-protein interaction is required to understand the specific functions of bacterial protein which might be influenced by the neighboring proteins. This helps in finding its significance in bacterial physiology. In the provided network, the protein PCS8235_RS05845 formed the highest nodes (8 interactions) with other HPs that might have critical cellular functions and the potential to go under unconventional protein secretions. Most of the HPs play a physiological role in the cytoplasmic region and have hydrophilic properties. Therefore, they can be considered as drug targets rather than vaccines [29, 67]. This cytoplasmic characteristic was evaluated by subcellular localization and physicochemical properties. PCS8235_RS05845 also has TRP structural motif along with virulent and pathogenic properties. Previous studies reported that some TRP-containing proteins are involved directly in virulence-associated functions, especially, in maintaining proper functions of the type III secretion system and class II chaperones [67]. The next highest interactive protein is PCS8235_RS02595 (7 interactions). Although this protein is uncharacterized still now, its interaction with other proteins can suggest its significant role in the network. PCS8235_RS02595 is interconnected to Pneumococcal vaccine antigen A (PCS8235_RS03660), cell division protein Zap A (PCS8235_RS00650), putative Anti-CRISPR (PCS8235_RS04815) and Regulator of chromosomal segregation (RocS) (PCS8235_RS03470) proteins. Zap A is a Z-ring-associated protein found in Bacillus subtilis [68]. This protein binds to tubulin-like protein FtsZ at cytokinesis and therefore takes part during bacterial cell division [69]. Putative Anti-CRISPR may inhibit the bacterial CRISPR-Cas system. RocS participates in chromosome segregation and nucleoid protection in S. pneumoniae [70].

The druggability of the essential HPs was checked to ensure the proteins’ susceptibility toward the active site of the small inhibitor molecules. In other studies, the drug targets in S. pneumoniae were explored; however, it has only screened the potential drug targets for the pathogen [71]. Here, only the compounds with both druggability and essentiality are chosen for further analysis. The 3D structure of the targeted proteins was then analyzed and the best model was selected for molecular docking. The basis of docking was to identify the relationship between the hypothetical proteins and available druggable ligands. A blind docking test for each of the whole proteins with expected drug candidates reveals satisfactory drug affinity towards the ligands. The strongest binding affinity could be noticed between PCS8235_RS04650 and Artenimol with the binding affinity of -9.5 kcal/mol. MD simulations also support that these druggable compounds will change the targeted bacterial proteins by impeding their general physiology. Therefore, this predictive essential protein could be an eligible drug target for the development of new pneumococcal therapeutics.

Conclusion

This study focused on the HPs that are found in the core genome of S. pneumoniae and have a major role in the survival, virulence, and pathogenesis of the pathogen. The majority of the observed proteins are mainly involved in bacterial transcription and translation while a high proportion of the proteins were recognized as essential proteins. Hence they make up good targets for broad spectrum anti-streptococcal drug development.

Supporting information

S1 Data. List of all 498 genes that constitute the core genome of Streptococcus sp.

https://doi.org/10.1371/journal.pone.0272945.s001

(XLSX)

S2 Data. The core proteome of Streptococcus strains.

https://doi.org/10.1371/journal.pone.0272945.s002

(XLSX)

S3 Data. Result of gene ontology analysis (biological processes, molecular functions, and cellular components) of the 28 hypothetical proteins.

https://doi.org/10.1371/journal.pone.0272945.s003

(XLSX)

S4 Data. Result of virulence and pathogenicity analysis.

https://doi.org/10.1371/journal.pone.0272945.s004

(XLSX)

S5 Data. Functional enrichment of the interconnected proteins generated by STRING webserver.

https://doi.org/10.1371/journal.pone.0272945.s005

(DOCX)

S6 Data. Ramachandran plot statistics of the hypothetical proteins PCS8235_RS02820 and PCS8235_RS04650.

https://doi.org/10.1371/journal.pone.0272945.s006

(XLSX)

Citation: Chowdhury ZM, Bhattacharjee A, Ahammad I, Hossain MU, Jaber AA, Rahman A, et al. (2022) Exploration of Streptococcus core genome to reveal druggable targets and novel therapeutics against S. pneumoniae. PLoS ONE 17(8): e0272945. https://doi.org/10.1371/journal.pone.0272945

About the Authors:

Zeshan Mahmud Chowdhury

Roles: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft

Affiliation: Bioinformatics Division, National Institute of Biotechnology, Dhaka, Bangladesh

https://orcid.org/0000-0002-6263-4960

Arittra Bhattacharjee

Roles: Conceptualization, Writing – review & editing

Affiliation: Bioinformatics Division, National Institute of Biotechnology, Dhaka, Bangladesh

Ishtiaque Ahammad

Roles: Project administration, Resources, Supervision, Writing – review & editing

Affiliation: Bioinformatics Division, National Institute of Biotechnology, Dhaka, Bangladesh

Mohammad Uzzal Hossain

Roles: Project administration, Supervision, Writing – review & editing

Affiliation: Bioinformatics Division, National Institute of Biotechnology, Dhaka, Bangladesh

Abdullah All Jaber

Roles: Data curation

Affiliation: Department of Biochemistry & Microbiology, North South University, Dhaka, Bangladesh

https://orcid.org/0000-0001-8964-7350

Anisur Rahman

Roles: Visualization

Affiliation: Bioinformatics Division, National Institute of Biotechnology, Dhaka, Bangladesh

https://orcid.org/0000-0001-8279-9810

Preonath Chondrow Dev

Roles: Visualization

Affiliation: Bioinformatics Division, National Institute of Biotechnology, Dhaka, Bangladesh

Md. Salimullah

Roles: Supervision, Writing – review & editing

Affiliation: Molecular Biotechnology Division, National Institute of Biotechnology, Dhaka, Bangladesh

https://orcid.org/0000-0001-9895-9233

Chaman Ara Keya

Roles: Project administration, Supervision, Writing – review & editing

E-mail: [email protected]

Affiliation: Department of Biochemistry & Microbiology, North South University, Dhaka, Bangladesh

References

1. Ferreira-Coimbra J, Sarda C, Rello J. Burden of Community-Acquired Pneumonia and Unmet Clinical Needs. Adv Ther. 2020;37: 1302–1318. pmid:32072494

2. Metlay JP, Waterer GW, Long AC, Anzueto A, Brozek J, Crothers K, et al. Diagnosis and Treatment of Adults with Community-acquired Pneumonia. An Official Clinical Practice Guideline of the American Thoracic Society and Infectious Diseases Society of America. Am J Respir Crit Care Med. 2019;200: e45–e67. pmid:31573350

3. Rudan I, O’Brien KL, Nair H, Liu L, Theodoratou E, Qazi S, et al. Epidemiology and etiology of childhood pneumonia in 2010: estimates of incidence, severe morbidity, mortality, underlying risk factors and causative pathogens for 192 countries. J Glob Health. 2013;3: 010401. pmid:23826505

4. Brooks LRK, Mias GI. Streptococcus pneumoniae’s Virulence and Host Immunity: Aging, Diagnostics, and Prevention. Front Immunol. 2018;9: 1366. pmid:29988379

5. GBD 2016 Lower Respiratory Infections Collaborators. Estimates of the global, regional, and national morbidity, mortality, and aetiologies of lower respiratory infections in 195 countries, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Infect Dis. 2018;18: 1191–1210. pmid:30243584

6. Morens DM, Taubenberger JK, Fauci AS. Predominant Role of Bacterial Pneumonia as a Cause of Death in Pandemic Influenza: Implications for Pandemic Influenza Preparedness. J Infect Dis. 2008;198: 962–970. pmid:18710327

7. Knight GM, Glover RE, McQuaid CF, Olaru ID, Gallandat K, Leclerc QJ, et al. Antimicrobial resistance and COVID-19: Intersections and implications. eLife. 2021;10: e64139. pmid:33588991

8. Torres A, Peetermans WE, Viegi G, Blasi F. Risk factors for community-acquired pneumonia in adults in Europe: a literature review. Thorax. 2013;68: 1057–1065. pmid:24130229

9. Müllerova H, Chigbo C, Hagan GW, Woodhead MA, Miravitlles M, Davis KJ, et al. The natural history of community-acquired pneumonia in COPD patients: a population database analysis. Respir Med. 2012;106: 1124–1133. pmid:22621820

10. Taylor SN, Sanders CV. Unusual manifestations of invasive pneumococcal infection. Am J Med. 1999;107: 12S–27S. pmid:10451005

11. de Benedictis FM, Kerem E, Chang AB, Colin AA, Zar HJ, Bush A. Complicated pneumonia in children. Lancet Lond Engl. 2020;396: 786–798. pmid:32919518

12. Yf T, Yh K. Necrotizing pneumonia: a rare complication of pneumonia requiring special consideration. Curr Opin Pulm Med. 2012;18. pmid:22388585

13. Torres A, Menéndez R, Wunderink RG. Bacterial Pneumonia and Lung Abscess. Murray Nadels Textb Respir Med. 2016; 557–582.e22.

14. Tan TQ. Pediatric invasive pneumococcal disease in the United States in the era of pneumococcal conjugate vaccines. Clin Microbiol Rev. 2012;25: 409–419. pmid:22763632

15. McIntosh EDG, Reinert RR. Global prevailing and emerging pediatric pneumococcal serotypes. Expert Rev Vaccines. 2011;10: 109–129. pmid:21162625

16. Wang Y, Li J, Wang Y, Gu W, Zhu F. Effectiveness and practical uses of 23-valent pneumococcal polysaccharide vaccine in healthy and special populations. Hum Vaccines Immunother. 2018;14: 1003–1012. pmid:29261406

17. Zhang M, Chen H, Wu F, Li Q, Lin Q, Cao H, et al. Heightened Willingness toward Pneumococcal Vaccination in the Elderly Population in Shenzhen, China: A Cross-Sectional Study during the COVID-19 Pandemic. Vaccines. 2021;9: 212. pmid:33802327

18. Miranda C, Silva V, Capita R, Alonso-Calleja C, Igrejas G, Poeta P. Implications of antibiotics use during the COVID-19 pandemic: present and future. J Antimicrob Chemother. 2020;75: 3413–3416. pmid:32830266

19. Cherazard R, Epstein M, Doan T-L, Salim T, Bharti S, Smith MA. Antimicrobial Resistant Streptococcus pneumoniae: Prevalence, Mechanisms, and Clinical Implications. Am J Ther. 2017;24: e361–e369. pmid:28430673

20. Larsson M, Nguyen HQ, Olson L, Tran TK, Nguyen TV, Nguyen CTK. Multi-drug resistance in Streptococcus pneumoniae among children in rural Vietnam more than doubled from 1999 to 2014. Acta Paediatr Oslo Nor 1992. 2021;110: 1916–1923. pmid:33544434

21. Liñares J, Ardanuy C, Pallares R, Fenoll A. Changes in antimicrobial resistance, serotypes and genotypes in Streptococcus pneumoniae over a 30-year period. Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis. 2010;16: 402–410. pmid:20132251

22. Chaguza C, Cornick JE, Everett DB. Mechanisms and impact of genetic recombination in the evolution of Streptococcus pneumoniae. Comput Struct Biotechnol J. 2015;13: 241–247. pmid:25904996

23. Straume D, Stamsås GA, Håvarstein LS. Natural transformation and genome evolution in Streptococcus pneumoniae. Infect Genet Evol J Mol Epidemiol Evol Genet Infect Dis. 2015;33: 371–380. pmid:25445643

24. Salvadori G, Junges R, Morrison DA, Petersen FC. Competence in Streptococcus pneumoniae and Close Commensal Relatives: Mechanisms and Implications. Front Cell Infect Microbiol. 2019;9. https://www.frontiersin.org/article/10.3389/fcimb.2019.00094 pmid:31001492

25. Fani F, Leprohon P, Légaré D, Ouellette M. Whole genome sequencing of penicillin-resistant Streptococcus pneumoniae reveals mutations in penicillin-binding proteins and in a putative iron permease. Genome Biol. 2011;12: R115. pmid:22108223

26. Schroeder MR, Stephens DS. Macrolide Resistance in Streptococcus pneumoniae. Front Cell Infect Microbiol. 2016;6. https://www.frontiersin.org/article/10.3389/fcimb.2016.00098 pmid:27709102

27. Jacoby GA. Mechanisms of resistance to quinolones. Clin Infect Dis Off Publ Infect Dis Soc Am. 2005;41 Suppl 2: S120–126. pmid:15942878

28. Mitchell AM, Mitchell TJ. Streptococcus pneumoniae: virulence factors and variation. Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis. 2010;16: 411–418. pmid:20132250

29. Rabbi MdF, Akter SA, Hasan MdJ, Amin A. In Silico Characterization of a Hypothetical Protein from Shigella dysenteriae ATCC 12039 Reveals a Pathogenesis-Related Protein of the Type-VI Secretion System. Bioinforma Biol Insights. 2021;15: 11779322211011140. pmid:33994781

30. Imam N, Alam A, Ali R, Siddiqui MF, Ali S, Malik MZ, et al. In silico characterization of hypothetical proteins from Orientia tsutsugamushi str. Karp uncovers virulence genes. Heliyon. 2019;5: e02734. pmid:31720472

31. Nan J, Brostromer E, Liu X-Y, Kristensen O, Su X-D. Bioinformatics and Structural Characterization of a Hypothetical Protein from Streptococcus mutans: Implication of Antibiotic Resistance. PLOS ONE. 2009;4: e7245. pmid:19798411

32. Cui Y, Chen X, Luo H, Fan Z, Luo J, He S, et al. BioCircos.js: an interactive Circos JavaScript library for biological data visualization on web applications. Bioinforma Oxf Engl. 2016;32: 1740–1742. pmid:26819473

33. Dieckmann MA, Beyvers S, Nkouamedjo-Fankep RC, Hanel PHG, Jelonek L, Blom J, et al. EDGAR3.0: comparative genomics and phylogenomics on a scalable infrastructure. Nucleic Acids Res. 2021;49: W185–W192. pmid:33988716

34. Araujo FA, Barh D, Silva A, Guimarães L, Ramos RTJ. GO FEAT: a rapid web-based functional annotation tool for genomic and transcriptomic data. Sci Rep. 2018;8: 1794. pmid:29379090

35. Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, et al. Protein Identification and Analysis Tools on the ExPASy Server. In: Walker JM, editor. The Proteomics Protocols Handbook. Totowa, NJ: Humana Press; 2005. pp. 571–607. https://doi.org/10.1385/1-59259-890-0:571

36. Yu C-S, Lin C-J, Hwang J-K. Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci Publ Protein Soc. 2004;13: 1402–1406. pmid:15096640

37. PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes—PMC. [cited 10 May 2022]. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2887053/

38. Damte D, Suh J-W, Lee S-J, Yohannes SB, Hossain MdA, Park S-C. Putative drug and vaccine target protein identification using comparative genomic analysis of KEGG annotated metabolic pathways of Mycoplasma hyopneumoniae. Genomics. 2013;102: 47–56. pmid:23628646

39. Bendtsen JD, Kiemer L, Fausbøll A, Brunak S. Non-classical protein secretion in bacteria. BMC Microbiol. 2005;5: 58. pmid:16212653

40. Zhao L, Poschmann G, Waldera-Lupa D, Rafiee N, Kollmann M, Stühler K. OutCyte: a novel tool for predicting unconventional protein secretion. Sci Rep. 2019;9: 19448. pmid:31857603

41. Zhang R, Ou H-Y, Zhang C-T. DEG: a database of essential genes. Nucleic Acids Res. 2004;32: D271–272. pmid:14681410

42. Garg A, Gupta D. VirulentPred: a SVM based prediction method for virulent proteins in bacterial pathogens. BMC Bioinformatics. 2008;9: 62. pmid:18226234

43. Gupta A, Kapil R, Dhakan DB, Sharma VK. MP3: A Software Tool for the Prediction of Pathogenic Proteins in Genomic and Metagenomic Data. PLoS ONE. 2014;9: e93907. pmid:24736651

44. Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011;39: D561–D568. pmid:21045058

45. Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, et al. DrugBank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 2008;36: D901–D906. pmid:18048412

46. Kim DE, Chivian D, Baker D. Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004;32: W526–531. pmid:15215442

47. Schwede T, Kopp J, Guex N, Peitsch MC. SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res. 2003;31: 3381–3385. pmid:12824332

48. Hollingsworth SA, Karplus PA. A fresh look at the Ramachandran plot and the occurrence of standard structures in proteins. Biomol Concepts. 2010;1: 271–283. pmid:21436958

49. Small-Molecule Library Screening by Docking with PyRx | Springer Nature Experiments. [cited 10 May 2022]. https://experiments.springernature.com/articles/10.1007/978-1-4939-2269-7_19

50. Kirkwood J, Hargreaves D, O’Keefe S, Wilson J. Using isoelectric point to determine the pH for initial protein crystallization trials. Bioinformatics. 2015;31: 1444–1451. pmid:25573921

51. Gill SC, von Hippel PH. Calculation of protein extinction coefficients from amino acid sequence data. Anal Biochem. 1989;182: 319–326. pmid:2610349

52. Gamage DG, Gunaratne A, Periyannan GR, Russell TG. Applicability of Instability Index for In vitro Protein Stability Prediction. Protein Pept Lett. 2019;26: 339–347. pmid:30816075

53. Dutta B, Banerjee A, Chakraborty P, Bandopadhyay R. In silico studies on bacterial xylanase enzyme: Structural and functional insight. J Genet Eng Biotechnol. 2018;16: 749–756. pmid:30733796

54. Davies J, Davies D. Origins and Evolution of Antibiotic Resistance. Microbiol Mol Biol Rev MMBR. 2010;74: 417–433. pmid:20805405

55. Sharma AK, Dhasmana N, Dubey N, Kumar N, Gangwal A, Gupta M, et al. Bacterial Virulence Factors: Secreted for Survival. Indian J Microbiol. 2017;57: 1–10. pmid:28148975

56. Rajendran L, Knölker H-J, Simons K. Subcellular targeting strategies for drug design and delivery. Nat Rev Drug Discov. 2010;9: 29–42. pmid:20043027

57. Assefa M, Tigabu A, Belachew T, Tessema B. Bacterial profile, antimicrobial susceptibility patterns, and associated factors of community-acquired pneumonia among adult patients in Gondar, Northwest Ethiopia: A cross-sectional study. PLOS ONE. 2022;17: e0262956. pmid:35104293

58. Chaguza C, Yang M, Cornick JE, du Plessis M, Gladstone RA, Kwambana-Adams BA, et al. Bacterial genome-wide association study of hyper-virulent pneumococcal serotype 1 identifies genetic variation associated with neurotropism. Commun Biol. 2020;3: 1–12.

59. Andam CP, Hanage WP. Mechanisms of genome evolution of Streptococcus. Infect Genet Evol J Mol Epidemiol Evol Genet Infect Dis. 2015;33: 334–342. pmid:25461843

60. Hiller NL, Janto B, Hogg JS, Boissy R, Yu S, Powell E, et al. Comparative genomic analyses of seventeen Streptococcus pneumoniae strains: Insights into the pneumococcal supragenome. J Bacteriol. 2007;189: 8186–8195. pmid:17675389

61. Krzyściak W, Jurczak A, Kościelniak D, Bystrowska B, Skalniak A. The virulence of Streptococcus mutans and the ability to form biofilms. Eur J Clin Microbiol Infect Dis Off Publ Eur Soc Clin Microbiol. 2014;33: 499–515. pmid:24154653

62. Koonin EV, Tatusov RL. Computer analysis of bacterial haloacid dehalogenases defines a large superfamily of hydrolases with diverse specificity. Application of an iterative approach to database search. J Mol Biol. 1994;244: 125–132. pmid:7966317

63. Kuznetsova E, Nocek B, Brown G, Makarova KS, Flick R, Wolf YI, et al. Functional Diversity of Haloacid Dehalogenase Superfamily Phosphatases from Saccharomyces cerevisiae: BIOCHEMICAL, STRUCTURAL, AND EVOLUTIONARY INSIGHTS. J Biol Chem. 2015;290: 18678–18698. pmid:26071590

64. Attaiech L, Granadel C, Claverys J-P, Martin B. RadC, a Misleading Name? J Bacteriol. 2008;190: 5729–5732. pmid:18556794

65. Carabetta VJ, Tanner AW, Greco TM, Defrancesco M, Cristea IM, Dubnau D. A complex of YlbF, YmcA and YaaT regulates sporulation, competence and biofilm formation by accelerating the phosphorylation of Spo0A. Mol Microbiol. 2013;88: 283–300. pmid:23490197

66. Peng C, Lin Y, Luo H, Gao F. A Comprehensive Overview of Online Resources to Identify and Predict Bacterial Essential Genes. Front Microbiol. 2017;8. https://www.frontiersin.org/article/10.3389/fmicb.2017.02331 pmid:29230204

67. Cerveny L, Straskova A, Dankova V, Hartlova A, Ceckova M, Staud F, et al. Tetratricopeptide Repeat Motifs in the World of Bacterial Pathogens: Role in Virulence Mechanisms. Infect Immun. 2013;81: 629–635. pmid:23264049

68. Buss J, Coltharp C, Huang T, Pohlmeyer C, Wang S-C, Hatem C, et al. In vivo organization of the FtsZ-ring by ZapA and ZapB revealed by quantitative super-resolution microscopy. Mol Microbiol. 2013;89: 1099–1120. pmid:23859153

69. Gueiros-Filho FJ, Losick R. A widely conserved bacterial cell division protein that promotes assembly of the tubulin-like protein FtsZ. Genes Dev. 2002;16: 2544–2556. pmid:12368265

70. Mercy C, Ducret A, Slager J, Lavergne J-P, Freton C, Nagarajan SN, et al. RocS drives chromosome segregation and nucleoid protection in Streptococcus pneumoniae. Nat Microbiol. 2019;4: 1661–1670. pmid:31182798

71. Nayak S, Pradhan D, Singh H, Reddy MS. Computational screening of potential drug targets for pathogens causing bacterial pneumonia. Microb Pathog. 2019;130: 271–282. pmid:30914386