Computational Binding Study Hints at Ecdysone 20-Mono-Oxygenase as the Hitherto Unknown Target for Ring C-Seco Limonoid-Type Insecticides.

The insecticidal property of ring C-seco limonoids has been discovered empirically and the target protein identified, but, to date, the molecular mechanism of action has not been described at the atomic scale. We elucidate on computational grounds whether nine C-seco limonoids present sufficiently high affinity to bind specifically with the putative target enzyme of the insects (ecdysone 20-monooxygenase). To this end, 3D models of ligands and the receptor target were generated and their interaction energies estimated by docking simulations. As a proof of concept, the tetrahydro-isoquinolinyl propenamide derivative QHC is the reference ligand bound to aldosterone synthase in the complex with PDB entry 4ZGX. It served as the 3D template for target modeling via homology. QHC was successfully docked back to its crystal pose in a one-digit nanomolar range. The reported experimental binding affinities span over the nanomolar to lower micromolar range. All nine limonoids were found with strong affinities in the range of -9 < ΔG < -13 kcal/mol. The molt hormone ecdysone showed a comparable ΔG energy of -12 kcal/mol, whereas -11 kcal/mol was the back docking result for the liganded crystal 4ZGX. In conclusion, the nine C-seco limonoids were strong binders on theoretical grounds in an activity range between a ten-fold lower to a ten-fold higher concentration level than insecticide ecdysone with its known target receptor. The comparable or even stronger binding hints at ecdysone 20-monooxygenase as their target biomolecule. Our assumption, however, is in need of future experimental confirmation before conclusions with certainty can be drawn about the true molecular mechanism of action for the C-seco limonoids under scrutiny.

Molecular Dynamics and Docking Simulations of Homologous RsmE Methyltransferases Hints at a General Mechanism for Substrate Release upon Uridine Methylation on 16S rRNA.

In this study, molecular dynamics (MD) and docking simulations were carried out on the crystal structure of Neisseria Gonorrhoeae RsmE aiming at free energy of binding estimation (ΔGbinding) of the methyl transfer substrate S-adenosylmethionine (SAM), as well as its homocysteine precursor S-adenosylhomocysteine (SAH). The mechanistic insight gained was generalized in view of existing homology to two other crystal structures of RsmE from Escherichia coli and Aquifex aeolicus. As a proof of concept, the crystal poses of SAM and SAH were reproduced reflecting a more general pattern of molecular interaction for bacterial RsmEs. Our results suggest that a distinct set of conserved residues on loop segments between ß12, α6, and Met169 are interacting with SAM and SAH across these bacterial methyltransferases. Comparing molecular movements over time (MD trajectories) between Neisseria gonorrhoeae RsmE alone or in the presence of SAH revealed a hitherto unknown gatekeeper mechanism by two isoleucine residues, Ile171 and Ile219. The proposed gating allows switching from an open to a closed state, mimicking a double latch lock. Additionally, two key residues, Arg221 and Thr222, were identified to assist the exit mechanism of SAH, which could not be observed in the crystal structures. To the best of our knowledge, this study describes for the first time a general catalytic mechanism of bacterial RsmE on theoretical ground.

Targeting the Human Influenza a Virus: The Methods, Limitations, and Pitfalls of Virtual Screening for Drug-like Candidates Including Scaffold Hopping and Compound Profiling.

In this study, we describe the input data and processing steps to find antiviral lead compounds by a virtual screen. Two-dimensional and three-dimensional filters were designed based on the X-ray crystallographic structures of viral neuraminidase co-crystallized with substrate sialic acid, substrate-like DANA, and four inhibitors (oseltamivir, zanamivir, laninamivir, and peramivir). As a result, ligand-receptor interactions were modeled, and those necessary for binding were utilized as screen filters. Prospective virtual screening (VS) was carried out in a virtual chemical library of over half a million small organic substances. Orderly filtered moieties were investigated based on 2D- and 3D-predicted binding fingerprints disregarding the "rule-of-five" for drug likeness, and followed by docking and ADMET profiling. Two-dimensional and three-dimensional screening were supervised after enriching the dataset with known reference drugs and decoys. All 2D, 3D, and 4D procedures were calibrated before execution, and were then validated. Presently, two top-ranked substances underwent successful patent filing. In addition, the study demonstrates how to work around reported VS pitfalls in detail.

