A spin-thermodynamic approach to characterize spin dynamics in TEMPO-based samples for dissolution DNP at 7 T field

The spin dynamics of dissolution DNP samples consisting of 4.5 M [¹³C]urea in a mixture of (1/1)_Vol glycerol/water using 4-Oxo-TEMPO as a radical was investigated. We analyzed the DNP dynamics as function of radical concentration at 7 T and 3.4 T static magnetic field as well as function of deuteration of the solvent matrix at the high field. The spin dynamics could be reproduced in all cases, at least qualitatively, by a thermodynamic model based on spin temperatures of the nuclear Zeeman baths and an electron non-Zeeman (dipolar) bath. We find, however, that at high field (7 T) and low radical concentrations (25 mM) the nuclear spins do not reach the same spin temperature indicating a weak coupling of the two baths. At higher radical concentrations, as well as for all radical concentrations at low field (3.4 T), the two nuclear Zeeman baths reach the same spin temperature within experimental errors. Additionally, the spin system was prepared with different initial conditions. For these cases, the thermodynamic model was able to predict the time evolution of the system well. While the DNP profiles do not give clear indications to a specific polarization transfer mechanism, at high field (7 T) increased coupling is seen. The EPR line shapes cannot clarify this in absence of ELDOR type experiments, nevertheless DNP profiles and dynamics under frequency-modulated microwave irradiation illustrate the expected increase in coupling between electrons with increasing radical concentration.

Direct hyperpolarization of micro- and nanodiamonds for bioimaging applications - Considerations on particle size, functionalization and polarization loss

Due to the inherently long relaxation time of ¹³C spins in diamond, the nuclear polarization enhancement obtained with dynamic nuclear polarization can be preserved for a time on the order of about one hour, opening up an opportunity to use diamonds as a new class of long-lived contrast agents. The present communication explores the feasibility of using ¹³C spins in directly hyperpolarized diamonds for MR imaging including considerations for potential in vivo applications.

Nanometer size silicon particles for hyperpolarized MRI

Hyperpolarized silicon particles have been shown to exhibit long spin-lattice relaxation times at room temperature, making them interesting as novel MRI probes. Demonstrations of hyperpolarized silicon particle imaging have focused on large micron size particles (average particle size (APS) = 2.2 μm) as they have, to date, demonstrated much larger polarizations than nanoparticles. We show that also much smaller silicon-29 particles (APS = 55 ± 12 nm) can be hyperpolarized with superior properties. A maximum polarization of 12.6% in the solid state is reported with a spin-lattice relaxation time of 42 min at room temperature thereby opening a new window for MRI applications.

Dissolution DNP using trityl radicals at 7 T field

Characterization of direct ¹³C DNP at 1.4 K and 7 T field using trityl radicals.

Conceptual and instrumental progress in dissolution DNP

We discuss conceptual and instrumental progress in dissolution DNP since its introduction in 2003. In our view there are three critical steps in the dissolution DNP process: (i) The achievable polarization level in a sample. (ii) The time required to build up the polarization. (iii) The transfer of the sample to the measurement system with minimum loss of polarization. In this review we describe in detail these steps and the different methodological and instrumental implementations, which have been proposed to optimize them.

Fast and slow orientational fluctuations in membranes

Fluorescence anisotropy measurements on isotropically distributed membranes yield the well-known orientational order parameter as a measure of orientational fluctuations of the fluorophore that are fast compared to the fluorescence lifetime. Measurements on oriented membranes provide four order parameters, two for fast orientational fluctuations and two for slow, from which their approximate angular distributions can be derived. This is exemplified for the fluctuations of lipids in membranes with and without protein. Steady-state fluorescence anisotropy experiments on oriented membranes are carried out, using diphenylhexatriene as a fluorescence probe for lipid fluctuations and melittin as a membrane protein. Protein molecules in a fluid lipid membrane restrict the fast lipid fluctuations and induce slow lipid fluctuations. The angular distribution of the slow fluctuations indicates that lipid molecules at the protein surface are tilted with respect to the membrane normal. The spatial coherence length of these fluctuations is of the order of 15 A.

Aggregation state of melittin in lipid vesicle membranes

