International Journal of Pharmaceutical and Phytopharmacological Research
ISSN (Print): 2250-1029
ISSN (Online): 2249-6084
Publish with eIJPPR Submission
2020   Volume 10   Issue 4

In-Silico Evaluations of Some α -Methylene-γ-Butyrolactone Analogues

 

Feten Beji1,2, Reem al Harbi1, Arshi Naqvi1,3*

1 Department of Chemistry, Taibah University, 30002 Al-Madinah Al-Munawarah, Saudi Arabia.

2 Laboratory of Heteroatom Organic Chemistry, University of Carthage, Faculty of Sciences of Bizerte, 7021-Jarzouna, Tunisia.

3 BioDiscovery-Solutions for future, Plot No 29, FF2, 2nd street, Pearl Astragal Apartment, Perumbakkam, off Medavakkam, Solinganallur Main Road, Chennai (Tamil Nadu)-600100 (India).


ABSTRACT
The efficacy along with proper pharmacokinetics and toxicity profiling is the major determinants for the present scenario of successful drug development. The toxicity prospect accompanying weak ADMET (absorption, distribution, metabolism, elimination, and toxicity) profile are the vital justifications of late costly jeopardy of drug development. To predict the ADMET idiosyncrasy, in-silico inspection of some α-methylene-γ-butyrolactones (2-7) was induced on the foundation of several physico-chemical criteria to predict their pharmacokinetic, pharmacodynamic, drug-likeness, bioactivity and molecular docking profile exploiting several computational tactics. 


Key Words: lactones, α-methylene-γ-butyrolactones, in-silico molecular docking, ADMET, Drug likeness, bioactivity.


INTRODUCTION

Pharmaceutical consultation is an integral part of the retail sale of drugs to the public [1]. The liposomal formulations are targeted to deliver the important drug combinations to the body [2-4]. Lactones are compounds that are extensively dispersed in nature and several of them have been elucidated to date with influential bioactivity against a variety of human cancer cell lines, bacteria, and fungi. They exhibit interesting and useful biological activities. There is an antibacterial and antifungal property of lactones described in the literature. Much research, which

appears understandable, refers to antitumor and cytotoxic activity lactones. The factual “hunt” for the presence of lactose in plants began at the end of the 60s of the preceding century [5, 6]

Sesquiterpene lactones were estimated as one of the largest classes of natural compounds [7]. The α-methylene-γ-lactone ring is a prime structural skeleton in several natural products, majorly the sesquiterpene lactones (Figures 1a and 1b) [8, 9]. It was gauged in 1985, that almost 10% of the familiar 30,000 natural products comprises of this α-methylene γ-lactone functionality [10].

 

 

 

 

https://www.mdpi.com/molecules/molecules-21-00130/article_deploy/html/images/molecules-21-00130-g001.png

Figure 1a: General structure of α-methylene-γ-lactones

Figure 1b: Some sesquiterpene lactones.

 

 

This class of compounds has gain attraction due to their unique biological traits and in various cases, this high potency is reported to be due to the existence of the electrophilic exocyclic enoate moiety, which is responsible for trapping nucleophilic residues existing in the active site of target enzymes. They have been reported as DNA polymerase inhibitors, cellular steroidal inhibitors, nuclear vitamin D receptor inhibitors, blockers of tumor necrosis factor-α production, etc [11, 12]. The potency of these drug candidates is due to their cytotoxic, anti-inflammatory, antiallergenic, phytotoxic, anti-tumor, and antimicrobial properties [13-16].

It is computationally and economically unfeasible to develop and sieve candidates with anti-microbial activity from among innumerable compounds. The development of Computer-aided drug designing (CADD) or in silico notion is a promising shortcut to resolve the cost and time issues [17]. The computational approach legalizes the calculations of the copious number of quantitative descriptors on the foundation of molecular structural particulars and is very effective in optimizing important details such as biological activity or toxicity. Meanwhile, virtual screenings are useful in providing additional guidance for the design and development of new potent anti-microbial agents [18, 19].

To evaluate our previously reported α-methylene-g-butyrolactone derivatives [20] as anti-microbial agents virtually we exploited in-silico computational tools to virtually screen for their pharmacodynamic, pharmacokinetic, drug-likeness, bioactivity, and molecular docking traits.

 

EXPERIMENTAL

