Search

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 (CpHMD(MSλ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 CpHMD(MSλD) provides a realistic pH-dependent model for membrane proteins.

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.

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.

Hereditary apolipoprotein A-I (apoA-I) amyloidosis (AApoAI) is a life-threatening incurable genetic disorder whose molecular underpinnings and the full spectrum of afflicted organs are unclear. We report a new form of AApoAI with amyloid deposition in multiple organs, including an unprecedented retinal amyloidosis. Genetic and proteomic analyses identified Glu34Lys apoA-I as the fibrillar protein causing the clinical manifestations. A life-saving combined hepatorenal transplantation was performed for one Glu34Lys carrier. To elucidate structural underpinnings for amyloidogenic properties of Glu34Lys, we generated its recombinant globular domain 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 stability and altered aromatic residue packing. Molecular dynamics simulations revealed local helical unfolding and suggested that transient opening of Trp72 induced mutation-dependent structural perturbations in a sensitive region, including the major amyloid hotspot residues 14-22. We posit that a shift from the “closed” to an “open” orientation of Trp72 modulates structural protection of amyloid hotspots, suggesting a previously unknown early step in protein misfolding.

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.

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.

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.

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.

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.

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.

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.

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.

Environmental conditions such as pH are considered to play profound roles in determining proteins' structures and functions. In recent years, several experimental and computational techniques have been developed to improve our understanding of the effects of pH on proteins. However, the implicit representation of solvents in the current constant pH molecular dynamics (CpHMD) methods can lead to undesirable side effects in calculating pKa values of titratable side chains. These effects can be even more pronounced for membrane proteins where the rigid nature of the current implicit membrane models causes unrealistic pKa shifts in favor of neutral states. Inspired by recently published pKa values for Lys residues in model transmembrane peptides (Gleason NJ et al. (2013) PNAS 110:1692-1695), here we explore a newly developed CpHMD technique based on an all-atom representation of solvent/lipid molecules and proteins for simple transmembrane peptide models with a titratable side chain. The pKa values of the side chains along with the pH dependent configuration of these peptides in the membrane bilayer are investigated and compared with the available experimental values. Motivated by the excellent agreement between the calculated pKa values and the experiment for model peptides, the purely explicit CpHMD is also applied to study the pKa values of titratable side chains of influenza A proton channel (M2) and its pH dependent structural change in the membrane bilayer.

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, the W8A mutant form, and an extended construct with flanking residues were examined. It was found that the monomeric form of wild type Ebola fusion peptide adopts coil-helix-coil structure with a short helix from residue 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.

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.