Interfacial tryptophan residues: a role for the cation-pi effect?

Molecular dynamics simulations support a preference of cyclotide kalata B1 for phosphatidylethanolamine phospholipids

Kalata B1 (kB1), a naturally occurring cyclotide has been shown experimentally to bind lipid membranes that contain phosphatidylethanolamine (PE) phospholipids. Here, molecular dynamics simulations were used to explore its interaction with two phospholipids, palmitoyloleoylphosphatidylethanolamine (POPE), palmitoyloleoylphosphatidylcholine (POPC), and a heterogeneous membrane comprising POPC/POPE (90:10), to understand the basis for the selectivity of kB1 towards PE phospholipids. The simulations showed that in the presence of only 10 % POPE lipid, kB1 forms a stable binding complex with membrane bilayers. An ionic interaction between the E7 carboxylate group of kB1 and the ammonium group of PE headgroups consistently initiates binding of kB1 to the membrane. Additionally, stable noncovalent interactions such as hydrogen bonding (E7, T8, V10, G11, T13 and N15), cation-π (W23), and CH-π (W23) interactions between specific residues of kB1 and the lipid membrane play an important role in stabilizing the binding. These findings are consistent with a structure-activity relationship study on kB1 where lysine mutagenesis on the bioactive and hydrophobic faces of the peptide abolished membrane-dependent bioactivities. In summary, our simulations suggest the importance of residue E7 (in the bioactive face) in enabling kB1 to recognize and bind selectively to PE-containing phospholipids bilayers through ionic and hydrogen bonding interactions, and of W23 (in the hydrophobic face) for the association and insertion of kB1 into the lipid bilayer through cation-π and CH-π interactions. Overall, this work enhances our understanding of the molecular basis of the membrane binding and bioactivity of this prototypic cyclotide.

Tryptophan- and arginine-rich antimicrobial peptides: Anti-infectives with great potential

With the increasing prevalence of multidrug resistant (MDR) bacteria, there is a need to design synthetic antimicrobial peptides (AMPs) that are effective and selective for bacteria, i.e. non-toxic to mammalian cells. One design strategy, namely the use of tryptophan- and arginine-rich AMPs, is rooted in the study of natural AMPs that are composed mainly of these amino acids, such as lactoferricin, tritrpticin, and puroindoline. A number of important studies on these AMPs by the Vogel group are reviewed here. More recent work on W-/R-rich peptides is also presented. The examples show that these peptides represent anti-infectives with great potential.

Sensitive Detection of Biomolecular Adsorption by a Low-Density Surfactant Layer Using Sum-Frequency Vibrational Spectroscopy

Small molecules or proteins interact with a biomembrane in various ways for molecular recognition, structure stabilization, and transmembrane signaling. In this study, 1,2-dipalmitoyl-3-trimethylammonium-propane (DPTAP), having a choline group, was used to investigate this interaction by using sum-frequency vibrational spectroscopy. The sum-frequency spectrum characteristic of a neat monolayer changed to that of a bare air/water interface at a larger molecular area of the DPTAP molecules due to local laser heating. Upon introduction of the aromatic molecules in the subphase at around 120 Å2 per molecule, the sum-frequency signal suddenly reappeared due to molecular adhesion, and this was utilized to probe the adsorption of the aromatic ring molecules in the water subphase to the choline headgroup of the DPTAP by cation-π interaction. The onset concentrations of this sum-frequency signal change allowed a comparison of the relative interaction strengths between different aromatic molecules. A zwitterionic surfactant molecule (DPPC) was found to interact weakly compared to the cationic DPTAP molecule.

Structure-activity relationships of mitochondria-targeted tetrapeptide pharmacological compounds

Mitochondria play a central role in metabolic homeostasis, and dysfunction of this organelle underpins the etiology of many heritable and aging-related diseases. Tetrapeptides with alternating cationic and aromatic residues such as SS-31 (elamipretide) show promise as therapeutic compounds for mitochondrial disorders. In this study, we conducted a quantitative structure-activity analysis of three alternative tetrapeptide analogs, benchmarked against SS-31, that differ with respect to aromatic side chain composition and sequence register. We present the first structural models for this class of compounds, obtained with Nuclear Magnetic Resonance (NMR) and molecular dynamics approaches, showing that all analogs except for SS-31 form compact reverse turn conformations in the membrane-bound state. All peptide analogs bound cardiolipin-containing membranes, yet they had significant differences in equilibrium binding behavior and membrane interactions. Notably, analogs had markedly different effects on membrane surface charge, supporting a mechanism in which modulation of membrane electrostatics is a key feature of their mechanism of action. The peptides had no strict requirement for side chain composition or sequence register to permeate cells and target mitochondria in mammalian cell culture assays. All four peptides were pharmacologically active in serum withdrawal cell stress models yet showed significant differences in their abilities to restore mitochondrial membrane potential, preserve ATP content, and promote cell survival. Within our peptide set, the analog containing tryptophan side chains, SPN10, had the strongest impact on most membrane properties and showed greatest efficacy in cell culture studies. Taken together, these results show that side chain composition and register influence the activity of these mitochondria-targeted peptides, helping provide a framework for the rational design of next-generation therapeutics with enhanced potency.

Tuning Electrostatic and Hydrophobic Surfaces of Aromatic Rings to Enhance Membrane Association and Cell Uptake of Peptides

