Chromophore-Protein Interactions Affecting the Polyene Twist and π-π* Energy Gap of the Retinal Chromophore in Schizorhodopsins

The properties of a prosthetic group are broadened by interactions with its neighboring residues in proteins. The retinal chromophore in rhodopsins absorbs light, undergoes structural changes, and drives functionally important structural changes in proteins during the photocycle. It is therefore crucial to understand how chromophore-protein interactions regulate the molecular structure and electronic state of chromophores in rhodopsins. Schizorhodopsin is a newly discovered subfamily of rhodopsins found in the genomes of Asgard archaea, which are extant prokaryotes closest to the last common ancestor of eukaryotes and of other microbial species. Here, we report the effects of a hydrogen bond between a retinal Schiff base and its counterion on the twist of the polyene chain and the color of the retinal chromophore. Correlations between spectral features revealed the unexpected fact that the twist of the polyene chain is reduced as the hydrogen bond becomes stronger, suggesting that the twist is caused by tight atomic contacts between the chromophore and nearby residues. In addition, the strength of the hydrogen bond is the primary factor affecting the color-tuning of the retinal chromophore in schizorhodopsins. The findings of this study are valuable for manipulating the molecular structure and electronic state of the chromophore by controlling chromophore-protein interactions.

Cis-Trans Reisomerization Preceding Reprotonation of the Retinal Chromophore Is Common to the Schizorhodopsin Family: A Simple and Rational Mechanism for Inward Proton Pumping

The creation of unidirectional ion transporters across membranes represents one of the greatest challenges in chemistry. Proton-pumping rhodopsins are composed of seven transmembrane helices with a retinal chromophore bound to a lysine side chain via a Schiff base linkage and provide valuable insights for designing such transporters. What makes these transporters particularly intriguing is the discovery of both outward and inward proton-pumping rhodopsins. Surprisingly, despite sharing identical overall structures and membrane topologies, these proteins facilitate proton transport in opposite directions, implying an underlying rational mechanism that can transport protons in different directions within similar protein structures. In this study, we unraveled this mechanism by examining the chromophore structures of deprotonated intermediates in schizorhodopsins, a recently discovered subfamily of inward proton-pumping rhodopsins, using time-resolved resonance Raman spectroscopy. The photocycle of schizorhodopsins revealed the cis-trans thermal isomerization that precedes reprotonation at the Schiff base of the retinal chromophore. Notably, this order has not been observed in other proton-pumping rhodopsins, but here, it was observed in all seven schizorhodopsins studied across the archaeal domain, strongly suggesting that cis-trans thermal isomerization preceding reprotonation is a universal feature of the schizorhodopsin family. Based on these findings, we propose a structural basis for the remarkable order of events crucial for facilitating inward proton transport. The mechanism underlying inward proton transport by schizorhodopsins is straightforward and rational. The insights obtained from this study hold great promise for the design of transmembrane unidirectional ion transporters.

Heme Pocket Structure and Its Functional Implications in an Ancestral Globin Protein

Proteins have undergone evolutionary processes to achieve optimal stability, increased functionality, and novel functions. Comparative analysis of existent and ancestral proteins provides insights into the factors that influence protein stability and function. Ancestral sequence reconstruction allows us to deduce the amino acid sequences of ancestral proteins. Here, we present the structural and functional characteristics of an ancestral protein, AncMH, reconstructed to be the last common ancestor of hemoglobins and myoglobins. Our findings reveal that AncMH harbors heme and that the heme binds oxygen. Furthermore, we demonstrate that the ferrous heme in AncMH is pentacoordinated, similar to that of human adult hemoglobin and horse myoglobin. A detailed comparison of the heme pocket structure indicates that the heme pocket in AncMH is more similar to that of hemoglobin than that of myoglobin. However, the autoxidation of AncMH is faster than that of both hemoglobin and myoglobin. Collectively, our results suggest that ancestral proteins of hemoglobins and myoglobins evolved in steps, including the hexa- to pentacoordination transition, followed by stabilization of the oxygen-bound form.

Locating dynamic contributions to allostery via determining rates of vibrational energy transfer

