Bacteriocin AS-48 binding to model membranes and pore formation as revealed by coarse-grained simulations

Bacteriocin AS-48 is a membrane-interacting peptide that acts as a broad-spectrum antimicrobial against Gram-positive and Gram-negative bacteria. Prior Nuclear Magnetic Resonance experiments and the high resolution crystal structure of AS-48 have suggested a mechanism for the molecular activity of AS-48 whereby the peptide undergoes transition from a water-soluble to a membrane-bound state upon membrane binding. To help interpret experimental results, we here simulate the molecular dynamics of this binding mechanism at the coarse-grained level. By simulating the self-assembly of the peptide, we predict induction by the bacteriocin of different pore types consistent with a “leaky slit” model.     

The intrinsically disordered N-terminus of the voltage-dependent anion channel

The voltage-dependent anion channel (VDAC) is a critical β-barrel membrane protein of the mitochondrial outer membrane, which regulates the transport of ions and ATP between mitochondria and the cytoplasm. In addition, VDAC plays a central role in the control of apoptosis and is therefore of great interest in both cancer and neurodegenerative diseases. Although not fully understood, it is presumed that the gating mechanism of VDAC is governed by its N-terminal region which, in the open state of the channel, exhibits an α-helical structure positioned midway inside the pore and strongly interacting with the β-barrel wall. In the present work, we performed molecular simulations with a recently developed force field for disordered systems to shed new light on known experimental results, showing that the N-terminus of VDAC is an intrinsically disordered region (IDR). First, simulation of the N-terminal segment as a free peptide highlighted its disordered nature and the importance of using an IDR-specific force field to properly sample its conformational landscape. Secondly, accelerated dynamics simulation of a double cysteine VDAC mutant under applied voltage revealed metastable low conducting states of the channel representative of closed states observed experimentally. Related structures were characterized by partial unfolding and rearrangement of the N-terminal tail, that led to steric hindrance of the pore. Our results indicate that the disordered properties of the N-terminus are crucial to properly account for the gating mechanism of VDAC.     

Dual Mechanism of Ion Permeation through VDAC Revealed with Inorganic Phosphate Ions and Phosphate Metabolites

In the exchange of metabolites and ions between the mitochondrion and the cytosol, the voltage-dependent anion channel (VDAC) is a key element, as it forms the major transport pathway for these compounds through the mitochondrial outer membrane. Numerous experimental studies have promoted the idea that VDAC acts as a regulator of essential mitochondrial functions. In this study, using a combination of molecular dynamics simulations, free-energy calculations, and electrophysiological measurements, we investigated the transport of ions through VDAC, with a focus on phosphate ions and metabolites. We showed that selectivity of VDAC towards small anions including monovalent phosphates arises from short-lived interactions with positively charged residues scattered throughout the pore. In dramatic contrast, permeation of divalent phosphate ions and phosphate metabolites (AMP and ATP) involves binding sites along a specific translocation pathway. This permeation mechanism offers an explanation for the decrease in VDAC conductance measured in the presence of ATP or AMP at physiological salt concentration. The binding sites occur at similar locations for the divalent phosphate ions, AMP and ATP, and contain identical basic residues. ATP features a marked affinity for a central region of the pore lined by two lysines and one arginine of the N-terminal helix. This cluster of residues together with a few other basic amino acids forms a “charged brush” which facilitates the passage of the anionic metabolites through the pore. All of this reveals that VDAC controls the transport of the inorganic phosphates and phosphate metabolites studied here through two different mechanisms.     

Effects of signal peptide and adenylate on the oligomerization and membrane binding of soluble SecA

SecA protein, a cytoplasmic ATPase, plays a central role in the secretion of signal peptide-containing proteins. Here, we examined effects of signal peptide and ATP on the oligomerization, conformational change, and membrane binding of SecA. The wild-type (WT) signal peptide from the ribose-binding protein inhibited ATP binding to soluble SecA and stimulated release of ATP already bound to the protein. The signal peptide enhanced the oligomerization of soluble SecA, while ATP induced dissociation of SecA oligomer. Analysis of SecA unfolding with urea or heat revealed that the WT signal peptide induces an open conformation of soluble SecA, while ATP increased the compactness of SecA. We further obtained evidences that the signal peptide-induced oligomerization and the formation of open structure enhance the membrane binding of SecA, whereas ATP inhibits the interaction of soluble SecA with membranes. On the other hand, the complex of membrane-bound SecA and signal peptide was shown to resume nucleotide-binding activity. From these results, we propose that the translocation components affect the degree of oligomerization of soluble SecA, thereby modulating the membrane binding of SecA in early translocation pathway. A possible sequential interaction of SecA with signal peptide, ATP, and cytoplasmic membrane is discussed.     

