F.A. Khan, A.U. Rehman and M.A. Abbasi et al. / Journal of Molecular Structure 1223 (2021) 129141
3
buffer (4 mM, pH 6.8) containing urea (100 mM) was added in
each well and plate was incubated at 30 °C for additional 15 min.
Activity of urease was determined by the measurement of the
amount of NH3 formed during urease catalysis as mentioned by
Weatherburn [29]. After that, each well of the plate was inoculated
with 45 μL of phenol reagent (1% w/v phenol + 0.005% w/vNa2-
nitroprusside) and 70 μL of alkali reagents (0.5% w/v sodium
hydroxide+ 0.1% sodium hypochlorite). Finally, after 50 min the in-
crease in absorbance was measured at 630 nm by the use of mi-
croplate readers (Molecular Devices., CA, USA). The entire reactions
were performed (in triplicate) till to achieve a final volume 200 μL.
Analysis of the results was carried out by the SoftMax Pro6.3 soft-
ware (Molecular Device, USA). Thiourea was used as standard by
applying this method. The percentage (%) inhibition was measured
by using the formula;
which is important for structure-based drug designing. These re-
sults were then correlated with the experimental values. Molecular
Operating Environment (MOE 2009–2010, Quebec, Canada) docking
software was used to perform such studies (MOE. 2010). Ligands as
well as the receptor protein (urease) preparation are the two basic
steps prior to perform the docking.
(і) Ligands preparation
The 3D structures of all the ligands were sketched through
ChemDraw Ultra 12.0 (Cambridge Soft-2001, Cambridge, USA)
(Cambridge Soft. 2001) and saved in the mol format. The com-
pounds scaffolds were 3D protonated by Protonate 3D Option
in MOE software. The energies of identified ligands were mini-
mized by applying the default parameters of energy minimization
algorithm already adjusted in MOE (gradients: 0.05, force field:
MMFF94x). The Database of compounds series was created in mdb
file format.
%age inhibition = 100-(ODtest well/ODcontrol) × 100.
(іі) Preparation of receptor protein
The 3D structure of urease (4UBP) having resolution of 1.55 A,
˚
2.4. Determination of kinetic parameters
was retrieved from Protein Data Bank (PDB). H2O molecules from
the receptor protein were extracted. Like ligands, Protonate 3D op-
tion was applied for the protonation of urease 3D structure. Simi-
larly, the energy of urease protein was also minimized by applying
the default parameters of energy minimization algorithm with gra-
dients of 0.05 and force field Amber99, and saved in pdb file for-
mat. Later on, docking of all ligands was carried out in the binding
pocket of urease by following the default parameters of MOE-Dock
Program. In order to increase the accuracy of protocol, re-docking
was repeated [31]. After complete docking of all the ligands, the
most excellent 2D as well as 3D interaction images were chosen to
elucidate their specific types of interactions and bond lengths re-
spectively. All 3-D figures were made by discovery studio visualizer
The IC50 inhibitory potential of compounds is the concentration
of the test compounds that blocked the urea (substrate) hydrolysis
values represent by 50% was measured by monitoring the effects of
their different concentrations (from high to low) in the assay. Such
calculation was carried out by using the EZ-Fit Enzyme Kinetics
Program (Perrella Scientific Inc., Amherst, USA). The ES represents
the complex of Jack bean urease and urea, while P stand for the
product obtained as a result of reaction. The Lineweaver-Burk plot
was then applied to identify the type of inhibition, while dissocia-
tion and inhibition constants (Ki) values were calculated by apply-
ing the Dixon plots and their secondary replots [30]. Non-linear re-
gression equation was applied to find out the Km, Ki, and Vmax val-
ues. At first, the values of 1/Vmaxapp were calculated at each inter-
section points of lines of every inhibitor concentration on y-axis of
the Lineweaver-Burk plot and then replotted against various con-
centrations of inhibitor. In the subsequent proceeding, the slope of
each line on the Lineweaver-Burk plot obtained as a result of in-
hibitor concentration was plotted against various concentrations of
inhibitor. Competitive inhibition can be classically represented as;
2.7. Cytotoxicity evaluation
The cytotoxicity of various potential urease inhibitors were
tested against neutrophil cell lines. The following steps were used
for the evaluation of cytotoxicity.
(і) Isolation of human neutrophils
Heparinized fresh venous blood was taken from the young and
healthy volunteers in a city blood collection center. Isolation of
neutrophils was carried out through Siddiqui et al., 1995 method
[33]. Accordingly, mixing of whole blood was made with Ficoll-
Paque (Pharmacia Biotech Amersham, Uppsala). When sedimenta-
tion occurred, the unnecessary red blood cells (RBCs) were layered
in a buffy coat way on a 3.0 mL cushion of Ficoll. Centrifugation
was then carried out for 30 min at a rate of 1500 rpm. Super-
natant was discarded, while pellets were collected. Furthermore, it
was mixed with 0.83% of hypotonic ammonium chloride solution
in order to lyse the RBCs. Again, the solution was centrifuged and
Modified Hank’s Solution (MHS) was used for the washing of col-
lected neutrophils. Later on, resuspension of neutrophil cells was
performed in the same solution at the rate of 1 × 107 cells per
mL.
Scheme 2. Classical representation of competitive inhibition.
2.5. Statistical analysis
Each experiment was performed in triplicate and results were
expressed as a mean of three. SoftMax Pro 6.3 Software (Molecular
Devices, CA, USA) was used for the analysis of these results. The
graph plotting was carried out through GraFit program. By using
the same program; we have obtained the values of correlation co-
efficients, intercepts, slopes, and their standard errors from linear
regression analysis. The relationship for each and every line in all
the graphs was more than 0.99.
(іі) Assay procedure for neutrophil based cytotoxicity
Isolated human neutrophils (1 × 107cells per mL) were treated
first with compounds under trial for thirty minutes and then with
0.25 mM WST-1(Dojindo Laboratories Kumamoto, Japan) on shak-
ing water bath at 37 °C [34]. Change in absorbance was calculated
after incubation for 3 h, at 450 nm using 96-well plate reader
(Spectra-MAX-340, Molecular Devices, CA, USA). Here the OD rep-
resents the mean of 5 investigational replicates. Percent (%) viabil-
ity of cells was measured as:
2.6. Molecular docking
%age viability of cells = {(OD test compound × 100/OD control)
-100}–100
Molecular docking was performed for prediction the orienta-
tion of protein active site residues and potential inhibitors binding,