Investigation of Plant Antimicrobial Peptides against Selected Pathogenic Bacterial Species Using a Peptide-Protein Docking Approach

Mustafa, Ghulam; Mehmood, Rizwan; Mahrosh, Hafiza Salaha; Mehmood, Khalid; Ahmed, Shakeel

doi:https://doi.org/10.1155/2022/1077814

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results and Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 1077814 | https://doi.org/10.1155/2022/1077814

Investigation of Plant Antimicrobial Peptides against Selected Pathogenic Bacterial Species Using a Peptide-Protein Docking Approach

Ghulam Mustafa,¹Rizwan Mehmood,¹Hafiza Salaha Mahrosh,²Khalid Mehmood,³and Shakeel Ahmed^4,5

Academic Editor: Abdelwahab Omri

Received19 May 2021

Revised25 Nov 2021

Accepted01 Mar 2022

Published21 Mar 2022

Abstract

Antimicrobial resistance is the key threat to global health due to high morbidity and mortality. The alteration of bacterial proteins, enzymatic degradation, and change of membrane permeability towards antimicrobial agents are the key mechanisms of antimicrobial resistance. Based on the current condition, there is an urgent clinical need to develop new drugs to treat these bacterial infections. In the current study, the binding patterns of selected antimicrobial peptides (AMPs) with different multidrug-resistant bacterial strains have been analyzed. Among ten selected AMPs in this study, napin and snakin-1 exhibited the best scores and binding patterns. Napin exhibited strong interactions with penicillin-binding protein 1a of Acinetobacter baumannii (with a binding score of -158.7 kcal/mol and ten hydrogen bonds), with glucose-1-phosphate thymidylyltransferase of Mycobacterium tuberculosis H37Rv (with a binding score of -107.8 kcal/mol and twelve hydrogen bonds), and with streptomycin 3-adenylyltransferase protein of Salmonella enterica (with a binding score of -84.2 kcal/mol and four hydrogen bonds). Similarly, snakin-1 showed strong interactions with oxygen-insensitive NADPH nitroreductase of Helicobacter pylori (with a binding score of -105.0 kcal/mol and thirteen hydrogen bonds) and with penicillin-binding protein 2a of methicillin-resistant Staphylococcus aureus (with a binding score of -103.8 kcal/mol and twenty-three hydrogen bonds). The docking results were further validated by molecular dynamics simulations. The results of this computational approach support the evidence of efficiency of these AMPs as potent inhibitors of these specific proteins of bacterial strains. However, further validations are required to fully evaluate the potential of selected AMPs as drug candidates against these resistant bacterial strains.

1. Introduction

The excessive or inappropriate use of antibiotics has pushed the world towards the postantibiotic era. The bacterial strains mastered their own antidote that has led to severe resistance against many antimicrobial agents [1]. As the multidrug-resistant strains evolve, the invasive infections in hospitals and communities increase day by day in a complex pattern. Many antibiotics are failing now because of the occurrence of bacterial resistance due to mutational changes in the bacterial cellular machinery. Therefore, there is a dire need for potent antimicrobial agents which are less toxic and more effective. From the past few years, antimicrobial peptides (AMPs) have drawn much attention due to their vast therapeutic properties with fewer side effects [2].

AMPs are the first line of defense in plant innate immunity to protect them against microbial infections [3]. In contrast, some bacterial species also produce AMPs to counter other bacterial species in the amensalism relationship to compete and kill other bacterial strains for the same ecological niche [4]. About 17 families of AMPs have been reported with different antibacterial, antifungal, and antiviral activities [5]. These include defensin (PR-12 family), hevein-like peptide, thionin (PR-13 family), knottin, α-hairpinin, lipid transfer protein (LTP; PR-14 family), and snakin [6].

AMPs are classified as broad-spectrum antimicrobial agents that can regulate the innate immune system of various organisms such as bacteria, protozoa, fungi, plants, insects, and animals [7]. AMPs exhibit considerable structural and functional diversity that upholds their antimicrobial activity, microbial cell selectivity, and immunomodulatory properties which make them suitable drug candidates for the development of new therapies [8]. After microbial invasion into the host, the AMPs undergo genetic modifications for their expression or rapid transcription [9]. Being conserved in nature, these natural AMPs show specific resistance against certain pathogens including bacteria, pathogenic fungi, viruses, and parasites that invade into the metabolic machinery of the cell [10]. AMPs have been reported from almost all organisms, but more specifically, plant-derived antimicrobial peptides have been reported with structural and functional diversity [11]. Different studies have highlighted the physiological and therapeutic importance of various AMPs against multidrug-resistant pathogens [12]. The antimicrobial activity of AMPs makes them interestingly special as potential drug candidates against various pathogens [2].

Considering the current increase of antibiotic resistance in bacterial species because of their evolution, the discovery of new and effective natural compounds which could be employed in the treatment of various bacterial infections with fewer or no side effects compared to present antibiotics has become important in order to ensure the health of our future generations [13]. The current century witnesses an extraordinary advancement in the field of drug discovery because of the advent of in silico approaches of bioinformatics, molecular docking, and molecular dynamics simulation [14]. The advantage of the docking approach over traditional drug discovery is the predictions of protein pairs which enhance our knowledge about biological pathways by examining protein-peptide, protein-protein, or protein-ligand complexes which provide insights into the mechanisms of novel interactions [15].

The bacterial strains A. baumannii, M. tuberculosis, H. pylori, MRSA, and S. enterica used in this study have been reported by the World Health Organization (WHO) as the most deadly bacteria due to their multidrug resistance [16]. Ten efficient AMPs were selected from the literature and docked against five selected deadly bacterial strains. The purpose of this study was to explore the antibacterial activities of the most reported AMPs against five selected multidrug-resistant bacterial strains. The reported AMPs in this study are assumed to be useful for drug discovery professionals to check their bactericidal and/or bacteriostatic potentials in controlling these bacterial species. To provide a more meaningful in vivo prediction of efficacy of these AMPs, it is also necessary to combine the results and information of this study with pharmacokinetic and pharmacodynamic data because the clinical outcome must be the ultimate guide for curing any infection.

