Molecular Mechanism for Gramicidin Dimerization and Dissociation in Bilayers of Different Thickness

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Peptide Power: Mechanistic Insights into the Effect of Mitochondria-Targeted Tetrapeptides on Membrane Electrostatics from Molecular Simulations

Mitochondrial dysfunction is implicated in nine of the ten leading causes of death in the US, yet there are no FDA-approved therapeutics to treat it. Synthetic mitochondria-targeted peptides (MTPs), including the lead compound SS-31, offer promise, as they have been shown to restore healthy mitochondrial function and treat a variety of common diseases. At the cellular level, research has shown that MTPs accumulate strongly at the inner mitochondrial membrane (IMM), slow energy sinks (e.g., proton leaks), and improve ATP production. Modulation of electrostatic fields around the IMM has been implicated as a key aspect in the mechanism of action (MoA) of these peptides; however, molecular and mechanistic details have remained elusive. In this study, we employed all-atom molecular dynamics simulations (MD) to investigate the interactions of four MTPs with lipid bilayers and calculate their effect on structural and electrostatic properties. In agreement with previous experimental findings, we observed the modulation of the membrane surface and dipole potentials by MTPs. The simulations reveal that the MTPs achieve a reduction in the dipole potential by acting to disorder both lipid head groups and water layers proximal to the bilayer surface. We also find that MTPs decrease the bilayer thickness and increase the membrane's capacitance. These changes suggest that MTPs may enhance how much potential energy can be stored across the IMM at a given transmembrane potential difference. The MTPs also displace cations away from the bilayer surface, modulating the surface potential and offering an alternative mechanism for how these MTPs reduce mitochondrial energy sinks like proton leaks and mitigate Ca2+ accumulation stress. In conclusion, this study highlights the therapeutic potential of MTPs and underlines how interactions of MTPs with lipid bilayers serve as a fundamental component of their MoA.

Screening for bilayer-active and likely cytotoxic molecules reveals bilayer-mediated regulation of cell function

A perennial problem encountered when using small molecules (drugs) to manipulate cell or protein function is to assess whether observed changes in function result from specific interactions with a desired target or from less specific off-target mechanisms. This is important in laboratory research as well as in drug development, where the goal is to identify molecules that are unlikely to be successful therapeutics early in the process, thereby avoiding costly mistakes. We pursued this challenge from the perspective that many bioactive molecules (drugs) are amphiphiles that alter lipid bilayer elastic properties, which may cause indiscriminate changes in membrane protein (and cell) function and, in turn, cytotoxicity. Such drug-induced changes in bilayer properties can be quantified as changes in the monomer↔dimer equilibrium for bilayer-spanning gramicidin channels. Using this approach, we tested whether molecules in the Pathogen Box (a library of 400 drugs and drug-like molecules with confirmed activity against tropical diseases released by Medicines for Malaria Venture to encourage the development of therapies for neglected tropical diseases) are bilayer modifiers. 32% of the molecules in the Pathogen Box were bilayer modifiers, defined as molecules that at 10 µM shifted the monomer↔dimer equilibrium toward the conducting dimers by at least 50%. Correlation analysis of the molecules' reported HepG2 cell cytotoxicity to bilayer-modifying potency, quantified as the shift in the gramicidin monomer↔dimer equilibrium, revealed that molecules producing <25% change in the equilibrium had significantly lower probability of being cytotoxic than molecules producing >50% change. Neither cytotoxicity nor bilayer-modifying potency (quantified as the shift in the gramicidin monomer↔dimer equilibrium) was well predicted by conventional physico-chemical descriptors (hydrophobicity, polar surface area, etc.). We conclude that drug-induced changes in lipid bilayer properties are robust predictors of the likelihood of membrane-mediated off-target effects, including cytotoxicity.

Protein-protein interface analysis of the non-ribosomal peptide synthetase peptidyl carrier protein and enzymatic domains

Non-ribosomal peptide synthetases (NRPSs) are attractive targets for biosynthetic pathway engineering due to their modular architecture and the therapeutic relevance of their products. With catalysis mediated by specific protein-protein interactions formed between the peptidyl carrier protein (PCP) and its partner enzymes, NRPS enzymology and control remains fertile ground for discovery. This review focuses on the recent efforts within structural biology by compiling high-resolution structural data that shed light into the various protein-protein interfaces formed between the PCP and its partner enzymes, including the phosphopantetheinyl transferase (PPTase), adenylation (A) domain, condensation (C) domain, thioesterase (TE) domain and other tailoring enzymes within the synthetase. Integrating our understanding of how the PCP recognizes partner proteins with the potential to use directed evolution and combinatorial biosynthetic methods will enhance future efforts in discovery and production of new bioactive compounds.

The coming of age of water channels for separation membranes: from biological to biomimetic to synthetic

