Inhibited complete folding of consecutive human telomeric G-quadruplexes.

Noncanonical DNA structures, termed G-quadruplexes, are present in human genomic DNA and are important elements in many DNA metabolic processes. Multiple sites in the human genome have G-rich DNA stretches able to support formation of several consecutive G-quadruplexes. One of those sites is the telomeric overhang region that has multiple repeats of TTAGGG and is tightly associated with both cancer and aging. We investigated the folding of consecutive G-quadruplexes in both potassium- and sodium-containing solutions using single-molecule FRET spectroscopy, circular dichroism, thermal melting and molecular dynamics simulations. Our observations show coexistence of partially and fully folded DNA, the latter consisting of consecutive G-quadruplexes. Following the folding process over hours in sodium-containing buffers revealed fast G-quadruplex folding but slow establishment of thermodynamic equilibrium. We find that full consecutive G-quadruplex formation is inhibited by the many DNA structures randomly nucleating on the DNA, some of which are off-path conformations that need to unfold to allow full folding. Our study allows describing consecutive G-quadruplex formation in both nonequilibrium and equilibrium conditions by a unified picture, where, due to the many possible DNA conformations, full folding with consecutive G-quadruplexes as beads on a string is not necessarily achieved.

Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins.

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 µs. Four combinations of protein and water force fields were tested: ff19SB/OPC, ff19SB/TIP4P-D, ff03CMAP/TIP4P-D, and a99SB-disp/TIP4P-disp, with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 µs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and 3 J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.

Lipid Modulation of a Class B GPCR: Elucidating the Modulatory Role of PI(4,5)P2 Lipids.

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) lipids have been shown to stabilize an active conformation of class A G-protein coupled receptors (GPCRs) through a conserved binding site, not present in class B GPCRs. For class B GPCRs, previous molecular dynamics (MD) simulation studies have shown PI(4,5)P2 interacting with the Glucagon receptor (GCGR), which constitutes an important target for diabetes and obesity therapeutics. In this work, we applied MD simulations supported by native mass spectrometry (nMS) to study lipid interactions with GCGR. We demonstrate how tail composition plays a role in modulating the binding of PI(4,5)P2 lipids to GCGR. Specifically, we find the PI(4,5)P2 lipids to have a higher affinity toward the inactive conformation of GCGR. Interestingly, we find that in contrast to class A GPCRs, PI(4,5)P2 appear to stabilize the inactive conformation of GCGR through a binding site conserved across class B GPCRs but absent in class A GPCRs. This suggests differences in the regulatory function of PI(4,5)P2 between class A and class B GPCRs.

Beneficent and Maleficent Effects of Cations on Bufadienolide Binding to Na+,K+-ATPase.

Kinetic properties and crystal structures of the Na+,K+-ATPase in complex with cardiotonic steroids (CTS) revealed significant differences between CTS subfamilies (Laursen et al.). Thus, we found beneficial effects of K+ on bufadienolide binding, which strongly contrasted with the well-known antagonism between K+ and cardenolides. In order to understand this peculiarity of bufalin interactions, we used docking and molecular dynamics simulations of the complexes involving Na+,K+-ATPase, bufadienolides (bufalin, cinobufagin), and ions (K+, Na+, Mg2+). The results revealed that bufadienolide binding is affected by (i) electrostatic attraction of the lactone ring by a cation and (ii) the ability of a cation to stabilize and "shape" the site constituted by transmembrane helices of the α-subunit (αM1-6). The latter effect was due to varying coordination patterns involving amino acid residues from helix bundles αM1-4 and αM5-10. Substituents on the steroid core of a bufadienolide add to and modify the cation effects. The above rationale is fully consistent with the ion effects on the kinetics of Na+,K+-ATPase/bufadienolide interactions.

General Protocol for Constructing Molecular Models of Nanodiscs.

Nanodisc technology is increasingly being applied for structural and biophysical studies of membrane proteins. In this work, we present a general protocol for constructing molecular models of nanodiscs for molecular dynamics simulations. The protocol is written in python and based on geometric equations, making it fast and easy to modify, enabling automation and customization of nanodiscs in silico. The novelty being the ability to construct any membrane scaffold protein (MSP) variant fast and easy given only an input sequence. We validated and tested the protocol by simulating seven different nanodiscs of various sizes and with different membrane scaffold proteins, both circularized and noncircularized. The structural and biophysical properties were analyzed and shown to be in good agreement with previously reported experimental data and simulation studies.