Allosteric Binding of MDMA to the Human Serotonin Transporter (hSERT) via Ensemble Binding Space Analysis with ΔG Calculations, Induced Fit Docking and Monte Carlo Simulations.

Despite the recent promising results of MDMA (3,4-methylenedioxy-methamphetamine) as a psychotherapeutic agent and its history of misuse, little is known about its molecular mode of action. MDMA enhances monoaminergic neurotransmission in the brain and its valuable psychoactive effects are associated to a dual action on the 5-HT transporter (SERT). This drug inhibits the reuptake of 5-HT (serotonin) and reverses its flow, acting as a substrate for the SERT, which possesses a central binding site (S1) for antidepressants as well as an allosteric (S2) one. Previously, we characterized the spatial binding requirements for MDMA at S1. Here, we propose a structure-based mechanistic model of MDMA occupation and translocation across both binding sites, applying ensemble binding space analyses, electrostatic complementarity, and Monte Carlo energy perturbation theory. Computed results were correlated with experimental data (r = 0.93 and 0.86 for S1 and S2, respectively). Simulations on all hSERT available structures with Gibbs free energy estimations (ΔG) revealed a favourable and pervasive dual binding mode for MDMA at S2, i.e., adopting either a 5-HT or an escitalopram-like orientation. Intermediate ligand conformations were identified within the allosteric site and between the two sites, outlining an internalization pathway for MDMA. Among the strongest and more frequent interactions were salt bridges with Glu494 and Asp328, a H-bond with Thr497, a π-π with Phe556, and a cation-π with Arg104. Similitudes and differences with the allosteric binding of 5-HT and antidepressants suggest that MDMA may have a distinctive chemotype. Thus, our models may provide a framework for future virtual screening studies and pharmaceutical design and to develop hSERT allosteric compounds with a unique psychoactive MDMA-like profile.

Computational Molecular Characterization of the Interaction of Acetylcholine and the NMDA Receptor to Explain the Direct Glycine-Competitive Potentiation of NMDA-Mediated Neuronal Currents.

The activation of N-methyl-d-aspartate receptor (NMDAR) is triggered by the closure of bilobed (D1 and D2) clamshell-like clefts upon binding glycine (Gly) and glutamate. There is evidence that cholinergic compounds modulate NMDAR-mediated currents via direct receptor-ligand interactions; however, molecular bases are unknown. Here, we first propose a mechanistic structure-based explanation for the observed ACh-induced submaximal potentiation of NMDA-elicited currents in striatal neurons by predicting competitive inhibition with Gly. Then, the model was validated, in principle, by confirming that the coapplication of Gly and ACh significantly reduces these neuronal currents. Finally, we delineate the interplay of ACh with the NMDAR by a combination of computational strategies. Crystallographic ACh-bound complexes were studied, revealing a similar ACh binding environment on the GluN1 subunit of the NMDAR. We illustrate how ACh can occupy X-ray monomeric open, dimeric "semiopen" cleft conformations obtained by molecular dynamics and a full-active cryo-EM NMDAR structure, explaining the suboptimal NMDAR electrophysiological activity under the "Venus Flytrap model". At an evolutionary biology level, the binding mode of ACh coincides with that of the homologous ornithine-bound periplasmic LAO binding protein complex. Our computed results indicate an analogous mechanism of action, inasmuch as ACh may stabilize the GluN1 subunit "semiclosed" conformations by inducing direct and indirect D1-to-D2 interdomain bonds. Additionally, an alternative binding site was detected, shared by the known NMDAR allosteric modulators. Experimental and computed results strongly suggest that ACh acts as a Gly-competitive, submaximal potentiating agent of the NMDAR, possibly constituting a novel chemotype for multitarget-directed drug development, e.g., to treat Alzheimer's, and it may lead to a new understanding of glutamatergic neurotransmission.