Ligand identification

The ligands to be screened in the present study i.e. α-methylene-g-butyrolactones analogs have been synthesized and characterized previously [20].

 

 

1

R= CH2, C2H4, C3H6, C4H8, C5H10, C6H12.

2-7

Figure 1: Structures of the selected compounds 1-7 for in-silico studies.

 

 

In-silico predictive studies

Effective and efficient integration of pharmacokinetics/pharmacodynamics (PK/PD) traits of

compounds in initial stages of development will not only assist with the compound selection but also guides the devising of systematic clinical development tactics [21]. SwissADME [22] and ProTox [23] were utilized to predict the PK/PD properties of the target ligands 1-7. For the prediction of the drug-like properties, drug-likeness is considered to be one of the qualitative ideas employed and it plays a pivotal role in demonstrating whether the studied drug examinees are alike the known drugs or not. The targeted ligands 1-7 were assessed for predicting their drug-likeness on the motive of five distinct filters namely Egan [24], Ghose [25], Muegge [26], Veber [27] and Lipinski [28] rules which were further escorted by bioavailability and drug-likeness scores using the Molsoft software (https://www.molsoft.com).

Also, these compounds were evaluated for in-silico bioactivity prediction using the software from Molinspiration Cheminformatics server, which screens the compounds against six important classes of drugs, i.e. G protein-coupled receptors ligand (GPCR ligands), ion channel blockers/modulators, kinase inhibitors, nuclear receptor ligands, protease inhibitors, and enzyme inhibitors.

In-silico molecular docking screenings were engaged to understand the interaction between the drug and the chosen receptor. While crafting the drug molecules to be effective antimicrobial agents, the supreme targets considered as a prime choice are those enzymes that are associated with the biosynthetic procedure of the microbe cell wall. To achieve this goal, the enzyme glucosamine-6-phosphate synthase (GlmS, GlcN-6-P synthase, L-glutamine: D-fructose-6P amido-transferase, EC 2.6.1.16) was contemplated for performing the docking studies [29]. The earlier mentioned enzyme is required both in the fungal and bacterial cell walls in the outset of the main building block viz. N-acetyl Glucosamine (the core amino sugar) [30, 31]. ACD/Labs-Chemsketch program was engaged in crafting 3D atomic coordinates of the ligands. The PDB ID, 2VF5 was considered for docking studies and it was obtained from the protein data bank (PDB) (Source: www.rcsb.org/pdb/). Utilization of the Dundee PRODRG2 server endeavored to do energy minimization [32]. Autodock4 program having a graphical user interface from “Auto-Dock Tools (ADT, 1.5.6)” was approached for the docking studies [33] and was implemented to recognize torsion angles in ligands, add the solvent model and allocate the Gasteiger charges to both the ligand and the protein. A 60x60x60 point grid was constructed representing X, Y, Z-axis. The values of RMSD (root mean square deviation) was utilized for cluster investigation. Furthermore, the cluster obtaining the lowest value for energy was appraised to be the most authentic solution. UCSF Chimera 1.11.2 and mcule, a web interface was employed for 3D visualization of the ligand-protein alliance.

 

RESULT AND DISCUSSION

The research work outlined here was escorted with the motive to foretell the physicochemical hallmarks, pharmacokinetic/ADME and toxicity traits, drug-likeness, and molecular docking screenings of our previously reported compounds 1-7 using in-silico computational tools.

Assorted physicochemical attributes like some rotatable bonds, count of specific atom types, lipophilicity, water solubility, and molecular refractivity were itemized. Topological Polar Surface Area (TPSA) is a fundamental physiochemical trait appraised for gauging drug transport characteristics. Calculations were attempted to forecast in-silico % absorption for all the aimed compounds by the reported formula (%ABS = 109-(0.345 X TPSA) [33]. All the targeted ligands revealed excellent in-silico % absorption with the highest being 99.93%. These physicochemical traits are given in Table 1.

 

 

Table 1: Physicochemical properties of the selected compounds 1-7.

Comp.

No.

Fraction Csp3a

No. of rotatable bonds

HBAb

HBDc

iLogPd

M.Re

Log Sf

TPSAg

In-silico % absorption

1

0.88

5

5

0

1.92

50.57

S

71.64

84.28

2

0.40

0

2

0

1.31

24.85

S

26.30

99.93

3

0.50

0

2

0

1.57

29.65

S

26.30

99.93

4

0.57

1

2

0

1.79

34.46

S

26.30

99.93

5

0.62

2

2

0

2.02

39.27

S

26.30

99.93

6

0.67

2

2

0

2.23

44.07

S

26.30

99.93

7

0.70

4

2

0

2.49

48.88

S

26.30

99.93

aThe ratio of sp3 hybridized carbons over the total carbon count of the molecule; bnumber of hydrogen bond acceptors; cnumber of hydrogen bond donors; d lipophilicity; eMolar refractivity, fWater solubility (SILICOS-IT [S=Soluble]); gtopological polar surface area (Å2).

 

 

Predictions of the pharmacokinetic/ADME and toxicity properties of the tested compounds 1-7 are given in Table 2. All the screened compounds demonstrated high gastrointestinal (GI) absorption and emerged as P-gp (p-glycoprotein) non-inhibitors. All the tested molecules were efficient to move through the blood-brain barrier (BBB) except compound 2. The forecast for the HIA (human gastrointestinal absorption), passive BBB permeation, and P-gp substrates is revealed jointly in the BOILED-Egg diagram as displayed in Figure 2. None of the tested ligands impede the Cytochrome P450 isomers and emerged to be low on skin permeability. Compound 1,2 and 7 are predicted to be non-toxic (class 6) in nature while the rest of the compounds fall in class 5, that it may be harmful if swallowed.

 

 

Table 2. Pharmacokinetic/ADME and toxicity properties of the selected compounds 1-7.

 

Comp.

No

Pharmacokinetic/ADME properties

Toxicity

GI

absa

BBB permeantb

P-gp substratec

CYP1A2 inhibitord

CYP2C19 inhibitore

CYP2C9 inhibitorf

CYP2D6 inhibitorg

CYP3A4 inhibitorh

Log Kpi

LD50 (mg/kg)

1

High

Yes

No

No

No

No

No

No

-7.45

25000

2

High

No

No

No

No

No

No

No

-6.37

11990

3

High

Yes

No

No

No

No

No

No

-6.29

3700

4

High

Yes

No

No

No

No

No

No

-6.07

5000

5

High

Yes

No

No

No

No

No

No

-5.76

5000

6

High

Yes

No

No

No

No

No

No

-5.67

5000

7

High

Yes

No

No

No

No

No

No

-5.17

8000

aGastro Intestinal absorption, bBlood Brain Barrier permeant, cP-glycoprotein substrate, d CYP1A2: Cytochrome P450 family 1 subfamily A member 2 (PDB:2HI4), e CYP2C19: Cytochrome P450 family 2 subfamily C member 19 (PDB:4GQS), fCYP2C9: Cytochrome P450 family 2 subfamily C member 9 (PDB:1OG2), g CYP2D6: Cytochrome P450 family 2 subfamily D member 6 (PDB:5TFT), h CYP3A4: Cytochrome P450 family 3 subfamilies A member 4 (PDB:4K9T), iSkin permeation in cm/s.

 

Figure 2: BOILED-Egg diagram of the selected compounds 1-7.

 

 

Drug likeness is explored as an essential element that furnishes the base for the candidates to be an influential drug aspirant. Several rules namely Lipinski, Ghose, Veber, Egan, and Muegge were contemplated to predict drug-likeness of the candidate molecules (1-7) to realize whether they can be potential drug candidates as per some acute benchmarks like molecular weight, LogP, number of hydrogen bond donors and acceptors. The number of breaches to the above-revealed rules jointly with bioavailability and drug-likeness scores is revealed in Table 3. None of the compounds 1-7 violated Lipinski, Veber, and Egan rule. All the compounds except 1 and 7 breached Ghose rule with 1-3 violations. All the compounds except for compound 1, discarded to be drug-like with one violation as per Muegge rule. All the tested compounds exhibited bioavailability scores around 0.55 and flaunted moderate drug-likeness scores ranged from -1.67 to -0.84. The bioavailability radar of the compounds 1-7 is showcased in figure 3. Forecast of bioactivity scores of all the tested ligands against six different protein structures namely GPCR ligand, ion channel modulator, a kinase inhibitor, nuclear receptor ligand, protease inhibitor, and enzyme inhibitor is given in Table 4. All these integers stipulate binding affinity of the targeted ligands (1-7) to the declared receptors and enzymes. Most of the tested compounds have displayed negative bioactivity while compounds 1 and 7 displayed affinities towards some of the selected protein structures. The prime in-silico conformation was endorsed by employing molecular docking studies of the compounds 1-7 on GlcN-6-P synthase (PDB ID: 2VF5) as the target protein receptor. To achieve this objective AutoDock 4.0 was appointed. The docking pose of the ligand 6 is manifested in Figure 4 and their docking scores are shown in Table 3. Docking scores range from -3.6 to -5.1.

 

 

Table 3. Drug likeness predictions and docking scores of the selected compounds 1-7.

Comp.

No.

Lipinski

violations

Ghose

violations

Veber

violations

Egan

violations

Muegge

violations

Bioavailability

Score

Drug Likeness

score

Docking

score

1

0

0

0

0

0

0.55

-0.88

-5.1

2

0

3

0

0

1

0.55

-1.67

-3.6

3

0

3

0

0

1

0.55

-1.35

-4.0

4

0

3

0

0

1

0.55

-1.11

-4.1

5

0

2

0

0

1

0.55

-1.08

-4.3

6

0

1

0

0

1

0.55

-0.84

-4.7

7

0

0

0

0

1

0.55

-0.96

-4.6

 

Figure 3: Bioavailability Radar of the tested compounds (1-7) The pink area represents the optimal range for each property (lipophilicity: XLOGP3 between − 0.7 and + 5.0, size: MW between 150 and 500 g/mol, polarity: TPSA between 20 and 130 (Å2)., solubility: log S not higher than 6, saturation: the fraction of carbons in the sp3 hybridization not less than 0.25, and flexibility: no more than 9 rotatable bonds. In this example, the compound is predicted not orally bioavailable, because too flexible and too polar.

 

Figure 4: Docking of compound 6 into the active site of GlcN-6-P synthase (PDB ID: 2VF5).

Table 4: Bioactivity of the selected compounds 1-7.

Comp.

No.

GPCR Ligand

Ion channel modulator

Kinase inhibitor

Nuclear receptor ligand

Protease inhibitor

Enzyme inhibitor

1

-0.55

0.24

-0.70

-0.71

-0.18

0.31

2

-3.22

-3.34

-3.51

-2.24

-3.43

-2.92

3

-2.94

-2.62

-2.96

-2.07

-3.26

-2.13

4

-1.89

-1.46

-2.00

-1.38

-2.24

-1.16

5

-0.76

-0.46

-0.79

-0.32

-1.10

-0.09

6

-0.73

-0.38

-0.76

-0.27

-0.86

-0.11

7

-0.53

-0.28

-0.54

-0.12

-0.80

0.06

 

 

CONCLUSION

The research work consolidated here is related to the in-silico forecasting of the physicochemical and ADMET features, drug-likeness, and molecular docking of some previously reported α-methylene-γ-lactone derivatives utilizing online available servers.  The anticipated study revealed that all the screened compounds 1-7 displayed excellent in-silico % absorption where the highest obtained value was 99.93%. The screened compounds have showcased high GI absorption and are forecasted to be CNS active candidates as they can easily cross the blood-brain barrier (BBB) except compound 2 which is considered to be CNS inactive. All the tested compounds emerged to be non-inhibitors of P-gp and Cytochrome P450 isomers and possess low skin permeability. Compounds 1,2 and 7 are non-toxic as their LD50 is more than 5000 mg/kg while the rest of the compounds emerged to, may be harmful if swallowed. All the screened compounds were like a drug according to Lipinski, Veber, and Egan rule while except compound 1 and 7, all the candidates breached Ghose rule. Also, most of the compounds except compound 1, rejected to be drug-like with 1 violation according to the Muegge rule. Their bioavailability score is 0.55 and their drug-likeness and docking scores range from -1.67 to -0.84 and -3.6 to -5.1 respectively. The compound 6 has emerged to be the best compound out of the α-methylene-γ-lactone derivatives in the in-silico predictions with a drug-likeness score of -0.84 and docking score of -4.7. Most of the tested compounds have displayed low bioactivity scores while compound 1 and 7 flaunted affinity towards some of the selected protein structures. In conclusion, these forecasted evaluations give the information about the physico-chemical as well as ADMET attributes along with illuminating the drug-likeness, bioactivity, and binding/affinity fashion on the target protein of the tested ligands thus aiding the lead for the drugs of the future with additional efficacy.

REFERENCES

 

  1. Sergeevna, S. M., Efimovna, L. E., Vladimirovna, A. I. Improvement of pharmaceutical consultation process in drugstores. Journal of Advanced Pharmacy Education & Research| Jan-Mar, 2020; 10(1), 136-142.
  2. MP, S. C., Prasobh, G. R., Jaghatha, T., Aswathy, B. S., Remya, S. B. An Overview on Liposomal Drug Delivery System. International Journal of Pharmaceutical and Phytopharmacological Research, 2019; 9(2), 61-68.
  3. Mozayyeni, N., Morsali, A., Bozorgmehr, M. R., Beyramabadi, S. A. Structural and mechanistic studies of γ-Fe2O3 nanoparticle as troxacitabine drug nanocarrier. Archives of Pharmacy Practice, 2019; 10(1), 31-37.
  4. Alzahrani, S., Alosaimi, M. E., Alamri, A. A., Alotaibi, M., Almatar, E. A., Almanea, B. A., ... Al-Shams, A. A. Association between Knowledge and Drug Adherence in Patients with Hypertension in Saudi Arabia. Archives of Pharmacy Practice, 2019; 10(3), 71-6.
  5. Kupchan S.M., Kelsey J.E., Maruyama M., Cassady J.M., Hemıngway J.C., Knox J.R. Tumor inhibitors. XLI. Structural elucidation of tumor-inhibitory sesquiterpene lactones from Eupatorium rotundifolium. J. Org. Chem., 1969; 34(12): 3876–3883.
  6. Hanson R.L., Lardy H.A., Kupchan S.M. Inhibition of phosphofructokinase by quinose methide and α -methylene lactone tumor inhibitors. Science, 1970; 168(3929): 378–380.
  7. Franck B. Key building blocks of natural product biosynthesis and their significance in chemistry and medicine. Angew. Chem. Int. Ed. Engl., 1979; 18(6): 429–439.
  8. Gach  K,  Dlugosz  A,  Janecka  A.  The role of oxidative stress in anticancer activity of sesquiterpene lactones.   Naunyn-Schmiedeberg's archives of pharmacology. 2015; 388: 477-86.
  9. Wu, Y.; Wang, D.; Gao, Y.; Feng, J.; Zhang, X. New α-Methylene-γ-Butyrolactone Derivatives as Potential Fungicidal Agents: Design, Synthesis, and Antifungal Activities. Molecules 2016, 21, 130.
  10. Hoffmann, H. M. R.; Rabe, J. Angew. Chem. Int. Ed. Engl. 1985, 24, 94–110.
  11. Konaklieva, M. I.; Plotkin, B. J. Mini Rev. Med. Chem. 2005, 5, 73–95.
  12. Tim G. Elford, Dennis G. Hall, Chapter 4 - Total Synthesis of Chinensiolide B: An Exercise in Functional Group Selectivity with the Help of a Fortuitous Side Product and a Good Nose. Strategies and Tactics in Organic Synthesis, 10, 2014, 79-112.
  13. Janecki, T.; Blaszczyk, E.; Studzian, K.; Janecka, A.; Krajewska, U.; Rozalski, M. J. Med. Chem. 2005, 48, 3516–3521.
  14. Duraıpandıyan V., Al-Harbı N.A., Ignacımuthu S., MuthukumaR C. Antimicrobial activity of sesquiterpene lactones isolated from traditional medicinal plant, Costus speciosus (Koen ex. Retz.) Sm. BMC Complement Altern. Med., 2012; 12(13): 12–23.
  15. Wu Y, Wang D, Gao Y, Feng J, Zhang X. New α-methylene-γ-butyrolactone derivatives as potential fungicidal agents: design, synthesis and antifungal activities. Molecules. 2016 Feb;21(2):130.
  16. Janecka A., Wyrebska A., Gach K., Fichna J., Janecki T., Natural and synthetic α-methylenelactones and α-methylenelactams with anticancer potential. Drug Discov. Today, 2012; 17, 561–572.
  17. Gao, Y.; Li, J.; Song, Z.; Song, J.; Shang, S.; Xiao, G.; Wan, Z.; Rao, X. Turning renewable resources into value-added products: Development of rosin-based insecticide candidates. Ind. Crop. Prod. 2015, 76, 660–671.
  18. García, I.; Fall, Y.; Gómez, G.; González-Díaz, H. First computational chemistry multi-target model for anti-Alzheimer, anti-parasitic, anti-fungi, and anti-bacterial activity of GSK-3 inhibitors in vitro, in vivo, and in different cellular lines. Mol. Divers. 2011, 15, 561–567
  19. Prado-Prado, F.J.; Borges, F.; Perez-Montoto, L.G.; González-Díaz, H. Multi-target spectral moment: QSAR for antifungal drugs vs. different fungi species. Eur. J. Med. Chem. 2009, 44, 4051–4056.
  20. Feten Beji, Rafaa Besbes, Hassen Amri. Synthesis of α-alkylidene-ß-ethoxycarbonyl Cyclopentanones and - γ -butyrolactones. Synthetic Comm. 30, 21 (2000) 3947-3954.
  21. Miller R, Ewy W, Corrigan BW, Ouellet D, Hermann D, Kowalski KG, Lockwood P, Koup JR, Donevan S, El-Kattan A, Li CS. How modeling and simulation have enhanced decision making in new drug development. Journal of pharmacokinetics and pharmacodynamics. 2005 Apr 1;32(2):185-97.
  22. 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, 42717.
  23. Banerjee P, Dehnbostel F.O, Preissner R. Prediction is a Balancing Act: Importance of Sampling Methods to Balance Sensitivity and Specificity of Predictive Models based on Imbalanced Chemical Data Sets. Frontiers in Chemistry. 2018; 6, 362.
  24. Egan, W.J., Merz, K.M., Baldwin, J.J. Prediction of Drug Absorption Using Multivariate Statistics. J. Med. Chem. 2000; 43:3867-3877.
  25. Ghose, A.K., Viswanadhan, V.N., Wendoloski, J.J. A knowledge-based approach in designing combinatorial or medicinal chemistry libraries for drug discovery. 1. A qualitative and quantitative characterization of known drug databases. J Comb. Chem. 1999; 1:55-68.
  26. Muegge, I., Heald, S.L., Brittelli, D. Simple selection criteria for drug-like chemical matter. J. Med. Chem. 2001; 44:1841-1846.
  27. Veber, D.F., Johnson, S.R., Cheng, H.Y., Smith, B.R., Wars, K.W., Kopple, K.D. Molecular properties that influence the oral bioavailability of drug candidates. J. Med. Chem. 2002; 45:2615-2623.
  28. Lipinski, C.A., Lombardo, F., Dominy, B.W., Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug. Deliv. Rev. 2001; 46:3-26.
  29. Chmara, H.; Andruszkiewicz, R.; Borowski, E. Inactivation of glucosamine-6-phosphate synthetase from Salmonella typhimurium LT2 SL 1027 by N-beta-fumarylcarboxyamido-l-2,3-diamino-propionic acid. Biochem. Biophys. Res. Commun. 1984, 120, 865-872.
  30. Marshall, N.J.; Andruszkiewicz, R.; Gupta, S.; Milewski, S.; Payne. J.W. Structure-activity relationships for a series of peptidomimetic antimicrobial prodrugs containing glutamine analogues. J. Antimicrob. Chemother. 2003, 51, 821-831.
  31. Lather, A.; Sharma, S.; Khatkar, A. Virtual Screening of Novel Glucosamine-6-Phosphate Synthase Inhibitors. Comb. Chem. High Throughput Screen. 2018; 21,182-193.
  32. Schüttelkopf, A.W.; Aalten, D.M. PRODRG: A tool for high-throughput crystallography of protein-ligand complexes. Acta. Cryst. 2004; 60, 1355–1363.
  33. Morris, G.M.; Goodsell, D.S.; Halliday, R.S.; Huey, R.; Hart, W.E.; Belew, K.R.; Olson, A.J. Automated docking using a Lamarckian genetic algorithm and empirical binding free energy function. J. Comput. Chem. 1998; 19, 1639-1662.
References


 

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

International Journal of Pharmaceutical and Phytopharmacological Research
© 2024 All rights reserved