Crystal structure of osmoporin OmpC from E. coli at 2.0 A

Porins form transmembrane pores in the outer membrane of Gram-negative bacteria with matrix porin OmpF and osmoporin OmpC from Escherichia coli being differentially expressed depending on environmental conditions. The three-dimensional structure of OmpC has been determined to 2.0 A resolution by X-ray crystallography. As expected from the high sequence similarity, OmpC adopts the OmpF-like 16-stranded hollow beta-barrel fold with three beta-barrels associated to form a tight trimer. Unlike in OmpF, the extracellular loops form a continuous wall at the perimeter of the vestibule common to the three pores, due to a 14-residues insertion in loop L4. The pore constriction and the periplasmic outlet are very similar to OmpF with 74% of the pore lining residues being conserved. Overall, only few ionizable residues are exchanged at the pore lining. The OmpC structure suggests that not pore size, but electrostatic pore potential and particular atomic details of the pore linings are the critical parameters that physiologically distinguish OmpC from OmpF.     

Intermolecular protein interactions in solutions of bovine lens beta L-crystallin. Results from 1/T1 nuclear magnetic relaxation dispersion profiles.

We report the magnetic field dependence of 1/T1 of solvent water protons and deuterons (nuclear magnetic relaxation dispersion, or NMRD, profiles) for solutions of steer lens beta L-crystallin. Such data allow the study of intermolecular protein interactions over a wide concentration range, here 1-34% vol/vol, by providing a measure of the rotational relaxation time of solute macromolecules. We conclude that, for approximately less than 5% protein, the solute particles are noncompact, with a rotationally averaged volume approximately three times that of a compact 60-kD sphere. (Earlier results for alpha-crystallin, approximately 1,000 kD, from optical and osmotic measurements (Vérétout and Tardieu, 1989. J. Mol. Biol. 205:713-728), show a similar, approximately twofold, effect). At intermediate concentrations, to approximately 20% protein, there is evidence for limited association or oligomerization, as found for the structurally related gamma II-crystallin (Koenig et al. 1990. Biophys. J. 57:461-469), to a limiting size about two-thirds that of alpha-crystallin. The difference in NMRD behavior of the three classes of crystallins is consonant with their differing osmotic properties (Vérétout and Tardieu. J. Mol. Biol. 1989, 205:713-728; Kenworthy, McIntosh, and Magid. Biophys. J. 1992. 61:A477; Tardieu et al. 1992. Eur. Biophys. J. 21:1-12). We indicate how the unusual structures and interactions of these three classes of proteins can be combined to optimize transparency and minimize colloid osmotic difficulties in eye lens.     

Dynamics and energetics of water permeation through the aquaporin channel

Structural properties of water inside bovine aquaporin-1 are investigated by molecular simulation. The calculations, which are based on the recently determined X-ray structure at 2.2 A resolution (Sui et al., Nature 2001;414:872-878), are carried out on one monomeric subunit immersed in a water-n-octane-water bilayer. Molecular dynamics (MD) simulations suggest that His182, a fully conserved residue in the channel pore, is protonated in the delta position. Furthermore, they reveal a highly ordered water structure in the channel, induced by the electrostatic properties of the protein. Multiple-steering MD simulations are used to calculate the free-energy of water diffusion. To the best of our knowledge, this represents the first free-energy calculation based on the new, high-resolution structure of the pore. The calculated barrier is 2.5 kcal/mol, and it is associated to water permeation through the Asn-Pro-Ala (NPA) region of the pore, where water molecules are only hydrogen-bonded with themselves. These findings are fully consistent with those based on the previous MD studies on the human protein (de Groot and Grubmüller, Science 2001;294:2353-2357).     

Localization of a voltage gate in connexin46 gap junction hemichannels.

Cysteine replacement mutagenesis has identified positions in the first transmembrane domain of connexins as contributors to the pore lining of gap junction hemichannels (Zhou et al. 1997. Biophys. J. 72:1946-1953). Oocytes expressing a mutant cx46 with a cysteine in position 35 exhibited a membrane conductance sensitive to the thiol reagent maleimidobutyryl biocytin (MBB). MBB irreversibly reduced the single-channel conductance by 80%. This reactive cysteine was used to probe the localization of a voltage gate that closes cx46 gap junction hemichannels at negative potentials. MBB was applied to the closed channel either from outside (whole cell) or from inside (excised membrane patches). After washout of the thiol reagent the channels were tested at potentials at which the channels open. After extracellular application of MBB to intact oocytes, the membrane conductance was unaffected. In contrast, channels treated with intracellular MBB were blocked. Thus the cysteine in position 35 of cx46 is accessible from inside but not from the outside while the channel is closed. These results suggest that the voltage gate, which may be identical to the "loop gate" (Trexler et al. 1996. Proc. Natl. Acad. Sci. USA. 93:5836-5841), is located extracellular to the 35 position. The voltage gate results in regional closure of the pore rather than closure along the entire pore length.     

