Ligand-induced Conformational Changes Research Articles

Functional protein dynamics in a crystal.

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

Abstract 1942 Investigating Ligand-Induced Local Conformational Changes of Fluorescently Labeled G-Quadruplex Structures

Abstract 1942 Investigating Ligand-Induced Local Conformational Changes of Fluorescently Labeled G-Quadruplex Structures

Structural basis for excitatory neuropeptide signaling

Rapid signaling between neurons is mediated by ligand-gated ion channels, cell-surface proteins with an extracellular ligand-binding domain and a membrane-spanning ion channel domain. The degenerin/epithelial sodium channel (DEG/ENaC) superfamily is diverse in terms of its gating stimuli, with some DEG/ENaCs gated by neuropeptides, and others gated by pH, mechanical force or enzymatic activity. The mechanism by which ligands bind to and activate DEG/ENaCs is poorly understood. Here we dissected the structural basis for neuropeptide-gated activity of a neuropeptide-gated DEG/ENaC, FMRFamide-gated sodium channel 1 (FaNaC1) from the annelid worm Malacoceros fuliginosus, using cryo-electron microscopy. Structures of FaNaC1 in the ligand-free resting state and in several ligand-bound states reveal the ligand-binding site and capture the ligand-induced conformational changes of channel gating, which we verified with complementary mutagenesis experiments. Our results illuminate channel gating in DEG/ENaCs and offer a structural template for experimental dissection of channel pharmacology and ion conduction.

Quantifying a light-induced energetic change in bacteriorhodopsin by force spectroscopy

Ligand-induced conformational changes are critical to the function of many membrane proteins and arise from numerous intramolecular interactions. In the photocycle of the model membrane protein bacteriorhodopsin (bR), absorption of a photon by retinal triggers a conformational cascade that results in pumping a proton across the cell membrane. While decades of spectroscopy and structural studies have probed this photocycle in intricate detail, changes in intramolecular energetics that underlie protein motions have remained elusive to experimental quantification. Here, we measured these energetics on the millisecond time scale using atomic-force-microscopy-based single-molecule force spectroscopy. Precisely, timed light pulses triggered the bR photocycle while we measured the equilibrium unfolding and refolding of the terminal 8-amino-acid region of bR's G-helix. These dynamics changed when the EF-helix pair moved ~9 Å away from this end of the G helix during the "open" portion of bR's photocycle. In ~60% of the data, we observed abrupt light-induced destabilization of 3.4 ± 0.3 kcal/mol, lasting 38 ± 3 ms. The kinetics and pH-dependence of this destabilization were consistent with prior measurements of bR's open phase. The frequency of light-induced destabilization increased with the duration of illumination and was dramatically reduced in the triple mutant (D96G/F171C/F219L) thought to trap bR in its open phase. In the other ~40% of the data, photoexcitation unexpectedly stabilized a longer-lived putative misfolded state. Through this work, we establish a general single-molecule force spectroscopy approach for measuring ligand-induced energetics and lifetimes in membrane proteins.

Structural Insights into the Activation of Human Aryl Hydrocarbon Receptor by the Environmental Contaminant Benzo[a]pyrene and Structurally Related Compounds

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor belonging to the bHLH/PAS protein family and responding to hundreds of natural and chemical substances. It is primarily involved in the defense against chemical insults and bacterial infections or in the adaptive immune response, but also in the development of pathological conditions ranging from inflammatory to neoplastic disorders. Despite its prominent roles in many (patho)physiological processes, the lack of high-resolution structural data has precluded for thirty years an in-depth understanding of the structural mechanisms underlying ligand-binding specificity, promiscuity and activation of AHR. We recently reported a cryogenic electron microscopy (cryo-EM) structure of human AHR bound to the natural ligand indirubin, the chaperone Hsp90 and the co-chaperone XAP2 that provided the first experimental visualization of its ligand-binding PAS-B domain. Here, we report a 2.75 Å resolution structure of the AHR complex bound to the environmental pollutant benzo[a]pyrene (B[a]P). The structure substantiates the existence of a bipartite PAS-B ligand-binding pocket with a geometrically constrained primary binding site controlling ligand binding specificity and affinity, and a secondary binding site contributing to the binding promiscuity of AHR. We also report a docking study of B[a]P congeners that validates the B[a]P-bound PAS-B structure as a suitable model for accurate computational ligand binding assessment. Finally, comparison of our agonist-bound complex with the recently reported structures of mouse and fruit fly AHR PAS-B in different activation states suggests a ligand-induced loop conformational change potentially involved in the regulation of AHR function.

