Accelerated molecular dynamics simulation analysis of MSI-594 in a lipid bilayer

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.

Membrane-Active All-Hydrocarbon-Stapled α-Helical Amphiphilic Tat Peptides: Broad-Spectrum Antibacterial Activity and Low Incidence of Drug Resistance

Multidrug resistance against conventional antibiotics has dramatically increased the difficulty of treatment and accelerated the need for novel antibacterial agents. The peptide Tat (47-57) is derived from the transactivating transcriptional activator of human immunodeficiency virus 1, which is well-known as a cell-penetrating peptide in mammalian cells. However, it is also reported that the Tat peptide (47-57) has antifungal activity. In this study, a series of membrane-active hydrocarbon-stapled α-helical amphiphilic peptides were synthesized and evaluated as antibacterial agents against Gram-positive and Gram-negative bacteria, including multidrug-resistant strains. The impact of hydrocarbon staple, the position of aromatic amino acid residue in the hydrophobic face, the various types of aromatic amino acids, and the hydrophobicity on bioactivity were also investigated and discussed in this study. Among those synthesized peptides, analogues P3 and P10 bearing a l-2-naphthylalanine (Φ) residue at the first position and a Tyr residue at the eighth position demonstrated the highest antimicrobial activity and negligible hemolytic toxicity. Notably, P3 and P10 showed obviously enhanced antimicrobial activity against multidrug-resistant bacteria, low drug resistance, high cell selectivity, extended half-life in plasma, and excellent performance against biofilm. The antibacterial mechanisms of P3 and P10 were also preliminarily investigated in this effort. In conclusion, P3 and P10 are promising antimicrobial alternatives for the treatment of the antimicrobial-resistance crisis.

Insights into the Adsorption Mechanisms of the Antimicrobial Peptide CIDEM-501 on Membrane Models

CIDEM-501 is a hybrid antimicrobial peptide rationally designed based on the structure of panusin and panulirin template peptides. The new peptide exhibits significant antibacterial activity against multidrug-resistant pathogens (MIC = 2-4 μM) while conserving no toxicity in human cell lines. We conducted molecular dynamics (MD) simulations using the CHARMM-36 force field to explore the CIDEM-501 adsorption mechanism with different membrane compositions. Several parameters that characterize these interactions were analyzed to elucidate individual residues' structural and thermodynamic contributions. The membrane models were constructed using CHARMM-GUI, mimicking the bacterial and eukaryotic phospholipid compositions. Molecular dynamics simulations were conducted over 500 ns, showing rapid and highly stable peptide adsorption to bacterial lipids components rather than the zwitterionic eucaryotic model membrane. A predominant peptide orientation was observed in all models dominated by an electric dipole. The peptide remained parallel to the membrane surface with the center loop oriented to the lipids. Our findings shed light on the antibacterial activity of CIDEM-501 on bacterial membranes and yield insights valuable for designing potent antimicrobial peptides targeting multi- and extreme drug-resistant bacteria.

Exploring Cyclic Peptide Nanotube Stability Across Diverse Lipid Bilayers and Unveiling Water Transport Dynamics

Cyclic Peptide Nanotubes (CPNTs) have emerged as compelling candidates for various applications, particularly as nanochannels within lipid bilayers. In this study, the stability of two CPNTs, namely 8 × [(Cys-Gly-Met-Gly)2] and 8 × [(Gly-Leu)4], are comprehensively investigated across different lipid bilayers, including 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a mixed model membrane (POPE/POPG), and a realistic yeast model membrane. The results demonstrate that both CPNTs maintain their tubular structures in all lipid bilayers, with [(Cys-Gly-Met-Gly)2] showing increased stability over an extended period in these lipid membranes. The insertion of CPNTs shows negligible impact on lipid bilayer properties, including area per lipid, volume per lipid, and bilayer thickness. The study demonstrates that the CPNT preserves its two-line water movement pattern within all the lipid membranes, reaffirming their potential as water channels. The MSD curves further reveal that the dynamics of water molecules inside the nanotube are similar for all the bilayer systems with minor differences that arise due to different lipid environments.

A novel synthetic peptide analog enhanced antibacterial activity of the frog-derived skin peptide wuchuanin-A1

In recent years, there has been a growing focus on the development of novel antibacterial compounds for clinical applications, such as antimicrobial peptide (AMP). Among the developed AMP, wuchuanin-A1, a coil-shaped bioactive peptide derived from Odorrana wuchuanensis frog skin, has been reported to exhibit antibacterial, antifungal, and antioxidant activity, but there are limited studies on its potential as an antibacterial agent. Therefore, this study aims to molecularly modify the sequence of wuchuanin-A1 to enhance its antibacterial properties. The interaction of both the native and analog peptide with bacterial inner membranes was initially assessed using computational methods. Specific amino acid substitutions were then used to enhance the modified peptide's antibacterial efficacy, followed by several preliminary tests to evaluate its activity. This study bridges the gap in exploring the potential of wuchuanin-A1 for antibacterial purposes, providing insights into the design of effective antimicrobial agents.Communicated by Ramaswamy H. Sarma.

