Biochimica et Biophysica Acta
|
|
- Dylan Johnson
- 5 years ago
- Views:
Transcription
1 Biochimica et Biophysica Acta 1818 (2012) Contents lists available at ScienceDirect Biochimica et Biophysica Acta journal homepage: Review Molecular simulations and solid-state NMR investigate dynamical structure in rhodopsin activation Blake Mertz a, Andrey V. Struts a,b, Scott E. Feller c, Michael F. Brown a,d, a Department of Chemistry and Biochemistry, University of Arizona, Tucson, AZ 85721, USA b Department of Medical Physics, St. Petersburg State Medical University, St. Petersburg , Russia c Department of Chemistry, Wabash College, Crawfordsville, IN 47933, USA d Department of Physics, University of Arizona, Tucson, AZ 85721, USA article info abstract Article history: Received 26 June 2011 Received in revised form 1 August 2011 Accepted 1 August 2011 Available online 8 August 2011 Keywords: G protein-coupled receptor Membrane Molecular dynamics Solid-state NMR Rhodopsin Vision Rhodopsin has served as the primary model for studying G protein-coupled receptors (GPCRs) the largest group in the human genome, and consequently a primary target for pharmaceutical development. Understanding the functions and activation mechanisms of GPCRs has proven to be extraordinarily difficult, as they are part of a complex signaling cascade and reside within the cell membrane. Although X-ray crystallography has recently solved several GPCR structures that may resemble the activated conformation, the dynamics and mechanism of rhodopsin activation continue to remain elusive. Notably solid-state 2 H NMR spectroscopy provides key information pertinent to how local dynamics of the retinal ligand change during rhodopsin activation. When combined with molecular mechanics simulations of proteolipid membranes, a new paradigm for the rhodopsin activation process emerges. Experiment and simulation both suggest that retinal isomerization initiates the rhodopsin photocascade to yield not a single activated structure, but rather an ensemble of activated conformational states. This article is part of a Special Issue entitled: Membrane protein structure and function Published by Elsevier B.V. Contents 1. Introduction Rhodopsin: interface between experiment and simulation Solid-state 2 H NMR line shapes are related to dynamical structure of retinal ligand Molecular simulations further characterize structural dynamics and fluctuations of membrane proteins Molecular dynamics calculations allow testing of specific counterion models for rhodopsin Solid-state 2 H NMR spectroscopy investigates the dynamics of cofactors bound to membrane proteins Relaxation in 2 H NMR spectroscopy provides key insights into retinal mobility Methyl groups probe site-specific local interactions in the inactive dark state Changes in retinylidene methyl dynamics occur in the preactive Meta I state and signaling Meta II state Analyticalrelaxationtheoryandnumericaldynamicssimulationscharacterizetheenergylandscapeofrhodopsinactivation Solid-state 2 H NMR highlights the functional motions of retinal methyl groups Molecular dynamics simulations expand the 2 H NMR results to longer time scales Quantum mechanical calculations address methyl dihedral terms in the retinal force field Larger polyene compounds are necessary to model the retinylidene methyl groups Hyperconjugation extends to methyl groups of the retinal chromophore implicated in rhodopsin activation Putting it all together: the rhodopsin activation mechanism Abbreviations: CHARMM, Chemistry at Harvard Macromolecular Mechanics; EPR, electron paramagnetic resonance; FTIR, Fourier transform infrared; GPCR, G protein-coupled receptor; MD, molecular dynamics; Meta I, metarhodopsin I; Meta II, metarhodopsin II; 2MBD, 2-methyl-butadiene; 3MHT, 3-methyl-hexatriene; MM, molecular mechanics; POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine; PSB, protonated Schiff base; RDC, residual dipolar coupling; RQC, residual quadrupolar coupling; SDPC, 1-stearoyl-2- docosahexaenoyl-sn-glycero-3-phosphocholine; SDPE, 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphoethanolamine; QM, quantum mechanics This article is part of a Special Issue entitled: Membrane protein structure and function. Corresponding author at: Department of Chemistry and Biochemistry, University of Arizona, Tucson, AZ 85721, USA. Tel.: ; fax: address: mfbrown@u.arizona.edu (M.F. Brown) /$ see front matter 2011 Published by Elsevier B.V. doi: /j.bbamem
2 242 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) Retinal methyl groups assist in driving multiscale rhodopsin activation Ionic locks hold the key to rhodopsin transducin interactions Moving forward Acknowledgements References Introduction Rhodopsin is the G protein-coupled receptor (GPCR) responsible for dim light vision and plays a prominent role in our understanding of biological signaling. Investigating GPCRs lies at the forefront of pharmaceutical research, as they represent almost one-half of current drug targets [1]. In fact, rhodopsin is the most well-studied GPCR, on account of its relative ease of procurement and its well-characterized spectroscopic intermediates [2,3]. Upon photon absorption, an 11- cis all-trans isomerization of the covalently bound inverse agonist retinal leads to initiation of the rhodopsin photocascade. This isomerization causes conformational changes within rhodopsin that allow interaction with the G protein transducin (G t ) in the cytosol. The signaling protein transducin acts as an intermediary between rhodopsin and its effector cgmp phosphodiesterase. Catalysis of GDP GTP exchange by rhodopsin allows transducin to activate the cgmp phosphodiesterase, yielding hydrolysis of cgmp and closing of cyclic nucleotide-gated ion channels in the rod cellular plasma membrane. This last step leads to hyperpolarization of the rod and generation of a visual nerve impulse. Although the above signaling cascade is a blueprint for GPCR function, it remains incompletely understood at present particularly for rhodopsin activation at the molecular level. In this regard, structural biology has contributed insights into the workings of this system that are broadly significant for membrane function. The first GPCR crystal structure was of rhodopsin in the dark state [4]. Subsequent crystal structures of rhodopsin [5 10], opsin [11,12], the β 1 - and β 2 -adrenergic receptors [13 17], and the adenosine A 2A receptor [18] have now been solved, offering the potential for further mechanistic insight. Nonetheless, despite the availability of a number of GPCR structures, X-ray crystallography cannot completely explain membrane protein function. No structure has been solved in a native lipid bilayer environment. Additionally, each structure is a snapshot of the protein, and is unable to fully reveal the dynamics. Recently, our understanding of rhodopsin activation has been further advanced by reports of X-ray structures for the active Meta II state yet, the results are complicated by opposite orientations for the retinal ligand obtained by different research groups [9,10]. Application of spectroscopic techniques such as nuclear magnetic resonance (NMR) are thus needed to provide structural data for trapped photointermediates in order to fill these essential gaps in knowledge of the activation process [19,20]. Below we summarize a solid-state 2 H NMR approach that is unique in that it can be used to investigate both structure and dynamics of retinal within the rhodopsin binding pocket [21]. The 2 H NMR method is mainly limited to the immediate vicinity of isotopic labels. However, computational simulations can add to the 2 H NMR data by examining the entire proteolipid system, thereby contributing to a more complete picture of the receptor activation process. Here we demonstrate how experiment and theory have a synergistic interplay that yields new insights into the activation mechanism of rhodopsin in a nativelike membrane environment. Solid-state NMR spectroscopy probes the changes in the local structure and dynamics of the retinal ligand of rhodopsin at a site-specific level. In addition, molecular dynamics simulations allow extension to longer time scales, thus enabling further investigation of the energy landscape of rhodopsin activation. Numerical simulations require accurate parameterization of the relevant molecular mechanics force fields, which has motivated us to carry out quantum mechanical calculations for retinal model compounds. Last, we show how the retinal methyl groups are implicated in rhodopsin activation, involving subsequent interaction with the signaling G protein transducin in visual excitation. 2. Rhodopsin: interface between experiment and simulation Rhodopsin, the canonical GPCR, is composed of seven transmembrane helices connected by a series of extracellular and cytoplasmic loops [3]. Unlike most other GPCRs the activating ligand, retinal, is covalently bound to Lys 296 on transmembrane helix 7 (H7). As noted above, upon absorption of a photon retinal undergoes an 11-cis alltrans isomerization [2], initiating a photocascade culminating with activation of the G protein, transducin (Fig. 1). This photocascade consists of a spectroscopically well-defined series of intermediates: dark state bathorhodopsin lumirhodopsin Meta I Meta II [22]. It is important to appreciate that rhodopsin is unable to fully activate transducin until reaching the Meta II state. Because the equilibrium between the Meta I and Meta II states strongly depends on conditions such as ph, temperature, humidity, and membrane bilayer composition, one can physically trap these intermediates by Fig. 1. Photocascade of rhodopsin underlies multiscale activation mechanism involving conformational changes with different time scales. Upon light absorption 11-cis to trans isomerization converts retinal from an inverse agonist to an agonist. Rhodopsin intermediates are designated as photorhodopsin, bathorhodopsin, blue-shifted intermediate (BSI), lumirhodopsin, metarhodopsin I, and metarhodopsin II. The Meta II state comprises an ensemble of Meta II a, Meta II b, and Meta II b H + substates (activated ensemble mechanism). Note that various photoproducts can be trapped under different conditions of temperature, ph, or lipid composition.
3 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) changing thermodynamic state variables. This stabilization provides a means for studying the preactivated Meta I state as well as the fully activated Meta II state of rhodopsin. Techniques which have been used successfully in this regard include magnetic resonance methods, such as spin-label electron paramagnetic resonance (EPR) [23], solidstate 13 C NMR [20,24 30], and 2 H NMR spectroscopy [19,31 35] Solid-state 2 H NMR line shapes are related to dynamical structure of retinal ligand Solid-state 2 H NMR spectroscopy uses site-directed 2 H-labeling to investigate the local structure and dynamics of both proteins [19,32,33] and lipid bilayers [36 41]. In the case of membrane proteins, 2 HNMR provides orientational restraints (bond orientations), and when coupled with interatomic distance restraints it gives structural data analogous to solution NMR [19]. Residual quadrupolar couplings correspond directly to segmental order parameters, and can be measured along with dynamical parameters (nuclear spin relaxation rates). A more detailed explanation of this method can be found elsewhere [34]. In our research, retinal is synthesized with 2 H labels on a specific methyl group (C1-, C5-, C9-, or C13-methyl), reincorporated into the binding pocket of opsin, and recombined with synthetic phospholipids to produce a nativelike environment [19,32,33]. Isotopic labeling of each methyl group allows one to examine various regions of the binding pocket, and thus to delineate specific interactions during the activation of rhodopsin Molecular simulations further characterize structural dynamics and fluctuations of membrane proteins An important aspect is that 2 H NMR line shape analysis can be used as a means for validating computational simulations, which in turn can extend and aid the interpretation of spectroscopic observables. A representative example of this approach to spectral line shape analysis can be found in Fig. 2. Experimental 2 H NMR spectra of the retinal methyl groups are related to bond orientations with respect to the membrane normal. The theoretical 2 H NMR spectra of aligned samples can be calculated using C\C 2 H 3 bond orientational distributions (relative to the membrane) as obtained from molecular dynamics (MD) simulations, thus allowing testing of different molecular models [32,33,42]. Subsequently it is possible to directly compare the simulated spectra to the experimental spectra, providing an excellent means to use simulation in a complementary way with experimental observables. 3. Molecular dynamics calculations allow testing of specific counterion models for rhodopsin One of the major shortcomings of the 2 H NMR approach as applied to the retinal cofactor [35] is that it does not provide structural details for the entire protein. However, MD simulations can be used in combination with 2 H NMR to provide knowledge of the dynamics of the entire proteolipid system. Since rhodopsin was the first published GPCR crystal structure, it has been the most extensively studied by computation. Namely, it has been used to examine various aspects of GPCR activation including overall structural changes [43,44], lipid effects [45,46], oligomerization [47,48], and lipid protein interactions [49,50], among others. Among the hotly debated aspects of the rhodopsin activation mechanism has been the deprotonation event of the protonated Schiff base (PSB) in the activated Meta II state (Fig. 2). In the originally proposed counterion-switch model [51], rhodopsin has a neutral binding pocket, where Glu 113 is deprotonated and acts as a counterion to the PSB. The protonated Glu 181 residue acts as a Brønsted acid and undergoes proton transfer to the conjugate base Glu 113 upon isomerization of the retinal ligand, thus accommodating deprotonation of the Schiff base in the active Meta II state (Fig. 2a). Based on more recent data, an alternative complex-counterion model has also been developed [52], Fig. 2. Molecular dynamics (MD) simulations of rhodopsin allow testing of models for retinal Schiff base counterion based on bond angle distributions from solid-state 2 H NMR spectroscopy. (a) Scheme 1 shows a counterion-switch model where the counterion of the retinylidene protonated Schiff base (PSB) changes from Glu 113 to Glu 181 upon photoisomerization. (b) Scheme 2 is an alternative complex-counterion model in which both Glu 113 and Glu 181 are counterions to the retinal PSB of rhodopsin. (c, d) Results of MD simulations are compared to experimental 2 H NMR spectra for retinal 2 H-labeled at the C5-, C9-, or C13-methyl groups in the Meta I state at 100 C in aligned membranes. The orientation of the average membrane normal to the magnetic field was θ=0, 45, and 90. Theoretical 2 H NMR spectra (continuous lines) are shown compared to experimental results for (c) the counterion-switch model and (d) the complex-counterion model. Note the complex-counterion model agrees best with experimental solid-state 2 H NMR spectra. Figure adapted with permission from Ref. [55]. assigning a negative charge to the binding pocket. Here both Glu 113 and Glu 181 are deprotonated, with Glu 113 accepting the proton from the PSB (Fig. 2b). Several computational studies have examined this problem using different approaches, including Poisson Boltzmann methods, classical MD simulations, and hybrid quantum mechanics/molecular mechanics (QM/MM) approaches [53,54], with results tending to support the complex-counterion mechanism. To definitely answer this question, we set out to simulate both counterion models using long time-scale classical MD simulations, and then compared the results to experimental 2 H NMR data (Fig. 2). Each system consisted of a rhodopsin molecule embedded in a ternary lipid bilayer (1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine
4 244 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) (SDPC)/1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphoethanolamine (SDPE)/cholesterol), and was simulated for over 1.5 μs following retinal isomerization to correspond to the Meta I state [55]. A case clearly supporting the complex-counterion model emerged from this line of research. Theoretical 2 H NMR spectra for the C5-, C9-, and C13-methyl groups of retinal were calculated from the final 500 ns of each MD trajectory, representative of the preactivated Meta I state (Figs. 2c,d). Calculations involved weighting the spectra with bond orientational distributions, and using a Monte Carlo method to compute results for a static uniaxial distribution with mosaic spread. Furthermore, calculations were performed for three sample alignments (θ=0, 45, and 90 ) with respect to the magnetic field B 0. Using different sample tilt angles allows for a more accurate and complete interpretation of the spectral data. Another factor to consider is the mosaic spread of the NMR samples: its contribution to the line shape is minimized close to 0, meaning a more accurate result is produced. Clearly, there is a discrepancy between the experimental spectra of the Meta I state and the counterion-switch simulation (Fig. 2c), especially for the C9- and C13-methyl groups. This poor agreement occurs most evidently for the 0 tilt. In contrast, the complex-counterion model shows excellent agreement with the 2 H NMR data for Meta I (Fig. 2d), supporting this mechanism under which PSB deprotonation occurs subsequently in the active Meta II state [55]. By directly computing experimental observables, in this case the solid-state 2 H NMR spectral line shape, one can begin to appreciate the importance and complementary nature of experimental and computational techniques. Successfully duplicating experimental data also allows one to interpret other aspects of simulation results, such as dihedral fluctuations in the retinal polyene chain. Lastly, one is able to draw parallels between rhodopsin and other GPCRs that are activated by diffusible ligands. 4. Solid-state 2 H NMR spectroscopy investigates the dynamics of cofactors bound to membrane proteins We have emphasized above that solid-state NMR spectroscopy is an invaluable tool to study the structure and dynamics of membrane proteins, inasmuch as the proteolipid system cannot be dissolved and investigated by solution NMR. In particular, solid-state NMR line shapes entail structural restraints involving residual dipolar couplings (RDCs) or residual quadrupolar couplings (RQCs) (the former are related to internuclear distances and the latter to bond angles). Moreover, NMR relaxation provides unique information on molecular motions, and reveals intra- and intermolecular interactions through dynamical parameters. Specific interactions are further investigated by combining relaxation data with structural data Relaxation in 2 H NMR spectroscopy provides key insights into retinal mobility The complementary nature of solid-state 2 H NMR spectral line shape and relaxation studies is illustrated in Fig. 3, which shows experimental results for rhodopsin in the dark state with retinal 2 H-labeled at the C5-, C9-, or C13-methyl groups. Data are included for rhodopsin in unoriented recombinant membranes comprising 1-palmitoyl-2- oleoyl-sn-glycero-3-phosphocholine (POPC). The RQCs obtained from the powder-type spectra (Figs. 3a c) show that the amplitude of the fluctuations of the retinal methyl groups describe the dynamic structure of retinal embedded in the binding pocket of rhodopsin. The average structure of retinal is characterized by the order parameter S C3 = 1 2 3cos2 β IM 1 of the various methyl groups. Here β IM is the angle between the instantaneous C\C 2 H 3 axis direction (I) and the average methyl group (M) orientation. The residual quadrupolar splittings for the C5-, C9-, and C13-methyl groups indicate motional narrowing, due to rapid spinning of the methyl groups about their C 3 axes. Comparing the residual (motionally averaged) quadrupolar Fig. 3. Site-specific 2 H NMR relaxation manifests functional dynamics of retinal methyl groups within the ligand-binding cavity of rhodopsin. (a c) Powder-type 2 H NMR spectra of rhodopsin containing 11-cis-retinal with 2 H-labeled C5-, C9- or C13-methyl groups in POPC bilayers (1:50 protein to lipid molar ratio). Solid-state 2 H NMR line shapes demonstrate rapid spinning of C\C 2 H 3 groups at temperatures down to at least 160 C. (d, e) Inversion-recovery 2 H NMR spectra for retinylidene C9- and C13- methyl groups of rhodopsin in aligned POPC membranes at 150 C. The tilt angle θ between the average membrane normal and magnetic field is equal to 0. (f) Inversionrecovery plots for C9- and C13-methyl groups at 150 C demonstrate site-specific variations in longitudinal Zeeman (T 1Z ) relaxation times. Figure adapted with permission from Ref. [35]. coupling χ Q with the value of χ Q due to fast spinning of the methyl group yields S C3 0.9, representing an amplitude of the off-axial fluctuations of Δβ IM 15. Note that this motion includes both methyl group fluctuations with respect to the unsaturated polyene, as well as any reorientations of the retinal ligand within the binding pocket of rhodopsin. Because the solid-state 2 H NMR spectra are in the fast motional limit, they do not distinguish among the various possible contributions to the methyl group dynamics. However, site-specific differences in retinal methyl mobility can be delineated from partiallyrelaxed 2 H NMR spectra, as well as from the inversion recovery curves used to establish the spin lattice relaxation times (Figs. 3d f). To further elucidate how the retinal dynamics are implicated in rhodopsin function, we measured 2 H relaxation times (Zeeman, T 1Z, and quadrupolar order, T 1Q ) for aligned or nonaligned (powder-type) samples (Fig. 4). Data were acquired in the temperature range from 30 to 160 C (close to values used for X-ray [56] and spin-label EPR [23] studies), which is well below the melting temperature (T M ) of the lipid bilayer. Working below T M freezes out the rotational diffusion and large-scale (collective) motions of rhodopsin, and allows us to focus on the internal (local) dynamics of retinal. This approach constitutes a rather different strategy from that typically employed in solution NMR of proteins [57]. Moreover, the lipid dependence of rhodopsin intermediates [58 61] was used to quantitatively trap the Meta I and Meta II photoproducts in frozen aligned rhodopsin films. Now in the partially-relaxed 2 H NMR spectra (Figs. 3d,e), the null point for inversion recovery is related to T 1Z by t=t 1Z ln2. As a result, by comparing the data for the C9- and C13-methyl groups, we can immediately see that the spin lattice relaxation rates differ by more
5 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) below approximately 120 C. (In Ref. [35] the value of 1/ω 13 ns is a misprint and should be corrected as indicated here.) The 2 H NMR relaxation results are summarized in Fig. 4, which shows Arrhenius-type plots of the T 1Z relaxation times versus inverse temperature in the dark, Meta I, and Meta II states, respectively. A minimum is expected in plots of ln T 1Z versus 1/T, since matching the power spectrum (i.e. spectral density) of the fluctuations to the energy level transitions is optimal when the correlation time is near the inverse resonance (Larmor) frequency. Consequently, shifting the minimum to lower temperature corresponds to faster motions and vice versa. Quantitative evaluation assumes transition state theory leading to Aexp( E a /RT) for the temperature dependence. Here the preexponential factor A represents the intrinsic rate of rotation in the absence of a barrier, i.e. E a bbrt corresponding to infinite temperature, where E a is the barrier height (activation energy) for the methyl rotation [34]. The width of the minimum in the ln T 1Z versus 1/T plot is inversely related to the E a activation barrier for the methyl dynamics. For temperatures significantly above the T 1Z minimum, the plot essentially corresponds to an Arrhenius-type law, where the slope is proportional to E a (see below) Methyl groups probe site-specific local interactions in the inactive dark state Fig. 4. Solid-state 2 H NMR relaxation uncovers striking changes in retinal mobility and interactions during rhodopsin activation. (a c) Spin lattice (T 1Z ) relaxation times of retinylidene methyl groups were measured in (a) the dark state, (b) the Meta I state, and (c) the Meta II state at temperatures ranging from 30 to 160 C. Methyl dynamics were analyzed using either an axial 3-fold jump model, or a continuous diffusion model with coefficients D (for methyl group rotation about the 3-fold (C 3 ) symmetry axis) and D (for methyl axis reorientation). In panels (a c) an axial 3-fold jump model or a continuous diffusion model with D =0 are represented by solid lines; the dashed lines are for a continuous diffusion model having restricted off-axial diffusion with D =D. In panel (b) theoretical fits for the C5-methyl group in Meta I assume either two conformers with different bond orientations and axial diffusion coefficients (solid line), or two different rotational diffusion constants (D D ) (dashed line). Figure adapted with permission from Ref. [35]. than an order of magnitude (Fig. 3f). This variation is caused by sitespecific differences in the intramolecular dynamics of the retinylidene chromophore. Since the RQCs for the retinal methyl groups are nearly equivalent, the relaxation differences must arise from the methyl spinning rates, as opposed to the amplitude of the off-axial motions (S C3 ). What becomes clear is that distinct variations exist in the internal mobility of the retinal methyl groups in the dark state [35]. The T 1Z relaxation times for the dark, Meta I, and Meta II states increase with temperature for each of the methyl groups, except for the C5- methyl at extremely low temperatures. This means that motional correlation times are less than 1/ω 2 ns except for the C5-methyl group, whose experimental rotational correlation time exceeded 2 ns Keeping in mind that 2 H NMR relaxation can provide us with sitespecific variations in retinal mobility, we are now able to examine the behavior of 2 H-labeled retinal bound to rhodopsin in the inactive dark state. Fig. 4a shows that ln T 1Z has a temperature dependence that is nearly linear for the C9- and C13-methyl groups, although with quite different slopes. By contrast, a curve with a minimum near 120 C is observed for the C5-methyl group. Hence, the individual methyl groups possess markedly different dynamical and relaxation behavior, with three nonequivalent environments in the dark state. Smaller T 1Z values and rotational diffusion coefficients are correlated to larger activation energy (E a ) values and vice versa, which may indicate a role of nonbonded interactions. Other differences originate from the preexponential factor: for instance, the T 1Z value for the C9-methyl at the higher temperature (Fig. 4a) is due to a relatively small preexponential factor, as opposed to the activation energy. These differences in the rotational dynamics and activation energies of the methyl groups observed in the dark state (and upon transition to the Meta I and Meta II states, see below) can be attributed to a combination of intra- and intermolecular interactions of retinal within the chromophore binding pocket. The dynamical parameters obtained from Fig. 4a indicate that at low temperatures (b0 C),theβ-ionone ring (C5-methyl group) is the least mobile in the dark state. In addition, the C13-methyl group of the polyene chain has intermediate mobility, and the C9-methyl in the middle of the retinylidene cofactor is most mobile. Despite different temperature dependencies and E a values, the mobility of the C9- and C13-methyl groups is almost the same near physiological temperature. By contrast, the T 1Z time for the C5-methyl and correspondingly its rate of motion are 3-fold smaller (see below). Another interesting observation from the relaxation data in the dark state (Fig. 4a) is that the C9-methyl activation energy is significantly less than for the C5- and C13-methyl groups. A possible explanation is that the C9-methyl group dynamics are strongly determined by 1 6 interactions (C9-methyl hydrogen to H7 and C9-methyl hydrogen to H11), which can be approximated by 3-fold rotational potentials. Although these potentials may be the same order of magnitude as 1 7 interactions (C13-methyl hydrogen to H10), they are shifted with respect to each other by 60, thus producing a shallow potential and activation energy. This hypothesis is further supported by quantum mechanics (QM) calculations (see below). As a result, we propose that nonbonded 1 6 intramolecular interactions are the primary reason for the low activation energy E a and also the relatively small preexponential factor for 3-fold jumps of the C9-methyl group. According to the X-ray
6 246 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) structure [5] and MD simulations [62], the dynamics of both the C9- and C13-methyl groups are most likely relatively unaffected by intermolecular interactions in the dark state; on average, these groups are 4 A away from amino acid side chains. It follows that the activation energy for the C13-methyl rotational diffusion is mainly determined by the nonbonded 1 7 interaction with hydrogen H10. The high activation energy for the C5-methyl group in the dark state stems from interactions with hydrogen H8 of the retinal ligand and the Glu 122 side chain. This proposal is consistent with a Fourier transform infrared (FTIR) study of desmethyl and acyclic analogs of retinal bound to rhodopsin [63],and is supported by close proximity of the C5-methyl group to Glu 122 in the X- ray crystal structure [5] Changes in retinylidene methyl dynamics occur in the preactive Meta I state and signaling Meta II state It is striking that compared to the dark state (Fig. 4a) significant differences occur in the relaxation behavior for the C5- and C9-methyl groups in the Meta I and Meta II states (Figs. 4b, c). The E a for the crucial C9-methyl group (3.7 kj mol 1 ) is markedly increased (by a factor of 2) in Meta I, up to a level similar to the C13-methyl. Concurrently, the E a value for the C13-methyl group in the Meta I state (4.7 kj mol 1 )is reduced by about one-third, making its T 1Z values nearly identical to those of the C9-methyl group. The decreased C9-methyl mobility following isomerization is most likely related to the more crowded environment, which results from the polyene chain all-trans conformation and nonbonded interactions within the binding pocket of rhodopsin. The lack of high-resolution Meta I structural data makes it difficult to identify which specific amino acids affect the C9-methyl mobility, which could potentially explain changes in its activation energy. As mentioned previously, the reduction of the activation energy for the C13-methyl group (from 7.3 to 4.7 kj mol 1 )makessense, because the 11-cis all-trans isomerization places the C9- and C13- methyl groups on the same side of the molecule (Fig. 4). Hence they undergo similar intramolecular interactions, and experience a similar mobility. Interestingly, changes are also observed in the T 1Z values for the C5-methyl group in Fig. 4, with an emergence of a second temperature-dependent minimum in the Meta I state, and subsequent reversion to a single minimum in the Meta II state. The significance of the T 1Z minimum is that the relaxation time at this temperature unequivocally defines the effective correlation time, as described below. Evidently the C5-methyl relaxation data in the Meta I state do not conform to the framework of a 3-fold jump model with a single E a value. The two minima may represent two distinct types of motion, e.g. methyl group rotation and β-ionone ring reorientation. An alternative interpretation would be the presence of two 6-s-cis conformers of the retinal, with positive and negative torsion angles for the C6\C7 bond [62]. Regardless of the model used to describe the C5-methyl dynamics, its behavior remains similar in both the dark and Meta I states. The rotational dynamics are dramatically slower, and the activation energy E a for the C5-methyl (N10 kj mol 1 ) is much larger compared to the C9- and C13-methyl groups (Fig. 4). Furthermore, the minimum of the C5-methyl T 1Z temperature dependence occurs at an elevated temperature versus the C9- and C13-methyls. This essentially means that the C5-methyl rotation is restricted in the dark and Meta I states. It has been proposed recently that the β-ionone ring and, in particular, the C5-methyl disrupts the hydrogen bonding network around Glu 122 in the Meta I Meta II transition, due to retinal motion toward helix H5 [63,64]. This movement of the β-ionone ring has been regarded as an essential part of the rhodopsin activation mechanism. In this respect, the existence of two minima in the T 1Z temperature dependence may indicate a more restricted environment for C5- methyl group rotation in the Meta I state versus the dark state. Another important finding of this work is that relatively small changes occur in the retinal dynamics upon the transition from the Meta I to the Meta II state. The E a values and rotational diffusion constants increase slightly for the C9-methyl group (from 3.7 to 4.3 kj mol 1 ) and decrease for the C13-methyl (from 4.7 to 2.8 kj mol 1 ) (Fig. 4c). The E a value for the C13-methyl group progressively decreases from the dark state to Meta I, and then to the Meta II state. This is the result of a stepwise decrease of the preexponential factor from state to state during the activation process of rhodopsin. For the C5-methyl group, the minimum observed in the T 1Z temperature dependence shifts to lower temperature on going from the dark to the Meta II state (E a =10.3 kj mol 1 ). This position approximately corresponds to the average of the two minima in the Meta I state. The minimum also broadens, which would indicate a superposition of two nearby minima. Most striking, the overall T 1Z temperature dependence for the C5-methyl in the dark, Meta I, and Meta II states remains the same (differences in E a barely exceed 20% for a given model of the rotational dynamics), implying a stable environment for the β-ionone ring. A similar conclusion has been reached based on 13 C chemical shift data [64,65]. 5. Analytical relaxation theory and numerical dynamics simulations characterize the energy landscape of rhodopsin activation As a rule there are two approaches to the quantitative analysis of NMR relaxation data, such as those summarized above. First, one can employ the general theory of nuclear spin relaxation to interpret the relaxation times in terms of the power spectral densities and correlation functions, using simplified models for the rotational dynamics in analytical closed form. Alternatively, we can forgo an exact mathematical description, and interpret the experimental relaxation times using numerical molecular dynamics simulations. Here the temporal evolution of the molecular positions and orientations under the influences of a force field encapsulates knowledge of intra-and intermolecular interactions pertinent to the systems under investigation. Both approaches are valid, because their characteristic strengths and weaknesses are different and tend to be complementary Solid-state 2 H NMR highlights the functional motions of retinal methyl groups According to theory, the 2 H NMR relaxation rates are determined by the power spectrum (spectral density) of the thermal fluctuations of the quadrupolar coupling tensor near the resonance frequency [66]. They are given by R 1Z =1= T 1Z = 3 4 π2 χ 2 Q ½ J 1 ðω 0 Þ +4J 2 ð2ω 0 ÞŠ ð1þ and R 1Q =1= T 1Q = 9 4 π2 χ 2 Q J 1 ðω 0 Þ; ð2þ where χ Q is the static quadrupolar coupling constant, and ω 0 is the resonance (Larmor) frequency. Here J m (mω 0 ) denotes the spectral density (m=1,2), which is a function of the molecular geometry, the mean-square amplitudes of the motion, and the rotational correlation times (τ r ). The latter in turn are related to the rate constant for 3-fold rotations (k) about the C 3 methyl axis by 1/τ r =3 k (r=1) [39]. An alternative approach is to use a continuous rotational diffusion model [66]. Within a strong collision approximation one obtains 1/ t r 6D +(D D )r 2, where the rotational diffusion coefficients D and D characterize methyl group axial rotation and off-axial fluctuations, correspondingly. The temperature dependence of either
7 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) the jump rate or rotational diffusion constants assumes an Arrhenius law, as described above (cf. [34] for detailed description). In principle, T 1Z measurements at different orientations for aligned samples or at different frequencies for powder-type samples allow one to identify specific motional mechanisms. An axial jump model predicts a 20% difference in the methyl relaxation depending on orientation [67], whereas for an axial diffusion model the methyl relaxation is angular independent [39,66]. In practice, for the case of powder-type samples our T 1Z measurements were limited to the θ=90 orientation of the methyl axis to the main magnetic field B 0, due to the relatively low signal-to-noise ratio. Aligned samples were investigated at the θ=0 orientation relative to B 0. There was no detectable difference in the relaxation times for various orientations. Consequently, the experimental error of the T 1Z measurements (about 10%) did not allow us to make unambiguous conclusions about the reorientation mechanism. Both a 3-fold jump [67] and a continuous diffusion [66] model were used to fit the temperature dependence of the T 1Z and T 1Q relaxation times (Fig. 4). For the dark, Meta I, and Meta II states of rhodopsin, the T 1Z data could be fit with a 0:1 D D ratio to within experimental error. However, simultaneous fitting of the T 1Z and T 1Q temperature dependencies could only be performed for D =0 (results not shown), meaning that methyl group off-axial fluctuations are much slower than spinning about its 3-fold symmetry axis. The preexponential factor and E a values obtained using transition state theory (see above) vary according to the model, yet the qualitative picture remains the same Molecular dynamics simulations expand the 2 H NMR results to longer time scales For rhodopsin proteolipid membranes, molecular dynamics (MD) simulations offer the ability to study both smaller (ns) and larger (μs) time-scale dynamics for the entire system, providing the basis for extending insights from 2 H NMR spectroscopy. Many system aspects can be examined: protein lipid interactions [49,50], behavior of the lipid bilayer (i.e., bilayer thickness, order parameters, lipid diffusion, etc.) [68], protein conformational changes leading to the activated state [43,62,69], and agonist binding effects [17]. However, two limitations need to be considered when modeling GPCRs: (i) the inaccessibility of longer time scales (rhodopsin reaches Meta II in milliseconds) and (ii) the use of an empirically derived molecular mechanics (MM) force field to approximate atomic interactions. Clearly, molecular mechanics force fields are intended to correctly reproduce experimental phenomena. Because they are empirically derived, however, MM force fields are continually improved as experimental data and theory become more robust. This applies especially to retinal, since our latest 2 H NMR research clearly shows that each methyl group has site-specific intra- and intermolecular interactions (Fig. 4). Addressing this issue is biologically important, because experiment has shown that the retinal methyl groups are essential to rhodopsin activation [63,70]. As a result, we have used ab initio quantum mechanics calculations to examine the applicability of MM force fields for describing the rotational behavior of methyl groups in appropriate retinal model compounds. 6. Quantum mechanical calculations address methyl dihedral terms in the retinal force field In molecular mechanics, the force field encapsulates all information about the intra- and intermolecular forces that act upon the molecular system of interest. Analysis of bond angle bending, stretching, and torsions as well as electronic structure requires knowledge of intramolecular potential energies that derive ultimately from either quantum mechanical calculations or experimental data. For molecules such as retinal, a current drawback is that they are too large to be accurately computed with high-level quantum theory. Calculations must be restricted to comparatively small fragments that approximate the molecular orbitals. Here, recent advances in computing power have motivated us to reexamine the torsional potential energy terms for retinal currently used with CHARMM and other molecular mechanics programs. Our general aim is to more accurately model experimental 2 H NMR relaxation data for rhodopsin in different functional states Larger polyene compounds are necessary to model the retinylidene methyl groups Comparison of our 2 H NMR results with functional data [63,70] has demonstrated that each of the retinylidene methyl groups contributes individually and uniquely to rhodopsin activation. This finding is in agreement with several FTIR studies conducted on retinal desmethyl analogs, which show that removal of any methyl group prevents rhodopsin from becoming fully activated [63,70]. In contrast, previous MD simulations of the rhodopsin proteolipid system have neglected the retinal methyl torsion force constants by setting them to zero [71] (i.e. the potential energy was assumed to be determined solely by nonbonded terms). The discrepancy between experiment and simulation led us to a closer examination of the methyl dihedral parameters in the retinal force field (Fig. 5). The majority of molecular mechanics force fields draw from a combination of experimental and theoretical data to create a set of parameters that are used to solve Newton's equations of motion. In the case of retinal methyl dihedrals, quantum mechanical calculations (in this case, MP2 level of theory with a 6 31 G** basis set) are employed to formulate the force field parameters. When examining retinal methyl rotations, steric interactions with the polyene chain figure prominently. Specifically, methyl groups in an all-trans polyene chain experience 1 6 interactions Fig. 5. Quantum mechanical theory reveals the influences of steric interactions on torsional potential energy surfaces underlying retinal methyl rotation. (a) Model compounds corresponding to fragments of retinal used for QM calculations. (b) Comparison of QM (circles) and MM (lines) energies plotted against methyl torsion angle (ϕ) for 2-methyl-butadiene (2MBD) (compound 1) and 3-methyl-hexatriene (3MHT) (compound 2). (c) Plot of 1 6 distance in 3MHT between methyl hydrogen and C1 vinyl hydrogen (closed) and C5 vinyl hydrogen (open). (d) Plot of QM energy as a function of methyl torsion angle in 2MBD (compound 1), propane (compound 3), and propene (compound 4). (e) Activation energies (E a ) for C5-, C9-, and C13-methyl groups obtained from 2 H NMR data for the dark, Meta I, and Meta II states of rhodopsin versus standard MM dihedral energy barrier in CHARMM force field. Figure adapted with permission from Ref. [72].
8 248 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) between the methyl hydrogens and those hydrogen atoms bound to neighboring vinyl groups. For 2-methyl-butadiene (2MBD) (compound 1), the largest polyene fragment currently used in the retinal force field distributed with CHARMM, only one such 1 6 interaction is present (Fig. 5a). By extending the polyene chain by two carbons, creating the fragment 3-methyl-hexatriene (3MHT) (compound 2), we now possess a methyl group with two 1 6 interactions. Adding the extra 1 6 interaction reduces the 3MHT methyl energy barrier in comparison to 2MBD by raising the energies of both the minimum energy conformation and the eclipsed state (Fig. 5b). For 3MHT the 1 6 interhydrogen distances as a function of torsion angle also demonstrate the effect of dual 1 6 interactions (Fig. 5c), showing that 3MHT is nearly symmetrical over the course of a methyl rotation. The distance between the methyl group and the C1 vinyl hydrogen (direction of single bond) is slightly higher at its point of closest approach than the corresponding distance to the C5 vinyl hydrogen (direction of double bond). Hence, the steric repulsion at the torsional energy maximum is slightly greater than the repulsion at the minimum. These results suggest that 3MHT is the smallest possible fragment to accurately model the retinal C9-methyl group. The standard parametrization procedure is to fit the QM results to the torsional energy term in the molecular mechanics force field. For the MM package CHARMM, this contribution takes the form U(ϕ)=K ϕ [1+cos(nϕ+δ)]. By retaining the current force field terms and using a new dihedral potential one can perform MM simulations of the polyene fragments 2MBD and 3MHT and compare the results to QM calculations (Fig. 5b). It can be seen that the new dihedral parameters accurately reproduce the QM rotational profile, and give a marked improvement in the retinal force field Hyperconjugation extends to methyl groups of the retinal chromophore implicated in rhodopsin activation One should note that polyene chains have unique covalent bond characteristics, due to their single bond-double bond alternation. This bond alternation creates hyperconjugation effects, which essentially imparts shared single and double bond properties to each bond in the polyene chain. In order to fully understand the steric effects on 2MBD and 3MHT, it was necessary to examine two smaller model compounds, propane (compound 3) and propene (compound 4), Fig. 5a. The QM results show that the methyl rotation barrier of a fully saturated hydrocarbon chain (propane) is not necessarily equivalent to the same barrier of an unsaturated chain (propene), shifting both the barrier location and height (Fig. 5d). In addition, the methyl rotational barrier of propene is almost half (8 kj mol 1 )thatofpropane(15kjmol 1 ), which stems from different orbital interactions in these two molecules. The most stable methyl conformation in propene is at ϕ=0,wherethe methyl C\H bond is cis to the vinyl C_C bond. The reason for this stability is hyperconjugation: the other two methyl hydrogens are connected by a vector perpendicular to the C_C bond axis, which involves stabilization through interaction of the filled π-like CH 2 orbitals of the methyl group and the unoccupied π* molecularorbitalofthe double bond. Propane lacks these hyperconjugative interactions, with the most stable conformer being the staggered orientation. In the context of 2MBD and 3MHT, the dramatically increased rotational barrier of 2MBD compared to propene (Fig. 5d) is due to steric repulsion between the methyl and vinyl hydrogen atoms, which is maximized for 2MBD when both hydrogen atoms are eclipsed. Because of the 3-fold symmetry of the methyl group, when one of the methyl hydrogens is eclipsed (ϕ=180 ), the other methyl hydrogens lie at 60. Thus maximum steric repulsion coincides with the highest energy state arising from orbital interactions, with the net effect of raising the energy barrier by 50% in 2MBD compared to propene (Fig. 5d). This steric effect would predict a rotational barrier for 3MHT that is slightly higher than propene, while QM calculations give the opposite (Figs. 5b,d). Most likely, this discrepancy stems from the altered π orbital characteristics when increasing polyene chain unsaturation from one double bond (propene) to three (3MHT). Yet another hyperconjugation effect, which was not considered in our original 2 H NMR analysis [21], butwhich came to light in our theoretical investigations [72], was that protonation of the Schiff base yields a lowering of the C13-methyl rotational energy barrier. This finding is significant due to the fact that Meta II involves deprotonation of the PSB, which in an MD simulation would involve physically removing the proton from the terminus of retinal. Currently these time scales are not attainable but with advances in computational power, and the availability of new crystal structures in putative Meta II states [9,10], they will be soon. These QM calculations, guided by insight from experiment, ultimately led to improvements in the retinal force field [72]. Direct comparisons between the 2 H NMR data and QM results cannot be made, since the experimental analysis is made in the protein environment and the theoretical analysis is made in vacuo. However, it is clear from both approaches that intramolecular steric clashes with hydrogens play a large role in determining the activation energy profiles for retinal methyl rotations. The effort put into the analysis of such a small molecular component of the rhodopsin system, viz. retinal, is important due to the crucial role it plays in rhodopsin activation. An improved MM force field will thus allow future investigations of effective barrier heights in the full protein environment, with subsequent quantitative comparisons and detailed interpretations of experimental observations. 7. Putting it all together: the rhodopsin activation mechanism Based on the structural [11,12,20,23,65,73 76] and dynamical [34,35,45,46,49,55,62] data, a multiscale activation mechanism has been proposed, whereby retinal isomerization initiates collective helix fluctuations in the Meta I Meta II equilibrium with multiple activated substates. The proposed structural changes upon rhodopsin activation are shown in Fig. 6. In accord with our MD simulations (ns time scale) [68], helices H1 H4 constitute the transmembrane core of the receptor, whereas the cytoplasmic loops evince significant motions [68]. Rhodopsin activation is due to movements of helices H5 and H6 together with H7. In an expanded view of the photopathway (Fig. 1), the early photointermediates (photorhodopsin, bathorhodopsin, BSI, lumirhodopsin) are followed by the major reactions [77]: Meta I Meta II a Meta II b +H 3 O + Meta II b H + ; here Meta II b and Meta II b H + are the active conformational substates. The entropy-stabilized Meta II b substate is in equilibrium with the Meta II b H + substate, whose formation is enthalpically downhill [78]. Using 2 H NMR, we are able to show how the fast time-scale motions of retinal (ps) correspond to the helical movements associated with rhodopsin activation [35] Retinal methyl groups assist in driving multiscale rhodopsin activation In the dark state, the low activation barrier of the C9-methyl group (E a =1.4 kj mol 1 ) indicates the absence of steric clashes within the ligand-binding cavity of rhodopsin (see Figs. 4a, 5e). The small increase in E a value for the C9-methyl in the Meta I state (from 1.4 to 3.7 kj mol 1 ) signifies the local environment changes following isomerization. Consequently, the 2 H NMR data suggest that leading up to and including the Meta I state [19], the C9-methyl acts effectively as a hinge point for retinal isomerization, causing reorientation of the C13-methyl and the C_NH + groups. Rotation of the C13-methyl group and the PSB occurs toward the C9-methyl group rather than vice versa. This conclusion is supported by the crystal structures of batho- and lumirhodopsin [8]. The reduced C13-methyl activation energy in Meta I suggests the second extracellular (E2) loop and the retinal ligand move apart. Hence, the C13- methyl rotation displaces the β 4 strand of the E2 loop toward the extracellular side, disrupting the hydrogen-bonding network connecting the extracellular ends of helices H4, H5, and H6 [20]. At the opposite end of retinal, the β-ionone ring is displaced toward the H3 H5 helical
9 B. Mertz et al. / Biochimica et Biophysica Acta 1818 (2012) cytoplasmic side. Receptor activation culminates with protonation of Glu 134 of the E(D)RY motif from the aqueous medium [78] (Fig. 6a). Our proposed mechanism [35] varies somewhat from the recent X-ray structure attributed to Meta II [9]. The Meta II X-ray structure was obtained from opsin crystals soaked with all-trans-retinal, showing the ligand is rotated by 180 about its long axis compared to rhodopsin in the dark state (Fig. 6d) [9]. The retinal location is in agreement with most distance constraints obtained from dipolarassisted rotational-resonance NMR spectroscopy [30]. However,the authors mention that some larger distance deviations may indicate the existence of different Meta II substates observed in the NMR and X-ray experiments. Moreover, the X-ray crystal structure of the E113Q/N2C/D282C triple mutant in the active state shows an orientation of the retinal ligand that is opposite from the retinalsoaked opsin structure [10]. It is possible that the rotated orientation of retinal in Meta II obtained by opsin regeneration may stem from the crystal soaking protocol. At this point, it cannot be concluded that a singular Meta II structure has been unequivocally established. On the other hand, various experimental structures of the active state support the possibility of multiple Meta II substates being present in the Meta I Meta II equilibrium [35,75,78]. Fig. 6. Rhodopsin activation is triggered by retinal isomerization and stabilized by rearrangement of ionic locks. (a) Structure of rhodopsin (PDB accession number 2X72) shows hydrogen-bonding networks and ionic locks that rearrange or form after isomerization of the covalently bound retinal ligand. Spheres represent retinal; the sticks indicate residues in specific networks; the smooth van der Waals surface shows a transducin peptide fragment; and the ribbons depict the rhodopsin polypeptide backbone. (b) Rhodopsin in the dark state (PDB 1GZM) is stabilized by an ionic lock due to nonbonded interactions among Glu 134, Arg 135, and Glu 247 residues. (c) Upon retinal isomerization the rearrangement of the ionic lock entails new interactions between Arg 135 Tyr 223 and Glu 247 Lys 231 Thr 251 (rhodopsin in putative Meta II state; PDB 2X72). (d) Comparison of ligand binding pocket for the dark state (PDB 1GZM; green sticks and cyan retinal), the E113Q putative Meta II state (PDB 2X72; purple sticks and yellow retinal), and for retinal-soaked opsin (PDB 3PQR; cyan sticks and orange retinal). Note that all-trans-retinal in the retinal-soaked opsin is rotated by 180 about its long axis either compared to 11-cis-retinal in the dark state or the E113Q mutant. interface, producing increased packing for the C5-methyl as indicated by the relatively high E a value in the Meta I state (Figs. 4b, 5e). Additionally, movement of the β-ionone ring rearranges the hydrogen-bonded network about Glu 122 and His 211 connecting helices H3 and H5 [30]. This releases allosteric constraints on helix H5, allowing it to rotate its cytoplasmic end toward helix H6. Retinal isomerization is thus responsible for transmission of energy against the E2 loop by the C13-methyl group, and toward helices H3 and H5 by the C5-methyl of the β-ionone ring. The ionic lock involving the retinylidene PSB on H7 with its complex counterion due to Glu 113 (H3) and Glu 181 (E2) is broken by internal proton transfer from the PSB to Glu 113 in Meta II [52]. Itfollowsthat destabilization of the ionic lock may occur through the rotation of the C_NH + group due to retinal isomerization Ionic locks hold the key to rhodopsin transducin interactions Notably, retinal isomerization and the displacement of its PSB end away from the E2 loop create a potential steric clash of the middle of the polyene chain with Trp 265 of H6. To avoid steric hindrance, Trp 265 shifts toward helix H5 and in the cytoplasmic direction, enabling the concerted tilt of helix H6 away from the H1 H4 helical core (Fig. 6a). In addition to the PSB, a second ionic lock due to the E(D)RY motif is implicated with rhodopsin activation (Fig. 6b). Displacement of helix H5 occurs concurrently with the tilt of helix H6, bringing Tyr 223 closer to Arg 135 and Lys 231 near to Glu 247 (Fig. 6c) [12]. The charge adduct due to Glu 134 Arg 135 of the E(D)RY sequence in helix H3 with Glu 247 of helix H6 rearranges in Meta II with Arg 135 Tyr 223 and Glu 247 Lys 231 Thr 251 interactions [9] (Fig. 6c). Thus a second ionic lock is broken, exposing transducin (G t ) recognition elements on the 8. Moving forward According to the above picture, receptor activation occurs as a dynamical equilibrium of substates a striking example of how conformational entropy plays a role in GPCR biology and pharmacology. On the other hand, full characterization of the rhodopsin activation mechanism remains incomplete. In particular, the temporal sequence of the conformational changes is not completely established, and it is unclear when or whether long-axis rotation of the retinal takes place in the binding pocket. This confusion stems from the fact that the threedimensional structures of metarhodopsin I and metarhodopsin II were obtained with moderately low resolution, yielding somewhat contradictory results. Consequently, these X-ray structures require further confirmation and refinement. A fundamental question will be to determine if there is a unique activated rhodopsin conformation, or whether receptor function entails a conformational ensemble triggered by photoisomerization of the retinal ligand. In this context, solid-state NMR and computational simulation continue to provide key insights, because they are able to examine the dynamics of rhodopsin activation. But what is exciting to us and other researchers is that development of new NMR methods and improved computational hardware and software expand the horizons for elucidating new aspects of membrane protein function. More than ten years ago the first crystal structure of rhodopsin was published and just over five years ago, microsecondtime scale simulations were considered the state of the art. Currently there are over a dozen solved GPCR structures, and microsecond time scales are becoming much more common in simulations. The synergistic interplay in which experiment and simulation are used to examine rhodopsin activation will continue to answer many of the key questions about GPCR function as the field evolves in the future. Acknowledgements Research support from the U. S. National Institutes of Health (EY to BM and EY and EY to MFB) and the U. S. National Science Foundation (MCB to SEF) is gratefully acknowledged. References [1] D. Filmore, It's a GPCR world, Mod. Drug Discovery 7 (2004) [2] J.W. Lewis, D.S. Kliger, Absorption spectroscopy in studies of visual pigments: spectral and kinetic characterization of intermediates, Methods Enzymol. 315 (2000)
Computational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College
Computational Studies of the Photoreceptor Rhodopsin Scott E. Feller Wabash College Rhodopsin Photocycle Dark-adapted Rhodopsin hn Isomerize retinal Photorhodopsin ~200 fs Bathorhodopsin Meta-II ms timescale
More informationStructural Analysis and Dynamics of Retinal Chromophore in Dark and Meta I States of Rhodopsin from 2 H NMR of Aligned Membranes
doi:10.1016/j.jmb.2007.03.046 J. Mol. Biol. (2007) 372, 50 66 Structural Analysis and Dynamics of Retinal Chromophore in Dark and Meta I States of Rhodopsin from 2 H NMR of Aligned Membranes Andrey V.
More informationRetinal Proteins (Rhodopsins) Vision, Bioenergetics, Phototaxis. Bacteriorhodopsin s Photocycle. Bacteriorhodopsin s Photocycle
Molecular chanisms of Photoactivation and Spectral Tuning in Retinal Proteins Emad Tajkhorshid Theoretical and Computational Biophysics Group Beckman Institute University of Illinois at Urbana-Champaign
More informationI690/B680 Structural Bioinformatics Spring Protein Structure Determination by NMR Spectroscopy
I690/B680 Structural Bioinformatics Spring 2006 Protein Structure Determination by NMR Spectroscopy Suggested Reading (1) Van Holde, Johnson, Ho. Principles of Physical Biochemistry, 2 nd Ed., Prentice
More informationVisual pigments. Neuroscience, Biochemistry Dr. Mamoun Ahram Third year, 2019
Visual pigments Neuroscience, Biochemistry Dr. Mamoun Ahram Third year, 2019 References Webvision: The Organization of the Retina and Visual System (http://www.ncbi.nlm.nih.gov/books/nbk11522/#a 127) The
More informationCalculate a rate given a species concentration change.
Kinetics Define a rate for a given process. Change in concentration of a reagent with time. A rate is always positive, and is usually referred to with only magnitude (i.e. no sign) Reaction rates can be
More informationSolid-State Deuterium NMR Spectroscopy of Rhodopsin
Solid-State Deuterium NMR Spectroscopy of Rhodopsin Suchithranga M.D.C. Perera, Xiaolin Xu, Trivikram R. Molugu, Andrey V. Struts, and Michael F. Brown Abstract Rhodopsin is a prototype for the large Family
More informationTimescales of Protein Dynamics
Timescales of Protein Dynamics From Henzler-Wildman and Kern, Nature 2007 Dynamics from NMR Show spies Amide Nitrogen Spies Report On Conformational Dynamics Amide Hydrogen Transverse Relaxation Ensemble
More informationTimescales of Protein Dynamics
Timescales of Protein Dynamics From Henzler-Wildman and Kern, Nature 2007 Summary of 1D Experiment time domain data Fourier Transform (FT) frequency domain data or Transverse Relaxation Ensemble of Nuclear
More informationSAM Teacher s Guide Protein Partnering and Function
SAM Teacher s Guide Protein Partnering and Function Overview Students explore protein molecules physical and chemical characteristics and learn that these unique characteristics enable other molecules
More informationComputational modeling of G-Protein Coupled Receptors (GPCRs) has recently become
Homology Modeling and Docking of Melatonin Receptors Andrew Kohlway, UMBC Jeffry D. Madura, Duquesne University 6/18/04 INTRODUCTION Computational modeling of G-Protein Coupled Receptors (GPCRs) has recently
More informationCD Basis Set of Spectra that is used is that derived from comparing the spectra of globular proteins whose secondary structures are known from X-ray
CD Basis Set of Spectra that is used is that derived from comparing the spectra of globular proteins whose secondary structures are known from X-ray crystallography An example of the use of CD Modeling
More informationBIOCHEMISTRY GUIDED NOTES - AP BIOLOGY-
BIOCHEMISTRY GUIDED NOTES - AP BIOLOGY- ELEMENTS AND COMPOUNDS - anything that has mass and takes up space. - cannot be broken down to other substances. - substance containing two or more different elements
More informationA molecular-dynamics study of the rhodopsin chromophore using ultrasoft pseudopotentials
Progress of Theoretical Physics Supplement No. 138, 2000 107 A molecular-dynamics study of the rhodopsin chromophore using ultrasoft pseudopotentials Minoru Sugihara, 1 ) Peter Entel, 1 Hendrik Meyer,
More informationThree-layer ONIOM studies of the dark state of Rhodopsin: the protonation state of Glu181
Accepted Manuscript Three-layer ONIOM studies of the dark state of Rhodopsin: the protonation state of Glu181 Katherine. F. Hall, Thom Vreven, Michael J. Frisch, Michael J. Bearpark PII: S0022-2836(08)00987-X
More informationBiophysical Chemistry: NMR Spectroscopy
Relaxation & Multidimensional Spectrocopy Vrije Universiteit Brussel 9th December 2011 Outline 1 Relaxation 2 Principles 3 Outline 1 Relaxation 2 Principles 3 Establishment of Thermal Equilibrium As previously
More informationComparison between Bacteriorhodopsin and Halorhodopsin. Halorhodopsin (HR) and Bacteriorhodopsin (BR) belong to a subfamily of
Comparison between Bacteriorhodopsin and Halorhodopsin Halorhodopsin (HR) and Bacteriorhodopsin (BR) belong to a subfamily of heptahelical membrane proteins, the archaeal rhodopsins. They are found in
More informationHigh Frequency Electron Paramagnetic Resonance Studies of Mn 12 Wheels
High Frequency Electron Paramagnetic Resonance Studies of Mn 12 Wheels Gage Redler and Stephen Hill Department of Physics, University of Florida Abstract High Frequency Electron Paramagnetic Resonance
More information4. Stereochemistry of Alkanes and Cycloalkanes
4. Stereochemistry of Alkanes and Cycloalkanes Based on McMurry s Organic Chemistry, 6 th edition, Chapter 4 2003 Ronald Kluger Department of Chemistry University of Toronto The Shapes of Molecules! The
More informationBMB/Bi/Ch 173 Winter 2018
BMB/Bi/Ch 173 Winter 2018 Homework Set 8.1 (100 Points) Assigned 2-27-18, due 3-6-18 by 10:30 a.m. TA: Rachael Kuintzle. Office hours: SFL 220, Friday 3/2 4:00-5:00pm and SFL 229, Monday 3/5 4:00-5:30pm.
More informationScalar (contact) vs dipolar (pseudocontact) contributions to isotropic shifts.
Scalar (contact) vs dipolar (pseudocontact) contributions to isotropic shifts. Types of paramagnetic species: organic radicals, and complexes of transition metals, lanthanides, and actinides. Simplest
More informationProteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability
Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Dr. Andrew Lee UNC School of Pharmacy (Div. Chemical Biology and Medicinal Chemistry) UNC Med
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationNMR in Structural Biology
NMR in Structural Biology Exercise session 2 1. a. List 3 NMR observables that report on structure. b. Also indicate whether the information they give is short/medium or long-range, or perhaps all three?
More informationRHODOPSIN: Insights from Recent Structural Studies
Annu. Rev. Biophys. Biomol. Struct. 2002. 31:443 84 DOI: 10.1146/annurev.biophys.31.082901.134348 Copyright c 2002 by Annual Reviews. All rights reserved RHODOPSIN: Insights from Recent Structural Studies
More informationIntroduction to Computational Structural Biology
Introduction to Computational Structural Biology Part I 1. Introduction The disciplinary character of Computational Structural Biology The mathematical background required and the topics covered Bibliography
More informationVisual pigments. Neuroscience, Biochemistry Dr. Mamoun Ahram Third year, 2015
Visual pigments Neuroscience, Biochemistry Dr. Mamoun Ahram Third year, 2015 References Photoreceptors and visual pigments Webvision: The Organization of the Retina and Visual System (http://www.ncbi.nlm.nih.gov/books/nbk11522/#a127)
More informationAnalysis of NMR Spectra Part 2
Analysis of NMR Spectra Part 2-1- Analysis of NMR Spectra Part 2 "Things should be made as simple as possible, but not any simpler." Albert Einstein 1.1 Review of Basic NMR Concepts NMR analysis is a complex
More informationSpin Dynamics Basics of Nuclear Magnetic Resonance. Malcolm H. Levitt
Spin Dynamics Basics of Nuclear Magnetic Resonance Second edition Malcolm H. Levitt The University of Southampton, UK John Wiley &. Sons, Ltd Preface xxi Preface to the First Edition xxiii Introduction
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationEXAM I COURSE TFY4310 MOLECULAR BIOPHYSICS December Suggested resolution
page 1 of 7 EXAM I COURSE TFY4310 MOLECULAR BIOPHYSICS December 2013 Suggested resolution Exercise 1. [total: 25 p] a) [t: 5 p] Describe the bonding [1.5 p] and the molecular orbitals [1.5 p] of the ethylene
More informationSolid-state NMR and proteins : basic concepts (a pictorial introduction) Barth van Rossum,
Solid-state NMR and proteins : basic concepts (a pictorial introduction) Barth van Rossum, 16.02.2009 Solid-state and solution NMR spectroscopy have many things in common Several concepts have been/will
More information12. Spectral diffusion
1. Spectral diffusion 1.1. Spectral diffusion, Two-Level Systems Until now, we have supposed that the optical transition frequency of each single molecule is a constant (except when we considered its variation
More informationModule 20: Applications of PMR in Structural Elucidation of Simple and Complex Compounds and 2-D NMR spectroscopy
Subject Chemistry Paper No and Title Module No and Title Module Tag Paper 12: Organic Spectroscopy Module 20: Applications of PMR in Structural Elucidation of Simple and Complex Compounds and 2-D NMR spectroscopy
More informationConjugated Systems, Orbital Symmetry and UV Spectroscopy
Conjugated Systems, Orbital Symmetry and UV Spectroscopy Introduction There are several possible arrangements for a molecule which contains two double bonds (diene): Isolated: (two or more single bonds
More informationPrinciples of Physical Biochemistry
Principles of Physical Biochemistry Kensal E. van Hold e W. Curtis Johnso n P. Shing Ho Preface x i PART 1 MACROMOLECULAR STRUCTURE AND DYNAMICS 1 1 Biological Macromolecules 2 1.1 General Principles
More informationMOLECULAR DYNAMICS SIMULATIONS ON STRUCTURAL CONFORMATIONS OF RHODOPSIN AND PRION PROTEINS
Ó³ Ÿ. 2005.. 2, º 6(129).. 87Ä97 Š [51-72:530.145]+[51-72:541.1] MOLECULAR DYNAMICS SIMULATIONS ON STRUCTURAL CONFORMATIONS OF RHODOPSIN AND PRION PROTEINS Kh. T. Kholmurodov 1 Joint Institute for Nuclear
More informationDynamic Structure of Retinylidene Ligand of Rhodopsin Probed by Molecular Simulations
doi:10.1016/j.jmb.2007.06.047 J. Mol. Biol. (2007) 372, 906 917 Dynamic Structure of Retinylidene Ligand of Rhodopsin Probed by Molecular Simulations Pick-Wei Lau 1, Alan Grossfield 2, Scott E. Feller
More informatione 2m e c I, (7.1) = g e β B I(I +1), (7.2) = erg/gauss. (7.3)
Chemistry 126 Molecular Spectra & Molecular Structure Week # 7 Electron Spin Resonance Spectroscopy, Supplement Like the hydrogen nucleus, an unpaired electron in a sample has a spin of I=1/2. The magnetic
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 09 Quantum Mechanics/Molecular Mechanics (QM/MM) Techniques Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville,
More informationSUPPLEMENTARY INFORMATION
Figure S1. Secondary structure of CAP (in the camp 2 -bound state) 10. α-helices are shown as cylinders and β- strands as arrows. Labeling of secondary structure is indicated. CDB, DBD and the hinge are
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationForce Fields for MD simulations
Force Fields for MD simulations Topology/parameter files Where do the numbers an MD code uses come from? ow to make topology files for ligands, cofactors, special amino acids, ow to obtain/develop missing
More informationAlkanes. Introduction
Introduction Alkanes Recall that alkanes are aliphatic hydrocarbons having C C and C H bonds. They can be categorized as acyclic or cyclic. Acyclic alkanes have the molecular formula C n H 2n+2 (where
More informationFerdowsi University of Mashhad
Spectroscopy in Inorganic Chemistry Nuclear Magnetic Resonance Spectroscopy spin deuterium 2 helium 3 The neutron has 2 quarks with a -e/3 charge and one quark with a +2e/3 charge resulting in a total
More informationIntroduction to" Protein Structure
Introduction to" Protein Structure Function, evolution & experimental methods Thomas Blicher, Center for Biological Sequence Analysis Learning Objectives Outline the basic levels of protein structure.
More informationAll measurement has a limit of precision and accuracy, and this must be taken into account when evaluating experimental results.
Chapter 11: Measurement and data processing and analysis 11.1 Uncertainty and error in measurement and results All measurement has a limit of precision and accuracy, and this must be taken into account
More informationA Combined Optical and EPR Spectroscopy Study: Azobenzene-Based Biradicals as Reversible Molecular Photoswitches
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2017 A Combined Optical and EPR Spectroscopy Study: Azobenzene-Based Biradicals as Reversible
More information17.1 Classes of Dienes
17.1 Classes of Dienes There are three categories for dienes: Cumulated: pi bonds are adjacent. Conjugated: pi bonds are separated by exactly ONE single bond. Isolated: pi bonds are separated by any distance
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
More informationT H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y. jgp
S u p p l e m e n ta l m at e r i a l jgp Lee et al., http://www.jgp.org/cgi/content/full/jgp.201411219/dc1 T H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y S u p p l e m e n ta l D I S C U S
More informationElectron spin resonance
Quick reference guide Introduction This is a model experiment for electron spin resonance, for clear demonstration of interaction between the magnetic moment of the electron spin with a superimposed direct
More informationNMR Spectroscopy. Guangjin Hou
NMR Spectroscopy Guangjin Hou 22-04-2009 NMR History 1 H NMR spectra of water H NMR spectra of water (First NMR Spectra on Water, 1946) 1 H NMR spectra ethanol (First bservation of the Chemical Shift,
More informationsensors ISSN
Sensors 2009, 9, 9741-9804; doi:10.3390/s91209741 OPEN ACCESS sensors ISSN 1424-8220 www.mdpi.com/journal/sensors Review Photoreactions and Structural Changes of Anabaena Sensory Rhodopsin Akira Kawanabe
More informationNMR BMB 173 Lecture 16, February
NMR The Structural Biology Continuum Today s lecture: NMR Lots of slides adapted from Levitt, Spin Dynamics; Creighton, Proteins; And Andy Rawlinson There are three types of particles in the universe Quarks
More informationT 1, T 2, NOE (reminder)
T 1, T 2, NOE (reminder) T 1 is the time constant for longitudinal relaxation - the process of re-establishing the Boltzmann distribution of the energy level populations of the system following perturbation
More informationSlow symmetric exchange
Slow symmetric exchange ϕ A k k B t A B There are three things you should notice compared with the Figure on the previous slide: 1) The lines are broader, 2) the intensities are reduced and 3) the peaks
More information7a. Structure Elucidation: IR and 13 C-NMR Spectroscopies (text , , 12.10)
2009, Department of Chemistry, The University of Western Ontario 7a.1 7a. Structure Elucidation: IR and 13 C-NMR Spectroscopies (text 11.1 11.5, 12.1 12.5, 12.10) A. Electromagnetic Radiation Energy is
More informationBenzene a remarkable compound. Chapter 14 Aromatic Compounds. Some proposed structures for C 6 H 6. Dimethyl substituted benzenes are called xylenes
Benzene a remarkable compound Chapter 14 Aromatic Compounds Discovered by Faraday 1825 Formula C 6 H 6 Highly unsaturated, but remarkably stable Whole new class of benzene derivatives called aromatic compounds
More informationSTEREOCHEMISTRY OF ALKANES AND CYCLOALKANES CONFORMATIONAL ISOMERS
STEREOCHEMISTRY OF ALKANES AND CYCLOALKANES CONFORMATIONAL ISOMERS 1 CONFORMATIONAL ISOMERS Stereochemistry concerned with the 3-D aspects of molecules Rotation is possible around C-C bonds in openchain
More informationSupplementary Figure 1
Supplementary Figure 1 a Tyr191 Trp265 4-5 Å Tyr268 ~14 Å 4.3 Å 16 17 19 20 4.4 Å 11 15 Lys296 N 8 10 12 14 7.3 Å 18 Tyr191 4-5 Å 18 3.8-4 Å 8 10 12 11 17 16 19 ~ Å ~4 Å Tyr268 5.5-6.3 Å 14 N Lys296 15
More informationMeasuring Spin-Lattice Relaxation Time
WJP, PHY381 (2009) Wabash Journal of Physics v4.0, p.1 Measuring Spin-Lattice Relaxation Time L.W. Lupinski, R. Paudel, and M.J. Madsen Department of Physics, Wabash College, Crawfordsville, IN 47933 (Dated:
More informationProtein Dynamics. The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron.
Protein Dynamics The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron. Below is myoglobin hydrated with 350 water molecules. Only a small
More informationCOPYRIGHTED MATERIAL. Production of Net Magnetization. Chapter 1
Chapter 1 Production of Net Magnetization Magnetic resonance (MR) is a measurement technique used to examine atoms and molecules. It is based on the interaction between an applied magnetic field and a
More information17.1 Classes of Dienes
W 2/1 Due: HW14, spec03 Due: n/a M 2/6 Lecture HW14 grading Lect17a Conjugated π systems Lecture quiz Lect17b Lab Lab02 Qual Analysis II (cont) 7-1 17.1 Classes of Dienes There are three categories for
More informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/4/1/eaau413/dc1 Supplementary Materials for Structure and dynamics conspire in the evolution of affinity between intrinsically disordered proteins Per Jemth*, Elin
More informationConformational Isomers. Isomers that differ as a result of sigma bond rotation of C-C bond in alkanes
Conformational Isomers Isomers that differ as a result of sigma bond Isomers that differ as a result of sigma bond rotation of C-C bond in alkanes Bond Rotation and Newman Projections As carbon-carbon
More informationBiochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationConjugated Dienes and Ultraviolet Spectroscopy
Conjugated Dienes and Ultraviolet Spectroscopy Key Words Conjugated Diene Resonance Structures Dienophiles Concerted Reaction Pericyclic Reaction Cycloaddition Reaction Bridged Bicyclic Compound Cyclic
More informationMagnetic Resonance Spectroscopy
INTRODUCTION TO Magnetic Resonance Spectroscopy ESR, NMR, NQR D. N. SATHYANARAYANA Formerly, Chairman Department of Inorganic and Physical Chemistry Indian Institute of Science, Bangalore % I.K. International
More informationIntroduction to Relaxation Theory James Keeler
EUROMAR Zürich, 24 Introduction to Relaxation Theory James Keeler University of Cambridge Department of Chemistry What is relaxation? Why might it be interesting? relaxation is the process which drives
More informationH B. θ = 90 o. Lecture notes Part 4: Spin-Spin Coupling. θ θ
Lecture notes Part 4: Spin-Spin Coupling F. olger Försterling October 4, 2011 So far, spins were regarded spins isolated from each other. owever, the magnetic moment of nuclear spins also have effect on
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature17991 Supplementary Discussion Structural comparison with E. coli EmrE The DMT superfamily includes a wide variety of transporters with 4-10 TM segments 1. Since the subfamilies of the
More informationLectures Spectroscopy. Fall 2012
Lectures 19-20 Spectroscopy Fall 2012 1 spectroscopic principles (Chem 1M/1N exps. #6 and #11) 4 1 spectroscopic excitations ( E = h = hc/ ; = c ) (nm) (sec -1 ) radiation technique molecular excitation
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
More informationCHAPTER 2. Structure and Reactivity: Acids and Bases, Polar and Nonpolar Molecules
CHAPTER 2 Structure and Reactivity: Acids and Bases, Polar and Nonpolar Molecules 2-1 Kinetics and Thermodynamics of Simple Chemical Processes Chemical thermodynamics: Is concerned with the extent that
More informationLectures Spectroscopy. Fall 2012
Lectures 19-20 Spectroscopy Fall 2012 1 spectroscopic principles (Chem 1M/1N exps. #6 and #11) 4 spectroscopic excitations ( E = h = hc/ ; = c ) (nm) (sec -1 ) radiation technique molecular excitation
More informationIntroduction to Comparative Protein Modeling. Chapter 4 Part I
Introduction to Comparative Protein Modeling Chapter 4 Part I 1 Information on Proteins Each modeling study depends on the quality of the known experimental data. Basis of the model Search in the literature
More informationTheory and Applications of Residual Dipolar Couplings in Biomolecular NMR
Theory and Applications of Residual Dipolar Couplings in Biomolecular NMR Residual Dipolar Couplings (RDC s) Relatively new technique ~ 1996 Nico Tjandra, Ad Bax- NIH, Jim Prestegard, UGA Combination of
More informationELECTRON PARAMAGNETIC RESONANCE
ELECTRON PARAMAGNETIC RESONANCE = MAGNETIC RESONANCE TECHNIQUE FOR STUDYING PARAMAGNETIC SYSTEMS i.e. SYSTEMS WITH AT LEAST ONE UNPAIRED ELECTRON Examples of paramagnetic systems Transition-metal complexes
More informationIntroduction to Alkenes and Alkynes
Introduction to Alkenes and Alkynes In an alkane, all covalent bonds between carbon were σ (σ bonds are defined as bonds where the electron density is symmetric about the internuclear axis) In an alkene,
More informationChapter 7. Nuclear Magnetic Resonance Spectroscopy
Chapter 7 Nuclear Magnetic Resonance Spectroscopy I. Introduction 1924, W. Pauli proposed that certain atomic nuclei have spin and magnetic moment and exposure to magnetic field would lead to energy level
More informationExcited States Calculations for Protonated PAHs
52 Chapter 3 Excited States Calculations for Protonated PAHs 3.1 Introduction Protonated PAHs are closed shell ions. Their electronic structure should therefore be similar to that of neutral PAHs, but
More information4. Circular Dichroism - Spectroscopy
4. Circular Dichroism - Spectroscopy The optical rotatory dispersion (ORD) and the circular dichroism (CD) are special variations of absorption spectroscopy in the UV and VIS region of the spectrum. The
More informationSupporting Information
Supporting Information Boehr et al. 10.1073/pnas.0914163107 SI Text Materials and Methods. R 2 relaxation dispersion experiments. 15 NR 2 relaxation dispersion data measured at 1 H Larmor frequencies of
More informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 2. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Functional Protein Dynamics Proteins are dynamic and
More informationInterpreting and evaluating biological NMR in the literature. Worksheet 1
Interpreting and evaluating biological NMR in the literature Worksheet 1 1D NMR spectra Application of RF pulses of specified lengths and frequencies can make certain nuclei detectable We can selectively
More information4. Organic photosynthetic reactions
4. rganic photosynthetic reactions 100 4.1 eactions of ethenes and aromatic compounds Photoreactivity of ethenes E Geometrical isomerization In π-π* excited states, there is effectively no π bond and so
More informationProton-Pump Mechanism in Retinal Schiff Base: On the molecular structure of the M-state
Proton-Pump Mechanism in Retinal Schiff Base: On the molecular structure of the M-state Ayan Datta and Swapan K. Pati * Theoretical Sciences Unit and Chemistry and Physics of Materials Unit, Jawaharlal
More informationhand and delocalization on the other, can be instructively exemplified and extended
Text Related to Segment 8.0 00 Claude E. Wintner The ideas developed up to this point, concerning stereochemistry on the one hand and delocalization on the other, can be instructively exemplified and extended
More informationTDDFT as a tool in biophysics
TDDFT as a tool in biophysics The primary event in vision Robert Send Universität Karlsruhe 09.09.08 Robert Send TDDFT as a tool in biophysics 09.09.08 1 / 28 Outline 1 Human vision 2 The methods 3 The
More informationLipid Regulated Intramolecular Conformational Dynamics of SNARE-Protein Ykt6
Supplementary Information for: Lipid Regulated Intramolecular Conformational Dynamics of SNARE-Protein Ykt6 Yawei Dai 1, 2, Markus Seeger 3, Jingwei Weng 4, Song Song 1, 2, Wenning Wang 4, Yan-Wen 1, 2,
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationAnalysis of Conformational Influences Within (2R,3R)-butanediol
Analysis of Conformational Influences Within (2R,3R)-butanediol Zoila M. Estrada Abstract (2,3)-Butanediol is a simple polyol with two different isomers: the meso isomer and the racemic isomer. We examined
More informationPhysical Background Of Nuclear Magnetic Resonance Spectroscopy
Physical Background Of Nuclear Magnetic Resonance Spectroscopy Michael McClellan Spring 2009 Department of Physics and Physical Oceanography University of North Carolina Wilmington What is Spectroscopy?
More informationModule 4 : Third order nonlinear optical processes. Lecture 28 : Inelastic Scattering Processes. Objectives
Module 4 : Third order nonlinear optical processes Lecture 28 : Inelastic Scattering Processes Objectives In this lecture you will learn the following Light scattering- elastic and inelastic-processes,
More informationChapter 2: Chemical Basis of Life
Chapter 2: Chemical Basis of Life Chemistry is the scientific study of the composition of matter and how composition changes. In order to understand human physiological processes, it is important to understand
More informationIntroduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th
Abstract When conductive, gramicidin monomers are linked by six hydrogen bonds. To understand the details of dissociation and how the channel transits from a state with 6H bonds to ones with 4H bonds or
More informationExperimental Techniques in Protein Structure Determination
Experimental Techniques in Protein Structure Determination Homayoun Valafar Department of Computer Science and Engineering, USC Two Main Experimental Methods X-Ray crystallography Nuclear Magnetic Resonance
More informationLS1a Fall 2014 Problem Set #2 Due Monday 10/6 at 6 pm in the drop boxes on the Science Center 2 nd Floor
LS1a Fall 2014 Problem Set #2 Due Monday 10/6 at 6 pm in the drop boxes on the Science Center 2 nd Floor Note: Adequate space is given for each answer. Questions that require a brief explanation should
More information