Name: Wpływ izotopowy porficenu w procesie podwójnego transferu protonów badany metodą Enhanced QM/MM (Wideo) | JoVE Journal
Uploaded: 2019-07-19T15:00:00.000+00:00
Duration: 5 min 51 s
Description: Poznaj metodę Six QM/MM do badania izotopowych efektów w podwójnym transferze protonów w porficenie. Protokół umożliwia analizę mechanizmów reakcji bez predefiniowania współrzędnych reakcji.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

playerTrialDesktop

Te badania są wspierane przez Narodowy Kluczowy Program Badawczo-Rozwojowy Chin (2017YFA0206801, 2018YFA0208600), Fundację Nauk Przyrodniczych prowincji Jiangsu oraz Narodową Fundację Nauk Przyrodniczych Chin (91645116). L.X jest specjalnie mianowanym profesorem Zhong-Wu na Uniwersytecie Technologicznym Jiangsu. Autorzy przyjmują do wiadomości sugestie dr Hao Hu i dr Mingjun Yanga.

1. Model budynku
<ol>
	<li>Zbuduj strukturę porficonu: Otwórz oprogramowanie GaussView, klikając dwukrotnie myszą. Następnie kliknij przycisk Fragment elementu w menu GaussView, aby wybrać potrzebne elementy. Skonstruuj porficen. Następnie kliknij przycisk Plik, aby zapisać jako plik pdb.</li>
	<li>Solvate model: Solvate porphycen in a cubic TIP3P21 water box o długości krawędzi 38 A, wydając polecenie w systemie operacyjnym Linux: genbox_d -cp prp-vac.pdb -cs spc216.gro -o solv.pdb -maxsol 1484 -box 3.8.</li>
	<li>Zbuduj porficon deuterowy: Wydaj następujące polecenie, aby wygenerować plik topologii: cns < ppi_solv.inp. Następnie otwórz prp-wat.psf za pomocą polecenia vi i zmień masę H1 z 1.00800 na 2.01600, aby zastąpić jeden wewnątrzcząsteczkowy atom wodoru porficenu deuterem, aby zbudować pojedynczy porficen podstawiony deuterem.</li>
	<li>Ustaw normalne parametry symulacji MD: Metoda wprowadzania scctb, całka 0,5 fs i odcięcie 12 w pliku wejściowym MD, otwierając go za pomocą polecenia vi. 
	UWAGA: Przyjmij odległość odcięcia 12 A do obliczania zarówno vdW, jak i oddziaływań elektrostatycznych. Symulacja cząsteczki porficenu za pomocą metody DFTB/MIO22. Ustaw krok czasu integracji na 0,5 fs dla symulacji MD. Utrzymuj temperaturę systemu symulacji na poziomie 300 K za pomocą termostatu Langevin. Następnie wykonaj symulacje za pomocą oprogramowania QM4D, postępując zgodnie z poniższymi krokami.</li>
</ol>
2. Wstępne siedzenia
<ol>
	<li>Ustaw parametry temperatury: Wprowadź templow 260, temphigh 1100 i ntemp 160 w pliku wejściowym. 
	UWAGA: Zakres temperatur od 260K do 1100K został rozłożony przez oprogramowanie QM4D do 160 punktów temperatury podczas symulacji MD. Pliki wejściowe szablonów są zawarte w plikach uzupełniających.</li>
	<li>Zainicjuj ustawienia wstępne: Ustaw typ działania 100 i krok 120 000 w pliku wejściowym. Następnie wydaj następujące polecenie: $PATH/qm4d $INPUTFILE > $OUTPUTFILE. 
	UWAGA: Całkowity krok wynosi 120 000, ale można go dostosować w zależności od konkretnych potrzeb>. Sugerowane parametry symulacji MD zapisywane są w $INPUTFILE. To samo polecenie jest również używane w następujących krokach symulacji opcji i produkcji, z odpowiednio zmodyfikowanym plikiem wejściowym.</li>
	<li>Obliczanie rozłożonych energii
	<ol>
		<li>Wyodrębnij zmiany energii: Podczas etapu wstępnego monitorowania energii każdego członu, aby obliczyć wartości średnie, jak pokazano na <silna klasa="xfig">Rysunek 1. Użyj polecenia grep Linux, aby wydobyć energię w następujący sposób: 
		grep 'SITS-ener0' $INPUTFILE | awk '{a+=$3; b+=4 dolary; c+=$5}END{print a/NR,b/NR,c/NR}.</li>
		<li>Zmodyfikuj średnie energie w pliku wejściowym MD: Oblicz średnie energie na podstawie danych wyjściowych wiersza poleceń powyżej i zmodyfikuj linię vshift0 -30801.95; vshift1 -26.88; vshift2 -13888.28 w pliku wejściowym z nowo wygenerowanymi średnimi. 
		UWAGA: Liczby -30801.85, -26.88 i -13888.28 są średnimi energiami w obecnym systemie modelowym. Zmodyfikuj wartości w oparciu o konkretne systemy.</li>
	</ol>
	</li>
