Channelrhodopsin unchained: Structure and mechanism of a light-gated cation channel

The new and vibrant field of optogenetics was founded by the seminal discovery of channelrhodopsin, the first light-gated cation channel. Despite the numerous applications that have revolutionised neurophysiology, the functional mechanism is far from understood on the molecular level. An arsenal of biophysical techniques has been established in the last decades of research on microbial rhodopsins. However, application of these techniques is hampered by the duration and the complexity of the photoreaction of channelrhodopsin compared with other microbial rhodopsins. A particular interest in resolving the molecular mechanism lies in the structural changes that lead to channel opening and closure. Here, we review the current structural and mechanistic knowledge that has been accomplished by integrating the static structure provided by X-ray crystallography and electron microscopy with time-resolved spectroscopic and electrophysiological techniques. The dynamical reactions of the chromophore are effectively coupled to structural changes of the protein, as shown by ultrafast spectroscopy. The hierarchical sequence of structural changes in the protein backbone that spans the time range from 10^− 12 s to 10^− 3 s prepares the channel to open and, consequently, cations can pass. Proton transfer reactions that are associated with channel gating have been resolved. In particular, glutamate 253 and aspartic acid 156 were identified as proton acceptor and donor to the retinal Schiff base. The reprotonation of the latter is the critical determinant for channel closure. The proton pathway that eventually leads to proton pumping is also discussed. This article is part of a Special Issue entitled: Retinal Proteins — You can teach an old dog new tricks.     

Interaction between S100P and the anti-allergy drug cromolyn

The S100P protein has been known to mediate cell proliferation by binding the receptor for advanced glycation end products (RAGE) to activate signaling pathways, such as the extracellular regulated kinase (ERK) and nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) pathways. S100P/RAGE signaling is involved in a variety of diseases, such as cancer, metastasis, and diabetes. Cromolyn is an anti-allergy drug that binds S100P to block the interaction between S100P and RAGE. In the present study, we characterized the properties of the binding between cromolyn and calcium-bound S100P using various biophysical techniques. The binding affinity for S100P and cromolyn was measured to be in the millimolar range by fluorescence spectroscopy. NMR-HSQC titration experiments and HADDOCK modeling was employed to determine the spatial structure of the proposed heterotetramer model of the S100P–cromolyn complex. Additional MD simulation results revealed the important properties in the complex stability and conformational flexibility of the S100P–cromolyn complex. This proposed model has provided an understanding of the molecular level interactions of S100P–cromolyn complex.     

MD simulations of the central pore of ryanodine receptors and sequence comparison with 2B protein from coxsackie virus

The regulation of intracellular Ca^2 + triggers a multitude of vital processes in biological cells. Ca^2 + permeable ryanodine receptors (RyRs) are the biggest known ion channels and play a key role in the regulation of intracellular calcium concentrations, particularly in muscle cells. In this study, we construct a computational model of the pore region of the skeletal RyR and perform molecular dynamics (MD) simulations. The dynamics and distribution of Ca^2 + around the luminal pore entry of the RyR suggest that Ca^2 + ions are channeled to the pore entry due to the arrangement of (acidic) amino acids at the extramembrane surface of the protein. This efficient mechanism of Ca^2 + supply is thought to be part of the mechanism of Ca^2 + conductance of RyRs. Viral myocarditis is predominantly caused by coxsackie viruses that induce the expression of the protein 2B which is known to affect intracellular Ca^2 + homeostasis in infected cells. From our sequence comparison, it is hypothesized, that modulation of RyR could be due to replacement of its transmembrane domains (TMDs) by those domains of the viral channel forming protein 2B of coxsackie virus. This article is part of a Special Issue entitled: Viral Membrane Proteins — Channels for Cellular Networking.     

Structural dynamics of the cell wall precursor lipid II in the presence and absence of the lantibiotic nisin

Representing a physiological “Achilles' heel”, the cell wall precursor lipid II (L_II) is a prime target for various classes of antibiotics. Over the years L_II-binding agents have been recognized as promising candidates and templates in the search for new antibacterial compounds to complement or replace existing drugs. To elucidate the molecular structural basis underlying L_II functional mechanism and to better understand if and how lantibiotic binding alters the molecular behavior of L_II, we performed molecular dynamics (MD) simulations of phospholipid membrane-embedded L_II in the absence and presence of the L_II-binding lantibiotic nisin. In a series of 2 × 4 independent, unbiased 100 ns MD simulations we sampled the conformational dynamics of nine L_II as well as nine L_II–nisin complexes embedded in an aqueous 150 mM NaCl/POPC phospholipid membrane environment. We found that nisin binding to L_II induces a reduction of L_II mobility and flexibility, an outward shift of the L_II pentapeptide, an inward movement of the L_II disaccharide section, and an overall deeper insertion of the L_II tail group into the membrane. The latter effect might indicate an initial step in adopting a stabilizing, scaffold-like structure in the process of nisin-induced membrane leakage. At the same time nisin conformation and L_II interaction remain similar to the 1WCO L_II–nisin NMR solution structure.     

