https://doi.org/10.5812/ijhls.121109

Discovery of a Potent Inhibitor to Overcome Carbapenem Resistance in Pseudomonas aeruginosa Strains via Inhibition of VIM-2 Metallo-β-lactamases

authors:

avatar Seyedeh Sara Moosavi ORCID 1 , avatar Mahboobeh Madani ORCID 1 , * , avatar Sako Mirzaie ORCID 2 , ** , avatar Nima Hosseini Jazani ORCID 3

Department of Microbiology, Falavarjan Branch, Islamic Azad University, Isfahan, Iran
Department of Biochemistry, Sanandaj Branch, Islamic Azad University, Sanandaj, Iran
Department of Molecular Microbiology, Nextgen Laboratories, Newport Beach, California, USA
Corresponding Authors:
Journal of Health Reports and Technology: Vol.8, issue 2; e121109
published online: April 2, 2022
article type: Research Article
received: November 13, 2021
revised: November 23, 2021
accepted: February 8, 2022

how to cite: Moosavi S S, Madani M, Mirzaie S, Hosseini Jazani N. Discovery of a Potent Inhibitor to Overcome Carbapenem Resistance in Pseudomonas aeruginosa Strains via Inhibition of VIM-2 Metallo-β-lactamases. J Health Rep Technol. 2022;8(2):e121109. https://doi.org/10.5812/ijhls.121109.

Abstract

Background:

Mutations in bacteria frequently occur that display a crucial need for new antimicrobial agents. Metallo-β-lactamases (MBLs) are a growing threat to maintaining the effectiveness of beta-lactam antibiotics. Resistance to beta-lactam antibiotics is one of the most common types of antibiotic resistance, which causes the ineffectiveness of antibiotics.

Objectives:

This study aimed to identify a novel inhibitor using molecular dynamics simulations to inhibit VIM-2 Metallo-β-lactamases and overcome carbapenem resistance in Pseudomonas aeruginosa strains.

Methods:

Computational biology tools were employed for this study, including molecular dynamics, binding free energy, virtual screening, and docking. Natural compounds were taken from the ZINC databank and prepared. At the next stage, the prepared compounds were screened based on docking energy in the active site of VIM-2 MBL by Schrödinger (Maestro) software, and better compounds were selected. Captopril was chosen as a positive control inhibitor for VIM-2 MBLs. Ultimately, molecular dynamics simulations were performed using GROMACS software, and outputs were analyzed.

Results:

Maestro software's screening results showed that ZINC00517765 was the best inhibitor with -12.29 kcal mol-1 docking energy. The ADME investigations revealed that ZINC00517765 had an appropriate range of pharmacokinetics, lipophilicity, and drug-likeness features as an inhibitor of VIM-2 MBL. Molecular dynamics outcomes explicated that VIM-2 MBL in the presence of ZINC00517765 had better stability during simulation. The results of the MM-PBSA study illustrated that ZINC00517765 with -72.29 kJ mol-1 binding free energy was more potent than Captopril with -23.39 kJ mol-1.

Conclusions:

This study showed that VIM-2 MBL in the presence of ZINC00517765 has suitable stability during simulation. Also, more hydrogen bonds and stronger binding free energy than Captopril confirm that ZINC00517765 is a proper candidate for further studies and laboratory investigation.

1. Background

The decreased effect of β-lactam antibiotics on Gram-negative pathogens is a growing clinical problem, leading to antibacterial drug resistance in organisms such as Stenotrophomonas maltophilia, Pseudomonas aeruginosa, Klebsiella pneumoniae, and Escherichia coli (1, 2). Each year, about 2.8 million people are diagnosed with bacterial infections in the United States who are resistant to one or two types of antibiotics, leading to 35,000 deaths (3). Bacteria become resistant to antibiotics by mutations, and bacterial resistance manifests a vital need for novel antimicrobial agents (4, 5). Most bacterial resistance is due to the overexpression of β-lactamases, which can hydrolyze and inactivate many β-lactam antibiotics (6, 7). Antibiotic resistance has increased the attention of the World Health Organization in recent years as an important challenge for humans (7). Although β-lactam antibiotics such as carbapenem, cephalosporins, monobactams, and penicillins are currently used, antibiotic resistance necessitates studies to find new compounds with better and broader activity (8-10).

