Determination of the binding mode for anti-inflammatory natural product xanthohumol with myeloid differentiation protein 2

Abstract

It is recognized that myeloid differentiation protein 2 (MD-2), a coreceptor of toll-like receptor 4 (TLR4) for innate immunity, plays an essential role in activation of the lipopolysaccharide signaling pathway. MD-2 is known as a neoteric and suitable therapeutical target. Therefore, there is great interest in the development of a potent MD-2 inhibitor for anti-inflammatory therapeutics. Several studies have reported that xanthohumol (XN), an anti-inflammatory natural product from hops and beer, can block the TLR4 signaling by binding to MD-2 directly. However, the interaction between MD-2 and XN remains unknown. Herein, our work aims at characterizing interactions between MD-2 and XN. Using a combination of experimental and theoretical modeling analysis, we found that XN can embed into the hydrophobic pocket of MD-2 and form two stable hydrogen bonds with residues ARG-90 and TYR-102 of MD-2. Moreover, we confirmed that ARG-90 and TYR-102 were two necessary residues during the recognition process of XN binding to MD-2. Results from this study identified the atomic interactions between the MD-2 and XN, which will contribute to future structural design of novel MD-2-targeting molecules for the treatment of inflammatory diseases.

Keywords:

• myeloid differentiation 2
• xanthohumol
• binding mode
• inflammation
• molecular dynamics simulation

Introduction

Inflammatory responses are triggered by pathogenic bacteria or toxin-mediated injury, which are mediated by pattern receptors expressed by immune cells.¹ However, the subsequent production of toxic proinflammatory molecules may result in undesirable outcomes, largely due to collateral tissue damage. Such uncontrolled inflammation plays an important role in the pathogenesis of diseases, such as cardiovascular diseases and sepsis.^2–4 For instance, the extensively studied immunostimulatory component of bacteria, lipopolysaccharide (LPS), causes systemic inflammation and sepsis.^5–7

LPS-induced inflammation is mediated by the toll-like receptor 4 (TLR4). In this pathway, LPS-binding protein captures LPS, transferring it to CD14 (cluster of differentiation 14)⁸ for delivery to a coreceptor, myeloid differentiation protein 2 (MD-2).⁹ TLR4 is a type-I transmembrane protein with a characteristic of horseshoe/solenoid-shaped ectodomain. MD-2, contains 160 amino acid residues, adopts a β-cup fold with two antiparallel β-sheets that forms a deep hydrophobic pocket for LPS recognition.^10,11 LPS directly binds to the hydrophobic pocket in MD-2 protein. Upon LPS binding, the TLR4/MD-2 complex undergoes dimerization and recruits its downstream adaptors.^9,10 Two major downstream adaptor molecules of MyD88 and TIR-domain-containing adapter-inducing interferon-β are triggered by the dimerization of two TLR4/MD-2 complexes, which then activates the intracellular proinflammatory cascades.¹² In general, MD-2 plays an indispensable role in the activation of the LPS-TLR4 signaling pathway and thus has been confirmed as a suitable target for the therapeutic inhibition of TLR4 signaling.^13,14

The study of MD-2, as a target for anti-inflammatory therapy, has drawn great attention and a variety of MD-2 inhibitors have been reported, among which lipid IVa is a popular lipid-like antagonist. Lipid IVa, an underacylated lipid A with four fatty acid chains from Escherichia coli, is an intermediate in the biosynthetic pathway of LPS and can inhibit LPS with four acyl chains invading the hydrophobic pocket of MD-2.¹⁵ In the meantime, some other natural and synthetic nonlipid compounds have also been tested. For example, JTT-705 and auranofin form covalent bonds with residue Cys-133 of MD-2, while compounds like taxanes and curcumin form noncovalent interactions embedded into the hydrophobic pocket of MD-2.^16–18

Xanthohumol (XN) (Figure 1A) is a prenylated chalcone-type flavonoid of Humulus lupulus L. (Cannabaceae) and it is also used to make commercial flavoring products; XN has a variety of pharmacological functions including anti-inflammation, prevention and treatment of diabetes and antioxidation.^19–21 Previous researchers suggested that XN’s anti-inflammatory activities and chemopreventive effects are attributed to blocking LPS binding to MD-2. Peluso et al²² reported that XN can suppress LPS-stimulated inflammatory signaling through direct competitive binding to the hydrophobic cavity of MD-2. However, the binding mechanism of XN in MD-2 pocket are still unclear, preventing effective therapeutic design of XN targeting the MD-2 for the treatment of inflammatory diseases.²² In the present study, we used both computational and experimental methods to explore the interaction between XN and MD-2. Our results indicated that XN can embed into the hydrophobic pocket of MD-2 and form two stable hydrogen bonds with residues ARG-90 and TYR-102.

