Analiza konformacyjna polega na badaniu zmian energii cząsteczki, na ogół związku organicznego, w zależności od jego geometrii. Ponieważ wszystkie geometrie muszą odpowiadać tej samej chemicznie cząsteczce (a nie np. jej izomerom), zmiennymi opisującymi powierzchnię energii są kąty torsyjne obrotu wokół wiązań, z wyłączeniem obrotu wokół wiązań podwójnych oraz nieskorelowanych obrotów wokół wiązań w pierścieniach. Przykładowo: kwas fumarowy i maleinowy, które różnią się położeniem grup karboksylowych względem osi wiązania C=C, są izomerami a nie konformerami.
Podobnie, tylko niektóre zestawy kątów obrotu wokół wiązań C-C w pierścieniu cykloheksanu odpowiadają konformacjom tego związku i aby zmienić np. konformację krzesłową w konformację łódkową należy jednocześnie odpowiednio zmienić trzy wszystkie kąty dwuścienne w pierścieniu. W przeciwnym razie nastąpi rozerwanie któregoś z wiązań w pierścieniu (patrz rysunek).
To, że długości wiązań i kąty walencyjne są w przybliżeniu stałe a duża część kątów torsyjnych w cząsteczce jest stała lub współzmienna z innymi kątami torsyjnymi eliminuje ogromną większość stopni swobody opisujących powierzchnię energii konformacyjnej. Tym niemniej, nawet dla stosunkowo małych cząsteczek takich jak oligopeptydy czy oligosacharydy pozostaje ich wystarczająco dużo aby znalezienie niskoenergetycznych (a przez to najbardziej prawdopodobnych) konformacji, które są identyfikowane z minimami energii, było zagadnieniem bardzo trudnym. Każda jednostka peptydowa lub cukrowa wymaga dwóch kątów torsyjnych (oznaczanych jako i ) do opisu, zatem już w przypadku pentapeptydu jakim jest np. naturalny opioid Met-enkefalina, liczba zmiennych wynosi 10.
Przykład i zadania z tej części ćwiczeń są ograniczone do prostych cząsteczek, których konformacje określa jeden kąt torsyjny oraz jedynie do obliczania energii wybranych konformacji. Dla zobrazowania zagadnienia poniżej jest przedstawiony dobrze zapewne znany z kursu chemii organicznej wykres energii cząsteczki butanu w funkcji kąta obrotu wokół centralnego wiązania C-C. Zaniedbujemy tutaj mniej istotne zmiany energii spowodowane zmianami kątów obrotu skrajnych grup CH wokół wiązań C-C.
Znajdowanie minimów energii konformacyjnej dla podanego już wyrażenia na energię potencjalną w funkcji kąta lub kątów torsyjnych było treścią tematu 2, natomiast w tej części ćwiczeń wyrażenia na energię należy utworzyć na podstawie ogólnego wzoru na energię konformacyjną w przybliżeniu mechaniki molekularnej oraz podanych w zadaniach parametrów. W zadaniach z obecnej części zakłada się, że długości wiązań i kąty walencyjne nie zmieniają się, zatem jedynymi wkładami do energii konformacyjnej są wkłady niewiążące (), elektrostatyczne () i torsyjne ().
Siła działająca na dany atom w cząsteczce lub układzie cząsteczek określa zmianę jego prędkości a dalej trajektorię ruchu. Dlatego obliczanie sił jest kluczowym punktem dynamiki molekularnej, która umożliwia symulacje ewolucji czasowej układów a przez to na określenie np. funkcjonlnie ważnych ruchów cząsteczki białka (np. cytochromu w czasie oddychania komórkowego) czy też ruchliwości poszczególnych fragmentów białka, co ma znaczenie np. w ocenie wpływu mutacji na kancerogenność. Obliczanie sił elektrostatycznych pozwala z kolei ocenić miejsca cząsteczki receptora najbardziej podatne na wiązanie z kandydatami na leki.
Aby wyprowadzić sposób obliczania sił na podstawie wyrażeń na energię w polu siłowym przypomnijmy sobie znany ze szkoły wzór na pracę:
gdzie jest pracą wykonaną przez siłę a jest drogą. Wykonana praca jest tym większa im większa jest siła i dłuższa droga.
Jeżeli praca jest wykonywana przez siły potencjalne (konserwatywne) takie, jak np. siły grawitacyjne czy elektrostatyczne, energia układu zmaleje. Swobodne spadanie ciał w polu grawitacyjnym ziemskim jest doskonałą ilustracją tego prawidła. Jeżeli ciało o masie spadnie z wysokośći to jego energia potencjalna będzie mniejsza o ( jest przyspieszeniem ziemskim) w stosunku do energii w położeniu początkowym. Zatem możemy napisać:
Zakładając, że ruch odbywa się wzdłuż osi i oznaczając przesunięcie oraz przechodząc do granicy dostajemy:
Siła działająca w kierunku jest zatem pochodną cząstkową energii potencjalnej względem . Podobnie, składowe siły w kierunkach współrzędnych i są pochodnymi cząstkowymi energi potencjalnej względem tych współrzędnych a cały wektor siły, , jest wektorem gradientu energii potencjalnej ze znakiem minus:
Powyższe wyrażenie pozwala, na podstawie wyrażenia na energię potencjalną, obliczyć siłę działającą na wybrany atom. Jeżeli chcemy obliczyć siły działające na wszystkie atomy, musimy odróżnić od siebie współrzędne poszczególnych atomów nadając im numery, np. . Wektory siły działającej na atom o numerze otrzymujemy wtedy licząc gradient w kierunku jego współrzędnych:
Ponieważ obliczanie pochodnych całego wyrażenia na energię nawet dla niewielkich układów jest uciążliwe wykorzystuje się to, że większość wkładów do energii, w szczególności elektrostatyczne i niewiążące ale także oddziaływania pochodzące od odkształcania wiązań, wyraża się w mechanice molekularnej poprzez sumę wkładów pochodzących od par atomów, z których każdy zależy jedynie od odległości między tymi atomami:
gdzie jest odległością między atomami oraz . Korzystając z wzoru na pochodną funkcji złożonej mamy dla współrzędnej i-tego atomu (analogiczne wyprowadzenie dla pozostałych):
Ponieważ
Zatem
Zatem aby obliczyć siłę działającą na atom należy dla każdego oddziałującego z nim atomu obliczyć pochodną tego wkładu energii oddziaływań względem odległości tych atomów, zmienić znak tej pochodnej, przemnożyć przez wektor kierunku od atomu do atomu (łatwo sprawdzić, że wektor ten ma długość 1) a następnie wektorowo zsumować otrzymane w ten sposób składowe siły. Zwróćmy uwagę, że z powyższego wzoru wynika III prawo Newtona: siła pochodząca od atomu a działająca na atom ( ) jest co do wielkości równa a przeciwnie skierowana do tej pochodzącej od atomu i działającej na atom ( ), ponieważ w drugim przypadku pochodna będzie mnożona przez wektor , który jest przeciwnie skierowany do wektora .
Powyższe rozważania ilustruje rysunek.
Jeżeli energię wyrażamy w kcal/mol a odległość w angströmach, to podstawiając takie jednostki do wzoru na siłę dostajemy kcal/(mol). Tymczasem jednostką siły jest newton. Poniżej jest podane wyprowadzenie współczynnika przeliczającego siłę na newtony:
Konformacja naprzeciwległa (I) | Konformacja naprzemianległa (II) | |||||
Atom | x [] | y [] | z [] | x [] | y [] | z [] |
C1 | 0,000000 | 0,000000 | 0,000000 | 0,000000 | 0,000000 | 0,000000 |
C2 | 0,000000 | 0,000000 | 1,460000 | 0,000000 | 0,000000 | 1,460000 |
Cl1 | 1,649206 | 0,000000 | -0,642667 | -1,649206 | 0,000000 | -0,642667 |
Cl2 | -0,824603 | -1,428254 | -0,642667 | 0,824603 | 1,428254 | -0,642667 |
Cl3 | -0,824603 | 1,428254 | -0,642667 | 0,824603 | -1,428254 | -0,642667 |
Cl4 | 1,649206 | 0,000000 | 2,102667 | 1,649206 | 0,000000 | 2,102667 |
Cl5 | -0,824603 | 1,428254 | 2,102667 | -0,824603 | 1,428254 | 2,102667 |
Cl6 | -0,824603 | -1,428254 | 2,102667 | -0,824603 | -1,428254 | 2,102667 |
Dodatkowo, z uwagi na symetrię obu konformacji mamy oraz . Zatem wystarczy policzyć energie oddziaływań pomiędzy ClCl i ClCl, po czym pierwszą z nich pomnożyć przez 3 a drugą przez 6 i zsumować.
Atomy Cl1 i Cl4 różnią się tylko współrzędną więc odległość między nimi jest równa bezwzględnej wartości różnicy tej współrzędnej.
Parametry potencjału Lennarda-Jonesa dla pary oddziałujących atomów chloru wynoszą:
Energia oddziaływań niewiążących i elektrostatycznych wynoszą odpowiednio:
Energia torsyjna wynosi
Zatem energia konformacji naprzeciwległej wynosi kcal/mol.
Wartość eksperymentalna bariery rotacji heksachloroetanu wg. Y. Morino, E. Hirota, J. Chem. Phys. 28, 185 (1958) wynosi kcal/mol. Jakie są możliwe przyczyny, że obliczona bariera rotacji jest prawie dwukrotnie wyższa?
0,96 | 1,36 | |||
O1 | C2 | O3 |
Energia odkształcenia każdego wiązania C=O jest dana wzorem ( oznacza długość wiązania):
ponieważ z rysunku wynika, że suma długości obu wiązań C=O wynosi 2,32 . Zatem energia wyraża się następujące przez x:
Siła działająca w kierunku współrzędnej wynosi:
Wstawiając dostajemy kcal/mol/= N (przeliczenie jednostek jest wyprowadzone w części teoretycznej). Siła działa w prawo ponieważ wiązanie O1=C2 jest krótsze a C2=O3 dłuższe niż nienaprężone wiązanie o długości równowagowej 1,16 .
Sposób podany powyżej jest najszybszym sposobem rozwiązania tego konkretnego zadania i dobrze ilustruje definicję siły jako wektora przeciwnego do gradientu energii potencjalnej. Jednak dla bardziej złożonych wyrażeń na energię znacznie lepiej jest najpierw obliczyć składowe siły pochodzące od oddziaływań danego atomu z poszczególnymi atomami (w tym przypadku są to wyłącznie oddziaływania wiążące opisane potencjałem harmonicznym) a następnie je zsumować, jak opisano w części teoretycznej. Ten ogólniejszy sposób jest zilustrowany poniżej.
Przy składowej siły pochodzącej od wiązania C2O1 występuje 1, ponieważ wektor ma kierunek osi , natomiast przy składowej siły pochodzącej od wiązania C2O3 występuje -1, ponieważ wektor jest skierowany przeciwnie do osi .
ponieważ długość wiązania Cl1H2 jest równa jego długości równowagowej. Na rysunku pokazane są również wektory jednostkowe skierowane do atomu na który działa siła od atomów działających na niego siłą, oraz . Ponieważ rozważamy tylko jeden wymiar (oś ), pierwszy wektor (od Cl1 do H2) jest równy 1 a drugi (od Cl3 do H2) -1.
Energia oddziaływania pomiędzy H2 a Cl3 wyraża się następująco (zauważmy, że kcal/mol, ); we wzorze , ponieważ występuje tylko jedna odległość:
Po obliczeniu pochodnej względem dostajemy:
Wstawiając do powyższego wyrażenia dostajemy kcal/mol/ = N. Siła jest skierowana w prawo co jest zgodne z przewidywaniem, że dodatnio naładowany atom H2 cząsteczki chlorowodoru jest przyciągany przez ujemnie naładowany jon chlorkowy. Jednak przy szacowaniu sił działających na atomy nie należy pochopnie uwzględniać tylko jednego rodzaju oddziaływń (w tym przypadku elektrostatycznych). Przykładem, że takie oszacowanie może prowadzić do mylnych wniosków jest zadanie 4.
Dwa jony sodowe i dwa jony chlorkowe w próżni są naprzemiennie umieszczone w wierzchołach kwadratu. Zakładając, że wszystkie odległości są większe niż suma promieni van der Waalsa poszczególnych par jonów, naszkicować siły działające na każdy jon oraz ich wypadkowe i przewidzieć zachowanie układu.
Ponieważ wszystkie odległości są większe niż suma promieni van der Waalsa oddziałujących jonów a układ znajduje się w próżni, oddziaływania elektrostatyczne będą zdecydowanie dominowały nad niewiążącymi. Wystarczy zatem uwzględnić jedynie siły elektrostatyczne. (Uwaga! Sytuacja mogła by się zmienić gdyby układ znajdował się w wodzie czy innym rozpuszczalniku o dużej stałej dielektrycznej.) Na każdy jon sodowy działają siły pochodzące od sąsiadujących z nim jonów chlorkowych oraz jonu sodowego leżącego na przekątnej kwadratu a sytuacja jest analogiczna dla jonów chlorkowych. Siły działające wzdłuż boków kwadratu (ciemnoniebieskie strzałki) są skierowane do przeciwjonu a ich wypadkowe do środka kwadratu wzdłuż jego przekątnej (jasnoszare strzałki). Jeżeli wartość każdej składowej oznaczymy to wartość wypadkowej wynosi . Z kolei siła działająca od leżącego w przeciwległym rogu jonu tego samego znaku jest skierowana od środka kwadratu (odpowiednio zielone i różowe strzałki) a jej wartość jest dwukrotnie mniejsza niż ponieważ
W takim razie Wypadkowa siła działająca na każdy jon (czarne strzałki) jest skierowana do środka kwadratu wzdłuż jego przekątnej. Zatem jony będą poruszały się po przekątnych kwadratu w kierunku jego środka, aż odpychanie elektronów związane z zakazem Pauliego (wszystkie jony są układami zamkniętopowłokowymi) zrównoważy przyciąganie elektrostatyczne i przyciąganie wynikające z oddziaływań dyspersyjnych. Gdyby jony były ładunkami punktowymi nie istniałoby położenie równowagi i wszystkie jony zbiegłyby się w jednym punkcie.
Konformacja prostopadła (I) | Konformacja naprzeciwległa (II) | |||||
Atom | x [] | y [] | z [] | x [] | y [] | z [] |
H1 | 0,000000 | -0,895669 | -0,316667 | 0,895669 | 0,000000 | -0,316667 |
O2 | 0,000000 | 0,000000 | 0,000000 | 0,000000 | 0,000000 | 0,000000 |
O3 | 0,000000 | 0,000000 | 1,480000 | 0,000000 | 0,000000 | 1,480000 |
H4 | 0,895670 | 0,000000 | 1,796667 | 0,895670 | 0,000000 | 1,796667 |
Długości wiązań i kąty walencyjne są takie same dla obu konformacji.
Energia torsyjna jest dana wzorem:
Ładunki na atomach tlenu i wodoru wynoszą odpowiednio e i e. Stałe potencjału Lennarda-Jonesa wynoszą kcal/mol, , kcal/mol, .
0,93 | 2,87 | |||
F1 | — | H2 | &dotfill#dotfill; | F3 |
-0,6e | +0,6e | -1,0e |