Determining rates of energy transfer across non-covalent contacts for different states of a protein can provide information about dynamic and associated entropy changes during transitions between states. We investigate the relationship between rates of energy transfer across polar and nonpolar contacts and contact dynamics for the β₂-adrenergic receptor, a rhodopsin-like G-protein coupled receptor, in an antagonist-bound inactive state and agonist-bound active state. From structures sampled during molecular dynamics (MD) simulations, we find the active state to have, on average, a lower packing density, corresponding to generally more flexibility and greater entropy than the inactive state. Energy exchange networks (EENs) are computed for the inactive and active states from the results of the MD simulations. From the EENs, changes in the rates of energy transfer across polar and nonpolar contacts are found for contacts that remain largely intact during activation. Change in dynamics of the contact, and entropy associated with the dynamics, can be estimated from the change in rates of energy transfer across the contacts. Measurement of change in the rates of energy transfer before and after the transition between states thereby provides information about dynamic contributions to activation and allostery.

Time-resolved spectroscopic mapping of vibrational energy flow in proteins: Understanding thermal diffusion at the nanoscale

Vibrational energy exchange between various degrees of freedom is critical to barrier-crossing processes in proteins. Hemeproteins are well suited for studying vibrational energy exchange in proteins because the heme group is an efficient photothermal converter. The released energy by heme following photoexcitation shows migration in a protein moiety on a picosecond timescale, which is observed using time-resolved ultraviolet resonance Raman spectroscopy. The anti-Stokes ultraviolet resonance Raman intensity of a tryptophan residue is an excellent probe for the vibrational energy in proteins, allowing the mapping of energy flow with the spatial resolution of a single amino acid residue. This Perspective provides an overview of studies on vibrational energy flow in proteins, including future perspectives for both methodologies and applications.

Energy Transport in Class B GPCRs: Role of Protein-Water Dynamics and Activation

We compute energy exchange networks (EENs) through glucagon-like peptide-1 receptor (GLP-1R), a class B G-protein-coupled receptor (GPCR), in inactive and two active states, one activated by a peptide ligand and the other by a small molecule agonist, from results of molecular dynamics simulations. The reorganized network upon activation contains contributions from structural as well as from dynamic changes and corresponding entropic contributions to the free energy of activation, which are estimated in terms of the change in rates of energy transfer across non-covalent contacts. The role of water in the EENs and in activation of GLP-1R is also investigated. The dynamics of water in contact with the central polar network of the transmembrane region is found to be significantly slower for both activated states compared to the inactive state. This result is consistent with the contribution of water molecules to activation of GLP-1R previously suggested and resembles water dynamics in parts of the transmembrane region found in earlier studies of rhodopsin-like GPCRs.

Dependence of Vibrational Energy Transfer on Distance in a Four-Helix Bundle Protein: Equidistant Increments with the Periodicity of α Helices

Vibrational energy exchanges between various degrees of freedom are critical to barrier-crossing processes in proteins. Heme proteins are highly suitable for studies of the vibrational energy exchanges in proteins. The migration of excess energy released by heme in a protein moiety can be observed using time-resolved anti-Stokes ultraviolet resonance Raman spectroscopy. The anti-Stokes resonance Raman intensity of a tryptophan residue is an excellent probe for the excess energy and the spatial resolution of a single amino acid residue can be achieved. Here, we studied dependence of vibrational energy transfer on the distance in cytochrome b₅₆₂, which is a heme-containing, four-helix bundle protein. The vibrational energy transfer from the heme group to a single tryptophan residue introduced by site-directed mutagenesis was examined for different heme-tryptophan distances by a quasi-constant length with the periodicity of α helices. Taken together with structural data obtained by molecular dynamics simulations, the energy transfer could be well described by the model of classical thermal diffusion, which suggests that continuum media provide a good approximation of the protein interior, of which the atomic packing density is very high.

FTIR and Raman Spectroscopy of Rhodopsins

