Search

ABC transporters are integral membrane proteins which couple the energy of ATP hydrolysis to the translocation of solutes across cell membranes. BtuCD is a approximately 1100-residue protein found in the inner membrane of Gram-negative bacteria which transports vitamin B12. Vitamin B12 is bound in the periplasm by BtuF, which delivers the solute to the periplasmic entrance of the transporter protein complex BtuCD. Molecular dynamics simulations of the BtuCD and BtuCDF complexes (in a lipid bilayer) and of the isolated BtuD and BtuF proteins (in water) have been used to explore the conformational dynamics of this complex transport system. Overall, seven simulations have been performed, with and without bound ATP, corresponding to a total simulation time of 0.1 micros. Binding of ATP drives closure of the nucleotide-binding domains (NBDs) in BtuD in a symmetrical fashion, but not in BtuCD. It seems that ATP constrains the flexibility of the NBDs in BtuCD such that their closure may only occur upon binding of BtuF to the complex. Upon introduction of BtuF, and concomitant with NBD association, one ATP-binding site displays a closure, while the opposite site remains relatively unchanged. This asymmetry may reflect an initial step in the "alternating hydrolysis" mechanism and is consistent with measurements of nucleotide-binding stoichiometries. Principal components analysis of the simulation of BtuCD reveals motions that are comparable to those suggested in current transport models.

