In silico induced effect of N-ε-lysine acetylation on microtubule stability and subsequent interaction of microtubule-associated proteins

Plant systems have been considered valuable models for addressing fundamental questions of microtubule (MT) organization due to their considerable practical utility. Protein acetylation is a very common protein modification, and therate of acetylation can be modulated in cells in different biological states, and these changes can be detected at a molecular level. Here, we focused on K40, K112, and K394 residues as putative acetylation sites, which were shown to exist in both plants and mammals. Such residual effect of acetylation causes critical but unclear effect on MT stability. In turn, it was shown that acetylation indirectly affects the probability of interaction with different MAPs (Microtubule-associated proteins). In a multiscale study using an all-atom force field to reproduce several lattice-forming elements found on the surface the microtubule, we assembled a fragment of a plant microtubule composed of nine tubulins and used it as a model object along with the existing human complex. Triplets of tubulins assembled in a lattice cell were then simulated for both human and plant protein complexes, using a coarse-grained force field. We then analyzed the trajectories and identified some critical deformations of the MAP interaction surface. The initial coordinates were used to investigate the structural scenario in which autophagy-related protein 8 (ATG8) was able to interact with the MT fragment.

Aqueous Molecular Dynamics for Understanding Glycosaminoglycan Recognition by Proteins

Glycosaminoglycans (GAGs) are biopolymers that exist in most organisms. GAGs are known to bind to hundreds of proteins and partake in multiple biological processes such as growth, morphogenesis, inflammation, infection, and others. Their intrinsic structural heterogeneity and conformational variability introduce major challenges in experimental studies. On the other hand, recent advances in force field development and computational technology have yielded phenomenal opportunity to study thousands of GAG sequences simultaneously to understand recognition of target protein(s). Here, we describe experimental setup for conventional molecular dynamics simulations of GAGs to position an experimental biologist favorably in performance, analysis and interpretation of stability, specificity, and conformational properties of GAGs, while also elucidating their interactions with amino acid residues of a protein at an atomistic level in presence of water.

Native-like membrane models of E. coli polar lipid extract shed light on the importance of lipid composition complexity

BACKGROUND

Lipid-protein interactions stabilize protein oligomers, shape their structure, and modulate their function. Whereas in vitro experiments already account for the functional importance of lipids by using natural lipid extracts, in silico methods lack behind by embedding proteins in single component lipid bilayers. However, to accurately complement in vitro experiments with molecular details at very high spatio-temporal resolution, molecular dynamics simulations have to be performed in natural(-like) lipid environments.

RESULTS

To enable more accurate MD simulations, we have prepared four membrane models of E. coli polar lipid extract, a typical model organism, each at all-atom (CHARMM36) and coarse-grained (Martini3) representations. These models contain all main lipid headgroup types of the E. coli inner membrane, i.e., phosphatidylethanolamines, phosphatidylglycerols, and cardiolipins, symmetrically distributed between the membrane leaflets. The lipid tail (un)saturation and propanylation stereochemistry represent the bacterial lipid tail composition of E. coli grown at 37∘C until 3/4 of the log growth phase. The comparison of the Simple three lipid component models to the complex 14-lipid component model Avanti over a broad range of physiologically relevant temperatures revealed that the balance of lipid tail unsaturation and propanylation in different positions and inclusion of lipid tails of various length maintain realistic values for lipid mobility, membrane area compressibility, lipid ordering, lipid volume and area, and the bilayer thickness. The only Simple model that was able to satisfactory reproduce most of the structural properties of the complex Avanti model showed worse agreement of the activation energy of basal water permeation with the here performed measurements. The Martini3 models reflect extremely well both experimental and atomistic behavior of the E. coli polar lipid extract membranes. Aquaporin-1 embedded in our native(-like) membranes causes partial lipid ordering and membrane thinning in its vicinity. Moreover, aquaporin-1 attracts and temporarily binds negatively charged lipids, mainly cardiolipins, with a distinct cardiolipin binding site in the crevice at the contact site between two monomers, most probably stabilizing the tetrameric protein assembly.

CONCLUSIONS

The here prepared and validated membrane models of E. coli polar lipids extract revealed that lipid tail complexity, in terms of double bond and cyclopropane location and varying lipid tail length, is key to stabilize membrane properties over a broad temperature range. In addition, they build a solid basis for manifold future simulation studies on more realistic lipid membranes bridging the gap between simulations and experiments.