Free energy determinants of peptide association with lipid bilayers

Gibberellin 2-Oxidases in Potato (Solanum tuberosum L.): Cloning, Characterization, In Silico Analysis and Molecular Docking

Gibberellins (GAs; tetracyclic di-terpenoid carboxylic acids) are endogenous plant growth regulators responsible for stimulating plant growth and development from seed germination to plant maturity. In potato (Solanum tuberosum L.), GA levels are known to be crucial in the complex process of tuberization. Gibberellin 2-oxidases (GA2oxs) inactivate bioactive GAs during stolon swelling and early stages of tuberization as evident from the predominant expression of a member of this gene family namely GA2ox1. We isolated and characterized a 1105-bp cDNA clone encoding a 340-aa GA2ox1 form, designated St-GA2ox1, using total RNA from growing tuber of a potato (Solanum tuberosum L.) cultivar, Kufri Chipsona-1 (KC-1) based on RT-PCR approach. A total of 26 GA2ox sequences were also retrieved from potato genome database and analysed. Multiple sequence alignment revealed sequence relatedness between the GA2oxs. Crucial protein motifs were identified. Phylogenetic analysis revealed the evolutionary relationships between the GA2oxs. Three-dimensional structure of St-GA2ox1 was predicted by using AlphaFold tool, validated by the predicted local-distance difference test and Ramachandran Plot. Structural analysis and molecular docking were carried out to identify domains, binding sites and affinity for the ligand. The STRING database and hydropathy analysis revealed the presence of a putative interaction site for other enzymes. Expression Atlas database and semi-quantitative RT-PCR revealed the expression patterns of various GA2ox forms in different potato organs. This comprehensive report would be useful in providing new insights into possible underlying mechanisms involved in tuber development, and could facilitate the targeted alteration of genes responsible to combat the stress and enhance tuber production.

Slow release of a synthetic auxin induces formation of adventitious roots in recalcitrant woody plants

Clonal propagation of plants by induction of adventitious roots (ARs) from stem cuttings is a requisite step in breeding programs. A major barrier exists for propagating valuable plants that naturally have low capacity to form ARs. Due to the central role of auxin in organogenesis, indole-3-butyric acid is often used as part of commercial rooting mixtures, yet many recalcitrant plants do not form ARs in response to this treatment. Here we describe the synthesis and screening of a focused library of synthetic auxin conjugates in Eucalyptus grandis cuttings and identify 4-chlorophenoxyacetic acid-L-tryptophan-OMe as a competent enhancer of adventitious rooting in a number of recalcitrant woody plants, including apple and argan. Comprehensive metabolic and functional analyses reveal that this activity is engendered by prolonged auxin signaling due to initial fast uptake and slow release and clearance of the free auxin 4-chlorophenoxyacetic acid. This work highlights the utility of a slow-release strategy for bioactive compounds for more effective plant growth regulation.

Computational identification of a multi-peptide vaccine candidate in E2 glycoprotein against diverse Hepatitis C virus genotypes

Hepatitis C Virus (HCV) is estimated to affect nearly 180 million people worldwide, culminating in ∼0.7 million yearly casualties. However, a safe vaccine against HCV is not yet available. This study endeavored to identify a multi-genotypic, multi-epitopic, safe, and globally competent HCV vaccine candidate. We employed a consensus epitope prediction strategy to identify multi-epitopic peptides in all known envelope glycoprotein (E2) sequences, belonging to diverse HCV genotypes. The obtained peptides were screened for toxicity, allergenicity, autoimmunity and antigenicity, resulting in two favorable peptides viz., P2 (VYCFTPSPVVVG) and P3 (YRLWHYPCTV). Evolutionary conservation analysis indicated that P2 and P3 are highly conserved, supporting their use as part of a designed multi-genotypic vaccine. Population coverage analysis revealed that P2 and P3 are likely to be presented by >89% Human Leukocyte Antigen (HLA) molecules from six geographical regions. Indeed, molecular docking predicted the physical binding of P2 and P3 to various representative HLAs. We designed a vaccine construct using these peptides and assessed its binding to toll-like receptor 4 (TLR-4) by molecular docking and simulation. Subsequent analysis by energy-based and machine learning tools predicted high binding affinity and pinpointed the key binding residues (i.e. hotspots) in P2 and P3. Also, a favorable immunogenic profile of the construct was predicted by immune simulations. We encourage the scientific community to validate our vaccine construct in vitro and in vivo.Communicated by Ramaswamy H. Sarma.

Molecular, in silico and expression studies on lipoxygenases (LOXs) in potato (Solanum tuberosum L.)

