figur 5
dipol-dipol interaktioner i glycin. Akser er vist i de forskudte lp’-lp’ vinkler . Energiplotter af (a) Lennard-Jones 12-6 potentialer i det reviderede sæt steriske sammenstød; (b) alle elektrostatiske interaktioner; (c)-(f) de individuelle dipol-dipol-interaktioner af glycin-rygraden (se figur 1a for backbone skematisk af dipolerne). Energiparametre blev taget fra CHARMM22. De lyse områder viser regioner med minimal energi.
dipol-dipol-interaktioner i glycin
det reviderede glycin-steriske kort forklarer ikke den diagonale form af regionerne Karr, al, karrp, karrpr og Karr. I det generiske Ramachandran-plot blev det fundet, at regionernes diagonale form kunne reproduceres ved hjælp af elektrostatiske dipol-dipol-interaktioner, men kun når dipol-dipol-interaktionerne blev betragtet individuelt. Den samlede elektrostatiske interaktion gengiver ikke det observerede Ramachandran-plot . Her bruger vi den samme tilgang til behandling af individuelle elektrostatiske dipol-dipol-interaktioner langs glycin-rygraden.
Vi beregner energikortet for KRP-KRP for de 4 dipol-dipol-interaktioner i glycin-rygradsinteraktionen: COi-1···CO, NH * * * NHi + 1, CO * * * NH og COi-1···NHi + 1 (Figur 5c-F). De elektrostatiske interaktioner beregnes med Lennard-Jones-potentialerne i de steriske sammenstød identificeret i afsnittet ovenfor. Vi finder ud af, at formerne for de forskellige regioner i glycine Ramachandran-plottet (figur 3A) gengives (figur 5). Co * * * NH-interaktionen producerer den diagonale al -, Kristian-og Kristian-region (figur 5E). NH * * * Nhi+1-interaktionen producerer også en diagonal al-og kurr-region (figur 5D). Regionen er symmetrisk i forhold til aL-regionen. COi-1···CO-interaktion producerer minima svarende til de regioner, hvor der er tale om et Prip-og prippr-niveau (figur 5c).
i det originale glycin steriske kort (figur 2A) er regionen nær (kurr, Kurr) = (-180 Kurr, 180 Kurr) forbudt på grund af et sterisk sammenstød mellem O og H. Alligevel har glycin densitet i denne region i det observerede Ramachandran-plot (figur 3a). Dette kan også ses i frekvensfordelingen af d (O···H) (figur 3c), hvor der er en top ved d (O···H) ~ 2,4 liter. Ved denne top er O-og H-atomerne i kontakt, da VDD-diameteren er 2,5 liter. Således overvinder den gunstige Co···HN-dipol-dipol-interaktion i glycin-regionen af glycin den steriske afstødning af O-og H-atomerne (figur 5e).
Pre-proline Ramachandran-plottet
Schimmel og Flory argumenterede i 1968 for, at pre – prolin – aminosyrer forud for prolin-har et særligt begrænset Ramchandran-plot sammenlignet med det generiske Ramachandran-plot . Dette blev endelig observeret i proteindatabasen af MacArthur og Thornton (figur 1b) .
der er tre hovedforskelle mellem Pre-proline Ramachandran-plottet og det generiske Ramachandran-plot. I pre-proline Ramachandran-plottet er der en stor ekskluderet vandret strimmel ved -40 liter < liter < 50 liter, som begrænser Al-og prisregioner. Al-regionen flyttes højere op. Disse to funktioner blev gengivet i Schimmel-Flory beregning og efterfølgende beregninger . Den tredje funktion er et lille ben af tæthed, der stikker ud under Kristian-regionen (figur 1b; lilla i figur 2C). Karplus kaldte dette for regionen , som er unik for pre-proline.
tidligere beregninger fokuserede ikke på de individuelle interaktioner og redegjorde ikke for kurr-regionen. Her identificerer vi de nøjagtige steriske sammenstød, der bestemmer pre-proline Ramachandran plot. Vi vil derefter analysere de interaktioner, der er ansvarlige for regionen.
steriske interaktioner i præprolin-rygraden
i præprolin, i stedet for en interaktion med NH-atomet i den efterfølgende generiske aminosyre, interagerer præprolinen med en CH2-gruppe af den efterfølgende prolin (figur 1b). CH2-gruppen udøver en meget større sterisk effekt på præ-prolin Ramachandran-plottet. MacArthur og Thornton foreslog, at den dominerende effekt skyldes N···C-Kurri+1 og C-Kurri··· – Kurri+1-steriske sammenstød. Her kan vi analysere effektiviteten af hvert sammenstød ved at analysere de statistiske fordelinger direkte.
vi betragter de Koler-Koler afhængige interaktioner, der involverer C-Koler, h-Koler1 og H-koler2-atomer i den efterfølgende prolin (figur 1b). For hver interaktion genererer vi konturplottet i LP-lp af VDD-diameterafstanden. Ved at sammenligne konturplottet med den observerede tæthed i pre-proline Ramachandran-plottet identificerer vi de interaktioner, der fremkalder det bedste match i grænserne (figur 6a, interaktionerne identificeres i figur 2C). Vi fandt ud af, at den del, der blev taget ud af den nederste venstre kur-region af den observerede tæthed, skyldes Oi-1···C Kursi+1 sterisk sammenstød. En anden begrænsning på Al-og Kurt-regionerne skyldes H * * * C-Kurri+1 sterisk sammenstød.
vi overvejer derefter de afhængige interaktioner. I præ – prolin-frekvensfordelingen fandt vi tre forskellige toppe (nederste figur 6B). Den mest venstreorienterede top ved kr. ~ -50 kr. svarer til kr.-regionen i PR.prolin. Vi fokusere på de to toppe i β-regionen 50° << 180° større peak, der er centreret om ψ ~ 150° svarer til ßS region af generiske Ramachandran plot. I det generiske Ramachandran-plot er denne kartsregion afgrænset af C-Karr···O og C-Karr···Ni+1 steriske sammenstød. I pre-proline svarer den mindre Top centreret om karrus ~ 70 karrus til karrus-regionen og forekommer i en region, der ville blive udelukket af C Karrus···o sterisk sammenstød. I stedet er den mindre Top afgrænset nedenfra af N···C Kursi+1 sterisk sammenstød. Dette kan ses ved at sammenligne kurvefordelingen med modelkurven for N···C Kursi+1 vs. kursist (midt i figur 6b).
figur 6
pre-Proline parametre. (A) Ramachandran-plottet. De stiplede linjer viser de steriske sammenstød, der definerer nogle af grænserne for de observerede tætheder (se figur 2C). (B) fordelingen af forskellige inter-atomare interaktioner som en funktion af Kross. De stiplede linjer viser grænsen for VDD-diameterne. Den solide grå linje giver modelkurven beregnet med ideel geometri. I bunden er frekvensfordelingen af prisvinklen.
Ved hjælp af parametre fra CHARMM22 beregner vi Lennard-Jones 12-6 potentialet på grund af de reviderede steriske sammenstød (figur 7a). Lennard-Jones-potentialer kan ikke redegøre for den kurdiske region.
Figur 7
Energiplotter i pre-prolin som en funktion af Kurt-Kurt. Energiplotter af (a) Lennard-Jones 12-6 potentialer i det reviderede sæt steriske sammenstød; COi-1···C-pucker + 1-dipol-dipol-interaktioner, når den efterfølgende prolinring er i (b) op-pucker og (c) ned-pucker. De lyse områder viser regioner med lav energi.
interaktioner, der stabiliserer pre-prolin-regionen
da regionen (lilla i figur 2b) bringer interaktionen med C-regionen···O I sterisk konflikt, skal der være en kompenserende interaktion, der stabiliserer regionen. Hvad er denne interaktion? For at forstå denne interaktion overvejer vi en analogi med kurr-regionen i det generiske Ramachandran-plot. I regionen er en forvrænget Koi-1···HNI+1 hydrogenbinding dannes, hvilket bringer Hi + 1-atomet i kontakt med Oi-1-atomet. På samme måde er Oi-1-atomet i pre-prolin i pre-prolin-regionen i pre-prolin i kontakt med H-Lyr1-og H-Lyr2-atomerne (se figur 4b; tabel 1), hvilket antyder, at COi-1-gruppen interagerer med C-Lyrh Lyri+1-gruppen af den efterfølgende prolin.
tabel 1 parametre for co···brintbindingen
kan C-gruppen med Ch-1 interagere med COi-1? En sådan interaktion ville falde ind under klassen af CH···o svag hydrogenbinding, en veldokumenteret interaktion i proteiner . Undersøgelser af CH * * * o svag hydrogenbinding bruger et afstandskriterium på d (H···O) < 2,8 liter . Der er lidt vinkelafhængighed fundet i Ch * * * O-bindingen omkring H-atomet, hvor der ofte anvendes et vinkelkriterium for ch Ohh > 90-liter. Dette er meget mere tilladt end geometrien af den kanoniske hydrogenbinding. I tabel 1 opregner vi hydrogenbindingsparametrene for COi-1···C-Kurt-Kurt+1-interaktion i Kurt-regionen. Som proline kan påtage sig to forskellige store konformationer, op og ned pucker, målinger af geometrien af COi-1···C-Kurt+1 interaktion skal også opdeles i forhold til op og ned pucker. Den observerede geometri af COi-1···C-Kursth-geometri+1 opfylder de geometriske kriterier for den svage hydrogenbinding (tabel 1).
som COi-1···C Kurth Kursi + 1 svag hydrogenbinding er en tæt kontakt, vi er nødt til at modellere interaktionen for at forstå dens afhængighed af Kurt-kurvvinklerne. Til modelleringen, vi overvejer strategier, der er blevet brugt til den analoge COi-1···HNI + 1 hydrogenbinding. COi-1···HNI + 1 hydrogenbinding er blevet modelleret i kvantemekaniske undersøgelser, hvor det viste sig, at det var den mindste energikonformation i vakuum . En enklere tilgang, der modellerede hydrogenbindingen med elektrostatiske dipol-dipol-interaktioner, finder også et minimum i larr-regionen .
Her modellerer vi COi-1···C-L-L-L + 1 svag hydrogenbinding som en elektrostatisk dipol-dipol-interaktion (se metoder). Hvordan modellerer vi C-Kurt+1-gruppen som en elektrostatisk dipol? Bhattacharyya og Chakrabarti fandt ud af, at af CH-grupperne i prolin danner C-Kurt-gruppen de mest ch···o hydrogenbindinger. C-Karatatomet sidder ved siden af det elektronudtagende n-atom og er således mere surt end de andre C-atomer. Derfor placerer vi en lille negativ delvis ladning på C-Cholatom. I vores model finder vi et energiminimum i kurr-regionen for både op-rynken (figur 7B) og ned-rynken (figur 7c). Vi konkluderer, at COi-1···C-Kursi+1h-kursi1+1 svag hydrogenbinding stabiliserer kursisregionen i pre-prolin.