Conformational Changes in the 5-HT3A Receptor Extracellular Domain Measured by Voltage-Clamp Fluorometry.

The 5-hydroxytryptamine (5-HT) type 3 receptor is a member of the cysteine (Cys)-loop receptor super family of ligand-gated ion channels in the nervous system and is a clinical target in a range of diseases. The 5-HT3 receptor mediates fast serotonergic neurotransmission by undergoing a series of conformational changes initiated by ligand binding that lead to the rapid opening of an intrinsic cation-selective channel. However, despite the availability of high-resolution structures of a mouse 5-HT3 receptor, many important aspects of the mechanistic basis of 5-HT3 receptor function and modulation by drugs remain poorly understood. In particular, there is little direct evidence for the specific conformational changes predicted to occur during ligand-gated channel activation and desensitization. In the present study, we used voltage-clamp fluorometry (VCF) to measure conformational changes in regions surrounding the orthosteric binding site of the human 5-HT3A (h5-HT3A) receptor during binding of 5-HT and different classes of 5-HT3 receptor ligands. VCF utilizes parallel measurements of receptor currents with photon emission from fluorescent reporter groups covalently attached to specific positions in the receptor structure. Reporter groups that are highly sensitive to the local molecular environment can, in real time, report conformational changes as changes in fluorescence that can be correlated with changes in receptor currents reporting the functional states of the channel. Within the loop C, D, and E regions that surround the orthosteric binding site in the h5-HT3A receptor, we identify positions that are amenable to tagging with an environmentally sensitive reporter group that reports robust fluorescence changes upon 5-HT binding and receptor activation. We use these reporter positions to characterize the effect of ligand binding on the local structure of the orthosteric binding site by agonists, competitive antagonists, and allosterically acting channel activators. We observed that loop C appears to show distinct fluorescence changes for ligands of the same class, while loop D reports similar fluorescence changes for all ligands binding at the orthosteric site. In contrast, the loop E reporter position shows distinct changes for agonists, antagonists, and allosteric compounds, suggesting the conformational changes in this region are specific to ligand function. Interpretation of these results within the framework of current models of 5-HT3 and Cys-loop mechanisms are used to expand the understanding of how ligand binding in Cys-loop receptors relates to channel gating. SIGNIFICANCE STATEMENT: The 5-HT3 receptor is an important ligand-gated ion channel and drug target in the central and peripheral nervous system. Determining how ligand binding induced conformational changes in the receptor is central for understanding the structural mechanisms underlying 5-HT3 receptor function. Here, we employ voltage-gated fluorometry to characterize conformational changes in the extracellular domain of the human 5-HT3 receptor to identify intrareceptor motions during binding of a range of 5-HT3 receptor agonists and antagonists.

Cholesterol binding to a conserved site modulates the conformation, pharmacology, and transport kinetics of the human serotonin transporter.

The serotonin transporter (SERT) is important for reuptake of the neurotransmitter serotonin from the synaptic cleft and is also the target of most antidepressants. It has previously been shown that cholesterol in the membrane bilayer affects the conformation of SERT. Although recent crystal structures have identified several potential cholesterol-binding sites, it is unclear whether any of these potential cholesterol sites are occupied by cholesterol and functionally relevant. In the present study, we focus on the conserved cholesterol site 1 (CHOL1) located in a hydrophobic groove between TM1a, TM5, and TM7. By molecular dynamics simulations, we demonstrate a strong binding of cholesterol to CHOL1 in a membrane bilayer environment. In biochemical experiments, we find that cholesterol depletion induces a more inward-facing conformation favoring substrate analog binding. Consistent with this, we find that mutations in CHOL1 with a negative impact on cholesterol binding induce a more inward-facing conformation, and, vice versa, mutations with a positive impact on cholesterol binding induce a more outward-facing conformation. This shift in transporter conformation dictated by the ability to bind cholesterol in CHOL1 affects the apparent substrate affinity, maximum transport velocity, and turnover rates. Taken together, we show that occupation of CHOL1 by cholesterol is of major importance in the transporter conformational equilibrium, which in turn dictates ligand potency and serotonin transport activity. Based on our findings, we propose a mechanistic model that incorporates the role of cholesterol binding to CHOL1 in the function of SERT.

