In silico approach for potential antimicrobial agents through antiviral, molecular docking, molecular dynamics, pharmacokinetic and bioactivity predictions of galactopyranoside derivatives

Abstract Carbohydrates derivatives played an important and vital role in a large part to design new potential drug against various pathogens, in medicinal chemistry. Also, monosaccharide derivatives are very useful in biological chemistry because monosaccharide binding agents were shown to possess antiviral, antibacterial, and antifungal activity towards HIV, Dengue coronaviruses, and many members of the family Enterobacteriaceae, including strains of Escherichia coli, Salmonella typhimurium, Klebsiella pneumoniae, Shigella flexneri, and Enterobacter aerogenes. The present study was therefore interesting in identifying the binding affinity of a series of monosaccharides previously synthesized through their physicochemical and pharmacokinetic properties. These monosaccharides are derivatives of methyl-β-D-galactopyranoside (MGP) coupled to the protein deacetoxycephalosporin C synthase (DAOCS) of Streptomyces clavuligerus (1RXF). Molecular docking study against DAOCS proteins (PDB ID: 1RXF) was performed assess the antimicrobial activities of these derivatives along with antiviral prediction. The results of molecular docking studies showed that derivatives 8 and 10 both presented the highest score (-7.2 kcal/mol), reflecting their good binding affinity toward the active pocket of the pathogen S. clavuligerus (1RXF). They are found to be potent antimicrobial agents. These findings were further validated in dynamics (200 ns) milieu and intermolecular stability was evaluated via highly acceptable Molecular Mechanics Poisson-Boltzmann Surface Area (MMPB)/Molecular Mechanics/Generalized Born Surface Area (GBSA) binding free energies estimation methods. The dynamics of the complexes are highly stable with the strengthening of docked interactions as highlighted by MMPB/GBSA method (net binding energy < −40 kcal/mol). Density functional theory (DFT) calculations were carried out on the derivatives of MGP in order to determine the partial atomic charges, electronic energies, enthalpies, entropies and polarizability, at the B3LYP/3-21G level. The in silico ADMET prediction carried out provided access to the bioactivity and toxicity parameters which established that most of the different ligands designed are safe to use as drugs. This in silico study provided an overall view of the possible pharmacological and therapeutic profile of MGP derivatives. It therefore, paved the way for the discovery of promising next-generation drugs against pathogenic microbes.