In the last two decades, Metallo-β-lactamases (MBLs) have been identified as major role-players in the mechanism of resistance to carbapenems, which are the last resort of β-lactam antibiotics for the treatment of infections caused by multi-drug-resistant pathogens (10, 11). β-lactamases (BLs) are arranged in two classes, including MBLs (metallo-BLs) containing a hydroxyl ion (as active nucleophile) and one or two zinc ions (in the binding site), and SBLs (serine-BLs) utilizing a catalytic serine for antibiotic hydrolysis (12, 13). MBLs are zinc-dependent enzymes that hydrolyze almost all β-lactams, such as cephalosporins and carbapenems (14). SBLs have Ambler classes (A, C, and D), but MBLs cover class B with three subgroups (B1, B2, and B3). Also, chromosomal and plasmid-mediated genes can encode BLs (7).

The unique structure of the active site, di-zinc coordination, and water bridge hydrolysis mechanism enable subclass B1 to inactivate a variety of substrates (15). The β-lactam ring of antibiotics can be hydrolyzed by SBLs based on the serine residue through an enzyme-acyl intermediate around the active site. On the other hand, MBLs can catalyze the hydrolysis of the β-lactam ring by utilizing the Zn ions (one or two) in the active site (16, 17). The most common clinical problems in MBLs are with the B1 subclass, including VIM (Verona integrin-encoded Metallo-β-lactamase), NDM (New Delhi Metallo-β-lactamase), and IMP (IMiPenamase) (18, 19). VIM is currently the most common among the B1 subclass of MBLs in P. aeruginosa, and the primary source of global prevalence is VIM-2 (20). Developmental studies have critical roles in finding new inhibitors to overcome MBLs and treat bacterial resistance (21). In parts of Asia, VIM exists in more than 99% of patients with MBL-positive multidrug-resistant strains (22).

Captopril ((D-3-mercapto-2-methylpropranoyl-L-proline) found by Bristol-Myers in Boston) has been used for several decades as a suitable drug for high blood pressure prevention. Captopril also can inhibit different MBLs of all subclasses (23, 24). In the present study, captopril was selected as a positive control for the inhibition of VIM-2 MBL. In addition, virtual screening was chosen to find strong and new inhibitors for time and cost-saving purposes. Also, molecular dynamics simulations were performed by the GROMACS package.

2. Objectives

The purpose of this investigation was to recognize a potent inhibitor for VIM-2 MBL and compare it with captopril as a reference inhibitor.

3. Methods

3.1. Molecular Docking

Natural compounds from the ZINC databank (25) in mol2 format containing 61953 structures were achieved. The VIM-2 MBL crystal structure (ID: 4C1E) in the PDB format was obtained from the RCSB (Protein Data Bank) (23). Both proteins and compounds were prepared, and Maestro 12.5 (26) was employed for the docking study. In addition, the ADME (absorption, distribution, metabolism, and excretion) calculation was performed using an online platform (27). Captopril was chosen as the positive control for docking and molecular dynamics (MD) simulation (23).

3.2. Molecular Dynamic