Revealing a Dual Role of Ganglioside Lipids in the Aggregation of Membrane-Associated Islet Amyloid Polypeptide.

Amyloid formation of the human islet amyloid polypeptide (hIAPP) correlates with a loss of insulin-producing beta cells in patients with type II diabetes mellitus. In this study, we investigated the binding of hIAPP to bilayers consisting of ganglioside lipids and dioleoylphosphatidylcholine (DOPC), which is a physiologically relevant lipid species for pancreatic beta cell-associated aggregation. The membrane interactions are studied computationally using a combination of coarse-grained, umbrella sampling, and atomistic molecular dynamics simulations. Herein, we demonstrate how the hIAPP peptides accumulate in the areas with a high content of ganglioside lipids. We have characterized two distinct binding modes of hIAPP on ganglioside-rich membranes, with both binding modes formed due to electrostatic interaction between the cationic peptides and the anionic ganglioside headgroup. We observed that binding in the ganglioside headgroup region induced conformational changes of the peptide towards an aggregation prone conformation, rich in ß-strands. In contrast, the binding of hIAPP near the ganglioside-enriched areas mobilizes the peptide, preventing it from conformational changes and potentially shields it from interactions with other peptides. This suggests a dual role of ganglioside lipids, affecting the aggregation of hIAPP by either accelerating or inhibiting amyloid formation depending on the membrane binding and the ganglioside concentration.

A direct interaction of cholesterol with the dopamine transporter prevents its out-to-inward transition.

Monoamine transporters (MATs) carry out neurotransmitter reuptake from the synaptic cleft, a key step in neurotransmission, which is targeted in the treatment of neurological disorders. Cholesterol (CHOL), a major component of the synaptic plasma membrane, has been shown to exhibit a modulatory effect on MATs. Recent crystal structures of the dopamine transporter (DAT) revealed the presence of two conserved CHOL-like molecules, suggesting a functional protein-CHOL direct interaction. Here, we present extensive atomistic molecular dynamics (MD) simulations of DAT in an outward-facing conformation. In the absence of bound CHOL, DAT undergoes structural changes reflecting early events of dopamine transport: transition to an inward-facing conformation. In contrast, in the presence of bound CHOL, these conformational changes are inhibited, seemingly by an immobilization of the intracellular interface of transmembrane helix 1a and 5 by CHOL. We also provide evidence, from coarse grain MD simulations that the CHOL sites observed in the DAT crystal structures are preserved in all human monoamine transporters (dopamine, serotonin and norepinephrine), suggesting that our findings might extend to the entire family.

Molecular Modeling Investigation of the Interaction between Humicola insolens Cutinase and SDS Surfactant Suggests a Mechanism for Enzyme Inactivation.

One of the largest commercial applications of enzymes and surfactants is as main components in modern detergents. The high concentration of surfactant compounds usually present in detergents can, however, negatively affect the enzymatic activity. To remedy this drawback, it is of great importance to characterize the interaction between the enzyme and the surfactant molecules at an atomistic resolution. The protein enzyme cutinase from the thermophilic and saprophytic fungus called Humicola insolens (HiC) is a promising candidate for use in detergents thanks to its hydrolase activity targeting mostly biopolyesters (e.g., cutin). HiC is, however, inhibited by low concentrations of sodium dodecyl sulfate (SDS), an ubiquitous surfactant. In this work, we investigate the interaction between HiC and SDS using molecular dynamics simulations. Simulations of HiC dissolved in different aqueous concentrations of SDS show the interaction between HiC and SDS monomers, as well as the formation and dynamics of SDS micelles on the surface of the enzyme. These results suggest a mechanism of cutinase inhibition by SDS, which involves the nucleation of aggregates of SDS molecules on hydrophobic patches on the cutinase surface. Notably, a primary binding site for monomeric SDS is identified near the active site of HiC constituting a possible nucleation point for micelles and leading to the blockage of the entrance to the enzymatic site. Detailed analysis of the simulations allow us to suggest a set of residues from the SDS binding site on HiC to probe as engineered mutations aimed at reducing SDS binding to HiC, thereby decreasing SDS inhibition of HiC.

