Importance of the broad regional interaction for spectral tuning in Natronobacterium pharaonis phoborhodopsin (sensory rhodopsin II)

Natronobacterium pharaonis phoborhodopsin (ppR; also called N. pharaonis sensory rhodopsin II, NpsRII) is a photophobic sensor in N. pharaonis, and has a shorter absorption maximum (lambdamax, 500 nm) than those of other archaeal retinal proteins (lambdamax, 560-590 nm) such as bacteriorhodopsin (bR). We constructed chimeric proteins between bR and ppR to investigate the long range interactions effecting the color regulation among archaeal retinal proteins. The lambdamax of B-DEFG/P-ABC was 545 nm, similar to that of bR expressed in Escherichia coli (lambdamax, 550 nm). B-DEFG/P-ABC means a chimera composed of helices D, E, F, and G of bR and helices A, B, and C of ppR. This indicates that the major factor(s) determining the difference in lambdamax between bR and ppR exist in helices DEFG. To specify the more minute regions for the color determination between bR and ppR, we constructed 15 chimeric proteins containing helices D, E, F, and G of bR. According to the absorption spectra of the various chimeric proteins, the interaction between helices D and E as well as the effect of the hydroxyl group around protonated Schiff base on helix G (Thr-204 for ppR and Ala-215 for bR) are the main factors for spectral tuning between bR and ppR.     

Thermal denaturing of bacteriorhodopsin by X-Ray scattering from oriented purple membranes

We present a temperature-dependent x-ray diffraction study of thin films of purple membranes (PMs) with the native membrane protein bacteriorhodopsin (BR). The high degree of alignment with respect to the silicon substrates allows for the application of modern interface-sensitive scattering techniques. Here we focus on the structural changes of BR in PMs at the thermal denaturing transition. A partial unfolding of the helices is observed rather than the complete unfolding process known from helix to coil transitions. While BR remains threaded into the lipid bilayer in the denatured state, changes in the short-range lateral structures are associated with the partial unfolding of the transmembrane helices.     

Stability of bacteriorhodopsin alpha-helices and loops analyzed by single-molecule force spectroscopy

The combination of high-resolution atomic force microscopy imaging and single-molecule force spectroscopy allows the identification, selection, and mechanical investigation of individual proteins. In a recent paper we had used this technique to unfold and extract single bacteriorhodopsins (BRs) from native purple membrane patches. We show that subsets of the unfolding spectra can be classified and grouped to reveal detailed insight into the individualism of the unfolding pathways. We have further developed this technique and analysis to report here on the influence of pH effects and local mutations on the stability of individual structural elements of BR against mechanical unfolding. We found that, although the seven transmembrane alpha-helices predominantly unfold in pairs, each of the helices may also unfold individually and in some cases even only partially. Additionally, intermittent states in the unfolding process were found, which are associated with the stretching of the extracellular loops connecting the alpha-helices. This suggests that polypeptide loops potentially act as a barrier to unfolding and contribute significantly to the structural stability of BR. Chemical removal of the Schiff base, the covalent linkage of the photoactive retinal to the helix G, resulted in a predominantly two-step unfolding of this helix. It is concluded that the covalent linkage of the retinal to helix G stabilizes the structure of BR. Trapping mutant D96N in the M state of the proton pumping photocycle did not affect the unfolding barriers of BR.     

Molecular force modulation spectroscopy revealing the dynamic response of single bacteriorhodopsins

Recent advances in atomic force microscopy allowed globular and membrane proteins to be mechanically unfolded on a single-molecule level. Presented is an extension to the existing force spectroscopy experiments. While unfolding single bacteriorhodopsins from native purple membranes, small oscillation amplitudes (6-9 nm) were supplied to the vertical displacement of the cantilever at a frequency of 3 kHz. The phase and amplitude response of the cantilever-protein system was converted to reveal the elastic (conservative) and viscous (dissipative) contributions to the unfolding process. The elastic response (stiffness) of the extended parts of the protein were in the range of a few tens pN/nm and could be well described by the derivative of the wormlike chain model. Discrete events in the viscous response coincided with the unfolding of single secondary structure elements and were in the range of 1 microNs/m. In addition, these force modulation spectroscopy experiments revealed novel mechanical unfolding intermediates of bacteriorhodopsin. We found that kinks result in a loss of unfolding cooperativity in transmembrane helices. Reconstructing force-distance spectra by the integration of amplitude-distance spectra verified their position, offering a novel approach to detect intermediates during the forced unfolding of single proteins.     

