Dynamics of GLP-1R peptide agonist engagement are correlated with kinetics of G protein activation

The glucagon-like peptide-1 receptor (GLP-1R) has broad physiological roles and is a validated target for treatment of metabolic disorders. Despite recent advances in GLP-1R structure elucidation, detailed mechanistic understanding of how different peptides generate profound differences in G protein-mediated signalling is still lacking. Here we combine cryo-electron microscopy, molecular dynamics simulations, receptor mutagenesis and pharmacological assays, to interrogate the mechanism and consequences of GLP-1R binding to four peptide agonists; glucagon-like peptide-1, oxyntomodulin, exendin-4 and exendin-P5. These data reveal that distinctions in peptide N-terminal interactions and dynamics with the GLP-1R transmembrane domain are reciprocally associated with differences in the allosteric coupling to G proteins. In particular, transient interactions with residues at the base of the binding cavity correlate with enhanced kinetics for G protein activation, providing a rationale for differences in G protein-mediated signalling efficacy from distinct agonists.

Tweezer-type binding cavity formed by the helical folding of a carbazole–pyridine oligomer

We have synthesised a new aromatic foldamer based on the carbazole–pyridine oligomers that adopt helical conformations via dipole-dipole interactions and π-stacking between two ethynyl bond-linked monomers. This foldamer scaffold has...

Crystal structure of a putative short-chain dehydrogenase/reductase from Paraburkholderia xenovorans

Paraburkholderia xenovorans degrades organic wastes, including polychlorinated biphenyls. The atomic structure of a putative dehydrogenase/reductase (SDR) from P. xenovorans (PxSDR) was determined in space group P21 at a resolution of 1.45 Å. PxSDR shares less than 37% sequence identity with any known structure and assembles as a prototypical SDR tetramer. As expected, there is some conformational flexibility and difference in the substrate-binding cavity, which explains the substrate specificity. Uniquely, the cofactor-binding cavity of PxSDR is not well conserved and differs from those of other SDRs. PxSDR has an additional seven amino acids that form an additional unique loop within the cofactor-binding cavity. Further studies are required to determine how these differences affect the enzymatic functions of the SDR.

Exploring the Retinal Binding Cavity of Archaerhodopsin-3 by Replacing the Retinal Chromophore With a Dimethyl Phenylated Derivative

Rhodopsins act as photoreceptors with their chromophore retinal (vitamin-A aldehyde) and they regulate light-dependent biological functions. Archaerhodopsin-3 (AR3) is an outward proton pump that has been widely utilized as a tool for optogenetics, a method for controlling cellular activity by light. To characterize the retinal binding cavity of AR3, we synthesized a dimethyl phenylated retinal derivative, (2E,4E,6E,8E)-9-(2,6-Dimethylphenyl)-3,7-dimethylnona-2,4,6,8-tetraenal (DMP-retinal). QM/MM calculations suggested that DMP-retinal can be incorporated into the opsin of AR3 (archaeopsin-3, AO3). Thus, we introduced DMP-retinal into AO3 to obtain the non-natural holoprotein (AO3-DMP) and compared some molecular properties with those of AO3 with the natural A1-retinal (AO3-A1) or AR3. Light-induced pH change measurements revealed that AO3-DMP maintained slow outward proton pumping. Noteworthy, AO3-DMP had several significant changes in its molecular properties compared with AO3-A1 as follows; 1) spectroscopic measurements revealed that the absorption maximum was shifted from 556 to 508 nm and QM/MM calculations showed that the blue-shift was due to the significant increase in the HOMO-LUMO energy gap of the chromophore with the contribution of some residues around the chromophore, 2) time-resolved spectroscopic measurements revealed the photocycling rate was significantly decreased, and 3) kinetical spectroscopic measurements revealed the sensitivity of the chromophore binding Schiff base to attack by hydroxylamine was significantly increased. The QM/MM calculations show that a cavity space is present at the aromatic ring moiety in the AO3-DMP structure whereas it is absent at the corresponding β-ionone ring moiety in the AO3-A1 structure. We discuss these alterations of the difference in interaction between the natural A1-retinal and the DMP-retinal with binding cavity residues.