2. Materials and Methods

2.1. Selection and Retrieval of Antimicrobial Peptides

On the basis of reported antimicrobial activities in the literature, ten AMPs were selected for analyses against deadly pathogenic bacterial species. The 3D structures of selected peptides were downloaded from Protein Data Bank (https://www.rcsb.org/) in .pdb format. The selected AMPs and their PDB IDs were as follows: napin (PDB ID: 1PNB), snakin-1 (PDB ID: 5E5Q), knot1 domain-containing protein (PDB ID: 7C31), Amaranthus caudatus-AMP2 (PDB ID: 1MMC), EcAMP1 (PDB ID: 2L2R), Nigellin-1.1 (PDB ID: 2NB2), plant defensin NsD7 (PDB ID: 5KK4), flower-specific gamma-thionin (PDB ID: 6DMZ), acyclotide ribe 31 (PDB ID: 7KPD), and antimicrobial peptide 1a (PDB ID: 2LB7).

2.2. Retrieval of Receptor Proteins of Pathogenic Bacterial Species

For the current study, the five most virulent bacterial strains were selected. The 3D structures of receptor proteins of selected bacterial strains were downloaded from PDB in .pdb format. The selected proteins with their PDB IDs used in this study as receptor proteins were as follows: penicillin-binding protein 1a of Acinetobacter baumannii (PDB ID: 3UDF), glucose-1-phosphate thymidylyltransferase of Mycobacterium tuberculosis (PDB ID: 6B5E), oxygen-insensitive NADPH nitroreductase of Helicobacter pylori (PDB ID: 3QDL), penicillin-binding protein 2a of methicillin-resistant Staphylococcus aureus (PDB ID: 1MWT), and streptomycin 3-adenylyltransferase of Salmonella enterica (PDB ID: 6FZB).

2.3. Protein-Protein Docking

The hosts produce specific immune responses upon any pathogenic attack. Therefore, to observe the binding patterns and interactions between selected AMPs and specific bacterial receptor proteins, HADDOCK v.2.4 was used to carry out the docking analysis [17]. The educational version of the PyMOL Molecular Graphics System was used to predict the active site of each receptor protein and to visualize and draw interactions between AMPs and receptor proteins [18]. Besides, PDBsum was used to validate the interactions of amino acid residues involved in the docked complex [19]. Chain A of each receptor protein was used in docking studies.

2.4. Molecular Dynamics Simulation

Molecular dynamics (MD) simulations were performed for 100 nanoseconds using Desmond, a package of Schrödinger LLC. The initial stages of protein and peptide complexes for molecular dynamics simulation were obtained from docking studies. Molecular docking studies provide a prediction of binding status in static conditions. Simulations were carried out to predict the binding status in the physiological environment. The protein-peptide complexes were preprocessed using Protein Preparation Wizard or Maestro, which also include optimization and minimization of complexes. All systems were prepared by the System Builder tool. A solvent model with an orthorhombic box was selected as TIP3P (Transferable Intermolecular Interaction Potential 3 Points). The OPLS_2005 force field was used in the simulation. The models were made neutral by adding counter ions where needed. To mimic the physiological conditions, 0.15 M salt (NaCl) was added. The NPT ensemble with 300 K temperature and 1 atm pressure was selected for complete simulation. The models were relaxed before the simulation. The trajectories were saved after every 100 ps for analysis, and the stability of simulations was evaluated by calculating the root mean square deviation (RMSD) of the protein and the ligand over time. The Desmond simulation trajectories were analyzed. RMSD, root mean square fluctuation (RMSF), and protein-ligand contacts were calculated from MD trajectory analysis.

3. Results and Discussion

Protein-protein interactions play vital roles in cellular activities as these hold the major tasks of the biological machinery. These complex (e.g., protein-peptide complex or protein-ligand complex) interactions help to sort out the mysterious signaling and pathways related to the functioning of living systems. To decode the possible interactions of AMPs against bacterial receptor proteins, computer-mediated molecular docking was performed. HADDOCK server v.2.4 was used to predict the possible protein-protein interactions between AMPs and pathogenic bacterial receptor proteins.

Medicinal plants have a long history in both the traditional and modern medicines in different communities across the world. The crude extract and herbal decoction of many plants have been reported with biologically active compounds that play a significant role in the treatment of many diseases [20]. Recent studies on natural flora have shown the presence of biologically active peptides and phytochemicals with reported therapeutic activities [21]. Proteins are the workhorse of cells with diverse cellular functions as these serve as messengers, modifiers, scaffolds, catalysts, and signal receptors. Protein interactions with other proteins, DNA, RNA, and peptides are responsible for various biological activities [22].

Among ten selected AMPs, only two peptides (i.e., napin and snakin-1) showed the best binding interactions and HADDOCK scores against selected pathogenic bacterial species. The peptide-protein complexes with the lowest binding energy were considered to be the most stable ones (Table 1). Napin showed strong interactions and hydrogen bonding with proteins of three bacterial strains (i.e., Acinetobacter baumannii, Mycobacterium tuberculosis, and Salmonella enterica) while snakin-1 binds with the interacting residues of two bacterial proteins (i.e., Helicobacter pylori and methicillin-resistant Staphylococcus aureus).

3.1. Interactions between AMPs and PBP1a of Acinetobacter baumannii

Acinetobacter baumannii is a Gram-negative bacterium and has been listed in the group of ESKAPE pathogens which are responsible for a variety of infections, most commonly respiratory and urinary tract infections. The lethality and prevalence of Acinetobacter have increased due to resistance against different antibiotics [23]. Over the last three decades, this bacterium has acquired resistance against antibiotics due to adaptational changes in enzymes and cellular proteins [24]. The mortality rate of Acinetobacter outbreaks was 50-65% as most of the patients died within 48 hours of hospitalization [25].

Penicillin-binding proteins are diverse bifunctional enzymes and classified as PBP1a and PBP1b that perform the assembly of the bacterial cell wall. Different β-lactam antibiotics disrupt the bacterial cell wall synthesis by covalently inactivating the penicillin-binding proteins. Access to these periplasmic targets can be helpful to inhibit bacterial activity and appeared as an emerging step in meeting the new challenges represented by multidrug-resistant bacteria [26]. In this study, napin with a HADDOCK score of -158.7 kcal/mol showed strong interactions with active amino acids of penicillin-binding protein 1a of A. baumannii. The protein is highlighted as grey color with highlighted red sticks as interactive amino acid residues in Figure 1. In a study, thirty-four compounds were tested using a molecular docking approach against PBP1a of Acinetobacter baumannii, and neogrifoline and 3,11-dioxolanosta-8,24(Z)-diene-26-oic acid exhibited the best results [13] and could be potential drug candidates against this bacterial species. In another study, Skariyachan et al. [27] used herbal-based ligands to predict receptor-ligand interactions by molecular docking. They revealed that the herbal ligand imipenem exhibited a binding energy of −5.9 kcal/mol when docked with PBP1a and suggested that the lead compound and the target could be used for structure-based drug designing against A. baumannii.

(a)

(b)

Figure 1

Protein-protein interactions between napin and penicillin-binding protein 1a of A. baumannii. (a) PBP1a is represented in grey color with red interacting residues, and napin is shown in blue color with green interacting residues. (b) All interacting residues between napin and PBP1a of A. baumannii complex; hydrogen bonds are shown in blue color, and salt bridges are represented by red-colored lines. Other representing properties of amino acids are represented by different colors (positive: blue, neutral: green, negative: red, aliphatic: grey, Pro&Gly: orange, aromatic: pink, and Cys: yellow).

3.2. Interactions between AMPs and RmlA of Mycobacterium tuberculosis

Mycobacterium tuberculosis is a facultative intracellular pathogen and causative agent of tuberculosis. It has remained one of the main causes of increased mortality and morbidity with approximately two million deaths worldwide [28]. Tuberculosis (TB) is known as a major threat to humanity for the past three decades due to the emergence of multidrug-resistant strains of this bacterium [29]. The cell wall of M. tuberculosis consists of three interconnected molecules including arabinogalactan, mycolic acids, and peptidoglycan which are responsible for bacterial cell viability [30]. The arabinogalactan is linked to the sixth position of the muramic acid residue of peptidoglycan via disaccharide linker α-l-rhamnosyl-(1 → 3)-α-d-N-acetylglucosaminosyl-1-phosphate [31].

The rhamnose is the precursor in the rhamnosyl biosynthesis pathway not only in M. tuberculosis but also in a wide range of bacterial species [28]. Therefore, rhamnose plays a crucial role in the attachment of arabinogalactan to peptidoglycan in the bacterial cell wall. Glucose-1-phosphate thymidylyltransferase (RmlA) serves as the leading enzyme of the rhamnose biosynthesis pathway and is therefore essential for the survival of M. tuberculosis [32]. The glucose-1-phosphate thymidylyltransferase (RmlA) also plays an essential role in bacterial cell wall viability, and therefore, RmlA could serve as a major target in the prevention of this infection. In this study, we focused on protein-protein docking between selected AMPs and glucose-1-phosphate thymidylyltransferase of M. tuberculosis to explore the potential of AMPs against RmlA. Among ten AMPs, napin with a score of -107.8 kcal/mol showed twelve hydrogen bonds with glucose-1-phosphate thymidylyltransferase of M. tuberculosis. The interactions have been shown in Figure 2. Mansuri et al. [32] docked two compounds (i.e., 6-[(2R,3S,5R)-5-[5-(2-aminoethyl)-2,4-dioxo-1,2,3,4-tetrahydropyrimidin-1-yl]-3-hydroxyoxolan-2-yl] hexanoic acid and 4-(2-{1-[(1S,3S,4S)-3-(5-carboxypentyl)-4-hydroxy-2-methylidenecyclopentyl]-2,4-dioxo-1,2,3,4-tetrahydropyrimidin-5-yl}ethyl)morpholin-4-ium) against glucose-1-phosphate thymidylyltransferase of M. tuberculosis and reported that these compounds could be used as competitive inhibitors.

(a)

(b)

Figure 2

Protein-protein interactions between napin and glucose-1-phosphate thymidylyltransferase of M. tuberculosis H37Rv. (a) Glucose-1-phosphate thymidylyltransferase is represented in yellow-orange color with deep salmon-colored interacting residues, and napin is shown in blue color with green interacting residues. (b) All interacting residues between napin and glucose-1-phosphate thymidylyltransferase protein of A. baumannii complex; hydrogen bonds are shown in blue color, and salt bridges are represented by red-colored lines. Other representing properties of amino acids are represented by different colors (positive: blue, neutral: green, negative: red, aliphatic: grey, Pro&Gly: orange, aromatic: pink, and Cys: yellow).

3.3. Interactions between AMPs and Streptomycin 3-Adenylyltransferase of Salmonella enterica

Salmonella enterica is a facultative Gram-negative intracellular bacterium that infects both animals and humans. Growing pieces of evidence have pointed towards more severity of infections due to the resistant strains of S. enterica [33]. The streptomycin and spectinomycin belong to aminocyclitol and aminoglycoside families of antibiotics, respectively, and bind to the bacterial ribosome and interfere with the protein biosynthesis.

Currently, the most reported resistance in S. enterica is due to the inactivation of these drugs by aminoglycoside-modifying enzymes such as aminoglycoside nucleotidyltransferases (ANTs). Aminoglycoside (3) (9) adenylyltransferase AadA from S. enterica is the member of the ANT(3)-Ia family that O-adenylates the streptomycin and spectinomycin at specific positions. Thus, AadA is the promising target to hinder the protein machinery of this bacterium [34]. In this study, napin showed the best binding patterns with streptomycin 3-adenylyltransferase protein of S. enterica with a score of -84.2 kcal/mol and only four hydrogen bonds (Figure 3). The HADDOCK server predicts and gives different numbers of hydrogen bonds for different protein-peptide complexes as it depends on how strongly the peptide is interacting with the receptor protein. In a study, Prabhu et al. [35] used molecular docking and molecular dynamics simulation approaches to predict structural, binding, and pharmacokinetic properties of different selected compounds. They docked these compounds against streptomycin 3-adenylyltransferase of Serratia marcescens and reported that the best five identified compounds could be used as potential drug entities to develop antipathogenic agents.

(a)

(b)

Figure 3

Protein-protein interactions between napin and streptomycin 3-adenylyltransferase protein of S. enterica. (a) Streptomycin 3-adenylyltransferase is represented in light blue color with orange-colored interacting residues, and napin is shown in blue color with green interacting residues. (b) All interacting residues between napin and streptomycin 3-adenylyltransferase protein of A. baumannii complex; hydrogen bonds are shown in blue color, and salt bridges are represented by red-colored lines. Other representing properties of amino acids are represented by different colors (positive: blue, neutral: green, negative: red, aliphatic: grey, Pro&Gly: orange, aromatic: pink, and Cys: yellow).

3.4. Interactions between AMPs and Oxygen-Insensitive NADPH Nitroreductase of Helicobacter pylori

Helicobacter pylori is a microaerophilic Gram-negative bacterium that colonizes the gastric mucosa. H. pylori is one of the world’s most common pathogens affecting about 50% of the world’s population. This bacterium affects the human gastrointestinal tract and causes upper gastrointestinal tract infections such as chronic gastritis, ulcerative colitis, gastrointestinal/mucosa-associated lymphoid tissue (MALT) lymphoma, duodenal ulcer, and gastric carcinoma [36]. The frequent consumption and high dosage of the metronidazole (MTZ) antibiotic are responsible for resistance in the bacterium with a negative impact on the treatment efficacy. The mutations in the rdxA gene (oxygen-insensitive) and frxA gene (flavin reductase) that encode NADPH nitroreductases have been associated with the onset of metronidazole resistance by H. pylori [37]. Snakin-1 with a score of -105.0 kcal/mol showed 13 hydrogen bonds with oxygen-insensitive NADPH nitroreductase of H. pylori. The interactions have been displayed in Figure 4. In a study, Mulimani et al. [38] used 100 compounds in a molecular docking study against oxygen-insensitive NADPH nitroreductase of H. pylori and reported that benzimidazoles and oxacillin exhibited the best results and therefore could be used as potential inhibitors.

(a)

(b)

Figure 4

Protein-protein interactions between snakin-1 and oxygen-insensitive NADPH nitroreductase of H. pylori. (a) NADPH nitroreductase is represented in yellow color with orange-colored interacting residues, and snakin-1 is shown in blue color with green interacting residues. (b) All interacting residues between snakin-1 and NADPH nitroreductase protein of A. baumannii complex; hydrogen bonds are shown in blue color, and salt bridges are represented by red-colored lines. Other representing properties of amino acids are represented by different colors (positive: blue, neutral: green, negative: red, aliphatic: grey, Pro&Gly: orange, aromatic: pink, and Cys: yellow).

3.5. Interactions between AMPs and PBP2a of MRSA

Staphylococcus aureus is a Gram-positive pathogen that is capable of spreading a wide spectrum of infections. Different strains of S. aureus have been isolated in the past which are resistant to multiple drugs and responsible for the severe outbreaks of infections worldwide. Among all the isolated strains, MRSA is notably known as a resistant strain due to the unique genetic element staphylococcal chromosomal cassette mec (SCCmec) that carries the mecA gene which encodes penicillin-binding protein 2a (PBP2a) [39]. PBP2a is the product of the mutant gene which plays a role as surrogate transpeptidase in the absence of other PBPs [40].

S. aureus (MRSA) manifests the resistance to methicillin and other β-lactams due to the production of PBP2a that exhibits transpeptidase activity for cell wall biosynthesis [39]. Hence, there is a dire need for new drugs to inhibit the biosynthesis of the cell wall of S. aureus. In this study, snakin-1 showed the best binding pattern with a score of -103.8 kcal/mol and exhibited twenty-three hydrogen bonds with PBP2a of S. aureus (MRSA). The detail has been discussed in cartoon representation with highlighted portions in Figure 5. In a study, Murugavel et al. [41] revealed the inhibitory activity of methyl (2E)-2-{[N-(2-formylphenyl) (4-methylbenzene)sulfonamido]methyl}-3-(4-chlorophenyl) prop-2-enoate (MFMSC) when docked against PBP-2X. Similarly, Levy et al. [42] docked the compound trans 2-(aminomethyl)-4-oxazol-5-yl-7-oxo-1,6-diazabicyclo[3.2.1]oct-3-en-6-yl] hydrogen sulfate (CPD4) against PBP2 of Escherichia coli and reported that CPD4 could be a potential scaffold for the development of active molecules which would be effective against a broad range of bacterial species.