Lipoxygenases (LOXs) namely 9-LOXs and 13-LOXs catalyse the oxygenation of polyunsaturated fatty acids to produce fatty acid hydroperoxides which are crucial in growth, development and stress responses in plants. Here, we isolated and characterized a 2723-bp cDNA encoding a distinct 861-aa 9-LOX form, designated StKCLX-1, using tuber total RNA from an Indian potato cultivar, Kufri Chipsona-1 through RT-PCR. A total of 17 LOX genes distributed in different chromosomes were identified and characterized in the potato genome. Multiple sequence alignment revealed highly conserved amino acids in the crucial domains, motifs and variable N-terminal regions between the LOX classes. A total of 36 LOXs from potato, tomato and Arabidopsis were used in phylogenetic analysis. A 3-D structure of StKCLX-1 was predicted by AlphaFold tool, validated through the predicted local-distance difference test (pLDDT) and Ramachandran Plot. Molecular docking predicted the nature of receptor-ligand interactions. STRING database was used to predict the protein-protein interactions. Expression patterns of the LOXs in the potato organs were examined by Expression Atlas and semi-quantitative RT-PCR. 9-LOX activity was noticed at early stages of tuberization, and significantly increased in the freshly-harvested mature tubers. This report would be useful in gaining insights into the structure-function relationships of the LOXs and corresponding multigene family-prerequisites for understanding tuber development in potato.

The mitochondria-targeted peptide SS-31 binds lipid bilayers and modulates surface electrostatics as a key component of its mechanism of action

Mitochondrial dysfunction underlies many heritable diseases, acquired pathologies, and aging-related declines in health. Szeto-Schiller (SS) peptides comprise a class of amphipathic tetrapeptides that are efficacious toward a wide array of mitochondrial disorders and are believed to target mitochondrial membranes because they are enriched in the anionic phospholipid cardiolipin (CL). However, little is known regarding how SS peptides interact with or alter the physical properties of lipid bilayers. In this study, using biophysical and computational approaches, we have analyzed the interactions of the lead compound SS-31 (elamipretide) with model and mitochondrial membranes. Our results show that this polybasic peptide partitions into the membrane interfacial region with an affinity and a lipid binding density that are directly related to surface charge. We found that SS-31 binding does not destabilize lamellar bilayers even at the highest binding concentrations; however, it did cause saturable alterations in lipid packing. Most notably, SS-31 modulated the surface electrostatics of both model and mitochondrial membranes. We propose nonexclusive mechanisms by which the tuning of surface charge could underpin the mitoprotective properties of SS-31, including alteration of the distribution of ions and basic proteins at the interface, and/or modulation of bilayer physical properties. As a proof of concept, we show that SS-31 alters divalent cation (calcium) distribution within the interfacial region and reduces the energetic burden of calcium stress in mitochondria. The mechanistic details of SS-31 revealed in this study will help inform the development of future compound variants with enhanced efficacy and bioavailability.

New Insights into the Structure-Function Relationship of the Endosomal-Type Na+, K+/H+ Antiporter NHX6 from Mulberry (Morus notabilis)

The endosomal-type Na⁺, K⁺/H⁺ antiporters (NHXs) play important roles in K⁺, vesicle pH homeostasis, and protein trafficking in plant. However, the structure governing ion transport mechanism and the key residues related to the structure–function of the endosomal-type NHXs remain unclear. Here, the structure-function relationship of the only endosomal-type NHX from mulberry, MnNHX6, was investigated by homology modeling, mutagenesis, and localization analyses in yeast. The ectopic expression of MnNHX6 in arabidopsis and Nhx1 mutant yeast can enhance their salt tolerance. MnNHX6’s three-dimensional structure, established by homology modeling, was supported by empirical, phylogenetic, and experimental data. Structure analysis showed that MnNHX6 contains unusual 13 transmembrane helices, but the structural core formed by TM5-TM12 assembly is conserved. Localization analysis showed that MnNHX6 has the same endosomal localization as yeast Nhx1/VPS44, and Arg402 is important for protein stability of MnNHX6. Mutagenesis analysis demonstrated MnNHX6 contains a conserved cation binding mechanism and a similar charge-compensated pattern as NHE1, but shares a different role in ion selectivity than the vacuolar-type NHXs. These results improve our understanding of the role played by the structure–function related key residues of the plant endosomal-type NHXs, and provide a basis for the ion transport mechanism study of endosomal-type NHXs.

Interplay between hydrophobicity and the positive-inside rule in determining membrane-protein topology

The energetics of membrane-protein interactions determine protein topology and structure: hydrophobicity drives the insertion of helical segments into the membrane, and positive charges orient the protein with respect to the membrane plane according to the positive-inside rule. Until recently, however, quantifying these contributions met with difficulty, precluding systematic analysis of the energetic basis for membrane-protein topology. We recently developed the dsTβL method, which uses deep sequencing and in vitro selection of segments inserted into the bacterial plasma membrane to infer insertion-energy profiles for each amino acid residue across the membrane, and quantified the insertion contribution from hydrophobicity and the positive-inside rule. Here, we present a topology-prediction algorithm called TopGraph, which is based on a sequence search for minimum dsTβL insertion energy. Whereas the average insertion energy assigned by previous experimental scales was positive (unfavorable), the average assigned by TopGraph in a nonredundant set is -6.9 kcal/mol. By quantifying contributions from both hydrophobicity and the positive-inside rule we further find that in about half of large membrane proteins polar segments are inserted into the membrane to position more positive charges in the cytoplasm, suggesting an interplay between these two energy contributions. Because membrane-embedded polar residues are crucial for substrate binding and conformational change, the results implicate the positive-inside rule in determining the architectures of membrane-protein functional sites. This insight may aid structure prediction, engineering, and design of membrane proteins. TopGraph is available online (topgraph.weizmann.ac.il).

