Molecular determinants of ASIC1 modulation by divalent cations

Acid-sensing ion channels (ASICs) are proton-gated cation channels widely expressed in the nervous system. ASIC gating is modulated by divalent cations as well as small molecules; however, the molecular determinants of gating modulation by divalent cations are not well understood. Previously, we identified two small molecules that bind to ASIC1a at a novel site in the acidic pocket and modulate ASIC1 gating in a manner broadly resembling divalent cations, raising the possibility that these small molecules may help to illuminate the molecular determinants of gating modulation by divalent cations. Here, we examined how these two groups of modulators might interact as well as mutational effects on ASIC1a gating and its modulation by divalent cations. Our results indicate that binding of divalent cations to an acidic pocket site plays a key role in gating modulation of the channel.

Discovery of Potent and Orally Bioavailable Pyridine N-Oxide-Based Factor XIa Inhibitors through Exploiting Nonclassical Interactions

Activated factor XI (FXIa) inhibitors are promising novel anticoagulants with low bleeding risk compared with current anticoagulants. The discovery of potent FXIa inhibitors with good oral bioavailability has been challenging. Herein, we describe our discovery effort, utilizing nonclassical interactions to improve potency, cellular permeability, and oral bioavailability by enhancing the binding while reducing polar atoms. Beginning with literature-inspired pyridine N-oxide-based FXIa inhibitor 1, the imidazole linker was first replaced with a pyrazole moiety to establish a polar C-H···water hydrogen-bonding interaction. Then, structure-based drug design was employed to modify lead molecule 2d in the P1' and P2' regions, with substituents interacting with key residues through various nonclassical interactions. As a result, a potent FXIa inhibitor 3f (K_i = 0.17 nM) was discovered. This compound demonstrated oral bioavailability in preclinical species (rat 36.4%, dog 80.5%, and monkey 43.0%) and displayed a dose-dependent antithrombotic effect in a rabbit arteriovenous shunt model of thrombosis.

Novel Reagent Space: Identifying Unorderable but Readily Synthesizable Building Blocks

Drug discovery building blocks available commercially or within an internal inventory cover a diverse range of chemical space and yet describe only a tiny fraction of all chemically feasible reagents. Vendors will eagerly provide tools to search the former; there is no straightforward method of mining the latter. We describe a procedure and use case in assembling chemical structures not available for purchase but that could likely be synthesized in one robust chemical transformation starting from readily available building blocks. Accessing this vast virtual chemical space dramatically increases our curated collection of reagents available for medicinal chemistry exploration and novel hit generation, almost tripling the number of those with 10 or fewer atoms.

Hit-to-lead optimization and discovery of a potent, and orally bioavailable G protein coupled receptor kinase 2 (GRK2) inhibitor

G-protein coupled receptor kinase 2 (GRK2), which is upregulated in the failing heart, appears to play a critical role in heart failure (HF) progression in part because enhanced GRK2 activity promotes dysfunction of β-adrenergic signaling and myocyte death. An orally bioavailable GRK2 inhibitor could offer unique therapeutic outcomes that cannot be attained by current heart failure treatments that directly target GPCRs or angiotensin-converting enzyme. Herein, we describe the discovery of a potent, selective, and orally bioavailable GRK2 inhibitor, 8h, through high-throughput screening, hit-to-lead optimization, structure-based design, molecular modelling, synthesis, and biological evaluation. In the cellular target engagement assays, 8h enhances isoproterenol-mediated cyclic adenosine 3',5'-monophosphate (cAMP) production in HEK293 cells overexpressing GRK2. Compound 8h was further evaluated in a human stem cell-derived cardiomyocyte (HSC-CM) contractility assay and potentiated isoproterenol-induced beating rate in HSC-CMs.

Tandem Suzuki Coupling/Intramolecular Oxetane Ring Opening to Form Polycyclic Ring Systems

A tandem one-pot reaction featuring a cross-coupling followed by an intramolecular oxetane ring opening by mild nucleophiles is reported. The overall transformation comprises a carbon-carbon bond formation along with a carbon-heteroatom bond construction providing diverse multicyclic ring systems with a pendant hydroxymethyl handle for further elaboration. This approach constitutes a convergent method for rapid access to various scaffolds. Furthermore, a comparison of computed low-energy conformers is presented to rationalize instances in which cyclization was not observed.

Current and Future Roles of Artificial Intelligence in Medicinal Chemistry Synthesis

