A Real-World Study in Advanced Non–Small Cell Lung Cancer with KRAS Mutations

BACKGROUND: KRAS gene mutations are well known as a key driver of advanced non–small cell lung cancer (NSCLC). The impact of KRAS-mutant subtypes on the survival benefit from salvage chemotherapy is controversial. Here, we present a real-world study in patients across China with advanced NSCLC with KRAS mutations using a website-based patient self-report system. METHODS: We identified a total of 75 patients diagnosed with KRAS-mutant (determined by molecular sequencing) advanced NSCLC between 2014/5/9 and 2019/5/30. KRAS mutation subtypes were divided into G12C and non-G12C groups for statistical analysis. The clinicopathological characteristics and treatment survival benefit in all patients with a KRAS mutation were evaluated. Programmed death-ligand 1 (PD-L1) expression data were collected from 30 patients in the same cohort. RESULTS: In this study, 23 patients with stage IIIB NSCLC and 52 patients with stage IV NSCLC were enrolled with 58 men and 17 women; the median age was 60 years (39–84). All patients received regular chemotherapy/radiotherapy/targeted therapy/immune therapy as per the disease condition. Four main KRAS mutation subtypes were detected: G12C (33%), G12V (19%), G12A (12%), and G12D (12%). Three predominant KRAS comutations were detected: TP53-KRAS (31%), EGFR-KRAS (11%), and STK11-KRAS (8%). Compared with the KRAS non-G12C mutation subtype, patients with the KRAS G12C mutation had potentially longer progression-free survival (PFS) after first-line chemotherapy (4.7 vs. 2.5 months, p < 0.05). Pemetrexed-based chemotherapy appeared to be superior to taxanes- and gemcitabine-based chemotherapies in all patients (PFS: 5.0 vs. 1.5 and 2.3 months, respectively, p > 0.05). Cox regression analysis showed that the KRAS G12C mutation and pemetrexed-based first-line chemotherapy were positive influencers for PFS after first-line (hazard ratios = 0.31 and 0.55, respectively, P < 0.05), but not second-line chemotherapies. CONCLUSION: The KRAS G12C mutation could be a predictive biomarker for better survival benefit from first-line chemotherapy in patients with advanced NSCLC and KRAS mutations. The first-line chemotherapy regimen could possibly influence the outcome in patients with KRAS mutations. Larger and prospective clinical trials are warranted to confirm our conclusions.     

CDK1 Is a Synthetic Lethal Target for KRAS Mutant Tumours

Activating KRAS mutations are found in approximately 20% of human cancers but no RAS-directed therapies are currently available. Here we describe a novel, robust, KRAS synthetic lethal interaction with the cyclin dependent kinase, CDK1. This was discovered using parallel siRNA screens in KRAS mutant and wild type colorectal isogenic tumour cells and subsequently validated in a genetically diverse panel of 26 colorectal and pancreatic tumour cell models. This established that the KRAS/CDK1 synthetic lethality applies in tumour cells with either amino acid position 12 (p.G12V, pG12D, p.G12S) or amino acid position 13 (p.G13D) KRAS mutations and can also be replicated in vivo in a xenograft model using a small molecule CDK1 inhibitor. Mechanistically, CDK1 inhibition caused a reduction in the S-phase fraction of KRAS mutant cells, an effect also characterised by modulation of Rb, a master control of the G₁/S checkpoint. Taken together, these observations suggest that the KRAS/CDK1 interaction is a robust synthetic lethal effect worthy of further investigation.     

Conformational analysis of trimeric maleimide substituted 1,5,9-triazacyclododecane HIV fusion scaffolds

