Wpływ szczelności węzłów na poziomie molekularnym

Znaczenie

węzły i splątania występują w białkach, DNA i polimerach syntetycznych i są wykorzystywane do tworzenia podstawy splecionych nanomateriałów. Zrozumienie, w jaki sposób napięcie splątań molekularnych wpływa na właściwości, ma kluczowe znaczenie dla przyszłego projektowania wiązanych, tkanych i splątanych cząsteczek i materiałów. Jednak podczas gdy istnieje wiele badań teoretycznych na temat takich systemów, jest bardzo niewiele badań eksperymentalnych. Badamy tutaj wpływ szczelności węzłów w zakresie właściwości fizycznych i chemicznych dla trzech węzłów wiązanych w nici o pętli zamkniętej 20, 23 i 26 nm. Stwierdzamy, że szczelność węzłów molekularnych znacząco wpływa na reaktywność, konformację i ekspresję chiralności.

Streszczenie

trzy 819 węzłów w nici pętli zamkniętej o różnych długościach (∼20, 23 i 26 nm) wykorzystano do eksperymentalnej oceny konsekwencji szczelności węzłów na poziomie molekularnym. Dzięki zastosowaniu 1H NMR, spektroskopii uporządkowanej dyfuzyjnie (DOSY), dichroizmu kołowego (CD), spektrometrii mas dysocjacji wywołanej kolizją (Cid-MS) i symulacji dynamiki molekularnej (MD) na węzłach o różnej wielkości, stwierdzamy, że na strukturę, dynamikę i reaktywność łańcuchów molekularnych dramatycznie wpływa szczelność wiązania. Napięcie splątania powoduje różnice konformacji, wzmacnia ekspresję topologicznej chiralności, osłabia wiązania kowalencyjne, hamuje zdarzenia dekompleksji i zmienia właściwości absorpcyjne. Zrozumienie skutków dokręcania nanoskali węzłów może z pożytkiem informować o projektowaniu wiązanych i splątanych materiałów molekularnych.

  • węzły molekularne
  • Chemia supramolekularna
  • topologia chemiczna

węzły znajdują się w niektórych białkach (1), liniowym i kołowym DNA (2) oraz polimerach o wystarczającej długości i elastyczności (3). Splątania nici wpływają na rozmiar cząsteczkowy (4), stabilność (5) i różne właściwości mechaniczne (6⇓-8), chociaż wiele ze zrozumienia, w jaki sposób i dlaczego pozostaje niejasne. Do tej pory wpływ różnych cech strukturalnych, takich jak liczba skrzyżowań węzłów (9), wić (10), stosunek krzyżowania szkieletu (BCR) (11) i globalny promień krzywizny (12), na właściwości był głównie badany przez symulacje (13 ⇓ ⇓ -16), a nie eksperyment (17, 18). Szczelność węzłów (19⇓-21) jest szczególnie łatwą do docenienia cechą, znaną z naszych codziennych doświadczeń w makroskopowym świecie, która może mieć również znaczący wpływ na poziomie molekularnym (4, 22⇓⇓⇓⇓-27). Przewiduje się, że długie polimery mają obszary wiązane, które są zarówno częste, jak i ciasne (25). Szczelność wiązania jest związana z różnicami w termostabilności splątanych białek (26) i uważa się, że ma to wpływ na wytrzymałość na rozciąganie (27). Trudno jednak ocenić wewnętrzne efekty wiązania, porównując struktury o raczej odmiennym składzie chemicznym. Syntetyczne węzły molekularne monodyzowane są idealnymi modelami, za pomocą których można ocenić wpływ dokręcania węzłów na właściwości fizyczne i chemiczne (17, 18). Tutaj zgłaszamy doświadczalnie określone różnice własności w zestawie trzech wiązanych cząsteczek, które różnią się tylko długością elastycznych regionów (łańcuchów alkilowych), które oddzielają bardziej sztywne sekcje (złożone z pierścieni aromatycznych) nici. Stosując strategię oplatania wcześniej używaną do złożenia (28) niezwykle ciasno zawiązanego węzła pętli 819 z 192 atomami (1) (29), łańcuchy zakończone alkenami użyte do zamknięcia wiązanej struktury zostały przedłużone bez zmiany w inny sposób wyniku syntezy węzła (Fig. 1). W wyniku tego powstał zestaw trzech 819 węzłów molekularnych, które różnią się jedynie długością łańcuchów alkilowych w pętli (1 ma długość 192 atomów, 2 ma długość 216 atomów, a 3 ma długość 240 atomów, co stanowi 25% różnicy w długości nici obejmującej tylko regiony elastyczne). Zbadaliśmy właściwości węzłów różnej wielkości za pomocą spektrometrii 1H NMR, spektrometrii masowej, spektroskopii UV i dichroizmu kołowego i wykorzystaliśmy badania obliczeniowe, aby wyjaśnić rolę rosnącej szczelności wiązania W zmianie zachowania fizycznego i chemicznego.

