Cheminformatics identification of phenolics as modulators of penicillin binding protein (PBP) 2x of Streptococcus pneumoniae towards interventive antibacterial therapy

Infections caused by multidrug-resistant Streptococcus pneumoniae remain the leading cause of pneumonia-related deaths in children < 5 years globally, and mutations in penicillin-binding protein (PBP) 2x have been identified as the major cause of resistance in the organism to beta-lactams. Thus, the development of new modulators with enhanced binding of PBP2x is highly encouraged. In this study, phenolics, due to their reported antibacterial activities, were screened against the active site of PBP2x using structure-based pharmacophore and molecular docking techniques, and the ability of the top-hit phenolics to inhibit the active and allosteric sites of PBP2x was refined through 120 ns molecular dynamic simulation. Except for gallocatechin gallate and lysidicichin, respectively, at the active and allosteric sites of PBP2x, the top-hit phenolics had higher negative binding free energy (ΔGbind) than amoxicillin [active site (−19.23 kcal/mol), allosteric site (−33.75 Kcal/mol)]. Although silicristin had the best broad-spectrum effects at the active (−38.41 kcal/mol) and allosteric (−50.54 kcal/mol) sites of PBP2x, the high thermodynamic entropy (4.90 Å) of the resulting complex might suggest the need for its possible structural refinement for enhanced potency. Interestingly, silicristin had a predicted synthetic feasibility score of < 5 and quantum calculations using the DFT B3LYP/6-31G+ (dp) revealed that silicristin is less stable and more reactive than amoxicillin. These findings point to the possible benefits of the top-hit phenolics, and most especially silicristin, in the direct and synergistic treatment of infections caused by S. pneumoniae. Accordingly, silicristin is currently the subject of further confirmatory in vitro research.


Introduction
The incidence of antibiotic-resistant Streptococcus pneumoniae has been on the increase since its first report in the 1960s [1].Globally, S. pneumoniae is the leading cause of community-acquired pneumonia, bacteremia, meningitis, and otitis media [2] and in children < 5 years of age, the pathogen remains the leading cause of pneumonia-related deaths [3][4][5].With more than half of S. pneumoniae isolates frequently showing reduced susceptibility to beta-lactam antibiotics [6,7], the most frequently administered antibiotics for pneumococcal infections, treatment of the diseases has become problematic, necessitating the need for the discovery and development of novel drug candidates with an alternative mode of binding of key therapeutic target of the organism [6].Penicillin-binding proteins (PBPs) are the typical target enzymes for beta-lactams and are engaged in the late phases of peptidoglycan biosynthesis, where their essential function involve transpeptidation of two muropeptides [8].Mutation of important PBPs is a resistance mechanism which has been described in various bacteria species [9,10].The mutation entails modifying key PBPs such that interaction with beta-lactams occurs at considerably higher antibiotic doses than with PBPs from susceptible strains, and therefore the biological activity of the antibiotics is significantly diminished [8].Thus, identifying improved modulators against critical PBPs of S. pneumoniae may aid in the battle against infections caused by multidrug resistant S. pneumoniae.S. pneumoniae possesses six PBPs, including PBP1a, PBP1b, and PB2a, which are Class A high-molecular-mass PBPs (hmm-PBPs), PBP2x and PBP2b (Class B hmm-PBPs), and PBP3, a D, D-carboxypeptidase low-molecularmass PBP (lmm-PBP) [11].Antibiotic resistance in S. pneumoniae has been linked to all six PBPs [9].However, only the class B hmm-PBPs (PBP2x and PBP2b) are key beta-lactam targets, and thus far, only PBP2x has been characterized as a major PBP target in clinical isolates since PBP2x genes that express low-affinity variants of S. pneumoniae can transfer beta-lactam resistance to susceptible strains [12].Also, resistance due to mutations in PBP 2b and 1a of S. pneumoniae has been shown to occur in strains that already have a low affinity PBP2x, suggesting PBP2x as the most essential PBP in S. pneumoniae that is responsible for resistance [12].Beta-lactams interact with PBP2x by acylating the active (actv) site serine residue located in the first of three conserved motifs (S*XXK, SXN, and KS/TG) of PBP2x, and mutations in this region have been attributed to resistance, which is the primary resistance mechanism for beta-lactam antibiotics in S. pneumoniae (Tsui et al., 2014).Interestingly, a PASTA domain located adjacent to the actv site (27 Å), whose function has been well-deciphered over the years, has in recent times been demonstrated to actv as an allosteric (allo) site in the generation of promising peptidoglycan and thus, when appropriately bound, could aid beta-lactam entry into the actv site [12,13].
Bernardo-Garca et al. [12] showed that the complexation of a pentapeptide strain at the PASTA domains prompts another pentapeptide strand in the same nascent peptidoglycan strand to enter the actv site for the turnover events.Thus, the PASTA domain of PBP2x could serve as a good therapeutic target for identifying drug candidates that could actv in synergy with beta-lactams [14], and in this study, both the actv site and the allo site of PBP2x were employed in the identification of drug candidates against infections caused by S. pneumoniae.
Plant and plant-derived metabolites have long been used in the treatment of human diseases.Interestinling the increased interest in plant-based antimicrobials in recent years, particularly phenolic compounds, has resulted in the identification of numerous phenolic candidates demonstrating promising pharmacological properties [15,16].
Phenolics ability to cause membrane damage, bind to virulence factors such as PBPs, and cause the suppression of bacteria biofilm formation, among other properties, has been implicated in its ability to inhibit both Grampositive and Gram-negative microorganisms [17,18].The characteristic hydroxyl groups attached to the benzene rings of phenolics have been shown to enhance strong protein binding while increasing affinity for cytoplasmic membranes [16].Alongside their direct antibacterial activity against microorganisms, several phenolics have been shown to work in synergy with conventional antibiotics, an approach that is encouraged in the treatment of infections caused by multi-drug-resistant bacteria [19].Several phenolics have been identified from various plant parts, with repositories in several databases [20].Repurposing phenolics could have the potential to reduce the time and costs to identify potent antimicrobials against PBP2x of S. pneumoniae [20].Thus, in this study, all the currently known phenolics were repurposed against PBP2x by compiling phenolics that could interact with the protein from the ZINCPharmer database using structural-based pharmacophore and molecular docking techniques.The limitations of virtual screening techniques relating to lack of energy refinement and assessment of stability were addressed using molecular dynamic simulation and mechanics/GB surface area (MMGBSA) approaches.This was carried out to identify potentially novel and actv antibacterial candidates that can treat or act in synergy with conventional antibiotics in the treatment of infections caused by S. pneumoniae.

