Search

Antimicrobial resistance presents a significant health care crisis. The mutation F98Y in Staphylococcus aureus dihydrofolate reductase (SaDHFR) confers resistance to the clinically important antifolate trimethoprim (TMP). Propargyl-linked antifolates (PLAs), next generation DHFR inhibitors, are much more resilient than TMP against this F98Y variant, yet this F98Y substitution still reduces efficacy of these agents. Surprisingly, differences in the enantiomeric configuration at the stereogenic center of PLAs influence the isomeric state of the NADPH cofactor. To understand the molecular basis of F98Y-mediated resistance and how PLAs' inhibition drives NADPH isomeric states, we used protein design algorithms in the osprey protein design software suite to analyze a comprehensive suite of structural, biophysical, biochemical, and computational data. Here, we present a model showing how F98Y SaDHFR exploits a different anomeric configuration of NADPH to evade certain PLAs' inhibition, while other PLAs remain unaffected by this resistance mechanism. Author summary: Antimicrobial resistance is a major healthcare crisis. While we were developing novel enzyme inhibitors to combat methicillin-resistant Staphylococcus aureus (MRSA), we found that the chirality of both inhibitor and cofactor can have a large influence on inhibitor potency. Our detailed study of enantiomeric propargyl-linked antifolates (PLAs) shows that the chiral state of inhibitors can affect the chiral state of the cofactor. Moreover, the bacterial enzyme target can exploit cooperative chirality to evade inhibitor binding. We call this phenomenon chiral evasion. Using crystal structures, biochemical assays, computational protein design algorithms, and statistical mechanics, a detailed mechanism for chiral evasion is proposed. While the concept that different enantiomers have different biology is well known, MRSA is unique: we do not know of any other cases where a single mutation (F98Y) flips the chirality preference for cofactor binding and induces stereospecificity for drug binding. Thus, we illuminate the effect of this clinically relevant resistance mutation on the obligate cofactor binding site. These new insights will be useful to develop more durable antibiotics that are resilient to resistance. [ABSTRACT FROM AUTHOR]

The K* algorithm provably approximates partition functions for a set of states (e.g., protein, ligand, and protein-ligand complex) to a user-specified accuracy ε. Often, reaching an ε-approximation for a particular set of partition functions takes a prohibitive amount of time and space. To alleviate some of this cost, we introduce two new algorithms into the osprey suite for protein design: fries, a Fast Removal of Inadequately Energied Sequences, and EWAK*, an Energy Window Approximation to K*. fries pre-processes the sequence space to limit a design to only the most stable, energetically favorable sequence possibilities. EWAK* then takes this pruned sequence space as input and, using a user-specified energy window, calculates K* scores using the lowest energy conformations. We expect fries/EWAK* to be most useful in cases where there are many unstable sequences in the design sequence space and when users are satisfied with enumerating the low-energy ensemble of conformations. In combination, these algorithms provably retain calculational accuracy while limiting the input sequence space and the conformations included in each partition function calculation to only the most energetically favorable, effectively reducing runtime while still enriching for desirable sequences. This combined approach led to significant speed-ups compared to the previous state-of-the-art multi-sequence algorithm, BBK*, while maintaining its efficiency and accuracy, which we show across 40 different protein systems and a total of 2,826 protein design problems. Additionally, as a proof of concept, we used these new algorithms to redesign the protein-protein interface (PPI) of the c-Raf-RBD:KRas complex. The Ras-binding domain of the protein kinase c-Raf (c-Raf-RBD) is the tightest known binder of KRas, a protein implicated in difficult-to-treat cancers. fries/EWAK* accurately retrospectively predicted the effect of 41 different sets of mutations in the PPI of the c-Raf-RBD:KRas complex. Notably, these mutations include mutations whose effect had previously been incorrectly predicted using other computational methods. Next, we used fries/EWAK* for prospective design and discovered a novel point mutation that improves binding of c-Raf-RBD to KRas in its active, GTP-bound state (KRasGTP). We combined this new mutation with two previously reported mutations (which were highly-ranked by osprey) to create a new variant of c-Raf-RBD, c-Raf-RBD(RKY). fries/EWAK* in osprey computationally predicted that this new variant binds even more tightly than the previous best-binding variant, c-Raf-RBD(RK). We measured the binding affinity of c-Raf-RBD(RKY) using a bio-layer interferometry (BLI) assay, and found that this new variant exhibits single-digit nanomolar affinity for KRasGTP, confirming the computational predictions made with fries/EWAK*. This new variant binds roughly five times more tightly than the previous best known binder and roughly 36 times more tightly than the design starting point (wild-type c-Raf-RBD). This study steps through the advancement and development of computational protein design by presenting theory, new algorithms, accurate retrospective designs, new prospective designs, and biochemical validation. Author summary: Computational structure-based protein design is an innovative tool for redesigning proteins to introduce a particular or novel function. One such function is improving the binding of one protein to another, which can increase our understanding of important protein systems. Herein we introduce two novel, provable algorithms, fries and EWAK*, for more efficient computational structure-based protein design as well as their application to the redesign of the c-Raf-RBD:KRas protein-protein interface. These new algorithms speed-up computational structure-based protein design while maintaining accurate calculations, allowing for larger, previously infeasible protein designs. Additionally, using fries and EWAK* within the osprey suite, we designed the tightest known binder of KRas, a heavily studied cancer target that interacts with a number of different proteins. This previously undiscovered variant of a KRas-binding domain, c-Raf-RBD, has potential to serve as a tool to further probe the protein-protein interface of KRas with its effectors and its discovery alone emphasizes the potential for more successful applications of computational structure-based protein design. [ABSTRACT FROM AUTHOR]

Endoplasmic reticulum calcium homeostasis is critical for cellular functions and is disrupted in diverse pathologies including neurodegeneration and cardiovascular disease. Owing to the high concentration of calcium within the ER, studying this subcellular compartment requires tools that are optimized for these conditions. To develop a single-fluorophore genetically encoded calcium indicator for this organelle, we targeted a low affinity variant of GCaMP3 to the ER lumen (GCaMPer (10.19)). A set of viral vectors was constructed to express GCaMPer in human neuroblastoma cells, rat primary cortical neurons, and human induced pluripotent stem cell-derived cardiomyocytes. We observed dynamic changes in GCaMPer (10.19) fluorescence in response to pharmacologic manipulations of the ER calcium store. Additionally, periodic calcium efflux from the ER was observed during spontaneous beating of cardiomyocytes. GCaMPer (10.19) has utility in imaging ER calcium in living cells and providing insight into luminal calcium dynamics under physiologic and pathologic states. [ABSTRACT FROM AUTHOR]