Solution structure of the ribosome recycling factor from Aquifex aeolicus

The solution structure of ribosome recycling factor (RRF) from hyperthermophilic bacterium, Aquifex aeolicus, was determined by heteronuclear multidimensional NMR spectroscopy. Fifteen structures were calculated using restraints derived from NOE, J-coupling, and T1/T2 anisotropies. The resulting structure has an overall L-shaped conformation with two domains and is similar to that of a tRNA molecule. The domain I (corresponding to the anticodon stem of tRNA) is a rigid three alpha-helix bundle. Being slightly different from usual coiled-coil arrangements, each helix of domain I is not twisted but straight and parallel to the main axis. The domain II (corresponding to the portion with the CCA end of tRNA) is an alpha/beta domain with an alpha-helix and two beta-sheets, that has some flexible regions. The backbone atomic root-mean-square deviation (rmsd) values of both domains were 0.7 A when calculated separately, which is smaller than that of the molecule as a whole (1.4 A). Measurement of 15N-[1H] NOE values show that the residues in the corner of the L-shaped molecule are undergoing fast internal motion. These results indicate that the joint region between two domains contributes to the fluctuation in the orientation of two domains. Thus, it was shown that RRF remains the tRNA mimicry in solution where it functions.     

Structure and binding mode of a ribosome recycling factor (RRF) from mesophilic bacterium

X-ray and NMR analyses on ribosome recycling factors (RRFs) from thermophilic bacteria showed that they display a tRNA-like L-shaped conformation consisting of two domains. Since then, it has been accepted that domain I, consisting of a three-helix bundle, corresponds to the anticodon arm of tRNA and domain II and a beta/alpha/beta sandwich structure, corresponds to the acceptor arm. In this study, we obtained a RRF from a mesophilic bacterium, Vibrio parahaemolyticus, by gene cloning and carried out an x-ray analysis on it at 2.2 A resolution. This RRF was shown to be active in an in vitro assay system using Escherichia coli polysomes and elongation factor G (EF-G). In contrast, the above-mentioned RRFs from thermophilic bacteria were inactive in such a system. Analysis of the relative orientations between the two domains in the structures of various RRFs, including this RRF from mesophilic bacterium, revealed that domain II rotates about the long axis of the helix bundle of domain I. To elucidate the ribosome binding site of RRF, the peptide fragment (RRF-DI) corresponding to domain I of RRF was expressed and characterized. RRF-DI is bound to 70 S ribosome and the 50 S subunit with an affinity similar to that of wild-type RRF. But it does not bind to the 30 S subunit. These findings caused us to reinvestigate the concept of the mimicry of RRF to tRNA and to propose a new model where domain I corresponds to the acceptor arm of tRNA and domain II corresponds to the anticodon arm. This is just the reverse of a model that is now widely accepted. However, the new model is in better agreement with published biological findings.     

Backbone dynamics of inactive, active, and effector-bound Cdc42Hs from measurements of (15)N relaxation parameters at multiple field strengths.

Cdc42Hs, a member of the Ras superfamily of GTP-binding proteins, initiates a cascade that begins with the activation of several kinases, including p21-activated kinase (PAK). We have previously determined the structure of Cdc42Hs and found that the regions involved in effector (Switch I) and regulator (Switch II) actions are partially disordered [Feltham, J. L., et al. (1997) Biochemistry 36, 8755-8766]. Recently, we used a 46-amino acid fragment of PAK (PBD46) to define the binding surface on Cdc42Hs, which includes the beta2 strand and a portion of Switch I [Guo, W., et al. (1998) Biochemistry 37, 14030-14037]. Here we describe the backbone dynamics of three constructs of [(15)N]Cdc42Hs (GDP-, GMPPCP-, and GMPPCP- and PBD46-bound) using (15)N-(1)H NMR measurements of T(1), T(1)(rho), and the steady-state NOE at three magnetic field strengths. Residue-specific values of the generalized order parameters (S(s)(2) and S(f)(2)), local correlation time (tau(e)), and exchange rate (R(ex)) were obtained using the Lipari-Szabo model-free formalism. Residues in Switch I were found to exhibit high-amplitude (low-order) motions on a nanosecond time scale, whereas those in Switch II experience low-amplitude motion on the nanosecond time scale and chemical (conformational) exchange on a millisecond time scale. The Insert region of Cdc42Hs-GDP exhibits high-order, nanosecond motions; the time scale of motion in the Insert is reduced in Cdc42Hs-GMPPCP and Cdc42Hs-PBD46. Overall, significant flexibility was observed mainly in the regions of Cdc42Hs that are involved in protein-protein interactions (Switch I, Switch II, and Insert), and flexibility was reduced upon interaction with a protein ligand. These results suggest that protein flexibility is important for high-affinity binding interactions.     