Comparison, Analysis, and Molecular Dynamics Simulations of Structures of a Viral Protein Modeled Using Various Computational Tools

The structural analysis of proteins is a major domain of biomedical research. Such analysis requires resolved three-dimensional structures of proteins. Advancements in computer technology have led to progress in biomedical research. In silico prediction and modeling approaches have facilitated the construction of protein structures, with or without structural templates. In this study, we used three neural network-based de novo modeling approaches-AlphaFold2 (AF2), Robetta-RoseTTAFold (Robetta), and transform-restrained Rosetta (trRosetta)-and two template-based tools-the Molecular Operating Environment (MOE) and iterative threading assembly refinement (I-TASSER)-to construct the structure of a viral capsid protein, hepatitis C virus core protein (HCVcp), whose structure have not been fully resolved by laboratory techniques. Templates with sufficient sequence identity for the homology modeling of complete HCVcp are currently unavailable. Therefore, we performed domain-based homology modeling for MOE simulations. The templates for each domain were obtained through sequence-based searches on NCBI and the Protein Data Bank. Then, the modeled domains were assembled to construct the complete structure of HCVcp. The full-length structure and two truncated forms modeled using various computational tools were compared. Molecular dynamics (MD) simulations were performed to refine the structures. The root mean square deviation of backbone atoms, root mean square fluctuation of Cα atoms, and radius of gyration were calculated to monitor structural changes and convergence in the simulations. The model quality was evaluated through ERRAT and phi-psi plot analysis. In terms of the initial prediction for protein modeling, Robetta and trRosetta outperformed AF2. Regarding template-based tools, MOE outperformed I-TASSER. MD simulations resulted in compactly folded protein structures, which were of good quality and theoretically accurate. Thus, the predicted structures of certain proteins must be refined to obtain reliable structural models. MD simulation is a promising tool for this purpose.

Conformation and Orientation of Antimicrobial Peptides MSI-594 and MSI-594A in a Lipid Membrane

There is significant interest in the development of antimicrobial compounds to overcome the increasing bacterial resistance to conventional antibiotics. Studies have shown that naturally occurring and de novo-designed antimicrobial peptides could be promising candidates. MSI-594 is a synthetic linear, cationic peptide that has been reported to exhibit a broad spectrum of antimicrobial activities. Investigation into how MSI-594 disrupts the cell membrane is important for better understanding the details of this antimicrobial peptide (AMP)'s action against bacterial cells. In this study, we used two different synthetic lipid bilayers: zwitterionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and anionic 7:3 POPC/1-palmitoyl-2-oleoyl-sn-glycero-3-phospho(1'-rac-glycerol) (POPG). Sum frequency generation (SFG) vibrational spectroscopy and attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR) were used to determine the orientations of MSI-594 and its analogue MSI-594A associated with zwitterionic POPC and anionic 7:3 POPC/POPG lipid bilayers. The simulated ATR-FTIR and SFG spectra using nuclear magnetic resonance (NMR)-determined structures were compared with experimental spectra to optimize the bent angle between the N- (1-11) and C- (12-24) termini helices and the membrane orientations of the helices; since the NMR structure of the peptide was determined from lipopolysaccharide (LPS) micelles, the optimization was needed to find the most suitable conformation and orientation in lipid bilayers. The reported experimental results indicate that the optimized MSI-594 helical hairpin structure adopts a complete lipid bilayer surface-bound orientation (denoted "face-on") in both POPC and 7:3 POPC/POPG lipid bilayers. The analogue peptide, MSI-584A, on the other hand, exhibited a larger bent angle between the N- (1-11) and C- (12-24) termini helices with the hydrophobic C-terminal helix inserted into the hydrophobic region of the bilayer (denoted "membrane-inserted") when interacting with both POPC and 7:3 POPC/POPG lipid bilayers. These experimental findings on the membrane orientations suggest that both peptides are likely to disrupt the cell membrane through the carpet mechanism.

The role of hydrophobic patches of de novo designed MSI-78 and VG16KRKP antimicrobial peptides on fragmenting model bilayer membranes

Antimicrobial peptides (AMPs) with cell membrane lysing capability are considered potential candidates for the development of the next generation of antibiotics. Designing novel AMPs requires an in-depth understanding of the mechanism of action of the peptides. In this work, we used various biophysical techniques including ³¹P solid-state NMR to examine the interaction of model membranes with amphipathic de novo-designed peptides. Two such peptides, MSI-78 and VG16KRKP, were designed with different hydrophobicity and positive charges. The model lipid membranes were constituted by mixing lipids of varying degrees of 'area per lipid' (APL), which directly affected the packing properties of the membrane. The observed emergence of the isotropic peak in ³¹P NMR spectra as a function of time is a consequence of the fragmentation of the membrane mediated by the peptide interaction. The factors such as the charges, overall hydrophilicity of the AMPs, as well as lipid membrane packing, contributed to the kinetics of membrane fragmentation. Furthermore, we anticipate the designed AMPs follow the carpet and toroidal pore mechanisms when lysing the cell membrane. This study highlights the significance of the effect of the overall charges and the hydrophobicity of the novel AMPs designed for antimicrobial activity.

