An approach to design, In silico predictions and molecular docking studies of 1,3,4-oxadiazolyl sulphonamides as possible inhibitors of bacterial Mur enzymes, the amino acid ligases in peptidoglycan synthesis

B. HASEENA BANU

Easy Links

Research article Open Access

An approach to design, In silico predictions and molecular docking studies of 1,3,4-oxadiazolyl sulphonamides as possible inhibitors of bacterial Mur enzymes, the amino acid ligases in peptidoglycan synthesis

B. HASEENA BANU*

Krishna Teja Pharmacy College, Affiliated to JNTUA, Chadalawada Nagar, Renigunta Road, Tirupati, AP, Indi

*Corresponding author: B. HASEENA BANU, Krishna Teja Pharmacy College, Affiliated to JNTUA, Chadalawada Nagar, Renigunta Road, Tirupati, AP, India, Tel: 9949640769; E-mail:

Received: 13 November, 2017; Accepted: 07 December, 2017; Published: 11 December, 2017

Citation: B. Haseena Banu (2017) An approach to design, In silico predictions and molecular docking studies of 1,3,4-oxadiazolyl sulphonamides as possible inhibitors of bacterial Mur enzymes, the amino acid ligases in peptidoglycan synthesis. J Adv Res Biotech 2(3): 1-12. DOI: http://dx.doi.org/10.15226/2475-4714/2/3/00126

Abstract

The extensive use of antibiotics in hospitals and community since their introduction into medical practice has created major evolutionary pressures in bacteria to develop various resistance mechanisms. This phenomenon has led to increased morbidity, mortality and health care costs. The search for new antibacterial agents directed towards novel targets became highly imperative. The biosynthetic pathway of cytoplasmic peptidoglycan precursor is currently gaining much interest as a target site for antibacterial therapy. Since the mid-1990s, many inhibitors of the Mur cytoplasmic enzymes have been reported, but none has yet led to the development of a clinically utilized therapeutic agent. In view of the potentiality of sulphonamide and 1,3,4-oxadiazole moieties, it was planned to incorporate both the scaffolds. Thus a series of 4-Amino- N-[(5- phenyl-1,3,4- oxadiazol-2-yl)methyl] substituted benzene-1- sulfonamides (1-10) were designed, predicted for molecular properties and docking studies were conducted on Mur enzymes. The presence of the two pharmacophoric scaffolds, oxadiazole and sulphonamide influenced the good bioactive scores. Docking studies further supported the biological activities. The results indicated the possible inhibitory activity towards the Mur enzymes. The results also indicate the selectivity towards Mur A, Mur C, & Mur E enzymes. Of all the compounds, compound 6 and compound 9 with isopropyl substitution exhibited the highest scores against all the three Mur enzymes. Compound 7 showed good binding interactions with Mur C and Mur E.

Abbreviations

Alr: Alanine Racemase

Ddl: D-Ala-D-Ala Ligases

GPCR: G-Protein-Coupled Receptor

IUPAC: The International Union of Pure And Applied Chemistry

Mur A Udp-N-Acetylglucosamine Enolpyruvyltransferase

Mur B Udp-N-Acetylenolpyruvylglucosamine Reductase

Mur C- N-Acetylmuramate L-Alanine Ligase

Mur D- N-Acetylmuramoyl-L-Alanine D-Glutamate Ligase

Mur E - N-Acetylmuramoyl-L-Alanyl-D-Glutamate Meso-Diaminopimelate Ligase

PASS: Prediction Of Activity Spectra For Biologically Active Substances.

PDB ID-Protein Data Bank-Identification

QSAR: Quantitative Structure Activity Relationship

Rtob: Rotatable Bonds

SAR: Structure Activity Relationship

SD File: Structure Data File

SMILES: Simplified Molecular-Input Line-Entry System

TPSA: Topological Polar Surface Area

UDP : Uridine Diphosphate

Introduction

The treatment of microbial infections still remains an important and challenging medical problem. Sulfonamides (sulfa drugs) were the first drugs largely employed and systematically used as preventive and chemotherapeutic agents against various diseases [1]. Over 30 drugs containing this functionality are in clinical use, including antihypertensive agent bosentan, antibacterial, antiprotozoal, antifungal, antiinflammatory, lipoxygenase inhibitor, nonpeptidic vasopressin receptor antagonists and translation initiation inhibitors [2-8].

