Small molecules enhance the potency of natural antimicrobial peptides.

The skin-associated microbiome plays an important role in general well-being and in a variety of treatable skin conditions. In this regard, endogenous antimicrobial peptides have both a direct and indirect role in determining the composition of the microbiota. We demonstrate here that certain small molecular species can amplify the antimicrobial potency of naturally occurring antimicrobial peptides. In this study, we have used niacinamide, a form of vitamin B3 naturally found in foods and widely used in cosmetic skincare products, and two of its structural analogs, to investigate their cooperativity with the human antimicrobial peptide LL37 on the bacterium Staphylococcus aureus. We observed a clear synergistic effect of niacinamide and, to some extent, N-methylnicotinamide, whereas isonicotinamide showed no significant cooperativity with LL37. Adaptively biased molecular dynamics simulations using simplified model membrane substrates and single peptides revealed that these molecules partition into the headgroup region of an anionic bilayer used to mimic the bacterial membrane. The simulated effects on the physical properties of the simulated model membrane are well correlated with experimental activity observed in real biological assays despite the simplicity of the model. In contrast, these molecules have little effect on zwitterionic bilayers that mimic a mammalian membrane. We conclude that niacinamide and N-methylnicotinamide can therefore potentiate the activity of host peptides by modulating the physical properties of the bacterial membrane, and to a lesser extent through direct interactions with the peptide. The level of cooperativity is strongly dependent on the detailed chemistry of the additive, suggesting an opportunity to fine-tune the behavior of host peptides.

Switching Cytolytic Nanopores into Antimicrobial Fractal Ruptures by a Single Side Chain Mutation.

Disruption of cell membranes is a fundamental host defense response found in virtually all forms of life. The molecular mechanisms vary but generally lead to energetically favored circular nanopores. Here, we report an elaborate fractal rupture pattern induced by a single side-chain mutation in ultrashort (8-11-mers) helical peptides, which otherwise form transmembrane pores. In contrast to known mechanisms, this mode of membrane disruption is restricted to the upper leaflet of the bilayer where it exhibits propagating fronts of peptide-lipid interfaces that are strikingly similar to viscous instabilities in fluid flow. The two distinct disruption modes, pores and fractal patterns, are both strongly antimicrobial, but only the fractal rupture is nonhemolytic. The results offer wide implications for elucidating differential membrane targeting phenomena defined at the nanoscale.

Modulation of Antimicrobial Peptide Potency in Stressed Lipid Bilayers.

It is shown that the tendency of an archetypal antimicrobial peptide to insert into and perforate a simple lipid bilayer is strongly modulated by tensile stress in the membrane. The results, obtained through molecular dynamics simulations, have been demonstrated with several lipid compositions and appear to be general, although quantitative details differ. The findings imply that the potency of antimicrobial peptides may not be a purely intrinsic chemical property and, instead, depends on the mechanical state of the target membrane.

Author Correction: Tuneable poration: host defense peptides as sequence probes for antimicrobial mechanisms.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Tuneable poration: host defense peptides as sequence probes for antimicrobial mechanisms.

The spread of antimicrobial resistance stimulates discovery strategies that place emphasis on mechanisms circumventing the drawbacks of traditional antibiotics and on agents that hit multiple targets. Host defense peptides (HDPs) are promising candidates in this regard. Here we demonstrate that a given HDP sequence intrinsically encodes for tuneable mechanisms of membrane disruption. Using an archetypal HDP (cecropin B) we show that subtle structural alterations convert antimicrobial mechanisms from native carpet-like scenarios to poration and non-porating membrane exfoliation. Such distinct mechanisms, studied using low- and high-resolution spectroscopy, nanoscale imaging and molecular dynamics simulations, all maintain strong antimicrobial effects, albeit with diminished activity against pathogens resistant to HDPs. The strategy offers an effective search paradigm for the sequence probing of discrete antimicrobial mechanisms within a single HDP.

Cooperative Modes of Action of Antimicrobial Peptides Characterized with Atomistic Simulations: A Study on Cecropin B.

Antimicrobial peptides (AMPs) are widely occurring host defense agents of interest as one route for addressing the growing problem of multidrug-resistant pathogens. Understanding the mechanisms behind their antipathogen activity is instrumental in designing new AMPs. Herein, we present an all-atom molecular dynamics and free energy study on cecropin B (CB) and its constituent domains. We find a cooperative mechanism in which CB inserts into an anionic model membrane with its amphipathic N-terminal segment, supported by the hydrophobic C-terminal segment of a second peptide. The two peptides interact via a Glu···Lys salt bridge and together sustain a pore in the membrane. Using a modified membrane composition, we demonstrate that when the lower leaflet is overall neutral, insertion of the cationic segment is retarded and thus this mode of action is membrane specific. The observed mode of action utilizes a flexible hinge, a common structural motif among AMPs, which allows CB to insert into the membrane using either or both termini. Data from both unbiased trajectories and enhanced sampling simulations indicate that a requirement for CB to be an effective AMP is the interaction of its hydrophobic C-terminal segment with the membrane. Simulations of these segments in isolation reveal their aggregation in the membrane and a different mechanism of supporting pore formation. Together, our results show the complex interaction of different structural motifs of AMPs and, in particular, a potential role for electronegative side chains in an overall cationic AMP.