Affixing N-terminal α-helix to the wall of the voltage-dependent anion channel does not prevent its voltage gating

The voltage-dependent anion channel (VDAC) governs the free exchange of ions and metabolites between the mitochondria and the rest of the cell. The three-dimensional structure of VDAC1 reveals a channel formed by 19 β-strands and an N-terminal α-helix located near the midpoint of the pore. The position of this α-helix causes a narrowing of the cavity, but ample space for metabolite passage remains. The participation of the N-terminus of VDAC1 in the voltage-gating process has been well established, but the molecular mechanism continues to be debated; however, the majority of models entail large conformational changes of this N-terminal segment. Here we report that the pore-lining N-terminal α-helix does not undergo independent structural rearrangements during channel gating. We engineered a double Cys mutant in murine VDAC1 that cross-links the α-helix to the wall of the β-barrel pore and reconstituted the modified protein into planar lipid bilayers. The modified murine VDAC1 exhibited typical voltage gating. These results suggest that the N-terminal α-helix is located inside the pore of VDAC in the open state and remains associated with β-strand 11 of the pore wall during voltage gating.     

Markov state modeling of membrane transport proteins

The flux of ions and molecules in and out of the cell is vital for maintaining the basis of various biological processes. The permeation of substrates across the cellular membrane is mediated through the function of specialized integral membrane proteins commonly known as membrane transporters. These proteins undergo a series of structural rearrangements that allow a primary substrate binding site to be accessed from either side of the membrane at a given time. Structural insights provided by experimentally resolved structures of membrane transporters have aided in the biophysical characterization of these important molecular drug targets. However, characterizing the transitions between conformational states remains challenging to achieve both experimentally and computationally. Though molecular dynamics simulations are a powerful approach to provide atomistic resolution of protein dynamics, a recurring challenge is its ability to efficiently obtain relevant timescales of large conformational transitions as exhibited in transporters. One approach to overcome this difficulty is to adaptively guide the simulation to favor exploration of the conformational landscape, otherwise known as adaptive sampling. Furthermore, such sampling is greatly benefited by the statistical analysis of Markov state models. Historically, the use of Markov state models has been effective in quantifying slow dynamics or long timescale behaviors such as protein folding. Here, we review recent implementations of adaptive sampling and Markov state models to not only address current limitations of molecular dynamics simulations, but to also highlight how Markov state modeling can be applied to investigate the structure–function mechanisms of large, complex membrane transporters.     

Role of glycosylation on the ensemble of conformations in the MUC1 immunodominant epitope

MUC1 is a membrane glycoprotein, which in adenocarninomas is overexpressed and exhibits truncated O‐glycosylation. Overexpression and altered glycosylation make MUC1 into a candidate for immunotherapy. Monoclonal antibodies directed against MUC1 frequently bind an immunodominant epitope that contains a single site for O‐glycosylation. Glycosylation with tumor carbohydrate antigens such as the Tn‐antigen (GalNAc‐O‐Ser/Thr) results in antibodies binding with higher affinity. One proposed model to explain the enhanced affinity of antibodies for the glycosylated antigen is that the addition of a carbohydrate alters the conformational properties, favoring a binding‐competent state. The conformational effects associated with Tn glycosylation of the MUC1 antigen was investigated using solution‐state NMR and molecular dynamics. NMR experiments revealed distinct substructures of the glycosylated MUC1 peptides compared with the unglycosylated peptide. Molecular dynamics simulations of the MUC1 glycopeptide and peptide revealed distinguishing differences in their conformational preferences. Furthermore, the glycopeptide displayed a smaller conformational sampling compared with the peptide, suggesting that the glycopeptide sampled a narrower conformational space and is less dynamic. A comparison of the computed ensemble of conformations assuming random distribution, NMR models, and molecular dynamics simulations indicated that the MUC1 glycopeptide and aglycosylated peptide sampled structurally distinctly ensembles and that these ensembles were different from that of the random coil. Together, these data support the hypothesis that that conformational pre‐selection could be an essential feature of these peptides that dictates the binding affinities to MUC1 specific antibodies.     