An analysis of the conformational preferences of three trimeric maleimide substituted 1,5,9-triazacyclododecane derivatives, proposed as cross linking reagents for HIV-1 fusion inhibitors, is presented. Exhaustive sampling was performed using the mixed Low Mode Monte Carlo conformational searching technique on the corresponding OPLS2005/GBSA(water) potential energy surface. Geometric structure, molecular length, and hydrogen bonding patterns of the compounds are analyzed. Global minimum energy structures were verified as minima using B3LYP/6-31G1 geometry optimization. All structures adopt a crown-like 12-membered ring conformation; however, the system with the shortest maleimide arms (1a) can also adopt alternative ring orientations. Overall, derivatives with longer maleimide arms were more flexible and resulted in ensembles with a larger number of low energy structures. Comparison with biological inhibition data indicates that there is very little relationship between molecular size and the ability of the scaffold to orient CD4M9 miniproteins for optimal inhibition; however hydrophobicity may play a role.     

Oncogenic K-Ras Binds to an Anionic Membrane in Two Distinct Orientations: A Molecular Dynamics Analysis

K-Ras is a membrane-associated GTPase that cycles between active and inactive conformational states to regulate a variety of cell signaling pathways. Somatic mutations in K-Ras are linked to 15–20% of all human tumors. K-Ras attaches to the inner leaflet of the plasma membrane via a farnesylated polybasic domain; however, the structural details of the complex remain poorly understood. Based on extensive (7.5 μs total) atomistic molecular dynamics simulations here we show that oncogenic mutant K-Ras interacts with a negatively charged lipid bilayer membrane in multiple orientations. Of these, two highly populated orientations account for ∼54% of the conformers whose catalytic domain directly interacts with the bilayer. In one of these orientation states, membrane binding involves helices 3 and 4 of the catalytic domain in addition to the farnesyl and polybasic motifs. In the other orientation, β-strands 1–3 and helix 2 on the opposite face of the catalytic domain contribute to membrane binding. Flexibility of the linker region was found to be important for the reorientation. The biological significance of these observations was evaluated by initial experiments in cells overexpressing mutant K-Ras as well as by an analysis of Ras-effector complex structures. The results suggest that only one of the two major orientation states is capable of effector binding. We propose that the different modes of membrane binding may be exploited in structure-based drug design efforts for cancer therapy.     

Effect of protein matrix on CP29 spectra and energy transfer pathways

We investigate energy transfer pathways between strongly coupled chlorophylls (Chls) in the CP29 (LHCII B4.1) antenna complex of Pisum sativum, including the possibility of higher energy states. We test for the environmental effects caused by the protein, membrane and solvent using a hybrid QM/MM approach. Classical molecular dynamics simulations of the full CP29 complex embedded in a DOPC membrane have been performed, followed by calculations of the time dependent DFT spectra of all Chls at several timesteps. The relative orientations of transition dipole moments (TDMs) were specifically analyzed, including and excluding the point charge field (PCF) of the surrounding environment.The PCF is found to drastically shift the spectra of specific Chls, while the majority of Chls is mostly unaffected. The net effect on the sum spectrum is however found to be negligible: The few strong changes in Chl spectra cancel each other due to being opposite in sign. We further find that the spectra of the Chls coordinating to water show a blue shift upon introduction of the environment. Conversely, the spectra of the Chls coordinating to glutamine show a red shift upon activation of the PCF.As the main influence of the PCF for tuning the couplings, we identify the energetic position of the individual chromophores. The fine-tuning, especially for states energetically above the Q_y state, is however controlled by the changes in the TDM orientations. We also find an indication for the PCF to steer potentially harmful high energy excitations away from the PSII core complex.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

KRAS G12V mutation upregulates PD-L1 expression via TGF-β/EMT signaling pathway in human non-small-cell lung cancer

