Probing the Druggablility on the Interface of the Protein-Protein Interaction and Its Allosteric Regulation Mechanism on the Drug Screening for the CXCR4 Homodimer

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

Antibacterial activities of plant-derived xanthones

The increasing threat to global health posed by antibiotic resistance remains a serious concern. This troublesome scenario has steered a need for the discovery and evaluation of novel antibacterial agents. Natural products are the main sources of antimicrobials used in clinical practice, serving as a rich reservoir for the discovery of new antibiotics. Pharmaceutical phenolics especially xanthones widely exist in the plant kingdom, and are important plant metabolites. They possess versatile biological activities, including antiviral, antibacterial, neurotrophic, and anticancer. In the present study, we focus on the antibacterial activities of phytoxanthones and summarize their structures and sources, categories and drug-likeness evaluations, and antibacterial activities. A total of 226 different plant xanthones are identified through the NETs screening, and most of them are distributed in Clusiaceae family. These phytoxanthones are divided into four groups according to the intrinsic structural properties, including the most common simple xanthones and the majority of biprenylated ones. Moreover, their physicochemical parameters are calculated and the structure-activity relationships are discussed as well. These results indicate that the biprenylated xanthone derivatives may be promising antibacterial candidates and that the natural products of plants may be a poorly understood repository for the discovery of novel antibacterial agents.

Molecular Mechanisms of Diverse Activation Stimulated by Different Biased Agonists for the β2-Adrenergic Receptor

β2AR is an important drug target protein involving many diseases. Biased drugs induce specific signaling and provide additional clinical utility to optimize β2AR-based therapies. However, the biased signaling mechanism has not been elucidated. Motivated by the issue, we chose four agonists with divergent bias (balanced agonist, G-protein-biased agonist, and β-arrestin-biased agonists) and utilized Gaussian accelerated molecular dynamics simulation coupled with a dynamic network to probe the molecular mechanisms of distinct biased activation induced by the structural differences between the four agonists. Our simulations reveal that the G-protein-biased agonist induces an open conformation with the outward shifts of TM6 and TM7 for the intracellular domain, which will be beneficial to couple G protein. In contrast, the β-arrestin-biased agonists regulate an occluded conformation with a slightly outward movement of TM6 and an inward shift of TM7, which should favor β-arrestin signaling. The balanced agonist does not induce an observable outward shift for TM6 but, along with a slight tilt for TM7, leads to an inactive-like conformation. In addition, our results reveal the first time that ICL3 presents specific conformations with different agonists. The G-protein-biased agonist drives ICL3 to open so that the G protein-binding pocket can be available, while the β-arrestin-biased agonists induce ICL3 to form a closed conformation with a stable local α-helix. MM/PBSA analysis further reveals that the hydroxyl groups in the resorcinol of the G-protein-biased agonist form strong interactions with Y5.38 and S5.42, thus preventing tilting of the TM5 extracellular end. The catechol of the balanced agonist and the β-arrestin-biased ones induces the rearrangement of two hydrophobic residues F6.52 and W6.48. However, different from the balanced agonist, the ethyl substituent of β-arrestin-biased agonists forms additional hydrophobic interactions with W6.48 and F6.51 after the rearrangement, which should contribute to the β-arrestin bias. The shortest pathway analysis further reveals that the three residues Y7.43, N7.45, and N7.49 are crucial for allosterically regulating G-protein-biased signaling, while the two residues W6.48 and F6.44 make an important contribution to regulate β-arrestin-biased signaling. For the balanced agonist NE, the allosteric regulation pathway simultaneously involves the residue associated with G-protein-biased signaling like S5.46 and the residues related to β-arrestin-biased signaling like W6.48 and F6.44, thus producing unbiased signaling. The observations could advance our understanding of the biased activation mechanism on class A GPCRs and provide a useful guideline for the design of biased drugs.

Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation

The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W318^7.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W293^6.48 and F289^6.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y336^7.53, F343^8.50, and D340^8.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L158^3.43 and V286^6.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.

Exploring the Molecular Basis of Substrate and Product Selectivities of Nocardicin Bifunctional Thioesterase

D-amino acid introduction in peptides can enrich their biological activities and pharmacological properties as potential drugs. This achievement of stereochemical inversion usually owes to an epimerase or racemase. Interestingly, a unique bifunctional thioesterase (NocTE), which is incorporated in nonribosomal peptide synthetase (NRPS) NocA-NocB assembly line for the biosynthesis of monocyclic β-lactam antibiotic nocardicin A, can control the generation of D-products with high stereochemical purity. However, the molecular basis of NocTE selectivity on substrates and products is still unclear. Herein, we constructed a series of systems with different peptides varying in stereochemistry, length, and composition to investigate the substrate selectivity. The studies on binding affinities and loading conformations elucidated the important roles of peptide length and β-lactam ring in substrate selectivity. Through energy decomposition and interaction analyses, some key residues involved in substrate selectivity were captured. On the other hand, natural product undergoing epimerization was found to be liberated from the active pocket more easily in comparison with its diastereomer (epi-nocardicin G), explaining the superiority of nocardicin G. These results provide detailed molecular insights into the exquisite control of substrate and product scopes for NocTE, and encourage to diversification of substrates and final products for NRPS assembly line. The molecular insights into substrate and product selectivities of unique bifunctional thioesterase NocTE were illustrated via several molecular simulations and free energy calculations, contributing to expanding substrate and product scopes of nonribosomal peptide synthetases.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.