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

Molecular-dynamics simulations of macromolecular diffraction, part I: Preparation of protein crystal simulations

Molecular-dynamics (MD) simulations of protein crystals enable the prediction of structural and dynamical features of both the protein and the solvent components of macromolecular crystals, which can be validated against diffraction data from X-ray crystallographic experiments. The simulations have been useful for studying and predicting both Bragg and diffuse scattering in protein crystallography; however, the preparation is not yet automated and includes choices and tradeoffs that can impact the results. Here we examine some of the intricacies and consequences of the choices involved in setting up MD simulations of protein crystals for the study of diffraction data, and provide a recipe for preparing the simulations, packaged in an accompanying Jupyter notebook. This article and the accompanying notebook are intended to serve as practical resources for researchers wishing to put these models to work.

Tris(hydroxymethyl)aminomethane Linker-Bearing Triazine-Based Triglucosides for Solubilization and Stabilization of Membrane Proteins

High-resolution membrane protein structures are essential for a fundamental understanding of the molecular basis of diverse cellular processes and for drug discovery. Detergents are widely used to extract membrane-spanning proteins from membranes and maintain them in a functional state for downstream characterization. Due to limited long-term stability of membrane proteins encapsulated in conventional detergents, development of novel agents is required to facilitate membrane protein structural study. In the current study, we designed and synthesized tris(hydroxymethyl)aminomethane linker-bearing triazine-based triglucosides (TTGs) for solubilization and stabilization of membrane proteins. When these glucoside detergents were evaluated for four membrane proteins including two G protein-coupled receptors, a few TTGs including TTG-C10 and TTG-C11 displayed markedly enhanced behaviors toward membrane protein stability relative to two maltoside detergents [DDM (n-dodecyl-β-d-maltoside) and LMNG (lauryl maltose neopentyl glycol)]. This is a notable feature of the TTGs as glucoside detergents tend to be inferior to maltoside detergents at stabilizing membrane proteins. The favorable behavior of the TTGs for membrane protein stability is likely due to the high hydrophobicity of the lipophilic groups, an optimal range of hydrophilic-lipophilic balance, and the absence of cis-trans isomerism.

The structural integrity of the membrane-embedded bacterial division complex FtsQBL studied with molecular dynamics simulations

The FtsQBL is an essential molecular complex sitting midway through bacterial divisome assembly. To visualize and understand its structure, and the consequences of its membrane anchorage, we produced a model of the E. coli complex using the deep-learning prediction utility, AlphaFold 2. The heterotrimeric model was inserted into a 3-lipid model membrane and subjected to a 500-ns atomistic molecular dynamics simulation. The model is superb in quality and captures most experimentally derived structural features, at both the secondary structure and the side-chain levels. The model consists of a uniquely interlocking module contributed by the C-terminal regions of all three proteins. The functionally important constriction control domain residues of FtsB and FtsL are located at a fixed vertical position of ∼43-49 Å from the membrane surface. While the periplasmic domains of all three proteins are well-defined and rigid, the single transmembrane helices of each are flexible and their collective twisting and bending contribute to most structural variations, according to principal component analysis. Considering FtsQ only, the protein is more flexible in its free state relative to its complexed state-with the biggest structural changes located at the elbow between the transmembrane helix and the α-domain. The disordered N-terminal domains of FtsQ and FtsL associate with the cytoplasmic surface of the inner membrane instead of freely venturing into the solvent. Contact network analysis highlighted the formation of the interlocking trimeric module in FtsQBL as playing a central role in mediating the overall structure of the complex.

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.

Flexible Loop in Carbohydrate-Binding Module 48 Allosterically Modulates Substrate Binding of the 1,4-α-Glucan Branching Enzyme

The 1,4-α-glucan branching enzyme (GBE, EC 2.4.1.18) catalyzes the formation of α-1,6 branching points in starch and plays a key role in synthesis. To obtain mechanistic insights into the catalytic action of the enzyme, we first determined the crystal structure of GBE from Rhodothermus obamensis STB05 (RoGBE) to a resolution of 2.39 Å (PDB ID: 6JOY). The structure consists of three domains: domain A, domain C, and the carbohydrate-binding module 48 (CBM48). An engineered truncated mutant lacking the CBM48 domain (ΔCBM48) showed significantly reduced ligand binding affinity and enzyme activity. Comparison of the structures of RoGBE with other GBEs showed that CBM48 of RoGBE had a longer flexible loop. Truncation of the flexible loops resulted in reduced binding affinity and activity, thereby substantiating the importance of the optimum loop structure for catalysis. In essence, our study shows that CBM48, especially the flexible loop, plays an important role in substrate binding and enzymatic activity of RoGBE. Further, based on the structural analysis, kinetics, and activity assays on wild type and mutants, as well as homology modeling, we proposed a mechanistic model (called the "lid model") to illustrate how the flexible loop triggers substrate binding, ultimately leading to catalysis.

