Next Article in Journal
High-Rate Anaerobic Digestion of Waste Activated Sludge by Integration of Electro-Fenton Process
Next Article in Special Issue
GPR6 Structural Insights: Homology Model Construction and Docking Studies
Previous Article in Journal
Characterization and Dimethyl Phthalate Flocculation Performance of the Cationic Polyacrylamide Flocculant P(AM-DMDAAC) Produced by Microwave-Assisted Synthesis
Previous Article in Special Issue
Considerations for Docking of Selective Angiotensin-Converting Enzyme Inhibitors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

TLR4-Targeting Therapeutics: Structural Basis and Computer-Aided Drug Discovery Approaches

Department of Molecular Science and Technology, Ajou University, Suwon 16499, Korea
*
Author to whom correspondence should be addressed.
Molecules 2020, 25(3), 627; https://doi.org/10.3390/molecules25030627
Submission received: 5 January 2020 / Revised: 28 January 2020 / Accepted: 29 January 2020 / Published: 31 January 2020
(This article belongs to the Special Issue Computational Methods for Drug Discovery and Design)

Abstract

:
The integration of computational techniques into drug development has led to a substantial increase in the knowledge of structural, chemical, and biological data. These techniques are useful for handling the big data generated by empirical and clinical studies. Over the last few years, computer-aided drug discovery methods such as virtual screening, pharmacophore modeling, quantitative structure-activity relationship analysis, and molecular docking have been employed by pharmaceutical companies and academic researchers for the development of pharmacologically active drugs. Toll-like receptors (TLRs) play a vital role in various inflammatory, autoimmune, and neurodegenerative disorders such as sepsis, rheumatoid arthritis, inflammatory bowel disease, Alzheimer’s disease, multiple sclerosis, cancer, and systemic lupus erythematosus. TLRs, particularly TLR4, have been identified as potential drug targets for the treatment of these diseases, and several relevant compounds are under preclinical and clinical evaluation. This review covers the reported computational studies and techniques that have provided insights into TLR4-targeting therapeutics. Furthermore, this article provides an overview of the computational methods that can benefit a broad audience in this field and help with the development of novel drugs for TLR-related disorders.

1. Introduction

Toll-like receptor 4 (TLR4) belongs to the pattern recognition receptor family, which plays a key role in the human defense mechanism and responds to invading pathogens with high selectivity and sensitivity [1,2]. TLR4 is sensitive to pathogen-associated molecular patterns (PAMPs) such as lipopolysaccharide (LPS) and lipo-oligosaccharide. Moreover, TLR4 recognizes PAMPs from fungi, viruses, and mycoplasmas [3]. In addition to PAMPs, TLR4 can be activated by certain endogenous ligands produced due to tissue injury and/or inflammation [3,4]. This receptor–ligand interaction initiates an intracellular signaling cascade that leads to the subsequent proinflammatory response [5]. Thus, due to the involvement of TLR4 in various pathological conditions, it is considered a potential therapeutic target for drug development. Many drugs targeting TLR4 are already in clinical trials, and the two Food and Drug Administration (FDA)-approved drugs named Bacillus Calmette–Guerin (a vaccine used against tuberculosis) and monophosphoryl lipid A (MPLA; used for the treatment of bladder and cervical cancer) are available on the market [6]. Nonetheless, there is a need for new drugs for life-threating diseases involving TLR4 dysregulation. Thus, the hunt for a potential drug targeting TLR4 is still in progress. Over the past decades, computer-aided drug discovery (CADD) has made tremendous progress in the development of drugs including small organic molecules, aptamers, peptides, and antibodies. CADD is a cost-effective method that has accelerated the drug development process by reducing the total timespan from target identification to hit discovery and optimization [7]. Generally, CADD is categorized into two main approaches: structure-based drug design (SBDD) [8] and ligand-based drug design (LBDD) [9,10,11]. In SBDD, the structural information on the target molecule is used to identify hotspot residues and key interactions that are crucial for their biological activity [7]. Such information can be helpful for designing new drug candidates to modulate the function of proteins involved in a disease. In the absence of protein three-dimensional (3D) structure, the LBDD method is employed, which takes into account the information on known ligands for the specific target. Structure-activity relationship (SAR) analysis is one of the applications of LBDD method. In SAR analysis, the chemical structure and biological activity of known ligands can be utilized to design new drugs or optimize known drugs to increase their efficacy [10]. The main steps in the drug discovery pipeline are target identification and validation followed by lead discovery and optimization (Figure 1). In next stage, the optimized hit is tested empirically, whereby toxicity and pharmacokinetic properties are evaluated. Afterwards, nontoxic and effective hits proceed to preclinical and clinical trials, and the most effective drug that passes the clinical trials is distributed in the market [7]. This review focuses on the activation mechanism and structural details of TLR4 and most importantly, the therapeutic importance of TLR4 along with the progress in drug development through computational techniques. Furthermore, light is shed on various computational studies that have addressed the importance of TLR4 by means of various state-of-the-art computational tools.

2. Structure and Interaction Interfaces of TLR4

TLR4 is a transmembrane receptor that belongs to the superfamily of leucine-rich repeat (LRR)-containing proteins. LRRs are 20–29 residues long sequence motifs present in various proteins and provide unique structural features for protein–protein interaction. The full-length structure of TLR4 is composed of three domains: an ectodomain, transmembrane helical domain, and intracellular Toll-interleukin-1 receptor (TIR) domain. Thus far, none of full-length crystal structures of TLR4 has been resolved, but the structures of individual domains have been reported. Crystal structure of the extracellular domain (ECD) of human TLR4 in complex with myeloid differentiation factor 2 (MD2) and LPS has been resolved [12]. The ECD consists of 608 amino acid (aa) residues and has a horseshoe-like shape comprising a concave and convex surface [13]. Parallel β-sheets collectively make up the concave surface of the ECD, while the convex surface consists of loops and 310 helices [14]. The ECD is responsible for the recognition of specific structural patterns in invading microbial molecules, where the binding of these molecules with the receptor leads to the activation and dimerization of TLR4 [5,15].
The ECD of human TLR4 (27–631 aa) is further classified into three subdomains: N-terminal domain, central domain, and C-terminal domain (Figure 2a). The N-terminal domain is a highly hydrophobic region that extends from residue 27 to residue 201 and contains LRR modules 1 through 6. The central domain contains a hypervariable region and is important for LPS recognition. It is composed of LRR modules 7 to 12. Lastly, the C-terminal domain contains LRR modules 13 to 22 [12,13]. β-Sheets of the C-terminal domain are 28% larger in radii as compared to the N-terminal domain. The N-terminal and central domains of TLR4 collectively constitute a primary contact surface comprising two evolutionarily and chemically distant regions, the A patch and B patch. These patches enable charge complementarity to its coreceptor MD2 for binding, forming a stable heterodimer [12,16]. The A patch is evolutionarily conserved and negatively charged and interacts with a positively charged MD2 surface. By contrast, patch B is positively charged and interacts with the negatively charged surface of MD2. Any mutation in the primary interface can disrupt TLR4–MD2 binding [13]. The transmembrane domain of human TLR4 (632–652 aa) forms a helical structure of 21 residues followed by loop and a juxtamembrane helix. The intracellular TIR domain is composed of 187 aa and is responsible for the initiation of the downstream signaling cascade [5,17].
MD2 has a β-cup fold structure composed of two antiparallel β-sheets separated from each other and forming a hydrophobic pocket, which is ideal for the binding of ligands such as LPS. The β-cup fold structure of MD2 is different from the immunoglobulin fold because it lacks a conserved disulfide bridge that connects the two β-sheets thus forming a large internal pocket. The internal surface of the pocket is made up of hydrophobic residues, and the outer region contains the positively charged residues that promote the LPS binding (Figure 2b) [13,18]. Electrostatic interactions, hydrogen bonds, and a few hydrophobic residues contribute to the binding of MD2 to TLR4.

3. The Signaling Mechanism of TLR4 and Its Role in Inflammatory Diseases

