Waiting
Login processing...

Trial ends in Request Full Access Tell Your Colleague About Jove
JoVE Journal Biochemistry
Click here for the English version

Biochemistry

Dechiffrera de strukturella effekterna av att aktivera EGFR somatiska mutationer med molekylär dynamik Simulering

Published: May 20, 2020 doi: 10.3791/61125
Automatic Translation
1, 2, 2,3,4, 1
1Structural Bioinformatics Laboratory, Biochemistry, Faculty of Science and Engineering, Åbo Akademi University, 2Medicity Research Laboratories and Institute of Biomedicine, University of Turku, 3Turku Bioscience Centre, University of Turku and Åbo Akademi University, 4Department of Oncology and Radiotherapy, University of Turku and Turku University Hospital

Summary

Målet med detta protokoll är att använda molekylärdynamiksimuleringar för att undersöka de dynamiska strukturella förändringar som sker på grund av aktiverande mutationer av EGFR-kinasproteinet.

Abstract

Ett flertal somatiska mutationer som förekommer i den epidermala tillväxtfaktorreceptorn (EGFR) familj (ErbB) av receptortyrosin kinaser (RTK) har rapporterats från cancerpatienter, även om relativt få har testats och visat sig orsaka funktionella förändringar i ErbBs. ErbB-receptorerna dimeriseras och aktiveras vid ligandbindning, och dynamiska konformationella förändringar av receptorerna är inneboende för induktion av nedströmssignalering. För två mutationer som visas experimentellt för att ändra EGFR-funktion, A702V och Δ746ELREA750 deletionsmutation, illustrerar vi i följande protokoll hur simuleringar av molekylär dynamik (MD) kan sondera (1) konformationsstabiliteten hos den muterade tyrosinkinasstrukturen i jämförelse med EGFR av vild typ; 2) strukturella konsekvenser och konformationsövergångar och deras förhållande till observerade funktionella förändringar, (3) effekter av mutationer på styrkan av bindande ATP samt för bindning mellan kinasdomänerna i den aktiverade asymmetriska dimer; och (4) effekterna av mutationerna på viktiga interaktioner inom EGFR-bindningsstället som är associerat med det aktiverade enzymet. Protokollet ger en detaljerad steg-för-steg-procedur samt vägledning som kan vara mer allmänt användbar för undersökning av proteinstrukturer med hjälp av MD-simuleringar som ett sätt att sondera strukturell dynamik och förhållandet till biologisk funktion.

Introduction

Den mänskliga epidermala tillväxtfaktorreceptorn (EGFR) familj (ErbB) av receptortyrosinkinaser (RTKs) omfattar fyra medlemmar - EGFR/ErbB1/HER1, ErbB2/HER2, ErbB3/HER3 och ErbB4/HER4. ErbB-receptorerna reglerar fundamentala cellulära processer som celltillväxt och -proliferation, differentiering, migrationoch överlevnad 1,2, och är därmed potenta proto-onkogener. Avvikande aktivitet av ErbB-receptorer, särskilt EGFR och ErbB2, har ofta förknippats med mänskliga cancerformer gör ErbB-receptorer viktiga mål för cancer therapeutics2,3.

Flera somatiska förändringar av ERBB-gener har rapporterats från mänskliga maligniteter3,4,5. De bäst karakteriserade exemplen inkluderar de återkommande, aktiverande punktmutationerna och korta in-frame-deletioner i EGFR-kinasdomänen i icke-småcellig lungcancer (NSCLC). Dessa EGFR mutationer utgör viktiga drivkrafter för cancertillväxt, och förutsäga känslighet för EGFR inriktning cancerläkemedel6,7,8. I de flesta cancerformer förekommer dock somatiska mutationer i EGFR utanför dessa återkommande "hotspots" och fördelas över hela 1210-resterspannet i receptorn. Faktum är att de flesta av restprodukter längs EGFR primära sekvensen har visat sig vara muterade i mänskliga cancer9. Bortsett från de få hotspots, den funktionella betydelsen av den stora majoriteten av cancer-associerade EGFR mutationer förblir okänd.

Den monomeriska strukturen av ErbBs består av en stor amino terminal extracellulära domän, följt av en enda transmembrane helix leder till den intracellulära tyrosin Kinas domän och C-terminal svans region som innehåller dockningsplatser för intracellulära signalering proteiner. Ligandbindningen utlöser en dramatisk konformationsförändring i den extracellulära domänen, vilket underlättar bildningen av receptorsusponerare genom att exponera dimeriseringsarmarna som symmetriskt korsar varandra och interagerar med sina aromatiska/hydrofoba ytor. Vid receptor dimerbildning tyrosinkinasdomänerna kommer i kontakt asymmetriskt (Figur 1), vilket resulterar i aktivering av kinaserna som fosforylera C-terminala svansar av receptormonomerna, och därefter vid aktivering av nedströms signalering10,11.

Figure 1
Bild 1: EGFR-dimerns struktur. EGFR dimerizes när de extracellulära domänerna binda tillväxtfaktor (EGF, epidermal tillväxtfaktor). Mottagarkinasdomänen aktiveras sedan genom asymmetrisk interaktion med aktivatorkinasdomänen, och C-terminalsvansarna autofosforyleras vid tyrosinrester (Modifierad från Tamirat et al.12). Vänligen klicka här för att visa en större version av denna figur.

På grund av de dynamiska strukturella rearrangements som uppstår under monomeren dimer övergångar, tillsammans med kinas aktivering som är associerad med bildandet av en asymmetrisk dimer, kan mutationer längs den fulla längden av receptorstrukturen potentiellt ha en effekt på receptorfunktionen. Här beskriver vi flera exempel från våra tidigare studier där modellering av mutationen och visualisering var tillräckliga för att förklara konsekvenserna för funktion.

Exempel 1: En rapporterad mutation, D595V i ErbB413, ledde till ökad ErbB4 dimerisering och fosforylering14. Visualisering av platsen för mutationen var en kritisk faktor för att förstå de observerade funktionella effekter: D595V inträffade vid symmetrisk crossover av dimeric armarna i ektodomän (Figur 2A). Armarna är till stor del aromatiska och hydrofoba, och byte av polära asparaginsyra genom valine skulle förväntas öka den "klibbiga" hydrofoba interaktioner, stabilisera dimer och därmed öka den tid då fosforylering sker14. Det var en överraskning i början att hitta aspartat i varje arm, men i efterhand kan man tänka sig det som en timing mekanism för aktivitet, där polarsyran sidokedjor minska affinitet och livslängd intakt dimer och därmed begränsa kinas-medierad fosforylering och signalering. Byte av valine skulle då ta bort denna skyddsåtgärd genom att ytterligare stabilisera ErbB4 dimer.

