Search

We applied a new large two membrane model of human urea transporter B (UT-B) with urea concentration gradient between intracellular and extracellular solution to explore the structure and dynamics of urea permeation. The first solved human UT-B X-ray crystal structure (PDB 6QD5) was studied using atomistic force field based large scale molecular dynamics (MD) simulations. The urea binding sites in UT-B were identified by the combination of first principles molecular dynamics docking and hybrid steered molecular dynamics (hSMD). In the simulations, urea molecules bound to the two binding sites although urea molecules did not go through the channels. The potential of mean force (PMF) of urea along the channel was calculated by hSMD to investigate the selectivity of the urea permeation across the channel pore. The simulations provide an opportunity to explore the binding sites spontaneously to help understand the urea permeation gating and transportation mechanism. A clearer analysis of the selectivity filter is shown in more detail. We also furthered the research about small organic molecules in the specific environment, such as channels in this work. [ABSTRACT FROM AUTHOR]

Dopamine transporter (DAT) and sigma‐1 receptor (σ1R) are potential therapeutic targets to reduce the psychostimulant effects induced by methamphetamine (METH). Interaction of σ1R with DAT could modulate the binding of METH, but the molecular basis of the association of the two transmembrane proteins and how their interactions mediate the binding of METH to DAT or σ1R remain unclear. Here, we characterize the protein–ligand and protein–protein interactions at a molecular level by various theoretical approaches. The present results show that METH adopts a different binding pose in the binding pocket of σ1R and is more likely to act as an agonist. The relatively lower binding affinity of METH to σ1R supports the role of antagonists as inhibitors that protect against METH‐induced effects. We demonstrate that σ1R could bind to Drosophila melanogaster DAT (dDAT) through interactions with either the transmembrane helix α12 or α5 of dDAT. Our results showed that the truncated σ1R displays stronger association with dDAT than the full‐length σ1R. Although different helix–helix interactions between σ1R and dDAT lead to distinct effects on the dynamics of individual protein, both associations attenuate the binding affinity of METH to dDAT, particularly in the interactions with the helix α5 of dDAT. Together, the present study provides the first computational investigation on the molecular mechanism of coupling METH binding and the association of σ1R with dDAT. [ABSTRACT FROM AUTHOR]

Cytochrome P450 2J2 (CYP2J2) is responsible for the epoxidation of endogenous arachidonic acid, and is involved in the metabolism of exogenous drugs. To date, no crystal structure of CYP2J2 is available, and the proposed structural basis for the substrate recognition and specificity in CYP2J2 varies with the structural models developed using different computational protocols. In this study, we developed a new structural model of CYP2J2, and explored its sensitivity to substrate binding by molecular dynamics simulations of the interactions with chemically similar fluorescent probes. Our results showed that the induced-fit binding of these probes led to the preferred active poses ready for the catalysis by CYP2J2. Divergent conformational dynamics of CYP2J2 due to the binding of each probe were observed. However, a stable hydrophobic clamp composed of residues I127, F310, A311, V380, and I487 was identified to restrict any substrate access to the active site of CYP2J2. Molecular docking of a series of compounds including amiodarone, astemizole, danazol, ebastine, ketoconazole, terfenadine, terfenadone, and arachidonic acid to CYP2J2 confirmed the role of those residues in determining substrate binding and specificity of CYP2J2. In addition to the flexibility of CYP2J2, the present work also identified other factors such as electrostatic potential in the vicinity of the active site, and substrate strain energy and property that have implications for the interpretation of CYP2J2 metabolism. [ABSTRACT FROM AUTHOR]