A release of LPS from invading bacteria is an early sign of infection and is capable of inducing a strong immunological response [16]. LPS is composed of three structural components: (I) a conserved hydrophobic lipid A component, (II) a branched hydrophilic core polysaccharide chain, and (III) a hydrophilic repeating oligosaccharide side chain [15]. The engagement of TLR4 by LPS has a complicated mechanism and involves several proteins including LPS-binding protein (LBP), cluster of differentiation 14 (CD14), and MD2. LBP extracts LPS from the bacterial membrane and then transfers it to CD14, which assists with the transfer of LPS to the TLR4–MD2 complex [15,19,20]. Upon LPS recognition, the TLR4–MD2 homo-heterodimer forms that triggers the downstream signaling cascade. This process is schematically presented in Figure 3. TLR4 has a distinctive capability of initiating two signaling cascades: the myeloid differentiation primary response protein 88 (MyD88)-dependent pathway and TIR domain-containing adaptor inducing interferon-β (TRIF)-dependent (MyD88-independent) pathway (Figure 4) [15]. In the former, TLR4-MD2 dimerization at the cell surface recruits adaptor proteins TIR domain-containing adaptor protein (TIRAP) and MyD88. MyD88 interacts with the IL-1R–associated kinase (IRAK) complex, which recruits tumor necrosis factor receptor–associated factor 6 (TRAF-6) [21,22,23]. It promotes the activation of transcription of nuclear factor-κB (NF-ĸB) genes and proinflammatory cytokines tumor necrosis factor-α (TNF-α) and interleukin-6 (IL-6) [24,25,26]. Endocytosis of TLR4-MD2 complex activates the MyD88-independent pathway in the endosomal compartment [27]. Instead of TIRAP, TRIF-related adaptor molecule (TRAM) recruits TRIF and switches TRAF3 on. TRAF3 is responsible for downstream induction of TANK-binding kinase 1 (TBK1) and IKKε leading to the activation of transcription factor IRF3 and type I interferon [28,29]. The MyD88-independent pathway stimulates the late-phase induction of a mitogen-activated protein kinase (MAPK) and NF-ĸB [30].
Generally, inflammation has a protective role in the affected area; however, an unnecessary host response to various exogenous or endogenous molecules can lead to serious health problems [31]. Clinical data shows that any dysregulation of TLR4-mediated signaling triggered by PAMPs and damage-associated molecular patterns (DAMPs) may cause autoimmune, cardiovascular, neurological, and infectious diseases [32,33,34,35,36]. TLR4 also recognizes DAMPs including high-mobility group protein 1 (HMGB1), heparin sulfate, hyaluronic acid, heat shock proteins, and oxidized phospholipids [3,4].
In addition to immune cells, tumor cells express TLR4, where it acts as a double-edged sword on tumorigenesis. Various lines of evidence have revealed that TLR4 engagement within a tumor microenvironment can promote either antitumor immunity or tumor progression [37,38,39,40,41]. It has been reported that TNF-α boosts intratumoral penetration of chemotherapeutic agents and acts as an adjuvant for chemotherapy. Moreover, as discussed elsewhere, activation of two distinct signaling pathways by TLR4 and TNF-α secretion makes TLR4 a promising target for cancer immunotherapy [42,43]. Nevertheless, a better understanding of the mechanism of action of TLR4 in cancer and in other diseases is necessary to develop potent and successful drug candidates.

4. The Role of Computational Techniques in TLR4-based Therapeutics

Computational techniques have been extensively used in drug development to find potent therapeutic compounds [44]. So far, no effective compound against TLR4-mediated diseases is available on the market. The computational strategies can help with the discovery of novel compounds with promising therapeutic effects. Various computational methods employed in the drug discovery process include virtual screening (VS), quantitative SAR analysis, pharmacophore-based screening, molecular docking, and molecular dynamics (MD) simulations.

4.1. Target Identification and Structure Determination

The initial step in the first stage of the drug development process is target identification [7]. To identify the druggable biological molecule responsible for specific disease is an important step followed by functional characterization of the target molecule [45]. The characteristics of a druggable target molecule are safe (non-toxic), efficacious, and meet clinical requirements. Various computational methods like panel docking [46], or data mining/machine learning [47] and experimental techniques such as expression profiling [48], cell-based assays, and functional screening, are available for identification and selection of potent target molecules. Target validation involves confirmation of role of target molecule in the disease and whether its activity can be modulated using a small molecule drug or biologic [49]. Once the potential target molecule is identified, next step is structure determination. 3D structure information of therapeutically important targets is crucial for the development of a potent drug. Many integrated structural biology techniques are available such as cryo-electron microscopy, nuclear magnetic resonance method, and X-ray crystallography but in the absence of solution structure, computational techniques are available. Computational methods include homology modeling, threading, and ab initio modeling. Homology modeling is more reliable approach among in silico methods as it considers the structural information of homologous proteins. The predicted structures are further validated using various validation tools such as Ramachandran plot [7].
As described above, TLR4 is a promising therapeutic target because of its involvement in the pathogenesis of various diseases. To date, 12 solution structures of human TLR4 are available in the Protein Data Bank (PDB) and are listed in Table 1. Most of the crystal structures are those of the ECD N terminus, only one crystal structure of the C terminus of the ectodomain is available, and a few structures of the full-length ECD have been reported too. Furthermore, two NMR structures of the transmembrane helix are resolved, but unfortunately full-length structure of TLR4 is not determined yet.
A study has been conducted to identify the interaction surface of the TLR4-TIR domain with the TIR domain of adaptor proteins MAL and TRAM as well as the TIR-TIR interaction for TLR4 dimerization [52]. Because a solution structure of the TIR domain is not available, a homology model was constructed computationally. The homology model was energy minimized and optimized on Web servers for further use [52]. Three putative binding interfaces (for MAL, TRAM, and TLR4 TIR) were identified by alanine scanning mutagenesis. These sites were confirmed via protein-protein docking. The two sites include residues of the BB loop, whereas a third binding site is present on the paradoxical surface. Furthermore, these results were confirmed by in vitro assays [52].
To predict the full-length structure of TLR4, including the ectodomain, transmembrane, and cytoplasmic domains, a homology modeling technique has been used [53]. Online modeling server SWISS-MODEL was utilized to construct a full-length model of the TLR4 protein. All three domains were modeled independently and were combined into the single model by establishing peptide bonds between the residues present at N and C termini of flanking domains. For this purpose, Discovery studio Visualizer was used. Next, these peptide bonds were optimized, and the full-length model was energy minimized in the Gromacs software [53].

4.2. Lead Identification and Optimization

Lead identification is the next step in stage I of the drug discovery process. For identification of novel hits, VS is a robust approach [54]. In VS, a database of millions of compounds is screened and docked into a target molecule to choose hits with better binding affinity. There are two types of VS, i.e., ligand-based VS and structure-based VS. Ligand-based VS includes SAR analysis, QSAR analysis, and pharmacophore-based screening of huge libraries of active compounds. The highly active scaffolds from these libraries are identified considering various factors, such as descriptors, a similarity index, and a consensus pharmacophore [55]. By contrast, in structure-based VS, the 3D structure of a target protein and its active binding site are analyzed. The druglike molecules are docked into the binding site of the protein via suitable docking algorithms. Then, on the basis of scoring functions, the binding affinity of the docked complexes is evaluated [7]. The top hits are next optimized by structural modifications. While the main scaffold is preserved, side chains can be replaced to increase stability and efficacy. This task can be accomplished by computational methods and synthesized experimentally. The top hits are then tested in vitro and in vivo for confirmation of the results.

4.2.1. Small-Molecule Inhibitors of TLR4: Lead Identification & Preclinical Studies

Two small-molecule compounds, C11 and C15, have been identified (using the pharmacophore-mapping approach) that act as protein-protein interaction inhibitors [56]. A pharmacophore query has been generated for both the primary and dimerization interface of TLR4-MD2 and has been used to screen the ZINC database [57] by VS and molecular docking. The docked compounds were analyzed in silico for absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties and SAR. The next step was protein-protein inhibitor profiling of the screened-out compounds, and two potent compounds, C11 and C15, were selected [56]. MD simulations were performed on both compounds to check the stability of the docked complex, which showed that C15 is more stable than C11 [56].
Joce et al. have performed screening against TLR4 and MD2 and identified two antagonists, one for each target protein [58]. T5342126 was identified as a TLR4 antagonist, whereas T6071187 was found to be an inhibitor of the interaction interface of TLR4-MD2. Their activity was tested in vitro in RAW cells and it was observed that both compounds—T5342126 (at 2 µM) and T6071187 (at 200 nM)—inhibit LPS-induced TLR4 signaling [58]. Recently, in another study, derivatives of T5342126 and carvedilol [59] were tested for TLR4 inhibition [60]. These compounds were subjected to SAR analyses, and thus, novel derivatives of these two compounds were designed and synthesized. By in vitro assays, one of the derivatives, 8a, was identified as the most potent inhibitor [60] and chosen for further studies. This compound significantly inhibits inflammatory cytokines (nitric oxide and TNF-α) and downregulates TLR4 expression. Molecular docking was performed to investigate the binding mode of 8a, thereby revealing that the benzyl ring of this compound binds with the hydrophobic pocket of TLR4, and the carbazole group occupies the MD2 pocket [60].
VS has been performed on a library of 130,000 compounds obtained from the ZINC database [57], and ZINC02157367 was selected as a hit according to binding affinity. The structure of this lead molecule was optimized next to design a novel and effective compound. Considering Lipinski’s rule of five and hydrophobic active-site residues of TLR4, 12 compounds were designed [61]. Flexible docking was performed to rank the ligands by their binding energies followed by MD simulations using the all-atom optimized potentials for liquid simulations (OPLSAA) force field in Gromacs to evaluate the binding mode. The MD simulations indicated the stability of the receptor-ligand complex, and compounds 24, 14, and 17 showed the highest binding energies. To confirm these in silico results, all 12 compounds were evaluated in experimental studies, and the binding affinity of the most effective compound (17) was confirmed by surface plasmon resonance analysis [61].
In another study, five lipid A analogs (FP13-17) were designed rationally, and their binding modes were analyzed against TLR4-MD2 through molecular docking [62]. The docking analyses revealed that all the compounds fit into the MD2 pocket in an antagonist manner. To confirm the results, the docked complexes were subjected to MD simulations, and the findings indicated that all the compounds were antagonists. Later, experimental analyses were conducted to check the binding and biological activity of these compounds, and the results corroborated the antagonistic behavior of all five compounds [62].

4.2.2. TLR4-Inhibitory Peptides: Lead Identification & Preclinical Studies