Hsc70-interacting HPD loop of the J domain of polyomavirus T antigens fluctuates in ps to ns and micros to ms

The backbone dynamics of the J domain from polyomavirus T antigens have been investigated using 15N NMR relaxation and molecular dynamics simulation. Model-free relaxation analysis revealed picosecond to nanosecond motions in the N terminus, the I-II loop, the C-terminal end of helix II through the HPD loop to the beginning of helix III, and the C-terminal end of helix III to the C terminus. The backbone dynamics of the HPD loop and termini are dominated by motions with moderately large amplitudes and correlation times of the order of a nanosecond or longer. Conformational exchange on the microsecond to millisecond timescale was identified in the HPD loop, the N and C termini, and the I-II loop. A 9.7ns MD trajectory manifested concerted swings of the HPD loop. Transitions between major and minor conformations of the HPD loop featured distinct patterns of change in backbone dihedral angles and hydrogen bonds. Fraying of the C-terminal end of helix II and the N-terminal end of helix III correlated with displacements of the HPD loop. Correlation of crankshaft motions of Gly46 and Gly47 with the collective motions of the HPD loop suggested an important role of the two glycine residues in the mobility of the loop. Fluctuations of the HPD loop correlated with relative reorientation of side-chains of Lys35 and Asp44 that interact with Hsc70.     

Modulating LOV domain photodynamics with a residue alteration outside the chromophore binding site

Phototropins, a class of light-activated protein kinases, are essential for several blue light responses in plants and algae, including phototropism. These proteins contain two internal light, oxygen, and voltage sensitive (LOV) domains, which bind flavin chromophores and undergo a reversible photochemical formation of a cysteinyl-flavin adduct as part of the light sensing process. While the photodynamic properties of such photosensory domains are dictated by interactions between the chromophore and surrounding protein, more distant residues can play a significant role as well. Here we explore the role of the Phe434 residue in the photosensory response of the second LOV domain of Avena sativa phototropin 1 (AsLOV2), a model photochemical system for these LOV domains. Phe434 is more than 6 ? from the FMN chromophore in AsLOV2; nevertheless, an F434Y point mutation is likely to change several structural features of the chromophore binding site, as we demonstrate using molecular dynamics simulations. Transient absorption signals spanning 15 decades in time were compared for wild-type AsLOV2 and the F434Y mutant, showing that the latter has significantly altered photodynamics, including (i) a faster intersystem crossing leading to triplet formation on a nanosecond time scale, (ii) biphasic formation of adduct-state kinetics on the microsecond time scale, and (iii) greatly accelerated ground-state recovery kinetics on a second time scale. We present mechanistic models that link these spectroscopic differences to changes in the configuration of the critical cysteine residue and in the chromophore's accessibility to solvent and oxygen according to MD trajectories and purging experiments. Taken together, these results demonstrate the importance of residues outside the chromophore-binding pocket in modulating LOV domain photodynamics.     

Dynamics of Ca2+-saturated calmodulin D129N mutant studied by multiple molecular dynamics simulations