Aromatic groups are key mediators of protein–membrane association at cell surfaces, contributing to hydrophobic effects and π‐membrane interactions. Here we show electrostatic and hydrophobic influences of aromatic ring substituents on membrane affinity and cell uptake of helical, cyclic and cell penetrating peptides. Hydrophobicity is important, but subtle changes in electrostatic surface potential, dipoles and polarizability also enhance association with phospholipid membranes and cell uptake. A combination of fluorine and sulfur substituents on an aromatic ring induces microdipoles that enhance cell uptake of 12‐residue peptide inhibitors of p53‐HDM2 interaction and of cell‐penetrating cyclic peptides. These aromatic motifs can be readily inserted into peptide sidechains to enhance their cell uptake.

The Influence of Short Motifs on the Anticancer Activity of HB43 Peptide

Despite the remarkable similarity in amino acid composition, many anticancer peptides (ACPs) display significant differences in terms of activity. This strongly suggests that particular relative dispositions of amino acids (motifs) play a role in the interaction with their biological target, which is often the cell membrane. To better verify this hypothesis, we intentionally modify HB43, an ACP active against a wide variety of cancers. Sequence alignment of related ACPs by ADAPTABLE web server highlighted the conserved motifs that could be at the origin of the activity. In this study, we show that changing the order of amino acids in such motifs results in a significant loss of activity against colon and breast cancer cell lines. On the contrary, amino acid substitution in key motifs may reinforce or weaken the activity, even when the alteration does not perturb the amphipathicity of the helix formed by HB43 on liposomes mimicking their surface. NMR and MD simulations with different membrane models (micelles, bicelles, and vesicles) indicate that the activity reflects the insertion capability in cancer-mimicking serine-exposing membranes, supported by the insertion of N-terminal phenylalanine in the FAK motif and the anchoring to the carboxylate of phosphatidylserine by means of arginine side chains.

Improving the Activity of Antimicrobial Peptides Against Aquatic Pathogen Bacteria by Amino Acid Substitutions and Changing the Ratio of Hydrophobic Residues

With the increasing number of drug-resistant bacteria, there is an urgent need for new antimicrobial agents, and antimicrobial peptides (AMPs), which exist in the human non-specific immune system, are one of the most promising candidates. It is an effective optimization strategy to modify antimicrobial peptides (AMPs) according to the distribution of amino acids and hydrophobic characteristics. The addition of bacterial pheromones to the N short peptide can increase the ability to recognize bacteria. In this study, we designed and synthesized AMP1–6 by amino acid substitution of mBjAMP1. Additionally, P-6, S-6, and L-6 were designed and synthesized by adding bacterial pheromones based on 1–6. Functional tests showed that the four AMPs had the ability to kill Gram-negative Vibrio anguillarum, Pseudomonas mendocina, and Vibrio parahaemolyticus, and Gram-positive Micrococcus luteus and Listeria monocytogenes. Additionally, all four AMPs induced permeabilization and depolarization of bacterial cell membranes and increased intracellular reactive oxygen species (ROS) levels. Importantly, they had little or no mammalian cytotoxicity. At the same time, 1–6 and L-6 protected the stability of intestinal flora in Sebastes schlegelii and increased the relative abundance of Lactobacillaceae. In summary, our results indicate that the designed AMPs have broad application prospects as a new type of polypeptide antimicrobial agent.

Unfolding of Tryptophanoctyl Ester and Elastic Deformation of Host Micelles via RR'3 N+ ⋅⋅⋅π Interaction: Conceivable Relevance to Wrapping Process of Receptor Mediated Endocytosis

The interfacial properties of the mixed amphiphiles are modified by a stronger cation-π interaction between the quaternary ammonium head group of CTAB and the π-face of TROE, compared to the tyrosine analogue (TYOE). This eventually triggers a morphology transition through elastic deformation of the spherical micelles of CTAB to cylindrical/wormlike micelles. The unfolding of TROE and the molecular interactions in the nanoenvironment have been recognized by NMR spectroscopy and the physical characteristics of the entangled wormlike micelles are investigated by high resolution transmission electron microscopy (HRTEM), whereas the complex fluidic feature is examined by dynamic rheological measurements. Morphology tuning of the soft nanoaggregates of zwitterionic dodecylphosphocholine by the tryptophan analogue via choline-π interaction has unique biological consequences and we consider the significance of such interactions in facilitating endocytosis of a virion/nano particle(NP) in terms of a quantitative model. The implication in future research on drug development strategies is discussed.

Tryptophan, an Amino-Acid Endowed with Unique Properties and Its Many Roles in Membrane Proteins

Tryptophan is an aromatic amino acid with unique physico-chemical properties. It is often encountered in membrane proteins, especially at the level of the water/bilayer interface. It plays a role in membrane protein stabilization, anchoring and orientation in lipid bilayers. It has a hydrophobic character but can also engage in many types of interactions, such as π–cation or hydrogen bonds. In this review, we give an overview of the role of tryptophan in membrane proteins and a more detailed description of the underlying noncovalent interactions it can engage in with membrane partners.

Differential interactions of the antimicrobial peptide, RQ18, with phospholipids and cholesterol modulate its selectivity for microorganism membranes

BACKGROUND

Antimicrobial peptides (AMPs) are molecules with potential application for the treatment of microorganism infections. We, herein, describe the structure, activity, and mechanism of action of RQ18, an α-helical AMP that displays antimicrobial activity against Gram-positive and Gram-negative bacteria, and yeasts from the Candida genus.

METHODS

A physicochemical-guided design assisted by computer tools was used to obtain our lead peptide candidate, named RQ18. This peptide was assayed against Gram-positive and Gram-negative bacteria, yeasts, and mammalian cells to determine its selectivity index. The secondary structure and the mechanism of action of RQ18 were investigated using circular dichroism, large unilamellar vesicles, and molecular dynamic simulations.

RESULTS

