MD simulations of spontaneous membrane protein/detergent micelle formation

The in vitro study of membrane proteins for the purpose of physicochemical analysis or structure determination often relies upon successful reconstitution into detergent micelles. Moreover, a number of biological processes such as membrane protein folding and transport rely on lipid interactions which may resemble the micellar environment. Little is known about the structures of these micelles or the processes which lead to their formation. We therefore present two 50 ns all-atom molecular dynamics simulations of spontaneous dodecylphosphocholine micelle formation around representatives of the two major families of membrane proteins, a small beta-barrel protein, OmpA, and a model alpha-helical protein, glycophorin A. Despite differences in protein architecture, we highlight common mechanistic pathways in micelle formation, which are consistent with experimental studies. We characterize the exponential kinetics of detergent-protein adsorption and suggest a simple model which may explain the aggregation process. We also compare the results with 25 and 50 ns simulations of preformed micelles containing the same proteins. We confirm that the end structures of the self-assembled micelles are similar to those from their preformed counterparts, with each micelle presenting a bilayerlike environment to the enclosed protein.     

Ultraviolet resonance Raman spectroscopy of folded and unfolded states of an integral membrane protein

The vibrational structure of native anchoring tryptophan (Trp) and tyrosine residues in an integral membrane protein, bacterial outer membrane protein A (OmpA), have been investigated using UV resonance Raman (UVRR) spectroscopy for the first time. Spectra of native OmpA, a single-Trp mutant, and a Trp-less mutant were recorded in folded and unfolded states, and reveal significant changes in tryptophan structure and local environment. Salient alterations upon folding include loss of hydrogen-bonding character of indole N1H, evidenced by a shift in W17 frequency from 874 and 878 cm(-1), and growth in hydrophobicity of the local tryptophan environment, supported by increase in the ratio I1361/I1340. In addition to these site-specific changes in a single tryptophan residue, modification of the vibrational structure of the remaining native tryptophan and tyrosine amino acids is also evident. Finally, the UVRR data presented here indicate that the structures of OmpA folded in vesicle and folded in detergent may differ, and provide important foundations for ongoing studies of membrane protein folding.     

Multidimensional umbrella sampling and replica‐exchange molecular dynamics simulations for structure prediction of transmembrane helix dimers

Structural information of a transmembrane (TM) helix dimer is useful in understanding molecular mechanisms of important biological phenomena such as signal transduction across the cell membrane. Here, we describe an umbrella sampling (US) scheme for predicting the structure of a TM helix dimer in implicit membrane using the interhelical crossing angle and the TM–TM relative rotation angles as the reaction coordinates. This scheme conducts an efficient conformational search on TM–TM contact interfaces, and its robustness is tested by predicting the structures of glycophorin A (GpA) and receptor tyrosine kinase EphA1 (EphA1) TM dimers. The nuclear magnetic resonance (NMR) structures of both proteins correspond to the global free‐energy minimum states in their free‐energy landscapes. In addition, using the landscape of GpA as a reference, we also examine the protocols of temperature replica‐exchange molecular dynamics (REMD) simulations for structure prediction of TM helix dimers in implicit membrane. A wide temperature range in REMD simulations, for example, 250–1000 K, is required to efficiently obtain a free‐energy landscape consistent with the US simulations. The interhelical crossing angle and the TM–TM relative rotation angles can be used as reaction coordinates in multidimensional US and be good measures for conformational sampling of REMD simulations. © 2013 Wiley Periodicals, Inc.     

Dramatic destabilization of transmembrane helix interactions by features of natural membrane environments

Membrane proteins have evolved to fold and function in a lipid bilayer, so it is generally assumed that their stability should be optimized in a natural membrane environment. Yet optimal stability is not always in accord with optimization of function, so evolutionary pressure, occurring in a complex membrane environment, may favor marginal stability. Here, we find that the transmembrane helix dimer, glycophorin A (GpATM), is actually much less stable in the heterogeneous environment of a natural membrane than it is in model membranes and even common detergents. The primary destabilizing factors are electrostatic interactions between charged lipids and charged GpATM side chains, and nonspecific competition from other membrane proteins. These effects overwhelm stabilizing contributions from lateral packing pressure and excluded volume. Our work illustrates how evolution can employ membrane composition to modulate protein stability.     