(a)

(b)

Figure 5

Protein-protein interactions between snakin-1 and penicillin-binding protein 2a from methicillin-resistant Staphylococcus aureus. (a) PBP2a is represented in pink color with grey-colored interacting residues, and snakin-1 is shown in blue color with green interacting residues. (b) All interacting residues between snakin-1 and PBP2a protein of A. baumannii complex; hydrogen bonds are shown in blue color, and salt bridges are represented by red-colored lines. Other representing properties of amino acids are represented by different colors (positive: blue, neutral: green, negative: red, aliphatic: grey, Pro&Gly: orange, aromatic: pink, and Cys: yellow).

3.6. Molecular Dynamics Simulation

Based on the best HADDOCK scores, two complexes (i.e., napin with PBP1a of A. baumannii and snakin-1 with oxygen-insensitive NADPH nitroreductase of H. pylori) were selected for MD simulation studies. The evolution of RMSD values in the course of time for the C-alpha atoms of the peptide-bound proteins is shown in Figures 6 and 7. The RMSD plot of the complex Napin-PBP indicated that the complex reached stability at 20 ns, and the RMSD plot of the complex Snakin-NADPH indicated that the complex reached stability at 5 ns. From then, an average RMSD value of 1.5 Å for Napin-PBP and 1.0 Å for Snakin-NADPH persists up to 100 ns during the simulation period, which is quite acceptable. Peptides fit to proteins, and RMSD values fluctuated within 2.5 Å after being stable. These indicate that the peptides remained stably bound to the binding sites of their respective receptors during the simulation period. However, there is more deviation in Napin-PBP compared to Snakin-NADPH.