Artificial intelligence and machine learning have demonstrated their potential role in predictive chemistry and synthetic planning of small molecules; there are at least a few reports of companies employing in silico synthetic planning into their overall approach to accessing target molecules. A data-driven synthesis planning program is one component being developed and evaluated by the Machine Learning for Pharmaceutical Discovery and Synthesis (MLPDS) consortium, comprising MIT and 13 chemical and pharmaceutical company members. Together, we wrote this perspective to share how we think predictive models can be integrated into medicinal chemistry synthesis workflows, how they are currently used within MLPDS member companies, and the outlook for this field.

Breaking the Glass Ceiling in Simulation and Modeling: Women in Pharmaceutical Discovery

The topic of gender equality within the United States workforce is receiving a great deal of attention. The field of chemistry is no exception and is increasingly focused on taking steps to achieve gender diversity within the chemistry workforce. Over the past several years, many computational chemistry groups within large pharmaceutical companies have realized growth in the number of women, and here we discuss the key factors that we believe have played a role in attracting and retaining the authors of this review as computational chemists in pharma. Furthermore, we combine our professional experiences in the context of how computational methodology and technology have evolved over the past decades and how that evolution has facilitated the inclusion of more women into the field. Our hope is to be a part of a solution and provide insight that will allow the chemistry workforce to continue to make steps forward in attaining gender diversity in the workplace.

6-Benzhydryl-4-amino-quinolin-2-ones as Potent Cannabinoid Type 1 (CB1) Receptor Inverse Agonists and Chemical Modifications for Peripheral Selectivity

A novel series of 6-benzhydryl-4-amino-quinolin-2-ones was discovered as cannabinoid type 1 receptor (CB₁R) inverse agonists based on the high-throughput screening hit, compound 1a. Structure-activity relationships were studied to improve in vitro/in vivo pharmacology and restrict distribution to the peripheral circulation. We adopted several strategies such as increasing topological polar surface area, incorporating discrete polyethylene glycol side chains, and targeting P-glycoprotein (P-gp) to minimize access to the brain. Compound 6a is a P-gp substrate and a potent and highly selective CB₁R inverse agonist, demonstrating excellent in vivo metabolic stability and a low brain to plasma ratio. However, brain receptor occupancy studies showed that compound 6a may accumulate in brain with repeat dosing. This was evidenced by compound 6a inhibiting food intake and inducing weight loss in diet-induced obese mice. Thus, a strategy based on P-gp efflux may not be adequate for peripheral restriction of the disclosed quinolinone series.

Enhancement of kinase selectivity in a potent class of arylamide FMS inhibitors

Structure-activity relationship (SAR) studies on a highly potent series of arylamide FMS inhibitors were carried out with the aim of improving FMS kinase selectivity, particularly over KIT. Potent compound 17r (FMS IC50 0.7 nM, FMS cell IC50 6.1 nM) was discovered that had good PK properties and a greater than fivefold improvement in selectivity for FMS over KIT kinase in a cellular assay relative to the previously reported clinical candidate 4. This improved selectivity was manifested in vivo by no observed decrease in circulating reticulocytes, a measure of bone safety, at the highest studied dose. Compound 17r was highly active in a mouse pharmacodynamic model and demonstrated disease-modifying effects in a dose-dependent manner in a strep cell wall-induced arthritis model of rheumatoid arthritis in rats.

Insights from structural analysis of cFMS/inhibitor complexes: common interactions via three structurally dissimilar scaffolds

A small-molecule drug discovery effort can benefit from having several chemical series. Where multiple series are not available, it is often the goal of a project to find novel scaffolds. Structural studies of ligand/protein complexes provide important information on the interactions driving binding. By generalizing these, it is possible to find molecules lacking in similarity in their connectivity yet retaining the ability to interact with the same target protein. Our studies on inhibitors of the cFMS tyrosine kinase provide a dramatic example of three different chemical series that make the same key interactions with the target protein. Collectively, these structural data provide a striking example of the pharmacophore hypothesis at work. In addition, they should prompt one to employ a broad approach when attempting scaffold hopping or any search for a novel series. It is clear that molecules that bind with similar interactions to a target need not possess 2-dimensional molecular similarity.

Crystal structure of a soluble form of human monoglyceride lipase in complex with an inhibitor at 1.35 Å resolution

A high-resolution structure of a ligand-bound, soluble form of human monoglyceride lipase (MGL) is presented. The structure highlights a novel conformation of the regulatory lid-domain present in the lipase family as well as the binding mode of a pharmaceutically relevant reversible inhibitor. Analysis of the structure lacking the inhibitor indicates that the closed conformation can accommodate the native substrate 2-arachidonoyl glycerol. A model is proposed in which MGL undergoes conformational and electrostatic changes during the catalytic cycle ultimately resulting in its dissociation from the membrane upon completion of the cycle. In addition, the study outlines a successful approach to transform membrane associated proteins, which tend to aggregate upon purification, into a monomeric and soluble form.

