How GPCR Phosphorylation Patterns Orchestrate Arrestin-Mediated Signaling.

Binding of arrestin to phosphorylated G-protein-coupled receptors (GPCRs) controls many aspects of cell signaling. The number and arrangement of phosphates may vary substantially for a given GPCR, and different phosphorylation patterns trigger different arrestin-mediated effects. Here, we determine how GPCR phosphorylation influences arrestin behavior by using atomic-level simulations and site-directed spectroscopy to reveal the effects of phosphorylation patterns on arrestin binding and conformation. We find that patterns favoring binding differ from those favoring activation-associated conformational change. Both binding and conformation depend more on arrangement of phosphates than on their total number, with phosphorylation at different positions sometimes exerting opposite effects. Phosphorylation patterns selectively favor a wide variety of arrestin conformations, differently affecting arrestin sites implicated in scaffolding distinct signaling proteins. We also reveal molecular mechanisms of these phenomena. Our work reveals the structural basis for the long-standing "barcode" hypothesis and has important implications for design of functionally selective GPCR-targeted drugs.

Structure of a Signaling Cannabinoid Receptor 1-G Protein Complex.

Cannabis elicits its mood-enhancing and analgesic effects through the cannabinoid receptor 1 (CB1), a G protein-coupled receptor (GPCR) that signals primarily through the adenylyl cyclase-inhibiting heterotrimeric G protein Gi. Activation of CB1-Gi signaling pathways holds potential for treating a number of neurological disorders and is thus crucial to understand the mechanism of Gi activation by CB1. Here, we present the structure of the CB1-Gi signaling complex bound to the highly potent agonist MDMB-Fubinaca (FUB), a recently emerged illicit synthetic cannabinoid infused in street drugs that have been associated with numerous overdoses and fatalities. The structure illustrates how FUB stabilizes the receptor in an active state to facilitate nucleotide exchange in Gi. The results compose the structural framework to explain CB1 activation by different classes of ligands and provide insights into the G protein coupling and selectivity mechanisms adopted by the receptor.

Structure and mechanism of the cation-chloride cotransporter NKCC1.

Cation-chloride cotransporters (CCCs) mediate the electroneutral transport of chloride, potassium and/or sodium across the membrane. They have critical roles in regulating cell volume, controlling ion absorption and secretion across epithelia, and maintaining intracellular chloride homeostasis. These transporters are primary targets for some of the most commonly prescribed drugs. Here we determined the cryo-electron microscopy structure of the Na-K-Cl cotransporter NKCC1, an extensively studied member of the CCC family, from Danio rerio. The structure defines the architecture of this protein family and reveals how cytosolic and transmembrane domains are strategically positioned for communication. Structural analyses, functional characterizations and computational studies reveal the ion-translocation pathway, ion-binding sites and key residues for transport activity. These results provide insights into ion selectivity, coupling and translocation, and establish a framework for understanding the physiological functions of CCCs and interpreting disease-related mutations.

Molecular mechanism of GPCR-mediated arrestin activation.

Despite intense interest in discovering drugs that cause G-protein-coupled receptors (GPCRs) to selectively stimulate or block arrestin signalling, the structural mechanism of receptor-mediated arrestin activation remains unclear1,2. Here we reveal this mechanism through extensive atomic-level simulations of arrestin. We find that the receptor's transmembrane core and cytoplasmic tail-which bind distinct surfaces on arrestin-can each independently stimulate arrestin activation. We confirm this unanticipated role of the receptor core, and the allosteric coupling between these distant surfaces of arrestin, using site-directed fluorescence spectroscopy. The effect of the receptor core on arrestin conformation is mediated primarily by interactions of the intracellular loops of the receptor with the arrestin body, rather than the marked finger-loop rearrangement that is observed upon receptor binding. In the absence of a receptor, arrestin frequently adopts active conformations when its own C-terminal tail is disengaged, which may explain why certain arrestins remain active long after receptor dissociation. Our results, which suggest that diverse receptor binding modes can activate arrestin, provide a structural foundation for the design of functionally selective ('biased') GPCR-targeted ligands with desired effects on arrestin signalling.