Prospecting in silico antibacterial activity of a peptide from trypsin inhibitor isolated from tamarind seed

Bacterial infections have become a global concern, stimulating the growing demand for natural and biologically safe therapeutic agents with antibacterial action. This study was evaluated the genotoxicity of the trypsin inhibitor isolated from tamarind seeds (TTI) and the antibacterial effect of TTI theoric model, number 56, and conformation number 287 (TTIp 56/287) and derived peptides in silico. TTI (0.3 and 0.6 mg.mL⁻¹) did not cause genotoxicity in cells (p > 0.05). In silico, a greater interaction of TTIp 56/287 with the Gram-positive membrane (GP) was observed, with an interaction potential energy (IPE) of −1094.97 kcal.mol⁻¹. In the TTIp 56/287-GP interaction, the Arginine, Threonine (Thr), and Lysine residues presented lower IPE. In molecular dynamics (MD), Peptidotrychyme59 (TVSQTPIDIPIGLPVR) showed an IPE of −518.08 kcal.mol⁻¹ with the membrane of GP bacteria, and the Thr and Arginine residues showed the greater IPE. The results highlight new perspectives on TTI and its derived peptides antibacterial activity.

Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations

With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure-activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide-membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP-membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.

Antimicrobial Peptide Mechanism Studied by Scattering-Guided Molecular Dynamics Simulation

In an effort to combat rising antimicrobial resistance, our labs have rationally designed cationic, helical, amphipathic antimicrobial peptides (AMPs) as alternatives to traditional antibiotics since AMPs incur bacterial resistance in weeks, rather than days. One highly positively charged AMP, WLBU2 (+13e), (RRWV RRVR RWVR RVVR VVRR WVRR), has been shown to be effective in killing both Gram-negative (G(-)) and Gram-positive (G(+)) bacteria by directly perturbing the bacterial membrane nonspecifically. Previously, we used two equilibrium experimental methods: synchrotron X-ray diffuse scattering (XDS) providing lipid membrane thickness and neutron reflectometry (NR) providing WLBU2 depth of penetration into three lipid model membranes (LMMs). The purpose of the present study is to use the results from the scattering experiments to guide molecular dynamics (MD) simulations to investigate the detailed biophysics of the interactions of WLBU2 with LMMs of Gram-negative outer and inner membranes, and Gram-positive cell membranes, to elucidate the mechanisms of bacterial killing. Instead of coarse-graining, backmapping, or simulating without bias for several microseconds, all-atom (AA) simulations were guided by the experimental results and then equilibrated for ∼0.5 μs. Multiple replicas of the inserted peptide were run to probe stability and reach a combined time of at least 1.2 μs for G(-) and also 2.0 μs for G(+). The simulations with experimental comparisons help rule out certain structures and orientations and propose the most likely set of structures, orientations, and effects on the membrane. The simulations revealed that water, phosphates, and ions enter the hydrocarbon core when WLBU2 is positioned there. For an inserted peptide, the three types of amino acids, arginine, tryptophan, and valine (R, W, V), are arranged with the 13 Rs extending from the hydrocarbon core to the phosphate group, Ws are located at the interface, and Vs are more centrally located. For a surface state, R, W, and V are positioned relative to the bilayer interface as expected from their hydrophobicities, with Rs closest to the phosphate group, Ws close to the interface, and Vs in between. G(-) and G(+) LMMs are thinned ∼1 Å by the addition of WLBU2. Our results suggest a dual anchoring mechanism for WLBU2 both in the headgroup and in the hydrocarbon region that promotes a defect region where water and ions can flow across the slightly thinned bacterial cell membrane.

Atomic-Resolution Structures and Mode of Action of Clinically Relevant Antimicrobial Peptides

Global rise of infections and deaths caused by drug-resistant bacterial pathogens are among the unmet medical needs. In an age of drying pipeline of novel antibiotics to treat bacterial infections, antimicrobial peptides (AMPs) are proven to be valid therapeutics modalities. Direct in vivo applications of many AMPs could be challenging; however, works are demonstrating encouraging results for some of them. In this review article, we discussed 3-D structures of potent AMPs e.g., polymyxin, thanatin, MSI, protegrin, OMPTA in complex with bacterial targets and their mode of actions. Studies on human peptide LL37 and de novo-designed peptides are also discussed. We have focused on AMPs which are effective against drug-resistant Gram-negative bacteria. Since treatment options for the infections caused by super bugs of Gram-negative bacteria are now extremely limited. We also summarize some of the pertinent challenges in the field of clinical trials of AMPs.