Biphenylsulfonyl-thiophene-carboxamidine inhibitors of the complement component C1s

Complement activation has been implicated in disease states such as hereditary angioedema, ischemia-reperfusion injury, acute respiratory distress syndrome, and acute transplant rejection. Even though the complement cascade provides several protein targets for potential therapeutic intervention only two complement inhibitors have been approved so far for clinical use including anti-C5 antibodies for the treatment of paroxysmal nocturnal hemoglobinuria and purified C1-esterase inhibitor replacement therapy for the control of hereditary angioedema flares. In the present study, optimization of potency and physicochemical properties of a series of thiophene amidine-based C1s inhibitors with potential utility as intravenous agents for the inhibition of the classical pathway of complement is described.

Processing of small molecule databases for automated docking

Virtual screening involves the mining of small molecule databases from various sources. The small molecule databases used in virtual screening are typically processed, from simple 2D representations, to maximize their information content and to optimize them for input to the particular virtual screening technology being used. Processing interprets or adds molecular information related to connectivity, stereochemistry, protonation, tautomers and conformation. For virtual screening with an automated docking protocol, a technique that relies on specific intermolecular atom-atom contacts for ranking molecules, it is expected that the pre-processing protocol can affect the results of the docking experiment. The possible effects of processing on docking results have not been extensively studied, and this topic has only recently emerged as a significant aspect of the docking-based virtual screening process. One recent report highlights significant effects of different processing procedures on docking enrichment, while another outlines a general ligand preparation strategy. Here we survey and comment on recent practice in the field.

A novel series of arylsulfonylthiophene-2-carboxamidine inhibitors of the complement component C1s

Inhibiting the classical pathway of complement activation by attenuating the proteolytic activity of the serine protease C1s is a potential strategy for the therapeutic intervention in disease states such as hereditary angioedema, ischemia-reperfusion injury, and acute transplant rejection. A series of arylsulfonylthiophene-2-carboxamidine inhibitors of C1s were synthesized and evaluated for C1s inhibitory activity. The most potent compound had a Ki of 10nM and >1000-fold selectivity over uPA, tPA, FX(a), thrombin, and plasmin.

Docking: successes and challenges

The state of the art of various computational aspects of docking-based virtual screening of database of small molecules is presented. The review encompasses the different search algorithms and the scoring functions used in docking methods and their applications to protein and nucleic acid drug targets. Recent progress made in the development and application of methods to include target flexibility are summarized. The fundamental issues and challenges involved in comparing various docking methods are discussed. Limitations of current technologies as well as future prospects are presented.

Non-peptidic αvβ3 antagonists containing indol-1-yl propionic acids

We describe the synthesis and structure/activity relationship of RGD mimetics that are potent inhibitors of the integrin alpha(v)beta3. Indol-1-yl propionic acids containing a variety of basic moieties at the 5-position, as well as substitutions alpha and beta to the carboxy terminus were synthesized and evaluated. Novel compounds with improved potency have been identified.

Comparison of automated docking programs as virtual screening tools

The performance of several commercially available docking programs is compared in the context of virtual screening. Five different protein targets are used, each with several known ligands. The simulated screening deck comprised 1000 molecules from a cleansed version of the MDL drug data report and 49 known ligands. For many of the known ligands, crystal structures of the relevant protein-ligand complexes were available. We attempted to run experiments with each docking method that were as similar as possible. For a given docking method, hit rates were improved versus what would be expected for random selection for most protein targets. However, the ability to prioritize known ligands on the basis of docking poses that resemble known crystal structures is both method- and target-dependent.

Decrypting the biochemical function of an essential gene from Streptococcus pneumoniae using ThermoFluor technology

The protein product of an essential gene of unknown function from Streptococcus pneumoniae was expressed and purified for screening in the ThermoFluor affinity screening assay. This assay can detect ligand binding to proteins of unknown function. The recombinant protein was found to be in a dimeric, native-like folded state and to unfold cooperatively. ThermoFluor was used to screen the protein against a library of 3000 compounds that were specifically selected to provide information about possible biological functions. The results of this screen identified pyridoxal phosphate and pyridoxamine phosphate as equilibrium binding ligands (K(d) approximately 50 pM, K(d) approximately 2.5 microM, respectively), consistent with an enzymatic cofactor function. Several nucleotides and nucleotide sugars were also identified as ligands of this protein. Sequence comparison with two enzymes of known structure but relatively low overall sequence homology established that several key residues directly involved in pyridoxal phosphate binding were strictly conserved. Screening a collection of generic drugs and natural products identified the antifungal compound canescin A as an irreversible covalent modifier of the enzyme. Our investigation of this protein indicates that its probable biological role is that of a nucleoside diphospho-keto-sugar aminotransferase, although the preferred keto-sugar substrate remains unknown. These experiments demonstrate the utility of a generic affinity-based ligand binding technology in decrypting possible biological functions of a protein, an approach that is both independent of and complementary to existing genomic and proteomic technologies.

