Though crucial for understanding the function of large biomolecular systems, locating the minimum free energy paths (MFEPs) between their key conformational states is far from trivial due to their high-dimensional nature. Most existing pathsearching methods require a static collective variable space as input, encoding intuition or prior knowledge of the transition mechanism. Such information is, however, hardly available a priori and expensive to validate. To alleviate this issue, we have previously introduced a Traveling-salesman based Automated Path Searching method (TAPS) and demonstrated its efficiency on simple peptide systems. Having implemented a parallel version of this method, here we assess the performance of TAPS on three realistic systems (tens to hundreds of residues) in explicit solvents. We show that TAPS successfully located the MFEP for the ground/excited state transition of the T4 lysozyme L99A variant, consistent with previous findings. TAPS also helped identifying the important role of the two polar contacts in directing the loop-in/loop-out transition of the mitogen-activated protein kinase kinase (MEK1), which explained previous mutant experiments. Remarkably, at a minimal cost of 126 ns sampling, TAPS revealed that the Ltn40/Ltn10 transition of lymphotactin needs no complete unfolding/ refolding of its β-sheets and that five polar contacts are sufficient to stabilize the various partially unfolded intermediates along the MFEP. These results present TAPS as a general and promising tool for studying the functional dynamics of complex biomolecular systems.

Patients with metastatic acral lentiginous melanoma (ALM) suffer worse outcomes relative to patients with other forms of cutaneous melanoma (CM), and do not benefit as well to approved melanoma therapies. Identification of cyclin-dependent kinase 4 and 6 (CDK4/6) pathway gene alterations in > 60% of ALMs has led to clinical trials of the CDK4/6 inhibitor (CDK4i/6i) palbociclib for ALM; however, median progression free survival with CDK4i/6i treatment was only 2.2 months, suggesting existence of resistance mechanisms. Therapy resistance in ALM remains poorly understood; here we report hyperactivation of MAPK signaling and elevated cyclin D1 expression are a unified mechanism of both intrinsic and acquired CDK4i/6i resistance. MEK and/or ERK inhibition increases CDK4i/6i efficacy in a patient-derived xenograft (PDX) model of ALM and promotes a defective DNA repair, cell cycle arrested and apoptotic program. Notably, gene alterations poorly correlate with protein expression of cell cycle proteins in ALM or efficacy of CDK4i/6i, urging additional strategies when stratifying patients for CDK4i/6i trial inclusion. Concurrent targeting of the MAPK pathway and CDK4/6 represents a new approach to improve outcomes for patients with advanced ALM.

Resistance to combination BRAF/MEK inhibitor (BRAFi/MEKi) therapy arises in nearly every patient with BRAFV600E/K melanoma, despite promising initial responses. Achieving cures in this expanding BRAFi/MEKi-resistant cohort represents one of the greatest challenges to the field; few experience additional durable benefit from immunotherapy and no alternative therapies exist. To better personalize therapy in cancer patients to address therapy relapse, umbrella trials have been initiated whereby genomic sequencing of a panel of potentially actionable targets guide therapy selection for patients; however, the superior efficacy of such approaches remains to be seen. We here test the robustness of the umbrella trial rationale by analyzing relationships between genomic status of a gene and the downstream consequences at the protein level of related pathway, which find poor relationships between mutations, copy number amplification, and protein level. To profile candidate therapeutic strategies that may offer clinical benefit in the context of acquired BRAFi/MEKi resistance, we established a repository of patient-derived xenograft models from heavily pretreated patients with resistance to BRAFi/MEKi and/or immunotherapy (R-PDX). With these R-PDXs, we executed in vivo compound repurposing screens using 11 FDA-approved agents from an NCI-portfolio with pan-RTK, non-RTK and/or PI3K-mTOR specificity. We identify dasatinib as capable of restoring BRAFi/MEKi antitumor efficacy in ~70% of R-PDX tested. A systems-biology analysis indicates elevated baseline protein expression of canonical drivers of therapy resistance (e.g., AXL, YAP, HSP70, phospho-AKT) as predictive of MAPKi/dasatinib sensitivity. We therefore propose that dasatinib-based MAPKi therapy may restore antitumor efficacy in patients that have relapsed to standard-of-care therapy by broadly targeting proteins critical in melanoma therapy escape. Further, we submit that this experimental PDX paradigm could potentially improve preclinical evaluation of therapeutic modalities and augment our ability to identify biomarker-defined patient subsets that may respond to a given clinical trial.

A major risk factor for glaucomatous optic neuropathy is the level of intraocular pressure (IOP), which can lead to retinal ganglion cell axon injury and cell death. The optic nerve has a rostral unmyelinated portion at the optic nerve head followed by a caudal myelinated region. The unmyelinated region is differentially susceptible to IOP–induced damage in rodent models and in human glaucoma. While several studies have analyzed gene expression changes in the mouse optic nerve following optic nerve injury, few were designed to consider the regional gene expression differences that exist between these distinct areas. We performed bulk RNA–sequencing on the retina and on separately micro–dissected unmyelinated and myelinated optic nerve regions from naïve C57BL/6 mice, mice after optic nerve crush, and mice with microbead–induced experimental glaucoma (total = 36). Gene expression patterns in the naïve unmyelinated optic nerve showed significant enrichment of the Wnt, Hippo, PI3K–Akt, and transforming growth factor β pathways, as well as extracellular matrix—receptor and cell membrane signaling pathways, compared to the myelinated optic nerve and retina. Gene expression changes induced by both injuries were more extensive in the myelinated optic nerve than the unmyelinated region, and greater after nerve crush than glaucoma. Changes three and fourteen days after injury largely subsided by six weeks. Gene markers of reactive astrocytes did not consistently differ between injury states. Overall, the transcriptomic phenotype of the mouse unmyelinated optic nerve was significantly different from immediately adjacent tissues, likely dominated by expression in astrocytes, whose junctional complexes are inherently important in responding to IOP elevation.