New Malonate-Derived Tetraglucoside Detergents for Membrane Protein Stability

Membrane proteins are widely studied in detergent micelles, a membrane-mimetic system formed by amphiphilic compounds. However, classical detergents have serious limitations in their utility, particularly for unstable proteins such as eukaryotic membrane proteins and membrane protein complexes, and thus, there is an unmet need for novel amphiphiles with enhanced ability to stabilize membrane proteins. Here, we developed a new class of malonate-derived detergents with four glucosides, designated malonate-derived tetra-glucosides (MTGs), and compared these new detergents with previously reported octyl glucose neopentyl glycol (OGNG) and n-dodecyl-β-d-maltoside (DDM). When tested with two G-protein coupled receptors (GPCRs) and three transporters, a couple of MTGs consistently conferred enhanced stability to all tested proteins compared to DDM and OGNG. As a result of favorable behaviors for a range of membrane proteins, these MTGs have substantial potential for membrane protein research. This study additionally provides a new detergent design principle based on the effect of a polar functional group (i.e., ether) on protein stability depending on its position in the detergent scaffold.

Role of Computational Methods in Going beyond X-ray Crystallography to Explore Protein Structure and Dynamics

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.

Prediction of the Closed Conformation and Insights into the Mechanism of the Membrane Enzyme LpxR

Covalent modification of outer membrane lipids of Gram-negative bacteria can impact the ability of the bacterium to develop resistance to antibiotics as well as modulating the immune response of the host. The enzyme LpxR from Salmonella typhimurium is known to deacylate lipopolysaccharide molecules of the outer membrane; however, the mechanism of action is unknown. Here, we employ molecular dynamics and Monte Carlo simulations to study the conformational dynamics and substrate binding of LpxR in representative outer membrane models as well as detergent micelles. We examine the roles of conserved residues and provide an understanding of how LpxR binds its substrate. Our simulations predict that the catalytic H122 must be Nε-protonated for a single water molecule to occupy the space between it and the scissile bond, with a free binding energy of -8.5 kcal mol^-1. Furthermore, simulations of the protein within a micelle enable us to predict the structure of the putative "closed" protein. Our results highlight the need for including dynamics, a representative environment, and the consideration of multiple tautomeric and rotameric states of key residues in mechanistic studies; static structures alone do not tell the full story.

Temperature dependence of protein-water interactions in a gated yeast aquaporin

Regulation of aquaporins is a key process of living organisms to counteract sudden osmotic changes. Aqy1, which is a water transporting aquaporin of the yeast Pichia pastoris, is suggested to be gated by chemo-mechanical stimuli as a protective regulatory-response against rapid freezing. Here, we tested the influence of temperature by determining the X-ray structure of Aqy1 at room temperature (RT) at 1.3 Å resolution, and by exploring the structural dynamics of Aqy1 during freezing through molecular dynamics simulations. At ambient temperature and in a lipid bilayer, Aqy1 adopts a closed conformation that is globally better described by the RT than by the low-temperature (LT) crystal structure. Locally, for the blocking-residue Tyr31 and the water molecules inside the pore, both LT and RT data sets are consistent with the positions observed in the simulations at room-temperature. Moreover, as the temperature was lowered, Tyr31 adopted a conformation that more effectively blocked the channel, and its motion was accompanied by a temperature-driven rearrangement of the water molecules inside the channel. We therefore speculate that temperature drives Aqy1 from a loosely- to a tightly-blocked state. This analysis provides high-resolution structural evidence of the influence of temperature on membrane-transport channels.

Dimerization of the EphA1 receptor tyrosine kinase transmembrane domain: Insights into the mechanism of receptor activation