Depth of the Steroid Core Location Determines the Mode of Na,K-ATPase Inhibition by Cardiotonic Steroids

Cardiotonic steroids (CTSs) are specific inhibitors of Na,K-ATPase (NKA). They induce diverse physiological effects and were investigated as potential drugs in heart diseases, hypertension, neuroinflammation, antiviral and cancer therapy. Here, we compared the inhibition mode and binding of CTSs, such as ouabain, digoxin and marinobufagenin to NKA from pig and rat kidneys, containing CTSs-sensitive (α1S) and -resistant (α1R) α1-subunit, respectively. Marinobufagenin in contrast to ouabain and digoxin interacted with α1S-NKA reversibly, and its binding constant was reduced due to the decrease in the deepening in the CTSs-binding site and a lower number of contacts between the site and the inhibitor. The formation of a hydrogen bond between Arg111 and Asp122 in α1R-NKA induced the reduction in CTSs’ steroid core deepening that led to the reversible inhibition of α1R-NKA by ouabain and digoxin and the absence of marinobufagenin’s effect on α1R-NKA activity. Our results elucidate that the difference in signaling, and cytotoxic effects of CTSs may be due to the distinction in the deepening of CTSs into the binding side that, in turn, is a result of a bent-in inhibitor steroid core (marinobufagenin in α1S-NKA) or the change of the width of CTSs-binding cavity (all CTSs in α1R-NKA).

DockStream: a docking wrapper to enhance de novo molecular design

Recently, we have released the de novo design platform REINVENT in version 2.0. This improved and extended iteration supports far more features and scoring function components, which allows bespoke and tailor-made protocols to maximize impact in small molecule drug discovery projects. A major obstacle of generative models is producing active compounds, in which predictive (QSAR) models have been applied to enrich target activity. However, QSAR models are inherently limited by their applicability domains. To overcome these limitations, we introduce a structure-based scoring component for REINVENT. DockStream is a flexible, stand-alone molecular docking wrapper that provides access to a collection of ligand embedders and docking backends. Using the benchmarking and analysis workflow provided in DockStream, execution and subsequent analysis of a variety of docking configurations can be automated. Docking algorithms vary greatly in performance depending on the target and the benchmarking and analysis workflow provides a streamlined solution to identifying productive docking configurations. We show that an informative docking configuration can inform the REINVENT agent to optimize towards improving docking scores using public data. With docking activated, REINVENT is able to retain key interactions in the binding site, discard molecules which do not fit the binding cavity, harness unused (sub-)pockets, and improve overall performance in the scaffold-hopping scenario. The code is freely available at https://github.com/MolecularAI/DockStream.

Molecular Dynamic Simulation Analysis on the Inclusion Complexation of Plumbagin with β-Cyclodextrin Derivatives in Aqueous Solution

Stable encapsulation of medically active compounds can lead to longer storage life and facilitate the slow-release mechanism. In this work, the dynamic and molecular interactions between plumbagin molecule with β-cyclodextrin (BCD) and its two derivatives, which are dimethyl-β-cyclodextrin (MBCD), and 2-O-monohydroxypropyl-β-cyclodextrin (HPBCD) were investigated. Molecular dynamics simulations (MD) with GLYCAM-06 and AMBER force fields were used to simulate the inclusion complex systems under storage temperature (4 °C) in an aqueous solution. The simulation results suggested that HPBCD is the best encapsulation agent to produce stable host–guest binding with plumbagin. Moreover, the observation of the plumbagin dynamic inside the binding cavity revealed that it tends to orient the methyl group toward the wider rim of HPBCD. Therefore, HPBCD is a decent candidate for the preservation of plumbagin with a promising longer storage life and presents the opportunity to facilitate the slow-release mechanism.

Fragment-Based Virtual Screening to Discover Potential Plasmodium PI4KIIIβ Ligands

