2822
K. S. Ghosh et al. / Bioorg. Med. Chem. 16 (2008) 2819–2828
in the presence of inhibitors to monitor their behavior.
From the H NMR titration curve of the histidine C-2
1
proton (Fig. 4A), the observed pKa values for His 12,
His 105, and His 119 are 5.88, 6.59, and 6.22, respec-
tively, for free RNase A. When compound 5 and 6 com-
plex with RNase A (Fig. 4B and C, respectively), the
pKa values shift to 6.01, 6.62 and 6.77 for compound 5
and 6.23, 6.53, and 6.39 for compound 6, respectively.
Perturbation of the pKa values of His 12 and His 119
indicate that compounds 5 and 6 go to the active site
of RNase A. We observe that the pKa values for His
105 remain practically unaltered (<1%) which is also
similar to the reported change of pKa for His 105 in case
of 20-CMP22 and 30-CMP.23 The pKa value for His 119
shifts by ꢀ9% in the presence of compound 5. The in-
crease in the pKa value indicates that the presence of
the negative charge center has a direct influence on the
environment of the imidazole ring of His 119. An in-
crease in pKa is also observed when the negatively
charged phosphate group interacts with the His residues
for complexes with both 20-CMP22 and 30-CMP.23 In the
case of phosphate esters, however, there is a substantial
change (ꢀ28%). This can be attributed to the formation
of the dianion that is capable of perturbing the environ-
ment to a higher degree than compound 5 or 6, where
only a negatively polarized charge center on the car-
bonyl oxygen is resident. Interestingly, the pKa of histi-
dines are unaffected in the cytidine–RNase A complex,22
where such interactions of the phosphate dianion with
the His residues are absent. Our results thus indicate
that the inhibition of RNase A by compound 5 is pri-
marily because of its interaction with His 119. In case
of compound 6, relatively minor changes are observed
for His 12 (ꢀ6%) and His 119 (ꢀ3%). The presence of
the phenyl ring in this case could interact with Phe 120
resulting in a smaller perturbation to the His residues
but nevertheless occupying the active site. The abstrac-
tion of the proton from 20-OH by His 12 (acting as base)
for cytidine phosphates is not also possible for com-
pound 5 or 6, which has no 20-OH that explains the
much lower DpKa observed for the 30-O-carboxy esters.
The preferred docking poses for both the compounds
are consistent with these changes in pKa (discussed
later).
Figure 5. FT-IR difference spectra of RNase A complexes with
compound 5 (–•–) and compound 6 (ÆÆÆ) at RNase A:ligand as 1:1
(left hand axis) with FT-IR spectra of free RNase A (—) (right hand
axis).
tents, respectively. Similarly for compound 6 at the same
protein ligand ratio, difference spectra reflect changes in
a-helix (1662 cmꢁ1), random structure (1643 cmꢁ1), and
b-sheet (1635 cmꢁ1) content. The increase in intensity
observed in the difference spectra between 1549 and
1548 cmꢁ1 in the amide II region for compounds 5
and 6, respectively, can be attributed to the interaction
of the compounds with the backbone of the protein.
The distinctive changes in the content of secondary
structural elements are more prominent in the deconvo-
luted spectra of RNase A and its complexes with com-
pounds 5 and 6 (Fig. 6).
CD spectroscopic studies of free RNase A and its com-
plexes with compounds 5 and 6 also indicate a perturba-
tion in the secondary structure (Fig. 7). The percentage
of secondary structural elements for RNase A and its
complexes with compounds 5 and 6 at 1:0.5 and 1:1 ra-
tios was determined using SELCON3 method in
DICHROWEB.27 We find that there is an increase in
the a-helix content (ꢀ9%) with a concomitant reduction
in random (ꢀ5%) and turn structures (ꢀ4%) for com-
plexation with compound 5 at the higher ratio. An in-
crease in a-helix content of RNase A has also been
observed with the nucleoside 3-azido-3-deoxythymidine
(AZT).28 In case of compound 6 at the higher ratio, min-
or increases in both the a-helix and b-sheet contents
(ꢀ6% and 4%, respectively) are indicative more of con-
formational adjustments on complex formation rather
than a major change in secondary structure. This was
compensated by a decrease in the content of random
structure by ꢀ10%. The stabilization of b-sheet confor-
mation of RNase A has also been observed in the inter-
action of RNase A with (ꢁ)-epigallocatechin.15
The effect of binding of these compounds on the overall
secondary structure of RNase A was also investigated
by FT-IR and CD spectroscopy. In FT-IR spectros-
copy, amide I bands at 1645–1650 cmꢁ1 (mainly C@O
stretching) and the amide II band at 1548–1560 cmꢁ1
(C–N stretching coupled with an N–H bending mode)
of proteins change on interaction with small molecules
and ligands.24 The fitted Gaussian curves in the amide
I region with spectral ranges from 1610 to 1632 cmꢁ1
,
1636 to 1644 cmꢁ1, 1650 to 1662 cmꢁ1, and 1665 to
1680 cmꢁ1 are attributed to the b-sheet, random coil,
a-helix and turn structures, respectively.25,26 Changes
in the peak positions of the amide I and II bands in
the FT-IR spectra of RNase A are observed with both
compounds. This is visible in the difference spectra
(Fig. 5). With compound 5 at a protein ligand ratio of
1:1 the peaks at 1663, 1643, and 1629 cmꢁ1 indicate
changes in a-helix, random structure, and b-sheet con-
Protein–ligand docking studies were performed to ob-
tain some insight into the amino acid residues of RNase
A involved in interactions with compounds 5 and 6.
Though the crystal structure of the complex can repre-
sent specific details of the interactions, general observa-
tions may be obtained from docking studies. Our earlier
studies13 with another synthetic inhibitor of RNase A
revealed that docking poses are similar to those ob-
˚
tained from structural studies (RMSD 0.47 A). Interest-