OmpT: molecular dynamics simulations of an outer membrane enzyme

Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT

The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.

A Practical Implicit Membrane Potential for NMR Structure Calculations of Membrane Proteins

The highly anisotropic environment of the lipid bilayer membrane imposes significant constraints on the structures and functions of membrane proteins. However, NMR structure calculations typically use a simple repulsive potential that neglects the effects of solvation and electrostatics, because explicit atomic representation of the solvent and lipid molecules is computationally expensive and impractical for routine NMR-restrained calculations that start from completely extended polypeptide templates. Here, we describe the extension of a previously described implicit solvation potential, eefxPot, to include a membrane model for NMR-restrained calculations of membrane protein structures in XPLOR-NIH. The key components of eefxPot are an energy term for solvation free energy that works together with other nonbonded energy functions, a dedicated force field for conformational and nonbonded protein interaction parameters, and a membrane function that modulates the solvation free energy and dielectric screening as a function of the atomic distance from the membrane center, relative to the membrane thickness. Initial results obtained for membrane proteins with structures determined experimentally in lipid bilayer membranes show that eefxPot affords significant improvements in structural quality, accuracy, and precision. Calculations with eefxPot are straightforward to implement and can be used to both fold and refine structures, as well as to run unrestrained molecular-dynamics simulations. The potential is entirely compatible with the full range of experimental restraints measured by various techniques. Overall, it provides a useful and practical way to calculate membrane protein structures in a physically realistic environment.

Atomistic molecular-dynamics simulations enable prediction of the arginine permeation pathway through OccD1/OprD from Pseudomonas aeruginosa

Pseudomonas aeruginosa is a Gram-negative bacterium that does not contain large, nonspecific porins in its outer membrane. Consequently, the outer membrane is highly impermeable to polar solutes and serves as a barrier against the penetration of antimicrobial agents. This is one of the reasons why such bacteria are intrinsically resistant to antibiotics. Polar molecules that permeate across the outer membrane do so through substrate-specific channels proteins. To design antibiotics that target substrate-channel proteins, it is essential to first identify the permeation pathways of their natural substrates. In P. aeruginosa, the largest family of substrate-specific proteins is the OccD (previously reported under the name OprD) family. Here, we employ equilibrium and steered molecular-dynamics simulations to study OccD1/OprD, the archetypical member of the OccD family. We study the permeation of arginine, one of the natural substrates of OccD1, through the protein. The combination of simulation methods allows us to predict the pathway taken by the amino acid, which is enabled by conformational rearrangements of the extracellular loops of the protein. Furthermore, we show that arginine adopts a specific orientation to form the molecular interactions that facilitate its passage through part of the protein. We predict a three-stage permeation process for arginine.

A conformational landscape for alginate secretion across the outer membrane of Pseudomonas aeruginosa

The exopolysaccharide alginate is an important component of biofilms produced by Pseudomonas aeruginosa, a major pathogen that contributes to the demise of cystic fibrosis patients. Alginate exits the cell via the outer membrane porin AlgE. X-ray structures of several AlgE crystal forms are reported here. Whilst all share a common β-barrel constitution, they differ in the degree to which loops L2 and T8 are ordered. L2 and T8 have been identified as an extracellular gate (E-gate) and a periplasmic gate (P-gate), respectively, that reside on either side of an alginate-selectivity pore located midway through AlgE. Passage of alginate across the membrane is proposed to be regulated by the sequential opening and closing of the two gates. In one crystal form, the selectivity pore contains a bound citrate. Because citrate mimics the uronate monomers of alginate, its location is taken to highlight a route through AlgE taken by alginate as it crosses the pore. Docking and molecular-dynamics simulations support and extend the proposed transport mechanism. Specifically, the P-gate and E-gate are flexible and move between open and closed states. Citrate can leave the selectivity pore bidirectionally. Alginate docks stably in a linear conformation through the open pore. To translate across the pore, a force is required that presumably is provided by the alginate-synthesis machinery. Accessing the open pore is facilitated by complex formation between AlgE and the periplasmic protein AlgK. Alginate can thread through a continuous pore in the complex, suggesting that AlgK pre-orients newly synthesized exopolysaccharide for delivery to AlgE.

