Characterization of the mechanism of bile salt hydrolase substrate specificity by experimental and computational analyses

Bile salt hydrolases (BSHs) are currently being investigated as target enzymes for metabolic regulators in humans and as growth promoters in farm animals. Understanding structural features underlying substrate specificity is necessary for inhibitor design. Here, we used a multidisciplinary workflow including mass spectrometry, mutagenesis, molecular dynamic simulations, machine learning, and crystallography to demonstrate substrate specificity in Lactobacillus salivarius BSH, the most abundant enzyme in human and farm animal intestines. We show the preference of substrates with a taurine head and a dehydroxylated sterol ring for hydrolysis. A regression model that correlates the relative rates of hydrolysis of various substrates in various enzyme mutants with the residue-substrate interaction energies guided the identification of structural determinants of substrate binding and specificity. In addition, we found T208 from another BSH protomer regulating the hydrolysis. The designed workflow can be used for fast and comprehensive characterization of enzymes with a broad range of substrates.

Intermolecular Interactions in G Protein-Coupled Receptor Allosteric Sites at the Membrane Interface from Molecular Dynamics Simulations and Quantum Chemical Calculations

Allosteric modulators are called promising candidates in G protein-coupled receptor (GPCR) drug development by displaying subtype selectivity and more specific receptor modulation. Among the allosteric sites known to date, cavities at the receptor-lipid interface represent an uncharacteristic binding location that raises many questions about the ligand interactions and stability, the binding site structure, and how all of these are affected by lipid molecules. In this work, we analyze interactions in the allosteric sites of the PAR2, C5aR1, and GCGR receptors in three lipid compositions using molecular dynamics simulations. In addition, we performed quantum chemical calculations involving the symmetry-adapted perturbation theory (SAPT) and the natural population analysis to quantify the strength of intermolecular interactions. We show that besides classical hydrogen bonds, weak polar interactions such as O-HC, O-Br, and long-range electrostatics with the backbone amides contribute to the stability of allosteric modulators at the receptor-lipid interface. The allosteric cavities are detectable in various membrane compositions. The availability of polar atoms for interactions in such cavities can be assessed by water molecules from simulations. Although ligand-lipid interactions are weak, lipid tails play a role in ligand binding pose stability and the size of allosteric cavities. We discuss physicochemical aspects of ligand binding at the receptor-lipid interface and suggest a compound library enriched by weak donor groups for ligand search in such sites.

Biphenyl scaffold for the design of NMDA-receptor negative modulators: molecular modeling, synthesis, and biological activity

NMDA (N-methyl-d-aspartate) receptor antagonists are promising tools for the treatment of a wide variety of central nervous system impairments including major depressive disorder. We present here the activity optimization process of a biphenyl-based NMDA negative allosteric modulator (NAM) guided by free energy calculations, which led to a 100 times activity improvement (IC50 = 50 nM) compared to a hit compound identified in virtual screening. Preliminary calculation results suggest a low affinity for the human ether-a-go-go-related gene ion channel (hERG), a high affinity for which was earlier one of the main obstacles for the development of first-generation NMDA-receptor negative allosteric modulators. The docking study and the molecular dynamics calculations suggest a completely different binding mode (ifenprodil-like) compared to another biaryl-based NMDA NAM EVT-101.

Probe Confined Dynamic Mapping for G Protein-Coupled Receptor Allosteric Site Prediction

Targeting G protein-coupled receptors (GPCRs) through allosteric sites offers advantages over orthosteric sites in identifying drugs with increased selectivity and potentially reduced side effects. In this study, we developed a probe confined dynamic mapping protocol that allows the prediction of allosteric sites at both the GPCR extracellular and intracellular sides, as well as at the receptor-lipid interface. The applied harmonic wall potential enhanced sampling of probe molecules in a selected area of a GPCR while preventing membrane distortion in molecular dynamics simulations. The specific probes derived from GPCR allosteric ligand structures performed better in allosteric site mapping compared to commonly used cosolvents. The M₂ muscarinic, β₂ adrenergic, and P₂Y₁ purinergic receptors were selected for the protocol's retrospective validation. The protocol was next validated prospectively to locate the binding site of [5-fluoro-4-(hydroxymethyl)-2-methoxyphenyl]-(4-fluoro-1H-indol-1-yl)methanone at the D₂ dopamine receptor, and subsequent mutagenesis confirmed the prediction. The protocol provides fast and efficient prediction of key amino acid residues surrounding allosteric sites in membrane proteins and facilitates the structure-based design of allosteric modulators.

A Facile Approach to Bis(isoxazoles), Promising Ligands of the AMPA Receptor

A convenient synthetic approach to novel functionalized bis(isoxazoles), the promising bivalent ligands of the AMPA receptor, was elaborated. It was based on the heterocyclization reactions of readily available electrophilic alkenes with the tetranitromethane-triethylamine complex. The structural diversity of the synthesized compounds was demonstrated. In the electrophysiological experiments using the patch clamp technique on Purkinje neurons, the compound 1,4-phenylenedi(methylene)bis(5-aminoisoxazole-3-carboxylate) was shown to be highly potent positive modulator of the AMPA receptor, potentiating kainate-induced currents up to 70% at 10⁻¹¹ M.