Biological Membrane-Penetrating Peptides: Computational Prediction and Applications

Peptides comprise a versatile class of biomolecules that present a unique chemical space with diverse physicochemical and structural properties. Some classes of peptides are able to naturally cross the biological membranes, such as cell membrane and blood-brain barrier (BBB). Cell-penetrating peptides (CPPs) and blood-brain barrier-penetrating peptides (B3PPs) have been explored by the biotechnological and pharmaceutical industries to develop new therapeutic molecules and carrier systems. The computational prediction of peptides’ penetration into biological membranes has been emerged as an interesting strategy due to their high throughput and low-cost screening of large chemical libraries. Structure- and sequence-based information of peptides, as well as atomistic biophysical models, have been explored in computer-assisted discovery strategies to classify and identify new structures with pharmacokinetic properties related to the translocation through biomembranes. Computational strategies to predict the permeability into biomembranes include cheminformatic filters, molecular dynamics simulations, artificial intelligence algorithms, and statistical models, and the choice of the most adequate method depends on the purposes of the computational investigation. Here, we exhibit and discuss some principles and applications of these computational methods widely used to predict the permeability of peptides into biomembranes, exhibiting some of their pharmaceutical and biotechnological applications.

Probing the Functional Interaction Interface of Lipopolysaccharide and Antimicrobial Peptides: A Solution-State NMR Perspective

Antimicrobial peptides (AMPs) have been a topic of substantial research as the next-generation antibiotics. They have been extensively studied for the selectivity and action against microbial membrane lipids in imparting their targeted functioning. To determine the effectivity of the peptides against the Gram-negative pathogens, it is imperative to elucidate their role in interacting with the lipopolysaccharide moieties. Lipopolysaccharide is a major component of the outer membrane of the Gram-negative bacteria. It serves to protect the bacteria as well as govern the functionality of several antibacterial agents. It can prevent the access of the agents into the inner membrane of the bacteria, thus rendering them inactive. Several techniques have been employed to study the interaction for better designing of peptides; NMR spectroscopy is one of the most widely used techniques in determining the interactive properties of peptides with LPS as it can provide the details in atomistic level. NMR spectroscopy provides information about the structural and conformational changes as well as the dynamics of the interactions.

Spontaneous transmembrane pore formation by short-chain synthetic peptide

Amphiphilic β-peptides, which are synthetically designed short-chain helical foldamers of β-amino acids, are established potent biomimetic alternatives of natural antimicrobial peptides. An intriguing question is how the distinct molecular architecture of these short-chain and rigid synthetic peptides translates to its potent membrane-disruption ability. Here, we address this question via a combination of all-atom and coarse-grained molecular dynamics simulations of the interaction of mixed phospholipid bilayer with an antimicrobial 10-residue globally amphiphilic helical β-peptide at a wide range of concentrations. The simulation demonstrates that multiple copies of this synthetic peptide, initially placed in aqueous solution, readily self-assemble and adsorb at membrane interface. Subsequently, beyond a threshold peptide/lipid ratio, the surface-adsorbed oligomeric aggregate moves inside the membrane and spontaneously forms stable water-filled transmembrane pores via a cooperative mechanism. The defects induced by these pores lead to the dislocation of interfacial lipid headgroups, membrane thinning, and substantial water leakage inside the hydrophobic core of the membrane. A molecular analysis reveals that despite having a short architecture, these synthetic peptides, once inside the membrane, would stretch themselves toward the distal leaflet in favor of potential contact with polar headgroups and interfacial water layer. The pore formed in coarse-grained simulation was found to be resilient upon structural refinement. Interestingly, the pore-inducing ability was found to be elusive in a non-globally amphiphilic sequence isomer of the same β-peptide, indicating strong sequence dependence. Taken together, this work puts forward key perspectives of membrane activity of minimally designed synthetic biomimetic oligomers relative to the natural antimicrobial peptides.

Translocation of Endo-Functionalized Molecular Tubes across Different Lipid Bilayers: Atomistic Molecular Dynamics Simulation Study

Various artificial receptors, such as calixarenes, cyclodextrins, cucurbit[n]urils, and their acyclic compounds, pliiar[n]arenes, deep cavitands, and molecular tweezers, can permeate the lipid membranes and they are used as drug carriers to improve the drug solubility, stability, and bioavailability. Inspired by these, we have employed atomistic molecular dynamics simulation to examine the effects of endo-functionalized molecular tubes or naphthotubes (host-1a and host-1b) on seven different types of model lipid bilayers and the permeation properties of these receptors through these model lipid bilayers. Lipid types include six model lipid bilayers (POPC, POPE, DOPC, POPG, DPPE, POPE/POPG) and one realistic membrane (Yeast). We observe that these receptors are spontaneously translocated toward these model lipid bilayer head regions and do not proceed further into these lipid bilayer tail regions (reside at the interface between lipid head and lipid tail region), except for the DPPE-containing systems. In the DPPE model lipid bilayer-containing systems (1a-dppe and 1b-dppe), receptor molecules are only adsorbed on the bilayer surface and reside at the interface between lipid head and water. This finding is also supported by the biased free-energy profiles of these translocation processes. Passive transport of these receptors may be possible through these model lipid bilayers (due to low barrier height), except for DPPE bilayer-containing systems (that have a very high energy barrier at the center). The results from these simulations provide insight into the biocompatibility of host-1a or host-1b in microscopic detail. Based on this work, more research is needed to fully comprehend the role of these synthesized receptors as a prospective drug carrier.