EphA1 is a receptor tyrosine kinase (RTK) that plays a key role in developmental processes, including guidance of the migration of axons and cells in the nervous system. EphA1, in common with other RTKs, contains an N-terminal extracellular domain, a single transmembrane (TM) α-helix, and a C-terminal intracellular kinase domain. The TM helix forms a dimer, as seen in recent NMR studies. We have modeled the EphA1 TM dimer using a multiscale approach combining coarse-grain (CG) and atomistic molecular dynamics (MD) simulations. The one-dimensional potential of mean force (PMF) for this system, based on interhelix separation, has been calculated using CG MD simulations. This provides a view of the free energy landscape for helix–helix interactions of the TM dimer in a lipid bilayer. The resulting PMF profiles suggest two states, consistent with a rotation-coupled activation mechanism. The more stable state corresponds to a right-handed helix dimer interacting via an N-terminal glycine zipper motif, consistent with a recent NMR structure (2K1K). A second metastable state corresponds to a structure in which the glycine zipper motif is not involved. Analysis of unrestrained CG MD simulations based on representative models from the PMF calculations or on the NMR structure reveals possible pathways of interconversion between these two states, involving helix rotations about their long axes. This suggests that the interaction of TM helices in EphA1 dimers may be intrinsically dynamic. This provides a potential mechanism for signaling whereby extracellular events drive a shift in the repopulation of the underlying TM helix dimer energy landscape.

Dodecyl maltoside protects membrane proteins in vacuo

Molecular dynamics simulations have been used to characterize the effects of transfer from aqueous solution to a vacuum to inform our understanding of mass spectrometry of membrane-protein-detergent complexes. We compared two membrane protein architectures (an α-helical bundle versus a β-barrel) and two different detergent types (phosphocholines versus an alkyl sugar) with respect to protein stability and detergent packing. The β-barrel membrane protein remained stable as a protein-detergent complex in vacuum. Zwitterionic detergents formed conformationally destabilizing interactions with an α-helical membrane protein after detergent micelle inversion driven by dehydration in vacuum. In contrast, a nonionic alkyl sugar detergent resisted micelle inversion, maintaining the solution-phase conformation of the protein. This helps to explain the relative stability of membrane proteins in the presence of alkyl sugar detergents such as dodecyl maltoside.

Packing interface energetics in different crystal forms of the λ Cro dimer

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site-directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson-Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ∼5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures.

Dangerous liaisons between detergents and membrane proteins. The case of mitochondrial uncoupling protein 2

The extraction of membrane proteins from their native environment by detergents is central to their biophysical characterization. Recent studies have emphasized that detergents may perturb the structure locally and modify the dynamics of membrane proteins. However, it remains challenging to determine whether these perturbations are negligible or could be responsible for misfolded conformations, altering the protein's function. In this work, we propose an original strategy combining functional studies and molecular simulations to address the physiological relevance of membrane protein structures obtained in the presence of detergents. We apply our strategy to a structure of isoform 2 of an uncoupling protein (UCP2) binding an inhibitor recently obtained in dodecylphosphocholine detergent micelles. Although this structure shares common traits with the ADP/ATP carrier, a member of the same protein family, its functional and biological significance remains to be addressed. In the present investigation, we demonstrate how dodecylphosphocholine severely alters the structure as well as the function of UCPs. The proposed original strategy opens new vistas for probing the physiological relevance of three-dimensional structures of membrane proteins obtained in non-native environments.

Detergent-mediated protein aggregation

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

The simulation approach to lipid-protein interactions

The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.

Conformational exchange in a membrane transport protein is altered in protein crystals

Successful macromolecular crystallography requires solution conditions that may alter the conformational sampling of a macromolecule. Here, site-directed spin labeling is used to examine a conformational equilibrium within BtuB, the Escherichia coli outer membrane transporter for vitamin B(12). Electron paramagnetic resonance (EPR) spectra from a spin label placed within the N-terminal energy coupling motif (Ton box) of BtuB indicate that this segment is in equilibrium between folded and unfolded forms. In bilayers, substrate binding shifts this equilibrium toward the unfolded form; however, EPR spectra from this same spin-labeled mutant indicate that this unfolding transition is blocked in protein crystals. Moreover, crystal structures of this spin-labeled mutant are consistent with the EPR result. When the free energy difference between substates is estimated from the EPR spectra, the crystal environment is found to alter this energy by 3 kcal/mol when compared to the bilayer state. Approximately half of this energy change is due to solutes or osmolytes in the crystallization buffer, and the remainder is contributed by the crystal lattice. These data provide a quantitative measure of how a conformational equilibrium in BtuB is modified in the crystal environment, and suggest that more-compact, less-hydrated substates will be favored in protein crystals.

