Hydrophobicity scales: a thermodynamic looking glass into lipid-protein interactions

From PDB files to protein features: a comparative analysis of PDB bind and STCRDAB datasets

Understanding protein structures is crucial for various bioinformatics research, including drug discovery, disease diagnosis, and evolutionary studies. Protein structure classification is a critical aspect of structural biology, where supervised machine learning algorithms classify structures based on data from databases such as Protein Data Bank (PDB). However, the challenge lies in designing numerical embeddings for protein structures without losing essential information. Although some effort has been made in the literature, researchers have not effectively and rigorously combined the structural and sequence-based features for efficient protein classification to the best of our knowledge. To this end, we propose numerical embeddings that extract relevant features for protein sequences fetched from PDB structures from popular datasets such as PDB Bind and STCRDAB. The features are physicochemical properties such as aromaticity, instability index, flexibility, Grand Average of Hydropathy (GRAVY), isoelectric point, charge at pH, secondary structure fracture, molar extinction coefficient, and molecular weight. We also incorporate scaling features for the sliding windows (e.g., k-mers), which include Kyte and Doolittle (KD) hydropathy scale, Eisenberg hydrophobicity scale, Hydrophilicity scale, Flexibility of the amino acids, and Hydropathy scale. Multiple-feature selection aims to improve the accuracy of protein classification models. The results showed that the selected features significantly improved the predictive performance of existing embeddings.

AAontology: An ontology of amino acid scales for interpretable machine learning

Amino acid scales are crucial for protein prediction tasks, many of them being curated in the AAindex database. Despite various clustering attempts to organize them and to better understand their relationships, these approaches lack the fine-grained classification necessary for satisfactory interpretability in many protein prediction problems. To address this issue, we developed AAontology-a two-level classification for 586 amino acid scales (mainly from AAindex) together with an in-depth analysis of their relations-using bag-of-word-based classification, clustering, and manual refinement over multiple iterations. AAontology organizes physicochemical scales into 8 categories and 67 subcategories, enhancing the interpretability of scale-based machine learning methods in protein bioinformatics. Thereby it enables researchers to gain a deeper biological insight. We anticipate that AAontology will be a building block to link amino acid properties with protein function and dysfunctions as well as aid informed decision-making in mutation analysis or protein drug design.

Main Factors Shaping Amino Acid Usage Across Evolution

The standard genetic code determines that in most species, including viruses, there are 20 amino acids that are coded by 61 codons, while the other three codons are stop triplets. Considering the whole proteome each species features its own amino acid frequencies, given the slow rate of change, closely related species display similar GC content and amino acids usage. In contrast, distantly related species display different amino acid frequencies. Furthermore, within certain multicellular species, as mammals, intragenomic differences in the usage of amino acids are evident. In this communication, we shall summarize some of the most prominent and well-established factors that determine the differences found in the amino acid usage, both across evolution and intragenomically.

Structural features of interfacial water predict the hydrophobicity of chemically heterogeneous surfaces

The hydrophobicity of an interface determines the magnitude of hydrophobic interactions that drive numerous biological and industrial processes. Chemically heterogeneous interfaces are abundant in these contexts; examples include the surfaces of proteins, functionalized nanomaterials, and polymeric materials. While the hydrophobicity of nonpolar solutes can be predicted and related to the structure of interfacial water molecules, predicting the hydrophobicity of chemically heterogeneous interfaces remains a challenge because of the complex, non-additive contributions to hydrophobicity that depend on the chemical identity and nanoscale spatial arrangements of polar and nonpolar groups. In this work, we utilize atomistic molecular dynamics simulations in conjunction with enhanced sampling and data-centric analysis techniques to quantitatively relate changes in interfacial water structure to the hydration free energy (a thermodynamically well-defined descriptor of hydrophobicity) of chemically heterogeneous interfaces. We analyze a large data set of 58 self-assembled monolayers (SAMs) composed of ligands with nonpolar and polar end groups of different chemical identity (amine, amide, and hydroxyl) in five mole fractions, two spatial patterns, and with scaled partial charges. We find that only five features of interfacial water structure are required to accurately predict hydration free energies. Examination of these features reveals mechanistic insights into the interfacial hydrogen bonding behaviors that distinguish different surface compositions and patterns. This analysis also identifies the probability of highly coordinated water structures as a unique signature of hydrophobicity. These insights provide a physical basis to understand the hydrophobicity of chemically heterogeneous interfaces and connect hydrophobicity to experimentally accessible perturbations of interfacial water structure.

Comparison of hydrophobicity scales for predicting biophysical properties of antibodies

While antibody-based therapeutics have grown to be one of the major classes of novel medicines, some antibody development candidates face significant challenges regarding expression levels, solubility, as well as stability and aggregation, under physiological and storage conditions. A major determinant of those properties is surface hydrophobicity, which promotes unspecific interactions and has repeatedly proven problematic in the development of novel antibody-based drugs. Multiple computational methods have been devised for in-silico prediction of antibody hydrophobicity, often using hydrophobicity scales to assign values to each amino acid. Those approaches are usually validated by their ability to rank potential therapeutic antibodies in terms of their experimental hydrophobicity. However, there is significant diversity both in the hydrophobicity scales and in the experimental methods, and consequently in the performance of in-silico methods to predict experimental results. In this work, we investigate hydrophobicity of monoclonal antibodies using hydrophobicity scales. We implement several scoring schemes based on the solvent-accessibility and the assigned hydrophobicity values, and compare the different scores and scales based on their ability to predict retention times from hydrophobic interaction chromatography. We provide an overview of the strengths and weaknesses of several commonly employed hydrophobicity scales, thereby improving the understanding of hydrophobicity in antibody development. Furthermore, we test several datasets, both publicly available and proprietary, and find that the diversity of the dataset affects the performance of hydrophobicity scores. We expect that this work will provide valuable guidelines for the optimization of biophysical properties in future drug discovery campaigns.

Amphiphilic Gold Nanoparticles: A Biomimetic Tool to Gain Mechanistic Insights into Peptide-Lipid Interactions