[Ligands of the AMPA-subtype glutamate receptors: mechanisms of action and novel chemotypes]

Ionotropic glutamate receptors of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) subtype play a key role in synaptic plasticity representing one of the mechanisms for learning and memory formation. They can also serve as targets for the development of novel classes of pharmaceuticals for the treatment or substantive correction of many serious neurodegenerative and psychoneurological disorders. The search and studies of various types of AMPA receptor ligands attract considerable attention from academic organizations and pharmaceutical companies all over the world. This review mainly focuses on recent advances in this field. The architecture and operational mechanism of the receptor as well as its major binding sites and ligand types are considered. Special attention is paid to the studies of mechanisms of action and novel chemotypes of AMPA receptor agonists and competitive antagonists, positive and negative allosteric modulators, auxiliary protein-dependent allosteric modulators, and ion channel blockers.

Bivalent AMPA receptor positive allosteric modulators of the bis(pyrimidine) series

The first example of a novel class of AMPA receptor positive allosteric modulators of the bis(pyrimidine) series having a hydroquinone linker has been obtained and showed a potency to increase kainate-induced currents at subnanomolar concentrations.

Chemical space exploration guided by deep neural networks

A parametric t-SNE approach based on deep feed-forward neural networks was applied to the chemical space visualization problem. It is able to retain more information than certain dimensionality reduction techniques used for this purpose (principal component analysis (PCA), multidimensional scaling (MDS)). The applicability of this method to some chemical space navigation tasks (activity cliffs and activity landscapes identification) is discussed. We created a simple web tool to illustrate our work (http://space.syntelly.com).

A parametric t-SNE approach based on deep feed-forward neural networks was applied to the chemical space visualization problem.

Cytokinin activity of N6-benzyladenine derivatives assayed by interaction with the receptors in planta, in vitro, and in silico

Biological effects of hormones in both plants and animals are based on high-affinity interaction with cognate receptors resulting in their activation. The signal of cytokinins, classical plant hormones, is perceived in Arabidopsis by three homologous membrane receptors: AHK2, AHK3, and CRE1/AHK4. To study the cytokinin-receptor interaction, we used 25 derivatives of potent cytokinin N⁶-benzyladenine (BA) with substituents in the purine heterocycle and/or in the side chain. The study was focused primarily on individual cytokinin receptors from Arabidopsis. The main in planta assay system was based on Arabidopsis double mutants retaining only one isoform of cytokinin receptors and harboring cytokinin-sensitive reporter gene. Classical cytokinin biotest with Amaranthus seedlings was used as an additional biotest. In parallel, the binding of ligands to individual cytokinin receptors was assessed in the in vitro test system. Quantitative comparison of results of different assays confirmed the partial similarity of ligand-binding properties of receptor isoforms. Substituents at positions 8 and 9 of adenine moiety, elongated linker up to 4 methylene units, and replacement of N⁶ by sulfur or oxygen have resulted in the suppression of cytokinin activity of the derivative toward all receptors. Introduction of a halogen into position 2 of adenine moiety, on the contrary, often increased the ligand activity, especially toward AHK3. Features both common and distinctive of cytokinin receptors in Arabidopsis and Amaranthus were revealed, highlighting species specificity of the cytokinin perception apparatus. Correlations between the extent to which a compound binds to a receptor in vitro and its ability to activate the same receptor in planta were evaluated for each AHK protein. Interaction patterns between individual receptors and ligands were rationalized by structure analysis and molecular docking in sensory modules of AHK receptors. The best correlation between docking scores and specific binding was observed for AHK3. In addition, receptor-specific ligands have been discovered with unique properties to predominantly activate or block distinct cytokinin receptors. These ligands are promising for practical application and as molecular tools in the study of the cytokinin perception by plant cells.

Molecular design of proneurogenic and neuroprotective compounds-allosteric NMDA receptor modulators

N-Methyl-D-aspartic acid (NMDA) receptor is a promising target for treatment of neurodegenerative diseases and other brain disorders as well as for designing proneurogenic compounds able to stimulate neurogenesis in adult brain. We analyzed the structure of the binding site of negative allosteric modulators in the amino-terminal domain of the NMDA receptor and identified possible modes of their binding as well as performed molecular design of new modulators that significantly differ from the known ones in structure and binding mode. In addition, we formed a focused library of chemical compounds with potential neuroprotective and proneurogenic properties, desirable set of pharmacokinetic properties, and low toxicity, which can be the basis for development of new-generation drugs.

On mechanism of allosteric modulation of NMDA receptor via amino-terminal domains

A possible mechanism of action of the allosteric modulators of NMDA (N-methyl-d-aspartate) receptors is proposed that involves the stabilization of the twisted closed-clamshell configuration of the amino-terminal domains of GluN1 and GluN2B subunits by negative modulators while positive modulators stabilize a roughly parallel tight arrangement of these domains. These respective motions may play an important role in the transition between the open-channel and closed-channel states of the receptor. In addition, some features of the negative modulator binding site found by means of the molecular dynamics study and pocket analysis can be used in the rational design of the allosteric NMDA receptor modulators.