Efficient unfolding pattern recognition in single molecule force spectroscopy data

Background

Single-molecule force spectroscopy (SMFS) is a technique that measures the force necessary to unfold a protein. SMFS experiments generate Force-Distance (F-D) curves. A statistical analysis of a set of F-D curves reveals different unfolding pathways. Information on protein structure, conformation, functional states, and inter- and intra-molecular interactions can be derived.

Results

In the present work, we propose a pattern recognition algorithm and apply our algorithm to datasets from SMFS experiments on the membrane protein bacterioRhodopsin (bR). We discuss the unfolding pathways found in bR, which are characterised by main peaks and side peaks. A main peak is the result of the pairwise unfolding of the transmembrane helices. In contrast, a side peak is an unfolding event in the alpha-helix or other secondary structural element. The algorithm is capable of detecting side peaks along with main peaks. Therefore, we can detect the individual unfolding pathway as the sequence of events labeled with their occurrences and co-occurrences special to bR's unfolding pathway. We find that side peaks do not co-occur with one another in curves as frequently as main peaks do, which may imply a synergistic effect occurring between helices. While main peaks co-occur as pairs in at least 50% of curves, the side peaks co-occur with one another in less than 10% of curves. Moreover, the algorithm runtime scales well as the dataset size increases.

Conclusions

Our algorithm satisfies the requirements of an automated methodology that combines high accuracy with efficiency in analyzing SMFS datasets. The algorithm tackles the force spectroscopy analysis bottleneck leading to more consistent and reproducible results.     

Bilayer mechanical properties regulate the transmembrane helix mobility and enzymatic state of CD39

CD39 can exist in at least two distinct functional states depending on the presence and intact membrane integration of its two transmembrane helices. In native membranes, the transmembrane helices undergo dynamic rotational motions that are required for enzymatic activity and are regulated by substrate binding. In this study, we show that bilayer mechanical properties regulate conversion between the two enzymatic functional states by modulating transmembrane helix dynamics. Alteration of membrane properties by insertion of cone-shaped or inverse cone-shaped amphiphiles or by cholesterol removal switches CD39 to the same enzymatic state that removal or solubilization of the transmembrane domains does. The same membrane alterations increase the propensity of both transmembrane helices to rotate within the packed structure, resulting in a structure with greater mobility but not an altered primary conformation. Membrane alteration also abolishes the ability of the substrate to stabilize the helices in their primary conformation, indicating a loss of coupling between substrate binding and transmembrane helix dynamics. Removal of either transmembrane helix mimics the effect of membrane alteration on the mobility and substrate sensitivity of the remaining helix, suggesting that the ends of the extracellular domain have intrinsic flexibility. We suggest that a mechanical bilayer property, potentially elasticity, regulates CD39 by altering the balance between the stability and flexibility of its transmembrane helices and, in turn, of its active site.     

Assembly of single bacteriorhodopsin trimers in bilayer nanodiscs

Nanodiscs, phospholipid bilayer assemblies of controlled size, were used to self-assemble bacteriorhodopsin (bR) into single trimers. Self-assembly at optimal bR to Nanodisc and phospholipid stoichiometry yielded particles containing three bR molecules. Analysis of solution small angle X-ray scattering indicated that bacteriorhodopsin is embedded in a discoidal phospholipid bilayer structure. Formation of trimers, as evidenced by visible circular dichroism of the retinal absorbance bands, is facilitated in Nanodiscs at a specific size threshold, suggesting that a critical bilayer area or amount of lipid is necessary to maintain a native oligomeric state. The lipid to bR ratio in the assembly process was also found to be an important factor in determining oligomerization state. These nanoscale bilayers offer the opportunity to understand and control the assembly of oligomeric integral membrane proteins critical to macromolecular recognition and cellular signaling.     

Electrostatic and steric interactions determine bacteriorhodopsin single-molecule biomechanics

