Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques

Ion-driven rotary membrane motors: From structure to function

Ion-driven membrane motors, essential across all domains of life, convert a gradient of ions across a membrane into rotational energy, facilitating diverse biological processes including ATP synthesis, substrate transport, and bacterial locomotion. Herein, we highlight recent structural advances in the understanding of two classes of ion-driven membrane motors: rotary ATPases and 5:2 motors. The recent structure of the human F-type ATP synthase is emphasised along with the gained structural insight into clinically relevant mutations. Furthermore, we highlight the diverse roles of 5:2 motors and recent mechanistic understanding gained through the resolution of ions in the structure of a sodium-driven motor, combining insights into potential unifying mechanisms of ion selectivity and rotational torque generation in the context of their function as part of complex biological systems.

Exploring pathological link between antimicrobial and amyloid peptides

Amyloid peptides (AMYs) and antimicrobial peptides (AMPs) are considered as the two distinct families of peptides, characterized by their unique sequences, structures, biological functions, and specific pathological targets. However, accumulating evidence has revealed intriguing pathological connections between these peptide families in the context of microbial infection and neurodegenerative diseases. Some AMYs and AMPs share certain structural and functional characteristics, including the ability to self-assemble, the presence of β-sheet-rich structures, and membrane-disrupting mechanisms. These shared features enable AMYs to possess antimicrobial activity and AMPs to acquire amyloidogenic properties. Despite limited studies on AMYs-AMPs systems, the cross-seeding phenomenon between AMYs and AMPs has emerged as a crucial factor in the bidirectional communication between the pathogenesis of neurodegenerative diseases and host defense against microbial infections. In this review, we examine recent developments in the potential interplay between AMYs and AMPs, as well as their pathological implications for both infectious and neurodegenerative diseases. By discussing the current progress and challenges in this emerging field, this account aims to inspire further research and investments to enhance our understanding of the intricate molecular crosstalk between AMYs and AMPs. This knowledge holds great promise for the development of innovative therapies to combat both microbial infections and neurodegenerative disorders.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Structure and Orientation of the SARS-Coronavirus-2 Spike Protein at Air-Water Interfaces

The SARS coronavirus 2 (SARS-CoV-2) spike protein is located at the outermost perimeter of the viral envelope and is the first component of the virus to make contact with surrounding interfaces. The stability of the spike protein when in contact with surfaces plays a deciding role for infection pathways and for the viability of the virus after surface contact. While cryo-EM structures of the spike protein have been solved with high resolution and structural studies in solution have provided information about the secondary and tertiary structures, only little is known about the folding when adsorbed to surfaces. We here report on the secondary structure and orientation of the S1 segment of the spike protein, which is often used as a model protein for in vitro studies of SARS-CoV-2, at the air-water interface using surface-sensitive vibrational sum-frequency generation (SFG) spectroscopy. The air-water interface plays an important role for SARS-CoV-2 when suspended in aerosol droplets, and it serves as a model system for hydrophobic surfaces in general. The SFG experiments show that the S1 segment of the spike protein remains folded at the air-water interface and predominantly binds in its monomeric state, while the combination of small-angle X-ray scattering and two-dimensional infrared spectroscopy measurements indicate that it forms hexamers with the same secondary structure in aqueous solution.

Membrane Activity of a DNA-Based Ion Channel Depends on the Stability of Its Double-Stranded Structure

DNA nanotechnology has emerged as a promising method for designing spontaneously inserting and fully controllable synthetic ion channels. However, both insertion efficiency and stability of existing DNA-based membrane channels leave much room for improvement. Here, we demonstrate an approach to overcoming the unfavorable DNA-lipid interactions that hinder the formation of a stable transmembrane pore. Our all-atom MD simulations and experiments show that the insertion-driving cholesterol modifications can cause fraying of terminal base pairs of nicked DNA constructs, distorting them when embedded in a lipid bilayer. Importantly, we show that DNA nanostructures with no backbone discontinuities form more stable conductive pores and insert into membranes with a higher efficiency than the equivalent nicked constructs. Moreover, lack of nicks allows design and maintenance of membrane-spanning helices in a tilted orientation within the lipid bilayer. Thus, reducing the conformational degrees of freedom of the DNA nanostructures enables better control over their function as synthetic ion channels.