Vibrational spectroscopy such as FTIR and Raman spectroscopy is a powerful, sensitive, and informative method for studying protein structural changes in rhodopsins during their functions. The usefulness has been historically proven for the study of bacteriorhodopsin and bovine rhodopsin before their structural determination of rhodopsins. We now have atomic structures of many animal and microbial rhodopsins, and it is now important to know the structural dynamics of rhodopsins for function. FTIR and Raman spectroscopy provides useful information for this aim. In this chapter, we introduce the methods of FTIR and resonance Raman spectroscopy applied to rhodopsins. These vibrational methods offer deeper understanding on the mechanism how rhodopsins change their structures for function.

Activation-Induced Reorganization of Energy Transport Networks in the β2 Adrenergic Receptor

We compute energy exchange networks (EENs) through the β₂ adrenergic receptor (β₂AR), a G-protein coupled receptor (GPCR), in inactive and active states, based on the results of molecular dynamics simulations of this membrane bound protein. We introduce a new definition for the reorganization of EENs upon activation that depends on the relative change in rates of energy transfer across noncovalent contacts throughout the protein. On the basis of the reorganized network that we obtain for β₂AR upon activation, we identify a branched pathway between the agonist binding site and the cytoplasmic region, where a G-protein binds to the receptor when activated. The pathway includes all of the motifs containing molecular switches previously identified as contributing to the allosteric transition of β₂AR upon agonist binding. EENs and their reorganization upon activation are compared with structure-based contact networks computed for the inactive and active states of β₂AR.

Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations

Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.

Energy flow and intersubunit signalling in GSAM: A non-equilibrium molecular dynamics study

Non-equilibrium molecular dynamics simulations of vibrational energy flow induced by the imposition of a thermal gradient have been performed on the μ2-dimeric enzyme glutamate-1-semialdehyde aminomutase (GSAM), the key enzyme in the biosynthesis of chlorophyll, in order to identify energy transport pathways and to elucidate their role as potential allosteric communication networks for coordinating functional dynamics, specifically the negative cooperativity observed in the motion of the two active site gating loops. Fully atomistic MD simulations of thermal diffusion were executed with a GROMACS simulation package on a fully solvated GSAM enzyme by heating various active site target ligands (initially, catalytic intermediates and cofactors) to 300K while holding the remainder of the protein and the solvent bath at 10K and monitoring the temperature T(t) of all the enzyme residues as a function of time over a 1ns observation window. Energy is observed to be deposited in a relatively small number of discrete chains of residues most of which contribute to specific structural or biochemical functionality. Thermal linkages between all thermally active chains were established by isolating a specific pair of chains and performing a thermal diffusion simulation on the pair, one held at 300K and the other at 10K, with the rest of the protein frozen in its initial atomic configuration and thus thermally unresponsive. Proceeding in this way, it was possible to map out multiple pathways of vibrational energy flow leading from one of the active sites through a network of contiguous residues, many of which were evolutionarily conserved and linked by hydrogen bonds, into the other active site and ultimately to the other gating loop.

Role of atomic contacts in vibrational energy transfer in myoglobin

Heme proteins are ideal systems to investigate vibrational energy flow at the atomic level. Upon photoexcitation, a large amount of excess vibrational energy is selectively deposited on heme due to extremely fast internal conversion. This excess energy is redistributed to the surrounding protein moiety and then to water. Vibrational energy flow in myoglobin (Mb) was examined using picosecond time-resolved anti-Stokes ultraviolet resonance Raman (UVRR) spectroscopy. We used the Trp residue directly contacting the heme group as a selective probe for vibrationally excited populations. Trp residues were placed at different position close to the heme by site-directed mutagenesis. This technique allows us to monitor the excess energy on residue-to-residue basis. Anti-Stokes UVRR measurements for Mb mutants suggest that the dominant channel for energy transfer in Mb is the pathway through atomic contacts between heme and nearby amino acid residues as well as that between the protein and solvent water. It is found that energy flow through proteins is analogous to collisional exchange processes in solutions.

Recent developments in the computational study of protein structural and vibrational energy dynamics