The literature is enriched with excellent antimicrobial activity for compounds containing the 2,5-substituted 1,3,4-oxadiazole core [9-11]. Oliveira and co-workers reported synthesis and anti-staphylococcal activity of 1, 3, 4-oxadiazolines against strains of S.aureus, resistant to methicillin and amino glycosides and that encode efflux proteins (multidrug drugs resistant—MDR more active than the standard drug chloramphenicol [12]. The antibacterial and antifungal activity of 2-(5-amino-1,3,4-oxadiazol-2-yl)-4-bromophenol and 5-(3,5-dibromophenyl)- 1,3,4-oxadiazol-2-amine were investigated against gram-positive bacteria, gram-negative bacteria and fungal species approximately equal to the standard drugs of treatment streptomycin and Griseofulvin [13]. 2,5-disubstituted 1,3,4 oxadiazoles were reported as potent antibacterial and antifungal agents [14](Figure 1).

Despite of the development of several new antimicrobial agents, their use has been clinically limited due to bacterial resistance and pharmacokinetic insufficiencies [15]. The aim of the research was to target the biosynthetic pathway of cytoplasmic peptidoglycan, which is currently gaining much interest as a target site for antibacterial therapy and provides good prospects for discovering selective bacterial inhibitors [16]. Therefore, it prompted us to design the new molecules in such a way that the basic pharmacophoric scaffolds were conjugated with each other, sulfonamide moiety at 2^nd position of 1,3,4-oxadiazol ring and phenyl ring at position 5, altered by varying substituents on phenyl ring. These compounds were predicted molecular properties and drug likeness scores, compared the effect of different substitutions in the pharmacophore on the drug likeness properties. Molecular docking studies of the compounds on Mur enzymes were conducted for ensuring the possibility of mechanism of action.

1, 3, 4-oxadiazolyl sulphonamides
Figure 1: Design of new compounds by conjugation of sulfonamide and 1,3, 4-oxadiazole moieties

Methodology

The structural analogue based drug design has been performed using MARVIN SKETCH, in which molecules in-silico drug likeness and molecular property prediction tool are important the new molecules, designed on the basis of SAR and pharmacophore study, were then imported into JME molecular editor to allow different properties to be calculated.

Preparation of Mol and Mol2 files

The chemical structures drawn were saved in the Mol and Mol2 formats using Chemaxon software (www.chemaxon.com).

Prediction of molecular properties from chemical structure

Chemicalize is a free chemical structure miner and web search engine developed and owned by ChemAxon. The main purpose of chemicalize.org is to identify chemical names (SMILES, traditional and IUPAC names) on websites and convert them to chemical structures. chemicalize.org provides other services such as structure based predictions, chemical search, and a “chemicalized” web search (www.chemicalize.org).

Molinspiration offers broad range of cheminformatics software tools supporting molecule manipulation and processing, including SMILES and SDfile conversion, normalization of molecules, generation of tautomers, molecule fragmentation, calculation of various molecular properties needed in QSAR, molecular modeling and drug design, high quality molecule depiction, molecular database tools supporting substructure and similarity searches. Our products support also fragment-based virtual screening, bioactivity prediction and data visualization. Molinspiration tools are written in Java, therefore can be used practically on any computer platform (www.molinspiration.com).

Molsoft offers software tools and services in lead discovery, modeling, chem-informatics, bioinformatics, and corporate data management; and forms partnerships with biotechnology and pharmaceutical companies. For many years Molsoft is building unique technologies for structure prediction and is increasing our understanding of spatial organization of biological molecules and their interactions with each other, their biological substrates and drug-like molecules at atomic level (www.molsoft.com).

Selection of targets

The synthesis of the pentapeptide chain is facilitated by enzymes of the Mur pathway and involve MurA-F the D-Ala-D-Ala ligases (Ddl) and alanine racemase (Alr). MurA and –B convert UDP-N-acetylglucosamine into UDP-N-acetylmuramic acid, the nucleotide sugar used by the ligases MurC-F to produce the pentapeptide chain which is later transferred across the cytoplasmic membrane to be incorporated in the growing cell wall [16]. In contrast to well established approaches exploiting the later stages of prostaglandin synthesis (the discovery and development of β -lactam and glycopeptide antibiotics) structure-based approaches were used for the design and synthesis of selective inhibitors of early stage Mur enzymes(Figure 2).

Molecular Docking Studies

Docking calculations have been applied in pharmaceutical research for nearly two decades. Computational approaches that ‘dock’ small molecules into the structures of macromolecular targets and ‘score’ their potential complementarily to binding sites are widely used in hit identification and lead optimization. Indeed, there are now a number of drugs whose development was heavily influenced by / or based on structure-based design and screening strategies [18].

The docking process involves the prediction of ligand conformation and orientation (or posing) within a targeted binding site. In general, there are two aims of docking studies: accurate structural modeling and correct prediction of activity by the identification of molecular features responsible for specific biological recognition, or the prediction of compound modifications that improve potency [19].

Figure 2: The cytoplasmic steps of peptidoglycan biosynthesis catalyzed by Mur enzymes

In order to rationalize the different biological results shown and bioactivity score predictions and PASS predicted values obtained for active compounds in each series, molecular docking studies were performed using SWISS dock software.

Results and Discussion

Design and nomenclature

Prompted by the medicinal importance of 1,3,4-oxadiazolyl and sulfonamide moieties, the structural analogues of oxadiazolyl sufonamides were designed in such a way that, the analogues have both biological active moieties.

The structural analogues have been designed in such a way that it will show more drug likeness score than the prototype molecule but having the same pharmacophore essential for the activity. A series of ten 4-amino-N-[(5-phenyl-1,3,4-oxadiazol-2-yl)methyl] substituted benzene-1-sulfonamides were designed and their chemical structures were drawn using Marvin sketch.

Different substituents such as electron donating (chloro, isopropyl, tert-butyl), electron releasing (hydroxyl, methoxy, nitro) etc were substituted on the benzene ring to study their effect on the predicted properties and also biological activities.

The list of designed analogues with drug likeness scores, and other molecular properties have been listed below. The position of substituent on benzene ring was modified in these designed molecules to get increased drug likeness score. Table 1 gives the information about different substituted compounds and table 2 indicates the nomenclature generated using soft ware chemicalize.org (Table 1, Table 2).

Table 1: Design of 1, 3, 4-oxadiazolyl sulphonamides

Compound	R
1.	H
2.	4-OH
3.	4-Cl
4.	4-NH₂
5.	4-OH,3-CH₃
6.	4-CH(CH₃)₂
7.	4-C(CH₃)₃
8.	4-N(CH₃)₂
9.	4-NO₂
10.	3,4,5-(OCH₃)₃

Table 2: Nomenclature of the designed 1,3,4-oxadiazolyl sulphonamides using Chemicalize.org

S.No	R	Nomenclature
1	H	IUPAC: 4-Amino-N-[(5-Phenyl-1,3,4-Oxadiazol-2-L)Methyl]Benzene-1-Sulfonamide SMILES: INCHI: 1S/C15H14N4O3S/C16-12-6-8-13(9-7-12)23(20,21)17-10-14-18-19-15(22-14)11-4-2-1-3-5-11/H1-9,17H,10,16H2 INCHI KEY: WOJFXQNTVVDWNV-UHFFFAOYSA-N
2	4-OH	IUPAC: 4-Amino-N-{[5-(4-Hydroxyphenyl)-1,3,4-Oxadiazol-2-Yl]Methyl}Benzene-1-Sulfonamide SMILES: NC1=CC=C(C=C1)S(=O)(=O)NCC1=NN=C(O1)C1=CC=C(O)C=C1 INCHI: 1S/C15H14N4O4S/C16-11-3-7-13(8-4-11)24(21,22)17-9-14-18-19-15(23-14)10-1-5-12(20)6-2-10/H1-8,17,20H,9,16H2 INCHI KEY: KVWVQFHGSYNLSC-UHFFFAOYSA-N IUPAC: 4-Amino-N-{[5-(4-Chlorophenyl)-1,3,4-Oxadiazol-2-Yl]Methyl}Benzene-1-Sulfonamide
3	4-Cl	SMILES: NC1=CC=C(C=C1)S(=O)(=O)NCC1=NN=C(O1)C1=CC=C(CL)C=C1 INCHI: 1S/C15H13CLN4O3S/C16-11-3-1-10(2-4-11)15-20-19-14(23-15)9-18-24(21,22)13-7-5-12(17)6-8-13/H1-8,18H,9,17H2 INCHI KEY: PPQSXSJNGVRRFB-UHFFFAOYSA-N
4	4-NH₂	IUPAC: 4-Amino-N-{[5-(4-Aminophenyl)-1,3,4-Oxadiazol-2-Yl]Methyl}Benzene-1-Sulfonamide SMILES: NC1=CC=C(C=C1)C1=NN=C(CNS(=O)(=O)C2=CC=C(N)C=C2)O1 INCHI: 1S/C15H15N5O3S/C16-11-3-1-10(2-4-11)15-20-19-14(23-15)9-18-24(21,22)13-7-5-12(17)6-8-13/H1-8,18H,9,16-17H2 INCHI KEY: CFKKLKAYSDHLBZ-UHFFFAOYSA-N
5	4-OH, 3-OCH₃	IUPAC: 4-Amino-N-{[5-(4-Hydroxy-3-Methoxyphenyl)-1,3,4-Oxadiazol-2-Yl]Methyl}Benzene-1-Sulfonamide SMILES: COC1=C(O)C=CC(=C1)C1=NN=C(CNS(=O)(=O)C2=CC=C(N)C=C2)O1 INCHI: 1S/C16H16N4O5S/C1-24-14-8-10(2-7-13(14)21)16-20-19-15(25-16)9-18-26(22,23)12-5-3-11(17)4-6-12/H2-8,18,21H,9,17H2,1H3 INCHI KEY: DQYQKXMOKGSEGM-UHFFFAOYSA-N
6	4-CH(CH₃)₂	IUPAC: 4-Amino-N-({5-[4-(Propan-2-Yl)Phenyl]-1,3,4-Oxadiazol-2-Yl}Methyl)Benzene-1-Sulfonamide SMILES: CC(C)C1=CC=C(C=C1)C1=NN=C(CNS(=O)(=O)C2=CC=C(N)C=C2)O1 INCHI: 1S/C18H20N4O3S/C1-12(2)13-3-5-14(6-4-13)18-22-21-17(25-18)11-20-26(23,24)16-9-7-15(19)8-10-16/H3-10,12,20H,11,19H2,1-2H3 INCHI KEY: XALHGEIOZWCCHJ-UHFFFAOYSA-N
7	4-C(CH₃)₃	IUPAC: 4-Amino-N-{[5-(4-Tert-Butylphenyl)-1,3,4-Oxadiazol-2-Yl]Methyl}Benzene-1-Sulfonamide SMILES: CC(C)(C)C1=CC=C(C=C1)C1=NN=C(CNS(=O)(=O)C2=CC=C(N)C=C2)O1 INCHI: 1S/C19H22N4O3S/C1-19(2,3)14-6-4-13(5-7-14)18-23-22-17(26-18)12-21-27(24,25)16-10-8-15(20)9-11-16/H4-11,21H,12,20H2,1-3H3 INCHI KEY: QVTZXCYLTOXWHY-UHFFFAOYSA-N
8	4-N(CH₃)₂	IUPAC: 4-Amino-N-({5-[4-(Dimethylamino)Phenyl]-1,3,4-Oxadiazol-2-Yl}Methyl)Benzene-1-Sulfonamide SMILES: CN(C)C1=CC=C(C=C1)C1=NN=C(CNS(=O)(=O)C2=CC=C(N)C=C2)O1 INCHI: 1S/C17H19N5O3S/C1-22(2)14-7-3-12(4-8-14)17-21-20-16(25-17)11-19-26(23,24)15-9-5-13(18)6-10-15/H3-10,19H,11,18H2,1-2H3 INCHI KEY: NKCXFZSYQJLWAT-UHFFFAOYSA-N
9	4-NO₂	IUPAC: 4-Amino-N-{[5-(4-Nitrophenyl)-1,3,4-Oxadiazol-2-Yl]Methyl}Benzene-1-Sulfonamide SMILES:NC1=CC=C(C=C1)S(=O)(=O)NCC1=NN=C(O1)C1=CC=C(C=C1)N(=O)=O INCHI: 1S/C15H13N5O5S/C16-11-3-7-13(8-4-11)26(23,24)17-9-14-18-19-15(25-14)10-1-5-12(6-2-10)20(21)22/H1-8,17H,9,16H2 INCHI KEY: CHTKFGIKJUONOV-UHFFFAOYSA-N
10	3,4,5-(OCH₃)₃	IUPAC: 4-Amino-N-{[5-(3,4,5-Trimethoxyphenyl)-1,3,4-Oxadiazol-2-Yl]Methyl}Benzene-1-Sulfonamide SMILES: COC1=CC(=CC(OC)=C1OC)C1=NN=C(CNS(=O)(=O)C2=CC=C(N)C=C2)O1 INCHI: 1S/C18H20N4O6S/C1-25-14-8-11(9-15(26-2)17(14)27-3)18-22-21-16(28-18)10-20-29(23,24)13-6-4-12(19)5-7-13/H4-9,20H,10,19H2,1-3H3 INCHI KEY: AZCWENOKCBERLG-UHFFFAOYSA-N

Molecular properties predictions

Lipinski’s Rule of Five describes molecular properties important for a drug’s pharmacokinetics in the human body, including their absorption, distribution, metabolism, and excretion [20]. Using Chemicalize, molinspiration and molsoft, molecular properties were predicted for the designed compounds. The results obtained revealed that all the compounds tested obeyed Lipinski rule of five with not more than 5 hydrogen bond donors (OH and NH groups), not more than 10 hydrogen bond acceptors (notably N and O), not more than 15 rotatable bonds (rotb), molecular weight under 500 g/mol and a partition coefficient, log P less than 5 except for compound 3&8 (Tables 3-6). All the compounds tested also passed the other Lipinski like filter such as bioavailability.

Absorption is defined as the process involved in getting a drug from its dosage form into the body and the ability to predict the percent oral absorption is primary goal in the design, optimization, and selection of potential candidates in the development of oral drugs. Topological Polar surface area (TPSA) is another key property that has been linked to bioavailability and it was found that passively absorbed molecules with a TPSA more than 140 are thought to have low oral availability [21]. TPSA obtained for the tested compounds were below 140 indicating their good oral bioavailability. Thus these results predicted the good drug likeliness, solubility, permeability and oral bioavailability of compounds. No compound violated the Lipinski rule of five (Table 3, Table 4, Table 5, Table 6).

Table 3: Molecular properties prediction of title compounds using Chemicalize.org

S.NO	R	Molecular Formula	Molecular weight	Log P	TPSA	SASA
1	H	C₁₅H₁₄N₄O₃S	330.362	0.78	111.11	436.55
2	4-OH	C₁₅H₁₄N₄O₄S	346.361	0.47	131.34	447.81
3	4-Cl	C₁₅H₁₃ClN₄O₃S	364.807	1.38	111.11	452.77
4	4-NH₂	C₁₅H₁₅N₅O₃S	345.376	-0.05	137.13	451.37
5	4-OH,3-OCH₃	C₁₆H₁₆N₄O₅S	76.387	0.32	140.57	495.57
6	4-CH(CH₃)₂	C₂₅H₂₃N₃O₅	372.441	2.02	111.11	528.98
7	4-C(CH₃)₃	C₂₆H₂₅N₃O₅	386.468	2.32	111.11	564.79
8	4-N(CH₃)₂	C₁₆H₁₆N₄O₅S	376.387	0.32	140.57	495.57
9	4-NO₂	C₁₅H₁₃M.N₅O₅S	375.359	0.72	156.93	476.76
10	3,4,5-(OCH₃)₃	C₁₈H₂₀N₄O₆S	420.44	0.3	138.8	580.24

Log P- Partition co efficient
TPSA - Topological Polar surface area
SASA- Solvent Accessible Surface Area

Table 4: Molecular properties prediction of title compounds using Chemicalize.org

S.NO	R	Lipinski rule of five	Bio- availability	Ghose filter	Lead likeness	Mugge Filter	Veber filter
1	H	Yes	Yes	Yes	Yes	Yes	Yes
2	4-OH	Yes	Yes	Yes	No	Yes	Yes
3	4-Cl	Yes	Yes	Yes	Yes	Yes	Yes
4	4-NH₂	Yes	Yes	Yes	Yes	Yes	Yes
5	4-OH,3-OCH₃	Yes	Yes	Yes	Yes	Yes	No
6	4-CH(CH₃)₂	Yes	Yes	Yes	Yes	Yes	Yes
7	4-C(CH₃)₃	Yes	Yes	No	No	Yes	Yes
8	4-N(CH₃)₂	Yes	Yes	Yes	Yes	Yes	Yes
9	4-NO₂	Yes	Yes	Yes	Yes	No	No
10	3,4,5-(OCH₃)₃	Yes	Yes	Yes	Yes	Yes	Yes

Table 5: Molecular properties prediction of title compounds using molinspiration.com

S.No	R	nON	nOHNH	N violations	nrotob	Volume	Mi Logp	n Atoms
1	H	7	3	0	5	272.04	0.903	23
2	4-OH	8	4	0	5	280.058	0.424	24
3	4-Cl	7	3	0	5	285.576	1.581	24
4	4-NH₂	8	5	0	5	283.329	0.021	24
5	4-OH,3-OCH₃	9	4	0	6	305.604	0.243	26
6	4-CH(CH₃)₂	7	3	0	6	321.99	2.416	26
7	4-C(CH₃)₃	7	3	0	6	338.227	2.61	27
8	4-N(CH₃)₂	8	3	0	6	317.946	1.006	26
9	4-NO₂	10	3	0	6	295.374	0.862	26
10	3,4,5-(OCH₃)₃	10	3	0	8	348.677	0.535	29

Table 6: Molecular properties prediction of title compounds using molsoft.com

S.No	R	No of HBA	No of HBD	Mol Log p	Mol Log s	Mol PSA (A2)	Mol volume (A3)	No of stereo centres
1	H	6	3	1.7	-3.71(in Log(moles/L) 64.84 (in mg/L)	92.98 A2	274.16 A3	0
2	4-OH	7	4	1.44	-3.56(in Log(moles/L)94.63 (in mg/L)	110.60 A2	284.70 A3	0
3	4-Cl	6	3	1.6	-3.86(in Log(moles/L)94.63 (in mg/L)	111.06A2	264.16 A3	0
4	4-NH₂	6	5	1.2	-3.88(inLog(moles/L)) 45.80 (in mg/L)	113.79 A2	281.58 A3	0
5	4-OH,3-OCH₃	8	4	1.4	-3.83(inLog(moles/L)) 55.73 (in mg/L)	117.16 A2	317.35 A3	0
6	4-CH(CH₃)₂	6	3	2.82	-5.18(inLog(moles/L)) 2.43 (in mg/L)	92.98 A2	327.14 A3	0
7	4-C(CH₃)₃	6	3	3.31	-5.84(inLog(moles/L)) 0.56 (in mg/L)	92.98 A2	358.69 A3	0
8	4-N(CH₃)₂	6	3	1.82	-4.28(inLog(moles/L)) 19.73 (in mg/L)	95.79 A2	323.64 A3	0
9	4-NO₂	8	3	1.37	-4.80(inLog(moles/L)) 6.00 (in mg/L)	126.37 A2	299.07 A3	0
10	3,4,5-(OCH₃)₃	9	3	1.73	-3.87(inLog(moles/L)) 56.23 (in mg/L)	115.96 A2	369.19 A3	0

Log P-Partition coefficient, HBA-Hydrogen bond acceptor, HBD-Hydrogen bond donor

Druglikeness and bioactivity scores

Prediction of bioactivity score for the most important drug targets (GPCR ligands, kinase inhibitors, ion channel modulators, nuclear receptors. Drug likeness may be defined as a complex balance of various molecular properties and structure features which determine whether particular molecule is similar to the known drugs. Table 7 and 8 lists the predicted values of selected parameters for compounds. Positive scores were obtained for the compounds indicating the potential quinazolinone nucleus and amino acid residues. The presence of the two pharmacophoric scaffolds, oxadiazole and sulphonamide influenced the good bioactive scores. Compounds 6 & 7 with electron donating substituents showed positive scores as protease inhibitor and enzyme inhibitors. The compounds showed positive drug likeness scores which is a good indication of being active compounds. Among the compounds 6 & 7 were found to be potent with good druglikeness scores nearer to 1 (Table 7, Table 8).

Table 7: Bioactivity prediction of title compounds using molinspiration.com

S.No	R	GPCR Ligand	Ion Channel Modulator	Kinase Inhibitor	Nuclear Receptor Ligand	Protease Inhibitor	Enzyme Inhibitor
1	H	-0.18	-0.4	-0.12	-0.58	0.06	-0.07
2	4-OH	-0.13	-0.34	-0.08	-0.4	0.09	-0.02
3	4-Cl	-0.17	-0.39	-0.14	-0.57	0.03	-0.11
4	4-NH₂	-0.17	-0.38	-0.12	-0.54	0.07	-0.07
5	4-OH,3-OCH₃	-0.18	-0.39	-0.07	-0.47	-0.02	-0.05
6	4-CH(CH₃)₂	-0.16	-0.37	-0.17	-0.45	0.07	0.07
7	4-C(CH₃)₃	-0.12	-0.29	-0.12	-0.38	0.07	0.06
8	4-N(CH₃)₂	-0.15	-0.37	-0.09	-0.48	-0.06	-0.08
9	4-NO₂	-0.3	-0.4	-0.26	-0.6	0.06	-0.17
10	3,4,5-(OCH₃)₃	-0.21	-0.4	-0.11	-0.56	-0.04	-0.11

Table 8: Druglikeness score prediction of title compounds using Molsoft.com

S.No	R	Druglikeness score
1	H	0.41
2	4-OH	0.24
3	4-Cl	0.33
4	4-NH₂	0. 50
5	4-OH,3-OCH₃	0.74
6	4-CH(CH₃)₂	0.93
7	4-C(CH₃)₃	0.83
8	4-N(CH₃)₂	0.48
9	4-NO₂	0.36
10	3,4,5-(OCH₃)₃	0.74

Docking studies:

The virtual screening of drug like molecules against a protein target is a common strategy used to identify novel inhibitors. The designed compounds exhibited good binding interactions with Mur enzymes, specifically with x-ray crystallographic structures of UDP-N[Mur A]PDB ID:1a2n, UDP-N[Mur C]PDB ID:1GQ7, UDP-N [Mur E] PDB ID:2XJA. MurA is the only cytoplasmic step inhibited by a clinically used antibacterial agent. Fosfomycin, a naturally occurring broad spectrum antibiotic, is the best known inhibitor of MurA. [22]. It has been the drug of choice for the treatment of paediatric gastrointestinal infections resulting from Shiga-like toxin-producing Escherichia coli (STEC). It is also among the first-line agents for the treatment of bacterial infections of the urinary tract, which is a common health problem, particularly in women [23]. Inhibition of MurA enzyme by fosfomycin is competitive, the antibiotic acting as an analogue of PEP and forming a covalent bond with the active cysteine residue of the enzyme (Figure 3).

Figure 3:Title compounds selected for docking studies

Mur C enzyme catalyzes the stepwise formation of the peptide moiety of the peptidoglycan disaccharide peptide monomer unit. MurC is responsible of the addition of the first residue (L-alanine) onto the nucleotide precursor UDP-MurNAc. The first phosphinate inhibitor of MurE was designed with structural features based on Tanner’s previously reported MurD phosphinate inhibitor [23]. Several analogues of diaminopimelic acid were also tested as substrates or inhibitors of MurE [22].

Table 9 indicates the docking scores, binding energies and compared with the standard. All the compounds exhibited good binding energies, of which the propyl and t-butyl substituted phenyl ring containing compounds showed stronger interactions.

The detailed interaction of active compounds with Mur A, Mur C and Mur E are given in tables 10, 11 and 12 respectively. Most of the interactions of the compounds with amino acids at active site of the enzymes were due to hydrogen bonds, van der Waals and hydrophobic attractions. Compounds 6 [ΔG=-9.8(Kcal/mol)] and compound 9 [ΔG= -9.5(Kcal/mol)] were found to have good interactions at active site of Mur A enzyme than the standard Mur A inhibitor, Fosfomycin [ΔG= -4.8(Kcal/mol)] (Figure 4). The nitro group increased its efficiency or conformation of compounds towards active interaction. Schonbrunn et al reported similar type of interaction of the fluorescence probe 8-anilino-1-naphthalene sulfonate (ANS) with the antibiotic target MurA [25].

Figure 4: Hypothetical binding interactions of selected compounds and hydrogen bonds as black lines with bond length in the UDP-N [Mur A] (PDB ID:1a2n) active site with amino acids(wire model) (a) Compound Fosfomycin (blue colored stick model) Carbon(light blue), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white) (b) compound 9 (light brown colored ball and stick model) Carbon(light brown), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white) (c) Superimposition of compounds at the active site.(d) superimposition of the compounds in hydrophobic pocket Carbon(light blue, light brown), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white)

Compound 6 [ΔG=-9.7(Kcal/mol)] and compound 7[ΔG= -9.6(Kcal/mol)] were found to have good interactions at active site of Mur C enzyme than the standard Mur C inhibitor, Diamino pimilic acid Phosphinate [ΔG= -5.2 (Kcal/mol)](Figure 5). The active site contains amino acids as reported by Spraggon et al and the compounds exhibited strong interactions [25]. Compounds 6 [ΔG=-9.8(Kcal/mol)] and compound 7 [ΔG= -9.7(Kcal/mol)] were found to have good interactions at active site of Mur E enzyme than the standard Mur E inhibitor, Diamino pimilic acid Phosphinate [ΔG= -5.2 (Kcal/mol)] (Figure 6).

Of all the compounds, compound 6 and compound 9 with isopropyl substitution exhibited the highest scores against all the three Mur enzymes. Compound 7 showed good binding interactions with Mur C and Mur E. The Mur E enzyme of Mur pathway of Mycobacterium tuberculosis is an attractive drug target as it is unique to bacteria and is absent in mammalian cells [26].

This indicates the high specificity and good interactions of the compounds at binding pockets of these Mur enzymes. Docking studies were also conducted on Alpha amylase PDB ID: 3BC9 and Dihydropteroate synthase PDB ID:1ad4, but the compounds exhibited only mild interactions. Thus the compounds were found to be selective towards Mur enzymes (Table 9, Table 10, Table11, Table 12) (Figure 4, Figure 5, Figure 6).

Figure 5: Hypothetical binding interactions of selected compounds and hydrogen bonds as black lines with bond length in the UDP-N[Mur C](PDB ID:1GQ7) active site with amino acids(wire model) (a) Compound phosphinate (violet colored ball and stick model) Carbon(violet), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white) (b) compound 6 (light pink colored ball and stick model) Carbon(light pink), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white) (c) compound 7 (green colored ball and stick model) (d) superimposition of the compounds in hydrophobic pocket,Carbon(violet, light pink), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white)