Bacteriorhodopsin (bR) is a haloarchaeal membrane protein that converts the energy of single photons into large structural changes to directionally pump protons across purple membrane. This is achieved by a complex combination of local dynamic interactions controlling bR biomechanics at the submolecular level, producing efficient amplification of the retinal photoisomerization. Using single molecule force spectroscopy at different salt concentrations, we show that tryptophan (Trp) residues use steric specific interactions to create a rigid scaffold in bR extracellular region and are responsible for the main unfolding barriers. This scaffold, which encloses the retinal, controls bR local mechanical properties and anchors the protein into the membrane. Furthermore, the stable Trp-based network allows ion binding to two specific sites on the extracellular loops (BC and FG), which are involved in proton release and lateral transport. In contrast, the cytoplasmic side of bR is mainly governed by relatively weak nonspecific electrostatic interactions that provide the flexibility necessary for large cytoplasmic structural rearrangements during the photocycle. The presence of an extracellular Trp-based network tightly enclosing the retinal seems common to most haloarchaeal rhodopsins, and could be relevant to their exceptional efficiency.     

Mechanically unfolding protein L using a laser-feedback-controlled cantilever

Force spectroscopy using the atomic force microscope (AFM) can yield important information on the strength and lifetimes of the folded states of single proteins and their complexes when they are loaded with force. For example, by mechanically unfolding concatenated proteins at different velocities, a dynamic force spectrum can be built up that allows reconstruction of the energy landscape that the protein traverses during unfolding. To characterize fully the unfolding landscape, however, it is necessary both to explore the entire force spectrum and to characterize each species populated during unfolding. In the conventional AFM apparatus, force is applied to the protein construct through a compliant cantilever. This limits the dynamic range of the force spectrum that can be probed, and the cantilever recoil after unfolding may mask the presence of metastable intermediates. Here, we describe to our knowledge a new technique—constant-deflection AFM—in which the compliance of the AFM cantilever is removed. Using this technique, we show that protein L exhibits a more complex unfolding energy landscape than previously detected using the conventional technique. This technique is also able to detect the presence of a refolding intermediate whose formation is otherwise prevented by cantilever recoil.     

Sequential unfolding of individual helices of bacterioopsin observed in molecular dynamics simulations of extraction from the purple membrane

Multiple molecular dynamics simulations of bacterioopsin pulling from its C-terminus show that its alpha-helices unfold individually. In the first metastable state observed in the simulations, helix G is unfolded at its C-terminal segment while the rest of helix G (residues 200-216) is folded and opposes resistance because of a salt-bridge network consisting of Asp-212 and Lys-216 on helix G and Arg-82 and Asp-85 on helix C. Helix G unfolds inside the bundle because the external force is applied to its C-terminal end in a direction perpendicular to the surface of the membrane. Inversely, helix F has to flip by 180 degrees to exit from the membrane because the applied force and the helical N-C axis point in opposite directions. At the highest peak of the force, which cannot be interpreted in single-molecule force spectroscopy experiments, helix F has a pronounced kink at Pro-186. Mutation of Pro-186 and/or the charged side chains mentioned above, which are involved in very favorable electrostatic interactions in the low-dielectric region of the membrane, are expected to reduce the highest peak of the force. Helices E and D unfold in a similar way to helices G and F, respectively. Hence, the force-distance profile and sequence of events during forced unfolding of bacterioopsin are influenced by the up-and-down topology of the seven-helix bundle. The sequential extraction of individual helices from the membrane suggests that the spontaneous (un)folding of bacterioopsin proceeds through metastable bundles of fewer than seven helices. The metastable states observed in the simulations provide atomic level evidence that corroborates the interpretation of very recent force spectroscopy experiments of bacteriorhodopsin refolding.     

Orientational preferences of neighboring helices can drive ER insertion of a marginally hydrophobic transmembrane helix

α-helical integral membrane proteins critically depend on the correct insertion of their transmembrane α helices into the lipid bilayer for proper folding, yet a surprisingly large fraction of the transmembrane α helices in multispanning integral membrane proteins are not sufficiently hydrophobic to insert into the target membrane by themselves. How can such marginally hydrophobic segments nevertheless form transmembrane helices in the folded structure? Here, we show that a transmembrane helix with a strong orientational preference (N(cyt)-C(lum) or N(lum)-C(cyt)) can both increase and decrease the hydrophobicity threshold for membrane insertion of a neighboring, marginally hydrophobic helix. This effect helps explain the "missing hydrophobicity" in polytopic membrane proteins.     

Unfolding pathway of the colicin E1 channel protein on a membrane surface