In the next step, two complexes of ligand-protein (MBL-ZINC00517765 and MBL-Captopril) were extracted from the results of the docking study for MD simulation. The ACPYPE server (28) was used to produce the topology parameters file for compounds. GROMACS software (Amber99sb force field) and TIP3P water model with cubic box were applied to simulate the systems (29, 30). The physiological condition (the salt ion environment (NaCl 0.15M)) was fixed, and the charge was neutralized for three systems by adding Na+ or Cl- (9 Na, 11 Na, 10 NA ions to free MBL, MBL-ZINC00517765, MBL-Captopril, respectively) to the solvent (31). The temperature and pressure were determined by employing NPT and NVT ensembles (1 bar with Parrinello-Rahman barostat and 310K using a v-rescale thermostat) (32). After simulation, the radius of gyration (Rg) and principal component analysis (PCA) were employed to investigate the flexibility and motion, and root mean squared fluctuation (RMSF) with root mean square deviation (RMSD) were used to estimate the fluctuations and H-bond (inter-molecular hydrogen bonds). The simulation time for all systems was 120 nanoseconds.

3.3. Binding Free Energy

Molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) is a reliable method for calculating the binding free energy (ΔGbind). In the present study, the GROMACS Tool (g_mmpbsa) was used to study binding free energy. This approach has been used in computational drug design as an appropriate technique for measuring binding free energy (33-35). The following equation expresses the binding free energy:

G of binding = G of complex – (G of protein + G of ligand)

Where the free energy of a state (L denotes ligand and P represents protein) was computed through processing ligand-receptor conformations (36-40).

4. Results and Discussion

4.1. Virtual Screening

The docking study by Maestro 12.5 software determined the power of the new compounds based on inhibitory potential in the active site of VIM-2 MBL. The top nine compounds based on docking energy including ZINC00517765, ZINC04023122, ZINC12659408, ZINC12889936, ZINC01531046, ZINC02108368, ZINC72326107, ZINC20113415, and ZINC10475480 are listed in Figure 1. Molecular docking results (docking energy) for top structures were between -12.29 kcal mol-1 and -10.12 kcal mol-1. Docking's study displayed that ZINC00517765 with -12.29 kcal mol-1 was stronger than the positive control (captopril) with -10.12 kcal mol-1 docking energy.

Docking energies of Captopril and the top nine compounds, as well as 2D structures, ZINC IDs, and molecular weight
Docking energies of Captopril and the top nine compounds, as well as 2D structures, ZINC IDs, and molecular weight
Figure 1.

Docking energies of Captopril and the top nine compounds, as well as 2D structures, ZINC IDs, and molecular weight

4.2. ADME Investigations

One of the major problems with the usefulness of new inhibitors is their safety inadequacy. Drug development usually encounters various unforeseen challenges in different phases. The potency of the top inhibitors was measured in terms of pharmacokinetics, lipophilicity, and drug-likeness impacts. Captopril and the nine compounds including ZINC00517765, ZINC04023122, ZINC12659408, ZINC12889936, ZINC01531046, ZINC02108368, ZINC72326107, ZINC20113415, and ZINC10475480 were analyzed. Table 1 shows that ZINC00517765 was in the acceptable region as a drug in terms of absorption, distribution, metabolism, and excretion. Also, captopril features were similar to those of ZINC00517765.

Table 1.

Pharmacokinetics, Lipophilicity, and Drug-likeness Features of the Top Nine Compounds Compared with Captopril

Molecule NameGI AbsorptionLog Kp (Skin Permeation) cm/sLipinskiBioavailability ScoreSynthetic AccessibilityBBB Permeant
ZINC00517765High-6.49Yes; 0 violation0.562.78No
ZINC04023122High-5.83Yes; 0 violation0.852.98Yes
ZINC12659408Low-6.75Yes; 0 violation0.565.04No
ZINC12889936Low-8.50Yes; 1 violation: NorO>100.114.21No
ZINC01531046High-6.63Yes; 0 violation0.851.47No
ZINC02108368Low-7.32Yes; 0 violation0.563.76No
ZINC72326107High-6.85Yes; 0 violation0.563.07No
ZINC20113415High-10.29Yes; 0 violation0.563.84No
ZINC10475480High-8.19Yes; 0 violation0.562.76No
CaptoprilHigh-7.38Yes; 0 violation0.562.47No