Structural Characterization and Ligand-Induced Conformational Changes of SenB, a Se-Glycosyltransferase Involved in Selenoneine Biosynthesis.

Selenium (Se) is an essential micronutrient that is found naturally in proteins, nucleic acids, and natural products. Unlike selenoproteins and selenonucleic acids, little is known about the structures of biosynthetic enzymes that incorporate Se into small molecules. Here, we report the X-ray crystal structure of SenB, the first known Se-glycosyltransferase that was recently found to be involved in the biosynthesis of the Se-containing metabolite selenoneine. SenB catalyzes C-Se bond formation using selenophosphate and an activated uridine diphosphate sugar as a Se and glycosyl donor, respectively, making it the first known selenosugar synthase and one of only four bona fide C-Se bond-forming enzymes discovered to date. Our crystal structure, determined to 2.25 Å resolution, reveals that SenB is a type B glycosyltransferase, displaying the prototypical fold with two globular Rossmann-like domains and a catalytic interdomain cleft. By employing complementary structural biology techniques, we find that SenB undergoes both local and global substrate-induced conformational changes, demonstrating a significant increase in α-helicity and a transition to a more compact conformation. Our results provide the first structure of SenB and set the stage for further biochemical characterization in the future.

Drug-Induced Conformational Dynamics of P-Glycoprotein Underlies the Transport of Camptothecin Analogs.

P-glycoprotein (Pgp) plays a pivotal role in drug bioavailability and multi-drug resistance development. Understanding the protein's activity and designing effective drugs require insight into the mechanisms underlying Pgp-mediated transport of xenobiotics. In this study, we investigated the drug-induced conformational changes in Pgp and adopted a conformationally-gated model to elucidate the Pgp-mediated transport of camptothecin analogs (CPTs). While Pgp displays a wide range of conformations, we simplified it into three model states: 'open-inward', 'open-outward', and 'intermediate'. Utilizing acrylamide quenching of Pgp fluorescence as a tool to examine the protein's tertiary structure, we observed that topotecan (TPT), SN-38, and irinotecan (IRT) induced distinct conformational shifts in the protein. TPT caused a substantial shift akin to AMPPNP, suggesting ATP-independent 'open-outward' conformation. IRT and SN-38 had relatively moderate effects on the conformation of Pgp. Experimental atomic force microscopy (AFM) imaging supports these findings. Further, the rate of ATPase hydrolysis was correlated with ligand-induced Pgp conformational changes. We hypothesize that the separation between the nucleotide-binding domains (NBDs) creates a conformational barrier for substrate transport. Substrates that reduce the conformational barrier, like TPT, are better transported. The affinity for ATP extracted from Pgp-mediated ATP hydrolysis kinetics curves for TPT was about 2-fold and 3-fold higher than SN-38 and IRT, respectively. On the contrary, the dissociation constants (KD) determined by fluorescence quenching for these drugs were not significantly different. Saturation transfer double difference (STDD) NMR of TPT and IRT with Pgp revealed that similar functional groups of the CPTs are accountable for Pgp-CPTs interactions. Efforts aimed at modifying these functional groups, guided by available structure-activity relationship data for CPTs and DNA-Topoisomerase-I complexes, could pave the way for the development of more potent next-generation CPTs.

Repurposing conformational changes in ANL superfamily enzymes to rapidly generate biosensors for organic and amino acids