Induced fit, ensemble binding space docking and Monte Carlo simulations of MDMA 'ecstasy' and 3D pharmacophore design of MDMA derivatives on the human serotonin transporter (hSERT).

The popular recreational drug MDMA (3,4-methylenedioxy-methamphetamine) has a documented potential as a psychopharmacological clinical and research tool. This is due to its unique ability to promote reprocessing of traumatic memories, empathetic and pro-social states. Although it is established that MDMA exerts its behavioural effects via the serotonin transporter (SERT), the ligand-protein molecular interplay remains elusive. In order to shed light on the binding of MDMA and its primary congeneric entactogens (MDA, MBDB and MDAI), we first combined induced fit with Monte Carlo simulations. The computed interaction energies of the models correlated well with experimental activities (adjR2 = 0.78). Then we carried out 'ensemble binding space docking' on trajectories generated by interpolation of experimentally derived structures of the hSERT from the outward-open, and the occluded, to the inward-open states. This approach revealed low-energy alternative binding modes, suggesting high occupancy of the central site, yet considerable MDMA mobility within it, favouring the paroxetine-like orientation. Finally, we designed a pharmacophore that may be used to recognise hSERT-mediated serotonin releasers and uptake inhibitors of diverse chemical structure, identifying their active conformations and interacting residues. We conclude that the conserved amine-Asp98 ionic and edge-to-face π-π interactions are crucial to the mode of action of MDMA on the hSERT, underscoring the contributions of Tyr95 and gating residues Phe341, Tyr176 and Phe335. Amenable to experimental testing, our modelling may aid the rational design of novel entactogenic compounds and contribute to the understanding of an action mechanism, common and typical of psychotropic agents.

Are vanadium complexes druggable against the main protease Mpro of SARS-CoV-2? - A computational approach.

In silico techniques helped explore the binding capacities of the SARS-CoV-2 main protease (Mpro) for a series of metalloorganic compounds. Along with small size vanadium complexes a vanadium-containing derivative of the peptide-like inhibitor N3 (N-[(5-methylisoxazol-3-yl)carbonyl]alanyl-l-valyl-N1-((1R,2Z)-4-(benzyloxy)-4-oxo-1-{[(3R)-2-oxopyrrolidin-3-yl] methyl }but-2-enyl)-l-leucinamide) was designed from the crystal structure with PDB entry code 6LU7. On theoretical grounds our consensus docking studies evaluated the binding affinities at the hitherto known binding site of Chymotrypsin-like protease (3CLpro) of SARS-CoV-2 for existing and designed vanadium complexes. This main virus protease (Mpro) has a Cys-His dyad at the catalytic site that is characteristic of metal-dependent or metal-inhibited hydrolases. Mpro was compared to the human protein-tyrosine phosphatase 1B (hPTP1B) with a comparable catalytic dyad. HPTP1B is a key regulator at an early stage in the signalling cascade of the insulin hormone for glucose uptake into cells. The vanadium-ligand binding site of hPTP1B is located in a larger groove on the surface of Mpro. Vanadium constitutes a well-known phosphate analogue. Hence, its study offers possibilities to design promising vanadium-containing binders to SARS-CoV-2. Given the favourable physicochemical properties of vanadium nuclei, such organic vanadium complexes could become drugs not only for pharmacotherapy but also diagnostic tools for early infection detection in patients. This work presents the in silico design of a potential lead vanadium compound. It was tested along with 20 other vanadium-containing complexes from the literature in a virtual screening test by docking to inhibit Mpro of SARS-CoV-2.

