Nudged Elastic Band Study on the N-Terminal Domain Conformational Pathways of Annexin A1 from a Buried State to an Exposed State

Coarse-grained modeling of annexin A2-induced microdomain formation on a vesicle

Annexin A2 (A2)-induced microdomain formation is a key step in biological processes such as Ca2+-mediated exocytosis in neuroendocrine cells. In this work, a total of 15 coarse-grained molecular dynamics simulations were performed on vesicle models having a diameter of approximately 250 Å for 15 μs each using the Martini2 force field. Five simulations were performed in the presence of 10 A2, 5 in the presence of A2 but absence of PIP2, and 5 simulations in the absence of A2 but presence of PIP2. Consistent results were generated among the simulations. A2-induced PIP2 microdomain formation was observed and shown to occur in three phases: A2-vesicle association, localized A2-induced PIP2 clustering, and A2 aggregation driving PIP2 microdomain formation. The relationship between A2 aggregation and PIP2 microdomain formation was quantitatively described using a novel method which calculated the variance among protein and lipid positions via the Fréchet mean. A large reduction in PIP2 variance was observed in the presence of A2 but not in its absence. This reduction in PIP2 variance was proportional to the reduction observed in A2 variance and demonstrates that the observed PIP2 microdomain formation is dependent upon A2 aggregation. The three-phase model of A2-induced microdomain formation generated in this work will serve as a valuable guide for further experimental studies and the development of novel A2 inhibitors. No microdomain formation was observed in the absence of A2 and minimal A2-membrane interaction was observed in the absence of PIP2.

Modeling the Annexin A1-S100A11 heterotetramer: a molecular dynamics investigation of structure and correlated motion

Annexin A1 (A1) has been shown to form a tetrameric complex (A1t) with S100A11 which is implicated in calcium homeostasis and EGFR pathways. In this work, a full-length model of the A1t was generated for the first time. Multiple molecular dynamics simulations were performed on the complete A1t model for several hundred nanoseconds each to assess the structure and dynamics of A1t. These simulations yielded three structures for the A1 N-terminus (ND) which were identified via principal component analysis. The orientations and interactions of the first 11 A1-ND residues for all three structures were conserved, and their binding modes were strikingly similar to those of the Annexin A2 N-terminus in the Annexin A2-p11 tetramer. In this study, we provided detailed atomistic information for the A1t. Strong interactions were identified within the A1t between the A1-ND and both S100A11 monomers. Residues M3, V4, S5, E6, L8, K9, W12, E15, and E18 of A1 were the strongest interactions between A1 and the S100A11 dimer. The different conformations of the A1t were attributed to the interaction between W12 of the A1-ND with M63 of S100A11 which caused a kink in the A1-ND. Cross-correlation analysis revealed strong correlated motion throughout the A1t. Strong positive correlation was observed between the ND and S100A11 in all simulations regardless of conformation. This work suggests that the stable binding of the first 11 residues of A1-ND to S100A11 is potentially a theme for Annexin-S100 complexes and that the flexibility of the A1-ND allows for multiple conformations of the A1t.Communicated by Ramaswamy H. Sarma.

Interactions and conformational changes of annexin A2/p11 heterotetramer models on a membrane: a molecular dynamics study

Ca2+-dependent membrane-binding by the Annexin A2/p11 heterotetramer (A2t) plays an important role in various biological processes including fibrinogen activation and exocytosis in neuroendocrine cells. Two models where A2t associates with a single membrane surface were generated and used to perform molecular dynamics simulations. The first model mimics initial A2t-membrane binding through both Annexin A2 (A2) subunits of A2t (TS model) while the second model mimics A2t-binding through a single A2 subunit (OS model). Conformational changes were summarized using principal component analysis (PCA), simulation snapshots, and distance plots from the simulations. The full TS model, including the p11 dimer, fully associates with the membrane adopting a stable structure with little conformational variation as evidence by PCA. The unassociated subunits of the OS model moved toward the membrane. The molecular mechanics/Generalized-Born surface area (MMGBSA) method was applied to investigate the energetics of the models. The MMGBSA results demonstrated that R63 of p11 was the primary contributor to the p11-membrane interaction. The TS model results were both consistent with those found in the literature and provide novel insights about the specific residues driving the A2t-membrane interaction. Additionally, it represents the most complete model of A2t on the membrane surface available.Communicated by Ramaswamy H. Sarma.