Channelrhodopsin C1C2: Photocycle kinetics and interactions near the central gate.

Channelrhodopsins (ChR) are light-sensitive cation channels used in optogenetics, a technique that applies light to control cells (e.g., neurons) that have been modified genetically to express those channels. Although mutations are known to affect pore kinetics, little is known about how mutations induce changes at the molecular scale. To address this issue, we first measured channel opening and closing rates of a ChR chimera (C1C2) and selected variants (N297D, N297V, and V125L). Then, we used atomistic simulations to correlate those rates with changes in pore structure, hydration, and chemical interactions among key gating residues of C1C2 in both closed and open states. Overall, the experimental results show that C1C2 and its mutants do not behave like ChR2 or its analogous variants, except V125L, making C1C2 a unique channel. Our atomistic simulations confirmed that opening of the channel and initial hydration of the gating regions between helices I, II, III, and VII of the channel occurs with 1) the presence of 13-cis retinal; 2) deprotonation of a glutamic acid gating residue, E129; and 3) subsequent weakening of the central gate hydrogen bond between the same glutamic acid E129 and asparagine N297 in the central region of the pore. Also, an aspartate (D292) is the unambiguous primary proton acceptor for the retinal Schiff base in the hydrated channel.

Atomistic Study of Intramolecular Interactions in the Closed-State Channelrhodopsin Chimera, C1C2.

Channelrhodopsins (ChR1 and ChR2) are light-activated ion channels that enable photomobility of microalgae from the genus Chlamydomonas. Despite common use of ChR2 in optogenetics for selective control and monitoring of individual neurons in living tissue, the protein structures remain unresolved. Instead, a crystal structure of the ChR chimera (C1C2), an engineered combination of helices I-V from ChR1, without its C-terminus, and helices VI-VII from ChR2, is used as a template for ChR2 structure prediction. Surprisingly few studies have focused in detail on the chimera. Here, we present atomistic molecular dynamics studies of the closed-state, non-conducting C1C2 structure and protonation states. A new and comprehensive characterization of interactions in the vicinity of the gating region of the pore, namely between residues E90, E123, D253, N258, and the protonated Schiff base (SBH), as well as nearby residues K93, T127, and C128, indicates that the equilibrated C1C2 structure with both E123 and D253 deprotonated closely resembles the available crystal structure. In agreement with experimental studies on C1C2, no direct or water-mediated hydrogen bonding between an aspartate and a cysteine (D156-O…S-C128) that would define a direct-current gate in C1C2 was observed in our simulations. Finally, we show that a single hydrogen bond between a glutamic acid (E90) and an asparagine (N258) residue suffices to keep the gate of C1C2 closed and to disable free water and ion passage through the putative pore, in contrast to the double bond proposed earlier for ChR2. We anticipate that this work will provide context for studies of both the gating process and water and ion transport in C1C2, and will spark interest in further experimental studies on the chimera.

Direct observation of bis(dicarbollyl)nickel conformers in solution by fluorescence spectroscopy: an approach to redox-controlled metallacarborane molecular motors.

As a continuation of work on metallacarborane-based molecular motors, the structures of substituted bis(dicarbollyl)nickel complexes in Ni(III) and Ni(IV) oxidation states were investigated in solution by fluorescence spectroscopy. Symmetrically positioned cage-linked pyrene molecules served as fluorescent probes to enable the observation of mixed meso-trans/dl-gauche (pyrene monomer fluorescence) and dl-cis/dl-gauche (intramolecular pyrene excimer fluorescence with residual monomer fluorescence) cage conformations of the nickelacarboranes in the Ni(III) and Ni(IV) oxidation states, respectively. The absence of energetically disfavored conformers in solution--dl-cis in the case of nickel(III) complexes and meso-trans in the case of nickel(IV)--was demonstrated based on spectroscopic data and conformer energy calculations in solution. The conformational persistence observed in solution indicates that bis(dicarbollyl)nickel complexes may provide attractive templates for building electrically driven and/or photodriven molecular motors.

Molecular Dynamics Simulations of Channelrhodopsin Chimera, C1C2.

Molecular dynamics (MD) simulations have been successfully used for modeling dynamic behavior of biologically relevant systems, such as ion channels in representative environments to decode protein structure-function relationships. Protocol presented here describes steps for generating input files and modeling a monomer of transmembrane cation channel, channelrhodopsin chimera (C1C2), in representative environment of 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) planar lipid bilayer, TIP3P water and ions (Na+ and Cl-) using molecular dynamics package NAMD, molecular graphics/analysis tool VMD, and other relevant tools. MD simulations of C1C2 were performed at 303.15 K and in constant particle number, isothermal-isobaric (NpT) ensemble. The results of modeling have helped understand how key interactions in the center of the C1C2 channel contribute to channel gating and subsequent solvent transport across the membrane.

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2.

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

Carborane-derived local anesthetics are isomer dependent.

Clinically there is a need for local anesthetics with a greater specificity of action on target cells and longer duration. We have synthesized a series of local anesthetic derivatives we call boronicaines in which the aromatic phenyl ring of lidocaine was replaced with ortho-, meta-, C,C'-dimethyl meta- and para-carborane clusters. The boronicaine derivatives were tested for their analgesic activity and compared with lidocaine using standard procedures in mice following a plantar injection. The compounds differed in their analgesic activity in the following order: ortho-carborane = C,C'-dimethyl meta-carborane > para-carborane > lidocaine > meta-carborane derivative. Both ortho-boronicaine and C,C'-dimethyl meta-boronicaine had longer durations of analgesia than lidocaine. Differences in analgesic efficacies are rationalized by variations in chemical structure and protein binding characteristics.