Membrane proteins that function as transporters or receptors must communicate with both sides of the lipid bilayer in which they sit. This long distance communication enables transporters to move protons or other ions and small molecules across the bilayer and receptors to transmit an external signal to the cell. Hydrogen bonds, hydrogen-bond networks, and lipid−protein interactions are essential for the motions and functioning of the membrane protein and, consequently, of outmost interest to structural biology and numerical simulations. We present here Bridge, an algorithm tailored for efficient analyses of hydrogen-bond networks in membrane transporter and receptor proteins. For channelrhodopsin, a membrane protein whose functioning involves proton-transfer reactions, Bridge identifies extensive networks of protein−water hydrogen bonds and an unanticipated network that can bridge transiently two proton donors across a distance of ∼20 Å. Graphs of the protein hydrogen bonds reveal rapid propagation of structural changes within hydrogen-bond networks of mutant transporters and identify protein groups potentially important for the proton transfer activity. The algorithm is made available as a plugin for PyMol.

Although the outward-directed proton transport across biological membranes is well studied and its importance for bioenergetics is clearly understood, inward-directed light-driven proton pumping by microbial rhodopsins has remained a mystery both physiologically and mechanistically. A new family of Antarctic rhodopsins, which is a subgroup within a novel class of schizorhodopsins reported recently, includes a member, denoted as AntR, which proved amenable to extensive characterization with experiments and computation. Phylogenetic analyses identify AntR as distinct from the well-studied microbial rhodopsins that function as outward-directed ion pumps, and bioinformatics sequence analyses reveal amino acid substitutions at conserved sites essential for outward proton pumping. Modeling and numerical simulations of AntR, combined with advanced analyses using the graph theory and centrality measures from social sciences, identify the dynamic three-dimensional network of hydrogen-bonded water molecules and amino acid residues that function as communication hubs in AntR. This network undergoes major rearrangement upon retinal isomerization, showing important changes in the connectivity of the active center, retinal Schiff base, to the opposing sides of the membrane, as required for proton transport. Numerical simulations and experimental studies of the photochemical cycle of AntR by spectroscopy and sitedirected mutagenesis allowed us to identify pathways that could conduct protons in the direction opposite to that commonly known for outward-directed pumps.

Graphical abstract

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Adenylyl cyclase 9 (AC9) is a membrane-bound enzyme that converts ATP into cAMP. The enzyme is weakly activated by forskolin, fully activated by the G protein Gαs subunit and is autoinhibited by the AC9 C-terminus. Although our recent structural studies of the AC9-Gαs complex provided the framework for understanding AC9 autoinhibition, the conformational changes that AC9 undergoes in response to activator binding remains poorly understood. Here, we present the cryo-EM structures of AC9 in several distinct states: (i) AC9 bound to a nucleotide inhibitor MANT-GTP, (ii) bound to an artificial activator (DARPin C4) and MANT-GTP, (iii) bound to DARPin C4 and a nucleotide analogue ATPαS, (iv) bound to Gαs and MANT-GTP. The artificial activator DARPin C4 partially activates AC9 by binding at a site that overlaps with the Gαs binding site. Together with the previously observed occluded and forskolin-bound conformations, structural comparisons of AC9 in the four conformations described here show that secondary structure rearrangements in the region surrounding the forskolin binding site are essential for AC9 activation.