Figure 6: Hypothetical binding interactions of selected compounds and hydrogen bonds as black lines with bond length in the UDP-N[Mur E]PDB ID:2XJA active site with amino acids(wire model) (a) compound 7 (silver colored ball and stick model) and phosphinate(light blue colored ball and stick model) Carbon(silver), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white) (b) and (c) compound 6 (cyan colored ball and stick model) Carbon(cyan), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white) (d) superimposition of the compounds in hydrophobic pocket Carbon(silver, cyan), Sulfur(yellow), Oxygen(red), Nitrogen(blue) and Hydrogen (white)

Table 9: Docking studies of 1,3,4-oxadiazolyl sulphonamides

S.No	R
		Targets and Estimated binding energy ΔG(Kcal/mol) a
		UDP-N [Mur A] PDB ID:1a2n	UDP-N [Mur C] PDB ID:1GQ7	UDP-N [Mur E] PDB ID:2XJA	Alpha amylase PDB ID:3BC9	Dihydro pteroate synthase PDB ID:1ad4
		UDP-N [Mur A] PDB ID:1a2n	UDP-N [Mur C] PDB ID:1GQ7	UDP-N [Mur E] PDB ID:2XJA
1	H	-9	-9.2	-9.3	-7.5	-6.6
2	4-OH	-9.1	-9.1	-9.2	-7.7	-6.4
3	4-Cl	-9.3	-9.1	-8.3	-7.5	-6.9
4	4-NH₂	-9.2	-9.1	-8.2	-7.6	-6.8
5	4-OH, 3-OCH₃	-9.1	-9.1	-9.3	-7.5	-7
6	4-CH(CH₃)₂	-9.8	-9.7	-9.8	-7.9	-6.8
7	4-C(CH₃)₃	-8.8	-9.6	-9.7	-7.8	-6.6
8	4-N(CH₃)₂	-9	-9	-9.5	-7.6	-6.4
9	4-NO₂	-9.5	-9.5	-9.5	-7.7	-6.9
10	3,4,5-(OCH₃)₃	-8.7	-9.1	-8.3	-7.3	-6.8
11	Standard	-4.8 Fosfomycin	-5.2 Diamino pimilic acid Phosphinate	-5.2 Diamino pimilic acid Phosphinate
11	Standard	-4.8 Fosfomycin	-5.2 Diamino pimilic acid Phosphinate	-5.2 Diamino pimilic acid Phosphinate		Sulfadiazine -4.7

