Structural insight into allosteric modulation of protease-activated receptor 2.

Protease-activated receptors (PARs) are a family of G-protein-coupled receptors (GPCRs) that are irreversibly activated by proteolytic cleavage of the N terminus, which unmasks a tethered peptide ligand that binds and activates the transmembrane receptor domain, eliciting a cellular cascade in response to inflammatory signals and other stimuli. PARs are implicated in a wide range of diseases, such as cancer and inflammation. PARs have been the subject of major pharmaceutical research efforts but the discovery of small-molecule antagonists that effectively bind them has proved challenging. The only marketed drug targeting a PAR is vorapaxar, a selective antagonist of PAR1 used to prevent thrombosis. The structure of PAR1 in complex with vorapaxar has been reported previously. Despite sequence homology across the PAR isoforms, discovery of PAR2 antagonists has been less successful, although GB88 has been described as a weak antagonist. Here we report crystal structures of PAR2 in complex with two distinct antagonists and a blocking antibody. The antagonist AZ8838 binds in a fully occluded pocket near the extracellular surface. Functional and binding studies reveal that AZ8838 exhibits slow binding kinetics, which is an attractive feature for a PAR2 antagonist competing against a tethered ligand. Antagonist AZ3451 binds to a remote allosteric site outside the helical bundle. We propose that antagonist binding prevents structural rearrangements required for receptor activation and signalling. We also show that a blocking antibody antigen-binding fragment binds to the extracellular surface of PAR2, preventing access of the tethered ligand to the peptide-binding site. These structures provide a basis for the development of selective PAR2 antagonists for a range of therapeutic uses.

Biosensor-Enabled Deconvolution of the Avidity-Induced Affinity Enhancement for the SARS-CoV-2 Spike Protein and ACE2 Interaction.

Avidity is an effective and frequent phenomenon employed by nature to achieve extremely high-affinity interactions. As more drug discovery efforts aim to disrupt protein-protein interactions, it is becoming increasingly common to encounter systems that utilize avidity effects and to study these systems using surface-based technologies, such as surface plasmon resonance (SPR) or biolayer interferometry. However, heterogeneity introduced from multivalent binding interactions complicates the analysis of the resulting sensorgram. A frequently applied practice is to fit the data based on a 1:1 binding model, and if the fit does not describe the data adequately, then the experimental setup is changed to favor a 1:1 binding interaction. This reductionistic approach is informative but not always biologically relevant. Therefore, we aimed to develop an SPR-based assay that would reduce the heterogeneity to enable the determination of the kinetic rate constants for multivalent binding interactions using the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein and the human receptor angiotensin-converting enzyme 2 (ACE2) as a model system. We employed a combinatorial approach to generate a sensor surface that could distinguish between monovalent and multivalent interactions. Using advanced data analysis algorithms to analyze the resulting sensorgrams, we found that controlling the surface heterogeneity enabled the deconvolution of the avidity-induced affinity enhancement for the SARS-CoV-2 spike protein and ACE2 interaction.

Fragment-assisted hit investigation involving integrated HTS and fragment screening: Application to the identification of phosphodiesterase 10A (PDE10A) inhibitors.