Functional peptides are now widely used in a myriad of biomedical and clinical contexts, from cancer therapy and tumor targeting to the treatment of bacterial and viral infections. Underlying this diverse range of applications are the non-specific interactions that can occur between peptides and cell membranes, which, in many contexts, result in spontaneous internalization of the peptide within cells by avoiding energy-driven endocytosis. For this to occur, the amphipathicity and surface structural flexibility of the peptides play a crucial role and can be regulated by the presence of specific molecular residues that give rise to precise molecular events. Nevertheless, most of the mechanistic details regulating the encounter between peptides and the membranes of bacterial or animal cells are still poorly understood, thus greatly limiting the biomimetic potential of these therapeutic molecules. In this arena, finely engineered nanomaterials—such as small amphiphilic gold nanoparticles (AuNPs) protected by a mixed thiol monolayer—can provide a powerful tool for mimicking and investigating the physicochemical processes underlying peptide-lipid interactions. Within this perspective, we present here a critical review of membrane effects induced by both amphiphilic AuNPs and well-known amphiphilic peptide families, such as cell-penetrating peptides and antimicrobial peptides. Our discussion is focused particularly on the effects provoked on widely studied model cell membranes, such as supported lipid bilayers and lipid vesicles. Remarkable similarities in the peptide or nanoparticle membrane behavior are critically analyzed. Overall, our work provides an overview of the use of amphiphilic AuNPs as a highly promising tailor-made model to decipher the molecular events behind non-specific peptide-lipid interactions and highlights the main affinities observed both theoretically and experimentally. The knowledge resulting from this biomimetic approach could pave the way for the design of synthetic peptides with tailored functionalities for next-generation biomedical applications, such as highly efficient intracellular delivery systems.

Sense the moment: A highly sensitive antimicrobial activity predictor based on hydrophobic moment

BACKGROUND

Computer-aided identification and design tools are indispensable for developing antimicrobial agents for controlling antibiotic-resistant bacteria. Antimicrobial peptides (AMPs) have aroused intense interest, since they have a broad spectrum of activity, and therefore, several systems for predicting antimicrobial peptides have been developed, using scalar physicochemical properties; however, regardless of the machine learning algorithm, these systems often fail in discriminating AMPs from their shuffled versions, leading to the need for new training methods to overcome this bias. Aiming to solve this bias, here we present "Sense the Moment", a prediction system capable of discriminating AMPs and shuffled versions.

METHODS

The system was trained using 776 entries: 388 from known AMPs and another 388 based on shuffled versions of known AMPs. Each entry contained the geometric average of three hydrophobic moments measured with different scales.

RESULTS

The model showed good accuracy (>80%) and excellent sensitivity (>90%) for AMP prediction, exceeding deep-learning-based methods.

CONCLUSION

Our results demonstrate the system's applicability, aiding in identifying and discarding non-AMPs, since the number of false negatives is lower than false positives.

GENERAL SIGNIFICANCE

The application of this model in virtual screening protocols for identifying and/or creating antimicrobial agents could aid in the identification of potential drugs to control pathogenic microorganisms and in solving the antibiotic resistance crisis.

AVAILABILITY

The system was implemented as a web application, available at .

Identifying nonadditive contributions to the hydrophobicity of chemically heterogeneous surfaces via dual-loop active learning

Hydrophobic interactions drive numerous biological and synthetic processes. The materials used in these processes often possess chemically heterogeneous surfaces that are characterized by diverse chemical groups positioned in close proximity at the nanoscale; examples include functionalized nanomaterials and biomolecules, such as proteins and peptides. Nonadditive contributions to the hydrophobicity of such surfaces depend on the chemical identities and spatial patterns of polar and nonpolar groups in ways that remain poorly understood. Here, we develop a dual-loop active learning framework that combines a fast reduced-accuracy method (a convolutional neural network) with a slow higher-accuracy method (molecular dynamics simulations with enhanced sampling) to efficiently predict the hydration free energy, a thermodynamic descriptor of hydrophobicity, for nearly 200 000 chemically heterogeneous self-assembled monolayers (SAMs). Analysis of this dataset reveals that SAMs with distinct polar groups exhibit substantial variations in hydrophobicity as a function of their composition and patterning, but the clustering of nonpolar groups is a common signature of highly hydrophobic patterns. Further molecular dynamics analysis relates such clustering to the perturbation of interfacial water structure. These results provide new insight into the influence of chemical heterogeneity on hydrophobicity via quantitative analysis of a large set of surfaces, enabled by the active learning approach.

In Silico Prediction of Steroids and Triterpenoids as Potential Regulators of Lipid Metabolism

This review focuses on a rare group of steroids and triterpenoids that share common properties as regulators of lipid metabolism. This group of compounds is divided by the type of chemical structure, and they represent: aromatic steroids, steroid phosphate esters, highly oxygenated steroids such as steroid endoperoxides and hydroperoxides, α,β-epoxy steroids, and secosteroids. In addition, subgroups of carbon-bridged steroids, neo steroids, miscellaneous steroids, as well as synthetic steroids containing heteroatoms S (epithio steroids), Se (selena steroids), Te (tellura steroids), and At (astatosteroids) were presented. Natural steroids and triterpenoids have been found and identified from various sources such as marine sponges, soft corals, starfish, and other marine invertebrates. In addition, this group of rare lipids is found in fungi, fungal endophytes, and plants. The pharmacological profile of the presented steroids and triterpenoids was determined using the well-known computer program PASS, which is currently available online for all interested scientists and pharmacologists and is currently used by research teams from more than 130 countries of the world. Our attention has been focused on the biological activities of steroids and triterpenoids associated with the regulation of cholesterol metabolism and related processes such as anti-hyperlipoproteinemic activity, as well as the treatment of atherosclerosis, lipoprotein disorders, or inhibitors of cholesterol synthesis. In addition, individual steroids and triterpenoids were identified that demonstrated rare or unique biological activities such as treating neurodegenerative diseases, Alzheimer's, and Parkinson's diseases with a high degree of certainty over 95 percent. For individual steroids or triterpenoids or a group of compounds, 3D drawings of their predicted biological activities are presented.

Cryo-EM Structure of Mechanosensitive Channel YnaI Using SMA2000: Challenges and Opportunities

