Structural dynamics of the mitochondrial ADP/ATP carrier revealed by molecular dynamics simulation studies

In Silico Analysis of the Structural Dynamics and Substrate Recognition Determinants of the Human Mitochondrial Carnitine/Acylcarnitine SLC25A20 Transporter

The Carnitine-Acylcarnitine Carrier is a member of the mitochondrial Solute Carrier Family 25 (SLC25), known as SLC25A20, involved in the electroneutral exchange of acylcarnitine and carnitine across the inner mitochondrial membrane. It acts as a master regulator of fatty acids β-oxidation and is known to be involved in neonatal pathologies and cancer. The transport mechanism, also known as "alternating access", involves a conformational transition in which the binding site is accessible from one side of the membrane or the other. In this study, through a combination of state-of-the-art modelling techniques, molecular dynamics, and molecular docking, the structural dynamics of SLC25A20 and the early substrates recognition step have been analyzed. The results obtained demonstrated a significant asymmetry in the conformational changes leading to the transition from the c- to the m-state, confirming previous observations on other homologous transporters. Moreover, analysis of the MD simulations' trajectories of the apo-protein in the two conformational states allowed for a better understanding of the role of SLC25A20 Asp231His and Ala281Val pathogenic mutations, which are at the basis of Carnitine-Acylcarnitine Translocase Deficiency. Finally, molecular docking coupled to molecular dynamics simulations lend support to the multi-step substrates recognition and translocation mechanism already hypothesized for the ADP/ATP carrier.

Mechanistic Insights into Multiple-step Transport of Mitochondrial ADP/ATP Carrier

The ADP/ATP carrier (AAC) is crucial for mitochondrial functions by importing ADP and exporting ATP across the inner mitochondrial membrane. However, the mechanism of highly specific ADP recognition and transport by AAC remains largely elusive. In this work, spontaneous ADP binding process to the ground c-state AAC was investigated through rigorous molecular dynamics simulations of over 31 microseconds in total. With improved simulation strategy, we have successfully identified a highly specific ADP binding site in the upper region of the cavity, and this site exhibits selectivity for ADP over ATP based on free-energy calculations. Sequence analyses on adenine nucleotide transporters also suggest that this subgroup uses the upper region of the cavity, rather than the previously proposed central binding site located at the bottom of the cavity to discriminate their substrates. Identification of the new site unveils the unusually high substrate specificity of AAC and explains the dependence of transport on the flexibility between anti and syn glycosidic conformers of ADP. Moreover, this new site together with the central site supports early biochemical findings. In light of these early findings, our simulations described a multi-step model in which the carrier uses different sites for substrate attraction, recognition and conformational transition. These results provide new insights into the transport mechanism of AAC and other adenine nucleotide transporters.

Investigating the Broad Matrix-Gate Network in the Mitochondrial ADP/ATP Carrier through Molecular Dynamics Simulations

The mitochondrial ADP/ATP carrier (AAC) exports ATP and imports ADP through alternating between cytosol-open (c-) and matrix-open (m-) states. The salt bridge networks near the matrix side (m-gate) and cytosol side (c-gate) are thought to be crucial for state transitions, yet our knowledge on these networks is still limited. In the current work, we focus on more conserved m-gate network in the c-state AAC. All-atom molecular dynamics (MD) simulations on a variety of mutants and the CATR-AAC complex have revealed that: (1) without involvement of other positive residues, the charged residues from the three Px[DE]xx[KR] motifs only are prone to form symmetrical inter-helical network; (2) R235 plays a determinant role for the asymmetry in m-gate network of AAC; (3) R235 significantly strengthens the interactions between H3 and H5; (4) R79 exhibits more significant impact on m-gate than R279; (5) CATR promotes symmetry in m-gate mainly through separating R234 from D231 and fixing R79; (6) vulnerability of the H2-H3 interface near matrix side could be functionally important. Our results provide new insights into the highly conserved yet variable m-gate network in the big mitochondrial carrier family.

Structural determinants of ligands recognition by the human mitochondrial basic amino acids transporter SLC25A29. Insights from molecular dynamics simulations of the c-state

