Intramolecular Ring-Opening Decomposition of Aryl Azetidines

The ring strain present in azetidines can lead to undesired stability issues. Herein, we described a series of N-substituted azetidines which undergo an acid-mediated intramolecular ring-opening decomposition via nucleophilic attack of a pendant amide group. Studies were conducted to understand the decomposition mechanism enabling the design of stable analogues.

PF-07059013: A Noncovalent Modulator of Hemoglobin for Treatment of Sickle Cell Disease

Sickle cell disease (SCD) is a genetic disorder caused by a single point mutation (β6 Glu → Val) on the β-chain of adult hemoglobin (HbA) that results in sickled hemoglobin (HbS). In the deoxygenated state, polymerization of HbS leads to sickling of red blood cells (RBC). Several downstream consequences of polymerization and RBC sickling include vaso-occlusion, hemolytic anemia, and stroke. We report the design of a noncovalent modulator of HbS, clinical candidate PF-07059013 (23). The seminal hit molecule was discovered by virtual screening and confirmed through a series of biochemical and biophysical studies. After a significant optimization effort, we arrived at 23, a compound that specifically binds to Hb with nanomolar affinity and displays strong partitioning into RBCs. In a 2-week multiple dose study using Townes SCD mice, 23 showed a 37.8% (±9.0%) reduction in sickling compared to vehicle treated mice. 23 (PF-07059013) has advanced to phase 1 clinical trials.

NMR spectroscopy: the swiss army knife of drug discovery

Nuclear magnetic resonance (NMR) spectroscopy has evolved into a powerful tool within drug discovery over the last two decades. While traditionally being used by medicinal chemists for small molecule structure elucidation, it can also be a valuable tool for the identification of small molecules that bind to drug targets, for the characterization of target-ligand interactions and for hit-to-lead optimization. Here, we describe how NMR spectroscopy is integrated into the Pfizer drug discovery pipeline and how we utilize this approach to identify and validate initial hits and generate leads.

Aminopyrazole Carboxamide Bruton's Tyrosine Kinase Inhibitors. Irreversible to Reversible Covalent Reactive Group Tuning

Potent covalent inhibitors of Bruton's tyrosine kinase (BTK) based on an aminopyrazole carboxamide scaffold have been identified. Compared to acrylamide-based covalent reactive groups leading to irreversible protein adducts, cyanamide-based reversible-covalent inhibitors provided the highest combined BTK potency and EGFR selectivity. The cyanamide covalent mechanism with BTK was confirmed through enzyme kinetic, NMR, MS, and X-ray crystallographic studies. The lead cyanamide-based inhibitors demonstrated excellent kinome selectivity and rat pharmacokinetic properties.

An Intracellular Allosteric Modulator Binding Pocket in SK2 Ion Channels Is Shared by Multiple Chemotypes

Small conductance potassium (SK) ion channels define neuronal firing rates by conducting the after-hyperpolarization current. They are key targets in developing therapies where neuronal firing rates are dysfunctional, such as in epilepsy, Parkinson's, and amyotrophic lateral sclerosis (ALS). Here, we characterize a binding pocket situated at the intracellular interface of SK2 and calmodulin, which we show to be shared by multiple small-molecule chemotypes. Crystallization of this complex revealed that riluzole (approved for ALS) and an analog of the anti-ataxic agent (4-chloro-phenyl)-[2-(3,5-dimethyl-pyrazol-1-yl)-pyrimidin-4-yl]-amine (CyPPA) bind to and allosterically modulate via this site. Solution-state nuclear magnetic resonance demonstrates that riluzole, NS309, and CyPPA analogs bind at this bipartite pocket. We demonstrate, by patch-clamp electrophysiology, that both classes of ligand interact with overlapping but distinct residues within this pocket. These data define a clinically important site, laying the foundations for further studies of the mechanism of action of riluzole and related molecules.

Discovery of Fragment-Derived Small Molecules for in Vivo Inhibition of Ketohexokinase (KHK)

Increased fructose consumption and its subsequent metabolism have been implicated in hepatic steatosis, dyslipidemia, obesity, and insulin resistance in humans. Since ketohexokinase (KHK) is the principal enzyme responsible for fructose metabolism, identification of a selective KHK inhibitor may help to further elucidate the effect of KHK inhibition on these metabolic disorders. Until now, studies on KHK inhibition with small molecules have been limited due to the lack of viable in vivo pharmacological tools. Herein we report the discovery of 12, a selective KHK inhibitor with potency and properties suitable for evaluating KHK inhibition in rat models. Key structural features interacting with KHK were discovered through fragment-based screening and subsequent optimization using structure-based drug design, and parallel medicinal chemistry led to the identification of pyridine 12.