4.3. Interactions of VIM-2 MBL with ZINC00517765 Before Simulation

As can be seen in Figure 2, MBL via Asn210 in its active site interacted with captopril by hydrogen bonding. Captopril also interacted with the MBL active site through Zn501 and Zn502. On the other hand, ZINC00517765 formed two hydrogen bonds by interaction with Asn210 from the MBL active site. Furthermore, Tyr67 and Trp87 in the MBL active site formed a pi-pi stacking interaction with ZINC00517765 and Arg205, and Zn501 and Zn502 interacted with the MBL active site, which confirmed that ZINC00517765 was more potent and formed more bonds with the MBL active site.

Interaction of ZINC00517765 and captopril with the VIM-2 MBL active site before simulation
Interaction of ZINC00517765 and captopril with the VIM-2 MBL active site before simulation
Figure 2.

Interaction of ZINC00517765 and captopril with the VIM-2 MBL active site before simulation

4.4. Stability Analysis

After MD simulation to evaluate the stability of the crystal structure and changes during the simulation, RMSD (Figure 3A), Rg (Figure 3B), RMSF (Figure 3C), Intermolecular H-bonds (Figure 3D), and PCA (Figure 3E) analyses were performed on MBL-ZINC00517765, MBL-Captopril, and free MBL. The RMSD results for free MBL, MBL-ZINC00517765, and MBL-Captopril were 2.18, 2.13, and 2.09 Å, respectively. The RMSD results showed that from 50 nm to the end of the simulation, all the systems obtained stable conditions, but the stability of the MBL-ZINC00517765 complex was better than those of the other two systems. Based on the Rg analysis, the protein fluctuated during the simulations, and these changes were fewer for the MBL-ZINC00517765 complex, indicating that this complex had better stability. The Rg values for all simulations were between 16.2 and 16.8 Å, but the Rg value for MBL-ZINC00517765 was between 16.2 and 16.7 Å. In addition, RMSF analysis was performed using the simulation results, which showed that free MBL in the residues Val36 with 2.21 Å, Glu38 with 2.20 Å, Asp63 with 3.62 Å, Tyr67 with 2.32 Å, Trp87 with 2.45 Å, His116 with 1.77 Å, and Val211 with 3.32 had the highest fluctuations. Also, Val41 with 2.20 Å, Arg127 with 1.73 Å, Arg141 with 2.75 Å, Cys198 with 1.81 Å, and Leu216 with 2.37 Å had the most fluctuations in the MBL-Captopril system. The RMSF analysis demonstrated that MBL-ZINC00517765 in the residues including Arg109 with 1.38 Å, Glu156 with 2.24 Å, Ser159 with 2.41 Å, Ser161 with 2.92 Å, and Glu232 with 1.98 Å had the most fluctuations. Inter-molecular H-bond analysis revealed that during simulation, ZINC00517765 formed more hydrogen bonds with the active site of MBL than did Captopril. The PCA displayed the movements of proteins. The PCA outcomes showed that MBL-ZINC00517765 had less motion than the other two systems and confirmed the results of other analyses, including Rg, intramolecular H-bonds, RMSF, and RMSD.

RMSD (A), Rg (B), RMSF (C), Intermolecular H-bonds (D), and PCA (E) values for MBL-ZINC00517765 (blue), MBL-Captopril (dark red), and free MBL (yellow)
RMSD (A), Rg (B), RMSF (C), Intermolecular H-bonds (D), and PCA (E) values for MBL-ZINC00517765 (blue), MBL-Captopril (dark red), and free MBL (yellow)
Figure 3.

RMSD (A), Rg (B), RMSF (C), Intermolecular H-bonds (D), and PCA (E) values for MBL-ZINC00517765 (blue), MBL-Captopril (dark red), and free MBL (yellow)

4.5. Free Binding Energy