In mitochondria, metabolic processes require the trafficking of solutes and organic molecules, such as amino acids. This task is accomplished by the Mitochondrial Carrier Family members (also known as SLC25), among which the SLC25A29 is responsible for the translocation of basic amino acids. In this regard, nitric oxide levels originated by the arginine mitochondrial catabolism have been shown to strongly affect cancer cells' metabolic status. Furthermore, the metabolic disease saccharopinuria has been linked to a mitochondrial dysregulation caused by a toxic intermediate of the lysine catabolism. In both cases, a reduction of the activity of SLC25A29 has been shown to ameliorate these pathological conditions. However, no detailed structural data are available on SLC25A29. In the present work, molecular modelling, docking and dynamics simulations have been employed to analyse the structural determinants of ligands recognition by SLC25A29 in the c-state. Results confirm and reinforce earlier predictions that Asn73, Arg160 and Glu161, and Arg257 represent the ligand contact points I, II, and III, respectively, and that Arg160, Trp204 and Arg257 form a stable interaction, likely critical for ligand binding and translocation. These results are discussed in view of the experimental data available for SLC25A29 and other homologous carriers of the same family.

Computational studies of the mitochondrial carrier family SLC25. Present status and future perspectives

The members of the mitochondrial carrier family, also known as solute carrier family 25 (SLC25), are transmembrane proteins involved in the translocation of a plethora of small molecules between the mitochondrial intermembrane space and the matrix. These transporters are characterized by three homologous domains structure and a transport mechanism that involves the transition between different conformations. Mutations in regions critical for these transporters’ function often cause several diseases, given the crucial role of these proteins in the mitochondrial homeostasis. Experimental studies can be problematic in the case of membrane proteins, in particular concerning the characterization of the structure–function relationships. For this reason, computational methods are often applied in order to develop new hypotheses or to support/explain experimental evidence. Here the computational analyses carried out on the SLC25 members are reviewed, describing the main techniques used and the outcome in terms of improved knowledge of the transport mechanism. Potential future applications on this protein family of more recent and advanced in silico methods are also suggested.

Sequence Features of Mitochondrial Transporter Protein Families

Mitochondrial carriers facilitate the transfer of small molecules across the inner mitochondrial membrane (IMM) to support mitochondrial function and core cellular processes. In addition to the classical SLC25 (solute carrier family 25) mitochondrial carriers, the past decade has led to the discovery of additional protein families with numerous members that exhibit IMM localization and transporter-like properties. These include mitochondrial pyruvate carriers, sideroflexins, and mitochondrial cation/H+ exchangers. These transport proteins were linked to vital physiological functions and disease. Their structures and transport mechanisms are, however, still largely unknown and understudied. Protein sequence analysis per se can often pinpoint hotspots that are of functional or structural importance. In this review, we summarize current knowledge about the sequence features of mitochondrial transporters with a special focus on the newly included SLC54, SLC55 and SLC56 families of the SLC solute carrier superfamily. Taking a step further, we combine sequence conservation analysis with transmembrane segment and secondary structure prediction methods to extract residue positions and sequence motifs that likely play a role in substrate binding, binding site gating or structural stability. We hope that our review will help guide future experimental efforts by the scientific community to unravel the transport mechanisms and structures of these novel mitochondrial carriers.

On the Detection and Functional Significance of the Protein-Protein Interactions of Mitochondrial Transport Proteins

Protein-protein assemblies are highly prevalent in all living cells. Considerable evidence has recently accumulated suggesting that particularly transient association/dissociation of proteins represent an important means of regulation of metabolism. This is true not only in the cytosol and organelle matrices, but also at membrane surfaces where, for example, receptor complexes, as well as those of key metabolic pathways, are common. Transporters also frequently come up in lists of interacting proteins, for example, binding proteins that catalyze the production of their substrates or that act as relays within signal transduction cascades. In this review, we provide an update of technologies that are used in the study of such interactions with mitochondrial transport proteins, highlighting the difficulties that arise in their use for membrane proteins and discussing our current understanding of the biological function of such interactions.

Molecular dynamics simulations on apo ADP/ATP carrier shed new lights on the featured motif of the mitochondrial carriers