aby ocenić wpływ szczelności węzłów na siłę wiązania kowalencyjnego, przeprowadzono eksperymenty z tandemową spektrometrią mas (MS-MS) na piku 2+ (Fig. 4 A, C i E) węzłów 1-3 aktywowanych przez dysocjację wywołaną kolizją (Cid). Spowodowało to mnożenie naładowanych jonów z utratą masy zgodną z nacięciem wiązania kowalencyjnego i późniejszym rozwinięciem wiązania wiązanego (Fig. 4 B, D i F). Po zwiększeniu znormalizowanej energii zderzenia w krokach 1-eV do momentu zerwania wiązania, luźniejsze węzły (2 i 3) wymagały znacznie większych energii (odpowiednio 35 I 55 eV) do rozpadu niż 1 (28 EV). Dla każdego rozcięcia węzła zawsze występuje w wiązaniu węgiel-tlen między centralnymi jednostkami bipirydyny i fenylo-bipirydyny, ale węzły wykazywały dwa różne typy wzoru fragmentacji (typ I, Fig. 4G i typ II, rys. 4H). W przypadku najciaśniejszego węzła 1, po fragmentacji (utracie cząsteczki wody; m / z = 1,681.75) następuje utrata sąsiedniej centralnej jednostki bipirydyny (m/z = 1,582.67; Fig. 4B, typ I). Natomiast najbardziej luźny węzeł 3 wymaga co najmniej 55 eV do fragmentacji i bezpośrednio tworzy dwa produkty rozpadu, z których każdy zawiera połowę węzła (m/z = 2,008.36 i 1,810.45; Fig. 4F, typ II). Oba rodzaje wzoru fragmentacji obserwuje się z węzłem 2 przy 35 eV (rys. 4D, m/z = 1 849,75 i 1 750,83 dla typu I oraz m / z = 1 840,75 i 1 642,83 dla typu II). Wyraźnie szczelność węzłów znacząco wpływa na siłę wiązania kowalencyjnego; im mocniej węzeł, tym łatwiej można złamać niektóre wiązania. Wyniki eksperymentów są zgodne z symulacjami, które sugerują, że ciasne wiązanie wymusza przyjęcie konformacji o napiętych długościach i kątach wiązania.

4.

