Blake Mertz

Retinal Ligand Mobility Explains Internal Hydration and Reconciles Active Rhodopsin Structures

Rhodopsin, the mammalian dim-light receptor, is one of the best-characterized G-protein-coupled receptors, a pharmaceutically important class of membrane proteins that has garnered a great deal of attention because of the recent availability of structural information. Yet the mechanism of rhodopsin activation is not fully understood. Here, we use microsecond-scale all-atom molecular dynamics simulations, validated by solid-state (2)H nuclear magnetic resonance spectroscopy, to understand the transition between the dark and metarhodopsin I (Meta I) states. Our analysis of these simulations reveals striking differences in ligand flexibility between the two states. Retinal is much more dynamic in Meta I, adopting an elongated conformation similar to that seen in the recent activelike crystal structures. Surprisingly, this elongation corresponds to both a dramatic influx of bulk water into the hydrophobic core of the protein and a concerted transition in the highly conserved Trp265(6.48) residue. In addition, enhanced ligand flexibility upon light activation provides an explanation for the different retinal orientations observed in X-ray crystal structures of active rhodopsin.

Computational analyses of the conformational itinerary along the reaction pathway of GH94 cellobiose phosphorylase.

GH94 cellobiose phosphorylase (CBP) catalyzes the phosphorolysis of cellobiose into alpha-D-glucose 1-phosphate (G1P) and D-glucose with inversion of anomeric configuration. The complex crystal structure of CBP from Cellvibrio gilvus had previously been determined; glycerol, glucose, and phosphate are bound to subsites -1, +1, and the anion binding site, respectively. We performed computational analyses to elucidate the conformational itinerary along the reaction pathway of this enzyme. autodock was used to dock cellobiose with its glycon glucosyl residue in various conformations and with its aglycon glucosyl residue in the low-energy 4C1 conformer. An oxocarbenium ion-like glucose molecule mimicking the transition state was also docked. Based on the clustering analysis, docked energies, and comparison with the crystallographic ligands, we conclude that the reaction proceeds from 1S3 as the pre-transition state conformer (Michaelis complex) via E3 as the transition state candidate to 4C1 as the G1P product conformer. The predicted reaction pathway of the inverting phosphorylase is similar to that proposed for the first-half glycosylation reaction of retaining cellulases, but is different from those for inverting cellulases. NAMD was used to simulate molecular dynamics of the enzyme. The 1S3 pre-transition state conformer is highly stable compared with other conformers, and a conformational change from 4C1 to 1,4B was observed.

Molecular simulations and solid-state NMR investigate dynamical structure in rhodopsin activation ☆

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.

Steric and electronic influences on the torsional energy landscape of retinal

We have performed quantum mechanical calculations for retinal model compounds to establish the rotational energy barriers for the C5-, C9-, and C13-methyl groups known to play an essential role in rhodopsin activation. Intraretinal steric interactions as well as electronic effects lower the rotational barriers of both the C9- and C13-methyl groups, consistent with experimental (2)H NMR data. Each retinal methyl group has a unique rotational behavior which must be treated individually. These results are highly relevant for the parameterization of molecular mechanics force fields which form the basis of molecular dynamics simulations of retinal proteins such as rhodopsin.

Analysis of functional divergence within two structurally related glycoside hydrolase families

Two glycoside hydrolase (GH) families were analyzed to detect the presence of functional divergence using the program DIVERGE. These two families, GH7 and GH16, each contain members related by amino acid sequence similarity, retaining hydrolytic mechanisms, and catalytic residue identity. GH7 and GH16 comprise GH Clan B, with a shared β‐jelly roll topology and mechanism. GH7 contains fungal cellobiohydrolases and endoglucanases and is divided into five main subfamilies, four of the former and one of the latter. Cluster comparisons between three of the cellobiohydrolase subfamilies and the endoglucanase subfamily identified specific amino acid residues that play a role in the functional divergence between the two enzyme types. GH16 contains subfamilies of bacterial agarases, xyloglucosyl transferases, 1,3‐β‐D‐glucanases, lichenases, and other enzymes with various substrate specificities and product profiles. Four cluster comparisons between these four main subfamilies again have identified amino acid residues involved in functional divergence between the subfamilies.

Novel Phosphotidylinositol 4,5-Bisphosphate Binding Sites on Focal Adhesion Kinase.

Focal adhesion kinase (FAK) is a protein tyrosine kinase that is ubiquitously expressed, recruited to focal adhesions, and engages in a variety of cellular signaling pathways. Diverse cellular responses, such as cell migration, proliferation, and survival, are regulated by FAK. Prior to activation, FAK adopts an autoinhibited conformation in which the FERM domain binds the kinase domain, blocking access to the activation loop and substrate binding site. Activation of FAK occurs through conformational change, and acidic phospholipids such as phosphatidylinositol 4,5-bisphosphate (PIP2) are known to facilitate this process. PIP2 binding alters the autoinhibited conformation of the FERM and kinase domains and subsequently exposes the activation loop to phosphorylation. However, the detailed molecular mechanism of PIP2 binding and its role in FAK activation remain unclear. In this study, we conducted coarse-grained molecular dynamics simulations to investigate the binding of FAK to PIP2. Our simulations identified novel areas of basic residues in the kinase domain of FAK that potentially undergo transient binding to PIP2 through electrostatic attractions. Our investigation provides a molecular picture of PIP2-initiated FAK activation and introduces promising new pathways for future studies of FAK regulation.