The ADP/ATP carrier (AAC) is a transporter responsible for the equal molar exchange of cytosolic ADP and ATP synthesized within mitochondrial matrix across the mitochondrial membrane. Its primary structure consists of three homologous repeats, and each repeat contains a conserved motif that is shared by all members of the mitochondrial carrier family (MCF). Although these MCF motif residues cluster together in the crystal structure of AAC, detailed analyses on the interactions among the motif residues are still limited. In the present study, all-atom molecular dynamics (MD) simulations of up to 10 μs have been carried out on AAC, and interactions and structural dynamics of the MCF motif residues have been specifically investigated. Our simulations have revealed: i) a very asymmetrical electrostatic network at the bottom of the pocket of apo AAC, ii) the asymmetrical interactions between the Pro kink region and the [YWF][KR] G motif in three repeats, iii) the role of the conserved Arg residues in stabilizing the C-ends of the odd-numbered helices, iv) the structural change of the [YWF][KR] G motif and its potential involvement in substrate translocation process. Our results highlight the asymmetry of the MCF residues in the three repeats, which might contribute to the ability of the carriers to transport the asymmetrical substrates. Our observations provide microscopic basis for further research on the translocation mechanism of mitochondrial carriers.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Mitochondrial ADP/ATP Carrier in Dodecylphosphocholine Binds Cardiolipins with Non-native Affinity

Biophysical investigation of membrane proteins generally requires their extraction from native sources using detergents, a step that can lead, possibly irreversibly, to protein denaturation. The propensity of dodecylphosphocholine (DPC), a detergent widely utilized in NMR studies of membrane proteins, to distort their structure has been the subject of much controversy. It has been recently proposed that the binding specificity of the yeast mitochondrial ADP/ATP carrier (yAAC3) toward cardiolipins is preserved in DPC, thereby suggesting that DPC is a suitable environment in which to study membrane proteins. In this communication, we used all-atom molecular dynamics simulations to investigate the specific binding of cardiolipins to yAAC3. Our data demonstrate that the interaction interface observed in a native-like environment differs markedly from that inferred from an NMR investigation in DPC, implying that in this detergent, the protein structure is distorted. We further investigated yAAC3 solubilized in DPC and in the milder dodecylmaltoside with thermal-shift assays. The loss of thermal transition observed in DPC confirms that the protein is no longer properly folded in this environment.

The switching mechanism of the mitochondrial ADP/ATP carrier explored by free-energy landscapes

The ADP/ATP carrier (AAC) of mitochondria has been an early example for elucidating the transport mechanism alternating between the external (c-) and internal (m-) states (M. Klingenberg, Biochim. Biophys. Acta 1778 (2008) 1978-2021). An atomic resolution crystal structure of AAC is available only for the c-state featuring a three repeat transmembrane domain structure. Modeling of transport mechanism remained hypothetical for want of an atomic structure of the m-state. Previous molecular dynamics studies simulated the binding of ADP or ATP to the AAC remaining in the c-state. Here, a full description of the AAC switching from the c- to the m-state is reported using well-tempered metadynamics simulations. Free-energy landscapes of the entire translocation from the c- to the m-state, based on the gyration radii of the c- and m-gates and of the center of mass, were generated. The simulations revealed three free-energy basins attributed to the c-, intermediate- and m-states separated by activation barriers. These simulations were performed with the empty and with the ADP- and ATP-loaded AAC as well as with the poorly transported AMP and guanine nucleotides, showing in the free energy landscapes that ADP and ATP lowered the activation free-energy barriers more than the other substrates. Upon binding AMP and guanine nucleotides a deeper free-energy level stabilized the intermediate-state of the AAC2 hampering the transition to the m-state. The structures of the substrate binding sites in the different states are described producing a full picture of the translocation events in the AAC.

Multi-scaled normal mode analysis method for dynamics simulation of protein-membrane complexes: A case study of potassium channel gating motion correlations

Membrane proteins play critically important roles in many cellular activities such as ions and small molecule transportation, signal recognition, and transduction. In order to fulfill their functions, these proteins must be placed in different membrane environments and a variety of protein-lipid interactions may affect the behavior of these proteins. One of the key effects of protein-lipid interactions is their ability to change the dynamics status of membrane proteins, thus adjusting their functions. Here, we present a multi-scaled normal mode analysis (mNMA) method to study the dynamics perturbation to the membrane proteins imposed by lipid bi-layer membrane fluctuations. In mNMA, channel proteins are simulated at all-atom level while the membrane is described with a coarse-grained model. mNMA calculations clearly show that channel gating motion can tightly couple with a variety of membrane deformations, including bending and twisting. We then examined bi-channel systems where two channels were separated with different distances. From mNMA calculations, we observed both positive and negative gating correlations between two neighboring channels, and the correlation has a maximum as the channel center-to-center distance is close to 2.5 times of their diameter. This distance is larger than recently found maximum attraction distance between two proteins embedded in membrane which is 1.5 times of the protein size, indicating that membrane fluctuation might impose collective motions among proteins within a larger area. The hybrid resolution feature in mNMA provides atomic dynamics information for key components in the system without costing much computer resource. We expect it to be a conventional simulation tool for ordinary laboratories to study the dynamics of very complicated biological assemblies. The source code is available upon request to the authors.

