Vibrational Energy Transfer from Heme through Atomic Contacts in Proteins

Iminothioindoxyl Donors with Exceptionally High Cross Section for Protein Vibrational Energy Transfer

Various protein functions are related to vibrational energy transfer (VET) as an important mechanism. The underlying transfer pathways can be experimentally followed by ultrafast Vis-pump/IR-probe spectroscopy with a donor-sensor pair of non-canonical amino acids (ncAAs) incorporated in a protein. However, so far only one donor ncAA, azulenylalanine (AzAla), exists, which suffers from a comparably low Vis extinction coefficient. Here, we introduce two novel donor ncAAs based on an iminothioindoxyl (ITI) chromophore. The dimethylamino-ITI (DMA-ITI) and julolidine-ITI (J-ITI) moieties overcome the limitation of AzAla with a 50 times higher Vis extinction coefficient. While ITI moieties are known for ultrafast photoswitching, DMA-ITI and J-ITI exclusively form a hot ground state on the sub-ps timescale instead, which is essential for their usage as vibrational energy donor. In VET measurements of donor-sensor dipeptides we investigate the performance of the new donors. We observe 20 times larger signals compared to the established AzAla donor, which opens unprecedented possibilities for the study of VET in proteins.

Noncovalent π Interactions in Mutated Aquomet-Myoglobin Proteins: A QM/MM and Local Vibrational Mode Study

Protein dynamics and function is strongly connected to the energy flow taking place. Myoglobin (Mb) and its mutations are ideal systems to study the process of vibrational energy transfer (VET) at the molecular level. Anti-Stokes ultraviolet resonance Raman studies using a tryptophan (Trp) probe, introduced at different Mb positions by amino acid replacement, have suggested that the amount of VET depends on the position of the Trp probe relative to the heme group. Inspired by this experimental work, we explored the strength of noncovalent π interactions, as well as covalent interactions for both the axial and distal ligands bound to iron in aquomet-Mb with the local vibrational mode analysis (LMA), originally developed by Konkoli and Cremer. Two sets of noncovalent interactions were investigated: (1) the interaction between the water ligand and Trp rings and (2) the interaction between the Trp and the porphyrin rings of the heme group. We assessed the strength of these noncovalent interactions via a special local mode force constant. Various Trp-modified water-bound ferric Mb proteins in the ground state were studied (6 in total) using gas-phase and QM/MM calculations followed by LMA. Our results disclose that VET is indeed dependent on the position of the Trp probe relative to the heme group but also on the tautomeric nature of distal histidine. They provide new guidelines on how to assess noncovalent π interactions in proteins utilizing LMA and how to use these data to explore VET, and more generally protein dynamics and function.

Course of Plasmodium infection studied using 2D-COS on human erythrocytes

BACKGROUND

The threat of malaria is still present in the world. Recognizing the type of parasite is important in determining a treatment plan. The golden routine involves microscopic diagnostics of Giemsa-stained thin blood smears, however, alternative methods are also constantly being sought, in order to gain an additional insight into the course of the disease. Spectroscopic methods, e.g., Raman spectroscopy, are becoming increasingly popular, due to the non-destructive nature of these techniques.

METHODS

The study included patients hospitalized for malaria caused by Plasmodium falciparum or Plasmodium vivax, in the Department of Infectious Diseases at the University Hospital in Krakow, Poland, as well as healthy volunteers. The aim of this study was to assess the possibility of using Raman spectroscopy and 2D correlation (2D-COS) spectroscopy in understanding the structural changes in erythrocytes depending on the type of attacking parasite. EPR spectroscopy and two-trace two-dimensional (2T2D) correlation was also used to examine the specificity of paramagnetic centres found in the infected human blood.

RESULTS