Leveraging nonstructural data to predict structures and affinities of protein-ligand complexes.

Over the past five decades, tremendous effort has been devoted to computational methods for predicting properties of ligands-i.e., molecules that bind macromolecular targets. Such methods, which are critical to rational drug design, fall into two categories: physics-based methods, which directly model ligand interactions with the target given the target's three-dimensional (3D) structure, and ligand-based methods, which predict ligand properties given experimental measurements for similar ligands. Here, we present a rigorous statistical framework to combine these two sources of information. We develop a method to predict a ligand's pose-the 3D structure of the ligand bound to its target-that leverages a widely available source of information: a list of other ligands that are known to bind the same target but for which no 3D structure is available. This combination of physics-based and ligand-based modeling improves pose prediction accuracy across all major families of drug targets. Using the same framework, we develop a method for virtual screening of drug candidates, which outperforms standard physics-based and ligand-based virtual screening methods. Our results suggest broad opportunities to improve prediction of various ligand properties by combining diverse sources of information through customized machine-learning approaches.

Delineating the Ligand-Receptor Interactions That Lead to Biased Signaling at the µ-Opioid Receptor.

Biased agonists, which selectively stimulate certain signaling pathways controlled by a G protein-coupled receptor (GPCR), hold great promise as drugs that maximize efficacy while minimizing dangerous side effects. Biased agonists of the µ-opioid receptor (µOR) are of particular interest as a means to achieve analgesia through G protein signaling without dose-limiting side effects such as respiratory depression and constipation. Rational structure-based design of biased agonists remains highly challenging, however, because the ligand-mediated interactions that are key to activation of each signaling pathway remain unclear. We identify several compounds for which the R- and S-enantiomers have distinct bias profiles at the µOR. These compounds serve as excellent comparative tools to study bias because the identical physicochemical properties of enantiomer pairs ensure that differences in bias profiles are due to differences in interactions with the µOR binding pocket. Atomic-level simulations of compounds at µOR indicate that R- and S-enantiomers adopt different poses that form distinct interactions with the binding pocket. A handful of specific interactions with highly conserved binding pocket residues appear to be responsible for substantial differences in arrestin recruitment between enantiomers. Our results offer guidance for rational design of biased agonists at µOR and possibly at related GPCRs.

Gi- and Gs-coupled GPCRs show different modes of G-protein binding.

More than two decades ago, the activation mechanism for the membrane-bound photoreceptor and prototypical G protein-coupled receptor (GPCR) rhodopsin was uncovered. Upon light-induced changes in ligand-receptor interaction, movement of specific transmembrane helices within the receptor opens a crevice at the cytoplasmic surface, allowing for coupling of heterotrimeric guanine nucleotide-binding proteins (G proteins). The general features of this activation mechanism are conserved across the GPCR superfamily. Nevertheless, GPCRs have selectivity for distinct G-protein family members, but the mechanism of selectivity remains elusive. Structures of GPCRs in complex with the stimulatory G protein, Gs, and an accessory nanobody to stabilize the complex have been reported, providing information on the intermolecular interactions. However, to reveal the structural selectivity filters, it will be necessary to determine GPCR-G protein structures involving other G-protein subtypes. In addition, it is important to obtain structures in the absence of a nanobody that may influence the structure. Here, we present a model for a rhodopsin-G protein complex derived from intermolecular distance constraints between the activated receptor and the inhibitory G protein, Gi, using electron paramagnetic resonance spectroscopy and spin-labeling methodologies. Molecular dynamics simulations demonstrated the overall stability of the modeled complex. In the rhodopsin-Gi complex, Gi engages rhodopsin in a manner distinct from previous GPCR-Gs structures, providing insight into specificity determinants.

Conformational selectivity in cytochrome P450 redox partner interactions.