RQ18 was not cytotoxic to human lung fibroblasts, peripheral blood mononuclear cells, red blood cells, or Vero cells at MIC values, exhibiting a high selectivity index. Circular dichroism analysis and molecular dynamic simulations revealed that RQ18 presents varying structural profiles in aqueous solution, TFE/water mixtures, SDS micelles, and lipid bilayers. The peptide was virtually unable to release carboxyfluorescein from large unilamellar vesicles composed of POPC/cholesterol, model that mimics the eukaryotic membrane, indicating that vesicles' net charges and the presence of cholesterol may be related with RQ18 selectivity for bacterial and fungal cell surfaces.

CONCLUSIONS

RQ18 was characterized as a membrane-active peptide with dual antibacterial and antifungal activities, without compromising mammalian cells viability, thus reinforcing its therapeutic application.

GENERAL SIGNIFICANCE

These results provide further insight into the complex process of AMPs interaction with biological membranes, in special with systems that mimic prokaryotic and eukaryotic cell surfaces.

Designing a leucine-rich antibacterial nonapeptide with potent activity against mupirocin-resistant MRSA via a structure-activity relationship study

Staphylococcus aureus is the main aetiological agent responsible for the majority of human skin infections. Of particular concern is the methicillin-resistant variety, commonly known as MRSA. The extensive use of the first-line topical antibiotic of choice, mupirocin, has inevitably resulted in the emergence of resistant strains, signalling an urgent need for the development of new antibacterials with new mechanisms of action. In this work, we describe how we designed a novel cationic nonapeptide, containing only leucine and two lysine residues, with potent anti-MRSA activity and a rapid bactericidal mode of action. Coupled to a favourable safety profile towards human skin fibroblasts, we believe nonapeptide 11 has high potential for further development as a mupirocin replacement candidate to treat skin infections caused by MRSA.

Structure and sequence features of mussel adhesive protein lead to its salt-tolerant adhesion ability

Mussels can strongly adhere to hydrophilic minerals in sea habitats by secreting adhesive proteins. The adhesion ability of these proteins is often attributed to the presence of Dopa derived from posttranslational modification of Tyr, whereas the contribution of structural feature is overlooked. It remains largely unknown how adhesive proteins overcome the surface-bound water layer to establish underwater adhesion. Here, we use molecular dynamics simulations to probe the conformations of adhesive protein Pvfp-5β and its salt-tolerant underwater adhesion on superhydrophilic mica. Dopa and positively charged basic residues form pairs, in this intrinsically disordered protein, and these residue pairs can lead to firm surface binding. Our simulations further suggest that the unmodified Tyr shows similar functions on surface adhesion by forming pairing structure with a positively charged residue. We confirm the presence of these residue pairs and verify the strong binding ability of unmodified proteins using nuclear magnetic resonance spectroscopy and lap shear tests.

Capturing Choline-Aromatics Cation-π Interactions in the MARTINI Force Field

Cation−π interactions play an important role in biomolecular recognition, including interactions between membrane phosphatidylcholine lipids and aromatic amino acids of peripheral proteins. While molecular mechanics coarse grain (CG) force fields are particularly well suited to simulate membrane proteins in general, they are not parameterized to explicitly reproduce cation−π interactions. We here propose a modification of the polarizable MARTINI coarse grain (CG) model enabling it to model membrane binding events of peripheral proteins whose aromatic amino acid interactions with choline headgroups are crucial for their membrane binding. For this purpose, we first collected and curated a dataset of eight peripheral proteins from different families. We find that the MARTINI CG model expectedly underestimates aromatics–choline interactions and is unable to reproduce membrane binding of the peripheral proteins in our dataset. Adjustments of the relevant interactions in the polarizable MARTINI force field yield significant improvements in the observed binding events. The orientation of each membrane-bound protein is comparable to reference data from all-atom simulations and experimental binding data. We also use negative controls to ensure that choline–aromatics interactions are not overestimated. We finally check that membrane properties, transmembrane proteins, and membrane translocation potential of mean force (PMF) of aromatic amino acid side-chain analogues are not affected by the new parameter set. This new version “MARTINI 2.3P” is a significant improvement over its predecessors and is suitable for modeling membrane proteins including peripheral membrane binding of peptides and proteins.

Improved Modeling of Cation-π and Anion-Ring Interactions Using the Drude Polarizable Empirical Force Field for Proteins

Cation-π interactions are noncovalent interactions between a π-electron system and a positively charged ion that are regarded as a strong noncovalent interaction and are ubiquitous in biological systems. Similarly, though less studied, anion-ring interactions are present in proteins along with in-plane interactions of anions with aromatic rings. As these interactions are between a polarizing ion and a polarizable π system, the accuracy of the treatment of these interactions in molecular dynamics (MD) simulations using additive force fields (FFs) may be limited. In the present work, to allow for a better description of ion-π interactions in proteins in the Drude-2013 protein polarizable FF, we systematically optimized the parameters for these interactions targeting model compound quantum mechanical (QM) interaction energies with atom pair-specific Lennard-Jones parameters along with virtual particles as selected ring centroids introduced to target the QM interaction energies and geometries. Subsequently, MD simulations were performed on a series of protein structures where ion-π pairs occur to evaluate the optimized parameters in the context of the Drude-2013 FF. The resulting FF leads to a significant improvement in reproducing the ion-π pair distances observed in experimental protein structures, as well as a smaller root-mean-square differences and fluctuations of the overall protein structures from experimental structures. Accordingly, the optimized Drude-2013 protein polarizable FF is suggested for use in MD simulations of proteins where cation-π and anion-ring interactions are critical. © 2019 Wiley Periodicals, Inc.

A comprehensive computational study of amino acid interactions in membrane proteins