Five Novel Non-Sialic Acid-Like Scaffolds Inhibit In Vitro H1N1 and H5N2 Neuraminidase Activity of Influenza a Virus.

Neuraminidase (NA) of influenza viruses enables the virus to access the cell membrane. It degrades the sialic acid contained in extracellular mucin. Later, it is responsible for releasing newly formed virions from the membrane of infected cells. Both processes become key functions within the viral cycle. Therefore, it is a therapeutic target for research of the new antiviral agents. Structure-activity relationships studies have revealed which are the important functional groups for the receptor-ligand interaction. Influenza virus type A NA activity was inhibited by five scaffolds without structural resemblance to sialic acid. Intending small organic compound repositioning along with drug repurposing, this study combined in silico simulations of ligand docking into the known binding site of NA, along with in vitro bioassays. The five proposed scaffolds are N-acetylphenylalanylmethionine, propanoic 3-[(2,5-dimethylphenyl) carbamoyl]-2-(piperazin-1-yl) acid, 3-(propylaminosulfonyl)-4-chlorobenzoic acid, ascorbic acid (vitamin C), and 4-(dipropylsulfamoyl) benzoic acid (probenecid). Their half maximal inhibitory concentration (IC50) was determined through fluorometry. An acidic reagent 2'-O-(4-methylumbelliferyl)-α-dN-acetylneuraminic acid (MUNANA) was used as substrate for viruses of human influenza H1N1 or avian influenza H5N2. Inhibition was observed in millimolar ranges in a concentration-dependent manner. The IC50 values of the five proposed scaffolds ranged from 6.4 to 73 mM. The values reflect a significant affinity difference with respect to the reference drug zanamivir (p < 0.001). Two compounds (N-acetyl dipeptide and 4-substituted benzoic acid) clearly showed competitive mechanisms, whereas ascorbic acid reflected non-competitive kinetics. The five small organic molecules constitute five different scaffolds with moderate NA affinities. They are proposed as lead compounds for developing new NA inhibitors which are not analogous to sialic acid.

Chemometric Models of Differential Amino Acids at the Navα and Navß Interface of Mammalian Sodium Channel Isoforms.

(1) Background: voltage-gated sodium channels (Navs) are integral membrane proteins that allow the sodium ion flux into the excitable cells and initiate the action potential. They comprise an α (Navα) subunit that forms the channel pore and are coupled to one or more auxiliary ß (Navß) subunits that modulate the gating to a variable extent. (2) Methods: after performing homology in silico modeling for all nine isoforms (Nav1.1α to Nav1.9α), the Navα and Navß protein-protein interaction (PPI) was analyzed chemometrically based on the primary and secondary structures as well as topological or spatial mapping. (3) Results: our findings reveal a unique isoform-specific correspondence between certain segments of the extracellular loops of the Navα subunits. Precisely, loop S5 in domain I forms part of the PPI and assists Navß1 or Navß3 on all nine mammalian isoforms. The implied molecular movements resemble macroscopic springs, all of which explains published voltage sensor effects on sodium channel fast inactivation in gating. (4) Conclusions: currently, the specific functions exerted by the Navß1 or Navß3 subunits on the modulation of Navα gating remain unknown. Our work determined functional interaction in the extracellular domains on theoretical grounds and we propose a schematic model of the gating mechanism of fast channel sodium current inactivation by educated guessing.

Dapsone is not a Pharmacodynamic Lead Compound for its Aryl Derivatives.