Figure 1 XN’s binding activity to MD-2 protein.

Notes: (A) The structure of XN. (B) SPR shows direct interactions between XN and the wild type of MD-2 protein. (C) ELISA assay showing XN inhibition of biotin-LPS binding to wild-type MD-2 protein. Data are mean values (± SEM) of at least three separate repeated experiments (**P<0.01).
Abbreviations: ELISA, enzyme-linked immunosorbent assay; LPS, lipopolysaccharide; MD-2, myeloid differentiation protein 2; SEM, standard error of the mean; SPR, surface plasmon resonance; XN, xanthohumol; K_a, association (‘on rate’); K_d, dissociation rates (‘off rate’); K_D, equilibrium dissociation constant (‘binding constant’); WT, wild type.

Materials and methods

Reagents

The mouse RAW 264.7 macrophages were purchased from ATCC (Manassas, VA, USA). The recombinant human MD-2 (rhMD-2, or rhMD-2 mutant) proteins were purchased from Biowit Technologies (Shenzhen, People’s Republic of China); XN, LPS (from Salmonella typhosa); and the anti-MD-2 antibody from eBioscience (San Diego, CA, USA). No ethical approval was required for the use of these cell lines.

Surface plasmon resonance analysis

The binding affinity of XN to rhMD-2 and rhMD-2 mutants was determined using a ProteOn XPR36 Protein Interaction Array system (Bio-Rad Laboratories, Hercules, CA, USA) with an HTE sensor chip (ProteOn™, #176-5033). Briefly, MD-2 protein (in acetate acid buffer pH 5.5) was loaded to the sensors that were activated with 10 mM NiSO₄. The XN samples (at 100, 50, 25, 12.5, and 6.25 µM) were prepared with running buffer (phosphate buffered saline [PBS], 0.1% sodium dodecyl sulfate, 5% dimethyl sulfoxide). Sensor and sample plates were placed on the instrument. The XN samples were then captured in flow cell 1, leaving the second flow cell as a blank. Five concentrations were injected simultaneously at a flow rate of 30 µM/min for 120 seconds of association phase, followed with 120 seconds of dissociation phase at 25°C. The final graphs were obtained by subtracting blank sensorgrams from the duplex or quadruplex sensorgrams. Data were analyzed with ProteOn manager software. K_D (equilibrium dissociation constant [‘binding constant’]) was calculated by global fitting of the kinetic data from various concentrations of XN using 1:1 Langmuir binding model.

Enzyme-linked immunosorbent assay

A 96-well microplate was coated overnight at 4°C with anti-human MD-2 antibody in 10 mM Tris–HCl buffer (pH 7.5). The plate was washed with PBS Tween-20 and blocked with 3% bovine serum albumin in high purity water for 1.5 hours at room temperature. rhMD-2 (4 µg/mL) in 10 mM Tris–HCl buffer (pH 7.5) was added to a precoated plate and incubated for 1.5 hours at room temperature. After washing with PBS Tween-20, biotin-labeled LPS (InvivoGen, San Diego, CA, USA) was incubated for 1 hour at room temperature. Compound XN (0.1 or 1 µM) was added to the biotin-LPS before addition to the wells. After further washing streptavidin-conjugated horseradish peroxidase (Beyotime, Shanghai, People’s Republic of China) was added to the wells for 1 hour at room temperature. The horseradish peroxidase activity was determined using TMB substrate solution (eBioscience, San Diego, CA, USA). The optical density of each well was measured at 450 nm. Data were expressed as mean values (± standard deviation) of three separate experiments, each performed in duplication.

Molecular docking

The binding pose of probe XN in MD-2’s binding site was predicted by the software Glide (Schrödinger, Inc., NY, USA). The crystal structure of MD-2 in complex with lipid IVa was downloaded from the protein data bank (PDB code: 2E59) and prepared using the Protein Preparation Wizard in Maestro (Schrödinger, Inc.).¹⁵ The process included removal of all nonbonded heteroatoms and water molecules. Hydrogen bonds were added and optimized to the structure. Other preparation steps involved removal of bad contents, optimization of bond lengths and creation of disulfide bonds. The XN was built by using Maestro and converted to a 3D structure from the 2D structure using LigPrep (Schrödinger, Inc.). The resulting structures were saved in Maestro format.

The crystal structure of MD-2 was used to generate grid with van der Waals radius scaling 1.00 and partial charge cut-off at 0.25. The scoring grid was generated by enclosing the residues 30 Å around lipid IVa in the binding site using Receptor Grid Generation (Schrödinger, Inc.), following the standard procedure. The receptor grid file and the prepared XN were docked using Glide standard precision, while the ligand sampling was set to be flexible. Glide score, an empirical docking scoring functions implemented by the Optimized Potentials for Liquid Simulations 2005 force field, was used to infer the affinity and further analyze the binding mode and molecular dynamics (MD) simulations.

MD simulations

MD simulations of XN/MD-2 complex system and apo MD-2 system were carried out, respectively. The initial XN/MD-2 complex came from the top ranked binding pose in the glide docking results. Before MD simulation, XN was extracted from the XN/MD-2 complex and the charge was prepared by gaussian09 and antechamber.^23–25 Molecular mechanics parameters from the ff99SB and GAFF force fields were assigned to the protein and the ligand, respectively, using the LEaP module of AMBER (Assisted Model Building with Energy Refinement) 11 software packages.^26–29 The MD-2 structure of apo MD-2 system was prepared by Protein Preparation Wizard of Maestro. All the water molecules and the specified ligand were deleted. The two systems were all solvated in a box of TIP3P water molecules with a hydration shell of 10 Å. In addition, an appropriate number of chloride ions were used to neutralize these two systems.

MD simulation was carried out with AMBER 11 software packages and the two systems were using the same protocol.^30,31 Prior to the MD productive simulation, we performed an equilibration protocol consisting of an initial minimization of the water box of 1,000 steps, 500 steps for the steepest descent and 500 steps in the conjugate gradient. Then, the TIP3P water box was heated at constant volume until 300 K using a time constant for the heat bath with a coupling time of 100 ps. Equilibration was at 300 K and constant pressure of 100 ps. Before the production of MD, the whole system was equilibrated at 100 ps at a constant pressure of 1 bar and turned on the Langevin temperature scaling with a collision frequency of 2 ps. The last step used 50 ns of MD simulation without any restraints. Furthermore, nonbonded interactions were cut off at 8.0 Å. Periodic boundary conditions were turned on in every step of the whole process.³² Particle mesh Ewald was also performed to deal with the long-range electrostatic interactions under periodic boundary conditions.³³ The SHAKE method was used to constrain hydrogen atoms and the time step was set to 2 fs.^34,35 The coordinates were saved every 10 ps for the subsequent analysis.

Results and discussion

Determination of the direct binding of XN to MD-2

Previous studies reported that MD-2 is a potential target for XN but the binding mode of XN to MD-2 is still unknown.²² In this study, we performed biochemical and computational experiments to determine the molecular interactions between MD-2 and XN. At the beginning, surface plasmon resonance (SPR) experiments were carried out to assess the binding affinity of XN to MD-2. As shown in Figure 1B, XN bound to the MD-2 protein in a dose-dependent manner with a relatively high affinity with a K_D value of 0.00046 M. Furthermore, using a biotin–streptavidin-based enzyme-linked immunosorbent assay (ELISA) system, we found that biotin-tagged LPS (biotin-LPS) bound rhMD-2 previously coated onto plates, while coincubation with XN significantly reduced binding of biotin-LPS with rhMD-2 (Figure 1C). These results suggest that XN has the same binding site as LPS with a higher binding affinity.

The binding mode of XN to MD-2

Molecular docking was used to analyze the binding mode of XN to MD-2. Docking results showed that XN could be inserted into the large hydrophobic binding pocket of MD-2 and easily occupy a large portion of the LPS binding site (Figure 2A and B), which was consistent with our previous experimental results. In addition, other important interactions to MD-2 were also observed such as hydrogen bonding to residues ARG-90, GLU-92, and TYR-102 (Figure 2C). As with previous study, XN and related prenylated flavonoids had the same interactions with the TYR-102 residue, underscoring the role of TYR-102 in XN’s binding to MD-2. Our results were similar with the previously reported study. Additionally, based on the molecular docking results, we selected a docked complex that interacted with TYR-102 with lowest binding energy.

Figure 2 Molecular docking analysis of XN to the activity cavity of MD-2.

Notes: (A) XN overlapped with LPS in the binding site of MD-2 together with TLR4. TLR4 (cyan), MD-2 (wheat), LPS (green sticks), XN (blue spheres). (B) Overlapped region of XN (blue sticks) and LPS (green sticks). (C) Molecular docking analysis of three-dimensional binding pose between XN and MD-2. MD-2 (wheat cartoon), XN (cyan sticks), hydrogen bonds (black dotted lines).
Abbreviations: LPS, lipopolysaccharide; MD-2, myeloid differentiation protein 2; TLR4, toll-like receptor 4; XN, xanthohumol.

The stability of the docked complex

The results from the above docking analysis lead us to explore the dynamic behavior of the XN/MD-2 complex. In order to confirm the stability of docked conformation and to probe the function of XN on MD-2, conventional MD simulations on XN/MD-2 complex and apo MD-2 were conducted. For the initial structure as reference, the root-mean square deviation (RMSD) was detected for all Cα atoms in the proteins and all heavy atoms of XN. From our study, the trajectories that we established were more stable after 10 ns, therefore, only the latter part of the trajectories were considered for further analysis. As shown in Figure 3A, after 10 ns of MD simulation, the RMSD values of Cα atoms for MD-2 and XN gradually stabilized to 3.0 and 0.4 Å, respectively. However, the apo MD-2 converged the RMSD values in the system at 3.5 Å after 10 ns simulation (Figure 3B). Therefore, the above results suggest that two systems became stable in the end. Then, we extracted the initial and last snapshots from MD trajectories of two systems to verify the conformation change of the two systems. As depicted in Figure 3C and D, the XN/MD-2 complex should be stable; the docked XN/MD-2 complex in the end snapshot of the MD simulation structure and the initial structure were examined to show that they overlapped with slight positional difference. But the apo MD-2 showed a significant deviation until the end of the simulation, resulting in a backbone RMSD of about 3.5 Å that revealed a greatly changed conformation. What’s more, to quantify the “openness” of the hydrophobic cavity of MD-2, we cited a collective distance (D_total) defined as the total distance between the center of mass of MD-2 and to each center of mass of β strand.³⁶ In Figure 4A, the time evolution of D_total for apo MD-2 as well as the values of D_total became smaller about 5 Å than from the beginning (Video S1). On the contrary, in Figure 4B, values of D_total of XN/MD-2 complex increased slightly, which indicate that XN can stabilize the XN/MD-2 complex (Video S2). Garate and Oostenbrink³⁶ also had previously reported the same result with the simulated apo MD-2 and showed that the cavity of MD-2 closes rapidly in the absence of any ligand, indicating a strong hydrophobic nature and the instability of apo MD-2. Therefore XN not only directly binds to MD-2, but also plays a crucial role in the hydrophobic cavity of the MD-2 conformation.

Figure 3 Backbone RMSDs are shown as a function of time for apo MD-2 and XN/MD-2 complex structures at 50 ns.

Notes: (A) Time evolution of the RMSD of MD-2 and XN are shown with black lines and red lines, respectively. (B) The time evolution of the RMSD of apo MD-2. (C) The XN/MD-2 complex system. The alignments of starting (wheat) and final (blue) structures. (D) Apo MD-2 system. The alignments of starting (wheat) and final (blue) structures.
Abbreviations: MD-2, myeloid differentiation protein 2; RMSD, root-mean square deviation; XN, xanthohumol.

Figure 4 The “openness” of apo MD-2 and XN/MD-2 complex.

Notes: (A) Total distance between center of mass of apo MD-2 and every center of mass of β strand. (B) Total distance between center of mass of XN/MD-2 complex and every center of mass of β strand.
Abbreviations: MD-2, myeloid differentiation protein 2; XN, xanthohumol.

The main driving force for XN/MD-2 binding is a combination of several strong hydrogen bonds, which plays an important role in the protein-ligand recognition process. An analysis of the hydrogen bonds was then carried out to explore the interactions between MD-2 and XN where the most frequent hydrogen bonding partner was considered. As shown in Figure 5A, the distance between atoms in the docked complex, which form hydrogen bonds, was monitored during the whole MD simulation process and plotted. As a result, vast profiles of hydrogen bonds were observed, which fluctuated from 2.5 to 4.5 Å with an average of 3.0–3.5 Å hydrogen bonds in the XN/MD-2 complex. Regarding the XN/MD-2 complex, amidogens of ARG-90 and the hydroxyl of TYR-120 formed hydrogen bonds with hydroxyl and methoxy groups of XN, respectively, but ARG-90 was more stable than TYR-120 (Figure 5B and C). Moreover, although the distance between the two mentioned hydrogen bonds fluctuated slightly, they stabilized slightly with a distance ~3–3.5 Å in 50 ns of MD simulation. The detailed distances of evolutions in the two hydrogen bonds are illustrated in Figure 5B and C. Lastly, it is unknown as to why the hydrogen bond between GLU-92 and XN disappeared rapidly compared to the docking complex. This may be related to an inadequate amount of attraction.

Figure 5 Key hydrogen interactions between XN and residues ARG-90 and TYR-102 of MD-2.

Notes: (A) XN (yellow sticks), residues AGR-90 and TYR-102 (white sticks covered with transparent surface), hydrogen bonds (black dotted lines); hydrogen atoms are omitted. (B) The distance for the hydrogen bonds between XN and the residue TYR-102 as a function of time. (C) The distances for the hydrogen bonds between XN and residue ARG-90 as a function of time.
Abbreviations: MD-2, myeloid differentiation protein 2; XN, xanthohumol.

The weak interaction between XN and MD-2^R90A/Y102A mutant

Based on the computational results in previous studies, we hypothesized that if residues ARG-90 and TYR-102 transform into alanine, the binding affinity may decline. To test our hypothesis, we analyzed the double-mutant protein MD-2^R90A/Y102A using SPR to determine whether these two residues would affect the binding affinity of XN to MD-2 (Figure 6A). We found that the K_D of XN bound to MD-2^R90A/Y102A was 0.939 M, which was >2,000 times higher than that of the wild-type MD-2. In addition, ELISA assay showed that XN was not able to prevent the binding of biotin-LPS with the MD-2^R90A/Y102A mutant (Figure 6B). Furthermore, the side chains in the mutation of MD-2 ^R90A/Y102A were unable to attract hydrogen bonds to XN, indicating that ARG-90 and TYR-102 were the key amino acids on MD-2 that bind XN. These two results further support the notion that the residues ARG-90 and TYR-102 are essential to the recognition process of XN binding to MD-2 proteins.

Figure 6 XN shows low binding affinity to the MD-2^R90A/Y102A mutant.

Notes: (A) SPR assay indicating loss of high affinity binding of XN with the R90A and Y102A mutants of MD-2. (B) ELISA assay of XN competition with biotin-LPS binding to the MD-2^R90A/Y102A mutant.
Abbreviations: ELISA, enzyme-linked immunosorbent assay; LPS, lipopolysaccharide; MD-2, myeloid differentiation protein 2; SPR, surface plasmon resonance; XN, xanthohumol; K_a, association (‘on rate’); K_d, dissociation rates (‘off rate’); K_D, equilibrium dissociation constant (‘binding constant’).

Conclusion

In this report we described the interaction mode between XN and MD-2 protein. SPR and ELISA experiments were performed to assess the binding affinity and the binding mode of XN to the MD-2 wild-type protein. We applied molecular docking and molecular dynamic simulations to model the dynamic behavior of molecules and macromolecular assemblies at an atomic level. Moreover, the MD-2 site-directed mutants were used to verify the residues of ARG-90 and TYR-102, both of which are crucial for XN to compete with LPS. We found that XN formed two stable hydrogen bonds to ARG-90 and TYR-102 of MD-2, stabilizing the complex. We put forward a strategy to increase hydrophobic property and add hydrogen bond donors or hydrogen bond acceptors for the purpose of forming hydrogen bonds with GLU-92, which was mentioned in the docking study. What’s more, we can retain some pharmacophores and employ scaffold hopping method, known as lead hopping, to discover novel potent MD-2 inhibitors.^37,38 In conclusion, by combining theoretical and experimental approaches, we characterized a structure-based drug design of a potential anti-inflammatory inhibitor by targeting the MD-2 protein in the TLR4 signaling pathway.

Acknowledgments

Financial support was provided by the National Natural Science Funding of China (21272179, 81473242, and 81570027), and Zhejiang Provincial Natural Science Foundation (LY14C050004).

Supplementary materials

Video S1 50 ns of apo MD-2 protein molecular dynamics simulation.

Note: The cavity of MD-2 closes rapidly in the absence of any ligand.
Abbreviation: MD-2, myeloid differentiation protein 2.

https://www.youtube.com/watch?v=xiehyaDcvYw

Video S2 50 ns of xanthohumol/MD-2 complex molecular dynamics simulation.

Notes: Xanthohumol was stable in the hydrophobic cavity of MD-2. MD-2 and xanthohumol were shown as cartoon and sticks, respectively. The sticks of upper left and lower right corner were TYR-102 and ARG-90.
Abbreviation: MD-2, myeloid differentiation protein 2.

https://www.youtube.com/watch?v=53AazFB5Kpo

Disclosure

The authors report no conflicts of interest in this work.