Therapeutic Peptides Are Preferentially Solubilized in Specific Microenvironments within PEG-PLGA Polymer Nanoparticles.

Polymeric nanoparticles are a highly promising drug delivery formulation. However, a lack of understanding of the molecular mechanisms that underlie their drug solubilization and controlled release capabilities has hindered the efficient clinical translation of such technologies. Polyethylene glycol-poly(lactic-co-glycolic) acid (PEG-PLGA) nanoparticles have been widely studied as cancer drug delivery vehicles. In this letter, we use unbiased coarse-grained molecular dynamics simulations to model the self-assembly of a PEG-PLGA nanoparticle and its solubulization of the anticancer peptide, EEK, with good agreement with previously reported experimental structural data. We applied unsupervised machine learning techniques to quantify the conformations that polymers adopt at various locations within the nanoparticle. We find that the local microenvironments formed by the various polymer conformations promote preferential EEK solubilization within specific regions of the NP. This demonstrates that these microenvironments are key in controlling drug storage locations within nanoparticles, supporting the rational design of nanoparticles for therapeutic applications.

Rapid Peptide Cyclization Inspired by the Modular Logic of Nonribosomal Peptide Synthetases.

Nonribosomal cyclic peptides (NRcPs) are structurally complex natural products and a vital pool of therapeutics, particularly antibiotics. Their structural diversity arises from the ability of the multidomain enzyme assembly lines, nonribosomal peptide synthetases (NRPSs), to utilize bespoke nonproteinogenic amino acids, modify the linear peptide during elongation, and catalyze an array of cyclization modes, e.g., head to tail, side chain to tail. The study and drug development of NRcPs are often limited by a lack of easy synthetic access to NRcPs and their analogues, with selective macrolactamization being a major bottleneck. Herein, we report a generally applicable chemical macrocyclization method of unprecedented speed and selectivity. Inspired by biosynthetic cyclization, it combines the deprotected linear biosynthetic precursor peptide sequence with a highly reactive C-terminus to produce NRcPs and analogues in minutes. The method was applied to several NRcPs of varying sequences, ring sizes, and cyclization modes including rufomycin, colistin, and gramicidin S with comparable success. We thus demonstrate that the linear order of modules in NRPS enzymes that determines peptide sequence encodes the key structural information to produce peptides conformationally biased toward macrocyclization. To fully exploit this conformational bias synthetically, a highly reactive C-terminal acyl azide is also required, alongside carefully balanced pH and solvent conditions. This allows for consistent, facile cyclization of exceptional speed, selectivity, and atom efficiency. This exciting macrolactamization method represents a new enabling technology for the biosynthetic study of NRcPs and their development as therapeutics.

A least-squares-fitting procedure for an efficient preclinical ranking of passive transport across the blood-brain barrier endothelium.

The treatment of various disorders of the central nervous system (CNS) is often impeded by the limited brain exposure of drugs, which is regulated by the human blood-brain barrier (BBB). The screening of lead compounds for CNS penetration is challenging due to the biochemical complexity of the BBB, while experimental determination of permeability is not feasible for all types of compounds. Here we present a novel method for rapid preclinical screening of libraries of compounds by utilizing advancements in computing hardware, with its foundation in transition-based counting of the flux. This method has been experimentally validated for in vitro permeabilities and provides atomic-level insights into transport mechanisms. Our approach only requires a single high-temperature simulation to rank a compound relative to a library, with a typical simulation time converging within 24 to 72 h. The method offers unbiased thermodynamic and kinetic information to interpret the passive transport of small-molecule drugs across the BBB.

Development of membrane-active peptide therapeutics in oncology.

Membrane-active peptides play an essential role in many living organisms and their immune systems and counter many infectious diseases. Many have dual or multiple mechanisms and can synergize with other molecules, like peptides, proteins, and small molecules. Although membrane-active peptides have been intensively studied in the past decades and more than 3500 sequences have been identified, only a few received approvals from the US Food and Drug Administration. In this review, we investigated all the peptide therapeutics that have entered the market or were subjected to preclinical and clinical studies to understand how they succeeded. With technological advancement (e.g., chemical modifications and pharmaceutical formulations) and a better understanding of the mechanism of action and the potential targets, we found at least five membrane-active peptide drugs that have entered preclinical/clinical phases and show promising results for cancer treatment. We summarized our findings in this review and provided insights into membrane-active anticancer peptide therapeutics.

Tuning of a Membrane-Perforating Antimicrobial Peptide to Selectively Target Membranes of Different Lipid Composition.