An alamethicin channel in a lipid bilayer: molecular dynamics simulations

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.     

The current-voltage relation of an aqueous pore in a lipid bilayer membrane

In a recent paper, Glaser et al. [1988) Biochim. Biophys. Acta 940, 275-287) use an inappropriate formula for the current-voltage relation of a pore to analyze the data from an electroporation experiment. They proceed to use their data analysis to support their contention that the voltage dependence of the membrane conductivity is caused entirely by the nonlinearity of the current-voltage relation of the pores, changes in pore size and number being unimportant. This contention is no longer justified. The formula that Glaser et al. use is valid for a one-dimensional problem, but they apply it to a pore that has a three-dimensional structure. Although the corresponding one-dimensional problem can be reduced to quadratures, only upper and lower bounds on the conductivity can be found for the three-dimensional problem without solving a partial differential equation.     

A new gating site in human aquaporin-4: Insights from molecular dynamics simulations

Aquaporin-4 (AQP4) is the predominant water channel in different organs and tissues. An alteration of its physiological functioning is responsible for several disorders of water regulation and, thus, is considered an attractive target with a promising therapeutic and diagnostic potential. Molecular dynamics (MD) simulations performed on the AQP4 tetramer embedded in a bilayer of lipid molecules allowed us to analyze the role of spontaneous fluctuations occurring inside the pore. Following the approach by Hashido et al. [Hashido M, Kidera A, Ikeguchi M (2007) Biophys J 93: 373–385], our analysis on 200 ns trajectory discloses three domains inside the pore as key elements for water permeation. Herein, we describe the gating mechanism associated with the well-known selectivity filter on the extracellular side of the pore and the crucial regulation ensured by the NPA motifs (asparagine, proline, alanine). Notably, on the cytoplasmic side, we find a putative gate formed by two residues, namely, a cysteine belonging to the loop D (C178) and a histidine from loop B (H95). We observed that the spontaneous reorientation of the imidazole ring of H95 acts as a molecular switch enabling H-bond interaction with C178. The occurrence of such local interaction seems to be responsible for the narrowing of the pore and thus of a remarkable decrease in water flux rate. Our results are in agreement with recent experimental observations and may represent a promising starting point to pave the way for the discovery of chemical modulators of AQP4 water permeability.     

Molecular dynamics simulation of sucrose- and trehalose-coated carboxy-myoglobin

We performed a room temperature molecular dynamics (MD) simulation on a system containing 1 carboxy-myoglobin (MbCO) molecule in a sucrose-water matrix of identical composition (89% [sucrose/(sucrose + water)] w/w) as for a previous trehalose-water-MbCO simulation (Cottone et al., Biophys J 2001;80:931-938). Results show that, as for trehalose, the amplitude of protein atomic mean-square fluctuations, on the nanosecond timescale, is reduced with respect to aqueous solutions also in sucrose. A detailed comparison as a function of residue number evidences mobility differences along the protein backbone, which can be related to a different efficacy in bioprotection. Different heme pocket structures are observed in the 2 systems. The joint distribution of the magnitude of the electric field at the CO oxygen atom and of the angle between the field and the CO unit vector shows a secondary maximum in sucrose, absent in trehalose. This can explain the CO stretching band profile (A substates distribution) differences evidenced by infrared spectroscopy in sucrose- and trehalose-coated MbCO (Giuffrida et al., J Phys Chem B 2004;108:15415-15421), and in particular the appearance of a further substate in sucrose. Analysis of hydrogen bonds at the protein-solvent interface shows that the fraction of water molecules shared between the protein and the sugar is lower in sucrose than in trehalose, in spite of a larger number of water molecules bound to the protein in the former system, thus indicating a lower protein-matrix coupling, as recently observed by Fourier transform infrared (FTIR) experiments (Giuffrida et al., J Phys Chem B 2004;108:15415-15421).     

Supramolecular structures of peptide assemblies in membranes by neutron off-plane scattering: method of analysis

