PubMed Health⌕ Search

Biomedical subjects

Peter J Bond

Publications and source records attributed to Peter J Bond.

At least 19 recordsLinked to original sources
PubMed2006

Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer.

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.

Cell Membrane↗
PubMed2006

Membrane simulations of OpcA: gating in the loops?

Mobility of extracellular loops may play an important role in the function of outer membrane proteins from Gram-negative bacteria. Molecular dynamics simulations of OpcA from Neisseria meningitidis, embedded in a lipid bilayer, have been used to explore the relationship between the crystal structure and dynamic function of this protein. The results suggest that the crystal environment may constrain the membrane protein structure in a nonphysiological state. The presence of lipids and physiological salt concentrations result in changes in the conformation of the extracellular loops of OpcA, leading to opening of a pore, and to modulation of the molecular surface implicated in recognition of proteoglycan. These changes may be related to the role of OpcA in pathogenesis via modulation of the conformation of a possible sialic acid binding site.

Bacterial Outer Membrane Proteins↗
PubMed2006

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.

Amino Acid Sequence↗
PubMed2006

Cysteine-scanning mutagenesis and disulfide mapping studies of the conserved domain of the twin-arginine translocase TatB component.

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.

Amino Acid Sequence↗
PubMed2006

MD simulations of Mistic: conformational stability in detergent micelles and water.

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.

Bacillus subtilis↗
PubMed2006

Membrane protein dynamics and detergent interactions within a crystal: a simulation study of OmpA.

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.

Bacterial Outer Membrane Proteins↗
PubMed2006

Modeling and simulations of a bacterial outer membrane protein: OprF from Pseudomonas aeruginosa.

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.

Bacteria↗
PubMed2006

Insertion and assembly of membrane proteins via simulation.

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.

Bacterial Outer Membrane Proteins↗
PubMed2005

Anionic phospholipid interactions with the potassium channel KcsA: simulation studies.

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.

Anions↗
PubMed2005

Molecular dynamics simulations of GlpF in a micelle vs in a bilayer: conformational dynamics of a membrane protein as a function of environment.

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.

Aquaporins↗
PubMed2005

Membrane protein structure quality in molecular dynamics simulation.

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.

Computer Simulation↗
PubMed2004

MD simulations of spontaneous membrane protein/detergent micelle formation.

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

Bacterial Outer Membrane Proteins↗
PubMed2004

Conformational sampling and dynamics of membrane proteins from 10-nanosecond computer simulations.

In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small beta-barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C(alpha)-atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extra-membranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations.

Computer Simulation↗
PubMed2004

Lipid-protein interactions of integral membrane proteins: a comparative simulation study.

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.

Bacterial Outer Membrane Proteins↗
PubMed2003

Lipid/protein interactions and the membrane/water interfacial region.

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

Bacterial Outer Membrane Proteins↗
PubMed2003

Membrane protein dynamics versus environment: simulations of OmpA in a micelle and in a bilayer.

The bacterial outer membrane protein OmpA is one of the few membrane proteins whose structure has been solved both by X-ray crystallography and by NMR. Crystals were obtained in the presence of detergent, and the NMR structure is of the protein in a detergent micelle. We have used 10 ns duration molecular dynamics simulations to compare the behaviour of OmpA in a detergent micelle and in a phospholipid bilayer. The dynamic fluctuations of the protein structure seem to be ca 1.5 times greater in the micelle environment than in the lipid bilayer. There are subtle differences between the nature of OmpA-detergent and OmpA-lipid interactions. As a consequence of the enhanced flexibility of the OmpA protein in the micellar environment, side-chain torsion angle changes are such as to lead to formation of a continuous pore through the centre of the OmpA molecule. This may explain the experimentally observed channel formation by OmpA.

Bacterial Outer Membrane Proteins↗
PubMed2002

OmpA: a pore or not a pore? Simulation and modeling studies.

The bacterial outer membrane protein OmpA is composed of an N-terminal 171-residue beta-barrel domain (OmpA(171)) that spans the bilayer and a periplasmic, C-terminal domain of unknown structure. OmpA has been suggested to primarily serve a structural role, as no continuous pore through the center of the barrel can be discerned in the crystal structure of OmpA(171). However, several groups have recorded ionic conductances for bilayer-reconstituted OmpA(171). To resolve this apparent paradox we have used molecular dynamics (MD) simulations on OmpA(171) to explore the conformational dynamics of the protein, in particular the possibility of transient formation of a central pore. A total of 19 ns of MD simulations of OmpA(171) have been run, and the results were analyzed in terms of 1) comparative behavior of OmpA(171) in different bilayer and bilayer-mimetic environments, 2) solvation states of OmpA(171), and 3) pore characteristics in different MD simulations. Significant mobility was observed for residues and water molecules within the beta-barrel. A simulation in which putative gate region side chains of the barrel interior were held in a non-native conformation led to an open pore, with a predicted conductance similar to experimental measurements. The OmpA(171) pore has been shown to be somewhat more dynamic than suggested by the crystal structure. A gating mechanism is proposed to explain its documented channel properties, involving a flickering isomerization of Arg138, forming alternate salt bridges with Glu52 (closed state) and Glu128 (open state).

Arginine↗