Diverse Scientific Benchmarks for Implicit Membrane Energy Functions

Energy functions are fundamental to biomolecular modeling. Their success depends on robust physical formalisms, efficient optimization, and high-resolution data for training and validation. Over the past 20 years, progress in each area has advanced soluble protein energy functions. Yet, energy functions for membrane proteins lag behind due to sparse and low-quality data, leading to overfit tools. To overcome this challenge, we assembled a suite of 12 tests on independent data sets varying in size, diversity, and resolution. The tests probe an energy function's ability to capture membrane protein orientation, stability, sequence, and structure. Here, we present the tests and use the franklin2019 energy function to demonstrate them. We then identify areas for energy function improvement and discuss potential future integration with machine-learning-based optimization methods. The tests are available through the Rosetta Benchmark Server (https://benchmark.graylab.jhu.edu/) and GitHub (https://github.com/rfalford12/Implicit-Membrane-Energy-Function-Benchmark).

Prediction of amphipathic helix-membrane interactions with Rosetta

Amphipathic helices have hydrophobic and hydrophilic/charged residues situated on opposite faces of the helix. They can anchor peripheral membrane proteins to the membrane, be attached to integral membrane proteins, or exist as independent peptides. Despite the widespread presence of membrane-interacting amphipathic helices, there is no computational tool within Rosetta to model their interactions with membranes. In order to address this need, we developed the AmphiScan protocol with PyRosetta, which runs a grid search to find the most favorable position of an amphipathic helix with respect to the membrane. The performance of the algorithm was tested in benchmarks with the RosettaMembrane, ref2015_memb, and franklin2019 score functions on six engineered and 44 naturally-occurring amphipathic helices using membrane coordinates from the OPM and PDBTM databases, OREMPRO server, and MD simulations for comparison. The AmphiScan protocol predicted the coordinates of amphipathic helices within less than 3Å of the reference structures and identified membrane-embedded residues with a Matthews Correlation Constant (MCC) of up to 0.57. Overall, AmphiScan stands as fast, accurate, and highly-customizable protocol that can be pipelined with other Rosetta and Python applications.

Amphipathic helices are important targets as antibacterial peptides and as domains of membrane proteins that play a role in sensing the membrane environment. Understanding how amphipathic helices interact with membrane enables us to design better peptides and understand how membrane proteins use them to interact with their environment. However, there is a limited number of tools available for the modeling of amphipathic helices in membranes. Implicit membrane models can be used for this purpose as simplistic representations of the membrane environment. In this work, we developed the AmphiScan protocol that can be used to predict membrane coordinates of amphipathic helices starting with a helix structure in an implicit membrane environment. We benchmarked the performance of AmphiScan on engineered LK peptides, naturally-occurring amphipathic helices, and hydrophobic and hydrophilic peptides. Our approach provides a reliable and customizable tool to model amphipathic helix–membrane interactions, and pose a platform for the screening of amphipathic helix properties in silico.

Protein Structure Prediction and Design in a Biologically Realistic Implicit Membrane

Protein design is a powerful tool for elucidating mechanisms of function and engineering new therapeutics and nanotechnologies. Although soluble protein design has advanced, membrane protein design remains challenging because of difficulties in modeling the lipid bilayer. In this work, we developed an implicit approach that captures the anisotropic structure, shape of water-filled pores, and nanoscale dimensions of membranes with different lipid compositions. The model improves performance in computational benchmarks against experimental targets, including prediction of protein orientations in the bilayer, ΔΔG calculations, native structure discrimination, and native sequence recovery. When applied to de novo protein design, this approach designs sequences with an amino acid distribution near the native amino acid distribution in membrane proteins, overcoming a critical flaw in previous membrane models that were prone to generating leucine-rich designs. Furthermore, the proteins designed in the new membrane model exhibit native-like features including interfacial aromatic side chains, hydrophobic lengths compatible with bilayer thickness, and polar pores. Our method advances high-resolution membrane protein structure prediction and design toward tackling key biological questions and engineering challenges.