Although clinical data suggest remarkable promise for targeting programmed cell death protein-1 (PD-1) and ligand (PD-L1) signaling in non-small-cell lung cancer (NSCLC), it is still largely undetermined which subtype of patients will be responsive to checkpoint blockade. In the present study, we explored whether PD-L1 was regulated by mutant Kirsten rat sarcoma viral oncogene homolog (KRAS), which is frequently mutated in NSCLC and results in poor prognosis and low survival rates. We verified that PD-L1 levels were dramatically increased in KRAS mutant cell lines, particularly in NCI-H441 cells with KRAS G12V mutation. Overexpression of KRAS G12V remarkably elevated PD-L1 messenger RNA and protein levels, while suppression of KRAS G12V led to decreased PD-L1 levels in NCI-H441 cells. Consistently, higher levels of PD-L1 were observed in KRAS-mutated tissues as well as tumor tissues-derived CD4⁺ and CD8⁺ T cells using a tumor xenograft in B-NDG mice. Mechanically, both in vitro and in vivo assays found that KRAS G12V upregulated PD-L1 via regulating the progression of epithelial-to-mesenchymal transition (EMT). Moreover, pembrolizumab activated the antitumor activity and decreased tumor growth with KRAS G12V mutated NSCLC. This study demonstrates that KRAS G12V mutation could induce PD-L1 expression and promote immune escape via transforming growth factor-β/EMT signaling pathway in KRAS-mutant NSCLC, providing a potential therapeutic approach for NSCLC harboring KRAS mutations.     

Influence of the orientation of hydroxyl groups on the puckering modes of furanoid rings

The influence of the orientations of HO-2 and HO-3 on the puckering of the ring in (2S,3R)-tetrahydrofuran-2,3-diol (1) and (3R)-tetrahydrofuran-3-ol (2) has been investigated by consistent force-field methods. In the E or W regions of the pseudorotation cycle, there are no local energy minima. The local minima in the N and S regions are affected by different orientations of the hydroxyl groups, with up to 12.2 kJ.mol^?1 for 1 and 6.9 kJ.mol^?1 for 2. If statistical weights of different rotamers are used to calculate N/S equilibrium constants, both conformational states are equally populated. However, if the orientations of the hydroxyl groups are restricted to simulate behaviour of these molecular fragments in DNA double-helices, the equilibrium shifts towards N. The results indicate that, in the vapor state, a unique one-dimensional N-S energy-barrier does not exist and that there is a spectrum of barriers depending on the orientation of the hydroxyl groups.     

Atomistic View of the Conformational Activation of Src Kinase Using the String Method with Swarms-of-Trajectories

The inactive-to-active conformational transition of the catalytic domain of human c-Src tyrosine kinase is characterized using the string method with swarms-of-trajectories with all-atom explicit solvent molecular dynamics simulations. The activation process occurs in two main steps in which the activation loop (A-loop) opens first, followed by the rotation of the αC helix. The computed potential of mean force energy along the activation pathway displays a local minimum, which allows the identification of an intermediate state. These results show that the string method with swarms-of-trajectories is an effective technique to characterize complex and slow conformational transitions in large biomolecular systems.     

Single-stranded DNA scanning and deamination by APOBEC3G cytidine deaminase at single molecule resolution

APOBEC3G (Apo3G) is a single-stranded (ss)DNA cytosine deaminase that eliminates HIV-1 infectivity by converting C → U in numerous small target motifs on the minus viral cDNA. Apo3G deaminates linear ssDNA in vitro with pronounced spatial asymmetry favoring the 3′ → 5′ direction. A similar polarity observed in vivo is believed responsible for initiating localized C → T mutational gradients that inactivate the virus. When compared with double-stranded (ds)DNA scanning enzymes, e.g. DNA glycosylases that excise rare aberrant bases, there is a paucity of mechanistic studies on ssDNA scanning enzymes. Here, we investigate ssDNA scanning and motif-targeting mechanisms for Apo3G using single molecule Förster resonance energy transfer. We address the specific issue of deamination asymmetry within the general context of ssDNA scanning mechanisms and show that Apo3G scanning trajectories, ssDNA contraction, and deamination efficiencies depend on motif sequence, location, and ionic strength. Notably, we observe the presence of bidirectional quasi-localized scanning of Apo3G occurring proximal to a 5′ hot motif, a motif-dependent DNA contraction greatest for 5′ hot > 3′ hot > 5′ cold motifs, and diminished mobility at low salt. We discuss the single molecule Förster resonance energy transfer data in terms of a model in which deamination polarity occurs as a consequence of Apo3G binding to ssDNA in two orientations, one that is catalytically favorable, with the other disfavorable.     