A novel series of potent and selective small molecule inhibitors of the complement component C1s

Activation of the classical pathway of complement has been implicated in disease states such as hereditary angioedema, ischemia-reperfusion injury and acute transplant rejection. The trypsin-like serine protease C1s represents a pivotal upstream point of control in the classical pathway of complement activation and is therefore likely to be a useful target in the therapeutic intervention of these disease states. A series of thiopheneamidine-based inhibitors of C1s has been optimized to give a 70 nM inhibitor that inhibits the classical pathway of complement activation in vitro.

Discovery and SAR of novel Naphthyridines as potent inhibitors of spleen tyrosine kinase (SYK)

The discovery of novel 5,7-disubstituted[1,6]naphthyridines as potent inhibitors of Spleen Tyrosine Kinase (SYK) is discussed. The SAR reveals the necessity for a 7-aryl group with preference towards para substitution and that this in combination with 5-aminoalkylamino substituents further improved the potency of the compounds. The initial SAR as well as a survey of the other positions is discussed in detail.

Synthesis of thiophene-2-carboxamidines containing 2-aminothiazoles and their biological evaluation as urokinase inhibitors

The serine protease urokinase (uPa) has been implicated in the progression of both breast and prostate cancer. Utilizing structure based design, the synthesis of a series of substituted 4-[2-amino-1,3-thiazolyl]-thiophene-2-carboxamidines is described. Further optimization of this series by substitution of the terminal amine yielded urokinase inhibitors with excellent activities.

ASP147 in the third transmembrane helix of the rat mu opioid receptor forms ion-pairing with morphine and naltrexone

We tested the hypotheses that the carboxylate side chain of Asp147 of the mu opioid receptor interacts with the protonated nitrogen of naltrexone and morphine and that this interaction is important for pharmacological properties of the two compounds. Mutation of Asp147 to Ala or Asn substantially reduced the affinity of naltrexone and the affinity, potency and efficacy of morphine, while the Glu mutant had similar properties as the wildtype, indicating the significant role of the carboxylate group of Asp147 in receptor binding and activation. This role could be due to its direct interaction with ligands or involvement in interhelical interactions. The unprotonated analogs of naltrexone and morphine, cyclopropylcarbonyl noroxymorphone (CPCNOM) and N-formylnormorphine (NFNM), respectively, were used to discriminate between these mechanisms. CPCNOM was much less potent as an antagonist and had substantially lower affinity for the mu receptor than naltrexone. Similarly, NFNM was unable to activate the mu receptor and had much lower affinity than morphine. These results indicate the importance of the protonated nitrogen. Notably, the D147A and D147N mutations did not appreciably affect the binding affinities of CPCNOM and NFNM. In addition, the D147E mutant had similar affinities for CPCNOM and NFNM as the D147A and D147N mu receptors. Thus, the carboxylate group of Asp147 is not important for binding of the two unprotonated compounds. These results indicate that the carboxylate group of Asp147 of the mu receptor interacts directly with the protonated nitrogen of naltrexone and morphine and this interaction is important for binding and receptor activation.

Design and synthesis of diaminopyrrolidinone inhibitors of human osteoclast cathepsin K

The structure-based design and synthesis of lactam-constrained azapeptide inhibitors of human cathepsin K are described. Enhanced stability to proteolytic cleavage over acyclic analogues is discussed.

Structure-based design of non-peptide, carbohydrazide-based cathepsin K inhibitors

Using binding models which were based on the X-ray crystal structure of an amino acid-based active site-spanning inhibitor complexed with cathepsin K, Cbz-leucine mimics have been developed, leading ultimately to the design of a potent cathepsin K inhibitor free of amino acid components. These mimics, which consist of alpha-substituted biphenylacetyl groups in place of Cbz-leucine moieties, effectively mimic all aspects of the Cbz-leucine moieties which are important for inhibitor binding. The predicted directions of binding for the inhibitors were confirmed by mass spectral analysis of their complexes with cathepsin K, which gave results consistent with acylation of the enzyme and loss of the acylhydrazine portion of the inhibitor which binds on the S' side of the active site. The binding models were found to be very predictive of relative inhibitor potency as well as direction of inhibitor binding. These results strengthen the validity of a strategy involving iterative cycles of structure-based design and inhibitor synthesis and evaluation for the discovery of non-peptide inhibitors.