Chirality-Dependent Adsorption between Amphipathic Peptide and POPC Membrane

The interactions between chiral molecules and cell membranes have attracted more and more attention in recent decades, due to their importance in molecular science and medical applications. It is observed that some peptides composed of different chiral amino acids may have distinct interactions with a membrane. How does the membrane exhibit a selective behavior related to the chirality of the peptides? Microscopically, the interactions between the peptides and the membrane are poorly understood. In this work, we study the interactions between an amphipathic peptide (C6) and POPC membrane with simulations. The kinetics and thermodynamics of peptide enantiomers during the adsorption to the membrane are characterized with direct simulations and umbrella sampling. It is observed that there are slow kinetics for the peptide composed of D-type amino acids. Along the observed pathways, the free energy landscapes are determined with umbrella sampling techniques. A free-energy barrier for the peptide composed of D-amino acids is observed, which is consistent with the kinetic observations. The results indicate the concurrent adsorption and rotation of the peptide helix. The local interactions between the peptides and the membrane are examined in detail, including the contact interactions between the peptides and the membrane, and the distributions of the lipids around the peptide. There are observable differences of the local interactions for the cases related to different peptide enantiomers. These results further demonstrate the importance of the rotation of peptide helix during the adsorption. More interestingly, all these kinetic differences between peptide enantiomers can be explained based on the conformations of the residue Trp and interactions between Trp and lipid molecules. These results give us a molecular understanding of the mechanism of the chirality-dependent peptide-membrane interactions, and may provide clues to designing systems which are sensitive to the chirality of membranes.

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

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

Engineering monolayer poration for rapid exfoliation of microbial membranes

A novel mechanism of monolayer poration leading to the rapid exfoliation and lysis of microbial membranes is reported.

The spread of bacterial resistance to traditional antibiotics continues to stimulate the search for alternative antimicrobial strategies. All forms of life, from bacteria to humans, are postulated to rely on a fundamental host defense mechanism, which exploits the formation of open pores in microbial phospholipid bilayers. Here we predict that transmembrane poration is not necessary for antimicrobial activity and reveal a distinct poration mechanism that targets the outer leaflet of phospholipid bilayers. Using a combination of molecular-scale and real-time imaging, spectroscopy and spectrometry approaches, we introduce a structural motif with a universal insertion mode in reconstituted membranes and live bacteria. We demonstrate that this motif rapidly assembles into monolayer pits that coalesce during progressive membrane exfoliation, leading to bacterial cell death within minutes. The findings offer a new physical basis for designing effective antibiotics.

Folding Simulations of the Transmembrane Helix of Virus Protein U in an Implicit Membrane Model

Vpu is an 81-amino-acid auxiliary membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). One of its roles is to amplify viral release by self-assembling in homo-oligomers to form functional water-filled pores enabling the flux of ions across the membrane. Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation. With a recently developed implicit membrane model and an efficient Monte Carlo (MC) algorithm using concerted backbone rotations, we simulate the folding of the transmembrane domain of Vpu at atomic resolution. The implicit membrane environment is based on the generalized Born theory and enables very long time scale events, such as folding to be observed using detailed all-atom representation of the protein. Such studies are currently computationally unfeasible with fully explicit lipid bilayer molecular dynamics simulations. The correct helical transmembrane structure of Vpu is predicted from extended conformations and remains stably inserted. Tilt and kink angles agree well with experimental estimates from NMR measurements. The experimentally observed change in tilt angle in membranes of varying hydrophobic width is accurately reproduced. The extensive simulation of a pentamer of the Vpu transmembrane domain in the implicit membrane gives results similar to the ones reported previously for fully explicit bilayer simulations.

Evaluating membrane affinity by integrating protein orientations