Recent developments in the computational study of energy transport in proteins are reviewed, including advances in both methodology and applications. The concept of energy exchange network (EEN) is discussed, and a recent calculation of EENs for the allosteric protein FixL is reviewed, which illustrates how residues and protein regions involved in the allosteric transition can be identified. Recent work has examined relations between EENs and protein dynamics as well as structure. We review some of the computational studies carried out on several proteins that explore connections between energy conductivity across polar contacts in proteins and between proteins and water and equilibrium dynamics of the contacts, and we discuss some of the recent experimental work that addresses this topic.

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Protein Dynamics Preceding Photoisomerization of the Retinal Chromophore in Bacteriorhodopsin Revealed by Deep-UV Femtosecond Stimulated Raman Spectroscopy

Bacteriorhodopsin is a prototypical photoreceptor protein that functions as a light-driven proton pump. The retinal chromophore of bacteriorhodopsin undergoes C₁₃═C₁₄ trans-to-cis isomerization upon photoexcitation, and it has been believed to be the first event that triggers the cascaded structural changes in bacteriorhodopsin. We investigated the protein dynamics of bacteriorhodopsin using deep-ultraviolet resonance femtosecond stimulated Raman spectroscopy. It was found that the stimulated Raman signals of tryptophan and tyrosine residues exhibit significant changes within 0.2 ps after photoexcitation while they do not noticeably change during the isomerization process. This result implies that the protein environment changes first, and its change is small during isomerization. The obtained femtosecond stimulated Raman data indicate that ultrafast change is induced in the protein part by the sudden creation of the large dipole of the excited-state chromophore, providing an environment that realizes efficient and selective isomerization.

Langmuir Nanoarchitectonics from Basic to Frontier

Methodology to combine nanotechnology and these organization processes has been proposed as a novel concept of nanoarchitectonics, which can fabricate functional materials with nanolevel units. As an instant nanoarchitectonics approach, confining systems within a two-dimensional plane to drastically reduce translational motion freedom can be regarded as one of the rational approaches. Supramolecular chemistry and nanofabrication and their related functions at the air-water interface with the concept of nanoarchitectonics would lead to the creation of a novel methodology of Langmuir nanoarchitectonics. In this feature article, we briefly summarize research efforts related to Langmuir nanoarchitectonics including the basics for anomalies in molecular interactions such as highly enhanced molecular recognition capabilities. It is also extended to frontiers including the fabrication of supramolecular receptors and two-dimensional patterns with subnanometer-scale structural regulation, manual control of molecular machines and receptors by hand-motion-like macroscopic actions, and the regulation of cell fates at nanoarchitected arrays of nanocarbon assemblies and at direct liquid interfaces.

Materials nanoarchitectonics at two-dimensional liquid interfaces

Much attention has been paid to the synthesis of low-dimensional materials from small units such as functional molecules. Bottom-up approaches to create new low-dimensional materials with various functional units can be realized with the emerging concept of nanoarchitectonics. In this review article, we overview recent research progresses on materials nanoarchitectonics at two-dimensional liquid interfaces, which are dimensionally restricted media with some freedoms of molecular motion. Specific characteristics of molecular interactions and functions at liquid interfaces are briefly explained in the first parts. The following sections overview several topics on materials nanoarchitectonics at liquid interfaces, such as the preparation of two-dimensional metal-organic frameworks and covalent organic frameworks, and the fabrication of low-dimensional and specifically structured nanocarbons and their assemblies at liquid-liquid interfaces. Finally, interfacial nanoarchitectonics of biomaterials including the regulation of orientation and differentiation of living cells are explained. In the recent examples described in this review, various materials such as molecular machines, molecular receptors, block-copolymer, DNA origami, nanocarbon, phages, and stem cells were assembled at liquid interfaces by using various useful techniques. This review overviews techniques such as conventional Langmuir-Blodgett method, vortex Langmuir-Blodgett method, liquid-liquid interfacial precipitation, instructed assembly, and layer-by-layer assembly to give low-dimensional materials including nanowires, nanowhiskers, nanosheets, cubic objects, molecular patterns, supramolecular polymers, metal-organic frameworks and covalent organic frameworks. The nanoarchitecture materials can be used for various applications such as molecular recognition, sensors, photodetectors, supercapacitors, supramolecular differentiation, enzyme reactors, cell differentiation control, and hemodialysis.