Mechanosensitive channels respond to mechanical forces exerted on the cell membrane and play vital roles in regulating the chemical equilibrium within cells and their environment. High-resolution structural information is required to understand the gating mechanisms of mechanosensitive channels. Protein-lipid interactions are essential for the structural and functional integrity of mechanosensitive channels, but detergents cannot maintain the crucial native lipid environment for purified mechanosensitive channels. Recently, detergent-free systems have emerged as alternatives for membrane protein structural biology. This report shows that while membrane-active polymer, SMA2000, could retain some native cell membrane lipids on the transmembrane domain of the mechanosensitive-like YnaI channel, the complete structure of the transmembrane domain of YnaI was not resolved. This reveals a significant limitation of SMA2000 or similar membrane-active copolymers. This limitation may come from the heterogeneity of the polymers and nonspecific interactions between the polymers and the relatively large hydrophobic pockets within the transmembrane domain of YnaI. However, this limitation offers development opportunities for detergent-free technology for challenging membrane proteins.

3D interaction homology: Hydropathic interaction environments of serine and cysteine are strikingly different and their roles adapt in membrane proteins

Atomic-resolution protein structural models are prerequisites for many downstream activities like structure-function studies or structure-based drug discovery. Unfortunately, this data is often unavailable for some of the most interesting and therapeutically important proteins. Thus, computational tools for building native-like structural models from less-than-ideal experimental data are needed. To this end, interaction homology exploits the character, strength and loci of the sets of interactions that define a structure. Each residue type has its own limited set of backbone angle-dependent interaction motifs, as defined by their environments. In this work, we characterize the interactions of serine, cysteine and S-bridged cysteine in terms of 3D hydropathic environment maps. As a result, we explore several intriguing questions. Are the environments different between the isosteric serine and cysteine residues? Do some environments promote the formation of cystine S-S bonds? With the increasing availability of structural data for water-insoluble membrane proteins, are there environmental differences for these residues between soluble and membrane proteins? The environments surrounding serine and cysteine residues are dramatically different: serine residues are about 50% solvent exposed, while cysteines are only 10% exposed; the latter are more involved in hydrophobic interactions although there are backbone angle-dependent differences. Our analysis suggests that one driving force for -S-S- bond formation is a rather substantial increase in burial and hydrophobic interactions in cystines. Serine and cysteine become less and more, respectively, solvent-exposed in membrane proteins. 3D hydropathic environment maps are an evolving structure analysis tool showing promise as elements in a new protein structure prediction paradigm.

Ceramide-1-phosphate transfer protein (CPTP) regulation by phosphoinositides

Ceramide-1-phosphate transfer proteins (CPTPs) are members of the glycolipid transfer protein (GLTP) superfamily that shuttle ceramide-1-phosphate (C1P) between membranes. CPTPs regulate cellular sphingolipid homeostasis in ways that impact programmed cell death and inflammation. CPTP downregulation specifically alters C1P levels in the plasma and trans-Golgi membranes, stimulating proinflammatory eicosanoid production and autophagy-dependent inflammasome-mediated cytokine release. However, the mechanisms used by CPTP to target the trans-Golgi and plasma membrane are not well understood. Here, we monitored C1P intervesicular transfer using fluorescence energy transfer (FRET) and showed that certain phosphoinositides (phosphatidylinositol 4,5 bisphosphate (PI-(4,5)P2) and phosphatidylinositol 4-phosphate (PI-4P)) increased CPTP transfer activity, whereas others (phosphatidylinositol 3-phosphate (PI-3P) and PI) did not. PIPs that stimulated CPTP did not stimulate GLTP, another superfamily member. Short-chain PI-(4,5)P2, which is soluble and does not remain membrane-embedded, failed to activate CPTP. CPTP stimulation by physiologically relevant PI-(4,5)P2 levels surpassed that of phosphatidylserine (PS), the only known non-PIP stimulator of CPTP, despite PI-(4,5)P2 increasing membrane equilibrium binding affinity less effectively than PS. Functional mapping of mutations that led to altered FRET lipid transfer and assessment of CPTP membrane interaction by surface plasmon resonance indicated that di-arginine motifs located in the α-6 helix and the α3-α4 helix regulatory loop of the membrane-interaction region serve as PI-(4,5)P2 headgroup-specific interaction sites. Haddock modeling revealed specific interactions involving the PI-(4,5)P2 headgroup that left the acyl chains oriented favorably for membrane embedding. We propose that PI-(4,5)P2 interaction sites enhance CPTP activity by serving as preferred membrane targeting/docking sites that favorably orient the protein for function.

Computational prediction of protein aggregation: Advances in proteomics, conformation-specific algorithms and biotechnological applications

Protein aggregation is a widespread phenomenon that stems from the establishment of non-native intermolecular contacts resulting in protein precipitation. Despite its deleterious impact on fitness, protein aggregation is a generic property of polypeptide chains, indissociable from protein structure and function. Protein aggregation is behind the onset of neurodegenerative disorders and one of the serious obstacles in the production of protein-based therapeutics. The development of computational tools opened a new avenue to rationalize this phenomenon, enabling prediction of the aggregation propensity of individual proteins as well as proteome-wide analysis. These studies spotted aggregation as a major force driving protein evolution. Actual algorithms work on both protein sequences and structures, some of them accounting also for conformational fluctuations around the native state and the protein microenvironment. This toolbox allows to delineate conformation-specific routines to assist in the identification of aggregation-prone regions and to guide the optimization of more soluble and stable biotherapeutics. Here we review how the advent of predictive tools has change the way we think and address protein aggregation.

pH-Dependent Aggregation in Intrinsically Disordered Proteins Is Determined by Charge and Lipophilicity

Protein aggregation is associated with an increasing number of human disorders and premature aging. Moreover, it is a central concern in the manufacturing of recombinant proteins for biotechnological and therapeutic applications. Nevertheless, the unique architecture of protein aggregates is also exploited by nature for functional purposes, from bacteria to humans. The relevance of this process in health and disease has boosted the interest in understanding and controlling aggregation, with the concomitant development of a myriad of algorithms aimed to predict aggregation propensities. However, most of these programs are blind to the protein environment and, in particular, to the influence of the pH. Here, we developed an empirical equation to model the pH-dependent aggregation of intrinsically disordered proteins (IDPs) based on the assumption that both the global protein charge and lipophilicity depend on the solution pH. Upon its parametrization with a model IDP, this simple phenomenological approach showed unprecedented accuracy in predicting the dependence of the aggregation of both pathogenic and functional amyloidogenic IDPs on the pH. The algorithm might be useful for diverse applications, from large-scale analysis of IDPs aggregation properties to the design of novel reversible nanofibrillar materials.