Table 10: Binding interactions of active compounds with UDP-N[Mur A]PDB ID:1a2n (His-Histidine, Pro-Proline, Ser-Serine, Thr-Threonine, Gly-Glycine, Val-Valine, Arg-Arginine, Leu-Leucine)

Compound	Hydrogen bond interactions	Other interactions
Fosfomycin	Phosphate(OH)-His_125A(~3.128 Ǻ)	van der Waals interactions of methyl group with Leu_124A and Arg_120A
Fosfomycin	Phosphate(OH)-Pro_123A(~2.155 Ǻ) Cyclic(O)-Ser-_162A(~2.289 Ǻ)
Compound 6	OxadiazoleN1-Gly _164A(~1.996 Ǻ), N-Ser_162A(~2.071 Ǻ)	Hydrophobic aromatic pockets-(His_125B, Pro_303A, Pro_293A) van der Waals interactions -(Gly_164A,Val_163A, Thr_304A)
Compound 9	NO₂- Thr_304A(~2.342 Ǻ), Oxadiazole N- Val_163A(~2.094Ǻ) Sulfonamide SO₂(O)- His_125A(~3.229 Ǻ) Sulfonamide SO₂(O)- Arg_91A(~2.454 Ǻ)	van der Waals interactions -(Gly_164A,Val_163A,Thr_304A) Hydrophobic aromatic pockets-(His_125B, Pro_303A,Pro_293A)