Supramolecular assembly of VDAC in native mitochondrial outer membranes

The voltage-dependent anion channel (VDAC) is the most abundant protein in the mitochondrial outer membrane (MOM). Due to its localization, VDAC is involved in a wide range of processes, such as passage of ATP out of mitochondria, and particularly plays a central role in apoptosis. Importantly, the assembly of VDAC provides interaction with a wide range of proteins, some implying oligomerization. However, many questions remain as to the VDAC structure, its supramolecular assembly, packing density, and oligomerization in the MOM is unknown. Here we report the so far highest resolution view of VDAC and its native supramolecular assembly. We have studied yeast MOM by high-resolution atomic force microscopy (AFM) in physiological buffer and found VDAC in two distinct types of membrane domains. We found regions where VDAC was packed at high density (approximately 80%), rendering the membrane a voltage-dependent molecular sieve. In other domains, VDAC has a low surface density (approximately 20%) and the pore assembly ranges from single molecules to groups of up to 20. We assume that these groups are mobile in the lipid bilayer and allow association and dissociation with the large assemblies. VDAC has no preferred oligomeric state and no long-range order was observed in densely packed domains. High-resolution topographs show an eye-shaped VDAC with 3.8 nm x 2.7 nm pore dimensions. Based on the observed VDAC structure and the pair correlation function (PCF) analysis of the domain architectures, we propose a simple model that could explain the phase behavior of VDAC, and illustrates the sensitivity of the molecular organization to conditions in the cell, and the possibility for modulation of its assembly. The implication of VDAC in cytochrome c release from the mitochondria during cell apoptosis has made it a target in cancer research.     

Effect of Membrane Thickness on Conformational Sampling of Phospholamban from Computer Simulations

The conformational sampling of monomeric, membrane-bound phospholamban is described from computer simulations. Phospholamban (PLB) plays a key role as a regulator of sarcoplasmic reticulum calcium ATPase. An implicit membrane model is used in conjunction with replica exchange molecular dynamics simulations to reach μs-ms timescales. The implicit membrane model was also used to study the effect of different membrane thicknesses by scaling the low-dielectric region. The conformational sampling with the membrane model mimicking dipalmitoylphosphatidylcholine bilayers is in good agreement overall with experimental measurements, but consists of a wide variety of different conformations including structures not described previously. The conformational ensemble shifts significantly in the presence of thinner or thicker membranes. This has implications for the structure and dynamics of PLB in physiological membranes and offers what we believe to be a new interpretation of previous experimental measurements of PLB in detergents and microsomal membrane.     

Concentration dependent ion selectivity in VDAC: a molecular dynamics simulation study

The voltage-dependent anion channel (VDAC) forms the major pore in the outer mitochondrial membrane. Its high conducting open state features a moderate anion selectivity. There is some evidence indicating that the electrophysiological properties of VDAC vary with the salt concentration. Using a theoretical approach the molecular basis for this concentration dependence was investigated. Molecular dynamics simulations and continuum electrostatic calculations performed on the mouse VDAC1 isoform clearly demonstrate that the distribution of fixed charges in the channel creates an electric field, which determines the anion preference of VDAC at low salt concentration. Increasing the salt concentration in the bulk results in a higher concentration of ions in the VDAC wide pore. This event induces a large electrostatic screening of the charged residues promoting a less anion selective channel. Residues that are responsible for the electrostatic pattern of the channel were identified using the molecular dynamics trajectories. Some of these residues are found to be conserved suggesting that ion permeation between different VDAC species occurs through a common mechanism. This inference is buttressed by electrophysiological experiments performed on bean VDAC32 protein akin to mouse VDAC.     

The function of complexes between the outer mitochondrial membrane pore (VDAC) and the adenine nucleotide translocase in regulation of energy metabolism and apoptosis

