Přehled pro střední školy a uchazeče o studium na Přírodovědecké fakultě Univerzity Palackého v Olomouci.
Vyzkoušejte si reálné simulace v prohlížeči.
Úvodní slovo
Předmluva
Když se řekne chemie, téměř každý si ihned představí bílý plášť, laboratoř a chemické sklo. Málokdo ale ví, že chemie se dá studovat i pomocí počítačových modelů a simulací. Právě představení základních principů a aplikací počítačové chemie je předmětem této učebnice.
Cílem tohoto textu je informovat středoškolské studenty o tomto atraktivním a stále se rozvíjejícím oboru a dát jim solidní základ při hledání vysokoškolských oborů. Tato učebnice vznikla v rámci příprav bakalářského oboru Počítačová chemie na Přírodovědecké fakultě Univerzity Palackého v Olomouci.
💡 Jak pracovat s učebnicí
Text je doplněn o živá online cvičení přímo v prohlížeči. Každá kapitola začíná motivací a propojením s tématy, která již znáte ze střední školy (Newtonovy zákony, atomové orbitaly, elektronové konfigurace).
Kapitola 1
Úvod: Chemie bez baněk, aneb proč vůbec počítat?
Jedná se o obor, ve kterém se pomocí počítačových modelů studují chemické vlastnosti a procesy. Pod pojmem model se schovávají dvě klíčové části:
Struktura studovaného systému: Musíme počítači přesně definovat, z jakých atomů se skládá a jakou má geometrii (prostorové souřadnice x, y, z).
Výpočetní metoda: Fyzikální rovnice a algoritmy, které ze struktury spočítají hledané vlastnosti (energii, spektra, reaktivitu).
Proč počítat v počítači místo zkumavky?
1. Rozlišení: V 18 ml vody (1 mol) je neuvěřitelných 6·10²³ molekul. Experiment vždy měří průměr celého tohoto obrovského souboru. V počítači můžeme izolovat jedinou molekulu a prozkoumat její chování s atomární přesností.
2. Podpora experimentu: Přechodové stavy chemických reakcí trvají zlomky femtosekund a „nedáte je do baňky“. Počítač je jediná cesta, jak tyto klíčové meziprodukty spatřit.
3. Průzkum a screening: Počítač dokáže otestovat miliony chemických struktur (např. kandidátů na léčiva či materiálů pro baterie) během několika hodin, aniž bychom museli syntetizovat gram toxického odpadu.
⚡ Interaktivní cvičení 1
Trojúhelník kompromisů ve výpočetní chemii
Vyzkoušejte pravidlo „2 ze 3“. Klikněte nebo táhněte kurzorem v trojúhelníku a zjistěte, jaká výpočetní metoda nejlépe odpovídá zvolenému poměru přesnosti, velikosti a rychlosti.
Doporučená metoda
Bod D: Víceškálové metody (QM/MM)
Kombinace metod nebo využití neuronových sítí k dosažení kvantové přesnosti za zlomek času.
Škála: Tisíce atomů s kvantovou přesností
📌 Zapamatuj si:
Model = geometrická struktura + výpočetní metoda.
Metodu vybíráme jako kompromis mezi přesností, velikostí systému a dostupným výpočetním časem.
Kapitola 2
Kvantová chemie: Svět elektronů a jader
Všechny chemické vlastnosti – od barvy látky až po sílu chemické vazby – jsou určeny chováním elektronů v molekule. Protože je elektron téměř 2000krát lehčí než proton, nechová se jako klasická kulička, ale vykazuje vlnově-částicový dualismus.
Prostor okolo jádra, kde se elektron s největší pravděpodobností nachází, nazýváme atomový orbital (s, p, d, f). Kombinací atomových orbitalů vznikají molekulové orbitaly (MO).
⚡ Interaktivní cvičení 2
Potenciální křivka molekuly H₂ a překryv orbitalů
Posouvejte jezdcem vzdálenost mezi dvěma jádry vodíku R(H–H). Sledujte, jak se mění energie systému podle Morseova potenciálu a kde leží energetické minimum (vazebná délka 0,74 Å).
0.74 Å
Energie E: -4.75 eV
🌌 Vyzkoušejte plnohodnotný 3D Simulátor atomových orbitalů
Interaktivní webový nástroj od prof. Karla Berky (PřF UPOL). Umožňuje 3D Monte Carlo vizualizaci orbitalů pro libovolná kvantová čísla n, l, m.
Když v počítači sestavíme model molekuly, počáteční polohy atomů obvykle neodpovídají skutečnému rovnovážnému stavu – vazby mohou být příliš natažené nebo atomy nepřirozeně blízko u sebe. Cílem optimalizace geometrie je nalézt takové prostorové uspořádání jader, které má nejnižší možnou energii (tzv. energetické minimum).
Podle Bornovy–Oppenheimerovy aproximace jsou těžká atomová jádra vůči lehkým elektronům v podstatě nehybná. Pro každou pevnou pozici jader nejprve vyřešíme elektronovou strukturu (kvantový SCF cyklus) a získáme celkovou energii systému. Závislost energie na souřadnicích všech atomů tvoří vícerozměrnou plochu potenciální energie (PES – Potential Energy Surface), po které se při optimalizaci pohybujeme jako v kopcovité krajině.
🔁 Cyklus optimalizace probíhá ve čtyřech krocích:
Výpočet energie (SCF): Vyřešení Schrödingerovy rovnice pro aktuální polohy jader.
Výpočet sil (gradientu ∇E): Zjištění sil Fi = −∇iE působících na jednotlivé atomy (směr největšího růstu/poklesu energie).
Posun jader: Optimalizační algoritmus posune atomy ve směru klesající energie.
Test konvergence: Pokud jsou síly i změna energie menší než zadaná tolerance (|∇E| < ε), geometrie konvergovala a máme stabilní molekulu! V opačném případě se cyklus opakuje.
⚡ Interaktivní cvičení 2B
Optimalizace geometrie a hledání minima na ploše PES
Klikněte kamkoliv do vrstevnicové mapy energie a zvolte optimalizační algoritmus. Sledujte, jak se algoritmus propracovává údolím k minimu (červený křížek) a jak se liší trajektorie Nejprudšího spádu, Sdružených gradientů a metody BFGS.
📊 Průběh konvergence
Iterace: 0
Energie E: -- kcal/mol
Síla |∇E| (RMS gradient): --
⚪ Klikněte do mapy pro start
Sdružené gradienty (CG): Využívá informaci z předchozího kroku k volbě ortogonálního směru. V zakřivených údolích netrpí křížovými oscilacemi jako nejprudší spád a k minimu dorazí mnohem rychleji.
Kapitola 3
Klasická mechanika: Molekula jako kuličky a pružinky
Když potřebujeme simulovat velké systémy (tisíce molekul vody, proteiny či membrány), kvantová mechanika by byla příliš pomalá. Sáhneme proto po klasické molekulové dynamice (MD).
Atomy aproximujeme jako hmotné kuličky s parciálním nábojem a vazby mezi nimi jako mechanické pružinky řídící se Hookeovým zákonem:
F = −k · (r − r0)
⚡ Interaktivní cvičení 3
Vibrační pohyby molekul a Infračervená (IR) spektroskopie
Zvolte normální vibrační mód a sledujte animaci kmitání atomů i odpovídající absorpční pík na infračerveném spektru.
H₂O: Deformační kmit | Frekvence: 1595 cm⁻¹
Mění se úhel vazby H-O-H. Atomy vodíku se synchronně přibližují a oddalují od sebe.
Kapitola 4
Hrubozrnné a víceškálové metody: Od atomů k biologickým systémům
Biologické procesy v buňce probíhají na obrovských škálách milionů atomů a v časových horizontech milisekund. V hrubozrnných modelech (Coarse-Grained / CG), jako je silové pole MARTINI, slučujeme skupiny 3–4 těžkých atomů do jediných zjednodušených interakčních koulí (tzv. beadů). Tím dramaticky snižujeme počet stupňů volnosti a můžeme simulovat obří buněčné motory, biologické membrány či virové obaly.
Pokud potřebujeme studovat enzymatickou reakci v obřím proteinovém komplexu, propojíme obě metody do víceškálového modelu QM/MM: samotné reaktivní centrum popíšeme exaktně kvantovou mechanikou (QM) a zbytek obřího proteinu s rozpouštědlem klasickou mechanikou (MM). Za vývoj víceškálových modelů byla v roce 2013 udělena Nobelova cena za chemii (Martin Karplus, Michael Levitt, Arieh Warshel).
⚡ Interaktivní cvičení 4
3D Mesoscale & Hrubozrnný prohlížeč (Mol*)
Vyzkoušejte si interaktivní 3D prohlížení obřích makromolekulárních komplexů a mesoscale struktur přímo v prohlížeči pomocí špičkového vědeckého nástroje Mol* (Molstar.org).
Bakteriální bičíkový motor (Motor Hook): Molekulární kloub a univerzální převodovka rotujícího bakteriálního bičíku složená ze stovek proteinových podjednotek. V hrubozrnném zobrazení vidíme dynamické uspořádání celého molekulárního stroje.
Chemická informatika: Molekuly jako data a návrh léčiv
Chemická informatika (chemoinformatika) propojuje chemii s datovou vědou a počítači. Místo zdlouhavých fyzikálních simulací reprezentuje molekuly jako data a umožňuje prohledávat knihovny se stovkami milionů sloučenin během zlomku vteřiny.
📝 Jak počítač kóduje molekulu?
SMILES: Textový zápis struktury (např. c1ccccc1 pro benzen, CCO pro ethanol).
InChIKey: 27místný kontrolní hash řetězec IUPAC pro vyhledávání v databázích.
Molekulární graf (MOL / SDF): Atomy jako vrcholy, kovalentní vazby jako hrany grafu.
💊 Lipinského pravidlo pěti (Rule of 5)
Předpovídá, zda má molekula šanci stát se účinným perorálním lékem (vstřebatelným trávicím traktem):
V moderní vědecké praxi a farmaceutickém průmyslu se k výpočtu deskriptorů využívají chemoinformatické knihovny RDKit (standard pro Python / C++), CDK (Chemistry Development Kit) a Open Babel. Tyto nástroje využívají exaktní matematické a empirické algoritmy:
Molekulová hmotnost (MW): Exaktní sumace atomových hmotností z grafu molekuly podle tabulek IUPAC.
Lipofilita (logP):Wildman–Crippen model rozkládá molekulu na 68 atomových typů a sčítá jejich experimentálně kalibrované hydrofobní příspěvky (SlogP).
Rozpustnost ve vodě (logS):Delaneyho ESOL model počítá rozpustnost ze vzorce: logS = 0.16 − 0.63·logP − 0.0062·MW + 0.066·RotB − 0.74·(Aromatický podíl).
Polární povrch (TPSA):Ertl–Rohde topologický algoritmus sčítá 2D tabulkové plochy atomů O, N a jejich navázaných H bez nutnosti generovat náročnou 3D konformaci.
H-donory / akceptory a aromatické kruhy: Hledání podstruktur pomocí chemických regulárních výrazů (SMARTS) a detekce cyklů (SSSR algoritmus).
Výzkumníci na Katedře fyzikální chemie PřF UPOL pracují v rámci infrastruktury ELIXIR CZ s celosvětovými databázemi ChEMBL (2,4+ milionu bioaktivních látek), PubChem (115+ milionů struktur) a ZINC (miliardy sloučenin pro virtuální screening).
⚡ Interaktivní cvičení 5
2D Editor molekul (Ketcher / SmilesDrawer) & kalkulátor vlastností
Nakreslete vlastní molekulu na plátno pomocí tlačítek a šablon, nebo vyberte předvolené léčivo. Editor v reálném čase generuje SMILES, vykresluje profesionální 2D strukturu (včetně aromatických jader) a počítá molekulární deskriptory podle standardů RDKit.
Návrh léčiv a dokování: Klíč, zámek a virtuální screening
Většina léků funguje tak, že se jako klíč zasune do specifické dutiny (aktivního místa) enzymu a zablokuje jeho nežádoucí činnost. Molekulární dokování je výpočetní metoda, která simuluje usazení malé molekuly do vazebné kapsy a vyhodnotí vazebnou afinitu (ΔG).
⚡ Interaktivní cvičení 6
Interaktivní dokovací simulátor: Klíč a zámek
Přetáhněte molekulu ligandu myší/dotykem do aktivního místa enzymu. Tlačítkem „Otočit ligand“ vyzkoušejte správnou orientaci pro vytvoření vodíkových můstků (červený akceptor, modrý donor, žlutý hydrofobní kruh).
Vazebná energie ΔG: -2.0 kcal/mol
❄️ Volný ligand mimo vazebnou kapsu.
Kapitola 7
Skládání bílkovin: Od sekvence k funkci (AI revoluce)
Po více než 50 let bylo předpovězení 3D prostorového tvaru bílkoviny pouze z její lineární sekvence aminokyselin jedním z největších nevyřešených problémů biologie (tzv. Levinthalův paradox).
V roce 2020 přišel zlom díky umělé inteligenci AlphaFold od společnosti Google DeepMind a modelu ESMFold. V roce 2024 byla za tento objev udělena Nobelova cena za chemii (Demis Hassabis, John Jumper, David Baker).
⚡ Interaktivní cvičení 7
3D Prohlížeč proteinů & AlphaFold DB (Mol*)
Vyberte protein z databáze a prohlédněte si jeho 3D strukturu přímo v interaktivním Mol* okně. U modelů z AlphaFold DB uvidíte zbarvení podle predikované spolehlivosti (pLDDT skóre). Můžete myší strukturu libovolně otáčet a přibližovat.
Cytochrom P450 3A4 (AlphaFold DB: P08684) (pLDDT: 94.8%) — Nejdůležitější lidský enzym zodpovědný za metabolickou přeměnu a odbourávání více než 50 % všech klinicky užívaných léčiv v játrech.
Kapitola 8
Studium reakcí: Reakční koordináta, přechodový stav a katalýza
Aby chemická reakce proběhla, musí reaktanty překonat energetickou bariéru – aktivační energii (Ea). Na vrcholu této bariéry leží přechodový stav (transition state, značeno ‡) – nestabilní uspořádání jader s nejvyšší energií na reakční dráze (sedlový bod 1. řádu na ploše potenciální energie).
⛰️ Přechodový stav (‡) & SN2 substituce
U bimolekulární nukleofilní substituce (např. F⁻ + CH₃Cl → CH₃F + Cl⁻) dochází k útoku nukleofilu zezadu a k Waldenově inverzi (obrácení deštníku).
V přechodovém stavu [F···CH₃···Cl]‡ má centrální uhlík 5 vazebných kontaktů v trigonálně bipyramidálním uspořádání.
🧪 Jak funguje katalyzátor?
Katalyzátor (enzym nebo kovový povrch) nemění celkovou termodynamickou rovnováhu (ΔH), ale nabízí alternativní reakční cestu s nižší aktivační energií (Ea).
Dle Arrheniovy rovnice \(k = A \cdot e^{-E_a / RT}\) snížení bariéry o pouhých 25 kJ/mol zrychlí reakci při pokojové teplotě více než 25 000×!
⚡ Interaktivní cvičení 8
Reakční koordináta, přechodový stav & vliv katalýzy
Posouvejte jezdec reakční koordináty (nebo spusťte animaci) a sledujte přeměnu reaktantů přes přechodový stav v grafu energie i v živé molekulární geometrii. Vyzkoušejte vliv katalyzátoru na snížení bariéry a zrychlení reakce.
ReaktantyProdukty
📈 Profil potenciální energie (PES)ΔE = 0.0 kJ/mol
⚛️ Molekulární geometrie podél reakční dráhyξ = 0 %
Aktuální fáze: Reaktanty
Aktivační bariéra (Eₐ): 42.0 kJ/mol
Reakční teplo (ΔH): -125.0 kJ/mol (Exotermní)
Zrychlení katalyzátorem: 1× (bez katalýzy)
Kapitola 9
Materiálové vlastnosti: Fotovoltaika, baterie a nanomateriály
Výpočetní metody hrají zásadní roli v moderní energetice. V Olomouci na Katedře fyzikální chemie PřF UPOL se vědci zabývají modelováním:
2D grafenových derivátů (fluorografen, grafénová kyselina) pro superkondenzátory.
Materiálů pro akumulátory (rychlejší transport lithných a sodných iontů bez rizika vznícení).
Perovskitových solárních článků s vysokou účinností přeměny světla na elektřinu.
Kapitola 10
Kým se stát? Co dělají výpočetní chemičky a chemici?
Práce výpočetního chemika spojuje kreativitu vědce a flexibilitu IT specialisty. Místo práce v digestoři trávíte den modelováním molekul, psaním skriptů v Pythonu a diskusemi nad 3D strukturami s experimentátory z celého světa.
Chcete se stát součástí tohoto fascinujícího světa? Podejte si přihlášku na studijní program Počítačová chemie na PřF UPOL a otevřete si dveře k moderní kariéře!