In a recent study [63], computational techniques were utilized to enhance the stability and activity of a TLR4-inhibitory peptide identified through phage display. The PIP2 peptide inhibits TLR4-induced proinflammatory cytokines. PIP2 was cyclized via introduction of a β-lactam bridge using in silico tools; this change improved the half-maximal inhibitory concentration (IC50) threefold. The optimized peptide was next tested in vitro and in vivo and caused a significant reduction of inflammation in a rat model of rheumatoid arthritis. MD simulations were performed to confirm the stability of both peptides and to understand the binding mechanism [63].
Using an ab initio computational approach, novel cyclic and linear peptides for TLR4 coreceptors CD14 and MD2 have been designed [64]. These peptides were designed based on 3D structures of the target protein in ab initio design software. Schrodinger and the PyMOL molecular graphic system [65] were used for structural visualization comparison. The missing loops were modeled by the homology modeling technique. Next, cyclic and linear peptides were designed with in silico methods and synthesized in the laboratory [64]. Experimental results indicated that the cyclic peptides are rigid as compared to the linear versions. Peptides with activity values in the micromolar range were able to trigger the TLR4-mediated signaling pathway and increase IL-1β production. These agonistic peptides are patented under patent number WO2017141248 [64].

4.2.3. Miscellaneous TLR4-Related Agents

Neoseptins, a new class of TLR4 agonists, have been discovered by Morin et al. [66]. Among 90,000 compounds, neoseptin-1 was identified, and SAR analyses were performed to improve its design, thus leading to the neoseptin-3 peptidomimetic. On the other hand, the activity of this peptidomimetic is species specific because it failed to show any effect on human TLR4; that is why this compound did not proceed to clinical development [66].
Aptamers are single-stranded DNA or RNA molecules having high specificity, binding affinity, and the best characteristics of antibodies and small molecules [67]. Two TLR4-blocking DNA aptamers, ApTLR#1R and ApTLR#4F, were reported recently [68]. These aptamers were picked form a random oligonucleotide pool and optimized to enhance the affinity. ApTLR#4F and its truncated form, ApTLR#4FT, were found to be highly effective and exerted long-lasting effects in patients with a stroke [68]. Aptamers are reported to be less toxic as compared to most of small-molecule drugs, thus holding promise for clinical use.
In another study, it was demonstrated that heparin-based nanoparticles bind to TLR4-MD2 and inhibit NF-ĸB activation in macrophages [69]. With the help of computational techniques, the binding mode of the nanoparticle and TLR4-MD2 complex was analyzed, and the interacting residues of TLR4 and MD2 were highlighted [69].
A machine-learning Cox proportional hazards regression model has been employed to analyze the influence of disease marker genes on the survival of patients with ovarian cancer. In this regression analysis involving 41 genes, the expression of five genes (TLR4, BSCL2, CDH1, ERBB2, and SCGB2A1) accounted for significant variation in the patients’ survival rate. This study suggests that TLR4 along with the other four genes is a potential drug target in ovarian cancer and that their expression is related to patient survival [70].

4.3. MD Simulations of TLR4

MD simulations were performed on the TLR4-MD2 complex with a bound natural ligand, ursolic acid (URA), which interferes the LPS binding [71]. URA is a lipophilic five ringed structure and plant-based natural compound. This study uncovered the important residues (Ile52, Leu54, Leu78, Ile80, Val82, Phe119, Phe121, Tyr131, and Cys133) for the interaction of URA and TLR4-MD2 on the basis of binding energy calculations and energy decomposition analysis. The binding mode of the inhibitor with the TLR4-MD2 complex was studied too. The diameter of the TLR4-MD2-URA complex after 150 ns MD simulations was estimated. The average diameter of the TLR4-MD2 complex was 4.43 nm, whereas in the presence of URA, the diameter diminished to an average value of 3.46 nm [71].
Recently, MD simulations were carried out to gain insight into the activation mechanism of TLR4-MD2 mouse protein structure. A 1.2 µs simulation was performed four times on four different complexes (TLR4-MD2 heterodimer, TLR4-MD2 homo-heterodimer, LPS-TLR4-MD2 homo-heterodimer, and neoseptin-3-TLR4-MD2 homo-heterodimer) to verify the stability of the complexes along with binding energy calculations [72]. The results showed stable interfaces and well-maintained structure of TLR4. Using molecular mechanics, Poisson–Boltzmann surface area (MM-PBSA) key residues were identified that play a crucial part in the dimerization and intracellular signaling of TLR4 [72].
Rhodobacter sphaeroides lipid A (RsLA)-induced TLR4-MD2 signaling has been studied by computational methods in different species (humans, horses, murine, and hamsters) [73]. MD simulations revealed that the RsLA backbone acquired an antagonist-like orientation in the murine and human TLR4-MD2 complex and inhibits downstream signaling. By contrast, it activates the TLR4 pathway by acquiring an agonist-like conformation in the hamster and horse complexes [73]. This dual behavior of LA is due to the binding orientation. During simulations, acyl chains in horse and hamster complexes folded back due to increased shift in the molecule. In addition, the spatial arrangement of G1cN1-G1cN2 in RsLA resembles lipid IVa in murine and human complexes. This structural change is responsible for the specific ligand behavior. Moreover, the stability of the Phe126 loop in MD2 was assessed, which is crucial for the activation of the TLR-MD2 complex. It was noted that this loop is stable in hamsters and horses as compared to murine and humans. These results provide a convincing explanation for the species-specific behavior of RsLA [73]. The importance of the Phe126 loop was reported in another study too [74]. It was demonstrated there that morphine and morphine-3-glucuronide each interacts with MD2 near this loop, thereby forming a complex, and its stability increases when it interacts with TLR4.
Recently, an interaction between HMGB1 and the TLR4-MD2 complex was studied by molecular docking and MD simulations [75]. In this study, crystal structure (PDB ID: 3FXI) of the TLR4-ECD was used. HMGB1 consist of 215 residues divided into two DNA-binding domains termed as A-box and B-box and a C-terminal domain. The cysteine residues present in the DNA-binding domains makes a disulfide bridge and induces structural changes [76]. Protein-protein docking of HMGB1, MD2, and TLR4 was performed on docking server ZDOCK [77]. MD simulations were conducted to characterize the behavior of full-length HMGB1, docked complexes of TLR4, and mutants by means of the OPLS force field. Mutagenesis and surface plasmon resonance analyses were conducted to study the interactions. Their data revealed that the N terminus of TLR4 binds to HMGB1 A-box but does not help with dimer formation, thereby preventing the launch of downstream signaling and HMGB1-induced inflammation. Meanwhile, the B-box fragment of HMGB1 promotes the TLR4 dimerization, which results in the activation of the downstream proinflammatory signaling cascade and cytokine production [75].

5. The Molecular Understanding of TLR4-targeting Drugs

Given the role of TLR4 in the pathogenesis of many diseases, several therapeutic modulators have been devised to regulate TLR4 expression, and a few of these compounds are currently being evaluated in clinical trials (Table 2). These compounds can be categorized as antibodies, small-molecule inhibitors, peptides, microRNAs, nanoparticles, lipid A analogs, and derivatives of natural products. Detailed information about TLR4 modulators can be found in the literature [1,60,78,79,80]. Most TLR4 agonists belong to the class of glycolipids that are LPS mimetic. LPS contains lipid A, which is a membrane-anchoring and biologically active moiety. It is composed of six fatty acid chains, two negatively charged phosphate esters, and a glucosamine disaccharide core [81]. Five of the six fatty acid chains are completely submerged into the hydrophobic pocket of the MD2 protein. Meanwhile, the other chain participates in a hydrophobic interaction with the other domains of the TLR4 molecule.
Based on the binding mode of LPS to the TLR4-MD2 heterodimer, several strategies have been implemented to chemically modify lipid A structure, e.g., the removal of phosphate groups or a change in the position, number, or length of carbon chains or in the glucosamine backbone. These modifications can switch the activity of lipid A from an agonist to antagonist mode. The compounds derived from lipid A structure are eritoran [82], aminoglycosides [83], MPLA [84], and lipid A mimetics (Figure 5). Eritoran is a well-known TLR4 antagonist that directly binds to the hydrophobic pocket of MD2 and blocks TLR4 dimerization [13]. Nevertheless, eritoran failed in a phase III clinical trial because it yielded poor outcomes in severe sepsis cases [85]. MPLA is a TLR4 agonist (extracted from Salmonella Minnesota), binds to TLR4-MD2, and has been used as a vaccine adjuvant against cervical cancer and hepatitis B [86]. Other well-known lipid A analogs include lipid IVa [87] and glucopyranosyl lipid adjuvant (GLA) [88]. Lipid IVa has a species-specific activity because it acts as an agonist in mice, whereas in humans, it behaves as an inhibitor [87]. GLA is a synthetic TLR4 agonist that is composed of a phosphate group, six acyl chains, and a disaccharide backbone. It is an effective adjuvant and enhances an immune response [89,90].
The inhibition of TLR4 by a small-molecule can be achieved by blocking (1) ligand-receptor interaction, (2) dimerization of the TLR4-MD2 complex, or (3) downstream signaling. Small-molecule antagonist TAK-242 (resatorvid) (Figure 5), a cyclohexane derivative binds to Cys747 in the intracellular domain of TLR4 and inhibits its interaction with adaptor molecules [91]. A cyclohexane ring present in TAK-242 interacts with the −SH2 group of cysteine of TIR domain. Nonetheless, this compound failed in phase III clinical trials because it did not lower the IL-6 level in sepsis or respiratory failure [92]. Various studies have shown that curcumin has anti-inflammatory properties and downregulates pro-inflammatory cytokines such as IL-1, -2, -6, -8, and TNF-α as well as MAPKs [93]. Curcumin (Figure 5) can inhibit both pathways of TLR4 i.e., MyD88 dependent and TRIF dependent [94]—and can fit into the MD2 pocket as confirmed by experimental and computational studies [95,96,97]. To enhance the stability of curcumin, its analogs have been synthesized, and among these analogs, L48H37 (Figure 5) features greater chemical stability and enhanced anti-inflammatory properties [98]. These compounds contain a conserved group 3-(4-hydroxyphenyl) acrylaldehyde, which directly binds to the MD2 pocket and blocks the LPS binding. Therapeutic applications of curcumin in a wide range of inflammatory diseases have been reported [99]; the same is true for metabolic syndrome [100], type II diabetes mellitus [101], systemic lupus erythematosus [102], melanoma [103], pulmonary diseases, and lung cancer [104]. Another well-known plant-derived inhibitor of TLR4 is sulforaphane (Figure 5), which inhibits TLR4 dimerization and has an anti-inflammatory effect [105]. It has been reported that sulforaphane inhibits LPS-induced secretion of cyclooxygenase 2, TNF-α, HMGB1, and inducible NO synthase in macrophages [106,107]. Sulforaphane binds to the cysteine residue located in the MD2 pocket at position 133, which performs an important function in the MD2-LPS interaction. Sulforaphane occupies the binding pocket and as a result disrupts the binding of LPS to the TLR4-MD2 complex [108].
In the last decade, developments in the field of therapeutic monoclonal antibodies were incredible (in terms of the number and impact) due to their high specificity and ability to treat a variety of medical conditions. NI-0101 is a humanized monoclonal antibody that interferes with TLR4 dimerization and eventually blocks its activation. NI-0101 has strong anti-inflammatory properties and suppresses cytokine production [109]. Currently, NI-0101 is being tested in phase II clinical trials against rheumatoid arthritis (Table 2).