Biosensors are powerful tools for detecting, real-time imaging, and quantifying molecules, but rapidly constructing diverse genetically encoded biosensors remains challenging. Here, we report a method to rapidly convert enzymes into genetically encoded circularly permuted fluorescent protein-based indicators to detect organic acids (GECFINDER). ANL superfamily enzymes undergo hinge-mediated ligand-coupling domain movement during catalysis. We introduce a circularly permuted fluorescent protein into enzymes hinges, converting ligand-induced conformational changes into significant fluorescence signal changes. We obtain 11 GECFINDERs for detecting phenylalanine, glutamic acid and other acids. GECFINDER-Phe3 and GECFINDER-Glu can efficiently and accurately quantify target molecules in biological samples in vitro. This method simplifies amino acid quantification without requiring complex equipment, potentially serving as point-of-care testing tools for clinical applications in low-resource environments. We also develop a GECFINDER-enabled droplet-based microfluidic high-throughput screening method for obtaining high-yield industrial strains. Our method provides a foundation for using enzymes as untapped blueprint resources for biosensor design, creation, and application.

Natural protein engineering in the Ω-loop: the role of Y221 in ceftazidime and ceftolozane resistance in Pseudomonas-derived cephalosporinase.

A wide variety of clinically observed single amino acid substitutions in the Ω-loop region have been associated with increased minimum inhibitory concentrations and resistance to ceftazidime (CAZ) and ceftolozane (TOL) in Pseudomonas-derived cephalosporinase and other class C β-lactamases. Herein, we demonstrate the naturally occurring tyrosine to histidine substitution of amino acid 221 (Y221H) in Pseudomonas-derived cephalosporinase (PDC) enables CAZ and TOL hydrolysis, leading to similar kinetic profiles (k cat = 2.3 ± 0.2 µM and 2.6 ± 0.1 µM, respectively). Mass spectrometry of PDC-3 establishes the formation of stable adducts consistent with the formation of an acyl enzyme complex, while spectra of E219K (a well-characterized, CAZ- and TOL-resistant comparator) and Y221H are consistent with more rapid turnover. Thermal denaturation experiments reveal decreased stability of the variants. Importantly, PDC-3, E219K, and Y221H are all inhibited by avibactam and the boronic acid transition state inhibitors (BATSIs) LP06 and S02030 with nanomolar IC50 values and the BATSIs stabilize all three enzymes. Crystal structures of PDC-3 and Y221H as apo enzymes and complexed with LP06 and S02030 (1.35-2.10 Å resolution) demonstrate ligand-induced conformational changes, including a significant shift in the position of the sidechain of residue 221 in Y221H (as predicted by enhanced sampling well-tempered metadynamics simulations) and extensive hydrogen bonding between the enzymes and BATSIs. The shift of residue 221 leads to the expansion of the active site pocket, and molecular docking suggests substrates orientate differently and make different intermolecular interactions in the enlarged active site compared to the wild-type enzyme.

Functional GPCR Expression in Eukaryotic LEXSY System

G protein-coupled receptors (GPCRs) form the largest superfamily of membrane proteins in the human genome, and represent one of the most important classes of drug targets. Their structural studies facilitate rational drug discovery. However, atomic structures of only about 20% of human GPCRs have been solved to date. Recombinant production of GPCRs for structural studies at a large scale is challenging due to their low expression levels and stability. Therefore, in this study, we explored the efficacy of the eukaryotic system LEXSY (Leishmania tarentolae) for GPCR production. We selected the human A2A adenosine receptor (A2AAR), as a model protein, expressed it in LEXSY, purified it, and compared with the same receptor produced in insect cells, which is the most popular expression system for structural studies of GPCRs. The A2AAR purified from both expression systems showed similar purity, stability, ligand-induced conformational changes and structural dynamics, with a remarkably higher protein yield in the case of LEXSY expression. Overall, our results suggest that LEXSY is a promising platform for large-scale production of GPCRs for structural studies.

The nuclear receptor LRH-1 discriminates between ligands using distinct allosteric signaling circuits.

