Time-averaged order parameter restraints in molecular dynamics simulations

A method is described that allows experimental \(S^2\) order parameters to be enforced as a time-averaged quantity in molecular dynamics simulations. The two parameters that characterize time-averaged restraining, the memory relaxation time and the weight of the restraining potential energy term in the potential energy function used in the simulation, are systematically investigated based on two model systems, a vector with one end restrained in space and a pentapeptide. For the latter it is shown that the backbone N–H order parameter of individual residues can be enforced such that the spatial fluctuations of quantities depending on atomic coordinates are not significantly perturbed. The applicability to realistic systems is illustrated for the B3 domain of protein G in aqueous solution.     

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.     

Convergence of molecular dynamics simulations of membrane proteins

The central question in evaluating almost any result from a molecular dynamics simulation is whether the calculation has converged. Unfortunately, assessing the ergodicity of a single trajectory is very difficult to do. In this work, we assess the sampling of molecular dynamics simulations of the membrane protein rhodopsin by comparing the results from 26 independent trajectories, each run for 100 ns. By examining principal components and cluster populations, we show that rhodopsin's fluctuations are not well described by 100 ns of dynamics, and that the sampling is not fully converged even for individual loops. The results serve as a reminder of the caution required when interpreting molecular dynamics simulations of macromolecules.     

Combined use of homo- and heteronuclear coupling constants as restraints in molecular dynamics simulations.

A penalty function for scalar coupling constants has been applied in molecular dynamics simulations as an experimental constraint. The function is based on the difference between the coupling constant calculated from the dihedral angle and the experimentally measured coupling constant. The method is illustrated on a model cyclic pentapeptide for which 3JHN-H alpha and 3JHN-C beta, both about the phi backbone dihedral angle, have been measured. The function is efficient in producing structures consistent with the scalar couplings, but removed from the conformation observed in solution. This arises from the lack of J restraints for the psi dihedral angle. Simulations with both nuclear Overhauser effect (NOE) and J-coupling restraints illustrates small but significant differences from simulations using only NOEs.     

Structure refinement of membrane proteins via molecular dynamics simulations

A refinement protocol based on physics‐based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge‐based or implicit membrane‐based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid‐facing residues. Scoring with knowledge‐based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane‐based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Automated builder and database of protein/membrane complexes for molecular dynamics simulations