Abstract The ubiquitous glucose transporter 1 (GLUT1) is physiologically and pathologically relevant in energy metabolism of the CNS, skeletal muscles, cancer cells etc. Extensive experiments on GLUT1 produced thorough understandings of its expressions, functions, and structures which were recently resolved to atomic accuracy. However, theoretical understandings are still controversial about how GLUT1 facilitates glucose diffusion across the cell membrane. Molecular dynamics (MD) simulations of the current literature have GLUT1 embedded in a symmetric bilayer of a single lipid type. They provide atomistic illustrations of the alternating access theory (AAT), but the simulation results are inconsistent with the undisputed experimental data of kinetics showing rapid transport of glucose at near-physiological temperatures, high Arrhenius activation barrier in zero-trans uptake, and large trans-acceleration at sub-physiological temperatures. In this research, we embedded GLUT1 in an asymmetric bilayer of multiple lipids to better mimic the erythrocyte membrane. We ran unbiased MD simulations at 37 °C and at 5 °C and found a new mechanism of glucose transport via GLUT1: The extracellular (EC) gate opened wide for EC glucopyranose at 37 °C and, only in the presence of intracellular (IC) glucose, at 5 °C. In the absence of IC glucose at 5 °C, the EC gate opened narrowly for acyclic glucose, gating out glucopyranose. This EC-gating mechanism is simpler than AAT and yet it well explains for the rapid glucose transport at near-physiological temperatures and large trans-acceleration at sub-physiological temperatures. It also explains why zero-trans uptake (involving the pyranose-to-aldehyde transformation) has an Arrhenius barrier ∼20 kcal/mol higher than the equilibrium exchange transport. Graphical abstract Image 1 Highlights • EC-gating activities found in unbiassed simulations of human GLUT1 in erythrocyte membrane. • EC-gating mechanism simpler than the alternating access mechanism for glucose transport. • EC-gate widely open for extracellular glucopyranose at physiological temperature. • EC-gating dependence on intracellular glucose at sub-physiological temperatures explaining full-range kinetics. [ABSTRACT FROM AUTHOR]

Glucose transporters (GLUTs), expressed in all types of human cells, are responsible for the uptake of sugars as the primary energy source for the normal functions of good cells and for the abnormal growth of cancer cells. The E. coli xylose permease (XylE), a homologue of human GLUTs, has been investigated more thoroughly than other major facilitator proteins in the current literature. In this paper, we present a molecular dynamics (MD) study of an all-atom model system to elucidate the atomistic details and the free-energy landscape along the path of binding a xylopyranose (XYP) from the extracellular space to the inside of the transporter protein XylE. From the MD simulations, the Gibbs free energy of binding was found to be − 4.4 kcal/mol in agreement with the experimental value of − 4.7 kcal/mol. The accuracy of our study is further shown in the computed hydration energy of XYP of − 14.6 kcal/mol in comparison with the experimental data of − 15.0 kcal/mol. Along the binding path, the Gibbs free energy of the XYP-XylE complex first rises from zero in the dissociated state to approximately 4 kcal/mol in the transition state (when XylE slightly increases its opening toward the extracellular side to accommodate XYP) before dropping down to − 9.0 kcal/mol in the bound state. These quantitative insights indicate the fast equilibration between the bound and the unbound states of XylE and XYP. They also serve as an atomistic-dynamic corroboration of the experimental conclusion that XylE is a high-affinity sugar transporter. [ABSTRACT FROM AUTHOR]

Amyloid-beta (Aβ) peptides, Aβ40 and the more neurotoxic Aβ42, have been the subject of many research efforts for Alzheimer's disease. In two recent independent investigations, the atomistic structure of Aβ42 fibril has been clearly established in the S-shaped conformation consisting of three β-sheets stabilized by salt bridges formed between the Lys28 sidechain and the C-terminus of Ala42. This structure distinctively differs from the long-known structure of Aβ40 in the β-hairpin shaped conformation consisting of two β-sheets. Recent in silico investigations based on all-atom models have reached closer agreement with the in vitro measurements of Aβ40 thermodynamics. In this study, we present an in silico investigation of Aβ42 thermodynamics. Using the established force field parameters in seven sets of all-atom simulations, we examined the stability of small Aβ42 oligomers in physiological saline. We computed the elongation affinity of the S-shaped Aβ42 fibril, reaching agreement with the experimental data. We also estimated the Arrhenius activation barrier along the elongation pathway (from the disordered conformation of a free Aβ42 peptide to its S-shaped conformation on a fibril) that amounts to about 16 kcal/mol, which is consistent with the experimental data. Based on these quantitative agreements, we conclude that aggregation of Aβ42 peptides into fibrils is thermodynamically slow without precipitation by extrinsic factors such as heparan sulfate proteoglycan and highlight the possibility to prevent Aβ42 aggregation by eliminating some precipitation factors or by increasing competitive agents to capture and transport free Aβ42 peptides from the cerebrospinal fluid. [ABSTRACT FROM AUTHOR]