</ol>
3. Opcja
<ol>
	<li>Inicjowanie opt-sits: Ustawia typ działania 0 w pliku wejściowym. Następnie zainicjuj program QM4D, wpisując polecenie, jak pokazano w kroku 2.2, aby rozpocząć krok optymalizacji.</li>
	<li>Monitorowanie zmian energii i wartości nk.
	<ol>
		<li>Wykreśl propagację energii z programem "grace" i upewnij się, że fluktuacja energii może pokryć najniższy i najwyższy koniec zakresu temperatur.</li>
		<li>Po optymalizacji zapisz końcowe wartości nk kroku opt-sits do nowego pliku, który jest nazwany jako nk.dat w tym protokole.</li>
	</ol>
	</li>
</ol>
4. Przeprowadzanie symulacji produkcyjnych
<ol>
	<li>Przygotuj plik wejściowy MD: Ustaw typ działania 1 w nowym pliku wejściowym, aby rozpocząć krok symulacji produkcji. Określa nazwę przechowywanego pliku nk jako nkfile nk.dat w pliku wejściowym. Liczba kroków czasowych została ustalona na 6 400 000 w obecnych systemach.</li>
	<li>Inicjowanie produkcyjnej symulacji MD: Wydaj następujące polecenie, aby uruchomić symulacje MD: $PATH/qm4d $INPUTFILE > $OUTPUTFILE. 
	UWAGA: Upewnij się, że wartości >nk mogą być odczytane przez oprogramowanie QM4D. Czas symulacji jest zależny od systemu, więc zmień krok symulacji w zależności od konkretnych wymagań. Wybierz odpowiednią liczbę kroków czasowych, aby zapewnić wystarczającą ilość czasu na symulację dla własnego systemu. Ten krok jest prawdopodobnie czasochłonny, więc zapisz pliki ponownego uruchomienia, aby uniknąć ponownego uruchamiania produkcji od samego początku po zakłóceniu.</li>
</ol>
5 Analiza danych
<ol>
	<li>Monitorowanie zmian odległości
	<ol>
		<li>Monitoruj proces tworzenia i zrywania wiązań w fazie produkcji, użyj polecenia grep, aby sprawdzić zmiany odległości H1-N1 i H1-N2 w czasie symulacji. Tę samą operację można przeprowadzić dla H2-N3 i H2-N4. Następnie należy wykreślić propagację odległości, korzystając ze skumulowanej wartości odległości podczas symulacji produkcyjnych.</li>
	</ol>
	</li>
	<li>Wyodrębnianie współrzędnych reakcji
	<ol>
		<li>Wyodrębnij współrzędne reakcji i warunki energetyczne z pliku wyjściowego produkcji wygenerowanego przez QM4D za pomocą polecenia grep: 
		grep 'Dyc 1' $OUTPUTFILE | awk '{print $5}' > odległość1; 
		grep 'ener0' $OUTPUTFILE > ener0.</li>
		<li>Uporządkuj dane w czterech kolumnach: q1, q2, U0 i U' (U0 i U' to normalna energia wyjściowa i energia ważona) i zapisz je w pliku danych w każdym przedziale czasowym.</li>
	</ol>
	</li>
	<li>Obliczanie darmowej energii
	<ol>
		<li>Oblicz energię swobodną, wydając następujące polecenie: 
		sits-pmf 300 $INPUTFILE PMF2 [hist_minx hist_maxx num_binsx] [hist_miny hist_maxy num_binsy] > $OUTPUTFILE. 
		UWAGA: sits-pmf to metoda analizy oparta na histogramie. [hist_minx hist_maxx num_binsx] określa zakres i liczbę przedziałów dla pierwszej współrzędnej reakcji. Drugie współrzędne reakcji można ustawić za pomocą [hist_miny hist_maxy num_binsy].</li>
		<li>Aby rzutować energię swobodną na dwuwymiarowy krajobraz, wpisz następujące polecenie: 
		sits-pmf 300 h1-2d.dat PMF2 -0,6 0,6 24 2,45 4,25 36 > sits-pmf.out. 
		UWAGA: Użyj łącznie 24 przedziałów i 36 pojemników, aby pokryć zmiany odległości w dwóch wybranych współrzędnych reakcji, odpowiednio q>1 i q2. Zapisz dane PMF 2D dla każdego wodoru/deuteru w pliku sits-pmf.out.</li>
	</ol>
	</li>
</ol>