Efficient molecular mechanics simulations of the folding,orientation, and assembly of peptides in lipid bilayers using an implicit atomic solvation model

Membrane proteins comprise a significant fraction of the proteomes of sequenced organisms and are the targets of approximately half of marketed drugs. However, in spite of their prevalence and biomedical importance, relatively few experimental structures are available due to technical challenges. Computational simulations can potentially address this deficit by providing structural models of membrane proteins. Solvation within the spatially heterogeneous membrane/solvent environment provides a major component of the energetics driving protein folding and association within the membrane. We have developed an implicit solvation model for membranes that is both computationally efficient and accurate enough to enable molecular mechanics predictions for the folding and association of peptides within the membrane. We derived the new atomic solvation model parameters using an unbiased fitting procedure to experimental data and have applied it to diverse problems in order to test its accuracy and to gain insight into membrane protein folding. First, we predicted the positions and orientations of peptides and complexes within the lipid bilayer and compared the simulation results with solid-state NMR structures. Additionally, we performed folding simulations for a series of host–guest peptides with varying propensities to form alpha helices in a hydrophobic environment and compared the structures with experimental measurements. We were also able to successfully predict the structures of amphipathic peptides as well as the structures for dimeric complexes of short hexapeptides that have experimentally characterized propensities to form beta sheets within the membrane. Finally, we compared calculated relative transfer energies with data from experiments measuring the effects of mutations on the free energies of translocon-mediated insertion of proteins into lipid bilayers and of combined folding and membrane insertion of a beta barrel protein.     

非天然相互作用对均聚多肽链螺旋形成动力学过程的影响的计算机模拟研究

以螺旋结构的形成过程为研究对象,基于粗粒化的格子模型和动态蒙特卡罗模拟方法,初步探讨了非天然氢键相互作用对均聚多肽链螺旋折叠动力学过程的影响.研究发现,非天然氢键的引入虽然延缓了其热力学转变的发生,但也从整体上降低了折叠动力学过程的能垒,在某一特定温度之下,反而可以提高折叠速率.对其折叠路径分布的分析表明,非天然氢键可以减少慢速折叠路径的发生,而后者是导致折叠时间增加的主要因素.另一方面,比较特定温度下多肽链链构象及螺旋片断随时间的演化进程,发现非天然氢键在一定程度上影响了天然氢键的形成以及天然态构象的稳定存在,同时也加快了其部分解折叠过程.这说明,非天然相互作用的存在有利于蛋白质构象的快速动态调整,从而行使其相应的生物功能.     

CHARMM-GUI Nanodisc Builder for modeling and simulation of various nanodisc systems

