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.

Ca2+ ions facilitate the organization of the Annexin A2/S100A10 heterotetramer.

Annexin A2 (A2) is a member of the Annexin family, which contains Ca2+ -regulated phospholipid-binding proteins. Annexins associate with S100 proteins to form heterotetramers. The A2/S100A10 heterotetramer (A2t) is the most extensively studied of these heterotetramers. It induces membrane microdomain formation, causes membrane budding, and facilitates proliferation of some cancers. In this work, the first molecular dynamics (MD) study on the complete A2t of 868 amino acids was performed. MD trajectories of more than 600 ns each were generated for complete A2t complexes with and without Ca2+ ions. The outward extension of membrane-binding residues A2-K279 and A2-K281 was shown to be inhibited in the absence of Ca2+ as they were captured by Ca2+ -binding residue D322. F-actin binding residue A2-D339 was observed to occupy either an exposed or buried state in the absence of Ca2+ , while it only occupied the buried state in the presence of Ca2+ . The observed motions of the A2t subunits are highly organized with a strongly correlated central region which is negatively correlated with the periphery of the complex. The central region contains the S100A10 (p11) dimer, A2-N, and A2-I, while the periphery contains A2-II, A2-III, and A2-IV. Novel interactions between A2 and p11 were identified. A2 residues outside of A2-N (K80, R77, E82, and R145) had strong interactions with p11. Residue R145 of A2 may have a significant effect on the dynamics of the system, with its interaction resulting in asymmetric motions of A2. The presented results provide novel insights to inform future experimental studies.

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

Mapping the conformational pathway of biomolecules is a great challenge because of the large size and complexity of biomolecules. The nudged elastic band (NEB) method has been applied to study the reaction pathways for both small organic molecules and small peptides of a few amino acids. In this work, for the first time, the NEB method was employed to study the conformational pathways of Annexin A1, a membrane-binding protein of 334 amino acids. The N-terminal domain conformational change from the buried state within the core domain to the exposed state outside the core domain is a vital step for Annexin A1 to interact with membranes or target proteins. In this work, multiple molecular dynamics simulations using the NEB method were performed to simulate the N-terminal domain conformational pathway of Annexin A1. Our results suggested that the N-terminal domain of Annexin A1 is removed from the repeat III of the core domain in a sliding motion. The loop region of repeat III covering the N-terminal helix in the buried state does not lift up for the N-terminal to swing out of the pocket; instead, the N-terminal pulls out from the bottom of the core domain. The N-terminal domain linker region (S27-N42) flexibility is critical for the N-terminal domain conformational changes. Our results also suggested a two-step folding process for the helix D in repeat III, M247-V250 folds first followed by the folding of L251-E254. The results demonstrated that the NEB method could be an effective tool for theoretical studies on conformational pathways of biomolecules.