<ol><li>Waluk, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Ground-+and+excited-state+tautomerism+in+porphycenes%5BTitle%5D)+AND+%22Accounts+of+Chemical+Research%22%5BJournal%5D)">Ground- and excited-state tautomerism in porphycenes.</a> Accounts of Chemical Research. 39 (12), 945-952 (2006).</li>
<li>Waluk, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Spectroscopy+and+tautomerization+studies+of+porphycenes%5BTitle%5D)+AND+%22Chemical+Reviews%22%5BJournal%5D)">Spectroscopy and tautomerization studies of porphycenes.</a> Chemical Reviews. 117 (4), 2447-2480 (2017).</li>
<li>Kozlowski, P. M., Zgierski, M. Z., Baker, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+inner-hydrogen+migration+and+ground-state+structure+of+porphycene%5BTitle%5D)+AND+%22The+Journal+of+Chemical+Physics%22%5BJournal%5D)">The inner-hydrogen migration and ground-state structure of porphycene.</a> The Journal of Chemical Physics. 109 (14), 5905-5913 (1998).</li>
<li>Yoshikawa, T., Sugawara, S., Takayanagi, T., Shiga, M., Tachikawa, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Theoretical+study+on+the+mechanism+of+double+proton+transfer+in+porphycene+by+path-integral+molecular+dynamics+simulations%5BTitle%5D)+AND+%22Chemical+Physics+Letters%22%5BJournal%5D)">Theoretical study on the mechanism of double proton transfer in porphycene by path-integral molecular dynamics simulations.</a> Chemical Physics Letters. 496 (1-3), 14-19 (2010).</li>
<li>Smedarchina, Z., Shibl, M. F., Kühn, O., Fernández-Ramos, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+tautomerization+dynamics+of+porphycene+and+its+isotopomers+-+concerted+versus+stepwise+mechanisms%5BTitle%5D)+AND+%22Chemical+Physics+Letters%22%5BJournal%5D)">The tautomerization dynamics of porphycene and its isotopomers - concerted versus stepwise mechanisms.</a> Chemical Physics Letters. 436 (4-6), 314-321 (2007).</li>
<li>Wolfsberg, M., Hook, W. A., Paneth, P., Rebelo, L. P. N. Isotope effects. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aWolfsberg%2C+M+%22The+Chemical%2C+Geological%2C+and+Bio+Sciences%22">The Chemical, Geological, and Bio Sciences</a>. , Springer. Dordrecht. (2009).</li>
<li>Pietrzak, M., Shibl, M. F., Bröring, M., Kühn, O., Limbach, H. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((1H%2F2H+NMR+studies+of+geometric+H%2FD+isotope+effects+on+the+coupled+hydrogen+bonds+in+porphycene+derivatives%5BTitle%5D)+AND+%22Journal+of+the+American+Chemical+Society%22%5BJournal%5D)">1H/2H NMR studies of geometric H/D isotope effects on the coupled hydrogen bonds in porphycene derivatives.</a> Journal of the American Chemical Society. 129 (2), 296-304 (2007).</li>
<li>Yoshikawa, T., Sugawara, S., Takayanagi, T., Shiga, M., Tachikawa, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Quantum+tautomerization+in+porphycene+and+its+isotopomers%3A+Path-integral+molecular+dynamics+simulations%5BTitle%5D)+AND+%22Chemical+Physics%22%5BJournal%5D)">Quantum tautomerization in porphycene and its isotopomers: Path-integral molecular dynamics simulations.</a> Chemical Physics. 394 (1), 46-51 (2012).</li>
<li>Ladenthin, J. N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Force-induced+tautomerization+in+a+single+molecule%5BTitle%5D)+AND+%22Nature+Chemistry%22%5BJournal%5D)">Force-induced tautomerization in a single molecule.</a> Nature Chemistry. 8 (10), 935-940 (2016).</li>
<li>Dellago, C., Bolhuis, P. G., Csajka, F. S., Chandler, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transition+path+sampling+and+the+calculation+of+rate+constants%5BTitle%5D)+AND+%22The+Journal+of+Chemical+Physics%22%5BJournal%5D)">Transition path sampling and the calculation of rate constants.</a> The Journal of Chemical Physics. 108 (5), 1964-1977 (1998).</li>
<li>Bolhuis, P. G., Chandler, D., Dellago, C., Geissler, P. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transition+path+sampling%3A+Throwing+ropes+over+rough+mountain+passes%2C+in+the+dark%5BTitle%5D)+AND+%22Annual+Review+of+Physical+Chemistry%22%5BJournal%5D)">Transition path sampling: Throwing ropes over rough mountain passes, in the dark.</a> Annual Review of Physical Chemistry. 53 (1), 291-318 (2002).</li>
<li>Torrie, G. M., Valleau, J. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Nonphysical+sampling+distributions+in+monte+carlo+free-energy+estimation%3A+Umbrella+sampling%5BTitle%5D)+AND+%22Journal+of+Computational+Physics%22%5BJournal%5D)">Nonphysical sampling distributions in monte carlo free-energy estimation: Umbrella sampling.</a> Journal of Computational Physics. 23 (2), 187-199 (1977).</li>
<li>Kästner, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Umbrella+sampling%5BTitle%5D)+AND+%22Wiley+Interdisciplinary+Reviews%3A+Computational+Molecular+Science%22%5BJournal%5D)">Umbrella sampling.</a> Wiley Interdisciplinary Reviews: Computational Molecular Science. 1 (6), 932-942 (2011).</li>
<li>Gao, Y. Q. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+integrate-over-temperature+approach+for+enhanced+sampling%5BTitle%5D)+AND+%22Journal+of+Chemical+Physics%22%5BJournal%5D)">An integrate-over-temperature approach for enhanced sampling.</a> Journal of Chemical Physics. 128 (6), 064105(2008).</li>
<li>Yang, L., Gao, Y. Q. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+selective+integrated+tempering+method%5BTitle%5D)+AND+%22Journal+of+Chemical+Physics%22%5BJournal%5D)">A selective integrated tempering method.</a> Journal of Chemical Physics. 131 (21), 214109(2009).</li>
<li>Yang, M., Yang, L., Gao, Y., Hu, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Combine+umbrella+sampling+with+integrated+tempering+method+for+efficient+and+accurate+calculation+of+free+energy+changes+of+complex+energy+surface%5BTitle%5D)+AND+%22The+Journal+of+Chemical+Physics%22%5BJournal%5D)">Combine umbrella sampling with integrated tempering method for efficient and accurate calculation of free energy changes of complex energy surface.</a> The Journal of Chemical Physics. 141 (4), 044108(2014).</li>
<li>Yang, Y. I., Zhang, J., Che, X., Yang, L., Gao, Y. Q. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Efficient+sampling+over+rough+energy+landscapes+with+high+barriers%3A+A+combination+of+metadynamics+with+integrated+tempering+sampling%5BTitle%5D)+AND+%22The+Journal+of+Chemical+Physics%22%5BJournal%5D)">Efficient sampling over rough energy landscapes with high barriers: A combination of metadynamics with integrated tempering sampling.</a> The Journal of Chemical Physics. 144 (9), 094105(2016).</li>
<li>Xie, L., Shen, L., Chen, Z. N., Yang, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Efficient+free+energy+calculations+by+combining+two+complementary+tempering+sampling+methods%5BTitle%5D)+AND+%22The+Journal+of+Chemical+Physics%22%5BJournal%5D)">Efficient free energy calculations by combining two complementary tempering sampling methods.</a> The Journal of Chemical Physics. 146 (2), 024103(2017).</li>
<li>Xie, L., Cheng, H., Fang, D., Chen, Z. N., Yang, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enhanced+QM%2FMM+sampling+for+free+energy+calculation+of+chemical+reactions%3A+A+case+study+of+double+proton+transfer%5BTitle%5D)+AND+%22The+Journal+of+Chemical+Physics%22%5BJournal%5D)">Enhanced QM/MM sampling for free energy calculation of chemical reactions: A case study of double proton transfer.</a> The Journal of Chemical Physics. 150 (4), 044111(2019).</li>
<li>Hu, X., Hu, H., Yang, W. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aHu%2C+X+%22QM4D%3A+an+integrated+and+versatile+quantum+mechanical%2Fmolecular+mechanical+simulation+package+%28http%3A%2F%2Fwww.qm4d.info%2F%29%22">QM4D: an integrated and versatile quantum mechanical/molecular mechanical simulation package (http://www.qm4d.info/)</a>. , (2016).</li>
<li>Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., Klein, M. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comparison+of+simple+potential+functions+for+simulating+liquid+water%5BTitle%5D)+AND+%22The+Journal+of+Chemical+Physics%22%5BJournal%5D)">Comparison of simple potential functions for simulating liquid water.</a> The Journal of Chemical Physics. 79 (2), 926-935 (1983).</li>
<li>Walewski, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Scc-dftb+energy+barriers+for+single+and+double+proton+transfer+processes+in+the+model+molecular+systems+malonaldehyde+and+porphycene%5BTitle%5D)+AND+%22International+Journal+of+Quantum+Chemistry%22%5BJournal%5D)">Scc-dftb energy barriers for single and double proton transfer processes in the model molecular systems malonaldehyde and porphycene.</a> International Journal of Quantum Chemistry. 106 (3), 636-640 (2006).</li>
<li>Bakowies, D., Thiel, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Hybrid+models+for+combined+quantum+mechanical+and+molecular+mechanical+approaches%5BTitle%5D)+AND+%22The+Journal+of+Physical+Chemistry%22%5BJournal%5D)">Hybrid models for combined quantum mechanical and molecular mechanical approaches.</a> The Journal of Physical Chemistry. 100 (25), 10580-10594 (1996).</li>
<li>Hu, H., Yang, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Development+and+application+of+ab+initio+QM%2FMM+methods+for+mechanistic+simulation+of+reactions+in+solution+and+in+enzymes%5BTitle%5D)+AND+%22Journal+of+Molecular+Structure%3A+THEOCHEM%22%5BJournal%5D)">Development and application of ab initio QM/MM methods for mechanistic simulation of reactions in solution and in enzymes.</a> Journal of Molecular Structure: THEOCHEM. 898 (1-3), 17-30 (2009).</li>
<li>Xie, L., Yang, M., Chen, Z. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Understanding+the+entropic+effect+in+chorismate+mutase+reaction+catalyzed+by+isochorismate-pyruvate+lyase+from+pseudomonas+aeruginosa+%28PchB%29%5BTitle%5D)+AND+%22Catalysis+Science%2C+Technology%22%5BJournal%5D)">Understanding the entropic effect in chorismate mutase reaction catalyzed by isochorismate-pyruvate lyase from pseudomonas aeruginosa (PchB).</a> Catalysis Science, Technology. 9 (4), 957-965 (2019).</li>
<li>Shibl, M. F., Pietrzak, M., Limbach, H. H., Kühn, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Geometric+H%2FD+isotope+effects+and+cooperativity+of+the+hydrogen+bonds+in+porphycene%5BTitle%5D)+AND+%22ChemPhysChem%22%5BJournal%5D)">Geometric H/D isotope effects and cooperativity of the hydrogen bonds in porphycene.</a> ChemPhysChem. 8 (2), 315-321 (2007).</li>
</ol>