In a previous paper (Yang et al., Biophys. J. 75:641-645, 1998), we showed a simple, efficient method of recording the diffraction patterns of supramolecular peptide assemblies in membranes where the samples were prepared in the form of oriented multilayers. Here we develop a method of analysis based on the diffraction theory of two-dimensional liquids. Gramicidin was used as a prototype model because its pore structure in membrane in known. At full hydration, the diffraction patterns of alamethicin and magainin are similar to gramicidin except in the scale of q (the momentum transfer of scattering), clearly indicating that both alamethicin and magainin form pores in membranes but of different sizes. When the hydration of the multilayer samples was decreased while the bilayers were still fluid, the in-plane positions of the membrane pores became correlated from one bilayer to the next. We believe that this is a new manifestation of the hydration force. The effect is most prominent in magainin patterns, which are used to demonstrate the method of analysis. When magainin samples were further dehydrated or cooled, the liquid-like diffraction turned into crystal-like patterns. This discovery points to the possibility of investigating the supramolecular structures with high-order diffraction.     

Pore formation by LamB of Escherichia coli in lipid bilayer membranes.

Lipid bilayer experiments were performed in the presence of different Escherichia coli LamB preparations. These LamB preparations formed two types of pores in the membranes. Large pores, which had a single-channel conductance of 2.7 nS and comprised about 1 to 6% of the total pores, were presumably contaminants which might have been induced together with LamB. LamB itself formed small pores with a single-channel conductance of 160 pS in 1 M KCl. These pores could be completely blocked by the addition of maltose and maltodextrins. Titration of the pore conductance with maltotriose suggested that there was a binding site inside the pores with a Ks of 2.5 X 10(-4) M for maltotriose. On the basis of our data we concluded that the structure of the LamB channels is quite different from the structures of the channels of general diffusion porins, such as OmpF and OmpC.     

Dynamics of heteropentameric nicotinic acetylcholine receptor: implications of the gating mechanism

The dynamics characteristics of the currently available structure of Torpedo nicotinic acetylcholine receptor (nAChR), including the extracellular, transmembrane, and intracellular domains (ICDs), were analyzed using the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We found that a symmetric quaternary twist motion, reported previously in the literature in a homopentameric receptor (Cheng et al. J Mol Biol 2006;355:310-324; Taly et al. Biophys J 2005;88:3954-3965), occurred also in the heteropentameric Torpedo nAChR. We believe, however, that the symmetric twist alone is not sufficient to explain a large body of experimental data indicating asymmetry and subunit nonequivalence during gating. Here we report our results supporting the hypothesis that a combination of symmetric and asymmetric motions opens the gate. We show that the asymmetric motion involves tilting of the TM2 helices. Furthermore, our study reveals three additional aspects of channel dynamics: (1) loop A serves as an allosteric mediator between the ligand binding loops and those at the domain interface, particularly the linker between TM2 and TM3; (2) the ICD can modulate the pore dynamics and thus should not be neglected in gating studies; and (3) the F loops, which are peculiarly longer and poorly-conserved in non-alpha-subunits, have important dynamical implications.     

Properties of chemically modified porin from Escherichia coli in lipid bilayer membranes

Purified porin OmpF from Escherichia coli outer membrane was chemically modified by acetylation and succinylation of amino groups and by amidation of the carboxyl groups. Native and chemically modified porins were incorporated into lipid bilayer membranes and the permeability properties of the pores were studied. Acetylation and succinylation of the porin trimers had almost no influence on the single channel conductance in the presence of small cations and anions and the cation selectivity remained essentially unchanged as compared with the native porin. Amidation had also only little influence on the single channel conductance and changed the pore conductance at maximum by less than 50%, whereas the cation selectivity of the porin is completely lost after amidation. The results suggest that the structure of the porin pore remains essentially unchanged after chemical modification of the pores and that their cation selectivity is caused by an excess of negatively charged groups inside the pore and/or on the surface of the protein. Furthermore, it seems very unlikely that the pore contains any positively charged group at neutral pH.     

"Opening" the ferritin pore for iron release by mutation of conserved amino acids at interhelix and loop sites.