BACKGROUND: The relatedness between the linear equations of thermodynamics and QSAR was studied thanks to the recently elucidated crystal structure complexes between sulfonamide pterin conjugates and dihydropteroate synthase (DHPS) together with a published set of thirty- six synthetic dapsone derivatives with their reported entropy-driven activity data. Only a few congeners were slightly better than dapsone. OBJECTIVE: Our study aimed at demonstrating the applicability of thermodynamic QSAR and to shed light on the mechanistic aspects of sulfone binding to DHPS. METHODS: To this end ligand docking to DHPS, quantum mechanical properties, 2D- and 3D-QSAR as well as Principle Component Analysis (PCA) were carried out. RESULTS: The short aryl substituents of the docked pterin-sulfa conjugates were outward oriented into the solvent space without interacting with target residues which explains why binding enthalpy (ΔH) did not correlate with potency. PCA revealed how chemically informative descriptors are evenly loaded on the first three PCs (interpreted as ΔG, ΔH and ΔS), while chemically cryptic ones reflected higher dimensional (complex) loadings. CONCLUSION: It is safe to utter that synthesis efforts to introduce short side chains for aryl derivatization of the dapsone scaffold have failed in the past. On theoretical grounds we provide computed evidence why dapsone is not a pharmacodynamic lead for drug profiling because enthalpic terms do not change significantly at the moment of ligand binding to target.

Exploring electrostatic patterns of human, murine, equine and canine TLR4/MD-2 receptors.

Electrostatic interactions between phosphate anions and Toll-like receptor 4 / Myeloid differentiation factor-2 (TLR4/MD-2) protein complexes of human, murine, equine and canine species were computed. Such knowledge can provide mechanistic information about recognising LPS-like ligands, since anionic phosphate groups belong to the structural features of LPS with their diphosphorylated diglucosamine backbone. Sequence composition analyses, electrostatic interaction potentials and docked energies as well as molecular dynamics studies evaluated the phosphate interactions within the triangular LPS binding site (wedge). According to electrostatic analyses, human, horse and dog wedges possess phosphate-binding sites with indistinct positive and negative charge distributions, but the murine wedge shows a unique strong negative net charge at the site where antagonists bind in other species (Pan). Docking of a phosphate mono-anion (probe) confirmed its repulsion at this Pan site, but the Pag site of the murine wedge attracted the probe. It is occupied by phosphate groups of agonists in other species (Pag). Molecular dynamics trajectories show a variable degree of random walk across the wedges, that is, not following electrostatic preferences (neither Pag nor Pan). In summary, two opposing electrostatic patterns exist -murine versus human, equine and canine species - all of which reflect the potential dual activity mode of under-acylated ligands such as lipid IVA.

Induced fit for cytochrome P450 3A4 based on molecular dynamics.

The present study aims at numerically describing to what extent substrate - enzyme complexes in solution may change over time as a natural process of conformational changes for a liganded enzyme in comparison to those movements which occur independently from substrate interaction, i.e. without a ligand. To this end, we selected structurally known pairs of liganded / unliganded CYP450 3A4 enzymes with different geometries hinting at induced fit events. We carried out molecular dynamics simulations (MD) comparing the trajectories in a "cross-over" protocol: (i) we added the ligand to the unliganded crystal form which should adopt geometries similar to the known geometry of the liganded crystal structure during MD, and - conversely - (ii) we removed the bound ligand form the known liganded complex to test if a geometry similar to the known unliganded (apo-) form can be adopted during MD. To compare continues changes we measured root means square deviations and frequencies. Results for case (i) hint at larger conformational changes required for accepting the substrate during its approach to final position - in contrast to case (ii) when mobility is fairly reduced by ligand binding (strain energy). In conclusion, a larger conformational sampling prior to ligand binding and the freezing-in (rigidity) of conformations for bound ligands can be interpreted as two conditions linked to induced-fit.

Synthesis, in vitro and in vivo giardicidal activity of nitrothiazole-NSAID chimeras displaying broad antiprotozoal spectrum.