Mh OR5, an insect olfactory receptor (OR), has an occluded binding site for the odorant eugenol in both the open and closed states of the ion channel. We used atomistic molecular dynamics simulation (MD) and steered molecular dynamics to examine possible tunnels to the odorant binding site from the protein surface. Four high probability tunnels were identified in the MD results. Surprisingly, three of the tunnels connect the ligand binding site to the lipid bilayer. We found sharp 30%–50% increases or decreases in tunnel bottleneck areas over 70 nsec MD trajectories, both in the ligand-bound and unliganded OR structures. Steered MD showed that eugenol follows the tunnels to the protein surface, and the potential of mean force is quantitatively consistent with the known affinity of eugenol for Mh OR5. We examined AlphaFold-generated models of 21 other insect ORs, and we found that 19 had odorant binding sites and tunnels in similar positions to Mh OR5. The possibility of a tunnel between the odorant binding site and the lipid bilayer in insect ORs suggests new experiments to test molecular mechanisms for insect odorant reception. [Display omitted] • The insect olfactory receptor Mh OR5 was analyzed by molecular dynamics simulation. • Tunnels connect occluded ligand-binding sites to the exterior and to the lipid bilayer. • Tunnel bottlenecks undergo large area fluctuations over 5–10 ns. • Steered molecular dynamics shows ligand can enter or exit Mh OR5 via the tunnels. [ABSTRACT FROM AUTHOR]

Glucose transporter GLUT1 is ubiquitously expressed in the human body from the red cells to the blood-brain barrier to the skeletal muscles. It is physiologically relevant to understand how GLUT1 facilitates diffusion of glucose across the cell membrane. It is also pathologically relevant because GLUT1 deficiency causes neurological disorders and anemia and because GLUT1 overexpression fuels the abnormal growth of cancer cells. This article presents a quantitative investigation of GLUT1 based on all-atom molecular-dynamics (MD) simulations of the transporter embedded in lipid bilayers of asymmetric inner-and-outer-leaflet lipid compositions, subject to asymmetric intra-and-extra-cellular environments. This is in contrast with the current literature of MD studies that have not considered both of the aforementioned asymmetries of the cell membrane. The equilibrium (unbiased) dynamics of GLUT1 shows that it can facilitate glucose diffusion across the cell membrane without undergoing large-scale conformational motions. The Gibbs free-energy profile, which is still lacking in the current literature of GLUT1, quantitatively characterizes the diffusion path of glucose from the periplasm, through an extracellular gate of GLUT1, on to the binding site, and off to the cytoplasm. This transport mechanism is validated by the experimental data that GLUT1 has low water-permeability, uptake-efflux symmetry, and 10 kcal/mol Arrhenius activation barrier around 37 °C. [Display omitted] • Model system of two membrane patches asymmetric in EC and IC salines • Model RBC membrane asymmetric in inner and outer leaflet lipid compositions • EC-gating mechanism of GLUT1 in quantitative agreement with in vitro data on glucose transport and water permeation [ABSTRACT FROM AUTHOR]

Human carbonic anhydrase II (hCAII) represents an ultimate example of the perfectly efficient metalloenzymes, which is capable of catalyzing the hydration of carbon dioxide with a rate approaching the diffusion controlled limit. Extensive experimental studies of this physiologically important metalloprotein have been done to elucidate the fundamentals of its enzymatic actions: what residues anchor the Zn 2+ (or another divalent cation) at the bottom of the binding pocket; how the relevant residues work concertedly with the divalent cation in the reversible conversions between CO 2 and HCO 3 – ; what are the protonation states of the relevant residues and acetazolamide, an inhibitor complexed with hCAII, etc. In this article, we present a detailed computational study on the basis of the all-atom CHARMM force field where Zn 2+ is represented with a simple model of divalent cation using the transferrable parameters available from the current literature. We compute the hydration free energy of Zn 2+ , the characteristics of hCAII-Zn 2+ complexation, and the absolute free energy of binding acetazolamide to the hCAII-Zn 2+ complex. In each of these three problems, our computed results agree with the experimental data within the known margin of error without making any case-by-case adjustments to the parameters. The quantitatively accurate insights we gain in this all-atom molecular dynamics study should be helpful in the search and design of more specific inhibitors of this and other carbonic anhydrases. [ABSTRACT FROM AUTHOR]

