The effect of membrane curvature on the conformation of antimicrobial peptides: implications for binding and the mechanism of action

Short cationic antimicrobial peptides (AMPs) are believed to act either by inducing transmembrane pores or disrupting membranes in a detergent-like manner. For example, the antimicrobial peptides aurein 1.2, citropin 1.1, maculatin 1.1 and caerin 1.1, despite being closely related, appear to act by fundamentally different mechanisms depending on their length. Using molecular dynamics simulations, the structural properties of these four peptides have been examined in solution as well as in a variety of membrane environments. It is shown that each of the peptides has a strong preference for binding to regions of high membrane curvature and that the structure of the peptides is dependent on the degree of local curvature. This suggests that the shorter peptides aurein 1.2 and citropin 1.1 act via a detergent-like mechanism because they can induce high local, but not long-range curvature, whereas the longer peptides maculatin 1.1 and caerin 1.1 require longer range curvature to fold and thus bind to and stabilize transmembrane pores

Full-length structural model of RET3 and SEC21 in COPI: identification of binding sites on the appendage for accessory protein recruitment motifs

COPI, a 600 kD heptameric complex (consisting of subunits α, β, γ, δ, ε, ζ, and β′) “coatomer,” assembles non-clathrin-coated vesicles and is responsible for intra-Golgi and Golgi-to-ER protein trafficking. Here, we report the three-dimensional structures of the entire sequences of yeast Sec21 (γ-COPI mammalian ortholog), yeast Ret3 (ζ-COPI mammalian ortholog), and the results of successive molecular dynamics investigations of the subunits and assembly based on a protein–protein docking experiment. The three-dimensional structures of the subunits in their complexes indicate the residues of the two subunits that impact on assembly, the conformations of Ret3 and Sec21, and their binding orientations in the complexed state. The structure of the appendage domain of Sec21, with its two subdomains—the platform and the β-sandwich, was investigated to explore its capacity to bind to accessory protein recruitment motifs. Our study shows that a binding site on the platform is capable of binding the Eps15 DPF and epsin DPW2 peptides, whereas the second site on the platform and the site on the β-sandwich subdomain were found to selectively bind to the amphiphysin FXDXF and epsin DPW1 peptides, respectively. Identifying the regions of both the platform and sandwich subdomains involved in binding each peptide motif clarifies the mechanism through which the appendage domain of Sec21 engages with the accessory proteins during the trafficking process of non-clathrin-coated vesicles

Hydrophilicity Matching – A Potential Prerequisite for the Formation of Protein-Protein Complexes in the Cell

A binding event between two proteins typically consists of a diffusional search of binding partners for one another, followed by a specific recognition of the compatible binding sites resulting in the formation of the complex. However, it is unclear how binding partners find each other in the context of the crowded, constantly fluctuating, and interaction-rich cellular environment. Here we examine the non-specific component of protein-protein interactions, which refers to those physicochemical properties of the binding partners that are independent of the exact details of their binding sites, but which can affect their localization or diffusional search for one another. We show that, for a large set of high-resolution experimental 3D structures of binary, transient protein complexes taken from the DOCKGROUND database, the binding partners display a surprising, statistically significant similarity in terms of their total hydration free energies normalized by a size-dependent variable. We hypothesize that colocalization of binding partners, even within individual cellular compartments such as the cytoplasm, may be influenced by their relative hydrophilicity, potentially in response to local hydrophilic gradients

Polarizable Water Model for the Coarse-Grained MARTINI Force Field

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models

Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the β1/β2, β3/β4, and β5/β6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK

An All-Atom Model of the Chromatin Fiber Containing Linker Histones Reveals a Versatile Structure Tuned by the Nucleosomal Repeat Length