fragmentacja 819 węzłów 1-3 w ramach eksperymentów CID-MS. A) CID-MS węzła 1. B) MS-MS 2+ (m / z = 1,690.75) od węzła 1, obserwuje się tylko fragmentację typu I. C) CID-MS węzła 2. D) MS-MS 2+ (m / z = 1859.00) od węzła 2 obserwuje się oba rodzaje fragmentacji. E) CID-MS węzła 3. F) MS-MS 2+ (m/z = 2,027.09) węzła 3, obserwuje się tylko fragmentację typu II. Nie zaobserwowano żadnych jonów fragmentarycznych odpowiadających makrocykliom, które mogłyby powstać z ogniw molekularnych (np. ogniwa katenanowego lub Solomona) przy fragmentacji jednego pierścienia (17) w zakresie energii zderzeń, co dodatkowo dowodzi, że obecne są tylko węzły. (G I H) dwie różne drogi fragmentacji (typ i, G; Typ II, H) dla prędkości 819 węzłów, odpowiadające wynikom MS-MS pokazanym w punktach B, D i E; dekolt wiązania jest oznaczony czerwoną falistą linią, pozostałości utracone w pierwszym kroku są zabarwione na fioletowo, a te utracone w drugim kroku są zabarwione na niebiesko.

Analiza indeksów wiązań Wiberga, zaimplementowana w pakiecie programów Gaussian09 (35) została przeprowadzona za pomocą funkcji Meta-korelacji hybrydowej Minnesota 2006 (M06-2X) (36) i zestawu baz Pople ’ A Z FUNKCJAMI polaryzacyjnymi, 6-31g(d), na globalnych minimalnych strukturach energetycznych podwójnie naładowanych 1 i 3 w celu obliczenia kolejności wiązania wiązania węgiel–tlen między centralnymi jednostkami bipirydyny i fenylo–bipirydyny (Ce–O–Cd, Fig. 1). Obliczenia dotyczące najciaśniejszego węzła, 1, pokazują, że jedno z tych wiązań węgiel–tlen jest znacznie słabsze od innych i jest bardziej prawdopodobne, że pęknie w eksperymentach CID-MS (dodatek SI, Fig. S38). Po odcięciu powstały końcowy atom tlenu może zreorganizować swoją strukturę elektroniczną, aby umożliwić utratę cząsteczki wody (rys. 4G). W najbardziej luźnym węźle, 3, wszystkie wiązania węgiel-tlen mają podobne rozkazy wiązania, które są wyższe niż słabe Wiązanie C-O W 1. Może to skutkować niemal równoczesnym rozszczepianiem się dwóch wiązań C-O przy wyższej energii (dodatek SI, rys. S38), w celu wytworzenia dwóch połówek 3 zaobserwowanych doświadczalnie (rys. 4H).

badano również wpływ szczelności węzłów na właściwości widmowe 1-3. Widma UV-vis węzłów 1 i 2 są podobne, z niewielkim zmniejszeniem absorpcji przy λmax = 306 nm dla 2. Jednak absorpcja najbardziej luźnego węzła, 3, znacznie się zmniejsza i wiąże się z przesunięciem ku czerwieni o 10 nm (dodatek SI, Fig. S25). Symulowane widma, wygenerowane na podstawie teoretycznych obliczeń funkcjonalnych gęstości zależnej od czasu (37) na poziomie teorii M06-2x/6–31g(d), pokazują, że elektroniczne widma węzłów są bardzo wrażliwe na dwa kąty skrętne, θ1 i θ2, w obrębie ugrupowań fenylo–bipirydynowych (Fig. Dodatek 3D i SI, Fig. S35-S37). Obserwowane przesunięcie ku czerwieni węzła 3 jest spowodowane zmiennością θ1 i θ2, ponieważ zwiększona elastyczność luźnego węzła pozwala chromoforom zaludniać bardziej sprzężone, spłaszczone konformacje.