Ferritin concentrates, stores, and detoxifies iron in most organisms. The iron is a solid, ferric oxide mineral (< or =4500 Fe) inside the protein shell. Eight pores are formed by subunit trimers of the 24 subunit protein. A role for the protein in controlling reduction and dissolution of the iron mineral was suggested in preliminary experiments [Takagi et al. (1998) J. Biol. Chem. 273, 18685-18688] with a proline/leucine substitution near the pore. Localized pore disorder in frog L134P crystals coincided with enhanced iron exit, triggered by reduction. In this report, nine additional substitutions of conserved amino acids near L134 were studied for effects on iron release. Alterations of a conserved hydrophobic pair, a conserved ion pair, and a loop at the ferritin pores all increased iron exit (3-30-fold). Protein assembly was unchanged, except for a slight decrease in volume (measured by gel filtration); ferroxidase activity was still in the millisecond range, but a small decrease indicates slight alteration of the channel from the pore to the oxidation site. The sensitivity of reductive iron exit rates to changes in conserved residues near the ferritin pores, associated with localized unfolding, suggests that the structure around the ferritin pores is a target for regulated protein unfolding and iron release.     

Transmembrane peptide NB of influenza B: a simulation, structure, and conductance study

The putative transmembrane segment of the ion channel forming peptide NB from influenza B was synthesized by standard solid-phase peptide synthesis. Insertion into the planar lipid bilayer revealed ion channel activity with conductance levels of 20, 61, 107, and 142 pS in a 0.5 M KCl buffer solution. In addition, levels at -100 mV show conductances of 251 and 413 pS. A linear current-voltage relation reveals a voltage-independent channel formation. In methanol and in vesicles the peptide appears to adopt an alpha-helical-like structure. Computational models of alpha-helix bundles using N = 4, 5, and 6 NB peptides per bundle revealed water-filled pores after 1 ns of MD simulation in a solvated lipid bilayer. Calculated conductance values [using HOLE (Smart et al. (1997) Biophys. J. 72, 1109-1126)] of ca. 20, 60, and 90 pS, respectively, suggested that the multiple conductance levels seen experimentally must correspond to different degrees of oligomerization of the peptide to form channels.     

Characterizing the secondary hydration shell on hydrated myoglobin, hemoglobin, and lysozyme powders by its vitrification behavior on cooling and its calorimetric glass-->liquid transition and crystallization behavior on reheating.

For hydrated metmyoglobin, methemoglobin, and lysozyme powders, the freezable water fraction of between approximately 0.3-0.4 g water/g protein up to approximately 0.7-0.8 g water/g protein has been fully vitrified by cooling at rates up to approximately 1500 K min-1 and the influence of cooling rate characterized by x-ray diffractograms. This vitreous but freezable water fraction started to crystallize at approximately 210 K to cubic ice and at approximately 240 K to hexagonal ice. Measurements by differential scanning calorimetry have shown that this vitreous but freezable water fraction undergoes, on reheating at a rate of 30 K min-1, a glass-->liquid transition with an onset temperature of between approximately 164 and approximately 174 K, with a width of between approximately 9 and approximately 16 degrees and an increase in heat capacity of between approximately 20 and approximately 40 J K-1 (mol of freezable water)-1 but that the glass transition disappears upon crystallization of the freezable water. These calorimetric features are similar to those of water imbibed in the pores of a synthetic hydrogel but very different from those of glassy bulk water. The difference to glassy bulk water's properties is attributed to hydrophilic interaction and H-bonding of the macromolecules' segments with the freezable water fraction, which thereby becomes dynamically modified. Abrupt increase in minimal or critical cooling rate necessary for complete vitrification is observed at approximately 0.7-0.8 g water/g protein, which is attributed to an abrupt increase of water's mobility, and it is remarkably close to the threshold value of water's mobility on a hydrated protein reported by Kimmich et al. (1990, Biophys. J. 58:1183). The hydration level of approximately 0.7-0.8 g water/g protein is approximately that necessary for completing the secondary hydration shell.     

Ions and counterions in a biological channel: a molecular dynamics simulation of OmpF porin from Escherichia coli in an explicit membrane with 1 M KCl aqueous salt solution