Modeling and Mutational Analysis of the Binding Mode for the Multimodal Antidepressant Drug Vortioxetine to the Human 5-HT3A Receptor.

5-Hydroxytryptamine3 (5-HT3) receptors are ligand-gated ion channels that mediate neurotransmission by serotonin in the central nervous system. Pharmacological inhibition of 5-HT3 receptor activity has therapeutic potential in several psychiatric diseases, including depression and anxiety. The recently approved multimodal antidepressant vortioxetine has potent inhibitory activity at 5-HT3 receptors. Vortioxetine has an inhibitory mechanism that differs from classic 5-HT3 receptor competitive antagonists despite being believed to bind in the same binding site. Specifically, vortioxetine shows partial agonist activity followed by persistent and insurmountable inhibition. We have investigated the binding mode of vortioxetine at the human 5-HT3A receptor through computational and in vitro experiments to provide insight into the molecular mechanisms behind the unique pharmacological profile of the drug. We find that vortioxetine binds in a manner different from currently known 5-HT3A orthosteric ligands. Specifically, while the binding pattern of vortioxetine mimics some aspects of both the setron class of competitive antagonists and 5-hydroxytryptamine (5-HT) with regards to interactions with residues of the aromatic box motif in the orthosteric binding site, vortioxetine also forms interactions with residues not previously described to be important for the binding of either setrons or 5-HT such as Val202 on Loop F. Our results expand the framework for understanding how orthosteric ligands drive 5-HT3 receptor function, which is of importance for the potential future development of novel classes of 5-HT3 receptor antagonists.

On the effect of mutations in bovine or camel chymosin on the thermodynamics of binding κ-caseins.

Bovine and camel chymosins are aspartic proteases that are used in dairy food manufacturing. Both enzymes catalyze proteolysis of a milk protein, κ-casein, which helps to initiate milk coagulation. Surprisingly, camel chymosin shows a 70% higher clotting activity than bovine chymosin for bovine milk, while exhibiting only 20% of the unspecific proteolytic activity. By contrast, bovine chymosin is a poor coagulant for camel milk. Although both enzymes are marketed commercially, the disparity in their catalytic activity is not yet well understood at a molecular level, due in part to a lack of atomistic resolution data about the chymosin-κ-casein complexes. Here, we report computational alanine scanning calculations of all four chymosin-κ-casein complexes, allowing us to elucidate the influence that individual residues have on binding thermodynamics. Of the 12 sequence differences in the binding sites of bovine and camel chymosin, eight are shown to be particularly important for understanding differences in the binding thermodynamics (Asp112Glu, Lys221Val, Gln242Arg, Gln278Lys. Glu290Asp, His292Asn, Gln294Glu, and Lys295Leu. Residue in bovine chymosin written first). The relative binding free energies of single-point mutants of chymosin are calculated using the molecular mechanics three dimensional reference interaction site model (MM-3DRISM). Visualization of the solvent density functions calculated by 3DRISM reveals the difference in solvation of the binding sites of chymosin mutants.

Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology.

This is a perspective article entitled "Frontiers in computational biophysics: understanding conformational dynamics of complex lipid mixtures relevant to biology" which is following a CECAM meeting with the same name.

A direct view of the complex multi-pathway folding of telomeric G-quadruplexes.

G-quadruplexes (G4s) are DNA secondary structures that are capable of forming and function in vivo The propensity of G4s to exhibit extreme polymorphism and complex dynamics is likely to influence their cellular function, yet a clear microscopic picture of their folding process is lacking. Here we employed single-molecule FRET microscopy to obtain a direct view of the folding and underlying conformational dynamics of G4s formed by the human telomeric sequence in potassium containing solutions. Our experiments allowed detecting several folded states that are populated in the course of G4 folding and determining their folding energetics and timescales. Combining the single-molecule data with molecular dynamics simulations enabled obtaining a structural description of the experimentally observed folded states. Our work thus provides a comprehensive thermodynamic and kinetic description of the folding of G4s that proceeds through a complex multi-route pathway, involving several marginally stable conformational states.

Identification of Key Interactions in the Initial Self-Assembly of Amylin in a Membrane Environment.