Fifteen independent 1-nsec MD simulations of fully solvated Ca(2+) saturated calmodulin (CaM) mutant D129N were performed from different initial conditions to provide a sufficient statistical basis to gauge the significance of observed dynamical properties. In all MD simulations the four Ca(2+) ions remained in their binding sites, and retained a single water ligand as observed in the crystal structure. The coordination of Ca(2+) ions in EF-hands I, II, and III was sevenfold. In EF-hand IV, which was perturbed by the mutation of a highly conserved Asp129, an anomalous eightfold Ca(2+) coordination was observed. The Ca(2+) binding loop in EF-hand II was observed to dynamically sample conformations related to the Ca(2+)-free form. Repeated MD simulations implicate two well-defined conformations of Ca(2+) binding loop II, whereas similar effect was not observed for loops I, III, and IV. In 8 out of 15 MD simulations Ca(2+) binding loop II adopted an alternative conformation in which the Thr62 >C=O group was displaced from the Ca(2+) coordination by a water molecule, resulting in the Ca(2+) ion ligated by two water molecules. The alternative conformation of the Ca(2+) binding loop II appears related to the "closed" state involved in conformational exchange previously detected by NMR in the N-terminal domain fragment of CaM and the C-terminal domain fragment of the mutant E140Q. MD simulations suggest that conformations involved in microsecond exchange exist partially preformed on the nanosecond time scale.     

Mechanistic investigation of domain specific unfolding of human serum albumin and the effect of sucrose

This study is devoted to understand the unfolding mechanism of a multidomain protein, human serum albumin (HSA), in absence and presence of the sucrose by steady‐state and time‐resolved fluorescence spectroscopy with domain specific marker molecules and is further being substantiated by molecular dynamics (MD) simulation. In water, the domain III of HSA found to unfold first followed by domains I and II as the concentration of GnHCl is increased in the medium. The sequential unfolding behavior of different domains of HSA remains same in presence of sucrose; however, a higher GnHCl concentration is required for unfolding, suggesting stabilizing effect of sucrose on HSA. Domain I is found to be most stabilized by sucrose. The stabilization of domain II is somewhat similar to domain I, but the effect of sucrose on domain III is found to be very small. MD simulation also predicted a similar behavior of sucrose on HSA. The stabilizing effect of sucrose is explained in terms of the entrapment of water molecules in between HSA surface and sucrose layer as well as direct interaction between HSA and sucrose.     

Exploring the flexibility of ribosome recycling factor using molecular dynamics

Ribosome recycling factor is proposed to be flexible, and that flexibility is believed to be important to its function. Here we use molecular dynamics to test the flexibility of Escherichia coli RRF (ecRRF) with and without decanoic acid bound to a hydrophobic pocket between domains 1 and 2, and Thermus thermophilus RRF (ttRRF) with and without a mutation in the hinge between domains 1 and 2. Our simulations show that the structure of ecRRF rapidly goes from having an interdomain angle of 124 degrees to an angle of 98 degrees independently of the presence of decanoic acid. The simulations also show that the presence or absence of decanoic acid leads to changes in ecRRF flexibility. Simulations of wild-type and mutant ttRRF (R32G) show that mutating Arg-32 to glycine decreases RRF flexibility. This was unexpected because the range of dihedral angles for arginine is limited relative to glycine. Furthermore, the interdomain angle of wild-type T. thermophilus goes from 81 degrees to 118 degrees whereas the R32G mutant remains very close to the crystallographic angle of 78 degrees . We propose that this difference accounts for the fact that mutant ttRRF complements an RRF deficient strain of E. coli whereas wild-type ttRRF does not. When the ensemble of RRF structures is modeled into the ribosomal crystal structure, a series of overlaps is found that corresponds with regions where conformational changes have been found in the cryoelectron microscopic structure of the RRF/ribosome complex, and in the crystal structure of a cocomplex of RRF with the 50S subunit. There are also overlaps with the P-site, suggesting that RRF flexibility plays a role in removing the deacylated P-site tRNA during termination of translation.     

Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations. Proteins 31:116–127, 1998. © 1998 Wiley-Liss, Inc.     

(1)H/(15)N heteronuclear NMR spectroscopy shows four dynamic domains for phospholamban reconstituted in dodecylphosphocholine micelles

