Impact of Cholesterol Concentration and Lipid Phase on Structure and Fluctuation of Amyloid Precursor Protein.

Elevated levels of cellular cholesterol have been identified as one factor contributing to the onset of Alzheimer's disease (AD). Specific interaction between cholesterol and the amyloid precursor protein (APP), investigated via NMR experiments and computational studies, has been proposed to play a critical role in the processing of APP by secretases and the biogenesis of amyloid-ß (Aß) protein. We present all-atom molecular dynamics simulations of the 40-residue congener of the C-terminal domain of APP, C9916-55 (C99), in cholesterol-enriched DMPC lipid bilayers. We investigated the effect of cholesterol concentration on the conformational ensemble of wild-type C99 and C99-cholesterol associations at the low pH of endosomal environments, at which residues E22 and D23 are neutral. C99 was also characterized in liquid ordered domains for Dutch (E22Q) and Iowa (D23N) Familial AD mutants at low pH and for the wild-type sequence using protonation states characteristic of neutral pH. Our results reproduce the equilibrium constant of past NMR characterizations of the C99-cholesterol interaction but are not consistent with the C99-cholesterol binding hypothesis. We find that the lifetimes of both DMPC and cholesterol complexed with C99 display a power-law distribution of residence lifetimes. Longer-lived C99-DMPC and C99-cholesterol complexes are primarily stabilized by salt bridges and hydrogen bonds of lysine amines to phosphate and hydroxyl groups. Nevertheless, specific interfaces for C99-cholesterol association which are not present for DMPC can be identified. Changes to C99-cholesterol interfaces are found to depend on C99 tilt angle and orientation of the juxtamembrane domain of C99 containing residues E22 and D23. These observations support a more nuanced view of the C99-cholesterol interaction than has previously been suggested. We propose that cholesterol modulates the conformation and activity of C99 and other small transmembrane proteins indirectly through induction of the liquid ordered phase and directly through hydrogen bonding. This suggests a critical role for membrane heterogeneity introduced by cholesterol in modulating the structural ensemble of C99 and the production of Aß.

M2 amphipathic helices facilitate pH-dependent conformational transition in influenza A virus.

The matrix-2 (M2) protein from influenza A virus is a tetrameric, integral transmembrane (TM) protein that plays a vital role in viral replication by proton flux into the virus. The His37 tetrad is a pH sensor in the center of the M2 TM helix that activates the channel in response to the low endosomal pH. M2 consists of different regions that are believed to be involved in membrane targeting, packaging, nucleocapsid binding, and proton transport. Although M2 has been the target of many experimental and theoretical studies that have led to significant insights into its structure and function under differing conditions, the main mechanism of proton transport, its conformational dynamics, and the role of the amphipathic helices (AHs) on proton conductance remain elusive. To this end, we have applied explicit solvent constant pH molecular dynamics using the multisite λ-dynamics approach (CpHMDMSλD) to investigate the buried ionizable residues comprehensively and to elucidate their effect on the conformational transition. Our model recapitulates the pH-dependent conformational transition of M2 from closed to open state when the AH domain is included in the M2 construct, revealing the role of the amphipathic helices on this transition and shedding light on the proton-transport mechanism. This work demonstrates the importance of including the amphipathic helices in future experimental and theoretical studies of ion channels. Finally, our work shows that explicit solvent CpHMDMSλD provides a realistic pH-dependent model for membrane proteins.

Exploring the impact of proteins on the line tension of a phase-separating ternary lipid mixture.

