Buried lysine, but not arginine, titrates and alters transmembrane helix tilt

l-Phenylalanine Partitioning Mechanisms in Model Biological Membranes

Time-resolved fluorescence spectroscopy in combination with differential scanning calorimetry (DSC) was used to study the chemical interactions that occur when l-phenylalanine is introduced to solutions containing phosphatidylcholine vesicles. Studies reported in this work address open questions about l-Phe's affinity for lipid vesicle bilayers, the effects of l-Phe partitioning on bilayer properties, l-Phe's solvation within a lipid bilayer, and the amount of l-Phe within that local solvation environment. DSC data show that l-Phe reduces the amount of heat necessary to melt saturated phosphatidylcholine bilayers from their gel to liquid-crystalline state but does not change the transition temperature (T_gel-lc). Time-resolved emission shows only a single l-Phe lifetime at low temperatures corresponding to l-Phe remaining solvated in aqueous solution. At temperatures close to T_gel-lc, a second, shorter lifetime appears that is assigned to l-Phe already embedded within the membrane that becomes hydrated as water starts to permeate the lipid bilayer. This new lifetime is attributed to a conformationally restricted rotamer in the bilayer's polar headgroup region and accounts for up to 30% of the emission amplitude. Results reported for dipalmitoylphosphatidylcholine (DPPC, 16:0) lipid vesicles prove to be general, with similar effects observed for dimyristoylphosphatidylcholine (DMPC, 14:0) and distearoylphosphatidylcholine (DSPC, 18:0) vesicles. Taken together, these results create a complete and compelling picture of how l-Phe associates with model biological membranes. Furthermore, this approach to examining amino acid partitioning into membranes and the resulting solvation forces points to new strategies for studying the structure and chemistry of membrane-soluble peptides and selected membrane proteins.

Improving the Therapeutic Index of Smp24, a Venom-Derived Antimicrobial Peptide: Increased Activity against Gram-Negative Bacteria

Antimicrobial peptides (AMPs) are naturally occurring compounds which possess a rapid killing mechanism and low resistance potential. Consequently, they are being viewed as potential alternatives to traditional antibiotics. One of the major factors limiting further development of AMPs is off-target toxicity. Enhancements to antimicrobial peptides which can maximise antimicrobial activity whilst reducing mammalian cytotoxicity would make these peptides more attractive as future pharmaceuticals. We have previously characterised Smp24, an AMP derived from the venom of the scorpion Scorpio maurus palmatus. This study sought to better understand the relationship between the structure, function and bacterial selectivity of this peptide by performing single amino acid substitutions. The antimicrobial, haemolytic and cytotoxic activity of modified Smp24 peptides was determined. The results of these investigations were compared with the activity of native Smp24 to determine which modifications produced enhanced therapeutic indices. The structure–function relationship of Smp24 was investigated by performing N-terminal, mid-chain and C-terminal amino acid substitutions and determining the effect that they had on the antimicrobial and cytotoxic activity of the peptide. Increased charge at the N-, mid- and C-termini of the peptide resulted in increased antimicrobial activity. Increased hydrophobicity at the N-terminus resulted in reduced haemolysis and cytotoxicity. Reduced antimicrobial, haemolytic and cytotoxic activity was observed by increased hydrophobicity at the mid-chain. Functional improvements have been made to modified peptides when compared with native Smp24, which has produced peptides with enhanced therapeutic indices.

The Dynamic Range of Acidity: Tracking Rules for the Unidirectional Penetration of Cellular Compartments

Labeled ammonium cations with pK_a∼7.4 accumulate in acidic organelles because they can be neutralized transiently to cross the membrane at cytosolic pH 7.2 but not at their internal pH<5.5. Retention in early endosomes with less acidic internal pH was achieved recently using weaker acids of up to pK_a 9.8. We report here that primary ammonium cations with higher pK_a 10.6, label early endosomes more efficiently. This maximized early endosome tracking coincides with increasing labeling of Golgi networks with similarly weak internal acidity. Guanidinium cations with pK_a 13.5 cannot cross the plasma membrane in monomeric form and label the plasma membrane with selectivity for vesicles embarking into endocytosis. Self‐assembled into micelles, guanidinium cations enter cells like arginine‐rich cell‐penetrating peptides and, driven by their membrane potential, penetrate mitochondria unidirectionally despite their high inner pH. The resulting tracking rules with an approximated dynamic range of pK_a change ∼3.5 are expected to be generally valid, thus enabling the design of chemistry tools for biology research in the broadest sense. From a practical point of view, most relevant are two complementary fluorescent flipper probes that can be used to image the mechanics at the very beginning of endocytosis.

H1 helix of colicin U causes phospholipid membrane permeation

In light of an increasing number of antibiotic-resistant bacterial strains, it is essential to understand an action imposed by various antimicrobial agents on bacteria at the molecular level. One of the leading mechanisms of killing bacteria is related to the alteration of their plasmatic membrane. We study bio-inspired peptides originating from natural antimicrobial proteins colicins, which can disrupt membranes of bacterial cells. Namely, we focus on the α-helix H1 of colicin U, produced by bacterium Shigella boydii, and compare it with analogous peptides derived from two different colicins. To address the behavior of the peptides in biological membranes, we employ a combination of molecular simulations and experiments. We use molecular dynamics simulations to show that all three peptides are stable in model zwitterionic and negatively charged phospholipid membranes. At the molecular level, their embedment leads to the formation of membrane defects, membrane permeation for water, and, for negatively charged lipids, membrane poration. These effects are caused by the presence of polar moieties in the considered peptides. Importantly, simulations demonstrate that even monomeric H1 peptides can form toroidal pores. At the macroscopic level, we employ experimental co-sedimentation and fluorescence leakage assays. We show that the H1 peptide of colicin U incorporates into phospholipid vesicles and disrupts their membranes, causing leakage, in agreement with the molecular simulations. These insights obtained for model systems seem important for understanding the mechanisms of antimicrobial action of natural bacteriocins and for future exploration of small bio-inspired peptides able to disrupt bacterial membranes.

Exploring the trimerization process of a transmembrane helix with an ionizable residue by molecular dynamics simulations: a case study of transmembrane domain 5 of LMP-1

The oligomerization of membrane proteins is an important biological process that plays a critical role in the initialization of membrane protein receptor signaling. Unveiling how transmembrane segments oligomerize is critical for understanding the mechanism of membrane receptor signaling activation. Owing to the complicated membrane environment and the extraordinary dynamic properties of the ionizable residues in the transmembrane segment, it is extremely challenging to thoroughly understand the oligomerization process of the transmembrane domain. In this study, transmembrane domain 5 (TMD5) of latent membrane protein-1 from Epstein-Barr virus was used as a prototype model to investigate the trimerization process of the transmembrane segment with ionizable residues. The trimerization process of TMD5 was rebuilt and investigated via conventional molecular dynamics simulations and constant-pH molecular dynamics simulations. When TMD5s approached each other, the tilting angles of the TMD5 monomer decreased. TMD5s formed stable trimers until two interacting sites (D150s and Q139s) along each transmembrane helix were created to lock the TMD5s. The pK_a values of D150 shifted toward neutral states in the membrane environment. When TMD5s were monomers, the pK_a shift of D150 was mainly influenced by its microenvironment in the lipid bilayer. When TMD5s were moving close to each other, protein-protein interactions became the main contributing factor for the pK_a shift of D150s. Overall, this work elucidates the behavior of the TMD5 helix and the pK_a shift of ionizable residue D150 in the process of TMD5 oligomerization. This study may provide insight into the development of agents for targeting the oligomerization of membrane proteins.

