D. E. Shaw Research

As for man-made machines, the functional cycles of enzymes and other proteins involve the coordinated motion of their parts. These motions, and the mechanical properties of proteins more generally, therefore link protein structure, function, and evolution. For many proteins, experimental access to protein mechanics remains elusive. I will describe our efforts to access these dynamics by electric field-stimulated X-ray crystallography and complementary methods. I will illustrate this approach using two examples: our work on the direct visualization of K+ ion permeation through an ion channel, and on the dissection of structural intermediates enzymatic cycle of dihydrofolate reductase.

Estimating free energy differences, an important problem in computational drug discovery and in a wide range of other application areas, commonly involves a computationally intensive process of sampling a family of high-dimensional probability distributions and a procedure for computing estimates based on those samples. The variance of the free energy estimate of interest typically depends strongly on how the total computational resources available for sampling are divided among the distributions, but determining an efficient allocation is difficult without sampling the distributions. Here we introduce the Times Square sampling algorithm, a novel on-the-fly estimation method that dynamically allocates resources in such a way as to significantly accelerate the estimation of free energies and other observables, while providing rigorous convergence guarantees for the estimators. We also show that it is possible, surprisingly, for on-the-fly free energy estimation to achieve lower asymptotic variance than the maximum-likelihood estimator MBAR, raising the prospect that on-the-fly estimation could reduce variance in a variety of other statistical applications. Supplementary materials for this article are available online.

Protein tyrosine phosphatase 1B (PTP1B) is a negative regulator of the insulin and leptin signaling pathways, making it a highly attractive target for the treatment of type II diabetes. For PTP1B to perform its enzymatic function, a loop referred to as the "WPD loop" must transition between open (catalytically incompetent) and closed (catalytically competent) conformations, which have both been resolved by X-ray crystallography. Although prior studies have established this transition as the rate-limiting step for catalysis, the transition mechanism for PTP1B and other PTPs has been unclear. Here we present an atomically detailed model of WPD loop transitions in PTP1B based on unbiased, long-timescale molecular dynamics simulations and weighted ensemble simulations. We found that a specific WPD loop region─the PDFG motif─acted as the key conformational switch, with structural changes to the motif being necessary and sufficient for transitions between long-lived open and closed states of the loop. Simulations starting from the closed state repeatedly visited open states of the loop that quickly closed again unless the infrequent conformational switching of the motif stabilized the open state. The functional importance of the PDFG motif is supported by the fact that it is well conserved across PTPs. Bioinformatic analysis shows that the PDFG motif is also conserved, and adopts two distinct conformations, in deiminases, and the related DFG motif is known to function as a conformational switch in many kinases, suggesting that PDFG-like motifs may control transitions between structurally distinct, long-lived conformational states in multiple protein families.

This study investigates the effect of water and CO2 exposure on the performance of epoxy-based coatings under conditions commonly found in pipeline applications. The permeability of fusion bonded epoxy (FBE) decreases as CO2 pressure increases and the presence of water facilitates gas transport through the coating. The latter is in conflict with theories of membrane selectivity and competitive transport in gas/vapor systems, in which water is the predominant permeant for its lower kinetic diameter and higher condensability. Our results show that this anomaly was due to the dynamic transformation of permeable channels in the coating structure. Microstructural characterization of FBE after exposure to CO2/H2O mixtures showed that carbonation of wollastonite filler particles results in a change in shape, volume, and chemical composition of these filler particles. The chemical reaction between wollastonite and CO2 in the presence of water led to a degradation that caused changes in porosity, permeability, and filler volume. To verify these findings for the coating in the presence of adhesion forces, electrochemical impedance spectroscopy (EIS) was also used. The tests conducted showed that exposure to a CO2/H2O mixture causes filler particles within the coating to undergo carbonation, leading to the formation of new microchannels within the epoxy network. Initially, the carbonation of fillers led to an increase in the pore resistance of the coating, which was attributed to the plugging of micropore channels on the coating surface. However, the subsequent decreasing trend of this parameter suggested that water infiltration into the coating had increased due to this degradation. The carbonation of wollastonite within the FBE coating leads to a volume change in the fillers, resulting in the creation of small voids or gaps in the epoxy matrix. This, in turn, facilitates the easier penetration of water and dissolved gas to the underlying substrate. Consequently, the accumulation of water inside coating increases the dielectric constant, resulting in a higher capacitance of the coating. It appears that high concentrations of CO2 in wet conditions, even at low pressures, can have negative impacts on the barrier performance of the coating.