Transmembrane proteins play a fundamental role in a wide series of biological processes but, despite their importance, they are less studied than globular proteins, essentially because their embedding in lipid membranes hampers their experimental characterization. In this paper, we improved our understanding of their structural stability through the development of new knowledge-based energy functions describing amino acid pair interactions that prevail in the transmembrane and extramembrane regions of membrane proteins. The comparison of these potentials and those derived from globular proteins yields an objective view of the relative strength of amino acid interactions in the different protein environments, and their role in protein stabilization. Separate potentials were also derived from α-helical and β-barrel transmembrane regions to investigate possible dissimilarities. We found that, in extramembrane regions, hydrophobic residues are less frequent but interactions between aromatic and aliphatic amino acids as well as aromatic-sulfur interactions contribute more to stability. In transmembrane regions, polar residues are less abundant but interactions between residues of equal or opposite charges or non-charged polar residues as well as anion-π interactions appear stronger. This shows indirectly the preference of the water and lipid molecules to interact with polar and hydrophobic residues, respectively. We applied these new energy functions to predict whether a residue is located in the trans- or extramembrane region, and obtained an AUC score of 83% in cross validation, which demonstrates their accuracy. As their application is, moreover, extremely fast, they are optimal instruments for membrane protein design and large-scale investigations of membrane protein stability.

Interfacial Aromatics Mediating Cation-π Interactions with Choline-Containing Lipids Can Contribute as Much to Peripheral Protein Affinity for Membranes as Aromatics Inserted below the Phosphates

Membrane-binding interfaces of peripheral proteins are restricted to a small part of their exposed surface, so the ability to engage in strong selective interactions with membrane lipids at various depths in the interface, both below and above the phosphates, is an advantage. Driven by their hydrophobicity, aromatic amino acids preferentially partition into membrane interfaces, often below the phosphates, yet enthalpically favorable interactions with the lipid headgroups, above the phosphate plane, are likely to further stabilize high interfacial positions. Using free-energy perturbation, we calculate the energetic cost of alanine substitution for 11 interfacial aromatic amino acids from 3 peripheral proteins. We show that the involvement in cation-π interactions with the headgroups (i) increases the ΔΔG_transfer as compared with insertion at the same depth without cation-π stabilization and (ii) can contribute at least as much as deeper insertion below the phosphates, highlighting the multiple roles of aromatics in peripheral membrane protein affinity.

Interaction of monomeric and self-assembled aromatic amino acids with model membranes: self-reproduction phenomena

The spontaneous formation of amyloid structures of proteins is responsible for several major human neurodegenerative diseases.

Cation-π Interactions between Methylated Ammonium Groups and Tryptophan in the CHARMM36 Additive Force Field

Cation-π interactions between tryptophan and choline or trimethylated lysines are vital for many biological processes. The performance of the additive CHARMM36 force field against target quantum mechanical data is shown to reproduce QM equilibrium geometries but required modified Lennard-Jones potentials to accurately reproduce the QM interaction energies. The modified parameter set allows accurate modeling, including free energies, of cation-π indole-choline and indole-trimethylated lysines interactions relevant for protein-ligand, protein-membrane, and protein-protein interfaces.

Augmentation of the antibacterial activities of Pt5-derived antimicrobial peptides (AMPs) by amino acid substitutions: Design of novel AMPs against MDR bacteria

The ever-growing concerns on multi-drug resistant (MDR) bacteria lead to urgent demands for novel antibiotics including antimicrobial peptides (AMPs). Pt5, a peptide consisting of the C-terminal 55 residues of zebrafish phosvitin, has been shown to function as an antibacterial agent. Here we used Pt5 as a template to design new AMPs by shortening the sequence and substituting with tryptophan (W) and lysine (K) at selected positions. Among the resultant Pt5-derived peptides, Pt5-1c showed the strongest antimicrobial activity against both Gram-negative and Gram-positive bacteria, including MDR bacteia, with the minimum inhibitory concentrations (MICs) ranging from 1.2 μM to 4.8 μM. Electron microscopic examination showed that Pt5-1c was able to kill the bacteria directly. ELISA revealed that Pt5-1c possessed high affinity to lipopolysaccharide (LPS), lipoteichoic acid (LTA) and peptidoglycan (PGN). Importantly, Pt5-1c was able to disrupt the bacterial membrane by a combined action of membrane depolarization and permeabilization, with little cytotoxicity to mammalian cells. Taken together, these findings suggest that Pt5-1c has considerable potential for future development as novel peptide antibiotics against MDR bacteria.

Structure-activity relationship studies of ultra-short peptides with potent activities against fluconazole-resistant Candida albicans

Vulvovaginal candidiasis (VVC) is a genital fungal infection afflicting approximately 75% of women globally and is primarily caused by the yeast Candida albicans. The extensive use of fluconazole, the first-line antifungal drug of choice, has led to the emergence of fluconazole-resistant C. albicans, creating a global clinical concern. This, coupled to the lack of new antifungal drugs entering the market over the past decade, has made it imperative for the introduction of new antifungal drug classes. Peptides with antifungal properties are deemed potential drug candidates due to their rapid membrane-disrupting mechanism of action. By specifically targeting and rapidly disrupting fungal membranes, they reduce the chances of resistance development and treatment duration. In a previous screening campaign involving an antimicrobial peptide library, we identified an octapeptide (IKIKIKIK-NH₂) with potent activity against C. albicans. Herein, we report a structure-activity relationship study on this peptide with the aim of designing a more potent peptide for further development. The lead peptide was then tested against a panel of fluconazole-resistant C. albicans, subjected to a fungicidal/static determination assay, a human dermal fibroblast viability assay and a homozygous profiling assay to gain insights into its mechanism of action and potential for further development as a topical antifungal agent.

