G Model
CCLET 4579 No. of Pages 4
X. Li et al. / Chinese Chemical Letters xxx (2018) xxx–xxx
3
Fig. 3. Fluorescence spectra of TAR RNA-ICR upon the addition of complexes 1a (A) and 2b (B). The arrows show the change upon the increasing amount of complex
concentration. Inset: plot of F0/F vs. [complex] for the titration of the complex to ICR-191. Plot of lg(F–F0)/F0 vs. lg[[Dt] – n(F–F0)/F [Rt]] for the titration of complexes 1a (C) and
2b (D) to ICR. Binding assay conditions: The fluorescence displacement experiments were performed in Na2HPO4-NaH2PO4 buffer solution (20 mmol/L, pH 7.0 at 25 ꢃC,
containing 20 mmol/L NaCl) with the concentrations of both ICR 191 and TAR RNA at 1 mmol/L. All experiments were performed on a fluorescence spectrophotometer Hitachi
F7000. The excitation wavelength for ligand ICR 191 was set at 411 nm, emission wavelength was set at 420 nm, scaning range 420-600 nm. The slit width was set at (5, 5) nm,
scanning speed 1200 nm/min and photomultiplier voltage was set as 600 V, response: 2.0 s.
Table 1
Binding sites n and binding constants Ka for complexes 1a, 2b and 3.
RNA
TAR
Compd.
n
Ka (ꢂ 104 L/mol)
R2
1a
2b
3
1a
2b
0.85
1.50
0.98
0.29
0.49
1.6
4.0
2.6
0.9
2.2
0.946
0.977
0.999
0.982
0.997
RRE
GCCAGAUUUGAGCCUGGGAGCUCUCUGGC (PDB ID: 1QD3, a se-
quence used in oligonucleotide study) obtained from the Protein
4.2 package to study the RNA-binding properties. The copper
parameters, a vdW radii of 0.96 Å and a vdW well depth of 0.01
kcal/mol, used in the docking calculation were taken from ref. [50].
The coordinates of complexes 1a and 2b were taken from their
crystal structures as a CIF file and converted to the PDB format
using Mercury software 3.10. The receptor (TAR RNA) and the
ligands (dinuclear complexes) files were prepared using AutoDock
Tools. The heteroatoms including water molecules were deleted
and polar hydrogen atoms and Kollman charges were added to the
receptor molecule. All other bonds were allowed to be rotatable. In
the docking analysis, the binding site was assigned across all of the
minor and major grooves of the RNA molecule, which was enclosed
in a box with the number of grid points in x ꢂ y ꢂ z directions,
70 ꢂ 70 ꢂ 70 and a grid spacing of 0.375 Å. Initially, AutoGrid was
run to generate the grid map of various atoms of the ligands and
receptor. After the completion of the grid map, AutoDock was run
by using autodock parameters as follows: GA population size, 150;
maximum number of energy evaluations 2,500,000; and the
number of generations 27,000. A total of 10 runs were carried out. A
maximum of 30 conformers were considered for each molecule,
and the root-mean-square (RMS) cluster tolerance was set to 2.0 Å
in each run. All calculations were performed on an Intel Core i5
based machine running Windows as the operating system. For
each of the docking cases, the lowest energy docked conformation,
according to the Autodock scoring function, was selected as the
binding mode. Visualization of the docked pose was done by using
PyMOL (The PyMOL Molecular Graphics System, Version 2.0,
Schrödinger, LLC) molecular graphics program [51].
Fig. 4. Visualization of the docked poses of complexes 1a (A) and 2b (B) with TAR
RNA. The wireframe model of TAR with the dimetallic complexes (ball and stick) and
showing the interaction.
2.894 Å, leading to an increase of the electrostatic energy and thus
the stronger binding affinity of complex 2b to TAR RNA, which
coincides with the results obtained from the experimental studies
(Table 2). While the AutoDock calculation indicated 2b did not
directly interact with TAR RNA, we predicted a hydrogen bond
network might be involved via water molecule(s) [52–54]. It
should be pointed out that the modeling studies of the interaction
only allow a comparison between the two potential binding
ligands. However, the above spectroscopy experiments together
with the molecular docking studies can provide a comprehensive
understanding of the RNA binding activities, proving that the
dinuclear metallic complexes can strongly recover the intrinsic
fluorescence of ICR through the static binding mechanism with the
ratio of the dicopperII complex being close to 1:1 while the
dinickelII being 1.5:1, and the most possible binding site is in the
proximity of stem area. Thus, the reason that compounds 1a and 2b
bind preferably with TAR may be ascribed to the longer distance
between upper stem and lower stem in RRE.
With the purpose of gaining more information concerning the
favorably ligands of 2b vs. 1a towards TAR RNA, the binding sites
and the binding poses of the two dimetallic complexes in the
docking model with 1QD3 are shown in Fig. 4, from which it can be
seen that the binding poses of the two dimetallic complexes are
different. While complex 1a (O62) can interact with the NH7 (A20)
(1.868 Å), 2b forms four different bindings: OP2 (U38)-N3, 2.878 Å;
OP1 (G36)-O3, 2.678 Å, OP1 (G36)-O1, 2.936 Å; OP2 (G21)-O2,
Please cite this article in press as: X. Li, et al., Selective recognition of HIV RNA by dinuclear metallic ligands, Chin. Chem. Lett. (2018), https://