Optimal pH 8.5 to 9 for the Hydrolysis of Vixotrigine and Other Basic Substrates of Carboxylesterase-1 in Human Liver Microsomes

Vixotrigine is a voltage- and use-dependent sodium channel blocker under investigation for the potential treatment of neuropathic pain. One of the major in vivo metabolic pathways of vixotrigine in humans is the hydrolysis of the carboxamide to form the carboxylic acid metabolite M14.The in vitro formation of M14 in human hepatocytes was inhibited by the carboxylesterase (CES) inhibitor Bis(4-nitrophenyl) phosphate in a concentration-dependent manner. The hydrolysis reaction was identified to be catalyzed by recombinant human CES1b.Initial observation of only trace level formation of M14 in human liver microsomes at pH 7.4 caused us to doubt the involvement of CES1, an enzyme localized at the endoplasmic reticulum and the dominant carboxylesterase in human liver. Further investigation has revealed that optimal pH for the hydrolysis of vixotrigine and two other basic substrates of CES1, methylphenidate and oseltamivir, in human liver microsomes was pH 8.5 to 9 which is higher than their respective pK_a(base), suggesting that neutral form of basic substrates is probably preferred for CES1 catalysis in liver microsomes.

BTK-inhibitor drug covalent binding to lysine in human serum albumin using LC-MS/MS

Irreversible Bruton's tyrosine kinase (BTK) inhibitor drugs are designed to bind covalently to a free-thiol cysteine in the BTK protein active site. However, these reactive drugs bind to off-target proteins as well. In this study, seven BTK-inhibitor drugs containing acrylamide warheads were incubated with human serum albumin (HSA) and analyzed using an LC-MS/MS peptide mapping approach to determine the amino acid sites of drug covalent binding. Significant adduction at the free-thiol cysteine of HSA was only observed for two of the drugs. However, significant adduction was observed for at least four lysine residues. This is just a small percentage of the 59 total lysine residues in HSA. These four lysine residues are likely partially buried, accessible to the drugs, and exist at least partially in a neutral state. The levels of adduction observed in the in-vitro experimental conditions are only indicative of a relative propensity for adduction with the individual lysine residues of HSA, and are not in-vivo predictions. Widespread off-target lysine binding could impact clearance and bioavailability for irreversible inhibitor drugs. However, the extent of the impact on clearance may be limited in comparison to conjugation with glutathione.

Illuminating Disorder Induced by Glu in a Stable Arg-Anchored Transmembrane Helix

Membrane proteins are vital for biological function and are complex to study. Even in model peptide-lipid systems, the combined influence or interaction of pairs of chemical groups still is not well understood. Disordered proteins, whether in solution or near lipid membranes, are an emerging paradigm for the initiation and control of biological function. The disorder can involve molecular orientation as well as molecular folding. This paper reports an astonishing induction of disorder when one Glu residue is introduced into a highly stable 23-residue transmembrane helix. The parent helix is anchored by a single Arg residue, tilted at a well-defined angle with respect to the DOPC bilayer normal and undergoes rapid cone precession. When Glu is introduced two residues away from Arg, with 200° (or 160°) radial separation, the helix properties change radically to exhibit a multiplicity of three or more disordered states. The helix characteristics have been monitored by deuterium (²H) NMR spectroscopy as functions of the pH and lipid bilayer composition. The disordered multistate behavior of the (Glu, Arg)-containing helix varies with the lipid bilayer thickness and pH. The results highlight a fundamental induction of protein multistate properties by a single Glu residue in a lipid membrane environment.

Thermodynamics-Based Molecular Modeling of α-Helices in Membranes and Micelles

The Folding of Membrane-Associated Peptides (FMAP) method was developed for modeling α-helix formation by linear peptides in micelles and lipid bilayers. FMAP 2.0 identifies locations of α-helices in the amino acid sequence, generates their three-dimensional models in planar bilayers or spherical micelles, and estimates their thermodynamic stabilities and tilt angles, depending on temperature and pH. The method was tested for 723 peptides (926 data points) experimentally studied in different environments and for 170 single-pass transmembrane (TM) proteins with available crystal structures. FMAP 2.0 detected more than 95% of experimentally observed α-helices with an average error in helix end determination of around 2, 3, 4, and 5 residues per helix for peptides in water, micelles, bilayers, and TM proteins, respectively. Helical and nonhelical residue states were predicted with an accuracy from 0.86 to 0.96, and the Matthews correlation coefficient was from 0.64 to 0.88 depending on the environment. Experimental micelle- and membrane-binding energies and tilt angles of peptides were reproduced with a root-mean-square deviation of around 2 kcal/mol and 7°, respectively. The TM and non-TM states of hydrophobic and pH-triggered α-helical peptides in various lipid bilayers were reproduced in more than 95% of cases. The FMAP 2.0 web server (https://membranome.org/fmap) is publicly available to explore the structural polymorphism of antimicrobial, cell-penetrating, fusion, and other membrane-binding peptides, which is important for understanding the mechanisms of their biological activities.

Lipid-Dependent Titration of Glutamic Acid at a Bilayer Membrane Interface

The ionization properties of protein side chains in lipid-bilayer membranes will differ from the canonical values of side chains exposed to an aqueous solution. While the propensities of positively charged side chains of His, Lys, and Arg to release a proton in lipid membranes have been rather well characterized, the propensity for a negatively charged Glu side chain to receive a proton and achieve the neutral state in a bilayer membrane has been less well characterized. Indeed, the ionization of the glutamic acid side chain has been predicted to depend on its depth of burial in a lipid membrane but has been difficult to verify experimentally. To address the issue, we incorporated an interfacial Glu residue at position 4 of a distinct 23-residue transmembrane helix and used ²H NMR to examine the helix properties as a function of pH. We observe that the helix tilt and azimuthal rotation vary little with pH, but the extent of helix unraveling near residues 3 and 4 changes as the Glu residue E4 titrates. Remarkably, the ²H quadrupolar splitting for the side chain of alanine A3 responds to pH with an apparent pK _a of 4.8 in 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 6.3 in 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC), but is unchanged up to pH 8.0 in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in the presence of residue E4. With bilayers composed of alkali-stable ether-linked lipids, the side chain of A3 responds to pH with an apparent pK _a of 11.0 in the ether analogue of DOPC. These results suggest that the depth dependence of Glu ionization in lipid-bilayer membranes may be steeper than previously predicted or envisioned.

1H/13C/15N triple-resonance experiments for structure determinaton of membrane proteins by oriented-sample NMR