Islet amyloid polypeptide, also known as amylin, forms aggregates that reduce the amount of insulin-producing cells in patients with type II diabetes mellitus. Much remains unknown about the process of aggregation and cytotoxicity, but it is known that certain cell membrane components can alter the rate of aggregation. Using atomistic molecular dynamics simulations combined with the highly mobile membrane mimetic model incorporating enhanced sampling of lipid diffusion, we investigate interaction of amylin peptides with the membrane components as well as the self-assembly of amylin. Consistent with experimental evidence, we find that an initial membrane-bound α-helical state folds into stable ß-sheet structures upon self-assembly. Our results suggest the following mechanism for the initial phase of amylin self-assembly. The peptides move around on the membrane with the positively charged N-terminus interacting with the negatively charged lipid headgroups. When the peptides start to interact, they partly unfold and break some of the contacts with the membrane. The initial interactions between the peptides are dominated by aromatic and hydrophobic interactions. Oligomers are formed showing both intra- and interpeptide ß-sheets, initially with interactions mainly in the C-terminal domain of the peptides. Decreasing the pH to 5.5 is known to inhibit amyloid formation. At low pH, His18 is protonated, adding a fourth positive charge at the peptide. With His18 protonated, no oligomerization is observed in the simulations. The additional charge gives a strong midpoint anchoring of the peptides to negatively charged membrane components, and the peptides experience additional interpeptide repulsion, thereby preventing interactions.

Reactive Center Loop Insertion in α-1-Antitrypsin Captured by Accelerated Molecular Dynamics Simulation.

Protease inhibition by metastable serine protease inhibitors (serpins) is mediated by one of the largest functional intradomain conformational changes known in biology. In this extensive structural rearrangement, protease-serpin complex formation triggers cleavage of the serpin reactive center loop (RCL), its subsequent insertion into central ß-sheet A, and covalent trapping of the target protease. In this study, we present the first detailed accelerated molecular dynamics simulation of the insertion of the fully cleaved RCL in α-1-antitrypsin (α1AT), the archetypal member of the family of human serpins. Our results reveal internal water pathways that allow the initial incorporation of side chains of RCL residues into the protein interior. We observed structural plasticity of the helix F (hF) element that blocks the RCL path in the native state, which is in excellent agreement with previous experimental reports. Furthermore, the simulation suggested a novel role of hF and the connected turn (thFs3A) as chaperones that support the insertion process by reducing the conformational space available to the RCL. Transient electrostatic interactions of RCL residues potentially fine-tune the serpin inhibitory activity. On the basis of our simulation, we generated the α1AT mutants K168E, E346K, and K168E/E346K and analyzed their inhibitory activity along with their intrinsic stability and heat-induced polymerization. Remarkably, the E346K mutation exhibited enhanced inhibitory activity along with an increased rate of premature structural collapse (polymerization), suggesting a significant role of E346 in the gatekeeping of the strain in the metastable native state.

Structural Basis for Simvastatin Competitive Antagonism of Complement Receptor 3.

The complement system is an important part of the innate immune response to infection but may also cause severe complications during inflammation. Small molecule antagonists to complement receptor 3 (CR3) have been widely sought, but a structural basis for their mode of action is not available. We report here on the structure of the human CR3 ligand-binding I domain in complex with simvastatin. Simvastatin targets the metal ion-dependent adhesion site of the open, ligand-binding conformation of the CR3 I domain by direct contact with the chelated Mg(2+) ion. Simvastatin antagonizes I domain binding to the complement fragments iC3b and C3d but not to intercellular adhesion molecule-1. By virtue of the I domain's wide distribution in binding kinetics to ligands, it was possible to identify ligand binding kinetics as discriminator for simvastatin antagonism. In static cellular experiments, 15-25 µm simvastatin reduced adhesion by K562 cells expressing recombinant CR3 and by primary human monocytes, with an endogenous expression of this receptor. Application of force to adhering monocytes potentiated the effects of simvastatin where only a 50-100 nm concentration of the drug reduced the adhesion by 20-40% compared with untreated cells. The ability of simvastatin to target CR3 in its ligand binding-activated conformation is a novel mechanism to explain the known anti-inflammatory effects of this compound, in particular because this CR3 conformation is found in pro-inflammatory environments. Our report points to new designs of CR3 antagonists and opens new perspectives and identifies druggable receptors from characterization of the ligand binding kinetics in the presence of antagonists.