Energetic interactions of a protein with lipid bilayers determine its propensity to reside in the membrane. Here we seek to evaluate the membrane interactions for EMAPII, a protein found to be released from the cell by unknown mechanisms, as well as several other proteins. Using a knowledge-based coarse-grained membrane potential, we calculate the free energy profiles for these proteins by integrating out the orientation degrees of freedom. Due to the invariance of energy under in-plane rotations about the membrane normal, the orientation space can be reduced to two dimensions and mapped onto the surface of a unit sphere, thus making visualization, sampling and integration more convenient. The integrated free energy profiles determine the relative probabilities along the membrane normal for the proteins regardless of their orientations, and display distinctive characteristics for membrane proteins and water-soluble proteins. The membrane interactions for EMAPII exhibit typical features of a water-soluble protein, with a high energetic barrier to enter or cross the membrane. Our results thus suggest that similar to the non-classical export of FGF1, the release of EMAPII would involve more complicated mechanisms than simple passive diffusion across the membrane.

Molecular understanding of a potential functional link between antimicrobial and amyloid peptides

Antimicrobial and amyloid peptides do not share common sequences, typical secondary structures, or normal biological activity but both the classes of peptides exhibit membrane-disruption ability to induce cell toxicity. Different membrane-disruption mechanisms have been proposed for antimicrobial and amyloid peptides, individually, some of which are not exclusive to either peptide type, implying that certain common principles may govern the folding and functions of different cytolytic peptides and associated membrane disruption mechanisms. Particularly, some antimicrobial and amyloid peptides have been identified to have dual complementary amyloid and antimicrobial properties, suggesting a potential functional link between amyloid and antimicrobial peptides. Given that some similar structural and membrane-disruption characteristics exist between the two classes of peptides, this review summarizes major findings, recent advances, and future challenges related to antimicrobial and amyloid peptides and strives to illustrate the similarities, differences, and relationships in the sequences, structures, and membrane interaction modes between amyloid and antimicrobial peptides, with a special focus on direct interactions of the peptides with the membranes. We hope that this review will stimulate further research at the interface of antimicrobial and amyloid peptides - which has been studied less intensively than either type of peptides - to decipher a possible link between both amyloid pathology and antimicrobial activity, which can guide drug design and peptide engineering to influence peptide-membrane interactions important in human health and diseases.

Transferring the PRIMO Coarse-Grained Force Field to the Membrane Environment: Simulations of Membrane Proteins and Helix-Helix Association

An extension of the recently developed PRIMO coarse-grained force field to membrane environments, PRIMO-M, is described. The membrane environment is modeled with the heterogeneous dielectric generalized Born (HDGB) methodology that simply replaces the standard generalized Born model in PRIMO without further parametrization. The resulting model was validated by comparing amino acid insertion free energy profiles and application in molecular dynamics simulations of membrane proteins and membrane-interacting peptides. Membrane proteins with 148–661 amino acids show stable root-mean-squared-deviations (RMSD) between 2 and 4 Å for most systems. Transmembrane helical peptides maintain helical shape and exhibit tilt angles in good agreement with experimental or other simulation data. The association of two glycophorin A (GpA) helices was simulated using replica exchange molecular dynamics simulations yielding the correct dimer structure with a crossing angle in agreement with previous studies. Finally, conformational sampling of the influenza fusion peptide also generates structures in agreement with previous studies. Overall, these findings suggest that PRIMO-M can be used to study membrane bound peptides and proteins and validates the transferable nature of the PRIMO coarse-grained force field.

Biological applications of classical electrostatics methods

Continuum electrostatics modeling of solvation based on the Poisson–Boltzmann (PB) equation has gained wide acceptance in biomolecular applications such as energetic analysis and structural visualization. Successful application of the PB solvent models requires careful calibration of the solvation parameters. Extensive testing and validation is also important to ensure accuracy in their applications. Limitation in the continuum modeling of solvation is also a known issue in certain biomolecular applications. Growing interest in membrane systems has further spurred developmental efforts to allow inclusion of membrane in the PB solvent models. Despite their past successes due to careful parameterization, algorithm development and parallel implementation, there is still much to be done to improve their transferability from the small molecular systems upon which they were developed and validated to complex macromolecular systems as advances in technology continue to push forward, providing ever greater computational resources to researchers to study more interesting biological systems of higher complexity.

More than the sum of its parts: coarse-grained peptide-lipid interactions from a simple cross-parametrization