We designed and synthesized five new 5-nitrothiazole-NSAID chimeras as analogues of nitazoxanide, using a DCC-activated amidation. Compounds 1-5 were tested in vitro against a panel of five protozoa: 2 amitochondriates (Giardia intestinalis, Trichomonas vaginalis) and 3 kinetoplastids (Leishmania mexicana, Leishmania amazonensis and Trypanosoma cruzi). All chimeras showed broad spectrum and potent antiprotozoal activities, with IC50 values ranging from the low micromolar to nanomolar order. Compounds 1-5 were even more active than metronidazole and nitazoxanide, two marketed first-line drugs against giardiasis. In particular, compound 4 (an indomethacin hybrid) was one of the most potent of the series, inhibiting G. intestinalis growth in vitro with an IC50 of 0.145µM. Compound 4 was 38-times more potent than metronidazole and 8-times more active than nitazoxanide. The in vivo giardicidal effect of 4 was evaluated in a CD-1 mouse model obtaining a median effective dose of 1.709µg/kg (3.53nmol/kg), a 321-fold and 1015-fold increase in effectiveness after intragastric administration over metronidazole and nitazoxanide, respectively. Compounds 1 and 3 (hybrids of ibuprofen and clofibric acid), showed potent giardicidal activities in the in vitro as well as in the in vivo assays after oral administration. Therefore, compounds 1-5 constitute promising drug candidates for further testing in experimental chemotherapy against giardiasis, trichomoniasis, leishmaniasis and even trypanosomiasis infections.

Characterization of specific allosteric effects of the Na+ channel ß1 subunit on the Nav1.4 isoform.

The mechanism of inactivation of mammalian voltage-gated Na+ channels involves transient interactions between intracellular domains resulting in direct pore occlusion by the IFM motif and concomitant extracellular interactions with the ß1 subunit. Navß1 subunits constitute single-pass transmembrane proteins that form protein-protein associations with pore-forming α subunits to allosterically modulate the Na+ influx into the cell during the action potential of every excitable cell in vertebrates. Here, we explored the role of the intracellular IFM motif of rNav1.4 (skeletal muscle isoform of the rat Na+ channel) on the α-ß1 functional interaction and showed for the first time that the modulation of ß1 is independent of the IFM motif. We found that: (1) Nav1.4 channels that lack the IFM inactivation particle can undergo a "C-type-like inactivation" albeit in an ultraslow gating mode; (2) ß1 can significantly accelerate the inactivation of Nav1.4 channels in the absence of the IFM motif. Previously, we identified two residues (T109 and N110) on the ß1 subunit that disrupt the α-ß1 allosteric modulation. We further characterized the electrophysiological effects of the double alanine substitution of these residues demonstrating that it decelerates inactivation and recovery from inactivation, abolishes the modulation of steady-state inactivation and induces a current rundown upon repetitive stimulation, thus causing a general loss of function. Our results contribute to delineating the process of the mammalian Na+ channel inactivation. These findings may be relevant to the design of pharmacological strategies, targeting ß subunits to treat pathologies associated to Na+ current dysfunction.

Antiprotozoal Nitazoxanide Derivatives: Synthesis, Bioassays and QSAR Study Combined with Docking for Mechanistic Insight.

In view of the serious health problems concerning infectious diseases in heavily populated areas, we followed the strategy of lead compound diversification to evaluate the near-by chemical space for new organic compounds. To this end, twenty derivatives of nitazoxanide (NTZ) were synthesized and tested for activity against Entamoeba histolytica parasites. To ensure drug-likeliness and activity relatedness of the new compounds, the synthetic work was assisted by a quantitative structure-activity relationships study (QSAR). Many of the inherent downsides - well-known to QSAR practitioners - we circumvented thanks to workarounds which we proposed in prior QSAR publication. To gain further mechanistic insight on a molecular level, ligand-enzyme docking simulations were carried out since NTZ is known to inhibit the protozoal pyruvate ferredoxin oxidoreductase (PFOR) enzyme as its biomolecular target.

Predicting a double mutant in the twilight zone of low homology modeling for the skeletal muscle voltage-gated sodium channel subunit beta-1 (Nav1.4 ß1).