The use of designed antimicrobial peptides as drugs has been impeded by the absence of simple sequence-structure-function relationships and design rules. The likely cause is that many of these peptides permeabilize membranes via highly disordered, heterogeneous mechanisms, forming aggregates without well-defined tertiary or secondary structure. We suggest that the combination of high-throughput library screening with atomistic computer simulations can successfully address this challenge by tuning a previously developed general pore-forming peptide into a selective pore-former for different lipid types. A library of 2916 peptides was designed based on the LDKA template. The library peptides were synthesized and screened using a high-throughput orthogonal vesicle leakage assay. Dyes of different sizes were entrapped inside vesicles with varying lipid composition to simultaneously screen for both pore size and affinity for negatively charged and neutral lipid membranes. From this screen, nine different LDKA variants that have unique activity were selected, sequenced, synthesized, and characterized. Despite the minor sequence changes, each of these peptides has unique functional properties, forming either small or large pores and being selective for either neutral or anionic lipid bilayers. Long-scale, unbiased atomistic molecular dynamics (MD) simulations directly reveal that rather than rigid, well-defined pores, these peptides can form a large repertoire of functional dynamic and heterogeneous aggregates, strongly affected by single mutations. Predicting the propensity to aggregate and assemble in a given environment from sequence alone holds the key to functional prediction of membrane permeabilization.

Membrane proteins bind lipids selectively to modulate their structure and function.

Previous studies have established that the folding, structure and function of membrane proteins are influenced by their lipid environments and that lipids can bind to specific sites, for example, in potassium channels. Fundamental questions remain however regarding the extent of membrane protein selectivity towards lipids. Here we report a mass spectrometry approach designed to determine the selectivity of lipid binding to membrane protein complexes. We investigate the mechanosensitive channel of large conductance (MscL) from Mycobacterium tuberculosis and aquaporin Z (AqpZ) and the ammonia channel (AmtB) from Escherichia coli, using ion mobility mass spectrometry (IM-MS), which reports gas-phase collision cross-sections. We demonstrate that folded conformations of membrane protein complexes can exist in the gas phase. By resolving lipid-bound states, we then rank bound lipids on the basis of their ability to resist gas phase unfolding and thereby stabilize membrane protein structure. Lipids bind non-selectively and with high avidity to MscL, all imparting comparable stability; however, the highest-ranking lipid is phosphatidylinositol phosphate, in line with its proposed functional role in mechanosensation. AqpZ is also stabilized by many lipids, with cardiolipin imparting the most significant resistance to unfolding. Subsequently, through functional assays we show that cardiolipin modulates AqpZ function. Similar experiments identify AmtB as being highly selective for phosphatidylglycerol, prompting us to obtain an X-ray structure in this lipid membrane-like environment. The 2.3 Å resolution structure, when compared with others obtained without lipid bound, reveals distinct conformational changes that re-position AmtB residues to interact with the lipid bilayer. Our results demonstrate that resistance to unfolding correlates with specific lipid-binding events, enabling a distinction to be made between lipids that merely bind from those that modulate membrane protein structure and/or function. We anticipate that these findings will be important not only for defining the selectivity of membrane proteins towards lipids, but also for understanding the role of lipids in modulating protein function or drug binding.

Simulation-Guided Rational de Novo Design of a Small Pore-Forming Antimicrobial Peptide.

In the age of failing small-molecule antibiotics, tapping the near-infinite structural and chemical repertoire of antimicrobial peptides (AMPs) offers one of the most promising routes toward developing next-generation antibacterial compounds. One of the key impediments en route is the lack of methodologies for systematic rational design and optimization of new AMPs. Here we present a new simulation-guided rational design approach and apply it to develop a potent new AMP. We show that unbiased atomic detail molecular dynamics (MD) simulations are able to predict structures formed by evolving peptide designs enabling structure-based rational fine-tuning of functional properties. Starting from a 14-residue poly leucine template we demonstrate the design of a minimalistic potent new AMP. Consisting of only four types of amino acids (LDKA), this peptide forms large pores in microbial membranes at very low peptide-to-lipid ratios (1:1000) and exhibits low micromolar activity against common Gram-positive and Gram-negative pathogenic bacteria. Remarkably, the four amino acids were sufficient to encode preferential poration of bacterial membranes with negligible damage to red blood cells at bactericidal concentrations. As the sequence is too short to span cellular membranes, pores are formed by stacking of channels in each bilayer leaflet.

Molecular Dynamics Simulations Are Redefining Our View of Peptides Interacting with Biological Membranes.

