A multidomain outer membrane protein from Pasteurella multocida: Modelling and simulation studies of PmOmpA

PmOmpA is a two-domain outer membrane protein from Pasteurella multocida. The N-terminal domain of PmOmpA is a homologue of the transmembrane β-barrel domain of OmpA from Escherichia coli, whilst the C-terminal domain of PmOmpA is a homologue of the extra-membrane Neisseria meningitidis RmpM C-terminal domain. This enables a model of a complete two domain PmOmpA to be constructed and its conformational dynamics explored via MD simulations of the protein embedded within two different phospholipid bilayers (DMPC and DMPE). The conformational stability of the transmembrane β-barrel is similar to that of a homology model of OprF from Pseudomonas aeruginosa in bilayer simulations. There is a degree of water penetration into the interior of the β-barrel, suggestive of a possible transmembrane pore. Although the PmOmpA model is stable over 20 ns simulations, retaining its secondary structure and fold integrity throughout, substantial flexibility is observed in a short linker region between the N- and the C-terminal domains. At low ionic strength, the C-terminal domain moves to interact electrostatically with the lipid bilayer headgroups. This study demonstrates that computational approaches may be applied to more complex, multi-domain outer membrane proteins, rather than just to transmembrane β-barrels, opening the possibility of in silico proteomics approaches to such proteins.     

Structure refinement of membrane proteins via molecular dynamics simulations

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

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.     

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.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Incompatibility of outer membrane proteins OmpA and OmpF of Escherichia coli with secretion in Bacillus subtilis: Fusions with secretable peptides

The secretion of the outer membrane proteins OmpA and OmpF of Escherichia coli has previously been found to be blocked at an early intracellular step, when these proteins were fused to a bacillar signal sequence and expressed in Bacillus subtilis. We have now fused these proteins to long secretable polypeptides, the amino-terminal portions of alpha-amylase or beta-lactamase. In spite of this, no secretion of the fusion proteins was detected in B. subtilis. With the exception of a small fraction of the beta-lactamase fusion, the proteins were cell-bound with uncleaved signal sequences. Protease accessibility indicated that the fusion proteins were not even partially exposed on the outer surface of the cytoplasmic membrane. Thus there was no change of the location compared to the OmpA or OmpF fused to the signal sequence only. We conclude that, like OmpA and OmpF, the fusion proteins fold into an export-incompatible conformation in B. subtilis before the start of translocation, which we postulate to be a late post-translational event.     

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.     

Outer membrane proteins: comparing X-ray and NMR structures by MD simulations in lipid bilayers

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 × 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the ‘quality’ of the starting structure seems mainly to influence the simulation stability of the transmembrane β-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal α-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a ∼10 ns timescale.     

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.     

Driven to near-experimental accuracy by refinement via molecular dynamics simulations

Protein model refinement has been an essential part of successful protein structure prediction. Molecular dynamics simulation-based refinement methods have shown consistent improvement of protein models. There had been progress in the extent of refinement for a few years since the idea of ensemble averaging of sampled conformations emerged. There was little progress in CASP12 because conformational sampling was not sufficiently diverse due to harmonic restraints. During CASP13, a new refinement method was tested that achieved significant improvements over CASP12. The new method intended to address previous bottlenecks in the refinement problem by introducing new features. Flat-bottom harmonic restraints replaced harmonic restraints, sampling was performed iteratively, and a new scoring function and selection criteria were used. The new protocol expanded conformational sampling at reduced computational costs. In addition to overall improvements, some models were refined significantly to near-experimental accuracy.     

The mechanism of proton translocation in respiratory complex I from molecular dynamics

Abstract