Critical Evaluation and Compilation of Physicochemical Determinants and Membrane Interactions of MMGP1 Antifungal Peptide

A growing issue of pathogen resistance to antibiotics has fostered the development of innovative approaches for novel drug development. Here, we report the physicochemical and biological properties of an antifungal peptide, MMGP1, based on computational analysis. Computation of physicochemical properties has revealed that the natural biological activities of MMGP1 are coordinated by its intrinsic properties such as net positive charge (+5.04), amphipathicity, high hydrophobicity, low hydrophobic moment, and higher isoelectric point (11.915). Prediction of aggregation hot spots in MMGP1 had revealed the presence of potentially aggregation-prone segments that can nucleate in vivo aggregation (on the membrane), whereas no aggregating regions were predicted for in vitro aggregation (in solutions) of MMGP1. This ability of MMGP1 to form oligomeric aggregates on membrane further substantiates its direct-cell penetrating potency. Monte Carlo simulation of the interactions of MMGP1 in the aqueous phase and different membrane environments revealed that increasing the proportion of acidic lipids on membrane had led to increase in the peptide helicity. Furthermore, the peptide adopts energetically favorable transmembrane configuration, by inserting peptide loop and helix termini into the membrane containing >60% of anionic lipids. The charged lipid-based insertion of MMGP1 into membrane might be responsible for the selectivity of peptide toward fungal cells. Additionally, MMGP1 possessed DNA-binding property. Computational docking has identified DNA-binding residues (TRP3, SER4, MET7, ARG8, PHE10, ALA11, GLY20, THR21, ARG22, MET23, TRP34, and LYS36) in MMGP1 crucial for its DNA-binding property. Furthermore, computational mutation analysis revealed that aromatic amino acids are crucial for in vivo aggregation, membrane insertion, and DNA-binding property of MMGP1. These data provide new insight into the molecular determinants of MMGP1 antifungal activity and also serves as the template for the design of novel peptide antibiotics.

Mutational scanning reveals the determinants of protein insertion and association energetics in the plasma membrane

Insertion of helix-forming segments into the membrane and their association determines the structure, function, and expression levels of all plasma membrane proteins. However, systematic and reliable quantification of membrane-protein energetics has been challenging. We developed a deep mutational scanning method to monitor the effects of hundreds of point mutations on helix insertion and self-association within the bacterial inner membrane. The assay quantifies insertion energetics for all natural amino acids at 27 positions across the membrane, revealing that the hydrophobicity of biological membranes is significantly higher than appreciated. We further quantitate the contributions to membrane-protein insertion from positively charged residues at the cytoplasm-membrane interface and reveal large and unanticipated differences among these residues. Finally, we derive comprehensive mutational landscapes in the membrane domains of Glycophorin A and the ErbB2 oncogene, and find that insertion and self-association are strongly coupled in receptor homodimers.

The role and significance of potential lipid-binding regions in the mitochondrial protein import motor: an in-depth in silico study

Over the last two decades, an impressive progress has been made in the identification of novel factors in the translocation machineries of the mitochondrial protein import and their possible roles. The role of lipids and possible protein–lipids interactions remains a relatively unexplored territory. Investigating the role of potential lipid-binding regions in the sub-units of the mitochondrial motor might help to shed some more light in our understanding of protein–lipid interactions mechanistically. Bioinformatics results seem to indicate multiple potential lipid-binding regions in each of the sub-units. The subsequent characterization of some of those regions in silico provides insight into the mechanistic functioning of this intriguing and essential part of the protein translocation machinery. Details about the way the regions interact with phospholipids were found by the use of Monte Carlo simulations. For example, Pam18 contains one possible transmembrane region and two tilted surface bound conformations upon interaction with phospholipids. The results demonstrate that the presented bioinformatics approach might be useful in an attempt to expand the knowledge of the possible role of protein–lipid interactions in the mitochondrial protein translocation process.

Structure, dynamics and implied gating mechanism of a human cyclic nucleotide-gated channel

Cyclic nucleotide-gated (CNG) ion channels are nonselective cation channels, essential for visual and olfactory sensory transduction. Although the channels include voltage-sensor domains (VSDs), their conductance is thought to be independent of the membrane potential, and their gating regulated by cytosolic cyclic nucleotide-binding domains. Mutations in these channels result in severe, degenerative retinal diseases, which remain untreatable. The lack of structural information on CNG channels has prevented mechanistic understanding of disease-causing mutations, precluded structure-based drug design, and hampered in silico investigation of the gating mechanism. To address this, we built a 3D model of the cone tetrameric CNG channel, based on homology to two distinct templates with known structures: the transmembrane (TM) domain of a bacterial channel, and the cyclic nucleotide-binding domain of the mouse HCN2 channel. Since the TM-domain template had low sequence-similarity to the TM domains of the CNG channels, and to reconcile conflicts between the two templates, we developed a novel, hybrid approach, combining homology modeling with evolutionary coupling constraints. Next, we used elastic network analysis of the model structure to investigate global motions of the channel and to elucidate its gating mechanism. We found the following: (i) In the main mode of motion, the TM and cytosolic domains counter-rotated around the membrane normal. We related this motion to gating, a proposition that is supported by previous experimental data, and by comparison to the known gating mechanism of the bacterial KirBac channel. (ii) The VSDs could facilitate gating (supplementing the pore gate), explaining their presence in such 'voltage-insensitive' channels. (iii) Our elastic network model analysis of the CNGA3 channel supports a modular model of allosteric gating, according to which protein domains are quasi-independent: they can move independently, but are coupled to each other allosterically.