Membrane Anchoring of α-Helical Proteins: Role of Tryptophan

The function of membrane proteins relies on a defined orientation of protein relative to lipid. In apparent correlation to protein anchoring, tryptophan residues are enriched in the lipid headgroup region. To characterize the thermodynamic and structural basis of this relationship in α-helical membrane proteins, we examined the role of three conserved tryptophans in the folding of the heterodimeric integrin αIIbβ3 transmembrane (TM) complex in phospholipid bicelles and mammalian membranes. In the homogenous lipid environment of bicelles, tryptophan was replaceable by residues of distinct polarities. The appropriate polarity was guided by the electrostatic potential of the tryptophan surrounding, suggesting that tryptophan can complement diverse environments by adjusting the orientation of its anisotropic side chain to achieve site-specific anchoring. As a sole membrane anchor, tryptophan made a contribution of 0.4 kcal/mol to TM complex stability in bicelles. In membranes, it proved more difficult to replace tryptophan even by tyrosine, indicating a superior capacity to interact with heterogeneous lipids of biological membranes. Interestingly, at intracellular TM helix ends, where integrin activation is initiated, sequence motifs that interact with lipids via opposing polarity patterns were found to restrict TM helix orientations beyond tryptophan anchoring. In contrast to bicelles, phenylalanine became the least accepted substitute in membranes, demonstrating an increased role of the hydrophobic effect. Altogether, our study implicates a wide amphiphilic range of tryptophan, membrane complexity, and the hydrophobic effect to be important factors in tryptophan membrane anchoring.

Antifungal peptides: a potential new class of antifungals for treating vulvovaginal candidiasis caused by fluconazole-resistant Candida albicans

Vulvovaginal candidiasis/candidosis is a common fungal infection afflicting approximately 75% of women globally caused primarily by the yeast Candida albicans. Fluconazole is widely regarded as the antifungal drug of choice since its introduction in 1990 due to its high oral bioavailability, convenient dosing regimen and favourable safety profile. However, its widespread use has led to the emergence of fluconazole-resistant C. albicans, posing a universal clinical concern. Coupled to the dearth of new antifungal drugs entering the market, it is imperative to introduce new drug classes to counter this threat. Antimicrobial peptides (AMPs) are potential candidates due to their membrane-disrupting mechanism of action. By specifically targeting fungal membranes and being rapidly fungicidal, they can reduce the chances of resistance development and treatment duration. Towards this goal, we conducted a head-to-head comparison of 61 short linear AMPs from the literature to identify the peptide with the most potent activity against fluconazole-resistant C. albicans. The 11-residue peptide, P11-6, was identified and assayed against a panel of clinical C. albicans isolates followed by fungicidal/static determination and a time-kill assay to gauge its potential for further drug development. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

Improving the Force Field Description of Tyrosine-Choline Cation-π Interactions: QM Investigation of Phenol-N(Me)4+ Interactions

Cation-π interactions between tyrosine amino acids and compounds containing N,N,N-trimethylethanolammonium (N(CH3)3) are involved in the recognition of histone tails by chromodomains and in the recognition of phosphatidylcholine (PC) phospholipids by membrane-binding proteins. Yet, the lack of explicit polarization or charge transfer effects in molecular mechanics force fields raises questions about the reliability of the representation of these interactions in biomolecular simulations. Here, we investigate the nature of phenol-tetramethylammonium (TMA) interactions using quantum mechanical (QM) calculations, which we also use to evaluate the accuracy of the additive CHARMM36 and Drude polarizable force fields in modeling tyrosine-choline interactions. We show that the potential energy surface (PES) obtained using SAPT2+/aug-cc-pVDZ compares well with the large basis-set CCSD(T) PES when TMA approaches the phenol ring perpendicularly. Furthermore, the SAPT energy decomposition reveals comparable contributions from electrostatics and dispersion in phenol-TMA interactions. We then compared the SAPT2+/aug-cc-pVDZ PES obtained along various approach directions to the corresponding PES obtained with CHARMM, and we show that the force field accurately reproduces the minimum distances while the interaction energies are underestimated. The use of the Drude polarizable force field significantly improves the interaction energies but decreases the agreement on distances at energy minima. The best agreement between force field and QM PES is obtained by modifying the Lennard-Jones terms for atom pairs involved in the phenol-TMA cation-π interactions. This is further shown to improve the correlation between the occupancy of tyrosine-choline cation-π interactions obtained from molecular dynamics simulations of a bilayer-bound bacterial phospholipase and experimental affinity data of the wild-type protein and selected mutants.

Highly Dynamic Anion-Quadrupole Networks in Proteins

The dynamics of anion-quadrupole (or anion-π) interactions formed between negatively charged (Asp/Glu) and aromatic (Phe) side chains are for the first time computationally characterized in RmlC (Protein Data Bank entry 1EP0 ), a homodimeric epimerase. Empirical force field-based molecular dynamics simulations predict anion-quadrupole pairs and triplets (anion-anion-π and anion-π-π) are formed by the protein during the simulated trajectory, which suggests that the anion-quadrupole interactions may provide a significant contribution to the overall stability of the protein, with an average of -1.6 kcal/mol per pair. Some anion-π interactions are predicted to form during the trajectory, extending the number of anion-quadrupole interactions beyond those predicted from crystal structure analysis. At the same time, some anion-π pairs observed in the crystal structure exhibit marginal stability. Overall, most anion-π interactions alternate between an "on" state, with significantly stabilizing energies, and an "off" state, with marginal or null stabilizing energies. The way proteins possibly compensate for transient loss of anion-quadrupole interactions is characterized in the RmlC aspartate 84-phenylalanine 112 anion-quadrupole pair observed in the crystal structure. A double-mutant cycle analysis of the thermal stability suggests a possible loss of anion-π interactions compensated by variations of hydration of the residues and formation of compensating electrostatic interactions. These results suggest that near-planar anion-quadrupole pairs can exist, sometimes transiently, which may play a role in maintaining the structural stability and function of the protein, in an otherwise very dynamic interplay of a nonbonded interaction network as well as solvent effects.

