{"id":2421,"date":"2023-09-19T16:11:55","date_gmt":"2023-09-19T20:11:55","guid":{"rendered":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/?page_id=2421"},"modified":"2025-02-14T12:51:05","modified_gmt":"2025-02-14T17:51:05","slug":"examples","status":"publish","type":"page","link":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/resources\/gamd\/examples\/","title":{"rendered":"Examples"},"content":{"rendered":"<h2>Examples of LiGaMD simulations<\/h2>\n<h3>Trypsin-Benzamidine Binding:<\/h3>\n<p class=\"body\" align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-2402\" src=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/LiGaMD-TOC-1024x442.jpg\" alt=\"Examples of LiGaMD simulations: Trypsin-Benzamidine Binding\" width=\"750\" height=\"324\" srcset=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/LiGaMD-TOC-1024x442.jpg 1024w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/LiGaMD-TOC-300x130.jpg 300w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/LiGaMD-TOC-768x332.jpg 768w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/LiGaMD-TOC-600x259.jpg 600w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/LiGaMD-TOC.jpg 1200w\" sizes=\"auto, (max-width: 750px) 100vw, 750px\" \/><\/p>\n<p class=\"body\" align=\"justify\">All-atom ligand Gaussian accelerated molecular dynamics (LiGaMD) simulations captured repeptitive binding of the benzamidine (BEN) ligand to the trypsin model protein, which enabled us to characterize the ligand binding free energy profiles, pathways and kinetic rate constants.<\/p>\n<h4 class=\"style-scope ytd-watch-metadata\"><strong>Simulation Movie: Trypsin-benzamidine Sim2<\/strong><\/h4>\n<div class=\"entry-content-asset\"><iframe loading=\"lazy\" title=\"Trypsin-benzamidine Sim2\" width=\"1140\" height=\"855\" src=\"https:\/\/www.youtube.com\/embed\/oYO42n49kL8?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe><\/div>\n<p>Ligand Gaussian accelerated Molecular Dynamcis (LiGaMD) simulation of repetitive dissociation and binding of a ligand molecule benzamidine to target protein trypsin.<\/p>\n<p>&nbsp;<\/p>\n<h4 class=\"style-scope ytd-watch-metadata\"><strong>Simulation Movie:\u00a0<\/strong><strong>Trypsin-benzamidine Sim5<\/strong><\/h4>\n<div class=\"entry-content-asset\"><iframe loading=\"lazy\" title=\"Trypsin-benzamidine Sim5\" width=\"1140\" height=\"855\" src=\"https:\/\/www.youtube.com\/embed\/6pWazFW7AhM?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe><\/div>\n<p>Ligand Gaussian accelerated Molecular Dynamcis (LiGaMD) simulation of repetitive dissociation and binding of a ligand molecule benzamidine to target protein trypsin.<\/p>\n<p class=\"body\" style=\"text-align: left\" align=\"justify\">Reference: Miao, Y.; A. Bhattarai; and J. Wang, <a href=\"https:\/\/doi.org\/10.1021\/acs.jctc.0c00395\">Ligand Gaussian accelerated molecular dynamics (LiGaMD): Characterization of ligand binding thermodynamics and kinetics.<\/a> <em>Journal of Chemical Theory and Computation<\/em>, <strong>16<\/strong>(9): 5526-5547.<\/p>\n<h2 style=\"text-align: left\" align=\"center\"><span class=\"title\"><strong>Examples of GaMD simulations<\/strong><\/span><\/h2>\n<h3 class=\"subtitle\"><strong>Protein folding: <\/strong><\/h3>\n<p class=\"body\" align=\"center\"><a href=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/chigFolding.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-2403\" src=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/chigFolding-300x94.jpg\" alt=\"chig Folding\" width=\"720\" height=\"226\" srcset=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/chigFolding-300x94.jpg 300w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/chigFolding-1024x322.jpg 1024w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/chigFolding-768x241.jpg 768w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/chigFolding-600x189.jpg 600w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/chigFolding.jpg 1200w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><\/a><\/p>\n<p class=\"body\" align=\"justify\">Folding of chignolin was captured in GaMD simulations: the unfolded (\u201cU\u201d), intermediate (\u201cI\u201d) and folded (\u201cF\u201d) states (blue) aligned to the PDB native structure (red). Residues including Tyr2, Asp3, Pro4, Glu5, Thr6, Thr8 and Trp9 are represented by sticks. Notably, Trp9 and Pro4 form hydrophobic interactions in the intermediate state.<\/p>\n<h4><strong>Simulation Movie: Chignolin Folding<\/strong><\/h4>\n<div class=\"entry-content-asset\"><iframe loading=\"lazy\" title=\"chignolin folding\" width=\"1140\" height=\"855\" src=\"https:\/\/www.youtube.com\/embed\/mHiHxCe_3bc?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe><\/div>\n<p class=\"body\" align=\"justify\">Reference: Miao, Y.; Feher, V. A.; McCammon, J. A., <a href=\"http:\/\/pubs.acs.org\/doi\/abs\/10.1021\/acs.jctc.5b00436\">Gaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling and Free Energy Calculation.<\/a> <em>J Chemical Theory and Computation,<\/em> <strong>11<\/strong>(8): 3584-3595.<\/p>\n<h3 class=\"subtitle\"><strong>Protein-ligand Binding: <\/strong><\/h3>\n<p class=\"body\" align=\"center\"><a href=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/T4l-ligandBinding-opt.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-2430\" src=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/T4l-ligandBinding-opt-300x271.jpg\" alt=\"T4l ligand Binding\" width=\"720\" height=\"650\" srcset=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/T4l-ligandBinding-opt-300x271.jpg 300w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/T4l-ligandBinding-opt-1024x924.jpg 1024w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/T4l-ligandBinding-opt-768x693.jpg 768w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/T4l-ligandBinding-opt-600x542.jpg 600w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/T4l-ligandBinding-opt.jpg 1200w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><\/a><\/p>\n<p class=\"body\" align=\"justify\">A pathway of benzene binding to the T4-lysozyme was observed during the GaMD simulation. The unbound (\u201cU\u201d), intermediate (\u201cI\u201d) and bound (\u201cB\u201d) poses of the protein-ligand complex (blue) with the protein C-terminal domain (residues 80-160) aligned to the PDB native structure (red). The protein and benzene are represented by ribbons and spheres, respectively, and they are colored by blue for the simulation structure while red for the PDB native structure, except that in the binding pathway the simulated benzene is represented by lines and colored by simulation time in a BWR color scale. Residues with heavy atoms found within 3 \u00c5 of benzene are represented by sticks.<\/p>\n<h4><strong>Simulation Movie: T4-lysozyme ligand binding<\/strong><\/h4>\n<div class=\"entry-content-asset\"><iframe loading=\"lazy\" title=\"T4-lysozyme ligand binding\" width=\"1140\" height=\"855\" src=\"https:\/\/www.youtube.com\/embed\/2HpfIRJOcdc?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe><\/div>\n<p style=\"text-align: left\" align=\"center\">GaMD simulation of ligand binding to T4 lysozyme<\/p>\n<p class=\"body\" style=\"text-align: left\" align=\"justify\">Reference: Miao, Y.; Feher, V. A.; McCammon, J. A., <a href=\"http:\/\/pubs.acs.org\/doi\/abs\/10.1021\/acs.jctc.5b00436\">Gaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling and Free Energy Calculation.<\/a> <em>J Chemical Theory and Computation,<\/em> <strong>11<\/strong>(8): 3584-3595.<\/p>\n<h3 class=\"subtitle\"><strong>Graded Activation and Free Energy Landscapes of a Muscarinic G-Protein-Coupled Receptor: <\/strong><\/h3>\n<p class=\"body\" align=\"justify\">G-protein\u2013coupled receptors (GPCRs) represent primary targets of about one-third of currently marketed drugs. The structure, dynamics, and function of GPCRs result from complex free energy landscapes. In this work, we have applied Gaussian accelerated molecular dynamics (GaMD) to study the ligand-dependent behavior of the M2 muscarinic GPCR. Extensive GaMD simulations have revealed distinct structural flexibility and free energy profiles that depict graded activation of the M2 receptor. We have captured both dissociation and binding of an orthosteric ligand in a single all-atom GPCR simulation. GaMD is well poised to study large biomolecules and ligand recognition for drug discovery.<\/p>\n<p class=\"body\" align=\"center\"><a href=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Graded-Activation-opt.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-2424\" src=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Graded-Activation-opt.jpg\" alt=\"GPCR Graded Activation\" width=\"720\" height=\"701\" srcset=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Graded-Activation-opt.jpg 1200w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Graded-Activation-opt-300x292.jpg 300w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Graded-Activation-opt-1024x998.jpg 1024w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Graded-Activation-opt-768x748.jpg 768w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Graded-Activation-opt-600x585.jpg 600w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><\/a><\/p>\n<p class=\"body\" align=\"justify\">The M2 receptor (ribbons) samples a large conformational space with significant structural rearrangements, especially for the TM6 helix. Binding of the inverse agonist QNB (green spheres) confines the receptor in the inactive state. Without the G protein or mimetic nanobody, the partial agonist ARC (yellow spheres) biases the M2 receptor to visit an intermediate state \u201cI1\u201d (orange ribbons). ARC is able to dissociate completely to the bulk solvent via the extracellular vestibule and rebinds to the receptor repeatedly during a 2030 ns GaMD simulation. In comparison, the full agonist IXO (red spheres) biases the receptor further, sampling both intermediate \u201cI1\u201d (orange ribbons) and \u201cI2\u201d (purple ribbons). IXO escapes out of the orthosteric pocket and visits the extracellular vestibule in one of the GaMD simulations. By adding the G protein mimetic nanobody (purple surface), the M2 receptor is stabilized in the fully active state (red ribbons) as bound by IXO or ARC, although ARC adopts two alternative conformations in the orthosteric pocket, ARC-P1 (yellow spheres) and ARC-P1&#8242; (cyan spheres).<\/p>\n<p class=\"body\" align=\"center\"><a href=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Ligand-Dissociation-and-Binding-opt.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-2425\" src=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Ligand-Dissociation-and-Binding-opt-195x300.jpg\" alt=\"GPCR Ligand Dissociation and Binding\" width=\"500\" height=\"767\" srcset=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Ligand-Dissociation-and-Binding-opt-195x300.jpg 195w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Ligand-Dissociation-and-Binding-opt-667x1024.jpg 667w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Ligand-Dissociation-and-Binding-opt-600x921.jpg 600w, https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-content\/uploads\/sites\/1385\/2023\/09\/GPCR-Ligand-Dissociation-and-Binding-opt.jpg 750w\" sizes=\"auto, (max-width: 500px) 100vw, 500px\" \/><\/a><\/p>\n<p class=\"body\" align=\"justify\">Pathways of partial agonist dissociation and binding observed in GaMD simulation. (A) Timecourse of the ARC?Asp1033.32 distance during 2030 ns simulation. Four dissociation and three binding events are labeled. (B-H) Schematic representations of the ligand pathways during (B) \u201cD1\u201d, (C) \u201cB1\u201d, (D) \u201cB2\u201d, (E) \u201cB3\u201d, (F) \u201cD2\u201d, (G) \u201cD3\u201d and (H) \u201cD4\u201d. The receptor is represented by blue ribbons and the ligand by sticks colored by the position along the membrane normal. (I) Ten lowest energy structural clusters of ARC that are labeled and colored in a GWR scale according to the PMF values.<\/p>\n<h4 class=\"style-scope ytd-watch-metadata\"><strong>Simulation Movie: M2-ARC dissociation and binding<\/strong><\/h4>\n<div class=\"entry-content-asset\"><iframe loading=\"lazy\" title=\"M2-ARC dissociation and binding\" width=\"1140\" height=\"855\" src=\"https:\/\/www.youtube.com\/embed\/971Ix2UaKRc?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe><\/div>\n<p class=\"body\" align=\"justify\">The dissociation and binding of partial agonist ARC observed in GaMD simulation of the M2-ARC complex system. Not only does ARC escape out of the receptor orthosteric pocket, but it also dissociates completely and rebinds multiple times to the receptor during the 2030 ns GaMD simulation.<\/p>\n<p class=\"body\" align=\"justify\">Reference: Y. Miao and J. A. McCammon (2016), <a href=\"http:\/\/www.pnas.org\/content\/early\/2016\/10\/06\/1614538113.abstract\" rel=\"noopener\">Graded activation and free energy landscapes of a muscarinic G protein-coupled receptor<\/a>. <em>Proc Natl Acad Sci U S A<\/em>, 113 (43): 12162\u201312167.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Examples of LiGaMD simulations Trypsin-Benzamidine Binding: All-atom ligand Gaussian accelerated molecular dynamics (LiGaMD) simulations captured repeptitive binding of the benzamidine (BEN) ligand to the trypsin model protein, which enabled us to characterize the ligand binding free energy profiles, pathways and kinetic rate constants. Simulation Movie: Trypsin-benzamidine Sim2 Ligand Gaussian accelerated Molecular Dynamcis (LiGaMD) simulation of &hellip; <a href=\"https:\/\/www.med.unc.edu\/pharm\/miaolab\/resources\/gamd\/examples\/\" aria-label=\"Read more about Examples\">Read more<\/a><\/p>\n","protected":false},"author":22429,"featured_media":0,"parent":2390,"menu_order":3,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_acf_changed":false,"layout":"","cellInformation":"","apiCallInformation":"","footnotes":"","_links_to":"","_links_to_target":""},"class_list":["post-2421","page","type-page","status-publish","hentry","odd"],"acf":[],"_links_to":[],"_links_to_target":[],"_links":{"self":[{"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/pages\/2421","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/users\/22429"}],"replies":[{"embeddable":true,"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/comments?post=2421"}],"version-history":[{"count":0,"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/pages\/2421\/revisions"}],"up":[{"embeddable":true,"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/pages\/2390"}],"wp:attachment":[{"href":"https:\/\/www.med.unc.edu\/pharm\/miaolab\/wp-json\/wp\/v2\/media?parent=2421"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}