Human single-chain urokinase is activated by the omptins PgtE of Salmonella enterica and Pla of Yersinia pestis despite mutations of active site residues

Fibrinolysis is important in cell migration and tightly regulated by specific inhibitors and activators; of the latter, urokinase (uPA) associates with enhancement of cell migration. Active uPA is formed through cleavage of the single-chain uPA (scuPA). The Salmonella enterica strain 14028R cleaved human scuPA at the peptide bond Lys158-Ile159, the site cleaved also by the physiological activator human plasmin. The cleavage led to activation of scuPA, while no cleavage or activation were detected with the mutant strain 14028R lacking the omptin protease PgtE. Complementation and expression studies confirmed the role of PgtE in scuPA activation. Similar cleavage and activation of scuPA were detected with recombinant Escherichia coli expressing the omptin genes pla from Yersinia pestis, ompT and ompP from E. coli, sopA from Shigella flexneri, and leo from Legionella pneumophila. For these omptins the activation of scuPA is the only shared function so far detected. Only poor cleavage and activation of scuPA were seen with YcoA of Y. pestis and YcoB of Yersinia pseudotuberculosis that are considered to be proteolytically inactive omptin variants. Point mutations of active site residues in Pla and PgtE had different effects on the proteolysis of plasminogen and of scuPA, indicating versatility in omptin proteolysis.

Reduced lateral mobility of lipids and proteins in crowded membranes

Coarse-grained molecular dynamics simulations of the E. coli outer membrane proteins FhuA, LamB, NanC, OmpA and OmpF in a POPE/POPG (3∶1) bilayer were performed to characterise the diffusive nature of each component of the membrane. At small observation times (<10 ns) particle vibrations dominate phospholipid diffusion elevating the calculated values from the longer time-scale bulk value (>50 ns) of 8.5×10⁻⁷ cm² s⁻¹. The phospholipid diffusion around each protein was found to vary based on distance from protein. An asymmetry in the diffusion of annular lipids in the inner and outer leaflets was observed and correlated with an asymmetry in charged residues in the vicinity of the inner and outer leaflet head-groups. Protein rotational and translational diffusion were also found to vary with observation time and were inversely correlated with the radius of gyration of the protein in the plane of the bilayer. As the concentration of protein within the bilayer was increased, the overall mobility of the membrane decreased reflected in reduced lipid diffusion coefficients for both lipid and protein components. The increase in protein concentration also resulted in a decrease in the anomalous diffusion exponent α of the lipid. Formation of extended clusters and networks of proteins led to compartmentalisation of lipids in extreme cases.

Biological membranes are selective barriers which control the entry/exit of molecules to/from the interior of a cell. They are composed of a lipid bilayer in which are embedded many membrane proteins. Whilst the individual components of membranes are relatively well characterised, the lateral organization and dynamics of the membrane remain less well understood. The lateral mobility of constituent membrane species affects many processes, including how quickly proteins complexes form and protein recruitment occurs, how quickly lipids can be modified/lysed, and the formation of disordered and ordered microdomains. Biological membranes can contain as much as 50% protein. The dynamics of these crowded environments differ greatly from the sparsely populated membranes often studied in silico or in vitro. We use molecular dynamics computer simulations to quantify how mobility within the membrane decreases as the protein concentration increases. We calculate a baseline diffusion of both lipids and selected bacterial outer membrane proteins in the simplest of systems, namely a single protein in a large lipid bilayer patch. In this case diffusion can be correlated with the size of the protein. We observe how proteins affect the mobility of adjacent lipids. As the protein concentration within our systems increases we show that diffusion of both the proteins and lipids is reduced.

Yersinia infection tools-characterization of structure and function of adhesins