Use of papain as a model for the structure-based design of cathepsin K inhibitors: crystal structures of two papain-inhibitor complexes demonstrate binding to S'-subsites

Papain has been used as a surrogate enzyme in a drug design effort to obtain potent and selective inhibitors of cathepsin K, a new member of the papain superfamily of cysteine proteases that is selectively and highly expressed in osteoclasts and is implicated in bone resorption. Here we report the crystal structures of two papain-inhibitor complexes and the rational design of novel cathepsin K inhibitors. Unlike previously known crystal structures of papain-inhibitor complexes, our papain structures show ligand binding extending deep within the S'-subsites. The two inhibitor complexes, carbobenzyloxyleucinyl-leucinyl-leucinal and carbobenzyloxy-L-leucinyl-L-leucinyl methoxymethyl ketone, were refined to 2.2- and 2.5-A resolution with R-factors of 0.190 and 0. 217, respectively. The S'-subsite interactions with the inhibitors are dominated by an aromatic-aromatic stacking and an oxygen-aromatic ring edge interaction. The knowledge of S'-subsite interactions led to a design strategy for an inhibitor spanning both subsites and yielded a novel, symmetric inhibitor selective for cathepsin K. Simultaneous exploitation of both S- and S'-sites provides a general strategy for the design of cysteine protease inhibitors having high specificity to their target enzymes.

Structure-based design of cathepsin K inhibitors containing a benzyloxy-substituted benzoyl peptidomimetic

Peptidomimetic cathepsin K inhibitors have been designed using binding models which were based on the X-ray crystal structure of an amino acid-based, active site-spanning inhibitor complexed with cathepsin K. These inhibitors, which contain a benzyloxybenzoyl group in place of a Cbz-leucine moiety, maintained good inhibitory potency relative to the amino acid-based inhibitor, and the binding models were found to be very predictive of relative inhibitor potency. The binding mode of one of the inhibitors was confirmed by X-ray crystallography, and the crystallographically determined structure is in close qualitative agreement with the initial binding model. These results strengthen the validity of a strategy involving iterative cycles of structure-based design, inhibitor synthesis and evaluation, and crystallographic structure determination for the discovery of peptidomimetic inhibitors.

Design of potent and selective human cathepsin K inhibitors that span the active site

Potent and selective active-site-spanning inhibitors have been designed for cathepsin K, a cysteine protease unique to osteoclasts. They act by mechanisms that involve tight binding intermediates, potentially on a hydrolytic pathway. X-ray crystallographic, MS, NMR spectroscopic, and kinetic studies of the mechanisms of inhibition indicate that different intermediates or transition states are being represented that are dependent on the conditions of measurement and the specific groups flanking the carbonyl in the inhibitor. The species observed crystallographically are most consistent with tetrahedral intermediates that may be close approximations of those that occur during substrate hydrolysis. Initial kinetic studies suggest the possibility of irreversible and reversible active-site modification. Representative inhibitors have demonstrated antiresorptive activity both in vitro and in vivo and therefore are promising leads for therapeutic agents for the treatment of osteoporosis. Expansion of these inhibitor concepts can be envisioned for the many other cysteine proteases implicated for therapeutic intervention.

Determination of the amino acid residue involved in [3H]beta-funaltrexamine covalent binding in the cloned rat mu-opioid receptor

We previously demonstrated that [3H]beta-funaltrexamine ([3H]beta-FNA) labeled the rat mu opioid receptor expressed in Chinese hamster ovary cells with high specificity, and [3H]beta-FNA-labeled receptors migrated as one broad band with a mass of 80 kDa. In this study, we determined the region and then the amino acid residue of the mu receptor involved in the covalent binding of [3H]beta-FNA. [3H]beta-FNA-labeled receptors were solubilized and purified to approximately 10% purity by immunoaffinity chromatography with antibodies against a C-terminal domain peptide. The site of covalent bond formation was determined to be within Ala206-Met243 by CNBr cleavage of partially purified labeled mu receptors and determinations of sizes of labeled receptor fragments. The amino acid residue of beta-FNA covalent incorporation was then determined by site-directed mutagenesis studies within this region. Mutation of Lys233 to Ala, Arg, His, and Leu completely eliminated covalent binding of [3H]beta-FNA, although these mutants bound beta-FNA with high affinity. Mutations of other amino acid residues did not affect covalent binding of [3H]beta-FNA. These results indicate that [3H]beta-FNA binds covalently to Lys233. Since [3H]beta-FNA is a rigid molecule, the information will be very useful for molecular modeling of interaction between morphinans and the mu receptor.