In the nucleus of eukaryotic cells, histone proteins organize the linear genome into a functional and hierarchical architecture. In this paper, we use the crystal structures of the nucleosome core particle, B-DNA and the globular domain of H5 linker histone to build the first all-atom model of compact chromatin fibers. In this 3D jigsaw puzzle, DNA bending is achieved by solving an inverse kinematics problem. Our model is based on recent electron microscopy measurements of reconstituted fiber dimensions. Strikingly, we find that the chromatin fiber containing linker histones is a polymorphic structure. We show that different fiber conformations are obtained by tuning the linker histone orientation at the nucleosomes entry/exit according to the nucleosomal repeat length. We propose that the observed in vivo quantization of nucleosomal repeat length could reflect nature's ability to use the DNA molecule's helical geometry in order to give chromatin versatile topological and mechanical properties

Structural Heterogeneity and Quantitative FRET Efficiency Distributions of Polyprolines through a Hybrid Atomistic Simulation and Monte Carlo Approach

Förster Resonance Energy Transfer (FRET) experiments probe molecular distances via distance dependent energy transfer from an excited donor dye to an acceptor dye. Single molecule experiments not only probe average distances, but also distance distributions or even fluctuations, and thus provide a powerful tool to study biomolecular structure and dynamics. However, the measured energy transfer efficiency depends not only on the distance between the dyes, but also on their mutual orientation, which is typically inaccessible to experiments. Thus, assumptions on the orientation distributions and averages are usually made, limiting the accuracy of the distance distributions extracted from FRET experiments. Here, we demonstrate that by combining single molecule FRET experiments with the mutual dye orientation statistics obtained from Molecular Dynamics (MD) simulations, improved estimates of distances and distributions are obtained. From the simulated time-dependent mutual orientations, FRET efficiencies are calculated and the full statistics of individual photon absorption, energy transfer, and photon emission events is obtained from subsequent Monte Carlo (MC) simulations of the FRET kinetics. All recorded emission events are collected to bursts from which efficiency distributions are calculated in close resemblance to the actual FRET experiment, taking shot noise fully into account. Using polyproline chains with attached Alexa 488 and Alexa 594 dyes as a test system, we demonstrate the feasibility of this approach by direct comparison to experimental data. We identified cis-isomers and different static local environments as sources of the experimentally observed heterogeneity. Reconstructions of distance distributions from experimental data at different levels of theory demonstrate how the respective underlying assumptions and approximations affect the obtained accuracy. Our results show that dye fluctuations obtained from MD simulations, combined with MC single photon kinetics, provide a versatile tool to improve the accuracy of distance distributions that can be extracted from measured single molecule FRET efficiencies

The role of interfacial lipids in stabilizing membrane protein oligomers

Oligomerization of membrane proteins in response to lipid binding has a critical role in many cell-signalling pathways1 but is often difficult to define2 or predict3. Here we report the development of a mass spectrometry platform to determine simultaneously the presence of interfacial lipids and oligomeric stability and to uncover how lipids act as key regulators of membrane-protein association. Evaluation of oligomeric strength for a dataset of 125 α-helical oligomeric membrane proteins reveals an absence of interfacial lipids in the mass spectra of 12 membrane proteins with high oligomeric stability. For the bacterial homologue of the eukaryotic biogenic transporters (LeuT4, one of the proteins with the lowest oligomeric stability), we found a precise cohort of lipids within the dimer interface. Delipidation, mutation of lipid-binding sites or expression in cardiolipin-deficient Escherichia coli abrogated dimer formation. Molecular dynamics simulation revealed that cardiolipin acts as a bidentate ligand, bridging across subunits. Subsequently, we show that for the Vibrio splendidus sugar transporter SemiSWEET5, another protein with low oligomeric stability, cardiolipin shifts the equilibrium from monomer to functional dimer. We hypothesized that lipids are essential for dimerization of the Na+/H+ antiporter NhaA from E. coli, which has the lowest oligomeric strength, but not for the substantially more stable homologous Thermus thermophilus protein NapA. We found that lipid binding is obligatory for dimerization of NhaA, whereas NapA has adapted to form an interface that is stable without lipids. Overall, by correlating interfacial strength with the presence of interfacial lipids, we provide a rationale for understanding the role of lipids in both transient and stable interactions within a range of α-helical membrane proteins, including G-protein-coupled receptors