węzeł 819 jest wewnętrznie chiralny ze względu na swoją topologię. Enancjomery 1-3 izolowano chiralną wysokosprawną chromatografią cieczową (HPLC) i analizowano spektroskopią CD (Fig. 5 i dodatek SI, Fig. S23 i S24). Każda para enancjomerów węzłów daje widmo CD o jednakowym, ale przeciwstawnym kształcie i znaku (38, 39). Gdy węzły stają się luźniejsze od 1 do 3, zmiana środowiska wokół chromoforów powoduje przesunięcie sygnału CD w widmach odpowiednich cząsteczek. Symulacje i eksperymenty 1H NMR sugerują, że bipirydyny i etery fenolowe są zamknięte w chiralnej konformacji w najciaśniejszym węźle, 1, który charakteryzuje się silnymi oddziaływaniami π–π. Jednak luźniejsze węzły mają różne konformacje, które Średnio się dają słabszą odpowiedź CD. Przesunięcie ku czerwieni jest konsekwencją bardziej sprzężonych, płaskich konformacji chromoforów w luźniejszych węzłach. Na dynamikę i konformację splątanych struktur wyraźnie znacząco wpływa szczelność, a ekspresję chiralności można przestrajać przez dokręcenie lub rozluźnienie splątania.

5.

widma CD każdego enancjomeru cząsteczkowego 819 węzłów 1, 2 i 3.

Analiza symulacji dynamiki molekularnej sugeruje, że trendy właściwości obserwowane dla 819 węzłów 1-3 nie są specyficzną konsekwencją specyficznych aromatycznych pierścieni oporowych ich struktur. Fluktuacje atomowe mierzą elastyczność (dużych) układów molekularnych. Są one związane ze współczynnikiem temperatury B (40), który można zmierzyć eksperymentalnie i jest zwykle zgłaszany wraz ze strukturą białka w banku danych o białkach. W symulacjach dynamiki molekularnej wielkość fluktuacji atomowych pierścieni aromatycznych, łańcuchów alifatycznych i całego węzła została uznana za funkcję szczelności( Tabela 1): im mocniej węzeł, tym mniejsze fluktuacje. Trendy są zasadniczo podobne dla każdego z trzech monitorowanych elementów strukturalnych (aromatycznych, alifatycznych i całej cząsteczki).

zobacz tę tabelę:

  • Zobacz inline
  • Zobacz popup
Tabela 1.

fluktuacje atomowe, Δ (Å2), interakcje fragmentu aromatycznego, Efrag (kcal/mol) i energia odkształcenia łańcucha alkilowego, Estrain (kcal/mol), z trzech węzłów (1-3)

struktury migawek dynamiki molekularnej, które z powodzeniem symulowały dane eksperymentalne, wykorzystano do zbadania interakcji między fragmentami aromatycznymi a energiami szczepu łańcuchów alkilowych w trzech węzłach (Tabela 1). Obliczenia wykonano na poziomie teorii M06-2x/6-31g(d). Dla oddziaływań między pierścieniami aromatycznymi wyznaczono korektę błędu superpozycji zestawu bazowego za pomocą metody counterpoise (41). W kolejności malejącej szczelności węzłów (1 do 2 do 3) Średnia energia interakcji między fragmentami aromatycznymi zmniejsza się od -48,3 do -40,1 do -23,6 kcal / mol (Tabela 1). Zmienność odzwierciedla większą dynamikę nici, gdy węzły stają się luźniejsze. W tej samej serii średnie Energie szczepów łańcuchów alkilowych zwiększają się z 53,8, 160,0 i 260.1 kcal/mol, ponieważ do każdego łańcucha dodaje się więcej grup metylenowych, zmuszonych przez topologię węzła do przyjęcia napiętych kątów dwuściennych i starć sterycznych CH–CH. Interakcje układania aromatycznego nie generują preferowanych konformacji w bardziej elastycznych węzłach, ale raczej tworzą się w celu złagodzenia napięcia, gdy szczelność wiązania zmusza cząsteczkę do zwiększenia zwartości. Układanie pierścieni aromatycznych jest konsekwencją szczelności węzła, a nie siłą napędową jego konformacji/struktury. Pierścienie aromatyczne w łańcuchach ligandu tris (2,2 ’ – bipirydyny) są niezbędne do syntezy węzłów, a w niektórych przypadkach (1H NMR, widma CD) dostarczają użytecznych sond do eksperymentalnego wykrywania zachowania. Gdyby π-stacking był dominującą interakcją, która determinowała ich konformację, węzły nie byłyby wysoce dynamicznymi układami, które pokazują eksperymenty i symulacje. Dlatego wydaje się prawdopodobne, że trendy obserwowane dla 1-3 nie są bezpośrednio związane z ich szczególnym składem molekularnym, ale ogólnie będą odzwierciedlać aspekty zachowania dotyczące szczelności węzłów w nanoskali.