Nanodiscs are discoidal protein–lipid complexes that have wide applications in membrane protein studies. Modeling and simulation of nanodiscs are challenging due to the absence of structures of many membrane scaffold proteins (MSPs) that wrap around the membrane bilayer. We have developed CHARMM-GUI Nanodisc Builder ( http://www.charmm-gui.org/input/nanodisc ) to facilitate the setup of nanodisc simulation systems by modeling the MSPs with defined size and known structural features. A total of 11 different nanodiscs with a diameter from 80 to 180 Å are made available in both the all-atom CHARMM and two coarse-grained (PACE and Martini) force fields. The usage of the Nanodisc Builder is demonstrated with various simulation systems. The structures and dynamics of proteins and lipids in these systems were analyzed, showing similar behaviors to those from previous all-atom and coarse-grained nanodisc simulations. We expect the Nanodisc Builder to be a convenient and reliable tool for modeling and simulation of nanodisc systems. © 2019 Wiley Periodicals, Inc.     

Self-assembly of Asymmetric Dimer Particles in Supported Copolymer Bilayer

Using self-consistent field and density functional theories, we investigate the self-assembly behavior of asymmetric dimer particles in a supported AB block copolymer bilayer. Asym-metric dimer particles are amphiphilic molecules composed by two different spheres. One prefers to A block of copolymers and the other likes B block when they are introduced into the copolymer bilayer. The two layer structure of the dimer particles is formed within the bilayer. Due to the presence of the substrate surface, the symmetry of the two leaflets of the bilayer is broken, which may lead to two different layer structures of dimer particles within each leaflet of the bilayer. With the increasing concentration of the asymmetric dimer particles, in-plane structure of the dimer particles undergoes sparse square, hexagonal, dense square, and cylindrical structures. In a further condensed packing, a bending cylindrical structurecomes into being. Here we verify that the entropic effect of copolymers, the enthalpy of the system and the steric repulsion of the dimer particles are three important factors determing the self-assembly of dimer particles within the supported copolymer bilayer.     

Insights into the recognition and association of transmembrane alpha-helices. The free energy of alpha-helix dimerization in glycophorin A

The free energy of alpha-helix dimerization of the transmembrane (TM) region of glycophorin A was estimated from a 125-ns molecular dynamics (MD) simulation in a membrane mimetic. The free energy profile was obtained by allowing the TM helical segments to diffuse reversibly along the reaction pathway. Partition of the potential of mean force into free energy components illuminates the critical steps of alpha-helix recognition and association. At large separations, the TM segments are pushed together by the solvent, allowing initial, but not necessarily native, interhelical interactions to occur. This early recognition stage precedes the formation of native contacts, which is accompanied by a tilt of the helices, characteristic of the dimeric structure. This step is primarily driven by the van der Waals helix-helix interactions. Free energy perturbation calculations of the L75A and I76A point mutations reveal a disruption in helix-helix association due to a loss of favorable dispersion interactions. Additional MD simulations of the native TM dimer and of a single alpha-helix confirm that, prior to association, individual alpha-helices are independently stable, in agreement with the "two-stage" model of integral membrane protein folding.     

Alamethicin channels in a membrane: molecular dynamics simulations

Alamethicin (Alm) is a 20 residue peptide which forms a kinked alpha-helix in membrane and membrane-mimetic environments. Ion channels formed by intramembraneous aggregates of Alm are thought to be formed by bundles of approximately parallel Alm helices surrounding a central bilayer pore. Different channel conductance levels correspond to different numbers of helices per bundle, ranging from N = 5 to N > 8. Calculation of the predicted pKA values of the ring of Glu18 sidechains at the C-terminal mouth of the pore suggests that at neutral pH most or all of these sidechains will remain protonated. Nanosecond molecular dynamics (MD) simulations of N = 5, 6, 7 and 8 bundles of Alm helices in a POPC bilayer have been run, corresponding to a total simulation time of 4 ns. These simulations explore the stability and conformational dynamics of these helix bundle channels when embedded in a full phospholipid bilayer in an aqueous environment. The structural and dynamic properties of water in these model channels are examined. As in earlier in vacuo simulations (J. Breed, R. Sankararamakrishnan, I. D. Kerr and M. S. P. Sansom, Biophys. J., 1996, 70, 1643) the dipole moments of water molecules within the pores are aligned antiparallel to the helix dipoles. This helps to contribute to the stability of the helix bundles.     

Implicit two-phase solvation model as a tool to assess conformation and energetics of proteins in membrane-mimetic media

 A heterogeneous implicit membrane-mimetic model is applied to simulations of membrane proteins. The model employs atomic solvation parameters for gas–water and gas–cyclohexane transfer. It is used to analyze structure, energetics, and orientation with respect to the bilayer of two polypeptides with different modes of membrane binding – hydrophobic segment of human glycophorin A (GpA) and cytotoxin II from Naja naja oxiana snake venom (CTX). The native state of GpA represents a transmembrane (TM) α helix, while CTX is a water-soluble protein, which is able to interact with the cell membrane. The conformational space of the polypeptides was explored in Monte Carlo simulations. The results show that the most stable conformers of GpA represent a TM α helix. They are additionally stabilized by an applied TM voltage. The results also show that CTX inserts with its three loops, does not cross the hydrophobic layer, and stays partially immersed in the membrane. This agrees well with the experimental data, thus confirming the validity of the solvation model. Received: 13 June 2000 / Accepted: 15 September 2000 / Published online: 19 January 2001     

OmpT在Tween-20胶束中折叠的单分子研究

The cell envelope of gram-negative bacteria consists of the outer membrane (OM), inner membrane (IM), and periplasm. The β-barrel outer membrane proteins (OMPs) embedded in the OM perform diverse and significant functions such as signaling, transporting, and proteolysis. The OMPs of gram-negative bacteria share similar folding pathways with that of mitochondria and chloroplasts. Therefore, the study of the OMP folding mechanism not only provides insights into antimicrobial drug design but also helps elucidate mitochondrial and chloroplast biogenesis. Most knowledge about OMP folding was obtained from ensemble experiments where OMPs were usually at micromolar concentrations and prone to aggregate, which is different from the physiological environment in the cells. Unlike ensemble techniques, single-molecule detection (SMD) can measure OMPs from nano- to picomolar concentrations and prevent aggregation. In this work, we investigated the folding of OmpT, one of the OMPs, in Tween-20 and n-dodecyl β-d-maltopyranoside (DDM) micelles by SMD. We prepared monodisperse OmpT and observed both unfolded and folded OmpT in Tween-20 and DDM micelles under different urea concentrations by single-molecule fluorescence resonance energy transfer (FRET). The folded OmpT in Tween-20 is structurally similar to the native OmpT folded in DDM but exhibits weaker resistance to urea. In contrast, OmpA barely folds and OmpC hardly folds in Tween-20 micelles. We confirmed that folded OmpT forms complexes with detergent micelles and estimated the number of bound Tween-20 and DDM molecules per OmpT by fluorescence correlation spectroscopy. We compared the enzymatic activity of OmpT folded in two detergents with a fluorescent peptide as substrate, and found that the folded form of OmpT in Tween-20 possesses weaker enzymatic activity than that in DDM. We also investigated the folding properties of OmpT, OmpA, and OmpC in the presence of the β-barrel assembly machine (BAM) complex. OmpT folds efficiently in liposome even without the BAM complex; OmpA only folds with the help of the BAM complex; and OmpC does not fold with or without the BAM complex. Based on the comparison of the folding of OmpT, OmpA, and OmpC in detergent micelles and in the presence of the BAM complex, we propose that OmpT has stronger folding tendency than OmpA and OmpC, which supports the idea that the exact role of the BAM complex is dependent on the distinct folding properties of individual OMPs. Since Tween-20 is a widely used reagent to block nonspecific adsorption in SMD experiments, our results also remind people to exercise caution to prevent possible wrong interpretations caused by the interaction between proteins and Tween-20.     

Cytoskeleton mediated effective elastic properties of model red blood cell membranes

The plasma membrane of human red blood cells consists of a lipid bilayer attached to a regular network of underlying cytoskeletal polymers. We model this system at a dynamic coarse-grained level, treating the bilayer as an elastic sheet and the cytoskeletal network as a series of phantom entropic springs. In contrast to prior simulation efforts, we explicitly account for dynamics of the cytoskeletal network, both via motion of the protein anchors that attach the cytoskeleton to the bilayer and through breaking and reconnection of individual cytoskeletal filaments. Simulation results are explained in the context of a simple mean field percolation model and comparison is made to experimental measurements of red blood cell fluctuation amplitudes.     

Dynamic control of protein folding pathway with a polymer of tunable hydrophobicity

While the knowledge of protein folding in a dilute solution is now well-advanced, little is known of the influence of surrounding conditions on the folding kinetics, in particular when the protein is in a dynamically responsive environment. Here we report a new procedure to control the pathways of protein folding by using a thermally responsive polymer that varies its hydrophobicity concomitant with the protein structural changes. The advantages of folding in a dynamic environment have been demonstrated first by Langevin dynamics simulations on the basis of coarse-grained models for both the protein and polymer and then by experiments for lysozyme refolding in the presence of poly(N-isopropylacrylamide-co-N-tert-butylacrylamide), a thermal responsive polymer that varies its hydrophobicity in response to temperature. The simulation suggests that decreasing the polymer hydrophobicity during the folding process may result in an optimized free-energy landscape that enhances both the folding yield and kinetics. The experiments affirm that an optimal folding condition can be identified when structural transitions of the protein collaborate with the polymer hydrophobicity tuned by variation of temperature.     

Formation of Vesicles from Amphiphilic Random Copolymers in Solution: A Dissipative Particle Dynamics Simulation Study

Five coarse-grained models were built for amphiphilic random copolymers. The self-assembly of amphiphilic random copolymers in selective solvent was investigated via dissipative particle dynamics simulations. The simulation results showed that the content of hydrophilic particles and the repulsive parameter between solvent and copolymer particles were two key factors of the vesicle formation. We report herein on how to control the self-assembled morphology evolution. The two mechanisms of vesicle formation from amphiphilic random copolymers are found through investigating the dynamic processes of vesicle formation, which is in accordance with the experiment and simulation results of amphiphilic block copolymer reported in the literature.      

Lipid/protein interactions and the membrane/water interfacial region

Despite the importance of lipid/protein interactions in the folding, assembly, stability, and function of membrane proteins, information at an atomic level on how such proteins interact with the lipids that surround them remains sparse. The dynamics and flexible nature of the protein/bilayer interaction make it difficult to study, for example, by crystallographic means. However, based on recent progress in molecular simulations of membranes it is possible to address this problem computationally. This communication reports one of the first attempts to use multiple ns molecular simulations to establish a qualitative picture of the intermolecular interactions between the lipids of a bilayer and two topologically different membrane proteins for which a high resolution (2 A or better) X-ray structure is available.     

The Self-Association of the KRAS4b Protein is Altered by Lipid-Bilayer Composition and Electrostatics

KRAS is a peripheral membrane protein that regulates multiple signaling pathways, and is mutated in ≈30 % of cancers. Transient self-association of KRAS is essential for activation of the downstream effector RAF and oncogenicity. The presence of anionic phosphatidylserine (PS) lipids in the membrane was shown to promote KRAS self-assembly, however, the structural mechanisms remain elusive. Here, we employed nanodisc bilayers of defined lipid compositions, and probed the impact of PS concentration on KRAS self-association. Paramagnetic NMR experiments demonstrated the existence of two transient dimer conformations involving alternate electrostatic contacts between R135 and either D153 or E168 on the “α4/5-α4/5” interface, and revealed that lipid composition and salt modulate their dynamic equilibrium. These dimer interfaces were validated by charge-reversal mutants. This plasticity demonstrates how the dynamic KRAS dimerization interface responds to the environment, and likely extends to the assembly of other signaling complexes on the membrane.     

Coarse-grained Dynamics Simulation in Polymer Systems: from Structures to Material Properties

Polymers are widely used in our daily life and industry because of their intrinsic characteristics, such as multi-functionality, low cost, light mass, ease of processability, and excellent chemical stability. Polymers have multiscale space-time properties, which are mainly reflected in the fact that the properties of polymer systems depend not only on chemical structure and molecular properties, but also to a large extent on the aggregation state of molecules, that is, phase structure and condensed state structure. Thanks to the continuous development of simulation methods and the rapid improvement of scientific computation, computer simulation has played an increasingly important role in investigating the structure and properties of polymer systems. Among them, coarse-grained dynamics simulations provide a powerful tool for studying the self-assembly structure and dynamic behavior of polymers, such as glass transition and entanglement dynamics. This review summarizes the coarse-grained models and methods in the dynamic simulations for polymers and their composite systems based on graphics processing unit(GPU) algorithms, and discusses the characteristics, applications, and advantages of different simulation methods. Based on recent studies in our group, the main progress of coarse-grained simulation methods in studying the structure, properties and physical mechanism of polymer materials is reviewed. It is anticipated to provide a reference for further development of coarse-grained simulation methods and software suitable for polymer research.