Among the seventeen species of the Gram-negative genus Yersinia, three have been shown to be virulent and pathogenic to humans and animals-Y. enterocolitica, Y. pseudotuberculosis, and Y. pestis. In order to be so, they are armoured with various factors that help them adhere to tissues and organelles, cross the cellular barrier and escape the immune system during host invasion. The group of proteins that mediate pathogen-host interactions constitute adhesins. Invasin, Ail, YadA, YadB, YadC, Pla, and pH 6 antigen belong to the most prominent and best-known Yersinia adhesins. They act at different times and stages of infection complementing each other by their ability to bind a variety of host molecules such as collagen, fibronectin, laminin, β1 integrins, and complement regulators. All the proteins are anchored in the bacterial outer membrane (OM), often forming rod-like or fimbrial-like structures that protrude to the extracellular milieu. Structural studies have shown that the anchor region forms a β-barrel composed of 8, 10, or 12 antiparallel β-strands. Depending on the protein, the extracellular part can be composed of several domains belonging to the immunoglobulin fold superfamily, or form a coiled-coil structure with globular head domain at the end, or just constitute several loops connecting individual β-strands in the β-barrel. Those extracellular regions define the activity of each adhesin. This review focuses on the structure and function of these important molecules, and their role in pathogenesis.

Structural basis for activation of an integral membrane protease by lipopolysaccharide

Omptins constitute a unique family of outer membrane proteases that are widespread in Enterobacteriaceae. The plasminogen activator (Pla) of Yersinia pestis is an omptin family member that is very important for development of both bubonic and pneumonic plague. The physiological function of Pla is to cleave (activate) human plasminogen to form the plasma protease plasmin. Uniquely, lipopolysaccharide (LPS) is essential for the catalytic activity of all omptins, including Pla. Why omptins require LPS for enzymatic activity is unknown. Here, we report the co-crystal structure of LPS-free Pla in complex with the activation loop peptide of human plasminogen, its natural substrate. The structure shows that in the absence of LPS, the peptide substrate binds deep within the active site groove and displaces the nucleophilic water molecule, providing an explanation for the dependence of omptins on LPS for enzymatic activity.

pH-dependent channel gating in connexin26 hemichannels involves conformational changes in N-terminus

Connexin (Cx) hemichannels controlling an exchange of ions and metabolites between the cytoplasm and extracellular milieu can be modulated by the variation of intracellular pH during physiological and pathological conditions. To address the mechanism by which the pH exerts its effect on hemichannels, we have performed two 100-ns molecular dynamics simulations of the Cx26 channel in both acidic and neutral states. The results show that: 1) transmembrane domains undergo clockwise motions around the pore axis under both acidic and neutral conditions, while extracellular segments keep stable. 2) Under neutral condition, Cx26 has a tightly closed configuration that occurs through the assembly of N-terminal helix (NTH) region. This shows a constriction formed by the interhelical interactions of Asp2 and Met1 from neighboring NTH, which shapes the narrowest segment (pore radius<2Å) of the pore, preventing the passage of ions from the extracellular side. This indicates that Asp2 may act as a channel gate. 3) Under the acidic condition, the constriction is relieved by the protonation of Asp2 causing interruption of interhelical interactions, Cx26 has a flexibly opening pore (pore radius>4.5Å) around NTH region, allowing the passage of chloride ions unimpeded by the side-chain Asp2. While in the extracellular part two chloride ions interact with the side-chain Lys41 from three subunits. Finally, we provide a plausible mechanism of pH-dependent gating of hemichannel that involves protonation of the aspartic residues, suggesting that the pH sensitivity of hemichannel permeability is a sophisticated mechanism for cell regulating ion permeation.

GPU-accelerated atom and dynamic bond visualization using hyperballs: a unified algorithm for balls, sticks, and hyperboloids