Table 11: Binding interactions of active compounds with UDP-N[Mur C]PDB ID:1GQ7(His-Histidine, Pro-Proline, Ser-Serine, Thr-Threonine, Gly-Glycine, Val-Valine, Arg-Arginine, Lys-Lysine, Tyr-Tyrosine)

Compound	Hydrogen bond interactions	Other interactions
Phosphinate	-	van der Waals interactions of methyl group with Thr _120B,Gly_123B,Gly_342B and Arg_321B, Aromatic interactions with His _286B
Compound 6	Sulfonamide NH₂- His _343B (~2.593 Ǻ)SO₂-Arg ₃₂₁(~1.975Ǻ), SO₂-Thr_126B(~2.2768 Ǻ)	Hydrophobic aromatic pockets-( Tyr_341B, His_122B,His_343B, His_286B) Van der Waals interactions -( Lys _124B, Gly_123B,Gly_342B,Gly_347B,Val _250B)
Compound 7	NH₂ sulfonamide-His _286B(~3.272 Ǻ), Thr _351B(~2.151Ǻ) Oxadiazole O-Arg_321B(~1.972 Ǻ), N-Lys _124B(~1.957Ǻ)	van der Waals interactions -( Lys _124B, Gly_123B, Gly_342B,Gly_347B,Val_250B, Thr_125B, Thr_126B) Hydrophobic aromatic pockets-(His_125B, Pro_303A,Pro_293A)

