Kathleen M. Frey

Predicting resistance mutations using protein design algorithms.

Drug resistance resulting from mutations to the target is an unfortunate common phenomenon that limits the lifetime of many of the most successful drugs. In contrast to the investigation of mutations after clinical exposure, it would be powerful to be able to incorporate strategies early in the development process to predict and overcome the effects of possible resistance mutations. Here we present a unique prospective application of an ensemble-based protein design algorithm, K∗, to predict potential resistance mutations in dihydrofolate reductase from Staphylococcus aureus using positive design to maintain catalytic function and negative design to interfere with binding of a lead inhibitor. Enzyme inhibition assays show that three of the four highly-ranked predicted mutants are active yet display lower affinity (18-, 9-, and 13-fold) for the inhibitor. A crystal structure of the top-ranked mutant enzyme validates the predicted conformations of the mutated residues and the structural basis of the loss of potency. The use of protein design algorithms to predict resistance mutations could be incorporated in a lead design strategy against any target that is susceptible to mutational resistance.

Crystal structures of wild-type and mutant methicillin-resistant Staphylococcus aureus dihydrofolate reductase reveal an alternate conformation of NADPH that may be linked to trimethoprim resistance.

Both hospital- and community-acquired Staphylococcus aureus infections have become major health concerns in terms of morbidity, suffering and cost. Trimethoprim-sulfamethoxazole (TMP-SMZ) is an alternative treatment for methicillin-resistant S. aureus (MRSA) infections. However, TMP-resistant strains have arisen with point mutations in dihydrofolate reductase (DHFR), the target for TMP. A single point mutation, F98Y, has been shown biochemically to confer the majority of this resistance to TMP. Using a structure-based approach, we have designed a series of novel propargyl-linked DHFR inhibitors that are active against several trimethoprim-resistant enzymes. We screened this series against wild-type and mutant (F98Y) S. aureus DHFR and found that several are active against both enzymes and specifically that the meta-biphenyl class of these inhibitors is the most potent. In order to understand the structural basis of this potency, we determined eight high-resolution crystal structures: four each of the wild-type and mutant DHFR enzymes bound to various propargyl-linked DHFR inhibitors. In addition to explaining the structure-activity relationships, several of the structures reveal a novel conformation for the cofactor, NADPH. In this new conformation that is predominantly associated with the mutant enzyme, the nicotinamide ring is displaced from its conserved location and three water molecules complete a network of hydrogen bonds between the nicotinamide ring and the protein. In this new position, NADPH has reduced interactions with the inhibitor. An equilibrium between the two conformations of NADPH, implied by their occupancies in the eight crystal structures, is influenced both by the ligand and the F98Y mutation. The mutation induced equilibrium between two NADPH-binding conformations may contribute to decrease TMP binding and thus may be responsible for TMP resistance.

Towards the Understanding of Resistance Mechanisms in Clinically Isolated Trimethoprim-resistant, Methicillin-resistant Staphylococcus aureus Dihydrofolate Reductase

Resistance to therapeutics such as trimethoprim-sulfamethoxazole has become an increasing problem in strains of methicillin-resistant Staphylococcus aureus (MRSA). Clinically isolated trimethoprim-resistant strains reveal a double mutation, H30N/F98Y, in dihydrofolate reductase (DHFR). In order to develop novel and effective therapeutics against these resistant strains, we evaluated a series of propargyl-linked antifolate lead compounds for inhibition of the mutant enzyme. For the propargyl-linked antifolates, the F98Y mutation generates minimal (between 1.2- and 6-fold) losses of affinity and the H30N mutation generates greater losses (between 2.4- and 48-fold). Conversely, trimethoprim affinity is largely diminished by the F98Y mutation (36-fold) and is not affected by the H30N mutation. In order to elucidate a mechanism of resistance, we determined a crystal structure of a complex of this double mutant with a lead propargyl-linked antifolate. This structure suggests a resistance mechanism consistent both for the propargyl-linked class of antifolates and for trimethoprim that is based on the loss of a conserved water-mediated hydrogen bond.