Antiporters of the mitochondrial carrier family

The eukaryotic transport protein family SLC25 consists of mitochondrial carriers (MCs) that are recognized on the sequence level by a threefold repeated and conserved signature motif. The majority of MCs characterized so far catalyzes strict exchanges of substrates across the mitochondrial inner membrane. The substrates are nucleotides, metabolic intermediates, and cofactors that are required in cytoplasmic and matrix metabolism. This review summarizes and discusses the current knowledge of the antiport mechanism(s) of MCs that has been deduced from determining transport characteristics and by analyzing structural, sequence, and mutagenesis data. The mode of transport varies among different MCs with respect to how the substrate translocation depends on the electrical and pH gradients across the mitochondrial inner membrane, for example, the ADP/ATP carrier is electrogenic (electrophoretic), the GTP/GDP carrier is dependent on the pH gradient, the aspartate/glutamate carrier is dependent on both, and the oxoglutarate/malate carrier is independent of them. The structure of the bovine ADP/ATP carrier consists of a six-transmembrane α-helix bundle with a pseudo-threefold symmetry and a closed matrix gate. By using this structure as a template in homology modeling, residues engaged in substrate binding and the formation of a cytoplasmic gate in MCs have been proposed. The functional importance of the residues of the binding site, the matrix, and the cytoplasmic gates is supported by transport activities of different MCs with single point mutations. Cumulative evidence has been used to postulate a general transport mechanism for MCs.

Single-nucleotide evolution quantifies the importance of each site along the structure of mitochondrial carriers

Mitochondrial carriers are membrane-embedded proteins consisting of a tripartite structure, a three-fold pseudo-symmetry, related sequences, and similar folding whose main function is to catalyze the transport of various metabolites, nucleotides, and coenzymes across the inner mitochondrial membrane. In this study, the evolutionary rate in vertebrates was screened at each of the approximately 50,000 nucleotides corresponding to the amino acids of the 53 human mitochondrial carriers. Using this information as a starting point, a scoring system was developed to quantify the evolutionary pressure acting on each site of the common mitochondrial carrier structure and estimate its functional or structural relevance. The degree of evolutionary selection varied greatly among all sites, but it was highly similar among the three symmetric positions in the tripartite structure, known as symmetry-related sites or triplets, suggesting that each triplet constitutes an evolutionary unit. Based on evolutionary selection, 111 structural sites (37 triplets) were found to be important. These sites play a key role in structure/function of mitochondrial carriers and are involved in either conformational changes (sites of the gates, proline-glycine levels, and aromatic belts) or in binding and specificity of the transported substrates (sites of the substrate-binding area in between the two gates). Furthermore, the evolutionary pressure analysis revealed that the matrix short helix sites underwent different degrees of selection with high inter-paralog variability. Evidence is presented that these sites form a new sequence motif in a subset of mitochondrial carriers, including the ADP/ATP translocator, and play a regulatory function by interacting with ligands and/or proteins of the mitochondrial matrix.

Mutations in the Mitochondrial Citrate Carrier SLC25A1 are Associated with Impaired Neuromuscular Transmission

Background and Objective

Congenital myasthenic syndromes are rare inherited disorders characterized by fatigable weakness caused by malfunction of the neuromuscular junction. We performed whole exome sequencing to unravel the genetic aetiology in an English sib pair with clinical features suggestive of congenital myasthenia.

Methods

We used homozygosity mapping and whole exome sequencing to identify the candidate gene variants. Mutant protein expression and function were assessed in vitro and a knockdown zebrafish model was generated to assess neuromuscular junction development.

Results