Fragment-based drug discovery has led to six approved drugs, but the small sizes of the chemical fragments used in such methods typically result in only weak interactions between the fragment and its target molecule, which makes it challenging to experimentally determine the three-dimensional poses fragments assume in the bound state. One computational approach that could help address this difficulty is long-timescale molecular dynamics (MD) simulations, which have been used in retrospective studies to recover experimentally known binding poses of fragments. Here, we present the results of long-timescale MD simulations that we used to prospectively discover binding poses for two series of fragments in allosteric pockets on a difficult and important pharmaceutical target, protein tyrosine phosphatase 1b (PTP1b). Our simulations reversibly sampled the fragment association and dissociation process. One of the binding pockets found in the simulations has not to our knowledge been previously observed with a bound fragment, and the other pocket adopted a very rare conformation. We subsequently obtained high-resolution crystal structures of members of each fragment series bound to PTP1b, and the experimentally observed poses confirmed the simulation results. To the best of our knowledge, our findings provide the first demonstration that MD simulations can be used prospectively to determine fragment binding poses to previously unidentified pockets.

<div>Abstract<p>Current clinical RAF inhibitors (RAFi) inhibit monomeric BRAF (mBRAF) but are less potent against dimeric BRAF (dBRAF). RAFi equipotent for mBRAF and dBRAF have been developed but are predicted to have lower therapeutic index. Here we identify a third class of RAFi that selectively inhibits dBRAF over mBRAF. Molecular dynamic simulations reveal restriction of the movement of the BRAF αC-helix as the basis of inhibitor selectivity. Combination of inhibitors based on their conformation selectivity (mBRAF- plus dBRAF-selective plus the most potent BRAF–MEK disruptor MEK inhibitor) promoted suppression of tumor growth in BRAF<sup>V600E</sup> therapy–resistant models. Strikingly, the triple combination showed no toxicities, whereas dBRAF-selective plus MEK inhibitor treatment caused weight loss in mice. Finally, the triple combination achieved durable response and improved clinical well-being in a patient with stage IV colorectal cancer. Thus, exploiting allosteric properties of RAF and MEK inhibitors enables the design of effective and well-tolerated therapies for BRAF<sup>V600E</sup> tumors.</p>Significance:<p>This work identifies a new class of RAFi that are selective for dBRAF over mBRAF and determines the basis of their selectivity. A rationally designed combination of RAF and MEK inhibitors based on their conformation selectivity achieved increased efficacy and a high therapeutic index when used to target BRAF<sup>V600E</sup> tumors.</p><p><i>See related commentary by Zhang and Bollag, p. 1620</i>.</p><p><i>This article is highlighted in the In This Issue feature, p. 1601</i></p></div>

Abstract Fragment-based drug discovery has led to six approved drugs, but the small size of the chemical fragments used in such methods typically results in only weak interactions between the fragment and its target molecule, which makes it challenging to experimentally determine the three-dimensional poses fragments assume in the bound state. One computational approach that could help address this difficulty is long-timescale molecular dynamics (MD) simulation, which has been used in retrospective studies to recover experimentally known binding poses of fragments. Here, we present the results of long-timescale MD simulations that we used to prospectively discover binding poses for two series of fragments in allosteric pockets on a difficult and important pharmaceutical target, protein-tyrosine phosphatase 1b (PTP1b). Our simulations reversibly sampled the fragment association and dissociation process. One of the binding pockets found in the simulations has not to our knowledge been previously observed with a bound fragment, and the other pocket adopted a very rare conformation. We subsequently obtained high-resolution crystal structures of members of each fragment series bound to PTP1b, and the experimentally observed poses confirmed the simulation results. To the best of our knowledge, our findings provide the first demonstration that MD simulations can be used prospectively to determine fragment binding poses to previously unidentified pockets.