Inhibiting complex IL-17A and IL-17RA interactions with a linear peptide

IL-17A is a pro-inflammatory cytokine that has been implicated in autoimmune and inflammatory diseases. Monoclonal antibodies inhibiting IL-17A signaling have demonstrated remarkable efficacy, but an oral therapy is still lacking. A high affinity IL-17A peptide antagonist (HAP) of 15 residues was identified through phage-display screening followed by saturation mutagenesis optimization and amino acid substitutions. HAP binds specifically to IL-17A and inhibits the interaction of the cytokine with its receptor, IL-17RA. Tested in primary human cells, HAP blocked the production of multiple inflammatory cytokines. Crystal structure studies revealed that two HAP molecules bind to one IL-17A dimer symmetrically. The N-terminal portions of HAP form a β-strand that inserts between two IL-17A monomers while the C-terminal section forms an α helix that directly blocks IL-17RA from binding to the same region of IL-17A. This mode of inhibition suggests opportunities for developing peptide antagonists against this challenging target.

Short Hydrophobic Peptides with Cyclic Constraints Are Potent Glucagon-like Peptide-1 Receptor (GLP-1R) Agonists

Cyclic constraints are incorporated into an 11-residue analogue of the N-terminus of glucagon-like peptide-1 (GLP-1) to investigate effects of structure on agonist activity. Cyclization through linking side chains of residues 2 and 5 or 5 and 9 produced agonists at nM concentrations in a cAMP assay. 2D NMR and CD spectra revealed an N-terminal β-turn and a C-terminal helix that differentially influenced affinity and agonist potency. These structures can inform development of small molecule agonists of the GLP-1 receptor to treat type 2 diabetes.

Structural basis for AMPK activation: natural and synthetic ligands regulate kinase activity from opposite poles by different molecular mechanisms

AMP-activated protein kinase (AMPK) is a principal metabolic regulator affecting growth and response to cellular stress. Comprised of catalytic and regulatory subunits, each present in multiple forms, AMPK is best described as a family of related enzymes. In recent years, AMPK has emerged as a desirable target for modulation of numerous diseases, yet clinical therapies remain elusive. Challenges result, in part, from an incomplete understanding of the structure and function of full-length heterotrimeric complexes. In this work, we provide the full-length structure of the widely expressed α1β1γ1 isoform of mammalian AMPK, along with detailed kinetic and biophysical characterization. We characterize binding of the broadly studied synthetic activator A769662 and its analogs. Our studies follow on the heels of the recent disclosure of the α2β1γ1 structure and provide insight into the distinct molecular mechanisms of AMPK regulation by AMP and A769662.

Optimizing glucokinase activator binding kinetics to lower in vivo hypoglycemia risk

GK activators represent a promising strategy for treatment of T2DM; however, drug candidates have failed in clinical trials due to narrow TI between efficacy and hypoglycemia.

Spirocyclic sulfamides as β-secretase 1 (BACE-1) inhibitors for the treatment of Alzheimer's disease: utilization of structure based drug design, WaterMap, and CNS penetration studies to identify centrally efficacious inhibitors

β-Secretase 1 (BACE-1) is an attractive therapeutic target for the treatment and prevention of Alzheimer's disease (AD). Herein, we describe the discovery of a novel class of BACE-1 inhibitors represented by sulfamide 14g, using a medicinal chemistry strategy to optimize central nervous system (CNS) penetration by minimizing hydrogen bond donors (HBDs) and reducing P-glycoprotein (P-gp) mediated efflux. We have also taken advantage of the combination of structure based drug design (SBDD) to guide the optimization of the sulfamide analogues and the in silico tool WaterMap to explain the observed SAR. Compound 14g is a potent inhibitor of BACE-1 with excellent permeability and a moderate P-gp liability. Administration of 14g to mice produced a significant, dose-dependent reduction in central Aβ(X-40) levels at a free drug exposure equivalent to the whole cell IC(50) (100 nM). Furthermore, studies of the P-gp knockout mouse provided evidence that efflux transporters affected the amount of Aβ lowering versus that observed in wild-type (WT) mouse at an equivalent dose.

Design of a multi-purpose fragment screening library using molecular complexity and orthogonal diversity metrics