The outer mitochondrial membrane pore (VDAC) changes its structure either voltage-dependently in artificial membranes or physiologically by interaction with the adenine nucleotide translocase (ANT) in the c-conformation. This interaction creates contact sites and leads in addition to a specific organisation of cytochrome c in the VDAC-ANT complexes. The VDAC structure that is specific for contact sites generates a signal at the surface for several proteins in the cytosol to bind with high capacity, such as hexokinase, glycerol kinase and Bax. If the VDAC binding site is not occupied by hexokinase, the VDAC-ANT complex has two critical qualities: firstly, Bax gets access to cytochrome c and secondly the ANT is set in its c-conformation that easily changes conformation into an unspecific channel (uniporter) causing permeability transition. Activity of bound hexokinase protects against both, it hinders Bax binding and employs the ANT as anti-porter. The octamer of mitochondrial creatine kinase binds to VDAC from the inner surface of the outer membrane. This firstly restrains interaction between VDAC and ANT and secondly changes the VDAC structure into low affinity for hexokinase and Bax. Cytochrome c in the creatine kinase complex will be differently organised, not accessible to Bax and the ANT is run as anti-porter by the active creatine kinase octamer. However, when, for example, free radicals cause dissociation of the octamer, VDAC interacts with the ANT with the same results as described above: Bax-dependent cytochrome c release and risk of permeability transition pore opening.     

The membrane-bound structure and topology of a human α-defensin indicate a dimer pore mechanism for membrane disruption

Defensins are cationic and disulfide-bonded host defense proteins of many animals that target microbial cell membranes. Elucidating the three-dimensional structure, dynamics, and topology of these proteins in phospholipid bilayers is important for understanding their mechanisms of action. Using solid-state nuclear magnetic resonance spectroscopy, we have now determined the conformation, dynamics, oligomeric state, and topology of a human α-defensin, HNP-1, in DMPC/DMPG bilayers. Two-dimensional correlation spectra show that membrane-bound HNP-1 exhibits a conformation similar to that of the water-soluble state, except for the turn connecting strands β2 and β3, whose side chains exhibit immobilization and conformational perturbation upon membrane binding. At high protein/lipid ratios, rapid (1)H spin diffusion from the lipid chains to the protein was observed, indicating that HNP-1 was well inserted into the hydrocarbon core of the bilayer. Arg Cζ-lipid (31)P distances indicate that only one of the four Arg residues forms tight hydrogen-bonded guanidinium-phosphate complexes. The protein is predominantly dimerized at high protein/lipid molar ratios, as shown by (19)F spin diffusion experiments. The presence of a small fraction of monomers and the shallower insertion at lower protein concentrations suggest that HNP-1 adopts concentration-dependent oligomerization and membrane-bound structure. These data strongly support a "dimer pore" topology of HNP-1 in which the polar top of the dimer lines an aqueous pore while the hydrophobic bottom faces the lipid chains. In this structure, R25 lies closest to the membrane surface among the four Arg residues. The pore does not have a high degree of lipid disorder, in contrast to the toroidal pores formed by protegrin-1, a two-stranded β-hairpin antimicrobial peptide. These results provide the first glimpse into the membrane-bound structure and mechanism of action of human α-defensins.     

Effect of Sec62 on the conformation of the Sec61 channel in yeast

Most eukaryotic secretory and membrane proteins are funneled by the Sec61 complex into the secretory pathway. Furthermore, some substrate peptides rely on two essential accessory proteins, Sec62 and Sec63, being present to assist with their translocation via the Sec61 channel in post-translational translocation. Cryo-electron microscopy (cryo-EM) recently succeeded in determining atomistic structures of unbound and signal sequence-engaged Sec complexes from Saccharomyces cerevisiae, involving the Sec61 channel and the proteins Sec62, Sec63, Sec71 and Sec72. In this study, we investigated the conformational effects of Sec62 on Sec61. Indeed, we observed in molecular dynamics simulations that the conformational dynamics of lateral gate, plug and pore region of Sec61 are altered by the presence/absence of Sec62. In molecular dynamics simulations that were started from the cryo-EM structures of Sec61 coordinated to Sec62 or of apo Sec61, we observed that the luminal side of the lateral gate gradually adopts a closed conformation similar to the apo state during unbound state simulations. In contrast, it adopts a wider conformation in the bound state. Furthermore, we demonstrate that the conformation of the active (substrate-bound) state of the Sec61 channel shifts toward an alternative conformation in the absence of the substrate. We suggest that the signal peptide holds/stabilizes the active state conformation of Sec61 during post-translational translocation. Thus, our study explains the effect of Sec62 on the conformation of the Sec61 channel and describes the conformational transitions of Sec61 channel.     