Figure 2
Figur 2: Placering av en ErbB4 aktiverande mutation och mutationer som producerar kinasdöda ErbB4. (A) D595 (aktivera D595V mutation) är belägen på aromatiska / hydrofoba dimeric armarna i ErbB4 ektodomän modell; vapenmedarbetaren på tillväxtfaktorbindning, (närliggande rester visas som pinnar). (B) I ErbB4, G802 (inaktivera G802dup mutation) hjälper till att bilda bindningsfickan runt adenin ringen av ATP och katalytisk D861 (inaktivera D861Y mutationen) binder både Mg2 + (visas inte) och γ-fosfat grupp ATP. Vänligen klicka här för att visa en större version av denna figur.

Exempel 2: Man kan förutse att somatiska mutationer som riktar sig mot ATP-bindande platsen för kinasdomänen skulle förändra eller eliminera enzymatisk aktivitet som leder till en nedsatt eller kinasdöd receptor oförmögen att signalera. Av nio rapporterade mutationer från patienter med bröst, mag, kolorektal eller NSCLC15, två av de nio mutationer när testade hade mycket minskat fosforylering verksamhet16: G802dup (G → GG) och D861Y. Båda inaktivera somatiska mutationer hittades inom ATP bindande platsen för tyrosin Kinas domän struktur (Figur 2B): flexibel glycin, dupliceras, skulle förändra adenin ringen webbplats och små asparaginsyra ersättas med den skrymmande tyrosin nära terminalen fosfater skulle fysiskt hindra Mg2 +-ATP från bindande. Men sedan ErbB4 kan bilda en heterodimer med ErbB2 - ErbB2 binder inte en tillväxt dela upp i faktorer och beror på anslutning med en ErbB som gör för att heterodimerize - Erbbenen(aktivet),-ErbB4(kinase-dead) heterodimer skulle stimulera cellproliferation via Erk/Akting-signalvägen ännu celler skulle inte differentiera på grund av Kinas-döda ErbB4 och brist på STAT5 utbildningsavsnittaktivering 16.

I nyare studier blev det uppenbart att de dynamiska rörelserna hos ErbBs var relevanta för att förstå effekterna av vissa mutanter på ErbB-funktion, särskilt mutationer som förekommer inom domänen tyrosinkinas. Domänen för tyrosinkinas består av en N-lob (huvudsakligen β-ark) och C-lobe (till stor del alfa spiralformad), som är åtskilda av det katalytiska platsen där ATP binder. N-loben omfattar αC-helixen och P-slingan, medan aktiveringen (A-loop) och katalytiska slingor förekommer i C-loben17,18,19. Crystal strukturer av tyrosin kinas domänen visade två inaktiva konformationer, majoriteten av strukturer har Src-liknande inaktiva tillstånd. I den aktiva konformationen pekar A-slingans katalytiska aspartat mot ATP-bindningsstället och αC-helixen är orienterad mot ATP-bindningsfickan ("αC-in"-konformation), som bildar en stark glutamat-lysinjon-par-interaktion.

Eftersom ErbBs och komponenten kinas domän är mycket dynamiska enheter, och särskilt för fall där effekterna av mutationer på funktion och biologisk aktivitet sannolikt kommer att vara tätt kopplade till konformational tillstånd av ErbBs, är det viktigt att bedöma mutationer med avseende på det intervall av dynamiska förändringar de skulle uppleva. Röntgenkristallstrukturer av ErbBs ger statiska ögonblicksbilder av 3D-strukturen, som kan eller inte kan vara relevanta för att förstå de dynamiska konsekvenserna av en mutation. För att sondera omfånget av dynamiska förändringar som motsvarar "energilandskapet" som finns tillgängligt för en tredimensionell (3D) struktur används simuleringar av molekyldynamik (MD) i stor utsträckning20. När det gäller mutationer som skulle leda till lokala konformationella förändringar inom tyrosin kinas domän eller stabilisering av ett komplex, simuleringar på storleksordningen 100 ns kan vara tillräcklig. Dock kräver större skala konformationella förändringar (t.ex. övergångar mellan de aktiva och inaktiva konformationer av kinas domänen) en längre simuleringstid - på storleksordningen mikrosekonderna21.

Med avseende på det protokoll som beskrivs nedan vi överväga två aktiverande mutationer inom tyrosin kinas domän (Bild 3). Båda mutationer är belägna inom kinas domänen på platser som upplever lokala konformationella förändringar som dikterar om kinas är aktiv eller inte, och därmed MD simuleringar tillämpades i båda instanserna. I det första fallet anser vi förändringar som direkt påverkar ATP bindningsstället och katalytiska maskiner i EGFR mottagare Kinas domän, särskilt undersöka konsekvenserna av en exon 19 radering mutation som är allmänt inblandad i NSCLC4,7. Δ746ELREA750-mutationen, som minskar längden på β3-αC-slingan som föregår αC-helixen - helixen som rör sig mot bindnings/aktiv-platsen på kinasaktivering och deltar i att bilda den kritiska elektrostatiska interaktionen mellan E762 av helix och K745 genom att placera lysin för interaktion med ATP - predisponerar domänen föraktivering 12. I det andra fallet anser vi A702V mutation av EGFR, visat sig vara en ny gain-of-function aktiverande mutation avslöjas av iScream plattform9 och identifieras i en NSCLC patient22. Alanin-702 på mottagarkinasdomänen ligger på juxtamembrane segment B vid gränssnittet för mottagaren och aktivorkinasdomänerna, i vilka detta asymmetriska kinas dimerkomplex och kinaskonformationella förändringar krävs föraktivering 9.

Figure 3
Bild 3: Den asymmetriska kinasdomänen dimer av EGFR. A702V mutationen skulle vara belägen vid det kritiska gränssnittet för aktivor och mottagare Kinas domäner, intill αC helix och nära isoleucin 941 av aktivator kinas. Konformationsförändringar framkallade genom bildande av den asymmetriska dimern leda till kinasaktivering. β3-αC-slingan som innehåller ELREA-sekvensen direkt föregår αC-helixen; under aktiveringen rör sig αC-helixen inåt mot ATP-bindningsstället. Vänligen klicka här för att visa en större version av denna figur.

Protocol

OBS: De detaljerade steg som vidtagits för att undersöka effekterna av ΔELREA och A702V-mutationen på EGFR-strukturen med hjälp av MD-simuleringar diskuteras på följande sätt:

1. Strukturberedning