Population shuffling between ground and high energy excited states

Stochastic processes powered by thermal energy lead to protein motions traversing time-scales from picoseconds to seconds. Fundamental to protein functionality is the utilization of these dynamics for tasks such as catalysis, folding, and allostery. A hierarchy of motion is hypothesized to connect and synergize fast and slow dynamics toward performing these essential activities. Population shuffling predicts a “top-down” temporal hierarchy, where slow time-scale conformational interconversion leads to a shuffling of the free energy landscape for fast time-scale events. Until now, population shuffling was only applied to interconverting ground states. Here, we extend the framework of population shuffling to be applicable for a system interconverting between low energy ground and high energy excited states, such as the SH3 domain mutants G48M and A39V/N53P/V55L from the Fyn tyrosine kinase, providing another tool for accessing the structural dynamics of high energy excited states. Our results indicate that the higher energy gauche− rotameric state for the leucine χ₂ dihedral angle contributes significantly to the distribution of rotameric states in both the major and minor forms of the SH3 domain. These findings are corroborated with unrestrained molecular dynamics (MD) simulations on both the major and minor states of the SH3 domain demonstrating high correlations between experimental and back-calculated leucine χ₂ rotameric populations. Taken together, we demonstrate how fast time-scale rotameric side-chain population distributions can be extracted from slow time-scale conformational exchange data further extending the scope and the applicability of the population shuffling model.     

Substrate‐driven conformational changes in ClC‐ec1 observed by fluorine NMR

The CLC ‘Cl⁻ channel'' family consists of both Cl⁻/H⁺ antiporters and Cl⁻ channels. Although CLC channels can undergo large, conformational changes involving cooperativity between the two protein subunits, it has been hypothesized that conformational changes in the antiporters may be limited to small movements localized near the Cl⁻ permeation pathway. However, to date few studies have directly addressed this issue, and therefore little is known about the molecular movements that underlie CLC-mediated antiport. The crystal structure of the Escherichia coli antiporter ClC-ec1 provides an invaluable molecular framework, but this static picture alone cannot depict the protein movements that must occur during ion transport. In this study we use fluorine nuclear magnetic resonance (NMR) to monitor substrate-induced conformational changes in ClC-ec1. Using mutational analysis, we show that substrate-dependent ¹⁹F spectral changes reflect functionally relevant protein movement occurring at the ClC-ec1 dimer interface. Our results show that conformational change in CLC antiporters is not restricted to the Cl⁻ permeation pathway and show the usefulness of ¹⁹F NMR for studying conformational changes in membrane proteins of known structure.     

p.G12C KRAS mutation prevalence in non-small cell lung cancer: Contribution from interregional variability and population substructures among Hispanics

BackgroundThe KRAS exon 2 p. G12C mutation in patients with lung adenocarcinoma has been increasing in relevance due to the development and effectiveness of new treatment medications. Studies around different populations indicate that regional variability between ethnic groups and ancestries could play an essential role in developing this molecular alteration within lung cancer.MethodsIn a prospective and retrospective cohort study on samples from lung adenocarcinoma from 1000 patients from different administrative regions in Colombia were tested for the KRAS p.G12C mutation. An analysis of STR populations markers was conducted to identify substructure contributions to mutation prevalence.ResultsIncluded were 979 patients with a national mean frequency for the KRAS exon 2 p.G12C mutation of 7.97% (95%CI 6.27–9.66%). Variation between regions was also identified with Antioquia reaching a positivity value of 12.7% (95%CI 9.1–16.3%) in contrast to other regions such as Bogota DC (Capital region) with 5.4% (2.7–8.2%) and Bolivar with 2.4% (95%CI 0–7.2%) (p-value = 0.00262). Furthermore, Short tandem repeat population substructures were found for eight markers that strongly yielded association with KRAS exon 2 p.G12C frequency reaching an adjusted R2 of 0.945 and a p-value of < 0.0001.ConclusionsWidespread identification of KRAS exon 2 p.G12C mutations, especially in cases where NGS is not easily achieved is feasible at a population based level that can characterize regional and national patterns of mutation status. Furthermore, this type of mutation prevalence follows a population substructure pattern that can be easily determined by population and ancestral markers such as STR.     

