Estimating the Lipophobic Contributions in Model Membranes

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

The impact of bilayer composition on the dimerization properties of the Slg1 stress sensor TMD from a multiscale analysis

The impact of mutual interactions between the transmembrane domains of membrane proteins and lipids on bilayer properties has gained major attraction. Most simulation studies of membranes rely on the Martini force field, which has proven extremely helpful in providing molecular insights into realistic systems. Accordingly, an evaluation of the accuracy of the Martini force field is crucial to be able to correctly interpret the reported data. In this study, we combine atomistic and coarse-grained Martini simulations to investigate the properties of transmembrane domains (TMDs) in a model yeast membrane. The results show that the TMD binding state (monomeric and dimeric with positive or negative crossing angle) and the membrane composition significantly influence the properties around the TMDs and change TMD-TMD and TMD-lipid affinities. Furthermore, ergosterol (ERG) exhibits a strong affinity to TMD dimers. Importantly, the right-handed TMD dimer configuration is stabilized via TMD-TMD contacts by the addition of asymmetric anionic phosphatidylserine (PS). The coarse-grained simulations corroborate many of these findings, with two notable exceptions: a systematic overestimation of TMD-ERG interaction and lack of stabilization of the right-handed TMD dimers with the addition of PS.

Role of Cholesterol in Transmembrane Dimerization of the ErbB2 Growth Factor Receptor

The association of ErbB2 growth factor receptors is critical for cell growth and potentiates tumor proliferation in several cancer types. An important aspect in ErbB2 association is the role of lipids such as cholesterol, especially since their metabolism is often reprogrammed in cancer cells. Here, we have coupled metadynamics with coarse-grain simulations to identify cholesterol effects in the transmembrane dimerization of ErbB2 receptors. Overall, cholesterol interactions are observed with the receptor that directly tunes the association energetics. Several dimer conformations are identified both in the presence and absence of cholesterol, although the dimer regime appears to be more favorable in the presence of cholesterol. We observe an overall modulation of the underlying energy profile and the symmetric active and inactive conformational states are not distinguished in the presence of cholesterol. We show that cholesterol binds to the receptor transmembrane domain at a site (CRAC motif) that overlaps with the dimer interface (SmXXXSm motif). The competition between the transmembrane interactions and cholesterol interactions decides the final conformational landscape. Our work is an important step toward characterizing cholesterol effects in ErbB2 membrane receptor function.

Lipid-dependent conformational landscape of the ErbB2 growth factor receptor dimers

Altered lipid metabolism has been linked to cancer development and progression. Several roles have been attributed to the increased saturation and length of lipid acyl tails observed in tumors, but its effect on signaling receptors is still emerging. In this work, we have analyzed the lipid dependence of the ErbB2 growth factor receptor dimerization that plays an important role in the pathogenesis of breast cancer. We have performed coarse-grain ensemble molecular dynamics simulations to comprehensively sample the ErbB2 monomer-dimer association. Our results indicate a dynamic dimer state with a complex conformational landscape that is modulated with increasing lipid tail length. We resolve the native N-terminal "active" and C-terminal "inactive" conformations in all membrane compositions. However, the relative population of the N-terminal and C-terminal conformers is dependent on length of the saturated lipid tails. In short-tail membranes, additional non-specific dimers are observed which are reduced or absent in long-tailed bilayers. Our results indicate that the relative population as well as the structure of the dimer state is modulated by membrane composition. We have correlated these differences to local perturbations of the membrane around the receptor. Our work is an important step in characterizing ErbB dimers in healthy and diseased states and emphasize the importance of sampling lipid dynamics in understanding receptor association.

Resolving the conformational dynamics of ErbB growth factor receptor dimers

The combinatorial dimerization of the ErbB growth factor receptors (ErbB¹- ErbB⁴) are critical for their function. Here, we have characterized the conformational dynamics of ErbB transmembrane homo-dimers and hetero-dimers by using a coarse-grain simulation framework. All dimers, except ErbB^4-4 and ErbB^1-4, exhibit at least two conformations. The reported NMR structures correspond to one of these conformations, representing the N-terminal active state in ErbB^1-1 (RH2), ErbB^2-2 (RH1) and ErbB^4-4 (RH) homo-dimers and the LH dimer in ErbB^3-3 homo-dimer, validating the computational approach. Further, we predict a right-handed ErbB^3-3 dimer conformer that warrants experimental testing. The five hetero-dimers that have not yet been experimentally resolved display prominent right-handed dimers associating by the SmXXXSm motif. Our results provide insights into the constitutive signaling of ErbB⁴ after cleavage of the extracellular region. The presence of the inactive-like dimer conformers leading to symmetric kinase domains gives clues on the autoinhibition of the receptor dimers. The dimer states characterized here represent an important step towards understanding the combinatorial cross associations in the ErbB family.