6. Conclusions

TLR4 plays a key role in the host defense mechanism and is capable of recognizing a variety of PAMPs and DAMPs. On the other hand, any dysregulation of the TLR4 pathway can initiate various diseases. Many TLR4-targeting drugs are being assessed in clinical trials, and a few of them have failed in phase III. CADD has emerged as a robust alternative approach to identify novel and potent TLR4 modulators. Computational techniques are cost effective, efficient and can be used to optimize a lead compound. Further, computational techniques such as protein modeling, protein-protein docking, and MD simulations have greatly contributed to the understanding of poorly studied crucial interfaces that may be helpful for drug discovery. Recent advances in the pharmaceutical industry are facilitating the development of safe and potent drugs.

Author Contributions

Conceptualization, Q.u.A., M.B., and S.C.; Investigation and Formal Analysis, Q.u.A. and M.B.; Visualization and Writing—Original Draft Preparation, Q.u.A., M.B., and S.C.; Supervision and Project Administration, S.C. All authors have read and agreed to the final version of the manuscript.

Funding

This research was funded by the National Research Foundation of Korea, grant numbers NRF-2019M3A9A8065098, 2019M3D1A1078940 and 2019R1A6A1A11051471.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, L.; Fu, W.; Zheng, L.; Wang, Y.; Liang, G. Recent progress in the discovery of myeloid differentiation 2 (MD2) modulators for inflammatory diseases. Drug Discov. Today 2018, 23, 1187–1202. [Google Scholar] [CrossRef]
  2. Kawai, T.; Akira, S. TLR signaling. Cell Death Differ. 2006, 13, 816–825. [Google Scholar] [CrossRef] [Green Version]
  3. Lucas, K.; Maes, M. Role of the Toll Like receptor (TLR) radical cycle in chronic inflammation: Possible treatments targeting the TLR4 pathway. Mol. Neurobiol. 2013, 48, 190–204. [Google Scholar] [CrossRef] [PubMed]
  4. Imai, Y.; Kuba, K.; Neely, G.G.; Yaghubian-Malhami, R.; Perkmann, T.; van Loo, G.; Ermolaeva, M.; Veldhuizen, R.; Leung, Y.H.; Wang, H.; et al. Identification of oxidative stress and Toll-like receptor 4 signaling as a key pathway of acute lung injury. Cell 2008, 133, 235–249. [Google Scholar] [CrossRef] [PubMed]
  5. Akira, S.; Takeda, K. Toll-like receptor signalling. Nat. Rev. Immunol. 2004, 4, 499–511. [Google Scholar] [CrossRef] [PubMed]
  6. Vacchelli, E.; Galluzzi, L.; Eggermont, A.; Fridman, W.H.; Galon, J.; Sautes-Fridman, C.; Tartour, E.; Zitvogel, L.; Kroemer, G. Trial watch: FDA-approved Toll-like receptor agonists for cancer therapy. Oncoimmunology 2012, 1, 894–907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Batool, M.; Ahmad, B.; Choi, S. A Structure-Based Drug Discovery Paradigm. Int. J. Mol. Sci. 2019, 20, 2783. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Kalyaanamoorthy, S.; Chen, Y.P. Structure-based drug design to augment hit discovery. Drug Discov. Today 2011, 16, 831–839. [Google Scholar] [CrossRef]
  9. Yu, W.; MacKerell, A.D., Jr. Computer-Aided Drug Design Methods. Methods Mol. Biol. 2017, 1520, 85–106. [Google Scholar]
  10. Ou-Yang, S.S.; Lu, J.Y.; Kong, X.Q.; Liang, Z.J.; Luo, C.; Jiang, H. Computational drug discovery. Acta Pharmacol. Sin. 2012, 33, 1131–1140. [Google Scholar] [CrossRef] [Green Version]
  11. Drwal, M.N.; Griffith, R. Combination of ligand- and structure-based methods in virtual screening. Drug Discov. Today Technol. 2013, 10, e395–e401. [Google Scholar] [CrossRef] [PubMed]
  12. Park, B.S.; Song, D.H.; Kim, H.M.; Choi, B.S.; Lee, H.; Lee, J.O. The structural basis of lipopolysaccharide recognition by the TLR4-MD-2 complex. Nature 2009, 458, 1191–1195. [Google Scholar] [CrossRef] [PubMed]
  13. Kim, H.M.; Park, B.S.; Kim, J.I.; Kim, S.E.; Lee, J.; Oh, S.C.; Enkhbayar, P.; Matsushima, N.; Lee, H.; Yoo, O.J.; et al. Crystal structure of the TLR4-MD-2 complex with bound endotoxin antagonist Eritoran. Cell 2007, 130, 906–917. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kobe, B.; Kajava, A.V. The leucine-rich repeat as a protein recognition motif. Curr. Opin. Struct. Biol. 2001, 11, 725–732. [Google Scholar] [CrossRef]
  15. Maeshima, N.; Fernandez, R.C. Recognition of lipid A variants by the TLR4-MD-2 receptor complex. Front. Cell Infect. Microbiol. 2013, 3, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Park, B.S.; Lee, J.O. Recognition of lipopolysaccharide pattern by TLR4 complexes. Exp. Mol. Med. 2013, 45, e66. [Google Scholar] [CrossRef] [Green Version]
  17. Mineev, K.S.; Goncharuk, S.A.; Goncharuk, M.V.; Volynsky, P.E.; Novikova, E.V.; Aresinev, A.S. Spatial structure of TLR4 transmembrane domain in bicelles provides the insight into the receptor activation mechanism. Sci. Rep. 2017, 7, 6864. [Google Scholar] [CrossRef] [Green Version]
  18. Ohto, U.; Fukase, K.; Miyake, K.; Satow, Y. Crystal structures of human MD-2 and its complex with antiendotoxic lipid IVa. Science 2007, 316, 1632–1634. [Google Scholar] [CrossRef]
  19. Lu, Y.C.; Yeh, W.C.; Ohashi, P.S. LPS/TLR4 signal transduction pathway. Cytokine 2008, 42, 145–151. [Google Scholar] [CrossRef]
  20. Gay, N.J.; Symmons, M.F.; Gangloff, M.; Bryant, C.E. Assembly and localization of Toll-like receptor signalling complexes. Nat. Rev. Immunol. 2014, 14, 546–558. [Google Scholar] [CrossRef]
  21. Ferrao, R.; Zhou, H.; Shan, Y.; Liu, Q.; Li, Q.; Shaw, D.E.; Li, X.; Wu, H. IRAK4 dimerization and trans-autophosphorylation are induced by Myddosome assembly. Mol. Cell 2014, 55, 891–903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Barton, G.M.; Medzhitov, R. Toll-like receptor signaling pathways. Science 2003, 300, 1524–1525. [Google Scholar] [CrossRef] [PubMed]
  23. Batool, M.; Anwar, M.A.; Choi, S. Toll-like receptors targeting technology for the treatment of lymphoma. Expert. Opin. Drug Discov. 2016, 11, 1047–1059. [Google Scholar] [CrossRef] [PubMed]
  24. Oosenbrug, T.; van de Graaff, M.J.; Ressing, M.E.; van Kasteren, S.I. Chemical Tools for Studying TLR Signaling Dynamics. Cell Chem. Biol. 2017, 24, 801–812. [Google Scholar] [CrossRef] [Green Version]
  25. Kagan, J.C.; Medzhitov, R. Phosphoinositide-mediated adaptor recruitment controls Toll-like receptor signaling. Cell 2006, 125, 943–955. [Google Scholar] [CrossRef] [Green Version]
  26. Takeuchi, O.; Akira, S. Pattern recognition receptors and inflammation. Cell 2010, 140, 805–820. [Google Scholar] [CrossRef] [Green Version]
  27. Kagan, J.C.; Su, T.; Horng, T.; Chow, A.; Akira, S.; Medzhitov, R. TRAM couples endocytosis of Toll-like receptor 4 to the induction of interferon-beta. Nat. Immunol. 2008, 9, 361–368. [Google Scholar] [CrossRef] [Green Version]
  28. Gomez, R.; Villalvilla, A.; Largo, R.; Gualillo, O.; Herrero-Beaumont, G. TLR4 signalling in osteoarthritis--finding targets for candidate DMOADs. Nat. Rev. Rheumatol. 2015, 11, 159–170. [Google Scholar] [CrossRef]
  29. Molteni, M.; Bosi, A.; Rossetti, C. Natural Products with Toll-Like Receptor 4 Antagonist Activity. Int. J. Inflam. 2018, 2018, 2859135. [Google Scholar] [CrossRef] [Green Version]
  30. Yamamoto, M.; Sato, S.; Hemmi, H.; Hoshino, K.; Kaisho, T.; Sanjo, H.; Takeuchi, O.; Sugiyama, M.; Okabe, M.; Takeda, K.; et al. Role of adaptor TRIF in the MyD88-independent toll-like receptor signaling pathway. Science 2003, 301, 640–643. [Google Scholar] [CrossRef]
  31. Seymour, C.W.; Liu, V.X.; Iwashyna, T.J.; Brunkhorst, F.M.; Rea, T.D.; Scherag, A.; Rubenfeld, G.; Kahn, J.M.; Shankar-Hari, M.; Singer, M.; et al. Assessment of Clinical Criteria for Sepsis: For the Third International Consensus Definitions for Sepsis and Septic Shock (Sepsis-3). JAMA 2016, 315, 762–774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Abdollahi-Roodsaz, S.; Joosten, L.A.; Roelofs, M.F.; Radstake, T.R.; Matera, G.; Popa, C.; van der Meer, J.W.; Netea, M.G.; van den Berg, W.B. Inhibition of Toll-like receptor 4 breaks the inflammatory loop in autoimmune destructive arthritis. Arthritis Rheum. 2007, 56, 2957–2967. [Google Scholar] [CrossRef]
  33. Erridge, C. The roles of Toll-like receptors in atherosclerosis. J. Innate. Immun 2009, 1, 340–349. [Google Scholar] [CrossRef] [PubMed]
  34. Panaro, M.A.; Lofrumento, D.D.; Saponaro, C.; De Nuccio, F.; Cianciulli, A.; Mitolo, V.; Nicolardi, G. Expression of TLR4 and CD14 in the central nervous system (CNS) in a MPTP mouse model of Parkinson’s-like disease. Immunopharmacol. Immunotoxicol. 2008, 30, 729–740. [Google Scholar] [CrossRef]
  35. De Paola, M.; Mariani, A.; Bigini, P.; Peviani, M.; Ferrara, G.; Molteni, M.; Gemma, S.; Veglianese, P.; Castellaneta, V.; Boldrin, V.; et al. Neuroprotective effects of toll-like receptor 4 antagonism in spinal cord cultures and in a mouse model of motor neuron degeneration. Mol. Med. 2012, 18, 971–981. [Google Scholar] [CrossRef] [PubMed]
  36. Tsujimoto, H.; Ono, S.; Efron, P.A.; Scumpia, P.O.; Moldawer, L.L.; Mochizuki, H. Role of Toll-like receptors in the development of sepsis. Shock 2008, 29, 315–321. [Google Scholar] [CrossRef] [Green Version]
  37. Dajon, M.; Iribarren, K.; Cremer, I. Toll-like receptor stimulation in cancer: A pro- and anti-tumor double-edged sword. Immunobiology 2017, 222, 89–100. [Google Scholar] [CrossRef] [PubMed]
  38. Molteni, M.; Marabella, D.; Orlandi, C.; Rossetti, C. Melanoma cell lines are responsive in vitro to lipopolysaccharide and express TLR-4. Cancer Lett. 2006, 235, 75–83. [Google Scholar] [CrossRef]
  39. Adams, S. Toll-like receptor agonists in cancer therapy. Immunotherapy 2009, 1, 949–964. [Google Scholar] [CrossRef] [Green Version]
  40. Drexler, S.K.; Foxwell, B.M. The role of toll-like receptors in chronic inflammation. Int. J. Biochem. Cell Biol. 2010, 42, 506–518. [Google Scholar] [CrossRef]
  41. Shetab Boushehri, M.A.; Lamprecht, A. TLR4-Based Immunotherapeutics in Cancer: A Review of the Achievements and Shortcomings. Mol. Pharm. 2018, 15, 4777–4800. [Google Scholar] [CrossRef] [PubMed]
  42. Garay, R.P.; Viens, P.; Bauer, J.; Normier, G.; Bardou, M.; Jeannin, J.F.; Chiavaroli, C. Cancer relapse under chemotherapy: Why TLR2/4 receptor agonists can help. Eur. J. Pharmacol. 2007, 563, 1–17. [Google Scholar] [CrossRef] [PubMed]
  43. Xiao, L.; Kim, J.; Lim, M.; Dai, B.; Yang, L.; Reed, S.G.; Baltimore, D.; Wang, P. A TLR4 agonist synergizes with dendritic cell-directed lentiviral vectors for inducing antigen-specific immune responses. Vaccine 2012, 30, 2570–2581. [Google Scholar] [CrossRef] [Green Version]
  44. Sliwoski, G.; Kothiwale, S.; Meiler, J.; Lowe, E.W., Jr. Computational methods in drug discovery. Pharmacol. Rev. 2014, 66, 334–395. [Google Scholar] [CrossRef] [Green Version]
  45. Katsila, T.; Spyroulias, G.A.; Patrinos, G.P.; Matsoukas, M.T. Computational approaches in target identification and drug discovery. Comput. Struct. Biotechnol. J. 2016, 14, 177–184. [Google Scholar] [CrossRef] [Green Version]
  46. Li, H.; Gao, Z.; Kang, L.; Zhang, H.; Yang, K.; Yu, K.; Luo, X.; Zhu, W.; Chen, K.; Shen, J.; et al. TarFisDock: A web server for identifying drug targets with docking approach. Nucleic Acids Res. 2006, 34, W219–W224. [Google Scholar] [CrossRef] [Green Version]
  47. Nidhi; Glick, M.; Davies, J.W.; Jenkins, J.L. Prediction of biological targets for compounds using multiple-category Bayesian models trained on chemogenomics databases. J. Chem. Inf. Model. 2006, 46, 1124–1133. [Google Scholar] [CrossRef]
  48. Smith, S.M.; Murray, D.W. An overview of microRNA methods: Expression profiling and target identification. Methods Mol. Biol. 2012, 823, 119–138. [Google Scholar]
  49. Hughes, J.P.; Rees, S.; Kalindjian, S.B.; Philpott, K.L. Principles of early drug discovery. Br. J. Pharmacol. 2011, 162, 1239–1249. [Google Scholar] [CrossRef] [Green Version]
  50. Han, J.; Kim, H.J.; Lee, S.C.; Hong, S.; Park, K.; Jeon, Y.H.; Kim, D.; Cheong, H.K.; Kim, H.S. Structure-based rational design of a Toll-like receptor 4 (TLR4) decoy receptor with high binding affinity for a target protein. PLoS ONE 2012, 7, e30929. [Google Scholar] [CrossRef]
  51. Ohto, U.; Yamakawa, N.; Akashi-Takamura, S.; Miyake, K.; Shimizu, T. Structural analyses of human Toll-like receptor 4 polymorphisms D299G and T399I. J. Biol. Chem. 2012, 287, 40611–40617. [Google Scholar] [CrossRef] [Green Version]
  52. Bovijn, C.; Ulrichts, P.; De Smet, A.S.; Catteeuw, D.; Beyaert, R.; Tavernier, J.; Peelman, F. Identification of interaction sites for dimerization and adapter recruitment in Toll/interleukin-1 receptor (TIR) domain of Toll-like receptor 4. J. Biol. Chem. 2012, 287, 4088–4098. [Google Scholar] [CrossRef] [Green Version]
  53. Patra, M.C.; Kwon, H.K.; Batool, M.; Choi, S. Computational Insight Into the Structural Organization of Full-Length Toll-Like Receptor 4 Dimer in a Model Phospholipid Bilayer. Front. Immunol. 2018, 9, 489. [Google Scholar] [CrossRef] [Green Version]
  54. Lavecchia, A.; Di Giovanni, C. Virtual screening strategies in drug discovery: A critical review. Curr. Med. Chem. 2013, 20, 2839–2860. [Google Scholar] [CrossRef]
  55. Pedretti, A.; Mazzolari, A.; Gervasoni, S.; Vistoli, G. Rescoring and Linearly Combining: A Highly Effective Consensus Strategy for Virtual Screening Campaigns. Int. J. Mol. Sci. 2019, 20, 2060. [Google Scholar] [CrossRef] [Green Version]
  56. Mishra, V.; Pathak, C. Structural insights into pharmacophore-assisted in silico identification of protein-protein interaction inhibitors for inhibition of human toll-like receptor 4 - myeloid differentiation factor-2 (hTLR4-MD-2) complex. J. Biomol. Struct. Dyn. 2019, 37, 1968–1991. [Google Scholar] [CrossRef]
  57. Sterling, T.; Irwin, J.J. ZINC 15--Ligand Discovery for Everyone. J. Chem. Inf. Model. 2015, 55, 2324–2337. [Google Scholar] [CrossRef]
  58. Joce, C.; Stahl, J.A.; Shridhar, M.; Hutchinson, M.R.; Watkins, L.R.; Fedichev, P.O.; Yin, H. Application of a novel in silico high-throughput screen to identify selective inhibitors for protein-protein interactions. Bioorg Med. Chem. Lett. 2010, 20, 5411–5413. [Google Scholar] [CrossRef] [Green Version]
  59. Amirshahrokhi, K.; Khalili, A.R. Carvedilol attenuates paraquat-induced lung injury by inhibition of proinflammatory cytokines, chemokine MCP-1, NF-kappaB activation and oxidative stress mediators. Cytokine 2016, 88, 144–153. [Google Scholar] [CrossRef]
  60. Xu, Y.; Chen, S.; Cao, Y.; Zhou, P.; Chen, Z.; Cheng, K. Discovery of novel small molecule TLR4 inhibitors as potent anti-inflammatory agents. Eur. J. Med. Chem. 2018, 154, 253–266. [Google Scholar] [CrossRef]
  61. Honegr, J.; Dolezal, R.; Malinak, D.; Benkova, M.; Soukup, O.; Almeida, J.; Franca, T.C.C.; Kuca, K.; Prymula, R. Rational Design of a New Class of Toll-Like Receptor 4 (TLR4) Tryptamine Related Agonists by Means of the Structure- and Ligand-Based Virtual Screening for Vaccine Adjuvant Discovery. Molecules 2018, 23, 102. [Google Scholar] [CrossRef] [Green Version]
  62. Cochet, F.; Facchini, F.A.; Zaffaroni, L.; Billod, J.M.; Coelho, H.; Holgado, A.; Braun, H.; Beyaert, R.; Jerala, R.; Jimenez-Barbero, J.; et al. Novel carboxylate-based glycolipids: TLR4 antagonism, MD-2 binding and self-assembly properties. Sci. Rep. 2019, 9, 919. [Google Scholar] [CrossRef] [Green Version]
  63. Achek, A.; Shah, M.; Seo, J.Y.; Kwon, H.K.; Gui, X.; Shin, H.J.; Cho, E.Y.; Lee, B.S.; Kim, D.J.; Lee, S.H.; et al. Linear and Rationally Designed Stapled Peptides Abrogate TLR4 Pathway and Relieve Inflammatory Symptoms in Rheumatoid Arthritis Rat Model. J. Med. Chem. 2019, 62, 6495–6511. [Google Scholar] [CrossRef]
  64. Michaeli, A.; Mezan, S.; Kuhbacher, A.; Finkelmeier, D.; Elias, M.; Zatsepin, M.; Reed, S.G.; Duthie, M.S.; Rupp, S.; Lerner, I.; et al. Computationally Designed Bispecific MD2/CD14 Binding Peptides Show TLR4 Agonist Activity. J. Immunol. 2018, 201, 3383–3391. [Google Scholar] [CrossRef]
  65. Lill, M.A.; Danielson, M.L. Computer-aided drug design platform using PyMOL. J. Comput. Aided. Mol. Des. 2011, 25, 13–19. [Google Scholar] [CrossRef]
  66. Morin, M.D.; Wang, Y.; Jones, B.T.; Su, L.; Surakattula, M.M.; Berger, M.; Huang, H.; Beutler, E.K.; Zhang, H.; Beutler, B.; et al. Discovery and Structure-Activity Relationships of the Neoseptins: A New Class of Toll-like Receptor-4 (TLR4) Agonists. J. Med. Chem. 2016, 59, 4812–4830. [Google Scholar] [CrossRef] [Green Version]
  67. Keefe, A.D.; Pai, S.; Ellington, A. Aptamers as therapeutics. Nat. Rev. Drug Discov. 2010, 9, 537–550. [Google Scholar] [CrossRef]
  68. Fernandez, G.; Moraga, A.; Cuartero, M.I.; Garcia-Culebras, A.; Pena-Martinez, C.; Pradillo, J.M.; Hernandez-Jimenez, M.; Sacristan, S.; Ayuso, M.I.; Gonzalo-Gobernado, R.; et al. TLR4-Binding DNA Aptamers Show a Protective Effect against Acute Stroke in Animal Models. Mol. Ther. 2018, 26, 2047–2059. [Google Scholar] [CrossRef] [Green Version]
  69. Babazada, H.; Yanamoto, S.; Hashida, M.; Yamashita, F. Binding and structure-kinetic relationship analysis of selective TLR4-targeted immunosuppressive self-assembling heparin nanoparticles. Int. J. Pharm. 2018, 552, 76–83. [Google Scholar] [CrossRef]
  70. Hossain, M.A.; Saiful Islam, S.M.; Quinn, J.M.W.; Huq, F.; Moni, M.A. Machine learning and bioinformatics models to identify gene expression patterns of ovarian cancer associated with disease progression and mortality. J. Biomed. Inform. 2019, 100, 103313. [Google Scholar] [CrossRef]
  71. Niu, X.; Yu, Y.; Guo, H.; Yang, Y.; Wang, G.; Sun, L.; Gao, Y.; Yu, Z.; Wang, H. Molecular modeling reveals the inhibition mechanism and binding mode of ursolic acid to TLR4-MD2. Comput. Theor. Chem. 2018, 1123, 73–78. [Google Scholar] [CrossRef]
  72. Tafazzol, A.; Duan, Y. Key residues in TLR4-MD2 tetramer formation identified by free energy simulations. PLoS Comput. Biol. 2019, 15, e1007228. [Google Scholar] [CrossRef]
  73. Anwar, M.A.; Panneerselvam, S.; Shah, M.; Choi, S. Insights into the species-specific TLR4 signaling mechanism in response to Rhodobacter sphaeroides lipid A detection. Sci. Rep. 2015, 5, 7657. [Google Scholar] [CrossRef] [Green Version]
  74. Shah, M.; Anwar, M.A.; Yesudhas, D.; Krishnan, J.; Choi, S. A structural insight into the negative effects of opioids in analgesia by modulating the TLR4 signaling: An in silico approach. Sci. Rep. 2016, 6, 39271. [Google Scholar] [CrossRef]
  75. Sun, S.; He, M.; VanPatten, S.; Al-Abed, Y. Mechanistic insights into high mobility group box-1 (HMGb1)-induced Toll-like receptor 4 (TLR4) dimer formation. J. Biomol. Struct. Dyn. 2019, 37, 3721–3730. [Google Scholar] [CrossRef]
  76. He, M.; Bianchi, M.E.; Coleman, T.R.; Tracey, K.J.; Al-Abed, Y. Exploring the biological functional mechanism of the HMGB1/TLR4/MD-2 complex by surface plasmon resonance. Mol. Med. 2018, 24, 21. [Google Scholar] [CrossRef]
  77. Pierce, B.G.; Wiehe, K.; Hwang, H.; Kim, B.H.; Vreven, T.; Weng, Z. ZDOCK server: Interactive docking prediction of protein-protein complexes and symmetric multimers. Bioinformatics 2014, 30, 1771–1773. [Google Scholar] [CrossRef]
  78. Zaffaroni, L.; Peri, F. Recent advances on Toll-like receptor 4 modulation: New therapeutic perspectives. Future Med. Chem. 2018, 10, 461–476. [Google Scholar] [CrossRef] [Green Version]
  79. Kuzmich, N.N.; Sivak, K.V.; Chubarev, V.N.; Porozov, Y.B.; Savateeva-Lyubimova, T.N.; Peri, F. TLR4 Signaling Pathway Modulators as Potential Therapeutics in Inflammation and Sepsis. Vaccines (Basel) 2017, 5, 34. [Google Scholar] [CrossRef] [Green Version]
  80. Billod, J.M.; Lacetera, A.; Guzman-Caldentey, J.; Martin-Santamaria, S. Computational Approaches to Toll-Like Receptor 4 Modulation. Molecules 2016, 21, 994. [Google Scholar] [CrossRef] [Green Version]
  81. Caroff, M.; Karibian, D. Structure of bacterial lipopolysaccharides. Carbohydr. Res. 2003, 338, 2431–2447. [Google Scholar] [CrossRef]
  82. Rossignol, D.P.; Lynn, M. Antagonism of in vivo and ex vivo response to endotoxin by E5564, a synthetic lipid A analogue. J. Endotoxin Res. 2002, 8, 483–488. [Google Scholar] [CrossRef] [Green Version]
  83. Bowen, W.S.; Minns, L.A.; Johnson, D.A.; Mitchell, T.C.; Hutton, M.M.; Evans, J.T. Selective TRIF-dependent signaling by a synthetic toll-like receptor 4 agonist. Sci. Signal. 2012, 5, ra13. [Google Scholar] [CrossRef] [Green Version]
  84. Mata-Haro, V.; Cekic, C.; Martin, M.; Chilton, P.M.; Casella, C.R.; Mitchell, T.C. The vaccine adjuvant monophosphoryl lipid A as a TRIF-biased agonist of TLR4. Science 2007, 316, 1628–1632. [Google Scholar] [CrossRef]
  85. Opal, S.M.; Laterre, P.F.; Francois, B.; LaRosa, S.P.; Angus, D.C.; Mira, J.P.; Wittebole, X.; Dugernier, T.; Perrotin, D.; Tidswell, M.; et al. Effect of eritoran, an antagonist of MD2-TLR4, on mortality in patients with severe sepsis: The ACCESS randomized trial. JAMA 2013, 309, 1154–1162. [Google Scholar] [CrossRef] [Green Version]
  86. Gnjatic, S.; Sawhney, N.B.; Bhardwaj, N. Toll-like receptor agonists: Are they good adjuvants? Cancer J. 2010, 16, 382–391. [Google Scholar] [CrossRef] [Green Version]
  87. Yoon, S.I.; Hong, M.; Han, G.W.; Wilson, I.A. Crystal structure of soluble MD-1 and its interaction with lipid IVa. Proc. Natl. Acad. Sci. USA 2010, 107, 10990–10995. [Google Scholar] [CrossRef] [Green Version]
  88. Arias, M.A.; Van Roey, G.A.; Tregoning, J.S.; Moutaftsi, M.; Coler, R.N.; Windish, H.P.; Reed, S.G.; Carter, D.; Shattock, R.J. Glucopyranosyl Lipid Adjuvant (GLA), a Synthetic TLR4 agonist, promotes potent systemic and mucosal responses to intranasal immunization with HIVgp140. PLoS ONE 2012, 7, e41144. [Google Scholar] [CrossRef] [Green Version]
  89. Pantel, A.; Cheong, C.; Dandamudi, D.; Shrestha, E.; Mehandru, S.; Brane, L.; Ruane, D.; Teixeira, A.; Bozzacco, L.; Steinman, R.M.; et al. A new synthetic TLR4 agonist, GLA, allows dendritic cells targeted with antigen to elicit Th1 T-cell immunity in vivo. Eur. J. Immunol 2012, 42, 101–109. [Google Scholar] [CrossRef]
  90. Hill, D.L.; Pierson, W.; Bolland, D.J.; Mkindi, C.; Carr, E.J.; Wang, J.; Houard, S.; Wingett, S.W.; Audran, R.; Wallin, E.F.; et al. The adjuvant GLA-SE promotes human Tfh cell expansion and emergence of public TCRbeta clonotypes. J. Exp. Med. 2019, 216, 1857–1873. [Google Scholar] [CrossRef]
  91. Matsunaga, N.; Tsuchimori, N.; Matsumoto, T.; Ii, M. TAK-242 (resatorvid), a small-molecule inhibitor of Toll-like receptor (TLR) 4 signaling, binds selectively to TLR4 and interferes with interactions between TLR4 and its adaptor molecules. Mol. Pharmacol. 2011, 79, 34–41. [Google Scholar] [CrossRef] [Green Version]
  92. Rice, T.W.; Wheeler, A.P.; Bernard, G.R.; Vincent, J.L.; Angus, D.C.; Aikawa, N.; Demeyer, I.; Sainati, S.; Amlot, N.; Cao, C.; et al. A randomized, double-blind, placebo-controlled trial of TAK-242 for the treatment of severe sepsis. Crit. Care. Med. 2010, 38, 1685–1694. [Google Scholar] [CrossRef]
  93. Goel, A.; Kunnumakkara, A.B.; Aggarwal, B.B. Curcumin as “Curecumin”: From kitchen to clinic. Biochem. Pharmacol. 2008, 75, 787–809. [Google Scholar] [CrossRef] [Green Version]
  94. Youn, H.S.; Saitoh, S.I.; Miyake, K.; Hwang, D.H. Inhibition of homodimerization of Toll-like receptor 4 by curcumin. Biochem. Pharmacol. 2006, 72, 62–69. [Google Scholar] [CrossRef]
  95. Gradisar, H.; Keber, M.M.; Pristovsek, P.; Jerala, R. MD-2 as the target of curcumin in the inhibition of response to LPS. J. Leukoc. Biol. 2007, 82, 968–974. [Google Scholar] [CrossRef]
  96. Zhu, H.T.; Bian, C.; Yuan, J.C.; Chu, W.H.; Xiang, X.; Chen, F.; Wang, C.S.; Feng, H.; Lin, J.K. Curcumin attenuates acute inflammatory injury by inhibiting the TLR4/MyD88/NF-kappaB signaling pathway in experimental traumatic brain injury. J. Neuroinflamm. 2014, 11, 59. [Google Scholar] [CrossRef] [Green Version]
  97. Guney Eskiler, G.; Deveci Ozkan, A.; Kaleli, S.; Bilir, C. Inhibition of TLR4/TRIF/IRF3 Signaling Pathway by Curcumin in Breast Cancer Cells. J. Pharm. Pharm. Sci. 2019, 22, 281–291. [Google Scholar] [CrossRef]
  98. Wang, Y.; Shan, X.; Dai, Y.; Jiang, L.; Chen, G.; Zhang, Y.; Wang, Z.; Dong, L.; Wu, J.; Guo, G.; et al. Curcumin Analog L48H37 Prevents Lipopolysaccharide-Induced TLR4 Signaling Pathway Activation and Sepsis via Targeting MD2. J. Pharmacol. Exp. Ther. 2015, 353, 539–550. [Google Scholar] [CrossRef] [Green Version]
  99. Shehzad, A.; Rehman, G.; Lee, Y.S. Curcumin in inflammatory diseases. Biofactors 2013, 39, 69–77. [Google Scholar] [CrossRef]
  100. Panahi, Y.; Khalili, N.; Hosseini, M.S.; Abbasinazari, M.; Sahebkar, A. Lipid-modifying effects of adjunctive therapy with curcuminoids-piperine combination in patients with metabolic syndrome: Results of a randomized controlled trial. Complement. Ther. Med. 2014, 22, 851–857. [Google Scholar] [CrossRef]
  101. Panahi, Y.; Khalili, N.; Sahebi, E.; Namazi, S.; Simental-Mendia, L.E.; Majeed, M.; Sahebkar, A. Effects of Curcuminoids Plus Piperine on Glycemic, Hepatic and Inflammatory Biomarkers in Patients with Type 2 Diabetes Mellitus: A Randomized Double-Blind Placebo-Controlled Trial. Drug Res. (Stuttg) 2018, 68, 403–409. [Google Scholar] [CrossRef]
  102. Momtazi-Borojeni, A.A.; Haftcheshmeh, S.M.; Esmaeili, S.A.; Johnston, T.P.; Abdollahi, E.; Sahebkar, A. Curcumin: A natural modulator of immune cells in systemic lupus erythematosus. Autoimmun. Rev. 2018, 17, 125–135. [Google Scholar] [CrossRef]
  103. Mirzaei, H.; Naseri, G.; Rezaee, R.; Mohammadi, M.; Banikazemi, Z.; Mirzaei, H.R.; Salehi, H.; Peyvandi, M.; Pawelek, J.M.; Sahebkar, A. Curcumin: A new candidate for melanoma therapy? Int. J. Cancer 2016, 139, 1683–1695. [Google Scholar] [CrossRef]
  104. Lelli, D.; Sahebkar, A.; Johnston, T.P.; Pedone, C. Curcumin use in pulmonary diseases: State of the art and future perspectives. Pharmacol. Res. 2017, 115, 133–148. [Google Scholar] [CrossRef]
  105. Youn, H.S.; Kim, Y.S.; Park, Z.Y.; Kim, S.Y.; Choi, N.Y.; Joung, S.M.; Seo, J.A.; Lim, K.M.; Kwak, M.K.; Hwang, D.H.; et al. Sulforaphane suppresses oligomerization of TLR4 in a thiol-dependent manner. J. Immunol. 2010, 184, 411–419. [Google Scholar] [CrossRef] [Green Version]
  106. Heiss, E.; Herhaus, C.; Klimo, K.; Bartsch, H.; Gerhauser, C. Nuclear factor kappa B is a molecular target for sulforaphane-mediated anti-inflammatory mechanisms. J. Biol. Chem. 2001, 276, 32008–32015. [Google Scholar] [CrossRef] [Green Version]
  107. Killeen, M.E.; Englert, J.A.; Stolz, D.B.; Song, M.; Han, Y.; Delude, R.L.; Kellum, J.A.; Fink, M.P. The phase 2 enzyme inducers ethacrynic acid, DL-sulforaphane, and oltipraz inhibit lipopolysaccharide-induced high-mobility group box 1 secretion by RAW 264.7 cells. J. Pharmacol. Exp. Ther. 2006, 316, 1070–1079. [Google Scholar] [CrossRef]
  108. Koo, J.E.; Park, Z.Y.; Kim, N.D.; Lee, J.Y. Sulforaphane inhibits the engagement of LPS with TLR4/MD2 complex by preferential binding to Cys133 in MD2. Biochem. Biophys. Res. Commun. 2013, 434, 600–605. [Google Scholar] [CrossRef]
  109. Monnet, E.; Lapeyre, G.; Poelgeest, E.V.; Jacqmin, P.; Graaf, K.; Reijers, J.; Moerland, M.; Burggraaf, J.; Min, C. Evidence of NI-0101 pharmacological activity, an anti-TLR4 antibody, in a randomized phase I dose escalation study in healthy volunteers receiving LPS. Clin. Pharmacol. Ther. 2017, 101, 200–208. [Google Scholar] [CrossRef]
