Setting up and optimization of membrane protein simulations.
In this paper we describe a method for setting up an atomistic simulation of a membrane protein in a hydrated lipid bilayer and report the effect of differing electrostatic parameters on the drift in the protein structure during the subsequent simulation. The method aims to generate a suitable cavity in the interior of a lipid bilayer, using the solvent-accessible surface of the protein as a template, during the course of a short steered molecular dynamics simulation of a solvated lipid membrane. This is achieved by a two-stage process: firstly, lipid molecules whose headgroups are inside a cylindrical volume equivalent to that defined by the protein surface are removed; then the protein-lipid interface is optimized by applying repulsive forces perpendicular to the protein surface, and of gradually increased magnitude, to the remaining lipid atoms inside the volume occupied by the protein surface until it is emptied. The protein itself may then be inserted. Using the bacterial membrane proteins KcsA and FhuA as test cases, we show how the method achieves the formation of a suitable cavity in the interior of a dimyristoylphosphatidylcholine lipid bilayer without perturbing the configuration of the non-interfacial regions of the previously equilibrated lipid bilayer, even in cases of membrane proteins with irregular geometrical shapes. In addition, we compare subsequent simulations in which the long-range electrostatic interactions are treated via either a cut-off or a particle-mesh Ewald method. The results show that the drift from the initial structure is less in the latter case, especially for KcsA and for the non-core secondary structural elements (i.e. surface loops) of both proteins.
Faraldo-Gómez
José D
JD
Laboratory of Molecular Biophysics, Department of Biochemistry, The Rex Richards Building, University of Oxford, South Parks Road, Oxford OX1 3QU, UK.
Smith
Graham R
GR
Sansom
Mark S P
MS
D001425:Bacterial Outer Membrane Proteins / Q000737:chemistry*; D001426:Bacterial Proteins / Q000737:chemistry*; D003196:Computer Graphics; D003198:Computer Simulation; D004563:Electrochemistry / Q000379:methods*; D029968:Escherichia coli Proteins / Q000737:chemistry*; D008051:Lipid Bilayers / Q000737:chemistry*; D046911:Macromolecular Substances; D008565:Membrane Proteins / Q000737:chemistry; D008567:Membranes, Artificial; D008958:Models, Molecular; D015221:Potassium Channels / Q000737:chemistry*; D011487:Protein Conformation; D011786:Quality Control; D011991:Receptors, Virus / Q000737:chemistry*; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D055672:Static Electricity; D014867:Water / Q000737:chemistry*
D001425:Bacterial Outer Membrane Proteins; D001426:Bacterial Proteins; D029968:Escherichia coli Proteins; C050900:FhuA protein, E coli; D008051:Lipid Bilayers; D046911:Macromolecular Substances; D008565:Membrane Proteins; D008567:Membranes, Artificial; D015221:Potassium Channels; D011991:Receptors, Virus; C107790:prokaryotic potassium channel; D014867:Water