The heme iron of cytochromes P450 must be reduced to bind and activate molecular oxygen for substrate oxidation. Reducing equivalents are derived from a redox partner, which requires the formation of a protein-protein complex. A subject of increasing discussion is the role that redox partner binding plays, if any, in favoring significant structural changes in the P450s that are required for activity. Many P450s now have been shown to experience large open and closed motions. Several structural and spectral studies indicate that the well-studied P450cam adopts the open conformation when its redox partner, putidaredoxin (Pdx), binds, whereas recent NMR studies indicate that this view is incorrect. Given the relevance of this discrepancy to P450 chemistry, it is important to determine whether Pdx favors the open or closed form of P450cam. Here, we have used both computational and experimental isothermal titration calorimetry studies that unequivocally show Pdx favors binding to the open form of P450cam. Analyses of molecular-dynamic trajectories also provide insights into intermediate conformational states that could be relevant to catalysis.

Insights into the Dynamics and Dissociation Mechanism of a Protein Redox Complex Using Molecular Dynamics.

Leishmania major peroxidase (LmP) is structurally and functionally similar to the well-studied yeast Cytochrome c peroxidase (CCP). A recent Brownian dynamics study showed that L. major Cytochrome c (LmCytc) associates with LmP by forming an initial complex with the N-terminal helix A of LmP, followed by a movement toward the electron transfer (ET) site observed in the LmP-LmCytc crystal structure. Critical to forming the active electron transfer complex is an intermolecular Arg-Asp ion pair at the center of the interface. If the dissociation reaction is effectively the reverse of the association reaction, then rupture of the Asp-Arg ion pair should be followed by movement of LmCytc back toward LmP helix A. To test this possibility, we have performed multiple molecular dynamics (MD) simulations of the LmP-LmCytc complex. In five separate simulations, LmCytc is observed to indeed move toward helix A, and in two of the simulations, the Asp-Arg ion pair breaks, which frees LmCytc to fully associate with the LmP helix A secondary binding site. These results support the "bind and crawl" or "velcro" mechanism of association, wherein LmCytc forms a nonspecific electrostatic complex with LmP helix A, followed by a "crawl" toward the ET-active site, where the Asp-Arg ion pair holds the LmCytc in position for rapid ET. These simulations also point to Tyr134LmP as being important in the association/dissociation reactions. Experimentally mutating Tyr134 to Phe was found to decrease Km by 3.6-fold, which is consistent with its predicted role in complex formation by MD simulations.

"Bind and Crawl" Association Mechanism of Leishmania major Peroxidase and Cytochrome c Revealed by Brownian and Molecular Dynamics Simulations.

Leishmania major, the parasitic causative agent of leishmaniasis, produces a heme peroxidase (LmP), which catalyzes the peroxidation of mitochondrial cytochrome c (LmCytc) for protection from reactive oxygen species produced by the host. The association of LmP and LmCytc, which is known from kinetics measurements to be very fast (∼10(8) M(-1) s(-1)), does not involve major conformational changes and has been suggested to be dominated by electrostatic interactions. We used Brownian dynamics simulations to investigate the mechanism of formation of the LmP-LmCytc complex. Our simulations confirm the importance of electrostatic interactions involving the negatively charged D211 residue at the LmP active site, and reveal a previously unrecognized role in complex formation for negatively charged residues in helix A of LmP. The crystal structure of the D211N mutant of LmP reported herein is essentially identical to that of wild-type LmP, reinforcing the notion that it is the loss of charge at the active site, and not a change in structure, that reduces the association rate of the D211N variant of LmP. The Brownian dynamics simulations further show that complex formation occurs via a "bind and crawl" mechanism, in which LmCytc first docks to a location on helix A that is far from the active site, forming an initial encounter complex, and then moves along helix A to the active site. An atomistic molecular dynamics simulation confirms the helix A binding site, and steady state activity assays and stopped-flow kinetics measurements confirm the role of helix A charges in the association mechanism.

Enzymatic Mechanism of Leishmania major Peroxidase and the Critical Role of Specific Ionic Interactions.