Ever since the first molecular mechanics computer simulations of biological molecules became possible, there has been the dream to study all complex biological phenomena in silico, simply bypassing the enormous experimental challenges and their associated costs. For this, two inherent requirements need to be met: First, the time scales achievable in simulations must reach up to the millisecond range and even longer. Second, the computational model must accurately reproduce what is measured experimentally. Despite some recent successes, the general consensus in the field to date has been that neither of these conditions have yet been met and that the dream will be realized, if at all, only in the distant future. In this Account, we show that this view is wrong; instead, we are actually in the middle of the in silico molecular dynamics (MD) revolution, which is reshaping how we think about protein function. The example explored in this Account is a recent advance in the field of membrane-active peptides (MAPs). MD simulations have succeeded in accurately capturing the process of peptide binding, folding, and partitioning into lipid bilayers as well as revealing how channels form spontaneously from polypeptide fragments and conduct ionic and other cargo across membranes, all at atomic resolution. These game-changing advances have been made possible by a combination of steadily advancing computational power, more efficient algorithms and techniques, clever accelerated sampling schemes, and thorough experimental verifications. The great advantage of MD is the spatial and temporal resolution, directly providing a molecular movie of a protein undergoing folding and cycling through a functional process. This is especially important for proteins with transitory functional states, such as pore-forming MAPs. Recent successes are demonstrated here for the large class of antimicrobial peptides (AMPs). These short peptides are an essential part of the nonadaptive immune system for many organisms, ubiquitous in nature, and of particular interest to the pharmaceutical industry in the age of rising bacterial resistance to conventional antibiotic treatments. Unlike integral membrane proteins, AMPs are sufficiently small to allow converged sampling with the unbiased high-temperature sampling methodology outlined here and are relatively easy to handle experimentally. At the same time, AMPs exhibit a wealth of complex and poorly understood interactions with lipid bilayers, which allow not only tuning and validation of the simulation methodology but also advancement of our knowledge of protein-lipid interactions at a fundamental level. Space constraints limit our discussion to AMPs, but the MD methodologies outlined here can be applied to all phenomena involving peptides in membranes, including cell-penetrating peptides, signaling peptides, viral channel forming peptides, and fusion peptides, as well as ab initio membrane protein folding and assembly. For these systems, the promise of MD simulations to predict the structure of channels and to provide complete-atomic-detail trajectories of the mechanistic processes underlying their biological functions appears to rapidly become a reality. The current challenge is to design joint experimental and computational benchmarks to verify and tune MD force fields. With this, MD will finally fulfill its promise to become an inexpensive, powerful, and easy-to-use tool providing atomic-detail insights to researchers as part of their investigations into membrane biophysics and beyond.

Role of the Interaction Motif in Maintaining the Open Gate of an Open Sodium Channel.

Voltage-gated sodium channels undergo transitions between open, closed, and inactivated states, enabling regulation of the translocation of sodium ions across membranes. A recently published crystal structure of the full-length prokaryotic NavMs crystal structure in the activated open conformation has revealed the presence of a novel motif consisting of an extensive network of salt bridges involving residues in the voltage sensor, S4-S5 linker, pore, and C-terminal domains. This motif has been proposed to be responsible for maintaining an open conformation that enables ion translocation through the channel. In this study, we have used long-time molecular dynamics calculations without artificial restraints to demonstrate that the interaction network of full-length NavMs indeed prevents a rapid collapse and closure of the gate, in marked difference to earlier studies of the pore-only construct in which the gate had to be restrained to remain open. Interestingly, a frequently discussed "hydrophobic gating" mechanism at nanoscopic level is also observed in our simulations, in which the discontinuous water wire close to the gate region leads to an energetic barrier for ion conduction. In addition, we demonstrate the effects of in silico mutations of several of the key residues in the motif on the open channel's stability and functioning, correlating them with existing functional studies on this channel and homologous disease-associated mutations in human sodium channels; we also examine the effects of truncating/removing the voltage sensor and C-terminal domains in maintaining an open gate.

Computed Free Energies of Peptide Insertion into Bilayers are Independent of Computational Method.

We show that the free energy of inserting hydrophobic peptides into lipid bilayer membranes from surface-aligned to transmembrane inserted states can be reliably calculated using atomistic models. We use two entirely different computational methods: high temperature spontaneous peptide insertion calculations as well as umbrella sampling potential-of-mean-force (PMF) calculations, both yielding the same energetic profiles. The insertion free energies were calculated using two different protein and lipid force fields (OPLS protein/united-atom lipids and CHARMM36 protein/all-atom lipids) and found to be independent of the simulation parameters. In addition, the free energy of insertion is found to be independent of temperature for both force fields. However, we find major difference in the partitioning kinetics between OPLS and CHARMM36, likely due to the difference in roughness of the underlying free energy surfaces. Our results demonstrate not only a reliable method to calculate insertion free energies for peptides, but also represent a rare case where equilibrium simulations and PMF calculations can be directly compared.