We have performed time-resolved fluorescence energy transfer measurements using melittin as donor and a modified melittin as acceptor. The melittin molecules were bound to fluid vesicle membranes of dimyristoylphosphatidylcholine. Analysis of the temporal decay of the energy transfer and of its variation with the donor and acceptor concentrations led to the conclusion that melittin in fluid membranes is usually monomeric. Only at the high melittin/lipid molar ratio of 1/200 and high ionic strength evidence for aggregation was obtained, the percentage of aggregated melittin molecules being of the order of 10%. The shortcomings of previous steady-state measurements of fluorescence energy transfer between melittin molecules are discussed.

Enhanced internal dynamics of a membrane transport protein during substrate translocation

Conformational changes are essential for the activity of many proteins. If, or how fast, internal fluctuations are related to slow conformational changes that mediate protein function is not understood. In this study, we measure internal fluctuations of the transport protein lactose permease in the presence and absence of substrate by tryptophan fluorescence spectroscopy. We demonstrate that nanosecond fluctuations of alpha-helices are enhanced when the enzyme transports substrate. This correlates with previously published kinetic data from transport measurements showing that millisecond conformational transitions of the substrate-loaded carrier are faster than those in the absence of substrate. These findings corroborate the hypothesis of the hierarchical model of protein dynamics that predicts that slow conformational transitions are based on fast, thermally activated internal motions.

Effects of ligand binding on the internal dynamics of maltose-binding protein

Ligand binding to proteins often causes large conformational changes. A typical example is maltose-binding protein (MBP), a member of the family of periplasmic binding proteins of Gram-negative bacteria. Upon binding of maltose, MBP undergoes a large structural change that closes the binding cleft, i.e. the distance between its two domains decreases. In contrast, binding of the larger, nonphysiological ligand beta-cyclodextrin does not result in closure of the binding cleft. We have investigated the dynamic properties of MBP in its different states using time-resolved tryptophan fluorescence anisotropy. We found that the 'empty' protein exhibits strong internal fluctuations that almost vanish upon ligand binding. The measured relaxation times corresponding to internal fluctuations can be interpreted as originating from two types of motion: wobbling of tryptophan side-chains relative to the protein backbone, and orientational fluctuations of entire domains. After binding of a ligand, domain motions are no longer detectable and the fluctuations of some of the tryptophan side-chains become rather restricted. This transformation into a more rigid state is observed upon binding of both ligands, maltose and the larger beta-cyclodextrin. The fluctuations of tryptophan side-chains in direct contact with the ligand, however, are affected in a slightly different way by the two ligands.

Adhesion forces of lipids in a phospholipid membrane studied by molecular dynamics simulations

Lipid adhesion forces can be measured using several experimental techniques, but none of these techniques provide insight on the atomic level. Therefore, we performed extensive nonequilibrium molecular dynamics simulations of a phospholipid membrane in the liquid-crystalline phase out of which individual lipid molecules were pulled. In our method, as an idealization of the experimental setups, we have simply attached a harmonic spring to one of the lipid headgroup atoms. Upon retraction of the spring, the force needed to drag the lipid out of the membrane is recorded. By simulating different retraction rates, we were able to investigate the high pull rate part of the dynamical spectrum of lipid adhesion forces. We find that the adhesion force increases along the unbinding path, until the point of rupture is reached. The maximum value of the adhesion force, the rupture force, decreases as the pull rate becomes slower, and eventually enters a friction-dominated regime. The computed bond lengths depend on the rate of rupture, and show some scatter due to the nonequilibrium nature of the experiment. On average, the bond length increases from approximately 1.7 nm to 2.3 nm as the rates go down. Conformational analyses elucidate the detailed mechanism of lipid-membrane bond rupture. We present results of over 15 ns of membrane simulations. Implications for the interpretation and understanding of experimental rupture data are discussed.

Interaction between HLA-DM and HLA-DR involves regions that undergo conformational changes at lysosomal pH

Antigenic peptide loading of major histocompatibility complex class II molecules is enhanced by lysosomal pH and catalyzed by the HLA-DM molecule. The physical mechanism behind the catalytic activity of DM was investigated by using time-resolved fluorescence anisotropy (TRFA) and fluorescence binding studies with the dye 8-anilino-1-naphthalenesulfonic acid (ANS). We demonstrate that the conformations of both HLA-DM and HLA-DR3, irrespective of the composition of bound peptide, are pH sensitive. Both complexes reversibly expose more nonpolar regions upon protonation. Interaction of DM with DR shields these hydrophobic domains from the aqueous environment, leading to stabilization of the DM and DR conformations. At lysosomal pH, the uncovering of additional hydrophobic patches leads to a more extensive DM-DR association. We propose that DM catalyzes class II peptide loading by stabilizing the low-pH conformation of DR, favoring peptide exchange. The DM-DR association involves a larger hydrophobic surface area with DR/class II-associated invariant chain peptides (CLIP) than with stable DR/peptide complexes, explaining the preferred association of DM with the former. The data support a release mechanism of DM from the DM-DR complex through reduction of the interactive surface, upon binding of class II molecules with antigenic peptide or upon neutralization of the DM-DR complex at the cell surface.

