메타다이나믹스
Metadynamics메타다이나믹스(MTD, METAD 또는 MetaD)는 컴퓨터 물리, 화학 및 생물학에서 사용되는 컴퓨터 시뮬레이션 방법입니다.이 값은 시스템의 에너지 지형 형태에 의해 에르고디시티가 방해받는 시스템의 자유 에너지 및 기타 상태 기능을 추정하는 데 사용됩니다.그것은 2002년[1] 알레산드로 라이오와 미켈레 파리넬로에 의해 처음 제안되었고, 보통 분자역학 시뮬레이션에 적용된다.MTD는 적응적으로 편향된 분자 역학,[2] 적응 반응 좌표력[3] 및 국소 고도 우산 [4]표본 추출과 같은 많은 최신 방법과 매우 유사하다.보다 최근에는 원래 메타다이나믹스와 성질이 좋은 메타다이나믹스가[5] 중요도 샘플링의 맥락에서 도출되었고 적응 바이어싱 전위 [6]설정의 특수한 경우로 나타났다.MTD는 Wang-Landau [7]샘플링과 관련이 있습니다.
서론
이 기술은 감압,[8] 터널링,[9] 타부 검색,[10] 국소 고도,[11] 구조 홍수,[12] Engkvist-Karlström[13] 및 Adaptive Biasing Force [14]방법을 포함한 많은 관련 방법을 기반으로 한다.
메타다이나믹스는 비공식적으로 "자유 에너지 웰을 계산 [15]모래로 채우는 것"으로 묘사되어 왔다.알고리즘에서는 시스템이 몇 가지 집합 변수(CV)로 기술될 수 있다고 가정합니다.시뮬레이션 중에 집합변수에 의해 결정된 공간에서의 시스템 위치를 계산하고 시스템의 실제 에너지 경관에 양의 가우스 전위를 가산한다.이렇게 하면 시스템이 이전 시점으로 돌아가는 것이 권장되지 않습니다.시뮬레이션의 진화 동안, 점점 더 많은 가우시안들이 요약되고, 따라서 시스템이 전체 에너지 지형을 탐색할 때까지 이전 단계로 돌아가는 것을 점점 더 주저하게 됩니다. 이 시점에서 수정된 자유 에너지는 집합 변수의 함수로서 일정하게 됩니다. 이것이 집합 변수의 이유입니다.심하게 변동하기 시작합니다.이 시점에서 에너지 지형은 모든 가우시안 합계의 반대로 회복될 수 있다.
정확도와 계산 비용 간의 비율을 최적화하기 위해 가우스 높이와 가우스 폭뿐만 아니라 두 가우스 함수의 추가 시간 간격도 조정됩니다.단순히 가우스의 크기를 변경함으로써 메타다이나믹스를 큰 가우시안(Gaussian)을 사용하여 에너지 경관의 대략적인 지도를 매우 빠르게 생성하거나 작은 가우시안([1]Gaussian)을 사용하여 보다 세밀한 설명을 위해 사용할 수 있습니다.일반적으로 성질이 좋은 메타다이나믹스를[5] 사용하여 가우스 크기를 적응적으로 변경합니다.또, 가우스폭은 적응형 가우스 메타나믹스로 [16]적응할 수 있다.
메타다이나믹스는 적응형 우산 표본 추출과 같은 방법에서 [1]탐색에 에너지 지형의 초기 추정치를 필요로 하지 않는다는 장점이 있다.그러나 복잡한 시뮬레이션을 위해 적절한 집합 변수를 선택하는 것은 간단한 일이 아닙니다.일반적으로 적절한 집합 변수 집합을 찾기 위해서는 몇 가지 시행이 필요하지만, 필수 좌표,[17] 스케치 [18]맵 및 비선형 데이터 중심 집합 [19]변수와 같은 몇 가지 자동 절차가 제안된다.
멀티 리플리카 어프로치
사용성과 병렬 성능을 개선하기 위해 독립적인 메타다이나믹 시뮬레이션(복제)을 함께 결합할 수 있다.이러한 방법에는 멀티 워커 MTD,[20] 병렬 감쇠 MTD,[21] 바이어스 교환 MTD [22]및 집합 가변 감쇠 MTD가 제안됩니다.[23]마지막 세 가지는 병렬 템퍼링 방법과 유사하며 복제 교환을 사용하여 샘플링을 개선합니다.일반적으로 복제품 교환에는 Metropolis-Hastings 알고리즘이 사용되지만 무한[24] 스와핑 및 Swa-Todo[25] 알고리즘이 복제품 [26]환율을 향상시킵니다.
고차원적 접근법
일반적인 (단일 리플리카) MTD 시뮬레이션에는 최대 3개의 CV가 포함될 수 있으며, 멀티 리플리카 접근방식을 사용하더라도 실제로는 8개의 CV를 초과하기 어렵습니다.이 제한은 가우스 함수(커넬)를 추가하여 구성된 바이어스 전위로부터 발생합니다.Kernel Density Estimator(KDE; 커널 밀도 추정기)의 특수한 경우입니다.일정한 KDE 정확도를 위해 필요한 커널 수는 차원 수에 따라 기하급수적으로 증가합니다.따라서 바이어스 전위의 동일한 정확도를 유지하기 위해 MTD 시뮬레이션 길이는 CV의 수에 따라 기하급수적으로 증가해야 합니다.또한 빠른 평가를 위한 바이어스 전위는 일반적으로 정규 [27]그리드를 사용하여 근사된다.그리드를 저장하는 데 필요한 메모리는 차원(CV)의 수와 함께 기하급수적으로 증가합니다.
메타다이나믹스의 고차원적 일반화는 NN2B이다.[28]이것은 두 가지 기계 학습 알고리즘, 즉 가장 가까운 이웃 밀도 추정기(NNDE)와 인공 신경 네트워크(ANN)를 기반으로 한다. NNDE는 KDE를 대체하여 짧은 편향 시뮬레이션에서 편향 전위의 업데이트를 추정하는 반면 ANN은 결과 편향 전위를 추정하는 데 사용된다.ANN은 고차원 함수의 메모리 효율적 표현으로, 역전파 [28][29]알고리즘에 의해 미분(바이어스 힘)이 효과적으로 계산됩니다.
적응 바이어스 전위를 위해 ANN을 이용하는 대체 방법은 추정을 [30]위해 평균 전위력을 사용한다.이 방법은 또한 적응 바이어싱 힘([31]ABF) 방법을 고차원적으로 일반화한 것입니다.또, 베이지안 [32]정규화를 이용해 ANN의 트레이닝을 개선해, ANN의 [30]앙상블을 트레이닝 하는 것으로 근사 오차를 추측할 수 있다.
알고리즘.
를 들어, { } {\ {i 1 의 위치에 있는 클래식N( 스타일N)의 시스템이 있다고 가정합니다. 좌표 3) {textstyle3입자 상호작용은 잠재적 VV ({ V V로 설명됩니다.잠재적 함수 형태(예: 고에너지 장벽으로 분리된 두 개의 국소 최소값)는 분자 역학 또는 몬테 카를로 방법을 사용한 에르고드 표본 추출을 방지한다.
오리지널 메타다이나믹스
MTD의 일반적인 개념은 샘플링된 상태의 재검색을 방지함으로써 시스템샘플링을 강화하는 것입니다.이는 시스템 H(\H)에 바이어스 V 바이어스(\text{를 증가시킴으로써 달성됩니다.
- + + 바이어스text{
바이어스 전위는 집합 변수 바이어스 바이어스 ( ){{의 함수이다. 집합 변수는 입자 위치 i의 함수이다 바이어스 전위는 속도(\에서 바이어스를 추가하여 지속적으로 업데이트됩니다. 여기서 s {는 시간 {\ t의 순간 집합 변수 값입니다.
- s= )= t( -→ t) { { \ { bias } ( \ { s } )= \ } ( \ vec { s } ) \ viag ( \{ } )
무한히 긴 시뮬레이션 sim {\에서 누적 바이어스 전위는 반대 기호(및 무관한 C {\C})를 가진 자유 에너지로 수렴됩니다.
효율적인 계산 구현을 위해 프로세스는 {\} 시간 간격( {\displaystyle \는 플로어 함수를 나타냄 {\ \delta} - 함수는 로컬화된 양의 커널 K(\ K로 대체됩니다.potential은 순간 집합 값 s j를 중심으로 하는 커널 함수의 합계가 됩니다
일반적으로 커널은 다차원 가우스 함수이며, 공분산 행렬에는 대각 비 0 요소만 있습니다.
- {\right}\
{\(\displaystyledisplaystyle\displaystyle 및displaystyle\는 priori로 결정되어 시뮬레이션 중에 일정하게 유지됩니다.
실행
아래는 분자역학(MD)에 기반한 MTD 베이스의 의사코드이며 여기서 { {\ 및{ v { {v는 N {\ N - system의 위치와 속도이다. 바이어스 는 n / { n /\ MD 마다 업데이트되며, 시스템 포스{ { {F {\에 기여합니다
{ {vec 및 { { 설정 {s}}}}}:=0 MD 단계: CV 계산: S :의{\n} MD 단계마다{{r 바이어스 잠재력:V (s) : = ( ) +- ( - 2s - s→ 2) { {bias s {c} {c} {} {c} {c} {c}:{s 원자력을 계산합니다.}}}} _ { { 및 {\ } {\style} {\}} } } } } } } }} }
자유 에너지 추정기
커널의 크기가 한정되어 있기 때문에 바이어스 전위가 평균값 주위에 변동합니다.바이어스 전위를 평균화함으로써 수렴 자유 에너지를 얻을 수 있다.평균은 집합 변수에 따른 움직임이 확산되는 부터 됩니다.
적용들
메타다이나믹스는 다음과 같은 연구에 사용되었습니다.
실장
뿌루마드
PLUMED는[39] 많은 MTD 알고리즘과 집합 변수를 구현하는 오픈 소스 라이브러리입니다.유연한 객체 지향[40][41] 설계로 여러 MD 프로그램(AMBER, GROMACS, LAMPS, NAMD, Quantum ESPLO, DL_POLY_4, CP2K 및 OpenMM)[42][43]과 인터페이스할 수 있습니다.
다른.
기타 MTD 실장은 Collective Variables Module(LAMPS, NAMD 및 GROMACS용), ORAC, CP2K [45]및 Desmond에 있습니다.
외부 링크
「 」를 참조해 주세요.
레퍼런스
- ^ a b c Laio, A.; Parrinello, M. (2002). "Escaping free-energy minima". Proceedings of the National Academy of Sciences of the United States of America. 99 (20): 12562–12566. arXiv:cond-mat/0208352. Bibcode:2002PNAS...9912562L. doi:10.1073/pnas.202427399. PMC 130499. PMID 12271136.
- ^ Babin, V.; Roland, C.; Sagui, C. (2008). "Stabilization of resonance states by an asymptotic Coulomb potential". J. Chem. Phys. 128 (2): 134101/1–134101/7. Bibcode:2008JChPh.128b4101A. doi:10.1063/1.2821102. PMID 18205437.
- ^ Barnett, C.B.; Naidoo, K.J. (2009). "Free Energies from Adaptive Reaction Coordinate Forces (FEARCF): An application to ring puckering". Mol. Phys. 107 (8): 1243–1250. Bibcode:2009MolPh.107.1243B. doi:10.1080/00268970902852608. S2CID 97930008.
- ^ Hansen, H.S.; Hünenberger, P.H. (2010). "Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water". J. Comput. Chem. 31 (1): 1–23. doi:10.1002/jcc.21253. PMID 19412904. S2CID 7367058.
- ^ a b Barducci, A.; Bussi, G.; Parrinello, M. (2008). "Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method". Physical Review Letters. 100 (2): 020603. arXiv:0803.3861. Bibcode:2008PhRvL.100b0603B. doi:10.1103/PhysRevLett.100.020603. PMID 18232845. S2CID 13690352.
- ^ Dickson, B.M. (2011). "Approaching a parameter-free metadynamics". Phys. Rev. E. 84 (3): 037701–037703. arXiv:1106.4994. Bibcode:2011PhRvE..84c7701D. doi:10.1103/PhysRevE.84.037701. PMID 22060542. S2CID 42243972.
- ^ Christoph Jungghans, Danny Perez 및 Thomas Vogel. "다음파 앙상블의 분자 역학: 왕-란다우 표본 추출의 동등성, 통계 온도 분자 역학 및 메타나믹스."화학이론과 계산 저널 10.5 (2014) : 1843-1847. doi : 10.1021/ct500077d
- ^ Crippen, Gordon M.; Scheraga, Harold A. (1969). "Minimization of polypeptide energy. 8. Application of the deflation technique to a dipeptide". Proceedings of the National Academy of Sciences. 64 (1): 42–49. Bibcode:1969PNAS...64...42C. doi:10.1073/pnas.64.1.42. PMC 286123. PMID 5263023.
- ^ Levy, A.V.; Montalvo, A. (1985). "The Tunneling Algorithm for the Global Minimization of Functions". SIAM J. Sci. Stat. Comput. 6: 15–29. doi:10.1137/0906002.
- ^ Glover, Fred (1989). "Tabu Search—Part I". ORSA Journal on Computing. 1 (3): 190–206. doi:10.1287/ijoc.1.3.190. S2CID 5617719.
- ^ Huber, T.; Torda, A.E.; van Gunsteren, W.F. (1994). "Local elevation: A method for improving the searching properties of molecular dynamics simulation". J. Comput.-Aided Mol. Des. 8 (6): 695–708. Bibcode:1994JCAMD...8..695H. CiteSeerX 10.1.1.65.9176. doi:10.1007/BF00124016. PMID 7738605. S2CID 15839136.
- ^ Grubmüller, H. (1995). "Predicting slow structural transitions in macromolecular systems: Conformational flooding". Phys. Rev. E. 52 (3): 2893–2906. Bibcode:1995PhRvE..52.2893G. doi:10.1103/PhysRevE.52.2893. hdl:11858/00-001M-0000-000E-CA15-8. PMID 9963736.
- ^ Engkvist, O.; Karlström, G. (1996). "A method to calculate the probability distribution for systems with large energy barriers". Chem. Phys. 213 (1): 63–76. Bibcode:1996CP....213...63E. doi:10.1016/S0301-0104(96)00247-9.
- ^ Darve, E.; Pohorille, A. (2001). "Calculating free energies using average force". J. Chem. Phys. 115 (20): 9169. Bibcode:2001JChPh.115.9169D. doi:10.1063/1.1410978. hdl:2060/20010090348. S2CID 5310339.
- ^ http://www.grs-sim.de/cms/upload/Carloni/Presentations/Marinelli.ppt[영구 데드링크]
- ^ Branduardi, Davide; Bussi, Giovanni; Parrinello, Michele (2012-06-04). "Metadynamics with Adaptive Gaussians". Journal of Chemical Theory and Computation. 8 (7): 2247–2254. arXiv:1205.4300. doi:10.1021/ct3002464. PMID 26588957. S2CID 20002793.
- ^ Spiwok, V.; Lipovová, P.; Králová, B. (2007). "Metadynamics in essential coordinates: free energy simulation of conformational changes". The Journal of Physical Chemistry B. 111 (12): 3073–3076. doi:10.1021/jp068587c. PMID 17388445.
- ^ Ceriotti, Michele; Tribello, Gareth A.; Parrinello, Michele (2013-02-22). "Demonstrating the Transferability and the Descriptive Power of Sketch-Map". Journal of Chemical Theory and Computation. 9 (3): 1521–1532. doi:10.1021/ct3010563. PMID 26587614.
- ^ Hashemian, Behrooz; Millán, Daniel; Arroyo, Marino (2013-12-07). "Modeling and enhanced sampling of molecular systems with smooth and nonlinear data-driven collective variables". The Journal of Chemical Physics. 139 (21): 214101. Bibcode:2013JChPh.139u4101H. doi:10.1063/1.4830403. hdl:2117/20940. ISSN 0021-9606. PMID 24320358.
- ^ Raiteri, Paolo; Laio, Alessandro; Gervasio, Francesco Luigi; Micheletti, Cristian; Parrinello, Michele (2005-10-28). "Efficient Reconstruction of Complex Free Energy Landscapes by Multiple Walkers Metadynamics †". The Journal of Physical Chemistry B. 110 (8): 3533–3539. doi:10.1021/jp054359r. PMID 16494409. S2CID 15595613.
- ^ Bussi, Giovanni; Gervasio, Francesco Luigi; Laio, Alessandro; Parrinello, Michele (October 2006). "Free-Energy Landscape for β Hairpin Folding from Combined Parallel Tempering and Metadynamics". Journal of the American Chemical Society. 128 (41): 13435–13441. doi:10.1021/ja062463w. PMID 17031956.
- ^ a b Piana, S.; Laio, A. (2007). "A bias-exchange approach to protein folding". The Journal of Physical Chemistry B. 111 (17): 4553–4559. doi:10.1021/jp067873l. hdl:20.500.11937/15651. PMID 17419610.
- ^ Gil-Ley, Alejandro; Bussi, Giovanni (2015-02-19). "Enhanced Conformational Sampling Using Replica Exchange with Collective-Variable Tempering". Journal of Chemical Theory and Computation. 11 (3): 1077–1085. arXiv:1502.02115. doi:10.1021/ct5009087. PMC 4364913. PMID 25838811.
- ^ Plattner, Nuria; Doll, J. D.; Dupuis, Paul; Wang, Hui; Liu, Yufei; Gubernatis, J. E. (2011-10-07). "An infinite swapping approach to the rare-event sampling problem". The Journal of Chemical Physics. 135 (13): 134111. arXiv:1106.6305. Bibcode:2011JChPh.135m4111P. doi:10.1063/1.3643325. ISSN 0021-9606. PMID 21992286. S2CID 40621592.
- ^ Suwa, Hidemaro (2010-01-01). "Markov Chain Monte Carlo Method without Detailed Balance". Physical Review Letters. 105 (12): 120603. arXiv:1007.2262. Bibcode:2010PhRvL.105l0603S. doi:10.1103/PhysRevLett.105.120603. PMID 20867621. S2CID 378333.
- ^ Galvelis, Raimondas; Sugita, Yuji (2015-07-15). "Replica state exchange metadynamics for improving the convergence of free energy estimates". Journal of Computational Chemistry. 36 (19): 1446–1455. doi:10.1002/jcc.23945. ISSN 1096-987X. PMID 25990969. S2CID 19101602.
- ^ "PLUMED: Metadynamics". plumed.github.io. Retrieved 2018-01-13.
- ^ a b Galvelis, Raimondas; Sugita, Yuji (2017-06-13). "Neural Network and Nearest Neighbor Algorithms for Enhancing Sampling of Molecular Dynamics". Journal of Chemical Theory and Computation. 13 (6): 2489–2500. doi:10.1021/acs.jctc.7b00188. ISSN 1549-9618. PMID 28437616.
- ^ Schneider, Elia; Dai, Luke; Topper, Robert Q.; Drechsel-Grau, Christof; Tuckerman, Mark E. (2017-10-11). "Stochastic Neural Network Approach for Learning High-Dimensional Free Energy Surfaces". Physical Review Letters. 119 (15): 150601. Bibcode:2017PhRvL.119o0601S. doi:10.1103/PhysRevLett.119.150601. PMID 29077427.
- ^ a b Zhang, Linfeng; Wang, Han; E, Weinan (2017-12-09). "Reinforced dynamics for enhanced sampling in large atomic and molecular systems. I. Basic Methodology". The Journal of Chemical Physics. 148 (12): 124113. arXiv:1712.03461. doi:10.1063/1.5019675. PMID 29604808. S2CID 4552400.
- ^ Comer, Jeffrey; Gumbart, James C.; Hénin, Jérôme; Lelièvre, Tony; Pohorille, Andrew; Chipot, Christophe (2015-01-22). "The Adaptive Biasing Force Method: Everything You Always Wanted To Know but Were Afraid To Ask". The Journal of Physical Chemistry B. 119 (3): 1129–1151. doi:10.1021/jp506633n. ISSN 1520-6106. PMC 4306294. PMID 25247823.
- ^ Sidky, Hythem; Whitmer, Jonathan K. (2017-12-07). "Learning Free Energy Landscapes Using Artificial Neural Networks". The Journal of Chemical Physics. 148 (10): 104111. arXiv:1712.02840. doi:10.1063/1.5018708. PMID 29544298. S2CID 3932640.
- ^ Ensing, B.; De Vivo, M.; Liu, Z.; Moore, P.; Klein, M. (2006). "Metadynamics as a tool for exploring free energy landscapes of chemical reactions". Accounts of Chemical Research. 39 (2): 73–81. doi:10.1021/ar040198i. PMID 16489726.
- ^ Gervasio, F.; Laio, A.; Parrinello, M. (2005). "Flexible docking in solution using metadynamics". Journal of the American Chemical Society. 127 (8): 2600–2607. doi:10.1021/ja0445950. PMID 15725015. S2CID 6304388.
- ^ Vargiu, A. V.; Ruggerone, P.; Magistrato, A.; Carloni, P. (2008). "Dissociation of minor groove binders from DNA: insights from metadynamics simulations". Nucleic Acids Research. 36 (18): 5910–5921. doi:10.1093/nar/gkn561. PMC 2566863. PMID 18801848.
- ^ Martoňák, R.; Laio, A.; Bernasconi, M.; Ceriani, C.; Raiteri, P.; Zipoli, F.; Parrinello, M. (2005). "Simulation of structural phase transitions by metadynamics". Zeitschrift für Kristallographie. 220 (5–6): 489. arXiv:cond-mat/0411559. Bibcode:2005ZK....220..489M. doi:10.1524/zkri.220.5.489.65078. S2CID 96851280.
- ^ Cruz, F.J.A.L.; de Pablo, J.J.; Mota, J.P.B. (2014), "Endohedral confinement of a DNA dodecamer onto pristine carbon nanotubes and the stability of the canonical B form", J. Chem. Phys., 140 (22): 225103, arXiv:1605.01317, Bibcode:2014JChPh.140v5103C, doi:10.1063/1.4881422, PMID 24929415, S2CID 15149133
- ^ Cruz, F.J.A.L.; Mota, J.P.B. (2016), "Conformational Thermodynamics of DNA Strands in Hydrophilic Nanopores", J. Phys. Chem. C, 120 (36): 20357–20367, doi:10.1021/acs.jpcc.6b06234
- ^ "PLUMED". www.plumed.org. Retrieved 2016-01-26.
- ^ Bonomi, Massimiliano; Branduardi, Davide; Bussi, Giovanni; Camilloni, Carlo; Provasi, Davide; Raiteri, Paolo; Donadio, Davide; Marinelli, Fabrizio; Pietrucci, Fabio (2009-10-01). "PLUMED: A portable plugin for free-energy calculations with molecular dynamics". Computer Physics Communications. 180 (10): 1961–1972. arXiv:0902.0874. Bibcode:2009CoPhC.180.1961B. doi:10.1016/j.cpc.2009.05.011. S2CID 4852774.
- ^ Tribello, Gareth A.; Bonomi, Massimiliano; Branduardi, Davide; Camilloni, Carlo; Bussi, Giovanni (2014-02-01). "PLUMED 2: New feathers for an old bird". Computer Physics Communications. 185 (2): 604–613. arXiv:1310.0980. Bibcode:2014CoPhC.185..604T. doi:10.1016/j.cpc.2013.09.018. S2CID 17904052.
- ^ "MD engines - PLUMED". www.plumed.org. Archived from the original on 2016-02-07. Retrieved 2016-01-26.
- ^ "howto:install_with_plumed [CP2K Open Source Molecular Dynamics ]". www.cp2k.org. Retrieved 2016-01-26.
- ^ Fiorin, Giacomo; Klein, Michael L.; Hénin, Jérôme (December 2013). "Using collective variables to drive molecular dynamics simulations". Molecular Physics. 111 (22–23): 3345–3362. Bibcode:2013MolPh.111.3345F. doi:10.1080/00268976.2013.813594. ISSN 0026-8976.
- ^ "Cp2K_Input / Motion / Free_Energy / Metadyn".