Leishmania major peroxidase (LmP) is very similar to the well-known yeast cytochrome c peroxidase (CcP). Both enzymes catalyze the peroxidation of cytochrome c. Like CcP, LmP reacts with H2O2 to form Compound I, which consists of a ferryl heme and a Trp radical, Fe(IV)═O;Trp(•+). Cytochrome c (Cytc) reduces the Trp radical to give Compound II, Fe(IV)═O;Trp, which is followed by an intramolecular electron transfer to give Fe(III)-OH;Trp(•+), and in the last step, Cytc reduces the Trp radical. In this study, we have used steady-state and single-turnover kinetics to improve our understanding of the overall mechanism of LmP catalysis. While the activity of CcP greatly increases with ionic strength, the kcat for LmP remains relatively constant at all ionic strengths tested. Therefore, unlike CcP, where dissociation of oxidized Cytc is limiting at low ionic strengths, association/dissociation reactions are not limiting at any ionic strength in LmP. We conclude that in LmP, the intramolecular electron transfer reaction, Fe(IV)═O;Trp to Fe(III)-OH;Trp(•+), is limiting at all ionic strengths. Unlike CcP, LmP depends on key intermolecular ion pairs to form the electron transfer competent complex. Mutating these sites causes the initial rate of association to decrease by 2 orders of magnitude and a substantial decrease in kcat. The drop in kcat is due to a switch in the rate-limiting step of the mutants from intramolecular electron transfer to the rate of association in forming the LmP-LmCytc complex. These studies show that while LmP and CcP form very similar complexes and exhibit similar activities, they substantially differ in how their activity changes as a function of ionic strength. This difference is primarily due to the heavy reliance of LmP on highly specific intermolecular ion pairs, while CcP relies mainly on nonpolar interactions.

Crystal structure of cindoxin, the P450cin redox partner.

The crystal structure of the flavin mononucleotide (FMN)-containing redox partner to P450cin, cindoxin (Cdx), has been determined to 1.3 Å resolution. The overall structure is similar to that of the FMN domain of human cytochrome P450 reductase. A Brownian dynamics-molecular dynamics docking method was used to produce a model of Cdx with its redox partner, P450cin. This Cdx-P450cin model highlights the potential importance of Cdx Tyr96 in bridging the FMN and heme cofactors as well P450cin Arg102 and Arg346. Each of the single-site Ala mutants exhibits ~10% of the wild-type activity, thus demonstrating the importance of these residues for binding and/or electron transfer. In the well-studied P450cam system, redox partner binding stabilizes the open low-spin conformation of P450cam and greatly decreases the stability of the oxy complex. In sharp contrast, Cdx does not shift P450cin to a low-spin state, although the stability of oxy-P450cin is decreased 10-fold in the presence of Cdx. This indicates that Cdx may have a modest effect on the open-closed equilibrium in P450cin compared to that in P450cam. It has been postulated that part of the effector role of Pdx on P450cam is to promote a significant structural change that makes available a proton relay network involving Asp251 required for O2 activation. The structure around the corresponding Asp in P450cin, Asp241, provides a possible structural reason for why P450cin is less dependent on its redox partner for functionally important structural changes.

Interactions of gaseous HNO3 and water with individual and mixed alkyl self-assembled monolayers at room temperature.

The major removal processes for gaseous nitric acid (HNO3) in the atmosphere are dry and wet deposition onto various surfaces. The surface in the boundary layer is often covered with organic films, but the interaction of gaseous HNO3 with them is not well understood. To better understand the factors controlling the uptake of gaseous nitric acid and its dissociation in organic films, studies were carried out using single component and mixtures of C8 and C18 alkyl self-assembled monolayers (SAMs) attached to a germanium (Ge) attenuated total reflectance (ATR) crystal upon which a thin layer of SiOx had been deposited. For comparison, diffuse reflectance infrared Fourier transform spectrometry (DRIFTS) studies were also carried out using a C18 SAM attached to the native oxide layer on the surface of silicon powder. These studies show that the alkyl chain length and order/disorder of the SAMs does not significantly affect the uptake or dissociation/recombination of molecular HNO3. Thus, independent of the nature of the SAM, molecular HNO3 is observed up to 70-90% relative humidity. After dissociation, molecular HNO3 is regenerated on all SAM surfaces when water is removed. Results of molecular dynamics simulations are consistent with experiments and show that defects and pores on the surfaces control the uptake, dissociation and recombination of molecular HNO3. Organic films on surfaces in the boundary layer will certainly be more irregular and less ordered than SAMs studied here, therefore undissociated HNO3 may be present on surfaces in the boundary layer to a greater extent than previously thought. The combination of this observation with the results of recent studies showing enhanced photolysis of nitric acid on surfaces suggests that renoxification of deposited nitric acid may need to be taken into account in atmospheric models.