Disentangling the low-energy states of the major light-harvesting complex of plants and their role in photoprotection

The ability to dissipate large fractions of their absorbed light energy as heat is a vital photoprotective function of the peripheral light-harvesting pigment–protein complexes in photosystem II of plants. The major component of this process, known as qE, is characterised by the appearance of low-energy (red-shifted) absorption and fluorescence bands. Although the appearance of these red states has been established, the molecular mechanism, their site and particularly their involvement in qE are strongly debated. Here, room-temperature single-molecule fluorescence spectroscopy was used to study the red emission states of the major plant light-harvesting complex (LHCII) in different environments, in particular conditions mimicking qE. It was found that most states correspond to peak emission at around 700 nm and are unrelated to energy dissipative states, though their frequency of occurrence increased under conditions that mimicked qE. Longer-wavelength emission appeared to be directly related to energy dissipative states, in particular emission beyond 770 nm. The ensemble average of the red emission bands shares many properties with those obtained from previous bulk in vitro and in vivo studies. We propose the existence of at least three excitation energy dissipating mechanisms in LHCII, each of which is associated with a different spectral signature and whose contribution to qE is determined by environmental control of protein conformational disorder. Emission at 700 nm is attributed to a conformational change in the Lut 2 domain, which is facilitated by the conformational change associated with the primary quenching mechanism involving Lut 1.     

Modeling of Hidden Structures Using Sparse Chemical Shift Data from NMR Relaxation Dispersion

NMR relaxation dispersion measurements report on conformational changes occurring on the μs-ms timescale. Chemical shift information derived from relaxation dispersion can be used to generate structural models of weakly populated alternative conformational states. Current methods to obtain such models rely on determining the signs of chemical shift changes between the conformational states, which are difficult to obtain in many situations. Here, we use a “sample and select” method to generate relevant structural models of alternative conformations of the C-terminal-associated region of Escherichia coli dihydrofolate reductase (DHFR), using only unsigned chemical shift changes for backbone amides and carbonyls (¹H, ¹⁵N, and ¹³C′). We find that CS-Rosetta sampling with unsigned chemical shift changes generates a diversity of structures that are sufficient to characterize a minor conformational state of the C-terminal region of DHFR. The excited state differs from the ground state by a change in secondary structure, consistent with previous predictions from chemical shift hypersurfaces and validated by the x-ray structure of a partially humanized mutant of E. coli DHFR (N23PP/G51PEKN). The results demonstrate that the combination of fragment modeling with sparse chemical shift data can determine the structure of an alternative conformation of DHFR sampled on the μs-ms timescale. Such methods will be useful for characterizing alternative states, which can potentially be used for in silico drug screening, as well as contributing to understanding the role of minor states in biology and molecular evolution.     

Cooperative Changes in Solvent Exposure Identify Cryptic Pockets,Switches, and Allosteric Coupling