Two-dimensional (2D) correlation spectroscopy facilitates the identification of the hidden relationship, allowing for the discrimination of Raman spectra obtained during the course of disease in human red blood cells, infected by P. falciparum or P. vivax. Synchronous cross-peaks indicate the processes taking place inside the erythrocyte during the export of the parasite protein towards the cell membrane. In contrast, moieties that generate asynchronous 2D cross-peaks are characteristic of the respective ligand-receptor domains. These changes observed during the course of the infection, have different dynamics for P. falciparum and P. vivax, as indicated by the asynchronous correlation cross-peaks. Two-trace two-dimensional (2T2D) spectroscopy, applied to EPR spectra of blood at the beginning of the infection, showed differences between P. falciparum and P. vivax.

CONCLUSIONS

A unique feature of 2D-COS is the ability to discriminate the collected Raman and EPR spectra. The changes observed during the course of a malaria infection have different dynamics for P. falciparum and P. vivax, indicated by the reverse sequence of events. For each type of parasite, a specific recycling process for iron was observed in the infected blood.

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.

Prediction of Allosteric Communication Pathways in Proteins

MOTIVATION

Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form.

RESULTS

Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Computational Study on the Thermal Conductivity of a Protein

Protein molecules are thermally fluctuating and tightly packed amino acid residues strongly interact with each other. Such interactions are characterized in terms of heat current at the atomic level. We calculated the thermal conductivity of a small globular protein, villin headpiece subdomain, based on the linear response theory using equilibrium molecular dynamics simulation. The value of its thermal conductivity was 0.3 ± 0.01 [W m^-1 K^-1], which is in good agreement with experimental and computational studies on the other proteins in the literature. Heat current along the main chain was dominated by local vibrations in the polypeptide bonds, with amide I, II, III, and A bands on the Fourier transform of the heat current autocorrelation function.

Energy Bilocalization Effect and the Emergence of Molecular Functions in Proteins

Proteins are among the most complex molecular structures, which have evolved to develop broad functions, such as energy conversion and transport, information storage and processing, communication, and regulation of chemical reactions. However, the mechanisms by which these dynamical entities coordinate themselves to perform biological tasks remain hotly debated. Here, a physical theory is presented to explain how functional dynamical behavior possibly emerge in complex/macro molecules, thanks to the effect that we term bilocalization of thermal vibrations. More specifically, our approach allows us to understand how structural irregularities lead to a partitioning of the energy of the vibrations into two distinct sets of molecular domains, corresponding to slow and fast motions. This shape-encoded spectral allocation, associated to the genetic sequence, provides a close access to a wide reservoir of dynamical patterns, and eventually allows the emergence of biological functions by natural selection. To illustrate our approach, the SPIKE protein structure of SARS-COV2 is considered.

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.

Through bonds or contacts? Mapping protein vibrational energy transfer using non-canonical amino acids

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.

Ultraviolet resonance Raman spectroscopy with a continuously tunable picosecond laser: Application to the supramolecular ligand guanidiniocarbonyl pyrrole (GCP)

We present a UVRR spectroscopy setup which is equipped with a picosecond pulsed laser excitation source continuously tunable in the 210-2600 nm wavelength range. This laser source is based on a three-stage optical parametric amplifier (OPA) pumped by a bandwidth-compressed second harmonic output of an amplified Yb:KGW laser. It provides <15 cm^-1 linewidth pulses below 270 nm, which is sufficient for resolving Raman lines of samples in condensed phase studies. For demonstrating the capability of this tunable setup for UVRR spectroscopy we present its application to the artificial ligand guanidiniocarbonyl pyrrole (GCP), a carboxylate binder used in peptide and protein recognition. A UVRR excitation study in the range 244-310 nm was performed for identifying the optimum laser excitation wavelength for UVRR spectroscopy of this ligand (λ_max = 298 nm) at submillimolar concentrations (400 µM) in aqueous solution. The optimum UVRR spectrum is observed for laser excitation with λ_exc = 266 nm. Only in the relatively narrow range of λ_exc = 266-275 nm UVRR spectra with a sufficiently high signal-to-noise ratio and without severe interference from autofluorescence (AF) were detectable. At longer excitation wavelengths the UVRR signal is masked by AF. At shorter excitation wavelengths the UVRR spectrum is sufficiently separated from the AF, but the resonance enhancement is not sufficient. The presented tunable UVRR setup provides the flexibility to also identify optimum conditions for other supramolecular ligands for peptide/protein recognition.