After simulation, the binding free energy was computed for the MBL-ZINC00517765 and MBL-Captopril complexes using MM-PBSA software. According to the RMSD chart, the binding free energy was assessed at the end of the simulation (with minimum fluctuations), and the results exhibited that MBL-ZINC00517765 with -72.29 KJ.mol-1 was stronger than MBL-Captopril with -23.39 KJ. mol-1. Also, the values of electrostatic energy and Van der Waals energy were better for MBL-ZINC00517765 than for MBL-Captopril (Table 2).

Table 2.

Molecular Energy Values of MBL-ZINC00517765 and MBL-Captopril Complexes

Energy (KJ.mol-1)MBL-CaptoprilMBL-ZINC00517765
ΔESASA-17.92 ± 0.27-20.93 ± 0.58
ΔEsolv253.69 ± 15.17275.57 ± 10.53
ΔEelect-173.77 ± 7.35-215.81 ± 8.36
ΔEvdw-85.39 ± 4.81-111.12 ± 8.85
ΔGbinding-23.39 ± 5.92-72.29 ± 10.74

5. Conclusions

According to the results of molecular dynamics, virtual screening, and docking studies, ZINC00517765 can inhibit MBL. The results showed that MBL in the presence of ZINC00517765 has stable conditions and is better than Captopril. Finally, further study on ZINC00517765 is recommended to overcome bacterial infections as a suitable inhibitor. We hope that the results of this study will help researchers to find a drug for MBL.