Mutants and molecular dockings reveal that the primary L-thyroxine binding site in human serum albumin is not the one which can cause familial dysalbuminemic hyperthyroxinemia.

BACKGROUND: Natural mutations of R218 in human serum albumin (HSA) result in an increased affinity for L-thyroxine and lead to the autosomal dominant condition of familial dysalbuminemic hyperthyroxinemia. METHODS: Binding was studied by equilibrium dialysis and computer modeling. RESULTS: Ten of 32 other isoforms tested had modified high-affinity hormone binding. L-thyroxine has been reported to bind to four sites (Tr) in HSA; Tr1 and Tr4 are placed in the N-terminal and C-terminal part of the protein, respectively. Site-directed mutagenesis gave new information about all the sites. CONCLUSIONS: It is widely assumed that Tr1 is the primary hormone site, and that this site, on a modified form, is responsible for the above syndrome, but the binding experiments with the genetic variants and displacement studies with marker ligands indicated that the primary site is Tr4. This new assignment of the high-affinity site was strongly supported by results of MM-PBSA analyses and by molecular docking performed on relaxed protein structure. However, dockings also revealed that mutating R218 for a smaller amino acid increases the affinity of Tr1 to such an extent that it can become the high-affinity site. GENERAL SIGNIFICANCE: Placing the high-affinity binding site (Tr4) and the one which can result in familial dysalbuminemic hyperthyroxinemia (Tr1) in two very different parts of HSA is not trivial, because in this way persons with and without the syndrome can have different types of interactions, and thereby complications, when given albumin-bound drugs. The molecular information is also useful when designing drugs based on L-thyroxine analogues.

Simulations of membrane-bound diglycosylated human prion protein reveal potential protective mechanisms against misfolding.

Prion diseases are associated with the misfolding of the prion protein (PrP) from its normal cellular form (PrPC ) to its infectious scrapie form (PrPSc ). Post-translational modifications in PrP in vivo can play an important role in modulating the process of misfolding. To gain more insight into the effects of post-translational modifications in PrP structure and dynamics and to test the hypothesis that such modifications can interact with the protein, we have performed molecular dynamics simulations of diglycosylated human PrPC bound to a lipid bilayer via a glycophosphatidylinositol anchor. Multiple simulations were performed at three different pH ranges to explore pH effects on structure and dynamics. In contrast to simulations of protein-only PrPC , no large effects were observed upon lowering the pH of the system. The protein tilted toward the membrane surface in all of the simulations and the putative PrPSc oligomerization sites became inaccessible, thereby offering a possible protective mechanism against PrPSc -induced misfolding of PrPC .

Conformational Dynamics of the Human Islet Amyloid Polypeptide in a Membrane Environment: Toward the Aggregation Prone Form.

Human islet amyloid polypeptide (hIAPP) is a 37-residue peptide hormone, which upon misfolding changes from the physiologically active monomer into pathological amyloid fibril aggregates in the pancreas of type 2 diabetes mellitus patients. During this process, the insulin-producing pancreatic ß-cells are damaged; however, the underlying mechanism of this mode of cytotoxicity remains elusive. It is known that anionic lipids accelerate amyloid fibril formation, implicating the importance of the cellular membrane in the process, and that a pH close to the level in the ß-cell secretory granules (pH 5.5) inhibits amyloid fibril formation. Using all-atom molecular dynamics simulations, we have investigated the membrane-associated monomer state of α-helical hIAPP, analyzed specific interactions of hIAPP with a mixed anionic-zwitterionic lipid membrane and examined the influence of pH on the structure and dynamics of hIAPP and its interaction with the membrane. We find that hIAPP primarily interacts with the membrane by forming favorable interactions between anionic lipids and the positively charged residues in the N-terminal part of the peptide. Rationalizing experimental findings, the simulations show that the N-terminal part of the peptide interacts with the membrane in the lipid headgroup region. At neutral pH, the C-terminal part of the peptide, which contains the residues that initiate fibril formation, displays a highly dynamic, unfolded state, which interacts with the membrane significantly less than the N-terminal part. Such an unfolded form can be proposed to contribute to the acceleration of fibril formation. At low pH, protonation of His18 mediates a stronger interaction of the C-terminal part with the membrane, resulting in the immobilization of the C-terminal part on the membrane surface that might constitute a mechanism by which low pH inhibits fibril formation.