The benefits of triple-resonance experiments for structure determination of macroscopically oriented membrane proteins by solid-state NMR are discussed. While double-resonance ¹H/¹⁵N experiments are effective for structure elucidation of alpha-helical domains, extension of the method of oriented samples to more complex topologies and assessing side-chain conformations necessitates further development of triple-resonance (¹H/¹³C/¹⁵N) NMR pulse sequences. Incorporating additional spectroscopic dimensions involving ¹³C spin-bearing nuclei, however, introduces essential complications arising from the wide frequency range of the ¹H-¹³C dipolar couplings and ¹³C CSA (>20 ​kHz), and the presence of the ¹³C-¹³C homonuclear dipole-dipole interactions. The recently reported ROULETTE-CAHA pulse sequence, in combination with the selective z-filtering, can be used to evolve the structurally informative ¹H-¹³C dipolar coupling arising from the aliphatic carbons while suppressing the signals from the carbonyl and methyl regions. Proton-mediated magnetization transfer under mismatched Hartman-Hahn conditions (MMHH) can be used to correlate ¹³C and ¹⁵N nuclei in such triple-resonance experiments for the subsequent ¹⁵N detection. The recently developed pulse sequences are illustrated for n-acetyl Leucine (NAL) single crystal and doubly labeled Pf1 coat protein reconstituted in magnetically aligned bicelles. An interesting observation is that in the case of ¹⁵N-labeled NAL measured at ¹³C natural abundance, the triple (¹H/¹³C/¹⁵N) MMHH scheme predominantly gives rise to long-range intermolecular magnetization transfers from ¹³C to ¹⁵N spins; whereas direct Hartmann-Hahn ¹³C/¹⁵N transfer is entirely intramolecular. The presented developments advance NMR of oriented samples for structure determination of membrane proteins and liquid crystals.

Targeting the transmembrane domain 5 of latent membrane protein 1 using small molecule modulators

Protein-protein interactions (PPIs) play a critical role in living cells and represent promising targets for the drug discovery and life sciences communities. However, lateral transmembrane PPIs are difficult targets for small-molecule inhibitor development given less structural information is known and fewer ligand discovery methods have been explored compared to soluble proteins. In this study, the interactions of the transmembrane domain 5 (TMD-5) of latent membrane protein 1 (LMP-1) of Epstein-Barr virus (EBV) were disrupted by pentamidine derivatives to curb the committed step of EBV infection. A pentamidine derivative 2 with a 7-atom di-amide linker had the best activity whilst switching the amide regiochemistry in the linker influenced membrane permeability and abolished anti TMD-5 activity. Molecular dynamics simulations were performed to understand the interaction between pentamidine derivatives and TMD-5, and to rationalise the observed structure-activity relationships. This study explicitly demonstrated that the interaction of small molecule with lipid should be considered alongside interaction with the protein target when designing small molecules targeting the PPIs of TMDs. In all, this study provides proof of concept for the rational design of small molecules targeting transmembrane PPIs.

In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides

Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.

Examination of pH dependency and orientation differences of membrane spanning alpha helices carrying a single or pair of buried histidine residues

We have employed the peptide framework of GWALP23 (acetyl-GGALWLALALALALALALWLAGA-amide) to examine the orientation, dynamics and pH dependence of peptides having buried single or pairs of histidine residues. When residue L8 is substituted to yield GWALP23-H⁸, acetyl-GGALWLAH⁸ALALALALALWLAGA-amide, the deuterium NMR spectra of ²H-labeled core alanine residues reveal a helix that occupies a single transmembrane orientation in DLPC, or in DMPC at low pH, yet shows multiple states at higher pH or in bilayers of DOPC. Moreover, a single histidine at position 8 or 16 in the GWALP23 framework is sensitive to pH. Titration points are observed near pH 3.5 for the deprotonation of H8 in lipid bilayers of DLPC or DMPC, and for H16 in DOPC. When residues L8 and L16 both are substituted to yield GWALP23-H^8,16, the ²H NMR spectra show, interestingly, no titration dependence from pH 2-8, yet bilayer thickness-dependent orientation differences. The helix with H8 and H16 is found to adopt a transmembrane orientation in thin bilayers of DLPC, a combination of transmembrane and surface orientations in DMPC, and then a complete transition to a surface bound orientation in the thicker DPoPC and DOPC lipid bilayers. In the surface orientations, alanine A7 no longer fits within the core helix. These results along with previous studies with different locations of histidine residues suggest that lipid hydrophobic thickness is a first determinant and pH a second determinant for the helical orientation, along with possible side-chain snorkeling, when the His residues are incorporated into the hydrophobic region of a lipid membrane-associated helix.

Structure and Dynamics of GPCRs in Lipid Membranes: Physical Principles and Experimental Approaches

Over the past decade, the vast amount of information generated through structural and biophysical studies of GPCRs has provided unprecedented mechanistic insight into the complex signalling behaviour of these receptors. With this recent information surge, it has also become increasingly apparent that in order to reproduce the various effects that lipids and membranes exert on the biological function for these allosteric receptors, in vitro studies of GPCRs need to be conducted under conditions that adequately approximate the native lipid bilayer environment. In the first part of this review, we assess some of the more general effects that a membrane environment exerts on lipid bilayer-embedded proteins such as GPCRs. This is then followed by the consideration of more specific effects, including stoichiometric interactions with specific lipid subtypes. In the final section, we survey a range of different membrane mimetics that are currently used for in vitro studies, with a focus on NMR applications.

Mechanism of protein cleavage at asparagine leading to protein-protein cross-links

Long-lived proteins (LLPs) are present in numerous tissues within the human body. With age, they deteriorate, often leading to the formation of irreversible modifications such as peptide bond cleavage and covalent cross-linking. Currently understanding of the mechanism of formation of these cross-links is limited. As part of an ongoing study, proteomics was used to characterise sites of novel covalent cross-linking in the human lens. In this process, Lys residues were found cross-linked to C-terminal aspartates that had been present in the original protein as Asn residues. Cross-links were identified in major lens proteins such as αA-crystallin, αB-crystallin and aquaporin 0. Quantification of the level of an AQP0/AQP0 cross-linked peptide showed increased cross-linking with age and in cataract lenses. Using model peptides, a mechanism of cross-link formation was elucidated that involves spontaneous peptide bond cleavage on the C-terminal side of Asn residues resulting in the formation of a C-terminal succinimide. This succinimide does not form cross-links, but can hydrolyse to a mixture of C-terminal Asn and C-terminal Asp amide peptides. The C-terminal Asp amide is unstable at neutral pH and decomposes to a succinic anhydride. If the side chain of Lys attacks the anhydride, a covalent cross-link will be formed. This multi-step mechanism represents a link between two spontaneous events: peptide bond cleavage at Asn and covalent cross-linking. Since Asn deamidation and cleavage are abundant age-related modifications in LLPs, this finding suggests that such susceptible Asn residues should also be considered as potential sites for spontaneous covalent cross-linking.

What Makes a Good Pore Former: A Study of Synthetic Melittin Derivatives