Pathway of detergent-mediated and peptide ligand-mediated refolding of heterodimeric class II major histocompatibility complex (MHC) molecules

We investigated the mechanism of refolding and reassembly of recombinant alpha and beta chains of the class II major histocompatibility molecules (MHC-II) HLA-DRB5*0101. Both chains were expressed in the cytosol of Escherichia coli, purified in urea and SDS, and reassembled to functional heterodimers by replacement of SDS by mild detergents, incubation in a redox-shuffling buffer and finally by oxidation and removal of detergent. Refolding was mediated by mild detergents and by peptide ligands. Early stages of structure formation were characterized by circular dichroism, fluorescence, and time-resolved fluorescence anisotropy decay (FAD) spectroscopies. We found that formation of secondary structure was detectable after replacement of SDS by mild detergents. At that stage the alpha and beta chains were still monomeric, the buffer was strongly reducing, and the folding intermediates did not yet interact with peptide ligands. Formation of folding intermediates capable of interacting with peptide ligands was detected after adjusting the redox potential with oxidized glutathione and incubation in mild detergents. We conclude that at that stage a tertiary structure close to the native structure is formed at least locally. The nature and concentration of detergent critically determined the refolding efficiency. We compared detergents with different carbohydrate headgroups, and with aliphatic chains ranging from C6 to C14 in length. For each of the detergents we observed a narrow concentration range for mediating refolding. Surprisingly, detergents with long aliphatic chains had to be used at higher concentrations than short-chain detergents, indicating that increasing the solubility of folding intermediates is not the only function of detergents during a refolding reaction. We discuss structure formation and interactions of detergents with stable folding intermediates. Understanding such interactions will help to develop rational strategies for refolding hydrophobic or oligomeric proteins.

Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature

Molecular dynamics simulations of 500 ps were performed on a system consisting of a bilayer of 64 molecules of the lipid dipalmitoylphosphatidylcholine and 23 water molecules per lipid at an isotropic pressure of 1 atm and 50 degrees C. Special attention was devoted to reproduce the correct density of the lipid, because this quantity is known experimentally with a precision better than 1%. For this purpose, the Lennard-Jones parameters of the hydrocarbon chains were adjusted by simulating a system consisting of 128 pentadecane molecules and varying the Lennard-Jones parameters until the experimental density and heat of vaporization were obtained. With these parameters the lipid density resulted in perfect agreement with the experimental density. The orientational order parameter of the hydrocarbon chains agreed perfectly well with the experimental values, which, because of its correlation with the area per lipid, makes it possible to give a proper estimate of the area per lipid of 0.61 +/- 0.01 nm2.

Investigation of the proton release channel of bacteriorhodopsin in different intermediates of the photo cycle. A molecular dynamics study

Molecular dynamics simulations on bacteriorhodopsin were performed starting from a conformation based on electron cryomicroscopy studies [Henderson, R., et al. (1990) J. Mol. Biol. 213, 899-929]. We examined the proton release channel in different intermediates of the bacteriorhodopsin photocycle. In the simulations of the ground state, two stable sets of conformations were observed differing in the distance of the guanidinium group of Arg82 to the Schiff base. The set of conformations in which Arg82 is located closer to the Schiff base has a lower potential energy and agrees better with experimental data than the other set. With both sets, we performed a series of simulations in which the chromophore was isomerized to different states using purposive and nonpurposive methods. The energetic consideration of the different states argues for the location of the guanidinium group of Arg82 close to the Schiff base. The results also show that no C13-C14, C14-C15 dicis conformation of the retinal occurs in the K/L-intermediate of the photocycle instead supporting the occurrence of C13-C14 cis in these intermediates. In a last series of simulations, we modeled the M-intermediate of the bacteriorhodopsin photocycle. Again, comparison to different experimental data indicates that Arg82 points toward the Schiff base. We conclude that the guanidinium group of Arg82 is located close to the Schiff base at a distance of approximately 4.5 A and stays there at least up to the M-intermediate of the photocycle.

