In meso crystal structure and docking simulations suggest an alternative proteoglycan binding site in the OpcA outer membrane adhesin

OpcA is an integral outer membrane adhesin protein from Neisseria meningitidis, the causative agent of meningococcal meningitis and septicemia. It binds to sialic acid (SA)-containing polysaccharides on the surface of epithelial cells. The crystal structure of OpcA showed that the protein adopts a 10-stranded beta-barrel structure, with five extensive loop regions on the extracellular side of the membrane. These form a crevice structure, lined with basic residues, which was hypothesized to act as the binding site for polysaccharide ligands. In the current study, a distinctly different OpcA structure has been obtained using crystals grown from a lipidic mesophase. Comparison of the two structures shows that the largest loop (L2), which closes over the end of the beta-barrel in the original crystal form, adopts a much more extended structure by reaching outward and away from the protein. The difference in conformation may be attributable to the absence of zinc ions from the crystallization conditions for the in meso crystal form: in the original structure, two zinc ions were bound to the external loops. Molecular dynamics (MD) simulations performed on the two OpcA models in a lipid bilayer environment demonstrated pronounced loop mobility. These observations support the view that the loop regions of OpcA are capable of a high degree of conformational flexibility. The original binding site for polysaccharide is not present in the in meso crystal form, and is disrupted during MD simulations. Docking analysis suggests a putative alternative location for the SA ligand in the new crystal form of OpcA.     

Membrane simulations of OpcA: gating in the loops?

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

The influenza fusion peptide adopts a flexible flat V conformation in membranes

Knowledge about the influenza fusion peptide (FP) membrane insertion mode is crucial for understanding its fusogenic mechanism. NMR and electron paramagnetic resonance experiments showed that in micelles, the FP inserted as a fixed-angle inverted V. In membranes, however, it was shown to insert as a straight α-helix (by molecular-dynamics simulations) and to adopt multiple kinked conformations (by solid-state NMR). In this work we performed explicit-solvent molecular-dynamics simulations of the influenza FP, and its F9A and W14A mutants, in POPC membranes. The Hα1 chemical shifts predicted from the molecular-dynamics structures are in excellent agreement with the experimental values obtained for the three peptides. The peptide orientation and conformations observed from the simulations lead to a flexible flat-V model in which the peptide lies almost flat on the membrane surface and alternates between kinked and straight-helix conformations.     

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 Å 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.     

Crystal structures of the Chromobacterium violaceumω-transaminase reveal major structural rearrangements upon binding of coenzyme PLP

The bacterial ω-transaminase from Chromobacterium violaceum (Cv-ωTA, EC2.6.1.18) catalyses industrially important transamination reactions by use of the coenzyme pyridoxal 5'-phosphate (PLP). Here, we present four crystal structures of Cv-ωTA: two in the apo form, one in the holo form and one in an intermediate state, at resolutions between 1.35 and 2.4 ?. The enzyme is a homodimer with a molecular mass of ～ 100 kDa. Each monomer has an active site at the dimeric interface that involves amino acid residues from both subunits. The apo-Cv-ωTA structure reveals unique 'relaxed' conformations of three critical loops involved in structuring the active site that have not previously been seen in a transaminase. Analysis of the four crystal structures reveals major structural rearrangements involving elements of the large and small domains of both monomers that reorganize the active site in the presence of PLP. The conformational change appears to be triggered by binding of the phosphate group of PLP. Furthermore, one of the apo structures shows a disordered 'roof?' over the PLP-binding site, whereas in the other apo form and the holo form the 'roof' is ordered. Comparison with other known transaminase crystal structures suggests that ordering of the 'roof' structure may be associated with substrate binding in Cv-ωTA and some other transaminases. DATABASE: The atomic coordinates and structure factors for the Chromobacterium violaceumω-transaminase crystal structures can be found in the RCSB Protein Data Bank (http://www.rcsb.org) under the accession codes 4A6U for the holoenzyme, 4A6R for the apo1 form, 4A6T for the apo2 form and 4A72 for the mixed form STRUCTURED DIGITAL ABSTRACT: ? -transaminases and -transaminases bind by dynamic light scattering (View interaction) ? -transaminase and -transaminase bind by x-ray crystallography (View interaction) ? -transaminase and -transaminase bind by x-ray crystallography (View interaction).     

Structure and interactions of RecA: plasticity revealed by molecular dynamics simulations