Pore formation by membrane-active peptides, naturally encountered in innate immunity and infection, could have important medical and technological applications. Recently, the well-studied lytic peptide melittin has formed the basis for the development of combinatorial libraries from which potent pore-forming peptides have been derived, optimized to work under different conditions. We investigate three such peptides, macrolittin70, which is most active at neutral pH; pHD15, which is active only at low pH; and MelP5_Δ6, which was rationally designed to be active at low pH but formed only small pores. There are three, six, and six acidic residues in macrolittin70, pHD15, and MelP5_Δ6, respectively. We perform multi-microsecond simulations in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) of hexamers of these peptides starting from transmembrane orientations at neutral pH (all residues at standard protonation), low pH (acidic residues and His protonated), and highly acidic environments in which C-termini are also protonated. Previous simulations of the parent peptides melittin and MelP5 in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) are repeated in POPC. We find that the most potent pore-forming peptides exhibit strong interpeptide interactions, including salt bridges, H-bonds, and polar interactions. Protonation of the C-terminus promotes helicity and pore size. The proximity of the peptides allows fewer lipid headgroups to line the pores than in previous simulations, making the pores intermediate between barrel stave and toroidal. Based on these structures and geometrical arguments, we attempt to rationalize the factors that under different conditions can increase or decrease pore stability and propose mutations that could be tested experimentally.

Comparing Interfacial Trp, Interfacial His and pH Dependence for the Anchoring of Tilted Transmembrane Helical Peptides

Charged and aromatic amino acid residues, being enriched toward the terminals of membrane-spanning helices in membrane proteins, help to stabilize particular transmembrane orientations. Among them, histidine is aromatic and can be positively charge at low pH. To enable investigations of the underlying protein-lipid interactions, we have examined the effects of single or pairs of interfacial histidine residues using the constructive low-dynamic GWALP23 (acetyl-GG²ALW⁵LALALALALALALW¹⁹LAG²²A-amide) peptide framework by incorporating individual or paired histidines at locations 2, 5, 19 or 22. Analysis of helix orientation by means of solid-state ²H NMR spectra of labeled alanine residues reveals marked differences with H^2,22 compared to W^2,22. Nevertheless, the properties of membrane-spanning H^2,22WALP23 helices show little pH dependence and are similar to those having Gly, Arg or Lys at positions 2 and 22. The presence of H5 or H19 influences the helix rotational preference but not the tilt magnitude. H5 affects the helical integrity, as residue 7 unwinds from the core helix; yet once again the helix orientation and dynamic properties show little sensitivity to pH. The overall results reveal that the detailed properties of transmembrane helices depend upon the precise locations of interfacial histidine residues.

Influence of interfacial tryptophan residues on an arginine-flanked transmembrane helix

The transmembrane helices of membrane proteins often are flanked by interfacial charged or aromatic residues that potentially help to anchor the membrane-spanning protein. For isolated single-span helices, the interfacial residues may be especially important for stabilizing particular tilted transmembrane orientations. The peptide RWALP23 (acetyl-GR²ALW(LA)₆LWLAR²²A-amide) has been employed to investigate the interplay between interfacial arginines and tryptophans. Here we replace the tryptophans of RWALP23 with A5 and A19, to investigate arginines alone with respect to helix fraying and orientation in varying lipid bilayers. Deuterated alanines incorporated into the central sequence allow the orientation and stability of the core helix to be assessed by means of solid -state ²H NMR in bilayers of DOPC, DMPC and DLPC. The helix tilt from the bilayer normal is found to increase slightly when R2 and R22 are present, and increases still further when the tryptophans W5 and W19 are replaced by alanines. The extent of helix dynamic averaging remains low in all cases. The preferred helix azimuthal rotation is essentially constant for all of the helices in each of the lipid membranes considered here. The alanines located outside of the core region of the peptide are sensitive to helical integrity. The new alanines, A5 and A19, therefore, provide new information about the length of the core helix and the onset of unraveling of the terminals. Residue A19 remains essentially on the central helix in each lipid membrane, while residues A3, A5 and A21 deviate from the core helix to an extent that depends on the membrane thickness. Differential unraveling of the two ends to expose peptide backbone groups for hydrogen bonding therefore acts together with specific interfacial side chains to stabilize a transmembrane helix.

NMR "Crystallography" for Uniformly (13 C, 15 N)-Labeled Oriented Membrane Proteins

In oriented-sample (OS) solid-state NMR of membrane proteins, the angular-dependent dipolar couplings and chemical shifts provide a direct input for structure calculations. However, so far only ¹ H-¹⁵ N dipolar couplings and ¹⁵ N chemical shifts have been routinely assessed in oriented ¹⁵ N-labeled samples. The main obstacle for extending this technique to membrane proteins of arbitrary topology has remained in the lack of additional experimental restraints. We have developed a new experimental triple-resonance NMR technique, which was applied to uniformly doubly (¹⁵ N, ¹³ C)-labeled Pf1 coat protein in magnetically aligned DMPC/DHPC bicelles. The previously inaccessible ¹ H_α -¹³ C_α dipolar couplings have been measured, which make it possible to determine the torsion angles between the peptide planes without assuming α-helical structure a priori. The fitting of three angular restraints per peptide plane and filtering by Rosetta scoring functions has yielded a consensus α-helical transmembrane structure for Pf1 protein.

Modelling of interactions between Aβ(25-35) peptide and phospholipid bilayers: effects of cholesterol and lipid saturation

Aggregation of amyloid beta (Aβ) peptides in neuronal membranes is a known promoter of Alzheimer’s disease. To gain insight into the molecular details of Aβ peptide aggregation and its effect on model neuronal membranes, we carried out molecular dynamics simulations of the Aβ(25–35) fragment of the amyloid precursor protein in phospholipid bilayers composed of either fully saturated or highly unsaturated lipids, in the presence or absence of cholesterol. It was found that the peptide does not penetrate through any of the considered membranes, but can reside in the headgroup region and upper part of the lipid tails showing a clear preference to a polyunsaturated cholesterol-free membrane. Due to the ordering and condensing effect upon addition of cholesterol, membranes become more rigid facilitating peptide aggregation on the surface. Except for the case of the cholesterol-free saturated lipid bilayer, the peptides have a small effect on the membrane structure and ordering. It was also found that the most “active” amino-acid for peptide–lipid and peptide–cholesterol interaction is methionine-35, followed by asparagine-27 and serine-26, which form hydrogen bonds between peptides and polar atoms of lipid headgroups. These amino acids are also primarily responsible for peptide aggregation. This work will be relevant for designing strategies to develop drugs to combat Alzheimer’s disease.

Molecular dynamics simulations of Aβ(25–35) peptides in phospholipid bilayers are carried out to investigate the effect of polyunsaturated lipids and cholesterol on aggregation of the peptides.

Side-Chain Protonation States of a Fluorescent Arginine

An arginine derivative with a fluorescent side-chain, Boc-Arg(Nap)-OH, was prepared by palladium(0)-catalyzed coupling of Boc-Arg-OH with a 4-bromonaphthalimide. The presence of the fluorophore lowers the pK_a of the side-chain guanidinium group by several orders of magnitude, to 9.0 (±0.1), allowing the derivative to access an electrically neutral protonation state that is not generally available to arginine itself. Computational modeling (DFT) predicts that protonation takes place at the side-chain C═N atom that bears the fluorophore. Calculated electronic absorptions for the protonated (356 nm) and neutral species (440 nm) are in good agreement with experiment. When irradiated with light, excited-state proton transfer (ESPT) occurs from cationic side-chains to suitably basic solvents, resulting in fluorescence emission from the neutral species. Arg(Nap) can be incorporated into peptides with sterically accessible N-termini using specially adapted conditions of solid-phase peptide synthesis.

Gating control of the cardiac sodium channel Nav1.5 by its β3-subunit involves distinct roles for a transmembrane glutamic acid and the extracellular domain