Molecular dynamics simulations of membrane proteins have provided deeper insights into their functions and interactions with surrounding environments at the atomic level. However, compared to solvation of globular proteins, building a realistic protein/membrane complex is still challenging and requires considerable experience with simulation software. Membrane Builder in the CHARMM-GUI website (http://www.charmm-gui.org) helps users to build such a complex system using a web browser with a graphical user interface. Through a generalized and automated building process including system size determination as well as generation of lipid bilayer, pore water, bulk water, and ions, a realistic membrane system with virtually any kinds and shapes of membrane proteins can be generated in 5 minutes to 2 hours depending on the system size. Default values that were elaborated and tested extensively are given in each step to provide reasonable options and starting points for both non-expert and expert users. The efficacy of Membrane Builder is illustrated by its applications to 12 transmembrane and 3 interfacial membrane proteins, whose fully equilibrated systems with three different types of lipid molecules (DMPC, DPPC, and POPC) and two types of system shapes (rectangular and hexagonal) are freely available on the CHARMM-GUI website. One of the most significant advantages of using the web environment is that, if a problem is found, users can go back and re-generate the whole system again before quitting the browser. Therefore, Membrane Builder provides the intuitive and easy way to build and simulate the biologically important membrane system.     

OmpT: molecular dynamics simulations of an outer membrane enzyme

Five molecular dynamics simulations (total duration >25 ns) have been performed on the Escherichia coli outer membrane protease OmpT embedded in a dimyristoylphosphatidylcholine lipid bilayer. Globally the protein is conformationally stable. Some degree of tilt of the beta-barrel is observed relative to the bilayer plane. The greatest degree of conformational flexibility is seen in the extracellular loops. A complex network of fluctuating H-bonds is formed between the active site residues, such that the Asp210-His212 interaction is maintained throughout, whereas His212 and Asp83 are often bridged by a water molecule. This supports a catalytic mechanism whereby Asp83 and His212 bind a water molecule that attacks the peptide carbonyl. A configuration yielded by docking calculations of OmpT simulation snapshots and a model substrate peptide Ala-Arg-Arg-Ala was used as the starting point for an extended Huckel calculation on the docked peptide. These placed the lowest unoccupied molecular orbital mainly on the carbon atom of the central C=O in the scissile peptide bond, thus favoring attack on the central peptide by the water held by residues Asp83 and His212. The trajectories of water molecules reveal exchange of waters between the intracellular face of the membrane and the interior of the barrel but no exchange at the extracellular mouth. This suggests that the pore-like region in the center of OmpT may enable access of water to the active site from below. The simulations appear to reveal the presence of specific lipid interaction sites on the surface of the OmpT barrel. This reveals the ability of extended MD simulations to provide meaningful information on protein-lipid interactions.     

Insights into buforin II membrane translocation from molecular dynamics simulations

Buforin II is a histone-derived antimicrobial peptide that readily translocates across lipid membranes without causing significant membrane permeabilization. Previous studies showed that mutating the sole proline of buforin II dramatically decreases its translocation. As well, researchers have proposed that the peptide crosses membranes in a cooperative manner by forming transient toroidal pores. This paper reports molecular dynamics simulations designed to investigate the structure of buforin II upon membrane entry and evaluate whether the peptide is able to form toroidal pore structures. These simulations showed a relationship between protein–lipid interactions and increased structural deformations of the buforin N-terminal region promoted by proline. Moreover, simulations with multiple peptides show how buforin II can embed deeply into membranes and potentially form toroidal pores. Together, these simulations provide structural insight into the translocation process for buforin II in addition to providing more general insight into the role proline can play in antimicrobial peptides.     

Comparative molecular dynamics simulations of amphotericin B-cholesterol/ergosterol membrane channels

Amphotericin B (AmB) is a very effective anti-fungal polyene macrolide antibiotic whose usage is limited by its toxicity. Lack of a complete understanding of AmB's molecular mechanism has impeded attempts to design less toxic AmB derivatives. The antibiotic is known to interact with sterols present in the cell membrane to form ion channels that disrupt membrane function. The slightly higher affinity of AmB toward ergosterol (dominant sterol in fungal cells) than cholesterol (mammalian sterol) is regarded as the most essential factor on which antifungal chemotherapy is based. To study these differences at the molecular level, two realistic model membrane channels containing molecules of AmB, sterol (cholesterol or ergosterol), phospholipid, and water were studied by molecular dynamics (MD) simulations. Comparative analysis of the simulation data revealed that the sterol type has noticeable effect on the properties of AmB membrane channels. In addition to having a larger size, the AmB channel in the ergosterol-containing membrane has a more pronounced pattern of intermolecular hydrogen bonds. The interaction between the antibiotic and ergosterol is more specific than between the antibiotic and cholesterol. These observed differences suggest that the channel in the ergosterol-containing membrane is more stable and, due to its larger size, would have a higher ion conductance. These observations are in agreement with experiments.     

Self-assembling of peptide/membrane complexes by atomistic molecular dynamics simulations

Model biological membranes consisting of peptide/lipid-bilayer complexes can nowadays be studied by classical molecular dynamics (MD) simulations at atomic detail. In most cases, the simulation starts with an assumed state of a peptide in a preformed bilayer, from which equilibrium configurations are difficult to obtain due to a relatively slow molecular diffusion. As an alternative, we propose an extension of reported work on the self-organization of unordered lipids into bilayers, consisting of including a peptide molecule in the initial random configuration to obtain a membrane-bound peptide simultaneous to the formation of the lipid bilayer. This strategy takes advantage of the fast reorganization of lipids, among themselves and around the peptide, in an aqueous environment. Model peptides of different hydrophobicity, CH3-CO-W2L18W2-NH2 (WL22) and CH3-CO-W2A18W2-NH2 (WA22), in dipalmitoyl-phosphatidylcholine (DPPC) are used as test cases. In the equilibrium states of the peptide/membrane complexes, achieved in time ranges of 50-100 ns, the two peptides behave as expected from experimental and theoretical studies. The strongly hydrophobic WL22 is inserted in a transmembrane configuration and the marginally apolar, alanine-based WA22 is found in two alternative states: transmembrane inserted or parallel to the membrane plane, embedded close to the bilayer interface, with similar stability. This shows that the spontaneous assembly of peptides and lipids is an unbiased and reliable strategy to produce and study models of equilibrated peptide/lipid complexes of unknown membrane-binding mode and topology.     

Conformational regions of Boc-Ala-Aib-Ala-OMe. Sampling with molecular dynamics simulations using time averaging of distance restraints.

The method of time averaging of distance restraints in molecular dynamics simulations is applied to Boc-Ala-Aib-Ala-OMe in order to demonstrate the improved sampling properties of this method compared to conventional distance restraining. Two conformational regions, beta-turn type II and gamma-turn, are seen during MD runs at a simulation temperature of 500 K, while in simulations with conventional distance restraining, no conformational transitions could be observed for temperatures up to 1000 K.     

On using time-averaging restraints in molecular dynamics simulation

Introducing experimental values as restraints into molecular dynamics (MD) simulations to bias the values of particular molecular properties, such as nuclear Overhauser effect intensities or distances, 3J coupling constants, chemical shifts or crystallographic structure factors, towards experimental values is a widely used structure refinement method. To account for the averaging of experimentally derived quantities inherent in the experimental techniques, time-averaging restraining methods may be used. In the case of structure refinement using 3J coupling constants from NMR experiments, time-averaging methods previously proposed can suffer from large artificially induced structural fluctuations. A modified time-averaged restraining potential energy function is proposed which overcomes this problem. The different possible approaches are compared using stochastic dynamics simulations of antamanide, a cyclic peptide of ten residues.     

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.     

An NMR database for simulations of membrane dynamics

Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental 13C-1H and 2H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.     

OmpA: Gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of ∼ 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.     

Creating realistic silica nanopores for molecular dynamics simulations

We present a detailed approach to create realistic silica pores for computer simulations especially molecular dynamics (MD) simulations. These pores are essential for all different kinds of simulations with liquids in silica confinements. Despite wide use of silica pores in simulations, a detailed documentation how to create these pores for simulations still lacks. This issue is of high significance because with the help of this paper every researcher can build own silica pores with desired geometries and is not stick to already existing pores. We discuss problems that might occur during the whole process and how to solve these problems. So far more than 3 different silica pores have been created with this method and used successfully as confinement material in MD simulations.     

Higher order multipole moments for molecular dynamics simulations

In conventional force fields, the electrostatic potential is represented by atom-centred point charges. This choice is in principle arbitrary, but technically convenient. Point charges can be understood as the first term of multipole expansions, which converge with an increasing number of terms towards the accurate representation of the molecular potential given by the electron density distribution. The use of multipole expansions can therefore improve the force field accuracy. Technically, the implementation of atomic multipoles is more involved than the use of point charges. Important points to consider are the orientation of the multipole moments during the trajectory, conformational dependence of the atomic moments and stability of the simulations which are discussed here.