Fragment Based Drug Discovery (FBDD) continues to advance as an efficient and alternative screening paradigm for the identification and optimization of novel chemical matter. To enable FBDD across a wide range of pharmaceutical targets, a fragment screening library is required to be chemically diverse and synthetically expandable to enable critical decision making for chemical follow-up and assessing new target druggability. In this manuscript, the Pfizer fragment library design strategy which utilized multiple and orthogonal metrics to incorporate structure, pharmacophore and pharmacological space diversity is described. Appropriate measures of molecular complexity were also employed to maximize the probability of detection of fragment hits using a variety of biophysical and biochemical screening methods. In addition, structural integrity, purity, solubility, fragment and analog availability as well as cost were important considerations in the selection process. Preliminary analysis of primary screening results for 13 targets using NMR Saturation Transfer Difference (STD) indicates the identification of uM-mM hits and the uniqueness of hits at weak binding affinities for these targets.

Synthesis of macrocyclic, potential protease inhibitors using a generic scaffold

A generic macrocyclic peptide structure 2 was designed as a potential inhibitor of a range of proteinases, by using as a basis for the design the known structures of a series of enzyme-inhibitor complexes. The macrocyclic nature of the target 2 was chosen so as to reduce the entropic advantage in the hydrolytic enzymatic step, and thereby to inhibit the function of the enzyme. The nature of the linking group was identified as a benzoxazole by molecular modeling, so as to preserve the recognized conformation of the peptide chain. The specificity of the potential inhibitor was tuned by variation of the P(1) group (by incorporating phenylalanine, aspartic acid, or lysine), to allow recognition by different enzyme classes. The targets were prepared from the bis-amino acid derivative 5, itself prepared using the Pd-catalyzed coupling of an organozinc reagent with the iodobenzothiazole 7 and subsequent macrocyclization of the open-chain derivatives 22-24 using HATU. None of the macrocylic compounds 25, 28-30, and 32 inhibited their target enzymes. NMR and MS studies on the interaction of macrocycle 29 and chymotrypsin established that compound 29 was in fact a substrate of the enzyme. This result indicated that while the design had been partially successful in identifying a compound that bound, the reduction in entropic advantage due to its macrocyclic nature was not sufficient to allow 29 to act as an inhibitor.

Cyclization of N-terminal S-carbamoylmethylcysteine causing loss of 17 Da from peptides and extra peaks in peptide maps

Enzymatic digests of proteins S-alkylated with iodoacetamide may contain peptides with N-terminal S-carbamoylmethylcysteine. These can be partly converted to a form with 17 Da lower mass and increased HPLC retention. Proof by synthesis supported by MS/MS and NMR spectroscopy was used to show that N-terminal S-carbamoylmethyl-L-cysteine can cyclize, losing NH3 to form an N-terminal residue of (R)-5-oxoperhydro-1,4-thiazine-3-carboxylic acid. The abbreviation Otc is proposed for the (R)-5-oxoperhydro-1,4-thiazine-3-carbonyl residue. The rate of cyclization is significant in 0.1 M NH4HCO3 at 37 degrees C, with the half-life of the acyclic form being 10-12 h for several peptides tested. This is similar to the rate at which N-terminal pyroglutamate forms from N-terminal glutamine.

Allosteric inhibition of fructose-1,6-bisphosphatase by anilinoquinazolines

Anilinoquinazolines currently of interest as inhibitors of tyrosine kinases have been found to be allosteric inhibitors of the enzyme fructose 1,6-bisphosphatase. These represent a new approach to inhibition of F16BPase and serve as leads for further drug design. Enzyme inhibition is achieved by binding at an unidentified allosteric site.

Topology of the pore-region of a K+ channel revealed by the NMR-derived structures of scorpion toxins

The architecture of the pore-region of a voltage-gated K+ channel, Kv1.3, was probed using four high affinity scorpion toxins as molecular calipers. We established the structural relatedness of these toxins by solving the structures of kaliotoxin and margatoxin and comparing them with the published structure of charybdotoxin; a homology model of noxiustoxin was then developed. Complementary mutagenesis of Kv1.3 and these toxins, combined with electrostatic compliance and thermodynamic mutant cycle analyses, allowed us to identify multiple toxin-channel interactions. Our analyses reveal the existence of a shallow vestibule at the external entrance to the pore. This vestibule is approximately 28-32 A wide at its outer margin, approximately 28-34 A wide at its base, and approximately 4-8 A deep. The pore is 9-14 A wide at its external entrance and tapers to a width of 4-5 A at a depth of approximately 5-7 A from the vestibule. This structural information should directly aid in developing topological models of the pores of related ion channels and facilitate therapeutic drug design.

A thermodynamic scale for the beta-sheet forming tendencies of the amino acids