The residue-wise RMSF values of proteins bound to their respective peptides are shown in Figures 8 and 9. The residues are showing higher peaks corresponding to loop regions, as identified from MD trajectories (Figure S1 and S2), or N- and C-terminal zones. Low RMSF values of binding site residues indicate the stability of the peptide binding with the protein.

Most of the important interactions of protein-peptide determined with MD are hydrogen bonds, as depicted in Figures 10 and 11.

In silico studies on pharmacological and nutraceutical applications of plant-derived compounds are more efficient and less laborious [43]. Molecular docking is an elaborative method which helps scientists to forecast the best binding patterns and interactions between desired compounds before the experimental and laboratory approaches [44]. Computational biology applies different coherent and integrated approaches to analyze and explore new features from large collections of biological data such as whole genomes and proteomes of eukaryotes and prokaryotes [45]. Different docking techniques such as ligand-based molecular docking, protein-protein docking, protein-peptide docking, and induced fit docking have greatly influenced the field of drug discovery. Computer-aided drug designing using molecular docking and molecular dynamics simulation approaches provides rapid screening of novel and potential drug candidates to predict drug-receptor interactions [46]. In recent years, plant-derived natural compounds have been proved as potent drug candidates and inhibitors of many pathogenic proteins that play crucial roles in the pathogenesis of various diseases [47]. Recently, using computational modeling, the molecular dynamics aspects of moxifloxacin-induced resistance in M. tuberculosis DNA gyrase A and C have been studied by Pandey et al. [48]. Similarly, Bera et al. [49] also used molecular docking and simulation approaches to study interactions of Echinocandin B with the multidrug resistance-associated protein family of ATP-binding transporter protein.

