Detergent-mediated protein aggregation.

Because detergents are commonly used to solvate membrane proteins for structural evaluation, much attention has been devoted to assessing the conformational bias imparted by detergent micelles in comparison to the native environment of the lipid bilayer. Here, we conduct six 500-ns simulations of a system with >600,000 atoms to investigate the spontaneous self assembly of dodecylphosphocholine detergent around multiple molecules of the integral membrane protein PagP. This detergent formed equatorial micelles in which acyl chains surround the protein's hydrophobic belt, confirming existing models of the detergent solvation of membrane proteins. In addition, unexpectedly, the extracellular and periplasmic apical surfaces of PagP interacted with the headgroups of detergents in other micelles 85 and 60% of the time, respectively, forming complexes that were stable for hundreds of nanoseconds. In some cases, an apical surface of one molecule of PagP interacted with an equatorial micelle surrounding another molecule of PagP. In other cases, the apical surfaces of two molecules of PagP simultaneously bound a neat detergent micelle. In these ways, detergents mediated the non-specific aggregation of folded PagP. These simulation results are consistent with dynamic light scattering experiments, which show that, at detergent concentrations ≥600 mM, PagP induces the formation of large scattering species that are likely to contain many copies of the PagP protein. Together, these simulation and experimental results point to a potentially generic mechanism of detergent-mediated protein aggregation.

Construction and validation of a homology model of the human voltage-gated proton channel hHV1.

The topological similarity of voltage-gated proton channels (H(V)1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether H(V)1s have a similar structure. We present the construction and validation of a homology model of the human H(V)1 (hH(V)1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hH(V)1 by exploiting the homology of hH(V)1 with VSDs of K(+) and Na(+) channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hH(V)1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp(112). Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hH(V)1. Each Arg residue in the S4 helix of hH(V)1 was replaced by His to test accessibility using Zn(2+) as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg(211), the third Arg residue in the S4 helix in hH(V)1, remains accessible to the internal solution and is located near the charge transfer center, Phe(150).

Structure of saposin A lipoprotein discs.

The saposins are small, membrane-active proteins that exist in both soluble and lipid-bound states. Saposin A has roles in sphingolipid catabolism and transport and is required for the breakdown of galactosylceramide by ß-galactosylceramidase. In the absence of lipid, saposin A adopts a closed monomeric apo conformation typical of this family. To study a lipid-bound state of this protein, we determined the crystal structure of saposin A in the presence of detergent to 1.9 Å resolution. The structure reveals two chains of saposin A in an open conformation encapsulating 40 internally bound detergent molecules organized in a highly ordered bilayer-like hydrophobic core. The complex provides a high-resolution view of a discoidal lipoprotein particle in which all of the internalized acyl chains are resolved. Saposin A lipoprotein discs exhibit limited selectivity with respect to the incorporated lipid, and can solubilize phospholipids, sphingolipids, and cholesterol into discrete, monodisperse particles with mass of approximately 27 kDa. These discs may be the smallest possible lipoprotein structures that are stabilized by lipid self-assembly.

DNA and lipid bilayers: self-assembly and insertion.

DNA-lipid complexes are of biomedical importance as delivery vectors for gene therapy. To gain insight into the interactions of DNA with zwitterionic and cationic (dimyristoyltrimethylammonium propane (DMTAP)) lipids, we have used coarse-grained molecular dynamics simulations to study the self-assembly of DPPC and DPPC/DMTAP lipid bilayers in the presence of a DNA dodecamer. We observed the spontaneous formation of lipid bilayers from initial systems containing randomly placed lipids, water-counterions and DNA. In both the DPPC and DPPC/DMTAP simulations, the DNA molecule is located at the water-lipid headgroup interface, lying approximately parallel to the plane of the bilayer. We have also calculated the potential of mean force for transferring a DNA dodecamer through a DPPC/DMTAP bilayer. A high energetic barrier to DNA insertion into the hydrophobic core of the bilayer is observed. The DNA adopts a transmembrane orientation only in this region. Local bilayer deformation in the vicinity of the DNA molecule is observed, largely as a result of the DNA-DMTAP headgroup attraction.

Conformational change in an MFS protein: MD simulations of LacY.

Molecular dynamics simulations of lactose permease (LacY) in a phospholipid bilayer reveal the conformational dynamics of the protein. In inhibitor-bound simulations (i.e., those closest to the X-ray structure) the protein was stable, showing little conformational change over a 50 ns timescale. Movement of the bound inhibitor, TDG, to an alternative binding mode was observed, so that it interacted predominantly with the N-terminal domain and with residue E269 from the C-terminal domain. In multiple ligand-free simulations, a degree of domain closure occurred. This switched LacY to a state with a central cavity closed at both the intracellular and periplasmic ends. This may resemble a possible intermediate in the transport mechanism. Domain closure occurs by a combination of rigid-body movements of domains and of intradomain motions of helices, especially TM4, TM5, TM10, and TM11. A degree of intrahelix flexibility appears to be important in the conformational change.

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.

Modeling, docking, and simulation of the major facilitator superfamily.

X-ray structures are known for three members of the Major Facilitator Superfamily (MFS) of membrane transporter proteins, thus enabling the use of homology modeling to extrapolate to other MFS members. However, before employing such models for, e.g., mutational or docking studies, it is essential to develop a measure of their quality. To aid development of such metrics, two disparate MFS members (NupG and GLUT1) have been modeled. In addition, control models were created with shuffled sequences, to mimic poor quality homology models. These models and the template crystal structures have been examined in terms of both static and dynamic indicators of structural quality. Comparison of the behavior of modeled structures with the crystal structures in molecular dynamics simulations provided a metric for model quality. Docking of the inhibitor forskolin to GLUT1 and to a control model revealed significant differences, indicating that we may identify accurate models despite low sequence identity between target sequences and templates.

KcsA closed and open: modelling and simulation studies.

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height approximately 5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of approximately 3.3 A Calpha RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid ( < 5 ns) exit of all three K+ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.

Ion channel gating: insights via molecular simulations.

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.