Solution structure of the transmembrane domain of the insulin receptor in detergent micelles

The insulin receptor (IR) binds insulin and plays important roles in glucose homeostasis by regulating the tyrosine kinase activity at its C-terminus. Its transmembrane domain (TMD) is shown to be important for transferring conformational changes induced by insulin across the cell membrane to regulate kinase activity. In this study, a construct IR_940–988 containing the TMD was expressed and purified for structural studies. Its solution structure in dodecylphosphocholine (DPC) micelles was determined. The sequence containing residues L962 to Y976 of the TMD of the IR in micelles adopts a well-defined helical structure with a kink formed by glycine and proline residues present at its N-terminus, which might be important for its function. Paramagnetic relaxation enhancement (PRE) and relaxation experimental results suggest that residues following the TMD are flexible and expose to aqueous solution. Although purified IR_940–988 in micelles existed mainly as a monomeric form verified by gel filtration and relaxation analysis, cross-linking study suggests that it may form a dimer or oligomers under micelle conditions.     

Characterization of the water defect at the HIV-1 gp41 membrane spanning domain in bilayers with and without cholesterol using molecular simulations

The membrane spanning domain (MSD) of human immunodeficiency virus 1 (HIV-1) envelope glycoprotein gp41 is important for fusion and infection. We used molecular dynamics (MD) simulations (3.4 μs total) to relate membrane and peptide properties that lead to water solvation of the α-helical gp41 MSD's midspan arginine in pure dipalmitoylphosphatidylcholine (DPPC) and in 50/50 DPPC/cholesterol membranes. We find that the midspan arginine is solvated by water that penetrates the inner leaflet, leading to a so-called water defect. The water defect is surprisingly robust across initial conditions and membrane compositions, but the presence of cholesterol modulates its behavior in several key ways. In the cholesterol-containing membranes, fluctuations in membrane thickness and water penetration depth are localized near the midspan arginine, and the MSD helices display a tightly regulated tilt angle. In the cholesterol-free membranes, thickness fluctuations are not as strongly correlated to the peptide position and tilt angles vary significantly depending on protein position relative to boundaries between domains of differing thickness. Cholesterol in an HIV-1 viral membrane is required for infection. Therefore, this work suggests that the colocalized water defect and membrane thickness fluctuations in cholesterol-containing viral membranes play an important role in fusion by bringing the membrane closer to a stability limit that must be crossed for fusion to occur.     

Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

Revisiting the bilayer structures of fluid phase phosphatidylglycerol lipids: Accounting for exchangeable hydrogens

We recently published two papers detailing the structures of fluid phase phosphatidylglycerol (PG) lipid bilayers (Ku?erka et al., 2012 J. Phys. Chem. B 116: 232–239; Pan et al., 2012 Biochim. Biophys. Acta Biomembr. 1818: 2135–2148), which were determined using the scattering density profile model. This hybrid experimental/computational technique utilizes molecular dynamics simulations to parse a lipid bilayer into components whose volume probabilities follow simple analytical functional forms. Given the appropriate scattering densities, these volume probabilities are then translated into neutron scattering length density (NSLD) and electron density (ED) profiles, which are used to jointly refine experimentally obtained small angle neutron and X-ray scattering data. However, accurate NSLD and ED profiles can only be obtained if the bilayer's chemical composition is known. Specifically, in the case of neutron scattering, the lipid's exchangeable hydrogens with aqueous D₂O must be accounted for, as they can have a measureable effect on the resultant lipid bilayer structures. This was not done in our above-mentioned papers. Here we report on the molecular structures of PG lipid bilayers by appropriately taking into account the exchangeable hydrogens. Analysis indicates that the temperature-averaged PG lipid areas decrease by 1.5 to 3.8 Å², depending on the lipid's acyl chain length and unsaturation, compared to PG areas when hydrogen exchange was not taken into account.     

Low-mass molecular dynamics simulation: A simple and generic technique to enhance configurational sampling

CLN025 is one of the smallest fast-folding proteins. Until now it has not been reported that CLN025 can autonomously fold to its native conformation in a classical, all-atom, and isothermal–isobaric molecular dynamics (MD) simulation. This article reports the autonomous and repeated folding of CLN025 from a fully extended backbone conformation to its native conformation in explicit solvent in multiple 500-ns MD simulations at 277 K and 1 atm with the first folding event occurring as early as 66.1 ns. These simulations were accomplished by using AMBER forcefield derivatives with atomic masses reduced by 10-fold on Apple Mac Pros. By contrast, no folding event was observed when the simulations were repeated using the original AMBER forcefields of FF12SB and FF14SB. The results demonstrate that low-mass MD simulation is a simple and generic technique to enhance configurational sampling. This technique may propel autonomous folding of a wide range of miniature proteins in classical, all-atom, and isothermal–isobaric MD simulations performed on commodity computers—an important step forward in quantitative biology.