The channel-forming domain of colicin E1 is composed of a soluble helical bundle which, upon membrane binding, unfolds to form an extended, two-dimensional helical net in the membrane interfacial layer. To characterize the pathway of unfolding of the protein and the structure of the surface-bound intermediate, the time-course of intra-protein distance changes and unfolding on a millisecond time-scale were determined from the kinetics of changes in the efficiency of fluorescence resonance energy transfer, and of the donor-acceptor overlap integral, between each of six individual tryptophan residues and a Cys-conjugated energy transfer acceptor (C509-AEDANS). Comparison of the rate constants revealed the following order of events associated with unfolding of the protein at the membrane surface: (A) movement of the hydrophobic core helices VIII-IX, coincident with a small change in Trp-Cys509 distances of the outer helices; (B) unfolding of surface helices in the helical bundle in the order: helix I, helices III, IV, VI, VII, and helix V; (C) a slow (time-scale, seconds) condensation of the surface-bound helices. The rate of protein unfolding events increased with increasing anionic lipid content. Unfolding did not occur below the lipid thermal phase transition, indicating that unfolding requires mobility in the interfacial layer. The structure of the two-dimensional membrane-bound intermediate in the steady-state was inferred to consist of a quasi-circular arrangement of eight helices embedded in the membrane interfacial layer and anchored by the hydrophobic helical hairpin. The pathway of unfolding of the colicin channel at the membrane surface, catalyzed by electrostatic and hydrophobic forces, is the first described for a membrane-active protein. It is proposed that the pathway and principles described for the colicin protein are relevant to membrane protein import.     

Revisiting Hydrophobic Mismatch with Free Energy Simulation Studies of Transmembrane Helix Tilt and Rotation

Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-μs umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to ∼10° is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation.     

Unfolding Simulations of Holomyoglobin from Four Mammals: Identification of Intermediates and β-Sheet Formation from Partially Unfolded States

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R _g ∼1.48–1.51 nm, helicity ∼75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0–3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2–11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations.     

Robust Driving Forces for Transmembrane Helix Packing

The packing structures of transmembrane helices are traditionally attributed to patterns in residues along the contact surface. In this view, besides keeping the helices confined in the membrane, the bilayer has only a minor effect on the helices structure. Here, we use two different approaches to show that the lipid environment has a crucial effect in determining the cross-angle distribution of packed helices. We analyzed structural data of a membrane proteins database. We show that the distribution of cross angles of helix pairs in this database is statistically indistinguishable from the cross-angle distribution of two noninteracting helices imbedded in the membrane. These results suggest that the cross angle is, to a large extent, determined by the tilt angle of the individual helices. We test this hypothesis using molecular simulations of a coarse-grained model that contains no specific residue interactions. These simulations reproduce the same cross-angle distribution as found in the database. As the tilt angle of a helix is dominated by hydrophobic mismatch between the protein and surrounding lipids, our results indicate that hydrophobic mismatch is the dominant factor guiding the transmembrane helix packing. Other short-range forces might then fine-tune the structure to its final configuration.     

Apolipoprotein A-I adopts a belt-like orientation in reconstituted high density lipoproteins

Apolipoprotein A-I (apoA-I) is the major protein associated with high density lipoprotein (HDL), and its plasma levels have been correlated with protection against atherosclerosis. Unfortunately, the structural basis of this phenomenon is not fully understood. Over 25 years of study have produced two general models of apoA-I structure in discoidal HDL complexes. The "belt" model states that the amphipathic helices of apoA-I are aligned perpendicular to the acyl chains of the lipid bilayer, whereas the "picket fence" model argues that the helices are aligned parallel with the acyl chains. To distinguish between the two models, various single tryptophan mutants of apoA-I were analyzed in reconstituted, discoidal HDL particles composed of phospholipids containing nitroxide spin labels at various positions along the acyl chain. We have previously used this technique to show that the orientation of helix 4 of apoA-I is most consistent with the belt model. In this study, we performed additional control experiments on helix 4, and we extended the results by performing the same analysis on the remaining 22-mer helices (helices 1, 2, 5, 6, 7, 8, and 10) of human apoA-I. For each helix, two different mutants were produced that each contained a probe Trp occurring two helical turns apart. In the belt model, the two Trp residues in each helix should exhibit maximal quenching at the same nitroxide group position on the lipid acyl chains. For the picket fence model, maximal quenching should occur at two different levels in the bilayer. The results show that the majority of the helices are in an orientation that is consistent with a belt model, because most Trp residues localized to a position about 5 A from the center of the bilayer. This study corroborates a belt hypothesis for the majority of the helices of apoA-I in phospholipid discs.     