Mimicking Natural Photosynthesis: Designing Ultrafast Photosensitized Electron Transfer into Multiheme Cytochrome Protein Nanowires

Efficient nanomaterials for artificial photosynthesis require fast and robust unidirectional electron transfer (ET) from photosensitizers through charge-separation and accumulation units to redox-active catalytic sites. We explored the ultrafast time-scale limits of photo-induced charge transfer between a Ru(II)tris(bipyridine) derivative photosensitizer and PpcA, a 3-heme c-type cytochrome serving as a nanoscale biological wire. Four covalent attachment sites (K28C, K29C, K52C, and G53C) were engineered in PpcA enabling site-specific covalent labeling with expected donor-acceptor (DA) distances of 4-8 Å. X-ray scattering results demonstrated that mutations and chemical labeling did not disrupt the structure of the proteins. Time-resolved spectroscopy revealed three orders of magnitude difference in charge transfer rates for the systems with otherwise similar DA distances and the same number of covalent bonds separating donors and acceptors. All-atom molecular dynamics simulations provided additional insight into the structure-function requirements for ultrafast charge transfer and the requirement of van der Waals contact between aromatic atoms of photosensitizers and hemes in order to observe sub-nanosecond ET. This work demonstrates opportunities to utilize multi-heme c-cytochromes as frameworks for designing ultrafast light-driven ET into charge-accumulating biohybrid model systems, and ultimately for mimicking the photosynthetic paradigm of efficiently coupling ultrafast, light-driven electron transfer chemistry to multi-step catalysis within small, experimentally versatile photosynthetic biohybrid assemblies.

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.

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.

Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study

To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids β-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.

Site-Resolved Observation of Vibrational Energy Transfer Using a Genetically Encoded Ultrafast Heater

Allosteric information transfer in proteins has been linked to distinct vibrational energy transfer (VET) pathways in a number of theoretical studies. Experimental evidence for such pathways, however, is sparse because site-selective injection of vibrational energy into a protein, that is, localized heating, is required for their investigation. Here, we solved this problem by the site-specific incorporation of the non-canonical amino acid β-(1-azulenyl)-l-alanine (AzAla) through genetic code expansion. As an exception to Kasha's rule, AzAla undergoes ultrafast internal conversion and heating after S₁ excitation while upon S₂ excitation, it serves as a fluorescent label. We equipped PDZ3, a protein interaction domain of postsynaptic density protein 95, with this ultrafast heater at two distinct positions. We indeed observed VET from the incorporation sites in the protein to a bound peptide ligand on the picosecond timescale by ultrafast IR spectroscopy. This approach based on genetically encoded AzAla paves the way for detailed studies of VET and its role in a wide range of proteins.

Thermodynamic Scaling of Interfering Hemoglobin Strain Field Waves

Hemoglobin (Hgb) forms tetramers (dimerized α-β dimers), which enhance its globular stability and may also facilitate small gas molecule transport, as shown by recent all-atom Newtonian solvated simulations. Hydropathic bioinformatic thermodynamic scaling enables close comparisons of hemoglobin dimers with myoglobin and neuroglobin, and reveals many nonlocal wave-like features of strained Hgb structures at the coarse-grained amino acid level. The thermodynamic analysis employs two hydropathic scales, one describing abrupt first-order unfolding transitions, the other continuous second-order transitions. Small molecule exchange at hemes is a first-order process. Wave-like collective tetrameric features appropriate to ligand absorption and release, seen in optical experiments (short times), are identified thermodynamically at long times. Strain fields localized near hemes interfere with extended strain fields associated with dimer interfacial misfit, resulting in novel wavelength dependent dimer correlation function Fano antiresonances.