We report the backbone dynamics of monomeric phospholamban in dodecylphosphocholine micelles using (1)H/(15)N heteronuclear NMR spectroscopy. Phospholamban is a 52-amino acid membrane protein that regulates Ca-ATPase in cardiac muscle. Phospholamban comprises three structural domains: a transmembrane domain from residues 22 to 52, a connecting loop from 17 to 21, and a cytoplasmic domain from 1 to 16 that is organized in an "L"-shaped structure where the transmembrane and the cytoplasmic domain form an angle of approximately 80 degrees (Zamoon et al., 2003; Mascioni et al., 2002). T(1), T(2), and (1)H/(15)N nuclear Overhauser effect values measured for the amide backbone resonances were interpreted using the model-free approach of Lipari and Szabo. The results point to the existence of four dynamic domains, revealing the overall plasticity of the cytoplasmic helix, the flexible loop, and part of the transmembrane domain (residues 22-30). In addition, using Carr-Purcell-Meiboom-Gill-based experiments, we have characterized phospholamban dynamics in the micros-ms timescale. We found that the majority of the residues in the cytoplasmic domain, the flexible loop, and the first ten residues of the transmembrane domain undergo dynamics in the micros-ms range, whereas minimal dynamics were detected for the transmembrane domain. Hydrogen/deuterium exchange factors measured at different temperatures support the existence of slow motion in both the loop and the cytoplasmic helix. We propose that these dynamic properties are critical factors in the biomolecular recognition of phospholamban by Ca-ATPase and other interacting proteins such as protein kinase A and protein phosphatase 1.     

NMR backbone dynamics of the human type I interferon binding subunit, a representative cytokine receptor

The antiviral and antiproliferative activities of type I interferons (IFNs) are mediated by a common receptor, and its second subunit (IFNAR2) exhibits nanomolar affinity to both IFNalpha and IFNbeta subtypes. We have previously determined the structure of the IFN-binding extracellular domain of IFNAR2 (IFNAR2-EC) using multidimensional NMR [Chill, J. H., Quadt, S. R., Levy, R., Schreiber, G. E., and Anglister, J. (2003) Structure 11, 791-802], showing it to comprise two fibronectin domains linked by a hinge. As the first cytokine receptor structure determined in the unliganded state and in solution, IFNAR2-EC offers an opportunity to characterize the dynamics of the cytokine receptor family and their correlation to biological function. Backbone dynamics of IFNAR2-EC were investigated using 15N relaxation at 11.74 and 18.79 T, and measurements of residual dipolar couplings (RDCs). Dynamics of the binding site distinguish between rigid structural domains, which stabilize the binding site conformation, and a more flexible binding interface which interacts with the ligand. Measurements of diffusional anisotropy and RDCs and model-free analysis all show that the backbone of the hinge interdomain region of IFNAR2-EC is rigid on the picosecond to nanosecond time scale. Signal transduction in cytokines receptors is initiated by ligand-induced juxtaposition of the two receptor subunits, triggering the mutual phosphorylation of kinases associated to their cytoplasmic domains. The rigidity of the hinge ensures correct positioning of the receptor subunits in the ternary signaling complex and modulates the interaction between kinases in the cytoplasm, thereby controlling the rate and efficiency of phosphorylation.     

Comparison of multiple molecular dynamics trajectories calculated for the drug-resistant HIV-1 integrase T66I/M154I catalytic domain

HIV-1 integrase (IN) is an essential enzyme for the viral replication and an interesting target for the design of new pharmaceuticals for multidrug therapy of AIDS. Single and multiple mutations of IN at residues T66, S153, or M154 confer degrees of resistance to several inhibitors that prevent the enzyme from performing its normal strand transfer activity. Four different conformations of IN were chosen from a prior molecular dynamics (MD) simulation on the modeled IN T66I/M154I catalytic core domain as starting points for additional MD studies. The aim of this article is to understand the dynamic features that may play roles in the catalytic activity of the double mutant enzyme in the absence of any inhibitor. Moreover, we want to verify the influence of using different starting points on the MD trajectories and associated dynamical properties. By comparison of the trajectories obtained from these MD simulations we have demonstrated that the starting point does not affect the conformational space explored by this protein and that the time of the simulation is long enough to achieve convergence for this system.     

Domain motions of the Mip protein from Legionella pneumophila