Table 12: Binding interactions of active compounds with UDP-N[Mur C]PDB ID:1GQ7(His-Histidine, Pro-Proline, Ser-Serine, Thr-Threonine, Gly-Glycine, Val-Valine, Arg-Arginine, Lys-Lysine, Tyr-Tyrosine)

Compound	Hydrogen bond interactions	Other interactions
Phosphinate	Phosphate (O)-Ser_123D(3.276)	Hydrophobic aromatic pockets(Tyr_311D), van der Waals interactions of methyl group with Ser_370D, Ser_123D, Leu _346D and Thr_127D
Compound 6	NH₂ of sulphonamide-Lys_125D (~3.536 Ǻ) , Thr_148D (~3.069Ǻ) SO₂ with Ser_123D (~3.291 Ǻ)	Hydrophobic aromatic pockets-(Tyr _311D and Tyr _361D) van der Waals interactions -(Ser _370D, Ser_123D,Leu_346D, Thr_127D and Arg_345D)
Compound 6	Oxadiazole N-Thr_127D (~3.276 Ǻ) , N-Asn_315D(~3.329 Ǻ)
Compound 7	Oxadiazole N-Ser_123D (~3.276 Ǻ) N-Asn_315D (~3.329 Ǻ) Sulfonamide(NH)-Arg 3_45D (~3.334 Ǻ)	van der Waals interactions -(Ser_370D, Ser_123D,Leu_346D, Thr_127D and Arg_345D ) Hydrophobic aromatic pockets-(Tyr_311D and Tyr_361D)