The auxiliary β3-subunit is an important functional regulator of the cardiac sodium channel Nav1.5, and some β3 mutations predispose individuals to cardiac arrhythmias. The β3-subunit uses its transmembrane α-helix and extracellular domain to bind to Nav1.5. Here, we investigated the role of an unusually located and highly conserved glutamic acid (Glu-176) within the β3 transmembrane region and its potential for functionally synergizing with the β3 extracellular domain (ECD). We substituted Glu-176 with lysine (E176K) in the WT β3-subunit and in a β3-subunit lacking the ECD. Patch-clamp experiments indicated that the E176K substitution does not affect the previously observed β3-dependent depolarizing shift of V_½ of steady-state inactivation but does attenuate the accelerated recovery from inactivation conferred by the WT β3-subunit. Removal of the β3-ECD abrogated both the depolarizing shift of steady-state inactivation and the accelerated recovery, irrespective of the presence or absence of the Glu-176 residue. We found that steady-state inactivation and recovery from inactivation involve movements of the S4 helices within the DIII and DIV voltage sensors in response to membrane potential changes. Voltage-clamp fluorometry revealed that the E176K substitution alters DIII voltage sensor dynamics without affecting DIV. In contrast, removal of the ECD significantly altered the dynamics of both DIII and DIV. These results imply distinct roles for the β3-Glu-176 residue and the β3-ECD in regulating the conformational changes of the voltage sensors that determine channel inactivation and recovery from inactivation.

Influence of Lipid Saturation, Hydrophobic Length and Cholesterol on Double-Arginine-Containing Helical Peptides in Bilayer Membranes

Membrane proteins are essential for many cell processes yet are more difficult to investigate than soluble proteins. Charged residues often contribute significantly to membrane protein function. Model peptides such as GWALP23 (acetyl-GGALW⁵ LAL⁸ LALALAL¹⁶ ALW¹⁹ LAGA-amide) can be used to characterize the influence of specific residues on transmembrane protein domains. We have substituted R8 and R16 in GWALP23 in place of L8 and L16, equidistant from the peptide center, and incorporated specific ² H-labeled alanine residues within the central sequence for detection by solid-state ² H NMR spectroscopy. The resulting pattern of [² H]Ala quadrupolar splitting (Δν_q ) magnitudes indicates the core helix for R^8,16 GWALP23 is significantly tilted to give a similar transmembrane orientation in thinner bilayers with either saturated C12:0 or C14:0 acyl chains (1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) or 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)) or unsaturated C16:1 Δ9 cis acyl chains. In bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC; C18:1 Δ9 cis) multiple orientations are indicated, whereas in longer, unsaturated 1,2-dieicosenoyl-sn-glycero-3-phosphocholine (DEiPC; C20:1 Δ11 cis) bilayers, the R^8,16 GWALP23 helix adopts primarily a surface orientation. The inclusion of 10-20 mol % cholesterol in DOPC bilayers drives more of the R^8,16 GWALP23 helix population to the membrane surface, thereby allowing both charged arginines access to the interfacial lipid head groups. The results suggest that hydrophobic thickness and cholesterol content are more important than lipid saturation for the arginine peptide dynamics and helix orientation in lipid membranes.

Breaking the Backbone: Central Arginine Residues Induce Membrane Exit and Helix Distortions within a Dynamic Membrane Peptide

Transmembrane domains of membrane proteins sometimes contain conserved charged or ionizable residues which may be essential for protein function and regulation. This work examines the molecular interactions of single Arg residues within a highly dynamic transmembrane peptide helix. To this end, we have modified the GW^4,20ALP23 (acetyl-GGAW⁴(AL)₇AW²⁰AGA-amide) model peptide framework to incorporate Arg residues near the center of the peptide. Peptide helix formation, orientation and dynamics were analyzed by means of solid-state NMR spectroscopy to monitor specific ²H- or ¹⁵N-labeled residues. GW^4,20ALP23 itself adopts a tilted orientation within lipid bilayer membranes. Nevertheless, the GW^4,20ALP23 helix exhibits moderate to high dynamic averaging of NMR observables, such as ²H quadrupolar splittings or ¹⁵N-¹H dipolar couplings, due to competition between the interfacial Trp residues on opposing helix faces. Here we examine how the helix dynamics are impacted by the introduction of a single Arg residue at position 12 or 14. Residue R14 restricts the helix to low dynamic averaging and a well-defined tilt that varies inversely with the lipid bilayer thickness. To compensate for the dominance of R14, the competing Trp residues cause partial unwinding of the helix at the C-terminal. By contrast, R¹²GW^4,20ALP23 exits the DOPC bilayer to an interfacial surface-bound location. Interestingly, multiple orientations are exhibited by a single residue, Ala-9. Quadrupolar splittings generated by ²H-labeled residues A3, A5, A7, and A9 do not fit to the α-helical quadrupolar wave plot defined by residues A11, A13, A15, A17, A19, and A21. The discontinuity at residue A9 implicates a helical swivel distortion and an apparent 3₁₀-helix involving the N-terminal residues preceding A11. These molecular features suggest that, while arginine residues are prominent factors controlling transmembrane helix dynamics, the influence of interfacial tryptophan residues cannot be ignored.

Universal principles of membrane protein assembly, composition and evolution

Structural diversity in α-helical membrane proteins (MP) arises from variations in helix-helix crossings and contacts that may bias amino acid usage. Here, we reveal systematic changes in transmembrane amino acid frequencies (f) as a function of the number of helices (n). For eukarya, breaks in f(n) trends of packing (Ala, Gly and Pro), polar, and hydrophobic residues identify different MP assembly principles for 2≤n≤7, 8≤n≤12 and n≥13. In bacteria, the first f break already occurs after n = 6 in correlation to an earlier n peak in MP size distribution and dominance of packing over polar interactions. In contrast to the later n brackets, the integration levels of helix bundles continuously increased in the first, most populous brackets indicating the formation of single structural units (domains). The larger first bracket of eukarya relates to a balance of polar and packing interactions that enlarges helix-helix combinatorial possibilities (MP diversity). Between the evolutionary old, packing and new, polar residues f anti-correlations extend over all biological taxa, broadly ordering them according to evolutionary history and allowing f estimates for the earliest forms of life. Next to evolutionary history, the amino acid composition of MP is determined by size (n), proteome diversity, and effective amino acid cost.

Validation of protein backbone structures calculated from NMR angular restraints using Rosetta