The homodimeric 45.6 kDa (total mass) Mip protein, a virulence factor from Legionella pneumophila, was investigated with solution NMR spectroscopy and molecular dynamics (MD) simulations. Two Mip monomers are dimerized via an N-terminal helix bundle that is connected via a long alpha-helix to a C-terminal FKBP domain in each subunit. More than 85% of the amino acids were identified in triple-resonance NMR spectra. (15)N relaxation analysis showed a bimodal distribution of R(1)/R(2) values, with the lower ratio in the N-terminal domain. Relaxation dispersion measurements confirmed that these reduced ratios did not originate from conformational exchange. Thus, two different correlation times (tau(c)) can be deduced, reflecting partly uncoupled motions of both domains. Relaxation data of a Mip(77)(-)(213) monomer mutant were similar to those observed in the dimer, corroborating that the FKBP domain, including part of the connecting helix, behaves as one dynamic entity. MD simulations (18 ns) of the Mip dimer also yielded two different correlation times for the two domains and thus confirm the independence of the domain motions. Principal component analysis of the dihedral space covariance matrix calculated from the MD trajectory suggests a flexible region in the long connecting helix that acts as a hinge between the two domains. Such motion provides a possible explanation of how Mip can bind to complex molecular components of the extracellular matrix and mediate alveolar damage and bacterial spread in the lung.     

NMR structure and dynamics of the second transmembrane domain of the neuronal acetylcholine receptor beta 2 subunit

Structure and backbone dynamics of a selectively [(15)N]Leu-labeled 28-residue segment of the extended second transmembrane domain (TM2e) of the human neuronal nicotinic acetylcholine receptor (nAChR) beta(2) subunit were studied by (1)H and (15)N solution-state NMR in dodecylphosphocholine micelles. The TM2e structure was determined on the basis of the nuclear Overhauser effects (NOEs) and the hydrogen bond restraints, which were inferred from the presence of H(alpha)(i)-H(N)(i+3), H(alpha)(i)-H(beta)(i+3), and H(alpha)(i)-H(N)(i+4) NOE connectivity and from the slow amide hydrogen exchange with D(2)O. The TM2e structure of the nAChR beta(2) subunit contains a helical region between T4 and K22. Backbone dynamics were calculated using the model-free approach based on the (15)N relaxation rate constants, R(1) and R(2), and on the (15)N-[(1)H] NOE. The data acquired at 9.4 and 14.1 T and calculations using different dynamic models demonstrated no conformational exchange and internal motions on the nanosecond time scale. The global tumbling time of TM2e in micelles was 14.4 +/- 0.2 ns; the NOE values were greater than 0.63 at 9.4 T, and the order parameter, S(2), was 0.83-0.96 for all (15)N-labeled leucine residues, suggesting a restricted internal motion. This is the first report of NMR structure and backbone dynamics of the second transmembrane domain of the human nAChR beta(2) subunit in a membrane-mimetic environment, providing the basis for subsequent studies of subunit interactions in the transmembrane domain complex of the neuronal nAChR.     

Polymorphism of the epidermal growth factor receptor extracellular ligand binding domain: the dimer interface depends on domain stabilization

Epidermal growth factor receptors (EGFRs) and their cytoplasmic tyrosine kinases play important roles in cell proliferation and signaling. The EGFR extracellular domain (sEGFR) forms a dimer upon the binding of ligands, such as epidermal growth factor (EGF) and transforming growth factor α (TGFα). In this study, multiple molecular dynamics (MD) simulations of the 2:2 EGF·sEGFR3-512 dimer and the 2:2 TGFα·sEGFR3-512 dimer were performed in solvent and crystal environments. The simulations of systems comprising up to half a million atoms reveal part of the structural dynamics of which sEGFR dimers are capable. The solvent simulations consistently exhibited a prominent conformational relaxation from the initial crystal structures on the nanosecond time scale, leading to symmetry breaking and more extensive contacts between the two sEGFR monomers. In the crystal control simulation, this symmetry breaking and compaction was largely suppressed by crystal packing contacts. The simulations also provided evidence that the disordered domain IV of sEGFR may act as a stabilizing spacer in the dimer. Thus, the simulations suggest that the sEGFR dimer can take diverse configurations in solvent environments. These biologically relevant conformations of the EGFR signal transduction network can be controlled by contacts among the structural domains of sEGFR and its ligands.     