The use of a long-lifetime component of tryptophan to detect slow orientational fluctuations of proteins

The membrane protein porin and a synthetic polypeptide of 21 hydrophobic residues were inserted into detergent micelles or lipid membranes, and the fluorescence of their single tryptophan residue was measured in the time-resolved and polarized mode. In all cases, the tryptophan fluorescence exhibits a long-lifetime component of about 20 ns. This long-lifetime component was exploited to detect slow orientational motions in the range of tens of nanoseconds via the anisotropy decay. For this purpose, the analysis of the anisotropy has to be extended to account for different orientations of the dipoles of the short- and long-lifetime components. This is demonstrated for porin and the polypeptide solubilized in micelles, in which the longest relaxation time reflects the rotational diffusion of the micelle. When the polypeptide is inserted into lipid membranes, it forms a membrane-spanning alpha-helix, and the slowest relaxation process is interpreted as reflecting orientational fluctuations of the helix.

Modeling of halorhodopsin and rhodopsin based on bacteriorhodopsin

Bacteriorhodopsin (BR), halorhodopsin (HR), and rhodopsin (Rh) all belong to the class of seven-helix membrane proteins. For BR, a structural model at atomic resolution is available; for HR, diffraction data are available only down to a resolution of 6 A in the membrane plane, and for Rh, down to 9 A. BR and HR are closely related proteins with a sequence homology of 34%, while Rh does not share any sequence homology with BR. An atomic model for HR is derived that is based on sequence alignment and the atomic model for BR and is improved by molecular dynamics simulations. The model structure obtained accounts well for the experimentally observed difference between HR and BR in the projection map, where HR exhibits a higher density in the region between helices D and E. The reason for this difference lies partially in the different side chains and partially in slightly different helix tilts. The scattering amplitudes and phases of the model structure are calculated and agree with the experimental data down to a resolution of about 8 A. If the helix positions are adopted from the projection map for HR and used as input in the model, this number improves to 7 A. Analogously, an atomic model for Rh is derived based on the atomic model for BR and subjected to molecular dynamics simulations. Optimal agreement with the experimental projection map for Rh is obtained when the entire model structure is rotated slightly about two axes in the membrane plane. The agreement with the experimental projection map is not as satisfactory as for HR, but the results indicate that even for a nonhomologous, but structurally related, protein such as Rh, an acceptable model structure can be derived from the structure of BR.

Proton transport across transient single-file water pores in a lipid membrane studied by molecular dynamics simulations

To test the hypothesis that water pores in a lipid membrane mediate the proton transport, molecular dynamic simulations of a phospholipid membrane, in which the formation of a water pore is induced, are reported. The probability density of such a pore in the membrane was obtained from the free energy of formation of the pore, which was computed from the average force needed to constrain the pore in the membrane. It was found that the free energy of a single file of water molecules spanning the bilayer is 108(+/-10) kJ/mol. From unconstrained molecular dynamic simulations it was further deduced that the nature of the pore is very transient, with a mean lifetime of a few picoseconds. The orientations of water molecules within the pore were also studied, and the spontaneous translocation of a turning defect was observed. The combined data allowed a permeability coefficient for proton permeation across the membrane to be computed, assuming that a suitable orientation of the water molecules in the pore allows protons to permeate the membrane relatively fast by means of a wirelike conductance mechanism. The computed value fits the experimental data only if it is assumed that the entry of the proton into the pore is not rate limiting.

Folding and membrane insertion of the trimeric beta-barrel protein OmpF

We have studied folding and membrane insertion of the porin OmpF and compared it to OmpA. Both are beta-barrel membrane proteins from the outer membrane of Escherichia coli, OmpF forming trimers and OmpA monomers. Each of them can be unfolded in solubilized form in a water/urea mixture. Refolding is initiated by dilution into a dispersion of lipid vesicles or lipid/detergent vesicles, whereupon OmpF and OmpA refold and insert into the membranes. Folding and insertion of the monomers proceed in a similar way for the two proteins, but native OmpF appears more slowly and with a lower yield than native OmpA because of trimerization of OmpF. The dependence of the yield of refolding, membrane insertion, and trimerization on pH, lipid concentration, and the presence of detergent was investigated. Trimerization of OmpF is shown to take place at or in the membrane and a membrane-inserted dimer is detected as an intermediate of this process.

Kinetics of folding and membrane insertion of a beta-barrel membrane protein