The alteration of bacterial proteins and membrane permeability are the major reasons for the rapid emergence of bacterial resistance towards antibiotics [50]. The level of antibiotic resistance in bacteria is rising very rapidly due to the emergence of new bacterial strains. Along with modifications in bacterial proteins and enzymes, the second leading cause of the emergence of antibiotic resistance is the behavioral changes and prescriptions of antibiotics [51]. With the continuous use of these drugs against bacterial infections without any valid prescription, the emergence and spread of bacterial resistance would worsen in the coming time [52]. Bacterial resistance towards antibiotics is the main cause of outbreaks of specific infections. The mutations and advanced modifications in the bacterial machinery are the leading causes of antimicrobial resistance. The modified enzymes and mutant genes make bacteria no longer responsive to specific antibiotics. According to WHO, these multidrug-resistant bacteria are known as superbugs because current antibiotics are inefficient against them [53]. On the basis of the current situation, there is an urgent need for new and novel antimicrobial drugs that could inhibit the bacterial machinery directly.

The antibiotics which are used to treat A. baumannii, M. tuberculosis, H. pylori, S. aureus, and S. enterica are becoming more limited. Little is known about the contributions of penicillin-binding protein 1a, glucose-1-phosphate thymidylyltransferase, oxygen-insensitive NADPH nitroreductase, penicillin-binding protein 2a, and streptomycin 3-adenylyltransferase which could serve as excellent targets for the development of new antibiotics. In the literature, using molecular docking and simulation studies, the compounds from plant sources have shown better binding with strong interactions, hydrogen bonding, binding energy, and other weak interactions with their normal targets compared to different antibacterial agents which have been used conventionally (e.g., imipenem, polymyxin E, and clinafloxacin) [27]. The aim of this study was therefore to explore the binding interactions between active amino acid residues of AMPs and five different bacterial receptor proteins. This study will help for further validation and exploration of new drugs to cure these infections. The results of this study have proved that these infections are no longer incurable but the treatments might be hidden in some other aspects.

4. Conclusion

The current study focuses on the computational prediction of antimicrobial peptides of plant sources against prioritized targets of different bacterial species. Molecular docking and molecular dynamics simulations suggested that among ten peptides, only napin and snakin-1 showed strong binding interactions and the best scores with the selected receptor proteins of bacterial strains. Napin showed noteworthy interactions with three bacterial strains including A. baumannii, M. tuberculosis, and S. enterica while snakin-1 revealed binding interactions with interacting residues of two bacterial proteins such as H. pylori and methicillin-resistant Staphylococcus aureus. Further, the molecular dynamics simulation studies also confirmed that the peptides napin and snakin-1 remained firmly bound to the binding sites of proteins PBP1a of A. baumannii and oxygen-insensitive NADPH nitroreductase of H. pylori, respectively, during the simulation period. The findings of this study will help for further advancement and development of new drugs from the natural flora. The aim of this study was to explore the interactive sites and interactions between AMPs and active amino acids of selected bacterial receptor proteins in order to inhibit and target them directly.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Authors’ Contributions

