Dynamic architecture of mycobacterial outer membranes revealed by all-atom simulations

  1. Turner P Brown
  2. Matthieu Chavent
  3. Wonpil Im  Is a corresponding author
  1. Department of Bioengineering, Lehigh University, United States
  2. Université de Toulouse, CNRS, Laboratoire de Microbiologie et de Génétique Moléculaires, CBI, France
  3. Department of Biological Sciences, Lehigh University, United States
6 figures and 2 additional files

Figures

Schematic of a mycobacterial outer membrane composed of mycolic acids (MAs), trehalose di- and mono-mycolate (TDM and TMM), phthiocerol dimycocerosate (PDIM), sulphoglycolipids (SGL), and pentacyl and diacyl trehalose (PAT and DAT).

This membrane is covalently attached to the peptidoglycan via the arabinogalactan.

Chemical structures of mycobacterial lipids and membranes snapshots.

(A) Chemical structure of α-mycolic acid (α-MA) with 78 carbons. (B–D). Initial three-dimensional (3D) conformations (left) and snapshots after 3 μs production of pure mycolic acid bilayers (right) with fully extended (MA_eU), semi-folded (MA_eZ), and fully folded (MA_W) conformations. Lipids in the upper leaflet are colored in light, medium, and dark blue, and lipids in the lower leaflet are colored in white. Oxygen atoms are shown as van der Waals spheres. (E–J) Chemical structures (left) and initial conformations (right) of the six outer leaflet lipids. (K) Snapshot after 3 μs production of the symmetric system containing outer leaflet lipids at 313 K (All_Lipids_313). (L) Snapshot after 3 μs production of the whole mycobacterial outer membrane with asymmetry of lipids at 313 K (Asym_313 system). Lipid colors in (K) and (L) match the colors of the initial 3D conformations.

Figure 3 with 1 supplement
Folding patterns and organization of mycomembrane inner leaflet symmetric bilayers.

(A, D, G) Overlaid structures of one mycolic acid (A MA_eU, B MA_eZ, and C MA_W) from every 10 frames of the final 500 ns production. Structures were aligned on 15 carbons at the top of the α-alkyl chain and 3 carbons at the top of the β-hydroxy chain: named C10 to C27 in the MA topologies. Oxygen atoms are displayed as red spheres. (B, E, H) Lipid dynamics at 313 K (solid lines) and 333 K (dashed lines). The three columns are terminal carbon density profiles, average lipid order parameters of carbons (Vermeer et al., 2007), and cyclopropane density profiles in 2 Å-wide slabs along the z-axis (i.e. the membrane normal) with the bilayer center at z=0. For carbon density profiles, lipids belonging to the lower leaflet are depicted in grey and those from the upper leaflet are depicted in blue. (C, F, I) Snapshots with varying initial compositions after 3 μs production. Lower leaflet lipids are shown in white, upper leaflet lipids are shown in light, medium, and dark blue. Oxygen atoms are displayed as red spheres. Terminal carbons and cyclopropane carbons are shown as spheres colored according to their leaflet.

Figure 3—figure supplement 1
Time series of inner leaflet symmetric bilayers with various starting conformations of α-mycolic acid.

(A) Membrane thickness. (B) Area per lipid.

Figure 4 with 1 supplement
Phase transition of pure fully-extended α-mycolic acid bilayers.

(A) Lateral diffusion coefficients from 313 K to 353 K. (B) Terminal carbon (C75) z-density profiles separated by leaflet along the temperature gradient. Upper leaflet lipids are shown with solid lines, and lower leaflet lipids shown with dashed lines. (C) Average order parameters of MA at varying temperatures. (D) Cyclopropane carbons (CC1, CC2) z-density profiles along the temperature gradient. (E) Snapshots along the temperature gradient, with lower leaflet lipids shown in white, upper leaflet lipids shown in colored lines, and oxygen atoms shown as spheres. Line and shaded area colors in (B, C, and D) correspond to those of upper leaflet lipids in E.

Figure 4—figure supplement 1
Time series of MA_eU bilayers at increasing temperature.

(A) Membrane thickness. (B) Area per lipid.

Figure 5 with 5 supplements
Lipid aggregation in mycomembrane outer leaflet symmetric bilayers.