We have studied the kinetics of folding and membrane insertion of the outer membrane protein OmpA of Escherichia coli. In the native structure, its membrane-inserted domain forms a beta-barrel. The protein was unfolded in solubilized form in water/urea, and refolding was induced by dilution of urea and simultaneous addition of lipid vesicles. Three transitions along the folding pathway could be distinguished. Their characteristic times lie below a second, in the range of minutes, and in the range of an hour. The fast process corresponds to the transition from the unfolded state in water/urea to a misfolded state in water, the moderately slow process to a transition from the misfolded state to a partially folded state in the membrane, and the slow process to the transition from the partially folded to the native state. The partially folded state in the membrane is interpreted as the analogue of the molten globule state of soluble proteins.

Lipid vesicle adsorption versus formation of planar bilayers on solid surfaces

The absorption and spreading behavior of lipid vesicles composed of either palmitoyloleoylphosphatidylcholine (POPC) or Escherichia coli lipid upon contact with a glass surface was examined by fluorescence measurements. Fluorescently labeled lipids were used to determine 1) the amount of lipid adsorbed at the surface, 2) the extent of fusion of the vesicles upon contact with the surface, 3) the ability of the adsorbed lipids to undergo lateral diffusion, and 4) the accessibility of the adsorbed lipids by external water soluble molecules. The results of these measurements indicate that POPC vesicles spread on the surface and form a supported planar bilayer, whereas E. coli lipid vesicles adsorb to the surface and form a supported vesicle layer. Supported planar bilayers were found to be permeable for small molecules, whereas supported vesicles were impermeable and thus represented immobilized, topologically separate compartments.

Structure and fluctuations of bacteriorhodopsin in the purple membrane: a molecular dynamics study

Molecular dynamics simulations on bacteriorhodopsin were performed starting from the model structure described by Henderson et al. The simulations were gradually improved by first treating a monomer in vacuum and then adding further monomers, lipids, and water to finally simulate a unit cell of the hexagonal lattice of the purple membrane containing a trimer and lipids and water on both sides. During all simulations, the protein structure moved away from the model structure to reach a root-mean-square (r.m.s.) deviation of 2 to 3 A. In the simulations with the trimer, the structures of the three monomers differed by about the same amount and averaging over them led to an average structure with a considerably smaller r.m.s. deviation. The best average structure obtained had an r.m.s. deviation from the model structure of 1.3 A. Fluctuations of the protein, the lipids, and water were analyzed in detail. As expected, the membrane-spanning helices of the protein fluctuate less than the peripheral loops. Unexpected, however, was the finding that the fluctuations of the protein are asymmetric with respect to the midplane of the membrane. The fluctuations of the loops and the ends of the helices on the inner side of the membrane are much stronger than on the outer side. This asymmetry is also reflected by the fluctuations for the lipids, the lipids of the inner leaflet fluctuating more strongly than those of the outer leaflet. The asymmetry was observed only in the presence of water on both sides of the membrane. On the average, nine water molecules were found inside the protein, most of them undergoing exchange with external water.

A nontransportable substrate for lactose permease

A substrate for lactose permease of Escherichia coli was synthesized that binds to the protein with a relatively high affinity, but is not transported to any detectable extent. This substrate, 6'-[(N-phenylalanylphenylalanyl)amino]hexyl 1-thio-beta-D-galactoside, is a peptide galactoside composed of a bulky aromatic dipeptide that is linked to galactose via an aminohexyl spacer. Binding of the peptide galactoside to lactose permease in cytoplasmic membranes was determined in a competition assay yielding a dissociation constant of 150 microM. Transport was measured by a counterflow assay using lipid vesicles with reconstituted lactose permease. An upper limit for the rate constant of transport was obtained as 0.02 s-1, 3 orders of magnitude smaller than the value for lactose.

Characterization of two membrane-bound forms of OmpA