Comparative molecular dynamics (MD) simulations enable us to explore the conformational dynamics of the active sites of distantly related enzymes. We have used the BioSimGrid (http://www.biosimgrid.org) database to facilitate such a comparison. Simulations of four enzymes were analyzed. These included three hydrolases and a transferase, namely acetylcholinesterase, outer-membrane phospholipase A, outer-membrane protease T, and PagP (an outer-membrane enzyme which transfers a palmitate chain from a phospholipid to lipid A). A set of 17 simulations were analyzed corresponding to a total of approximately 0.1 micros simulation time. A simple metric for active-site integrity was used to demonstrate the existence of clusters of dynamic conformational behaviour of the active sites. Small (i.e. within a cluster) fluctuations appear to be related to the function of an enzymatically active site. Larger fluctuations (i.e. between clusters) correlate with transitions between catalytically active and inactive states. Overall, these results demonstrate the potential of a comparative MD approach to analysis of enzyme function. This approach could be extended to a wider range of enzymes using current high throughput MD simulation and database methods.

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.

The precise nature of redox modulation of N-methyl-d-aspartate (NMDA) receptors is still unclear, although it is thought to be related to the formation and breaking of disulfide bonds. Recent structural data demonstrated the way in which disulfide bonds in the ligand-binding core of the NR1 subunit are arranged. However, the structures were not able to reconcile existing experimental data that examined the effects of mutating these cysteine residues. We have used molecular dynamics (MD) simulations of a series of in silico mutations to try and address this in terms of the current structure of the NR1 ligand-binding domain. A double mutation that removes the disulfide bridge between C744 and C798 gives rise to greater interlobe mobility which was predicted from the crystal structure information but, unexpectedly, also appears to predispose the receptor toward greater flexibility in the hinge region. Removal of the disulfide bond between C454 and C420 did not show any appreciable difference from the "wild-type" simulation, suggesting that removal of this would not change receptor properties, which is in agreement with experimental findings. Furthermore, the position of the C454 side chain could be characterized into discrete rotamers, which may reflect the observation of alternative density in the crystal structure for this residue. Simulations in which two of the disulfide bonds are removed via mutations to alanine (C420A and C436A) resulted in a tendency of the protein to adopt a partially closed conformation.

Coarse-grained molecular dynamics simulations are used to explore the interaction with a phospholipid bilayer of the voltage sensor (VS) domain and the S4 helix from the archaebacterial voltage-gated potassium (Kv) channel KvAP. Multiple 2-mus self-assembly simulations reveal that the isolated S4 helix may adopt either interfacial or transmembrane (TM) locations with approximately equal probability. In the TM state, the insertion of the voltage-sensing region of S4 is facilitated via local bilayer deformation that, combined with side chain "snorkeling," enables its Arg side chains to interact with lipid headgroups and water. Multiple 0.2-mus self-assembly simulations of the VS domain are also performed, along with simulations of MscL and KcsA, to permit comparison with more "canonical" integral membrane protein structures. All three stably adopt a TM orientation within a bilayer. For MscL and KcsA, there is no significant bilayer deformation. In contrast, for the VS, there is considerable local deformation, which is again primarily due to the lipid-exposed S4. It is shown that for both the VS and isolated S4 helix, the positively charged side chains of S4 are accommodated within the membrane through a combination of stabilizing interactions with lipid glycerol and headgroup regions, water, and anionic side chains. Our results support the possibility that bilayer deformation around key gating charge residues in Kv channels may result in "focusing" of the electrostatic field, and indicate that, when considering competing models of voltage-sensing, it is essential to consider the dynamics and structure of not only the protein but also of the local lipid environment.

The nature of voltage sensing by voltage-activated ion channels is a key problem in membrane protein structural biology. The way in which the voltage-sensor (VS) domain interacts with its membrane environment remains unclear. In particular, the known structures of Kv channels do not readily explain how a positively charged S4 helix is able to stably span a lipid bilayer. Extended (2 x 50 ns) molecular dynamics simulations of the high-resolution structure of the isolated VS domain from the archaebacterial potassium channel KvAP, embedded in zwitterionic and in anionic lipid bilayers, have been used to explore VS/lipid interactions at atomic resolution. The simulations reveal penetration of water into the center of the VS and bilayer. Furthermore, there is significant local deformation of the lipid bilayer by interactions between lipid phosphate groups and arginine side chains of S4. As a consequence of this, the electrostatic field is "focused" across the center of the bilayer.

The glycophorin helix dimer is a paradigm for the exploration of helix-helix interactions in integral membrane proteins. Two NMR structures of the dimer are known, one in a detergent micelle and one in a lipid bilayer. Multiple (4 x 50 ns) molecular dynamics simulations starting from each of the two NMR structures, with each structure in either a dodecyl phosphocholine (DPC) micelle or a dimyristoyl phosphatidylcholine (DMPC) bilayer, have been used to explore the conformational dynamics of the helix dimer. Analysis of the helix-helix interaction, mediated by the GxxxG sequence motif, suggests convergence of the simulations to a common model. This is closer to the NMR structure determined in a bilayer than to micelle structure. The stable dimer interface in the final simulation model is characterized by (i) Gly/Gly packing and (ii) Thr/Thr interhelix H-bonds. These results demonstrate the ability of extended molecular dynamics simulations in a lipid bilayer environment to refine membrane protein structures or models derived from experimental data obtained in protein/detergent micelles.

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.

The cytoplasmic membrane protein TatB is an essential component of the Escherichia coli twin-arginine (Tat) protein translocation pathway. Together with the TatC component it forms a complex that functions as a membrane receptor for substrate proteins. Structural predictions suggest that TatB is anchored to the membrane via an N-terminal transmembrane alpha-helix that precedes an amphipathic alpha-helical section of the protein. From truncation analysis it is known that both these regions of the protein are essential for function. Here we construct 31 unique cysteine substitutions in the first 42 residues of TatB. Each of the substitutions results in a TatB protein that is competent to support Tat-dependent protein translocation. Oxidant-induced disulfide cross-linking shows that both the N-terminal and amphipathic helices form contacts with at least one other TatB protomer. For the transmembrane helix these contacts are localized to one face of the helix. Molecular modeling and molecular dynamics simulations provide insight into the possible structural basis of the transmembrane helix interactions. Using variants with double cysteine substitutions in the transmembrane helix, we were able to detect cross-links between up to five TatB molecules. Protein purification showed that species containing at least four cross-linked TatB molecules are found in correctly assembled TatBC complexes. Our results suggest that the transmembrane helices of TatB protomers are in the center rather than the periphery of the TatBC complex.

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane beta-barrel containing a central alpha-helix. The role of this alpha-helix, and of the conformational dynamics of the beta-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane beta-barrel is stable even in the absence of the alpha-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a 'plug'. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.

Contemporary structural biology has an increased emphasis on high-throughput methods. Biomolecular simulations can add value to structural biology via the provision of dynamic information. However, at present there are no agreed measures for the quality of biomolecular simulation data. In this Letter, we suggest suitable measures for the quality assurance of molecular dynamics simulations of biomolecules. These measures are designed to be simple, fast, and general. Reporting of these measures in simulation papers should become an expected practice, analogous to the reporting of comparable quality measures in protein crystallography. We wish to solicit views and suggestions from the simulation community on methods to obtain reliability measures from molecular-dynamics trajectories. In a database which provides access to previously obtained simulations [Formula: see text] for example BioSimGrid ( http://www.biosimgrid.org/ ) [Formula: see text] the user needs to be confident that the simulation trajectory is suitable for further investigation. This can be provided by the simulation quality measures which a user would examine prior to more extensive analyses.

SecY is the central channel protein of the SecYEbeta translocon, the structure of which has been determined by X-ray diffraction. Extended (15 ns) MD simulations of the isolated SecY protein in a phospholipid bilayer have been performed to explore the relationship between protein flexibility and the mechanisms of channel gating. In particular, principal components analysis of the simulation trajectory has been used to probe the intrinsic flexibility of the isolated SecY protein in the absence of the gamma-subunit (SecE) clamp. Analysis and visualization of the principal eigenvectors support a "plug and clamshell" model of SecY channel gating. The simulation results also indicate that hydrophobic gating at the central pore ring prevents leakage of water and ions through the channel in the absence of a translocating peptide.

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.

Mistic is an unusual membrane protein from Bacillus subtilis. It appears to fold and insert autonomously into a lipid bilayer and has been suggested as a tool that aids the targeting of eukaryotic membrane proteins to bacterial membranes. The NMR structure of Mistic in detergent (LDAO) micelles has revealed it to be a four alpha-helix bundle. From a structural perspective, Mistic does not resemble other membrane proteins. Its external surface is not very hydrophobic, and standard methods do not predict any of its helices to be in the transmembrane orientation. Molecular dynamics simulations (simulation times approximately 30 ns) in water and in detergent micelles have been used to explore the conformational stability of Mistic as a function of its environment. In water, the protein is stable, exhibiting no significant change in fold on a 30 ns time scale. In contrast, in three simulations in detergent micelles, the partial unfolding of Mistic occurred, whereby the H4 helix drifted away from the H1-H3 core. This was due to the penetration of detergent molecules between H4 and the remainder of the protein. This is unlike the behavior of several other membrane proteins, both alpha-helix bundles and beta-barrels, in comparable detergent micelle simulations. The unfolding of H4 from the H1-H3 core of Mistic could be partially reversed by a simulation in which the detergent molecules were removed, and the unfolded protein was simulated in water. These results suggest that Mistic may not be a stable integrated membrane protein but rather that it may undergo a conformational change upon interaction with a membrane or membrane-like environment.

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the 'Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid and serotonin. Cryo-electron microscopy has yielded a three-dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height about 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united-atom force field for the protein instead of an all-atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.

Combining structure determinations from nuclear magnetic resonance (NMR) data and molecular dynamics simulations (MD) under the same environmental conditions revealed a startling asymmetry in the intrinsic conformational stability of secondary structure in the transmembrane domain of lactose permease (LacY). Eleven fragments, corresponding to transmembrane segments (TMs) of LacY, were synthesized, and their secondary structure in solution was determined by NMR. Eight of the TMs contained significant regions of helical structure. MD simulations, both in DMSO and in a DMPC bilayer, showed sites of local stability of helical structure in these TMs, punctuated by regions of conformational instability, in substantial agreement with the NMR data. Mapping the stable regions onto the crystal structure of LacY reveals a marked asymmetry, contrasting with the pseudosymmetry in the static structure: the secondary structure in the C-terminal half is more stable than in the N-terminal half. The relative stability of secondary structure is likely exploited in the transport mechanism of LacY. Residues supporting proton conduction are in more stable regions of secondary structure, while residues key to substrate binding are found in considerably unstable regions of secondary structure.

Molecular dynamics (MD) simulations are used to explore the dynamics of a membrane protein in its crystal environment. A 50-ns-duration simulation (at a temperature of 300 K) is performed for the crystallographic unit cell of the bacterial outer membrane protein OmpA. The unit cell contains four protein molecules, plus detergent molecules and water. An excellent correlation between simulated and experimental values of crystallographic B factors is observed. Effectively, 0.2 micros of protein trajectories are obtained, allowing a critical assessment of simulation quality. Some deficiency in conformational sampling is demonstrated, but averaging over multiple trajectories improves this limitation. The previously undescribed structure and dynamics of detergent molecules in a unit cell are reported here, providing insight into the interactions important in the formation and stabilization of the crystalline environment at room temperature. In particular, we show that at room temperature the detergent molecules form a dynamic, extended micellar structure spreading over adjacent OmpA monomers within the crystal.

The pH-sensitive renal potassium channel Kir1.1 is important for K+ homeostasis. Disruption of the pH-sensing mechanism causes type II Bartter syndrome. The pH sensor is thought to be an anomalously titrated lysine residue (K80) that interacts with two arginine residues as part of an 'RKR triad'. We show that a Kir1.1 orthologue from Fugu rubripes lacks this lysine and yet is still highly pH sensitive, indicating that K80 is not the H+ sensor. Instead, K80 functionally interacts with A177 on transmembrane domain 2 at the 'helix-bundle crossing' and controls the ability of pH-dependent conformational changes to induce pore closure. Although not required for pH inhibition, K80 is indispensable for the coupling of pH gating to the extracellular K+ concentration, explaining its conservation in most Kir1.1 orthologues. Furthermore, we demonstrate that instead of interacting with K80, the RKR arginine residues form highly conserved inter- and intra-subunit interactions that are important for Kir channel gating and influence pH sensitivity indirectly.

The mechanism of partial agonism at N-methyl-D-aspartate receptors is an unresolved issue, especially with respect to the role of protein dynamics. We have performed multiple molecular dynamics simulations (7 x 20 ns) to examine the behavior of the ligand-binding core of the NR1 subunit with a series of ligands. Our results show that water plays an important role in stabilizing different conformations of the core and how a closed cleft conformation of the protein might be stabilized in the absence of ligands. In the case of ligand-bound simulations with both full and partial agonists, we observed that ligands within the binding cleft may undergo distinct conformational changes, without grossly influencing the degree of cleft closure within the ligand-binding domain. In agreement with recently published crystallographic data, we also observe similar changes in backbone torsions corresponding to the hinge region between the two lobes for the partial agonist, D-cycloserine. This observation rationalizes the classification of D-cycloserine as a partial agonist and should provide a basis with which to predict partial agonism in this class of receptor by analyzing the behavior of these torsions with other potential ligands.

OprF is a major outer membrane protein from Pseudomonas aeruginosa, a homolog of OmpA from Escherichia coli. The N-terminal domains of both proteins have been demonstrated to form low conductance channels in lipid bilayers. Homology models, consisting of an eight-stranded beta-barrel, of the N-terminal domain OprF have been constructed based on the crystal structure of the corresponding domain from E. coli OmpA. OprF homology models have been evaluated via a set (6 x 10 ns) of simulations of the beta-barrel embedded within a solvated dimyristoyl-phosphatidylcholine (DMPC) bilayer. The conformational stability of the models is similar to that of the crystal structure of OmpA in comparable simulations. There is a degree of water penetration into the pore-like center of the OprF barrel. The presence of an acidic/basic (E8/K121) side-chain interaction within the OprF barrel may form a "gate" able to close/open a central pore. Lipid-protein interactions within the simulations were analyzed and revealed that aromatic side-chains (Trp, Tyr) of OprF interact with lipid headgroups. Overall, the behavior of the OprF model in simulations supports the suggestion that this molecule is comparable to OmpA. The simulations help to explain the mechanism of formation of low conductance pores within the outer membrane., (2006 Wiley-Liss, Inc.)

Ion channel-forming peptides enable us to study the conformational dynamics of a transmembrane helix as a function of sequence and environment. Molecular dynamics simulations are used to study the conformation and dynamics of three 22-residue peptides derived from the second transmembrane domain of the glycine receptor (NK4-M2GlyR-p22). Simulations are performed on the peptide in four different environments: trifluoroethanol/water; SDS micelles; DPC micelles; and a DMPC bilayer. A hierarchy of alpha-helix stabilization between the different environments is observed such that TFE/water < micelles < bilayers. Local clustering of trifluoroethanol molecules around the peptide appears to help stabilize an alpha-helical conformation. Single (S22W) and double (S22W,T19R) substitutions at the C-terminus of NK4-M2GlyR-p22 help to stabilize a helical conformation in the micelle and bilayer environments. This correlates with the ability of the W22 and R19 side chains to form H-bonds with the headgroups of lipid or detergent molecules. This study provides a first atomic resolution comparison of the structure and dynamics of NK4-M2GlyR-p22 peptides in membrane and membrane-mimetic environments, paralleling NMR and functional studies of these peptides.

Analysis of the crystal structures of the intact voltage-sensitive potassium channel KvAP (from Aeropyrum pernix) and Kv1.2 (from rat brain), along with the isolated voltage sensor (VS) domain from KvAP, raises the question of the exact nature of the voltage-sensing conformational change that triggers activation of Kv and related voltage-gated channels. Molecular dynamics simulations of the isolated VS of KvAP in a detergent micelle environment at two different temperatures (300 K and 368 K) have been used to probe the intrinsic flexibility of this domain on a tens-of-nanoseconds timescale. The VS contains a positively charged (S4) helix which is packed against a more hydrophobic S3 helix. The simulations at elevated temperature reveal an intrinsic flexibility/conformational instability of the S3a region (i.e., the C-terminus of the S3 helix). It is also evident that the S4 helix undergoes hinge bending and swiveling about its central I130 residue. The conformational instability of the S3a region facilitates the motion of the N-terminal segment of S4 (i.e., S4a). These simulations thus support a gating model in which, in response to depolarization, an S3b-S4a "paddle" may move relative to the rest of the VS domain. The flexible S3a region may in turn act to help restore the paddle to its initial conformation upon repolarization.

Interactions of lipids are central to the folding and stability of membrane proteins. Coarse-grained molecular dynamics simulations have been used to reveal the mechanisms of self-assembly of protein/membrane and protein/detergent complexes for representatives of two classes of membrane protein, namely, glycophorin (a simple alpha-helical bundle) and OmpA (a beta-barrel). The accuracy of the coarse-grained simulations is established via comparison with the equivalent atomistic simulations of self-assembly of protein/detergent micelles. The simulation of OmpA/bilayer self-assembly reveals how a folded outer membrane protein can be inserted in a bilayer. The glycophorin/bilayer simulation supports the two-state model of membrane folding, in which transmembrane helix insertion precedes dimer self-assembly within a bilayer. The simulations also suggest that a dynamic equilibrium exists between the glycophorin helix monomer and dimer within a bilayer. The simulated glycophorin helix dimer is remarkably close in structure to that revealed by NMR. Thus, coarse-grained methods may help to define mechanisms of membrane protein (re)folding and will prove suitable for simulation of larger scale dynamic rearrangements of biological membranes.

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG >> PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.

Ionotropic glutamate receptors mediate fast synaptic transmission in the mammalian central nervous system and play an important role in many different functions, including memory and learning. They have also been implicated in a variety of neuropathologies and as such have generated widespread interest in their structure and function. Molecular Dynamics simulations (5 x 20 ns) of the ligand-binding core of the GluR2 glutamate receptor were performed. Through simulations of both wild type and the L650T mutant, we show that the degree of protein flexibility can be correlated with the extent to which the binding cleft is open. In agreement with recent experiments, the simulations of kainate with the wild-type construct show a slight increase in beta-sheet content that we are able to localize to two specific regions. During one simulation, the protein made a transition from an open-cleft conformation to a closed-cleft conformation. This closed cleft conformation closely resembles the closed-cleft crystal structure, thus indicating a potential pathway for conformational change associated with receptor activation. Analysis of the binding pocket suggests that partial agonists possess a greater degree of flexibility within the pocket that may help to explain why they are less efficient at opening the channel than full agonists. Examination of water molecules surrounding the ligands reveals that mobility in distinct subsites can be a discriminator between full and partial agonism and will be an important consideration in the design of drugs against these receptors.

ATP-sensitive potassium (K(ATP)) channels conduct potassium ions across cell membranes and thereby couple cellular energy metabolism to membrane electrical activity. Here, we report the heterologous expression and purification of a functionally active K(ATP) channel complex composed of pore-forming Kir6.2 and regulatory SUR1 subunits, and determination of its structure at 18 A resolution by single-particle electron microscopy. The purified channel shows ATP-ase activity similar to that of ATP-binding cassette proteins related to SUR1, and supports Rb(+) fluxes when reconstituted into liposomes. It has a compact structure, with four SUR1 subunits embracing a central Kir6.2 tetramer in both transmembrane and cytosolic domains. A cleft between adjacent SUR1s provides a route by which ATP may access its binding site on Kir6.2. The nucleotide-binding domains of adjacent SUR1 appear to interact, and form a large docking platform for cytosolic proteins. The structure, in combination with molecular modelling, suggests how SUR1 interacts with Kir6.2.

Proteins possessing the same fold may undergo similar motions, particularly if these motions involve large conformational transitions. The increasing amounts of structural data provide a useful starting point with which to test this hypothesis. We have performed a total of 0.29 micros of molecular dynamics across a series of proteins within the same fold family (periplasmic binding proteinlike) in order to address to what extent similarity of motion exists. Analysis of the local conformational space on these timescales (10-20 ns) revealed that the behavior of the proteins could be readily distinguished between an apo-state and a ligand-bound state. Moreover, analysis of the root-mean-square fluctuations reveals that the presence of the ligand exerts a stabilizing effect on the protein, with similar motions occurring, but with reduced magnitude. Furthermore, the conformational space in the presence of the ligand appears to be dictated by sequence but not by the type of ligand present. In contrast, apo-simulations showed considerable overlap of conformational space across the fold as a result of their ability to undergo larger fluctuations. Indeed, we observed several transitions from different simulations between states corresponding to the closed-cleft and open-cleft forms of the fold, with the predominant motions being conserved across the different proteins. Thus, large-scale conformational changes do indeed appear to be conserved across this fold architecture, but smaller conformational motions appear to reflect the differences in sequence and local fold., (Proteins 2005. 2005 Wiley-Liss, Inc.)

The effect of a short loop connecting two transmembrane alpha-helices was studied using molecular dynamics simulations. Helices F and G from bacteriorhodopsin and two corresponding polyalanine helices were embedded in octane and POPC membranes in a transmembrane configuration both with and without the inter-helical loop. The results indicate that the membrane environment and the sequence of the loop are more influential on the dynamics and structure of the motif than the presence of a loop as such, at least for the time-scales investigated. The four residues in the FG loop are stabilized by four hydrogen bonds. These hydrogen bonds are not present in the polyalanine loop, causing it to be more flexible than the FG loop. This effect was observed independently of the protein environment, stressing the importance of the sequence. The structural analysis indicates that the loop has weak stabilizing properties in all environments. The stabilization due to the presence of the loop was strongest in a simulation of the FG fragment in a membrane-mimetic octane slab. In the simulations of the helix-loop-helix motif embedded in an explicit lipid bilayer model, the lipid bilayer interface compensates to a large extent for the absence of the loop.

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

Molecular dynamics simulations may be used to probe the interactions of membrane proteins with lipids and with detergents at atomic resolution. Examples of such simulations for ion channels and for bacterial outer membrane proteins are described. Comparison of simulations of KcsA (an alpha-helical bundle) and OmpA (a beta-barrel) reveals the importance of two classes of side chains in stabilizing interactions with the head groups of lipid molecules: (i) tryptophan and tyrosine; and (ii) arginine and lysine. Arginine residues interacting with lipid phosphate groups play an important role in stabilizing the voltage-sensor domain of the KvAP channel within a bilayer. Simulations of the bacterial potassium channel KcsA reveal specific interactions of phosphatidylglycerol with an acidic lipid-binding site at the interface between adjacent protein monomers. A combination of molecular modelling and simulation reveals a potential phosphatidylinositol 4,5-bisphosphate-binding site on the surface of Kir6.2.

Interactions between membrane proteins and detergents are important in biophysical and structural studies and are also biologically relevant in the context of folding and transport. Despite a paucity of high-resolution data on protein-detergent interactions, novel methods and increased computational power enable simulations to provide a means of understanding such interactions in detail. Simulations have been used to compare the effect of lipid or detergent on the structure and dynamics of membrane proteins. Moreover, some of the longest and most complex simulations to date have been used to observe the spontaneous formation of membrane protein-detergent micelles. Common mechanistic steps in the micelle self-assembly process were identified for both alpha-helical and beta-barrel membrane proteins, and a simple kinetic mechanism was proposed. Recently, simplified (i.e. coarse-grained) models have been utilized to follow long timescale transitions in membrane protein-detergent assemblies.

Our goal was to assess the relationship between membrane protein quality, output from protein quality checkers and output from molecular dynamics (MD) simulations. Membrane transport proteins are essential for a wide range of cellular processes. Structural features of integral membrane proteins are still under-explored due to experimental limitations in structure determination. Computational techniques can be used to exploit biochemical and medium resolution structural data, as well as sequence homology to known structures, and enable us to explore the structure-function relationships in several transmembrane proteins. The quality of the models produced is vitally important to obtain reliable predictions. An examination of the relationship between model stability in molecular dynamics (MD) simulations derived from RMSD (root mean squared deviation) and structure quality assessment from various protein quality checkers was undertaken. The results were compared to membrane protein structures, solved at various resolution, by either X-ray or electron diffraction techniques. The checking programs could predict the potential success of MD in making functional conclusions. MD stability was shown to be a good indicator for the quality of structures. The quality was also shown to be dependent on the resolution at which the structures were determined.

The recently determined structure of a mammalian voltage-gated potassium channel has important implications for our understanding of voltage-sensing and gating mechanisms in channels. It is also the first crystal structure of an overexpressed eukaryotic membrane protein.

Biomolecular computer simulations are now widely used not only in an academic setting to understand the fundamental role of molecular dynamics on biological function, but also in the industrial context to assist in drug design. In this paper, two applications of Grid computing to this area will be outlined. The first, involving the coupling of distributed computing resources to dedicated Beowulf clusters, is targeted at simulating protein conformational change using the Replica Exchange methodology. In the second, the rationale and design of a database of biomolecular simulation trajectories is described. Both applications illustrate the increasingly important role modern computational methods are playing in the life sciences.

Transport by ABC proteins requires a cycle of ATP-driven conformational changes of the nucleotide binding domains (NBDs). We compare three molecular dynamics simulations of dimeric MJ0796: with ATP was present at both NBDs; with ATP at one NBD but ADP at the other; and without any bound ATP. In the simulation with ATP present at both NBDs, the dimeric protein interacts with the nucleotides in a symmetrical manner. However, if ADP is present at one binding site then both NBD-NBD and protein-ATP interactions are enhanced at the opposite site.

The prediction of transmembrane (TM) helices plays an important role in the study of membrane proteins, given the relatively small number (approximately 0.5% of the PDB) of high-resolution structures for such proteins. We used two datasets (one redundant and one non-redundant) of high-resolution structures of membrane proteins to evaluate and analyse TM helix prediction. The redundant (non-redundant) dataset contains structure of 434 (268) TM helices, from 112 (73) polypeptide chains. Of the 434 helices in the dataset, 20 may be classified as 'half-TM' as they are too short to span a lipid bilayer. We compared 13 TM helix prediction methods, evaluating each method using per segment, per residue and termini scores. Four methods consistently performed well: SPLIT4, TMHMM2, HMMTOP2 and TMAP. However, even the best methods were in error by, on average, about two turns of helix at the TM helix termini. The best and worst case predictions for individual proteins were analysed. In particular, the performance of the various methods and of a consensus prediction method, were compared for a number of proteins (e.g. SecY, ClC, KvAP) containing half-TM helices. The difficulties of predicting half-TM helices suggests that current prediction methods successfully embody the two-state model of membrane protein folding, but do not accommodate a third stage in which, e.g., short helices and re-entrant loops fold within a bundle of stable TM helices.

ATP-sensitive potassium (K(ATP)) channels are found in a wide variety of cell types where they couple cell metabolism to electrical activity. In glucose-sensing tissues, these channels respond to fluctuating changes in blood glucose concentration, but in other tissues they are activated only under ischemic conditions or in response to hormonal stimulation. Although K(ATP) channels in different tissues have different regulatory subunits, in almost all cases (except vascular smooth muscle) the pore-forming subunit is the inwardly rectifying K(+) channel Kir6.2. This article reviews recent studies of Kir6.2, focussing on the relation between channel structure and function, and on naturally occurring mutations in Kir6.2 that lead to human disease. New insights into the location of the ATP-binding site, the permeation pathway for K(+), and the gating of the pore provided by homology modelling are discussed in relation to functional studies. Gain-of-function mutations in Kir6.2 cause permanent neonatal diabetes mellitus (PNDM) by reducing the ATP sensitivity of the K(ATP) channel and increasing the K(ATP) current, which is predicted to inhibit beta-cell electrical activity and insulin secretion. Mutations at specific residues, that cause a greater decrease in ATP sensitivity, are associated with additional neurological symptoms. The molecular mechanism underlying the differences in ATP sensitivity produced by these two classes of mutations is discussed. We speculate on how some mutations lead to neurological disease and why no obvious cardiac symptoms are observed. We also consider the implications of these studies for type-2 diabetes.

Distributions of each amino acid in the trans-membrane domain were calculated as a function of the membrane normal using all currently available alpha-helical membrane protein structures with resolutions better than 4 A. The results were compared with previous sequence- and structure-based analyses. Calculation of the average hydrophobicity along the membrane normal demonstrated that the protein surface in the membrane domain is in fact much more hydrophobic than the protein core. While hydrophobic residues dominate the membrane domain, the interfacial regions of membrane proteins were found to be abundant in the small residues glycine, alanine, and serine, consistent with previous studies on membrane protein packing. Charged residues displayed nonsymmetric distributions with a preference for the intracellular interface. This effect was more prominent for Arg and Lys resulting in a direct confirmation of the positive inside rule. Potentials of mean force along the membrane normal were derived for each amino acid by fitting Gaussian functions to the residue distributions. The individual potentials agree well with experimental and theoretical considerations. The resulting implicit membrane potential was tested on various membrane proteins as well as single trans-membrane alpha-helices. All membrane proteins were found to be at an energy minimum when correctly inserted into the membrane. For alpha-helices both interfacial (i.e. surface bound) and inserted configurations were found to correspond to energy minima. The results demonstrate that the use of trans-membrane amino acid distributions to derive an implicit membrane representation yields meaningful residue potentials.

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.

Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10 ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The beta-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with "dimer-of-dimers" motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.

Inwardly rectifying potassium (Kir) channels control cell membrane K+ fluxes and electrical signalling in diverse cell types. Heterozygous mutations in the human Kir6.2 gene (KCNJ11), the pore-forming subunit of the ATP-sensitive (K(ATP)) channel, cause permanent neonatal diabetes mellitus. However, the I296L mutation also results in developmental delay, muscle weakness and epilepsy. We investigated the functional effects of the I296L mutation by expressing wild-type or mutant Kir6.2/SUR1 channels in Xenopus oocytes. The mutation caused a marked increase in resting whole-cell K(ATP) currents by reducing channel inhibition by ATP, in both homomeric and simulated heterozygous states. Kinetic analysis showed that the mutation impaired ATP sensitivity indirectly, by stabilizing the open state of the channel and possibly also by means of an allosteric effect on ATP binding and/or transduction. The results implicate a new region in Kir-channel gating and suggest that disease severity is correlated with the extent of reduction in ATP sensitivity.

The structure of a homopentameric alpha7 nicotinic acetylcholine receptor is modelled by combining structural information from two sources: the X-ray structure of a water soluble acetylcholine binding protein from Lymnea stagnalis, and the electron microscopy derived structure of the transmembrane domain of the Torpedo nicotinic receptor. The alpha7 nicotinic receptor model is generated by simultaneously optimising: (i) chain connectivity, (ii) avoidance of stereochemically unfavourable contacts, and (iii) contact between the beta1-beta2 and M2-M3 loops that have been suggested to be involved in transmission of conformational change between the extracellular and transmembrane domains. A Gaussian network model was used to predict patterns of residue mobility in the alpha7 model. The results of these calculations suggested a flexibility gradient along the transmembrane domain, with the extracellular end of the domain more flexible that the intracellular end. Poisson-Boltzmann (PB) energy calculations and atomistic (molecular dynamics) simulations were used to estimate the free energy profile of a Na+ ion as a function of position along the axis of the pore-lining M2 helix bundle of the transmembrane domain. Both types of calculation suggested a significant energy barrier to exist in the centre of the (closed) pore, consistent with a "hydrophobic gating" model. Estimations of the PB energy profile as a function of ionic strength suggest a role of the extracellular domain in determining the cation selectivity of the alpha7 nicotinic receptor. These studies illustrate how molecular models of members of the nicotinic receptor superfamily of channels may be used to study structure-function relationships.

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

ATP-sensitive potassium (KATP) channels couple cell metabolism to electrical activity by regulating K+ flux across the plasma membrane. Channel closure is mediated by ATP, which binds to the pore-forming subunit (Kir6.2). Here we use homology modelling and ligand docking to construct a model of the Kir6.2 tetramer and identify the ATP-binding site. The model is consistent with a large amount of functional data and was further tested by mutagenesis. Ligand binding occurs at the interface between two subunits. The phosphate tail of ATP interacts with R201 and K185 in the C-terminus of one subunit, and with R50 in the N-terminus of another; the N6 atom of the adenine ring interacts with E179 and R301 in the same subunit. Mutation of residues lining the binding pocket reduced ATP-dependent channel inhibition. The model also suggests that interactions between the C-terminus of one subunit and the 'slide helix' of the adjacent subunit may be involved in ATP-dependent gating. Consistent with a role in gating, mutations in the slide helix bias the intrinsic channel conformation towards the open state.

Background: Nicotinic acetylcholine receptors (nAChRs) are pentameric proteins that are important drug targets for a variety of diseases including Alzheimer's, schizophrenia and various forms of epilepsy. One of the most intensively studied nAChR subunits in recent years has been alpha7. This subunit can form functional homomeric pentamers (alpha7)5, which can make interpretation of physiological and structural data much simpler. The growing amount of structural, pharmacological and physiological data for these receptors indicates the need for a dedicated and accurate database to provide a means to access this information in a coherent manner., Description: A7DB http://www.lgics.org/a7db/ is a new relational database of manually curated experimental physiological data associated with the alpha7 nAChR. It aims to store as much of the pharmacology, physiology and structural data pertaining to the alpha7 nAChR. The data is accessed via web interface that allows a user to search the data in multiple ways: 1) a simple text query 2) an incremental query builder 3) an interactive query builder and 4) a file-based uploadable query. It currently holds more than 460 separately reported experiments on over 85 mutations., Conclusions: A7DB will be a useful tool to molecular biologists and bioinformaticians not only working on the alpha7 receptor family of proteins but also in the more general context of nicotinic receptor modelling. Furthermore it sets a precedent for expansion with the inclusion of all nicotinic receptor families and eventually all cys-loop receptor families.

Octyl glucoside (OG) is a detergent widely employed in structural and functional studies of membrane proteins. To better understand the nature of protein-OG interactions, molecular dynamics simulations (duration 10 ns) have been used to explore an alpha-helical membrane protein, GlpF, in OG micelles and in DMPC bilayers. Greater conformational drift of the extramembraneous protein loops, from the initial X-ray structure, is seen for the GlpF-OG simulations than for the GlpF-DMPC simulation. The mobility of the transmembrane alpha-helices is approximately 1.3x higher in the GlpF-OG than the GlpF-DMPC simulations. The detergent is seen to form an irregular torus around the protein. The presence of the protein leads to a small perturbation in the behavior of the alkyl chains in the OG micelle, namely an approximately 15% increase in the trans-gauche(-)-gauche(+) transition time. Aromatic side chains (Trp, Tyr) and basic side chains (Arg, Lys) play an important role in both protein-detergent (OG) and protein-lipid (DMPC) interactions.