Proteins are dynamic molecules that undergo conformational changes to a broad spectrum of different excited states. Unfortunately, the small populations of these states make it difficult to determine their structures or functional implications. Computer simulations are an increasingly powerful means to identify and characterize functionally relevant excited states. However, this advance has uncovered a further challenge: it can be extremely difficult to identify the most salient features of large simulation data sets. We reasoned that many functionally relevant conformational changes are likely to involve large, cooperative changes to the surfaces that are available to interact with potential binding partners. To examine this hypothesis, we introduce a method that returns a prioritized list of potentially functional conformational changes by segmenting protein structures into clusters of residues that undergo cooperative changes in their solvent exposure, along with the hierarchy of interactions between these groups. We term these groups exposons to distinguish them from other types of clusters that arise in this analysis and others. We demonstrate, using three different model systems, that this method identifies experimentally validated and functionally relevant conformational changes, including conformational switches, allosteric coupling, and cryptic pockets. Our results suggest that key functional sites are hubs in the network of exposons. As a further test of the predictive power of this approach, we apply it to discover cryptic allosteric sites in two different β-lactamase enzymes that are widespread sources of antibiotic resistance. Experimental tests confirm our predictions for both systems. Importantly, we provide the first evidence, to our knowledge, for a cryptic allosteric site in CTX-M-9 β-lactamase. Experimentally testing this prediction did not require any mutations and revealed that this site exerts the most potent allosteric control over activity of any pockets found in β-lactamases to date. Discovery of a similar pocket that was previously overlooked in the well-studied TEM-1 β-lactamase demonstrates the utility of exposons.     

Conformational selection in G-proteins: lessons from Ras and Rho

The induced fit model has traditionally been invoked to describe the activating conformational change of the monomeric G-proteins, such as Ras and Rho. With this scheme, the presence or absence of the γ-phosphate of GTP leads to an instantaneous switch in conformation. Here we describe atomistic molecular simulations that demonstrate that both Ras and Rho superfamily members harbor an intrinsic susceptibility to sample multiple conformational states in the absence of nucleotide ligand. By comparing the distribution of conformers in the presence and absence of nucleotide, we show that conformational selection is the dominant mechanism by which Ras and Rho undergo nucleotide-dependent conformational changes. Furthermore, the pattern of correlated motions revealed by these simulations predicts a preserved allosteric coupling of the nucleotide-binding site with the membrane interacting C-terminus in both Rho and Ras.     

Thermodynamic dissection of large‐scale domain motions coupled with ligand binding of enzyme I

Domain motions are central to the biological functions of many proteins. The energetics of the motions, however, is often difficult to characterize when motions are coupled with the ligand binding. Here, we determined the thermodynamic parameters of individual domain motions and ligand binding of enzyme I (EI) using strategic domain-deletion mutants that selectively removed particular motions. Upon ligand binding, EI employs two large-scale domain motions, the hinge motion and the swivel motion, to switch between conformational states of distinct domain−domain orientations. Calorimetric analysis of the EI mutants separated the free energy changes of the binding and motions, demonstrating that the unfavorable hinge motion (ΔG = 1.5 kcal mol⁻¹) was driven by the favorable swivel motion (ΔG = −5.2 kcal mol⁻¹). The large free energy differences could be explained by the physicochemical nature of the domain interfaces associated with the motions; the hinge motion employed much narrower interface than the swivel motion without any hydrogen bonds or salt bridges. The small heat capacity further suggested that the packing of the domain interfaces associated with the hinge motion was less compact than that commonly observed in proteins. Lastly, thermodynamic analysis of phosphorylated EI suggests that the domain motions are regulated by the ligand binding and the phosphorylation states. Taken together, the thermodynamic dissection approach illustrates how multiple motions and ligand binding are energetically connected during the functional cycle of EI.     

Membrane-bound KRAS approximates an entropic ensemble of configurations

KRAS4B is a membrane-anchored signaling protein and a primary target in cancer research. Predictions from molecular dynamics simulations that have previously shaped our mechanistic understanding of KRAS signaling disagree with recent experimental results from neutron reflectometry, NMR, and thermodynamic binding studies. To gain insight into these discrepancies, we compare this body of biophysical data to back-calculated experimental results from a series of molecular simulations that implement different subsets of molecular interactions. Our results show that KRAS4B approximates an entropic ensemble of configurations at model membranes containing 30% phosphatidylserine lipids, which is not significantly shaped by interactions between the globular G-domain of KRAS4B and the lipid membrane. These findings revise our understanding of KRAS signaling and promote a model in which the protein samples the accessible conformational space in a near-uniform manner while being available to bind to effector proteins.