The insertion of the outer membrane protein A (OmpA) into lipid bilayers was studied by limited proteolysis, polarized Fourier transform infrared (FTIR) spectroscopy, and fluorescence spectroscopy. In the native state, OmpA is thought to form a barrel of eight antiparallel beta-strands. For the present study, it was isolated in an unfolded form, purified, and exposed to performed vesicles of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC), dimyristoylphosphatidylcholine (DMPC), dipalmitoylphosphatidylcholine (DPPC), and three phospholipids that were brominated in different positions of their sn-2 chains (4,5-BrPC, 9,10-BrPC, and 11,12-BrPC). Limited proteolysis revealed two membrane-bound forms of OmpA, namely an "adsorbed" (35 kDa) and an "inserted" (30 kDa) form [Surrey, T., & Jähnig, F. (1992) Proc. Natl. Acad. Sci. U.S.A. 89, 7457-7461]. Which form was found after membrane binding and refolding depended on the lipids used and on the temperature. Polarized attenuated total reflection (ATR)-FTIR spectra were recorded with OmpA bound to germanium-supported bilayers in both forms. The position of the amide I' band indicated quite large fractions of beta-structure of OmpA in both membrane-bound forms (35-45% in the adsorbed form and 45-55% in the inserted form). Measurements of the linear dichroism of the amide I' bands in the inserted form are consistent with an antiparallel beta-barrel in which the strands are inclined at about 36 degrees from the membrane normal. The average angle of the beta-strands to the bilayer normal is likely larger in the 35 kDa form than in the inserted form.(ABSTRACT TRUNCATED AT 250 WORDS)

A long lifetime component in the tryptophan fluorescence of some proteins

The tryptophan fluorescence of two membrane proteins (outer membrane protein A and lactose permease), a 21-residue hydrophobic peptide, three soluble proteins (rat serum albumin, ribonuclease T1, and azurin), and N-acetyltryptophanamide (NATA) was investigated by time-resolved measurements extended over 65 ns. A long lifetime component with a characteristic time of 25 ns and an amplitude below 1% was found for outer membrane protein A, lactose permease, the peptide in lipid membranes, and azurin in water, but not for rat serum albumin, ribonuclease T1, and NATA in water. When outer membrane protein A was dissolved and unfolded in guanidinum hydrochloride, the long lifetime component disappeared. Hence, a hydrophobic environment seems to be a necessary requirement for the long lifetime component to be present. However, NATA dissolved in butanol does not exhibit the long lifetime component, while the peptide dissolved in the same solvent under conditions which preserve its helical structure does show the long lifetime. Thus, a regular secondary structure for the polypeptide chain to which the tryptophan residue belongs seems to be a second necessary requirement for the long lifetime component to be present. The long lifetime component may therefore be seen in the context of protein substates.

A synthetic analogue of melittin aggregates in large oligomers

An analogue of melittin synthesized in the group of E. T. Kaiser (DeGrado, W. F., F. J. Keźdy, and E. T. Kaiser. 1981. J. Am. Chem. Soc. 103:679-681) was investigated by Raman spectroscopy and fluorescence anisotropy decay. In water, the analogue is completely alpha-helical and aggregates in large oligomers of about 50 monomers. In vesicle membranes, it undergoes orientational fluctuations similar to melittin. The most significant difference from melittin, therefore, is the formation of straight helixes and their aggregation in large oligomers in water. We interpret this as a consequence of the lacking proline residue in the analogue. We, furthermore, hypothesize that the increased tendency for aggregation causes the increased hemolytic activity of the analogue.

Refolding and oriented insertion of a membrane protein into a lipid bilayer

We have studied the refolding and membrane insertion of the outer membrane protein OmpA of Escherichia coli. The protein was extracted from its native membrane by sonication in the presence of urea and dissolved in the urea/water mixture in unfolded form. In this form it was purified. Upon addition of preformed lipid vesicles, the protein spontaneously refolded and inserted into the vesicle membranes. The vesicles had to be small and the lipids had to be in the fluid state. The insertion occurred in an oriented manner.

Modeling of the structure of bacteriorhodopsin. A molecular dynamics study

Secondary structure predictions for membrane proteins are relatively reliable and permit the construction of model structures that may serve as initial conformations for molecular dynamics simulations. This might provide a scheme to predict the three-dimensional structures of membrane proteins. The feasibility of such an approach is tested for bacteriorhodopsin. We were not able to fully predict the kidney-shaped structure of bacteriorhodopsin. However, features compatible with this structure developed in a simulation starting from a circular arrangement of the seven predicted helices. When instead we started from the kidney shape, assigning the seven predicted helices in different ways to those on the structure, we could distinguish between the different assignments on the basis of energy and tilt of the helices. In this way we could select the correct assignment from a few others. For the correct assignment, the helices spontaneously adopted a tilt that agrees remarkably well with the experimental model structure derived by others. The root-mean-square deviation between our best molecular dynamics structure and the experimental model structure is 3.8 A, caused mainly by deviations in the internal degrees of freedom of the helices.