A paradigm for drug discovery using a conformation from the crystal structure of a presentation scaffold

We describe a structural validation of the use of presentation scaffolds for control and elucidation of bioactive conformations of peptides. The protein REI-RGD34--produced by inserting the sequence RIPRGDMP into the CDR1 loop region of the immunoglobulin VL domain REI--strongly inhibits fibrinogen binding to the integrins alpha IIb beta 3 and alpha V beta 3. In the X-ray crystal structure of their protein at 2.4 A resolution, the RGD-containing loop exhibits defined electron density that is consistent with models for the bioactive conformations of ligands of these receptors based on previous small-molecule studies. Furthermore, a search of a small-molecule database with conformational information derived from the structure of REI-RGD34 identified constrained peptides and peptidomimetics known to be antagonists of the platelet receptor alpha IIb beta 3.

A check on rational drug design: crystal structure of a complex of human immunodeficiency virus type 1 protease with a novel gamma-turn mimetic inhibitor

We have previously reported (Newlander et al., J. Med. Chem. 1993, 36, 2321-2331) the design of human immunodeficiency virus type 1 (HIV-1) protease inhibitors incorporating C7 mimetics that lock three amino acid residues of a peptide sequence into a gamma-turn. The design of one such compound, SB203238, was based on X-ray structures of reduced amide aspartyl protease inhibitors. It incorporates a gamma-turn mimetic in the P2-P1' position, where the carbonyl of the C7 ring is replaced with an sp3 methylene group yielding a constrained reduced amide. It shows competitive inhibition with Ki = 430 nM at pH 6.0. The three-dimensional structure of SB203238 bound to the active site of HIV-1 protease has been determined at 2.3 A resolution by X-ray diffraction and refined to a crystallographic R-factor (R = sigma magnitude of Fo magnitude of - magnitude of Fc magnitude of /sigma magnitude of Fo magnitude of, where Fo and Fc are the observed and calculated structure factor amplitudes, respectively) of 0.177. The inhibitor lies in an extended conformation in the active site; however, because of the constrained geometry of the C7 ring, it maintains fewer hydrogen bonds with the protein than in most other HIV-1 protease-inhibitor complexes. More importantly, the inhibitor binds to the enzyme differently than predicted in its design, by binding with the P2-P1' alpha-carbon atoms shifted by approximately one-half a residue toward the N-terminus from their presumed positions. This study illustrates the importance of structural information in an approach to rational drug design.

An orally bioavailable HIV-1 protease inhibitor containing an imidazole-derived peptide bond replacement: crystallographic and pharmacokinetic analysis

(2R,4S,5S,1'S)-2-Phenylmethyl-4-hydroxy-5-(tert-butoxycarbonyl) amino-6-phenylhexanoyl-N-(1'-imidazo-2-yl)-2'-methylpropanamide (compound 2) is a tripeptide analogue inhibitor of HIV-1 protease in which a C-terminal imidazole substituent constitutes an isoelectronic, structural mimic of a carboxamide group. Compound 2 is a potent inhibitor of the protease (K(i) = 18 nM) and inhibits HIV-1 acute infectivity of CD4+ T-lymphocytes (IC50 = 570 nM). Crystallographic analysis of an HIV-1 protease-compound 2 complex demonstrates that the nitrogen atoms of the imidazole ring assume the same hydrogen-bonding interactions with the protease as amide linkages in other peptide analogue inhibitors. The sole substitution of the C-terminal carboxamide of a hydroxyethylene-containing tripeptide analogue with an imidazole group imparts greatly improved pharmacokinetic and oral bioavailability properties on the compound compared to its carboxamide-containing homologue (compound 1). While the oral bioavailability of compound 1 in rats was negligible, compound 2 displayed oral bioavailabilities of 30% and 14%, respectively, in rats and monkeys.

Rational design, synthesis, and crystallographic analysis of a hydroxyethylene-based HIV-1 protease inhibitor containing a heterocyclic P1'--P2' amide bond isostere