wnioski

seria molekularnych 819 węzłów różniących się jedynie długością łańcuchów alkilowych łączących sztywne regiony aromatyczne umożliwia badanie wpływu szczelności węzłów na szereg właściwości fizycznych i chemicznych, umożliwiając racjonalizację obserwacji eksperymentalnych za pomocą symulacji obliczeniowych. Szczelność węzłów cząsteczkowych 819 wpływa na reaktywność (szybkość demetlacji wiązanych ligandów i energie zderzeń wymagane do zerwania wiązania), konformację (sztywne regiony preferujące znajdować się poza obszarem splątania w luźniejszych strukturach) oraz ekspresję topologicznej chiralności (przejawiającej się w widmach CD i 1H NMR). Wyniki te dostarczają pewnych eksperymentalnych dowodów na skutki dokręcania węzłów na poziomie molekularnym, co powinno okazać się przydatne w zrozumieniu roli wiązania w splątanych polimerach oraz w projektowaniu przyszłych wiązanych i splecionych nanomateriałów (42 ⇓ ⇓ -45).

materiały i metody

synteza.

węzły molekularne 2 i 3 zostały przygotowane przez modyfikację strategii oplatania stosowanej wcześniej (29) do złożenia węzła 1 (rys. 1). Poddanie każdego metalizowanego węzła ((PF6) 7Cl, (PF6)7CL lub (PF6)7CL) do 1:1 roztwór MeCN/Naohaq (1 M) w temperaturze 80 °C, a następnie chromatografia wykluczająca wielkość, dawały odpowiednie zdemetalowane węzły 1-3, które charakteryzowały się NMR i spektrometrią mas. Enancjomery 1-3 izolowano chiralnym HPLC i analizowano spektroskopią CD.

Symulacje Molekularno-Dynamiczne.

symulacje dynamiki molekularnej przeprowadzono przy użyciu pakietu programów AMBER 12.0. Węzły zostały sparametryzowane za pomocą ogólnego pola siłowego bursztynu i standardowej procedury potencjału elektrostatycznego wykonanej w celu przypisania ładunków atomom za pomocą Antechamberu. Struktura krystalograficzna metalizowanego węzła 8+ została użyta jako struktura wyjściowa dla 1 oraz jako wzór dla 2 i 3.

podziękowania

dziękujemy China 1000 talents Plan, East China Normal University, Engineering and Physical Sciences Research Council (EP/P027067/1) i European Research Council, Advanced Grant 339019 za finansowanie, a University of Manchester do President ’ s Doctoral Scholar Award (do L. Z.). D. A. L. jest chińskim profesorem 1000 talentów „Topnotch Talent” i profesorem Royal Society Research.

Przypisy

  • ↵1 do kogo można kierować korespondencję. E-mail: francesco. zerbetto{at}unibo.it albo david.leigh{at}manchester.ac.uk.
  • wkład autora: L. Z., F. Z., and D. A. L. designed research; L. Z., J.-F. L., A. A., M. C., and F. Z. performed research; A. A., M. C., and F. Z. analyzed data; and L. Z., J.-F. L., A. A., M. C., F. Z., and D. A. L. wrote the paper.

  • autorzy nie deklarują konfliktu interesów.

  • Ten artykuł jest bezpośrednim przekazem PNAS.

  • Ten artykuł zawiera dodatkowe informacje online na www.pnas.org/lookup/suppl/doi:10.1073/pnas.1815570116/-/DCSupplemental.

opublikowany na licencji PNAS.