Solvation Free Energy as a Measure of Hydrophobicity: Application to Serine Protease Binding Interfaces

Solvation and hydrophobicity play a key role in a variety of biological mechanisms. In substrate binding, but also in structure-based drug design, the thermodynamic properties of water molecules surrounding a given protein are of high interest. One of the main algorithms devised in recent years to quantify thermodynamic properties of water is the grid inhomogeneous solvation theory (GIST), which calculates these features on a grid surrounding the protein. Despite the inherent advantages of GIST, the computational demand is a major drawback, as calculations for larger systems can take days or even weeks. Here, we present a GPU accelerated version of the GIST algorithm, which facilitates efficient estimates of solvation free energy even of large biomolecular interfaces. Furthermore, we show that GIST can be used as a reliable tool to evaluate protein surface hydrophobicity. We apply the approach on a set of nine different proteases calculating localized solvation free energies on the surface of the binding interfaces as a measure of their hydrophobicity. We find a compelling agreement with the hydrophobicity of their substrates, i.e., peptides, binding into the binding cleft, and thus our approach provides a reliable description of hydrophobicity characteristics of these biological interfaces.

GeTFEP: A general transfer free energy profile of transmembrane proteins

Free energy of transferring amino acid side-chains from aqueous environment into lipid bilayers, known as transfer free energy (TFE), provides important information on the thermodynamic stability of membrane proteins. In this study, we derived a TFE profile named General Transfer Free Energy Profile (GeTFEP) based on computation of the TFEs of 58 β-barrel membrane proteins (βMPs). The GeTFEP agrees well with experimentally measured and computationally derived TFEs. Analysis based on the GeTFEP shows that residues in different regions of the transmembrane (TM) segments of βMPs have different roles during the membrane insertion process. Results further reveal the importance of the sequence pattern of TM strands in stabilizing βMPs in the membrane environment. In addition, we show that GeTFEP can be used to predict the positioning and the orientation of βMPs in the membrane. We also show that GeTFEP can be used to identify structurally or functionally important amino acid residue sites of βMPs. Furthermore, the TM segments of α-helical membrane proteins can be accurately predicted with GeTFEP, suggesting that the GeTFEP is of general applicability in studying membrane protein.

Hydrophobic Characteristic Is Energetically Preferred for Cysteine in a Model Membrane Protein

The naturally occurring amino acid cysteine has often been implicated with a crucial role in maintaining protein structure and stability. An intriguing duality in the intrinsic hydrophobicity of the cysteine side chain is that it exhibits both polar as well as hydrophobic characteristics. Here, we have utilized a cysteine-scanning mutational strategy on the transmembrane β-barrel PagP to examine the membrane depth-dependent energetic contribution of the free cysteine side chain (thiolate) versus the parent residue at an experimental pH of 9.5 in phosphatidylcholine vesicles. We find that introduction of cysteine causes destabilization at several of the 26 lipid-facing sites of PagP that we mutated in this study. The destabilization is minimal (0.5-1.5 kcal/mol) when the mutation is toward the bilayer midplane, whereas it is higher in magnitude (3.0-5.0 kcal/mol) near the bilayer interface. These observations suggest that cysteine forms more favorable interactions with the hydrophobic lipid core as compared to the amphiphilic water-lipid interface. The destabilizing effect is more pronounced when cysteine replaces the interfacial aromatics, which are known to participate in tertiary interaction networks in transmembrane β-barrels. Our observations from experiments involving the introduction of cysteine at the bilayer midplane further strengthen previous views that the free cysteine side chain does possess strongly apolar characteristics. Additionally, the free energy changes observed upon cysteine incorporation show a depth-dependent correlation with the estimated energetic cost of partitioning derived from reported hydrophobicity scales. Our results and observations from the thermodynamic analysis of the PagP barrel may explain why cysteine, despite possessing a polar sulfhydryl group, tends to behave as a hydrophobic (rather than polar) residue in folded protein structures.

Development of a Structure-Based, pH-Dependent Lipophilicity Scale of Amino Acids from Continuum Solvation Calculations

Lipophilicity is a fundamental property to characterize the structure and function of proteins, motivating the development of lipophilicity scales. We report a versatile strategy to derive a pH-adapted scale that relies on theoretical estimates of distribution coefficients from conformational ensembles of amino acids. This is accomplished by using an accurately parametrized version of the IEFPCM/MST continuum solvation model as an effective way to describe the partitioning between n-octanol and water, in conjunction with a formalism that combines partition coefficients of neutral and ionic species of residues and the corresponding p K_a values of ionizable groups. Two weighting schemes are considered to derive solvent-like and protein-like scales, which have been calibrated by comparison with other experimental scales developed in different chemical/biological environments and pH conditions as well as by examining properties such as the retention time of small peptides and the recognition of antigenic peptides. A straightforward extension to nonstandard residues is enabled by this efficient methodological strategy.

Coarse-grained model of titrating peptides interacting with lipid bilayers

Molecular-level computer simulations of peptide aggregation, translocation, and protonation at and in biomembranes are impeded by the large time and length scales involved. We present a computationally efficient, coarse-grained, and solvent-free model for the interaction between lipid bilayers and peptides. The model combines an accurate description of mechanical membrane properties with a new granular representation of the dielectric mismatch between lipids and the aqueous phase. All-atom force fields can be easily mapped onto the coarse-grained model, and parameters for coarse-grained monopeptides accurately extrapolate to membrane permeation free energies for the corresponding dipeptides and tripeptides. Acid-base equilibria of titratable amino acid residues are further studied using a constant-pH ensemble, capturing protonation state changes upon membrane translocation. Important differences between histidine, lysine, and arginine are observed, which are in good agreement with experimental observations.

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.

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.

Computational membrane biophysics: From ion channel interactions with drugs to cellular function