Screening of phenolics at the actv site of PBP2x using structure-based pharmacophore and molecular docking
Following the evaluation of the binding positioning and interaction of the 1555 phenolics obtained via structurebased pharmacophore at the actv site of PBP2x using the molecular docking approach, the top 20 phenolics had higher docking scores than all the beta-lactam antibiotics used as standards.The top 20 phenolics had docking scores ranging between -8.4 and -9.9 kcal/mol, with ZINC03793048 having the highest docking score, while the beta-lactams had docking scores between -7.1 and -8.0 kcal/mol, with amoxicillin having the highest docking score (Figure 1), and hence it was selected as a standard for further analysis.This observation suggests the probable better interactions and orientation of the top-hit compounds at the actv site of PBP2x relative to the standards.However, as virtual screening often produces inactive molecules in the preclinical and clinical stages of drug development [21], the identified top 20 compounds were further screened cheminformatically to identify the top hit phenolics that will be pharmacokinetically-friendly, have drug-like properties, and have a synthetic score that can enable structural alterations for enhanced druggability (Figure 1).Using these approaches, the top hit phenolics obtained and the standards passed the Lipinski rule of 5 (LRo5), which serves as a pioneering physicochemical filter, relating the physicochemical parameters of compounds to their pharmacokinetic properties and oral bioavailability [22].The LRo5 limits the molecular weight, hydrogen donor, hydrogen acceptor, and octanol coefficient of a bioactive drug to 500 g/mol, 5, 10, and 5, respectively (Figure 1).Also, the top-hit phenolics and the standards had a synthetic score of < 5 (Figure 2) that would encourage structural alteration to improve potency [23].GCG (ZINC03870413) had the lowest synthetic score (4.20) among the top-hit phenolics, with a value that compares well with amoxicillin (4.17), suggesting that the compound will be the most susceptible to structural modification for enhanced druggability.Other predicted pharmacokinetic properties (Table 1), such as the bioavailability scores, solubility, and gastrointestinal tract absorption of the top-hit phenolics, all point to the drug-like properties of the compounds [24], as they all have characteristics that compare well with amoxicillin.
The top-hit phenolics toxicity profiles suggest that they are non-inhibitors of most of the cytochrome P450 (CYP), an important isoenzyme in drug metabolism [25].GCG and amoxicillin are non-inhibitors of all the isoenzymes, suggesting that they would not result in drug-drug interactions (DDI) when co-administered with phase I-and phase II-dependent drugs.However, other top hit phenolics inhibit CYP3A4, suggesting their potential ability to cause DDI with drugs that are metabolized by the isoenzyme.This observation pinpoints to the need for further structural alteration to improve the toxicity profile of the compounds, more so as some of the compounds [silicristin, lysidicichin, and ECMG] were also predicted to be active for immunotoxicity (Table 1).However, ECBE, GCG, and amoxicillin are predicted to be non-active for all the toxicity endpoints.Remarkably, all the top-hit phenolics belong to either toxicity class 4 or 5 and can be deployed as medications [22,24].Generally, the top-hit phenolics in this study have pharmacokinetic, synthetic, and drug-like properties that suggest them as potential drug candidates that can be structurally modified for enhanced potency.