Eleven independent simulations, each involving three consecutive molecules in the RecA filament, carried out on the protein from Mycobacterium tuberculosis, Mycobacterium smegmatis and Escherichia coli and their Adenosine triphosphate (ATP) complexes, provide valuable information which is complementary to that obtained from crystal structures, in addition to confirming the robust common structural framework within which RecA molecules from different eubacteria function. Functionally important loops, which are largely disordered in crystal structures, appear to adopt in each simulation subsets of conformations from larger ensembles. The simulations indicate the possibility of additional interactions involving the P-loop which remains largely invariant. The phosphate tail of the ATP is firmly anchored on the loop while the nucleoside moiety exhibits substantial structural variability. The most important consequence of ATP binding is the movement of the ‘switch’ residue. The relevant simulations indicate the feasibility of a second nucleotide binding site, but the pathway between adjacent molecules in the filament involving the two nucleotide binding sites appears to be possible only in the mycobacterial proteins.     

Substrate-induced transmembrane signaling in the cobalamin transporter BtuB

The outer membranes of Gram-negative bacteria possess transport proteins essential for uptake of scarce nutrients. In TonB-dependent transporters, a conserved sequence of seven residues, the Ton box, faces the periplasm and interacts with the inner membrane TonB protein to energize an active transport cycle. A critical mechanistic step is the structural change in the Ton box of the transporter upon substrate binding; this essential transmembrane signaling event increases the affinity of the transporter for TonB and enables active transport to proceed. We have solved crystal structures of BtuB, the outer membrane cobalamin transporter from Escherichia coli, in the absence and presence of cyanocobalamin (vitamin B(12)). In these structures, the Ton box is ordered and undergoes a conformational change in the presence of bound substrate. Calcium has been implicated as a necessary factor for the high-affinity binding (K(d) approximately 0.3 nM) of cyanocobalamin to BtuB. We observe two bound calcium ions that order three extracellular loops of BtuB, thus providing a direct (and unusual) structural role for calcium.     

Recognition of saccharides by the OpcA, OpaD, and OpaB outer membrane proteins from Neisseria meningitidis

The adhesion of the pathogen Neisseria meningitidis to host cell surface proteoglycan, mediated by the integral outer membrane proteins OpcA and Opa, plays an important part in the processes of colonization and invasion by the bacterium. The precise specificities of the OpcA and Opa proteins are, however, unknown. Here we use a fluorescence-based binding assay to show that both proteins bind to mono- and disaccharides with high affinity. Binding of saccharides caused a quench in the intrinsic fluorescence emission of both proteins, and mutation of selected Tyr residues within the external loop regions caused a substantial decrease in fluorescence. We suggest that the intrinsic fluorescence arises from resonance energy transfer from Tyr to Trp residues in the beta-barrel portion of the structure. OpcA bound sialic acid with a Kd of 0.31 microM and was shown to be specific for pyranose saccharides. The binding specificities of two different Opa proteins were compared; unlike OpcA, neither protein bound to monosaccharides, but both bound to maltose, lactose, and sialic acid-containing oligosaccharides, with Kd values in the micromolar range. OpaB had a 10-fold higher affinity for sialic acid-containing ligands than OpaD as a result of the mutation Y165V, which was shown to restore this specificity to OpaD. Finally, the OpcA- and Opa-dependent adhesion of meningococci to epithelial cells was shown to be partially inhibited by exogenously added sialic acid and maltose. The results show that OpcA and the Opa proteins can be thought of as outer membrane lectins and that simple saccharides can modulate their recognition of complex proteoglycan receptors.     

Structural flexibility in Trypanosoma brucei enolase revealed by X-ray crystallography and molecular dynamics

Enolase is a validated drug target in Trypanosoma brucei. To better characterize its properties and guide drug design efforts, we have determined six new crystal structures of the enzyme, in various ligation states and conformations, and have carried out complementary molecular dynamics simulations. The results show a striking structural diversity of loops near the catalytic site, for which variation can be interpreted as distinct modes of conformational variability that are explored during the molecular dynamics simulations. Our results show that sulfate may, unexpectedly, induce full closure of catalytic site loops whereas, conversely, binding of inhibitor phosphonoacetohydroxamate may leave open a tunnel from the catalytic site to protein surface offering possibilities for drug development. We also present the first complex of enolase with a novel inhibitor 2-fluoro-2-phosphonoacetohydroxamate. The molecular dynamics results further encourage efforts to design irreversible species-specific inhibitors: they reveal that a parasite enzyme-specific lysine may approach the catalytic site more closely than crystal structures suggest and also cast light on the issue of accessibility of parasite enzyme-specific cysteines to chemically modifying reagents. One of the new sulfate structures contains a novel metal-binding site IV within the catalytic site cleft.     