Water channels are one of the key pillars driving the development of next-generation desalination and water treatment membranes. Over the past two decades, the rise of nanotechnology has brought together an abundance of multifunctional nanochannels that are poised to reinvent separation membranes with performances exceeding those of state-of-the-art polymeric membranes within the water-energy nexus. Today, these water nanochannels can be broadly categorized into biological, biomimetic and synthetic, owing to their different natures, physicochemical properties and methods for membrane nanoarchitectonics. Furthermore, against the backdrop of different separation mechanisms, different types of nanochannel exhibit unique merits and limitations, which determine their usability and suitability for different membrane designs. Herein, this review outlines the progress of a comprehensive amount of nanochannels, which include aquaporins, pillar[5]arenes, I-quartets, different types of nanotubes and their porins, graphene-based materials, metal- and covalent-organic frameworks, porous organic cages, MoS₂, and MXenes, offering a comparative glimpse into where their potential lies. First, we map out the background by looking into the evolution of nanochannels over the years, before discussing their latest developments by focusing on the key physicochemical and intrinsic transport properties of these channels from the chemistry standpoint. Next, we put into perspective the fabrication methods that can nanoarchitecture water channels into high-performance nanochannel-enabled membranes, focusing especially on the distinct differences of each type of nanochannel and how they can be leveraged to unlock the as-promised high water transport potential in current mainstream membrane designs. Lastly, we critically evaluate recent findings to provide a holistic qualitative assessment of the nanochannels with respect to the attributes that are most strongly valued in membrane engineering, before discussing upcoming challenges to share our perspectives with researchers for pathing future directions in this coming of age of water channels.

How to Combat Gram-Negative Bacteria Using Antimicrobial Peptides: A Challenge or an Unattainable Goal?

Antimicrobial peptides (AMPs) represent a promising and effective alternative for combating pathogens, having some advantages compared to conventional antibiotics. However, AMPs must also contend with complex and specialised Gram-negative bacteria envelops. The variety of lipopolysaccharide and phospholipid composition in Gram-negative bacteria strains and species are decisive characteristics regarding their susceptibility or resistance to AMPs. Such biological and structural barriers have created delays in tuning AMPs to deal with Gram-negative bacteria. This becomes even more acute because little is known about the interaction AMP–Gram-negative bacteria and/or AMPs’ physicochemical characteristics, which could lead to obtaining selective molecules against Gram-negative bacteria. As a consequence, available AMPs usually have highly associated haemolytic and/or cytotoxic activity. Only one AMP has so far been FDA approved and another two are currently in clinical trials against Gram-negative bacteria. Such a pessimistic panorama suggests that efforts should be concentrated on the search for new molecules, designs and strategies for combating infection caused by this type of microorganism. This review has therefore been aimed at describing the currently available AMPs for combating Gram-negative bacteria, exploring the characteristics of these bacteria’s cell envelop hampering the development of new AMPs, and offers a perspective regarding the challenges for designing new AMPs against Gram-negative bacteria.

Mechanisms underlying drug-mediated regulation of membrane protein function

The hydrophobic coupling between membrane proteins and their host lipid bilayer provides a mechanism by which bilayer-modifying drugs may alter protein function. Drug regulation of membrane protein function thus may be mediated by both direct interactions with the protein and drug-induced alterations of bilayer properties, in which the latter will alter the energetics of protein conformational changes. To tease apart these mechanisms, we examine how the prototypical, proton-gated bacterial potassium channel KcsA is regulated by bilayer-modifying drugs using a fluorescence-based approach to quantify changes in both KcsA function and lipid bilayer properties (using gramicidin channels as probes). All tested drugs inhibited KcsA activity, and the changes in the different gating steps varied with bilayer thickness, suggesting a coupling to the bilayer. Examining the correlations between changes in KcsA gating steps and bilayer properties reveals that drug-induced regulation of membrane protein function indeed involves bilayer-mediated mechanisms. Both direct, either specific or nonspecific, binding and bilayer-mediated mechanisms therefore are likely to be important whenever there is overlap between the concentration ranges at which a drug alters membrane protein function and bilayer properties. Because changes in bilayer properties will impact many diverse membrane proteins, they may cause indiscriminate changes in protein function.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

Molecular Dynamics for Antimicrobial Peptide Discovery

Although antimicrobial resistance is an increasingly significant public health concern, there have only been two new classes of antibiotics approved for human use since the 1960s. Understanding the mechanisms of action of antibiotics is critical for novel antibiotic discovery, but novel approaches are needed that do not exclusively rely on experiments. Molecular dynamics simulation is a computational tool that uses simple models of the atoms in a system to discover nanoscale insights into the dynamic relationship between mechanism and biological function. Such insights can lay the framework for elucidating the mechanism of action and optimizing antibiotic templates. Antimicrobial peptides represent a promising solution to escalating antimicrobial resistance, given their lesser tendency to induce resistance than that of small-molecule antibiotics. Simulations of these agents have already revealed how they interact with bacterial membranes and the underlying physiochemical features directing their structure and function. In this minireview, we discuss how traditional molecular dynamics simulation works and its role and potential for the development of new antibiotic candidates with an emphasis on antimicrobial peptides.