Position and orientational preferences of drug-like compounds in lipid membranes: a computational and NMR approach

Permeation of drugs across lipid bilayers is a key factor in dictating how effective they will be. In vivo, the issue is compounded by the presence of drug-exporter proteins such as P-glycoprotein. However, despite intense effort, exactly what controls permeation and susceptibility to export is still poorly understood. In this work we examine two well-studied drugs for which interaction with P-glycoprotein has been studied before: amitriptyline, a known substrate and clozapine, which is not a substrate. Extensive MD simulations, including potential of mean force (PMF) profiles of the compounds in all possible protonation states, reveal that the preferred location of the compounds in different bilayers in different protonation states is remarkably similar. For both molecules in charged states, there is a substantial barrier to crossing the bilayer. Clozapine however, shows an energetic barrier to movement across the bilayer even in a protonation state that results in an uncharged molecule. For amitriptyline there is only a very small barrier of approximately 1.3 kcal mol(-1). Further analysis revealed that the conformational and orientational behavior of the two compounds was also similar, with the sidechain interacting with the lipid headgroups. This effect was much stronger if the sidechain was charged (protonated). These interactions with lipid bilayers were confirmed by NMR ROESY experiments. The results are discussed in terms of their potential interactions with export proteins like P-glycoprotein.

Gramicidin A disassembles large conductive clusters of its lysine-substituted derivatives in lipid membranes

N-terminally substituted lysine derivatives of gramicidin A (gA), [Lys1]gA and [Lys3]gA, but not glutamate- or aspartate-substituted peptides have been previously shown to cause the leakage of carboxyfluorescein from liposomes. Here, the leakage induction was also observed for [Arg1]gA and [Arg3]gA, while [His1]gA and [His3]gA were inactive at neutral pH. The Lys3-containing analogue with all tryptophans replaced by isoleucines did not induce liposome leakage, similar to gA. This suggests that the presence of both tryptophans and N-terminal cationic residues is critical for pore formation. Remarkably, the addition of gA blocked the leakage induced by [Lys3]gA. By examining with fluorescence correlation spectroscopy the peptide-induced leakage of fluorescent markers from liposomes, we estimated the diameter of pores responsible for the leakage to be about 1.6 nm. Transmission electron cryo-microscopy imaging of liposomes with [Lys3]gA showed that the liposomal membranes contained high electron density particles with a size of about 40 Å, suggesting the formation of peptide clusters. No such clusterization was observed in liposomes incorporating gA or a mixture of gA with [Lys3]gA. Three-dimensional reconstruction of the clusters was compatible with their pentameric arrangement. Based on experimental data and computational modeling, we suggest that the large pore formed by [Lys3]gA represents a barrel-stave oligomeric cluster formed by antiparallel double-stranded helical dimers (DH). In a tentative model, the pentamer of dimers may be stabilized by aromatic Trp-Trp and cation-π Trp-Lys interactions between the neighboring DHs. The inhibiting effect of gA on the [Lys3]gA-induced leakage can be attributed to breaking of cation-π interactions, which prevents peptide clusterization and pore formation.

Impact of cation-π interactions on the cell voltage of carbon nanotube-based Li batteries

Carbon nanotube (CNT)-based Li batteries have attracted wide attention because of their high capacity, high cyclability and high energy density and are believed to be one of the most promising electrochemical energy storage systems. In CNT-based Li batteries, the main interaction between the Li(+) ions and the CNT is the cation-π interaction. However, up to now, it is still not clear how this interaction affects the storage characteristics of CNT-based Li batteries. Here, using density functional theory (DFT) calculations, we report a highly favorable impact of cation-π interactions on the cell voltage of CNT-based Li batteries. Considering both Li(+)-π interaction and Li-π interaction, we show that cell voltage enhances with the increase of the CNT diameter. In addition, when the Li(+) ion adsorbs on the external wall, the cell voltage is larger than that when it adsorbs on the internal wall. This suggests that CNTs with a large diameter and a low array density are more advantageous to enhance storage performance of CNT-based Li batteries. Compared with Li(+) ions on the (4,4) CNT internal wall, the cell voltage of Li(+) on the (10,10) CNT external wall is 0.55 V higher, which indicates an improvement of about 38%. These results will be helpful for the design of more efficient CNT-based Li batteries.

Archetypal tryptophan-rich antimicrobial peptides: properties and applications

Drug-resistant microorganisms ('superbugs') present a serious challenge to the success of antimicrobial treatments. Subsequently, there is a crucial need for novel bio-control agents. Many antimicrobial peptides (AMPs) show a broad-spectrum activity against bacteria, fungi or viruses and are strong candidates to complement or substitute current antimicrobial agents. Some AMPs are also effective against protozoa or cancer cells. The tryptophan (Trp)-rich peptides (TRPs) are a subset of AMPs that display potent antimicrobial activity, credited to the unique biochemical properties of tryptophan that allow it to insert into biological membranes. Further, many Trp-rich AMPs cross bacterial membranes without compromising their integrity and act intracellularly, suggesting interactions with nucleic acids and enzymes. In this work, we overview some archetypal TRPs derived from natural sources, i.e., indolicidin, tritrpticin and lactoferricin, summarising their biochemical properties, structures, antimicrobial activities, mechanistic studies and potential applications.