Abstract The non-receptor protein tyrosine phosphatase SHP2 (PTPN11) plays an important role in the regulation of RAS/MAPK signal transduction downstream of growth factor receptor activation. Loss of SHP2 activity suppresses tumor cell growth, making SHP2 a potential target for cancer therapy. Here we report the discovery of GDC-1971 (formerly RLY-1971), a highly potent, selective, and orally bioavailable small-molecule SHP2 inhibitor that stabilizes SHP2 in a closed, auto-inhibited conformation. GDC-1971 potently inhibits both wild-type SHP2 (IC50 <1nM) and the E76K activating mutant (IC50 <250nM) in biochemical assays. In standard 2-dimensional and anchorage-independent growth conditions, GDC-1971 inhibits cellular proliferation in models harboring receptor tyrosine kinases (RTKs), SHP2, NF1, KRAS, or BRAF mutations in a dose-dependent manner. GDC-1971 potently inhibits the proliferation of cellular models harboring KRAS G12C or G12A mutations (median IC50 <80 nM) compared to models harboring other KRAS G12, G13 or Q61 mutations (median IC50 >1 uM), indicating a link between KRAS GTP hydrolysis and SHP2 dependency. Despite this trend, some non-KRAS G12C or G12A cell lines harboring other KRAS mutations responded to GDC-1971 in vitro, suggesting some heterogeneity in RTK/SHP2 signaling dependence in subsets of other KRAS mutants. In vivo, GDC-1971 demonstrates dose-dependent RAS/MAPK pathway inhibition and induces significant tumor-growth inhibition in human xenograft models harboring EGFR and KRAS alterations at continuous daily doses that are well tolerated. Given the reported role of SHP2 as a critical mediator of resistance to targeted therapies, we assessed the activity of GDC-1971 combinations in multiple contexts. We observed increased suppression of the MAPK signaling cascade and anti-proliferation synergy when combining GDC-1971 with EGFR, ALK, and KRAS G12C inhibitors in vitro. The observed in vitro synergy translated to dramatic anti-tumor growth effects in vivo. GDC-1971 in combination with the KRAS G12C covalent inhibitor GDC-6036 resulted in significant regressions at doses well below those required for single agent activity in a KRAS G12C-mutant NSCLC xenograft model. In rodent and dog toxicology studies, GDC-1971 is well tolerated at exposures above those required to induce regression in xenograft models. The biochemical and cellular potency and favorable pharmaceutical properties of GDC-1971 support the further clinical development in RTK/MAPK pathway altered tumors using continuous daily dosing alone and in combination with other targeted agents, including the KRAS-G12C inhibitor GDC-6036 (clinical trial NCT04449874). Citation Format: Bret Williams, Alexander Taylor, Olivia Orozco, Christopher Owen, Elizabeth Kelley, Andre Lescarbeau, Kelley Shortsleeves, Randy Kipp, Vy Nguyen, Erin Brophy, Jeremy Wilbur, Yong Tang, David Lanzetta, Nigel Waters, Sherri Smith, Fabrizio Giordanetto, Paul Maragakis, Jack Greismann, Lindsay Willmore, Eric Therrien, Yang Xiao, Marie Evangelista, Luca Gerosa, Eva Lin, Mark Merchant, Alfonso Arrazate, Emily Chan, Pablo Sáenz-López Larrocha, Stefan Chun, Thomas Hunsaker, Gauri Deshmukh, Christine M. Bowman, David E. Shaw, Mark Murcko, Mahesh Padval, W Patrick Walters, James Watters, Donald A. Bergstrom. Discovery and characterization of the potent, allosteric SHP2 inhibitor GDC-1971 for the treatment of RTK/RAS driven tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3327.