The rapid development of experimental and computational techniques has changed fundamentally our understanding of cellular-membrane transport. The advent of powerful computers and refined force-fields for proteins, ions, and lipids has expanded the applicability of Molecular Dynamics (MD) simulations. A myriad of cellular responses is modulated through the binding of endogenous and exogenous ligands (e.g. neurotransmitters and drugs, respectively) to ion channels. Deciphering the thermodynamics and kinetics of the ligand binding processes to these membrane proteins is at the heart of modern drug development. The ever-increasing computational power has already provided insightful data on the thermodynamics and kinetics of drug-target interactions, free energies of solvation, and partitioning into lipid bilayers for drugs. This review aims to provide a brief summary about modeling approaches to map out crucial binding pathways with intermediate conformations and free-energy surfaces for drug-ion channel binding mechanisms that are responsible for multiple effects on cellular functions. We will discuss post-processing analysis of simulation-generated data, which are then transformed to kinetic models to better understand the molecular underpinning of the experimental observables under the influence of drugs or mutations in ion channels. This review highlights crucial mathematical frameworks and perspectives on bridging different well-established computational techniques to connect the dynamics and timescales from all-atom MD and free energy simulations of ion channels to the physiology of action potentials in cellular models. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Characterizing Residue-Bilayer Interactions Using Gramicidin A as a Scaffold and Tryptophan Substitutions as Probes

Previous experiments have shown that the lifetime of a gramicidin A dimer channel (which forms from two nonconducting monomers) in a lipid bilayer is modulated by mutations of the tryptophan (Trp) residues at the bilayer-water interface. We explore this further using extensive molecular dynamics simulations of various gA dimer and monomer mutants at the Trp positions in phosphatidylcholine bilayers with different tail lengths. gA interactions with the surrounding bilayer are strongly modulated by mutating these Trp residues. There are three principal effects: eliminating residue hydrogen bonding ability (i.e., reducing the channel-monolayer coupling strength) reduces the extent of the bilayer deformation caused by the assembled dimeric channel; a residue's size and geometry affects its orientation, leading to different hydrogen bonding partners; and increasing a residue's hydrophobicity increases the depth of gA monomer insertion relative to the bilayer center, thereby increasing the lipid bending frustration.

Energetics of side-chain partitioning of β-signal residues in unassisted folding of a transmembrane β-barrel protein

The free energy of water-to-interface amino acid partitioning is a major contributing factor in membrane protein folding and stability. The interface residues at the C terminus of transmembrane β-barrels form the β-signal motif required for assisted β-barrel assembly in vivo but are believed to be less important for β-barrel assembly in vitro. Here, we experimentally measured the thermodynamic contribution of all 20 amino acids at the β-signal motif to the unassisted folding of the model β-barrel protein PagP. We obtained the partitioning free energy for all 20 amino acids at the lipid-facing interface (ΔΔG⁰_w,i(φ)) and the protein-facing interface (ΔΔG⁰_w,i(π)) residues and found that hydrophobic amino acids are most favorably transferred to the lipid-facing interface, whereas charged and polar groups display the highest partitioning energy. Furthermore, the change in non-polar surface area correlated directly with the partitioning free energy for the lipid-facing residue and inversely with the protein-facing residue. We also demonstrate that the interface residues of the β-signal motif are vital for in vitro barrel assembly, because they exhibit a side chain–specific energetic contribution determined by the change in nonpolar accessible surface. We further establish that folding cooperativity and hydrophobic collapse are balanced at the membrane interface for optimal stability of the PagP β-barrel scaffold. We conclude that the PagP C-terminal β-signal motif influences the folding cooperativity and stability of the folded β-barrel and that the thermodynamic contributions of the lipid- and protein-facing residues in the transmembrane protein β-signal motif depend on the nature of the amino acid side chain.

Transmembrane helices containing a charged arginine are thermodynamically stable

Hydrophobic amino acids are abundant in transmembrane (TM) helices of membrane proteins. Charged residues are sparse, apparently due to the unfavorable energetic cost of partitioning charges into nonpolar phases. Nevertheless, conserved arginine residues within TM helices regulate vital functions, such as ion channel voltage gating and integrin receptor inactivation. The energetic cost of arginine in various positions along hydrophobic helices has been controversial. Potential of mean force (PMF) calculations from atomistic molecular dynamics simulations predict very large energetic penalties, while in vitro experiments with Sec61 translocons indicate much smaller penalties, even for arginine in the center of hydrophobic TM helices. Resolution of this conflict has proved difficult, because the in vitro assay utilizes the complex Sec61 translocon, while the PMF calculations rely on the choice of simulation system and reaction coordinate. Here we present the results of computational and experimental studies that permit direct comparison with the Sec61 translocon results. We find that the Sec61 translocon mediates less efficient membrane insertion of Arg-containing TM helices compared with our computational and experimental bilayer-insertion results. In the simulations, a combination of arginine snorkeling, bilayer deformation, and peptide tilting is sufficient to lower the penalty of Arg insertion to an extent such that a hydrophobic TM helix with a central Arg residue readily inserts into a model membrane. Less favorable insertion by the translocon may be due to the decreased fluidity of the endoplasmic reticulum (ER) membrane compared with pure palmitoyloleoyl-phosphocholine (POPC). Nevertheless, our results provide an explanation for the differences between PMF- and experiment-based penalties for Arg burial.

Structurally detailed coarse-grained model for Sec-facilitated co-translational protein translocation and membrane integration

We present a coarse-grained simulation model that is capable of simulating the minute-timescale dynamics of protein translocation and membrane integration via the Sec translocon, while retaining sufficient chemical and structural detail to capture many of the sequence-specific interactions that drive these processes. The model includes accurate geometric representations of the ribosome and Sec translocon, obtained directly from experimental structures, and interactions parameterized from nearly 200 μs of residue-based coarse-grained molecular dynamics simulations. A protocol for mapping amino-acid sequences to coarse-grained beads enables the direct simulation of trajectories for the co-translational insertion of arbitrary polypeptide sequences into the Sec translocon. The model reproduces experimentally observed features of membrane protein integration, including the efficiency with which polypeptide domains integrate into the membrane, the variation in integration efficiency upon single amino-acid mutations, and the orientation of transmembrane domains. The central advantage of the model is that it connects sequence-level protein features to biological observables and timescales, enabling direct simulation for the mechanistic analysis of co-translational integration and for the engineering of membrane proteins with enhanced membrane integration efficiency.

Phosphatidylserine Stimulates Ceramide 1-Phosphate (C1P) Intermembrane Transfer by C1P Transfer Proteins