MD Simulations and FRET Reveal an Environment-Sensitive Conformational Plasticity of Importin-β

The nuclear pore complex mediates nucleocytoplasmic transport of macromolecules in eukaryotic cells. Transport through the pore is restricted by a hydrophobic selectivity filter comprising disordered phenylalanine-glycine-rich repeats of nuclear pore proteins. Exchange through the pore requires specialized transport receptors, called exportins and importins, that interact with cargo proteins in a RanGTP-dependent manner. These receptors are highly flexible superhelical structures composed of HEAT-repeat motifs that adopt various degrees of extension in crystal structures. Here, we performed molecular-dynamics simulations using crystal structures of Importin-β in its free form or in complex with nuclear localization signal peptides as the starting conformation. Our simulations predicted that initially compact structures would adopt extended conformations in hydrophilic buffers, while contracted conformations would dominate in more hydrophobic solutions, mimicking the environment of the nuclear pore. We confirmed this experimentally by Förster resonance energy transfer experiments using dual-fluorophore-labeled Importin-β. These observations explain seemingly contradictory crystal structures and suggest a possible mechanism for cargo protection during passage of the nuclear pore. Such hydrophobic switching may be a general principle for environmental control of protein function.     

Solutes alter the conformation of the ligand binding loops in outer membrane transporters

The binding and recognition of ligands by bacterial outer membrane transport proteins is mediated in part by interactions made through their extracellular loops. Here, site-directed spin labeling (SDSL) and electron paramagnetic resonance (EPR) spectroscopy were used to examine the effect of stabilizing solutes on the extracellular loops in BtuB, the vitamin B12 transporter, and FecA, the ferric citrate transporter. EPR spectra from the extracellular loops of FecA and BtuB arise from dynamic backbone segments, and distance measurements made by double electron-electron resonance indicate that the second extracellular loop in BtuB samples a wide range of conformations. These conformations are dramatically restricted upon substrate binding. In addition, the EPR spectra from nitroxide labels attached to the extracellular loops in BtuB and FecA are highly sensitive to solutes, and at every site examined the motion of the label is significantly reduced in the presence of stabilizing osmolytes, such as polyethylene glycols. For the second extracellular loop in BtuB, the solute-induced structural changes are small, but they are sufficient to bring spin-labeled side chains into tertiary contact with other portions of the protein. The spectroscopic changes seen by SDSL suggest that high concentrations of stabilizing solutes, such as those used to generate membrane protein crystals, result in a more compact and ordered state of the protein than is seen under more physiological conditions.     

Study of structural dynamics of ligand-activated membrane receptors by means of principal component analysis

The structural dynamics of three different ligand-activated G-protein coupled receptors (GPCRs) and the photoreactive receptor rhodopsin from mammals were comparatively studied. As a result, diagrams demonstrating the main structural differences between the studied membrane receptors were obtained. These diagrams represent the projection of the crystal structures of rhodopsin photointermediates and ligand-activated receptors onto the plane defined by the principal components. Thus, we were able to associate the activation process of the receptors with large-scale movements of their individual transmembrane (TM) domains. In addition, the dynamics of extracellular loops of ligand-activated receptors responsible for recognition and initial binding of ligands was studied. Based on these results, two parameters of functionally significant structural dynamics of membrane receptors can be thoroughly analyzed simultaneously — movements of individual TM helices and of extracellular loops.     

Membrane Interaction of Bound Ligands Contributes to the Negative Binding Cooperativity of the EGF Receptor

The epidermal growth factor receptor (EGFR) plays a key role in regulating cell proliferation, migration, and differentiation, and aberrant EGFR signaling is implicated in a variety of cancers. EGFR signaling is triggered by extracellular ligand binding, which promotes EGFR dimerization and activation. Ligand-binding measurements are consistent with a negatively cooperative model in which the ligand-binding affinity at either binding site in an EGFR dimer is weaker when the other site is occupied by a ligand. This cooperativity is widely believed to be central to the effects of ligand concentration on EGFR-mediated intracellular signaling. Although the extracellular portion of the human EGFR dimer has been resolved crystallographically, the crystal structures do not reveal the structural origin of this negative cooperativity, which has remained unclear. Here we report the results of molecular dynamics simulations suggesting that asymmetrical interactions of the two binding sites with the membrane may be responsible (perhaps along with other factors) for this negative cooperativity. In particular, in our simulations the extracellular domains of an EGFR dimer spontaneously lay down on the membrane in an orientation in which favorable membrane contacts were made with one of the bound ligands, but could not be made with the other. Similar interactions were observed when EGFR was glycosylated, as it is in vivo.     