Modeling Aromatic Liquids: Toluene, Phenol, and Pyridine

Aromatic groups are now acknowledged to play an important role in many systems of interest. However, existing molecular mechanics methods provide a poor representation of these groups. In a previous paper, we have shown that the molecular mechanics treatment of benzene can be improved by the incorporation of an explicit representation of the aromatic π electrons. Here, we develop this concept further, developing charge-separation models for toluene, phenol, and pyridine. Monte Carlo simulations are used to parametrize the models, via the reproduction of experimental thermodynamic data, and our models are shown to outperform an existing atom-centered model. The models are then used to make predictions about the structures of the liquids at the molecular level and are tested further through their application to the modeling of gas-phase dimers and cation-π interactions.

Identification, Functional Characterization, and Evolution of Terpene Synthases from a Basal Dicot

Bay laurel (Laurus nobilis) is an agriculturally and economically important dioecious tree in the basal dicot family Lauraceae used in food and drugs and in the cosmetics industry. Bay leaves, with their abundant monoterpenes and sesquiterpenes, are used to impart flavor and aroma to food, and have also drawn attention in recent years because of their potential pharmaceutical applications. To identify terpene synthases (TPSs) involved in the production of these volatile terpenes, we performed RNA sequencing to profile the transcriptome of L. nobilis leaves. Bioinformatic analysis led to the identification of eight TPS complementary DNAs. We characterized the enzymes encoded by three of these complementary DNAs: a monoterpene synthase that belongs to the TPS-b clade catalyzes the formation of mostly 1,8-cineole; a sesquiterpene synthase belonging to the TPS-a clade catalyzes the formation of mainly cadinenes; and a diterpene synthase of the TPS-e/f clade catalyzes the formation of geranyllinalool. Comparison of the sequences of these three TPSs indicated that the TPS-a and TPS-b clades of the TPS gene family evolved early in the evolution of the angiosperm lineage, and that geranyllinalool synthase activity is the likely ancestral function in angiosperms of genes belonging to an ancient TPS-e/f subclade that diverged from the kaurene synthase gene lineages before the split of angiosperms and gymnosperms.

Discovery of an ultra-short linear antibacterial tetrapeptide with anti-MRSA activity from a structure-activity relationship study

The overuse and misuse of antibiotics has resulted in the emergence of drug-resistant pathogenic bacteria, including meticillin-resistant Staphylococcus aureus (MRSA), the primary pathogen responsible for human skin and soft-tissue infections. Antibacterial peptides are known to kill bacteria by rapidly disrupting their membranes and are deemed plausible alternatives to conventional antibiotics. One advantage of their membrane-targeting mode of action is that bacteria are unlikely to develop resistance as changing their cell membrane structure and morphology would likely involve extensive genetic mutations. However, major concerns in using peptides as antibacterial drugs include their instability towards plasma proteases, toxicity towards human cells due to their membrane-targeting mode of action and high manufacturing cost. These concerns can be mitigated by developing peptides as topical agents, by the judicial selection of amino acids and developing very short peptides respectively. In this preliminary report, we reveal a linear, non-hemolytic tetrapeptide with rapid bactericidal activity against MRSA developed from a structure-activity relationship study based on the antimicrobial hexapeptide WRWRWR-NH2. Our finding opens promising avenues for the development of ultra-short antibacterials to treat multidrug-resistant MRSA skin and soft tissue infections.

Ortho-7 bound to the active-site gorge of free and OP-conjugated acetylcholinesterase: cation-π interactions

Despite the immense importance of cation-π interactions prevailing in bispyridinium drug acetylcholinesterase (AChE) complexes, a precise description of cation-π interactions at molecular level has remained elusive. Here, we consider a bispyridinium drug, namely, ortho-7 in three different structures of AChE, with and without complexation with organophosphorus (OP) compounds for detailed investigation using all atom molecular dynamics simulation. By quantum mechanical calculations, Y72, W86, Y124, W286, Y337, and Y341 aromatic residues of the enzyme are investigated for possible cation-π interactions with ortho-7. The cation-π interactions in each of the protein-drug complexes are studied using distance, angle, a suitable functional form of them, and electrostatic criteria. The variation of cation-π functional is remarkably consistent with that of the Columbic variation. It is clearly observed that cation-π interactions for some of the residues in the catalytic active site (CAS) and peripheral anionic site (PAS) of the enzyme are either enhanced or reduced based on the nature of OP conjugation (i.e., nerve gas, tabun or pesticide, fenamiphos) when compared with the OP-free enzyme. The strength of cation-π interaction is strongly dependent on the type OP conjugation. The effect of conjugation at CAS is also seen to influence the cation-π interaction at the PAS region. The variation of cation-π interactions on the type of conjugating OP compounds might be suggestive of a reason as to why wide spectrum drug against any OP poisoning is yet to arrive in the market.

Aromatic Interactions in Organocatalyst Design: Augmenting Selectivity Reversal in Iminium Ion Activation

Substituting N-methylpyrrole for N-methyindole in secondary-amine-catalysed Friedel-Crafts reactions leads to a curious erosion of enantioselectivity. In extreme cases, this substrate dependence can lead to an inversion in the sense of enantioinduction. Indeed, these closely similar transformations require two structurally distinct catalysts to obtain comparable selectivities. Herein a focussed molecular editing study is disclosed to illuminate the structural features responsible for this disparity, and thus identify lead catalyst structures to further exploit this selectivity reversal. Key to effective catalyst re-engineering was delineating the non-covalent interactions that manifest themselves in conformation. Herein we disclose preliminary validation that intermolecular aromatic (CH-π and cation-π) interactions between the incipient iminium cation and the indole ring system is key to rationalising selectivity reversal. This is absent in the N-methylpyrrole alkylation, thus forming the basis of two competing enantio-induction pathways. A simple L-valine catalyst has been developed that significantly augments this interaction.