Protein grabs a ligand by extending anchor residues: molecular simulation for Ca2+ binding to calmodulin loop

The structural difference in proteins between unbound and bound forms directly suggests the importance of the conformational plasticity of proteins. However, pathways that connect two-end structures and how they are coupled to the binding reaction are not well understood at atomic resolution. Here, we analyzed the free-energy landscape, explicitly taking into account coupling between binding and conformational change by performing atomistic molecular dynamics simulations for Ca2+ binding to a calmodulin loop. Using the AMBER force field with explicit water solvent, we conducted umbrella sampling for the free-energy surface and steered molecular dynamics for the pathway search. We found that, at an early stage of binding, some key residue side chains extend their "arms" to catch Ca2+ and, after catching, they carry the Ca2+ to the center of the binding pocket. This grabbing motion resulted in smooth and stepwise exchange in coordination partners of Ca2+ from water oxygen to atoms in the calmodulin loop. The key residue that first caught the ion was one of the two acidic residues, which are highly conserved. In the pathway simulations, different pathways were observed between binding and dissociation reactions: The former was more diverse than the latter.     

Conformational Dynamics of a Ligand-Free Adenylate Kinase

Adenylate kinase (AdK) is a phosphoryl-transfer enzyme with important physiological functions. Based on a ligand-free open structure and a ligand-bound closed structure solved by crystallography, here we use molecular dynamics simulations to examine the stability and dynamics of AdK conformations in the absence of ligands. We first perform multiple simulations starting from the open or the closed structure, and observe their free evolutions during a simulation time of 100 or 200 nanoseconds. In all seven simulations starting from the open structure, AdK remained stable near the initial conformation. The eight simulations initiated from the closed structure, in contrast, exhibited large variation in the subsequent evolutions, with most (seven) undergoing large-scale spontaneous conformational changes and approaching or reaching the open state. To characterize the thermodynamics of the transition, we propose and apply a new sampling method that employs a series of restrained simulations to calculate a one-dimensional free energy along a curved pathway in the high-dimensional conformational space. Our calculated free energy profile features a single minimum at the open conformation, and indicates that the closed state, with a high (∼13 kcal/mol) free energy, is not metastable, consistent with the observed behaviors of the unrestrained simulations. Collectively, our simulations suggest that it is energetically unfavorable for the ligand-free AdK to access the closed conformation, and imply that ligand binding may precede the closure of the enzyme.     

VDAC and peripheral channelling complexes in health and disease

VDAC changes its structure either voltage dependent in artificial membranes or physiologically by interaction with the c conformation of the adenine nucleotide translocator (ANT). This interaction creates contact sites and leads to a specific organisation of cytochrome c in the VDAC ANT complexes. The VDAC structure specific for contact sites thus generates a signal at the surface for several proteins in the cytosol to bind with high affinity such as hexokinase, glycerolkinase and Bax. If the VDAC binding site is not occupied by hexokinase, the VDAC ANT complex has two critical qualities: firstly, external Bax gets access to the cytochrome c and secondly the ANT stays in the c conformation that easily changes the structure to an unspecific uni-porter causing permeability transition. Activity of bound hexokinase protects against both, it hinders Bax binding and employs the ANT as specific anti-porter. The octamer of mitochondrial creatine kinase binds to VDAC from the inner surface of the outer membrane. This firstly hinders direct interaction between VDAC and ANT and secondly changes porin structure into low affinity for hexokinase and external Bax. Cytochrome c in the creatine kinase complex will be differently organised not accessible to Bax and the ANT is run as anti-porter by the active octamer. However, when free radicals cause dissociation of the octamer, VDAC interacts with the ANT with the same results as described above: Bax dependent cytochrome c release and risk of permeability transition pore opening.     

Motions of the SecA protein motor bound to signal peptide: Insights from molecular dynamics simulations