Energy components refinement of the top-hit phenolics against the actv site of PBP2x
The molecular docking study has been shown to be a poor predictor of binding affinity and the identification of active inhibitors due to its scoring functions inability to account for the dynamics of the atoms of the ligand and the protein involved [21].Thus, in this study, energy refinement of the docking score was done using the MMGBSA approach [26].Using this approach, the ΔGbind, which gives the sum of all the intermolecular interactions that are present between a ligand and a target during a simulation, was calculated.Several energy components, such as the ΔE elec , ΔE gas , ΔE vdw and ΔG sol contributed to the total ΔGbind of a ligand to a receptor.
In this study, the ΔG sol , which estimates the the electrostatic solvation of the solute [27], had the highest contribution to the total ΔGbind observed for all the investigated complexes (Table 2).Three of the top-hit phenolics [ECMG (-42.18kcal/mol), silicristin (-38.41 kcal/mol), and lysidicichin (-29.85 kcal/mol)] had higher ΔGbind than amoxicillin, with ECMG having the highest negative value (Table 2).This observation, is not only consistent with previous studies on the antibacterial properties of phenolics [17], but indicative of the enhanced potential of the three phenolics as modulators of PPBP2x, with ECMG being the most promising compound.
ECBE (-19.72 kcal/mol) had ΔGbind scores that compared with amoxicillin (-19.23 kcal/mol), while the value observed for GCG (-14.93 kcal/mol) was significantly lower relative to amoxicillin (Table 2).This observation is suggestive of the lesser ability of ECBE and most especially GCG to modulate PBP2x of S. pneumoniae, thus contradicting the observed docking scores, where the two compounds had higher scores than amoxicillin.This finding substantiates the report of Cerón-Carrasco [21] that virtual screening required further energy refinement in the selection of active molecules.Overall, the findings from the thermodynamic ΔGbind of this study demonstrated the benefit of ECMG, silicristin, and lysidicichin over amoxicillin in the management/treatment of infections caused by multidrug resistant S. pneumoniae.To understand the structural stability of PBP2x and its complexes after ligand binding, the RMSD was assessed.

166
The RMSD is a function of the time-dependent deviation of a complex structure from its apo structure, with lower 167 RMSD values indicating better complex stability [28].In this study, the RMSD seems to equilibrate before 10 ns 168 for all the investigated systems, while each system subsequently took a more diverse route, which appears to 169 converge at 85 ns of simulation, thus impacting the different average RMSD values observed for each system 170 (Figure 2).The bound systems with average values ranging between 3.31 and 4.90 Å, generally had a reduced

171
RMSD value compared to the apo-PBP2x (6.32 Å) with lysidicichin and ECMG having the lowest values, which 172 compared well with amoxicillin (Table 3).The reduced RMSD value of PBP2x following binding of ligands 173 generally denotes enhanced thermodynamic stability of the complexes [29].This observation agrees with the 174 report of Sainsbury et al. [30], where experimental beta-lactam binding of PBP3 increases the thermodynamic  191 The ROG estimates the degree of folding with great folding meaning better thermodynamic organization and 192 sometimes also stability in a complex, a phenomenon related to protein (topology, size, and length) and ligand (molecular weight) parameters and their point of contact [31].Like the RMSD, following ligand binding of PBP2x, a marginal reduction in the ROG values was observed with the bound systems relative to the apo-PBP2x (Table 3, Figure 3).The reduced ROG value following ligand binding of PBP2x could be suggestive of the increased structural compactness in the bound systems due to folding [32].The bound system had average ROG values ranging between 29.15 and 30.01 Å with lysidicichin having the lowest value (Table 3).Lysidicichin having the lowest ROG value agrees with the RMSD findings of this study, where the compound also had the lowest RMSD value, and thus suggests the advantage of the compound as a potential PBP2x inhibitor.However, with a variance of 0.86 Å (Table 3) among the bound systems regarding the average ROG, it can be inferred that the compactness of PBP2x following ligand binding had little impact on ranking the top-hit phenolics.Root mean square fluctuation evaluates the mobility of residues following ligand binding of PBP2x [33].In this study, higher swaying in amino acid residues between 25 and 175 and a lesser fluctuation between residues 250 and 425, and 550 and 625 were noted in all the systems, with the degree of swaying varying among systems (Figure 4a), thus affecting the average RMSF.Interestingly, Ser337, which is the catalytic residue in PBP2x [10], occurred where the fluctuation was minimal (Figure 4b).Ser337 had fluctuations ranging between 0.87 Å and 1.29 Å in all the system with the apo-PBP2x having the highest fluctuation of Ser337 (Table 4).This observation suggests that ligand binding generally reduced the fluctuation of Ser337, indicating an active involvement of Ser337 in the intramolecular binding of other residues of PBP2x or the intermolecular binding of the ligands.The reduced effect of ligand binding on the Ser337 of PBP2x is in agreement with our previous report, where the gatekeeper residue (Tyr446) of PBP2a had reduced fluctuation following binding of phenolics [34].Silicristin (0.87 Å) and ECMG (0.88 Å) with the highest negative ΔGbind against PBP2x also had the lowest fluctuation of Ser337, pinpointing the importance of reduced fluctuation of the catalytic residue in the observed ΔGbind in the study.However, except for the ECMG-PBP2x complex (1.64 Å) , the fluctuation of residues in the whole protein was greatly enhanced following the top-hit phenolics and amoxicillin binding of PBP2x (Table 3).The apo-PBP2x (2.16 Å) (Table 3), suggesting the advantage of ECBE, lysidicichin, and GCG as potential inhibitors of PBP2x.
Overall, the RMSF findings of this study have demonstrated the importance of reduced fluctuation of Ser337 in the observed ΔGbind noted in this study and that ECMG among the investigated ligands caused enhanced structural organization in PBP2x, which could have contributed to its higher negative ΔGbind.Silicristin, on the other hand, despite having a high ΔGbind, might require further structural adjustment to increase its potency due to its high RMSF value, and this observation is also noted by the RMSD finding of this study.
12 239  Like the ROG, the SASA also analyzes protein folding and surface area changes during a simulation, with higher