Design and Optimization of Selectivity-Tunable Toll-like Receptor 7/8 Agonists as Novel Antibody-Drug Conjugate Payloads.

Toll-like receptors 7 and 8 are involved in modulating the adaptive and innate immune responses, and their activation has shown promise as a therapeutic strategy in the field of immuno-oncology. While systemic exposure to TLR7/8 agonists can result in poor tolerance, combination therapies and targeted delivery through antibody-drug conjugates (ADCs) can help mitigate adverse effects. Described herein is the identification of a novel and potent series of pyrazolopyrimidine-based TLR7/8 agonists with tunable receptor selectivity. Representative agonists from this series were successfully able to induce the production of various proinflammatory cytokines and chemokines from human peripheral blood mononuclear cells. Anti-HER2-25 and anti-HER2-26 ADCs made from this class of payloads demonstrated mechanism-based activation of TLR7/8 in a THP1/N87 coculture system.

P450cin active site water: implications for substrate binding and solvent accessibility.

In P450cin, Tyr81, Asp241, Asn242, two water molecules, and the substrate participate in a complex H-bonded network. The role of this H-bonded network in substrate binding and catalysis has been probed by crystallography, spectroscopy, kinetics, isothermal titration calorimetry (ITC), and molecular dynamics. For the Y81F mutant, the substrate binds about 20-fold more weakly and Vmax decreases by about 30% in comparison to WT. The enhanced susceptibility of the heme to H2O2-mediated destruction in Y81F suggests that this mutant favors the open, low-spin conformational state. Asn242 H-bonds directly with the substrate, and replacing this residue with Ala results in water taking the place of the missing Asn side chain. This mutant exhibits a 70% decrease in activity. Crystal structures and molecular dynamics simulations of substrate-bound complexes show that the solvent has more ready access to the active site, especially for the N242A mutant. This accounts for about a 64% uncoupling of electron transfer from substrate hydroxylation. These data indicate the importance of the interconnected water network on substrate binding and on the open/closed conformational equilibrium, which are both critically important for maintaining high-coupling efficiency.

Discovery and multimerization of cross-reactive single-domain antibodies against SARS-like viruses to enhance potency and address emerging SARS-CoV-2 variants.

Coronaviruses have been the causative agent of three epidemics and pandemics in the past two decades, including the ongoing COVID-19 pandemic. A broadly-neutralizing coronavirus therapeutic is desirable not only to prevent and treat COVID-19, but also to provide protection for high-risk populations against future emergent coronaviruses. As all coronaviruses use spike proteins on the viral surface to enter the host cells, and these spike proteins share sequence and structural homology, we set out to discover cross-reactive biologic agents targeting the spike protein to block viral entry. Through llama immunization campaigns, we have identified single domain antibodies (VHHs) that are cross-reactive against multiple emergent coronaviruses (SARS-CoV, SARS-CoV-2, and MERS). Importantly, a number of these antibodies show sub-nanomolar potency towards all SARS-like viruses including emergent CoV-2 variants. We identified nine distinct epitopes on the spike protein targeted by these VHHs. Further, by engineering VHHs targeting distinct, conserved epitopes into multi-valent formats, we significantly enhanced their neutralization potencies compared to the corresponding VHH cocktails. We believe this approach is ideally suited to address both emerging SARS-CoV-2 variants during the current pandemic as well as potential future pandemics caused by SARS-like coronaviruses.