Binding of Zn(II) to Tropomyosin Receptor Kinase A in Complex with Its Cognate Nerve Growth Factor: Insights from Molecular Simulation and in Vitro Essays.

The binding of the human nerve growth factor (NGF) protein to tropomyosin receptor kinase A (TrkA) is associated with Alzhemeir's development. Owing to the large presence of zinc(II) ions in the synaptic compartments, the zinc ions might be bound to the complex in vivo. Here, we have identified a putative zinc binding site using a combination of computations and experiments. First, we have predicted structural features of the NGF/TrkA complex in an aqueous solution by molecular simulation. Metadynamics free energy calculations suggest that these are very similar to those in the X-ray structure. Here, the "crab" structure of the NGF shape binds tightly to two TrkA "pincers". Transient conformations of the complex include both more extended and more closed conformations. Interestingly, the latter features facial histidines (His60 and His61) among the N-terminal D1-D3 domains, each of which is a potential binding region for biometals. This suggests the presence of a four-His Zn binding site connecting the two chains. To address this issue, we investigated the binding of a D1-D3 domains' peptide mimic by stability constant and nuclear magnetic resonance measurements, complemented by density functional theory-based calculations. Taken together, these establish unambiguously a four-His coordination of the metal ion in the model systems, supporting the presence of our postulated binding site in the NGF/TrkA complex.

Designing the Sniper: Improving Targeted Human Cytolytic Fusion Proteins for Anti-Cancer Therapy via Molecular Simulation.

Targeted human cytolytic fusion proteins (hCFPs) are humanized immunotoxins for selective treatment of different diseases including cancer. They are composed of a ligand specifically binding to target cells genetically linked to a human apoptosis-inducing enzyme. hCFPs target cancer cells via an antibody or derivative (scFv) specifically binding to e.g., tumor associated antigens (TAAs). After internalization and translocation of the enzyme from endocytosed endosomes, the human enzymes introduced into the cytosol are efficiently inducing apoptosis. Under in vivo conditions such enzymes are subject to tight regulation by native inhibitors in order to prevent inappropriate induction of cell death in healthy cells. Tumor cells are known to upregulate these inhibitors as a survival mechanism resulting in escape of malignant cells from elimination by immune effector cells. Cytosolic inhibitors of Granzyme B and Angiogenin (Serpin P9 and RNH1, respectively), reduce the efficacy of hCFPs with these enzymes as effector domains, requiring detrimentally high doses in order to saturate inhibitor binding and rescue cytolytic activity. Variants of Granzyme B and Angiogenin might feature reduced affinity for their respective inhibitors, while retaining or even enhancing their catalytic activity. A powerful tool to design hCFPs mutants with improved potency is given by in silico methods. These include molecular dynamics (MD) simulations and enhanced sampling methods (ESM). MD and ESM allow predicting the enzyme-protein inhibitor binding stability and the associated conformational changes, provided that structural information is available. Such "high-resolution" detailed description enables the elucidation of interaction domains and the identification of sites where particular point mutations may modify those interactions. This review discusses recent advances in the use of MD and ESM for hCFP development from the viewpoints of scientists involved in both fields.

EGFR oligomerization organizes kinase-active dimers into competent signalling platforms.

Epidermal growth factor receptor (EGFR) signalling is activated by ligand-induced receptor dimerization. Notably, ligand binding also induces EGFR oligomerization, but the structures and functions of the oligomers are poorly understood. Here, we use fluorophore localization imaging with photobleaching to probe the structure of EGFR oligomers. We find that at physiological epidermal growth factor (EGF) concentrations, EGFR assembles into oligomers, as indicated by pairwise distances of receptor-bound fluorophore-conjugated EGF ligands. The pairwise ligand distances correspond well with the predictions of our structural model of the oligomers constructed from molecular dynamics simulations. The model suggests that oligomerization is mediated extracellularly by unoccupied ligand-binding sites and that oligomerization organizes kinase-active dimers in ways optimal for auto-phosphorylation in trans between neighbouring dimers. We argue that ligand-induced oligomerization is essential to the regulation of EGFR signalling.

Design of human granzyme B variants resistant to serpin B9.