Conformational transition of an alpha-helix studied by molecular dynamics

Molecular dynamics simulations were performed on a 20-residue polyalanine helix and a spontaneous transition from a kinked to a straight conformation was observed. The kinetics of the transition was analyzed within the framework of the Kramers model for chemical reactions and within a random walk model. The Kramers model which is based on diffusion along a one-dimensional reaction pathway and the crossing of an energy barrier was found to be inadequate. Instead, a random walk model based on diffusion in the high-dimensional phase space of the system was found to be compatible with the data. The high dimensionality of the phase space permits the system to circumvent high energy barriers and diffuse rapidly at about constant energy, but decelerates the reaction since in the labyrinth of pathways the transition state is reached rarely.

The opacity proteins of Neisseria gonorrhoeae strain MS11 are encoded by a family of 11 complete genes

Variants of Neisseria gonorrhoeae MS11 show distinct colony morphologies because of the expression of a class of surface components called opacity (Opa, PII) proteins. Southern analyses combined with molecular cloning of genomic DNA from a single variant of MS11 has identified 11 opa genes contained in separate loci. These opa genes code for distinct opacity proteins which are distinguishable at their variable domains. The opa gene analyses were also extended to divergent variants of MS11. These studies have shown that, during in vitro and in vivo culture, 10 of the 11 opa genes did not undergo significant change in their primary sequence. However, in these variants, one gene (opaE) underwent non-reciprocal inter-opa recombinations to generate newer Opa variants. Phylogenic analysis of the opa gene sequences suggests that the opa gene family have evolved by a combination of gene duplication, gene replacement and partial inter-opa recombination events.

Collective vibrations of an alpha-helix. A molecular dynamics study

The internal dynamics of a 20-residue polyalanine helix was investigated by molecular dynamics simulations. Special attention was paid to the collective vibrations of the helix backbone. The stretch and bend vibrations could be assigned unambiguously to oscillations with periods of 1.4 and 4.3 ps, respectively. The influence of the environment on the dynamics of the collective vibrations was studied by coupling the helix to a heat bath and by adding water molecules. In the presence of water, the stretch vibration becomes more strongly dampened, but still exists as a vibration , while the bend vibration becomes overdamped and degenerates into a relaxation process. The results are compared with available experimental data.

Biosensors based on membrane transport proteins

We propose a novel class of biosensors based on membrane bound receptors or transport proteins as the sensing element. The protein is incorporated in a planar lipid bilayer which covers the transducer. The transducer may detect an electric current, a voltage, or a change in fluorescence. A prototype lactose sensor is presented which consists of a quartz slide covered by a lipid membrane containing the protein lactose permease from Escherichia coli. This protein is a lactose/H+ cotransporter, hence lactose in the external medium initiates lactose/H+ cotransport across the lipid membrane. This leads to a rise in proton concentration in the small volume between the lipid membrane and the quartz surface which can be detected by a pH-sensitive fluorescence dye.

Refolding of an integral membrane protein. OmpA of Escherichia coli

OmpA is an integral membrane protein from the outer membrane of Escherichia coli. Purified, lipopolysaccharide-free OmpA was denatured by boiling in sodium dodecyl sulfate (SDS). Refolding was then induced by replacement of SDS with the nonionic detergent octylglucoside. The structure of both the denatured and refolded protein were investigated by SDS-gel electrophoresis, protease digestion, Raman and fluorescence spectroscopy. Refolded OmpA could be reconstituted into membranes of the synthetic lipid dimyristoylphosphatidylcholine. Thus, lipopolysaccharide is neither necessary for proper folding of OmpA nor for its insertion into lipid membranes. Based on this result, models for sorting of OmpA into the outer membrane of E. coli are discussed.

Structure predictions of membrane proteins are not that bad

A simple generalization of the well-known hydrophobicity analysis is sufficient to predict membrane-spanning amphiphilic alpha-helices and beta-strands. The use of this method is illustrated for bacteriorhodopsin, OmpA protein and the photosynthetic reaction centre.

Restoration of membrane incorporation of an Escherichia coli outer membrane protein (OmpA) defective in membrane insertion