Structure and dynamics of parallel beta-sheets, hydrophobic core, and loops in Alzheimer's A beta fibrils

We explore the relative contributions of different structural elements to the stability of Abeta fibrils by molecular-dynamics simulations performed over a broad range of temperatures (298 K to 398 K). Our fibril structures are based on solid-state nuclear magnetic resonance experiments of Abeta(1-40) peptides, with sheets of parallel beta-strands connected by loops and stabilized by interior salt bridges. We consider models with different interpeptide interfaces, and different staggering of the N- and C-terminal beta-strands along the fibril axis. Multiple 10-20 ns molecular-dynamics simulations show that fibril segments with 12 peptides are stable at ambient temperature. The different models converge toward an interdigitated side-chain packing, and present water channels solvating the interior D23/K28 salt bridges. At elevated temperatures, we observe the early phases of fibril dissociation as a loss of order in the hydrophilic loops connecting the two beta-strands, and in the solvent-exposed N-terminal beta-sheets. As the most dramatic structural change, we observe collective sliding of the N- and C-terminal beta-sheets on top of each other. The interior C-terminal beta-sheets in the hydrophobic core remain largely intact, indicating that their formation and stability is crucial to the dissociation/elongation and stability of Abeta fibrils.     

Interactions of the C2 domain of human factor V with a model membrane

Activated coagulation Factor V is an important cofactor of the coagulation cascade that catalyzes the formation of the prothrombinase complex on the surface of membranes rich in phosphatidyl-L-serine (PS). Here we report molecular dynamics simulations of the two crystallographic structures (the open and closed conformations) of domain C2 of coagulation Factor V (FaVC2). The calculations were performed in water (1.5 ns for each conformation) and in the presence of a neutral phospholipid bilayer model (POPE; 10 ns for each conformation) in order to describe the dynamics of the free (plasma circulating) and membrane bound forms of FaVC2. Water simulations confirmed the hypothesis that the plasma circulating form is in the closed conformation. In contrast, the membrane simulations showed that both conformations are energetically compatible with membrane binding. We have investigated the mechanism, the dynamics, and the energetics of the binding process. Our data are consistent with published estimates of the immersion depth of the aromatic residues (W26 and W27), and with mutagenesis studies involving specific residues located on the spikes at the bottom of the FaVC2 structure. Electrostatic interactions between the phospholipid head groups and hydrophilic residues at the bottom of the structure play a key role in the binding process by creating a large number of hydrogen bonds that anchor the protein to the membrane. The simulations identified a stable phospholipid binding pocket reminiscent of a previously suggested PS interaction site. Our structural data could contribute to the design of potential inhibitors able to disrupt membrane association.     

Catalytic transitions in the human MDR1 P-glycoprotein drug binding sites

Multidrug resistance proteins that belong to the ATP-binding cassette family like the human P-glycoprotein (ABCB1 or Pgp) are responsible for many failed cancer and antiviral chemotherapies because these membrane transporters remove the chemotherapeutics from the targeted cells. Understanding the details of the catalytic mechanism of Pgp is therefore critical to the development of inhibitors that might overcome these resistances. In this work, targeted molecular dynamics techniques were used to elucidate catalytically relevant structures of Pgp. Crystal structures of homologues in four different conformations were used as intermediate targets in the dynamics simulations. Transitions from conformations that were wide open to the cytoplasm to transition state conformations that were wide open to the extracellular space were studied. Twenty-six nonredundant transitional protein structures were identified from these targeted molecular dynamics simulations using evolutionary structure analyses. Coupled movement of nucleotide binding domains (NBDs) and transmembrane domains (TMDs) that form the drug binding cavities were observed. Pronounced twisting of the NBDs as they approached each other as well as the quantification of a dramatic opening of the TMDs to the extracellular space as the ATP hydrolysis transition state was reached were observed. Docking interactions of 21 known transport ligands or inhibitors were analyzed with each of the 26 transitional structures. Many of the docking results obtained here were validated by previously published biochemical determinations. As the ATP hydrolysis transition state was approached, drug docking in the extracellular half of the transmembrane domains seemed to be destabilized as transport ligand exit gates opened to the extracellular space.     