References

  • 1.

    Bush K. Proliferation and significance of clinically relevant beta-lactamases. Ann N Y Acad Sci. 2013;1277:84-90. [PubMed ID: 23346859]. https://doi.org/10.1111/nyas.12023.

  • 2.

    Walsh TR. Emerging carbapenemases: a global perspective. Int J Antimicrob Agents. 2010;36 Suppl 3:S8-14. [PubMed ID: 21129630]. https://doi.org/10.1016/S0924-8579(10)70004-2.

  • 3.

    Thomas CA, Cheng Z, Yang K, Hellwarth E, Yurkiewicz CJ, Baxter FM, et al. Probing the mechanisms of inhibition for various inhibitors of metallo-beta-lactamases VIM-2 and NDM-1. J Inorg Biochem. 2020;210:111123. [PubMed ID: 32622213]. https://doi.org/10.1016/j.jinorgbio.2020.111123.

  • 4.

    Zhang J, Wang S, Wei Q, Guo Q, Bai Y, Yang S, et al. Synthesis and biological evaluation of Aspergillomarasmine A derivatives as novel NDM-1 inhibitor to overcome antibiotics resistance. Bioorg Med Chem. 2017;25(19):5133-41. [PubMed ID: 28784300]. https://doi.org/10.1016/j.bmc.2017.07.025.

  • 5.

    Darbani R, Farshadfar C, Tavana S, Saljoughi H, Zonouri SS. Identification of DNA gyrase Subunit a Mutations Associated with Ciprofloxacin Resistance in Staphylococcus aureus Isolated from Nasal Infection in Kurdistan-Iran. J Mol Biol Res. 2017;7(1):186. https://doi.org/10.5539/jmbr.v7n1p186.

  • 6.

    Frere JM. Beta-lactamases and bacterial resistance to antibiotics. Mol Microbiol. 1995;16(3):385-95. [PubMed ID: 7565100]. https://doi.org/10.1111/j.1365-2958.1995.tb02404.x.

  • 7.

    Gianquinto E, Tondi D, D'Arrigo G, Lazzarato L, Spyrakis F. Can We Exploit beta-Lactamases Intrinsic Dynamics for Designing More Effective Inhibitors? Antibiotics (Basel). 2020;9(11):833. [PubMed ID: 33233339]. [PubMed Central ID: PMC7700307]. https://doi.org/10.3390/antibiotics9110833.

  • 8.

    Walsh C. Molecular mechanisms that confer antibacterial drug resistance. Nature. 2000;406(6797):775-81. [PubMed ID: 10963607]. https://doi.org/10.1038/35021219.

  • 9.

    Page MG. Cephalosporins in clinical development. Expert Opin Investig Drugs. 2004;13(8):973-85. [PubMed ID: 15268635]. https://doi.org/10.1517/13543784.13.8.973.

  • 10.

    Meini MR, Llarrull LI, Vila AJ. Overcoming differences: The catalytic mechanism of metallo-beta-lactamases. FEBS Lett. 2015;589(22):3419-32. [PubMed ID: 26297824]. [PubMed Central ID: PMC4640939]. https://doi.org/10.1016/j.febslet.2015.08.015.

  • 11.

    Walsh TR, Toleman MA, Poirel L, Nordmann P. Metallo-beta-lactamases: the quiet before the storm? Clin Microbiol Rev. 2005;18(2):306-25. [PubMed ID: 15831827]. [PubMed Central ID: PMC1082798]. https://doi.org/10.1128/CMR.18.2.306-325.2005.

  • 12.

    Livermore DM. beta-Lactamases in laboratory and clinical resistance. Clin Microbiol Rev. 1995;8(4):557-84. [PubMed ID: 8665470]. [PubMed Central ID: PMC172876]. https://doi.org/10.1128/CMR.8.4.557.

  • 13.

    Bush K, Jacoby GA, Medeiros AA. A functional classification scheme for beta-lactamases and its correlation with molecular structure. Antimicrob Agents Chemother. 1995;39(6):1211-33. [PubMed ID: 7574506]. [PubMed Central ID: PMC162717]. https://doi.org/10.1128/AAC.39.6.1211.

  • 14.

    Hinchliffe P, Tanner CA, Krismanich AP, Labbe G, Goodfellow VJ, Marrone L, et al. Structural and Kinetic Studies of the Potent Inhibition of Metallo-beta-lactamases by 6-Phosphonomethylpyridine-2-carboxylates. Biochemistry. 2018;57(12):1880-92. [PubMed ID: 29485857]. [PubMed Central ID: PMC6007964]. https://doi.org/10.1021/acs.biochem.7b01299.

  • 15.

    Shen B, Yu Y, Chen H, Cao X, Lao X, Fang Y, et al. Inhibitor discovery of full-length New Delhi metallo-beta-lactamase-1 (NDM-1). PLoS One. 2013;8(5). e62955. [PubMed ID: 23675445]. [PubMed Central ID: PMC3652859]. https://doi.org/10.1371/journal.pone.0062955.

  • 16.

    Fisher JF, Meroueh SO, Mobashery S. Bacterial resistance to beta-lactam antibiotics: compelling opportunism, compelling opportunity. Chem Rev. 2005;105(2):395-424. [PubMed ID: 15700950]. https://doi.org/10.1021/cr030102i.

  • 17.

    Wang X, Yang Y, Gao Y, Niu X. Discovery of the Novel Inhibitor Against New Delhi Metallo-beta-Lactamase Based on Virtual Screening and Molecular Modelling. Int J Mol Sci. 2020;21(10). [PubMed ID: 32443639]. [PubMed Central ID: PMC7279046]. https://doi.org/10.3390/ijms21103567.

  • 18.

    Docquier JD, Mangani S. An update on beta-lactamase inhibitor discovery and development. Drug Resist Updat. 2018;36:13-29. [PubMed ID: 29499835]. https://doi.org/10.1016/j.drup.2017.11.002.

  • 19.

    Aitha M, Marts AR, Bergstrom A, Moller AJ, Moritz L, Turner L, et al. Biochemical, mechanistic, and spectroscopic characterization of metallo-beta-lactamase VIM-2. Biochemistry. 2014;53(46):7321-31. [PubMed ID: 25356958]. [PubMed Central ID: PMC4245990]. https://doi.org/10.1021/bi500916y.

  • 20.

    Persoon MC, Voor In't Holt AF, Wielders CCH, Gommers D, Vos MC, Severin JA. Mortality associated with carbapenem-susceptible and Verona Integron-encoded Metallo-beta-lactamase-positive Pseudomonas aeruginosa bacteremia. Antimicrob Resist Infect Control. 2020;9(1):1-8. [PubMed ID: 32014058]. [PubMed Central ID: PMC6998810]. https://doi.org/10.1186/s13756-020-0682-4.

  • 21.

    Farina D, Spyrakis F, Venturelli A, Cross S, Tondi D, Costi MP. The inhibition of extended spectrum beta-lactamases: hits and leads. Curr Med Chem. 2014;21(12):1405-34. [PubMed ID: 24180276]. https://doi.org/10.2174/09298673113206660323.

  • 22.

    Edelstein MV, Skleenova EN, Shevchenko OV, D'Souza J W, Tapalski DV, Azizov IS, et al. Spread of extensively resistant VIM-2-positive ST235 Pseudomonas aeruginosa in Belarus, Kazakhstan, and Russia: a longitudinal epidemiological and clinical study. Lancet Infect Dis. 2013;13(10):867-76. [PubMed ID: 23845533]. https://doi.org/10.1016/S1473-3099(13)70168-3.

  • 23.

    Brem J, van Berkel SS, Zollman D, Lee SY, Gileadi O, McHugh PJ, et al. Structural Basis of Metallo-beta-Lactamase Inhibition by Captopril Stereoisomers. Antimicrob Agents Chemother. 2016;60(1):142-50. [PubMed ID: 26482303]. [PubMed Central ID: PMC4704194]. https://doi.org/10.1128/AAC.01335-15.

  • 24.

    Vidt DG, Bravo EL, Fouad FM. Medical intelligence drug therapy: captopril. N Engl J Med. 1982;306(4):214-9. [PubMed ID: 7033784]. https://doi.org/10.1056/NEJM198201283060405.

  • 25.

    Sterling T, Irwin JJ. ZINC 15--Ligand Discovery for Everyone. J Chem Inf Model. 2015;55(11):2324-37. [PubMed ID: 26479676]. [PubMed Central ID: PMC4658288]. https://doi.org/10.1021/acs.jcim.5b00559.

  • 26.

    Beard H, Cholleti A, Pearlman D, Sherman W, Loving KA. Applying physics-based scoring to calculate free energies of binding for single amino acid mutations in protein-protein complexes. PLoS One. 2013;8(12). e82849. [PubMed ID: 24340062]. [PubMed Central ID: PMC3858304]. https://doi.org/10.1371/journal.pone.0082849.

  • 27.

    Daina A, Michielin O, Zoete V. SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci Rep. 2017;7:1-13. [PubMed ID: 28256516]. [PubMed Central ID: PMC5335600]. https://doi.org/10.1038/srep42717.

  • 28.

    Sousa da Silva AW, Vranken WF. ACPYPE - AnteChamber PYthon Parser interfacE. BMC Res Notes. 2012;5(1):1-8. [PubMed ID: 22824207]. [PubMed Central ID: PMC3461484]. https://doi.org/10.1186/1756-0500-5-367.

  • 29.

    Ardalan N, Mirzaie S, Sepahi AA, Khavari-Nejad RA. Novel mutant of Escherichia coli asparaginase II to reduction of the glutaminase activity in treatment of acute lymphocytic leukemia by molecular dynamics simulations and QM-MM studies. Med Hypotheses. 2018;112:7-17. [PubMed ID: 29447943]. https://doi.org/10.1016/j.mehy.2018.01.004.

  • 30.

    Ardalan N, Akhavan Sepahi A, Khavari-Nejad RA. Development of Escherichia coli asparaginase II for the Treatment of Acute Lymphocytic Leukemia: In Silico Reduction of asparaginase II Side Effects by a Novel Mutant (V27F). Asian Pac J Cancer Prev. 2021;22(4):1137-47. [PubMed ID: 33906306]. [PubMed Central ID: PMC8325130]. https://doi.org/10.31557/APJCP.2021.22.4.1137.

  • 31.

    Yang B, Lin SJ, Ren JY, Liu T, Wang YM, Li CM, et al. Molecular Docking and Molecular Dynamics (MD) Simulation of Human Anti-Complement Factor H (CFH) Antibody Ab42 and CFH Polypeptide. Int J Mol Sci. 2019;20(10). [PubMed ID: 31130605]. [PubMed Central ID: PMC6566401]. https://doi.org/10.3390/ijms20102568.

  • 32.

    Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J Chem Theory Comput. 2008;4(3):435-47. [PubMed ID: 26620784]. https://doi.org/10.1021/ct700301q.

  • 33.

    Koushki EH, Abolghasemi S, Mollica A, Aghaeepoor M, Moosavi SS, Farshadfar C, et al. Structure-based virtual screening, molecular docking and dynamics studies of natural product and classical inhibitors against human dihydrofolate reductase. Netw Model Anal Health Inform Bioinform. 2020;9(1):1-21. https://doi.org/10.1007/s13721-020-00244-9.

  • 34.

    Farshadfar C, Mollica A, Rafii F, Noorbakhsh A, Nikzad M, Seyedi SH, et al. Novel potential inhibitor discovery against tyrosyl-tRNA synthetase from Staphylococcus aureus by virtual screening, molecular dynamics, MMPBSA and QMMM simulations. Mol Simul. 2020;46(7):507-20. https://doi.org/10.1080/08927022.2020.1726911.

  • 35.

    Noorbakhsh A, Askandar RH, Alhagh MS, Farshadfar C, Seyedi SH, Ahmadizad M, et al. Prevention of SARS-CoV-2 Proliferation with a Novel and Potent Main Protease Inhibitor by Docking, ADMET, MM-PBSA, and Molecular Dynamics Simulation. J Comput Biol Chem. 2021;20(3):305-22. https://doi.org/10.1142/s2737416521500149.

  • 36.

    Wang C, Greene D, Xiao L, Qi R, Luo R. Recent Developments and Applications of the MMPBSA Method. Front Mol Biosci. 2017;4:87. [PubMed ID: 29367919]. [PubMed Central ID: PMC5768160]. https://doi.org/10.3389/fmolb.2017.00087.

  • 37.

    Kumari R, Kumar R, Open Source Drug Discovery C, Lynn A. g_mmpbsa--a GROMACS tool for high-throughput MM-PBSA calculations. J Chem Inf Model. 2014;54(7):1951-62. [PubMed ID: 24850022]. https://doi.org/10.1021/ci500020m.

  • 38.

    Dash R, Choi HJ, Moon IS. Mechanistic insights into the deleterious roles of Nasu-Hakola disease associated TREM2 variants. Sci Rep. 2020;10(1):1-14. [PubMed ID: 32107424]. [PubMed Central ID: PMC7046722]. https://doi.org/10.1038/s41598-020-60561-x.

  • 39.

    Noorbakhsh A, Hosseininezhadian Koushki E, Farshadfar C, Ardalan N. Designing a natural inhibitor against human kynurenine aminotransferase type II and a comparison with PF-04859989: a computational effort against schizophrenia. J Biomol Struct Dyn. 2021:1-14. [PubMed ID: 33645449]. https://doi.org/10.1080/07391102.2021.1893817.

  • 40.

    Shayan S, Jamaran S, Askandar RH, Rahimi A, Elahi A, Farshadfar C, et al. The SARS-Cov-2 Proliferation Blocked by a Novel and Potent Main Protease Inhibitor via Computer-aided Drug Design. Iran J Pharm Res. 2021;20(3):399-418.