Conclusion

The search for new antibacterial agents directed towards novel targets became highly imperative. Inspired by the literature survey about the potentiality of oxadiazole and sulphonamide moieties derivatives, a series of 4-Amino-N-[(5-phenyl-1,3,4-oxadiazol-2-yl) methyl] substituted benzene-1-sulfonamides were designed. In silico evaluation including molecular properties like bioavailability, Lipinski’s rule of five reflected the oral bioavailability of these compounds. Good druglikeness and bioactivity scores further indicated the probable potentiality of these compounds as future drugs. The presence of the two pharmacophoric scaffolds, oxadiazole and sulphonamide influenced the good bioactive scores. Compounds 6 & 7 with electron donating substituents showed positive scores as protease inhibitor and enzyme inhibitors.

Docking studies further supported the biological activities. The results indicated the possible inhibitory activity towards the Mur enzymes. The results also indicate the selectivity towards Mur A, Mur C, & Mur E enzymes. Of all the compounds, compound 6 and compound 9 with isopropyl substitution exhibited the highest scores against all the three Mur enzymes. Compound 7 showed good binding interactions with Mur C and Mur E. Further research in this approach may be extended to decrease the in vitro evaluation on these enzymes. In addition, the synthesis of these compounds was also observed to be feasible. Thus this study will be promising in finding novel antimicrobial agents with the Mur enzyme inhibitory activity (or) peptidoglycan synthesis inhibitory activity and therapeutic potential as novel anti-TB drugs.

Acknowledgement

I sincerely thank the management of Krishna Teja Pharmacy College, Tirupati for providing necessary facilities for the above research work.

ReferencesTop