Introduction
Carbohydrates are the most available biomolecules belonging to organic compounds class. They come from living organism's kingdom. Around the world, scientists are investing in medicinal research to find more effective and safer antimicrobial agents for the treatment of diseases caused by pathogenic organisms. For a long time, carbohydrates remained an attractive research subject for scientists due to the important role played in the functioning of biological systems, namely the growth and proliferation of cells, communication between cells and the immune response (Nogueira, Parmanhan, Farias, & Corrêa, 2009). They are also the source of metabolic energy supply and responsible for the fine-tuning of cell-cell interactions and many other important biological processes (Sears & Wong, 1996;Werz, 2007). Other advantages of carbohydrates also pointed out are the numerous and diverse biological activities which they exert against pathogenic germs. In addition, recent advances in the treatment of infectious diseases have established that the use of carbohydrates could contain the effects of resistant races of microbes. In this case, the deep understanding of the antimicrobial activity mechanism is a better way which undoubtedly leads to synthesizing new antimicrobial compounds with improved and effective activity.
Moreover, selective acylation as well as evaluation of microbial activities of carbohydrates revealed that in many cases the combination of two or more hetero-aromatic rings and acyl groups enhanced the biological activity of the parent nucleus (Kawsar, Faruk, Rahman, Fujii, & Ozeki, 2014;Kawsar et al., 2013;Misbah et al., 2020). Monosaccharide derivatives are known for their broad spectrum of biological activity against both Gram negative and positive organisms such as Escherichia coli, Bacillus subtilis, Salmonella typhi and Staphylococcus aureus (Islam et al., 2019;Kawsar et al., 2013;2014;Misbah et al., 2020). In a recent study, some monosaccharide analogs were shown to exhibit remarkable inhibitory activity against cancer cell (Kawsar, Hosen, Fujii, Ozeki, & Thermochemical, 2020). Modifications to the hydroxyl group level in the nucleoside and monosaccharide structure have been shown very useful in the design of potent antiviral (M. Maowa, Asraful, et al., 2021;Yasmin et al., 2021a) and antimicrobial agents (Kawsar & Hossain, 2020;. Thus, taking into account the observations and results mentioned above, we have chosen to substitute the hydrogen atom of the hydroxyl group of MGP with acyl substituents (including aromatic and hetero-aromatic group). Quantum chemistry methods are then used to determine some of electronic and thermochemical characteristics of the various resulting structures. Likewise, molecular docking studies were carried out to investigate the possible binding mechanism of MGP derivatives with DAOCS protein S. clavuligerus (PDB: 1RXF). The docking predictions on the binding mode of the compounds and interactions as not much reliable were thoroughly validated by dynamics stability in molecular dynamics simulations. Also, intermolecular strength of the complexes was confirmed through MMPB/GBSA analysis. On the other hand, and pharmacokinetic study using in silico ADMET prediction methods in order to demonstrate their bioavailability and toxicity, on the other hand. In general, the objective is to have, for MGP and derivatives, an in-depth knowledge of the in silico computational and pharmacokinetic properties.

Geometry optimization
In computational chemistry, quantum mechanical methods are widely used to calculate thermal, molecular orbital and molecular electrostatic properties (Berman et al., 2000). Geometry optimization of all the synthesized analogs is carried out using Gaussian 09 program (Frisch, 2009). Density functional theory (DFT) with Beck's (B) (Becke, 1988) three parameter hybrid model, Lee, Yang and Parr's (LYP) (Lee et al., 1988) correlation functional with 3-21 G basis set are employed to optimize and to predict their thermal and physicochemical properties. Parameters such as enthalpy, entropy, heat capacity, total energy, partial atomic charge and polarizability were calculated for all studied compounds.

Antiviral activities determination
Antiviral molecules (AVMs) are a special class of antimicrobial drugs used for the treatment of viral infections. They react by inhibiting the growth of the viral pathogen(s) inside the host cell. For the calculation of antiviral activity, we used the AVCPred program (Qureshi, Kaur, & Kumar, 2017) which displayed the inhibition percentage.

Preparation of protein and molecular docking
The molecular docking is one of the best and suitable approaches to identify the drug interactions with the protein DAOCS of Streptomyces clavuligerus. At this level, there is the blind docking method which consists of searching the whole surface of the protein molecule for binding sites.
The 3 D structure of DAOCS (PDB ID. 1RXF) ( Figure  1) was collected in pdb format from the protein data bank (Helen M Berman et al., 2000). All hetero-atoms and water molecules were took away by using PyMol (version 1.3) software packages (DeLano, 2002). Swiss-Pdb viewer software (version 4.1.0) was employed for energy minimization of the protein (Guex & Peitsch, 1997) . Then optimized drugs were subjected from molecular docking study against S. clavuligerus protein (1RXF). In fine, molecular docking simulation was rendered by PyRx software (version 0.8) (Dallakyan & Olson, 2015) considering the protein as the macromolecule and the drug as the ligand. AutodockVina was used for docking analysis, and AutoDock Tools (ADT) of the MGL software package was used to convert pdb format into a Pdbqt format to input protein and ligands. The size of grid box in AutoDockVina was kept at 57.96, 48.57, and 47.22 Å for X, Y, Z directions respectively. Both macromolecule and ligand structures in Pdbqt format after the docking were used under Accelrys Discovery Studio (version 4.1) program to explore and visualize the result of docking and the search for non-bonding interactions between ligands and amino acid residues of receptor protein (Version, 2017). PDBsum program was also used to check the validation of the protein (PDB: 1RXF) following the scheme of Ramachandran (Figure 1). It revealed that 93.20% of residues in the allowed region and no residues were missed.

Validation protocol of docking
The biological validation of docking outcomes were investigated by extracting the co-crystallized ligand of the fungal protein (PDB ID: 1RXF) and re-docked it maintaining the same position. The ground state energy pose attained on re-docking and the co-crystallized ligands were superimposed, and its root mean square deviation (RMSD) score was calculated between these two superimposed ligands poses. To validate the docking procedure, The RMSD score range must be within the 2 Å (Onodera, Satou, & Hirota, 2007;Warren et al., 2006).

Molecular dynamics simulations
The selected complexes of compound 8 and 10 were prepared through LEAP module of AMBER20 (Case et al., 2020). First, the complexes were submerged in TIP3PBOX allowing a 12.0 Å padding distance between the protein and box edge. Subsequently, sufficient number of counter ions were used to neutralize charge of each system. The ANTECHAMBER program (Wang, Wang, Kollman, & Case, 2001) of the AMBER was employed to parameterize the compounds via GAFF force field (Wang, Wolf, Caldwell, Kollman, & Case, 2004), while FF14SB force field (Maier et al., 2015) was applied to the receptor protein. All atoms simulations of the complexes were performed after hydrogen atoms position optimization, water molecules optimization and entire systems minimization. The simulations were carried out with time scale of 2 fs and cutoff radius of 8.0 Å for non-bounded interactions. The long-range electrostatic interactions are calculated using particle-mesh Ewald method (Petersen, 1995) and hydrogen bonds were constrained via SHAKE algorithm (Kr€ autler, van Gunsteren, & H€ unenberger, 2001). Systems heating was achieved in gradual phase from 0 to 300 K. Temperature control was accomplished through Langevin algorithm (Eisenmann & Stillfjord, 2021). Systems were then equilibrated for 100 ps at 1 atm pressure and 300 K. The same set of conditions were applied in 100 ns of simulations using PMEMD.CUDA from AMBER20. The simulations trajectories were investigated by CPPTRAJ module (Roe & Cheatham, 2013). The systems intermolecular docked conformation stability was examined by performing carbon alpha root mean square deviation (RMSD) (Kirchmair, Markt, Distinto, Wolber, & Langer, 2008) and root mean square fluctuation (RMSF) (Sahu & Pattanayak, 2019).

MMPB/GBSA binding free energies
The intermolecular interactions energy based on simulation trajectories was estimated using MMPB/ GBSA method (Wang et al., 2019) by considering MMPBSA.py module (Miller et al., 2012) of the AMBER20. The binding energy analysis takes into account explicitly computed parameters involving gas phase energy, conformation entropy energy and solvation energy. Details of the MMPB/GBSA methodology applied herein are presented in Haq, Abro, Raza, Liedl, & Azam, 2017(Haq et al., 2017. In total, 500 frames where the simulation trajectories were selected at regular time scale and used for the net binding free energy estimation.

Physicochemical pharmacokinetic analysis
The admetSAR 2.0 (http://lmmd.ecust.edu.cn/admet-sar2/about) server was used to evaluate the pharmacokinetics and toxicity parameters related to drug absorption, metabolism and toxicity of the parent drug and its new designed analogs (Cheng et al., 2012). Using structure similarity search methods, ADMETSAR predicts the latest and most comprehensive manually accurate data for diverse chemicals associated with known ADME/T profiles. Molinspiration program (https://www.molinspiration.com/cgi-bin/properties) (Ali & Alam Khan, 2013) was employed to analyze the drug-likeness properties (Lipinski, Lombardo, Dominy, & Feeney, 1997) of lead compounds.

Results and discussion
In the present study, eleven derivatives of MGP obtained by proceeding by structural modification with aliphatic and aromatic chains (2-10) ( Table 1) were considered. Geometry optimization and antiviral activity prediction were performed to highlight the mode of their antimicrobial behavior. Observed activities were then rationalized by calculating parameters such as thermodynamic parameters, partial atomic charge, through molecular docking, with the combination in silico pharmacokinetic and druglikeness properties.

Structural designation of designed MGP derivatives
Atomic identification and structural variation of substituted MGP derivatives are displayed in Table 1. Different aliphatic (acetyl, butyryl, palmitoyl, myristoyl and steroyl) and aromatic (trityl, 3-chlorobenzoyl, cinnamoyl, and p-toluenesulfonyl) groups were subjected to substitute the hydrogen atom of MGP hydroxyl (-OH) group in order to observe the variation of biological activities.

FTIR spectrum
One of the most commonly employed methods that have been applied to detect organic and polymeric elements is Fourier transform infrared analysis (FTIR). This approach employs infrared light to inspect test materials and assess their chemical characteristics through using the FTIR analysis method. The FTIR system transmits infrared light of approximately 10,000 to 100 cm À1 through a material, with some energy consumed and some escaping through. Signals from analyzers normally range from 4000 cm À1 to 400 cm À1 which reflect the molecular fingerprint of the sample. The FTIR spectrum of MGP derivatives 2-10 represent in Figure 2.

Antiviral activity prediction
In order to check the antiviral mode of action of inserted aliphatic and aromatic groups, prediction of antiviral activities of studied carbohydrate derivatives (2-10) have been decided. Then predicted antiviral activities are compared to that of Azidothymidine (AZT, antiviral drug) using an antiviral application (Table 2) (Qureshi et al., 2017). The predicted antiviral activities revealed that the derivatives of MGP (2-10) possess potential antiviral efficacy when compared to their parent molecule. The aliphatic derivatives (2-6) and aromatic derivatives (7-10) exhibited promising scores in comparison with standard drugs Azidothymidine (AZT).

Physicochemical analysis
A simple alteration of chemical structure significantly influences the structural properties including thermal and molecular orbital properties. Spontaneity of a reaction and stability of a product can be known from calculating its free energy and enthalpy values respectively (Cohen & Benson, 1993). Highly negative values are more eligible for gaining thermal stability. In drug design, hydrogen bond formation and non-bonded interactions are also influenced by dipole moment.
Electronic energy is a significant criterion to represent the interaction of binding partners, where negative value favors spontaneous binding and interaction. In the present study, the calculated properties values of all the MGP derivatives are greater negative than those of the parent MGP. Hence, these observations indicated that the attachment of ester group could improve the interaction and binding of these molecules with different microbial enzymes. As shown in Table 3, physicochemical data revealed that, among all the synthesized derivatives, only compound (6) with the highest molecular weight showed the highest value for polarizability (524.901a.u), Heat capacity (309.642 cal/mol-kelvin) and Entropy (425.076 cal/mol-kelvin). All the derivatives with aromatic ring as substituent displayed higher score comparatively to the aliphatic one. Finally, this discussion proves that the substitution carried out at the level of hydroxyl (-OH) groups of MGP increases significantly its thermodynamic properties which indicate the inherent stability of the synthesized derivatives. Moreover, presence of bulky acylating groups also suggesting the possible improvement of these physicochemical parameters.

Atomic partial charge analysis: NBO
The partial atomic charge is essential for such purposes in molecular computations as a simplified representation of global charge distribution in a molecule and predicting its conformational behavior. NBO analysis is carried out in order to know the charge distribution and the intrinsic property of the interactions in the designed structure. Polarity of chemical bonds often influences the structure and reactivity of a molecule (Heinz & Suter, 2004). The molecular dipole moment is a vector, which does not clearly define the polarity of the molecule. Different methods have been proposed for assigning partial charges to the interacting atoms within a molecule. The approaches of Mulliken Population Analysis and Natural Bond Orbital were employed to compute the partial charges of all drug interacting atoms (Gross, Seybold, & Hadad, 2002). They are the most used population analysis methods. Dipole moment and molecular polarizability are related to atomic charges (Mulliken, 1955). Here, all the hydrogen atoms showed positive charge in both methods and other electronegative atom (O), negative charge in both methods as expected (Figure 3). In derivative 5 (C-1, C-2 and C-25) showed greater positive charges. These greater positive charges are found due to the presence of higher electronegative element like oxygen (O-9, O-21, and O-33), and hydrogen H-7 which exhibit the higher positive value than the other hydrogen because of oxygen atom of hydroxyl group. Similarly, the carbon atoms identified in the derivatives 6 (C-1, C-6, and C-25), 7 (C-1, C-25), and 8 (C-5, C-14) displayed positive charge in both methods due to the presence of the oxygen atom of carbonyl group. In all derivatives, the carbon (C-25) showed the highest positive charge in all the other drugs. This fact is due to the presence of greater electronegative atom (N, Cl and O). As derivatives (2-6) contain aliphatic substituents and particularly derivatives 5 and 6 possess the longest carbon chain. It was found that these two compounds showed results that are improved. Finally, we depicted derivative 5 and 6 ( Figure 3) from both aliphatic and aromatic series, as these compounds exhibited clear graphical view while the rest of the derivatives gave slightly scribble view.

Validation study
The identification of the molecular docking algorithms to confirm the conformation of the protein-bound ligand, performed by re-docking of the co-crystallized ligand pose which revealed for the validation effectiveness of the docking results. The superimposed pose between the docked ligand conformation and the co-crystallized ligand conformation revealed the RMSD value is < 2 Å. The complex was then observed to interact with the all similar amino acid residues compared to the ones reported in the current study. The large symmetric structure   can be interacted in the active site of protein during docking, as the case in this study; the RMSD score would be at a very high level. Again, the small molecules can attain low RMSD score easily even when imposed randomly. Some reported works have found a new strategy for the efficacy of docking poses on the basis of visual realization. The outcomes of this visual observation showed the same interactions as in the experimental binding mode, as found in docking.

Molecular docking simulation
Molecular docking is a type of bioinformatics modeling which deals with the interaction of two or more molecules to give the stable adduct. In Table 4 and Table S1 the derivatives of MGP (2-5, 8-10) showed excellent docking performance and the chosen proteins were discovered to communicate with significant amino acids. We docked totally 9 complexes having better antimicrobial activities for one protein with eight MGP derivatives. We choose the deacetoxycephalosporin c synthase (DAOCS) protein (PDB: 1RXF) from S. clavuligerus bacteria and all designed derivatives were docked into the same binding pocket of DAOCS using similar optimized docking conditions to investigate their binding mode. The outcomes of the docking analysis showed that obtained binding affinities are ranging from À5.2 to À7.2 kcal/mol for all derivatives and parent compound. As shown in Table 4, all MGP derivatives (2-6 and 8-10) displayed higher binding affinity than parent MGP. These results indicated that insertion of a long carbon chain/aromatic ring allow the molecule to improve binding affinity, while addition of hetero groups like -Cl and -CH 3 made some fluctuations in binding affinities. However, modification with halogenated aromatic rings increased the binding affinity. The docked pose clearly showed that the drugs molecules bind within the active site of the DAOCS macromolecular structure (Figure 4). Comparatively the aromatic analogs displayed better binding score than aliphatic analogs. The interactions between inhibitors and bordering residues of DAOCS are illustrated in the 2 D schematics. These pictures were obtained by importing them from docking output to the Discovery Studio Visualizer and lig-plot ( Figures 5 and 6 and S1-S2). Obtained results show  that amino acids participate in the pattern of interactions between the ligand and protein with an important contribution to the total energy of interaction.
Among aliphatic derivatives, compounds 6 displays binding score which is fluctuating more than other aliphatic analogs with new binding sites. Aromatic derivative (8-10) displayed higher fluctuating binding score (-7.2, À6.9 and À7.2 kcal/mol) compared to aliphatic derivatives. These derivatives are also interacting with catalytic active site of the bacterial protein through hydrogen bond (Arg160, Arg162, Arg306, Asn304, Ile305, Pro71, Tyr184 Asp284 Thr308, and Phe283), hydrophobic bond (Arg307, Arg74, Ile192, His243, Val262 and Lys239) and electrostatic bond (Arg307, Arg303). Pi-cation  and pi-sigma interactions are found only at the level of aromatic derivatives. So, these outcomes show clearly that the fact of having high electron density, aromatic substituents can easily increase the binding skill as well as the antibacterial power of MGP derivatives. According to Phe264, all the derivatives displayed the maximum p-interactions with the His243, His183, Arg306And Arg307 denoting the tight binding with the active site. Reports suggest that phenylalanine and histidine residues are considered as the principal component of the p-alkyl, p-sigma and p-cation. They are responsible for the accessibility of small molecules to the active site. Binding energy and binding mode were improved in derivatives (2-5 and 8-10) because of significant hydrogen bonding. It was observed that the alterations of the -OH group in MGP exalted the p-p interactions with the amino acid chain on the binding site, whereas their polarity improvement resulted in the formation of hydrogen bond interactions.
The most significant H-bonds were obtained for the derivative (3 and 9), forming with identical residues such as Arg306, Arg160, Tyr184, Tyr302, Asn304, Phe283, Asp284, Thr72 and Ile305. Hydrogen-bonds execute a vital function in shaping the specificity of ligand binding with the receptor, drug design in chemical and biological processes, molecular recognition and biological activity (Perlstein, 2001).
Hydrogen bond surface and hydrophobic surface of derivative (3) is represented in Figure 7. We realize that the analyzed MGP derivatives bind within the active site of the DAOCS protein of S. clavuligerus, which is inevitable to prevent the protein from its infections. The calculated binding affinities varied from À5.2 to À7.2 kcal/mol which suggesting that the molecules can spontaneously interact within the binding site of bacterial protein. The docking findings therefore indicate that the derivatives 8 and 10 may possess excellent effectiveness for S. clavuligerus.

Molecular dynamics simulations studies
Molecular dynamics simulations were conducted in order to decipher selected compounds docked conformation stability with the proteins as well as to determine the strength of interactions between the molecules. It can be seen in Figure 8, that both selected docked complexes are depicting stable dynamics behavior: the DAOCS-compound 8 complex is relatively more stable than that of compound 10. The latter system experienced some deviations between times 10 ns to 65 ns however, received considerable stable towards simulation end time. RMSD for both   systems are around 2.5 Å. The RMSD was further extended to validate stability of systems (Figure 8). DAOCS-compound 8 complex was unveiled to gain better stability compared to DAOCS-compound 10 in the last 100 ns. Nevertheless, both systems revealed good in-termolecular binding conformation and interaction throughout the length of simulation time. Next, RMSF was calculated, which concluded that majority of the receptor residues are within acceptable flexibility range with RMSF value < 2 Å (Figure 9). The terminal regions of the proteins showed some fluctuations though that were to minor to effect overall systems stability.

MMPB/GBSA binding free energy analysis
The MMPB/GBSA analysis was performed using snapshots extracted from simulation trajectories to get validation on the docking predictions. Likewise, the docking studies the binding free energy analysis by MMPB/GBSA method also reported both the complexes as excellent stable complexes as demonstrated by negative net binding free energy that is À47.5 kcal/mol for DAOCS-compound 8 complex in MM/GBSA and À41.49 kcal/mol in MM/PBSA. While for DAOCS-compound 10 complex, the net MM/ GBSA binding free energy is À46.71 kcal/mol and MM/PBSA binding free energy is À39.24 kcal/mol. These energy values are four times more than that predicted by the docking study which may be the fact that as multiple dynamics conformations of the complexes are selected compared to just one docked conformation in docking studies. Also, it indicates that as simulation proceeds the binding affinity of the compounds get more strengthen. Details of the different energy parameters contribution to the net binding energy values is given in Table 5.

Pharmacokinetic and drug-likeness properties
Molecular parameters like membrane permeability and bioavailability of lead molecules depend on some fundamental features of molecules like partition coefficient (logP), molecular weight (MW), and number of hydrogen bond acceptors/donors that is associated with Lipinski ''rule of five''. Table 6 exhibits derivatives (2-4) are in agreement with the Lipinski rule of five indicating the good bioavailability of MGP derivatives, whereas derivatives (5-10) don't satisfy half of Lipinski's rules.
Derivatives with low molecular weight are easily absorbed, diffused and transported as compared to high molecular weight ones. As molecular weight increases except certain limit, the bulkiness of the molecules is also increasing. TPSA (Topological Polar Surface Area) is very important physicochemical property of molecule that provides the information about polarity of compounds. This parameter was screened for analyzing drug transport properties. Polar surface area is the sum of all polar atoms mainly oxygen and nitrogen including attached hydrogen. The drug-likeness score of lead derivatives is predicted with combination of GPCR, ion channel modulator, kinase inhibitor, nuclear receptor ligands, protease inhibitor, and enzyme inhibitor, which has been employed to identify the efficacy of molecules. When the bioactivity score is larger that means that the probability of the specific molecule to be active is higher. If bioactivity score of molecules is greater than 0.00, that molecule is qualified to be promising biological activities. When the scores are ranging from À0.50 to 0.00, they are considered to be moderately active. Finally, if bioactivity scores are less than À0.50 it is presumed to be inactive. The bioactivity scores of all designed MGP derivatives are displayed in Table 7. The obtained values of drug likeness scores exhibited that derivatives (2-4 and 8-10) followed the promising drug-likeness.
The bioactivity score gives the information about the binding cascade of the drugs that is applied for the development of a new functional drug with increased binding selectivity profile and less undesirable effects. All the selected MGP derivatives were evaluated to toxicity profile and given in Table 8 above.

Calculation of QSAR and pIC 50
The QSAR study is a computational modeling pathway for investigating correlations among the structural

Conclusions
In this investigation, the inherent characteristic of stability and biochemical behavior of MGP and its derivatives are studied for suggesting new antimicrobial compounds with improved activity. Important properties of biological chemistry and chemical reactivity such as antiviral activity, atomic partial charges and physicochemical properties were evaluated. They indicated that these derivatives are good candidate molecules for the development of drugs. All the designed MGP derivatives have improved physicochemical score than parent MGP. It is resulting that modified compounds are more reactive than parent drug. Insertion of various aliphatic and aromatic groups in MGP structure can significantly improve their biological activity. Antiviral prediction of MGP derivatives 2-10 revealed their promising potential activity against various pathogens. These observations were rationalized by molecular docking that pointed out promising antimicrobial effectiveness of MGP derivatives. Many of these analogs showed remarkable binding interactions and binding energy with DAOCS protein of S. clavuligerus. MGP derivatives (3-4, 8, and 10) have shown in silico potent ability to fight against S. clavuligerus. The selected complexes of derivatives 8 and 10 are also dynamically stable and interaction between the compound and receptor protein is dominated by van der Waals and electrostatic interactions. Finally, these derivatives were analyzed for their pharmacokinetic and drug-likeness properties which allowed understanding that the combination of toxicity prediction, bioactivity prediction, and drug-likeness had promising results. These promising results are obtained because most of designed molecules have improved kinetic parameters and all rules of drug-likeness are satisfied. In terms of biological activity, the results obtained were interesting result. Finally, this research may be useful to understand the chemical, thermal, physicochemical, biological and pharmacokinetic properties of MGP derivatives.

Funding
This work was supported by the "Ministry of Science and Technology (MOST)", Government of Bangladesh