Atomistic Characterization of Gramicidin Channel Formation

We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA-gA stacking orientations that were stabilized by four or two hydrogen bonds (4HB or 2HB). When examining the temporal evolution of the dimerization, we found that two bilayer-inserted gA subunits can form the 6HB dimer directly, with no discernible intermediate states, as well as through paths that involve the 2HB and 4HB dimers.

Deciphering the Structure of the Gramicidin A Channel in the Presence of AOT Reverse Micelles in Pentane Using Molecular Dynamics Simulations

Structural studies of proteins and, in particular, integral membrane proteins (IMPs) using solution NMR spectroscopy approaches are challenging due to not only their inherent structural complexities but also the fact that they need to be solubilized in biomimetic environments (such as micelles), which enhances the slow molecular reorientation. To deal with these difficulties and increase the effective rate of molecular reorientation, the encapsulation of IMPs in the aqueous core of the reverse micelle (RM) dissolved in a low-viscosity solvent has been proven to be a viable approach. However, the effect of the reverse micelle (RM) environment on the IMP structure and function is little known. To gain insight into these aspects, this article presents a series of atomistic unconstrained molecular dynamics (MD) of a model ion channel (gramicidin A, gA) with RMs formed with anionic surfactant diacyl chain bis(2-ethylhexyl) sodium succinate (AOT) in pentane at a water-to-surfactant molar ratio (W₀) of 6. The simulations were carried out with different protocols and starting conditions for a total of 2.4 μs and were compared with other MDs used with the gA channel inserted in models of the SDS micelle or the DMPC membrane. We show here that in the presence of AOT RMs the gA dimer did not look like the "dumbbell-like" model anticipated by experiments, where the C-terminal parts of the gA are capped with two RMs and the rest of the dimer is protected from the oil solvent by the AOT acyl chains. In contrast, the MD simulations reveal that the AOT, Na⁺, and water formed two well-defined and elongated RMs attached to the C-terminal ends of the gA dimer, while the rest is in direct contact with the pentane. The initial β^6.3 secondary structure of the gA is well conserved and filled with 6-9 waters, as in SDS micelles or the DMPC membrane. Finally, the water movement inside the gA is strongly affected by the presence of RMs at each extremity, and no passage of water molecules through the gA channel is observed even after a long simulation period, whereas the opposite was found for gA in SDS and DMPC.

Assessing the Perturbing Effects of Drugs on Lipid Bilayers Using Gramicidin Channel-Based In Silico and In Vitro Assays

Partitioning of bioactive molecules, including drugs, into cell membranes may produce indiscriminate changes in membrane protein function. As a guide to safe drug development, it therefore becomes important to be able to predict the bilayer-perturbing potency of hydrophobic/amphiphilic drugs candidates. Toward this end, we exploited gramicidin channels as molecular force probes and developed in silico and in vitro assays to measure drugs' bilayer-modifying potency. We examined eight drug-like molecules that were found to enhance or suppress gramicidin channel function in a thick 1,2-dierucoyl-sn-glycero-3-phosphocholine (DC22:1PC) but not in thin 1,2-dioleoyl-sn-glycero-3-phosphocholine (DC18:1PC) lipid bilayer. The mechanism underlying this difference was attributable to the changes in gramicidin dimerization free energy by drug-induced perturbations of lipid bilayer physical properties and bilayer-gramicidin interactions. The combined in silico and in vitro approaches, which allow for predicting the perturbing effects of drug candidates on membrane protein function, have implications for preclinical drug safety assessment.

Shorter Antibacterial Peptide Having High Selectivity for E. coli Membranes and Low Potential for Inducing Resistance

Antimicrobial peptides (AMPs) have been recognised as a significant therapeutic option for mitigating resistant microbial infections. It has been found recently that Plasmodium falciparum-derived, 20 residue long, peptide 35409 had antibacterial and haemolytic activity, making it an AMP having reduced selectivity, and suggesting that it should be studied more extensively for obtaining new AMPs having activity solely targeting the bacterial membrane. Peptide 35409 was thus used as template for producing short synthetic peptides (<20 residues long) and evaluating their biological activity and relevant physicochemical characteristics for therapeutic use. Four of the sixteen short peptides evaluated here had activity against E. coli without any associated haemolytic effects. The 35409-1 derivative (17 residues long) had the best therapeutic characteristics as it had high selectivity for bacterial cells, stability in the presence of human sera, activity against E. coli multiresistant clinical isolates and was shorter than the original sequence. It had a powerful membranolytic effect and low potential for inducing resistance in bacteria. This peptide’s characteristics highlighted its potential as an alternative for combating infection caused by E. coli multiresistant bacteria and/or for designing new AMPs.