Aggregation of Lysozyme in the Presence of a Mixed Bilayer of POPC and POPG

Understanding the molecular mechanisms by which amyloidogenic proteins interact with membranes is a challenging task. Amyloid accumulates from many human diseases have been observed to contain membrane lipids. In this work, coarse-grained molecular dynamics simulations have been used to inspect hen egg white lysozyme (HEWL) aggregation and membrane association in the presence of a pure POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) bilayer and a POPC and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol) mixed bilayer. It was observed that, in both cases, two HEWLs formed aggregates. In the presence of a mixed bilayer, after aggregation, the aggregated system started to interact with the membrane. It has been found that one of the lysozymes which came closer to the mixed bilayer unfolded more. The process of the initial insertion of an aggregated system in the mixed bilayer has been analyzed. The structural rearrangements of the protein and lipids were analyzed as well along the course of the simulation. Although with a pure POPC bilayer, aggregation was observed, the aggregated system moved away from the membrane. We believe that our study will provide considerable insights into lysozyme aggregation in the presence of a membrane environment.

Computational Methods and Tools in Antimicrobial Peptide Research

The evolution of antibiotic-resistant bacteria is an ongoing and troubling development that has increased the number of diseases and infections that risk going untreated. There is an urgent need to develop alternative strategies and treatments to address this issue. One class of molecules that is attracting significant interest is that of antimicrobial peptides (AMPs). Their design and development has been aided considerably by the applications of molecular models, and we review these here. These methods include the use of tools to explore the relationships between their structures, dynamics, and functions and the increasing application of machine learning and molecular dynamics simulations. This review compiles resources such as AMP databases, AMP-related web servers, and commonly used techniques, together aimed at aiding researchers in the area toward complementing experimental studies with computational approaches.

Quality of force fields and sampling methods in simulating pepX peptides: a case study for intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are a group of proteins that lack well-defined structures under native conditions and carry out crucial physiological functions in various biochemical pathways. Due to the heterogeneous nature of IDPs, molecular dynamics simulations have been extensively adopted to investigate the conformational ensembles and dynamic properties of these proteins. However, their accuracy remains limited by the development of force fields and sampling algorithms. Here, we evaluated the quality of both force fields and enhanced sampling algorithms based on five short pepX peptides. Our results show that the more extended conformational ensembles sampled by the AMOEBA polarizable force field present a higher ability to reproduce experimental NMR observables than AMBER and CHARMM classical force fields. Moreover, a better agreement with experiments is achieved in the simulation of IaMD (integrated accelerated molecular dynamics) than in aMD (accelerated molecular dynamics). The results together indicate that the combination of AMOEBA force field and IaMD enhanced sampling might be a better choice for simulating IDPs. This work may provide important clues for developments and applications of force fields and enhanced sampling methods in future simulations of IDPs.

Molecular Dynamics Simulation of the Interaction of Two Linear Battacin Analogs with Model Gram-Positive and Gram-Negative Bacterial Cell Membranes

Antimicrobial peptides (AMPs) are a potential solution to the increasing threat of antibiotic resistance, but successful design of active but nontoxic AMPs requires understanding their mechanism of action. Molecular dynamics (MD) simulations can provide atomic-level information regarding how AMPs interact with the cell membrane. Here, we have used MD simulations to study two linear analogs of battacin, a naturally occurring cyclic, lipidated, nonribosomal AMP. Like battacin, these analogs are active against Gram-negative multidrug resistant and Gram-positive bacteria, but they are less toxic than battacin. Our simulations show that this activity depends upon a combination of positively charged and hydrophobic moieties. Favorable interactions with negatively charged membrane lipid head groups drive association with the membrane and insertion of hydrophobic residues, and the N-terminal lipid anchors the peptides to the membrane surface. Both effects are required for stable membrane binding.

Aggregation State of Synergistic Antimicrobial Peptides

By integrating various simulation and experimental techniques, we discovered that antimicrobial peptides (AMPs) may achieve synergy at an optimal concentration and ratio, which can be caused by aggregation of the synergistic peptides. On multiple time and length scales, our studies obtain novel evidence of how peptide coaggregation in solution can affect the disruption of membranes by synergistic AMPs. Our findings provide crucial details about the complex molecular origins of AMP synergy, which will help guide the future development of synergistic AMPs as well as applications of anti-infective peptide cocktail therapies.

