scholarly journals Inverse design of cholesterol attracting transmembrane helices reveals a paradoxical role of hydrophobic length

2021 ◽  
Author(s):  
Jeroen Methorst ◽  
Niek van Hilten ◽  
Herre Jelger Risselada

The occurrence of linear cholesterol-recognition motifs in alpha-helical transmembrane domains has long been debated. Here, we demonstrate the ability of a genetic algorithm guided by coarse-grained molecular dynamics simulations---a method coined evolutionary molecular dynamics (evo-MD)---to directly resolve the sequence which maximally attracts/sorts cholesterol within a single-pass alpha-helical transmembrane domain (TMDs). We illustrate that the evolutionary landscape of cholesterol attraction in membrane proteins is characterized by a sharp, well-defined global optimum. Surprisingly, this optimal solution features an unusual short hydrophobic block, consisting of typically only eight short chain hydrophobic amino acids, surrounded by three successive lysines. Owing to the membrane thickening effect of cholesterol, cholesterol-enriched ordered phases favor TMDs characterized by a long rather than a short hydrophobic length. However, this short hydrophobic pattern evidently offers a pronounced net advantage for the binding of free cholesterol in both coarse-grained and atomistic simulations. Attraction is mediated by the unique ability of cholesterol to snorkel within the hydrophobic core of the membrane and thereby shield deeply located lysines from the unfavorable hydrophobic surrounding. Since this mechanism of attraction is of a thermodynamic nature and is not based on molecular shape specificity, a large diversity of sub-optimal cholesterol attracting sequences can exist. The puzzling sequence variability of proposed linear cholesterol-recognition motifs is thus consistent with sub-optimal, unspecific binding of cholesterol. Importantly, since evo-MD uniquely enables the targeted design of recognition motifs for distinct fluid lipid membranes, we foresee wide applications for evo-MD in the biological and biomedical fields.

2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Md. Iqbal Mahmood ◽  
Hiroshi Noguchi ◽  
Kei-ichi Okazaki

Abstract F-Bin/Amphiphysin/Rvs (F-BAR) domain proteins play essential roles in biological processes that involve membrane remodelling, such as endocytosis and exocytosis. It has been shown that such proteins transform the lipid membrane into tubes. Notably, Pacsin1 from the Pacsin/Syndapin subfamily has the ability to transform the membrane into various morphologies: striated tubes, featureless wide and thin tubes, and pearling vesicles. The molecular mechanism of this interesting ability remains elusive. In this study, we performed all-atom (AA) and coarse-grained (CG) molecular dynamics simulations to investigate the curvature induction and sensing mechanisms of Pacsin1 on a membrane. From AA simulations, we show that Pacsin1 has internal structural flexibility. In CG simulations with parameters tuned from the AA simulations, spontaneous assembly of two Pacsin1 dimers through lateral interaction is observed. Based on the complex structure, we show that the regularly assembled Pacsin1 dimers bend a tensionless membrane. We also show that a single Pacsin1 dimer senses the membrane curvature, binding to a buckled membrane with a preferred curvature. These results provide molecular insights into polymorphic membrane remodelling.


Soft Matter ◽  
2018 ◽  
Vol 14 (24) ◽  
pp. 5019-5030 ◽  
Author(s):  
Eric J. Spangler ◽  
P. B. Sunil Kumar ◽  
Mohamed Laradji

The self-assembly of spherical nanoparticles, resulting from their adhesion on tensionless lipid membranes, is investigated through molecular dynamics simulations of a coarse-grained implicit-solvent model for self-assembled lipid membranes.


2020 ◽  
Author(s):  
Florencia Klein ◽  
Daniela Cáceres-Rojas ◽  
Monica Carrasco ◽  
Juan Carlos Tapia ◽  
Julio Caballero ◽  
...  

<p>Although molecular dynamics simulations allow for the study of interactions among virtually all biomolecular entities, metal ions still pose significant challenges to achieve an accurate structural and dynamical description of many biological assemblies. This is particularly the case for coarse-grained (CG) models. Although the reduced computational cost of CG methods often makes them the technique of choice for the study of large biomolecular systems, the parameterization of metal ions is still very crude or simply not available for the vast majority of CG- force fields. Here, we show that incorporating statistical data retrieved from the Protein Data Bank (PDB) to set specific Lennard-Jones interactions can produce structurally accurate CG molecular dynamics simulations. Using this simple approach, we provide a set of interaction parameters for Calcium, Magnesium, and Zinc ions, which cover more than 80% of the metal-bound structures reported on the PDB. Simulations performed using the SIRAH force field on several proteins and DNA systems show that using the present approach it is possible to obtain non-bonded interaction parameters that obviate the use of topological constraints. </p>


2021 ◽  
Vol 22 (14) ◽  
pp. 7375
Author(s):  
Julie Ledoux ◽  
Alain Trouvé ◽  
Luba Tchertanov

The kinase insert domain (KID) of RTK KIT is the key recruitment region for downstream signalling proteins. KID, studied by molecular dynamics simulations as a cleaved polypeptide and as a native domain fused to KIT, showed intrinsic disorder represented by a set of heterogeneous conformations. The accurate atomistic models showed that the helical fold of KID is mainly sequence dependent. However, the reduced fold of the native KID suggests that its folding is allosterically controlled by the kinase domain. The tertiary structure of KID represents a compact array of highly variable α- and 310-helices linked by flexible loops playing a principal role in the conformational diversity. The helically folded KID retains a collapsed globule-like shape due to non-covalent interactions associated in a ternary hydrophobic core. The free energy landscapes constructed from first principles—the size, the measure of the average distance between the conformations, the amount of helices and the solvent-accessible surface area—describe the KID disorder through a collection of minima (wells), providing a direct evaluation of conformational ensembles. We found that the cleaved KID simulated with restricted N- and C-ends better reproduces the native KID than the isolated polypeptide. We suggest that a cyclic, generic KID would be best suited for future studies of KID f post-transduction effects.


Sign in / Sign up

Export Citation Format

Share Document