Nuclear receptors (NRs) are transcription factors that regulate essential biological processes in response to cognate ligands. An important part of NR function involves ligand-induced conformational changes that recruit coregulator proteins to the activation function surface (AFS), ~15 Å away from the ligand binding pocket. Ligands must communicate with the AFS to recruit appropriate coregulators and elicit different transcriptional outcomes, but this communication is poorly understood. These studies illuminate allosteric communication networks underlying activation of liver receptor homolog-1 (LRH-1), a NR that regulates development, metabolism, cancer progression and intestinal inflammation. Using >100 microseconds of all-atom molecular dynamics simulations involving 74 LRH-1 complexes, we identify distinct signaling circuits used by active and inactive ligands for AFS communication. Inactive ligands communicate via strong, coordinated motions along paths through the receptor to the AFS. Activating ligands disrupt the "inactive" circuit and induce connectivity with a second allosteric site. Ligand-contacting residues in helix 7 help mediate the switch between circuits, suggesting new avenues for developing LRH-1-targeted therapeutics. We also elucidate aspects of coregulator signaling, showing that localized, destabilizing fluctuations are induced by inappropriate ligand-coregulator pairings. These studies have uncovered novel features of LRH-1 allostery, and the quantitative approach used to analyze many simulations provides a framework to study allosteric signaling in other receptors. This article is protected by copyright. All rights reserved.

Capturing heterogeneous conformers of cobalamin riboswitch by cryo-EM.

RNA conformational heterogeneity often hampers its high-resolution structure determination, especially for large and flexible RNAs devoid of stabilizing proteins or ligands. The adenosylcobalamin riboswitch exhibits heterogeneous conformations under 1 mM Mg2+ concentration and ligand binding reduces conformational flexibility. Among all conformers, we determined one apo (5.3 Å) and four holo cryo-electron microscopy structures (overall 3.0-3.5 Å, binding pocket 2.9-3.2 Å). The holo dimers exhibit global motions of helical twisting and bending around the dimer interface. A backbone comparison of the apo and holo states reveals a large structural difference in the P6 extension position. The central strand of the binding pocket, junction 6/3, changes from an 'S'- to a 'U'-shaped conformation to accommodate ligand. Furthermore, the binding pocket can partially form under 1 mM Mg2+ and fully form under 10 mM Mg2+ within the bound-like structure in the absence of ligand. Our results not only demonstrate the stabilizing ligand-induced conformational changes in and around the binding pocket but may also provide further insight into the role of the P6 extension in ligand binding and selectivity.

Highly Sensitive and Specific Detection of Influenza A Viruses Using Bimolecular Fluorescence Complementation (BiFC) Reporter System.

In this study, we developed a highly sensitive and specific bimolecular fluorescence complementation (BiFC)-based influenza A virus (IAV)-sensing system by combining a galactose/glucose-binding protein (GGBP) with an N-terminal large domain (YN1-172) and a C-terminal small domain (YC173-239) made up of enhanced yellow fluorescence protein (eYFP). The GGBP-based BiFC reporter exhibits the fluorescence reconstitution as a result of conformational changes in GGBP when lactose, which was derived from 6'-silalyllactose and used as a substrate for neuraminidase (NA), binds to GGBP in the presence of IAV. The system showed a linear dynamic range extending from 1 × 100 to 1 × 107 TCID50/mL, and it had a detection limit of 1.1 × 100 TCID50/mL for IAV (H1N1), demonstrating ultra-high sensitivity. Our system exhibited fluorescence intensity enhancements in the presence of IAV, while it displayed weak fluorescence signals when exposed to NA-deficient viruses, such as RSV A, RSV B, adenovirus and rhinovirus, thereby indicating selective responses for IAV detection. Overall, our system provides a simple, highly sensitive and specific IAV detection platform based on BiFC that is capable of detecting ligand-induced protein conformational changes, obviating the need for virus culture or RNA extraction processes.

Linking folding dynamics and function of SAM/SAH riboswitches at the single molecule level.

Riboswitches are regulatory elements found in bacterial mRNAs that control downstream gene expression through ligand-induced conformational changes. Here, we used single-molecule FRET to map the conformational landscape of the translational SAM/SAH riboswitch and probe how co-transcriptional ligand-induced conformational changes affect its translation regulation function. Riboswitch folding is highly heterogeneous, suggesting a rugged conformational landscape that allows for sampling of the ligand-bound conformation even in the absence of ligand. The addition of ligand shifts the landscape, favoring the ligand-bound conformation. Mutation studies identified a key structural element, the pseudoknot helix, that is crucial for determining ligand-free conformations and their ligand responsiveness. We also investigated ribosomal binding site accessibility under two scenarios: pre-folding and co-transcriptional folding. The regulatory function of the SAM/SAH riboswitch involves kinetically favoring ligand binding, but co-transcriptional folding reduces this preference with a less compact initial conformation that exposes the Shine-Dalgarno sequence and takes min to redistribute to more compact conformations of the pre-folded riboswitch. Such slow equilibration decreases the effective ligand affinity. Overall, our study provides a deeper understanding of the complex folding process and how the riboswitch adapts its folding pattern in response to ligand, modulates ribosome accessibilityand the role of co-transcriptional folding in these processes.