Ray casting on graphics processing units (GPUs) opens new possibilities for molecular visualization. We describe the implementation and calculation of diverse molecular representations such as licorice, ball-and-stick, space-filling van der Waals spheres, and approximated solvent-accessible surfaces using GPUs. We introduce HyperBalls, an improved ball-and-stick representation replacing tubes, linking the atom spheres by hyperboloids that can smoothly connect them. This type of depiction is particularly useful to represent dynamic phenomena, such as the evolution of noncovalent bonds. It is furthermore well suited to represent coarse-grained models and spring networks. All these representations can be defined by a single general algebraic equation that is adapted for the ray-casting technique and is well suited for execution on the GPU. Using GPU capabilities, this implementation can routinely, accurately, and interactively render molecules ranging from a few atoms up to huge macromolecular assemblies with more than 500,000 particles. In simple cases, based only on spheres, we have been able to display up to two million atoms smoothly.

An active site water network in the plasminogen activator pla from Yersinia pestis

The plasminogen activator Pla from Yersinia pestis is an outer membrane protease (omptin) that is important for the virulence of plague. Here, we present the high-resolution crystal structure of wild-type, enzymatically active Pla at 1.9 A. The structure shows a water molecule located between active site residues D84 and H208, which likely corresponds to the nucleophilic water. A number of other water molecules are present in the active site, linking residues important for enzymatic activity. The R211 sidechain in loop L4 is close to the nucleophilic water and possibly involved in the stabilization of the oxyanion intermediate. Subtle conformational changes of H208 result from the binding of lipopolysaccharide to the outside of the barrel, explaining the unusual dependence of omptins on lipopolysaccharide for activity. The Pla structure suggests a model for the interaction with plasminogen substrate and provides a more detailed understanding of the catalytic mechanism of omptin proteases.

Orientation of the Escherichia coli outer membrane protein OmpX in phospholipid bilayer membranes determined by solid-State NMR

The solid-state NMR orientation-dependent frequencies measured for membrane proteins in macroscopically oriented lipid bilayers provide precise orientation restraints for structure determination in membranes. Here we show that this information can also be used to supplement crystallographic structural data to establish the orientation of a membrane protein in the membrane. This is achieved by incorporating a few orientation restraints, measured for the Escherichia coli outer membrane protein OmpX in magnetically oriented lipid bilayers (bicelles), in a simulated annealing calculation with the coordinates of the OmpX crystal structure. The (1)H-(15)N dipolar couplings measured for the seven Phe residues of OmpX in oriented bilayers can be assigned by back-calculation of the NMR spectrum from the crystal structure and are sufficient to establish the three-dimensional orientation of the protein in the membrane, while the (15)N chemical shifts provide a measure of cross-validation for the analysis. In C14 lipid bilayers, OmpX adopts a transmembrane orientation with a 7 degrees tilt of its beta-barrel axis relative to the membrane normal, matching the hydrophobic thickness of the barrel with that of the membrane.

Complex molecular assemblies at hand via interactive simulations

Studying complex molecular assemblies interactively is becoming an increasingly appealing approach to molecular modeling. Here we focus on interactive molecular dynamics (IMD) as a textbook example for interactive simulation methods. Such simulations can be useful in exploring and generating hypotheses about the structural and mechanical aspects of biomolecular interactions. For the first time, we carry out low-resolution coarse-grain IMD simulations. Such simplified modeling methods currently appear to be more suitable for interactive experiments and represent a well-balanced compromise between an important gain in computational speed versus a moderate loss in modeling accuracy compared to higher resolution all-atom simulations. This is particularly useful for initial exploration and hypothesis development for rare molecular interaction events. We evaluate which applications are currently feasible using molecular assemblies from 1900 to over 300,000 particles. Three biochemical systems are discussed: the guanylate kinase (GK) enzyme, the outer membrane protease T and the soluble N-ethylmaleimide-sensitive factor attachment protein receptors complex involved in membrane fusion. We induce large conformational changes, carry out interactive docking experiments, probe lipid-protein interactions and are able to sense the mechanical properties of a molecular model. Furthermore, such interactive simulations facilitate exploration of modeling parameters for method improvement. For the purpose of these simulations, we have developed a freely available software library called MDDriver. It uses the IMD protocol from NAMD and facilitates the implementation and application of interactive simulations. With MDDriver it becomes very easy to render any particle-based molecular simulation engine interactive. Here we use its implementation in the Gromacs software as an example.