We identified a novel homozygous missense mutation in the SLC25A1 gene, encoding the mitochondrial citrate carrier. Mutant SLC25A1 showed abnormal carrier function. SLC25A1 has recently been linked to a severe, often lethal clinical phenotype. Our patients had a milder phenotype presenting primarily as a neuromuscular (NMJ) junction defect. Of note, a previously reported patient with different compound heterozygous missense mutations of SLC25A1 has since been shown to suffer from a neuromuscular transmission defect. Using knockdown of SLC25A1 expression in zebrafish, we were able to mirror the human disease in terms of variable brain, eye and cardiac involvement. Importantly, we show clear abnormalities in the neuromuscular junction, regardless of the severity of the phenotype.

Conclusions

Based on the axonal outgrowth defects seen in SLC25A1 knockdown zebrafish, we hypothesize that the neuromuscular junction impairment may be related to pre-synaptic nerve terminal abnormalities. Our findings highlight the complex machinery required to ensure efficient neuromuscular function, beyond the proteomes exclusive to the neuromuscular synapse.

The transmembrane prolines of the mitochondrial ADP/ATP carrier are involved in nucleotide binding and transport and its biogenesis

The mitochondrial ADP/ATP carrier (Ancp) is a paradigm of the mitochondrial carrier family, which allows cross-talk between mitochondria, where cell energy is mainly produced, and cytosol, where cell energy is mainly consumed. The members of this family share numerous structural and functional characteristics. Resolution of the atomic structure of the bovine Ancp, in a complex with one of its specific inhibitors, revealed interesting features and suggested the involvement of some particular residues in the movements of the protein to perform translocation of nucleotides from one side of the membrane to the other. They correspond to three prolines located in the odd-numbered transmembrane helices (TMH), Pro-27, Pro-132, and Pro-229. The corresponding residues of the yeast Ancp (Pro-43, Ser-147, and Pro-247) were mutated into alanine or leucine, one at a time and analysis of the various mutants evidenced a crucial role of Pro-43 and Pro-247 during nucleotide transport. Beside, replacement of Ser-147 with proline does not inactivate Ancp and this is discussed in view of the conservation of the three prolines at equivalent positions in the Ancp sequences. These prolines belong to the signature sequences of the mitochondrial carriers and we propose they play a dual role in the mitochondrial ADP/ATP carrier function and biogenesis. Unexpectedly their mutations cause more general effects on mitochondrial biogenesis and morphology, as evidenced by measurements of respiratory rates, cytochrome contents, and also clearly highlighted by fluorescence microscopy.

Mapping multiple potential ATP binding sites on the matrix side of the bovine ADP/ATP carrier by the combined use of MD simulation and docking

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier-AAC-was crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). The protein consists of a six-transmembrane helix bundle that defines the nucleotide translocation pathway, which is closed towards the matrix side due to sharp kinks in the odd-numbered helices. In this paper, we describe the interaction between the matrix side of the AAC transporter and the ATP(4-) molecule using carrier structures obtained through classical molecular dynamics simulation (MD) and a protein-ligand docking procedure. Fifteen structures were extracted from a previously published MD trajectory through clustering analysis, and 50 docking runs were carried out for each carrier conformation, for a total of 750 runs ("MD docking"). The results were compared to those from 750 docking runs performed on the X-ray structure ("X docking"). The docking procedure indicated the presence of a single interaction site in the X-ray structure that was conserved in the structures extracted from the MD trajectory. MD docking showed the presence of a second binding site that was not found in the X docking. The interaction strategy between the AAC transporter and the ATP(4-) molecule was analyzed by investigating the composition and 3D arrangement of the interaction pockets, together with the orientations of the substrate inside them. A relationship between sequence repeats and the ATP(4-) binding sites in the AAC carrier structure is proposed.

Tetrahydrofolate recognition by the mitochondrial folate transporter