The effects of allosteric and competitive inhibitors on ZIKV protease conformational dynamics explored through smFRET, nanoDSF, DSF, and 19F NMR

ABSTRACT: Zika and dengue viruses cause mosquito-borne diseases of high epidemic relevance. The viral NS2B-NS3 proteases play crucial roles in the pathogen replication cycle and are validated drug targets. They can adopt at least two conformations depending on the position of the NS2B cofactor. Recently, we reported ligand-induced conformational changes of dengue virus NS2B-NS3 protease by single-molecule Förster resonance energy transfer (smFRET). Here, we investigated the conformational dynamics of the homologous Zika virus protease through an integrated methodological approach combining smFRET, thermal shift assays (DSF and nanoDSF) and 19F NMR spectroscopy. Our results show that allosteric inhibitors favor the open conformation and competitive inhibitors stabilize the closed conformation of the Zika virus protease.

Crystal structure of Escherichia coli thiamine pyrophosphate-sensing riboswitch in the apo state

The thiamine pyrophosphate (TPP)-sensing riboswitch is one of the earliest discovered and most widespread riboswitches. Numerous structural studies have been reported for this riboswitch bound with various ligands. However, the ligand-free (apo) structure remains unknown. Here, we report a 3.1 Å resolution crystal structure of Escherichia coli TPP riboswitch in the apo state, which exhibits an extended, Y-shaped conformation further supported by small-angle X-ray scattering data and driven molecular dynamics simulations. The loss of ligand interactions results in helical uncoiling of P5 and disruption of the key tertiary interaction between the sensory domains. Opening of the aptamer propagates to the gene-regulatory P1 helix and generates the key conformational flexibility needed for the switching behavior. Much of the ligand-binding site at the three-way junction is unaltered, thereby maintaining a partially preformed pocket. Together, these results paint a dynamic picture of the ligand-induced conformational changes in TPP riboswitches that confer conditional gene regulation.

Analysis of the Binding Mechanism of Bioactive Coumarins with Ovalbumin: Further Investigation into the Inhibitory Effects toward Protein Fibrillation

Two bioactive coumarin derivatives, 7-hydroxycoumarin (7-HC) and 4-methyl-7-hydroxycoumarin (4-Me-7-HC), were used to investigate the molecular interactions with the main food allergen from egg white, ovalbumin (OVA), along with their inhibitory effects on OVA aggregation. The mechanism of binding interaction was elucidated by multispectroscopic and computational methods. Circular dichroism (CD) and Fourier transform infrared (FT-IR) experiments confirmed the ligand-induced conformational changes of OVA. The fluorescence spectroscopic experiments demonstrated that the quenching mechanism was an unusual static quenching method with binding constants (Kb) in the range of 104 M–1, indicating a moderate nature of binding between OVA and coumarins. Thioflavin T (ThT) and Congo red (CR) binding assays confirmed the production of OVA fibrils, and the aggregation was shown to be inhibited by the coumarin derivatives in vitro. Thermodynamic parameters for OVA-7-HC/4-Me-7-HC interactions showed positive ΔH and ΔS values, indicating hydrophobic forces predominated in the binding processes and a negative ΔG value, implying the spontaneity of the complex formation, which was then further confirmed by computational studies.

The START domain potentiates HD-ZIPIII transcriptional activity.