Human granzyme B (hGB) is a serine protease involved in immune-mediated apoptosis. Its cytotoxicity makes it potentially applicable in cancer therapy. However, the effectiveness of hGB can be hampered by the cytosolic expression of a natural protein inhibitor, human Serpin B9 (hSB9). Here, we used computational approaches to identify hGB mutations that can affect its binding to hSB9 without significantly decreasing its catalytic efficiency. Alanine-scanning calculations allowed us to identify residues of hGB important for the interaction with hSB9. Some variants were selected, and molecular dynamic simulations on the mutated hGB in complex with hSB9 in aqueous solution were carried out to investigate the effect of these variants on the stability of the complex. The R28K, R201A, and R201K mutants significantly destabilized the interaction of the protein with hSB9. Consistently, all of these variants also retained their activity in the presence of the Serpin B9 inhibitor in subsequent in vitro assays of wild-type and mutated hGB. In particular, the activity of R201K hGB with and without Serpin B9 is very similar to that of the wild-type protein. Hence, R201K hGB emerges as a promising species for antitumoral therapy applications.

Structural role of compensatory amino acid replacements in the α-synuclein protein.

A subset of familial Parkinson's disease (PD) cases is associated with the presence of disease-causing point mutations in human α-synuclein [huAS(wt)], including A53T. Surprisingly, the human neurotoxic amino acid 53T is present in non-primate, wild-type sequences of α-synucleins, including that expressed by mice [mAS(wt)]. Because huAS(A53T) causes neurodegeneration when expressed in rodents, the amino acid changes between the wild-type human protein [huAS(wt)] and mAS(wt) might act as intramolecular suppressors of A53T toxicity in the mouse protein, restoring its physiological structure and function. The lack of structural information for mAS(wt) in aqueous solution has prompted us to conduct a comparative molecular dynamics study of huAS(wt), huAS(A53T), and mAS(wt) in water at 300 K. The calculations are based on an ensemble of nuclear magnetic resonance-derived huAS(wt) structures. huAS(A53T) turns out to be more flexible and less compact than huAS(wt). Its central (NAC) region, involved in fibril formation by the protein, is more solvent-exposed than that of the wild-type protein, in agreement with nuclear magnetic resonance data. The compactness of mAS(wt) is similar to that of the human protein. In addition, its NAC region is less solvent-exposed and more rigid than that of huAS(A53T). All of these features may be caused by an increase in the level of intramolecular interactions on passing from huAS(A53T) to mAS(wt). We conclude that the presence of "compensatory replacements" in the mouse protein causes a significant change in the protein relative to huAS(A53T), restoring features not too dissimilar to those of the human protein.

Modulation of alpha-synuclein aggregation by dopamine analogs.

The action of dopamine on the aggregation of the unstructured alpha-synuclein (alpha-syn) protein may be linked to the pathogenesis of Parkinson's disease. Dopamine and its oxidation derivatives may inhibit alpha-syn aggregation by non-covalent binding. Exploiting this fact, we applied an integrated computational and experimental approach to find alternative ligands that might modulate the fibrillization of alpha-syn. Ligands structurally and electrostatically similar to dopamine were screened from an established library. Five analogs were selected for in vitro experimentation from the similarity ranked list of analogs. Molecular dynamics simulations showed they were, like dopamine, binding non-covalently to alpha-syn and, although much weaker than dopamine, they shared some of its binding properties. In vitro fibrillization assays were performed on these five dopamine analogs. Consistent with our predictions, analyses by atomic force and transmission electron microscopy revealed that all of the selected ligands affected the aggregation process, albeit to a varying and lesser extent than dopamine, used as the control ligand. The in silico/in vitro approach presented here emerges as a possible strategy for identifying ligands interfering with such a complex process as the fibrillization of an unstructured protein.

Computational Study of the DNA-Binding Protein Helicobacter pylori NikR: The Role of Ni(2.).

An integrated approach, combining atomistic molecular dynamics simulations, coarse-grained models, and solution NMR, was used to characterize the internal dynamics of HpNikR, a Ni-dependent transcription factor. Specifically, these methods were used to ascertain how the presence of bound Ni(2+) ions affects the stability of the known open, cis, and trans forms observed in the crystal structures of this protein as well as their interconversion capability. The consensus picture emerging from all the collected data hints at the interconversion of NikR among the three types of conformations, regardless of the content of bound Ni(2+). On the basis of atomistic and coarse-grained simulations, we deduce that the interconversion capability is particularly effective between the cis and the open forms and appreciably less so between the trans conformer and the other two forms. The presence of the bound Ni(2+) ions does, however, affect significantly the degree of the correlations on the two DNA-binding domains of NikR, which is significantly suppressed as compared to the apo form. Overall, the findings suggest that the binding of HpNikR to DNA occurs through a sophisticated multistep process involving both a conformational selection and an induced fit.