If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Activation of Human Toll-like Receptor 4 (TLR4)·Myeloid Differentiation Factor 2 (MD-2) by Hypoacylated Lipopolysaccharide from a Clinical Isolate of Burkholderia cenocepacia*
Departments of Chemical Sciences, University of Naples Federico II, Naples 80126, ItalyDepartment of Biotechnology, National Institute of Chemistry, Ljubljana 1000, Slovenia
Department of Chemistry and Biochemistry, Universidad CEU San Pablo, Boadilla del Monte, Madrid 28668, SpainDepartment of Biopharmaceutics and Pharmacodynamics, Medical University of Gdańsk, Gdańsk 80-416, Poland
Department of Biotechnology, National Institute of Chemistry, Ljubljana 1000, SloveniaCentre of Excellence NMR - Future Innovation for Sustainable Technologies, Ljubljana 1000, Slovenia
Department of Biotechnology, National Institute of Chemistry, Ljubljana 1000, SloveniaCentre of Excellence NMR - Future Innovation for Sustainable Technologies, Ljubljana 1000, Slovenia
Department of Microbiology and Immunology, University of Western Ontario, London N6A 5C1, CanadaCentre for Infection and Immunity, Queen's University Belfast, Belfast BT9 7AE, United Kingdom
To whom correspondence may be addressed: Dept. of Chemical and Physical Biology, Centre for Biological Research, CIB-CSIC, 28040, Madrid, Spain. Tel.: 34-0918373112; Fax: 34-915360432.
* This work was supported in part by Spanish Ministry of Economy and Competitiveness Grant CTQ2011-22724, Universidad CEU San Pablo Grants PC14/2011 and PC13/2012, the Slovenian Research Agency, Centre of Excellence NMR - Future Innovation for Sustainable Technologies, European structural funds, Cystic Fibrosis Canada, European Cooperation in Science and Technology Actions BM1003 and CM1102, and Marie Curie Initial Training Networks (ITN) GLYCOPHARM PITN-GA-2012-317297 and H2020-MSCA-ITN-2014-ETN TOLLerant. The authors declare that they have no conflict of interest with the content of this article. 1 Supported by an Airbus Military research contract. 2 Supported by a Higher Education Funding Council of England senior lecturership and a Financial and Fiscal Commission grant.
Lung infection by Burkholderia species, in particular Burkholderia cenocepacia, accelerates tissue damage and increases post-lung transplant mortality in cystic fibrosis patients. Host-microbe interplay largely depends on interactions between pathogen-specific molecules and innate immune receptors such as Toll-like receptor 4 (TLR4), which recognizes the lipid A moiety of the bacterial lipopolysaccharide (LPS). The human TLR4·myeloid differentiation factor 2 (MD-2) LPS receptor complex is strongly activated by hexa-acylated lipid A and poorly activated by underacylated lipid A. Here, we report that B. cenocepacia LPS strongly activates human TLR4·MD-2 despite its lipid A having only five acyl chains. Furthermore, we show that aminoarabinose residues in lipid A contribute to TLR4-lipid A interactions, and experiments in a mouse model of LPS-induced endotoxic shock confirmed the proinflammatory potential of B. cenocepacia penta-acylated lipid A. Molecular modeling combined with mutagenesis of TLR4-MD-2 interactive surfaces suggests that longer acyl chains and the aminoarabinose residues in the B. cenocepacia lipid A allow exposure of the fifth acyl chain on the surface of MD-2 enabling interactions with TLR4 and its dimerization. Our results provide a molecular model for activation of the human TLR4·MD-2 complex by penta-acylated lipid A explaining the ability of hypoacylated B. cenocepacia LPS to promote proinflammatory responses associated with the severe pathogenicity of this opportunistic bacterium.
A central theme in innate immunity involves recognition of conserved microbial molecules (pathogen-associated molecular patterns) by surface receptors expressed on phagocytic cells (
). Pathogen recognition triggers cell signaling cascades, leading to activation of transcription factors such as nuclear factor-κB (NF-κB) and interferon regulatory factors, which in turn stimulate production of inflammatory cytokines (e.g. TNF-α, IL-1β, and type I interferons) (
). Lipopolysaccharide (LPS), the major component of the Gram-negative bacterial outer membrane, elicits potent innate immune responses through interactions with a receptor complex composed of Toll-like receptor 4 (TLR4)
). Depending on the amount and chemical nature of LPS released from the pathogen, TLR4·MD-2 recognition stimulates a protective immune response or leads to uncontrolled inflammation associated with high mortality (
). Lipid A, composed of an acylated glucosamine disaccharide backbone, is the LPS moiety recognized by the TLR4·MD-2 receptor complex. Lipid A bioactivity depends on its chemical structure. The number and distribution of acyl chains and the presence of the phosphate groups in the diglucosamine backbone determine the agonistic and antagonistic activities of lipid A (
). Most enteric bacteria such as Escherichia coli produce hexa-acylated bisphosphorylated lipid A, which has the highest cytokine-inducing capacity in mammals. In contrast, tetra-acylated lipid A (such as lipid IVA) and most penta-acylated lipid A forms lack activity on human cells (
FIGURE 1Structure of B. cenocepacia LPS inner core and lipid A.B. cenocepacia lipid A is heterogeneous, being composed of a mixture of penta- and tetra-acylated species (
Lipid A interacts with a hydrophobic pocket formed by two anti-parallel β-sheets of MD-2 mediating the dimerization and activation of the TLR4·MD-2 complex (
). The MD-2 hydrophobic pocket can accommodate up to five acyl chains. For agonistic hexa-acylated E. coli lipid A, five acyl chains are buried within the pocket, whereas the sixth chain lies on a channel of the MD-2 surface, building a hydrophobic region and the dimerization interface required for interaction with the TLR4 partner (referred to here as TLR4*). This arrangement enables hydrophobic interactions bridging the TLR4·MD-2·LPS heterodimer and promoting the juxtaposition of the intracellular domains, leading to activation of signal transduction (
). By contrast, antagonist lipid IVA binds to human MD-2 with all four acyl chains completely buried in the hydrophobic pocket precluding dimerization and subsequent activation (
Bacteria regulate the degree of lipid A acylation in response to environmental conditions. For example, Pseudomonas aeruginosa strains possess penta-acylated lipid A but can produce hexa-acylated lipid A during infection in cystic fibrosis (CF) patients (
), gaining the ability to elicit stronger inflammation. Intriguingly, the non-CF pathogen Porphyromonas gingivalis produces a high heterogeneous LPS lipid A whose penta-acylated isoform potently activates the NF-κB pathway in human gingival fibroblasts in a similar manner to the E. coli LPS (
Tetra- and Penta-acylated lipid A structures of Porphyromonas gingivalis LPS differentially activate TLR4-mediated NF-κB signal transduction cascade and immuno-inflammatory response in human gingival fibroblasts.
). Therefore, the molecular mechanisms underlying TLR4 dimerization and activation cannot always be predicted based on the E. coli lipid A-TLR4·MD-2 paradigm.
Burkholderia cenocepacia is an opportunistic Gram-negative bacterium causing serious infections in CF patients (
). Understanding how B. cenocepacia infection elicits inflammation is crucial to finding new ways to improve treatment of infected CF patients. The structure of the B. cenocepacia LPS (LPSBC) was described previously (
). The lipid A moiety of LPSBC consists of a mixture of penta-acylated and tetra-acylated, bisphosphorylated diglucosamine backbone with two 4-amino-4-deoxy-l-arabinose (l-Ara4N) residues linked by phosphodiester linkages (see Fig. 1) (
). The acyl chains are made of 3-(R)-hydroxyhexadecanoic acid C16:0 (3-OH), 3-(R)-hydroxytetradecanoic acid C14:0 (3-OH), and hydroxytetradecanoic acid C14:0 (see Fig. 1) (
). The structural features of LPSBC suggest that the molecule would be a very poor TLR4·MD-2 complex agonist compared with the prototypic enterobacterial LPS. In contrast, several studies show that LPSBC as well as biological and synthetic lipid A analogues is highly proinflammatory (
Lipopolysaccharide (LPS) from Burkholderia cepacia is more active than LPS from Pseudomonas aeruginosa and Stenotrophomonas maltophilia in stimulating tumor necrosis factor alpha from human monocytes.
The complete structure and pro-inflammatory activity of the lipooligosaccharide of the highly epidemic and virulent gram-negative bacterium Burkholderia cenocepacia ET-12 (strain J2315).
Chemical synthesis of Burkholderia lipid A modified with glycosyl phosphodiester-linked 4-amino-4-deoxy-β-L-arabinose and its immunomodulatory potential.
). Therefore, the current model for binding/activation of human TLR4·MD-2 complex cannot explain why LPSBC is proinflammatory. In this work, we investigated the molecular basis of the mechanism of LPSBC recognition and report for the first time that acyl chains longer than those in enterobacterial lipid A together with the l-Ara4N residues allow B. cenocepacia lipid A (LABC) to fit into the binding pocket of MD-2 in a manner that promotes TLR4 dimerization, leading to activation of inflammatory responses in cellular and animal models.
Experimental Procedures
LPS Extraction and Purification
LPS was prepared from B. cenocepacia strains MH71 and MH75 (
). Both strains are isogenic derivatives of the K56-2 clinical isolate and contain a deletion of the wbiF gene that eliminates O-antigen production. MH75 also has a deletion removing genes involved in UDP-l-Ara4N synthesis and therefore cannot produce LPS with l-Ara4N residues (LPSBCΔAra). Both MH71 and MH75 also carry a suppressor mutation in the LPS transport gene lptG that allows for the transport of LPS devoid of l-Ara4N to the bacterial outer membrane (
). The P. aeruginosa RP73 clinical isolate was obtained from a chronically infected CF patient and kindly provided by Prof. Burkhard Tümmler (Klinische Forschergruppe, Medizinische Hochschule Hannover, Hannover, Germany) (
). Strains were plated on trypticase soy agar plates and cultured in trypticase soy broth at 37 °C. For large scale LPS purification, bacteria were treated with hot phenol/water (
). After extensive dialysis against distilled water, the extracted phases were subjected to enzymatic digestions to remove nucleic acids and protein contaminants. Water and phenol fractions were analyzed by 13.5% SDS-PAGE and silver staining (
). The LPS fraction was exclusively found in the water phase.
Isolation of Lipid A
Lipid A was obtained by hydrolysis of the LPS with 100 mm sodium acetate buffer, pH 4.4 (100 °C for 3 h). The solution was extracted three times with CHCl3/MeOH/H2O (100:100:30, v/v/v) and centrifuged (4 °C at 5000 × g for 15 min). The organic phase contained lipid A, and the water phase contained the core oligosaccharide. The former was further purified through several washes with distilled water and then lyophilized. The lipid A and core oligosaccharide structures of strains K56-2, MH71, and MH75 were characterized elsewhere (Fig. 1) (
Biosynthesis and structure of the Burkholderia cenocepacia K56-2 lipopolysaccharide core oligosaccharide: truncation of the core oligosaccharide leads to increased binding and sensitivity to polymyxin B.
Animal studies were conducted according to protocols approved by the San Raffaele Scientific Institute (Milan, Italy) Institutional Animal Care and Use Committee and adhered strictly to the Italian Ministry of Health guidelines for the use and care of experimental animals. Efforts were made to minimize the number of animals used and their suffering. C57BL/6 mice (20–22 g; male; Charles River) were challenged via intraperitoneal injection with 300 μg/mouse LPS from E. coli and P. aeruginosa RP73, LPSBC, and LPSBCΔAra. Control mice were challenged with sterile saline solution. Five hours after treatment, mice were sacrificed by CO2 administration, and blood was collected from heart puncture. The blood clot was left at room temperature for 15–30 min. The clot was removed by centrifugation at 2000 × g for 10 min at 4 °C. Sera were collected and stored at −20 °C. TNF-α concentration in sera was determined by ELISA (R&D Systems) according to the manufacturer's instructions using antibody pairs and recombinant standards from R&D Systems.
Human Embryonic Kidney (HEK) 293 Cell Activation and Luciferase Reporter Assays
Expression plasmids containing sequences of human TLR4 and MD-2 and the pELAM-1 firefly luciferase plasmid were a gift from Dr. C. Kirschning (Technical University of Munich, Munich, Germany). Expression plasmid containing sequence of mouse TLR4 was purchased from InvivoGen. Expression plasmid for mouse MD-2 was a gift from Dr. Y. Nagai (University of Tokyo, Tokyo, Japan). The Renilla luciferase phRL-TK plasmid was purchased from Promega. Recombinant MD-2 genes were cloned into pEF-BOS with FLAG and His tags on the C terminus. Recombinant TLR4 genes were cloned into pUNO with a C-terminal HA tag. Transfection reagent JetPEI was purchased from Polyplus Transfection (France) and was used according to the manufacturer's instructions. The HEK293 cells were provided by Dr. J. Chow (Eisai Research Institute, Andover, MA) and grown in DMEM supplemented with 10% FBS. The MD-2 mutants were made using the QuikChange site-directed mutagenesis kit (Stratagene) according to the manufacturer's instructions. All plasmids were sequenced to confirm the appropriate mutations. For the NF-κB-luciferase reporter assay, HEK293 cells were seeded in 96-well plates at 3 × 104 cells/well and incubated overnight in a humidified atmosphere (5% CO2) at 37 °C. The next day when cells were 60–80% confluent, they were co-transfected with the plasmids for MD-2 (10 ng), TLR4 (1 ng), NF-κB-dependent luciferase (50 ng), and constitutive Renilla luciferase (10 ng) using JetPEI transfection reagent (all amounts are in ng/well). Cells were stimulated 6 h after transfection with endotoxin preparations. Cells were lysed after 16 h of stimulation in 1× reporter assay lysis buffer (Promega) and analyzed for reporter gene activities using the Dual-Luciferase reporter assay system. Relative luciferase units were calculated by normalizing the luciferase activity of each sample with the constitutive Renilla activity measured within the same sample.
Statistical Analysis
Results were expressed as mean ± S.D. from experiments done in triplicate. Statistical calculations and tests in vitro and in vivo were performed using Student's t test considering p ≤ 0.05 as the limit of statistical significance and a p value of 0.001 as extremely significant.
Molecular Modeling: Building and Geometry Optimization of LPSBC and LPSBCΔAra
Three-dimensional coordinates were built by Maestro (
). Molecular mechanics optimization (universal force field), semiempirical calculations (AM1), and density functional theory (B3LYP/6-31G*) were subsequently applied using Gaussian 03 (
). Ammonium and carboxylic acid groups were considered as ionized. Conformational analysis on DS1, DS2, and DS3 disaccharides was performed with MacroModel (
). The parameters were as follows: MM3* force field and water solvent dielectric constant of 1.0. The charges from the force field were used with an extended cutoff. The best conformers were selected, and from them, molecular dynamics simulations (MDSs) with implicit water and MM3* as the force field were performed using Schrödinger Maestro 9.3 and Impact 5.8 (
) (MM3* force field; dielectric constant, 80.0; number of MDS steps, 100; time step (ps), 0.001).
The full three-dimensional structures of LABCΔAra, LABC, and LPSBC were built with Maestro using optimized disaccharide scaffolds obtained from the conformational analysis. Fatty acid chains, DS6, and aminoarabinoses were optimized separately using Gaussian 03 (B3LYP/6-31G*). The starting conformation for lipid A optimization was obtained from E. coli lipid A from Protein Data Bank code 3fxi. The full structure of LPSBC was built by putting together core region disaccharides and lipid A with five aliphatic chains. Final MDS of the geometry was performed with implicit water by means of Impact (MM3* force field; dielectric constant, 80.0; number of MDS steps, 100; time step (ps), 0.001). Three-dimensional structures of lipid A molecules for validation were extracted from their corresponding Protein Data Bank files (E. coli lipid A was obtained from Protein Data Bank code 3fxi, and Re chemotype of E. coli lipid A was obtained from Protein Data Bank code 3vq2), refined with the help of Maestro, and finally submitted to MDS with Impact (MM3* force field and implicit water).
Molecular Modeling
For docking studies with MD-2 and TLR4 proteins, three-dimensional coordinates from MD-2 protein were obtained from Protein Data Bank codes 3fxi (human MD-2) and 3vq2 (murine MD-2) and refined and minimized with the Protein Preparation Wizard module of Maestro using the AMBER force field (
). The TLR4 structure was obtained from Protein Data Bank code 3fxi and treated following the same procedure after removing water molecules and any other non-standard residue.
Two different docking methodologies were used for docking studies of LABCΔAra and LABC on MD-2: AutoDock and AutoDock Vina. Although AutoDock is a widely used semiempirical docking method (
). For human MD-2 (from Protein Data Bank code 3fxi), the Autogrid grid point spacing was set at 0.375 Å; center coordinates of the grid box were 29.00, −7.00, and 17.875 (x, y, z); and number of grid points in xyz was 41, 53, and 83. For murine MD-2 (from Protein Data Bank code 3vq2), the Autogrid grid point spacing was set at 0.375 Å and center coordinates of the grid box were −27.50, −15.50, and 22.00 (x, y, z), leading to 75 × 40 × 60 (x, y, z) grid points. The best result from each docking job with LABC was used as the starting geometry for subsequent docking calculations. Different combinations of allowed rotatable bonds were also considered for the ligands. Docking calculations with AutoDock were performed using a genetic algorithm (number of individuals in population, 150; maximum number of energy evaluations, 2,500,000–5,000,000; maximum number of generations, 27,000; number of top individuals to survive to next generation, 1; rate of gene mutation, 0.02; rate of crossover, 0.8; window size, 10; α parameter of Cauchy distribution, 0.0; β parameter of Cauchy distribution, 1.0). When docking LPSBC, flexible docking was also performed considering Asp-101, Glu-120, and Glu-122 as flexible residues. Docking calculations with AutoDock Vina were also performed. Coordinates and dimensions of grid boxes, starting geometries, and general methodology were the same as for AutoDock. When docking LPSBC, flexible docking was also performed considering Asp-101 as a flexible residue. Three-dimensional structures of the docked complexes were optimized by MDS with Impact (implicit water and AMBER* force field).
For docking experiments studies of the B. cenocepacia LPS core on human TLR4, we used AutoDock Vina 1.1.2. Torsional bonds from the ligand were allowed to rotate. Ammonium groups were considered as ionized. The grid box was built on TLR4 from Protein Data Bank code 3fxi with grid center at 26.00, −22.00, and 11.50 (x, y, z), number of grid points of 41 × 41 × 43 (x, y, z); and spacing of 0.375 Å. Final optimization of the three-dimensional structure of the docked complex was carried out by MDS with Impact (implicit water and AMBER* force field).
The full complex of human TLR4·MD-2 with B. cenocepacia LPS core was built by merging docked MD-2·LPS core complex with docked LPS inner core·TLR4 using Protein Data Bank code 3fxi as a template. The resulting structure was optimized by MDS with implicit water (AMBER* force field). Coupling of two TLR4·MD-2·LPSBC complexes was finally performed also using Protein Data Bank code 3fxi as a template. The full three-dimensional structure of the dimer complex was submitted to MDS with implicit water and AMBER* force field (charges from, force field; cutoff, extended; method, Polak-Ribière conjugate gradient; maximum iterations, 500; converge on, gradient; convergence threshold, 0.05). This final structure for the TLR4·MD-2·LPSBC complex was used to generate the mutant TLR4·MD-2·LPSBC complex (D294A, R322A, S415A*, and S416A*) by changing the corresponding side chains to Ala side chain. The resulting structure was submitted to MDS with implicit water and AMBER* force field. Energy analysis was performed by means of the MM-ISMSA method (
LPS and Lipid A from B. cenocepacia Activate the TLR4·MD-2 Complex and Have Proinflammatory Activity in Vivo
TLR4·MD-2 activation by LPSBC and LABC was examined using the NF-κB reporter luciferase assay in transiently transfected HEK293 and HEK293T cells co-expressing murine (mMD-2) or human MD-2 (hMD-2) and murine (mTLR4) or human TLR4 (hTLR4). Initial results revealed that LPSBC and LABC activate both human and mouse TLR4·MD-2 complexes (FIGURE 2, FIGURE 3, A and B). Because LABC possesses longer acyl chains than those in E. coli lipid A and constitutively carries one or two l-Ara4N residues, to further understand the structure-function relationship of LPSBC, we utilized LPS from a mutant strain lacking the ability to produce l-Ara4N (LPSBCΔAra). NF-κB reporter luciferase assays demonstrated that at 5 ng/ml LPSBC and LPSBCΔAra induced activation of TLR4·MD-2 complexes (Fig. 2). Similar results were obtained in dose-response experiments with purified LABC and LABCΔAra (at 5, 10, and 50 ng/ml) in cells co-expressing mTLR4·mMD-2 and hTLR4·hMD-2 (Fig. 3, A and B, respectively), which also demonstrated consistent agonistic activity at the lowest concentration of 5 ng/ml. However, NF-κB activation by mouse (Fig. 3A) and human (Fig. 3B) complexes was significantly lower with LABCΔAra, suggesting that the l-Ara4N modification of the lipid A plays a role in bioactivity of the LPS molecule.
FIGURE 2LPSBC activates the murine and human TLR4·MD-2 complexes. NF-κB activation upon stimulation of HEK293 mTLR4·mMD-2 and hTLR4·hMD-2 after 6 h with 5 ng/ml LPSBC and LPSBCΔAra. Stimulation for 6 h with E. coli LPS was used as the control. The data are pooled from three independent experiments done in triplicate. Bars indicate S.D.; significance was calculated in comparison with stimulation with E. coli LPS (*, p < 0.05; **, p < 0.01). Curly brackets indicate significance calculated comparing LPSBC and LPSBCΔAra (*, p < 0.05; **, p < 0.01). NS, no stimulation.
FIGURE 3Lipid A activation of murine and human TLR4·MD-2 complexes.A, NF-κB activation upon stimulation of HEK293 mTLR4·mMD-2 after 6 h with LABC and LABCΔAra. Stimulations for 6 h with E. coli LPS, hexa-acylated lipid A, and lipid IVA were used as controls. The data are pooled from three independent experiments done in triplicate. Bars indicate S.D.; significance was calculated in comparison with stimulation with hexa-acylated E. coli lipid A (*, p < 0.05; **, p < 0.01; ***, p < 0.001); curly brackets indicate significance calculated comparing LPSBC and LPSBCΔAra (**p < 0.01). B, NF-κB activation upon stimulation of HEK293 hTLR4·hMD-2 after 6 h with LABC and LABCΔAra. Stimulations for 6 h with E. coli LPS, hexa-acylated lipid A, and lipid IVA were used as controls. The data are pooled from three independent experiments done in triplicate. Bars indicate S.D.; significance was calculated in comparison with stimulation with hexa-acylated E. coli lipid A (*, p < 0.05; ***, p < 0.001); curly brackets indicate significance calculated comparing LPSBC and LPSBCΔAra (*, p < 0.05). NS, no stimulation.
Comparisons with NF-κB activation elicited by synthetic E. coli-type hexa-acylated lipid A revealed significantly less activation by LPSBC and LABCΔAra (all at 50 ng/ml) (Fig. 3, A and B). Lower activation by LPSBC and LPSBCΔAra could depend on weaker binding to MD-2 or less efficient activation of TLR4. Therefore, we examined whether LPSBC and LPSBCΔAra could interfere with TLR4 signaling elicited by the E. coli LPS by a competition assay in which HEK293 cells transfected with hTLR4·hMD-2 were preincubated with LPSBC or LPSBCΔAra for 1 h and then restimulated with E. coli LPS for 4 h (Fig. 4). Preincubation of cells with 10 ng/ml LPSBC followed by stimulation by 1 ng/ml E. coli LPS revealed weaker NF-κB activation than that with E. coli LPS alone (p < 0.05; Fig. 4). Interestingly, a 2-fold reduction in NF-κB activation was evident in the case of preincubation with 10 ng/ml LPSBCΔAra (p < 0.01; Fig. 4). Furthermore, a significant difference was found in activation by 10 ng/ml E. coli LPS alone compared with that obtained after preincubation of cells with 100 ng/ml LPSBC and LPSBCΔAra (p < 0.01; Fig. 4). In contrast, this inhibitory effect was reduced by increasing the concentration of the preincubated LPSBC and LPSBCΔAra (100 ng/ml) added to the 1 ng/ml E. coli LPS; such effect is probably due to the contribution of a weak activation of MD-2·TLR4 by LPSBC and LPSBCΔAra. Collectively, these results suggest that LPSBCΔAra efficiently binds to TLR4·MD-2 but is a weaker agonist compared with the E. coli LPS and LPSBC.
FIGURE 4B. cenocepacia LPSBCΔAra effects on the E. coli LPS agonist activity. The potential antagonist activity of LPSBC and LPSBCΔAra on hexa-acylated E. coli LPS was assayed. NF-κB activation upon stimulation of HEK293 hTLR4 after 1 h with LPSBC (1, 10, and 100 ng/ml) and LPSBCΔAra (1, 10, and 100 ng/ml) followed by exposure to E. coli LPS (1 and 10 ng/ml) for 4 h is shown. The data are pooled from three independent experiments done in triplicate. Bars indicate S.D.; significance was calculated in comparison with stimulation with E. coli LPS (1 and 10 ng/ml) (*, p < 0.05; **, p < 0.01). NS, no stimulation.
The proinflammatory activity of LPS from B. cenocepacia was assessed using an established in vivo model of endotoxic shock. C57Bl/6 mice were challenged by intraperitoneal injection of LPS from E. coli and P. aeruginosa RP73 (a CF strain expressing penta- and tetra-acylated lipid A (
)), LPSBC, and LPSBCΔAra. After 5 h of treatment, TNF-α levels in sera were measured to evaluate proinflammatory and endotoxic potential. Compared with the saline solution control, penta-acylated LPSBC induced robust TNF-α production (LPSBCversus control, p < 0.01; Fig. 5), whereas LPSBCΔAra induced 3-fold less systemic release of TNF-α in sera (LPSBCΔAraversus LPSBC, p < 0.01; Fig. 5), suggesting a role for l-Ara4N residues in the endotoxic LPS response in agreement with the in vitro results (Fig. 2).
FIGURE 5Proinflammatory and endotoxic potential of LPSBC in C57Bl/6 mice. From three to five mice per group were challenged via intraperitoneal injection with 300 μg/mouse LPS from E. coli and P. aeruginosa RP73, LPSBC, and LPSBCΔAra. TNF-α levels in sera were quantified after 5 h of treatment. Treatment with sterile saline solution was used as control (Ctrl). The data are pooled from two independent experiments. Bars indicate mean ± S.D. Statistical analysis was made for pair wise comparisons (*, p < 0.05; **, p < 0.01).
Identification of MD-2 Amino Acid Residues Involved in the Interaction with LPSBC
To further investigate the molecular recognition of LPSBC and LPSBCΔAra by TLR4·MD-2, we constructed targeted hMD-2 and mMD-2 mutants that were co-expressed with hTLR4 or mTLR4 in transfected HEK293 cells. We focused on residues at the dimerization interface of the MD-2 protein (
). In particular, we examined Val-82 located in proximity to the pocket entrance at a loop between β-strands 5 and 6. This region contains several conserved, solvent-exposed hydrophobic residues that contribute to crucial hydrophobic interactions in the activated receptor complex (
). Val-82 was replaced by phenylalanine (V82F) to augment the hydrophobic interactions of the MD-2·lipid A complex with TLR4. In comparison with HEK293 cells transfected with parental MD-2, hTLR4·hMD-2 V82F led to increased activation upon stimulation with LPSBC and LPSBCΔAra, whereas E. coli LPS activation was not significantly modified (Fig. 6). This agrees with the notion that hexa-acylated lipid A provides stronger hydrophobic interactions with TLR4, and therefore increased hydrophobicity at position 82 augments the activation by penta-acylated LPS. No significant differences were detected between stimulation with LPSBC and LABCΔAra on hTLR4·hMD-2 V82F. Hydrophobic interactions between mTLR4·mMD-2 and LA are expected to be similar to those involving the hTLR4·hMD-2 complex because nearly all the hydrophobic residues are conserved (
). However, the electrostatic interactions that contribute to receptor selectivity differ between mMD-2 and hMD-2 proteins. Particularly, hMD-2 has a positively charged lysine at position 122 near the pocket entrance, whereas mMD-2 has a negatively charged glutamic acid (
) and are important for the species-specific differences in the recognition of lipid IVA between hMD-2 and mMD-2. We therefore introduced amino acid replacements at positions 122 and 125 in both MD-2 homologues to test whether the substitutions could affect LA binding and influence the selectivity of TLR4 activation by the LABC variants. The E122K replacement in mMD-2 introducing the equivalent lysine of hMD-2 had little or no effect on the activation induced by LABC, LABCΔAra, or the hexa-acylated lipid A. In contrast, this replacement significantly decreased activation by lipid IVA (p < 0.001; Fig. 7) likely due to removal of the repulsive forces toward its bisphosphorylated disaccharide backbone. These repulsive interactions have been proposed to be important for positioning of lipid IVA in an appropriate orientation for TLR4 dimerization (
). Thus, the unchanged LABC/LABCΔAra responsiveness suggests that the single residue at position 122 of the murine MD-2 protein is not critical for LABC/LABCΔAra signaling.
FIGURE 6B. cenocepacia LPS acylation pattern is responsible for the NF-κB activation. NF-κB activation upon stimulation of HEK293 hTLR4·hMD-2 after 6 h with LPSBC and LPSBCΔAra is shown. Stimulation for 6 h with E. coli LPS was used as the control. The same protocol was used to stimulate HEK293 hTLR4·hMD-2 V82F cells. The data are pooled from three independent experiments done in triplicate. Bars indicate S.D.; significance was calculated in comparison with E. coli LPS (*, p < 0.05; **, p < 0.01; ***, p < 0.001). NS, no stimulation.
FIGURE 7Lys-122 and Lys-125 in the B. cenocepacia LPS signaling on murine and human TLR4·MD-2 complexes. An NF-κB-luciferase reporter assay stimulating HEK293 hTLR4·hMD-2, HEK293 mTLR4·mMD-2, and mTLR4·hMD-2 after 6 h with LPSBC and LPSBCΔAra was carried out. Stimulations for 6 h with E. coli LPS and with lipid IVA were used as controls. The same protocol was used to stimulate HEK293 mTLR4·mMD-2 E1222K and HEK293 mTLR4·mMD-2 E122K/L125K cells. The data are pooled from two independent experiments done in triplicate. Bars indicate S.D.; statistical analysis was calculated for pairwise comparisons (*, p < 0.05; ***, p < 0.001). NS, no stimulation.
Stimulation of cells expressing the double E122K/L125K mMD-2 replacement (Fig. 7) showed substantially increased activity of the mutant MD-2 regardless of the LPS variant (hexa-acylated lipid A versus hexa-acylated lipid A, p < 0.05; LABCversus LABC, p < 0.001; LABCΔAraversus LABCΔAra, p < 0.001, respectively, on mTLR4·mMD-2 and mTLR4·mMD-2 E122K L125K; Fig. 7) except with lipid IVA whose activity was strongly impaired (p < 0.001; Fig. 7). Conversely, hMD-2 in combination with mTLR4 exhibited significantly stronger activation by LPSBCΔAra (p < 0.05; Fig. 7), demonstrating that the l-Ara4N affects interaction with the TLR4 ectodomain. Moreover, the hMD-2 K122E mutant had an additive effect in combination with mTLR4 on the activation induced by LABCΔAra, showing an increment of the NF-κB activity (LABCΔAra mTLR4·hMD-2 K122E versus LABCΔAra hTLR4·hMD-2, p < 0.001; LABCΔAra mTLR4·hMD-2 K122E versus LABCΔAra mTLR4·hMD-2, p < 0.05; Fig. 7). This probably reflects the presence of repulsive forces toward the phosphate groups of the lipid A that might facilitate the orientation toward the Cys-95/Cys-105 loop (
Together, these experiments using hMD-2 and mMD-2 mutants support the notion that the peculiar structure of the LABC plays a key role in its interaction with MD-2 as well as in the dimerization process. Our data suggest that the l-Ara4N molecules on the diglucosamine backbone influence the binding of lipid A to the TLR4·MD-2 complex, although their absence is not sufficient to impair signaling.
Molecular Modeling of LPSBC Binding to TLR4·MD-2
Our previous results suggested that different MD-2 binding affinities drive the agonist activity of LPS. To provide a model for the interactions of LPSBC and LPSBCΔAra, docking calculations were undertaken, and the two ligands were docked into hMD-2 and mMD-2 proteins. Validation of the docking protocols was performed with AutoDock and AutoDock Vina by docking E. coli lipid A into the MD-2 protein from their corresponding crystallographic structures (Protein Data Bank codes 3fxi (human) and 3vq2 (murine)). Values of root mean square deviation and predicted free energy of binding values indicated an excellent performance of both programs in predicting the crystallographic binding pose for both hMD-2 and mMD-2 proteins (Table 1).
TABLE 1Theoretical free energy of binding for the docking calculations with AutoDock 4.2 and AutoDock Vina (in bold)
LABC and LABCΔAra were predicted to bind both MD-2 homologues with binding poses that agreed with the crystallographic binding poses for E. coli lipid A (superimposition with E. coli lipid A in FIGURE 8, FIGURE 9). Docking poses consisted of four fatty acid chains deeply immersed into the hydrophobic pocket of MD-2, establishing van der Waals and CH-π interactions with the side chains of most of the lipophilic residues of the MD-2 pocket, mainly consisting of aliphatic Leu, Ile, and Val residues; aromatic Phe (numbering 76, 104, 119, 121, 126, 147, and 151); and Tyr (numbering 102 and 131) residues. All docked poses were found to participate in hydrophobic interactions with most of these residues. Higher efficiency in establishing lipophilic interactions can be deduced given the length of these four fatty acid chains. The fifth chain is placed into the groove or channel defined by Phe-126, Leu-87, Val-82, and Arg-90. This channel has been identified in the x-ray structures of the complex of TLR4·MD-2 with E. coli lipid A (
) as the allocation site for one lipid chain, allowing completion of the hydrophobic interface required for dimerization with the second TLR4·MD-2·ligand partner (referred to here as TLR4*·MD-2*). This precludes formation of the activated TLR4·MD-2·ligand multimer. Also Phe-126 has been proposed as a switch controlling the agonist/antagonist conformation of MD-2 because Phe-126 mutation prevents dimerization and abolishes downstream signaling (
FIGURE 8Predicted binding mode of the LPSBC core to hTLR4·MD-2. The computational model from docking followed by MDS is shown. LPSBC docked to MD-2 is shown in Corey-Pauling-Koltun colors. Superimposed fatty acids chains from E. coli lipid A are shown in different colors. Some representative residues from the MD-2 binding site are represented with carbon atoms in green.
FIGURE 9Predicted binding mode of the LPSBC core to hTLR4·MD-2. Details of LPSBC (Corey-Pauling-Koltun colors) docked to MD-2 are shown. Superimposed fatty acids chains from E. coli lipid A are shown with carbon atoms in blue.
Therefore, the model predicts that the greater length of the five fatty acid chains compared with E. coli lipid A places the fifth chain outside the MD-2 pocket, thus building the dimerization interface in the Phe-126, Leu-87, Val-82, and Arg-90 groove, analogously to the sixth fatty acid of E. coli lipid A, whereas the other four chains remain inside the MD-2 pocket, providing stabilizing interactions.
Furthermore, the disaccharide scaffolds of LABC and LABCΔAra together with the phosphate groups are predicted to be in the outer region and to establish electrostatic interactions with the polar residues that define the rim of the MD-2 pocket (Fig. 10). The higher predicted affinity for LPSBC could be explained in terms of additional anchorage points arising from the l-Ara4N residues. In particular, these l-Ara4N molecules can establish additional interactions with MD-2. The hydroxyl (OH) group from Tyr-102 establishes bridged H-bonding with the OH at position 2 of l-Ara4N and the OH from the acyl chain (Fig. 10), whereas OH at position 3 from the same l-Ara4N establishes a hydrogen bond with the Ser-118 OH group. In the case of l-Ara4N-2, the OH group from one of the acyl chains establishes a hydrogen bond with the ys-122 CO group, whereas OH-3 and ammonium groups, although of the H-bond distance range (4.5 Å), are well oriented toward carbonyl groups from Lys-122 and Gly-123 and may enable a putative H-bond networking at this site. These interactions could also be identified in the complex with mMD-2 (data not shown). This predicted binding mode of LABC could thus explain the unexpected significant level of activation in the biological assays. The lack of some lipophilic interactions arising from the absent sixth fatty acid chain in comparison with E. coli lipid A could also be counterbalanced by these extra polar interactions from the l-Ara4N residues.
FIGURE 10Predicted binding mode of the LPSBC core to hTLR4·MD-2. Detail of some H-bond interactions between the LPSBC inner core and TLR4·MD-2 involving Tyr-102 is shown. LPSBC is represented in Corey-Pauling-Koltun colors, and MD-2 protein is represented in green.
Molecular Model of the hTLR4·MD-2 dimer in Complex with LPSBC
To develop a structural model for the TLR4·MD-2 dimer in complex with LPSBC, we first docked the inner core oligosaccharide moiety of LPSBC (Fig. 1) on hTLR4 using as a guide the region where E. coli LPS inner core binds. The best binding pose was selected (theoretical binding free energy of −2.8 kcal mol−1 by AutoDock Vina) taking also into account a proper orientation for the building of the full complex. The docked pose superimposed with E. coli LPS core from Protein Data Bank code 3fxi is shown in Fig. 11. The hTLR4·core complex was then assembled to the MD-2·LABC complex (best result from AutoDock Vina calculations), leading to a full hTLR4·MD-2·LPSBC complex. The receptor dimer was built using Protein Data Bank code 3fxi as a template. The full three-dimensional structure of the dimer complex (Fig. 12) was subjected to MDS, which did not show meaningful differences with the docked complexes (data not shown). The three-dimensional model of the full complex supports the key interactions for the molecular recognition of LPSBC that were described above and also confirms that the l-Ara4N residues (including the one present in the LPS core region; Fig. 1) play a fundamental role in the complex formation as they participate in H-bonding at the dimerization interfaces (Fig. 13). First, both ammonium groups from l-Ara4N-2 and -3 (Fig. 13) (this latter is attached to a d-glycero-d-talo-2-octulosonic acid in the core moiety; Fig. 1) establish hydrogen bonds with the TLR4 Arg-264 and Asp-294 carboxylate groups. Second, the ammonium group of l-Ara4N-1 is close to the Asp-395 carboxylate and the Ser-416 OH group, which belong to the opposite TLR4*. Third, the OH-3 from l-Ara4N-1 establishes an H-bond with the Ser-415 CO group of TLR4* (Fig. 13). The three-dimensional models from docking could suggest higher activation for LPSBCversus LPSBCΔAra arising from the extra anchorage points from arabinoses. This agrees with the higher endotoxic potential of LPSBC in the biological assays. Based on this TLR4·MD-2·LPSBC model, we built a new model with the following four mutations: D294A, R322A, S415A*, and S416A*. These mutations involve the main residues establishing interactions with the ammonium groups from the three Ara4N residues and the OH group from the terminal glucose (Fig. 13). (i) Asp-294 side chain establishes polar interactions reinforced by hydrogen bonds with the ammonium group from Ara4N-3, and is in the proximity of the ammonium group from Ara4N-2, so this residue can be considered as a main anchorage point for TLR4·MD-2·LPSBC; Asp-294 was mutated to Ala to abolish this interaction. (ii) Ser-415* and Ser-416* are placed in the partner TLR4* and provide an anchorage point for the ammonium group from Ara4N-1 through polar interactions with the side chain of Ser-415* and the backbone CO groups; both Ser residues were mutated to Ala. (iii) Although not directly in contact with any of the Ara4N moieties, the Arg-322 side chain establishes a hydrogen bond with the OH group at position 3 from the terminal glucose, and we mutated it to compare its contribution to the global energy of the system; this interaction is absent in the complex of TLR4·MD-2 with E. coli LPS (Protein Data Bank code 3fxi) and could be considered as a distinctive interaction for LPSBC core. MDS and energy calculations of both systems (TLR4·MD-2·LPSBC and mutant TLR4·MD-2·LPSBC) have shown an important difference in the global energy: the wild type complex (TLR4·MD-2·LPSBC model) is around 32 kcal mol−1 more stable than the mutant counterpart mainly due to the coulombic term, pointing to the presence of important polar interactions (Table 2). Analysis of the contributions of each residue to the ligand binding energy reveals that this energy difference is mainly due to the lack of the interactions involving the mutated residues. (i) Interaction with Asp-294 is the strongest interaction in the TLR4·MD-2·LPSBC system (around 22 kcal mol−1) and is absent in the mutant. (ii) Also the interaction with Arg-322 is missing; this is one of the main interactions in the TLR4·MD-2·LPSBC system (around 9.6 kcal mol−1). (iii) The lower coulombic contribution to the interaction energy in the mutant TLR4·MD-2·LPSBC complex is due to the absence of the interaction between the Ser-415* side chain and the ammonium group. These results give reasonable evidence that the Ara4N ammonium groups provide additional anchorage interactions accounting for the final stability of the TLR4·MD-2·LPSBC complex. Abolishing these interactions leads to a less stable complex, suggesting that they are crucial for the binding.
FIGURE 11Superimposition of best dock results of LPSBC core. The docking was obtained using AutoDock Vina (depicted in red) and E. coli LPS core from Protein Data Bank code 3fxi (depicted in green). TLR4 is not shown for the sake of clarity.
FIGURE 12Predicted binding mode of the LPSBC core to hTLR4·MD-2. The computational model from docking followed by MDS is shown. Left panel, three-dimensional model of the dimer of LPSBC core in complex with TLR4·MD-2. Right panel, detail of the dimer of LPSBC inner core in complex with TLR4·MD-2.
FIGURE 13Predicted binding mode of the LPSBC core to hTLR4·MD-2. Details of some interactions between the LPSBC inner core and TLR4·MD-2 from the computational model are shown.
TABLE 2Analysis of the per residue contributions to the total ligand binding energy for both complexes TLR4·MD-2·LPSBC and the mutated TLR4·MD-2·LPSBC (D294A, R322A, S415A*, and S416A*)
) because the sixth acyl chain protrudes from the MD-2 binding pocket, bridging TLR4·MD-2 complex dimerization. On the contrary, underacylated lipid A molecules are not (or are poorly) sensed by human TLR4·MD-2 as they are potentially accommodated within the MD-2 pocket. Here, we demonstrate a role for l-Ara4N residues in lipid A that explains why the B. cenocepacia LPS (naturally consisting of a mixture of penta- and tetra-acylated forms) acts as a strong TLR4·MD-2 agonist. Our results agree with a recent report by Hollaus et al. (
Chemical synthesis of Burkholderia lipid A modified with glycosyl phosphodiester-linked 4-amino-4-deoxy-β-L-arabinose and its immunomodulatory potential.
) demonstrating that a synthetic Burkholderia lipid A substituted with l-Ara4N exclusively at the anomeric phosphate acts as a potent human TLR4·MD-2 agonist.
Our in vitro and in vivo results explain and expand the above conclusions and show that B. cenocepacia penta-acylated lipid A (both LABC and LABCΔAra) elicits an inflammatory response activating the TLR4·MD-2 complex, suggesting that the model of the receptor complex activation should also include other structural components. Furthermore, molecular modeling experiments indicated that the increased length of the two amide-linked acyl chains of the B. cenocepacia lipid A moiety (3-(R)-hydroxyhexadecanoic acid C16:0 (3-OH) versus 3-(R)-hydroxytetradecanoic acid C14:0 (3-OH) from the E. coli lipid A) compensates for the lack of one fatty acid chain filling the hydrophobic binding pocket of MD-2 protein, allowing the placement into the channel of the fifth acyl chain responsible for TLR4 dimerization. The replacement in E. coli lipid A of the secondary dodecanoic acid C12:0 with a hexadecanoic acid C16:0 results in a weaker LPS agonist (
Expression of a Porphyromonas gingivalis lipid A palmitylacyltransferase in Escherichia coli yields a chimeric lipid A with altered ability to stimulate interleukin-8 secretion.
). However, our data indicate that not only the increased length of the acyl chain but also its position is important for agonist activity.
Furthermore, our mutagenesis studies using two B. cenocepacia lipid A forms, LABC and LABCΔAra, uncovered differential responsiveness on transfected HEK293 cells, indicating an important role also for the l-Ara4N residues in the TLR4·LABC binding process. Particularly, the significantly different immunostimulatory activity of LABC compared with LABCΔAra (Fig. 7) observed in the presence of the human MD-2 K122E mutation might be related to the occurrence on LABC of the l-Ara4N residues, which could mask the negatively charged phosphate groups and then reduce the repulsive forces introduced with the K122E mutation. In contrast, the absence of l-Ara4N in LABCΔAra might explain its increased agonistic activity when tested on HEK293 cells transfected with mTLR4·hMD-2 K122E because the repulsive forces introduced might facilitate the orientation toward the Cys-95/Cys-105 loop (
). Therefore, it is reasonable to assume that the l-Ara4N modification with its positively charged ammonium group could favor the electrostatic interactions allowing receptor/lipid A binding. Intriguingly, recent work by Maeshima et al. (
) revealed several key charged amino acid residues in TLR4 and MD-2 mediating host-specific responses to the glucosamine-modified penta-acylated lipid A from Bordetella pertussis and its unmodified counterpart by human and mouse TLR4·MD-2 complexes. This further underscores the importance of positively charged residues decorating penta-acylated lipid A and interacting with human TLR4·MD-2 charged amino acid residues to promote complex activation (
Our molecular modeling experiments further supported the notion that l-Ara4N residues provide additional polar interactions affecting the LABC binding to the TLR4·MD-2 and contribute to anchoring the lipid A into the receptor complex. The network of hydrogen bonds and polar interactions contributing to anchoring the lipid A into the TLR4·MD-2 involves not only the sugar residues from the core oligosaccharide (for example, hydrogen bonds involving Arg-322; Fig. 13) but also the three l-Ara4N residues and the partner TLR4*. As stated above, mutagenesis studies demonstrate that several Lys residues involved in LPS binding play indispensable roles through the polar interactions with the phosphate groups (
). In our dimer model, Lys-122 is involved in binding the lipid A moiety of the LPSBC through polar interactions with the core; as for V82F, its role in building the groove that accommodates the acyl chain from lipid A, completing the interaction surface required for dimerization, has already been mentioned. The replacement of Val-82 by phenylalanine involves changes of van der Waals interactions into CH-π interactions, which may favor the binding of this chain, which is longer than the corresponding chain on E. coli LPS. Thus, the structural peculiarities of B. cenocepacia, i.e. its acylation pattern and the presence of the l-Ara4N residues in the lipid A region, exert a synergistic effect in the activation and dimerization of the LPS receptor, opening new insights into the comprehension of the molecular recognition of LPS by TLR4·MD-2.
Author Contributions
A. M. designed the research. F. D. L. executed the chemistry and immunochemistry experiments. Ł. K. and S. M.-S. performed the molecular dynamics experiments. F. D. L., A. M., Ł. K., S. M.-S., M. A. V., A. O., N. I. L., C. G., R. L., M. P., M. A. H., A. D. S., A. S., R. J., and A. B. contributed with analytical tools, reagents, or particular experiments. F. D. L., M. A. V., S. M.-S., and A. M. analyzed data and wrote the paper. F. D. L., A. M., Ł. K., S. M.-S., M. A. V., A. O., N. I. L., C. G., R. L., M. P., M. A. H., A. D. S., A. S., R. J., and A. B. analyzed the results and approved the final version of the manuscript.
Tetra- and Penta-acylated lipid A structures of Porphyromonas gingivalis LPS differentially activate TLR4-mediated NF-κB signal transduction cascade and immuno-inflammatory response in human gingival fibroblasts.
Lipopolysaccharide (LPS) from Burkholderia cepacia is more active than LPS from Pseudomonas aeruginosa and Stenotrophomonas maltophilia in stimulating tumor necrosis factor alpha from human monocytes.
The complete structure and pro-inflammatory activity of the lipooligosaccharide of the highly epidemic and virulent gram-negative bacterium Burkholderia cenocepacia ET-12 (strain J2315).
Chemical synthesis of Burkholderia lipid A modified with glycosyl phosphodiester-linked 4-amino-4-deoxy-β-L-arabinose and its immunomodulatory potential.
Biosynthesis and structure of the Burkholderia cenocepacia K56-2 lipopolysaccharide core oligosaccharide: truncation of the core oligosaccharide leads to increased binding and sensitivity to polymyxin B.
Expression of a Porphyromonas gingivalis lipid A palmitylacyltransferase in Escherichia coli yields a chimeric lipid A with altered ability to stimulate interleukin-8 secretion.
44901-4&id=gr1.jpg){kind=link}
44901-4&id=gr2.jpg){kind=link}
44901-4&id=gr3.jpg){kind=link}
44901-4&id=gr4.jpg){kind=link}
44901-4&id=gr5.jpg){kind=link}
44901-4&id=gr6.jpg){kind=link}
44901-4&id=gr7.jpg){kind=link}
44901-4&id=gr8.jpg){kind=link}
44901-4&id=gr9.jpg){kind=link}
44901-4&id=gr10.jpg){kind=link}
44901-4&id=gr11.jpg){kind=link}
44901-4&id=gr12.jpg){kind=link}
44901-4&id=gr13.jpg){kind=link}
44901-4&id=grt2.jpg){kind=link}