Interfacial systems are at the core of fascinating phenomena in many disciplines, such as biochemistry, soft-matter physics, and food science. However, the parametrization of accurate, reliable, and consistent coarse-grained (CG) models for systems at interfaces remains a challenging endeavor. In the present work, we explore to what extent two independently developed solvent-free CG models of peptides and lipids--of different mapping schemes, parametrization methods, target functions, and validation criteria--can be combined by only tuning the cross-interactions. Our results show that the cross-parametrization can reproduce a number of structural properties of membrane peptides (for example, tilt and hydrophobic mismatch), in agreement with existing peptide-lipid CG force fields. We find encouraging results for two challenging biophysical problems: (i) membrane pore formation mediated by the cooperative action of several antimicrobial peptides, and (ii) the insertion and folding of the helix-forming peptide WALP23 in the membrane.

Probing the orientational distribution of dyes in membranes through multiphoton microscopy

Numerous dyes are available or under development for probing the structural and functional properties of biological membranes. Exogenous chromophores adopt a range of orientations when bound to membranes, which have a drastic effect on their biophysical behavior. Here, we present a method that employs optical anisotropy data from three polarization-imaging techniques to establish the distribution of orientations adopted by molecules in monolayers and bilayers. The resulting probability density functions, which contain the preferred molecular tilt μ and distribution breadth γ, are more informative than an average tilt angle [φ]. We describe a methodology for the extraction of anisotropy data through an image-processing technology that decreases the error in polarization measurements by about a factor of four. We use this technique to compare di-4-ANEPPS and di-8-ANEPPS, both dipolar dyes, using data from polarized 1-photon, 2-photon fluorescence and second-harmonic generation imaging. We find that di-8-ANEPPS has a lower tilt but the same distributional width. We find the distribution of tilts taken by di-4-ANEPPS in two phospholipid membrane models: giant unilamellar vesicles and water-in-oil droplet monolayers. Both models result in similar distribution functions with average tilts of 52° and 47°, respectively.

Multiscale molecular dynamics simulations of membrane proteins

The time and length scales accessible by biomolecular simulations continue to increase. This is in part due to improvements in algorithms and computing performance, but is also the result of the emergence of coarse-grained (CG) potentials, which complement and extend the information obtainable from fully detailed models. CG methods have already proven successful for a range of applications that benefit from the ability to rapidly simulate spontaneous self-assembly within a lipid membrane environment, including the insertion and/or oligomerization of a range of "toy models," transmembrane peptides, and single- and multi-domain proteins. While these simplified approaches sacrifice atomistic level detail, it is now straightforward to "reverse map" from CG to atomistic descriptions, providing a strategy to assemble membrane proteins within a lipid environment, prior to all-atom simulation. Moreover, recent developments have been made in "dual resolution" techniques, allowing different molecules in the system to be modeled with atomistic or CG resolution simultaneously.

Solvation models and computational prediction of orientations of peptides and proteins in membranes

Membrane-associated peptides and proteins function in the highly heterogeneous environment of the lipid bilayer whose physico-chemical properties change non-monotonically along the bilayer normal. To simulate insertion of peptides and proteins into membranes and correctly reproduce the energetics of this process, an adequate solvation model and physically realistic representation of the lipid bilayer should be employed. We present a brief overview of the existing solvation models and their application for prediction of binding affinities and orientations of proteins in membranes. Particular emphasis is placed on the recently proposed PPM method, the corresponding web server, and the OPM database that were designed for positioning in membranes of integral and peripheral proteins with known three-dimensional structures.

Knowledge-based potential for positioning membrane-associated structures and assessing residue-specific energetic contributions