The separation of lipid mixtures into thermodynamically stable phase-separated domains is dependent on lipid composition, temperature, and system size. Using molecular dynamics simulations, the line tension between thermodynamically stable lipid domains formed from ternary mixtures of di-C16:0 PC:di-C18:2 PC:cholesterol at 40:40:20 mol. % ratio was investigated via two theoretical approaches. The line tension was found to be 3.1 ± 0.2 pN by capillary wave theory and 4.7 ± 3.7 pN by pressure tensor anisotropy approaches for coarse-grained models based on the Martini force field. Using an all-atom model of the lipid membrane based on the CHARMM36 force field, the line tension was found to be 3.6 ± 0.9 pN using capillary wave theory and 1.8 ± 2.2 pN using pressure anisotropy approaches. The discrepancy between estimates of the line tension based on capillary wave theory and pressure tensor anisotropy methods is discussed. Inclusion of protein in Martini membrane lipid mixtures was found to reduce the line tension by 25%-35% as calculated by the capillary wave theory approach. To further understand and predict the behavior of proteins in phase-separated membranes, we have formulated an analytical Flory-Huggins model and parameterized it against the simulation results. Taken together these results suggest a general role for proteins in reducing the thermodynamic cost associated with domain formation in lipid mixtures and quantifies the thermodynamic driving force promoting the association of proteins to domain interfaces.

Molecular Insights into Human Hereditary Apolipoprotein A-I Amyloidosis Caused by the Glu34Lys Mutation.

Hereditary apolipoprotein A-I (apoA-I) amyloidosis is a life-threatening incurable genetic disorder whose molecular underpinnings are unclear. In this disease, variant apoA-I, the major structural and functional protein of high-density lipoprotein, is released in a free form, undergoes an α-helix to intermolecular cross-ß-sheet conversion along with a proteolytic cleavage, and is deposited as amyloid fibrils in various organs, which can cause organ damage and death. Glu34Lys is the only known charge inversion mutation in apoA-I that causes human amyloidosis. To elucidate the structural underpinnings of the amyloidogenic behavior of Glu34Lys apoA-I, we generated its recombinant globular N-terminal domain (residues 1-184) and compared the conformation and dynamics of its lipid-free form with those of two other naturally occurring apoA-I variants, Phe71Tyr (amyloidogenic) and Leu159Arg (non-amyloidogenic). All variants showed reduced structural stability and altered aromatic residue packing. The greatest decrease in stability was observed in the non-amyloidogenic variant, suggesting that amyloid formation is driven by local structural perturbations at sensitive sites. Molecular dynamics simulations revealed local helical unfolding and suggested that transient opening of the Trp72 side chain induced mutation-dependent structural perturbations in a sensitive region, including the major amyloid hot spot residues Leu14-Leu22. We posit that a shift from the "closed" to the "open" orientation of the Trp72 side chain modulates structural protection of amyloid hot spots, suggesting a previously unknown early step in the protein misfolding pathway.

Critical size dependence of domain formation observed in coarse-grained simulations of bilayers composed of ternary lipid mixtures.

Model cellular membranes are known to form micro- and macroscale lipid domains dependent on molecular composition. The formation of macroscopic lipid domains by lipid mixtures has been the subject of many simulation investigations. We present a critical study of system size impact on lipid domain phase separation into liquid-ordered and liquid-disordered macroscale domains in ternary lipid mixtures. In the popular di-C16:0 PC:di-C18:2 PC:cholesterol at 35:35:30 ratio mixture, we find systems with a minimum of 1480 lipids to be necessary for the formation of macroscopic phase separated domains and systems of 10 000 lipids to achieve structurally converged conformations similar to the thermodynamic limit. To understand these results and predict the behavior of any mixture forming two phases, we develop and investigate an analytical Flory-Huggins model which is recursively validated using simulation and experimental data. We find that micro- and macroscale domains can coexist in ternary mixtures. Additionally, we analyze the distributions of specific lipid-lipid interactions in each phase, characterizing domain structures proposed based on past experimental studies. These findings offer guidance in selecting appropriate system sizes for the study of phase separations and provide new insights into the nature of domain structure for a popular ternary lipid mixture.

Exploring the structure and stability of cholesterol dimer formation in multicomponent lipid bilayers.