Computational and experimental characterization of dVHL establish a Drosophila model of VHL syndrome

The von Hippel-Lindau (VHL) cancer syndrome is associated with mutations in the VHL gene. The pVHL protein is involved in response to changes in oxygen availability as part of an E3-ligase that targets the Hypoxia-Inducible Factor for degradation. pVHL has a molten globule configuration with marginal thermodynamic stability. The cancer-associated mutations further destabilize it. The Drosophila homolog, dVHL, has relatively low sequence similarity to pVHL, and is also involved in regulating HIF1-α. Using in silico, in vitro and in vivo approaches we demonstrate high similarity between the structure and function of dVHL and pVHL. These proteins have a similar fold, secondary and tertiary structures, as well as thermodynamic stability. Key functional residues in dVHL are evolutionary conserved. This structural homology underlies functional similarity of both proteins, evident by their ability to bind their reciprocal partner proteins, and by the observation that transgenic pVHL can fully maintain normal dVHL-HIF1-α downstream pathways in flies. This novel transgenic Drosophila model is thus useful for studying the VHL syndrome, and for testing drug candidates to treat it.

Functional evaluation of autism-associated mutations in NHE9

NHE9 (SLC9A9) is an endosomal cation/proton antiporter with orthologues in yeast and bacteria. Rare, missense substitutions in NHE9 are genetically linked with autism but have not been functionally evaluated. Here we use evolutionary conservation analysis to build a model structure of NHE9 based on the crystal structure of bacterial NhaA and use it to screen autism-associated variants in the human population first by phenotype complementation in yeast, followed by functional analysis in primary cortical astrocytes from mouse. NHE9-GFP localizes to recycling endosomes, where it significantly alkalinizes luminal pH, elevates uptake of transferrin and the neurotransmitter glutamate, and stabilizes surface expression of transferrin receptor and GLAST transporter. In contrast, autism-associated variants L236S, S438P and V176I lack function in astrocytes. Thus, we establish a neurobiological cell model of a candidate gene in autism. Loss-of-function mutations in NHE9 may contribute to autistic phenotype by modulating synaptic membrane protein expression and neurotransmitter clearance.

Monte Carlo simulations of peptide-membrane interactions with the MCPep web server