Lipid headgroups alter huntingtin aggregation on membranes

Huntington's Disease is a fatal neurodegenerative disorder caused by expansion of a glutamine repeat region (polyQ) beyond a critical threshold within exon1 of the huntingtin protein (htt). As a consequence of polyQ expansion, htt associates into a variety of aggregate species that are thought to underlie cellular toxicity. Within cells, htt associates with numerous membranous organelles and surfaces that exert influence on the aggregation process. In particular, the first 17 amino acids at the N-terminus of htt (Nt17) serve as a lipid-binding domain that is intrinsically disordered in bulk solution but adopts an amphipathic α-helical structure upon binding membranes. Beyond this, Nt17 is implicated in initiating htt fibrillization. As the interaction between Nt17 and lipid membranes is likely influenced by lipid properties, the impact of lipid headgroups on htt-exon1 aggregation, membrane activity, and the ability to form protein:lipid complexes was determined. Htt-exon1 with a disease-length polyQ domain (46Q) was exposed to lipid vesicles comprised of lipids with either zwitterionic (POPC and POPE) or anionic (POPG and POPS) headgroups. With zwitterionic head groups, large lipid to peptide ratios were required to have a statistically significant impact on htt aggregation. Anionic lipids enhanced htt fibrillization, even at low lipid:protein ratios, and this was accompanied by changes in aggregate morphology. Despite the larger impact of anionic lipids, htt-exon1(46Q) was more membrane active with zwitterionic lipid systems. The ability of Nt17 to form complexes with lipids was also mediated by lipid headgroups as zwitterionic ionic lipids more readily associated with multimeric forms of Nt17 in comparison with anionic lipids. Collectively, these results highlight the complexity of htt/membrane interactions and the resulting impact on the aggregation process.

Plantaricin NC8 αβ exerts potent antimicrobial activity against Staphylococcus spp. and enhances the effects of antibiotics

The use of conventional antibiotics has substantial clinical efficacy, however these vital antimicrobial agents are becoming less effective due to the dramatic increase in antibiotic-resistant bacteria. Novel approaches to combat bacterial infections are urgently needed and bacteriocins represent a promising alternative. In this study, the activities of the two-peptide bacteriocin PLNC8 αβ were investigated against different Staphylococcus spp. The peptide sequences of PLNC8 α and β were modified, either through truncation or replacement of all L-amino acids with D-amino acids. Both L- and D-PLNC8 αβ caused rapid disruption of lipid membrane integrity and were effective against both susceptible and antibiotic resistant strains. The D-enantiomer was stable against proteolytic degradation by trypsin compared to the L-enantiomer. Of the truncated peptides, β1–22, β7–34 and β1–20 retained an inhibitory activity. The peptides diffused rapidly (2 min) through the bacterial cell wall and permeabilized the cell membrane, causing swelling with a disorganized peptidoglycan layer. Interestingly, sub-MIC concentrations of PLNC8 αβ substantially enhanced the effects of different antibiotics in an additive or synergistic manner. This study shows that PLNC8 αβ is active against Staphylococcus spp. and may be developed as adjuvant in combination therapy to potentiate the effects of antibiotics and reduce their overall use.

Dimerization energetics of the G-protein coupled bile acid receptor TGR5 from all-atom simulations

We describe the first extensive energetic evaluation of GPCR dimerization on the atomistic level by means of potential of mean force (PMF) computations and implicit solvent/implicit membrane end-point free energy calculations (MM-PBSA approach). Free energies of association computed from the PMFs show that the formation of both the 1/8 and 4/5 interface is energetically favorable for TGR5, the first GPCR known to be activated by hydrophobic bile acids and neurosteroids. Furthermore, formation of the 1/8 interface is favored over that of the 4/5 interface. Both results are in line with our previous FRET experiments in live cells. Differences in lipid-protein interactions are identified to contribute to the observed differences in free energies of association. A per-residue decomposition of the MM-PBSA effective binding energy reveals hot spot residues specific for both interfaces that form clusters. This knowledge may be used to guide the design of dimerization inhibitors or perform mutational studies to explore physiological consequences of distorted TGR5 association. Finally, we characterized the role of Y111, located in the conserved (D/E)RY motif, as a facilitator of TGR5 interactions. The types of computations performed here should be transferable to other transmembrane proteins that form dimers or higher oligomers as long as good structural models of the dimeric or oligomeric states are available. Such computations may help to overcome current restrictions due to an imperfect energetic representation of protein association at the coarse-grained level. © 2019 Wiley Periodicals, Inc.

Liposomes with pH responsive 'on and off' switches for targeted and intracellular delivery of antibiotics