The molecular structure modeling of the ß1 subunit of the skeletal muscle voltage-gated sodium channel (Nav1.4) was carried out in the twilight zone of very low homology. Structural significance can per se be confounded with random sequence similarities. Hence, we combined (i) not automated computational modeling of weakly homologous 3D templates, some with interfaces to analogous structures to the pore-bearing Nav1.4 α subunit with (ii) site-directed mutagenesis (SDM), as well as (iii) electrophysiological experiments to study the structure and function of the ß1 subunit. Despite the distant phylogenic relationships, we found a 3D-template to identify two adjacent amino acids leading to the long-awaited loss of function (inactivation) of Nav1.4 channels. This mutant type (T109A, N110A, herein called TANA) was expressed and tested on cells of hamster ovary (CHO). The present electrophysiological results showed that the double alanine substitution TANA disrupted channel inactivation as if the ß1 subunit would not be in complex with the α subunit. Exhaustive and unbiased sampling of "all ß proteins" (Ig-like, Ig) resulted in a plethora of 3D templates which were compared to the target secondary structure prediction. The location of TANA was made possible thanks to another "all ß protein" structure in complex with an irreversible bound protein as well as a reversible protein-protein interface (our "Rosetta Stone" effect). This finding coincides with our electrophysiological data (disrupted ß1-like voltage dependence) and it is safe to utter that the Nav1.4 α/ß1 interface is likely to be of reversible nature.

Identification of Navß1 residues involved in the modulation of the sodium channel Nav1.4.

Voltage-gated sodium channels (VGSCs) are heteromeric protein complexes that initiate action potentials in excitable cells. The voltage-gated sodium channel accessory subunit, Navß1, allosterically modulates the α subunit pore structure upon binding. To date, the molecular determinants of the interface remain unknown. We made use of sequence, knowledge and structure-based methods to identify residues critical to the association of the α and ß1 Nav1.4 subunits. The Navß1 point mutant C43A disrupted the modulation of voltage dependence of activation and inactivation and delayed the peak current decay, the recovery from inactivation, and induced a use-dependent decay upon depolarisation at 1 Hz. The Navß1 mutant R89A selectively delayed channel inactivation and recovery from inactivation and had no effect on voltage dependence or repetitive depolarisations. Navß1 mutants Y32A and G33M selectively modified the half voltage of inactivation without altering the kinetics. Despite low sequence identity, highly conserved structural elements were identified. Our models were consistent with published data and may help relate pathologies associated with VGSCs to the Navß1 subunit.

Reviewing and identifying amino acids of human, murine, canine and equine TLR4 / MD-2 receptor complexes conferring endotoxic innate immunity activation by LPS/lipid A, or antagonistic effects by Eritoran, in contrast to species-dependent modulation by lipid IVa.

There is literature evidence gathered throughout the last two decades reflecting unexpected species differences concerning the immune response to lipid IVa which provides the opportunity to gain more detailed insight by the molecular modeling approach described in this study. Lipid IVa is a tetra-acylated precursor of lipid A in the biosynthesis of lipopolysaccharide (LPS) in Gram-negative bacteria. Lipid A of the prototypic E. coli-type is a hexa-acylated structure that acts as an agonist in all tested mammalian species by innate immunorecognition via the Toll-like receptor 4 (TLR4)/myeloid differentiation factor 2 (MD-2) receptor complex. In contrast, lipid IVa is proinflammatory in mouse cells (agonism) but it remains inactive to human macrophages and even antagonizes the action of potent agonists like E. coli-type lipid A. This particular ambivalent activity profile of lipid IVa has been confirmed in other mammalian species: in equine cells Lipid IVa also acts in a weak agonistic manner, whereas being inactive and antagonizing the lipid A-induced activation of canine TLR4/MD-2. Intriguingly, the respective TLR4 amino acid sequences of the latter species are more identical to the human (67%, 68%) than to the murine (62%, 58%) ortholog. In order to address the unpaired activity-sequence dualism for human, murine, canine and equine species regarding the activity of lipid IVa as compared to LPS and lipid A and, we review the literature and computationally pinpoint the differential biological effects of lipid IVa versus LPS and lipid A to specific amino acid residues. In contrast to lipid IVa the structurally related synthetic compound Eritoran (E5564) acts consistently in an antagonistic manner in these mammalian species and serves as a reference ligand for molecular modeling in this study. The combined evaluation of data sets provided by prior studies and in silico homology mapping of differential residues of TLR4/MD-2 complexes lends detailed insight into the driving forces of the characteristic binding modes of the lipid A domain in LPS and the precursor structure lipid IVa to the receptor complex in individual mammalian species.