The rational design and synthesis of a highly potent inhibitor of HIV-1 protease have been accomplished. The inhibitor, SB 206343, is based on a model derived from the structure of the MVT-101/HIV-1 protease complex and contains a 4(5)-acylimidazole ring as an isosteric replacement for the P1'--P2' amide bond. It is a competitive inhibitor with an apparent inhibition constant of 0.6 nM at pH 6.0. The three-dimensional structure of SB 206343 bound in the active site of HIV-1 protease has been determined at 2.3 A resolution by X-ray diffraction techniques and refined to a crystallographic discrepancy factor, R (= sigma parallel Fo magnitude of/Fc parallel/sigma magnitude of), of 0.194. The inhibitor is held in the enzyme by a set of hydrophobic and polar interactions. N-3 of the imidazole ring participates in a novel hydrogen-bonding interaction with the bound water molecule, demonstrating the effectiveness of the imidazole ring as an isosteric replacement for the P1'--P2' amide bond in hydroxyethylene-based HIV-1 protease inhibitors. Also present are hydrogen-bonding interactions between N-1 of the imidazole ring and the carbonyl of Gly-127 as well as between the imidazole acyl carbonyl oxygen and the amide nitrogen of Asp-129, exemplifying the peptidomimetic nature of the 4(5)-acylimidazole isostere. All of these interactions are in qualitative agreement with those predicted by the model.

A shape- and chemistry-based docking method and its use in the design of HIV-1 protease inhibitors

The program DOCK [1,2] has been used successfully to identify molecules which will bind to a specified receptor [3]. The original method ranks molecules based on their shape complementarity to the receptor site and relies on the chemist to bring the appropriate electrostatic or hydrogen bond properties into the molecular skeletons obtained in the search. This is useful when screening a small database of compounds, where it is not likely that molecules with both the correct shape and electrostatic properties will be found. As large databases are more likely to have redundant molecular shapes with a variety of functionality (e.g., members of a congeneric series), it would be useful to have a method which identifies molecules with both the correct shape and functionality. To this end we have modified the DOCK 1.0 method to target user-specified atom types to selected positions in the receptor site. The target sites can be chosen based on structural evidence, calculation or inspection. Targeted-DOCK improves the ability of the DOCK method to find the crystallographically determined binding mode of a ligand. Additionally, targeted-DOCK searches a database of small molecules at 100-1000 times the rate of DOCK 1.0, allowing more molecules to be screened and more sophisticated scoring schemes to be employed. Targeted-DOCK has been used successfully in the design of a novel non-peptide inhibitor of HIV-1 protease.

Strong inhibition of fibrinogen binding to platelet receptor alpha IIb beta 3 by RGD sequences installed into a presentation scaffold

In order to probe the structural constraints on binding of RGD sequences to the platelet receptor alpha IIb beta 3 we have used recombinant DNA techniques to install the RGD sequence into 'presentation scaffolds', small proteins of known 3-D structure chosen to present guest sequences in constrained orientations. Using Escherichia coli expression systems we made sequence variants in which loop residues of the immunoglobulin VL domain REI and of human interleukin-1 beta were replaced (without changing polypeptide length) by the RGD sequence at positions predicted, based on small molecule studies, to orient the RGD moiety into an active conformation. These variants do not compete for fibrinogen binding to alpha IIb beta 3 up to almost 1 mM concentration. Unfolded or proteolytically fragmented forms of these same proteins do compete, however, showing that the RGD sequences in the mutants must be prohibited from binding by constraints imposed by scaffold structure. To suppress the effects of such structural constraints we constructed two sequence variants in which RGD-containing sequences 42-57 or 44-55 from the snake venom platelet antagonist kistrin were inserted (this increasing the length of the loop) into the third complementarity determining loop of REI. Both of these variants compete strongly for fibrinogen binding with IC50s in the nM range. These results, plus data on kistrin-related peptides also presented here, suggest that the molecular scaffold REI is capable of providing to an installed sequence a structural context and conformation beneficial to binding. The results also suggest that in order to bind well to alpha IIb beta 3, RGD sequences in protein ligands must either project significantly from the surface of the scaffold and/or retain a degree of conformational flexibility within the scaffold. Molecular scaffolds like REI should prove useful in the elucidation of structure-function relationships and the discovery of new active sequences, and may also serve as the basis for novel therapeutic agents.

Inhibition of human immunodeficiency virus-1 protease by a C2-symmetric phosphinate. Synthesis and crystallographic analysis

The human immunodeficiency virus type 1 (HIV-1) protease is a potential target of acquired immune deficiency syndrome (AIDS) therapy. A highly potent, perfectly symmetrical phosphinate inhibitor of this enzyme, SB204144, has been synthesized. It is a competitive inhibitor of HIV-1 protease, with an apparent inhibition constant of 2.8 nM at pH 6.0. The three-dimensional structure of SB204144 bound to the enzyme has been determined at 2.3-A resolution by X-ray diffraction techniques and refined to a crystallographic discrepancy factor, R (= sigma parallel F(o) magnitude to - Fc parallel/sigma magnitude of F(o)), of 0.178. The inhibitor is held in the enzyme active site by a set of hydrophobic and hydrophilic interactions, including an interaction between Arg8 and the center of the terminal benzene rings of the inhibitor. The phosphinate establishes a novel interaction with the two catalytic aspartates; each oxygen of the central phosphinic acid moiety interacts with a single oxygen of one aspartic acid, establishing a very short (2.2-2.4 A) oxygen-oxygen contact. As with the structures of penicillopepsin bound to phosphinate and phosphonate inhibitors [Fraser, M. E., Strynadka, N. C., Bartlett, P. A., Hanson, J. E., & James, M. N. (1992) Biochemistry 31, 5201-14], we interpret this short distance and the stereochemical environment of each pair of oxygens in terms of a hydrogen bond that has a symmetric single-well potential energy curve with the proton located midway between the two atoms.(ABSTRACT TRUNCATED AT 250 WORDS)