Synthetic and Crystallographic Studies of a New Inhibitor Series Targeting Bacillus anthracis Dihydrofolate Reductase

Bacillus anthracis, the causative agent of anthrax, poses a significant biodefense danger. Serious limitations in approved therapeutics and the generation of resistance have produced a compelling need for new therapeutic agents against this organism. Bacillus anthracis is known to be insensitive to the clinically used antifolate, trimethoprim, because of a lack of potency against the dihydrofolate reductase enzyme. Herein, we describe a novel lead series of B. anthracis dihydrofolate reductase inhibitors characterized by an extended trimethoprim-like scaffold. The best lead compound adds only 22 Da to the molecular weight and is 82-fold more potent than trimethoprim. An X-ray crystal structure of this lead compound bound to B. anthracis dihydrofolate reductase in the presence of NADPH was determined to 2.25 A resolution. The structure reveals several features that can be exploited for further development of this lead series.

Structure-Based Approach to the Development of Potent and Selective Inhibitors of Dihydrofolate Reductase from Cryptosporidium

Cryptosporidiosis is an emerging infectious disease that can be life-threatening in an immune-compromised individual and causes gastrointestinal distress lasting up to 2 weeks in an immune-competent individual. There are few therapeutics available for effectively treating this disease. We have been exploring dihydrofolate reductase (DHFR) as a potential target in Cryptosporidium. On the basis of the structure of the DHFR enzyme from C. hominis, we have developed a novel scaffold that led to the discovery of potent (38 nM) and efficient inhibitors of this enzyme. Recently, we have advanced these inhibitors to the next stage of development. Using the structures of both the protozoal and human enzymes, we have developed inhibitors with nanomolar potency (1.1 nM) against the pathogenic enzyme and high levels (1273-fold) of selectivity over the human enzyme.

Crystal Structures of HIV-1 Reverse Transcriptase with Picomolar Inhibitors Reveal Key Interactions for Drug Design

X-ray crystal structures at 2.9 Å resolution are reported for two complexes of catechol diethers with HIV-1 reverse transcriptase. The results help elucidate the structural origins of the extreme antiviral activity of the compounds. The possibility of halogen bonding between the inhibitors and Pro95 is addressed. Structural analysis reveals key interactions with conserved residues P95 and W229 of importance for design of inhibitors with high potency and favorable resistance profiles.

Picomolar Inhibitors of HIV Reverse Transcriptase Featuring Bicyclic Replacement of a Cyanovinylphenyl Group

Members of the catechol diether class are highly potent non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs). The most active compounds yield EC50 values below 0.5 nM in assays using human T-cells infected by wild-type HIV-1. However, these compounds such as rilpivirine, the most recently FDA-approved NNRTI, bear a cyanovinylphenyl (CVP) group. This is an uncommon substructure in drugs that gives reactivity concerns. In the present work, computer simulations were used to design bicyclic replacements for the CVP group. The predicted viability of a 2-cyanoindolizinyl alternative was confirmed experimentally and provided compounds with 0.4 nM activity against the wild-type virus. The compounds also performed well with EC50 values of 10 nM against the challenging HIV-1 variant that contains the Lys103Asn/Tyr181Cys double mutation in the RT enzyme. Indolyl and benzofuranyl analogues were also investigated; the most potent compounds in these cases have EC50 values toward wild-type HIV-1 near 10 nM and high-nanomolar activities toward the double-variant. The structural expectations from the modeling were much enhanced by obtaining an X-ray crystal structure at 2.88 Å resolution for the complex of the parent 2-cyanoindolizine 10b and HIV-1 RT. The aqueous solubilities of the most potent indolizine analogues were also measured to be ~40 μg/mL, which is similar to that for the approved drug efavirenz and ~1000-fold greater than for rilpivirine.