OBS: För att studera de strukturella effekterna av ΔELREA-mutationen bereds vildtyps- och mutantformer av APO-aktiva, ATP-bundna aktiva och apo inaktiva EGFR-monomerstrukturer enligt följande.

  1. Öppna visualiseringsprogrammet Chimera23 (https://www.cgl.ucsf.edu/chimera/) för att förbereda den aktiva EGFR-kinasstrukturen av vildtyp. I Arkiv-menyn klicka på alternativet Hämta efter ID och välj ProteinDatabankens (PDB24) databas och ange PDB-kod 2GS225 (2.8 Å resolution). PDB är ett förråd för 3D-strukturer som löses av olika experimentella tekniker, inklusive röntgenkristallografi, kärnmagnetisk resonansspektroskopi, kryo-elektronmikroskopi och neutrondiffraktion.
  2. Bygg de saknade strukturella elementen i 2GS2 genom att ta dessa segment från EGFR:s PDB-strukturer 1M1426 (2,6 Å) och 3W2S27 (1,9 Å). För att göra det, öppna 1M14 och 3W2S och överlagra dem på 2GS2 med hjälp av alternativet MatchMaker i jämförelsemenyn Verktyg → Struktur.
    1. Beskär ut de segment som ska läggas till från 1M14 och 3W2S. Välj slutatomer av restmaterialet före mellanrummen i 2GS2 och de atomer som skall läggas till från 1M14 och 3W2S (för detaljerad information om de tillagda segmenten, se tabell 1). På kommandoraden typ bond sel och hit enter. Den här strukturen är mallstrukturen.
  3. Använd mallstrukturen från steg 1.2 för att konstruera ΔELREA-borttagningsmutantformen av EGFR-kinasdomänen. Producera FAST-formatsekvensen av ΔELREA EGFR genom att spara sekvensen för mallstrukturen (Favorite → Sequence → File → save as) och sedan ta bort ELREA-sekvensen vid restnummer 746-750.
    1. Öppna ΔELREA EGFR-sekvensen i Chimera och anpassa den till sekvensen för 2GS2-mallstrukturen med hjälp av Sequence-menyn. I justeringsfönstret välj alternativet Struktur → Modeller (homology)28.
    2. På popup-fönstret ange 2GS2 sammansatta strukturen som mall och mutanten sekvensen som frågan som ska modelleras. Tryck sedan på OK. Välj en mutantmodell bland de resulterande modellerna, baserat på zDOPE-poängen (vanligtvis den lägsta poängen) och visuell inspektion.
  4. För att förbereda DEN ATP-bundna EGFR-kinasstrukturen av wild-typ, använd PDB-strukturen2ITX 29 (2,98 Å) som huvudstruktur. Bygga segment som saknas (se tabell 1) med hjälp av strukturerna 2GS625 (2.6 Å) och 3W2S som följer på förfarandet i steg 1.2. Konvertera liganden ANP i den resulterande strukturen till ATP genom att öppna PDB-filen i en textredigerare och ändra N3B kväveatom av ANP till en syreatom.
  5. Öppna strukturen i steg 1.4 i Chimera enligt vad som anges i steg 1.1. Lägg till en magnesiumjon till denna struktur från PDB-strukturen 2ITN29 (2,47 Å) för att uppnå en liknande positionering för Mg2+-jonen.
  6. Med den struktur som är resultatet av steg 1.5, modellera ΔELREA-mutantformen efter steg 1.3.
Apo aktiva EGFR Apo inaktiv EGFR ATP-bunden aktiv EGFR
Huvudstruktur 2GS2 2GS7 2ITX
Strukturer som används för att bygga saknade slingor 1M14 (723-725) 3W2S (958-984) 2GS6 (862-865)
3W2S (967-981) 4HJO (848-850) 3W2S (990-1001)

Tabell 1: Strukturer som används för att bygga sammansatta modeller av APO aktiva, apo inaktiva och ATP-bundna aktiva strukturer. Saknade regioner (aminosyran spänner i parentes) i det rektor strukturerar konstruerades från de listade strukturerar.

  1. För att förbereda den inaktiva EGFR-kinasstrukturen av vildtyp, öppna PDB-strukturen 2GS725 (2.6 Å) som i steg 1.1 och stryka de bundna liganderna och kristallografiska vattnen. Lägg till segmenten som saknas i 2GS7 (se tabell 1) från strukturerna 3W2S och 4HJO30 (2.75 Å) med hjälp av förfarandet i steg 1.2. Baserat på den slutliga inaktiva EGFR-strukturen, förbered mutantmodellen med hjälp av proceduren i steg 1.3.
    OBS: För A702V-mutationsstudien studeras EGFR-asymmetriska dimerstrukturen då mutationen är belägen vid den juxtamembrane B-segmentet av kinasdomänen som utgör en stor del av dimer-gränssnittet. De vilda och A702V mutanta EGFR-strukturerna bereds enligt följande:
  2. Den osymmetriska dimerstrukturen av vildtyp är konstruerad av PDB-strukturen 2GS2, som till en början visas i den monomeriska formen. För att konvertera till den biologiska monteringen som innehåller aktivatorn och mottagarkinaserna i det asymmetriska arrangemanget, öppna 2GS2 i Chimera som i (1.1) och utföra symmetriberäkningar genom att klicka på menyn Verktyg → Högre-Ordningsstruktur → Unit Cell. Markera 2GS2-strukturen och ange Gör kopior. Slutligen väljer du och sparar en enda asymmetrisk dimer från de flera kopior av dimer som resulterar från symmetriåtgärderna.
  3. Använda den osymmetriska EGFR-strukturen från steg 1.8, bygg A702V-mutanten genom att ersätta alanin 702 med valine med hjälp av alternativet Tools → Structure editing → Rotamers i Chimera.
    OBS: Sammantaget förbereds sex monomeriska och två dimeriska EGFR-strukturer för ΔELREA- respektive A702V-mutationsstudierna. Varje struktur bearbetas därefter för simulering med hjälp av guiden för beredning av proteiner i Maestro-programmet31 och MD-simuleringar görs med Amber-programmet32.
  4. Öppna strukturen i Maestro genom att använda alternativet → för importstruktur. Klicka sedan på knappen protein förberedelse guiden och välj följande: tillsätt väteatomer, bygga saknas sida-kedja atomer, bestämma protonation stater joniserbara rester vid pH 7.0 med hjälp av PROPKA, optimera orienteringen av asparagin, glutamin och histidin rester för väte bindning, och slutligen minimera strukturen.

2. Systeminställning

  1. Öppna språngprogrammet som ingår i amber-programpaketet. Importera ff14SB kraftfältet33 (källa leaprc.protein.ff14SB) och TIP3P vattenmolekyler34 (källa leaprc.water.tip3p). För de ATP-bundna systemen importera även parametrar för ATP35 (loadamberparams frcmod.phosphate, loadamberprep ATP.prep). Sedan, ladda strukturen (mol = loadpdb structure.pdb).
  2. Solvate strukturen i en oktaedral låda med explicit TIP3P vattenmolekyler som sträcker sig 10 Å i alla riktningar från ytatomer av proteinet (solvateoct mol TIP3PBOX 10,0).
  3. Kontrollera det byggda systemet (Kontrollera mol) och neutralisera det genom att lägga till nödvändiga joner (addions mol Na+ 0). För att tillräckligt modellera biomolekylära system, tillsätt ytterligare Na+/Cl- atomer till simuleringsrutan för att få systemets saltkoncentration till 0,15 M (addions mol Na+ X, addions mol Cl- X), där X ersätts med resultatet av: önskad saltkoncentration * antal vattenmolekyler * volym per vattenmolekyl * Avogadros antal.
  4. Generera och spara systemets topologi och koordinatfiler, som fungerar som ingångar för den efterföljande produktionssimuleringen (saveamberparm mol X.prmtop X.inpcrd).

3. Simulering av molekylär dynamik

  1. Använda Amber, initialt föremål simuleringssystemet till 5000 cykler av brantaste härkomst och konjugat gradient energi minimering att kringgå ogynnsamma konfigurationer. Genomför minimiseringen i flera steg, gradvis sänka fasthållningsanordningen som appliceras på soluteatomer från 25 kcal mol-1 Å-2 till 0 kcal mol-1 Å-2.
    1. I indatafilen för minimering, min.in, justera maxcyklen variabeln för den totala minimization cykeln (maxcyc = 5000) och ncyc att uppge antalet cykler för den brantaste nedstigning algoritm. Använd den restraint_wt variabeln för att tillämpa fasthållningskraften på soluteatomer som parametern restraintmask specificeras av. Kör sedan minimeringen enligt följande:
      $AMBERHOME/bin/sander -O -i min.in -o min.out -p X.prmtop -c X.inpcrd -r min.rst -ref X.inpcrd
      OBS: Strategin och de faktiska parametrar som används kan variera efter egna önskemål. Detaljer och vägledning finns från Amber manual och hemsida (https://ambermd.org/index.php)
  2. Värm upp systemet för 100 ps från 0 K till 300 K inställning en 10 kcal mol-1 Å-2 återhållsamhet på solute atomer. För att göra det, ställ in tempi = 0,0, temp0 = 300,0, dt = 0,002 ps, nstlim = 50000 och restraint_wt = 10 i heat.in-indatafilen. Utför uppvärmningen med följande kommando:
    $AMBERHOME/bin/sander -O -i heat.in -o heat.out -p X.prmtop -c min.rst -r heat.rst -x heat.mdcrd -ref min.rst
  3. Jämvikt systemet för 900 ps under en NPT ensemble; konstant antal atomer, temperatur (temp0 = 300,0) och tryck (ntp = 1), styra den med Berendsen-metoden (ntt = 1). Ställ in en 9 Å avstånd cutoff (cut = 9,0) för långväga elektrostatiska interaktioner. Sänk gradvis soluteatomens fasthållningsanordning till 0,1 kcal mol-1 Å-2 (restraint_wt = 0,1). Kör indatafilen för jämvikt equil.in beskriver ovanstående parametrar på följande sätt:
    $AMBERHOME/bin/sander -O -i equil.in -o equil.out -p X.prmtop -c heat.rst -r equil.rst -x equil.mdcrd -ref heat.rst
  4. Färdigställa jämvikten med en ohämmad 5 ns-simulering (set dt = 0.002 ps, ntslim = 2500000).
    $AMBERHOME/bin/sander -O -i equil_final.in -o equil_final.out -p X.prmtop -c equil.rst -r equil_final.rst -x equil_final.mdcrd -ref equil.rst
  5. Kontrollera att systemet har ekvilit genom att undersöka temperatur, tryck, densitet och energivärden.
    $AMBERHOME/bin/process_mdout.perl heat.out equil.out equil_final.out
    xmgrace sammanfattning. TEMP/DENSITET/ETOT/EPTOT/EKTOT
  6. Utför produktionssimuleringen för 100 ns (set dt = 0.002 ps, ntslim = 50000000 i prod.in) och spara konformationer var 10 ps (ntwx = 5000).
    $AMBERHOME/bin/sander -O -i prod.in -o prod.out -p X.prmtop -c equil_final.rst -r prod.rst -x prod.mdrcd -ref equil_final.rst

4. Analys

  1. Okulärbesiktning
    1. Visualisera konformationerna som samplas under simuleringarna av vilda och muterade EGFR-kinaser genom att öppna X.prmtop-filerna med bärnstensfärgad topologi och motsvarande prod.mdcrd-banfiler i VMD36. Med hjälp av bekväma sekundära strukturrepresentationer, analysera den övergripande strukturella dynamiken hos proteinerna från den registrerade banan. Visa specifika interaktioner mellan atomer/rester av intresse, såsom den katalytiskt väsentliga K745 - E762 saltbryggan.
    2. Alternativt, spara flera konformationer samplas under simuleringen i PDB-format och öppna dem med hjälp av Chimera programmet. Överlagra strukturerna på den inledande eller medianstrukturen med hjälp av alternativet MatchMaker. Visa den initiala/medianstrukturen i fast sken och resten av de justerade strukturerna i blekt vitt. Detta tillvägagångssätt gör det möjligt för en att visualisera de inspelade strukturella rörelser med mer klarhet.
      OBS: Förslag till effektiva representationer och bearbetning av konformationella ensembler från MDS finns i Melvin et al.37.
  2. RMSD- och RMSF-analys
    1. Compute rot-medelkvadratavvikelse (RMSD) och rot-genomsnittlig kvadrat fluctuation (RMSF) beräkningar med Cpptraj program38 att analysera den globala stabiliteten av proteinerna och undersöka flexibiliteten i de olika strukturella enheter. I rmsd.in- rmsf.in-indatafiler anger du stamnätsatomerna (för RMSD) och Cα-atomer (för RMSF) för den inledande strukturen som referens för RMS-montering. rmsf.in I filerna rmsd/rmsf.in importera de bärnstensfärgade topologifilerna (parm X.promtop) och motsvarande banfiler (trajin prod.mdcrd). Kör sedan kommandot Cpptraj -i rmsd/rmsf.in. Rita utdata för analys.
    2. Alternativt kan du anpassa konformationella ensembler och färga varje rester baserat på Cα atom RMSD. För att göra det, öppna konformationerna i Chimera och anpassa dem med Matchmaker alternativ.
      1. Gå till Verktyg → Avbildning → Rendera efter attribut. Välj Rester av den konformationella ensemblen och Cα RMSD som attribut och klicka på OK. Kedjespåret av konformationerna kommer då att färgas från blå → vita →, respektive reflekterande regioner av hög, medelhög och låg strukturell stabilitet.
  3. Analys av vätebindning
    1. Analysera vätebindningsinteraktionen mellan ATP och vildtyp/ΔELREA EGFRs. Förbered ett Cpptraj-skript, hbond.in, för att utföra denna uppgift. Definiera ett vätebindning med ett givare-acceptoravstånd på mindre än eller lika med 3,5 Å och en obligationsvinkel av större än eller lika med 135°. Ange analys endast för intermolekylära vätebindningar med variabeln nointramol d.v.s. vätebindningar mellan ATP och EGFR (hbond All nointramol dist 3.5 out nhb.agr avgout avghb.dat). Kör skriptet som Cpptraj -i hbond.in.
    2. Använd det här skriptet för att bedöma intramolekylära interaktioner, till exempel mellan restprodukterna K745 och E762, som är viktiga rester för EGFR-kinasaktivitet. För att göra det, specificera K745 som vätebindningsgivare och E762 som vätebindningsreceptoren i hbond.in och kör skriptet därefter.
  4. Övervakningsavstånd mellan atomer
    1. Mät avståndet mellan K745 och E762 genom att öppna banorna för vildtyp och ΔELREA apo EGFRs i VMD. Välj Cδ av Glu762 och Nz av Lys745 genom att klicka på Mus → etikett → obligation. Övervaka avståndet under simuleringen genom att rita en graf med → etiketter → bindnings→ graf.
  5. Gratis energiberäkningar
    1. För att beräkna de uppskattade bindande fria energierna mellan ATP och wild-type/ΔELREA EGFRs, och mellan aktivator och mottagarkinaser av VILD-typ/A702V EGFRs, använd den molekylära mekanik generaliserad Born surface area (MM-GBSA) modul39 som finns i AMBER-paketet. Ställ in ATP som ligand och EGFR som receptor i ΔELREA-studien. I A702V-studien anger du mottagarkinasen som liganden och aktivatornskinas som receptor.
    2. Först förbereda liganden, receptorn och ligand-receptor komplexa PDB filer separat i språng programmet inställning av PBRadii värdet till mbondi2. För PDB-filerna sparar du gasfasen bärnstenstopologi (.prmtop) och koordinatfiler (.inpcrd).
    3. Sedan, i mmgbsa-indatafilen, mmgbsa.in, ställ in igb = 2, saltcon = 0,1. Utför bindningsenergiberäkningar med hjälp av simuleringarnas banor, de förberedda receptorn/ligand-bärnstensfilerna och parametrarna i mmgbsa.in med MMPBSA.py-skriptet som finns i Amber enligt följande:
      $AMBERHOME/bin/MMPBSA.py -O -i mmgbsa.in -o mmgbsa.dat -sp X.prmtop -cp complex.prmtop -rp receptor.prmtop -lp ligand.prmtop -y prod.mdcrd -eo output.csv
    4. Analysera utdata, output.csv, genom att rita grafer.

Representative Results

Det beskrivna protokollet användes för att studera de strukturella effekterna av mutationerna ΔELREA och A702V på EGFR-kinasstrukturen. En tillämpning av protokollet var att undersöka effekten av mutationerna på lokala strukturella/conformational stabilitet genom att beräkna RMSD och RMSF värden från MD simuleringar. Eftersom A702V mutationen är belägen vid den juxtamemembrane B segmentet, rmsd av detta segment av mottagaren kinase i förhållande till startstrukturen beräknades för både vild-typ och A702V EGFRs. Resultatet (Figur 4A) visade att juxtramembrane B-segmentet av mutanten har ökat konformationsstabiliteten under 100 ns-simuleringen (genomsnittlig RMSD 0,7 Å - 95% konfidensintervall (CI) 0,009) jämfört med den wild-type EGFR-kinasdomänen (genomsnittlig RMSD 1.1 Å - 95% CI 0.01). Detta är mycket sannolikt resultatet av de snävare hydrofoba interaktioner vid dimer gränssnittet på grund av ersättning av alanin 702 (metylgruppen sidokedjan) till en skrymmande hydrofoba rester, valin (isopropyl grupp sidokedjan), vilket leder till ökad hydrofoba interaktioner av V702 på mottagaren kinas domän med isoleucin 941 av aktivator kinas domänen.

ΔELREA-mutationen är belägen vid β3-αC-slingan, intill den funktionellt kritiska αC-helixen; konformationen av αC-helixen är nyckeln till förskjutningar mellan egfr-kinas aktiva och inaktiva tillstånd. αC-helix konformationsstabiliteten i aktivt tillstånd bedömdes genom att undersöka RMSFEN över Cα-atomerna av rester inom helixen under MD-simuleringar (Figur 4B): övergripande, det finns lägre fluktuationer i mutanten (genomsnittlig RMSF 1.1 Å - 95% CI 0.4) i jämförelse med vildtypen (genomsnittlig RMSF 1.5 Å - 95% CI 0.57); med den högsta skillnaden i fluktuationer som registrerats för N-terminalens rester. De samplade konformationerna respektive superposerade på medianstrukturen för domänen med vildtypkinas och hos ΔELREA-kinasdomänen stöder också dessa resultat (Figur 4C): både vild- och ΔELREA-kinasdomänerna har övergripande liknande stabilitet för de superposerade konformationerna utom för β3-αC-slingan och αC-helixen, som är klart stabilare i ΔELREA EGFR. Dessa resultat visar att strykning av ELREA sekvensen hindrar förflyttning av den aktiva staten αC helix, därav hindra och därmed stabilisera den aktiva konformation. Vidare, eftersom αC-helixen utgör en del av det asymmetriska dimer-gränssnittet, skulle begränsningarna på αC-helixen i mutanten med stor tro stabilisera den asymmetriska dimer, vilket förlänger varaktigheten av det aktiverade tillståndet.

En annan tillämpning av protokollet är att undersöka beteendet hos viktiga intra- och intermolekylära interaktioner som äger rum under simuleringen. Således kan samspelet mellan K745 och E762, som är grundläggande för EGFR-enzymatisk aktivitet, analyserades för både den aktiva formen vildtyp och ΔELREA EGFR-kinas genom att mäta den procentuella beläggningen av vätebindningar som bildades mellan de två resternas polaratomer i sidokätting under MD-simuleringarna (Figur 5A): denna elektrostatiska interaktion bildades oftare i ΔELREA-kinasdomänen i jämförelse med domänen med vildtypskinas, på grund av den mer stabila αC-hecken som rymmer E762. Samspelet mellan Mg2+-ATP och domänerna för vildtyp och ΔELREA EGFR-kinas (Figur 5B) under simuleringen bedömdes också (Figur 5C): antalet vätebindningar var större för 5ΔELREA (genomsnittsvärde 4,0 - 95 % KI 0,03) än för EGFR av vildtyp (genomsnittsvärde 3,2 - 95 % KI 0,04). Ytterligare analys av vätebindningarna visade att K745 interagerar oftare med fosfatgrupperna av ATP i ELΔREA EGFR, som är kopplad till den stabilare K745-E762-interaktionen som noterats vid simuleringen av ΔELREA-mutanta EGFR-kinasdomänen.

MD simuleringar som beskrivs i protokollet är också användbara för att bedöma den relativa fria energin av bindning för protein-protein och protein-ligand interaktioner. Bindningsenergierna mellan aktivator- och mottagarkinasdomänerna av vildtyps- och A702V-EGFRs, och mellan ATP och domänerna för vilda och ΔELREA-mutanta EGFR-kinaser, var beräknat från molekylära mekaniker generaliserade Födda ytbehandlar område (MMGBSA) beräkningar (Figurera 6A): A702V-mutanten producerade en lägre genomsnittlig ΔGröra värderar (genomsnitt ΔGröra = -76 kcal/mol - 95% CI 0,47), som representerar mer gynnsamma dimer interaktioner, i motsats till den wild-type EGFR domän (genomsnitt ΔGbinda = -61 kcal/mol - 95% CI 0,61). Denna observation är förenlig med det mer stabila juxtamembrane B-segmentet och snävare dimer-gränssnitt på grund av ökade hydrofoba interaktioner som observerats för A702V EGFR-kinasdomänen. När det gäller ATP-bindning till domänerna för vildtyp och ΔELREA EGFR-kinaser (figur 6B) förutsäger MMGBSA-beräkningar starkare ATP-bindning med ΔELREA-mutanten (genomsnittlig ΔG-bindning -57 kcal/mol - 95% CI 0,43) i jämförelse med EGFR av vildtyp (genomsnittlig ΔG-bindning -48 kcal/mol - 95% CI 0,33).bind bind Detta utfall ligger i linje med det större antal vätebindningar som registrerats mellan ATP och ΔELREA EGFR (figur 5C) i jämförelse med domänen av vildtyp.

Protokollet kan också användas för att undersöka konformationsförändringar som observerats under en simulering. I den aktuella studien studerades effekterna av ΔELREA-mutationen på den inaktiva EGFR-konformationen genom visuell inspektion och superpositionering av de samplade konformationerna från simuleringen. Analysen avslöjade en inåtgående rörelse av αC-helixen i ΔELREA EGFR-kinasdomänen (Figur 7A), en strukturförändring som förväntas under övergången till det aktiva tillståndet. Däremot behöll αC-spiralen av inaktiva EGFR av vildtyp sin ursprungliga konformation (Figur 7B). Således MD simuleringar stödja förslaget att borttagning mutationen, visas experimentellt att öka Kinas verksamhet40,41, främjar en konformationell förskjutning från inaktiva kinas mot det aktiva staten.

Figure 4
Figur 4: Vildtyps- och mutant konformationsstabilitet för den aktiva EGFR-kinasdomänen under MD-simuleringar. (A) RMSD (ryggradsatomer) över den juxtamembrane B segmentet av wild-type (blå) och A702V (röd) mottagare kinas domän. (B) RMSF (Cα atomer) över rester av αC-helixen: vildtyp (blå) och ELΔREA (guld). (C) Superposed samplad konformationer av vildtyp (vänster) och ΔELREA (höger) EGFR-kinasdomän; kedjespår som är färgade baserat på RMSD (Cα atomer) av varje rest med avseende på medianstrukturen. Färgning varierar från blått till vitt till rött, som representerar regioner med hög till låg konformationsstabilitet. Observera att de "fria" N-terminalregionerna i den isolerade kinasdomänen, färgad röd, inte skulle uppvisa denna nivå av rörlighet i den intakta EGFR-strukturen. Siffror anpassade från Chakroborty et al.9 (Figur 4A återges med tillstånd av Journal of Biological Chemistry) och Tamirat et al.12. Vänligen klicka här för att visa en större version av denna figur.

Figure 5
Figur 5: Nyckelfunktioner som ses i den aktiva mottagarkinasen under MD-simuleringar: saltbryggan K745-E762, αC-helixen och interaktioner med ATP. (A) Procentuell beläggning av växelverkan K745-E762 under simuleringen av domänerna för vildtyp (blå) och ΔELREA (guld) EGFR-kinas. (B) Rester av vild-typ och ΔELREA mutant interagerar med ATP (pinnar). Mg2+ (grön) koordinater med ATP och D855. (C) Antalet vätebindningar som bildas av ATP med både domänerna för ΔELREA EGFR-kinas under MD-simuleringar. Figur från Tamirat et al.12. Vänligen klicka här för att visa en större version av denna figur.

Figure 6
Bild 6: Lägre relativa fria energier av bindning observeras för de muterade kinasdomänerna under simuleringarna. (A) Bindningsenergier beräknade för interaktionen mellan aktivator- och mottagarkinasdomänerna av vildtypsdomänerna (blå) och A702V (röd) EGFRs. (B) ΔG-bindning av ATP till wild-type (blå) och ΔELREA (guld) EGFR-kinasdomäner.bind Siffror anpassade från Chakroborty et al.9 (Figur 6A återges med tillstånd av Journal of Biological Chemistry) och Tamirat et al.12. Vänligen klicka här för att visa en större version av denna figur.

Figure 7
Figur 7: Superposerade konformationer från domänen med vildtyp och ΔELREA inaktiv EGFR-kinas. Konformation av αC-helix- och A-loopspiralen av (A) vild-typ (medianstruktur i blått) och (B) ΔELREA EGFRs (guld). Andra samplade konformationer, blekt vitt; initiala strukturer före MD-simuleringarna, rosa. Figur från Tamirat et al.12. Vänligen klicka här för att visa en större version av denna figur.

Discussion

Det protokoll som beskrivs i denna studie fokuserar på att använda molekylär dynamik simuleringar för att undersöka lokala och globala strukturella förändringar som uppstår från aktivera somatiska mutationer av EGFR kinas domän. Även om röntgenkristallstrukturer av vildtyp och mutanta EGFR ger ovärderlig strukturell insikt, skildrar de en eller några statiska framställningar. Inneboende till den biologiska funktionen av ErbBs är dock de nödvändiga övergångarna mellan den enzymatisk inaktiva och aktiva tyrosin kinas, åberopa dynamiska förändringar i både struktur och intramolekylära interaktioner mellan kinas monomerer. MD-simuleringar utfördes således för att undersöka den dynamiska karaktären hos EGFR tyrosinkinasdomänen, inklusive vildtypsstrukturen, den introducerade ΔELREA-deletionsmutationen och A702V-mutationen. Dessa simuleringar var framgångsrika i att klargöra den troliga rollen av dessa mutationer i strukturerna och hur deras effekter på konformation av tyrosin kinas domän skulle leda till experimentellt observerade ökningar i EGFR kinas verksamhet.

Ett avgörande steg i detta protokoll är användningen av en relevant struktur för att bedöma effekten av mutationen. Ett sätt att välja en relevant simuleringsingångsstruktur är att visualisera placeringen av mutationen i den statiska 3D-strukturen och undersöka dess möjliga påverkan med avseende på de angränsande aminosyrorna och strukturenheterna. I denna studie, till exempel, eftersom A702V EGFR mutationen ligger vid juxtamembrane B segmentet som bildar den asymmetriska dimer gränssnitt, användningen av dimer struktur för simuleringen i motsats till monomer är kritisk. Användningen av en monomerisk struktur skulle ha utsatt juxtamembrane B-segmentet av mottagaren kinase till lösningsmedlet, beröva den från den stabiliserande interaktioner, förstärkt av mutationen till en större hydrofoba rester och interaktioner med isoleucin 941 från C-lob resterna av aktivator kinas. Det är dessutom anmärkningsvärt att den 3D-struktur som koordinaterna representerar i en fil för preliminära budgetöverföringar inte nödvändigtvis motsvarar den biologiskt relevanta struktur som bör användas för studier. Till exempel, med strukturen av ErbB4, PDB-kod 3BCE, motsvarar PDB-koordinaterna en trimer, men detta beror på kristallkontakter (få kontakter mellan monomererna ses vid visualisering av denna struktur). Matriser inom PDB-filen kan användas (t.ex. inom Chimera) för att rekonstruera de kristallografiskt relaterade strukturerna, som kan visualiseras för att identifiera kedjor som motsvarar den biologiskt relevanta 3D-strukturen enligt vad som redovisas i den ursprungligapublikationen 42. Ett annat viktigt steg i protokollet är att korrekt förbereda simuleringen inmatningsstrukturen, såsom att bygga saknade aminosyror i olika sling regioner, och särskilt där ligger i närheten av mutationen. Även om ett flertal EGFR-strukturer av vildtyp finns i det preliminära budgettalet, finns endast ett begränsat antal muterade EGFR-strukturer tillgängliga. Följaktligen måste de muterade strukturerna också modelleras; för en enda restmutation som A702V användes Chimera för att mutera resthalten; medan, för ΔELREA-deletionsmutationen, Modeller användes.

De olika parametrar som utnyttjas i simuleringens indatafiler - till exempel antalet minimeringscykler, uppvärmning av systemet till önskad temperatur på en gång eller istället värmer långsamt genom flera mellanliggande temperaturer, tidsperioden för jämvikten och för produktionssimuleringarna - kan modifieras utifrån studiemolekylen, syftet med arbetet och ens egna preferenser. Samtidigt som man utför MD-simuleringar, är det också vanligt att komma över fel som kan uppstå från indatafilerna, frågor som rör simuleringsprogramvaran i bruk eller till och med ett användarfel. Därav, Det är mycket viktigt att förstå källan till felen genom att noggrant undersöka eventuella felmeddelanden. De flesta simuleringsprogram har en e-postlista där användarna kan ställa frågor till programutvecklarna och till andra användare genom vilka de flesta problem kan lösas. Dessutom ger användarhandböcker betydande hjälp för att förstå detaljerna i simuleringsprotokollet, inklusive antaganden och begränsningar. Även om MD-simulering är ett viktigt verktyg för att utforska molekylers dynamiska egenskaper, kom ihåg att beräkningsresultaten måste utvärderas noggrant tillsammans med andra informationskällor för att bedöma deras giltighet. När så är möjligt, arbeta tillsammans med forskare som är experter på de proteiner som studeras, särskilt där relevanta våt-lab experimentella studier görs, som tjänar till att ge resultat för strukturell tolkning samt att föreslå experiment som kan göras baserat på strukturella observationer för att testa hypoteser.

I denna studie var protokollet effektivt när det gäller att undersöka de dynamiska strukturella effekterna av mutationerna ΔELREA och A702V på EGFR-kinasstrukturerna. Simuleringarna visade att ΔELREA håller fast den funktionellt väsentliga αC-helixen och främjar en konformationsförskjutning från den inaktiva kinasen till en stabiliserad aktiv kinas. Simuleringsresultaten stöds oberoende av läkemedelsresponsdata som visade effekterna av tyrosinkinashämmare på lungcancercelllinjer med ΔELREA-deletionsmutationen och EGFR av vildtyp, där större hämning av läkemedel som erkänner den aktiva kinasreformationen rapporterades för ΔELREA än för EGFR12 av vildtyp. Med A702V-mutationen indikerar MD-simuleringarna, i jämförelse med vildtyp, ökad stabilisering av activator-mottagare kinase-gränssnittet samt högre affinitet av aktivatorn och mottagarenskinn för varandra, tillsammans stödjande underhåll av den aktiverade konformationen av EGFR-kinasen. A702V mutationen, som ligger på den juxtamembrane B-segmentet av mottagaren kinase, skulle öka hydrofoba interaktioner med aktivator kinas, fungerar för att förlänga varaktigheten av det aktiverade tillståndet. A702V-mutationen stöder cellernas överlevnad i avsaknad av tillväxtfaktor och identifierades i en in vitro-screening för EGFR-mutationer9.

Disclosures

Författarna har inget att avslöja.

Acknowledgments

Denna forskning finansieras genom bidrag till M.S.J från Finlands Akademi (308317, 320005), Sigrid Juselius Stiftelse och Tor, Joe och Pentti Borg minnesfond, samt till K.E. från Finlands Akademi (274728, 316796), Stiftelsen Cancer i Finland, och Åbo Universitetscentralsjukhus. M.Z.T. finansieras av Åbo Akademis doktorandnätverk för informations- och strukturbiologi. Vi tackar CSC IT Center for Science för datorresurserna och Dr. Jukka Lehtonen för IT-stödet inom biocenter Finlands bioinformatiknätverk; och Biocenter Finland strukturbiologi infrastruktur nätverk.

Materials

Name Company Catalog Number Comments
Amber software University of California, San Francisco Version 2018 Executable
Chimera program Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco Version 1.13.1 Executable
EGFR struture files The Protein Data Bank 3D coordinates of EGFR structures
Maestro Schrödinger LLC Version 2018-3 Executable
Modeller program The Andrej Šali Lab, Departments of Biopharmaceutical Sciences and Pharmaceutical Chemistry, University of California San Francisco Included in the Chimera program
VMD software Theoretical and Computational Biophysics Group, University of Illinois at Urbana-Champaign Version 1.9.3 Executable

DOWNLOAD MATERIALS LIST

References

  1. Yarden, Y., Sliwkowski, M. X. Untangling the ErbB signalling network. Nature Reviews Molecular Cell Biology. 2, 127-137 (2001).
  2. Lemmon, M. A., Schlessinger, J., Ferguson, K. M. The EGFR family: not so prototypical receptor tyrosine kinases. Cold Spring Harbor Perspectives in Biology. 6, a020768 (2014).
  3. Arteaga, C. L., Engelman, J. A. ERBB receptors: From oncogene discovery to basic science to mechanism-based cancer therapeutics. Cancer Cell. 2, 282-303 (2014).
  4. Mishra, R., Hanker, A. B., Garrett, J. T. Genomic alterations of ERBB receptors in cancer: Clinical implications. Oncotarget. 8, 114371-114392 (2017).
  5. Cerami, E., et al. cBioPortal for Cancer Genomics. , https://www.cbioportal.org (2020).
  6. Lynch, T. J., et al. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. New England Journal of Medicine. 350 (21), 2129-2139 (2004).
  7. Paez, J. G., et al. EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 304 (5676), 1497-1500 (2004).
  8. Pao, W., et al. EGF receptor gene mutations are common in lung cancers from "never smokers" and are associated with sensitivity of tumors to gefitinib and erlotinib. Proceedings of the National Academy of Sciences U.S.A. 101 (36), 13306-13311 (2004).
  9. Chakroborty, D., et al. Unbiased in vitro screen for activating EGFR mutations. Journal of Biological Chemistry. 294 (24), 9377-9389 (2019).
  10. Leahy, D. J. Structure and Function of the Epidermal Growth Factor (EGF/ErbB) Family of Receptors. Advances in Protein Chemistry. 68, 1-27 (2004).
  11. Roskoski, R. ErbB/HER protein-tyrosine kinases: Structures and small molecule inhibitors. Pharmacological Research. 87, 42-59 (2014).
  12. Tamirat, M. Z., Koivu, M., Elenius, K., Johnson, M. S. Structural characterization of EGFR exon 19 deletion mutation using molecular dynamics simulation. PLoS ONE. 14 (9), e0222814 (2019).
  13. Ding, L., et al. Somatic mutations affect key pathways in lung adenocarcinoma. Nature. 455, 1069-1075 (2008).
  14. Kurppa, K. J., Denessiouk, K., Johnson, M. S., Elenius, K. Activating somatic ERBB4 mutations in non small-cell lung cancer. Oncogene. 35 (10), 1283-1291 (2016).
  15. Soung, Y. H., et al. Somatic mutations of the ERBB4 kinase domain in human cancers. International Journal of Cancer. 118, 1426-1429 (2006).
  16. Tvorogov, D., et al. Somatic mutations of ERBB4: selective loss-of-function phenotype affecting signal transduction pathways in cancer. Journal of Biological Chemistry. 284, 5582-5591 (2009).
  17. Hubbard, S. R., Till, J. H. Protein tyrosine kinase structure and function. Annual Review of Biochemistry. 69 (1), 373-398 (2000).
  18. Huse, M., Kuriyan, J. The conformational plasticity of protein kinases. Cell. 109 (3), 275-282 (2002).
  19. Jura, N., et al. Catalytic control in the EGF receptor and its connection to general kinase regulatory mechanisms. Molecular Cell. 42, 9-22 (2011).
  20. Karplus, M., Kuriyan, M., J, Molecular dynamics and protein function. Proceedings of the National Academy of Sciences U.S.A. 102 (19), 6679-6685 (2005).
  21. Shan, Y., Arkhipov, A., Kim, E. T., Pan, A. C., Shaw, D. E. Transitions to catalytically inactive conformations in EGFR kinase. Proceedings of the National Academy of Sciences U.S.A. 110 (18), 7270-7275 (2013).
  22. Reckamp, K. L., et al. A phase I trial to determine the optimal biological dose of celecoxib when combined with erlotinib in advanced non-small cell lung cancer. Clinical Cancer Research. 12 (11 Pt 1), 3381-3388 (2006).
  23. Pettersen, E. F., et al. UCSF Chimera-a visualization system for exploratory research and analysis. Journal of Computational Chemistry. 25 (13), 1605-1612 (2004).
  24. Berman, H. M., et al. The Protein Data Bank. Nucleic Acids Research. 28 (1), 235-242 (2000).
  25. Zhang, X., Gureasko, J., Shen, K., Cole, P. A., Kuriyan, J. An Allosteric Mechanism for Activation of the Kinase Domain of Epidermal Growth Factor Receptor. Cell. 125 (6), 1137-1149 (2006).
  26. Stamos, J., Sliwkowski, M. X., Eigenbrot, C. Structure of the epidermal growth factor receptor kinase domain alone and in complex with a 4-anilinoquinazoline inhibitor. Journal of Biological Chemistry. 277 (48), 46265-46272 (2002).
  27. Sogabe, S., et al. Structure-Based Approach for the Discovery of Pyrrolo[3,2-d]pyrimidine-Based EGFR T790M/L858R Mutant Inhibitors. ACS Medicinal Chemistry Letters. 4 (2), 201-205 (2013).
  28. Sali, A., Blundell, T. L. Comparative protein modelling by satisfaction of spatial restraints. Journal of Molecular Biology. 234 (3), 779-815 (1993).
  29. Yun, C. H., et al. Structures of lung cancer-derived EGFR mutants and inhibitor complexes: mechanism of activation and insights into differential inhibitor sensitivity. Cancer Cell. 11 (3), 217-227 (2007).
  30. Park, J. H., Liu, Y., Lemmon, M. A., Radhakrishnan, R. Erlotinib binds both inactive and active conformations of the EGFR tyrosine kinase domain. Biochemical Journal. 448 (3), 417-423 (2012).
  31. Schrödinger LLC. Release 2018-3: Maestro. , https://www.schrodinger.com/maestro (2018).
  32. Case, D. A., et al. AMBER 2018. , University of California. San Francisco. (2018).
  33. Maier, J. A., Martinez, C., Kasavajhala, K., Wickstrom, L., Hauser, K. E., Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. Journal of Chemical Theory and Computation. 11 (8), 3696-3713 (2015).
  34. Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., Klein, M. L. Comparison of simple potential functions for simulating liquid water. Journal of Chemical Physics. 79 (2), 926-935 (1983).
  35. Meagher, K. L., Redman, L. T., Carlson, H. A. Development of polyphosphate parameters for use with the AMBER force field. Journal of Computational Chemistry. 24 (9), 1016-1025 (2003).
  36. Humphrey, W., Dalke, A., Schulten, K. VMD: Visual molecular dynamics. Journal of Molecular Graphics. 14 (1), 33-38 (1996).
  37. Melvin, R. L., Salsbury, F. R. Visualizing ensembles in structural biology. Journal of Molecular Graphics and Modelling. 67, 44-53 (2016).
  38. Roe, D. R., Cheatham, T. E. PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data. Journal of Chemical Theory and Computation. 9 (7), 3084-3095 (2013).
  39. Miller, B. R., et al. MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. Journal of Chemical Theory and Computation. 8 (9), 3314-3321 (2012).
  40. Guha, U., et al. Comparisons of tyrosine phosphorylated proteins in cells expressing lung cancer-specific alleles of EGFR and KRAS. Proceedings of the National Academy of Sciences U.S.A. 105 (37), 14112-14117 (2008).
  41. Furuyama, K., et al. Sensitivity and kinase activity of epidermal growth factor receptor (EGFR) exon 19 and others to EGFR-tyrosine kinase inhibitors. Cancer Science. 104 (5), 584-589 (2013).
  42. Qiu, C., et al. Mechanism of Activation and Inhibition of the HER4/ErbB4 Kinase. Structure. 6 (3), 460-467 (2008).

Tags

Denna månad i JoVE Molekylär dynamik simuleringar Proteinstruktur analyser Epidermal tillväxtfaktor receptor familj Somatiska mutationer Cancer
Dechiffrera de strukturella effekterna av att aktivera EGFR somatiska mutationer med molekylär dynamik Simulering
Play Video
PDF DOI DOWNLOAD MATERIALS LIST

Cite this Article

Tamirat, M. Z., Kurppa, K. J.,… More

Tamirat, M. Z., Kurppa, K. J., Elenius, K., Johnson, M. S. Deciphering the Structural Effects of Activating EGFR Somatic Mutations with Molecular Dynamics Simulation. J. Vis. Exp. (159), e61125, doi:10.3791/61125 (2020).

Less
Copy Citation Download Citation Reprints and Permissions
View Video

Get cutting-edge science videos from JoVE sent straight to your inbox every month.

Waiting X
Simple Hit Counter