Autorzy nie mają nic do ujawnienia.

Strukturę porficenu przedstawiono na rysunku 1. Do opisu reakcji chemicznych zachodzących w wodzie wykorzystano elektrostatyczny potencjał hybrydowy QM/MM metodą SITS 23,24. Transfer protonów zachodzi w obrębie porfirynu3, a zatem porficen jest ustawiony jako region QM, a woda przypominająca jest ustawiona jako obszar MM. W tym przypadku przyjęliśmy DFTB/MIO jako naszą metodę zarządzania jakością do obróbki po...

Proces przeniesienia protonów w porfycenach posiada potencjalne zastosowania w opracowywaniu przełączników molekularnych, tranzystorów i urządzeń do przechowywania informacji1,2. W szczególności tautomeryzacja w porfycenach poprzez proces podwójnego przeniesienia protonów wzbudziła szerokie zainteresowanie w dziedzinach spektroskopii i fotofizyki2. Wewnętrzne atomy wodoru w porfycenie mogą migrować z jednego izomeru trans do drugiego równoważnego izomeru trans poprzez proces podwójnego przeniesienia protonów, co przedstawiono na Rysunku 1. Zaproponowano dwa mechanizmy procesu podwójnego przeniesienia protonów: mechanizm skoordynowany oraz mechanizm stopniowy3,4. W skoordynowanym procesie podwójnego przeniesienia protonów oba atomy protonów przemieszczają się do stanu przejściowego synchronicznie w sposób symetryczny, podczas gdy w procesie stopniowym jeden proton kończy przeniesienie przed drugim protonem. Dwa atomy wodoru mogą przenosić się jednocześnie lub stopniowo, w zależności od siły korelacji między dwoma atomami wodoru5.
Podstawienie izotopowe było wykorzystywane do wykrywania właściwości strukturalnych cząsteczek oraz stałych szybkości kinetyki reakcji6. Pojedyncza substytucja deuterem wewnętrznego wodoru w porfirycenie prowadzi do asymetrycznego kształtu cząsteczki. Wiązanie wodorowe może ulec wydłużeniu lub skróceniu ze względu na różnicę mas między atomami wodoru i deuteru. Podstawienie izotopowe wprowadza zaburzenie w szkielecie porfirycenu. Pojawia się zatem pytanie, czy asymetryczna struktura wpłynie na proces przeniesienia protonu. Limbach i współpracownicy donieśli, że zastąpienie wodoru deuterem spowoduje skrócenie obu wiązań wodorowych, a sprzężenie kooperatywne dwóch wiązań wodorowych w porfirycenie może sprzyjać mechanizmowi skoordynowanemu7, podczas gdy Yoshikawa stwierdził, że deuteracja sprawi, iż mechanizm stopniowy będzie w większym stopniu przyczyniał się do reakcji niż mechanizm skoordynowany8. Opracowano techniki eksperymentalne, takie jak spektroskopia sił, aby uchwycić szczegóły tautomeryzacji w pojedynczej cząsteczce porfirycenu9. Jednakże, ze względu na przejściowy charakter tego procesu, eksperymentalne określenie szczegółów atomowych przeniesienia protonu wciąż stanowi wyzwanie.
Obliczenia teoretyczne i symulacje mogą służyć jako narzędzia komplementarne w wyjaśnianiu mechanizmów reakcji przeniesienia protonu. Spośród różnych metod teoretycznych symulacje dynamiki molekularnej (MD) pozwalają na monitorowanie ruchów dynamicznych każdego atomu i są szeroko stosowane do ujawniania złożonych mechanizmów w reakcjach chemicznych i enzymatycznych. Jednak konwencjonalne symulacje MD mają tendencję do występowania problemu niewystarczającego próbkowania, szczególnie gdy w analizowanym procesie istnieje wysoka bariera energetyczna. W związku z tym opracowano metody próbkowania wzmocnionego, które obejmują próbkowanie ścieżki przejścia10,11, próbkowanie parasolowe (US)12,13 oraz zintegrowane próbkowanie temperowania (ITS)14,15. Połączenie różnych metod próbkowania wzmocnionego może dodatkowo zwiększyć wydajność próbkowania16,17,18. Aby wykorzystać algorytmy próbkowania wzmocnionego w symulowaniu reakcji chemicznych, zaimplementowaliśmy niedawno metodę selektywnego zintegrowanego próbkowania temperowania (SITS) z potencjałami mechaniki kwantowej i mechaniki molekularnej (QM/MM)19. Zaproponowana metoda SITS-QM/MM łączy zalety obu podejść: metoda SITS przyspiesza próbkowanie i pozwala na eksplorację wszystkich możliwych kanałów reakcji bez wcześniejszej wiedzy o mechanizmie reakcji, a QM/MM zapewnia dokładniejszy opis procesów tworzenia i zrywania wiązań, których nie można symulować wyłącznie metodami MM. Zaimplementowane podejście SITS-QM/MM pozwoliło z sukcesem odkryć skoordynowany podwójny transfer protonów oraz nieskorelowany i skorelowany stopniowy mechanizm podwójnego transferu protonów w różnych układach, bez uprzedniego definiowania współrzędnych reakcji19. W przypadku porfyceny zgłoszono charakter stopniowego, lecz skorelowanego przeniesienia protonów19. W naszym badaniu hybrydowa metoda SITS-QM/MM została wykorzystana do zbadania efektu izotopowego w porfycenie; poniżej znajdują się szczegółowe opisy algorytmu i protokołu naszej metody.
Wdrożyliśmy metodę SITS z zastosowaniem hybrydowych potencjałów QM/MM. Potencjał efektywny SITS zdefiniowano tak, aby obejmował energię potencjalną w różnych temperaturach z czynnikami wagowymi nk, co pozwala na pokrycie szerszych zakresów temperatur,
<img alt="Wzór mechaniki statystycznej, U'(R), minimalizacja entropii, równanie do analizy rozkładu energii." src="/files/ftp_upload/60040/60040equ01.jpg" data-altupdatedat="1788282455000" data-alt="Równanie 1">
gdzie N to liczba członów kanonicznych, βk to temperatura odwrócona, a nk to odpowiadający czynnik wagowy dla każdego składnika kanonicznego. UE(R) oraz UN(R) reprezentują odpowiednio człony wzmocnione i niewzmocnione w SITS i są zdefiniowane jako,
<img alt="Równania potencjału elektrostatycznego; UE(R)=US+1/2Use, UN(R)=Ue+1/2Use; wzór matematyczny." src="/files/ftp_upload/60040/60040equ02.jpg" data-altupdatedat="1788282455000" data-alt="Równanie 2">
Us, Use oraz Ue to energia potencjalna podsystemu, oddziaływanie między podsystemem a otoczeniem oraz energia potencjalna otoczenia. Potencjał QM/MM jest wyrażony jako hybrydowa suma trzech komponentów,
<img alt="Równanie dla rozdziału energii potencjalnej: $ U = U_{qm} + U_{qm/mm} + U_{mm} $." src="/files/ftp_upload/60040/60040equ03.jpg" data-altupdatedat="1788282455000" data-alt="Równanie 3">
gdzie Uqm, Uqm/mm oraz Umm oznaczają odpowiednio składnik energii wewnętrznej podsystemu QM, energię oddziaływania między regionami QM i MM oraz energię oddziaływań wewnątrz podsystemu MM. Składnik Uqm/mm można dalej podzielić na trzy komponenty, obejmujące energię oddziaływań elektrostatycznych, van der Waalsa oraz kowalencyjnych między atomami QM i MM,
<img alt="Równanie mechaniki kwantowej/mechaniki molekularnej (QM/MM); komponenty energii; chemia teoretyczna." src="/files/ftp_upload/60040/60040equ04.jpg" data-altupdatedat="1788282455000" data-alt="Równanie 4">
W SITS przypisujemy <img alt="Wzór na energię potencjalną w mechanice kwantowej $U_{qm}$." src="/files/ftp_upload/60040/60040equ05.jpg" data-altupdatedat="1788282455000" data-alt="Equation 5"> oraz <img alt="Mechanika kwantowa i mechanika molekularna, symbol U^ele_qm/mm, stosowany w chemii obliczeniowej." src="/files/ftp_upload/60040/60040equ06.jpg" data-altupdatedat="1788282455000" data-alt="Equation 6"> do jednego składnika Us,
<img alt="Sformułowanie mechaniki kwantowej, równania: Us=Uqm+Uqme/qm/mm, schemat analizy oddziaływań cząsteczkowych." src="/files/ftp_upload/60040/60040equ07.jpg" data-altupdatedat="1788282455000" data-alt="Równanie 7">
Pełny potencjał układu został następnie rozłożony na energię podukładu Us, energię oddziaływania między podukładem a otoczeniem Use oraz energię otoczenia Ue. Przykładowo, w układzie opisanym w niniejszej pracy podukładem jest porfyrzyna, a otoczeniem woda.
Profil PMF wzdłuż zmiennej zbiorczej τ(R) wyznacza się jako:
<img alt="Równanie potencjału średniej siły, PMF(τ), schemat dla analizy mechaniki statystycznej." src="/files/ftp_upload/60040/60040equ08.jpg" data-altupdatedat="1788282455000" data-alt="Równanie 8">
Ogólnie stosowane współrzędne reakcji dla każdego przeniesienia wodoru w układzie N1−H1··· N2 to q1 = (r1−r2)/2 oraz q2 = r1 + r2, gdzie r1 jest odległością N1-H1, a r2 jest odległością H1-N2.
Metoda ta została zaimplementowana w pakiecie do symulacji MD w modelu QM/MM QM4D20. Pełny kod źródłowy oraz dokumentacja są dostępne pod adresem: http://www.qm4d.info Simply.
Zazwyczaj symulacje SITS-QM/MM MD obejmują cztery etapy: stan przedrównoważny (pre-sits), optymalizację nk (opt-sits), symulację produkcyjną oraz analizę danych.