Plasmodium species that cause malaria, a disease responsible for about half a million deaths per annum despite concerted efforts to combat it. The causative agent depends on type III beta phosphatidylinositol 4-kinase (PPI4K) during the development of merozoite. PPI4K is the only clinically validated Plasmodium kinase so far and its inhibitors are effective both in vitro and in vivo. In this work, a small library of ~22 000 fragments was virtually screened using PPI4K homology model to discover potential ligands of the enzyme. 16 virtual hits were selected based on ≤ -9.0 kcal/mol binding energy cut off and were subjected to similarity and substructure searching after they had passed PAINS screening. The derivatives obtained showed improved binding energies, which ranged from -10.00 to -13.80 kcal/mol. Moreover, the topmost ranking compound 31, with interesting drug-like quality was stable within the protein’s binding cavity during the 10 ns molecular dynamics simulation period. In addition, analysis of its binding pose revealed some unique binding interactions with PPI4K active site residues as the basis for the observed improved binding affinity. Overall, compound 31 appears to be a viable starting point for the development of PPI4K inhibitors with antimalarial activity.

POTENTIAL COUMARIN THIOSEMICARBAZONE HYBRIDS AS BRCA-1 MIMETICS FOR ER POSITIVE BREAST CANCER THERAPY: AN IN-SILICO APPROACH

The goal of this study is to find a novel class of BRCA-1 mimics for Estrogen Receptor α (ERα) breast cancer that works differently from conventional anti-estrogens. Breast Cancer Susceptibility Protein-1 (BRCA-1) is a protein was discovered to bind to ERα and decrease its function by a direct contact between regions inside BRCA1's amino terminus and the carboxy terminus of ERα. A novel class of hybrids with coumate and thiosemicarbazone scaffolds was created with the premise of developing small compounds that imitate the function of BRCA-1 to down regulate the ERα and inhibit the tumor activity of breast cancer cells. Using Schrodinger 2020-2, ADMET and in silico molecular docking tests of the proposed hybrids were performed on the BRCA-1 binding cavity of ERα. TSCO-XIV and TSCO-III are developed hybrids that have high docking scores and good binding interactions with important residues.

Specific Engineered G Protein Coupling to Histamine Receptors Revealed from Cellular Assay Experiments and Accelerated Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are targets of extracellular stimuli and hence occupy a key position in drug discovery. By specific and not yet fully elucidated coupling profiles with α subunits of distinct G protein families, they regulate cellular responses. The histamine H2 and H4 receptors (H2R and H4R) are prominent members of Gs- and Gi-coupled GPCRs. Nevertheless, promiscuous G protein and selective Gi signaling have been reported for the H2R and H4R, respectively, the molecular mechanism of which remained unclear. Using a combination of cellular experimental assays and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigated the coupling profiles of the H2R and H4R to engineered mini-G proteins (mG). We obtained coupling profiles of the mGs, mGsi, or mGsq proteins to the H2R and H4R from the mini-G protein recruitment assays using HEK293T cells. Compared to H2R–mGs expressing cells, histamine responses were weaker (pEC50, Emax) for H2R–mGsi and –mGsq. By contrast, the H4R selectively bound to mGsi. Similarly, in all-atom GaMD simulations, we observed a preferential binding of H2R to mGs and H4R to mGsi revealed by the structural flexibility and free energy landscapes of the complexes. Although the mG α5 helices were consistently located within the HR binding cavity, alternative binding orientations were detected in the complexes. Due to the specific residue interactions, all mG α5 helices of the H2R complexes adopted the Gs-like orientation toward the receptor transmembrane (TM) 6 domain, whereas in H4R complexes, only mGsi was in the Gi-like orientation toward TM2, which was in agreement with Gs- and Gi-coupled GPCRs structures resolved by X-ray/cryo-EM. These cellular and molecular insights support (patho)physiological profiles of the histamine receptors, especially the hitherto little studied H2R function in the brain, as well as of the pharmacological potential of H4R selective drugs.