Amphipathic solvation of indole: implications for the role of tryptophan in membrane proteins

The microscopic structure of the tryptophan side chain, indole, in an amphiphilic environment has been investigated using a combination of neutron diffraction measurements and simulations in solution. The results show that indole is preferentially solvated by hydrogen bonding interactions between water and alcohol -OH groups rather than the interaction being dominated by indole-methyl interactions. This has implications for understanding how tryptophan interacts with the amphipathic membrane environment to anchor proteins into membranes, where the results here suggest that the benzene ring of tryptophan interacts directly with the interfacial water at the membrane surface rather than being buried into the hydrophobic regions of the membrane bilayer.

Identification of Tyr residues that enhance folate substrate binding and constrain oscillation of the proton-coupled folate transporter (PCFT-SLC46A1)

The proton-coupled folate transporter (PCFT) mediates intestinal folate absorption and transport of folates across the choroid plexus. This study focuses on the role of Tyr residues in PCFT function. The substituted Cys-accessibility method identified four Tyr residues (Y291, Y362, Y315, and Y414) that are accessible to the extracellular compartment; three of these (Y291, Y362, and Y315) are located within or near the folate binding pocket. When the Tyr residues were replaced with Cys or Ala, these mutants showed similar (up to 6-fold) increases in influx Vmax and Kt/Ki for [(3)H]methotrexate and [(3)H]pemetrexed. When the Tyr residues were replaced with Phe, these changes were moderated or absent. When Y315A PCFT was used as representative of the mutants and [(3)H]pemetrexed as the transport substrate, this substitution did not increase the efflux rate constant. Furthermore, neither influx nor efflux mediated by Y315A PCFT was transstimulated by the presence of substrate in the opposite compartment; however, substantial bidirectional transstimulation of transport was mediated by wild-type PCFT. This resulted in a threefold greater efflux rate constant for cells that express wild-type PCFT than for cells that express Y315 PCFT under exchange conditions. These data suggest that these Tyr residues, possibly through their rigid side chains, secure the carrier in a high-affinity state for its folate substrates. However, this may be achieved at the expense of constraining the carrier's mobility, thereby decreasing the rate at which the protein oscillates between its conformational states. The Vmax generated by these Tyr mutants may be so rapid that further augmentation during transstimulation may not be possible.

Two homologous neutrophil serine proteases bind to POPC vesicles with different affinities: When aromatic amino acids matter

Neutrophil serine proteases Proteinase 3 (PR3) and human neutrophil elastase (HNE) are homologous antibiotic serine proteases of the polymorphonuclear neutrophils. Despite sharing a 56% sequence identity they have been shown to have different functions and localizations in the neutrophils. In particular, and in contrast to HNE, PR3 has been detected at the outer leaflet of the plasma membrane and its membrane expression is a risk factor in a number of chronic inflammatory diseases. Although a plethora of studies performed in various cell-based assays have been reported, the mechanism by which PR3, and possibly HNE bind to simple membrane models remains unclear. We used surface plasmon resonance (SPR) experiments to measure and compare the affinity of PR3 and HNE for large unilamellar vesicles composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). We also conducted 500-nanosecond long molecular dynamics simulations of each enzyme at the surface of a POPC bilayer to map the interactions between proteins and lipids and rationalize the difference in affinity observed in the SPR experiment. We find that PR3 binds strongly to POPC large unilamellar vesicles (Kd=9.2×10(-7)M) thanks to the insertion of three phenylalanines, one tryptophan and one leucine beyond the phosphate groups of the POPC lipids. HNE binds in a significantly weaker manner (Kd>10(-5)M) making mostly electrostatic interactions via lysines and arginines and inserting only one leucine between the hydrophobic lipid tails. Our results support the early reports that PR3, unlike HNE, is able to directly and strongly anchor directly to the neutrophil membrane.

Importance of indole N-H hydrogen bonding in the organization and dynamics of gramicidin channels

The linear ion channel peptide gramicidin represents an excellent model for exploring the principles underlying membrane protein structure and function, especially with respect to tryptophan residues. The tryptophan residues in gramicidin channels are crucial for the structure and function of the channel. In order to test the importance of indole hydrogen bonding for the biophysical properties of gramicidin channels, we monitored the effect of N-methylation of gramicidin tryptophans, using a combination of steady state and time-resolved fluorescence approaches along with circular dichroism spectroscopy. We show here that in the absence of the hydrogen bonding ability of tryptophans, tetramethyltryptophan gramicidin (TM-gramicidin) is unable to maintain the single stranded, head-to-head dimeric channel conformation in membranes. Our results show that TM-gramicidin displays a red-shifted fluorescence emission maximum, lower red edge excitation shift (REES), and higher fluorescence intensity and lifetime, consistent with its nonchannel conformation. This is in agreement with the measured location (average depth) of the 1-methyltryptophans in TM-gramicidin using the parallax method. These results bring out the usefulness of 1-methyltryptophan as a fluorescent tool to examine the hydrogen bonding ability of tryptophans in proteins and peptides. We conclude that changes in the hydrogen bonding ability of tryptophans, along with coupled changes in peptide backbone structure induce the loss of single stranded β(6.3) helical dimer conformation. These results agree with earlier results from size-exclusion chromatography and single-channel measurements for TM-gramicidin, and confirm the importance of indole hydrogen bonding for the conformation and function of ion channels and membrane proteins.