GM and HSM conceived and planned the experiments; HSM and RM carried out the molecular docking experiments; KM and SA carried out molecular docking simulation studies and drafted the manuscript; GM supervised the project and proofread the article. All authors discussed the results and commented on the manuscript.

Acknowledgments

The authors would like to gratefully acknowledge the Department of Biochemistry, Government College University Faisalabad, for providing space and facilities to accomplish this study.

Supplementary Materials

Figure S1: protein secondary structure element distribution by the residue index throughout the protein structure (Napin-PBP1a). Red columns indicate α-helices, and blue columns indicate β-strands. Figure S2: protein secondary structure element distribution by the residue index throughout the protein structure (Snakin-NADPH). Red columns indicate α-helices, and blue columns indicate β-strands. (Supplementary Materials). (Supplementary Materials)

References

N. Kumar, D. Sood, R. Tomar, and R. Chandra, “Antimicrobial peptide designing and optimization employing large-scale flexibility analysis of protein-peptide fragments,” ACS Omega, vol. 4, no. 25, pp. 21370–21380, 2019.
View at: Publisher Site | Google Scholar
M. Mahlapuu, J. Håkansson, L. Ringstad, and C. Björn, “Antimicrobial peptides: an emerging category of therapeutic agents,” Frontiers in Cellular and Infection Microbiology, vol. 6, p. 194, 2016.
View at: Publisher Site | Google Scholar
C. E. Salas, J. A. Badillo-Corona, G. Ramírez-Sotelo, and C. Oliver-Salvador, “Biologically active and antimicrobial peptides from plants,” BioMed research international, vol., vol. 2015, pp. 1–11, 2015.
View at: Publisher Site | Google Scholar
M. Hassan, M. Kjos, I. Nes, D. Diep, and F. Lotfipour, “Natural antimicrobial peptides from bacteria: characteristics and potential applications to fight against antibiotic resistance,” Journal of Applied Microbiology, vol. 113, no. 4, pp. 723–736, 2012.
View at: Publisher Site | Google Scholar
A. Iqbal, R. S. Khan, K. Shehryar et al., “Antimicrobial peptides as effective tools for enhanced disease resistance in plants,” Plant Cell, Tissue and Organ Culture (PCTOC), vol. 139, no. 1, pp. 1–15, 2019.
View at: Publisher Site | Google Scholar
J. P. Tam, S. Wang, K. H. Wong, and W. L. Tan, “Antimicrobial peptides from plants,” Pharmaceuticals, vol. 8, no. 4, pp. 711–757, 2015.
View at: Publisher Site | Google Scholar
H.-Y. Yi, M. Chowdhury, Y.-D. Huang, and X.-Q. Yu, “Insect antimicrobial peptides and their applications,” Applied Microbiology and Biotechnology, vol. 98, no. 13, pp. 5807–5822, 2014.
View at: Publisher Site | Google Scholar
Y. Qin, Z. D. Qin, J. Chen et al., “From antimicrobial to anticancer peptides: the transformation of peptides,” Recent Patents on Anti-Cancer Drug Discovery, vol. 14, no. 1, pp. 70–84, 2019.
View at: Publisher Site | Google Scholar
R. Nawrot, J. Barylski, G. Nowicki, J. Broniarczyk, W. Buchwald, and A. Goździcka-Józefiak, “Plant antimicrobial peptides,” Folia Microbiologica, vol. 59, no. 3, pp. 181–196, 2014.
View at: Publisher Site | Google Scholar
O. Silva, C. De La Fuente-Núñez, E. Haney et al., “An anti-infective synthetic peptide with dual antimicrobial and immunomodulatory activities,” Scientific Reports, vol. 6, no. 1, pp. 1–11, 2016.
View at: Publisher Site | Google Scholar
C. D. Fjell, J. A. Hiss, R. E. Hancock, and G. Schneider, “Designing antimicrobial peptides: form follows function,” Nature Reviews Drug Discovery, vol. 11, no. 1, pp. 37–51, 2012.
View at: Publisher Site | Google Scholar
V. Cunsolo, R. Schicchi, M. Chiaramonte et al., “Identification of new antimicrobial peptides from Mediterranean medical plant Charybdis pancration (Steinh.) Speta,” Antibiotics, vol. 9, no. 11, p. 747, 2020.
View at: Publisher Site | Google Scholar
M. J. Alves, H. J. Froufe, A. F. Costa et al., “Docking studies in target proteins involved in antibacterial action mechanisms: extending the knowledge on standard antibiotics to antimicrobial mushroom compounds,” Molecules, vol. 19, no. 2, pp. 1672–1684, 2014.
View at: Publisher Site | Google Scholar
R. Arif, S. Ahmad, G. Mustafa et al., “Molecular docking and simulation studies of antidiabetic agents devised from hypoglycemic polypeptide-P of Momordica charantia,” BioMed research international, vol. 2021, Article ID 5561129, 2021.
View at: Publisher Site | Google Scholar
Y. Matsuzaki, M. Ohue, N. Uchikoga, and Y. Akiyama, “Protein-protein interaction network prediction by using rigid-body docking tools: application to bacterial chemotaxis,” Protein and Peptide Letters, vol. 21, no. 8, pp. 790–798, 2014.
View at: Publisher Site | Google Scholar
S. R. Shrivastava, P. S. Shrivastava, and J. Ramasamy, “World health organization releases global priority list of antibiotic-resistant bacteria to guide research, discovery, and development of new antibiotics,” Journal of Medical Society, vol. 32, no. 1, pp. 76-77, 2012.
View at: Publisher Site | Google Scholar
G. Van Zundert, J. Rodrigues, M. Trellet et al., “The HADDOCK2.2 web server: user-friendly integrative modeling of biomolecular complexes,” Journal of Molecular Biology, vol. 428, no. 4, pp. 720–725, 2016.
View at: Publisher Site | Google Scholar
W. L. DeLano, “PyMOL: an open-source molecular graphics tool,” CCP4 Newsletter on protein crystallography, vol. 40, no. 1, pp. 82–92, 2002.
View at: Google Scholar
R. A. Laskowski, “PDBsum new things,” Nucleic acids research, vol. 37, Supplement 1, no. Database, pp. D355–D359, 2009.
View at: Publisher Site | Google Scholar
G. Mustafa, R. Arif, A. Atta, S. Sharif, and A. Jamil, “Bioactive compounds from medicinal plants and their importance in drug discovery in Pakistan,” vol. 1, no. 1, pp. 17–26, 2017.
View at: Publisher Site | Google Scholar
G. Mustafa, S. Ahmed, N. Ahmed, and A. Jamil, “Phytochemical and antibacterial activity of some unexplored medicinal plants of Cholistan desert,” Pakistan Journal of Botany, vol. 48, no. 5, pp. 2057–2062, 2016.
View at: Google Scholar
A. Alekseenko, M. Ignatov, G. Jones, M. Sabitova, and D. Kozakov, “Protein–protein and protein–peptide docking with ClusPro server,” Protein Structure Prediction, Springer, pp. 157–174, 2020.
View at: Google Scholar
V. Tiwari, M. Tiwari, and V. Solanki, “Polyvinylpyrrolidone-capped silver nanoparticle inhibits infection of carbapenem-resistant strain of Acinetobacter baumannii in the human pulmonary epithelial cell,” Frontiers in Immunology, vol. 8, p. 973, 2017.
View at: Publisher Site | Google Scholar
S. Zeidler and V. Müller, “Coping with low water activities and osmotic stress in Acinetobacter baumannii: significance, current status and perspectives,” Environmental Microbiology, vol. 21, no. 7, pp. 2212–2230, 2019.
View at: Publisher Site | Google Scholar
S. Skariyachan, A. G. Muddebihalkar, V. Badrinath et al., “Natural epiestriol-16 act as potential lead molecule against prospective molecular targets of multidrug resistant Acinetobacter baumannii -insight from _in silico_ modelling and _in vitro_ investigations,” Genetics and Evolution, vol. 82, p. 104314, 2020.
View at: Publisher Site | Google Scholar
S. Han, N. Caspers, R. P. Zaniewski et al., “Distinctive attributes of β-lactam target proteins in Acinetobacter baumannii relevant to development of new antibiotics,” Journal of the American Chemical Society, vol. 133, no. 50, pp. 20536–20545, 2011.
View at: Publisher Site | Google Scholar
S. Skariyachan, M. Manjunath, and N. Bachappanavar, “Screening of potential lead molecules against prioritised targets of multi-drug-resistant-Acinetobacter baumannii–insights from molecular docking, molecular dynamic simulations and in vitro assays,” Journal of Biomolecular Structure and Dynamics, vol. 37, no. 5, pp. 1146–1169, 2019.
View at: Publisher Site | Google Scholar
H. A. Brown, J. B. Thoden, P. A. Tipton, and H. M. Holden, “The structure of glucose-1-phosphate thymidylyltransferase from Mycobacterium tuberculosis reveals the location of an essential magnesium ion in the RmlA-type enzymes,” Protein Science, vol. 27, no. 2, pp. 441–450, 2018.
View at: Publisher Site | Google Scholar
G. Delogu, M. Sali, and G. Fadda, “The biology of Mycobacterium tuberculosis infection,” Mediterranean journal of hematology and infectious diseases, vol. 5, no. 1, p. e2013070, 2013.
View at: Publisher Site | Google Scholar
L. Alderwick, H. Birch, A. K. Mishra, L. Eggeling, and G. Besra, “Structure, function and biosynthesis of theMycobacterium tuberculosiscell wall: arabinogalactan and lipoarabinomannan assembly with a view to discovering new drug targets,” Biochemical Society Transactions, vol. 35, no. 5, pp. 1325–1328, 2007.
View at: Publisher Site | Google Scholar
S. Li, J. Kang, W. Yu et al., “Identification of M. tuberculosis Rv3441c and M. smegmatis MSMEG_1556 and essentiality of M. smegmatis MSMEG_1556,” PLoS One, vol. 7, no. 8, article e42769, 2012.
View at: Publisher Site | Google Scholar
R. Mansuri, M. Y. Ansari, J. Singh et al., “Computational elucidation of structural basis for ligand binding with Mycobacterium tuberculosis glucose-1-phosphate thymidylyltransferase (RmlA),” Current Pharmaceutical Biotechnology, vol. 17, no. 12, pp. 1089–1099, 2016.
View at: Publisher Site | Google Scholar
M. Silbergleit, A. A. Vasquez, C. J. Miller, J. Sun, and I. Kato, “Oral and intestinal bacterial exotoxins: potential linked to carcinogenesis,” Progress in Molecular Biology and Translational Science, vol. 171, pp. 131–193, 2020.
View at: Publisher Site | Google Scholar
A. L. Stern, S. E. Van der Verren, J. Näsvall, H. Gutiérrez-de-Terán, and M. Selmer, “Structural mechanism of AadA, a dual-specificity aminoglycoside adenylyltransferase from _Salmonella enterica_,” Journal of Biological Chemistry, vol. 293, no. 29, pp. 11481–11490, 2018.
View at: Publisher Site | Google Scholar
D. Prabhu, R. Vidhyavathi, and J. Jeyakanthan, “Computational identification of potent inhibitors for streptomycin -adenylyltransferase of _Serratia marcescens_,” Microbial Pathogenesis, vol. 103, pp. 94–106, 2017.
View at: Publisher Site | Google Scholar
I. Thung, H. Aramin, V. Vavinskaya et al., “Review article: the global emergence ofHelicobacter pyloriantibiotic resistance,” Alimentary Pharmacology & Therapeutics, vol. 43, no. 4, pp. 514–533, 2016.
View at: Publisher Site | Google Scholar
B. Marques, M. M. Donato, O. Cardoso, C. Luxo, A. Martinho, and N. Almeida, “Study of _rdxA_ and _frxA_ genes mutations in metronidazole-resistant and -susceptible _Helicobacter pylori_ clinical isolates from the central region of Portugal,” Journal of global antimicrobial resistance, vol. 17, pp. 300–304, 2019.
View at: Publisher Site | Google Scholar
N. D. Mulimani, B. B. Udapudi, and D. A. Yaraguppi, “In-silico drug design approach for identification of antibacterial drug leads for treatment of gastric lymphoma,” International Journal, vol. 2, no. 1, 2013.
View at: Google Scholar
M.-A. W. Shalaby, E. M. Dokla, R. A. Serya, and K. A. Abouzid, “Penicillin binding protein 2a: an overview and a medicinal chemistry perspective,” European Journal of Medicinal Chemistry, vol. 199, article 112312, 2020.
View at: Publisher Site | Google Scholar
S. B. Mohamed, T. A. Adlan, N. A. Khalafalla et al., “Proteomics and docking study targeting penicillin-binding protein and penicillin-binding protein2a of methicillin-resistant Staphylococcus aureus strain SO-1977 isolated from Sudan,” Evolutionary Bioinformatics, vol. 15, 2019.
View at: Google Scholar
S. Murugavel, D. Kannan, and M. Bakthadoss, “Experimental and computational approaches of a novel methyl (2 _E_ )-2-{[N-(2-formylphenyl)(4-methylbenzene)sulfonamido]methyl}-3-(4-chlorophenyl)prop-2-enoate: a potential antimicrobial agent and an inhibition of penicillin-binding protein,” Journal of Molecular Structure, vol. 1115, pp. 33–54, 2016.
View at: Publisher Site | Google Scholar
N. Levy, J.-M. Bruneau, E. Le Rouzic et al., “Structural basis for E. coli penicillin binding protein (PBP) 2 inhibition, a platform for drug design,” Journal of Medicinal Chemistry, vol. 62, no. 9, pp. 4742–4754, 2019.
View at: Publisher Site | Google Scholar
G. Mustafa, M. Majid, A. Ghaffar, M. Yameen, H. A. Samad, and H. S. Mahrosh, “Screening and molecular docking of selected phytochemicals against NS5B polymerase of hepatitis C virus,” Pakistan Journal of Pharmaceutical Sciences, vol. 33, no. 5, pp. 2317–2322, 2020.
View at: Google Scholar
H. S. Mahrosh and G. Mustafa, “An in silico approach to target RNA-dependent RNA polymerase of COVID-19 with naturally occurring phytochemicals,” Environment, Development and Sustainability, pp. 1–14, 2021.
View at: Publisher Site | Google Scholar
A. Mushtaq, G. Mustafa, T. M. Ansari, M. A. Shad, J. Cruz-Reyes, and A. Jamil, “Antiviral activity of hexapeptides derived from conserved regions of bacterial proteases against HCV NS3 protease,” Pakistan Journal of Pharmaceutical Sciences, vol. 34, no. 1, pp. 215–223, 2021.
View at: Google Scholar
G. Mustafa, H. S. Mahrosh, M. Salman et al., “Identification of peptides as novel inhibitors to target IFN-γ, IL-3, and TNF-α in systemic lupus erythematosus,” BioMed Research International, vol. 2021, Article ID 1124055, 11 pages, 2021.
View at: Publisher Site | Google Scholar
G. Mustafa, H. S. Mahrosh, and R. Arif, “In silico characterization of growth differentiation factors as inhibitors of TNF-alpha and IL-6 in immune-mediated inflammatory disease rheumatoid arthritis,” BioMed research international, vol. 2021, Article ID 5538535, 11 pages, 2021.
View at: Publisher Site | Google Scholar
B. Pandey, S. Grover, C. Tyagi et al., “Dynamics of fluoroquinolones induced resistance in DNA gyrase of Mycobacterium tuberculosis,” Journal of Biomolecular Structure and Dynamics, vol. 36, no. 2, pp. 362–375, 2018.
View at: Publisher Site | Google Scholar
K. Bera, P. Rani, G. Kishor, S. Agarwal, A. Kumar, and D. V. Singh, “Structural elucidation of transmembrane domain zero (TMD0) of EcdL: a multidrug resistance-associated protein (MRP) family of ATP-binding cassette transporter protein revealed by atomistic simulation,” Journal of Biomolecular Structure and Dynamics, vol. 36, no. 11, pp. 2938–2950, 2018.
View at: Publisher Site | Google Scholar
C. L. Ventola, “The antibiotic resistance crisis: part 1: causes and threats,” Pharmacy and therapeutics, vol. 40, no. 4, pp. 277–283, 2015.
View at: Google Scholar
M. U. Qamar, G. Mustafa, U. Qaisar et al., “Molecular detection of blaNDM and blaVIM in clinically isolated multi-drug resistant Escherichia coli in Pakistan,” Pakistan Journal of Pharmaceutical Sciences, vol. 32, no. S5, pp. 2305–2309, 2019.
View at: Google Scholar
M. Frieri, K. Kumar, and A. Boutin, “Antibiotic resistance,” Journal of Infection and Public Health, vol. 10, no. 4, pp. 369–378, 2017.
View at: Publisher Site | Google Scholar
A. Chakraborty, K. Poira, D. Saha, C. Halder, and S. Das, “Multidrug-resistant bacteria with activated and diversified MDR genes in Kolkata water: Ganga Action Plan and heterogeneous phyto-antibiotics tackling superbug spread in India,” American Journal of Drug Delivery and Therapeutics, vol. 5, no. 1, pp. 1–9, 2018.
View at: Google Scholar

Copyright

Copyright © 2022 Ghulam Mustafa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1321

Downloads

953

Citations