(A, C) Snapshots of All_Lipids systems after 3 μs production at 313 K and 333 K, respectively (see Supplementary file 1 for system name, composition, and basic statistics). Lipid colors match the coloring from Figure 2. Water and ions are omitted for clarity. (B, D) Headgroup z position time series of PDIM (purple) and PAT (green) at 313 K and 333 K, respectively. Each line is a separate lipid molecule. (E) Lateral diffusion coefficients from All_Lipids system for each lipid type. Solid bars are from the 313 K system, and hatched bars are from the 333 K system. Error bars are the standard errors across three replicas. (F) Comparison of average deuterium order parameters of carbons along the z-axis for PDIM (light purple) and PAT (light green) at 313 K (solid lines) and 333 K (dashed lines). Shaded regions are the standard errors across three replicas. All analysis is averaged over the last 500 ns production.

Figure 5—figure supplement 1
Headgroup Z position time series for all lipid types and all outer leaflet symmetric systems at 313 K.
Figure 5—figure supplement 2
Headgroup z position time series for all lipid types and all outer leaflet symmetric systems at 333 K.
Figure 5—figure supplement 3
Outer leaflet lipid dynamics.

(A, B). Snapshots of All_Lipids system after 3 μs at 313 K and 333 K, respectively. Lipid headgroups are shown as van der Waals spheres, and the rest of each lipid is shown with the QuickSurf drawing method. (C, D) Headgroup z position time series of all lipid types at 313 K and 333 K, respectively. (E) Average deuterium order parameters of carbons along the z-axis for all lipid types. Solid lines are 313 K and dashed lines are 333 K. Scd was averaged over the final 500 ns production.

Figure 5—figure supplement 4
Scd vs carbon number for outer leaflet symmetric bilayers.

All results shown are for systems simulated at 313 K.

Figure 5—figure supplement 5
Time series of outer leaflet symmetric systems at 313 K and 333 K.

(A, C) Membrane thickness. (B, D) Area per lipid.

Figure 6 with 1 supplement
Vertical heterogeneity and lipid dynamics in asymmetric mycomembranes.

(A, C) Snapshots after 3 μs production at 313 K and 333 K, respectively. Lipid colors match the coloring from Figure 2, and lipids are shown with the QuickSurf drawing method. (B, D) Comparison of average deuterium order parameters of carbons along the z-axis for MA_eU (light blue), PDIM (light purple), and PAT (light green) at 313 K and 333 K, respectively. Solid lines are from the symmetric systems and dashed lines are from the asymmetric systems. (E, G) Full lipid z-density profiles for each lipid at 313 K and 333 K, respectively. Shaded regions are the standard errors across three replicas. (F, H) Headgroup z position time series of MA_eU (light blue), PDIM (light purple), and PAT (light green) at 313 K and 333 K, respectively. Each line is a separate lipid molecule. All analyses are averaged over the last 500 ns of production.

Figure 6—figure supplement 1
Scd and headgroup time series for asymmetric systems.

(A, B). Snapshots after 3 μs of production of asymmetric bilayers at 313 K and 333 K, respectively. (C, D) Comparison of the average deuterium order parameter for each lipid type along the z-axis for symmetric (solid lines) and asymmetric (dashed lines) bilayers at 313 K and 333 K, respectively. (E, F) Headgroup z position time series for each lipid type in the asymmetric systems at 31 3 K and 333 K, respectively.

Additional files

Supplementary file 1

Composition and structural properties of simulated membrane systems.

System names indicate lipid composition and simulation temperature (K). Lipid composition is given as the number of molecules per leaflet in the order: MA_eU, MA_sZ, MA_W, PDIM, TDM, TMM, DAT, PAT, and SGL. Systems labeled ‘both’ are symmetric bilayers; asymmetric systems (Asym_313 and Asym_333) are reported with inner and outer leaflets listed separately. System size corresponds to the lateral box dimension (Å), membrane thickness to the bilayer thickness (Å), and APL to the area per lipid (Ų). Reported structural properties for asymmetric systems correspond to the full bilayer.

https://cdn.elifesciences.org/articles/108644/elife-108644-supp1-v1.docx
MDAR checklist
https://cdn.elifesciences.org/articles/108644/elife-108644-mdarchecklist1-v1.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Turner P Brown
  2. Matthieu Chavent
  3. Wonpil Im
(2026)
Dynamic architecture of mycobacterial outer membranes revealed by all-atom simulations
eLife 14:RP108644.
https://doi.org/10.7554/eLife.108644.3