SecA is an essential part of the Sec pathway for protein secretion in bacteria. In this pathway, SecA interacts with the N-terminal fragment of the secretory protein – the signal peptide, and couples binding and hydrolysis of adenosine triphosphate with movement of the secretory protein across the SecY protein translocon. How interactions with the signal peptide alter the conformational dynamics and long-distance conformational couplings of SecA is a key open question that we address here with molecular dynamics techniques. Analyses of protein motions indicate that the signal peptide alters SecA dynamics not only at the site where this peptide binds, but also at a nucleotide-binding domain. Hydrogen bond clusters contribute to the long-distance propagation of changes in SecA dynamics.     

Insight into the interactive residues between two domains of human somatic Angiotensin-converting enzyme and Angiotensin II by MM-PBSA calculation and steered molecular dynamics simulation

Angiotensin-converting enzyme (ACE), a membrane-bound zinc metallopeptidase, catalyzes the formation of Angiotensin-II (AngII) and the deactivation of bradykinin in the renin–angiotensin-aldosterone and kallikrein–kinin systems. As a hydrolysis product of ACE, AngII is regarded as an inhibitor and displays stronger competitive inhibition in the C-domain than the N-domain of ACE. However, the AngII binding differences between the two domains and the mechanisms behind AngII dissociation from the C-domain are rarely explored. In this work, molecular docking, Molecular Mechanics/Poisson–Boltzmann Surface Area calculation, and steered molecular dynamics (SMD) are applied to explore the structures and interactions in the binding or unbinding of AngII with the two domains of human somatic ACE. Calculated free energy values suggest that the C-domain–AngII complex is more stable than the N-domain–AngII complex, consistent with available experimental data. SMD simulation results imply that electrostatic interaction is dominant in the dissociation of AngII from the C-domain. Moreover, Gln106, Asp121, Glu123, and Tyr213 may be the key residues in the unbinding pathway of AngII. The simulation results in our work provide insights into the interactions between the two domains of ACE and its natural peptide inhibitor AngII at a molecular level. Moreover, the results provide theoretical clues for the design of new inhibitors.     

Toroidal pores formed by antimicrobial peptides show significant disorder

A large variety of antimicrobial peptides have been shown to act, at least in vitro, by poration of the lipid membrane. The nanometre size of these pores, however, complicates their structural characterization by experimental techniques. Here we use molecular dynamics simulations, to study the interaction of a specific class of antimicrobial peptides, melittin, with a dipalmitoylphosphatidylcholine bilayer in atomic detail. We show that transmembrane pores spontaneously form above a critical peptide to lipid ratio. The lipid molecules bend inwards to form a toroidally shaped pore but with only one or two peptides lining the pore. This is in strong contrast to the traditional models of toroidal pores in which the peptides are assumed to adopt a transmembrane orientation. We find that peptide aggregation, either prior or after binding to the membrane surface, is a prerequisite to pore formation. The presence of a stable helical secondary structure of the peptide, however is not. Furthermore, results obtained with modified peptides point to the importance of electrostatic interactions in the poration process. Removing the charges of the basic amino-acid residues of melittin prevents pore formation. It was also found that in the absence of counter ions pores not only form more rapidly but lead to membrane rupture. The rupture process occurs via a novel recursive poration pathway, which we coin the Droste mechanism.     

Investigation of pathways for the low-pH conformational transition in influenza hemagglutinin

Targeted molecular dynamics simulations were used to study the conformational transition of influenza hemagglutinin (HA) from the native conformation to putative fusogenic or postfusion conformations populated at low pH. Three pathways for this conformational change were considered. Complete dissociation of the globular domains of HA was observed in one pathway, whereas smaller rearrangements were observed in the other two. The fusion peptides became exposed and moved toward the target membrane, although occasional movement toward the viral membrane was also observed. The effective energy profiles along the paths show multiple barriers. The final low-pH structures, which are consistent with available experimental data, are comparable in effective energy to native HA. As a control, the uncleaved precursor HA0 was also forced along the same pathway. In this case both the final energy and the energy barrier were much higher than in the cleaved protein. This study suggests that 1) as proposed, the native conformation is the global minimum energy conformation for the uncleaved precursor but a metastable state for cleaved HA; 2) the spring-loaded conformational change is energetically plausible in full-length HA; and 3) complete globular domain dissociation is not necessary for extension of the coiled coil and fusion peptide exposure, but the model with complete dissociation has lower energy.