The results of a study to measure the beta-sheet forming propensities of the 20 naturally occurring amino acids are presented. The protein host for these studies is the 56 amino acid B1 domain of staphylococcal IgG binding protein G [Fahnestock, S.R., Alexander, P., Nagle, J., & Filpula, D. (1986) J. Bacteriol. 167, 870-880]. This protein was selected because it exhibits a reversible two-state thermal denaturation transition and its structure is known at high resolution. A suitable guest position in the protein was identified, and its neighboring environment was modified to minimize the potential for artifactual interactions. All 20 amino acids were individually substituted at the guest site, and their effect on the protein's thermal stability was determined. NMR was used to verify the structural integrity of several of the proteins with different amino acid substitutions at the guest site. The results of these studies provide a thermodynamic scale for the relative beta-sheet forming propensities of the amino acids that shows a clear correlation with the beta-sheet preferences derived from statistical surveys of proteins of known structure.

The CD58 (LFA-3) binding site is a localized and highly charged surface area on the AGFCC'C" face of the human CD2 adhesion domain

Using site-directed mutagenesis in conjunction with NMR structural data on the adhesion domain of human CD2, we have defined the binding region for CD58. Previous structural studies of rat and human CD2 indicate that this adhesion domain is immunoglobulin-like. Here we report that the CD58 binding site is a well-circumscribed, charged surface area covering approximately 770 A2 on the AGFCC'C" face of the CD2 beta barrel. This site contains beta-strand residues in the carboxyl-terminal half of the F strand (including Lys-82 and Tyr-86), the top of the C strand (Asp-32 and Lys-34), and the C' strand (Gln-46), which are all solvent exposed. In addition, several exposed residues on the FG loop (Gly-90, Lys-91, Asn-92, and Val-93), the CC' loop (Lys-41 and Lys-43), and the C'C" loop (Arg-48 and Lys-51) form this site. In contrast, neither residues on the more peripheral G and C" strands of the same CD2 surface nor residues on B, E, and D strands of the opposite face are involved in CD58 binding. This CD58 binding site is predicted to lie most distal to the T-lymphocyte surface membrane, with ready access to CD58 on the surface of the opposing antigen-presenting cell.

1H resonance assignments and secondary structure of the 13.6 kDa glycosylated adhesion domain of human CD2

Human CD2, a glycosylated transmembrane receptor found on all T-lymphocytes, plays a key role in facilitating cellular adhesion between T-cells and target cells or antigen-presenting cells by binding to its counter receptor CD58 (LFA-3) present on the surface of those cells. All CD2 adhesion functions are localized within the amino-terminal 105-residue domain, which contains a single high mannose N-glycan required for maintaining both the conformational stability and CD58 binding properties of the glycoprotein. In order to better understand the structural basis for CD2-CD58-mediated adhesion and the critical role of the carbohydrate moiety in maintaining the functional stability of the molecule, we have determined the secondary structure of the N-glycosylated adhesion domain of human CD2 (hu-sCD2(105)) using NMR spectroscopy. Most of the 1H resonance assignments have been obtained from 1H-1H homonuclear 2D NMR spectra, which were further extended by applying 1H-15N heteronuclear 2D experiments on a hu-sCD2(105) sample selectively labeled with [15N]lysine. Thus, 98% of all backbone 1H resonances and over 80% of all side chain 1H resonances have been assigned. An overall topology characteristic of an immunoglobulin variable domain is observed, which consists of two beta-sheets comprised of three (residues 16-20, 67-71, and 60-63) and five (residues 94-103, 80-86, 32-37, 45-47, and 53-55) antiparallel beta-strands, respectively, with a hydrophobic core sandwiched between them. A ninth beta-strand (residues 7-12) makes parallel contacts to the carboxy-terminal beta-strand. NOEs between the N-linked glycan and the protein have tentatively been identified.

Structure of the glycosylated adhesion domain of human T lymphocyte glycoprotein CD2

BACKGROUND

CD2, a T-cell specific surface glycoprotein, is critically important for mediating adherence of T cells to antigen-presenting cells or target cells. Domain 1 of human CD2 is responsible for cell adhesion, binding to CD58 (LFA-3) expressed on the cell to which the T cell binds. Human CD2 domain 1 requires N-linked carbohydrate to maintain its native conformation and ability to bind CD58. In contrast, rat CD2 does not require N-linked carbohydrate, and binds to a different ligand, CD48.

RESULTS

The three-dimensional structure of the glycosylated form of domain 1 of human CD2 has been determined by NMR spectroscopy. The overall structure resembles the typical beta-barrel of an immunoglobulin variable domain. Nuclear Overhauser enhancement contacts between the protein and the N-linked glycan have been tentatively identified.

CONCLUSION

