Molecular Modeling of Chemoreceptor:Ligand Interactions

Unpacking Alternative Features of the Bacterial Chemotaxis System

The bacterial chemotaxis system is one of the best-understood cellular pathways and serves as the model for signal transduction systems. Most chemotaxis research has been conducted with transmembrane chemotaxis systems from Escherichia coli and has established paradigms of the system that were thought to be universal. However, emerging research has revealed that many bacteria possess alternative features of their chemotaxis system, demonstrating that these systems are likely more complex than previously assumed. Here, we compare the canonical chemotaxis system of E. coli with systems that diverge in supramolecular architecture, sensory mechanisms, and protein composition. The alternative features have likely evolved to accommodate chemical specificities of natural niches and cell morphologies. Collectively, these studies demonstrate that bacterial chemotaxis systems are a rapidly expanding field that offers many new opportunities to explore this exceedingly diverse system.

Molecular Mechanism for Attractant Signaling to DHMA by E. coli Tsr

The attractant chemotaxis response of Escherichia coli to norepinephrine requires that it be converted to 3,4-dihydroxymandelic acid (DHMA) by the monoamine oxidase TynA and the aromatic aldehyde dehydrogenase FeaB. DHMA is sensed by the serine chemoreceptor Tsr, and the attractant response requires that at least one subunit of the periplasmic domain of the Tsr homodimer (pTsr) has an intact serine-binding site. DHMA that is generated in vivo by E. coli is expected to be a racemic mixture of the (R) and (S) enantiomers, so it has been unclear whether one or both chiral forms are active. Here, we used a combination of state-of-the-art tools in molecular docking and simulations, including an in-house simulation-based docking protocol, to investigate the binding properties of (R)-DHMA and (S)-DHMA to E. coli pTsr. Our studies computationally predicted that (R)-DHMA should promote a stronger attractant response than (S)-DHMA because of a consistently greater-magnitude piston-like pushdown of the pTsr α-helix 4 toward the membrane upon binding of (R)-DHMA than upon binding of (S)-DHMA. This displacement is caused primarily by interaction of DHMA with Tsr residue Thr156, which has been shown by genetic studies to be critical for the attractant response to L-serine and DHMA. These findings led us to separate the two chiral species and test their effectiveness as chemoattractants. Both the tethered cell and motility migration coefficient assays validated the prediction that (R)-DHMA is a stronger attractant than (S)-DHMA. Our study demonstrates that refined computational docking and simulation studies combined with experiments can be used to investigate situations in which subtle differences between ligands may lead to diverse chemotactic responses.

Isoflavones as Ah Receptor Agonists in Colon-Derived Cell Lines: Structure-Activity Relationships

Many of the protective responses observed for flavonoids in the gastrointestinal track resemble aryl hydrocarbon receptor (AhR)-mediated effects. Therefore, we examined the structure-activity relationships of isoflavones and isomeric flavone and flavanones as AhR ligands on the basis of their induction of CYP1A1, CYP1B1, and UGT1A1 gene expression in colon cancer Caco2 cells and young adult mouse colonocyte (YAMC) cells. Caco2 cells were significantly more Ah-responsive than YAMC cells, and this was due, in part, to flavonoid-induced cytotoxicity in the latter cell lines. The structure-activity relationships for the flavonoids were complex and both response and cell context specific; however, there was significant variability in the AhR activities of the isomeric substituted isoflavones and flavones. For example, 4',5,7-trihydroxyisoflavone (genistein) was AhR-inactive whereas 4',5,7-trihydroxyflavone (apigenin) induced CYP1A1, CYP1B1, and UGT1A1 in Caco2 cells. In contrast, both 5,7-dihydroxy-4-methoxy substituted isoflavone (biochanin A) and flavone (acacetin) induced all three AhR-responsive genes; 4',5,7-trimethoxyisoflavone was a potent AhR agonist, and the isomeric flavone was AhR-inactive. These results coupled with simulation studies modeling flavonoid interaction within the AhR binding pocket demonstrate that the orientation of the substituted phenyl ring at C-2 (flavones) or C-3 (isoflavones) on the common 4-H-chromen-4-one ring strongly influences the activities of isoflavones and flavones as AhR agonists.

Activation of COUP-TFI by a Novel Diindolylmethane Derivative

Chicken ovalbumin upstream promoter-transcription factor I (COUP-TFI) is an orphan receptor and member of the nuclear receptor superfamily. Among a series of methylene substituted diindolylmethanes (C-DIMs) containing substituted phenyl and heteroaromatic groups, we identified 1,1-bis(3'-indolyl)-1-(4-pyridyl)-methane (DIM-C-Pyr-4) as an activator of COUP-TFI. Structure activity studies with structurally diverse heteroaromatic C-DIMs showed that the pyridyl substituted compound was active and the 4-pyridyl substituent was more potent than the 2- or 3-pyridyl analogs in transactivation assays in breast cancer cells. The DIM-C-Pyr-4 activated chimeric GAL4-COUP-TFI constructs containing full length, C- or N-terminal deletions, and transactivation was inhibited by phosphatidylinositol-3-kinase and protein kinase A inhibitors. However, DIM-C-Pyr-4 also induced transactivation and interactions of COUP-TFI and steroid receptor coactivators-1 and -2 in mammalian two-hybrid assays, and ligand-induced interactions of the C-terminal region of COUP-TFI were not affected by kinase inhibitors. We also showed that DIM-C-Pyr-4 activated COUP-TFI-dependent early growth response 1 (Egr-1) expression and this response primarily involved COUP-TFI interactions with Sp3 and to a lesser extent Sp1 bound to the proximal region of the Egr-1 promoter. Modeling studies showed interactions of DIM-C-Pyr-4 within the ligand binding domain of COUP-TFI. This report is the first to identify a COUP-TFI agonist and demonstrate activation of COUP-TFI-dependent Egr-1 expression.

Virtual Screening of Chemical Compounds for Discovery of Complement C3 Ligands

The complement system is our first line of defense against foreign pathogens, but when it is not properly regulated, complement is implicated in the pathology of several autoimmune and inflammatory disorders. Compstatin is a peptidic complement inhibitor that acts by blocking the cleavage of complement protein C3 to the proinflammatory fragment C3a and opsonin fragment C3b. In this study, we aim to identify druglike small-molecule complement inhibitors with physicochemical, geometric, and binding properties similar to those of compstatin. We employed two approaches using various high-throughput virtual screening methods, which incorporate molecular dynamics (MD) simulations, pharmacophore model design, energy calculations, and molecular docking and scoring. We have generated a library of 274 chemical compounds with computationally predicted binding affinities for the compstatin binding site of C3. We have tested subsets of these chemical compounds experimentally for complement inhibitory activity, using hemolytic assays, and for binding affinity, using microscale thermophoresis. As a result, although none of the compounds showed inhibitory activity, compound 29 was identified to exhibit weak competitive binding against a potent compstatin analogue, therefore validating our computational approaches. Additional docking and MD simulation studies suggest that compound 29 interacts with C3 residues, which have been shown to be important in binding of compstatin to the C3c fragment of C3. Compound 29 is amenable to physicochemical optimization to acquire inhibitory properties. Additionally, it is possible that some of the untested compounds will demonstrate binding and inhibition in future experimental studies.