For 40 years, the existence and possible functional importance of cholesterol dimer formation has been discussed. Due to challenges associated with structural studies of membrane lipids, there has as yet been no direct experimental verification of the existence and relevance of the cholesterol dimer. Building on recent advances in lipid force fields for molecular simulation, in this work the structure and stability of the cholesterol dimer is characterized in POPC bilayers in absence and presence of sphingomyelin. The cholesterol dimer structural ensemble is found to consist of sub-states that reflect, but also differ from, previously proposed dimer structures. While face-to-face dimer structures predominate, no evidence is found for the existence of tail-to-tail dimers in POPC lipid bilayers. Near stoichiometric complex formation of cholesterol with sphingomyelin is found to effect cholesterol dimer structure without impacting population. Comparison with NMR-derived order parameters provide validation for the simulation model employed and conclusions drawn related to the structure and stability of cholesterol dimers in multicomponent lipid bilayers. © 2016 Wiley Periodicals, Inc.

Specific Binding of Cholesterol to C99 Domain of Amyloid Precursor Protein Depends Critically on Charge State of Protein.

Recent NMR chemical shift measurements of the 99 residue C-terminal fragment of amyloid precursor protein (APP-C99) in the presence of cholesterol provide evidence of binary complex formation between C99 and cholesterol in membrane mimetic environments. It has also been observed that the production of Aß protein is enhanced under conditions of high cholesterol concentration. In this study, we investigated the impact of the charge state of C99 on the structure and stability of the C99-cholesterol complex. We observed that the binding of C99 to cholesterol depends critically on the charge state of Glu 693 (E22) and Asp 694 (D23). Evaluation of the pKa values of the Asp and Glu side chains suggests that these residues may be predominantly neutral in existing experimental observations of a stable C99-cholesterol complex at lower pH (characteristic of the endosomal environment), while binding is destabilized near neutral pH (characteristic of the cytoplasm). These observations suggest that specific binding of cholesterol to C99 is a sensitive function of the pH encountered in vivo, with key E22 and D23 residues serving as a "pH switch" controlling C99-cholesterol binding.

Probing Site-Specific Structural Information of Peptides at Model Membrane Interface In Situ.

Isotope labeling is a powerful technique to probe detailed structures of biological molecules with a variety of analytical methods such as NMR and vibrational spectroscopies. It is important to obtain molecular structural information on biological molecules at interfaces such as cell membranes, but it is challenging to use the isotope labeling method to study interfacial biomolecules. Here, by individually (13)C═(16)O labeling ten residues of a peptide, Ovispirin-1, we have demonstrated for the first time that a site-specific environment of membrane associated peptide can be probed by the submonolayer surface sensitive sum frequency generation (SFG) vibrational spectroscopy in situ. With the peptide associated with a single lipid bilayer, the sinusoidal trend of the SFG line width and peak-center frequency suggests that the peptide is located at the interface beneath the lipid headgroup region. The constructive interferences between the isotope labeled peaks and the main peptide amide I peak contributed by the unlabeled components were used to determine the membrane orientation of the peptide. From the SFG spectral peak-center frequency, line width, and polarization dependence of the isotope labeled units, we deduced structural information on individual units of the peptide associated with a model cell membrane. We also performed molecular dynamics (MD) simulations to understand peptide-membrane interactions. The physical pictures described by simulation agree well with the SFG experimental result. This research demonstrates the feasibility and power of using isotope labeling SFG to probe molecular structures of interfacial biological molecules in situ in real time.

Membrane environment modulates the pKa values of transmembrane helices.

In this work, we apply the recently developed constant pH molecular dynamics technique to study protonation equilibria of titratable side chains in the context of simple transmembrane (TM) helices and explore the effect of pH on their configurations in membrane bilayers. We observe that, despite a significant shift toward neutral states, considerable population of different side chains stay in the charged state that give rise to pKa values around 9.6 for Asp and Glu and 4.5 to 6 for His and Lys side chains, respectively. These charged states are highly stabilized by favorable interactions between head groups, water molecules, and the charged side chains that are facilitated by substantial changes in the configuration of the peptides. The pH dependent configurations and the measured pKa values are in good agreement with relatively recent solid state NMR measurements. Our results presented here demonstrate that all-atom constant pH molecular dynamics can be applied to membrane proteins and peptides to obtain reliable pKa values and pH dependent behavior for these systems.