pH responsive drug delivery systems are one of the new strategies to address the spread of bacterial resistance to currently used antibiotics. The aim of this study was to formulate liposomes with 'On' and 'Off'' pH responsive switches for infection site targeting. The vancomycin (VCM) loaded liposomes had sizes below 100 nm, at pH 7.4. The QL-liposomes had a negative zeta potential at pH 7.4 that switched to a positive charge at acidic pH. VCM release from the liposome was quicker at pH 6 than pH 7.4. The OA-QL-liposome showed 4-fold lower MIC at pH 7.4 and 8- and 16-fold lower at pH 6.0 against both MSSA and MRSA compared to the bare drug. OA-QL liposome had a 1266.67- and 704.33-fold reduction in the intracellular infection for TPH-1 macrophage and HEK293 cells respectively. In vivo studies showed that the amount of MRSA recovered from mice treated with formulations was 189.67 and 6.33-fold lower than the untreated and bare VCM treated mice respectively. MD simulation of the QL lipid with the phosphatidylcholine membrane (POPC) showed spontaneous binding of the lipid to the bilayer membrane both electrostatic and Van der Waals interactions contributed to the binding. These studies demonstrated that the 'On' and 'Off' pH responsive liposomes enhanced the activity targeted and intracellular delivery VCM.

Accelerated Molecular Dynamics Simulation for Helical Proteins Folding in Explicit Water

In this study, we examined the folding processes of eight helical proteins (2I9M, TC5B, 1WN8, 1V4Z, 1HO2, 1HLL, 2KFE, and 1YYB) at room temperature using the explicit solvent model under the AMBER14SB force field with the accelerated molecular dynamics (AMD) and traditional molecular dynamics (MD), respectively. We analyzed and compared the simulation results obtained by these two methods based on several aspects, such as root mean square deviation (RMSD), native contacts, cluster analysis, folding snapshots, free energy landscape, and the evolution of the radius of gyration, which showed that these eight proteins were successfully and consistently folded into the corresponding native structures by AMD simulations carried out at room temperature. In addition, the folding occurred in the range of 40~180 ns after starting from the linear structures of the eight proteins at 300 K. By contrast, these stable folding structures were not found when the traditional molecular dynamics (MD) simulation was used. At the same time, the influence of high temperatures (350, 400, and 450 K) is also further investigated. Study found that the simulation efficiency of AMD is higher than that of MD simulations, regardless of the temperature. Of these temperatures, 300 K is the most suitable temperature for protein folding for all systems. To further investigate the efficiency of AMD, another trajectory was simulated for eight proteins with the same linear structure but different random seeds at 300 K. Both AMD trajectories reached the correct folded structures. Our result clearly shows that AMD simulation are a highly efficient and reliable method for the study of protein folding.

Enhanced Silkworm Cecropin B Antimicrobial Activity against Pseudomonas aeruginosa from Single Amino Acid Variation

Pseudomonas aeruginosa is an opportunistic bacterial pathogen causing severe infections in hospitalized and immunosuppressed patients, particularly individuals affected by cystic fibrosis. Several clinically isolated P. aeruginosa strains were found to be resistant to three or more antimicrobial classes indicating the importance of identifying new antimicrobials active against this pathogen. Here, we characterized the antimicrobial activity and the action mechanisms against P. aeruginosa of two natural isoforms of the antimicrobial peptide cecropin B, both isolated from the silkworm Bombyx mori. These cecropin B isoforms differ in a single amino acid substitution within the active portion of the peptide, so that the glutamic acid of the E53 CecB variant is replaced by a glutamine in the Q53 CecB isoform. Both peptides showed a high antimicrobial and membranolytic activity against P. aeruginosa, with Q53 CecB displaying greater activity compared with the E53 CecB isoform. Biophysical analyses, live-cell NMR, and molecular-dynamic-simulation studies indicated that both peptides might act as membrane-interacting elements, which can disrupt outer-membrane organization, facilitating their translocation toward the inner membrane of the bacterial cell. Our data also suggest that the amino acid variation of the Q53 CecB isoform represents a critical factor in stabilizing the hydrophobic segment that interacts with the bacterial membrane, determining the highest antimicrobial activity of the whole peptide. Its high stability to pH and temperature variations, tolerance to high salt concentrations, and low toxicity against human cells make Q53 CecB a promising candidate in the development of CecB-derived compounds against P. aeruginosa.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Potent and Broad-Spectrum Antimicrobial Activity of Analogs from the Scorpion Peptide Stigmurin

Scorpion venom constitutes a rich source of biologically active compounds with high potential for therapeutic and biotechnological applications that can be used as prototypes for the design of new drugs. The aim of this study was to characterize the structural conformation, evaluate the antimicrobial activity, and gain insight into the possible action mechanism underlying it, for two new analog peptides of the scorpion peptide Stigmurin, named StigA25 and StigA31. The amino acid substitutions in the native sequence for lysine residues resulted in peptides with higher positive net charge and hydrophobicity, with an increase in the theoretical helical content. StigA25 and StigA31 showed the capacity to modify their structural conformation according to the environment, and were stable to pH and temperature variation—results similar to the native peptide. Both analog peptides demonstrated broad-spectrum antimicrobial activity in vitro, showing an effect superior to that of the native peptide, being non-hemolytic at the biologically active concentrations. Therefore, this study demonstrates the therapeutic potential of the analog peptides from Stigmurin and the promising approach of rational drug design based on scorpion venom peptide to obtain new anti-infective agents.