Hexamerization of Anti-SARS CoV IgG1 Antibodies Improves Neutralization Capacity.

The SARS-CoV-2 pandemic and particularly the emerging variants have deepened the need for widely available therapeutic options. We have demonstrated that hexamer-enhancing mutations in the Fc region of anti-SARS-CoV IgG antibodies lead to a noticeable improvement in IC50 in both pseudo and live virus neutralization assay compared to parental molecules. We also show that hexamer-enhancing mutants improve C1q binding to target surface. To our knowledge, this is the first time this format has been explored for application in viral neutralization and the studies provide proof-of-concept for the use of hexamer-enhanced IgG1 molecules as potential anti-viral therapeutics.

Functionally selective signaling and broad metabolic benefits by novel insulin receptor partial agonists.

Insulin analogs have been developed to treat diabetes with focus primarily on improving the time action profile without affecting ligand-receptor interaction or functional selectivity. As a result, inherent liabilities (e.g. hypoglycemia) of injectable insulin continue to limit the true therapeutic potential of related agents. Insulin dimers were synthesized to investigate whether partial agonism of the insulin receptor (IR) tyrosine kinase is achievable, and to explore the potential for tissue-selective systemic insulin pharmacology. The insulin dimers induced distinct IR conformational changes compared to native monomeric insulin and substrate phosphorylation assays demonstrated partial agonism. Structurally distinct dimers with differences in conjugation sites and linkers were prepared to deliver desirable IR partial agonist (IRPA). Systemic infusions of a B29-B29 dimer in vivo revealed sharp differences compared to native insulin. Suppression of hepatic glucose production and lipolysis were like that attained with regular insulin, albeit with a distinctly shallower dose-response. In contrast, there was highly attenuated stimulation of glucose uptake into muscle. Mechanistic studies indicated that IRPAs exploit tissue differences in receptor density and have additional distinctions pertaining to drug clearance and distribution. The hepato-adipose selective action of IRPAs is a potentially safer approach for treatment of diabetes.

Discovery of Insulin Receptor Partial Agonists MK-5160 and MK-1092 as Novel Basal Insulins with Potential to Improve Therapeutic Index.

We have identified a series of novel insulin receptor partial agonists (IRPAs) with a potential to mitigate the risk of hypoglycemia associated with the use of insulin as an antidiabetic treatment. These molecules were designed as dimers of native insulin connected via chemical linkers of variable lengths with optional capping groups at the N-terminals of insulin chains. Depending on the structure, the maximal activation level (%Max) varied in the range of ∼20-70% of native insulin, and EC50 values remained in sub-nM range. Studies in minipig and dog demonstrated that IRPAs had sufficient efficacy to normalize plasma glucose levels in diabetes, while providing reduction of hypoglycemia risk. IRPAs had a prolonged duration of action, potentially making them suitable for once-daily dosing. Two lead compounds with %Max values of 30 and 40% relative to native insulin were selected for follow up studies in the clinic.

Structures of active-state orexin receptor 2 rationalize peptide and small-molecule agonist recognition and receptor activation.

Narcolepsy type 1 (NT1) is a chronic neurological disorder that impairs the brain's ability to control sleep-wake cycles. Current therapies are limited to the management of symptoms with modest effectiveness and substantial adverse effects. Agonists of the orexin receptor 2 (OX2R) have shown promise as novel therapeutics that directly target the pathophysiology of the disease. However, identification of drug-like OX2R agonists has proven difficult. Here we report cryo-electron microscopy structures of active-state OX2R bound to an endogenous peptide agonist and a small-molecule agonist. The extended carboxy-terminal segment of the peptide reaches into the core of OX2R to stabilize an active conformation, while the small-molecule agonist binds deep inside the orthosteric pocket, making similar key interactions. Comparison with antagonist-bound OX2R suggests a molecular mechanism that rationalizes both receptor activation and inhibition. Our results enable structure-based discovery of therapeutic orexin agonists for the treatment of NT1 and other hypersomnia disorders.