Nanometer-sized gold particles (AuNPs) are of peculiar interest because their behaviors in an aqueous solution are sensitive to changes in environmental factors including the size and shape of the solute ions. In order to determine these important characteristics, we performed all-atom molecular dynamics simulations on the icosahedral Au144 nanoparticles each coated with a homogeneous set of 60 thiolates (4-mercaptobenzoate, pMBA) in eight aqueous solutions having ions of varying sizes and shapes (Na+, K+, tetramethylamonium cation TMA+, tris-ammonium cation TRS+, Cl-, and OH-). For each solution, we computed the reversible work (potential of mean of force) to bring two nanoparticles together as a function of their separation distance. We found that the behavior of pMBA protected Au144 nanoparticles can be readily modulated by tuning their aqueous environmental factors (pH and solute ion combinations). We examined the atomistic details on how the sizes and shapes of solute ions quantitatively factor in the definitive characteristics of nanoparticle-environment and nanoparticle-nanoparticle interactions. We predict that tuning the concentrations of nonspherical composite ions such as TRS+ in an aqueous solution of AuNPs be an effective means to modulate the aggregation propensity desired in biomedical and other applications of small charged nanoparticles. [ABSTRACT FROM AUTHOR]

Plasmodium falciparum aquaporin (PfAQP) is a multifunctional channel protein in the plasma membrane of the malarial parasite that causes the most severe form of malaria infecting more than a million people a year. This channel protein facilitates transport of water and several solutes across the cell membrane. In order to better elucidate the fundamental interactions between PfAQP and its permeants and among the permeants, I conducted over three microseconds in silico experiments of atomistic models of the PfAQP–membrane system to obtain the free-energy profiles of five permeants (erythritol, water, glycerol, urea, and ammonia) throughout the amphipathic conducting pore of PfAQP. The profiles are analyzed in light of and shown to be consistent with the existent in vitro data. The binding affinities are computed using the free-energy profiles and the permeant fluctuations inside the channel. On this basis, it is predicted that erythritol, a permeant of PfAQP itself having a deep ditch in its permeation passageway, inhibits PfAQP's functions of transporting water and other solutes with an IC 50 in the range of high nanomolars. This leads to the possibility that erythritol, a sweetener generally considered safe, may inhibit or kill the malarial parasite in vivo without causing undesired side effects. Experimental studies are hereby called for to directly test this theoretical prediction of erythritol strongly inhibiting PfAQP in vitro and possibly inhibiting P. falciparum in vivo . [ABSTRACT FROM AUTHOR]

In order to elucidate how phosphatidylserine (PS6) interacts with AQP5 in a cell membrane, we developed a hybrid steered molecular dynamics (hSMD) method that involved: (1) Simultaneously steering two centers of mass of two selected segments of the ligand, and (2) equilibrating the ligand-protein complex with and without biasing the system. Validating hSMD, we first studied vascular endothelial growth factor receptor 1 (VEGFR1) in complex with N-(4-Chlorophenyl)-2-((pyridin-4-ylmethyl)amino)benzamide (8ST), for which the binding energy is known fromin vitroexperiments. In this study, our computed binding energy well agreed with the experimental value. Knowing the accuracy of this hSMD method, we applied it to the AQP5-lipid-bilayer system to answer an outstanding question relevant to AQP5’s physiological function: Will the PS6, a lipid having a single long hydrocarbon tail that was found in the central pore of the AQP5 tetramer crystal, actually bind to and inhibit AQP5’s central pore under near-physiological conditions, namely, when AQP5 tetramer is embedded in a lipid bilayer? We found,in silico, using the CHARMM 36 force field, that binding PS6 to AQP5 was a factor of 3 million weaker than “binding” it in the lipid bilayer. This suggests that AQP5’s central pore will not be inhibited by PS6 or a similar lipid in a physiological environment. [ABSTRACT FROM PUBLISHER]

Dopamine transporter mediates the neurotransmitter dopamine homeostasis in a sodium-dependent manner. The transport process involves an alternating access of a substrate to the extracellular and intracellular spaces, which is associated with different conformational states of the transporter. However, the underlying mechanism of modulation of the state transition remains elusive. Here we present a computational simulation study of human dopamine transporter to explore its two end states (outward-facing open and inward-facing open) that have not been determined experimentally. We show that the full-length transporter may tend to adopt the inward-facing open state in its free state. The binding of an amphetamine may not trap the transporter in the outward-facing open state with increasing length of the N-terminal. Furthermore, we identify distinct patterns in the interaction networks between the N-terminal and the intracellular region that could stabilize the state of the transporter, independent of substrate binding and phosphorylation. Our results reveal the essential role of the N-terminal dynamics in modulating the functional states of the dopamine transporter, providing molecular insights into the coupling of conformational transition and substrate passage in neurotransmitter transporters. When the N-terminal interacts with the intracellular region in an open manner, hDAT exhibits an outward-open conformation; whereas when the N-terminal collapses in the cytosol or interacts with the intracellular region in a closed manner, hDAT displays an inward-open conformation. [Display omitted] • The flexible N-terminal could modulate the conformational states of human dopamine transporter • The N-terminal Phosphorylation could increase the inward-open states of human dopamine transporter. • The collapsed N-terminal results in an inward-open state of human dopamine transporter. • The human dopamine transporter may share the same substrate transport path as human serotonin transporter. [ABSTRACT FROM AUTHOR]