Three-dimensional mapping of differential amino acids of human, murine, canine and equine TLR4/MD-2 receptor complexes conferring endotoxic activation by lipid A, antagonism by Eritoran and species-dependent activities of Lipid IVA in the mammalian LPS sensor system.

A literature review concerning the unexpected species differences of the vertebrate innate immune response to lipid IVA was published in CSBJ prior to the present computational study to address the unpaired activity-sequence correlation of prototypic E. coli -type lipid A and its precursor lipid IVA regarding human, murine, equine and canine species. To this end, their sequences and structures of hitherto known Toll-like receptor 4 (TLR4) and myeloid differentiation factor 2 (MD-2) complexes were aligned and their differential side chain patterns studied. If required due to the lack of the corresponding X-ray crystallographic data, three-dimensional models of TLR4/MD-2/ligand complexes were generated using mono and dimeric crystal structures as templates and in silico docking of the prototypic ligands lipid A, lipid IVA and Eritoran. All differential amino acids were mapped to pinpoint species dependency on an atomic scale, i.e. the possible concert of mechanistically relevant side chains. In its most abstract and general form the three-dimensional (3D-) models devise a triangular interface or "wedge" where molecular interactions between TLR4, MD-2 and ligand itself take place. This study identifies two areas in the wedge related to either agonism or antagonism reflecting why ligands like lipid IVA can possess a species dependent dual activity. Lipid IVA represents an imperfect (underacylated and backbone-flipped), low affinity ligand of mammalian TLR4/MD-2 complexes. Its specific but weak antagonistic activity in the human system is in particular due to the loss of phosphate attraction in the wedge-shaped region conferred by nonhomologous residue changes when compared to crystal and modeled structures of the corresponding murine and equine TLR4/MD-2 complexes. The counter-TLR4/MD-2 unit was also taken into account since agonist-mediated dimerization in a defined m-shaped complex composed of two TLR4/MD-2/agonist subunits triggers intracellular signaling during the innate immune response to bacterial endotoxin exposure.

Analgesic effect of leaf extract from Ageratina glabrata in the hot plate test

Ageratina glabrata (Kunth) R.M. King & H. Rob., Asteraceae (syn. Eupatorium glabratum Kunth) is widely distributed throughout Mexico and popularly known as "chamizo blanco" and "hierba del golpe" for its traditional use as external analgesic remedy. Though glabrata species has been chemically studied, there are no experimentally asserted reports about possible analgesic effects which can be inferred from its genus Ageratina. To fill the gap, we evaluated A. glabrata extracts in an animal model of nociception exploiting thermal stimuli. NMR and mass analyses identified a new thymol derivative, 10-benzoiloxy-6,8,9-trihydroxy-thymol isobutyrate (1), which was computationally converted into a ring-closed structure to explain interaction with the COX-2 enzyme in a ligand-receptor docking study. The resulting docked pose is in line with reported crystal complexes of COX-2 with chromene ligands. Based on the present results of dichloromethane extracts from its dried leaves, it is safe to utter that the plant possesses analgesic effects in animal tests which are mediated through inhibition of COX-2 enzyme.