The complex hydrophobic and hydrophilic milieus of membrane-associated proteins pose experimental and theoretical challenges to their understanding. Here, we produce a nonredundant database to compute knowledge-based asymmetric cross-membrane potentials from the per-residue distributions of C(β), C(γ) and functional group atoms. We predict transmembrane and peripherally associated regions from genomic sequence and position peptides and protein structures relative to the bilayer (available at http://www.degradolab.org/ez). The pseudo-energy topological landscapes underscore positional stability and functional mechanisms demonstrated here for antimicrobial peptides, transmembrane proteins, and viral fusion proteins. Moreover, experimental effects of point mutations on the relative ratio changes of dual-topology proteins are quantitatively reproduced. The functional group potential and the membrane-exposed residues display the largest energetic changes enabling to detect native-like structures from decoys. Hence, focusing on the uniqueness of membrane-associated proteins and peptides, we quantitatively parameterize their cross-membrane propensity, thus facilitating structural refinement, characterization, prediction, and design.

Beyond modeling: all-atom olfactory receptor model simulations

Olfactory receptors (ORs) are a type of GTP-binding protein-coupled receptor (GPCR). These receptors are responsible for mediating the sense of smell through their interaction with odor ligands. OR-odorant interactions marks the first step in the process that leads to olfaction. Computational studies on model OR structures can generate focused and novel hypotheses for further bench investigation by providing a view of these interactions at the molecular level beyond inferences that are drawn merely from static docking. Here we have shown the specific advantages of simulating the dynamic environment associated with OR-odorant interactions. We present a rigorous protocol which ranges from the creation of a computationally derived model of an olfactory receptor to simulating the interactions between an OR and an odorant molecule. Given the ubiquitous occurrence of GPCRs in the membranes of cells, we anticipate that our OR-developed methodology will serve as a model for the computational structural biology of all GPCRs.

Effects of clinically relevant MPL mutations in the transmembrane domain revealed at the atomic level through computational modeling

Background

Mutations in the thrombopoietin receptor (MPL) may activate relevant pathways and lead to chronic myeloproliferative neoplasms (MPNs). The mechanisms of MPL activation remain elusive because of a lack of experimental structures. Modern computational biology techniques were utilized to explore the mechanisms of MPL protein activation due to various mutations.

Results

Transmembrane (TM) domain predictions, homology modeling, ab initio protein structure prediction, and molecular dynamics (MD) simulations were used to build structural dynamic models of wild-type and four clinically observed mutants of MPL. The simulation results suggest that S505 and W515 are important in keeping the TM domain in its correct position within the membrane. Mutations at either of these two positions cause movement of the TM domain, altering the conformation of the nearby intracellular domain in unexpected ways, and may cause the unwanted constitutive activation of MPL's kinase partner, JAK2.

Conclusions

Our findings represent the first full-scale molecular dynamics simulations of the wild-type and clinically observed mutants of the MPL protein, a critical element of the MPL-JAK2-STAT signaling pathway. In contrast to usual explanations for the activation mechanism that are based on the relative translational movement between rigid domains of MPL, our results suggest that mutations within the TM region could result in conformational changes including tilt and rotation (azimuthal) angles along the membrane axis. Such changes may significantly alter the conformation of the adjacent and intrinsically flexible intracellular domain. Hence, caution should be exercised when interpreting experimental evidence based on rigid models of cytokine receptors or similar systems.

Determinants of membrane activity from mutational analysis of the HIV fusion peptide

We have synthesized a small library of 38 variants of the 23-residue fusion peptide domain found at the N-terminus of gp41 glycoprotein of HIV. This hydrophobic, glycine-rich sequence is critical for viral infectivity and is thought to be central in the membrane fusion of viral envelope with the host membrane. There has been extensive discussion regarding the origin of fusogenicity in this viral fusion sequence. Our library of fusion peptide variants was designed to address the biophysical importance of secondary structure, peptide flexibility, glycine content, and placement. We assayed each peptide for its ability to induce lipid mixing and membrane permeablization in synthetic vesicles. We find that the viral fusion peptide may be greatly simplified while retaining fusogenic function and minimizing membrane-permeablizing function; to the best of our knowledge, this is the first attempt to optimize fusogenic function of the HIV fusion peptide through sequence variation. Our data show that many flexible, linear, minimally hydrophobic peptides may achieve the biophysical function of fusion; glycine does not appear to be essential. These findings will be useful in the design of synthetic fusogens for cellular delivery.

A comparative study on the ability of two implicit solvent lipid models to predict transmembrane helix tilt angles

Free-energy profiles describing the relative orientation of membrane proteins along predefined coordinates can be efficiently calculated by means of umbrella simulations. Such simulations generate reliable orientational distributions but are difficult to converge because of the very long equilibration times of the solvent and the lipid bilayer in explicit representation. Two implicit lipid membrane models are here applied in combination with the umbrella sampling strategy to the simulation of the transmembrane (TM) helical segment from virus protein U (Vpu). The models are used to study both orientation and energetics of this α-helical peptide as a function of hydrophobic mismatch. We observe that increasing the degree of positive hydrophobic mismatch increased the tilt angle of Vpu. These findings agree well with experimental data and as such validate the solvation models used in this study.

[Prediction of the spatial structure of proteins: emphasis on membrane targets]

Knowledge of the spatial structure of proteins is a prerequisite for both awareness of their functional mechanisms and the framework for rational drug discovery and design. Meanwhile, direct structural determination is often hampered or impractical due to the complexity, expensiveness, and limited capabilities of experimental techniques. These issues are especially pronounced for integral membrane proteins. On numerous occasions, the theoretical prediction of protein structures may facilitate the process by exploiting physical or empirical principles. This paper surveys modern techniques for the prediction of the spatial structure of proteins using computer algorithms, and the main emphasis is placed on the most "complex" targets - membrane proteins (MPs). The first part of the review describes de novo methods based on empirical physical principles; in the second part, a comparative modeling philosophy, which accounts for the structure of related proteins, is described. Special focus is made regarding pharmacologically relevant classes of G-coupled receptors, receptor tyrosine ki-nases, and other MPs. Algorithms for the assessment of the models quality and potential fields of application of computer models are discussed.

Construction of an implicit membrane environment for the lattice Monte Carlo simulation of transmembrane protein

Due to the complexity of biological membrane, computer simulation of transmembrane protein's folding is challenging. In this paper, an implicit biological membrane environment has been constructed in lattice space, in which the lipid chains and water molecules were represented by the unoccupied lattice sites. The biological membrane was characterized with three features: stronger hydrogen bonding interaction, membrane lateral pressure, and lipophobicity index for the amino acid residues. In addition to the hydrocarbon core spanning region and the water solution, the lipid interface has also been represented in this implicit membrane environment, which was proved to be effective for the transmembrane protein's folding. The associated Monte Carlo simulations have been performed for SARS-CoV E protein and M2 protein segment (residues 18–60) of influenza A virus. It was found that the coil–helix transition of the transmembrane segment occurred earlier than the coil–globule transition of the two terminal domains. The folding process and final orientation of the amphipathic helical block in water solution are obviously influenced by its corresponding hydrophobicity/lipophobicity. Therefore, this implicit membrane environment, though in lattice space, can make an elaborate balance between different driving forces for the membrane protein's folding, thus offering a potential means for the simulation of transmembrane protein oligomers in feasible time.

Viruses: incredible nanomachines. New advances with filamentous phages

During recent decades, bacteriophages have been at the cutting edge of new developments in molecular biology, biophysics, and, more recently, bionanotechnology. In particular filamentous viruses, for example bacteriophage M13, have a virion architecture that enables precision building of ordered and defect-free two and three-dimensional structures on a nanometre scale. This could not have been possible without detailed knowledge of coat protein structure and dynamics during the virus reproduction cycle. The results of the spectroscopic studies conducted in our group compellingly demonstrate a critical role of membrane embedment of the protein both during infectious entry of the virus into the host cell and during assembly of the new virion in the host membrane. The protein is effectively embedded in the membrane by a strong C-terminal interfacial anchor, which together with a simple tilt mechanism and a subtle structural adjustment of the extreme end of its N terminus provides favourable thermodynamical association of the protein in the lipid bilayer. This basic physicochemical rule cannot be violated and any new bionanotechnology that will emerge from bacteriophage M13 should take this into account.

Peptide nanopores and lipid bilayers: interactions by coarse-grained molecular-dynamics simulations

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (alpha-helical and beta-barrel), and the seven different bilayer systems range in thickness from approximately 28 to approximately 43 A. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.

Inter-helical interactions in membrane proteins: analysis based on the local backbone geometry and the side chain interactions

The availability of a significant number of the structures of helical membrane proteins has prompted us to investigate the mode of helix-helix packing. In the present study, we have considered a dataset of alpha-helical membrane proteins representing structures solved from all the known superfamilies. We have described the geometry of all the helical residues in terms of local coordinate axis at the backbone level. Significant inter-helical interactions have been considered as contacts by weighing the number of atom-atom contacts, including all the side-chain atoms. Such a definition of local axis and the contact criterion has allowed us to investigate the inter-helical interaction in a systematic and quantitative manner. We show that a single parameter (designated as alpha), which is derived from the parameters representing the mutual orientation of local axes, is able to accurately capture the details of helix-helix interaction. The analysis has been carried out by dividing the dataset into parallel, anti-parallel, and perpendicular orientation of helices. The study indicates that a specific range of alpha value is preferred for interactions among the anti-parallel helices. Such a preference is also seen among interacting residues of parallel helices, however to a lesser extent. No such preference is seen in the case of perpendicular helices, the contacts that arise mainly due to the interaction of surface helices with the end of the trans-membrane helices. The study supports the prevailing view that the anti-parallel helices are well packed. However, the interactions between helices of parallel orientation are non-trivial. The packing in alpha-helical membrane proteins, which is systematically and rigorously investigated in this study, may prove to be useful in modeling of helical membrane proteins.

Surface activity of amphiphilic helical beta-peptides from molecular dynamics simulation

The surface activity of beta-peptides is investigated using molecular simulations. The type and display of hydrophobic and hydrophilic groups on helical beta-peptides is varied systematically. Peptides with 2/3 hydrophobic groups are found to be surface active, and to adopt an orientation parallel to the air-water interface. For select beta-peptides, we also determine the potential of mean force required to bring a peptide to the air-water interface. Facially amphiphilic helices with 2/3 hydrophobic groups are found to exhibit the lowest free energy of adsorption. The adsorption process is driven by a favorable energetic term and opposed by negative entropic changes. The temperature dependence of adsorption is also investigated; facially amphiphilic helices are found to adopt orientations that are largely independent of temperature, while nonfacially amphiphilic helices sample a broader range of interfacial orientations at elevated temperatures. The thermodynamics of adsorption of beta-peptides is compared to that of 1-octanol, a well-known surfactant, and ovispirin, a naturally occurring antimicrobial peptide. It is found that the essential difference lies in the sign of the entropy of adsorption, which is negative for beta- and alpha-peptides and positive for traditional surfactants such as octanol.

Free-energy profiles of membrane insertion of the M2 transmembrane peptide from influenza A virus

The insertion of the M2 transmembrane peptide from influenza A virus into a membrane has been studied with molecular-dynamics simulations. This system is modeled by an atomically detailed peptide interacting with a continuum representation of a membrane bilayer in aqueous solution. We performed replica-exchange molecular-dynamics simulations with umbrella-sampling techniques to characterize the probability distribution and conformation preference of the peptide in the solution, at the membrane interface, and in the membrane. The minimum in the calculated free-energy surface of peptide insertion corresponds to a fully inserted, helical peptide spanning the membrane. The free-energy profile also shows that there is a significant barrier for the peptide to enter into this minimum in a nonhelical conformation. The sequence of the peptide is such that hydrophilic amino acid residues at the ends of the otherwise primarily hydrophobic peptide create a trapped, U-shaped conformation with the hydrophilic residues associated with the aqueous phase and the hydrophobic residues embedded in the membrane. Analysis of the free energy shows that the barrier to insertion is largely enthalpic in nature, whereas the membrane-spanning global minimum is favored by entropy.

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of alpha-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 x 5 micros) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2x rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of approximately 16 degrees. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.

Coarse-grained MD simulations of membrane protein-bilayer self-assembly

Complete determination of a membrane protein structure requires knowledge of the protein position within the lipid bilayer. As the number of determined structures of membrane proteins increases so does the need for computational methods which predict their position in the lipid bilayer. Here we present a coarse-grained molecular dynamics approach to lipid bilayer self-assembly around membrane proteins. We demonstrate that this method can be used to predict accurately the protein position in the bilayer for membrane proteins with a range of different sizes and architectures.

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.

A generalized born implicit-membrane representation compared to experimental insertion free energies

An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.