Hansch C, Sammes PG, Taylor JB. Comprehensive Medicinal Chemistry. 1990; Vol. 2, Pergamon Press: Oxford, Chap. 7.1.
Kanda Y, Kawanishi Y, Oda K, Sakata T, Mihara S, Asakura K, et al. Synthesis and structure-activity relationships of potent and orally active sulfonamide ETB selective antagonists. Bioorg. & Med. Chem. 2001;9(4):897-907.
Stokes SS, Albert R, Buurman Ed T, Andrews B, Shapiro AB, Green OM, McKenzie, AR, Otterbein LR. Inhibitors of the acetyltransferase domain of N-acetylglucosamine-1-phosphate-uridylyl transferase/ glucosamine-1-phosphate acetyl transferase (GlmU). Part 2: Optimization of physical properties leading to antibacterial aryl sulfonamides. Bioorg. & Med. Chem. Lett, 2012; 22: 7019. (PMID:23099094)
Rahavi Ezabadi, Camoutsis C, Zoumpoulakis P, Geronikaki A, Sokovic M, Glamocilija J, et al. Sulfonamide-1,2,4-triazole derivatives as antifungal and antibacterial agents: Synthesis, biological evaluation, lipophilicity, and conformational studies. Bioorg. & Med. Chem, 2007;16(3):1150-1161. Doi: 10.1016/j.bmc.2007.10.082.
Chibale K, Haupt H, Kendrick H, Yardley V, Saravanamuthu A, Fairlamb A H, et al. Antiprotozoal and cytotoxicity evaluation of sulfonamide and urea analogues of quinacrine. Bioorg. & Med. Chem. Lett, 2001;11(19):2655-2657.
Roush WR, Cheng JM, KnappReed B, Alvarez Hernandez A, McKerrow JH, Hansell E, et al. Evaluation of Sulfonamide Derivatives of Dagenan Chloride as Lipoxygenase and α-Glucosidase Inhibitors. Bioorg. Med. Chem. Lett. 2001;11:2759-2762.
Serradeil-Le Gal C. An overview of SR121463, a selective non-peptide vasopressin V2 receptor antagonist. Cardiovascular Drug Rev. 2001,19(3):201-14.
Natarajan A Guo, Harbinski F, Fan YH, Chen H, Luus L, Dierck J, Aktas H, Chorev M, Halperin JA. Novel arylsulfoanilide−oxindole hybrid as an anticancer agent that inhibits translation initiation. J. Med. Chem. 2004;47(21):4979-4982. doi:10.1021/jm0496234.
Brain CT, Paul JM, Loong Y, Oakley PJ. Novel procedure for the synthesis of 1,3,4-oxadiazoles from 1,2-diacylhydrazines using polymer-supported Burgess reagent under microwave conditions. Tetrahedron Lett. 1999;40:3275-3278. doi: 10.1016/S0040-4039(99)00382-2
Chandrakantha B, Prakash Shetty, Vijesh N, Nishitha I, Arun M Isloor. Synthesis, characterization and biological activity of some new 1,3,4-oxadiazole bearing 2-flouro-4- methoxy phenyl moiety. Eur J Med Chem. 2010;45:1206–1210. doi: 10.1016/j.ejmech.2009.11.046
Sangshetti JN, Chabukswar AR, Shinde DB. Microwave assisted one pot synthesis of some Novel 2,5-disubstituted 1,3,4-oxadiazoles as antifungal agents. Bioorg. Med Chem Lett. 2011;21(1):444–448. Doi :10.1016/j.bmcl.2010.10.120
de Oliveira CS, Lira BF, Barbosa-Filho JM, Lorenzo JG, de Athayde-Filho PF. Synthetic approaches and pharmacological Activity of 1,3,4-Oxadiazoles: A Review of the Literature from 2000–2012. Molecules. 2012;17(9):10192-10231. Doi:10.3390/molecules170910192
Kumar B, Kumar A, Beheraand AK, Raj V. Latest Update on Pharmacological Activities of 1,3,4-Oxadiazole Derivatives. J Cell Sci Ther, 2016;7:233. doi:10.41 72/2157-7013.1000233
Om Prakash, Manoj Kumar, Rajesh Kumar, Chetan Sharma, Aneja KR. Hypervalent Iodine(III) mediated synthesis of novel unsymmetrical 2,5-disubstituted 1,3,4 oxadiazoles as antibacterial and antifungal agents. Eur J Med Chem. 2010;45:4252-4257. Doi:10.1016/j.ejmech.2010.06.023
Sudhamani H, Thaslim Basha Sk, Venkateswarlu N. Synthesis and characterization of new thiourea and urea derivatives of 6-fluoro-3-(piperidin-4-yl) benzo[d] isoxazole: In vitro Antimicrobial and Antioxidant activity, J Chem Sci. 2015;127(10):1739-1746. doi:10.1007/s12039-015-0935-6
Chopra I, Storey C, Falla TJ, Pearce JH. Antibiotics, peptidoglycan synthesis and genomics: the chlamydial anomaly revisited. Microbiology. 1998;144: 2637-2678. Doi:10.1099/00221287-144-10-2673.
Danial N N, Korsmeyer S J. Cell death: critical control points, Cell. 2004;116(2):205-219.
Douglas B Kitchen, Helene Decornez, John R Furr, Jurgen Bajorath. Docking and scoring in virtual screening for drug discovery: methods and applications. Nat Rev Drug Discov. 2004;3: 935-949. doi:10.1038/nrd1549
Brooijmans N, Kuntz I D. Molecular recognition and docking algorithms. Ann Rev Biophys Biolmol. Struct. 2003;32:335–373. Doi: 10.1146/annurev. biophys.32. 110601. 142532.
Lipinski CA, LombardoF, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug. Del Rev.1997;23:3–25. Doi: 10.1016/S0169-409X(96)00423-1
Lalitha P, Sivakamasundari S. Calculation of molecular lipophilicity and drug likeness for few heterocycles. Oriental J Chem. 2010;26:135-141
Auger G, van Heijenoort J, Vederas JC, and Blanot D. Effect of analogues of diaminopimelic acid on the meso-diaminopimelate-adding enzyme from Escherichia coli. FEBS Lett. 1996;391:171–174. Doi: 10.1016/S0014-5793(96)90569-4
Nicolle LE. Urinary tract infection: traditional pharmacologic therapies. Am J Med. 2002;113(1A): 35S–44S. doi: 10.1067/mda.2003
Zeng B, Wong KK, Pompliano DL, Reddy S, Tanner ME. A phosphinate inhibitor of the mesodiaminopimelic acid-adding enzyme (MurE) of peptidoglycan biosynthesis. J Org Chem. 1998;63(26):10081–10086.
Schonbrunn E, Eschenburg S, Luger K, Kabsch W, Amrhein N. Structural basis for the interaction of the fluorescence probe 8-anilino-1-naphthalene sulfonate (ANS) with the antibiotic target MurA. Proc Natl Acad Sci USA. 2000;97: 6345–6349. doi: 10.1073/pnas.120120397
Spraggon G1, Schwarzenbacher R, Kreusch A, Lee CC, Abdubek P, Ambing E, et al. Crystal structure of an UDP-N-acetylmuramate-alanine ligase MurC (TM0231) from Thermotoga maritima at 2.3 A ° resolution. Proteins, 2004; 55(4):1078–1081. Doi: 10.1002/prot.20034
Strancar K, Boniface A, Blanot D, Gobec S. Phosphinate inhibitors of UDP-N-acetylmuramoyl-L-alanyl-D-glutamate: L-lysine ligase (MurE). Arch Pharm Chem Life Sci. 2007;340(3):127–134. Doi: 10.1002/ardp.200600191
Chemaxon available at http:// www. chemaxon. com/ products/marvin/marvinsketch. [last accessed on 2016 April 14]. Marvin Sketch 16.11.21, Developed by Chemaxon worldwide (USA, Europe, Asia) Pune, India, Advent Informatics, , info@adventinformatics.com
Dogne JM, Supuran CT, Pratico DJ. Adverse Cardiovascular Effects of the Coxibs. J Med Chem. 2005;48: 2251-2257. Doi: 10.1021/jm0402059
Entrl P, Muhlbacher J, Rohde B, Selzer P. Web-based cheminformatics and molecular property prediction tools supporting drug design and development at Novartis. SAR QSAR Environ Res, 2003;14(15-6):321-328.Doi:10.1080/ 1062936031000167391
Molinspiration Available at http://www.molinspiration.com [last accessed on 2016 Feb 8]. Founded in 1986 as a spin-off of Bratislava University, Nova ulica, SK-900 26 Slovensky Grob, Slovak Republic, used Molinspiration Fragment-based Virtual Screening Engine v2016.02
Molsoft Available at http://molsoft.com/mprop/ http://molsoft.com/mprop/ [last accessed on 2016 Oct 4]. Version 3.8-5,Developed by Molsoft L.L.C. 11199 Sorrento Valley Road, S209 San Diego CA 92121, info@molsoft.com
Rathish IG, Javed K, Ahmad S, Bano S, Alam MS, Pillai KK, Singh S, Bagchi V. Bioorg. Synthesis and antiinflammatory activity of some new 1,3,5-trisubstituted pyrazolines bearing benzene sulfonamide. Med Chem Lett. 2009;19(1):255-258. Doi: 10.1016/j.bmcl.2011.05.061
Swiss Dock available at http://swissdock.vital-it.ch /docking #.[last accessed on 2016 Jun 25]. a free protein ligand docking web service powered by EADock DSS by the Molecular Modeling group of the Swiss Institute of Bioinformatics
Van Heijenoort J. Recent advances in the formation of the bacterial peptidoglycan monomer subunit. Nat Prod Rep. 2001; 18(5):503–519. Doi: 10.1039/A804532A