Membrane protein dynamics in different environments: simulation study of the outer membrane protein X in a lipid bilayer and in a micelle

The bacterial outer membrane protein OmpX from Escherichia coli has been investigated by molecular dynamics simulations when embedded in a phospholipid bilayer and as a protein-micelle aggregate. The resulting simulation trajectories were analysed in terms of structural and dynamic properties of the membrane protein. In agreement with experimental observations, highest relative stability was found for the β-barrel region that is embedded in the lipophilic phase, whereas an extracellular protruding β-sheet, which is a unique structural feature of OmpX that supposedly plays an important role in cell adhesion and invasion, shows larger structure fluctuations. Additionally, we investigated water permeation into the core of the β-barrel protein, which contains a tight salt-bridge and hydrogen-bond network, so that extensive water flux is unlikely. Differences between the bilayer and the micellar system were observed in the length of the barrel and its position inside the lipid environment, and in the protein interactions with the hydrophilic part of the lipids near the lipid/water interface. Those variations suggest that micelles and other detergent environments might not offer a wholly membrane-like milieu to promote adoption of the physiological conformational state by OmpX.

APBSmem: a graphical interface for electrostatic calculations at the membrane

Electrostatic forces are one of the primary determinants of molecular interactions. They help guide the folding of proteins, increase the binding of one protein to another and facilitate protein-DNA and protein-ligand binding. A popular method for computing the electrostatic properties of biological systems is to numerically solve the Poisson-Boltzmann (PB) equation, and there are several easy-to-use software packages available that solve the PB equation for soluble proteins. Here we present a freely available program, called APBSmem, for carrying out these calculations in the presence of a membrane. The Adaptive Poisson-Boltzmann Solver (APBS) is used as a back-end for solving the PB equation, and a Java-based graphical user interface (GUI) coordinates a set of routines that introduce the influence of the membrane, determine its placement relative to the protein, and set the membrane potential. The software Jmol is embedded in the GUI to visualize the protein inserted in the membrane before the calculation and the electrostatic potential after completing the computation. We expect that the ease with which the GUI allows one to carry out these calculations will make this software a useful resource for experimenters and computational researchers alike. Three examples of membrane protein electrostatic calculations are carried out to illustrate how to use APBSmem and to highlight the different quantities of interest that can be calculated.

Solution and crystal molecular dynamics simulation study of m4-cyanovirin-N mutants complexed with di-mannose

Cyanovirin-N (CVN) is a highly potent anti-HIV carbohydrate-binding agent that establishes its microbicide activity through interaction with mannose-rich glycoprotein gp120 on the virion surface. The m4-CVN and P51G-m4-CVN mutants represent simple models for studying the high-affinity binding site, B(M). A recently determined 1.35 A high-resolution structure of P51G-m4-CVN provided details on the di-mannose binding mechanism, and suggested that the Arg-76 and Glu-41 residues are critical components of high mannose specificity and affinity. We performed molecular-dynamics simulations in solution and a crystal environment to study the role of Arg-76. Network analysis and clustering were used to characterize the dynamics of Arg-76. The results of our explicit solvent solution and crystal simulations showed a significant correlation with conformations of Arg-76 proposed from x-ray crystallographic studies. However, the crystal simulation showed that the crystal environment strongly biases conformational sampling of the Arg-76 residue. The solution simulations demonstrated no conformational preferences for Arg-76, which would support its critical role as the residue that locks the ligand in the bound state. Instead, a comparative analysis of trajectories from >50 ns of simulation for two mutants revealed the existence of a very stable eight-hydrogen-bond network between the di-mannose ligand and predominantly main-chain atoms. This network may play a key role in the specific recognition and strong binding of mannose oligomers in CVN and its homologs.

Model structure of APOBEC3C reveals a binding pocket modulating ribonucleic acid interaction required for encapsidation

Human APOBEC3 (A3) proteins form part of the intrinsic immunity to retroviruses. Carrying 1 or 2 copies of a cytidine deaminase motif, A3s act by deamination of retroviral genomes during reverse transcription. HIV-1 overcomes this inhibition by the Vif protein, which prevents incorporation of A3 into virions. In this study we modeled and probed the structure of APOBEC3C (A3C), a single-domain A3 with strong antilentiviral activity. The 3-dimensional protein model was used to predict the effect of mutations on antiviral activity, which was tested in a Deltavif simian immunodeficiency virus (SIV) reporter virus assay. We found that A3C activity requires protein dimerization for antiviral activity against SIV. Furthermore, by using a structure-based algorithm for automated pocket extraction, we detected a putative substrate binding pocket of A3C distal from the zinc-coordinating deaminase motif. Mutations in this region diminished antiviral activity by excluding A3C from virions. We found evidence that the small 5.8S RNA specifically binds to this locus and mediates incorporation of A3C into virus particles.