Abstract: Among aquaglyceroporins that transport both water and glycerol across the cell membrane, Escherichia coli glycerol uptake facilitator (GlpF) is the most thoroughly studied. However, one question remains: Does glycerol modulate water permeation? This study answers this fundamental question by determining the three-dimensional potential of mean force of glycerol along the permeation path through GlpF's conducting pore. There is a deep well near the Asn-Pro-Ala (NPA) motifs (6.5kcal/mol below the bulk level) and a barrier near the selectivity filter (10.1kcal/mol above the well bottom). This profile owes its existence to GlpF's perfect steric arrangement: The glycerol–protein van der Waals interactions are attractive near the NPA but repulsive elsewhere in the conducting pore. In light of the single-file nature of waters and glycerols lining up in GlpF's amphipathic pore, it leads to the following conclusion: Glycerol modulates water permeation in the μM range. At mM concentrations, GlpF is glycerol-saturated and a glycerol residing in the well occludes the conducting pore. Therefore, water permeation is fully correlated to glycerol dissociation that has an Arrhenius activation barrier of 6.5kcal/mol. Validation of this theory is based on the existent in vitro data, some of which have not been given the proper attention they deserved: The Arrhenius activation barriers were found to be 7kcal/mol for water permeation and 9.6kcal/mol for glycerol permeation; The presence of up to 100mM glycerol did not affect the kinetics of water transport with very low permeability, in apparent contradiction with the existent theories that predicted high permeability (0M glycerol). [Copyright &y& Elsevier]

Abstract: Plasmodium falciparum aquaglyceroporin (PfAQP) is a multifunctional membrane protein in the plasma membrane of P. falciparum, the parasite that causes the most severe form of malaria. The current literature has established the science of PfAQP’s structure, functions, and hydrogen-bonding interactions but left unanswered the following fundamental question: does glycerol modulate water permeation through aquaglyceroporin that conducts both glycerol and water? This paper provides an affirmative answer to this question of essential importance to the protein’s functions. On the basis of the chemical-potential profile of glycerol from the extracellular bulk region, throughout PfAQP’s conducting channel, to the cytoplasmic bulk region, this study shows the existence of a bound state of glycerol inside aquaglyceroporin’s permeation pore, from which the dissociation constant is approximately 14μM. A glycerol molecule occupying the bound state occludes the conducting pore through which permeating molecules line up in single file by hydrogen-bonding with one another and with the luminal residues of aquaglyceroporin. In this way, glycerol inhibits permeation of water and other permeants through aquaglyceroporin. The biological implications of this theory are discussed and shown to agree with the existent in vitro data. It turns out that the structure of aquaglyceroporin is perfect for the van der Waals interactions between the protein and glycerol to cause the existence of the bound state deep inside the conducting pore and, thus to play an unexpected but significant role in aquaglyceroporin’s functions. [Copyright &y& Elsevier]

Abstract For its fundamental relevance, transport of water and glycerol across the erythrocyte membrane has long been investigated before and after the discovery of aquaporins (AQPs), the membrane proteins responsible for water and glycerol transport. AQP1 is abundantly expressed in the human erythrocyte for maintaining its hydrohomeostasis where AQP3 is also expressed (at a level ~30-folds lower than AQP1) facilitating glycerol transport. This research is focused on two of the remaining questions: How permeable is AQP3 to water? What is the glycerol-AQP3 affinity under near-physiological conditions? Through atomistic modelling and large-scale simulations, we found that AQP3 is two to three times more permeable to water than AQP1 and that the glycerol-AQP3 affinity is approximately 500/M. Using these computed values along with the data from the latest literature on AQP1 and on erythrocyte proteomics, we estimated the water and glycerol transport rates across the membrane of an entire erythrocyte. We used these rates to predict the time courses of erythrocyte swelling-shrinking in response to inward and outward osmotic gradients. Experimentally, we monitored the time course of human erythrocytes when subject to an osmotic or glycerol gradient with light scattering in a stopped-flow spectrometer. We observed close agreement between the experimentally measured and the computationally predicted time courses of erythrocytes, which corroborated our computational conclusions on the AQP3 water-permeability and the glycerol-AQP3 affinity. Graphical abstract Unlabelled Image Highlights • Human Aquaporin 3 is functionally similar to its bacterium and parasite analogues. • Aquaporin 3, a glycerol channel, is more permeable than Aquaporin 1. • Aquaporin 3 has an affinity of 500/M for glycerol. [ABSTRACT FROM AUTHOR]

