{"title":"Free energy, rates, and mechanism of transmembrane dimerization in lipid bilayers from dynamically unbiased molecular dynamics simulations","authors":"Emil Jackel, Gianmarco Lazzeri, Roberto Covino","doi":"arxiv-2408.01407","DOIUrl":null,"url":null,"abstract":"The assembly of proteins in membranes plays a key role in many crucial\ncellular pathways. Despite their importance, characterizing transmembrane\nassembly remains challenging for experiments and simulations. Equilibrium\nmolecular dynamics simulations do not cover the time scales required to sample\nthe typical transmembrane assembly. Hence, most studies rely on enhanced\nsampling schemes that steer the dynamics of transmembrane proteins along a\ncollective variable that should encode all slow degrees of freedom. However,\ngiven the complexity of the condensed-phase lipid environment, this is far from\ntrivial, with the consequence that free energy profiles of dimerization can be\npoorly converged. Here, we introduce an alternative approach, which relies only\non simulating short, dynamically unbiased trajectory segments, avoiding using\ncollective variables or biasing forces. By merging all trajectories, we obtain\nfree energy profiles, rates, and mechanisms of transmembrane dimerization with\nthe same set of simulations. We showcase our algorithm by sampling the\nspontaneous association and dissociation of a transmembrane protein in a lipid\nbilayer, the popular coarse-grained Martini force field. Our algorithm\nrepresents a promising way to investigate assembly processes in biologically\nrelevant membranes, overcoming some of the challenges of conventional methods.","PeriodicalId":501022,"journal":{"name":"arXiv - QuanBio - Biomolecules","volume":"58 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - QuanBio - Biomolecules","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2408.01407","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
The assembly of proteins in membranes plays a key role in many crucial
cellular pathways. Despite their importance, characterizing transmembrane
assembly remains challenging for experiments and simulations. Equilibrium
molecular dynamics simulations do not cover the time scales required to sample
the typical transmembrane assembly. Hence, most studies rely on enhanced
sampling schemes that steer the dynamics of transmembrane proteins along a
collective variable that should encode all slow degrees of freedom. However,
given the complexity of the condensed-phase lipid environment, this is far from
trivial, with the consequence that free energy profiles of dimerization can be
poorly converged. Here, we introduce an alternative approach, which relies only
on simulating short, dynamically unbiased trajectory segments, avoiding using
collective variables or biasing forces. By merging all trajectories, we obtain
free energy profiles, rates, and mechanisms of transmembrane dimerization with
the same set of simulations. We showcase our algorithm by sampling the
spontaneous association and dissociation of a transmembrane protein in a lipid
bilayer, the popular coarse-grained Martini force field. Our algorithm
represents a promising way to investigate assembly processes in biologically
relevant membranes, overcoming some of the challenges of conventional methods.