Correction to: Computed Free Energies of Peptide Insertion into Bilayers are Independent of Computational Method.

The original version of the article unfortunately contained an error in NIH support grant number RO1-GM74639 in the Acknowledgements section. The correct grant number is RO1-GM74637. This has been corrected with this erratum.

Transmembrane helices containing a charged arginine are thermodynamically stable.

Hydrophobic amino acids are abundant in transmembrane (TM) helices of membrane proteins. Charged residues are sparse, apparently due to the unfavorable energetic cost of partitioning charges into nonpolar phases. Nevertheless, conserved arginine residues within TM helices regulate vital functions, such as ion channel voltage gating and integrin receptor inactivation. The energetic cost of arginine in various positions along hydrophobic helices has been controversial. Potential of mean force (PMF) calculations from atomistic molecular dynamics simulations predict very large energetic penalties, while in vitro experiments with Sec61 translocons indicate much smaller penalties, even for arginine in the center of hydrophobic TM helices. Resolution of this conflict has proved difficult, because the in vitro assay utilizes the complex Sec61 translocon, while the PMF calculations rely on the choice of simulation system and reaction coordinate. Here we present the results of computational and experimental studies that permit direct comparison with the Sec61 translocon results. We find that the Sec61 translocon mediates less efficient membrane insertion of Arg-containing TM helices compared with our computational and experimental bilayer-insertion results. In the simulations, a combination of arginine snorkeling, bilayer deformation, and peptide tilting is sufficient to lower the penalty of Arg insertion to an extent such that a hydrophobic TM helix with a central Arg residue readily inserts into a model membrane. Less favorable insertion by the translocon may be due to the decreased fluidity of the endoplasmic reticulum (ER) membrane compared with pure palmitoyloleoyl-phosphocholine (POPC). Nevertheless, our results provide an explanation for the differences between PMF- and experiment-based penalties for Arg burial.

Mechanisms of Membrane Pore Formation by Amyloidogenic Peptides in Amyotrophic Lateral Sclerosis.

Using unbiased atomic-detailed molecular dynamics simulations, the C-terminal fragments of TDP-43 are observed to aggregate and form disordered-toroidal pores in a lipid bilayer. Cytotoxicity of TDP-43 may be inferred from the observation that the membrane pores catalyze lipid flip-flop between bilayer leaflets and conduct water at high rates.

Molecular dynamics of ion transport through the open conformation of a bacterial voltage-gated sodium channel.

The crystal structure of the open conformation of a bacterial voltage-gated sodium channel pore from Magnetococcus sp. (NaVMs) has provided the basis for a molecular dynamics study defining the channel's full ion translocation pathway and conductance process, selectivity, electrophysiological characteristics, and ion-binding sites. Microsecond molecular dynamics simulations permitted a complete time-course characterization of the protein in a membrane system, capturing the plethora of conductance events and revealing a complex mixture of single and multi-ion phenomena with decoupled rapid bidirectional water transport. The simulations suggest specific localization sites for the sodium ions, which correspond with experimentally determined electron density found in the selectivity filter of the crystal structure. These studies have also allowed us to identify the ion conductance mechanism and its relation to water movement for the NavMs channel pore and to make realistic predictions of its conductance properties. The calculated single-channel conductance and selectivity ratio correspond closely with the electrophysiology measurements of the NavMs channel expressed in HEK 293 cells. The ion translocation process seen in this voltage-gated sodium channel is clearly different from that exhibited by members of the closely related family of voltage-gated potassium channels and also differs considerably from existing proposals for the conductance process in sodium channels. These studies simulate sodium channel conductance based on an experimentally determined structure of a sodium channel pore that has a completely open transmembrane pathway and activation gate.

Lipid headgroups modulate membrane insertion of pHLIP peptide.

The pH low insertion peptide (pHLIP) is an important tool for drug delivery and visualization of acidic tissues produced by various maladies, including cancer, inflammation, and ischemia. Numerous studies indicate that pHLIP exists in three states: unfolded and soluble in water at neutral pH (State I), unfolded and bound to the surface of a phosphatidylcholine membrane at neutral pH (State II), and inserted across the membrane as an α-helix at low pH (State III). Here we report how changes in lipid composition modulate this insertion scheme. First, the presence of either anionic lipids, cholesterol, or phosphoethanolamine eliminates membrane binding at neutral pH (State II). Second, the apparent pKa for the insertion transition (State I → State III) is increased with increasing content of anionic lipids, suggesting that electrostatic interactions in the interfacial region modulate protonation of acidic residues of pHLIP responsible for transbilayer insertion. These findings indicate a possibility for triggering protonation-coupled conformational switching in proteins at membrane interfaces through changes in lipid composition.