Measuring or computing the single-channel permeability of aquaporins/aquaglyceroporins (AQPs) has long been a challenge. The measured values scatter over an order of magnitude but the corresponding Arrhenius activation energies converge in the current literature. Osmotic flux through an AQP was simulated as water current forced through the channel by kilobar hydraulic pressure or theoretically approximated as single-file diffusion. In this paper, we report large scale simulations of osmotic current under sub M gradient through three AQPs (water channels AQP4 and AQP5 and glycerol-water channel GlpF) using the mature particle mesh Ewald technique (PME) for which the established force fields have been optimized with known accuracy. These simulations were implemented with hybrid periodic boundary conditions devised to avoid the artifactitious mixing across the membrane in a regular PME simulation. The computed single-channel permeabilities at 5 °C and 25 °C are in agreement with recently refined experiments on GlpF. The Arrhenius activation energies extracted from our simulations for all the three AQPs agree with the in vitro measurements. The single-file diffusion approximations from our large-scale simulations are consistent with the current literature on smaller systems. From these unambiguous agreements among the in vitro and in silico studies, we observe the quantitative accuracy of the all-atom force fields of the current literature for water-channel biology. We also observe that AQP4, that is particularly rich in the central nervous system, is more efficient in water conduction and more temperature-sensitive than other water-only channels (excluding glycerol channels that also conduct water when not inhibited by glycerol). [ABSTRACT FROM AUTHOR]

Computing the ligand-protein binding affinity (or the Gibbs free energy) with chemical accuracy has long been a challenge for which many methods/approaches have been developed and refined with various successful applications. False positives and, even more harmful, false negatives have been and still are a common occurrence in practical applications. Inevitable in all approaches are the errors in the force field parameters we obtain from quantum mechanical computation and/or empirical fittings for the intra- and inter-molecular interactions. These errors propagate to the final results of the computed binding affinities even if we were able to perfectly implement the statistical mechanics of all the processes relevant to a given problem. And they are actually amplified to various degrees even in the mature, sophisticated computational approaches. In particular, the free energy perturbation (alchemical) approaches amplify the errors in the force field parameters because they rely on extracting the small differences between similarly large numbers. In this paper, we develop a hybrid steered molecular dynamics (hSMD) approach to the difficult binding problems of a ligand buried deep inside a protein. Sampling the transition along a physical (not alchemical) dissociation path of opening up the binding cavity---pulling out the ligand---closing back the cavity, we can avoid the problem of error amplifications by not relying on small differences between similar numbers. We tested this new form of hSMD on retinol inside cellular retinol-binding protein 1 and three cases of a ligand (a benzylacetate, a 2-nitrothiophene, and a benzene) inside a T4 lysozyme L99A/M102Q(H) double mutant. In all cases, we obtained binding free energies in close agreement with the experimentally measured values. This indicates that the force field parameters we employed are accurate and that hSMD (a brute force, unsophisticated approach) is free from the problem of error amplification suffered by many sophisticated approaches in the literature. [ABSTRACT FROM AUTHOR]

Abstract: To probe structural changes that occur when a membrane protein is transferred from lipid bilayers to SDS micelles, a fragment of bacteriorhodopsin containing transmembrane helical segments A and B was studied by fluorescence spectroscopy, molecular dynamics (MD) simulation, and stopped flow kinetics. In lipid bilayers, Förster resonance energy transfer (FRET) was observed between tyrosine 57 on helix B and tryptophans 10 and 12 on helix A. FRET efficiency decreased substantially when the peptide was transferred to SDS. MD simulation showed no evidence for significant disruption of helix–helix interactions in SDS micelles. However, a cluster of water molecules was observed to form a hydrogen-bonded network with the phenolic hydroxyl group of tyrosine 57, which probably causes the disappearance of tyrosine-to-tryptophan FRET in SDS. The tryptophan quantum yield decreased in SDS, and the change occurred at nearly the same rate as membrane solubilization. The results provide a clear example of the importance of corroborating distance changes inferred from FRET by using complementary methods. [Copyright &y& Elsevier]