Protein design algorithms predict viable resistance to an experimental antifolate.

Significance Computationally predicting drug resistance mutations early in the discovery phase would be an important breakthrough in drug development. The most meaningful predictions of target mutations will show reduced affinity for the drug while maintaining viability in the complex context of a cell. Here, the protein design algorithm K* in Osprey was used to predict a single-nucleotide polymorphism in the target dihydrofolate reductase that confers resistance to an experimental antifolate in the preclinical discovery phase. Excitingly, the mutation was also selected in bacteria under antifolate pressure, confirming the prediction of a viable molecular response to external stress. Methods to accurately predict potential drug target mutations in response to early-stage leads could drive the design of more resilient first generation drug candidates. In this study, a structure-based protein design algorithm (K* in the OSPREY suite) was used to prospectively identify single-nucleotide polymorphisms that confer resistance to an experimental inhibitor effective against dihydrofolate reductase (DHFR) from Staphylococcus aureus. Four of the top-ranked mutations in DHFR were found to be catalytically competent and resistant to the inhibitor. Selection of resistant bacteria in vitro reveals that two of the predicted mutations arise in the background of a compensatory mutation. Using enzyme kinetics, microbiology, and crystal structures of the complexes, we determined the fitness of the mutant enzymes and strains, the structural basis of resistance, and the compensatory relationship of the mutations. To our knowledge, this work illustrates the first application of protein design algorithms to prospectively predict viable resistance mutations that arise in bacteria under antibiotic pressure.

Structure-Based Evaluation of Non-nucleoside Inhibitors with Improved Potency and Solubility That Target HIV Reverse Transcriptase Variants.

The development of novel non-nucleoside inhibitors (NNRTIs) with activity against variants of HIV reverse transcriptase (RT) is crucial for overcoming treatment failure. The NNRTIs bind in an allosteric pocket in RT ∼10 Å away from the active site. Earlier analogues of the catechol diether compound series have picomolar activity against HIV strains with wild-type RT but lose potency against variants with single Y181C and double K103N/Y181C mutations. As guided by structure-based and computational studies, removal of the 5-Cl substitution of compound 1 on the catechol aryl ring system led to a new analogue compound 2 that maintains greater potency against Y181C and K103N/Y181C variants and better solubility (510 μg/mL). Crystal structures were determined for wild-type, Y181C, and K103N/Y181C RT in complex with both compounds 1 and 2 to understand the structural basis for these findings. Comparison of the structures reveals that the Y181C mutation destabilizes the binding mode of compound 1 and disrupts the interactions with residues in the pocket. Compound 2 maintains the same conformation in wild-type and mutant structures, in addition to several interactions with the NNRTI binding pocket. Comparison of the six crystal structures will assist in the understanding of compound binding modes and future optimization of the catechol diether series.

Toward New Therapeutics for Skin and Soft Tissue Infections: Propargyl-Linked Antifolates Are Potent Inhibitors of MRSA and Streptococcus pyogenes

Hospital- and community-acquired, complicated skin and soft tissue infections, often attributed to Staphylococcus aureus and Streptococcus pyogenes, present a significant health burden that is associated with increased health care costs and mortality. As these two species are difficult to discern on diagnosis and are associated with differential profiles of drug resistance, the development of an efficacious antibacterial agent that targets both organisms is a high priority. Herein we describe a structure-based drug development effort that has produced highly potent inhibitors of dihydrofolate reductase from both species. Optimized propargyl-linked antifolates containing a key pyridyl substituent display antibacterial activity against both methicillin-resistant S. aureus and S. pyogenes at MIC values below 0.1 µg/mL and minimal cytotoxicity against mammalian cells. Further evaluation against a panel of clinical isolates shows good efficacy against a range of important phenotypes such as hospital- and community-acquired strains as well as strains resistant to vancomycin.