A 5 ns all-atom molecular dynamics trajectory of Escherichia coli OmpF porin embedded in an explicit dimyristoyl-phosphatidylcholine (DMPC) bilayer bathed by a 1 M [KCl] aqueous salt solution is generated to explore the microscopic details of the mechanism of ion permeation. The atomic model includes the OmpF trimer, 124 DMPC, 13470 water molecules as well as 231 K+ and 201 Cl-, for a total of 70,693 atoms. The structural and dynamical results are in excellent agreement with the X-ray data. The global root-mean-square deviation of the backbone atoms relative to the X-ray structure is 1.4 A. A cluster of three fully charged arginine (Arg42, Arg82, and Arg132) facing two acidic residues (Asp113 and Glu117) on L3 in the narrowest part of the aqueous pore is observed to be very stable in the crystallographic conformation. In this region of the pore, the water molecules are markedly oriented perpendicular to the channel axis due to the strong transversal electrostatic field arising from those residues. On average the size of the pore is smaller during the simulation than in the X-ray structure, undergoing small fluctuations. No large movements of loop L3 leading to a gating of the pore are observed. Remarkably, it is observed that K+ and Cl- follow two well-separated average pathways spanning over nearly 40 A along the axis of the pore. In the center of the monomer, the two screw-like pathways have a left-handed twist, undergoing a counter-clockwise rotation of 180 degrees from the extracellular vestibule to the pore periplasmic side. In the pore, the dynamical diffusion constants of the ions are reduced by about 50% relative to their value in bulk solvent. Analysis of ion solvation across the channel reveals that the contributions from the water and the protein are complementary, keeping the total solvation number of both ions nearly constant. Unsurprisingly, K+ have a higher propensity to occupy the aqueous pore than Cl-, consistent with the cation selectivity of the channel. However, further analysis suggests that ion-ion pairs play an important role. In particular, it is observed that the passage of Cl- occurs only in the presence of K+ counterions, and isolated K+ can move through the channel and permeate on their own. The presence of K+ in the pore screens the negative electrostatic potential arising from OmpF to help the translocation of Cl- by formation of ion pairs.     

FTIR spectroscopic study of the dynamics of conformational substates in hydrated carbonyl-myoglobin films via temperature dependence of the CO stretching band parameters.

Two hydrated carbonyl myoglobin (MbCO) films, one containing (0.30 g water)/(g MbCO) from MbCO solution in water at pH 5.5 and the other (0.32 g water)/(gMbCO) from 0.1 M potassium phosphate buffer solution at pH 6.8, were studied by FTIR spectroscopy from 293 K to 78 K at selected temperatures on cooling and reheating. Above approximately 180 K the general trend in temperature dependence of half-bandwidths, peak maxima, and band area ratios of the A1 and A3 conformer bands is similar to those reported by Ansari et al. (1987. Biophys. J. 26:337) for MbCO in 75% glycerol/water solution, but abrupt changes in slopes at approximately 180-200 K and freezing-in of conformer populations, which could be taken as indicator for glass transition of the solvent or the protein, are absent for the hydrated MbCO films. This is interpreted in terms of an exceptionally broad distribution of relaxation times, and is in accord with conclusions from recent calorimetric annealing studies of hydrated protein powders (Sartor et al. 1994. Biophys. J. 66:249). Exchange between the three A conformers does not stop at approximately 180-200 K but occurs over the whole temperature region studied. These results are then discussed with respect to MbCO's behavior in the glass-->liquid transition region of glass-forming solvents, and it is concluded that, in analogy to the behavior of low-molecular-weight compounds with a distribution of rapidly interconverting conformers, freezing-in of MbCO's A conformer populations by the solvent should not be mistaken for a glass transition of MbCO.     

A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels

A computational algorithm based on Grand Canonical Monte Carlo (GCMC) and Brownian Dynamics (BD) is described to simulate the movement of ions in membrane channels. The proposed algorithm, GCMC/BD, allows the simulation of ion channels with a realistic implementation of boundary conditions of concentration and transmembrane potential. The method is consistent with a statistical mechanical formulation of the equilibrium properties of ion channels (; Biophys. J. 77:139-153). The GCMC/BD algorithm is illustrated with simulations of simple test systems and of the OmpF porin of Escherichia coli. The approach provides a framework for simulating ion permeation in the context of detailed microscopic models.     

A multiscale approach to the elastic moduli of biomembrane networks

We develop equilibrium fluctuation formulae for the isothermal elastic moduli of discrete biomembrane models at different scales. We account for the coupling of large stretching and bending strains of triangulated network models endowed with harmonic and dihedral angle potentials, on the basis of the discrete-continuum approach presented in Schmidt and Fraternali (J Mech Phys Solids 60:172-180, 2012). We test the proposed equilibrium fluctuation formulae with reference to a coarse-grained molecular dynamics model of the red blood cell (RBC) membrane (Marcelli et al. in Biophys J 89:2473-2480, 2005; Hale et al. in Soft Matter 5:3603-3606, 2009), employing a local maximum-entropy regularization of the fluctuating configurations (Fraternali et al. in J Comput Phys 231:528-540, 2012). We obtain information about membrane stiffening/softening due to stretching, curvature, and microscopic undulations of the RBC model. We detect local dependence of the elastic moduli over the RBC membrane, establishing comparisons between the present theory and different approaches available in the literature.