Genetic models for studying localized cell suicide that halt the spread of pathogen infection and immune response activation in plants include Arabidopsis accelerated-cell-death 11 mutant (acd11). In this mutant, sphingolipid homeostasis is disrupted via depletion of ACD11, a lipid transfer protein that is specific for ceramide 1-phosphate (C1P) and phyto-C1P. The C1P binding site in ACD11 and in human ceramide-1-phosphate transfer protein (CPTP) is surrounded by cationic residues. Here, we investigated the functional regulation of ACD11 and CPTP by anionic phosphoglycerides and found that 1-palmitoyl-2-oleoyl-phosphatidic acid or 1-palmitoyl-2-oleoyl-phosphatidylglycerol (≤15 mol %) in C1P source vesicles depressed C1P intermembrane transfer. By contrast, replacement with 1-palmitoyl-2-oleoyl-phosphatidylserine stimulated C1P transfer by ACD11 and CPTP. Notably, "soluble" phosphatidylserine (dihexanoyl-phosphatidylserine) failed to stimulate C1P transfer. Also, none of the anionic phosphoglycerides affected transfer action by human glycolipid lipid transfer protein (GLTP), which is glycolipid-specific and has few cationic residues near its glycolipid binding site. These findings provide the first evidence for a potential phosphoglyceride headgroup-specific regulatory interaction site(s) existing on the surface of any GLTP-fold and delineate new differences between GLTP superfamily members that are specific for C1P versus glycolipid.

Hydrophobic Clusters Raise the Threshold Hydrophilicity for Insertion of Transmembrane Sequences in Vivo

Insertion of a nascent membrane protein segment by the translocon channel into the bilayer is naturally promoted by high segmental hydrophobicity, but its selection as a transmembrane (TM) segment is complicated by the diverse environments (aqueous vs lipidic) the protein encounters and by the fact that most TM segments contain a substantial amount (∼30%) of polar residues, as required for protein structural stabilization and/or function. To examine the contributions of these factors systematically, we designed and synthesized a peptide library consisting of pairs of compositionally identical, but sequentially different, peptides with 19-residue core sequences varying (i) in Leu positioning (with five or seven Leu residues clustered into a contiguous "block" in the middle of the segment or "scrambled" throughout the sequence) and (ii) in Ser content (0-6 residues). The library was analyzed by a combination of biophysical and biological techniques, including HPLC retention times, circular dichroism measurements of helicity in micelle and phospholipid bilayer media, and relative blue shifts in Trp fluorescence maxima, as well as by the extent of membrane insertion in a translocon-mediated assay using microsomal membranes from dog pancreas endoplasmic reticulum. We found that local blocks of high hydrophobicity heighten the translocon's propensity to insert moderately hydrophilic sequences, until a "threshold hydrophilicity" is surpassed whereby segments no longer insert even in the presence of Leu blocks. This study codifies the prerequisites of apolar/polar content and residue positioning that define nascent TM segments, illustrates the accuracy in their prediction, and highlights how a single disease-causing mutation can tip the balance toward anomalous translocation/insertion.

Enthalpic and Entropic Contributions to Hydrophobicity

Hydrophobic hydration plays a key role in a vast variety of biological processes, ranging from the formation of cells to protein folding and ligand binding. Hydrophobicity scales simplify the complex process of hydration by assigning a value describing the averaged hydrophobic character to each amino acid. Previously published scales were not able to calculate the enthalpic and entropic contributions to the hydrophobicity directly. We present a new method, based on Molecular Dynamics simulations and Grid Inhomogeneous Solvation Theory, that calculates hydrophobicity from enthalpic and entropic contributions. Instead of deriving these quantities from the temperature dependence of the free energy of hydration or as residual of the free energy and the enthalpy, we directly obtain these values from the phase space occupied by water molecules. Additionally, our method is able to identify regions with specific enthalpic and entropic properties, allowing to identify so-called "unhappy water" molecules, which are characterized by weak enthalpic interactions and unfavorable entropic constraints.

Decrypting protein insertion through the translocon with free-energy calculations

Protein insertion into a membrane is a complex process involving numerous players. The most prominent of these players is the Sec translocon complex, a conserved protein-conducting channel present in the cytoplasmic membrane of bacteria and the membrane of the endoplasmic reticulum in eukaryotes. The last decade has seen tremendous leaps forward in our understanding of how insertion is managed by the translocon and its partners, coming from atomic-detailed structures, innovative experiments, and well-designed simulations. In this review, we discuss how experiments and simulations, hand-in-hand, teased out the secrets of the translocon-facilitated membrane insertion process. In particular, we focus on the role of free-energy calculations in elucidating membrane insertion. Amazingly, despite all its apparent complexity, protein insertion into membranes is primarily driven by simple thermodynamic and kinetic principles. This article is part of a Special Issue entitled: Membrane proteins edited by J.C. Gumbart and Sergei Noskov.

Integration of transmembrane domains is regulated by their downstream sequences

The Sec61 translocon catalyzes translocation of proteins into the endoplasmic reticulum and the lateral integration of transmembrane segments into the lipid bilayer. Integration is mediated by the hydrophobicity of a polypeptide segment consistent with thermodynamic equilibration between the translocon and the lipid membrane. Integration efficiency of a generic series of increasingly hydrophobic sequences (H-segments) was found to diverge significantly in different reporter constructs as a function of the ∼100 residues carboxyterminal of the H-segments. The hydrophobicity threshold of integration was considerably lowered by insertion of generic ∼20-residue peptides either made of flexible glycine-serine repeats, containing multiple negative charges, or consisting of an oligo-proline stretch. A highly flexible, 100-residue glycine-serine stretch maximally enhanced this effect. The apparent free energy of integration was found to be changed by more than 3 kcal/mol with the downstream sequences tested. The C-terminal sequences could also be shown to affect integration of natural mildly hydrophobic sequences. The results suggest that the conformation of the nascent polypeptide in the protected cavity between ribosome and translocon significantly influences the release of the H-segment into the bilayer.

Cholesterol enhances surface water diffusion of phospholipid bilayers