<table><tbody><tr><td>system operacyjny</td><td>CentOS Linux release 6.0</td><td></td><td></td></tr><tr><td>Oprogramowanie QM4D</td><td>http://www.qm4d.info/</td><td></td> <td>program własny</td></tr><tr><td>Komputer stacjonarny</td><td>HP</td><td></td><td></td></tr></tbody></table>

wpływ izotopowy porficenu w procesie podwójnego transferu protonów badany metodą enhanced qm/mm

Isotopic Effect in Double Proton Transfer Process of Porphycene Investigated by Enhanced QM/MM Method

Pojedyncze podstawienie deuteru w porficenie prowadzi do asymetrycznej geometrii molekularnej, co może wpływać na proces podwójnej transferu protonów w cząsteczce porficenu. W tym badaniu zastosowaliśmy ulepszoną metodę QM/MM o nazwie SITS-QM/MM w celu zbadania wpływu izotopów wodoru/deuteru (H/D) na transfer podwójnych protonów w porficenie. Zmiany odległości w symulacjach dynamiki molekularnej SITS-QM/MM sugerowały, że porficen podstawiony deuterem przyjął stopniowy mechanizm podwójnego transferu protonów. Analiza strukturalna i przesunięcia energii swobodnej procesu transferu podwójnych protonów wykazały, że asymetryczna substytucja izotopowa subtelnie ściska kowalencyjne wiązania wodorowe i może zmienić pierwotne położenie stanu przejściowego.