Molecular dynamics simulations of a bacterial autotransporter: NalP fromNeisseria meningitidis

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.

Structural biology of membrane-intrinsic beta-barrel enzymes: sentinels of the bacterial outer membrane

The outer membranes of Gram-negative bacteria are replete with integral membrane proteins that exhibit antiparallel beta-barrel structures, but very few of these proteins function as enzymes. In Escherichia coli, only three beta-barrel enzymes are known to exist in the outer membrane; these are the phospholipase OMPLA, the protease OmpT, and the phospholipidColon, two colonslipid A palmitoyltransferase PagP, all of which have been characterized at the structural level. Structural details have also emerged for the outer membrane beta-barrel enzyme PagL, a lipid A 3-O-deacylase from Pseudomonas aeruginosa. Lipid A can be further modified in the outer membrane by two beta-barrel enzymes of unknown structure; namely, the Salmonella enterica 3'-acyloxyacyl hydrolase LpxR, and the Rhizobium leguminosarum oxidase LpxQ, which employs O(2) to convert the proximal glucosamine unit of lipid A into 2-aminogluconate. Structural biology now indicates how beta-barrel enzymes can function as sentinels that remain dormant when the outer membrane permeability barrier is intact. Host immune defenses and antibiotics that perturb this barrier can directly trigger beta-barrel enzymes in the outer membrane. The ensuing adaptive responses occur instantaneously and rapidly outpace other signal transduction mechanisms that similarly function to restore the outer membrane permeability barrier.

pH-dependent interdomain tethers of CD1b regulate its antigen capture

As CD1 proteins recycle between the cell surface and endosomes, they show altered receptiveness to lipid antigen loading. We hypothesized that changes in proton concentration encountered within distinct endosomal compartments influence the charge state of residues near the entrance to the CD1 groove and thereby control antigen loading. Molecular dynamic models identified flexible areas of the CD1b heavy chain in the superior and lateral walls of the A' pocket. In these same areas, residues that carry charge in a pH-dependent manner (D60, E62) were found to tether the rigid alpha1 helix to flexible areas of the alpha2 helix and the 50-60 loop. After disruption of these tethers with acid pH or mutation, we observed increased association and dissociation of lipids with CD1b and preferential presentation of antigens with bulky lipid tails. We propose that ionic tethers act as molecular switches that respond to pH fluxes during endosomal recycling and regulate the conformation of the CD1 heavy chain to control the size and rate of antigens captured.

Microseconds dynamics simulations of the outer-membrane protease T

Conformational fluctuations of enzymes may play an important role for substrate recognition and/or catalysis, as it has been suggested in the case of the protease enzymatic superfamily. Unfortunately, theoretically addressing this issue is a problem of formidable complexity, as the number of the involved degrees of freedom is enormous: indeed, the biological function of a protein depends, in principle, on all its atoms and on the surrounding water molecules. Here we investigated a membrane protease enzyme, the OmpT from Escherichia coli, by a hybrid molecular mechanics/coarse-grained approach, in which the active site is treated with the GROMOS force field, whereas the protein scaffold is described with a Go-model. The method has been previously tested against results obtained with all-atom simulations. Our results show that the large-scale motions and fluctuations of the electric field in the microsecond timescale may impact on the biological function and suggest that OmpT employs the same catalytic strategy as aspartic proteases. Such a conclusion cannot be drawn within the 10- to 100-ns timescale typical of current molecular dynamics simulations. In addition, our studies provide a structural explanation for the drop in the catalytic activity of two known mutants (S99A and H212A), suggesting that the coarse-grained approach is a fast and reliable tool for providing structure/function relationships for both wild-type OmpT and mutants.

NMR structural studies of the bacterial outer membrane protein OmpX in oriented lipid bilayer membranes