The MCPep server (http://bental.tau.ac.il/MCPep/) is designed for non-experts wishing to perform Monte Carlo (MC) simulations of helical peptides in association with lipid membranes. MCPep is a web implementation of a previously developed MC simulation model. The model has been tested on a variety of peptides and protein fragments. The simulations successfully reproduced available empirical data and provided new molecular insights, such as the preferred locations of peptides in the membrane and the contribution of individual amino acids to membrane association. MCPep simulates the peptide in the aqueous phase and membrane environments, both described implicitly. In the former, the peptide is subjected solely to internal conformational changes, and in the latter, each MC cycle includes additional external rigid body rotational and translational motions to allow the peptide to change its location in the membrane. The server can explore the interaction of helical peptides of any amino-acid composition with membranes of various lipid compositions. Given the peptide’s sequence or structure and the natural width and surface charge of the membrane, MCPep reports the main determinants of peptide–membrane interactions, e.g. average location and orientation in the membrane, free energy of membrane association and the peptide’s helical content. Snapshots of example simulations are also provided.

Protein stability: a single recorded mutation aids in predicting the effects of other mutations in the same amino acid site

Motivation: Accurate prediction of protein stability is important for understanding the molecular underpinnings of diseases and for the design of new proteins. We introduce a novel approach for the prediction of changes in protein stability that arise from a single-site amino acid substitution; the approach uses available data on mutations occurring in the same position and in other positions. Our algorithm, named Pro-Maya (Protein Mutant stAbilitY Analyzer), combines a collaborative filtering baseline model, Random Forests regression and a diverse set of features. Pro-Maya predicts the stability free energy difference of mutant versus wild type, denoted as ΔΔG.

Results: We evaluated our algorithm extensively using cross-validation on two previously utilized datasets of single amino acid mutations and a (third) validation set. The results indicate that using known ΔΔG values of mutations at the query position improves the accuracy of ΔΔG predictions for other mutations in that position. The accuracy of our predictions in such cases significantly surpasses that of similar methods, achieving, e.g. a Pearson's correlation coefficient of 0.79 and a root mean square error of 0.96 on the validation set. Because Pro-Maya uses a diverse set of features, including predictions using two other methods, it also performs slightly better than other methods in the absence of additional experimental data on the query positions.

Availability: Pro-Maya is freely available via web server at http://bental.tau.ac.il/ProMaya.

Contact:nirb@tauex.tau.ac.il; wolf@cs.tau.ac.il

Supplementary Information:Supplementary data are available at Bioinformatics online.

Model-guided mutagenesis drives functional studies of human NHA2, implicated in hypertension

Human NHA2 is a poorly characterized Na(+)/H(+) antiporter recently implicated in essential hypertension. We used a range of computational tools and evolutionary conservation analysis to build and validate a three-dimensional model of NHA2 based on the crystal structure of a distantly related bacterial transporter, NhaA. The model guided mutagenic evaluation of transport function, ion selectivity, and pH dependence of NHA2 by phenotype screening in yeast. We describe a cluster of essential, highly conserved titratable residues located in an assembly region made of two discontinuous helices of inverted topology, each interrupted by an extended chain. Whereas in NhaA, oppositely charged residues compensate for partial dipoles generated within this assembly, in NHA2, polar but uncharged residues suffice. Our findings led to a model for transport mechanism that was compared to the well-known electroneutral NHE1 and electrogenic NhaA subtypes. This study establishes NHA2 as a prototype for the poorly understood, yet ubiquitous, CPA2 antiporter family recently recognized in plants and metazoans and illustrates a structure-driven approach to derive functional information on a newly discovered transporter.

Energetics of hydrophobic matching in lipid-protein interactions

Lipid chain length modulates the activity of transmembrane proteins by mismatch between the hydrophobic span of the protein and that of the lipid membrane. Relative binding affinities of lipids with different chain lengths are used to estimate the excess free energy of lipid-protein interaction that arises from hydrophobic mismatch. For a wide range of integral proteins and peptides, the energy cost is much less than the elastic penalty of fully stretching or compressing the lipid chains to achieve complete hydrophobic matching. The chain length dependences of the free energies of lipid association are described by a model that combines elastic chain extension with a free energy term that depends linearly on the extent of residual mismatch. The excess free energy densities involved lie in the region of 0.5-2.0 k(B)T x nm(-2). Values of this size could arise from exposure of hydrophobic groups to polar portions of the lipid or protein, but not directly to water, or alternatively from changes in tilt of the transmembrane helices that are energetically comparable to those activating mechanosensitive channels. The influence of hydrophobic mismatch on dimerization of transmembrane helices and their transfer between lipid vesicles, and on shifts in chain-melting transitions of lipid bilayers by incorporated proteins, is analyzed by using the same thermodynamic model. Segmental order parameters confirm that elastic lipid chain distortions are insufficient to compensate fully for the mismatch, but the dependence on chain length with tryptophan-anchored peptides requires that the free energy density of hydrophobic mismatch should increase with increasing extent of mismatch.

On the thermodynamic stability of a charged arginine side chain in a transmembrane helix

Biological membranes consist of bilayer arrangements of lipids forming a hydrophobic core that presents a physical barrier to all polar and charged molecules. This long-held notion has recently been challenged by biological translocon-based experiments that report small apparent free energies to insert charged side chains near the center of a transmembrane (TM) helix. We have carried out fully atomistic simulations to provide the free-energy profile for moving a TM helix containing a protonated Arg side chain across a lipid bilayer. Our results reveal the fundamental thermodynamics governing the stability of charged side chains in membranes and the microscopic interactions involved. Despite local membrane deformations, where large amounts of water and lipid head groups are pulled into the bilayer to interact with Arg, the free-energy barrier is 17 kcal/mol. We provide a rationale for the differences in our microscopic free energies and cell biological experiments using free-energy calculations that indicate that a protonated Arg at the central residue of a TM helix of the Leader peptidase might reside close to the interface and not at the membrane center. Our findings have implications for the gating mechanisms of voltage-gated ion channels, suggesting that movements of protonated Arg residues through the membrane will be prohibited.

An amino acid "transmembrane tendency" scale that approaches the theoretical limit to accuracy for prediction of transmembrane helices: relationship to biological hydrophobicity

Hydrophobicity analyses applied to databases of soluble and transmembrane (TM) proteins of known structure were used to resolve total genomic hydrophobicity profiles into (helical) TM sequences and mainly "subhydrophobic" soluble components. This information was used to define a refined "hydrophobicity"-type TM sequence prediction scale that should approach the theoretical limit of accuracy. The refinement procedure involved adjusting scale values to eliminate differences between the average amino acid composition of populations TM and soluble sequences of equal hydrophobicity, a required property of a scale having maximum accuracy. Application of this procedure to different hydrophobicity scales caused them to collapse to essentially a single TM tendency scale. As expected, when different scales were compared, the TM tendency scale was the most accurate at predicting TM sequences. It was especially highly correlated (r = 0.95) to the biological hydrophobicity scale, derived experimentally from the percent TM conformation formed by artificial sequences passing though the translocon. It was also found that resolution of total genomic sequence data into TM and soluble components could be used to define the percent probability that a sequence with a specific hydrophobicity value forms a TM segment. Application of the TM tendency scale to whole genomic data revealed an overlap of TM and soluble sequences in the "semihydrophobic" range. This raises the possibility that a significant number of proteins have sequences that can switch between TM and non-TM states. Such proteins may exist in moonlighting forms having properties very different from those of the predominant conformation.

Investigations of antimicrobial peptides in planar film systems

Planar systems--monolayers and films--constitute a useful platform for studying membrane-active peptides. Here, we summarize varied approaches for studying peptide organization and peptide-lipid interactions at the air/water interface, and focus on three representative antimicrobial membrane--associated peptides-alamethicin, gramicidin, and valinomycin. Experimental data, specifically surface pressure/area isotherms and Brewster angle microscopy images, provided information on peptide association and the effects of the lipid monolayers on peptide surface organization. In general, film analysis emphasized the effects of lipid layers in promoting peptide association and aggregation at the air/water interface. Importantly, the data demonstrated that in many cases peptide domains are phase-separated within the phospholipid monolayers, suggesting that this behavior contributes to the biological actions of membrane-active antimicrobial peptides.

Has the code for protein translocation been broken?

Polypeptides chains are segregated by the translocon channel into secreted or membrane-inserted proteins. Recent reports claim that an in vivo system has been used to break the "amino acid code" used by translocons to make the determination of protein type (i.e. secreted or membrane-inserted). However, the experimental setup used in these studies could have confused the derivation of this code, in particular for polar amino acids. These residues are likely to undergo stabilizing interactions with other protein components in the experiment, shielding them from direct contact with the inhospitable membrane. Hence, it is our view that the "code" for protein translocation has not yet been deciphered and that further experiments are required for teasing apart the various energetic factors contributing to protein translocation.

Computer simulation of partitioning of ten pentapeptides Ace-WLXLL at the cyclohexane/water and phospholipid/water interfaces

BACKGROUND

Peptide-membrane interactions play a key role in the binding, partitioning and folding of membrane proteins, the activity of antimicrobial and fusion peptides, and a number of other processes. To gain a better understanding of the thermodynamics of such interactions, White and Wimley created an interfacial hydrophobicity scale based of the transfer free energy from water to octanol or lipid bilayers of a series of synthetic peptapeptides (Ace-WLXLL, with X being any of the twenty natural amino acids) (White and Wimley (1996) Nat. Struct. Biol. 3, 842-848). In this study, we performed molecular dynamics simulations of a representative set of ten of these peptides (X = D, K, R, N, A, T, S, I, F and W) in two membrane mimetic interfaces: water-cyclohexane (10 ns) and a fully solvated dioleoylphosphatidylcholine (DOPC) bilayer (50 ns) using both constant pressure and constant area ensembles. We focus on partitioning of the ten peptides at the cyclohexane/water and lipid/water interfaces.

RESULTS

The peptides rapidly equilibrate (< 2 ns) and partition at the cyclohexane/water interface. The X3 guest residue assumes average orientations that depend on the nature of the side chain. At the DOPC/water interface, dynamics is much slower and convergence is difficult to achieve on a 50 ns timescale. Nonetheless, all peptides partition to the lipid/water interface with distributions with widths of 1-2 nm. The peptides assume a broad range of side chain and backbone orientations and have only a small effect on the area of the unit cell. On average, hydrophobic guest residues partition deeper into the hydrophobic core than hydrophilic residues. In some cases the peptides penetrate sufficiently deep to somewhat affect the distribution of the C=C double bond in DOPC. The relative distribution of the X3 guest residue compared to W1 and L5 is similar in the water/cyclohexane and water/lipid simulations. Snapshots show mostly extended backbone conformations in both environments. There is little difference between simulations at a constant area of 0.66 nm2 and simulations at constant pressure that approximately yield the same average area of 0.66 nm2.

CONCLUSION

These peptides were designed to assume extended conformations, which is confirmed by the simulations. The distribution of the X3 side chain depends on its nature, and can be determined from molecular dynamics simulations. The time scale of peptide motion at a phospholipids-water interface is too long to directly calculate the experimentally measured hydrophobicity scale to test and improve the simulation parameters. This should be possible at the water/cyclohexane interface and likely will become feasible in the future for the phospholipids/water case.

Melittin-induced bilayer leakage depends on lipid material properties: evidence for toroidal pores

The membrane-lytic peptide melittin has previously been shown to form pores in lipid bilayers that have been described in terms of two different structural models. In the "barrel stave" model the bilayer remains more or less flat, with the peptides penetrating across the bilayer hydrocarbon region and aggregating to form a pore, whereas in the "toroidal pore" melittin induces defects in the bilayer such that the bilayer bends sharply inward to form a pore lined by both peptides and lipid headgroups. Here we test these models by measuring both the free energy of melittin transfer (DeltaG degrees ) and melittin-induced leakage as a function of bilayer elastic (material) properties that determine the energetics of bilayer bending, including the area compressibility modulus (K(a)), bilayer bending modulus (k(c)), and monolayer spontaneous curvature (R(o)). The addition of cholesterol to phosphatidylcholine (PC) bilayers, which increases K(a) and k(c), decreases both DeltaG degrees and the melittin-induced vesicle leakage. In contrast, the addition to PC bilayers of molecules with either positive R(o), such as lysoPC, or negative R(o), such as dioleoylglycerol, has little effect on DeltaG degrees , but produces large changes in melittin-induced leakage, from 86% for 8:2 PC/lysoPC to 18% for 8:2 PC/dioleoylglycerol. We observe linear relationships between melittin-induced leakage and both K(a) and 1/R(o)(2). However, in contrast to what would be expected for a barrel stave model, there is no correlation between observed leakage and bilayer hydrocarbon thickness. All of these results demonstrate the importance of bilayer material properties on melittin-induced leakage and indicate that the melittin-induced pores are defects in the bilayer lined in part by lipid molecules.

A Calpha model for the transmembrane alpha helices of gap junction intercellular channels

Gap junction channels connect the cytoplasms of apposed cells via an intercellular conduit formed by the end-to-end docking of two hexameric hemichannels called connexons. We used electron cryomicroscopy to derive a three-dimensional density map at 5.7 angstroms in-plane and 19.8 angstroms vertical resolution, allowing us to identify the positions and tilt angles for the 24 alpha helices within each hemichannel. The four hydrophobic segments in connexin sequences were assigned to the alpha helices in the map based on biochemical and phylogenetic data. Analyses of evolutionary conservation and compensatory mutations in connexin evolution identified the packing interfaces between the helices. The final model, which specifies the coordinates of Calpha atoms in the transmembrane domain, provides a structural basis for understanding the different physiological effects of almost 30 mutations and polymorphisms in terms of structural deformations at the interfaces between helices, revealing an intimate connection between molecular structure and disease.

An automatic method for predicting transmembrane protein structures using cryo-EM and evolutionary data

The transmembrane (TM) domains of many integral membrane proteins are composed of alpha-helix bundles. Structure determination at high resolution (<4 A) of TM domains is still exceedingly difficult experimentally. Hence, some TM-protein structures have only been solved at intermediate (5-10 A) or low (>10 A) resolutions using, for example, cryo-electron microscopy (cryo-EM). These structures reveal the packing arrangement of the TM domain, but cannot be used to determine the positions of individual amino acids. The observation that typically, the lipid-exposed faces of TM proteins are evolutionarily more variable and less charged than their core provides a simple rule for orienting their constituent helices. Based on this rule, we developed score functions and automated methods for orienting TM helices, for which locations and tilt angles have been determined using, e.g., cryo-EM data. The method was parameterized with the aim of retrieving the native structure of bacteriorhodopsin among near- and far-from-native templates. It was then tested on proteins that differ from bacteriorhodopsin in their sequences, architectures, and functions, such as the acetylcholine receptor and rhodopsin. The predicted structures were within 1.5-3.5 A from the native state in all cases. We conclude that the computational method can be used in conjunction with cryo-EM data to obtain approximate model structures of TM domains of proteins for which a sufficiently heterogeneous set of homologs is available. We also show that in those proteins in which relatively short loops connect neighboring helices, the scoring functions can discriminate between near- and far-from-native conformations even without the constraints imposed on helix locations and tilt angles that are derived from cryo-EM.

Interactions of the M2delta segment of the acetylcholine receptor with lipid bilayers: a continuum-solvent model study

M2delta, one of the transmembrane segments of the nicotinic acetylcholine receptor, is a 23-amino-acid peptide, frequently used as a model for peptide-membrane interactions. In this and the companion article we describe studies of M2delta-membrane interactions, using two different computational approaches. In the present work, we used continuum-solvent model calculations to investigate key thermodynamic aspects of its interactions with lipid bilayers. M2delta was represented in atomic detail and the bilayer was represented as a hydrophobic slab embedded in a structureless aqueous phase. Our calculations show that the transmembrane orientation is the most favorable orientation of the peptide in the bilayer, in good agreement with both experimental and computational data. Moreover, our calculations produced the free energy of association of M2delta with the lipid bilayer, which, to our knowledge, has not been reported to date. The calculations included 10 structures of M2delta, determined by nuclear magnetic resonance in dodecylphosphocholine micelles. All the structures were found to be stable inside the lipid bilayer, although their water-to-membrane transfer free energies differed by as much as 12 kT. Although most of the structures were roughly linear, a single structure had a kink in its central region. Interestingly, this structure was found to be the most stable inside the lipid bilayer, in agreement with molecular dynamics simulations of the peptide and with the recently determined structure of the intact receptor. Our analysis showed that the kink reduced the polarity of the peptide in its central region by allowing the electrostatic masking of the Gln13 side chain in that area. Our calculations also showed a tendency for the membrane to deform in response to peptide insertion, as has been previously found for the membrane-active peptides alamethicin and gramicidin. The results are compared to Monte Carlo simulations of the peptide-membrane system, as presented in the accompanying article.

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.

Lipid-protein interactions in biological membranes: a structural perspective

Lipid molecules bound to membrane proteins are resolved in some high-resolution structures of membrane proteins. An analysis of these structures provides a framework within which to analyse the nature of lipid-protein interactions within membranes. Membrane proteins are surrounded by a shell or annulus of lipid molecules, equivalent to the solvent layer surrounding a water-soluble protein. The lipid bilayer extends right up to the membrane protein, with a uniform thickness around the protein. The surface of a membrane protein contains many shallow grooves and protrusions to which the fatty acyl chains of the surrounding lipids conform to provide tight packing into the membrane. An individual lipid molecule will remain in the annular shell around a protein for only a short period of time. Binding to the annular shell shows relatively little structural specificity. As well as the annular lipid, there is evidence for other lipid molecules bound between the transmembrane alpha-helices of the protein; these lipids are referred to as non-annular lipids. The average thickness of the hydrophobic domain of a membrane protein is about 29 A, with a few proteins having significantly smaller or greater thicknesses than the average. Hydrophobic mismatch between a membrane protein and the surrounding lipid bilayer generally leads to only small changes in membrane thickness. Possible adaptations in the protein to minimise mismatch include tilting of the helices and rotation of side chains at the ends of the helices. Packing of transmembrane alpha-helices is dependent on the chain length of the surrounding phospholipids. The function of membrane proteins is dependent on the thickness of the surrounding lipid bilayer, sometimes on the presence of specific, usually anionic, phospholipids, and sometimes on the phase of the phospholipid.

Stability of an ion channel in lipid bilayers: implicit solvent model calculations with gramicidin

Gramicidin is a helical peptide, 15 residues in length, which dimerizes to form ion-conducting channels in lipid bilayers. Here we report calculations of its free energy of transfer from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the desolvation free energy were calculated using implicit solvent models, in which gramicidin was described in atomic detail and the hydrocarbon region of the membrane was described as a slab of hydrophobic medium embedded in water. The free energy penalties from the lipid perturbation and membrane deformation effects, and the entropy loss associated with gramicidin immobilization in the bilayer, were estimated from a statistical thermodynamic model of the bilayer. The calculations were carried out using two classes of experimentally observed conformations: a head-to-head dimer of two single-stranded (SS) beta-helices and a double-stranded (DS) intertwined double helix. The calculations showed that gramicidin is likely to partition into the bilayer in all of these conformations. However, the SS conformation was found to be significantly more stable than the DS in the bilayer, in agreement with most of the experimental data. We tested numerous transmembrane and surface orientations of gramicidin in bilayers of various widths. Our calculations indicate that the most favorable orientation is transmembrane, which is indeed to be expected from a channel-forming peptide. The calculations demonstrate that gramicidin insertion into the membrane is likely to involve a significant deformation of the bilayer to match the hydrophobic width of the peptide (22 A), again in good agreement with experimental data. Interestingly, deformation of the bilayer was induced by all of the gramicidin conformations.

Continuum solvent model studies of the interactions of an anticonvulsant drug with a lipid bilayer

Valproic acid (VPA) is a short, branched fatty acid with broad-spectrum anticonvulsant activity. It has been suggested that VPA acts directly on the plasma membrane. We calculated the free energy of interaction of VPA with a model lipid bilayer using simulated annealing and the continuum solvent model. Our calculations indicate that VPA is likely to partition into the bilayer both in its neutral and charged forms, as expected from such an amphipathic molecule. The calculations also show that VPA may migrate (flip-flop) across the membrane; according to our (theoretical) study, the most likely flip-flop path at neutral pH involves protonation of VPA pending its insertion into the lipid bilayer and deprotonation upon departure from the other side of the bilayer. Recently, the flip-flop of long fatty acids across lipid bilayers was studied using fluorescence and NMR spectroscopies. However, the measured value of the flip-flop rate appears to depend on the method used in these studies. Our calculated value of the flip-flop rate constant, 20/s, agrees with some of these studies. The limitations of the model and the implications of the study for VPA and other fatty acids are discussed.

Interactions of cholesterol with lipid bilayers: the preferred configuration and fluctuations

The free energy difference associated with the transfer of a single cholesterol molecule from the aqueous phase into a lipid bilayer depends on its final location, namely on its insertion depth and orientation within the bilayer. We calculated desolvation and lipid bilayer perturbation contributions to the water-to-membrane transfer free energy, thus allowing us to determine the most favorable location of cholesterol in the membrane and the extent of fluctuations around it. The electrostatic and nonpolar contributions to the solvation free energy were calculated using continuum solvent models. Lipid layer perturbations, resulting from both conformational restrictions of the lipid chains in the vicinity of the (rigid) cholesterol backbone and from cholesterol-induced elastic deformations, were calculated using a simple director model and elasticity theory, respectively. As expected from the amphipathic nature of cholesterol and in agreement with the available experimental data, our results show that at the energetically favorable state, cholesterol's hydrophobic core is buried within the hydrocarbon region of the bilayer. At this state, cholesterol spans approximately one leaflet of the membrane, with its OH group protruding into the polar (headgroup) region of the bilayer, thus avoiding an electrostatic desolvation penalty. We found that the transfer of cholesterol into a membrane is mainly driven by the favorable nonpolar contributions to the solvation free energy, whereas only a small opposing contribution is caused by conformational restrictions of the lipid chains. Our calculations also predict a strong tendency of the lipid layer to elastically respond to (thermally excited) vertical fluctuations of cholesterol so as to fully match the hydrophobic height of the solute. However, orientational fluctuations of cholesterol were found to be accompanied by both an elastic adjustment of the surrounding lipids and by a partial exposure of the hydrophobic cholesterol backbone to the polar (headgroup) environment. Our calculations of the molecular order parameter, which reflects the extent of orientational fluctuations of cholesterol in the membrane, are in good agreement with available experimental data.