A mitochondrial carrier family (MCF) of transport proteins facilitates the transfer of charged small molecules across the inner mitochondrial membrane. The human genome has ∼50 genes corresponding to members of this family. All MCF proteins contain three repeats of a characteristic and conserved PX(D/E)XX(K/R) motif thought to be central to the mechanism of these transporters. The mammalian mitochondrial folate transporter (MFT) is one of a few MCF members, known as the P(I/L)W subfamily, that have evolved a tryptophan residue in place of the (D/E) in the second conserved motif; the function of this substitution (Trp-142) is unclear. Molecular dynamics simulations of the MFT in its explicit membrane environment identified this tryptophan, as well as several other residues lining the transport cavity, to be involved in a series of sequential interactions that coordinated the movement of the tetrahydrofolate substrate within the transport cavity. We probed the function of these residues by mutagenesis. The mutation of every residue identified by molecular dynamics to interact with tetrahydrofolate during simulated transit into the aqueous channel severely impaired folate transport. Mutation of the subfamily-defining tryptophan residue in the MFT to match the MCF consensus at this position (W142D) was incompatible with cell survival. These studies indicate that MFT Trp-142, as well as other residues lining the transporter interior, coordinate tetrahydrofolate descent and positioning of the substrate in the transporter basin. Overall, we identified residues in the walls and at the base of the transport cavity that are involved in substrate recognition by the MFT.

Conformational dynamics of the bovine mitochondrial ADP/ATP carrier isoform 1 revealed by hydrogen/deuterium exchange coupled to mass spectrometry

The mitochondrial adenine nucleotide carrier (Ancp) catalyzes the transport of ADP and ATP across the mitochondrial inner membrane, thus playing an essential role in cellular energy metabolism. During the transport mechanism the carrier switches between two different conformations that can be blocked by two toxins: carboxyatractyloside (CATR) and bongkrekic acid. Therefore, our understanding of the nucleotide transport mechanism can be improved by analyzing structural differences of the individual inhibited states. We have solved the three-dimensional structure of bovine carrier isoform 1 (bAnc1p) in a complex with CATR, but the structure of the carrier-bongkrekic acid complex, and thus, the detailed mechanism of transport remains unknown. Improvements in sample processing in the hydrogen/deuterium exchange technique coupled to mass spectrometry (HDX-MS) have allowed us to gain novel insights into the conformational changes undergone by bAnc1p. This paper describes the first study of bAnc1p using HDX-MS. Results obtained with the CATR-bAnc1p complex were fully in agreement with published results, thus, validating our approach. On the other hand, the HDX kinetics of the two complexes displays marked differences. The bongkrekic acid-bAnc1p complex exhibits greater accessibility to the solvent on the matrix side, whereas the CATR-bAnc1p complex is more accessible on the intermembrane side. These results are discussed with respect to the structural and biochemical data available on Ancp.

Structural approaches of the mitochondrial carrier family

The transport of solutes across the inner mitochondrial membrane is highly selective and necessitates membrane proteins mainly from the mitochondrial carrier family (MCF). These carriers are required for the transport of a variety of metabolites implicated in all the important processes occurring within the mitochondrial matrix. Due to its high abundance, the ADP/ATP carrier (AAC) is the member of the family that was studied most. It is the first mitochondrial carrier for which a high-resolution X-ray structure is known. The carrier was crystallized in the presence of a strong inhibitor, the carboxyatractyloside (CATR). The structure gives an insight not only into the overall fold of mitochondrial carriers in general but also into atomic details of the AAC in a conformation that is open toward the intermembrane space (IMS). Molecular dynamics simulations indicate the first events occurring to the carrier after the binding of ADP. A careful analysis of the primary sequences of all the carriers in light with the structure highlights properties of the protein that are related to the substrate.

Computational studies of H5N1 influenza virus resistance to oseltamivir

Influenza A (H5N1) virus is one of the world's greatest pandemic threats. Neuraminidase (NA) inhibitors, oseltamivir and zanamivir, prevent the spread of influenza, but drug-resistant viruses have reduced their effectiveness. Resistance depends on the binding properties of NA-drug complexes. Key residue mutations within the active site of NA glycoproteins diminish binding, thereby resulting in drug resistance. We performed molecular simulations and calculations to characterize the mechanisms of H5N1 influenza virus resistance to oseltamivir and predict potential drug-resistant mutations. We examined two resistant NA mutations, H274Y and N294S, and one non-drug-resistant mutation, E119G. Six-nanosecond unrestrained molecular dynamic simulations with explicit solvent were performed using NA-oseltamivir complexes containing either NA wild-type H5N1 virus or a variant. MM_PBSA techniques were then used to rank the binding free energies of these complexes. Detailed analyses indicated that conformational change of E276 in the Pocket 1 region of NA is a key source of drug resistance in the H274Y mutant but not in the N294S mutant.

Cardiolipin and mitochondrial carriers