W niniejszym protokole zbadano wpływ pojedynczej substytucji deuteru na proces podwójnego przeniesienia protonu w porfirycenie (Rycina 1). Energię potencjalną podsystemu QM oraz wody podczas etapu przedrównoważnego i optymalizacji sprawdzono, aby upewnić się, że zakres energii został rozszerzony do szerszego przedziału (Rycina 2). Reprezentatywne zmiany odległości i kątów (Rycina 3 i Rycina 4) oraz rzut...

Protokół wykorzystujący ulepszoną metodę QM/MM do zbadania wpływu izotopowego na proces podwójnego transferu protonów w porficenie jest przedstawiony tutaj.

Wpływ izotopowy porficenu w procesie podwójnego transferu protonów badany metodą Enhanced QM/MM

The single deuterium substitution in porphycene leads to an asymmetric molecular geometry, which may affect the double proton transfer process in the ...

Ulepszona metoda QM/MM pozwala przezwyciężyć proste problemy w normalnej symulacji QM/MM. Metoda Six QM/MM przyspiesza wybór montażu dla regionu QM i może uchwycić ścieżki reakcji chemicznej, co skutkuje zdefiniowaniem współrzędnych reakcji. Dzięki temu protokołowi udało nam się uchwycić szlaki reakcji chemicznej podwójnej wymiany protonów i ujawnić wpływ substytucji deuteru na mechanizm transferu prowansencji w wodzie.Protokół może być wykorzystany do zbadania substytucji halogenów lub deuteru w identyfikacji ciepła w odkrywaniu leków. Główną zaletą metody Six QM / MM jest to, że nie potrzebujemy jej do zdefiniowania współrzędnej reakcji solanki ani wprowadzenia urządzenia do szlaku reakcji chemicznej podczas badania mechanizmu reakcji. pozwalają nam zidentyfikować możliwe ścieżki reakcji, które pochodzą z React.Metoda, która może być stosowana i rozszerzona do metody zarządzania jakością wysokiego poziomu i może stać się ważnym narzędziem do badania mechanizmu reakcji chemicznej w roztworze. Aby rozpocząć tę procedurę, należy zainicjować ustawienia wstępne, ustawiając typ pracy jako 100, temp0 jako 300, templow jako 260, temphigh jako 1300 i krok jako 120 000 w pliku wejściowym.