Molecular dynamics of individual alpha-helices of bacteriorhodopsin in dimyristol phosphatidylocholine. I. Structure and dynamics.

Understanding the role of the lipid bilayer in membrane protein structure and dynamics is needed for tertiary structure determination methods. However, the molecular details are not well understood. Molecular dynamics computer calculations can provide insight into these molecular details of protein:lipid interactions. This paper reports on 10 simulations of individual alpha-helices in explicit lipid bilayers. The 10 helices were selected from the bacteriorhodopsin structure as representative alpha-helical membrane folding components. The bilayer is constructed of dimyristoyl phosphatidylcholine molecules. The only major difference between simulations is the primary sequence of the alpha-helix. The results show dramatic differences in motional behavior between alpha-helices. For example, helix A has much smaller root-mean-squared deviations than does helix D. This can be understood in terms of the presence of aromatic residues at the interface for helix A that are not present in helix D. Additional motions are possible for the helices that contain proline side chains relative to other amino acids. The results thus provide insight into the types of motion and the average structures possible for helices within the bilayer setting and demonstrate the strength of molecular simulations in providing molecular details that are not directly visualized in experiments.     

Tuning the Photocycle Kinetics of Bacteriorhodopsin in Lipid Nanodiscs

Monodisperse lipid nanodiscs are particularly suitable for characterizing membrane protein in near-native environment. To study the lipid-composition dependence of photocycle kinetics of bacteriorhodopsin (bR), transient absorption spectroscopy was utilized to monitor the evolution of the photocycle intermediates of bR reconstituted in nanodiscs composed of different ratios of the zwitterionic lipid (DMPC, dimyristoyl phosphatidylcholine; DOPC, dioleoyl phosphatidylcholine) to the negatively charged lipid (DOPG, dioleoyl phosphatidylglycerol; DMPG, dimyristoyl phosphatidylglycerol). The characterization of ion-exchange chromatography showed that the negative surface charge of nanodiscs increased as the content of DOPG or DMPG was increased. The steady-state absorption contours of the light-adapted monomeric bR in nanodiscs composed of different lipid ratios exhibited highly similar absorption features of the retinal moiety at 560 nm, referring to the conservation of the tertiary structure of bR in nanodiscs of different lipid compositions. In addition, transient absorption contours showed that the photocycle kinetics of bR was significantly retarded and the transient populations of intermediates N and O were decreased as the content of DMPG or DOPG was reduced. This observation could be attributed to the negatively charged lipid heads of DMPG and DOPG, exhibiting similar proton relay capability as the native phosphatidylglycerol (PG) analog lipids in the purple membrane. In this work, we not only demonstrated the usefulness of nanodiscs as a membrane-mimicking system, but also showed that the surrounding lipids play a crucial role in altering the biological functions, e.g., the ion translocation kinetics of the transmembrane proteins.     

Interactions of alpha-helices with lipid bilayers: a review of simulation studies

Membrane proteins, of which the majority seem to contain one or more alpha-helix, constitute approx. 30% of most genomes. A complete understanding of the nature of helix/bilayer interactions is necessary for an understanding of the structural principles underlying membrane proteins. This review describes computer simulation studies of helix/bilayer interactions. Key experimental studies of the interactions of alpha-helices and lipid bilayers are briefly reviewed. Surface associated helices are found in some membrane-bound enzymes (e.g. prostaglandin synthase), and as stages in the mechanisms of antimicrobial peptides and of pore-forming bacterial toxins. Transmembrane alpha-helices are found in most integral membrane proteins, and also in channels formed by amphipathic peptides or by bacterial toxins. Mean field simulations, in which the lipid bilayer is approximated as a hydrophobic continuum, have been used in studies of membrane-active peptides (e.g. alamethicin, melittin, magainin and dermaseptin) and of simple membrane proteins (e.g. phage Pf1 coat protein). All atom molecular dynamics simulations of fully solvated bilayers with transmembrane helices have been applied to: the constituent helices of bacteriorhodopsin; peptide-16 (a simple model TM helix); and a number of pore-lining helices from ion channels. Surface associated helices (e.g. melittin and dermaseptin) have been simulated, as have alpha-helical bundles such as bacteriorhodopsin and alamethicin. From comparison of the results from the two classes of simulation, it emerges that a major theoretical challenge is to exploit the results of all atom simulations in order to improve the mean field approach.