Hamiltonian Mapping Revisited: Calibrating Minimalist Models to Capture Molecular Recognition by Intrinsically Disordered Proteins.

Molecular recognition by intrinsically disordered proteins (IDPs) plays a central role in many critical cellular processes. Toward achieving detailed mechanistic understanding of IDP-target interactions, here we employ the "Hamiltonian mapping" methodology, which is rooted in the weighted histogram analysis method (WHAM), for the fast and efficient calibration of structure-based models in studies of IDPs. By performing reference simulations on a given Hamiltonian, we illustrate for two model IDPs how this method can extrapolate thermodynamic behavior under a range of modified Hamiltonians, in this case representing changes in the binding affinity (Kd) of the system. Given sufficient conformational sampling in a single trajectory, Hamiltonian mapping accurately reproduces Kd values from direct simulation. This method may be generally applied to systems beyond IDPs in force field optimization and in describing changes in thermodynamic behavior as a function of external conditions for connection with experiment.

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.

Dynamic Heterogeneous Dielectric Generalized Born (DHDGB): An implicit membrane model with a dynamically varying bilayer thickness.

An extension to the heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is presented to allow dynamic membrane deformations in response to membrane-inserted biomolecules during molecular dynamic simulations. The flexible membrane is implemented through additional degrees of freedom that represent the membrane deformation at the contact points of a membrane-inserted solute with the membrane. The extra degrees of freedom determine the dielectric and non-polar solvation free energy profiles that are used to obtain the solvation free energy in the presence of the membrane and are used to calculate membrane deformation free energies according to an elastic membrane model. With the dynamic HDGB (DHDGB) model the membrane is able to deform in response to the insertion of charged molecules thereby avoiding the overestimation of insertion free energies with static membrane models. The DHDGB model also allows the membrane to respond to the insertion of membrane-spanning solutes with hydrophobic mismatch. The model is tested with the membrane insertion of amino acid side chain analogs, arginine-containing helices, the WALP23 peptide, and the gramicidin A channel.

Effect of flanking residues on the conformational sampling of the internal fusion peptide from Ebola virus.

Fusion peptides mediate viral and host-cell membrane fusion during viral entry. The monomeric form of the internal fusion peptide from Ebola virus was studied in membrane bilayer and water environments with computer simulations using replica exchange sampling and an implicit solvent description of the environment. Wild-type Ebola fusion peptide (EFP), the W8A mutant form, and an extended construct with flanking residues were examined. It was found that the monomeric form of wild-type EFP adopts coil-helix-coil structure with a short helix from residues 8 to 11 mostly sampling orientations parallel to the membrane surface. W8A mutation disrupts the helicity in the N-terminal region of the peptide and leads to a preference for slightly oblique orientation relative to the membrane surface. The addition of flanking residues also alters the fusion peptide conformation with either a helix-break-helix structure or extended N and C-termini and reduced membrane insertion. In water, the fusion peptide is found to adopt structures with low helicity.

Conformational sampling of influenza fusion peptide in membrane bilayers as a function of termini and protonation states.

Influenza fusion peptide is critical for mediating the fusion of viral and host cell membranes during viral entry. The interaction of monomeric influenza fusion peptide with membranes is studied with replica exchange molecular dynamics simulations using a new implicit membrane model to effectively reach microsecond to millisecond time scales. The conformational sampling of the fusion peptide was studied as a function of different N- and C-termini, including an experimental construct with an additional C-terminal tag, as well as a function of protonation of acidic residues. It is found that the influenza fusion peptide mostly adopts helical structures with a pronounced kink at residues 11-13 with both N-terminal and C-terminal helices oriented mostly parallel to the membrane surface. A charged C-terminus and the presence of a charge C-terminal tag significantly alters the conformational sampling of the fusion peptide and results in more diverse conformational ensembles that include obliquely inserted N-terminal peptide structures. Protonation of acidic residues also affects the conformational sampling, however, based on pK(a) shift estimates the overall effect of pH = 5 on the conformational sampling of the influenza fusion peptide appears to be only minor.