Figure 1. The drug discovery and development pipeline. The initial step in drug discovery is target identification; the aim is to find a potential target via different approaches like genomics and proteomics. The next step is target validation, which classifies the molecular target and evaluates whether it is suitable for drug development by various in vivo and in vitro methods. Various techniques such as high-throughput screening (HTS), ligand-based drug design (LBDD), pharmacophore modeling, molecular docking, and structure-based drug design (SBDD) are employed to identify hit compounds. These hit molecules are then evaluated and modified to improve the activity of the lead compound (e.g., by quantitative SAR [QSAR] analysis). Prior to human clinical trials, preclinical studies are carried out that provide detailed knowledge on the toxicity and appropriate dose of drugs. Based on these findings, it is decided whether the drug is ready to be tested in humans. Clinical trials are the final stage of the drug development process. Once clinical trials are over, the successful drug is approved by the FDA and becomes available on the market for clinical use.
Figure 1. The drug discovery and development pipeline. The initial step in drug discovery is target identification; the aim is to find a potential target via different approaches like genomics and proteomics. The next step is target validation, which classifies the molecular target and evaluates whether it is suitable for drug development by various in vivo and in vitro methods. Various techniques such as high-throughput screening (HTS), ligand-based drug design (LBDD), pharmacophore modeling, molecular docking, and structure-based drug design (SBDD) are employed to identify hit compounds. These hit molecules are then evaluated and modified to improve the activity of the lead compound (e.g., by quantitative SAR [QSAR] analysis). Prior to human clinical trials, preclinical studies are carried out that provide detailed knowledge on the toxicity and appropriate dose of drugs. Based on these findings, it is decided whether the drug is ready to be tested in humans. Clinical trials are the final stage of the drug development process. Once clinical trials are over, the successful drug is approved by the FDA and becomes available on the market for clinical use.
Molecules 25 00627 g001
Figure 2. TLR4-MD2-LPS complex (PDB ID: 3FXI). (a) The ECD is divided into three subdomains: N-terminal, central, and C-terminal domains. The evolutionary conserved patches A and B resides in the N-terminal domain and central domain, respectively. (b) Key residues involved in the binding of LPS with TLR4-MD2 complex. TLR4 residue R241 is shown in green color and MD2 residues K58, Y102, S118, S120, K122, and G123 shown in magenta color. LPS is represented by red color.
Figure 2. TLR4-MD2-LPS complex (PDB ID: 3FXI). (a) The ECD is divided into three subdomains: N-terminal, central, and C-terminal domains. The evolutionary conserved patches A and B resides in the N-terminal domain and central domain, respectively. (b) Key residues involved in the binding of LPS with TLR4-MD2 complex. TLR4 residue R241 is shown in green color and MD2 residues K58, Y102, S118, S120, K122, and G123 shown in magenta color. LPS is represented by red color.
Molecules 25 00627 g002
Figure 3. The binding mechanism of LPS for TLR4 activation. A series of events takes place prior to LPS recognition by the TLR4-MD2 complex. (1) LBP extracts LPS from the bacterial membrane and (2) transfers it to membrane anchored CD14, where it binds to the hydrophobic pocket located at the N terminus and forms a monomeric complex. (3) CD14 facilitates LPS transfer to the TLR4-MD2 complex where it initiated the intracellular pathway.
Figure 3. The binding mechanism of LPS for TLR4 activation. A series of events takes place prior to LPS recognition by the TLR4-MD2 complex. (1) LBP extracts LPS from the bacterial membrane and (2) transfers it to membrane anchored CD14, where it binds to the hydrophobic pocket located at the N terminus and forms a monomeric complex. (3) CD14 facilitates LPS transfer to the TLR4-MD2 complex where it initiated the intracellular pathway.
Molecules 25 00627 g003
Figure 4. Illustration of TLR4 signaling pathway. Recognition of LPS activates MyD88-dependent and TRIF-dependent pathways. In MyD88-dependent pathway, recruitment of MyD88 by TIRAP initiates the interaction of IRAKs and TRAF6, that leads to the activation of transcription factors. Endocytosis of TLR4 initiates TRIF-dependent signaling, which involves recruitment of TRIF and TRAM, that leads to subsequent induction of TBK1 and IKKε, as well as activation of transcription factor IRF3. AP1, activated protein 1; ERK, extracellular-regulated kinase, IRAK, interleukin receptor-associated kinase; IKK, inhibitor of ĸ B kinase; IRF3, interferon response factor 3; JNK, c-Jun N-terminal kinase; MAPK, mitogen-activated protein kinase; MyD88, myeloid differentiation primary response protein 88; NF-ĸB, nuclear factor ĸB; p38, protein 38; TAK1, transforming growth factor β-activated kinase 1; TBK1, TANK-binding kinase 1; TIRAM, TRIF-related adaptor molecule; TIRAP, TIR domain-containing adaptor protein; TRAF, tumor necrosis factor receptor-associated factor.
Figure 4. Illustration of TLR4 signaling pathway. Recognition of LPS activates MyD88-dependent and TRIF-dependent pathways. In MyD88-dependent pathway, recruitment of MyD88 by TIRAP initiates the interaction of IRAKs and TRAF6, that leads to the activation of transcription factors. Endocytosis of TLR4 initiates TRIF-dependent signaling, which involves recruitment of TRIF and TRAM, that leads to subsequent induction of TBK1 and IKKε, as well as activation of transcription factor IRF3. AP1, activated protein 1; ERK, extracellular-regulated kinase, IRAK, interleukin receptor-associated kinase; IKK, inhibitor of ĸ B kinase; IRF3, interferon response factor 3; JNK, c-Jun N-terminal kinase; MAPK, mitogen-activated protein kinase; MyD88, myeloid differentiation primary response protein 88; NF-ĸB, nuclear factor ĸB; p38, protein 38; TAK1, transforming growth factor β-activated kinase 1; TBK1, TANK-binding kinase 1; TIRAM, TRIF-related adaptor molecule; TIRAP, TIR domain-containing adaptor protein; TRAF, tumor necrosis factor receptor-associated factor.
Molecules 25 00627 g004
Figure 5. Two-dimensional structures of important TLR4 ligands.
Figure 5. Two-dimensional structures of important TLR4 ligands.
Molecules 25 00627 g005
Table 1. 3D structures of human TLR4 domains available in PDB.
Table 1. 3D structures of human TLR4 domains available in PDB.
PDB IDMethodResolutionLengthDomainLigandMutationsRef.
2Z62X-ray1.70 Å27–228TLR4-ECD/MD2Eritoran-[13]
2Z63X-ray2.00 Å27–527TLR4-ECD/MD2Eritoran-[13]
2Z65X-ray2.70 Å27–228TLR4-ECD/MD2Eritoran-[13]
2Z66X-ray1.90 Å381–627TLR4-ECD/MD2Eritoran-[13]
3FXIX-ray3.10 Å27–631TLR4-ECD/MD2LPS-[12]
3UL7X-ray2.37 Å28–226TLR4-ECD-F63W[50]
3UL8X-ray2.50 Å27–228ECD-V134L[50]
3UL9X-ray2.45 Å28–228ECD-M141E[50]
3ULAX-ray3.60 Å27–228ECD-F63W[50]
4G8AX-ray2.40 Å23–629TLR4-ECD/MD2LPSD299G/T399I[51]
5NAMNMR-623–670TM Domain--[17]
5NAONMR-623–657TM Domain--[17]
Table 2. TLR4-targeting ligands currently in clinical trials.
Table 2. TLR4-targeting ligands currently in clinical trials.
CompoundClassDisease ModelActivityPhaseIdentifier/Status
EritoranGlycolipidInsulin sensitivityAntagonistIINCT02321111/Completed 3
GLA-SE 1GlycolipidStage III & IV soft tissue sarcomaAgonistINCT02180698/Completed
GLA (Mart-1 Antigen)GlycolipidStage IIA, IIB, IIC, IIIA, IIIB, IIIC & IV melanomaAgonistINCT02320305/Active
GLA-SEGlycolipidMetastatic colorectal cancer AgonistINCT03982121/Active
GSK1795091GlycolipidCancer, neoplasmsagonistINCT02798978/Completed
NI-0101AntibodyRheumatoid arthritisantagonistIINCT03241108/Completed
G100GlycolipidFollicular low-grade non-Hodgkin’s lymphomaAgonistINCT02501473/Active
CX-01PolysaccharideMyelodysplastic syndrome, AML 2AntagonistINCT02995655/Completed
GLA-SEGlycolipidMerkel cell carcinomaAgonistINCT02035657/Completed
G100GlycolipidLymphomaAgonistIINCT03742804/Not yet recruiting
GLA-SEGlycolipidSchistosomiasisAgonistIINCT03041766/Completed
GLAGlycolipidHookworm infectionAgonistINCT01717950/Completed
JKB-121Small moleculeNonalcoholic steatohepatitisAntagonistIINCT02442687/Completed
IbudilastSmall moleculeMethamphetamine use disorderAntagonistIINCT01860807/Completed
1 GLA: glucopyranosyl lipid A; 2 AML: acute myeloid leukemia; 3 NCT: national clinical trial.

Share and Cite

MDPI and ACS Style

Ain, Q.u.; Batool, M.; Choi, S. TLR4-Targeting Therapeutics: Structural Basis and Computer-Aided Drug Discovery Approaches. Molecules 2020, 25, 627. https://doi.org/10.3390/molecules25030627

AMA Style

Ain Qu, Batool M, Choi S. TLR4-Targeting Therapeutics: Structural Basis and Computer-Aided Drug Discovery Approaches. Molecules. 2020; 25(3):627. https://doi.org/10.3390/molecules25030627

Chicago/Turabian Style

Ain, Qurat ul, Maria Batool, and Sangdun Choi. 2020. "TLR4-Targeting Therapeutics: Structural Basis and Computer-Aided Drug Discovery Approaches" Molecules 25, no. 3: 627. https://doi.org/10.3390/molecules25030627

Article Metrics

Back to TopTop