Proteorhodopsin Activation Is Modulated by Dynamic Changes in Internal Hydration.

Proteorhodopsin, a member of the microbial rhodopsin family, is a seven-transmembrane α-helical protein that functions as a light-driven proton pump. Understanding the proton-pumping mechanism of proteorhodopsin requires intimate knowledge of the proton transfer pathway via complex hydrogen-bonding networks formed by amino acid residues and internal water molecules. Here we conducted a series of microsecond time scale molecular dynamics simulations on both the dark state and the initial photoactivated state of blue proteorhodopsin to reveal the structural basis for proton transfer with respect to protein internal hydration. A complex series of dynamic hydrogen-bonding networks involving water molecules exists, facilitated by water channels and hydration sites within proteorhodopsin. High levels of hydration were discovered at each proton transfer site-the retinal binding pocket and proton uptake and release sites-underscoring the critical participation of water molecules in the proton-pumping mechanism. Water-bridged interactions and local water channels were also observed and can potentially mediate long-distance proton transfer between each site. The most significant phenomenon is after isomerization of retinal, an increase in water flux occurs that connects the proton release group, a conserved arginine residue, and the retinal binding pocket. Our results provide a detailed description of the internal hydration of the early photointermediates in the proteorhodopsin photocycle under alkaline pH conditions. These results lay the fundamental groundwork for understanding the intimate role that hydration plays in the structure-function relationship underlying the proteorhodopsin proton-pumping mechanism, as well as providing context for the relationship of hydration in proteorhodopsin to other microbial retinal proteins.

Retinal Flip in Rhodopsin Activation

Rhodopsin is a well-characterized structural model of a G protein-coupled receptor. Photoisomerization of the covalently bound retinal triggers activation. Surprisingly, the x-ray crystal structure of the active Meta-II state has a 180° rotation about the long-axis of the retinal polyene chain. Unbiased microsecond-timescale all-atom molecular dynamics simulations show that the retinal cofactor can flip back to the orientation observed in the inactive state of rhodopsin under conditions favoring the Meta-I state. Our results provide, to our knowledge, the first evidence from molecular dynamics simulations showing how rotation of the retinal ligand within its binding pocket can occur in the activation mechanism of rhodopsin.

Molecular Dynamics Simulations as a Tool for Accurate Determination of Surfactant Micelle Properties

Molecular dynamics (MD) simulations were used to characterize the equilibrium size, shape, hydration, and self-assembly of dodecylphosphocholine (DPC) and dodecyl-β-D-maltoside (DDM) micelles. We show that DPC molecules self-assemble to form micelles with sizes within the range reported in the experimental literature. The equilibrium shape of DPC and DDM micelles as well as associated micellar radii are in agreement with small-angle X-ray scattering (SAXS) experiments and theoretical packing parameters. In addition, we show that hydration of the micelle interior is limited; however, flexibility of the acyl chains leads to dynamic encounters with the solvated outer shell of the micelle, providing an explanation for long-standing differences in models of micelle hydration. Altogether, our results provide fundamental understanding of physical characteristics of micelles that can be utilized to study other types of detergents and proteomicelle complexes.

Theoretical Evidence for Multiple Charge Transfer Pathways in Bacteriorhodopsin

The development of molecular-scale junctions utilizing biomolecules is a challenging field that requires intimate knowledge of the relationship between molecular structure and conductance characteristics. One of the key parameters to understanding conductance efficiency is the charge mobility, which strongly influences the response time of electronic devices. The charge mobility of bacteriorhodopsin (bR), a membrane protein that has been studied experimentally in detail, was theoretically investigated using extended Marcus-Hush theory. Charge mobilities of 1.3 × 10(-2) and 9.7 × 10(-4) cm(2)/(V s) for hole and electron transfer, respectively, were determined. The computed electron mobility is comparable to experimentally measured values (9 × 10(-4) cm(2)/(V s)). Interestingly, the pathways for hole and electron hopping were very distinct from each other, utilizing different transmembrane helices to traverse the protein. In particular, only the electron transfer pathway involved the retinal chromophore, indicating that the efficiency of charge transfer is directly affected by the tertiary arrangement of proteins. Our results provide a template for obtaining the molecular and electronic-level details that can reveal fundamental insights into experimental studies on protein electron transport and inform efficient design of biomolecular-based junctions on the nanoscale.