Respiratory complex I, the biggest enzyme of respiratory chain, plays a key role in energy production by the mitochondrial respiratory chain and has been implicated in many human neurodegenerative diseases. Recently, the crystal structure of respiratory complex I is reported. We perform 50?ns molecular dynamics simulations on the membrane domain of respiratory complex I under two hypothetical states (oxidized state and reduced state). We find that the density of water molecules in the trans-membrane domain under reduced state is bigger than that under oxidized state. The connecting elements (helix HL and β-hairpins-helix element) fluctuate stronger under reduced state than that under oxidized state, causing more internal water molecules and facilitating the proton conduction. The conformational changes of helix HL and the crucial charged residue Glu in TM5 play key roles in the mechanism of proton translocation. Our results illustrate the dynamic behavior and the potential mechanism of respiratory complex I, which provides the structural basis for drug design of respiratory complex I.     

The transmembrane domain of Neu in a lipid bilayer: molecular dynamics simulations

The results of full-atom molecular dynamics simulations of the transmembrane domains (TMDs) of both native, and Glu⁶⁶⁴-mutant (either protonated or unprotonated) Neu in an explicit fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer are presented. For the native TMD peptide, a 10.05 ns trajectory was collected, while for the mutant TMD peptides 5.05 ns trajectories were collected for each. The peptides in all three simulations display stable predominantly -helical hydrogen bonding throughout the trajectories. The only significant exception occurs near the C-terminal end of the native and unprotonated mutant TMDs just outside the level of the lipid headgroups, where -helical hydrogen bonding develops, introducing a kink in the backbone structure. However, there is no indication of the formation of a bulge within the hydrophobic region of either native or mutant peptides. Over the course of the simulation of the mutant peptide, it is found that a significant number of water molecules penetrate the hydrophobic region of the surrounding lipid molecules, effectively hydrating Glu⁶⁶⁴. If the energy cost of such water penetration is significant enough, this may be a factor in the enhanced dimerization affinity of Glu⁶⁶⁴-mutant Neu.     

A noncanonical binding site of chloramphenicol revealed via molecular dynamics simulations

Chloramphenicol, an antibiotic belonging to the family of amphenicols, is an inhibitor of translation. On the basis of X–ray structural analysis of the binding of chloramphenicol to free bacterial ribosomes, the chloramphenicol action mechanism that consists in preventing the binding of aminoacyl-tRNA to the A–site of the large subunit of the ribosome was adopted. However, the known structures of chloramphenicol complexes with bacterial ribosomes poorly explain the results of the experiments on the chemical modification of 23S rRNA, the resistance to chloramphenicol caused by mutations in 23S rRNA and, which is particularly important, the selectivity of chloramphenicol in suppression of translation, depending on the amino acid sequence of the nascent peptide. In the present study the putative structure of the chloramphenicol complex with a bacterial ribosome in the A,A/P,P–state has been obtained by molecular dynamics simulations methods. The proposed structure of the complex allows us to explain the results of biochemical studies of the interaction of chloramphenicol with the bacterial ribosome.     

Membrane topology and assembly of the outer membrane protein OmpA of Escherichia coli K12

The 325-residue outer membrane protein OmpA of Escherichia coli has been proposed to consist of a membrane-embedded moiety (residues 1 to about 170) and a C-terminal periplasmic region. The former is thought to comprise eight transmembrane segments in the form of antiparallel -strands, forming an amphiphilic connected by exposed turns. Several questions concerning this model were addressed. Thus no experimental evidence had been presented for the turns at the inner leaflet of the membrane and it was not known whether or not the periplasmic part of the polypeptide plays a role in the process of membrane incorporation. Oligonucleotides encoding trypsin cleavage sites were inserted at the predicted turn sites of the ompA gene and it was shown that the encoded proteins indeed become accessible to trypsin at the modified sites. Together with previous results, these data also show that the turns on both sides of the membrane do not possess specifically topogenic information. In two cases one of the two expected tryptic fragments was lost and could be detected at low concentration in only one case. Therefore, bilateral proteolytic digestion of outer membranes can cause loss of -strands and does not necessarily produce a reliable picture of protein topology. When ompA genes were constructed coding for proteins ending at residue 228 or 274, the membrane assembly of these proteins was shown to be partially defective with about 20% of the proteins not being assembled. No such defect was observed when, following the introduction of a premature stop codon, a truncated protein was produced ending with residue 171. It is concluded that (1) the proposed -barrel structure is essentially correct and (2) the periplasmic part of OmpA does not play an active role in, but can, when present in mutant form, interfere with membrane assembly.     

CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

Sorting signal of Escherichia coli OmpA is modified by oligo-(R)-3-hydroxybutyrate

Escherichia coli outer membrane protein A (OmpA) is a well-established model for the study of membrane assembly. Previous studies have shown that the essential sequence for outer membrane localization, known as the sorting signal, is contained in a segment of the eighth β-strand, residues 163-171. Sequential digestion of OmpA, purified from outer membranes or inclusion bodies with cyanogen bromide and Staphylococcus aureus GluC, yielded peptides 162-174(LSLGVSYRFGQGE). Western blot and chemical assays indicated that the peptide was covalently modified by oligo-(R)-3-hydroxybutyrate (cOHB), a flexible, amphipathic oligoester. MALDI/MS was consistent with modification of peptides 162-174 by up to ten R-3-hydroxybutyrate (HB) residues. Western blot analysis of mutants of the peptide, using anti-OHB IgG, indicated that cOHB modification was not inhibited by the single mutations S163G, S167G, Y168F, R169N or R169D; however, cOHB was not detected on peptides containing the double mutations S163G:S167G S163G:V166G, L162G:S167G, and L164G:S167G. MALDI/MS/MS of double mutant S163G:S167G confirmed the absence of cOHB-modification. The results suggest that cOHB may be attached to one or both serines, and point to the importance of the flanking hydrophobic residues. Modification by cOHB may play a role in outer membrane targeting and assembly of OmpA.     

Setting up and running molecular dynamics simulations of membrane proteins

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

A prokaryotic glutamate receptor: homology modelling and molecular dynamics simulations of GluR0

GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.     

Active site modeling in copper azurin molecular dynamics simulations

Active site modeling in molecular dynamics simulations is investigated for the reduced state of copper azurin. Five simulation runs (5 ns each) were performed at room temperature to study the consequences of a mixed electrostatic/constrained modeling for the coordination between the metal and the polypeptide chain, using for the ligand residues a set of charges that is modified with respect to the apo form of the protein by the presence of the copper ion.The results show that the different charge values do not lead to relevant effects on the geometry of the active site of the protein, as long as bond distance constraints are used for all the five ligand atoms. The distance constraint on the O atom of Gly45 can be removed without altering the active site geometry. The coordination between Cu and the other axial ligand Met121 is outlined as being flexible. Differences are found between the bonds of the copper ion with the two apparently equivalent N¹ atoms of His46 and His117.The overall findings are discussed in connection with the issue of determining a model for the active site of azurin suitable to be used in molecular dynamics simulations under unfolding conditions. Figure Model of azurin active site. Copper ligand residues are cut off at C position except Gly45, for which the portion of backbone connecting it to His46 is shown. Only polar H atoms are shown. All atoms are in standard colors (Cu in violet), and the five ligands are labeled     

Omp52 is a growth-phase-regulated outer membrane protein of Leptospira santarosai serovar Shermani

We report the expression and characterization of the omp52 gene of Leptospira santarosai serovar Shermani strain CCF that is isolated in Taiwan. omp52 was identified among pathogenic leptospires but not among non-pathogenic leptospires by using suppression subtractive hybridization in our previous study. With an open reading frame of 1371 bp that encodes 456 amino acids and a predicted molecular mass of 52.6 kDa, Omp52 was shown to be an outer membrane protein containing a C-terminal OmpA consensus domain and exposed on the cell surface. Furthermore, Omp52 increases dramatically during the stationary phase, indicating that the expression of Omp52 is environmentally regulated. By using immunoblotting analysis, we proved that Omp52 was expressed in human patients infected with leptospires. These observations suggest that Omp52 may play roles in the interaction of host cells and pathogens during infection.