Elucidating the physical effect of cholesterol (Chol) on biological membranes is necessary towards rationalizing their structural and functional role in cell membranes. One of the debated questions is the role of hydration water in Chol-embedding lipid membranes, for which only little direct experimental data are available. Here, we study the hydration dynamics in a series of Chol-rich and depleted bilayer systems using an approach termed (1)H Overhauser dynamic nuclear polarization (ODNP) NMR relaxometry that enables the sensitive and selective determination of water diffusion within 5-10 Å of a nitroxide-based spin label, positioned off the surface of the polar headgroups or within the nonpolar core of lipid membranes. The Chol-rich membrane systems were prepared from mixtures of Chol, dipalmitoyl phosphatidylcholine and/or dioctadecyl phosphatidylcholine lipid that are known to form liquid-ordered, raft-like, domains. Our data reveal that the translational diffusion of local water on the surface and within the hydrocarbon volume of the bilayer is significantly altered, but in opposite directions: accelerated on the membrane surface and dramatically slowed in the bilayer interior with increasing Chol content. Electron paramagnetic resonance (EPR) lineshape analysis shows looser packing of lipid headgroups and concurrently tighter packing in the bilayer core with increasing Chol content, with the effects peaking at lipid compositions reported to form lipid rafts. The complementary capability of ODNP and EPR to site-specifically probe the hydration dynamics and lipid ordering in lipid membrane systems extends the current understanding of how Chol may regulate biological processes. One possible role of Chol is the facilitation of interactions between biological constituents and the lipid membrane through the weakening or disruption of strong hydrogen-bond networks of the surface hydration layers that otherwise exert stronger repulsive forces, as reflected in faster surface water diffusivity. Another is the concurrent tightening of lipid packing that reduces passive, possibly unwanted, diffusion of ions and water across the bilayer.

Plant lipid environment and membrane enzymes: the case of the plasma membrane H+-ATPase

Several lipid classes constitute the universal matrix of the biological membranes. With their amphipathic nature, lipids not only build the continuous barrier that confers identity to every cell and organelle, but they are also active actors that modulate the activity of the proteins immersed in the lipid bilayer. The plasma membrane H(+)-ATPase, an enzyme from plant cells, is an excellent example of a transmembrane protein whose activity is influenced by the hydrophilic compartments at both sides of the membrane and by the hydrophobic domains of the lipid bilayer. As a result, an extensive documentation of the effect of numerous amphiphiles in the enzyme activity can be found. Detergents, membrane glycerolipids, and sterols can produce activation or inhibition of the enzyme activity. In some cases, these effects are associated with the lipids of the membrane bulk, but in others, a direct interaction of the lipid with the protein is involved. This review gives an account of reports related to the action of the membrane lipids on the H(+)-ATPase activity.

Hydrophobic blocks facilitate lipid compatibility and translocon recognition of transmembrane protein sequences

Biophysical hydrophobicity scales suggest that partitioning of a protein segment from an aqueous phase into a membrane is governed by its perceived segmental hydrophobicity but do not establish specifically (i) how the segment is identified in vivo for translocon-mediated insertion or (ii) whether the destination lipid bilayer is biochemically receptive to the inserted sequence. To examine the congruence between these dual requirements, we designed and synthesized a library of Lys-tagged peptides of a core length sufficient to span a bilayer but with varying patterns of sequence, each composed of nine Leu residues, nine Ser residues, and one (central) Trp residue. We found that peptides containing contiguous Leu residues (Leu-block peptides, e.g., LLLLLLLLLWSSSSSSSSS), in comparison to those containing discontinuous stretches of Leu residues (non-Leu-block peptides, e.g., SLSLLSLSSWSLLSLSLLS), displayed greater helicity (circular dichroism spectroscopy), traveled slower during sodium dodecyl sulfate–polyacrylamide gel electrophoresis, had longer reverse phase high-performance liquid chromatography retention times on a C-18 column, and were helical when reconstituted into 1-palmitoyl-2-oleoylglycero-3-phosphocholine liposomes, each observation indicating superior lipid compatibility when a Leu-block is present. These parameters were largely paralleled in a biological membrane insertion assay using microsomal membranes from dog pancreas endoplasmic reticulum, where we found only the Leu-block sequences successfully inserted; intriguingly, an amphipathic peptide (SLLSSLLSSWLLSSLLSSL; Leu face, Ser face) with biophysical properties similar to those of Leu-block peptides failed to insert. Our overall results identify local sequence lipid compatibility rather than average hydrophobicity as a principal determinant of transmembrane segment potential, while demonstrating that further subtleties of hydrophobic and helical patterning, such as circumferential hydrophobicity in Leu-block segments, promote translocon-mediated insertion.

Peptides with the same composition, hydrophobicity, and hydrophobic moment bind to phospholipid bilayers with different affinities

We investigated the dependence of membrane binding on amino acid sequence for a series of amphipathic peptides derived from δ-lysin. δ-Lysin is a 26 amino acid, N-terminally formylated, hemolytic peptide that forms an amphipathic α-helix bound at membrane–water interfaces. A shortened peptide, lysette, was derived from δ-lysin by deletion of the four N-terminal amino acid residues. Five variants of lysette were synthesized by altering the amino acid sequence such that the overall hydrophobic moment remained essentially the same for all peptides. Peptide–lipid equilibrium dissociation constants and helicities of peptides bound to zwitterionic lipid vesicles were determined by stopped-flow fluorescence and circular dichroism. We found that binding to phosphatidylcholine bilayers was a function of the helicity of the bound peptide alone and independent of the a priori hydrophobic moment or the ability to form intramolecular salt bridges. Molecular dynamics (MD) simulations on two of the peptides suggest that sequence determines the insertion depth into the bilayer. The location of the two aspartate residues at the C-terminus of lysette-2 leads to a loss of helical content in the simulations, which correlates with faster desorption from the bilayer as compared to lysette. We also found a systematic deviation of the experimentally determined dissociation constant and that predicted by the Wimley–White interfacial hydrophobicity scale. The reason for the discrepancy remains unresolved but appears to correlate with a predominance of isoleucine over leucine residues in the lysette family of peptides.

Atomistic simulations of pore formation and closure in lipid bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.

The importance of membrane defects-lessons from simulations