Expression in Escherichia coli, purification and characterization of Thermoanaerobacter tengcongensis ribosome recycling factor

A very promising approach to understanding the mechanism of protein thermostability is to investigate the structure-function relationship of homologous proteins with different thermostabilities. Ribosome recycling factor (RRF), which is an essential factor for protein synthesis in bacteria, may be a good candidate for such study. In this report, a ribosome recycling factor from Thermoanaerobacter tengcongensis was expressed and characterized. This protein contains 184 residues, shows 51.4% identity to that of Escherichia coli RRF, and has very strong antigenic cross-reactivity with antibody to E. coli RRF. In vivo activity assay shows that weak residual activity may remain in TteRRF in E. coli cells. Circular dichroism spectral analysis shows that TteRRF has a very similar secondary structure to that of E. coli RRF, implying that they have similar tertiary structures. However, their thermostabilities are significantly different. To find which domain of RRF is mainly responsible for maintaining stability, TteDI/EcoDII and EcoDI/TteDII RRF chimeras were created. Their domain I and domain II are from E. coli and T. tengcongensis RRFs, respectively. The results of GdnHCl and heat induced denaturation of the chimeric RRFs suggest that the domain I plays a major role in maintaining the stability of the RRF molecule.     

Molecular dynamics simulations as a tool for improving protein stability

Haloalkane dehalogenase (DhlA) was used as a model protein to explore the possibility to use molecular dynamics (MD) simulations as a tool to identify flexible regions in proteins that can serve as a target for stability enhancement by introduction of a disulfide bond. DhlA consists of two domains: an alpha/beta-hydrolase fold main domain and a cap domain composed of five alpha-helices. MD simulations of DhlA showed high mobility in a helix-loop-helix region in the cap domain, involving residues 184-211. A disulfide cross-link was engineered between residue 201 of this flexible region and residue 16 of the main domain. The mutant enzyme showed substantial changes in both thermal and urea denaturation. The oxidized form of the mutant enzyme showed an increase of the apparent transition temperature from 47.5 to 52.5 degrees C, whereas the T(m,app) of the reduced mutant decreased by more than 8 degrees C compared to the wild-type enzyme. Urea denaturation results showed a similar trend. Measurement of the kinetic stability showed that the introduction of the disulfide bond caused a decrease in activation free energy of unfolding of 0.43 kcal mol(-1) compared to the wild-type enzyme and also indicated that the helix-loop-helix region was involved early in the unfolding process. The results show that MD simulations are capable of identifying mobile protein domains that can successfully be used as a target for stability enhancement by the introduction of a disulfide cross-link.     

Solid-state NMR investigation of the dynamics of the soluble and membrane-bound colicin Ia channel-forming domain.

Solid-state NMR spectroscopy was employed to study the molecular dynamics of the colicin Ia channel domain in the soluble and membrane-bound states. In the soluble state, the protein executes small-amplitude librations (with root-mean-square angular fluctuations of 0-10 degrees ) in the backbone and larger-amplitude motions (16-17 degrees ) in the side chains. Upon membrane binding, the motional amplitudes increase significantly for both the backbone (12-16 degrees ) and side chains (23-29 degrees ), as manifested by the reduction in the C-H and H-H dipolar couplings and (15)N chemical shift anisotropy. These motions occur not only on the pico- to nanosecond time scales, but also on the microsecond time scale, as revealed by the (1)H rotating-frame spin-lattice relaxation times. Average motional correlation times of 0.8 and 1.2 micros were extracted for the soluble and membrane-bound states, respectively. In comparison, both forms of the colicin Ia channel domain are completely immobile on the millisecond scale. These results indicate that the colicin Ia channel domain has enhanced conformational mobility in the lipid bilayer compared to the soluble state. This membrane-induced mobility increase is consistent with the loss of tertiary structure of the protein in the membrane, which was previously suggested by the extended helical array model [Zakharov et al. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 4282-4287]. An extended structure would also facilitate protein interactions with the mobile lipids and thus increase the protein internal motions. We speculate that the large mobility of the membrane-bound colicin Ia channel domain is a prerequisite for channel opening in the presence of a voltage gradient.