246
SASA values indicating a rise in volume [35].The SASA plot had a more stable fluctuation throughout the 247 simulation for all the systems (Figure 5).Averagely, the bound system had slightly higher SASA values ranging 248 between 28065.41Å and 28612.65 Å compared to the apo-PBP2x, with lysidicichin having the lowest value

249
(Table 3).This observation highlighted the slightly higher surface area of PBP2x following ligand binding, which 250 is in contrast with the findings of Mousavi et al. [33], who reported that the interactions of ligands with proteins 251 can cause protein folding.The increase in surface area following ligand binding also did not agree with the ROG findings of this study.This observation might not necessarily mean less stability of the complexes, as a smaller SASA variance of 800 Å exists between the bound system and the apo-PBP2x (Table 3).However, lysidicichin, among the investigated compounds, has the lowest SASA value, which agrees with the ROG and RMSD findings of this study and thus further corroborates the advantage of the compound as PBP2x inhibitor.Overall, with a smaller SASA variance of 547 Å among the bound systems, it can be concluded that the surface area folding of PBP2x following ligand binding had little impact on ranking the top-hit phenolics.Hydrogen bonds and distances are important for the structural stability of proteins and can hence be evaluated during simulation to understand the influence of ligand binding on protein stability [36].A stable fluctuation in the number of intramolecular hydrogen bonds was observed with bonded and apo-PBP2x during the simulation (Figure 6a).This observation, while being consistent with other thermodynamic metrics of this study, further indicates the thermodynamic orderliness of the systems [37].Following ligand binding of PBP2x, a higher number of hydrogen bonds in PBP2x was observed with the bound system (ranging between 329.78 and 333.98) relative to the apo-PBP2x (Table 3), with ECMG having the highest number of hydrogen bonds.The increase in the hydrogen bond of the bound system relative to the apo-PBP2x is indicative of some intermolecular hydrogen bonds contributed by ligand binding of the protein, with ECMG having the highest intermolecular contribution of hydrogen bonds during the simulation.This observation could have contributed to ECMG having the highest negative ΔGbind against PBP2x and the higher negative ΔGbind observed with silicristin, ECBE, lysidicichin, and ECMG relative to amoxicillin.The relatively similar pattern of fluctuations in the distance of hydrogen bonds of PBP2x [with an average of 2.85 Å for each of the systems (Table 3)] following the top-hit phenolics binding of PBP2x, which continues to shorten as the simulation progresses (Figure 6b) also further emphasized the orderliness in the protein following the compounds binding at the actv site, an observation that agrees with other thermodynamic metrics of this study.