Based on our results, we propose a model showing how the N-linked glycan might be positioned in the human CD2 domain 1 structure. The model provides an explanation for the observed instability of deglycosylated human CD2, and allows residues that are important for CD58 binding to be differentiated from those affecting conformational stability via interactions with the glycan.

Effects of the presence of an aldehydic abasic site on the thermal stability and rates of helix opening and closing of duplex DNA

The presence of an abasic site in duplex DNA lowers the thermodynamic stability, as monitored by the optical melting temperature, and decreases the rate of imino proton exchange with water, by about an order of magnitude, as monitored by direct measurement of both the exchange lifetimes and the imino proton T1S. The exchange lifetimes of the imino protons with water as a function of base catalyst concentration were analyzed to determine the origin of the effect of the abasic site on imino exchange lifetimes. Analysis of the results showed that the helix opening rate is not significantly changed by the presence of an abasic site. The differences in exchange lifetimes are attributed to a faster helix closing rate in the presence of an abasic site. The faster rate of helix closing may be an important contribution to the stability of abasic sites in duplex DNA to base-catalyzed elimination reaction. It is noted that duplex DNAs containing analogues of the aldehydic abasic site apparently do not exhibit these exchange lifetime effects.

Toward a dynamical structure of DNA: comparison of theoretical and experimental NOE intensities

Comparisons of experimental and calculated interproton nuclear Overhauser effect (NOE) buildup curves for duplex d(CGCGAATTCGCG)2 have been made. The calculated NOEs are based on molecular dynamics simulations including counterions and water and on the single-structure canonical A, B, and crystal forms. The calculated NOE effects include consideration of the motions of individual interproton vectors and the anisotropic tumbling of the DNA. The effects due to inclusion of anisotropic tumbling are much larger than those due to the local motion, and both improve the agreement between calculated and experimental results. The predictions based on the dynamical models agree significantly better with experiment than those based on either of the canonical forms or the crystal structure.

Characterization of conformational features of DNA heteroduplexes containing aldehydic abasic sites

The DNA duplexes shown below, with D indicating deoxyribose aldehyde absic sites and numbering from 5' to 3', have been investigated by NMR. The 31P and 31P-1H correlation data indicate [formula: see text] that the backbones of these duplex DNAs are regular. One- and two-dimensional 1H NMR data indicate that the duplexes are right-handed and B-form. Conformational changes due to the presence of the abasic site extend to the two base pairs adjacent to the lesion site with the local conformation of the DNA being dependent on whether the abasic site is in the alpha or beta configuration. The aromatic base of residue A17 in the position opposite the abasic site is predominantly stacked in the helix as is G17 in the analogous sample. Imino lifetimes of the AT base pairs are much longer in samples with an abasic site than in those containing a Watson-Crick base pair. The conformational and dynamical properties of the duplex DNAs containing the naturally occurring aldehyde abasic site are different from those of duplex DNAs containing a variety of analogues of the abasic site.

Molecular dynamics of an in vacuo model of duplex d(CGCGAATTCGCG) in the B-form based on the amber 3.0 force field

The characteristics of 100 ps of molecular dynamics (MD) on the DNA dodecamer d(CGCGAATTCGCG) at 300 K are described and investigated. The simulation is based on an in vacuo model of the oligomer and the AMBER 3.0 force field configured in the manner of Singh, U. C., S. J. Weiner, and P. A. Kollman, (1985, Proc. Natl. Acad. Sci. USA. 82:755-759). The analysis of the results was carried out using the "curves, dials, and windows" procedure (Ravishanker, G., S. Swaminathan, D. L. Beveridge, R. Lavery, and H. Sklenar. 1989. J. Biomol. Struct. Dyn. 6:669-699). The results indicate this dynamical model to be a provisionally stable double helix which lies at approximately 3.2 A rms deviation from the canonical B-form. There is, however, a persistent nonplanarity in the base pair orientations which resemble that observed in canonical A-DNA. The major groove width is seen to narrow during the course of the simulation and the minor groove expands, contravariant to the alterations in groove width seen in the crystal structure of the native dodecamer (Drew, H. R., R. M. Wing, T. Takano, C. Broka, S. Tanaka, I. Itakura, and R. E. Dickerson, 1981. Proc. Natl. Acad. Sci. USA. 78:2179-2183). The propeller twist in the bases, the sequence dependence of the base pair roll and aspects of bending in the helix axis are in some degree of agreement with the crystal structure. The patterns in DNA bending are observed to follow Zhurkin theory (Zhurkin, V. B. 1985. J. Biomol. Struct. Dyn. 2:785-804.). The relationship between the dynamical model and structure in solution is discussed.