Absorption and folding of melittin onto lipid bilayer membranes via unbiased atomic detail microsecond molecular dynamics simulation.

Unbiased molecular simulation is a powerful tool to study the atomic details driving functional structural changes or folding pathways of highly fluid systems, which present great challenges experimentally. Here we apply unbiased long-timescale molecular dynamics simulation to study the ab initio folding and partitioning of melittin, a template amphiphilic membrane active peptide. The simulations reveal that the peptide binds strongly to the lipid bilayer in an unstructured configuration. Interfacial folding results in a localized bilayer deformation. Akin to purely hydrophobic transmembrane segments the surface bound native helical conformer is highly resistant against thermal denaturation. Circular dichroism spectroscopy experiments confirm the strong binding and thermostability of the peptide. The study highlights the utility of molecular dynamics simulations for studying transient mechanisms in fluid lipid bilayer systems. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Computational modeling of membrane proteins.

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as ß-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.

Peptide Folding in Translocon-Like Pores.

The cellular translocon, present in all three domains of life, is one of the most versatile and important biological nanopores. This complex molecular apparatus is directly responsible for the secretion of globular proteins across membranes as well as the insertion of integral membrane proteins into lipid bilayers. Recently determined structures of the archaean SecY translocon reveal an hour-glass-shaped pore, which accommodates the nascent peptide chain during translocation. While these structures provide important insights into ribosome binding to the translocon, threading of the nascent chain into the channel, and lateral gate opening for releasing the folded helical peptide into the membrane bilayer, the exact folding pathway of the peptide inside the protein-conducting channel during translocation and prior to the lateral release into the bilayer remains elusive. In the present study, we use molecular dynamics simulations to investigate atomic resolution peptide folding in hour-glass-shaped pore models that are based on the SecY translocon channel structure. The theoretical setup allows systematic variation of key determinants of folding, in particular the degree of confinement of the peptide and the hydration level of the pore. A 27-residue hydrophobic peptide was studied that is preferentially inserted into membranes by the translocon. Our results show that both pore diameter as well as channel hydration are important determinants for folding efficiency and helical stability of the peptide, therefore providing important insights into translocon gating and lateral peptide partitioning.

Unraveling the Molecular Dance: Insights into TREM2/DAP12 Complex Formation in Alzheimer's Disease through Molecular Dynamics Simulations.

Alzheimer's disease (AD) is a widespread neurodegenerative condition affecting millions globally. Recent research has implicated variants of the triggering receptor expressed in myeloid cells 2 (TREM2) as risk factors for AD. TREM2, an immunomodulatory receptor on microglial surfaces, plays a pivotal role in regulating microglial activation by association with DNAX-activation protein 12 (DAP12). Despite its significance, the mechanism underlying the formation of the complex between the transmembrane domains (TMDs) of TREM2 and DAP12 remains unclear. This study employs multiscale molecular dynamics (MD) simulations to investigate three TMD complex models, including two derived from experiments and one generated by AlphaFold2. Conducted within a lipid membrane consisting of an 80:20 mixture of phosphatidylcholine (POPC) and cholesterol, our analysis reveals hydrogen-bonding interactions between K26 of TREM2 and D16 of DAP12 in all three models, consistent with previous experimental findings. Our results elucidate the different spatial conformations observed in the models and offer insights into the structure of the TREM2/DAP12 TMD complex. Furthermore, we elucidate the role of charged residues in the assembly structure of the complex within the lipid membrane. These findings enhance our understanding of the molecular mechanism governing TREM2/DAP12 complex formation, providing a foundation for designing novel therapeutic strategies to address AD and other neurodegenerative diseases.

Conformational states of melittin at a bilayer interface.

The distribution of peptide conformations in the membrane interface is central to partitioning energetics. Molecular-dynamics simulations enable characterization of in-membrane structural dynamics. Here, we describe melittin partitioning into dioleoylphosphatidylcholine lipids using CHARMM and OPLS force fields. Although the OPLS simulation failed to reproduce experimental results, the CHARMM simulation reported was consistent with experiments. The CHARMM simulation showed melittin to be represented by a narrow distribution of folding states in the membrane interface.