Molecular docking of top hit phenolics at the allo site of PBP2x of S. pneumoniae
The PASTA domain of PBP2x has been demonstrated to act as an allo site, which, when properly occupied, could facilitate entering the actv site of PBP2x and thus prompt the acylation of Ser337 of the actv site [12].To investigate the ability of the top-ranking phenolics binding the allo site of PBP2x, molecular docking was carried out to assess the fitness of the phenolics at the site.The top-hit phenolics had higher negative docking scores at the allo site, ranging between -9.1 and -9.9 kcal/mol, than amoxicillin (-8.8 kcal/mol), with silicristin having the highest negative docking scores (Table 5).This observation might be indicative of the better fitness of the top-hit phenolics relative to amoxicillin at the allo site of PBP2x, with silicristin showing the better fitness.The docking scores observed at the allo site are comparable with those observed at the actv site, which also range between -9.2 and -9.9 kcal/mol, with silicristin having the highest negative docking scores at both sites.However, as molecular docking approach only serves as a screening process for obtaining the pose, orientation, and interactions of a ligand at the binding pocket of a protein, MMGBSA and molecular dynamic studies were further carried out to assess the ΔGbind, thermodynamic stability, and fluctuation of the actv site Ser337 following binding of the PBP2x allo site by the top-hit phenolics.5, Table S2).This observation is suggestive of the better affinity of the top-hit phenolics for the allo site of PBP2x, with silicristin being the most promising as a probable inhibitor at the site, and interestingly, this finding agrees with the docking findings of this study.The observed high affinity of silicristin at both the actv and allo sites of PBP2x points to the probable broad-spectrum activity of the compound as an inhibitor that can treat or act in synergy with conventional antibiotics in the treatment of infections caused by S. pneumoniae.This observation is consistent with previous studies where phenolics show good synergistic properties with conventional antibiotics [14 19].ECBE and GCG also had good ΔGbinds at the allo site of PBP2x, indicating their probable capability to act in synergy with conventional beta-lactams ( moderately higher negative ΔGbinds than amoxicillin (Table 5).Analysis of the thermodynamic stability suggests stable trajectories of the top-hit phenolics at the allo site of PBP2x with RMSD values that are less than 3 Å limit [38] (ranging between 1.74 and 2.19 Å) and negligible ROG and SASA variances of 0.23 Å and 1042 Å respectively, between the bonded systems and the apo-PBP2x (Figure S2, Table S3).These observations are revealing of the thermodynamic compatibility of the top-hit phenolics at the allo site of PBP2x.While two different experimental structures of PBP2x were used for the actv and allo site evaluations, the low RMSD values and smaller ROG variance following binding at the allo site relative to the actv site could be pinpointing the higher thermodynamic stability of the top-hit phenolics at the allo site.

Thermodynamic fluctuation of the actv site Ser337 following simulation of the top-ranked phenolics at the allo site of PBP2x of S. pneumoniae
The RMSF fluctuation of PBP2x following the top-hit phenolics and amoxicillin binding of PBP2x ranges between 9.61 and 15.79 Å with silicristin having the lowest value (Table 5, Figure 7a).Unlike the actv site, this observation implies a significantly higher fluctuation of the residues of PBP2x following the top-hit phenolics and amoxicillin binding at the allo site.However, it is worth acknowledging that the different experimental crystal structures of PBP2x used in this study for the actv and allo sites might have played a major role in the differences observed.Substantiating this observation is the fact that except for ECMG (15 Å), the top-hit phenolics and amoxicillin binding of PBP2x all resulted in a reduced fluctuation of residues relative to the apo-PBP2x (Table 5, Figure 7a).This observation implies that the residues were able to form stable intra-and intermolecular bonds [36] during the MD simulation of the investigated compounds at the allo site of PBP2x.Silicristin having the lowest RMSF value at 7.86 Å (Table 5) suggests the advantage of the compound as a probable inhibitor of the site, an observation that is consistent with the thermodynamic ΔGbind findings of this study.However, it is surprising that while silicristin had the lowest RMSF value in the whole protein, the compound instigated the highest fluctuation of Ser337 (Table 5, Figure 7b), a critical residue at the actv site PBP2x [10].The highest fluctuation of Ser337, in silicristin-PBP2x complex is coherent with the highest negative ΔGbind of silicristin (-50.54 kcal/mol) at the allo site of PBP2x.Similarly, ECBE (-45.82 kcal/mol) and GCG (-44.21 kcal/mol), with the second and third highest negative ΔGbind at the allo site, had higher fluctuations of Ser337 than other investigated compounds at 16.59 Å and 16.60 Å respectively (Table 5, Figure 7b).While other investigated compounds with moderately high negative ΔGbind at the allo site also have a higher fluctuation of Ser337 than apo-PBP2x (Table 5, Figure 7b).All these observations are suggestive of the high instability of Ser337 of the actv site following ligand binding of the allo site, with the ΔGbind of the ligands at the allo site partially correlating the fluctuation of the critical amino acid.Another conclusion from this is that there is less involvement of the critical amino acid in the intramolecular binding of other important residues at the actv site, which might mean less obstruction and thus encourage the entry of beta-lactam.This observation is consistent with the earlier report of Bernardo-Garca et al. [12], who demonstrated that the complexation of a pentapeptide strain at the PASTA domains, prompts another pentapeptide strand in the same nascent peptidoglycan strand to enter the actv site of PBP2x.Another study by Mahasenan et al. [39] also showed experimentally that the displacement of adjacent residues that can enter intramolecular interactions with Tyr446 (the critical amino acid in PBP2a) after allosteric binding of PBP2a with cefepime, which has also been demonstrated computationally by Aribisala and Sabiu [34], that appropriate occupancy of the allo site of PBP2a by phenolics caused high fluctuation of Tyr446.Thus, the high fluctuation of Ser337 following ligand binding to the allo site of PBP2x is consistent with previous report.
This can be further substantiated when the fluctuation of the critical amino acid is compared to the fluctuations when the top-hit phenolic and amoxicillin bind at the actv site.Binding of the top-hit phenolics at the actv site caused significant lesser fluctuation of Ser337 (between 0.87 Å and 1.19 Å) relative to the apo-PBP2x (1.27 Å) (Table 4, Figure 4b) while binding at the allo site caused significant higher fluctuation (between 14.03 and 18.75 Å) relative to the apo-PBP2x (11.44 Å) (Table 5, Figure 7b).This observation is revealing of the inter-or intramolecular involvement of Ser337 following ligand binding at the actv site, while binding at the allo site, on the other hand, causes lesser involvement of the residue in the intramolecular binding or obstruction of the actv site.pneumoniae over a 120 ns simulation period.

Interactions plots analysis of the top-hit phenolics at the actv and allo site of PBP2x of S. pneumoniae after 120 ns MD simulation
To understand how the nature of interactions impacted the ΔGbind and stability of the top-hit phenolics against PBP2x, the plots of interactions for each system after the 120 ns simulation were analyzed.The binding of drugs to biological targets has been demonstrated to rely on specific non-covalent interactions between partner molecules, including electrostatic interactions such as hydrogen bonds, van der Waals interactions, and hydrophobic effects [40].In this study, the number, nature, and lengths of interactions between partner molecules (ligands and PBP2x) were noted to have an impact on the ΔGbind and stability observed in this study.The plots of interactions of the reference standards and phenolics with the highest ΔGbind at both the actv and allo sites of  which had the highest negative ΔGbind against PBP2x interacted with the highest number of interactions at 15, three of which are hydrogen bond contacts with the lowest average distance at 2.84 Å with residues Glu295, Asn294, and Asp472 (Table 6, Figure 8).The high number of interactions and the significantly lower length of the hydrogen bond contacts could have drastically contributed to the observed high negative ΔGbind of ECMG against PBP2x.Silicristin, which has the second highest negative ΔGbind against PBP2x, had the highest number of interactions at 16, five of which are hydrogen bond contacts with an average length of 3.16 Å with residues Ser465, Tyr503, Glu295, Asn294 and Ser254 (Table 6, Figure S3).Compared to ECMG-PBP2x, even though silicristin-PBP2x had the highest number of interactions, including hydrogen bond contacts, the observed higher average distance of hydrogen bond contact in the silicristin-PBP2x complex could have contributed to its lower negative ΔGbind and stability relative to ECMG-PBP2x complex.The importance of hydrogen bond length observed in this study is consistent with the findings of Du et al. [36], who showed the relevance of reduced interaction length in the higher pull between two intramolecular or intermolecular molecules.Reduced bond length will facilitate the atoms of two molecules being held together more tightly, resulting in greater affinity and stability [36].Similarly, lysidicichin, with a moderately higher negative ΔGbind relative to amoxicillin (Table 6, Figure 9), and other top hit phenolics had 10 interactions, 4 of which were hydrogen bond contacts with an average distance of 2.98 Å with Ser504 (2), Gln499, and Thr443 (Table 6, Figure S3).The higher number and reduced average distance of hydrogen bond contacts of lysidicichin against PBP2x could have impacted the higher negative ΔGbind relative to amoxicillin-PBP2x, ECBE-PBP2x, and GCG-PBP2x, thus further highlighting the importance of the number and bond length in ΔGbind determination.ECBE and GCG against PBP2x only had 3 and 4 total interactions with PBP2x with no hydrogen bond contacts (Table 6, Figure S3).This observation could have immensely contributed to the significantly lower ΔGbind observed with the ECBE-PBP2x and GCG-PBP2x relative to other studied complexes in this study.At the allo site, a considerable similar number of interactions, ranging from 13 to 16 total interactions, were observed between the investigated compounds and the allo site of PBP2x, with GCG having the highest interaction (Table 6).Silicristin (Figure 10), ECBE (Figure S4), and amoxicillin (Figure 11) had the same number of interactions at 15.The high number of interactions in the silicristin, GCG, and ECBE complex with PBP2x could have impacted the higher negative ΔGbind of the compounds at the allo site of PBP2x.Despite having a similar number of total interactions with silicristin and ECBE, amoxicillin has a considerably lower negative ΔGbind, which could be partly attributable to the reduced number of hydrogen bond contacts in the amoxicillin-PBP2x complex (Figure 11) relative to silicristin (Figure 10) and ECBE (Table 6, Figure S4).The catalytic residues (such as Ser337, Tyr568, Pro335, Ser396, Trp374 and Ser395 at the actv site and Arg426, Gly422, Phe421, Asn417, and Val423 at the allo site) present in the molecular docking interactions of the investigated compounds with PBP2x were not observed in the interaction plots of both sites after the 120 ns MD simulation.Similar observation was noted by Al-Karmalawy et al. [41], following simulation of alacepril with hACE2 and might have been due to the dynamic nature of simulation employed.
Nevertheless, the lesser fluctuation of Ser337 (the catalytic residue in PBP2x) compared to the RMSF of the whole protein means that the top-ranked phenolics and amoxicillin during the simulation reside most at the actv region of the protein, as lesser fluctuation is tantamount to strong intermolecular bonding [42].

Molecular orbital properties of best performing phenolics at the actv and allo sites of PBP2x
Understanding the HOMO and LUMO characteristics of a compound could provide crucial chemical insights relevant to the actual use of the compound as a medicinal molecule [43,44].The energy gap (LUMO-HOMO) of a compound provides information about its stability; A larger energy gap means high stability, which can prevent ligand binding to a protein [43,44].In this study, silicristin (4.0510 eV) and ECMG (4.5388 eV) had lower energy gaps than amoxicillin (5.3234 eV) (Table 7).This observation points to the high instability of silicristin and ECMG relative to amoxicillin.The higher softness and lower hardness of silicristin and ECMG relative to amoxicillin (Table 7), might suggest better reactivity of silicristin and ECMG compounds compared to amoxicillin. .According to Pearson [45], higher chemical softness and lower chemical hardness signify higher compound reactivity.The high instability and better reactivity of silicristin and ECMG relative to amoxicillin could have positively impacted their better interactions at the actv and allo sites of PBP2x.Furthermore, the ability of an atom or a functional group of a compound to attract electrons could be assessed by evaluating the electronegativity and chemical potential.In this study, the higher electronegativity and higher negative chemical potential of silicristin and ECMG relative to amoxicillin (Table 7) means better capability of the atoms and functional groups of silicristin and ECMG to attracts electron at the binding pockets of PBP2x and hence their better affinity for PBP2x than amoxicillin in this study.This observation is also further supported by the higher electrophilicity index of silicristin and ECMG compared to amoxicillin (Table 7), which indicates the better tendency of these compound to acquire electron density relative to amoxicillin [46].3.0210

Conclusion
The findings from this study revealed the diverse capability of phenolics in the binding of the actv and allo sites of PBP2x, with ECMG, silicristin, and lysidicichin having advantages in the binding of the PBP2x actv site.On the other hand, ECBE, GCG and silicristin had better affinity at the allo site.The higher affinity of silicristin at both the actv and allo sites of PBP2x showed that the compound had the best broad-spectrum effects and thus had an advantage over other top hit phenolics and amoxicillin as inhibitors of PBP2x.However, relative to other top hit phenolics, silicristin elicited higher thermodynamic entropy following binding at the actv site.This discovery might point to the necessity for further structural modifications to the molecule for stability, which could greatly boost its affinity for PBP2x.Allosteric modulation of Ser337, the critical amino acid at the actv site of PBP2x, was noted as ECBE, GCG, and silicristin with higher negative ΔGbind at the allo site instigated more fluctuation of Ser337 than amoxicillin, lysidicichin, and ECMG with moderately high affinity.This observation could imply that higher affinity of the top-hits at the allo site of PBP2x causes less obstruction of the actv site due to the lesser involvement of Ser337 in intramolecular binding of adjacent residues.Comparably, binding of the top hit phenolics and amoxicillin at the actv site of PBP2x causes less fluctuation of Ser337, suggesting that the residue could be involved in inter-and intra-molecular interaction.These findings hint to the potential advantages of tophit phenolics, particularly silicristin, in the direct and synergistic treatment of S. pneumoniae infections with conventional beta-lactam medicines.More importantly, the top-hits phenolics had drug-like pharmacokinetics and toxicity profile that potentiate them as probable drug candidates.However, further in vitro and in vivo validation of the activities of the top-hits phenolics is highly suggested, and work is ongoing in this regard.

1. PBP2x acquisition, optimization, and identification of actv and allo sites
The experimental crystal structures of PBP2x from S. pneumoniae (1QMF for the actv site evaluation and 5OIZ for the allo site) were obtained from the protein data bank (PDB) (https://www.rcsb.org)[10,12].The crystal structures were optimized by eliminating water molecules and nonstandard amino acids using the UCSF Chimera v1.15 software tool [47].Afterword, the x-y-z coordinates and amino acid residues at the actv {center [x (116.[48] and thereafter validated from literature [10,12].demonstrated that the top-hit phenolics, amoxicillin, and the docked native inhibitors of 1QMF (actv) and 5OIZ (allosteric) had relative binding orientation at both sites of the protein with RMSD values of 1.5 Å (actv site) and 1.4 Å (allo site).Furthermore, all the ligands interacted with the key amino acid residue in the actv (Ser337) (Figure 12a) [10] and allo (Arg426) (Figure 12b) [12] sites of PBP2x.These findings supported the validity of the docking scores obtained in the study.

Pharmacokinetic and physicochemical prediction of phenolics
Using the simplified molecular input line entry system (SMILES)of each structure, the pharmacokinetics, physicochemical, and drug-like properties of the top twenty phenolics were predicted using the SwissADME web tools (http://swissadme.ch/index.php)and afterwards validated using the Molinspiration toolkit (https://www.molinspiration.com/cgi-bin/properties).On the other hand, the Protox II webserver (https://toxnew.charite.de/protox_II/)was employed in the prediction of the toxicity properties of the compounds.Coupled with the docking scores, the cheminformatics information obtained from these analyses were used in ranking the top hit phenolics.

Chemical fingerprinting of phenolics
their plots were made using Origin v 6.0.Similarly, the binding free energy (ΔGbind) was determined using an average of 120,000 snapshots collected from a 120 ns MD simulation trajectory using the equations below, utilizing the MMGBSA technique [26] ΔGbind Where Egas = the gas-phase energy, E int = the internal energy, E ele = coulomb energy E vdw = the van der Waals energy, Gsol = solvation-free energy from polar state, G GB = solvation−free energy from polar state non−polar states, S = total entropy and T= temperature.
The interaction plots of the top hit phenolics and amoxicillin against PBP2x at both actv and allo sites of the protein were visualized and analyzed with Liplot v2.2.8 (Wallace et al., 1995).

Quantum chemical calculations
The purpose of this investigation was to predict the molecular characteristics of the best phenolics at the actv and allo sites of PBP2x.The most generally used DFT/B3LYP/631G/+ (d, p) basis set [43,44]

Figure 1 :
Figure 1: Cheminformatics and docking scores ranking of the top-hit phenolics.The same colour and cluster signifies more similarity.Except for ZINC72320533, all top twenty phenolics had the same colour meaning they are more similar than the beta-lactams with different colour.Compounds of the same clusters with zero similarity scores are conformers and in ranking the top-hit, only one (compounds that fulfilled the Lo5 and has the highest docking score) from a set of conformers are chosen.The top-hit phenolics all had a synthetic sore of < 5 and possess the resorcinol structure as a common moiety.

Figure 2 :
Figure 2: Relative root means square deviation plots of alpha-carbon, amoxicillin, and top-hit phenolics at the

Figure 3 :
Figure 3: Relative radius of gyration plots of alpha−carbon, amoxicillin and top-hit phenolics at the actv site of PBP2x of S. pneumoniae over a 120 ns simulation period

Figure 4 :
Figure 4: Relative root means square fluctuation plots of alpha−carbon, amoxicillin and top-hit phenolics at the

Figure 5 :
Figure 5: Relative solvent−accessible surface area plots of alpha−carbon, amoxicillin and top-hit phenolics at the actv site of PBP2x of S. pneumoniae over a 120 ns simulation period

Figure 6 :
Figure 6: Time progression of (a) number of H-bonds, and (b) distance formed in PBP2x following ligand binding at the actv site during the simulation.

Figure 7 :
Figure 7: Relative root means square fluctuation plots of alpha−carbon, amoxicillin, and top-hit phenolics at the allo site of PBP2x (a) and the critical residue at the actv site of PBP2x following binding at the allo site (b) of S.

PBP2x are presented in Figures 8 -
11, while those of the other top-hit phenolics are shown in FiguresS3 and S4 for the actv and allo sites, respectively.At the actv site, except for silicristin, epicatechin 3-O-(3-O-methylgallate),

Figure 12 :
Figure 12: Docking protocol validation via superimposition of the docked top-hit phenolics with the experimental native inhibitor at the actv (a) and allosteric (b) sites of 1QMF and 5OIZ respectively.The docked experimental inhibitor (blue), reference standard (green) and top-hit phenolics [silicristin (magneta), ECBE (cyan), lysidicichin (red), ECMG (orange) and GCG (yellow)] all had an RMSD of 1.5 Å and 1.4 Å from the experimental native inhibitor at the actv and allo site of 1QMF and 5OIZ, respectively, suggesting partial binding orientation of the docked compounds and the experimental native inhibitors.

Table 2 .
Thermodynamic Energy (kcal/mol) of the top-five phenolics against the actv site of PBP2x of S.
Keyword: ΔE vdW means van der Waals energy, ΔGbind means total binding free energy, ΔE gas means gas phase 162 free energy, ΔE elec means electrostatic energy, while ΔG solv means solvation free energy 163 164 2.3 Thermodynamic information of PBP2x following binding at the actv site 165

Table 3 .
Thermodynamics data of the top-five phenolics after 120 ns MD simulation at the actv site of PBP2x of 175 stability of PBP3 through local conformational changes at the actv site.The fact that lysidicichin and ECMG had 176 RMSD values that compared well with amoxicillin following binding of PBP2x indicates their advantage as 177 inhibitors of PBP2x over the other top-hit phenolics.This observation partially agrees with the ΔGbind findings178of this study, where ECMG had the highest negative ΔGbind.However, silicristin, which also had higher ΔGbind 179 against PBP2x, had the highest RMSD value at 4.90 Å (Table3) and thus denotes that further structural adjustment 180 of the compound for stability might improve potency.Overall, the RMSD finding of the study showed that ligand 181 binding of PBP2x increases stability when compared to the apo-PBP2x and that while lysidicichin and ECMG are had an RMSF of 1.76 Å, while the bound systems had average RMSF values ranging between 1.64 Å and 2.62 Å with only ECMG-PBP2x having a lesser value than the apo-PBP2x.The lesser RMSF value of ECMG-PBP2x suggests enhanced structural organization with the binding domain of PBP2x due to strong intra and intermolecular binding and could have contributed to the higher negative ΔGbind observed in ECMG-PBP2x relative to other investigated complexes in this study.Like the findings on the RMSD, the silicritin-PBP2x complex had the highest fluctuation of residues in the whole protein, suggesting the need for structural adjustment of the compound for stability.Other top-hit phenolic complexes with PBP2x, such as ECBE-PBP2x (1.83 Å), lysidicichin-PBP2x (1.76 Å), GCG-PBP2x (1.85 Å) had relatively lower RMSF values than amoxicillin-PBP2x

Table 4 :
PBP2x actv site Ser337 fluctuation (Å) after 120 ns simulation of the top-hit phenolics at the actv site

Table 5 :
Molecular docking, thermodynamics ΔGbind, whole RMSF and RMSF values of Ser337 following binding of top hit phenolics at the allo site of PBP2x of S. pneumoniae.

Table 5
). ECMG and lysidicichin, on the other hand, had

Table 6 :
Interactions plots analysis of the top-hit phenolics at the actv and allo site of PBP2x of S. pneumoniae after 120 ns MD simulation.

Table 7 :
The CDFT of the best performing phenolics at the actv and allo site of PBP2x using DFT calculated by