The mechanism of sorting, to the outer membrane, of the 325-residue Escherichia coli protein OmpA has been investigated. It is thought to traverse the membrane eight times in antiparallel beta-strands, forming an amphiphilic beta-barrel which encompasses residues 1 to about 170; the COOH-terminal moiety is periplasmic. A mutant, carrying the substitutions Leu164----Pro and Val166----Asp within the last beta-strand (residues 160-170), has been described which was unable to assemble in the membrane (Klose, M., MacIntyre, S., Schwarz, H., and Henning, U. (1988) J. Biol. Chem. 263, 13297-13302). Linkers were inserted between the codons for residues 164 and 165 of the mutant protein. Of 13 different genes recovered, five encoded proteins which had regained the ability to assemble in the membrane. The properties of the mutant proteins, together with a structure prediction method, indicate the following rules for the final beta-strand to be compatible with, or possibly initiate, membrane insertion: (i) it must be amphiphilic or hydrophobic while its primary structure as such is fairly unimportant, (ii) it must extend over at least 9 residues, and (iii) it must not contain a proline residue around its center. One of the genes recovered coded for OmpA up to residue 164 and then followed by 10 linker-encoded residues. This 174-residue polypeptide was assembled in the membrane but did not, in contrast to all other proteins, expose sites sensitive to trypsin at the inner face of the membrane. This behavior agrees perfectly well with the OmpA model.

Internal dynamics of lactose permease

The transport protein lactose permease was reconstituted in vesicles of dimyristoylphosphatidylcholine, and the internal dynamics were studied by measuring the fluorescence anisotropy decay of the tryptophan residues and of a covalently bound pyrene label. For the tryptophans three relaxation processes and for the pyrene two relaxation processes with relaxation times in the nanosecond range were observed. The slowest process, of approximately 50 ns, is assigned to orientational fluctuations of membrane-spanning helices. When the temperature is decreased below the lipid-phase transition, this relaxation process is slowed down and restricted in amplitude. Because the transport rate is known to also decrease below the phase transition, this observation suggests a coupling between internal dynamics and transport. This coupling is analyzed on the basis of the Kramers relation for chemical reactions.

Fast measurement of galactoside transport by lactose permease

Lactose permease of Escherichia coli was reconstituted into vesicles of dimyristoylphosphatidylcholine, and the rate of galactoside counterflow was measured in the millisecond time range. The turnover number and the half-saturation constant for transport agree with the values known for cells. This result demonstrates that lactose permease is the sole protein necessary for galactoside transport. Furthermore, lactose permease seems not to require a high level of negatively charged lipids or a certain degree of unsaturation of the lipid hydrocarbon chains. However, the lipids must be in the fluid state, because the transport rate drastically decreases below the lipid ordered fluid phase transition.

Secondary structure of the promastigote surface protease of Leishmania

By Raman spectroscopic analysis we have determined the secondary structure of the promastigote surface protease, named PSP or gp63, of Leishmania major. It consist of nearly 50% antiparallel beta-strand, and less than 20% alpha-helix. These results are contrasted with the predominantly alpha-helical VSGs of the African trypanosomes and the alpha-helical metalloprotease thermolysin. The PSP of Leishmania thus represents a novel kind of membrane-anchored protease.

Dynamics of melittin in water and membranes as determined by fluorescence anisotropy decay

Fluorescence anisotropy decay measurements were performed on melittin in water and in membranes of dimyristoylphosphatidylcholine. The fluorescence of the single tryptophan residue of melittin and of a pyrene label attached to melittin was detected. In water, the slowest relaxation process in the anisotropy decay occurs with a relaxation time of 1.5 or 5.5 ns in the case of low or high ionic strength and corresponds to rotational diffusion of monomeric or tetrameric melittin. Superimposed on this slow process are fast processes in the subnanosecond range reflecting fluctuations of the fluorophores relative to the polypeptide backbone. In membranes, the fast relaxation processes are not much altered. A slow process with a relaxation time of 35 ns is observed and assigned to orientational fluctuations of the melittin helices in membranes.

The structure of a membrane-spanning polypeptide studied by molecular dynamics

We have performed a molecular dynamics simulation of a 46-residue segment of glycophorin which includes the hydrophobic membrane-spanning region of this protein. The presence of a membrane and of water is taken into account in a continuum approximation which makes use of phenomenological hydrophobic energies. The initial alpha-helical conformation and the membrane incorporation of the hydrophobic segment remain stable for the length of the simulation which is 100 ps. Moreover, when the hydrophobic segment is partially shifted out of the membrane, it moves back into the membrane. Superimposed on these deterministic effects one also observes thermal fluctuations in the form of bending and tilting of the membrane-spanning helix.