Characterizing and predicting the functional and conformational diversity of seven-transmembrane proteins

The activation of seven-transmembrane receptors (7TMRs) allows cells to sense their environment and convert extracellular signals (like hormone binding) into intracellular signals (through G protein-coupled and/or β arrestin-coupled pathways). A single 7TMR is capable of transducing a wide spectrum of physiological responses inside a cell by coupling to these pathways. This intracellular pleiotropic action is enabled by multiple conformations exhibited by these receptors. Developments in membrane protein structure determination technologies have led to a rapid increase in crystal structures for many 7TMRs. Majority of these receptors have been crystallized in their inactive conformation and, for some, one of the many active conformations has also been crystallized. Given the topological constraints of a lipid bilayer that results in a single fold of seven almost parallel TM helices connected by mostly unstructured loops, these structures exhibit a diversity of conformations not only across the receptors but also across the different functional forms for receptors with structures for one of the functionally active conformations. Here we present a method to characterize this conformational diversity in terms of transmembrane helix topology (TMHTOP) parameters and how to use these helix orientation parameters to predict functionally-distinct multiple conformations for these receptors. The TMHTOP parameters enable a quantification of the structural changes that underlie 7TMR activation and also sheds a unique mechanistic light on the pleiotropic nature of these receptors. It provides a common language to describe the 7TMR activation mechanisms as well as differences across many receptors in terms of visually intuitive structural parameters. Protein structure prediction methods can use these parameters to describe 7TMR conformational ensembles, which coupled to experimental data can be used to develop testable hypotheses for the structural basis of 7TMR functions.     

Crystallographic structure of SurA, a molecular chaperone that facilitates folding of outer membrane porins

The SurA protein facilitates correct folding of outer membrane proteins in gram-negative bacteria. The sequence of Escherichia coli SurA presents four segments, two of which are peptidyl-prolyl isomerases (PPIases); the crystal structure reveals an asymmetric dumbbell, in which the amino-terminal, carboxy-terminal, and first PPIase segments of the sequence form a core structural module, and the second PPIase segment is a satellite domain tethered approximately 30 A from this module. The core module, which is implicated in membrane protein folding, has a novel fold that includes an extended crevice. Crystal contacts show that peptides bind within the crevice, suggesting a model for chaperone activity whereby segments of polypeptide may be repetitively sequestered and released during the membrane protein-folding process.     

G 蛋白偶联受体寡聚体结构研究进展

G蛋白偶联受体(GPCR)是细胞膜上最大的一类受体,其通过构象变化激活下游G蛋白从而介导细胞响应多种来自内源和外界环境中的信号。自GPCR被发现以来,研究者就一直在努力解析GPCR的构象,x射线晶体衍射技术和GPCR蛋白质结晶技术的发展使得越来越多的GPCR单体在静息状态,以及与不同配体甚至G蛋白结合的晶体结构被成功解析。另一方面,FRET和电子显微技术的运用得到了GPCR二聚化和多聚化的多方面证据。本文将结合近年来该领域的进展,对GPCR寡聚体的结构和构象变化予以系统的综述,这些成果为研究GPCR的功能机制及其特异性的靶点药物开发提供了重要的基础。     

Anti-insulin antibody structure and conformation. II. Molecular dynamics with explicit solvent.

Molecular dynamics at 300 K was used as a conformation searching tool to analyze a knowledge-based structure prediction of an anti-insulin antibody. Solvation effects were modeled by packing water molecules around the antigen binding loops. Some loops underwent backbone and side-chain conformational changes during the 95-ps equilibration, and most of these new, lower potential energy conformations were stable during the subsequent 200-ps simulation. Alterations to the model include changes in the intraloop, main-chain hydrogen bonding network of loop H3, and adjustments of Tyr and Lys side chains of H3 induced by hydrogen bonding to water molecules. The structures observed during molecular dynamics support the conclusion of the previous paper that hydrogen bonding will play the dominant role in antibody-insulin recognition. Determination of the structure of the antibody by x-ray crystallography is currently being pursued to provide an experimental test of these results. The simulation appears to improve the model, but longer simulations at higher temperatures should be performed.