Multidimensional solid-state NMR spectra of oriented membrane proteins can be used to infer the backbone torsion angles and hence the overall protein fold by measuring dipolar couplings and chemical shift anisotropies, which depend on the orientation of each peptide plane with respect to the external magnetic field. However, multiple peptide plane orientations can be consistent with a given set of angular restraints. This ambiguity is further exacerbated by experimental uncertainty in obtaining and interpreting such restraints. The previously developed algorithms for structure calculations using angular restraints typically involve a sequential walkthrough along the backbone to find the torsion angles between the consecutive peptide plane orientations that are consistent with the experimental data. This method is sensitive to experimental uncertainty in interpreting the peak positions of as low as ± 10 Hz, often yielding high structural RMSDs for the calculated structures. Here we present a significantly improved version of the algorithm which includes the fitting of several peptide planes at once in order to prevent propagation of error along the backbone. In addition, a protocol has been devised for filtering the structural solutions using Rosetta scoring functions in order to find the structures that both fit the spectrum and satisfy bioinformatics restraints. The robustness of the new algorithm has been tested using synthetic angular restraints generated from the known structures for two proteins: a soluble protein 2gb1 (56 residues), chosen for its diverse secondary structure elements, i.e. an alpha-helix and two beta-sheets, and a membrane protein 4a2n, from which the first two transmembrane helices (having a total of 64 residues) have been used. Extensive simulations have been performed by varying the number of fitted planes, experimental error, and the number of NMR dimensions. It has been found that simultaneously fitting two peptide planes always shifted the distribution of the calculated structures toward lower structural RMSD values as compared to fitting a single torsion-angle pair. For each protein, irrespective of the simulation parameters, Rosetta was able to distinguish the most plausible structures, often having structural RMSDs lower than 2 Å with respect to the original structure. This study establishes a framework for de-novo protein structure prediction using a combination of solid-state NMR angular restraints and bioinformatics.

Functional Analysis of Conserved Transmembrane Charged Residues and a Yeast Specific Extracellular Loop of the Plasma Membrane Na+/H+ Antiporter of Schizosaccharomyces pombe

The Na⁺/H⁺ exchanger of the plasma membrane of S. pombe (SpNHE1) removes excess intracellular sodium in exchange for an extracellular proton. We examined the functional role of acidic amino acids of a yeast specific periplasmic extracellular loop 6 (EL6) and of Glu⁷⁴ and Arg⁷⁷ of transmembrane segment 3. Glu⁷⁴ and Arg⁷⁷ are conserved in yeast species while Glu⁷⁴ is conserved throughout various phyla. The mutation E74A caused a minor effect, while mutation R77A had a larger effect on the ability of SpNHE1 to confer salt tolerance. Mutation of both residues to Ala or Glu also eliminated the ability to confer salt tolerance. Arg³⁴¹ and Arg³⁴² were also necessary for SpNHE1 transport in S. pombe. Deletion of 3 out of 4 acidic residues (Asp³⁸⁹, Glu³⁹⁰, Glu³⁹², Glu³⁹⁷) of EL6 did not greatly affect SpNHE1 function while deletion of all did. Replacement of EL6 with a segment from the plant Na⁺/H⁺ exchanger SOS1 also did not affect function. We suggest that EL6 forms part of a cation coordination sphere, attracting cations for transport but that the region is not highly specific for the location of acidic charges. Overall, we identified a number of polar amino acids important in SpNHE1 function.

Molecular Dynamics Investigation into the Effect of Zinc(II) on the Structure and Membrane Interactions of the Antimicrobial Peptide Clavanin A

Clavanin A (ClavA) is an antimicrobial peptide (AMP) whose antimicrobial activity is enhanced in the presence of Zn(II) ions. The antimicrobial activity of ClavA has been shown to increase 16-fold in the presence of Zn(II) ions. In this study, we investigate the potential sources of this enhancement, namely, the effect of Zn(II) binding on the helical conformation of ClavA and on the ClavA interaction with a model for gram-negative bacterial membranes. In addition, we investigate the effect of Zn(II) on the membrane mechanical properties. We employed all-atom equilibrium molecular dynamics simulations initiated from both fully helical and random coil structures of ClavA. We observe that Zn(II) can stabilize an existing helical conformation in the Zn(II)-binding region, but we do not observe induction of helical conformations in systems initiated in random coil configurations. Zn(II) binding to ClavA provides more favorable electrostatics for membrane association in the C-terminal region. This is evidenced by longer and stronger C-terminal-lipid interactions. Zn(II) is also capable of modulating the membrane properties in a manner which favors ClavA insertion and the potential for enhanced translocation into the cell. This work provides insights into the role of divalent metal cations in the antimicrobial activity of ClavA. This information can be used for the development of synthetic AMPs containing motifs that can bind metals (metalloAMPs) for therapeutic and medical purposes.

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.

Actionable Activating Oncogenic ERBB2/HER2 Transmembrane and Juxtamembrane Domain Mutations

Deregulated HER2 is a target of many approved cancer drugs. We analyzed 111,176 patient tumors and identified recurrent mutations in HER2 transmembrane domain (TMD) and juxtamembrane domain (JMD) that include G660D, R678Q, E693K, and Q709L. Using a saturation mutagenesis screen and testing of patient-derived mutations we found several activating TMD and JMD mutations. Structural modeling and analysis showed that the TMD/JMD mutations function by improving the active dimer interface or stabilizing an activating conformation. Further, we found that HER2 G660D employed asymmetric kinase dimerization for activation and signaling. Importantly, anti-HER2 antibodies and small-molecule kinase inhibitors blocked the activity of TMD/JMD mutants. Consistent with this, a G660D germline mutant lung cancer patient showed remarkable clinical response to HER2 blockade.

In silico insights on diverse interacting partners and phosphorylation sites of respiratory burst oxidase homolog (Rbohs) gene families from Arabidopsis and rice

BACKGROUND

NADPH oxidase (Nox) is a critical enzyme involved in the generation of apoplastic superoxide (O₂^-), a type of reactive oxygen species (ROS) and hence regulate a wide range of biological functions in many organisms. Plant Noxes are the homologs of the catalytic subunit from mammalian NADPH oxidases and are known as respiratory burst oxidase homologs (Rbohs). Previous studies have highlighted their versatile roles in tackling different kind of stresses and in plant growth and development. In the current study, potential interacting partners and phosphorylation sites were predicted for Rboh proteins from two model species (10 Rbohs from Arabidopsis thaliana and 9 from Oryza sativa japonica). The present work is the first step towards in silico prediction of interacting partners and phosphorylation sites for Rboh proteins from two plant species.

RESULTS

In this work, an extensive range of potential partners (unique and common), leading to diverse functions were revealed from interaction networks and gene ontology classifications, where majority of AtRbohs and OsRbohs play role in stress-related activities, followed by cellular development. Further, 68 and 38 potential phosphorylation sites were identified in AtRbohs and OsRbohs, respectively. Their distribution, location and kinase specificities were also predicted and correlated with experimental data as well as verified with the other EF-hand containing proteins within both genomes.

CONCLUSIONS

Analysis of regulatory mechanisms including interaction with diverse partners and post-translational modifications like phosphorylation have provided insights regarding functional multiplicity of Rbohs. The bioinformatics-based workflow in the current study can be used to get insights for interacting partners and phosphorylation sites from Rbohs of other plant species.

Membrane-Induced p Ka Shifts in wt-pHLIP and Its L16H Variant