The beta-barrels found in the outer membranes of prokaryotic and eukaryotic organisms constitute an important functional class of proteins. Here we present solid-state NMR spectra of the bacterial outer membrane protein OmpX in oriented lipid bilayer membranes. We show that OmpX is folded in both glass-supported oriented lipid bilayers and in lipid bicelles that can be magnetically oriented with the membrane plane parallel or perpendicular to the direction of the magnetic field. The presence of resolved peaks in these spectra demonstrates that OmpX undergoes rotational diffusion around an axis perpendicular to the membrane surface. A tightly hydrogen-bonded domain of OmpX resists exchange with D2O for days and is assigned to the transmembrane beta-barrel, while peaks at isotropic resonance frequencies that disappear rapidly in D2O are assigned to the extracellular and periplasmic loops. The two-dimensional 1H/15N separated local field spectra of OmpX have several resolved peaks, and agree well with the spectra calculated from the crystal structure of OmpX rotated with the barrel axis nearly parallel (5 degrees tilt) to the direction of the magnetic field. The data indicate that it will be possible to obtain site-specific resonance assignments and to determine the structure, tilt, and rotation of OmpX in membranes using the solid-state NMR methods that are currently being applied to alpha-helical membrane 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 x 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 beta-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 alpha-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 approximately 10 ns timescale.

Three hydrolases and a transferase: comparative analysis of active-site dynamics via the BioSimGrid database

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.

Positioning of proteins in membranes: a computational approach

A new computational approach has been developed to determine the spatial arrangement of proteins in membranes by minimizing their transfer energies from water to the lipid bilayer. The membrane hydrocarbon core was approximated as a planar slab of adjustable thickness with decadiene-like interior and interfacial polarity profiles derived from published EPR studies. Applicability and accuracy of the method was verified for a set of 24 transmembrane proteins whose orientations in membranes have been studied by spin-labeling, chemical modification, fluorescence, ATR FTIR, NMR, cryo-microscopy, and neutron diffraction. Subsequently, the optimal rotational and translational positions were calculated for 109 transmembrane, five integral monotopic and 27 peripheral protein complexes with known 3D structures. This method can reliably distinguish transmembrane and integral monotopic proteins from water-soluble proteins based on their transfer energies and membrane penetration depths. The accuracies of calculated hydrophobic thicknesses and tilt angles were approximately 1 A and 2 degrees, respectively, judging from their deviations in different crystal forms of the same proteins. The hydrophobic thicknesses of transmembrane proteins ranged from 21.1 to 43.8 A depending on the type of biological membrane, while their tilt angles with respect to the bilayer normal varied from zero in symmetric complexes to 26 degrees in asymmetric structures. Calculated hydrophobic boundaries of proteins are located approximately 5 A lower than lipid phosphates and correspond to the zero membrane depth parameter of spin-labeled residues. Coordinates of all studied proteins with their membrane boundaries can be found in the Orientations of Proteins in Membranes (OPM) database:http://opm.phar.umich.edu/.

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.

Molecular dynamics simulations of proteins in lipid bilayers

With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.

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.

Molecular simulations and lipid–protein interactions: potassium channels and other membrane proteins

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 α-helical bundle) and OmpA (a β-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.

Orientation of β-Barrel Proteins OmpA and FhuA in Lipid Membranes. Chain Length Dependence from Infrared Dichroism†

The outer-membrane proteins OmpA and FhuA of Escherichia coli are monomeric beta-barrels of widely differing size. Polarized attenuated total reflection infrared spectroscopy has been used to determine the orientation of the beta-barrels in phosphatidylcholine host matrices of different lipid chain lengths. The linear dichroism of the amide I band from OmpA and FhuA in hydrated membranes generally increases with increasing chain length from diC(12:0) to diC(17:0) phosphatidylcholine, in both the fluid and gel phases. Measurements of the amide I and amide II dichroism from dry samples are used to deduce the strand tilt (beta = 46 degrees for OmpA and beta = 44.5 degrees for FhuA). These values are then used to deduce the order parameters, P(2)(cos alpha), of the beta-barrels from the amide I dichroic ratios of the hydrated membranes. The orientational ordering of the beta-barrels and their assembly in the membrane are discussed in terms of hydrophobic matching with the lipid chains.

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.