isotopic-effect-double-proton-transfer-process-porphycene

Research

JoVE Journal

Chemistry

Chemia

Ten film opisuje technikę spektroskopii magnetycznego rezonansu jądrowego do badania interakcji białko-białko pomiędzy 15znakowanymi N dzikimi i zmutowanymi białkami envoplakiny a nieznakowanym białkiem wimentyny. Udana interakcja między envoplakiną typu dzikiego a wimentyną prowadzi do rozległego poszerzenia linii i zaniku pików w widmach NMR, podczas gdy brak interakcji między zmutowaną envoplakiną a wimentyną skutkuje dobrze rozdzielonymi pikami w widmach NMR.

Magnetyczny rezonans jądrowy do badania oddziaływań białko-białko na poziomie atomowym

Kluczowe etapy funkcjonowania białka, w szczególności zmiany konformacyjne szkieletu i reakcje transferu protonów, często odbywają się w skali od mikrosekund do milisekund. Te dynamiczne procesy można badać za pomocą spektroskopii w podczerwieni z transformacją Fouriera ze skanowaniem czasowo-rozdzielczym, w szczególności w przypadku białek, których funkcja jest wyzwalana przez światło.

Transfer protonów i dynamika konformacji białek w białkach światłoczułych za pomocą spektroskopii w podczerwieni z transformacją Fouriera ze skanowaniem czasowym i transformacją Fouriera

