Keras/TensorFlow in Drug Design for Immunity Disorders.

Homeostasis of the host immune system is regulated by white blood cells with a variety of cell surface receptors for cytokines. Chemotactic cytokines (chemokines) activate their receptors to evoke the chemotaxis of immune cells in homeostatic migrations or inflammatory conditions towards inflamed tissue or pathogens. Dysregulation of the immune system leading to disorders such as allergies, autoimmune diseases, or cancer requires efficient, fast-acting drugs to minimize the long-term effects of chronic inflammation. Here, we performed structure-based virtual screening (SBVS) assisted by the Keras/TensorFlow neural network (NN) to find novel compound scaffolds acting on three chemokine receptors: CCR2, CCR3, and one CXC receptor, CXCR3. Keras/TensorFlow NN was used here not as a typically used binary classifier but as an efficient multi-class classifier that can discard not only inactive compounds but also low- or medium-activity compounds. Several compounds proposed by SBVS and NN were tested in 100 ns all-atom molecular dynamics simulations to confirm their binding affinity. To improve the basic binding affinity of the compounds, new chemical modifications were proposed. The modified compounds were compared with known antagonists of these three chemokine receptors. Known CXCR3 compounds were among the top predicted compounds; thus, the benefits of using Keras/TensorFlow in drug discovery have been shown in addition to structure-based approaches. Furthermore, we showed that Keras/TensorFlow NN can accurately predict the receptor subtype selectivity of compounds, for which SBVS often fails. We cross-tested chemokine receptor datasets retrieved from ChEMBL and curated datasets for cannabinoid receptors. The NN model trained on the cannabinoid receptor datasets retrieved from ChEMBL was the most accurate in the receptor subtype selectivity prediction. Among NN models trained on the chemokine receptor datasets, the CXCR3 model showed the highest accuracy in differentiating the receptor subtype for a given compound dataset.

Chemokine Receptors-Structure-Based Virtual Screening Assisted by Machine Learning.

Chemokines modulate the immune response by regulating the migration of immune cells. They are also known to participate in such processes as cell-cell adhesion, allograft rejection, and angiogenesis. Chemokines interact with two different subfamilies of G protein-coupled receptors: conventional chemokine receptors and atypical chemokine receptors. Here, we focused on the former one which has been linked to many inflammatory diseases, including: multiple sclerosis, asthma, nephritis, and rheumatoid arthritis. Available crystal and cryo-EM structures and homology models of six chemokine receptors (CCR1 to CCR6) were described and tested in terms of their usefulness in structure-based drug design. As a result of structure-based virtual screening for CCR2 and CCR3, several new active compounds were proposed. Known inhibitors of CCR1 to CCR6, acquired from ChEMBL, were used as training sets for two machine learning algorithms in ligand-based drug design. Performance of LightGBM was compared with a sequential Keras/TensorFlow model of neural network for these diverse datasets. A combination of structure-based virtual screening with machine learning allowed to propose several active ligands for CCR2 and CCR3 with two distinct compounds predicted as CCR3 actives by all three tested methods: Glide, Keras/TensorFlow NN, and LightGBM. In addition, the performance of these three methods in the prediction of the CCR2/CCR3 receptor subtype selectivity was assessed.

Computational and experimental approaches to probe GPCR activation and signaling.

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.

Helix 8 in chemotactic receptors of the complement system.

Host response to infection involves the activation of the complement system leading to the production of anaphylatoxins C3a and C5a. Complement factor C5a exerts its effect through the activation of C5aR1, chemotactic receptor 1, and triggers the G protein-coupled signaling cascade. Orthosteric and allosteric antagonists of C5aR1 are a novel strategy for anti-inflammatory therapies. Here, we discuss recent crystal structures of inactive C5aR1 in terms of an inverted orientation of helix H8, unobserved in other GPCR structures. An analysis of mutual interactions of subunits in the C5aR1-G protein complex has provided new insights into the activation mechanism of this distinct receptor. By comparing two C5aR receptors C5aR1 and C5aR2 we explained differences between their signaling pathways on the molecular level. By means of molecular dynamics we explained why C5aR2 cannot transduce signal through the G protein pathway but instead recruits beta-arrestin. A comparison of microsecond MD trajectories started from active and inactive C5aR1 receptor conformations has provided insights into details of local and global changes in the transmembrane domain induced by interactions with the Gα subunit and explained the impact of inverted H8 on the C5aR1 activation.