Members of the mitochondrial carrier family interact with cardiolipin (CL) as evident from a variety of functional and structural effects. CL stabilises carrier proteins on isolation with detergents, with the P(i) carrier as the prime example. CL is required for transport in reconstituted vesicles, prime examples are the P(i)- and ADP/ATP carrier (AAC). CL binds to the AAC in a graded manner; 6 CL/AAC dimer bind tightly as measured on the (31)P NMR time scale. 2 additional CL/dimer bind reversibly and a fast exchanging envelope of phospholipids includes CL as measured on the ESR time scale. In the crystal structure of the CAT-AAC complex 3 CL bind to the periphery of the AAC in a three-fold pseudo-symmetry. The binding of CL is implicated to contribute lowering the high transition energy barriers in the AAC. Para-functions of the AAC, as in the mitochondrial pore transition (MPT) and in cell death are linked to the CL binding of the AAC. Ca(++) or oxidants can sequester or destroy AAC bound CL, rendering AAC labile, allowing pore formation and degradation. Thus AAC, by being vital for energy transfer, constitutes an Achilles heel in the eukaryotic cell. AAC together with CL is also engaged in respiratory supercomplexes. Different from AAC the similarly structured uncoupling protein (UCP1) has no tightly bound CL, but CL addition lowers affinity of the inhibitory nucleotide binding that may contribute to the physiological regulation of the uncoupling activity by ATP.

Binding of ADP in the mitochondrial ADP/ATP carrier is driven by an electrostatic funnel

The ADP/ATP carrier (AAC) is a membrane protein of paramount importance for the energy-fueling function of the mitochondria, transporting ADP from the intermembrane space to the matrix and ATP in the opposite direction. On the basis of the high-resolution, 2.2-A structure of the bovine carrier, a total of 0.53 micros of classical molecular dynamics simulations were conducted in a realistic membrane environment to decipher the early events of ADP (3-) translocation across the inner membrane of the mitochondria. Examination of apo-AAC underscores the impermeable nature of the carrier, impeding passive transport of permeants toward the matrix. The electrostatic funnel illuminated from three-dimensional mapping of the electrostatic potential forms a privileged passageway anticipated to drive the diphosphate nucleotide rapidly toward the bottom of the internal cavity. This conjecture is verified in the light of repeated, independent numerical experiments, whereby the permeant is dropped near the mouth of the mitochondrial carrier. Systematic association of ADP (3-) to the crevice of the AAC, an early event in its transport across the inner membrane, is accompanied by the formation of an intricate network of noncovalent bonds. Simulations relying on the use of an adaptive biasing force reveal for the first time that the proposed binding site corresponds to a minimum of the free energy landscape delineating the translocation of ADP (3-) in the carrier. The present work paves the way to the design of novel nucleotides and new experiments aimed at unveiling key structural features in the chronology of ADP/ATP transport across the mitochondrial membrane.

Electrostatic funneling of substrate in mitochondrial inner membrane carriers

Exchange of ATP and ADP across mitochondrial membrane replenishes the cytoplasm with newly synthesized ATP and provides the mitochondria with the substrate ADP for oxidative phosphorylation. The sole means of this exchange is the mitochondrial ADP/ATP carrier (AAC), a membrane protein that is suggested to cycle between two conformationally distinct states, cytosolic-open (c-state) and matrix-open (m-state), thereby shuttling nucleotides across the inner mitochondrial membrane. However, the c-state is the only structurally resolved state, and the binding site of ADP remains elusive. Here, we present approximately 0.3 mus of all-atom MD simulations of the c-state revealing rapid, spontaneous binding of ADP to deeply positioned binding sites within the AAC lumen. To our knowledge, a complete ligand-binding event has heretofore not been described in full atomic detail in unbiased simulations. The identified ADP-bound state and additional simulations shed light on key structural elements and the initial steps involved in conversion to the m-state. Electrostatic analysis of trajectories reveals the presence of an unusually strong positive electrostatic potential in the lumen of AAC that appears to be the main driving force for the observed spontaneous binding of ADP. We provide evidence that the positive electrostatic potential is likely a common attribute among the entire family of mitochondrial carriers. In addition to playing a key role in substrate recruitment and translocation, the electropositivity of mitochondrial carriers might also be critical for their binding to the negatively charged environment of the inner mitochondrial membrane.