<title>Background</title> <p>The main aim of our molecular modelling investigations was to identify natural compounds for their ability to bind to the NF-kappaB p50 as a strategy to identify molecules exhibiting inhibitory activity on the molecular interactions of the transcription factor with its target DNA sequence. p50–p65 heterodimer is the predominant NF-kappaB complex in T-cells regulating HIV-1 infection and recent studies have shown that p50 unit of NF-kappaB is the one that mainly interacts with HIV-1 LTR [##REF##11160127##1##,##REF##10212261##2##]. The specific protein residues involved in DNA binding to the HIV-1 LTR NF-kappaB sites (sequence 5'-GGGACTTTCCC-3') have been identified [##REF##9450761##3##,##REF##7830764##4##]. Structurally different inhibitors of the NF-kappaB/DNA interactions with a rather low binding constant (in the range of 30 μM and 500 μM) are reported in the literature [##REF##10976530##5##, ####REF##15788163##6##, ##REF##15546742##7####15546742##7##]. Recently, some molecular modelling studies have predicted possible binding mode of the inhibitors molecules to the DNA binding region of subunit p50, starting from the crystallographic structure of the NF-kappaB homodimer [##REF##15788163##6##, ####REF##15546742##7##, ##REF##16759101##8##, ##REF##15124931##9####15124931##9##].</p> <p>In particular, Sharma et al. [##REF##16759101##8##] in an effort to rationalize the results obtained from EMSA studies on a set of aurintricarboxylic acid analogues, employed docking studies to explain the structure activity relationships observed within this class. To the best of our knowledge, nowadays the identification of new lead compounds for NF-kappaB inhibition through virtual screening of structures libraries is not yet reported in literature. In this paper, we present docking studies on a series of natural compounds previously identified within medicinal plant extracts by us, into NF-kappaB p50 protein target. After evaluation through electrophoretic mobility shift assays (EMSA), we obtained a fairly good agreement between experimental data and molecular modelling identification of bioactive and inactive compounds.</p>
<title>Methods</title> <title>Docking studies</title> <title>Ligands data and preparation</title> <p>The database of 27 natural structures used in our molecular docking studies, were derived from different medicinal plant extracts (Figure ##FIG##0##1##) as prepared in our laboratory. A dataset of 12 active compounds used as references molecules were collected from four publications [##REF##15788163##6##, ####REF##15546742##7##, ##REF##16759101##8##, ##REF##15124931##9####15124931##9##] reported by one laboratory (Figure ##FIG##1##2##). Ten of these inhibitors (<bold>1i-8i</bold>, <bold>11i </bold>and <bold>12i</bold>) were employed in starting docking studies (protocol 1) and in the Standard Similarity Scoring for subsequently docking simulations.</p> <p>Two inhibitory molecules (<bold>9i </bold>and <bold>10i</bold>) were used as test set in all docking simulations. The three-dimensional models of all the molecules under investigation were built by assembling fragments from the SYBYL 7.0 software package standard library [##UREF##0##10##]. Resulting geometries were optimized and molecular charges were assigned by a semi empirical molecular orbital calculations using the AM1 Hamiltonian [##UREF##1##11##] (module MOPAC implemented in SYBYL).</p> <title>Proteins data and preparation</title> <p>The three dimensional structure of the complex NF-kappaB-DNA [##REF##7830764##4##] was retrieved from the Protein Data Bank (PDB code: <ext-link ext-link-type="pdb" xlink:href="1NFK">1NFK</ext-link>). The cocrystallized DNA macromolecule was removed from the structure. p50 dimer and p50 monomers (chains A and B) were selected for the docking simulations and prepared using the graphical interface Maestro [##UREF##2##12##]. All water molecules were removed, the hydrogen atoms were added to the proteins and all atom force field (OPSL-2001) charges and atom types were assigned. Preparation and refinement were done running ProteinPrep job on the structure in a standard procedure. Minimizations were performed until the average root mean square deviation of non-hydrogen atoms reached 0.3 Å.</p> <title>Docking Simulations</title> <p>All molecules of plant extracts (<bold>1–27</bold>) and the known inhibitors (<bold>1i-12i</bold>) under study were docked in to the binding site of the receptor (PDB ID: <ext-link ext-link-type="pdb" xlink:href="1NFK">1NFK</ext-link>) using Glide (Grid-Based Ligand Docking With Energetics) software from Schrodinger [##REF##15027866##13##,##REF##15027865##14##]. Grids were prepared for each proteins with the exact same center and the size of the bounding box set on 30 Å. The coordinates of the enclosing box (x = -1,1958 Å; y = 9.0149 Å; z = 19,7598 Å) were defined starting from the set of active site residues involved in hydrogen bonds with the NF-kappaB recognition site of DNA (Arg54, Arg56, Tyr57, Cys59, Lys241, Gln306 and Thr143) and optimised including the double strands DNA helices volume by visual inspection. The Glide algorithm is based on a systematic search of positions, orientations, and conformations of the ligand in the receptor binding site using funnel type approach. The search begins with a rough positioning and scoring phase that significantly limits the search space and reduces the number of poses to be selected for minimization on the precomputed OPLS-2001 van der Waals and electrostatic grids for the protein. The 5–10 lowest-energy poses obtained from this stage are subjected to Monte Carlo simulations and the minimized poses accepted are then rescored using the GlideScore function, which is a more sophisticated version of ChemScore [##REF##9385547##15##]. This force field include additional terms accounting for solvation and repulsive interactions. In order to provide a better correlation between good poses and good scores, Glide Extra-Precision (XP) Mode was subsequently used on the conformations selected from Glide Standard Precision (SP) mode. The atom-pair superimposition of p50 chains A and B, prepared as described above, gave a minimum RMSD of 2,303 Å (heavy atoms). Considering the clear dependence of the docking accuracy of ligands on the protein structure, docking simulations were carried out with the same protocol on both A and B, considered as two slightly different conformations of the same structure.</p> <p>Unfortunately, complexes of NF-kappaB cocrystallized with inhibitors has not been solved. Therefore, a common self-docking procedure to evaluate the accuracy of the docking protocol adopted was not practicable. In order to overcame this situation, two structurally similar active compounds (<bold>9i </bold>and <bold>10i</bold>) were used as test set and docked into the DNA binding site of the protein. Moreover in following docking jobs, atom pair similarity (AP) scoring (Similscore) facility as implemented in Glide, was incorporated in GlideScore (G-score), based on the assumption that closely related chemical structure should share similar biological activity and physiochemical property [##UREF##3##16##]. Similscore can have a value between 0 and 1 as implemented in Glide. The adjusting G-score value is illustrated here below:</p> <p>1. if 0.0 ≤ SimilScore &lt; 0.3 → G-score = G-score+6.0</p> <p>2. if 0.3 ≤ SimilScore &lt; 0.7 → G-score = G-score+(0.7-Similscore)/(0.7-0.3)*6.0</p> <p>3. if 0.7 ≤ SimilScore &lt; 1.0 → G-score = G-score+0.0</p> <p>All inhibitors molecules, except for <bold>9i </bold>and <bold>10i</bold>, were used just as reference structures for AP similarity method.</p> <p>Based on the best final GlideScore ranking, the similarity docking procedure for subsequently docking simulations on p50 subunits was chosen.</p> <title>Preparation of nuclear extracts</title> <p>Nuclear extracts were prepared as described [##REF##12446679##18##]. Cell were washed twice with PBS and detached by trypsinization. After homogenization with Dounce B homogeneizer, nuclear proteins were obtained and protein concentration was determined using Bio-Rad protein assay. Nuclear extracts were brought to a concentration of 0.5 μg/μl for Electrophoretic Mobility Shift Assay (EMSA) experiments [##REF##12446679##18##].</p> <title>Electrophoretic Mobility Shift Assay (EMSA)</title> <p>EMSA was perfomed as previously described [##REF##12446679##18##, ####REF##12167479##19##, ##REF##17466942##20####17466942##20##]. Briefly, double-stranded synthetic oligodeoxynucleotides mimicking the NF-κB binding site present in the promoter of the IL-8 gene (IL-8 NF-κB, sense: 5'-AAT CGT GGA ATT TCC TCT-3') have been employed. Oligodeoxynucleotides were labeled with γ<sup>32</sup>-P-ATP using 10 Units of T4-polynucleotide-kinase (MBI Fermentas) in 500 mM Tris-HCl, pH 7.6, 100 mM MgCl<sub>2</sub>, 50 mM DTT, 1 mM spermidine, 1 mM EDTA in the presence of 50 μCi γ<sup>32</sup>-P-ATP) in a volume of 20 μl for 45 minutes at 37°C. Reaction was brought to 150 mM NaCl and 150 ng complementary oligodeoxynucleotide was added. Reaction temperature was increased to 100°C for 5 minutes and left diminishing to room temperature overnight. Nuclear extracts from IB3-1 cells or purified NF-κB p50 dimer protein (Promega) were used at the specified concentrations and poly(dI:dC) (1 mg per reaction) was also added to abolish nonspecific binding [##REF##18258920##21##]. After 5 min binding at room temperature, the samples were run at constant voltage (200 V) under low ionic strength conditions (0.25× TBE buffer: 22 mM Tris-borate, 0.4 mM EDTA) on 6% polyacrylamide gels. Gels were dried and subjected to standard autoradiographic procedures.</p> <title>Cell cultures and infection with <italic>Pseudomonas aeruginosa</italic></title> <p>IB3-1 cells have been obtained from LGC Promochem [##REF##17466942##20##]. Cells have been grown in LHC-8 basal medium (Biofluids), supplemented with 5% FBS in the absence of gentamycin [##REF##18258920##21##]. All culture flasks and plates have been coated with a solution containing 35 g/ml bovine collagen (Becton-Dickinson), 1 g/ml bovine serum albumin (Sigma) and 1 g/ml human fibronectin (Becton-Dickinson) as described. P. aeruginosa, PAO1 strain, was grown in trypticase soy broth (TSB) or agar (TSA) (Difco). Bacteria colonies from o/n cultures on TSA plates were grown in 20 ml TSB at 37°C. IB3-1 cells were infected with ranging doses of PAO1 at 37°C for 4 hours.</p> <title>Quantitation of transcripts of inflammatory genes</title> <p>Total RNA from IB3-1 cells was isolated using High Pure RNA Isolation Kit (Roche. Mannheim. Germany) [##REF##18258920##21##]. Total RNA (2.5 μg) was reverse-transcribed to cDNA using the High Capacity cDNA Archive Kit and random primers (Applied Biosystems) in a 100-μl reaction. The cDNA (2 μl) was then amplified for 50 PCR cycles using the Platinum<sup>® </sup>SYBR<sup>® </sup>Green qPCR SuperMix-UDG (Invitrogen) in an ABI Prism 5700 sequence detection system (Applied Biosystems). The real-time PCR reactions were performed in duplicates for both target and normalizer genes. Primer sequences for detection of IL-8 mRNA were GACCACACTGCGCCAACA (IL-8 forward) and GCTCTCTTCCATCAGAAAGTTACATAATTT(IL-8 reverse). Primer sets were purchased from Sigma-Genosys (The Woodlands. TX). Results were collected with Sequence Detection Software (version 1.3; Applied Biosystems). Relative quantification of gene expression was performed using the comparative threshold (C<sub>T</sub>) method as described by the manufacturer (Applied Biosystems User Bulletin 2). Changes in mRNA expression level were calculated following normalization to calibrator gene [##REF##12738678##22##]. The ratios obtained following normalization are expressed as -fold change over untreated samples.</p>
<title>Results and discussion</title> <title>Docking analysis</title> <p>The docking results for all the reference inhibitory compounds and the natural compounds under study are reported in Table ##TAB##0##1## (references compounds) and in Table ##TAB##1##2## (natural compounds). As shown in the 2.3 Å crystal structure, the DNA/p50 complex is formed by one DNA molecule and two p50 proteins each one consisting of two distinct domains connected by a10-residue linker. Both domains and the segments that connects them, form a sequence-specific DNA-binding surface by contributing 5 loops per subunit that fill the entire major groove of the DNA. The specific interactions that stabilized the NF-kappaB/DNA complex, occur over 10-bp forming the kB recognition site. Unlike many dimeric protein-DNA complexes, many residues of both subunits make specific base contacts in a non-contiguous cooperative network. The plasticity of centre region of the interface carry to the lack of symmetry exhibited by the interactions of Lys 241 from the linker segment, and Lys 272 and Arg 305 from the dimerization region with the symmetrical target site [##REF##7830764##4##]. In the subsequent experimental EMSA studies, a recombinant p50 protein that probably forms a monomer-dimer mixture in binding buffer solution will be used. On the base of structural and experimental assumptions as above mentioned, p50 dimer and monomer were employed as protein target in our molecular modelling investigation.</p> <p>In order to evaluate the impact of the introduction of the similarity penalty in the docking algorithm on the results, the positions of <bold>9i </bold>and <bold>10i </bold>used as test set in the final GlideScore ranking were compared in the two different procedures (Table ##TAB##1##2##).</p> <p>Known active compound <bold>10i </bold>was ranked at the top positions in both procedures, but only the introduction of the similarity parameter in the scoring function significantly increased the efficiency in <bold>9i </bold>ranking (Table ##TAB##1##2##). In fact the difference in glide-score among these two inhibitors with a similar inhibitory potency (500 μM) [##REF##15124931##9##] was very small (ΔGlideScore = 0.10). For each selected ligand the pose with best E-Model score (combination of energy grid score, GlideScore, and the internal strain of the ligand) was used for in-deph interaction analysis. Compound <bold>21 </bold>clearly showed highest score in respect to docked plant extracts (Table ##TAB##1##2##) outranking the known inhibitors at physiological pH in docking simulation to the dimer.</p> <p>Docked compounds <bold>1–27</bold>, <bold>9i </bold>and <bold>10i</bold>, occupied a region of the binding surface creates by the spatial relationship between the N-terminal domain of p50 subunit and the 10 residues long linker loop (Figure ##FIG##2##3##). Molecules <bold>21</bold>, <bold>9i </bold>and <bold>10i </bold>(Figure ##FIG##3##4##) were located in a small cleft surrounded by several polar amino acids (i.e. Tyr57, His109, His141, Tyr143 Lys144, Lys145, Ser208, Asp239, Lys241 and Ser208) and the highest score poses were superimposable with minimum RMSD of 1.36 Å for compounds <bold>9i </bold>and <bold>21</bold>. The RMSD was calculated by superimposing the following atoms pairs: heteroatoms involved in hydrogen bonding with the same residues of the protein (<bold>9i</bold>.O8 and <bold>21</bold>.O2'; <bold>9i</bold>.O7 and <bold>21</bold>.O1a) (Figure ##FIG##3##4B##) and the centroid of aromatic system of coumarin structure with the centroid of benzene ring of <bold>21</bold>. These compounds showed slightly different binding modes in p50 (chain A), p50 (chain B) and p50-p50 targets. Here we reported the highest score poses obtained from docking protocol including the similarity function. H-bond interactions between OH groups of coumarin structures (OH of benzene ring in <bold>21</bold>) and both NH of His141 and the carboxylic group of Asp239 showed to be important for ligands binding.</p> <p>Moreover, OH groups of coumarin moiety (carboxylate group in <bold>21</bold>) made an additional hydrogen bond with CO of the backbone and protonated NH<sub>3 </sub>group of Lys241 (Table ##TAB##2##3##). It is important to note that Lys241 could be involved in the stability of the DNA-binding conformation of the protein. In fact, as discussed above, this residue is situated in the flexible linker segment and interacts with Lys 272 and Arg 305 from the dimerization domain. Finally, carbonyl group of <bold>10i </bold>engage another H-bond with NH of the backbone of the Leu207. Compound <bold>21 </bold>showed the same binding mode of active ligands in the monomer configuration of the target, with the only difference of a stronger interaction of carboxylate group with Lys241 (Table ##TAB##2##3## and Figures ##FIG##4##5A##, ##FIG##5##6A##).</p> <p>Interestingly, the best pose of compound <bold>21 </bold>occupied a region formed by residues of both p50 units (chain A and chain B) of NF-kappaB dimer: Lys 145 and Thr143 of chain A and Tyr57, Lys144, Lys145, Glu60, Cys59, Thr143, Lys146 of chain B. In particular, the OH group of the ligand engages a hydrogen bond with the sidechain of Thr143 (chain B), and the carboxylate group forms a salt bridge stabilized by two hydrogen bonds with the side chain of Lys 145 (chain B). Moreover the phenyl structure of compound <bold>21 </bold>could be involved in a weak π-π stacking interaction with the aromatic moiety of Tyr57 (chain B) (centroid-centroid distance: 4.93 Å), a residue specific for kB DNA sequence 5'-GGGATTTCC-3', present in different cellular genes including HIV-LTR. Of course, further dynamics simulation on the protein-ligand complex should be necessary to validate this hypothesis. In addition, the amino group of Lys145 of the opposite p50 unit (chain A) could form an additional π-cation interaction with the aromatic group of <bold>21 </bold>(N<sup>+</sup>-H and centroid of benzene ring distance: 3.87 Å) (Figures ##FIG##4##5B##, ##FIG##5##6B##). These bridge structures are likely to reinforce the anchoring of this molecule to the DNA binding region of the dimer, and might account for the slight better G-score of <bold>21 </bold>in respect to the monomer configuration of the receptor. Moreover, all the residues of the protein involved in molecular interactions with molecule <bold>21</bold>, form hydrogen bonds also with DNA.</p> <p>All compounds with higher GlideScore and E-Model score clearly showed the ability to make a maximum number of hydrogen bonding, according with the result as previously reported on a flexible docking studies of known inhibitors <bold>9i </bold>and <bold>10i </bold>[##REF##15124931##9##], even if reported residues involved in binding interaction were different. The highest ranking poses of <bold>21</bold>, <bold>9i </bold>and <bold>10i </bold>formed 3–4 hydrogen bonding with the target protein, whereas molecules in medium positions in docking ranking not more than 2. According, structures not involved in hydrogen bonding were ranked in the last positions (Table ##TAB##1##2##). In particular, compound <bold>5 </bold>with a GlideScore &lt; 0 in similarity protocol lost the ability both to occupy the same positions of active ligands and to form hydrogen bonding with the protein (not shown). In house experimental data were in good agreement with the molecular modelling findings. In accordance with docking results, <bold>21 </bold>and <bold>5 </bold>showed to be active and inactive respectively in further EMSA experimental studies.</p> <title>Effects of compound 21 on NF-kappaB/DNA interactions</title> <p>The effects of compound <bold>21 </bold>on NF-kappaB interactions were first studied by electrophoretic mobility shift assay (EMSA) as described elsewhere [##REF##12446679##18##, ####REF##12167479##19##, ##REF##17466942##20##, ##REF##18258920##21####18258920##21##]. It is indeed well accepted that molecules binding NF-kappaB might retain inhibitory activity on molecular interaction between NF-kappaB and DNA [##REF##18258920##21##]. Accordingly, we performed EMSA in the presence of increasing amounts of compound <bold>21</bold>. In addition, compounds <bold>5 </bold>was used as possible negative control. This compound, indeed, is expected, from the docking analysis (Table ##TAB##1##2##), to be less active. In addition, extracts from <italic>Cupressus pyramidalis </italic>were also used, since this extract does not inhibit NF-kappaB/DNA interactions (data not shown). Finally, the known inhibitory compound <bold>9i </bold>was used as reference molecule. The results of the gel retardation analysis are shown in Figure ##FIG##6##7## and clearly demonstrate that compound <bold>21 </bold>inhibit the molecular interactions between nuclear factors (Figure ##FIG##6##7B##) or isolated NF-kappaB p50 (Figure ##FIG##6##7C## and ##FIG##6##7D##) and a target double stranded oligonucleotide mimicking the NF-kappaB binding sites. This effect was similar to that exhibited by the reference compound <bold>9i</bold>. Interestingly, compound <bold>5 </bold>and extracts from <italic>C. pyramidalis </italic>were found to be inactive (Figure ##FIG##6##7A##), fully in agreement with the docking data summarized in Table ##TAB##1##2##.</p> <title>Biological effects of compound 21: inhibition of <italic>Pseudomonas aeruginosa </italic>mediated increase of IL-8 mRNA</title> <p>Several experimental model system are available for biological validation of molecules inhibiting NF-kappaB function. In a recent paper we report that decoy oligonucleotides targeting NF-kappaB are powerful inhibitors of <italic>Pseudomoas aeruginosa </italic>mediated induction of IL-8 in cystic fibrosis IB3-1 cells [##REF##18258920##21##]. Besides the importance of these data for the theoretical point of view, our results are of great interest for the practical point of view, suggesting this treatment as a possible strategy for the therapy of inflammation associated with cystic fibrosis.</p> <p>When the effects of <italic>P. aeruginosa </italic>infection on the expression of pro-inflammatory genes of IB3-1 cells infected for 4 hours are analysed, the results shown in Figure ##FIG##7##8A## are obtained. In this preliminary experiment the content of RNAs coding for several pro-inflammatory proteins was analysed by RT-PCR. The results obtained indicate that IL-8 mRNA sequences sharply increases following PAO1 infection by 40 folds (range 12–92) in respect to basal level of uninfected cells, assumed to be 1. In addition to IL-8 mRNA, other genes induced by PAO1 are GRO-γ (15 folds, range 9–36), IL-6 (18 folds, range 5–44), IL-1β (5 folds, range 2–8), ICAM-1 (4 folds, range 2–8). On the contrary, very low increase of accumulation of IP-10, RANTES, MIP-1α, TNF-α, IFN-β, TGF-β and IL-10 mRNA was observed under these experimental conditions.</p> <p>Since NF-kappaB is one of the most important transcription factors regulating the expression of IL-8 gene [23] and the data reported in Figure ##FIG##6##7## demonstrate that compound <bold>21 </bold>inhibit NF-kappaB/DNA interactions, we tested the activity of this compound in inhibiting the expression of IL-8 gene in IB3-1 cells infected with PAO1. Cells were exposed to different concentrations of compound <bold>21 </bold>and then infected with 150 cfu/cell of PAO1. After 4 hours, cells were harvested, RNA extracted and quantitative RT-PCR analysis performed. The results obtained (Figure ##FIG##7##8B##) demonstrate that compound <bold>21 </bold>is a strong inhibitor of PAO-1 induced accumulation of IL-8 mRNA.</p>
<title>Results and discussion</title> <title>Docking analysis</title> <p>The docking results for all the reference inhibitory compounds and the natural compounds under study are reported in Table ##TAB##0##1## (references compounds) and in Table ##TAB##1##2## (natural compounds). As shown in the 2.3 Å crystal structure, the DNA/p50 complex is formed by one DNA molecule and two p50 proteins each one consisting of two distinct domains connected by a10-residue linker. Both domains and the segments that connects them, form a sequence-specific DNA-binding surface by contributing 5 loops per subunit that fill the entire major groove of the DNA. The specific interactions that stabilized the NF-kappaB/DNA complex, occur over 10-bp forming the kB recognition site. Unlike many dimeric protein-DNA complexes, many residues of both subunits make specific base contacts in a non-contiguous cooperative network. The plasticity of centre region of the interface carry to the lack of symmetry exhibited by the interactions of Lys 241 from the linker segment, and Lys 272 and Arg 305 from the dimerization region with the symmetrical target site [##REF##7830764##4##]. In the subsequent experimental EMSA studies, a recombinant p50 protein that probably forms a monomer-dimer mixture in binding buffer solution will be used. On the base of structural and experimental assumptions as above mentioned, p50 dimer and monomer were employed as protein target in our molecular modelling investigation.</p> <p>In order to evaluate the impact of the introduction of the similarity penalty in the docking algorithm on the results, the positions of <bold>9i </bold>and <bold>10i </bold>used as test set in the final GlideScore ranking were compared in the two different procedures (Table ##TAB##1##2##).</p> <p>Known active compound <bold>10i </bold>was ranked at the top positions in both procedures, but only the introduction of the similarity parameter in the scoring function significantly increased the efficiency in <bold>9i </bold>ranking (Table ##TAB##1##2##). In fact the difference in glide-score among these two inhibitors with a similar inhibitory potency (500 μM) [##REF##15124931##9##] was very small (ΔGlideScore = 0.10). For each selected ligand the pose with best E-Model score (combination of energy grid score, GlideScore, and the internal strain of the ligand) was used for in-deph interaction analysis. Compound <bold>21 </bold>clearly showed highest score in respect to docked plant extracts (Table ##TAB##1##2##) outranking the known inhibitors at physiological pH in docking simulation to the dimer.</p> <p>Docked compounds <bold>1–27</bold>, <bold>9i </bold>and <bold>10i</bold>, occupied a region of the binding surface creates by the spatial relationship between the N-terminal domain of p50 subunit and the 10 residues long linker loop (Figure ##FIG##2##3##). Molecules <bold>21</bold>, <bold>9i </bold>and <bold>10i </bold>(Figure ##FIG##3##4##) were located in a small cleft surrounded by several polar amino acids (i.e. Tyr57, His109, His141, Tyr143 Lys144, Lys145, Ser208, Asp239, Lys241 and Ser208) and the highest score poses were superimposable with minimum RMSD of 1.36 Å for compounds <bold>9i </bold>and <bold>21</bold>. The RMSD was calculated by superimposing the following atoms pairs: heteroatoms involved in hydrogen bonding with the same residues of the protein (<bold>9i</bold>.O8 and <bold>21</bold>.O2'; <bold>9i</bold>.O7 and <bold>21</bold>.O1a) (Figure ##FIG##3##4B##) and the centroid of aromatic system of coumarin structure with the centroid of benzene ring of <bold>21</bold>. These compounds showed slightly different binding modes in p50 (chain A), p50 (chain B) and p50-p50 targets. Here we reported the highest score poses obtained from docking protocol including the similarity function. H-bond interactions between OH groups of coumarin structures (OH of benzene ring in <bold>21</bold>) and both NH of His141 and the carboxylic group of Asp239 showed to be important for ligands binding.</p> <p>Moreover, OH groups of coumarin moiety (carboxylate group in <bold>21</bold>) made an additional hydrogen bond with CO of the backbone and protonated NH<sub>3 </sub>group of Lys241 (Table ##TAB##2##3##). It is important to note that Lys241 could be involved in the stability of the DNA-binding conformation of the protein. In fact, as discussed above, this residue is situated in the flexible linker segment and interacts with Lys 272 and Arg 305 from the dimerization domain. Finally, carbonyl group of <bold>10i </bold>engage another H-bond with NH of the backbone of the Leu207. Compound <bold>21 </bold>showed the same binding mode of active ligands in the monomer configuration of the target, with the only difference of a stronger interaction of carboxylate group with Lys241 (Table ##TAB##2##3## and Figures ##FIG##4##5A##, ##FIG##5##6A##).</p> <p>Interestingly, the best pose of compound <bold>21 </bold>occupied a region formed by residues of both p50 units (chain A and chain B) of NF-kappaB dimer: Lys 145 and Thr143 of chain A and Tyr57, Lys144, Lys145, Glu60, Cys59, Thr143, Lys146 of chain B. In particular, the OH group of the ligand engages a hydrogen bond with the sidechain of Thr143 (chain B), and the carboxylate group forms a salt bridge stabilized by two hydrogen bonds with the side chain of Lys 145 (chain B). Moreover the phenyl structure of compound <bold>21 </bold>could be involved in a weak π-π stacking interaction with the aromatic moiety of Tyr57 (chain B) (centroid-centroid distance: 4.93 Å), a residue specific for kB DNA sequence 5'-GGGATTTCC-3', present in different cellular genes including HIV-LTR. Of course, further dynamics simulation on the protein-ligand complex should be necessary to validate this hypothesis. In addition, the amino group of Lys145 of the opposite p50 unit (chain A) could form an additional π-cation interaction with the aromatic group of <bold>21 </bold>(N<sup>+</sup>-H and centroid of benzene ring distance: 3.87 Å) (Figures ##FIG##4##5B##, ##FIG##5##6B##). These bridge structures are likely to reinforce the anchoring of this molecule to the DNA binding region of the dimer, and might account for the slight better G-score of <bold>21 </bold>in respect to the monomer configuration of the receptor. Moreover, all the residues of the protein involved in molecular interactions with molecule <bold>21</bold>, form hydrogen bonds also with DNA.</p> <p>All compounds with higher GlideScore and E-Model score clearly showed the ability to make a maximum number of hydrogen bonding, according with the result as previously reported on a flexible docking studies of known inhibitors <bold>9i </bold>and <bold>10i </bold>[##REF##15124931##9##], even if reported residues involved in binding interaction were different. The highest ranking poses of <bold>21</bold>, <bold>9i </bold>and <bold>10i </bold>formed 3–4 hydrogen bonding with the target protein, whereas molecules in medium positions in docking ranking not more than 2. According, structures not involved in hydrogen bonding were ranked in the last positions (Table ##TAB##1##2##). In particular, compound <bold>5 </bold>with a GlideScore &lt; 0 in similarity protocol lost the ability both to occupy the same positions of active ligands and to form hydrogen bonding with the protein (not shown). In house experimental data were in good agreement with the molecular modelling findings. In accordance with docking results, <bold>21 </bold>and <bold>5 </bold>showed to be active and inactive respectively in further EMSA experimental studies.</p> <title>Effects of compound 21 on NF-kappaB/DNA interactions</title> <p>The effects of compound <bold>21 </bold>on NF-kappaB interactions were first studied by electrophoretic mobility shift assay (EMSA) as described elsewhere [##REF##12446679##18##, ####REF##12167479##19##, ##REF##17466942##20##, ##REF##18258920##21####18258920##21##]. It is indeed well accepted that molecules binding NF-kappaB might retain inhibitory activity on molecular interaction between NF-kappaB and DNA [##REF##18258920##21##]. Accordingly, we performed EMSA in the presence of increasing amounts of compound <bold>21</bold>. In addition, compounds <bold>5 </bold>was used as possible negative control. This compound, indeed, is expected, from the docking analysis (Table ##TAB##1##2##), to be less active. In addition, extracts from <italic>Cupressus pyramidalis </italic>were also used, since this extract does not inhibit NF-kappaB/DNA interactions (data not shown). Finally, the known inhibitory compound <bold>9i </bold>was used as reference molecule. The results of the gel retardation analysis are shown in Figure ##FIG##6##7## and clearly demonstrate that compound <bold>21 </bold>inhibit the molecular interactions between nuclear factors (Figure ##FIG##6##7B##) or isolated NF-kappaB p50 (Figure ##FIG##6##7C## and ##FIG##6##7D##) and a target double stranded oligonucleotide mimicking the NF-kappaB binding sites. This effect was similar to that exhibited by the reference compound <bold>9i</bold>. Interestingly, compound <bold>5 </bold>and extracts from <italic>C. pyramidalis </italic>were found to be inactive (Figure ##FIG##6##7A##), fully in agreement with the docking data summarized in Table ##TAB##1##2##.</p> <title>Biological effects of compound 21: inhibition of <italic>Pseudomonas aeruginosa </italic>mediated increase of IL-8 mRNA</title> <p>Several experimental model system are available for biological validation of molecules inhibiting NF-kappaB function. In a recent paper we report that decoy oligonucleotides targeting NF-kappaB are powerful inhibitors of <italic>Pseudomoas aeruginosa </italic>mediated induction of IL-8 in cystic fibrosis IB3-1 cells [##REF##18258920##21##]. Besides the importance of these data for the theoretical point of view, our results are of great interest for the practical point of view, suggesting this treatment as a possible strategy for the therapy of inflammation associated with cystic fibrosis.</p> <p>When the effects of <italic>P. aeruginosa </italic>infection on the expression of pro-inflammatory genes of IB3-1 cells infected for 4 hours are analysed, the results shown in Figure ##FIG##7##8A## are obtained. In this preliminary experiment the content of RNAs coding for several pro-inflammatory proteins was analysed by RT-PCR. The results obtained indicate that IL-8 mRNA sequences sharply increases following PAO1 infection by 40 folds (range 12–92) in respect to basal level of uninfected cells, assumed to be 1. In addition to IL-8 mRNA, other genes induced by PAO1 are GRO-γ (15 folds, range 9–36), IL-6 (18 folds, range 5–44), IL-1β (5 folds, range 2–8), ICAM-1 (4 folds, range 2–8). On the contrary, very low increase of accumulation of IP-10, RANTES, MIP-1α, TNF-α, IFN-β, TGF-β and IL-10 mRNA was observed under these experimental conditions.</p> <p>Since NF-kappaB is one of the most important transcription factors regulating the expression of IL-8 gene [23] and the data reported in Figure ##FIG##6##7## demonstrate that compound <bold>21 </bold>inhibit NF-kappaB/DNA interactions, we tested the activity of this compound in inhibiting the expression of IL-8 gene in IB3-1 cells infected with PAO1. Cells were exposed to different concentrations of compound <bold>21 </bold>and then infected with 150 cfu/cell of PAO1. After 4 hours, cells were harvested, RNA extracted and quantitative RT-PCR analysis performed. The results obtained (Figure ##FIG##7##8B##) demonstrate that compound <bold>21 </bold>is a strong inhibitor of PAO-1 induced accumulation of IL-8 mRNA.</p>
<title>Conclusion</title> <p>In the present work, we carried out docking studies on the dataset of 27 molecules found in different plant extracts to NF-kappaB-p50, with the purpose of developing a docking protocol fit for the target under study, eventually applicable for more time-consuming virtual screening of larger database of compounds.</p> <p>Usually, docking to protein structures that do not have a ligand present, as in the case of NF-kappaB, dramatically reduces the expected performance of structure-based methods. Therefore, the use of NF-kappaB as a target for virtual docking of natural compounds is not feasible. To overcome such a limitation, we proposed to enhance the simple docking procedure by means of a sort of combined target- and ligand-based drug design approach. Advantages of this combination strategy, based on a similarity parameter for the identification of weak binding chemical entities, are illustrated in this work with the discovery of a new lead compound for NF-kappaB. In this respect, this paper represents the first example of successfully individuation of a potential lead compound through molecular docking simulations on a NF-kappaB target. At the same time, information derived from this structure and its different binding modes, could carry through further lead optimization to more potent NF-kappaB inhibitors.</p> <p>In order to validate the approach, biochemical analyses based on EMSA were performed on compound <bold>21</bold>; the results obtained sustain the concept that the docking performance is predictive of a biochemical activity (Figure ##FIG##6##7##).</p> <p>Our results are of interest also from the practical point of view. The transcription factor NF-kappaB is indeed a very interesting target molecules in the design on anti-tumor, anti-inflammatory, pro-apoptotic drugs.</p> <p>In order to validate this last hypothesis, we have employed human cystic fibrosis IB3-1 tracheal epithelial cells. We have elsewhere reported that these cells activate, upon exposure to the bacterium <italic>Pseudomonas aeruginosa </italic>(the PAO-1 strain), the expression of several pro-inflammatory genes, including those coding interleukin-6 (IL-6) and interleukin-8 (IL-8). As supported by several groups, the expression of IL-8 is dependent from NF-kappaB activation. Accordingly, decoy molecules targeting NF-kappaB are strong inhibitors of the IL-8 expression. Therefore, PAO-1 infected IB3-1 cells are a very interesting model system to screen for IL-8 inhibitors. The results of our experiments, in agreement with both docking and EMSA data, demonstrate that compound <bold>21 </bold>is a strong inhibitor of IL-8 and should be considered of interest for modulation of the expression of this gene.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>The transcription factor NF-kappaB is a very interesting target molecule for the design on anti-tumor, anti-inflammatory and pro-apoptotic drugs. However, the application of the widely-used molecular docking computational method for the virtual screening of chemical libraries on NF-kappaB is not yet reported in literature. Docking studies on a dataset of 27 molecules from extracts of two different medicinal plants to NF-kappaB-p50 were performed with the purpose of developing a docking protocol fit for the target under study.</p> <title>Results</title> <p>We enhanced the simple docking procedure by means of a sort of combined target- and ligand-based drug design approach. Advantages of this combination strategy, based on a similarity parameter for the identification of weak binding chemical entities, are illustrated in this work with the discovery of a new lead compound for NF-kappaB. Further biochemical analyses based on EMSA were performed and biological effects were tested on the compound exhibiting the best docking score. All experimental analysis were in fairly good agreement with molecular modeling findings.</p> <title>Conclusion</title> <p>The results obtained sustain the concept that the docking performance is predictive of a biochemical activity. In this respect, this paper represents the first example of successfully individuation through molecular docking simulations of a promising lead compound for the inhibition of NF-kappaB-p50 biological activity and modulation of the expression of the NF-kB regulated IL8 gene.</p>
<title>Authors' contributions</title> <p>LP carried out all the bioinformatic procedures and the docking experiments. EF participated to the EMSA assays; MB purified the nuclear factors for EMSA analysis; I.M. performed semi-quantitative RT-PCR analysis; VB performed the treatment of IB3-1 cells with selected compounds; MCD performed infection IB3-1 cells with P. aeruginosa; EN performed quantitative RT-PCR analysis of IB-8 mRNA; GC was the responsible of the conception, design, analysis and interpretation of the data on cystic fibrosis cell lines; RG was the responsible of the coordination of the project and of the drafting of the manuscript. All authors read and approved the final manuscript.</p>
<title>Acknowledgements</title> <p>RG is granted by Fondazione CARIPARO, AIRC, Telethon and by MUR COFIN-2005. VB is fellow of the "Fondazione Cariverona", FQ and EN are fellows of the "Azienda Ospedaliera di Verona". This work was supported also by grants from the Italian Cystic Fibrosis Research Foundation (grants # 15/2004 and # 13/2007 to RG and GC) and Fondazione Cariverona – Bando 2005 – Malattie rare e della povertà (to GC).</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Structures of compounds found in <italic>Cupressus pyramidalis and Aegle marmelos </italic>extracts and used for docking simulations.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>NF-kappaB/DNA binding inhibitors used for atom-pair similarity scoring in docking.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Stereoview of compounds <bold>1–27</bold>, <bold>9i </bold>and <bold>10i </bold>docked in to DNA binding region of the NF-kappaB p50 monomer chain A. The macromolecule is highlighted in green.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>Superimposition of the docked poses of inhibitors <bold>9i</bold>, <bold>10i </bold>and compound <bold>21</bold>. A. the DNA binding site of NF-kappaB p50 (monomer chain A) is highlighted in green; B. the ligand atoms involved in hydrogen bonding are labeled. Compounds <bold>9i</bold>, (shown in purple), <bold>10i </bold>(shown in blue) and <bold>21 </bold>are illustrated in stick representations.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p>Binding modes of <bold>21 </bold>docked in to the active site of NF-kappaB p50: A. monomer and B. homodimer (chain A, shown in green). The residues involved in the interaction with the ligand are shown; the hydrogen bonding and the relative distances are indicated in purple.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>Poses of docked compound <bold>21 </bold>in to the DNA binding region of NF-kappaB p50. A. monomer and B. homodimer (chain A, shown in grey). The inhibitory activity of <bold>21 </bold>may be due to its ability to form a stable complex with the active conformation of the dimer and/or blocking the interaction of DNA with the monomer filling the protein binding site. The DNA was obtained from the crystal structure of the homodimer NF-kappaB (pdb code:<ext-link ext-link-type="pdb" xlink:href="1NFK">1NFK</ext-link>). The surface of the protein is represented in wire frame, the ligand and DNA are highlighted in VDW and stick representation respectively.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p>EMSA analysis. NF-kappaB p50 has been incubated for 15 min in the presence of increasing amounts of compounds <bold>5 </bold>(right side of panel A), extracts from <italic>C. pyramidalis </italic>(left side of panel A), compound <bold>21 </bold>(B-D) and reference compound <bold>9i </bold>(right side of panel C). After the incubation of compounds (or <italic>C. pyramidalis </italic>extracts) to NF-kappaB p50, a further 15 min incubation step was carried on with a <sup>32</sup>P-labelled double stranded oligonucleotide carrying the NF-kappaB binding sites. <italic>C. pyramidalis </italic>extracts were employed at a concentration ranging from 15 to 1000 μg/reaction as mentioned. Comparison between compound <bold>21</bold>, compound <bold>5</bold>, and compound <bold>9i </bold>was performed at 0.01, 0.1, 1 an 10 mM concentrations. Concentration-dependency of the effects displayed by compound <bold>21 </bold>was further analysed with intermediate dosages in the experiment depicted in panel D. The NF-kappaB/DNA complex were identified after polyacrylamide gel electrophoresis. Arrows indicate NF-kappaB/DNA complexes. Asterisks indicate free target DNA.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p>A. Effects of PAO-1 infection of cystic fibrosis IB3-1 cells on the expression of the indicated mRNA. Cells were infected with PAO-1 for 4 hours and then the mRNA analysed by RT-PCR. For RT-PCR analysis the PCR primers have been described in Bezzerri et al. [##REF##18258920##21##] B. Effects of compound 21 on the accumulation of IL-8 mRNA. Cells were incubated for 24 hours in the presence of the indicated concentrations of compound 21, then infected with 150 cfu/cell of PAO-1. Accumulation of IL-8 mRNA was determined by RT-PCR analysis [##REF##18258920##21##].</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Ranking of the poses of references inhibitory molecules (1i-8i and 11i-12i) in the target NF-kappaB p50 both as dimer (p50-p50) and as monomers (p50 A and B).</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center" colspan="6"><bold>Docking protocol 1</bold></td></tr></thead><tbody><tr><td align="center" colspan="2"><bold><italic>p50-p50</italic></bold></td><td align="center" colspan="2"><bold><italic>p50 (A)</italic></bold></td><td align="center" colspan="2"><bold><italic>p50 (B)</italic></bold></td></tr><tr><td colspan="2"><hr/></td><td colspan="2"><hr/></td><td colspan="2"><hr/></td></tr><tr><td align="center"><bold>Cpd.</bold></td><td align="center"><bold>G-Score</bold></td><td align="center"><bold>Cpd.</bold></td><td align="center"><bold>G-Score</bold></td><td align="center"><bold>Cpd</bold></td><td align="center"><bold>G-Score</bold></td></tr><tr><td colspan="6"><hr/></td></tr><tr><td align="center"><bold>8i</bold></td><td align="center">-5.88</td><td align="center"><bold>8i</bold></td><td align="center">-6.06</td><td align="center"><bold>1i</bold></td><td align="center">-6.08</td></tr><tr><td align="center"><bold>7i</bold></td><td align="center">-5.88</td><td align="center"><bold>7i</bold></td><td align="center">-5.83</td><td align="center"><bold>3i</bold></td><td align="center">-5.35</td></tr><tr><td align="center"><bold>2i</bold></td><td align="center">-5.84</td><td align="center"><bold>11i</bold></td><td align="center">-5.20</td><td align="center"><bold>7i</bold></td><td align="center">-5.06</td></tr><tr><td align="center"><bold>3i</bold></td><td align="center">-5.78</td><td align="center"><bold>1i</bold></td><td align="center">-5.07</td><td align="center"><bold>2i</bold></td><td align="center">-4.58</td></tr><tr><td align="center"><bold>11i</bold></td><td align="center">-3.87</td><td align="center"><bold>2i</bold></td><td align="center">-2.79</td><td align="center"><bold>11i</bold></td><td align="center">-3.59</td></tr><tr><td align="center"><bold>1i</bold></td><td align="center">-3.57</td><td align="center"><bold>3i</bold></td><td align="center">-2.77</td><td align="center"><bold>4i</bold></td><td align="center">-1.82</td></tr><tr><td align="center"><bold>4i</bold></td><td align="center">-0.76</td><td align="center"><bold>4i</bold></td><td align="center">&gt; 0</td><td align="center"><bold>8i</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>12i</bold></td><td align="center">&gt; 0</td><td align="center"><bold>12i</bold></td><td align="center">&gt; 0</td><td align="center"><bold>12i</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>5i</bold></td><td align="center">-</td><td align="center"><bold>5i</bold></td><td align="center">&gt; 0</td><td align="center"><bold>5i</bold></td><td align="center">-</td></tr><tr><td align="center"><bold>6i</bold></td><td align="center">-</td><td align="center"><bold>6i</bold></td><td align="center">&gt; 0</td><td align="center"><bold>6i</bold></td><td align="center">-</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Ranking of the poses of natural compounds and test set inhibitors (9i and 10i) in the target NF-kappaB p50 both as dimer (p50-p50) and as monomers (p50 A and B). In the docking protocol 1 the similarity scoring algorithm is not used.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center" colspan="2"><bold>Docking protocol 1</bold></td><td align="center" colspan="6"><bold>Docking protocol 2</bold></td></tr><tr><td colspan="2"><hr/></td><td colspan="6"><hr/></td></tr><tr><td align="center" colspan="2"><bold><italic>p50-p50</italic></bold></td><td align="center" colspan="2"><bold><italic>p50-p50</italic></bold></td><td align="center" colspan="2"><bold><italic>p50 (A)</italic></bold></td><td align="center" colspan="2"><bold><italic>p50 (B)</italic></bold></td></tr><tr><td colspan="2"><hr/></td><td colspan="2"><hr/></td><td colspan="2"><hr/></td><td colspan="2"><hr/></td></tr><tr><td align="center"><bold>Cpd.</bold></td><td align="center"><bold>G-score</bold></td><td/><td align="center"><bold>G-Score</bold></td><td align="center"><bold>Cpd.</bold></td><td align="center"><bold>G-Score</bold></td><td align="center"><bold>Cpd.</bold></td><td align="center"><bold>G-Score</bold></td></tr></thead><tbody><tr><td align="center"><bold><italic>10i</italic></bold></td><td align="center"><italic>-5.61</italic></td><td align="center"><bold><italic>21</italic></bold></td><td align="center"><italic>-2.34</italic></td><td align="center"><bold><italic>9i</italic></bold></td><td align="center"><italic>-3.81</italic></td><td align="center"><bold><italic>9i</italic></bold></td><td align="center"><italic>-2.13</italic></td></tr><tr><td align="center"><bold>15</bold></td><td align="center">-5.50</td><td align="center"><bold><italic>9i</italic></bold></td><td align="center"><italic>-2.29</italic></td><td align="center"><bold><italic>10i</italic></bold></td><td align="center"><italic>-2.27</italic></td><td align="center"><bold><italic>10i</italic></bold></td><td align="center"><italic>-1.43</italic></td></tr><tr><td align="center"><bold>25</bold></td><td align="center">-5.24</td><td align="center"><bold><italic>10i</italic></bold></td><td align="center"><italic>-2.19</italic></td><td align="center"><bold><italic>21</italic></bold></td><td align="center"><italic>-1.40</italic></td><td align="center"><bold><italic>21</italic></bold></td><td align="center"><italic>-0.32</italic></td></tr><tr><td align="center"><bold>23</bold></td><td align="center">-4.94</td><td align="center"><bold>18</bold></td><td align="center">-0.50</td><td align="center"><bold>18</bold></td><td align="center">-0.07</td><td align="center"><bold>18</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold><italic>21</italic></bold></td><td align="center"><italic>-4.91</italic></td><td align="center"><bold>23</bold></td><td align="center">-0.05</td><td align="center"><bold>10</bold></td><td align="center">&gt; 0</td><td align="center"><bold>15</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>24</bold></td><td align="center">-4.74</td><td align="center"><bold>20</bold></td><td align="center">&gt; 0</td><td align="center"><bold>23</bold></td><td align="center">&gt; 0</td><td align="center"><bold>23</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold><italic>9i</italic></bold></td><td align="center"><italic>-4.44</italic></td><td align="center"><bold>15</bold></td><td align="center">&gt; 0</td><td align="center"><bold>7</bold></td><td align="center">&gt; 0</td><td align="center"><bold>24</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>10</bold></td><td align="center">-4.40</td><td align="center"><bold>25</bold></td><td align="center">&gt; 0</td><td align="center"><bold>8</bold></td><td align="center">&gt; 0</td><td align="center"><bold>10</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>20</bold></td><td align="center">-4.37</td><td align="center"><bold>24</bold></td><td align="center">&gt; 0</td><td align="center"><bold>20</bold></td><td align="center">&gt; 0</td><td align="center"><bold>23</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>7</bold></td><td align="center">-4.28</td><td align="center"><bold>10</bold></td><td align="center">&gt; 0</td><td align="center"><bold>15</bold></td><td align="center">&gt; 0</td><td align="center"><bold>20</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>22</bold></td><td align="center">-4.13</td><td align="center"><bold>16</bold></td><td align="center">&gt; 0</td><td align="center"><bold>25</bold></td><td align="center">&gt; 0</td><td align="center"><bold>16</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>16</bold></td><td align="center">-3.48</td><td align="center"><bold>7</bold></td><td align="center">&gt; 0</td><td align="center"><bold>24</bold></td><td align="center">&gt; 0</td><td align="center"><bold>7</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>19</bold></td><td align="center">-3.71</td><td align="center"><bold>22</bold></td><td align="center">&gt; 0</td><td align="center"><bold>13</bold></td><td align="center">&gt; 0</td><td align="center"><bold>8</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>2</bold></td><td align="center">-3.39</td><td align="center"><bold>19</bold></td><td align="center">&gt; 0</td><td align="center"><bold>12</bold></td><td align="center">&gt; 0</td><td align="center"><bold>2</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>27</bold></td><td align="center">-3.09</td><td align="center"><bold>4</bold></td><td align="center">&gt; 0</td><td align="center"><bold>16</bold></td><td align="center">&gt; 0</td><td align="center"><bold>22</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>5</bold></td><td align="center">-3.05</td><td align="center"><bold>2</bold></td><td align="center">&gt; 0</td><td align="center"><bold>6</bold></td><td align="center">&gt; 0</td><td align="center"><bold>4</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>26</bold></td><td align="center">-2.89</td><td align="center"><bold>5</bold></td><td align="center">&gt; 0</td><td align="center"><bold>9</bold></td><td align="center">&gt; 0</td><td align="center"><bold>1</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>18</bold></td><td align="center">-2.69</td><td align="center"><bold>27</bold></td><td align="center">&gt; 0</td><td align="center"><bold>11</bold></td><td align="center">&gt; 0</td><td align="center"><bold>19</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>14</bold></td><td align="center">-2.47</td><td align="center"><bold>26</bold></td><td align="center">&gt; 0</td><td align="center"><bold>2</bold></td><td align="center">&gt; 0</td><td align="center"><bold>14</bold></td><td align="center">&gt; 0</td></tr><tr><td align="center"><bold>4</bold></td><td align="center">-2.00</td><td align="center"><bold>14</bold></td><td align="center">&gt; 0</td><td align="center"><bold>3</bold></td><td align="center">&gt; 0</td><td align="center"><bold>5</bold></td><td align="center">&gt; 0</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>intramolecular hydrogen bonds of the docked poses of 9i, 10i and 21 with the involved residues of the DNA binding region of NF-kappaB (see Figure 5 for ligand atom labels). The interatomic distances in Angstroms are shown.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Residue interaction</td><td align="center">Ligand atom</td><td align="center">Distance (Å)</td></tr></thead><tbody><tr><td align="center"><bold>(His141) </bold>-Nε-H::O</td><td align="center">O5 <bold>(9i)</bold></td><td align="center">2.12</td></tr><tr><td/><td align="center">O8 <bold>(10i)</bold></td><td align="center">1.88</td></tr><tr><td/><td align="center">O2' <bold>(21)</bold></td><td align="center">1.90</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="center"><bold>(Asp239) </bold>-COO::H</td><td align="center">H5 <bold>(9i)</bold></td><td align="center">1.91</td></tr><tr><td/><td align="center">H8 <bold>(10i)</bold></td><td align="center">2.09</td></tr><tr><td/><td align="center">H2' <bold>(21)</bold></td><td align="center">2.00</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="center"><bold>(Lys241) </bold>-CO::H</td><td align="center">H7 <bold>(10i)</bold></td><td align="center">1.96</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="center"><bold>(Lys241) </bold>-N<sup>+</sup>-H::O</td><td align="center">O7 <bold>(9i)</bold></td><td align="center">2.60</td></tr><tr><td/><td align="center">O1a <bold>(21)</bold></td><td align="center">2.17</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="center"><bold>(Leu207) </bold>-N-H::O</td><td align="center">O1 <bold>(10i)</bold></td><td align="center">2.18</td></tr></tbody></table></table-wrap>
<graphic xlink:href="1472-6807-8-38-1"/> <graphic xlink:href="1472-6807-8-38-2"/> <graphic xlink:href="1472-6807-8-38-3"/> <graphic xlink:href="1472-6807-8-38-4"/> <graphic xlink:href="1472-6807-8-38-5"/> <graphic xlink:href="1472-6807-8-38-6"/> <graphic xlink:href="1472-6807-8-38-7"/> <graphic xlink:href="1472-6807-8-38-8"/>
[{"article-title": ["TRIPOS"], "source": ["S Associates, St Louis, MO"], "year": ["1993"], "volume": ["version 7.0"]}, {"surname": ["Dewar", "Zoebish", "Healy"], "given-names": ["M", "G", "E"], "article-title": ["AM1: A new general purpose quantum mechanical molecular model"], "source": ["J Am Chem Soc"], "year": ["1985"], "volume": ["107"], "fpage": ["3902"], "lpage": ["3909"], "pub-id": ["10.1021/ja00299a024"]}, {"article-title": ["Maestro 7.0"], "source": ["Portland OR"], "year": ["2004"], "publisher-name": ["Schrodinger, LLC"]}, {"surname": ["Sheridan", "Miller", "Underwood", "Kearsley"], "given-names": ["RP", "MD", "DJ", "SK"], "article-title": ["Chemical Similarity Using Geometric Atom Pair Descriptors"], "source": ["J Chem Inf Comput Sci"], "year": ["1996"], "volume": ["36"], "fpage": ["128"], "lpage": ["36"], "pub-id": ["10.1021/ci950275b"]}, {"surname": ["Mohamadi", "Richards", "Guida", "Liscamp", "Lipton", "Caufield", "Chang", "Hendrickson", "Still"], "given-names": ["F", "NG", "W", "R", "M", "C", "G", "T", "W"], "article-title": ["Macromodels: An integrated software system for modeling or-ganic and bioorganic molecules using molecular mechanics"], "source": ["J Med Chem"], "year": ["1990"], "volume": ["11"], "fpage": ["440"], "lpage": ["467"]}]
BMC Struct Biol. 2008 Sep 3; 8:38