A symmetric inhibitor binds HIV-1 protease asymmetrically

Potential advantages of C2-symmetric inhibitors designed for the symmetric HIV-1 protease include high selectivity, potency, stability, and bioavailability. Pseudo-C2-symmetric monools and C2-symmetric diols, containing central hydroxymethylene and (R,R)-dihydroxyethylene moieties flanked by a variety of hydrophobic P1/P1' side chains, were studied as HIV-1 protease inhibitors. The monools and diols were synthesized in 8-10 steps from D-(+)-arabitol and D-(+)-mannitol, respectively. Monools with ethyl or isobutyl P1/P1' side chains were weak inhibitors of recombinant HIV-1 protease (Ki > 10 microM), while benzyl P1/P1' side chains afforded a moderately potent inhibitor (apparent Ki = 230 nM). Diols were 100-10,000x more potent than analogous monools, and a wider range of P1/P1' side chains led to potent inhibition. Both classes of compounds exhibited lower apparent Ki values under high-salt conditions. Surprisingly, monool and diol HIV-1 protease inhibitors were potent inhibitors of porcine pepsin, a prototypical asymmetric monomeric aspartic protease. These results were evaluated in the context of the pseudosymmetric structure of monomeric aspartic proteases and their evolutionary kinship with the retroviral proteases. The X-ray crystal structure of HIV-1 protease complexed with a symmetric diol was determined at 2.6 A. Contrary to expectations, the diol binds the protease asymmetrically and exhibits 2-fold disorder in the electron density map. Molecular dynamics simulations were conducted beginning with asymmetric and symmetric HIV-1 protease/inhibitor model complexes. A more stable trajectory resulted from the asymmetric complex, in agreement with the observed asymmetric binding mode. A simple four-point model was used to argue more generally that van der Waals and electrostatic force fields can commonly lead to an asymmetric association between symmetric molecules.

Structure-based design of nonpeptide inhibitors specific for the human immunodeficiency virus 1 protease

By using a structure-based computer-assisted search, we have found a butyrophenone derivative that is a selective inhibitor of the human immunodeficiency virus 1 (HIV-1) protease. The computer program creates a negative image of the active site cavity using the crystal structure of the HIV-1 protease. This image was compared for steric complementarity with 10,000 molecules of the Cambridge Crystallographic Database. One of the most interesting candidates identified was bromperidol. Haloperidol, a closely related compound and known antipsychotic agent, was chosen for testing. Haloperidol inhibits the HIV-1 and HIV-2 proteases in a concentration-dependent fashion with a Ki of approximately 100 microM. It is highly selective, having little inhibitory effect on pepsin activity and no effect on renin at concentrations as high as 5 mM. The hydroxy derivative of haloperidol has a similar effect on HIV-1 protease but a lower potency against the HIV-2 enzyme. Both haloperidol and its hydroxy derivative showed activity against maturation of viral polypeptides in a cell assay system. Although this discovery holds promise for the generation of nonpeptide protease inhibitors, we caution that the serum concentrations of haloperidol in normal use as an antipsychotic agent are less than 10 ng/ml (0.03 microM). Thus, concentrations required to inhibit the HIV-1 protease are greater than 1000 times higher than the concentrations normally used. Haloperidol is highly toxic at elevated doses and can be life-threatening. Haloperidol is not useful as a treatment for AIDS but may be a useful lead compound for the development of an antiviral pharmaceutical.

Computerized selection of potential DNA binding compounds

Using a general shape-search docking algorithm, potential DNA minor groove binding compounds were selected from a subset from the Cambridge Crystallographic Database consisting of almost 10,000 molecules. The crystal structure of the DNA dodecamer as observed in the d(CGCGAATTCGCG)2.netropsin complex served as the target receptor. Surprisingly, the highest scoring compound turned out to be CC-1065, a potent anti-tumour agent. Netropsin itself was number 6 on the list of highest scoring compounds. A number of the top-10 scoring compounds may serve as a source of inspiration for further drug design.