Manuskrypt opisuje metodę quasi-rezonansowej spektroskopii fluorescencyjnej wspomaganej fononem, która obejmuje zarówno spektroskopię o ograniczonej rozdzielczości laserowej, jak i spektroskopię fotoluminescencji (PL). Metoda ta wykorzystuje fonony optyczne do zapewnienia widm o ograniczonej szerokości linii struktur półprzewodnikowych podobnych do atomów w domenie energetycznej. Metodę tę można również łatwo zrealizować za pomocą pojedynczej konfiguracji spektrometru spektroskopii optycznej.

Spektroskopia fluorescencyjna quasi-rezonansowa wspomagana fononami o wysokiej rozdzielczości

Przedstawiono protokół monitorowania in situ dyfuzji cząsteczek gości w porowatych ośrodkach za pomocą obrazowania elektronowego rezonansu paramagnetycznego (EPR).

In Situ Monitorowanie dyfuzji cząsteczek gościnnych w ośrodkach porowatych za pomocą obrazowania elektronowego rezonansu paramagnetycznego

Przedstawiamy protokół do badania ultraszybkich koherencji wibracyjnych w wieloatomowych kationach rodników, które skutkują dysocjacją molekularną.

Pomiar ultraszybkich koherencji wibracyjnych w wieloatomowych kationach rodnikowych z jonizacją adiabatyczną w silnym polu

Wydajność kwantowa fotoizomeryzacji jest podstawową właściwością fotofizyczną, która powinna być dokładnie określona w badaniach nad nowo opracowanymi fotoprzełącznikami. W tym miejscu opisujemy zestaw procedur pomiaru wydajności kwantowej fotoizomeryzacji hydrazonu fotochromowego jako modelowego bistabilnego fotoprzełącznika.

Wyznaczanie wydajności kwantowej fotoizomeryzacji fotoprzełącznika hydrazonowego

Tutaj używamy stabilizatora polimerowego do przygotowania zawiesin struktur metaloorganicznych (MOF), które wykazują znacznie zmniejszone rozproszenie w swoich widmach absorpcji w stanie podstawowym i przejściowym. W przypadku tych zawieszeń MOF, protokół zawiera różne wytyczne do spektroskopowego charakteryzowania MOF w celu uzyskania danych możliwych do interpretacji.

Przewodnik techniczny dotyczący wykonywania pomiarów spektroskopowych na strukturach metaloorganicznych

Prezentujemy technikę pomiaru z dużą precyzją proporcji izotopów w narządach myszy.

Precyzyjne pomiary izotopowe stosowane do narządów myszy

W tym artykule opisujemy międzyfazową syntezę sprzężonych mikroporowatych polimerów (CMP) na podłożach protektorowych oraz rozpuszczanie substratu do przygotowania wolnostojących nanomembran CMP. Ponadto opiszemy, w jaki sposób kruche nanomembrany mogą być przenoszone na inne podłoża.

Synteza i transfer warstwa po warstwie wolnostojących mikroporowatych nanomembran polimerowych

Dwuwymiarowy materiał modelowy z dyskotycznego fosforanu cyrkonu został opracowany. Nieorganiczny kryształ o strukturze blaszkowej zsyntetyzowano metodami hydrotermalnymi, refluksowymi i mikrofalowymi. Podczas eksfoliacji cząsteczkami organicznymi warstwowe kryształy mogą zostać przekształcone w monowarstwy, a nematyczna faza ciekłokrystaliczna została utworzona przy wystarczającym stężeniu monowarstw.

Synteza i złuszczanie dyskotycznych fosforanów cyrkonu w celu otrzymania koloidalnych ciekłych kryształów

Poznaj metodę Six QM/MM do badania izotopowych efektów w podwójnym transferze protonów w porficenie. Protokół umożliwia analizę mechanizmów reakcji bez predefiniowania współrzędnych reakcji.

Wpływ izotopowy porficenu w procesie podwójnego transferu protonów badany metodą Enhanced QM/MM

Przeglądaj więcej filmów