If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, ChinaInstitute of Genetics, and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, ChinaInstitute of Genetics, and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, ChinaInstitute of Genetics, and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, ChinaInstitute of Genetics, and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
Cardiolipin (CL) has been shown to play a crucial role in regulating the function of proteins in the inner mitochondrial membrane. As the most abundant protein of the inner mitochondrial membrane, the ADP/ATP carrier (AAC) has long been the model of choice to study CL-protein interactions, and specifically bound CLs have been identified in a variety of crystal structures of AAC. However, how CL binding affects the structural dynamics of AAC in atomic detail remains largely elusive. Here we compared all-atom molecular dynamics simulations on bovine AAC1 in lipid bilayers with and without CLs. Our results show that on the current microsecond simulation time scale: 1) CL binding does not significantly affect overall stability of the carrier or structural symmetry at the matrix-gate level; 2) pocket volumes of the carrier and interactions involved in the matrix-gate network become more heterogeneous in parallel simulations with membranes containing CLs; 3) CL binding consistently strengthens backbone hydrogen bonds within helix H2 near the matrix side; and 4) CLs play a consistent stabilizing role on the domain 1-2 interface through binding with the R30:R71:R151 stacking structure and fixing the M2 loop in a defined conformation. CL is necessary for the formation of this stacking structure, and this structure in turn forms a very stable CL binding site. Such a delicate equilibrium suggests the strictly conserved R30:R71:R151stacking structure of AACs could function as a switch under regulation of CLs. Taken together, these results shed new light on the CL-mediated modulation of AAC function.
The inner mitochondrial membrane (IMM) exhibits a variety of unique properties: forming cristae leads to very large surface area; proteins are densely packed within it; and the membrane contains featured phospholipids, cardiolipins (CLs). IMM is crowded with membrane proteins involved in energy production and substrate transportation, including oxidative phosphorylation (OXPHOS) complexes and mitochondrial carriers. Some studies estimate the protein/lipid ratio reaches to as high as 3:1 (
). Because IMM is densely packed with OXPHOS complexes, there is little potential to improve OXPHOS capacity through increasing density of OXPHOS proteins. One important strategy is to increase the surface area of IMM through forming highly folded cristae structure. It is reported that CLs might play an important role on cristae biogenesis of mitochondria (
). CLs are unique lipids preserved during evolution from aerobic bacteria to mitochondria, and in mitochondria, CLs are mainly distributed within the inner leaflet of IMM where they are generated (
). The structure of CL is featured with a small hydroxyl head group and four hydrophobic acyl tails (Fig. 1A), and it usually carries two negative net charges in physiological conditions. Its special conical shape makes CLs have potential to relax curvature strain in cristae. Moreover, CLs have been reported to stabilize and regulate the function of a bunch of mitochondrial proteins such as cytochrome c, ADP/ATP carrier (AAC), and ATPase (
Fig. 1Structures of cardiolipin and the ADP/ATP carrier. A: Chemical and tertiary structure of CL. B: The tripartite secondary structure of AAC. Domains are shown in distinct colors with domain 1 in yellow, domain 2 in orange, and domain 3 in red. C: Side view (top) and matrix view (bottom) of the tertiary structure of AAC with three specifically bound CLs. CLs are highlighted in blue stick-ball model. AAC, ADP/ATP carrier; CL, cardiolipin.
AAC is the most abundant protein in IMM and plays a central role in OXPHOS. It is responsible for importing ADP into the mitochondrial matrix and exporting ATP out of the mitochondria for living cells in a 1:1 exchange. The exchange of ADP and ATP across the IMM is completed through AAC alternating from the cytosol-open (c-) or ADP-waiting conformation to the matrix-open (m-) or ATP-waiting conformation via the transition state (
). Benefitting from the high expression levels and the presence of state-specific inhibitors, AAC has long been paradigm of the big mitochondrial carrier family (MCF) (
), and this family functions to mediate a variety of metabolites across IMM. Currently 53 MCF members have been identified in human genome, and they form the biggest solute carrier (SLC) subfamily which is also referred to as SLC25. Like other members of MCF, AAC is featured with three homologous repeats (Fig. 1B), each containing about 100 amino acids with highly conserved MCF motif Px[DE]xx[KR]xRxQ-(matrix loop and matrix helix)-[DE]Gxxxx[YWF][KR]G near the matrix side (
), with the cysteine-containing matrix short helices hosting the sequence motif -50-QYKGxxDCxRK-60 (sequence numbering referred to the bovine AAC1 sequence) (
Specific cardiolipin binding interferes with labeling of sulfhydryl residues in the adenosine diphosphate/adenosine triphosphate carrier protein from beef heart mitochondria.
ADP/ATP carrier protein from beef heart mitochondria has high amounts of tightly bound cardiolipin, as revealed by phosphorus-31 nuclear magnetic resonance.
). In these structures, six transmembrane (TM) helices form a bundle surrounding a cone-shaped cavity opening toward the intermembrane space in c-state and toward the matrix in m-state. Each odd-numbered helix is connected to the following even-numbered helix at the matrix side through a flexible loop, a short amphipathic helix and a linker helix, which forms the basic structural scaffold for one repeat domain of AAC (Fig. 1B, C). In c-state structure, three repeated homologous domains of AAC show three-fold pseudosymmetry, and at each domain–domain interface, there is one specifically bound CL (Fig. 1C). Based on the static crystal structures, it was proposed that CLs bind with the exposed positively charged amide groups and helix dipoles at the N-terminal ends of even-numbered helices and short matrix helices (
). However, molecular dynamics (MD) simulations have shown that positively charged and polar residues also play an important role in binding with CLs (
Structural determinants of ligands recognition by the human mitochondrial basic amino acids transporter SLC25A29. Insights from molecular dynamics simulations of the c-state.
). The difference from crystal structures and MD simulations highlights the importance of introducing structural dynamics in exploring lipid–protein interactions, which is also supported by studies from other groups (
Lipid-loving ANTs: molecular simulations of cardiolipin interactions and the organization of the adenine nucleotide translocase in model mitochondrial membranes.
). Due to varied distributions of the positively charged and polar residues in the three specific CL binding sites, our previous work reported that the three sites are highly asymmetric, and the asymmetric CL binding behavior at the three sites is quite consistent with the asymmetry observed in the matrix side of the m-state structure (
). However, how CL binding affects the structural dynamics of AAC in atomic detail remains largely elusive.
In the current work, three 3-μs all-atom MD simulations on c-state apo AAC with CLs (CL-1, 2, 3) are compared to the three 3-μs simulations without CLs (PC-1, 2, 3). We analyzed the impact of CLs on both macroscopic parameters such as pocket volume of AAC and microscopic parameters such as interactions among MCF motif residues. Significantly, we have revealed a strictly conserved stacking structure R30:R71:R151 at the domain 1-2 interface, and we speculate that this stacking structure functions as a structural switch under regulation of CLs.
Materials and methods
Initial preparation of the simulation systems
The starting protein coordinates were obtained from the CATR-inhibited crystal structure of AAC (PDB entry: 1okc) (
). The termini residues of AAC were manipulated into a charged state by protonating the N terminus and deprotonating the C terminus. A pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane and a mixed membrane composed of POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamine, and CL were built with the Membrane Builder modules of VMD (
), respectively. To build the AAC-POPC system, a hole was made within the POPC lipid bilayer, and the carrier alone was embedded. To build the AAC-CL system, the carrier bound with three CLs was inserted manually into the mixed bilayer with the same method. Then, the system was solvated with TIP3P (
) water molecules, and 100 mmol/L Na+ and Cl− were added to the system in addition to the neutralized chloride ions. The AAC-POPC system contains 70,769 atoms, including one AAC molecule, 219 POPC lipids, 12,226 water molecules, 23 Na+, and 42 Cl-. The AAC-CL system contains 82,971 molecules, including one AAC molecule, 110 POPCs, 88 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamines, 22 CLs, 15,740 water molecules, 53 Na+, and 28 Cl-. The detailed steps for systems setup were reported in our previous papers (
Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations.
). The four-round energy minimizations were performed to remove unfavorable contacts and prevent the AAC-lipid system from collapsing at the subsequent MD simulations. Positional restraints were sequentially applied on all heavy atoms, main-chain atoms, and Cα atoms of AAC using a force constant of 1000 kcal/mol/Å2, and the last round was minimization without any positional restraint. Following energy minimizations, a series of equilibration steps were carried out for each system. In the first step, the system temperature was heated from 50 K to 310 K with positional restraints applied on heavy atoms of the protein, and the time step was set to 1 fs. In the subsequent three steps, the lipid bilayer and the solvent were fully equilibrated with positional restraints on all heavy atoms, main-chain atoms, and Cα atoms of AAC, respectively, and the time step was increased to 2 fs. After equilibration, three independent 3-μs production simulations were run in parallel for each system. The coordinates of each system were saved every 10 ps. The system temperature and pressure were constantly maintained (310 K, 1 bar) using the v-rescale method (
). The pdb files for the initial protein structure and structures after equilibrations are provided as a supplementary compressed file (pdb_structures.zip).
Data analyses and visualization
The programs in GROMACS 4.5.5 were utilized for most trajectory analyses. Specifically, root mean square deviation, root mean square fluctuations (RMSF), and spatial density function were calculated with the g_rms, g_rmsf, and g_spatial programs, respectively. Principal component analysis was calculated using g_covar and g_anaeig programs in Gromacs (
). Occupancies of salt bridges and H-bonds (hydrogen bonds) were all calculated based on the last 2 μs of each trajectory, with a distance cut-off of 0.33 nm. Time evolutions of the pocket volumes of AAC along the simulation trajectory were calculated with Epock (
Effect of CL on structural dynamics of c-state AAC in general
Structures of apo c-state AAC were well maintained in all six simulations with or without CLs in the lipid bilayer. Root mean square deviations of all Cα atoms of AAC during the simulations do not show obvious differences in the two conditions (Fig. 2A), indicating that the overall stability of AAC is not significantly affected by CL binding. Then we did RMSF calculations to further analyze the effect of CL binding on the structural flexibility of different regions in AAC. Of special interest, the M2 loop exhibits extreme stability in all three trajectories with CLs, while in all the other three trajectories without CLs the M2 loop becomes highly dynamic (Figs. 2B, D and S1). This suggests that CL dramatically stabilizes the M2 loop, but for the rest regions of the carrier, CL binding does not show obvious impact. In both conditions, the cytoplasmic sides of TM helices show much higher RMSF values than the matrix side, suggesting the cytoplasmic sides are more flexible (Fig. 2B). This is consistent with the fact that c-state structure closes at the matrix side and opens wide toward the cytoplasmic side. Moreover, flexibility at the cytoplasmic side is not symmetrically distributed among different domains. More specifically, the C2 loop and the cytoplasmic ends of H1 and H6 exhibit the highest RMSF values, while the values of the C1 loop is much lower than that of the C2 loop (Fig. 2B) which is consistent with previous simulations on shorter timescale (
). Principal component analysis also shows that major movements of AAC are observed in the region around the C2 loop and the cytoplasmic sides of H1 and H6, while the region around the C1 loop including cytoplasmic ends of H2 and H3 are relatively stable in both conditions with and without CLs (Fig. 2C, D). Our recent work had identified that the second basic patch (K91K95R187), N115 and the tyrosine ladder (Y186Y190Y194) including F191 form the specific ADP binding site (
). Thermodynamic stability of this region below the C1 loop could be crucial for efficient substrate recognition.
Fig. 2Effect of CL on the structural dynamics of AAC in general. A: Time evolutions of the root mean square deviation (RMSD) values of AAC Cα atoms during 3- s simulations in absence (red lines) and in presence (blue lines) of CLs. B: Time evolutions of root mean square fluctuations (RMSF) values of AAC backbone atoms in presence (red lines) and in absence of CLs (black lines). C–D: Principal component analysis on AAC trajectories. The motion of AAC C atoms is projected along the first eigenvector, with amplitude of the motion following the color scale from red to blue in simulations with CLs and from yellow to blue in simulations without CLs. AAC, ADP/ATP carrier; CL, cardiolipin.
). CATR is much bigger than the natural substrate ADP, and therefore it was proposed that the CATR-bound AAC could adopt an extreme configuration and open wider than the apo c-state which is the start of the translocation cycle (
). Our calculations show that in PC-1, PC-2, and CL-1, the pocket volume decreases first and then remains relatively stable till the end of the simulations (Fig. 3A, B and D). The significant decrease of the pocket volume at the beginning hundreds of nanoseconds indicates that the crystal conformation of AAC with CATR removed is actually unstable in the lipid bilayer. Similar results were also obtained in previous MD simulations on much shorter time scale (
). In PC-3, the pocket volume of AAC shows a decreasing trend throughout the simulation (Fig. 3C), and we reported similar results with an alternative method in our previous work (
). In CL-2 and CL-3, the pocket volumes fluctuate heavily during the 3-μs simulations. This result indicates that the cytoplasmic side of c-state AAC could be more flexible in presence of CL. However, we can not exclude the possibility that MD simulations with CL may need much longer time to get more converged results.
Fig. 3Time evolutions of pocket volume of c-state apo AAC in simulations without (A–C) and with (D–F) CLs. The smoothed average is shown in black lines, and the blue dashed lines highlight the volume calculated for the crystal structure of AAC. AAC, ADP/ATP carrier; CL, cardiolipin.
CL stabilizes the M2 loop through the conserved R30: R79: R151 stacking structure
Although the three parallel trajectories in presence of CL show heterogeneity based on the above pocket volume analysis (Fig. 3B), they consistently show the stabilized M2 loop (Fig. 2A). The stabilizing effect of CL on the M2 loop is also illustrated through superposed snapshots taken from the trajectories with and without CLs (Fig. 4A, B). The M2 loop connects the matrix end of H3 with the N terminus of the matrix helix h34 within domain 2. In the crystal structure of bovine AAC1, guanidinium group of R151 in the M2 loop stacks with guanidinium group of R71 of [YWF][KR]G motif in domain 1 and therefore attaches the M2 loop to domain 1. We observe that R71 also stacks with R30, the capping arginine of Px[DE]xx[KR]xR motif in domain 1. Different from the capping arginine residues R139 and R236 in the other two domains, R30 moves four residues ahead from its canonical position in sequence and therefore moves about one helical turn upward in structure (
). This movement makes the guanidinium group of R30 much closer to that of R71, which facilitates their stacking interaction. In this way, R151, R71, and R30 form a nice stacking structure (Fig. 4C), and similar stacking structure is commonly found in high-quality structures in PDB database (
Fig. 4The featured stacking structure on the domain 1-2 interface. A–B: Superposed conformations of 50 snapshots of domain 2 of AAC with and without CLs are extracted from MD simulation per 60 ns. The M2 loop is highlighted by the circle (red dashed line). C: The featured R30:R71:R151 stacking structure. R30, R71 and R151 are shown in light blue spheres, and the bound CL is shown in yellow sticks. The M2 loop is highlighted in salmon. D: Electrostatic interactions that stabilize the R30:R71:R151 stacking structure in the simulations with CL. The interactions are shown in black dashed lines. E: Dissociation of the stacking structure in the simulations without CL. The electrostatic interactions with the dissociated R30, R70, and R151 residues are shown in black dashed lines. AAC, ADP/ATP carrier; CL, cardiolipin; MD, molecular dynamics.
The electrostatic repulsion between these negatively charged guanidinium groups is balanced by interactions of these stacking residues with environment. In the crystal structure, R151 is also stabilized by side chain of D143 and backbone carbonyl of F153. R71 binds with backbone carbonyl oxygen atoms of D143 and G145. R30 forms hydrogen bonds with side chain of S68 and Q64 and backbone carbonyl oxygen atoms of A141 and D143. However, these interactions are not strong enough, and the R30:R71:R151 stacking structure is not maintained in all the simulations without CLs (Fig. 4E), which leads to thermodynamic flexibility of the M2 loop. In contrast, in presence of CLs, the stacking structure is well maintained throughout the simulations, and the M2 loop is extremely stable (Fig. 4A, C and D).
In our previous work, the stacked R71 and R151 bind to CL with high occupancies, and they are critical residues of the specific CL binding site I at the domain 1-2 interface (
). Therefore, these simulation results suggest that CL is necessary for the formation of the R30:R71:R151 stacking structure, and this stacking structure in turn forms a very stable CL binding site. Such delicate equilibrium prompts us to propose that the strictly conserved R30:R71:R151 stacking motif could serve as a structural switch under regulation of CLs. The simulation results also demonstrate that through stabilizing the R30:R71:R151 stacking structure, CLs play a crucial role in fixing the M2 loop in a pretty defined conformation.
CL binding stabilizes the domain 1-2 interface
Specific CL binding sites are located at domain-domain interfaces, involving positively charged residues, N termini of even-numbered and short matrix helices (
). To investigate the effect of CL binding on three domain-domain interfaces, we calculated distances between N terminus of even number helices and N terminus of matrix helices in the next domain with and without CLs. More specifically, minimum distances between C atoms of residues in triplet 75 (A75, V178, S275) and C atoms of corresponding residues in triplet 53 (L156, T253, I53) were calculated respectively for each simulation trajectory (Fig. 5A–C). Here, the triplet refers to three symmetric amino acids in the homologous tripartite structure, numbered according to the residue number in the first repeat of bovine AAC1 (
). Frequency distributions of these distances suggest that CL binding consistently show a stabilization effect on the domain 1-2 interface (Fig. 5A). For the other two interfaces, big variations among different trajectories either with or without CLs indicate that CL binding does not have significant effect (Fig. 5B, C). The stabilization effect of CL binding on the domain 1-2 interface is also consistent with the spatial distribution function results in which higher density of CLs is observed at this interface compared to those at the other two interfaces (Fig. 5D–F). These results imply that CL is of special significance in stabilizing the domain 1-2 interface.
Fig. 5Impact of CL binding on domain-domain interfaces. A–C: Distribution of minimum distances between N terminus of even-numbered helices and N terminus of matrix helices in MD simulations with and without CLs. D–F: Isosurface of CLs density in the three CL binding sites of AAC. AAC, ADP/ATP carrier; CL, cardiolipin; MD, molecular dynamics.
In the crystal structures of CATR-bound AAC, charged residues from the Px[DE]xx[KR] motifs form a symmetrical and cyclic interhelical salt bridge network on the matrix side (m-gate) (
). Based on simulations without CLs, we reported that after removing CATR, the charged residues from three Px[DE]xx[KR] motifs mainly form intrahelical salt bridges, and together with polar residues and positive residues K79, R235, and R279, form a broad asymmetric m-gate network (
). Here, we aim to evaluate how presence of CL may affect the salt bridges and H-bonds involved in the broad m-gate network.
Calculations show that in contrast to the simulations without CL, three parallel simulations with CLs exhibit tremendous heterogeneity for a majority of interactions involved in the m-gate network (Figs. 6, S2 and S3, supplemental Table S1). For example, occupancies of the interhelical salt bridge D134:R234 in CL-1, CL-2, and CL-3 are 96%, 58%, and 6% respectively; for the intrahelical D231:R234, the occupancies are 10%, 42%, and 94% in CL-1, CL-2, and CL-3 respectively; and for D134:R79, the occupancies are 4%, 91%, and 73% in these three trajectories. The tremendous heterogeneity might reflect the increased flexibility of AAC in presence of CL. However, we can not exclude the possibility that the simulations with CLs may need longer time to get fully relaxed.
Fig. 6Effect of CL binding on the m-gate network. Scatter plots for the occupancies of the salt bridges and H-bonds involved in the broad m-gate network. Each data point of group CL (n=3) or group PC (n=3) is signified as a red or blue solid circle. Average occupancies are represented as longer horizontal lines, and the standard deviations are shown as shorter lines. CL, cardiolipin.
Worthy of mentioning, despite the heterogeneity for many m-gate interactions in simulations with CLs, the E29:R279 salt bridge is extremely stable with a high occupancy of almost 100% in each trajectory either with or without CL. This result is consistent with our previous finding that R279 has little impact on other interactions of the m-gate network (
Structural determinants of ligands recognition by the human mitochondrial basic amino acids transporter SLC25A29. Insights from molecular dynamics simulations of the c-state.
Effect of CL binding on other MCF structural elements
In addition to the m-gate, we also evaluate the effect of CL binding on interactions among other MCF structural elements including proline kink, capping arginine, [DE]G, and [YWF][KR]G motifs (supplemental Table S2 and supplemental Fig. S4). Results show that the electrostatic interactions between the [YWF][KR]G motif and Pro kink region are very strong and not affected by CL binding. The interactions between capping arginines and termini of even-numbered helices are different among different domains. R30 binds carbonyl oxygen of A141 with high occupancy in each trajectory, either with or without CLs. R139 binds carbonyl oxygen of M238 with lower occupancy in simulations without CLs than in simulations with CLs. One possible reason is that without CLs, the M2 loop is very flexible and therefore E152 may dynamically bind with R139 and deteriorate the R139:M238(O) interaction. Heterogeneity of the R236:Q36(O) interactions in parallel simulations either with or without CLs may just result from thermodynamic fluctuations of the carrier. The occupancies of the Q64:R30 interactions are quite different in parallel simulations with CLs, and based on our observations, it is mainly caused by different conformations of S68 which may bind with either G65 backbone, R30 or Q64. The M3 loop is intervened between D167 and R139, and based on our simulation trajectories, the interaction between R139 and the M3 loop is generally stronger than that between D167 and the M3 loop, and CL binding does not directly affect these interactions. The E264:R236 salt bridge is extremely stable and not affected by CL binding, and this salt bridge was proposed to be important in AAC transport (
). In the previous work, we have reported that the [YWF][KR]G motif lost its featured -turn structure in one of the three trajectories without presence of CLs (
). While this featured structure was perfectly maintained in all the three simulations with CLs, suggesting CLs help stabilize the -turn structure of the [YWF][KR]G motif (
). In a previous study, the metric ψ was introduced to evaluate the deviation degree of backbone structural symmetry of triplet from the ideal C3 symmetry (
). Specifically, ψ was calculated through averaging the deviations of the three internal angles of the triangle formed by three Cα atoms of a triplet from the internal angle of an equilateral triangle (60°). Higher ψ values suggest a higher deviation from C3 symmetry for the triplet residues. In that work, the authors reported that presence of CLs slightly improves the structural symmetry of the triplets in H2, H4, and H6 of m-state AAC around m-gate level (
). Following the same method, our recent work on wild-type and mutant c-state AAC showed that the R79A mutation drastically deteriorates structural symmetry of the triples in H1, H3, and H5 near the m-gate level (
). Here, we use this method to investigate effect of CL binding on the structural symmetry of c-state AAC around m-gate level. The heat maps of time evolutions of ψ values show no significant changes during each simulation, and no big differences were observed among different trajectories (Fig. 7A). In addition to the terminal triplet 37 (V37, A142, M239), the triplet 33 (L33, T138, R235) consistently show higher ψ values compared to other triplet planes in both conditions, highlighting the crucial role of the strictly conserved R235 in introducing structural asymmetry to c-state AAC near the m-gate level. To compare the differences in the two conditions in a quantitative way, we calculated the averaged ψ values over all three 3-μs trajectories with or without CLs. The results show that all triplets show no significant differences (P > 0.05) (Fig. 7B). These results suggest that presence of CL does not significantly affect structural symmetry of c-state AAC around the m-gate level on the current simulation time scale.
Fig. 7Structural symmetry of c-state AAC around the m-gate level in simulations with and without CLs. A: Time evolutions of ψ values in c-state AAC with (left panel) and without (right panel) CLs. Gradient colors from blue to red are used to indicate the deviation from three-fold symmetry of each triplet. B: The average ψ values calculated over trajectories with (orange) and without (black) CLs. The ψ values are presented as mean ± SEM. To evaluate the differences of ψ values in the two groups (with CLs and without CLs), a two independent-samples t test was performed, with a two-tailed P-value less than 0.05 considered to be significant. The P-values of all the triplets are bigger than 0.05 and therefore are not labeled in the figure. AAC, ADP/ATP carrier; CL, cardiolipin.
Effect of CL binding on backbone H-bonds within TM helices
Flexibility within the TM helices is of special significance for the normal function of membrane proteins. To evaluate the effect of CL binding on the flexibility of each TM helix, we calculated the average distance of each backbone H-bond (Ni-Oi-4) within the TM helices over the simulations with and without CLs. The increased Ni-Oi-4 distance is usually caused by structural kinks or increased flexibility within the TM helix. Results show that CL binding does not affect the average distances of backbone H-bonds within the odd-numbered helices H1, H3, and H5 (Fig. 8A, C and E), and within each of these helices, the conserved kink proline (P27, P132, P229) leads to a sharp peak for the Ni-Oi-4 distances in the region nearby (from triplet 26 to triplet 29). Compared to H1 and H5, H3 is different in that a secondary peak consistently appears around F129 in all the six trajectories (Fig. 8C), and this could be caused by the bent H2 that connects tightly with the H3 through the short C1 loop.
Fig. 8Effect of CL binding on backbone H-bonds within TM helices. (A–F) The average distance between the backbone N atom of residue i and the backbone O atom of neighboring fourth residue (i-4) within TM helices of AAC in simulations with (red lines) and without (black lines) CLs. AAC, ADP/ATP carrier; CL, cardiolipin; TM, transmembrane. The arrow in figure B highlights the only proline residue P82 within the transmembrane region besides the three conserved kink prolines.
The three even-numbered helices H2, H4, and H6 are quite different from each other in their backbone H-bonds and effect upon CL binding. H2 is special in its hydrophilic property and bent backbone (
). The bent backbone of H2 is demonstrated by two sharp peaks in the backbone Ni-Oi-4 distances in the region from R79 to Q84 and in the region from K91 to Y94 (Fig. 8B). Parallel simulations show consistent effect of CL binding on the backbone Ni-Oi-4 distances of H2 (Fig. 8B). CL binding consistently decreases the Ni-Oi-4 distances in the region from R79 to P82 (Fig. 8B), and the decreased amplitude is maintained relative stable throughout the simulations (supplemental Fig. S5). Although the average distances of F88(O)-D92(N) and A89(O)-93(N) are also decreased upon CL binding (Fig. 8B), time evolutions show that the decrease is discrete (supplemental Fig. S5), indicating that the decrease is caused by transient conformational changes within H2 backbone. Moreover, CL binding increases the average distances of K91(O)-D95(N) and K93(O)-97(N) (Fig. 8B), and the time evolutions also show that the increase is discrete (supplemental Fig. S5). These results indicate that CL binding stabilizes the matrix side of H2 but makes the cytosolic side of H2 more dynamic.
For H4, the low average Ni-Oi-4 distances throughout the TM helix indicates that H4 conserves its perfect α-helical structure during the simulations, and this structure is not affected by CL binding (Fig. 8D). For H6, whether in presence of CL or not, increased backbone Ni-Oi-4 distances are observed in many regions of H6, and parallel trajectories show tremendous heterogeneity (Fig. 8F). This indicates the dynamic property of H6 in both conditions. In fact, instability of H6 and its conformational changes during the simulations without CLs have been described in detail in our previous work (
In the current work, we have investigated the effect of CL binding on the structure and dynamics of c-state AAC. Because of the dynamic property of AAC in apo state, parallel trajectories show considerable heterogeneity either with or without CLs, which makes it hard to get simple conclusions on the effect of CLs. However, we do observe that CLs play a consistent stabilizing role on the interface between domain 1 and domain 2. Our newly published results on wild-type and mutant AAC indicate the vulnerability of domain 1-2 interface near the matrix side without CLs (
). The current work shows that CLs could stabilize this interface through binding with the R30:R71:R151 stacking structure and therefore fix the M2 loop in a pretty defined conformation. Meanwhile, CLs stabilize the matrix part of the dynamic H2 and introduce more heterogeneity among parallel simulations for the m-gate network.
Previous studies have suggested that the M2 loop undergoes dramatic conformational changes during state transitions (
Importance of loops of mitochondrial ADP/ATP carrier for its transport activity deduced from reactivities of its cysteine residues with the sulfhydryl reagent eosin-5-maleimide.
). Crystal structures of AAC in the two states also show significant difference in the conformations of the M2 loop. Our MD simulations on c-state apo bovine AAC1 in the current work have shown that CL plays a critical role in stabilizing the crystal conformation of the M2 loop through binding with the R30:R71:R151 stacking structure. Without CLs, the R30:R71:R151 stacking structure is disrupted, and the M2 loop loses its crystal conformation and becomes very flexible in the simulations. All the three Arg residues within the stacking structure are strictly conserved throughout different species of AAC (Fig. 9A), and they form a highly asymmetric structural element at the matrix side of the carrier. Among the capping arginine residues from the three homologous domains, R30 is the only one that moves one helical turn upward in structure (
). It forms an extremely stable interdomain interaction with A141 but unstable intradomain interaction with Q64. Prior studies have noted the importance of these interactions. These interactions were reported to participate in protecting the dynamic interface between domains 1 and 2 (
). R71 is the only strictly conserved middle residue of the three [YWF][KR]G motifs (Fig. 9A). The side chain of R71 is exposed to the lipid bilayer, while its backbone oxygen is oriented toward the interior of the AAC pocket. In all crystal structures of AAC, R71 binds with the guanidine group of the R137 through its backbone oxygen and is probably the only surface residue that has a chance to involve in the m-gate electrostatic network. In the three matrix loops, R151 is the only strictly conserved basic residue located just before the 50-QYKGxxDCxRK-60 motif (
) (Fig. 9A). At the matrix side, the three domains close tightly in the c-state AAC, while domain 1 separates from the other two domains in the m-state. This implies that the R30:R71:R151 stacking structure at domain 1-2 interface needs to dissociate and reform repetitively during state transitions.
Fig. 9The solvent exposed β-turn structure of the [YWF][KR]G motif in domain 1 is protected by the R30:R71:R151 stacking structure. A: Sequence logo presentations on the neighboring residues of A26, R30, R71, and R151 and their counterpart residues in homologous domains. Forty-three sequences of various AAC subtypes from totally 24 species were used for alignment. B: Structural illustration for the R30:R71:R151 stacking structure (magenta spheres) and β-turn structure of the [YWF][KR]G motif (shown in red sticks) that is exposed to the cavity solvent （highlighted with red arrow） due to the small residue A26 (green spheres). Solvent is shown in blue surface, and the stably bound CL is shown in yellow sticks. AAC, ADP/ATP carrier; CL, cardiolipin.
Previous analysis based on 70 thousand crystal structures suggests that stacking arginines are potentially controlled or switched by charged membrane lipids and may play an important role in protein function (
). CLs are featured negatively charged lipids within IMM, and in our previous work on c-state AAC, we highlighted the dynamic property of CLs and their potential involvement in adapting to the conformational changes during state transitions (
). Our results also show that in the presence of CLs, the R30:R71:R151 stacking structure is stably maintained throughout our simulations (Fig. 3D), and meanwhile, this stabilized stacking structure forms a good CLs binding site, and both R71 and R151 exhibit very high CLs binding occupancies in previous atomic and coarse grain MD simulations (
Lipid-loving ANTs: molecular simulations of cardiolipin interactions and the organization of the adenine nucleotide translocase in model mitochondrial membranes.
). Without presence of CLs, the stacking dissociates in all three trajectories (Fig. 3E). Therefore, based on the dynamic property of CLs binding and sensitivity of the R30:R71:R151 stacking structure to CLs, we infer that CLs play an important role in the dynamic dissociation and reformation of the stacking structure.
As described in previous work, the small alanine (A26) instead of a tyrosine appears before P27, and this makes the β-turn structure of the [YWF][KR]G motif in domain 1 faces directly to the cavity solvent (Fig. 9B) and therefore contributes to the vulnerability of the domain 1-2 interface. This conserved vulnerable interface is protected by the conserved R30:R71:R151 stacking structure that is under regulation of CLs (Fig. 9B), and this leads us to infer that CLs may play a regulation role on the function of the carrier through this stacking structure.
Here we propose a model to describe the relationship between CLs, R30:R71:R151 stacking structure and state transition (Fig. 10). In the ground c-state, the R30:R71:R151 stacking structure at the vulnerable domain 1-2 interface is stabilized by CLs binding. One typical binding mode is that one phosphate group of CL binds simultaneously with R71 and R151, with the other phosphate group binding preferentially with N terminus of H2, and the other typical binding mode is that one phosphate group binds with R71, and the other phosphate group binds with both R151 and N terminus of H2 (
). Upon ADP binding, the original CL binding mode with the R30:R71:R151 stacking structure can easily be disturbed. This will lead to the dissociation of the stacking structure, and further facilitates opening of the vulnerable domain 1-2 interface. Therefore, this model explains, at least partially, the activation role of CLs on the transport of AAC (
Fig. 10A model for how CLs regulate the c- to m- state transition of AAC through the featured arginine stacking structure. The linked orange circles represent the two phosphate groups of CLs. Residues of the R30:R71:R151 stacking structure are shown as spheres. Domain 1 and 2 are colored in light blue and pink, respectively. The N atoms in the side chain of stacking arginines are colored in dark blue. AAC, ADP/ATP carrier; CL, cardiolipin.
). However, it was reported that the stabilization effect from the inhibitor CATR is much more significant, and little is known about the stabilization effect of CL on apo AAC (
). Well-tempered metadynamics simulations starting from c-state AAC embedded in the pure POPC lipid bilayer also showed that c-state is energetically more favorable than m-state (
ADP/ATP mitochondrial carrier MD simulations to shed light on the structural-dynamical events that, after an additional mutation, restore the function in a pathological single mutant.
Adenine nucleotide translocase is acetylated in vivo in human muscle: Modeling predicts a decreased ADP affinity and altered control of oxidative phosphorylation.
), and the c-state structure was well maintained in all these simulations. Nevertheless, a recent work demonstrated that in absence of CL, yeast AAC2 is mostly in the m-state (
), implying that only in presence of CLs c-state is more stable than m-state. In a previous model, the excess positive charges near the matrix side of AAC facilitate the opening of the m-gate when the bound CLs are sequestered by Ca2+ (
), which also implies that c-state becomes less stable without presence of CLs. Therefore, relationship between CL and AAC conformations seems rather complicated and merits further investigations. The current work is limited in that we only evaluate the effect of CL binding on structural dynamics of c-state AAC, and we believe that further work on m-state and state transition pathways combined with free energy calculations may help to better understand the regulation mechanism of CL on AAC conformations and functions.
Data availability
Data are available upon reasonable request to the corresponding author.
The authors declare that they have no conflicts of interest with the contents of this article.
Acknowledgments
MD simulations were carried out at National Supercomputer Center in Tianjin, and the calculations were performed on TianHe-1(A).
Author contribution
Q. Y. and S. Y. investigation; B. M., S. Y., and Q. Y. data curation; Q. Y. and S. Y. writing-original draft; X. C. conceptualization; X. C. supervision; X. C. funding acquisition; X. C. writing-review & editing.
Funding and additional information
The research was supported by the National Natural Science Foundation of China (Grant No. 32171241 to X. C.) and the Natural Science Foundation of Zhejiang Province (Grant No. LY18C050002 to X. C.).
Specific cardiolipin binding interferes with labeling of sulfhydryl residues in the adenosine diphosphate/adenosine triphosphate carrier protein from beef heart mitochondria.
ADP/ATP carrier protein from beef heart mitochondria has high amounts of tightly bound cardiolipin, as revealed by phosphorus-31 nuclear magnetic resonance.
Structural determinants of ligands recognition by the human mitochondrial basic amino acids transporter SLC25A29. Insights from molecular dynamics simulations of the c-state.
Lipid-loving ANTs: molecular simulations of cardiolipin interactions and the organization of the adenine nucleotide translocase in model mitochondrial membranes.
Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations.
Importance of loops of mitochondrial ADP/ATP carrier for its transport activity deduced from reactivities of its cysteine residues with the sulfhydryl reagent eosin-5-maleimide.
ADP/ATP mitochondrial carrier MD simulations to shed light on the structural-dynamical events that, after an additional mutation, restore the function in a pathological single mutant.
Adenine nucleotide translocase is acetylated in vivo in human muscle: Modeling predicts a decreased ADP affinity and altered control of oxidative phosphorylation.
00060-8&id=gr1.jpg){kind=link}
00060-8&id=gr2.jpg){kind=link}
00060-8&id=gr3.jpg){kind=link}
00060-8&id=gr4.jpg){kind=link}
00060-8&id=gr5.jpg){kind=link}
00060-8&id=gr6.jpg){kind=link}
00060-8&id=gr7.jpg){kind=link}
00060-8&id=gr8.jpg){kind=link}
00060-8&id=gr9.jpg){kind=link}
00060-8&id=gr10.jpg){kind=link}