The pH (low) insertion peptides (pHLIPs) is a family of peptides that are able to insert into a lipid bilayer at acidic pH. The molecular mechanism of pHLIPs insertion, folding, and stability in the membrane at low pH is based on multiple protonation events, which are challenging to study at the molecular level. More specifically, the relation between the experimental p K of insertion (p K^exp) of pHLIPs and the p K_a of the key residues is yet to be clarified. We carried out a computational study, complemented with new experimental data, and established the influence of (de)protonation of titrable residues on the stability of the peptide membrane-inserted state. Constant-pH molecular dynamics simulations were employed to calculate the p K_a values of these residues along the membrane normal. In the wt-pHLIP, we identified Asp14 as the key residue for the stability of the membrane-inserted state, and its p K_a value is strongly correlated with the experimental p K^exp measured in thermodynamics studies. Also, in order to narrow down the pH range at which pHLIP is stable in the membrane, we designed a new pHLIP variant, L16H, where Leu in the 16th position was replaced by a titrable His residue. Our results showed that the L16H variant undergoes two transitions. The calculated p K_a and experimentally observed p K^exp values are in good agreement. Two distinct p K^exp values delimit a pH range where the L16H peptide is stably inserted in the membrane, while, outside this range, the membrane-inserted state is destabilized and the peptide exits from the bilayer. pHLIP peptides have been successfully used to target cancer cells for the delivery of diagnostics and therapeutic agents to acidic tumors. The fine-tuning of the stability of the pHLIP inserted state and its restriction to a narrow well-defined pH range might allow the design of new peptides, able to discriminate between tissues with different extracellular pH values.

Counting charges on membrane-bound peptides

Quantifying the number of charges on peptides bound to interfaces requires reliable estimates of (i) surface coverage and (ii) surface charge, both of which are notoriously difficult parameters to obtain, especially at solid/water interfaces. Here, we report the thermodynamics and electrostatics governing the interactions of l-lysine and l-arginine octamers (Lys₈ and Arg₈) with supported lipid bilayers prepared.

Quantifying the number of charges on peptides bound to interfaces requires reliable estimates of (i) surface coverage and (ii) surface charge, both of which are notoriously difficult parameters to obtain, especially at solid/water interfaces. Here, we report the thermodynamics and electrostatics governing the interactions of l-lysine and l-arginine octamers (Lys₈ and Arg₈) with supported lipid bilayers prepared from a 9 : 1 mixture of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (sodium salt) (DMPG) from second harmonic generation (SHG) spectroscopy, quartz crystal microbalance with dissipation monitoring (QCM-D) and nanoplasmonic sensing (NPS) mass measurements, and atomistic simulations. The combined SHG/QCM-D/NPS approach provides interfacial charge density estimates from mean field theory for the attached peptides that are smaller by a factor of approximately two (0.12 ± 0.03 C m^–2 for Lys₈ and 0.10 ± 0.02 C m^–2 for Arg₈) relative to poly-l-lysine and poly-l-arginine. These results, along with atomistic simulations, indicate that the surface charge density of the supported lipid bilayer is neutralized by the attached cationic peptides. Moreover, the number of charges associated with each attached peptide is commensurate with those found in solution; that is, Lys₈ and Arg₈ are fully ionized when attached to the bilayer. Computer simulations indicate Lys₈ is more likely than Arg₈ to “stand-up” on the surface, interacting with lipid headgroups through one or two sidechains while Arg₈ is more likely to assume a “buried” conformation, interacting with the bilayer through up to six sidechains. Analysis of electrostatic potential and charge distribution from atomistic simulations suggests that the Gouy–Chapman model, which is widely used for mapping surface potential to surface charge, is semi-quantitatively valid; despite considerable orientational preference of interfacial water, the apparent dielectric constant for the interfacial solvent is about 30, due to the thermal fluctuation of the lipid–water interface.

In situ examination of a charged amino acid-induced structural change in lipid bilayers by sum frequency generation vibrational spectroscopy

The interactions between amino acids (AAs) and membranes represent various short-range and long-range interactions for biological phenomena; however, they are still poorly understood. In this study, we used cationic lysine and arginine as AA models, and systematically investigated the interactions between charged AAs and lipid bilayers using sum frequency generation vibrational spectroscopy (SFG-VS) in situ and in real time. The AA-induced dynamic structural changes of the lipid bilayer were experimentally monitored using the spectral features of CD₂, CD₃, the lipid head phosphate, and carbonyl groups in real time. Time-dependent SFG changes in the structure of the lipid bilayer provide direct evidence for the different interactions of lysine and arginine with the membrane. It was found that the discrepancy between lysine and arginine in binding with the lipid bilayer is due to the nature of the terminal functional groups. Arginine exhibits a more drastic impact on the membrane than lysine. SFG responses of the acyl chains, phosphate groups, and carbonyl groups provide evidence that the interaction between AAs and the membrane most likely follows an electrostatics and hydrogen bond-induced defect model. This work presents an exemplary method for comprehensive investigations of interactions between membranes and other functionally significant substances.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Refining amino acid hydrophobicity for dynamics simulation of membrane proteins

Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.

Effects of hydration on the protonation state of a lysine analog crossing a phospholipid bilayer – insights from molecular dynamics and free-energy calculations

Protonation states of amino acids crossing lipid bilayers from multidimensional free energy surfaces.

Fluorescamine Labeling for Assessment of Protein Conformational Change and Binding Affinity in Protein-Nanoparticle Interaction

Protein adsorption alters the "biological identity" of nanoparticles (NPs) and could affect how biosystems respond to invading NPs. Study of protein-NP interaction can help understand how the physicochemical properties of NPs impact the interaction and thus potentially guide the design of safer and more effective NPs for biomedical or other applications. Binding affinity between proteins and NPs and the occurrence of protein conformational change upon binding to NPs are two important aspects to be learned, but few methods are currently available to assess both simultaneously in a simple way. Herein, we demonstrated that the fluorescamine labeling method developed by our group not only could reveal protein conformational change upon adsorption to NPs, owing to its capability to label the primary amines exposed on protein surface, but also could be applied to measure the binding affinity. By screening the interaction between a large number of proteins and four types of NPs, the present study also revealed that protein adsorption onto NPs could be strongly affected by structure flexibility. The proteins with high structure flexibility experienced high degrees of conformation change when binding to the polystyrene NPs, which could potentially influence protein function. Overall, we demonstrate that our assay is a quick, simple, and high-throughput tool to reveal potential impacts on protein activity and evaluate the strength of protein-NP binding.

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.

Anhydrous Monoalkylguanidines in Aprotic and Nonpolar Solvents: Models for Deprotonated Arginine Side Chains in Membrane Environments

In this study, the synthesis of crystalline dodecylguanidine free base and its spectroscopic characterization in nonpolar environments are described. IR as well as ¹H and ¹⁵N NMR spectra of the free base dissolved in aprotic solvents are substantially different from the previously reported spectra of arginine, or other monoalkylguanidinium compounds, at high hydroxide concentrations. The current results provide improved modeling for the spectroscopic signals that would be expected from a deprotonated arginine in a nonpolar environment. On the basis of our spectra of the authentic dodecylguanidine free base, addition of large amounts of aqueous hydroxide to arginine or other monoalklyguanidinium salts does not deprotonate them. Instead, hydroxide addition leads to the formation of a guanidinium hydroxide complex, with a dissociation constant near ∼500 mM that accounts for the established arginine pK value of ∼13.7. We also report a method for synthesizing a compound containing both phenol and free-base guanidine groups, linked by a dodecyl chain that should be generalizable to other hydrocarbon linkers. Such alkyl-guanidine and phenolyl-alkyl-guanidine compounds can serve as small-molecule models for the conserved arginine-tyrosine groupings that have been observed in crystallographic structures of both microbial rhodopsins and G-protein-coupled receptors.

A membrane-inserted structural model of the yeast mitofusin Fzo1