Fragment-based drug design (FBDD) relies on direct elaboration of fragment hits and typically requires high resolution structural information to guide optimization. In fragment-assisted drug discovery (FADD), fragments provide information to guide selection and design but do not serve as starting points for elaboration. We describe FADD and high-throughput screening (HTS) campaign strategies conducted in parallel against PDE10A where fragment hit co-crystallography was not available. The fragment screen led to prioritized fragment hits (IC50's ∼500µM), which were used to generate a hypothetical core scaffold. Application of this scaffold as a filter to HTS output afforded a 4µM hit, which, after preparation of a small number of analogs, was elaborated into a 16nM lead. This approach highlights the strength of FADD, as fragment methods were applied despite the absence of co-crystallographical information to efficiently identify a lead compound for further optimization.

Gremlin1 preferentially binds to bone morphogenetic protein-2 (BMP-2) and BMP-4 over BMP-7.

Gremlin (Grem1) is a member of the DAN family of secreted bone morphogenetic protein (BMP) antagonists. Bone morphogenetic protein-7 (BMP-7) mediates protective effects during renal fibrosis associated with diabetes and other renal diseases. The pathogenic mechanism of Grem1 during diabetic nephropathy (DN) has been suggested to be binding and inhibition of BMP-7. However, the precise interactions between Grem1, BMP-7 and other BMPs have not been accurately defined. In the present study, we show the affinity of Grem1 for BMP-7 is lower than that of BMP-2 and BMP-4, using a combination of surface plasmon resonance and cell culture techniques. Using kidney proximal tubule cells and HEK (human embryonic kidney)-293 cell Smad1/5/8 phosphorylation and BMP-dependent gene expression as readouts, Grem1 consistently demonstrated a higher affinity for BMP-2>BMP-4>BMP-7. Cell-associated Grem1 did not inhibit BMP-2- or BMP-4-mediated signalling, suggesting that Grem1-BMP-2 binding occurred in solution, preventing BMP receptor activation. These data suggest that Grem1 preferentially binds to BMP-2 and this may be the dominant complex in a disease situation where levels of Grem1 and BMPs are elevated.

Crystal structure of microsomal prostaglandin E2 synthase provides insight into diversity in the MAPEG superfamily.

Prostaglandin E2 (PGE2) is a key mediator in inflammatory response. The main source of inducible PGE2, microsomal PGE2 synthase-1 (mPGES-1), has emerged as an interesting drug target for treatment of pain. To support inhibitor design, we have determined the crystal structure of human mPGES-1 to 1.2 Å resolution. The structure reveals three well-defined active site cavities within the membrane-spanning region in each monomer interface of the trimeric structure. An important determinant of the active site cavity is a small cytosolic domain inserted between transmembrane helices I and II. This extra domain is not observed in other structures of proteins within the MAPEG (Membrane-Associated Proteins involved in Eicosanoid and Glutathione metabolism) superfamily but is likely to be present also in microsomal GST-1 based on sequence similarity. An unexpected feature of the structure is a 16-Å-deep cone-shaped cavity extending from the cytosolic side into the membrane-spanning region. We suggest a potential role for this cavity in substrate access. Based on the structure of the active site, we propose a catalytic mechanism in which serine 127 plays a key role. We have also determined the structure of mPGES-1 in complex with a glutathione-based analog, providing insight into mPGES-1 flexibility and potential for structure-based drug design.

Drug discovery at the single molecule level: inhibition-in-solution assay of membrane-reconstituted ß-secretase using single-molecule imaging.

Inhibition-in-solution assays (ISA) employing surface-based biosensors such as surface plasmon resonance (SPR) are an effective screening approach in drug discovery. However, analysis of potent binders remains a significant hurdle due to limited sensitivity and accompanied depletion of the inhibiting compounds due to high protein concentrations needed for detectable binding signals. To overcome this limitation, we explored a microscopy-based single-molecule ISA compatible with liposome-reconstituted membrane proteins. Using a set of validated small molecule inhibitors against ß-secretase 1 (BACE1), the assay was benchmarked with respect to sensitivity and dynamic range against SPR. We demonstrate that the dynamic range of measurable affinities is greatly extended by more than 2 orders of magnitude as compared to SPR, thus facilitating measurements of highly potent (Kd < nM) compounds.

Preserved transmembrane protein mobility in polymer-supported lipid bilayers derived from cell membranes.

Supported lipid bilayers (SLBs) have contributed invaluable information about the physiochemical properties of cell membranes, but their compositional simplicity often limits the level of knowledge that can be gained about the structure and function of transmembrane proteins in their native environment. Herein, we demonstrate a generic protocol for producing polymer-supported lipid bilayers on glass surfaces that contain essentially all naturally occurring cell-membrane components of a cell line while still retaining transmembrane protein mobility and activity. This was achieved by merging vesicles made from synthetic lipids (PEGylated lipids and POPC lipids) with native cell-membrane vesicles to generate hybrid vesicles which readily rupture into a continuous polymer-supported lipid bilayer. To investigate the properties of these complex hybrid SLBs and particularly the behavior of their integral membrane-proteins, we used total internal reflection fluorescence imaging to study a transmembrane protease, ß-secretase 1 (BACE1), whose ectoplasmic and cytoplasmic domains could both be specifically targeted with fluorescent reporters. By selectively probing the two different orientations of BACE1 in the resulting hybrid SLBs, the role of the PEG-cushion on transmembrane protein lateral mobility was investigated. The results reveal the necessity of having the PEGylated lipids present during vesicle adsorption to prevent immobilization of transmembrane proteins with protruding domains. The proteolytic activity of BACE1 was unadulterated by the sonication process used to merge the synthetic and native membrane vesicles; importantly it was also conserved in the SLB. The presented strategy could thus serve both fundamental studies of membrane biophysics and the production of surface-based bioanalytical sensor platforms.

Equilibrium-fluctuation analysis for interaction studies between natural ligands and single G protein-coupled receptors in native lipid vesicles.

G protein-coupled receptors (GPCRs) constitute the most versatile family of cell-membrane receptors and have been increasingly identified as important mediators of many physiological functions. They also belong to one of the most central drug target classes, but current screening technologies are limited by the requirements of overexpression and stabilization of GPCRs. This calls for sensitivity-increased detection strategies preferably meeting single-molecule detection limits. This challenge is here addressed by employing total internal reflection fluorescence microscopy to characterize the interaction kinetics between CXCR3, a GPCR involved in inflammatory responses, and two of its chemokine ligands, CXCL10 and CXCL11. Fluorescence labeling of the lipid membrane, rather than the membrane protein itself, of GPCR-containing native vesicles, and immobilization of the corresponding ligand on the surface, enabled determination of the interaction kinetics using single-molecule equilibrium-fluctuation analysis. With a limit of detection of GPCR-containing vesicles in the low picomolar concentration regime, the results demonstrate the possibility to use inhibition in solution screening of high affinity ligands/drug candidates, which due to target-binding depletion of the inhibiting compounds is demanding using assays with more moderate detection limits.

Characterization of a small molecule inhibitor of plasminogen activator inhibitor type 1 that accelerates the transition into the latent conformation.

A novel class of small molecule inhibitors for plasminogen activator inhibitor type 1 (PAI-1), represented by AZ3976, was identified in a high throughput screening campaign. AZ3976 displayed an IC(50) value of 26 µm in an enzymatic chromogenic assay. In a plasma clot lysis assay, the compound was active with an IC(50) of 16 µm. Surprisingly, AZ3976 did not bind to active PAI-1 but bound to latent PAI-1 with a K(D) of 0.29 µm at 35 °C and a binding stoichiometry of 0.94, as measured by isothermal calorimetry. Reversible binding was confirmed by surface plasmon resonance direct binding experiments. The x-ray structure of AZ3976 in complex with latent PAI-1 was determined at 2.4 Å resolution. The inhibitor was bound in the flexible joint region with the entrance to the cavity located between α-helix D and ß-strand 2A. A set of surface plasmon resonance experiments revealed that AZ3976 inhibited PAI-1 by enhancing the latency transition of active PAI-1. Because AZ3976 only had measurable affinity for latent PAI-1, we propose that its mechanism of inhibition is based on binding to a small fraction in equilibrium with active PAI-1, a latent-like prelatent form, from which latent PAI-1 is then generated more rapidly. This mode of action, with induced accelerated latency transition of active PAI-1 may, together with supporting x-ray data, provide improved opportunities for small molecule drug design in the hunt for therapeutically useful PAI-1 inhibitors.

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations.

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Robust Prediction of Relative Binding Energies for Protein-Protein Complex Mutations Using Free Energy Perturbation Calculations.

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how Protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Identification of structural-kinetic and structural-thermodynamic relationships for thrombin inhibitors.

To improve our understanding of drug-target interactions, we explored the effect of introducing substituted amine residues with increased chain length in the P3 residue of the thrombin inhibitor melagatran. Inhibition, kinetic, and thermodynamic data obtained via stopped-flow spectroscopy (SF), isothermal microcalorimetry (ITC), and surface plasmon resonance (SPR) biosensor analysis were interpreted with the help of X-ray crystal structures of the enzyme-inhibitor complexes. The association rate became faster when the lipophilicity of the inhibitors was increased. This was coupled to an increased enthalpic component and a corresponding decreased entropic component. The dissociation rates were reduced with an increase in chain length, with only a smaller increase and a decrease in the enthalpic and entropic components, respectively. Overall, the affinity increased with an increase in chain length, with similar changes in the enthalpic and entropic components. ITC analysis confirmed the equilibrium data from SPR analysis, showing that the interaction of melagatran was the most enthalpy-driven interaction. Structural analysis of the thrombin-inhibitor complex showed that the orientation of the P1 and P2 parts of the molecules was very similar, but that there were significant differences in the interaction between the terminal part of the P3 side chain and the binding pocket. A combination of charge repulsion, H-bonds, and hydrophobic interactions could be used to explain the observed kinetic and thermodynamic profiles for the ligands. In conclusion, changes in the structure of a lead compound can have significant effects on its interaction with the target that translate directly into kinetic and thermodynamic effects. In contrast to what may be intuitively expected, hydrogen bond formation and breakage are not necessarily reflected in enthalpy gains and losses, respectively.

An enhanced biophysical screening strategy to investigate the affinity of ASOs for their target RNA.

The recent and rapid increase in the discovery of new RNA therapeutics has created the perfect terrain to explore an increasing number of novel targets. In particular, antisense oligonucleotides (ASOs) have long held the promise of an accelerated and effective drug design compared to other RNA-based therapeutics. Although ASOs in silico design has advanced distinctively in the past years, especially thanks to the several predictive frameworks for RNA folding, it is somehow limited by the wide approximation of calculating sequence affinity based on RNA-RNA/DNA sequences. None of the ASO modifications are taken into consideration, losing hybridization information particularly fundamental to ASOs that elicit their function through RNase H1-mediated mechanisms. Here we present an inexpensive and enhanced biophysical screening strategy to investigate the affinity of ASOs for their target RNA using several biophysical techniques such as high throughput differential scanning fluorimetry (DSF), circular dichroism (CD), isothermal calorimetry (ITC), surface plasmon resonance (SPR) and small-angle X-ray scattering (SAXS).

Direct cell extraction of membrane proteins for structure-function analysis.

Membrane proteins are the largest group of therapeutic targets in a variety of disease areas and yet, they remain particularly difficult to investigate. We have developed a novel one-step approach for the incorporation of membrane proteins directly from cells into lipid Salipro nanoparticles. Here, with the pannexin1 channel as a case study, we demonstrate the applicability of this method for structure-function analysis using SPR and cryo-EM.

Synthesis and structure-activity relationship of 4-(1,3-benzothiazol-2-yl)-thiophene-2-sulfonamides as cyclin-dependent kinase 5 (cdk5)/p25 inhibitors.

4-(1,3-Benzothiazol-2-yl)thiophene-2-sulfonamide (4a) was found to be a moderately potent inhibitor of cyclin-dependent kinase 5 (cdk5) from a HTS screen. The synthesis and SAR around this hit is described. The X-ray coordinates of ligand 4a with cdk5 are also reported, showing an unusual binding mode to the hinge region via a water molecule.

Application of optical biosensors in small-molecule screening activities.

The last two decades have seen remarkable progress and improvements in optical biosensor systems such that those are currently seen as an important and value-adding component of modern drug screening activities. In particular the introduction of microplate-based biosensor systems holds the promise to match the required throughput without compromising on data quality thus representing a sought-after complement to traditional fluidic systems. This article aims to highlight the application of the two most prominent optical biosensor technologies, namely surface plasmon resonance (SPR) and optical waveguide grating (OWG), in small-molecule screening and will present, review and discuss the advantages and disadvantages of different assay formats on these platforms. A particular focus will be on the specific advantages of the inhibition in solution assay (ISA) format in contrast to traditional direct binding assays (DBA). Furthermore we will discuss different application areas for both fluidic as well as plate-based biosensor systems by considering the individual strength of the platforms.

Importance of Binding Site Hydration and Flexibility Revealed When Optimizing a Macrocyclic Inhibitor of the Keap1-Nrf2 Protein-Protein Interaction.

Upregulation of the transcription factor Nrf2 by inhibition of the interaction with its negative regulator Keap1 constitutes an opportunity for the treatment of disease caused by oxidative stress. We report a structurally unique series of nanomolar Keap1 inhibitors obtained from a natural product-derived macrocyclic lead. Initial exploration of the structure-activity relationship of the lead, followed by structure-guided optimization, resulted in a 100-fold improvement in inhibitory potency. The macrocyclic core of the nanomolar inhibitors positions three pharmacophore units for productive interactions with key residues of Keap1, including R415, R483, and Y572. Ligand optimization resulted in the displacement of a coordinated water molecule from the Keap1 binding site and a significantly altered thermodynamic profile. In addition, minor reorganizations of R415 and R483 were accompanied by major differences in affinity between ligands. This study therefore indicates the importance of accounting both for the hydration and flexibility of the Keap1 binding site when designing high-affinity ligands.

Mining Natural Products for Macrocycles to Drug Difficult Targets.

Lead generation for difficult-to-drug targets that have large, featureless, and highly lipophilic or highly polar and/or flexible binding sites is highly challenging. Here, we describe how cores of macrocyclic natural products can serve as a high-quality in silico screening library that provides leads for difficult-to-drug targets. Two iterative rounds of docking of a carefully selected set of natural-product-derived cores led to the discovery of an uncharged macrocyclic inhibitor of the Keap1-Nrf2 protein-protein interaction, a particularly challenging target due to its highly polar binding site. The inhibitor displays cellular efficacy and is well-positioned for further optimization based on the structure of its complex with Keap1 and synthetic access. We believe that our work will spur interest in using macrocyclic cores for in silico-based lead generation and also inspire the design of future macrocycle screening collections.

Regenerable Biosensors for Small-Molecule Kinetic Characterization Using SPR.

A key activity in small-molecule drug discovery is the characterization of compound-target interactions. Surface plasmon resonance (SPR) is a flexible technique for this purpose, with a wide affinity range (micromoles to picomoles), low protein requirements, and the ability to characterize the kinetics of compound binding. However, a key requirement of SPR is the immobilization of the target protein to the surface of the sensor chip. The most commonly used immobilization techniques (covalent immobilization, streptavidin-biotin) are irreversible in nature, which can afford excellent baseline stability but impose limitations throughput for slowly dissociating compounds or unstable targets. Reversible immobilization (e.g., His-tag-Ni-NTA) is possible but typically precludes accurate quantification of slow dissociation kinetics due to baseline drift.Here we present our investigation of three immobilization strategies (dual-His-tagged target protein, His-tagged streptavidin, and switchavidin) that combine the robustness of irreversible immobilization with the flexibility of reversible immobilization. Each has its own advantages and limitations, and while a universal immobilization procedure remains to be found, these strategies add to the immobilization toolbox that enables previously out-of-scope applications. Such applications are highlighted in two examples that greatly increased throughput for the kinetic characterization of potent kinase inhibitors and kinetic profiling of covalent inhibitors.

Affinity-based, biophysical methods to detect and analyze ligand binding to recombinant proteins: matching high information content with high throughput.

Affinity-based technologies have become impactful tools to detect, monitor and characterize molecular interactions using recombinant target proteins. This can aid the understanding of biological function by revealing mechanistic details, and even more importantly, enables the identification of new improved ligands that can modulate the biological activity of those targets in a desired fashion. The selection of the appropriate technology is a key step in that process, as each one of the currently available technologies offers a characteristic type of biophysical information about the ligand-binding event. Alongside the indisputable advantages of each of those technologies they naturally display diverse restrictions that are quite frequently related to the target system to be studied but also to the affinity, solubility and molecular size of the ligands. This paper discusses some of the theoretical and experimental aspects of the most common affinity-based methods, what type of information can be gained from each one of those approaches, and what requirements as well as limitations are expected from working with recombinant proteins on those platforms and how those can be optimally addressed.