The defects and pores within lipid membranes are scientifically interesting and have a number of biological applications. Although lipid bilayers are extremely thin hydrophobic barriers, just ∼3 nm thick, they include diverse chemistry and have complex structures. Bilayers are soft and dynamic, and as a result, they can bend and deform in response to different stimuli by means of structural changes in their component lipids. Though defects occur within these structures, their transience and small size have made it difficult to characterize them. However, with recent advances in computer power and computational modeling techniques, researchers can now use simulations as a powerful tool to probe the mechanism and energies of defect and pore formation in a number of situations. In this Account, we present results from our detailed molecular dynamics computer simulations of hydrophilic pores and related defects in lipid bilayers at an atomistic level. Electroporation can be used to increase the permeability of cellular membranes, with potential therapeutic applications. Atomistic simulations of electroporation have illustrated the molecular details of this process, including the importance of water dipole interactions at the water-membrane interface. Characterization of the lipid-protein interactions provides an important tool for understanding transmembrane protein structure and thermodynamic stability. Atomistic simulations give a detailed picture of the free energies of model peptides and side chains in lipid membranes; the energetic cost of defect formation strongly influences the energies of interactions between lipids and polar and charged residues. Many antimicrobial peptides form hydrophilic pores in lipid membranes, killing bacteria or cancer cells. On the basis of simulation data, at least some of these peptides form defects and pores near the center of the bilayer, with a common disordered structure where hydrated headgroups form an approximately toroidal shape. The localization and trafficking of lipids supports general membrane structure and a number of important signaling cascades, such as those involving ceramide, diacylglycerol, and cholesterol. Atomistic simulations have determined the rates and free energies of lipid flip-flop. During the flip-flop of most phosphatidylcholine lipids, a hydrophilic pore forms when the headgroup moves near the center of the bilayer. Simulations have provided novel insight into many features of defects and pores in lipid membranes. Simulation data from very different systems and models show how water penetration and defect formation can determine the free energies of many membrane processes. Bilayers can deform and allow transient defects and pores when exposed to a diverse range of stimuli. Future work will explore many aspects of membrane defects with increased resolution and scope, including the study of more complex lipid mixtures, membrane domains, and large-scale membrane remodeling. Such studies will examine processes including vesicle budding and fusion, non-bilayer lipid phases, and interactions between lipid bilayers and other biomolecules. Simulations provide information that complements experimental studies, allowing microscopic insight into experimental observations and suggesting novel hypotheses and experiments. These studies should enable a deeper understanding of the role of lipid bilayers in cellular biology and support the development of future lipid-based biotechnology.

Antibacterial surface treatment for orthopaedic implants

It is expected that the projected increased usage of implantable devices in medicine will result in a natural rise in the number of infections related to these cases. Some patients are unable to autonomously prevent formation of biofilm on implant surfaces. Suppression of the local peri-implant immune response is an important contributory factor. Substantial avascular scar tissue encountered during revision joint replacement surgery places these cases at an especially high risk of periprosthetic joint infection. A critical pathogenic event in the process of biofilm formation is bacterial adhesion. Prevention of biomaterial-associated infections should be concurrently focused on at least two targets: inhibition of biofilm formation and minimizing local immune response suppression. Current knowledge of antimicrobial surface treatments suitable for prevention of prosthetic joint infection is reviewed. Several surface treatment modalities have been proposed. Minimizing bacterial adhesion, biofilm formation inhibition, and bactericidal approaches are discussed. The ultimate anti-infective surface should be “smart” and responsive to even the lowest bacterial load. While research in this field is promising, there appears to be a great discrepancy between proposed and clinically implemented strategies, and there is urgent need for translational science focusing on this topic.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Membrane-active peptides: binding, translocation, and flux in lipid vesicles

Recently, new and improved methods have been developed to measure translocation of membrane-active peptides (antimicrobial, cytolytic, and amphipathic cell-penetrating peptides) across lipid bilayer membranes. The hypothesis that translocation of membrane-active peptides across a lipid bilayer is determined by the Gibbs energy of insertion of the peptide into the bilayer is re-examined in the light of new experimental tests. The original hypothesis and its motivation are first revisited, examining some of the specific predictions that it generated, followed by the results of the initial tests. Translocation is understood as requiring two previous steps: binding and insertion in the membrane. The problem of peptide binding to membranes, its prediction, measurement, and calculation are addressed. Particular attention is given to understanding the reason for the need for amphipathic structures in the function of membrane-active peptides. Insertion into the membrane is then examined. Hydrophobicity scales are compared, and their influence on calculations is discussed. The relation between translocation and graded or all-or-none peptide-induced flux from or into lipid vesicles is also considered. Finally, the most recent work on translocation is examined, both experimental and from molecular dynamics simulations. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex biomolecular systems. The generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (T-REMD), is one of the most powerful methods to explore structures of biomolecules such as proteins, nucleic acids, carbohydrates, and also of lipid membranes. T-REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, and surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We test the method on two biological membrane systems: a fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer and a WALP23-POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane system. During these simulations, a random walk in surface tension space is realized. Large-scale lateral deformation (shrinking and stretching) of the membranes takes place in all of the replicas without collapse of the lipid bilayer structure. There is accelerated lateral diffusion of DPPC lipid molecules compared with conventional MD simulation, and a much wider range of tilt angle of the WALP23 peptide is sampled due to large deformation of the POPC lipid bilayer and through peptide-lipid interactions. Our method could be applicable to a wide variety of biological membrane systems.

Functional asymmetry within the Sec61p translocon

The Sec61 translocon forms a pore to translocate polypeptide sequences across the membrane and offers a lateral gate for membrane integration of hydrophobic (H) segments. A central constriction of six apolar residues has been shown to form a seal, but also to determine the hydrophobicity threshold for membrane integration: Mutation of these residues in yeast Sec61p to glycines, serines, aspartates, or lysines lowered the hydrophobicity required for integration; mutation to alanines increased it. Whereas four leucines distributed in an oligo-alanine H segment were sufficient for 50% integration, we now find four leucines in the N-terminal half of the H segment to produce significantly more integration than in the C-terminal half, suggesting functional asymmetry within the translocon. Scanning a cluster of three leucines through an oligo-alanine H segment showed high integration levels, except around the position matching that of the hydrophobic constriction in the pore where integration was strongly reduced. Both asymmetry and the position effect of H-segment integration disappeared upon mutation of the constriction residues to glycines or serines, demonstrating that hydrophobicity at this position within the translocon is responsible for the phenomenon. Asymmetry was largely retained, however, when constriction residues were replaced by alanines. These results reflect on the integration mechanism of transmembrane domains and show that membrane insertion of H segments strongly depends not only on their intrinsic hydrophobicity but also on the local conditions in the translocon interior. Thus, the contribution of hydrophobic residues in the H segment is not simply additive and displays cooperativeness depending on their relative position.