Mitofusins are large transmembrane GTPases of the dynamin-related protein family, and are required for the tethering and fusion of mitochondrial outer membranes. Their full-length structures remain unknown, which is a limiting factor in the study of outer membrane fusion. We investigated the structure and dynamics of the yeast mitofusin Fzo1 through a hybrid computational and experimental approach, combining molecular modelling and all-atom molecular dynamics simulations in a lipid bilayer with site-directed mutagenesis and in vivo functional assays. The predicted architecture of Fzo1 improves upon the current domain annotation, with a precise description of the helical spans linked by flexible hinges, which are likely of functional significance. In vivo site-directed mutagenesis validates salient aspects of this model, notably, the long-distance contacts and residues participating in hinges. GDP is predicted to interact with Fzo1 through the G1 and G4 motifs of the GTPase domain. The model reveals structural determinants critical for protein function, including regions that may be involved in GTPase domain-dependent rearrangements.

Membrane Active Antimicrobial Peptides: Translating Mechanistic Insights to Design

Antimicrobial peptides (AMPs) are promising next generation antibiotics that hold great potential for combating bacterial resistance. AMPs can be both bacteriostatic and bactericidal, induce rapid killing and display a lower propensity to develop resistance than do conventional antibiotics. Despite significant progress in the past 30 years, no peptide antibiotic has reached the clinic yet. Poor understanding of the action mechanisms and lack of rational design principles have been the two major obstacles that have slowed progress. Technological developments are now enabling multidisciplinary approaches including molecular dynamics simulations combined with biophysics and microbiology toward providing valuable insights into the interactions of AMPs with membranes at atomic level. This has led to increasingly robust models of the mechanisms of action of AMPs and has begun to contribute meaningfully toward the discovery of new AMPs. This review discusses the detailed action mechanisms that have been put forward, with detailed atomistic insights into how the AMPs interact with bacterial membranes. The review further discusses how this knowledge is exploited toward developing design principles for novel AMPs. Finally, the current status, associated challenges, and future directions for the development of AMP therapeutics are discussed.

Sensitivity enhancement for membrane proteins reconstituted in parallel and perpendicular oriented bicelles obtained by using repetitive cross-polarization and membrane-incorporated free radicals

Multidimensional separated local-field and spin-exchange experiments employed by oriented-sample solid-state NMR are essential for structure determination and spectroscopic assignment of membrane proteins reconstituted in macroscopically aligned lipid bilayers. However, these experiments typically require a large number of scans in order to establish interspin correlations. Here we have shown that a combination of optimized repetitive cross polarization (REP-CP) and membrane-embedded free radicals allows one to enhance the signal-to-noise ratio by factors 2.4-3.0 in the case of Pf1 coat protein reconstituted in magnetically aligned bicelles with their normals being either parallel or perpendicular to the main magnetic field. Notably, spectral resolution is not affected at the 2:1 radical-to-protein ratio. Spectroscopic assignment of Pf1 coat protein in the parallel bicelles has been established as an illustration of the method. The proposed methodology will advance applications of oriented-sample NMR technique when applied to samples containing smaller quantities of proteins and three-dimensional experiments.

Influence of glutamic acid residues and pH on the properties of transmembrane helices

Negatively charged side chains are important for the function of particular ion channels and certain other membrane proteins. To investigate the influence of single glutamic acid side chains on helices that span lipid-bilayer membranes, we have employed GWALP23 (acetyl-GGALW⁵LALALALALALALW¹⁹LAGA-amide) as a favorable host peptide framework. We substituted individual Leu residues with Glu residues (L12E or L14E or L16E) and incorporated specific ²H-labeled alanine residues within the core helical region or near the ends of the sequence. Solid-state ²H NMR spectra reveal little change for the core labels in GWALP23-E12, -E14 and -E16 over a pH range of 4 to 12.5, with the spectra being broader for samples in DOPC compared to DLPC bilayers. The spectra for samples with deuterium labels near the helix ends on alanines 3 and 21 show modest pH-dependent changes in the extent of unwinding of the helix terminals in DLPC and DOPC bilayers. The combined results indicate minor overall responses of these transmembrane helices to changes in pH, with the most buried residue E12 showing no pH dependence. While the Glu residues E14 and E16 may have high pK_a values in the lipid bilayer environment, it is also possible that a paucity of helix response is masking the pK_a values. Interestingly, when E16 is present, spectral changes at high pH report significant local unwinding of the core helix. Our results are consistent with the expectation that buried carboxyl groups aggressively hold their protons and/or waters of hydration.

Influence of High pH and Cholesterol on Single Arginine-Containing Transmembrane Peptide Helices

An essential component of mammalian cells, cholesterol exerts significant influence on the physical properties of the cell membrane and in turn its constituents, including membrane proteins. Although sparse, polar amino acid residues are highly conserved in membrane proteins and play pivotal roles in determining specific structural and functional properties. To improve our understanding of particular polar residues in the membrane environment, we have examined two specific “guest” Arg residues within a well-defined and deuterium-labeled “host” framework provided by the transmembrane helical peptide GWALP23 (acetyl-GGALWLALALALALALALWLAGA-amide). Solid-state ²H nuclear magnetic resonance (NMR) spectra from aligned bilayer membrane samples effectively report changes in the host helix properties because of the incorporation of the guest residues. The focus of this work is two-pronged. First, GWALP23-R14 was examined over a pH range of 2–13 in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) ester- or ether-linked bilayer membranes. Our findings indicate that the Arg guanidinium side chain remains charged over this entire range, in agreement with numerous molecular dynamics simulations. Second, GWALP23-R12 and GWALP23-R14 peptides were characterized in DOPC bilayers with varying cholesterol content. Our findings suggest that 10 or 20% cholesterol content has minimal impact on the orientation of the R14 peptide. Although the NMR signals are broader and weaker in the presence of 20% cholesterol, the deuterium quadrupolar splittings for [²H]Ala residues in GWALP23-R14 change very little. Conversely, cholesterol appears to modulate the multistate behavior of GWALP23-R12 and to favor a major interfacial state for the helix, bound at the bilayer surface. These results indicate a conditional sensitivity of a complex multistate transmembrane Arg-containing peptide helix to the presence of cholesterol.

Specific Binding of Cholesterol to C99 Domain of Amyloid Precursor Protein Depends Critically on Charge State of Protein

Recent NMR chemical shift measurements of the 99 residue C-terminal fragment of amyloid precursor protein (APP-C99) in the presence of cholesterol provide evidence of binary complex formation between C99 and cholesterol in membrane mimetic environments. It has also been observed that the production of Aβ protein is enhanced under conditions of high cholesterol concentration. In this study, we investigated the impact of the charge state of C99 on the structure and stability of the C99-cholesterol complex. We observed that the binding of C99 to cholesterol depends critically on the charge state of Glu 693 (E22) and Asp 694 (D23). Evaluation of the pKa values of the Asp and Glu side chains suggests that these residues may be predominantly neutral in existing experimental observations of a stable C99-cholesterol complex at lower pH (characteristic of the endosomal environment), while binding is destabilized near neutral pH (characteristic of the cytoplasm). These observations suggest that specific binding of cholesterol to C99 is a sensitive function of the pH encountered in vivo, with key E22 and D23 residues serving as a "pH switch" controlling C99-cholesterol binding.