Molecular Mechanism of Uptake of Cationic Photoantimicrobial Phthalocyanine across Bacterial Membranes Revealed by Molecular Dynamics Simulations

Phthalocyanines are aromatic macrocyclic compounds, which are structurally related to porphyrins. In clinical practice, phthalocyanines are used in fluorescence imaging and photodynamic therapy of cancer and noncancer lesions. Certain forms of the substituted polycationic metallophthalocyanines have been previously shown to be active in photodynamic inactivation of both Gram-negative and Gram-positive bacteria; one of them is zinc octakis(cholinyl)phthalocyanine (ZnPcChol⁸⁺). However, the molecular details of how these compounds translocate across bacterial membranes still remain unclear. In the present work, we have developed a coarse-grained (CG) molecular model of ZnPcChol⁸⁺ within the framework of the popular MARTINI CG force field. The obtained model was used to probe the solvation behavior of phthalocyanine molecules, which agreed with experimental results. Subsequently, it was used to investigate the molecular details of interactions between phthalocyanines and membranes of various compositions. The results demonstrate that ZnPcChol⁸⁺ has high affinity to both the inner and the outer model membranes of Gram-negative bacteria, although this species does not show noticeable affinity to the 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphatidylcholine membrane. Furthermore, we found out that the process of ZnPcChol⁸⁺ penetration toward the center of the outer bacterial membrane is energetically favorable and leads to its overall disturbance and formation of the aqueous pore. Such intramembrane localization of ZnPcChol⁸⁺ suggests their twofold cytotoxic effect on bacterial cells: (1) via induction of lipid peroxidation by enhanced production of reactive oxygen species (i.e., photodynamic toxicity); (2) via rendering the bacterial membrane more permeable for additional Pc molecules as well as other compounds. We also found that the kinetics of penetration depends on the presence of phospholipid defects in the lipopolysaccharide leaflet of the outer membrane and the type of counterions, which stabilize it. Thus, the results of our simulations provide a detailed molecular view of ZnPcChol⁸⁺ "self-promoted uptake", the pathway previously proposed for some small molecules crossing the outer bacterial membrane.

Molecular Insights into the Membrane Affinities of Model Hydrophobes

Membrane-active antibiotics are of great interest in fighting bacterial resistance. α-Mangostin is a membrane-active molecule, but there are no details of its mechanism of action at the atomistic level. We have employed free-energy simulations and microsecond-long conventional molecular dynamics simulations to study the mode of interaction of α-mangostin with a model bacterial membrane and compare it with the mechanisms of three hydrophobic molecules (ciprofloxacin, xanthone, and tetracycline). We find that α-mangostin is thermodynamically more favored to insert into the membrane compared to the other three molecules. Apart from tetracycline, which is largely hydrophilic, the other three molecules aggregate in water; however, only α-mangostin can penetrate into the lipid tail region of the membrane. When it reaches a high concentration in the lipid tail region, α-mangostin can form tubular clusters that span the two head group regions of the membrane, resulting in a large number of water translocations along the transmembrane aggregates. Structure-activity relationship analysis revealed two structural properties that characterize α-mangostin, namely, the two isoprenyl groups and the polar groups present in the aromatic rings, which result in "disruptive amphiphilicity" and hence its excellent membrane activity.

Insight into the interactions, residue snorkeling, and membrane disordering potency of a single antimicrobial peptide into different lipid bilayers

Pardaxin, with a bend-helix-bend-helix structure, is a membrane-active antimicrobial peptide that its membrane activity depends on the lipid bilayer composition. Herein, all-atom molecular dynamics (MD) simulations were performed to provide further molecular insight into the interactions, structural dynamics, orientation behavior, and cationic residues snorkeling of pardaxin in the DMPC, DPPC, POPC, POPG, POPG/POPE (3:1), and POPG/POPE (1:3) lipid bilayers. The results showed that the C-terminal helix of the peptide was maintained in all six types of the model-bilayers and pardaxin was tilted into the DMPC, DPPC, and POPG/POPE mixed bilayers more than the POPC and POPG bilayers. As well as, the structure of zwitterionic membranes was more affected by the peptide than the anionic bilayers. Taken together, the study demonstrated that the cationic residues of pardaxin snorkeled toward the interface of lipid bilayers and all phenylalanine residues of the peptide played important roles in the peptide-membrane interactions. We hope that this work will provide a better understanding of the interactions of antimicrobial peptides with the membranes.