Insights into protein-polysorbate interactions analysed by means of isothermal titration and differential scanning calorimetry

Therapeutic proteins formulated as liquid solutions at high protein concentration are very sensitive to chemical and physical degradation. Especially avoiding the formation of protein aggregates is very crucial for product quality. In order to stabilize the colloidal properties of protein therapeutics various excipient are used. Especially the detergents polysorbate 20 and 80 are common. However, the mechanism upon which the detergents protect the protein from aggregation is not really known. The present study investigates the interaction of polysorbate 20 and 80 with different proteins: lysozyme, bovine serum albumin (BSA) and an immunoglobulin. The interaction and binding of the detergents to the proteins is investigated by isothermal titration calorimetry (ITC). From ITC the thermodynamic parameters (DeltaH: change in enthalpy, DeltaS: entropy and DeltaG: free energy) upon binding are derived as well as the binding constant K (a). The thermal stability of the proteins in the presence of the detergent is assessed by differential scanning calorimetry (DSC). The results show that both detergents bind to BSA with K (a) between 8 and 12 x 10(3) M(-1) with DeltaH -50 to -60 kJ/mol (25 degrees C). One to two detergent molecules bind to BSA. The presence of both detergents induces a weak stabilisation of the thermal denaturation properties of BSA. However, the interaction of polysorbate 20 and 80 with lysozyme and the immunoglobulin is quite negligible. The presence of the detergents up to a concentration of 2 mM has no impact on the heat capacity curve neither a destabilisation nor a stabilisation of the native conformation is observed.

Conformational dynamics of the mitochondrial ADP/ATP carrier: a simulation study

The mitochondrial ADP/ATP carrier is a six helix bundle membrane transport protein, which couples the exit of ATP from the mitochondrial matrix to the entry of ADP. Extended (4x20 ns) molecular dynamics simulations of the carrier, in the presence and absence of bound inhibitor (carboxyatractyloside), have been used to explore the conformational dynamics of the protein in a lipid bilayer environment, in the presence and absence of the carboxyatractyloside inhibitor. The dynamic flexibility (measured as conformational drift and fluctuations) of the protein is reduced in the presence of bound inhibitor. Proline residues in transmembrane helices H1, H3 and H5 appear to form dynamic hinges. Fluctuations in inter-helix salt bridges are also observed over the time course of the simulations. Inhibitor-protein and lipid-protein interactions have been characterised in some detail. Overall, the simulations support a transport mechanism in which flexibility about the proline hinges enables a transition between a 'closed' and an 'open' pore-like state of the carrier protein.

Structure refinement of the OpcA adhesin using molecular dynamics

OpcA from Neisseria meningitidis, the causative agent of meningococcal meningitis and septicemia, is an integral outer membrane protein that facilitates meningococcal adhesion through binding the proteoglycan receptors of susceptible cells. Two structures of OpcA have been determined by x-ray diffraction to 2 A resolution, revealing dramatically different conformations in the extracellular loops--the protein domain implicated in proteoglycan binding. In the first structure, a positively charged crevice formed by loops 1 and 2 was identified as the site for binding proteoglycans, whereas in the second structure the crevice was not evident as loops 1 and 2 adopted different conformations. To reconcile these results, molecular-dynamics simulations were carried out on both structures embedded in a solvated lipid bilayer membrane. Free of crystal contacts and crystallization agents, the loops were observed to undergo large structural transformations, suggesting that the conformation of the loops in either x-ray structure is affected by crystallization. Subsequent simulations of both structures in their crystal lattices confirmed this conclusion. Based on our molecular-dynamics trajectories, we propose a model for OpcA that combines stable structural features of the available x-ray structures. In this model, all five extracellular loops of OpcA have stable secondary structures. The loops form a funnel that leads to the base of the beta-barrel and that includes Tyr-169 on its exposed surface, which has been implicated in proteoglycan binding.

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.

OmpA: gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of approximately 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.