The CLASS III HOMEODOMAIN-LEUCINE ZIPPER (HD-ZIPIII) transcription factors (TFs) were repeatedly deployed over 725 million years of evolution to regulate central developmental innovations. The START domain of this pivotal class of developmental regulators was recognized over 20 years ago, but its putative ligands and functional contributions remain unknown. Here, we demonstrate that the START domain promotes HD-ZIPIII TF homodimerization and increases transcriptional potency. Effects on transcriptional output can be ported onto heterologous TFs, consistent with principles of evolution via domain capture. We also show the START domain binds several species of phospholipids, and that mutations in conserved residues perturbing ligand binding and/or its downstream conformational readout abolish HD-ZIPIII DNA-binding competence. Our data present a model in which the START domain potentiates transcriptional activity and uses ligand-induced conformational change to render HD-ZIPIII dimers competent to bind DNA. These findings resolve a long-standing mystery in plant development and highlight the flexible and diverse regulatory potential coded within this widely distributed evolutionary module.

Dissecting the Thermodynamics of ATP Binding to GroEL One Nucleotide at a Time.

Variable-temperature electrospray ionization (vT-ESI) native mass spectrometry (nMS) is used to determine the thermodynamics for stepwise binding of up to 14 ATP molecules to the 801 kDa GroEL tetradecamer chaperonin complex. Detailed analysis reveals strong enthalpy-entropy compensation (EEC) for the ATP binding events leading to formation of GroEL-ATP7 and GroEL-ATP14 complexes. The observed variations in EEC and stepwise free energy changes of specific ATP binding are consistent with the well-established nested cooperativity model describing GroEL-ATP interactions, viz., intraring positive cooperativity and inter-ring negative cooperativity (Dyachenko A.; Proc. Natl. Acad. Sci. U.S.A.2013, 110, 7235-7239). Entropy-driven ATP binding is to be expected for ligand-induced conformational changes of the GroEL tetradecamer, though the magnitude of the entropy change suggests that reorganization of GroEL-hydrating water molecules and/or expulsion of water from the GroEL cavity may also play key roles. The capability for determining complete thermodynamic signatures (ΔG, ΔH, and -TΔS) for individual ligand binding reactions for the large, nearly megadalton GroEL complex expands our fundamental view of chaperonin functional chemistry. Moreover, this work and related studies of protein-ligand interactions illustrate important new capabilities of vT-ESI-nMS for thermodynamic studies of protein interactions with ligands and other molecules such as proteins and drugs.

An in silico comparative study of curcumin and 2-deoxyuridine nucleoside derivatives: Reveals the role of angiogenin in ER stress-induced apoptosis signaling.

Angiogenin (ANG) protein plays a crucial role in angiogenesis, neovascularization, and cancer metastasis in NSCLC (non-small cell lung cancer) via non-coding tiRNA. It protects the cell under ER (endoplasmic reticulum) stress-induced apoptosis through the translational reprogramming process. Although B82 (Curcumin derivatives) induces ER stress-induced apoptosis, its mechanism of action was not studied. Therefore, it was hypothesized that the ribonucleolytic activity of ANG may be regulated by B82, resulting in modulated ER stress signaling for apoptosis. Hence, we designed and proposed a synthesis scheme for RNA-based anti-angiogenic derivatives of 2-deoxyuridine nucleoside forming peptide bond with amino acids like serine (Ser-3) and para-hydroxy-phenyl glycine (Normtyr-1) and compared B82 with them to know the binding affinity with ANG, anti-angiogenic potential, and its probable mechanism of anti-RNase activity through MD simulation study. Therefore, using Gromos96 43a1 and 43a2 force fields, MD simulation was performed to investigate binding affinity, ligand-induced molecular surface area change, conformational change, and dynamics of catalytic site residues to predict ligand binding to ANG in this study. The obtained binding free energy (∆Gbind ) result showed the total average ∆Gbind as -113.480 ± 1.682 (Normtyr-1) > -53.038 ± 33.069 (B82) > -27.909 ± 16.438 (Ser-3)kJ/mole specify role of B82 in regulating ER stress signaling induced apoptosis through ANG ribonucleolytic activity inhibition, suitability of 43a2 force fields and methodology in ligand screening. It shows the crucial role of Leu115 and His13 residue involvement in total ∆Gbind contribution. Hence, based on the MD result, novel conformation of catalytic residues, and ∆Gbind , a promising combination candidate could be proposed for metastatic NSCLC therapy.