Instantaneous ion configurations in the K+ ion channel selectivity filter revealed by 2D IR spectroscopy

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Cholesterol modulates the structural dynamics of the paddle motif loop of KvAP voltage sensor

KvAP is a prokaryotic Kv channel, which has been widely used as a model system to understand voltage- and lipid-dependent gating mechanisms. In phospholipid membranes, the KvAP-VSD adopts the activated/'Up' conformation, whereas the presence of non-phospholipids in membranes favours the structural transition to resting/'Down' state. The S3b-S4 paddle motif loop of KvAP-VSD is functionally important as this participates in protein-protein interactions and is the target for animal toxins. In this study, we have monitored the modulatory role of cholesterol - the physiologically-relevant non-phospholipid - on the organization and dynamics of the S3b-S4 loop of the isolated KvAP-VSD in membranes by site-directed fluorescence approaches using the environmental sensitivity of 7-nitrobenz-2-oxa-1,3-diazol-4-yl-ethylenediamine (NBD) fluorescence. Our results show that cholesterol alters the dynamic nature (rotational and hydration dynamics) of S3b-S4 loop in a segmental fashion, i.e., the residues 110 to 114 and 115 to 117 behave differently in the presence of cholesterol, which is accompanied by considerable change in conformational heterogeneity. Further, quantitative depth measurements using the parallax quenching method reveal that the sensor loop is located at the shallow interfacial region of cholesterol-containing membranes, suggesting that the sensor loop organization is not directly correlated with S4 helix movement. Our results clearly show that cholesterol-induced changes in bilayer properties may not be the predominant factor for the sensor loop's altered structural dynamics, but can be attributed to the conformational change of the KvAP-VSD in cholesterol-containing membranes. Overall, these results are relevant for gating mechanisms, particularly the lipid-dependent gating, of Kv channels in membranes.

Insight into the K channel's selectivity from binding of K+, Na+ and water to N-methylacetamide

In potassium channels that conduct K+ selectively over Na+, which sites are occupied by K+ or water and the mechanism of selectivity are unresolved questions. The combination of the energetics and the constraints imposed by the protein structure yield the selective permeation and occupancy. To gain insight into the combination of structure and energetics, we performed density functional theory (DFT) calculations of multiple N-methyl acetamide (NMA) ligands binding to K+ and Na+, relative to hydrated K+ and Na+. NMA is an analogue of the amino acid backbone and provides the carbonyl binding to the ions that occurs in most binding sites of the K+ channel. Unconstrained optimal structures are obtained through geometry optimization calculations of the NMA ligand binding. The complexes formed by 8 NMA binding to the cations have the O atoms positioned in nearly identical locations as the O atoms in the selectivity filter. The transfer free energies between bulk water and K+ or Na+ bound to 8 NMA are almost identical, implying there is no selectivity by a single site. For water optimized with 8 NMA, binding is weak and O atoms are not positioned as in the K+ channel selectivity filter, suggesting that the ions are much more favored than water. Optimal structures of 8 NMA binding with two cations (K+ or Na+) are stable and have lower binding free energy than the optimal structures with just one cation. However, in the Na+ case, the optimal structure deforms and does not match the K+ channel; that is, two bound Na+ are destabilizing. In contrast, the two K+ structure is stabilized and the selectivity free energy favors K+. Overall, this study shows that binding site occupancy and the mechanism for K+ selectivity involves multiple K+ binding in multiple neighboring layers or sites of the K+ channel selectivity filter.

Artificial Intelligence-based Amide-II Infrared Spectroscopy Simulation for Monitoring Protein Hydrogen Bonding Dynamics

The structurally sensitive amide II infrared (IR) bands of proteins provide valuable information about the hydrogen bonding of protein secondary structures, which is crucial for understanding protein dynamics and associated functions. However, deciphering protein structures from experimental amide II spectra relies on time-consuming quantum chemical calculations on tens of thousands of representative configurations in solvent water. Currently, the accurate simulation of amide II spectra for whole proteins remains a challenge. Here, we present a machine learning (ML)-based protocol designed to efficiently simulate the amide II IR spectra of various proteins with an accuracy comparable to experimental results. This protocol stands out as a cost-effective and efficient alternative for studying protein dynamics, including the identification of secondary structures and monitoring the dynamics of protein hydrogen bonding under different pH conditions and during protein folding process. Our method provides a valuable tool in the field of protein research, focusing on the study of dynamic properties of proteins, especially those related to hydrogen bonding, using amide II IR spectroscopy.

Water inside the Selectivity Filter of a K+ Ion Channel: Structural Heterogeneity, Picosecond Dynamics, and Hydrogen Bonding

Water inside biological ion channels regulates the key properties of these proteins, such as selectivity, ion conductance, and gating. In this article, we measure the picosecond spectral diffusion of amide I vibrations of an isotope-labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100-2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope-labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D line shapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent or nonadjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations was observed on a picosecond timescale. These dynamics are in stark contrast with liquid water, which remains highly dynamic even in the presence of ions at high concentrations.

Femtosecond Electron-Transfer Dynamics across the D2O/Cs+/Cu(111) Interface: The Impact of Hydrogen Bonding

Hydrogen bonding is essential in electron-transfer processes at water-electrode interfaces. We study the impact of the H-bonding of water as a solvent molecule on real-time electron-transfer dynamics across a Cs+-Cu(111) ion-metal interface using femtosecond time-resolved two-photon photoelectron spectroscopy. We distinguish in the formed water-alkali aggregates two regimes below and above two water molecules per ion. Upon crossing the boundary of these regimes, the lifetime of the excess electron localized transiently at the Cs+ ion increases from 40 to 60 fs, which indicates a reduced alkali-metal interaction. Furthermore, the energy transferred to a dynamic structural rearrangement due to hydration is reduced from 0.3 to 0.2 eV concomitantly. These effects are a consequence of H-bonding in the water-water interaction and the beginning formation of a nanoscale water network. This finding is supported by real-space imaging of the solvatomers and vibrational frequency shifts of the OH stretching and bending modes calculated for these specific interfaces.

Water inside the selectivity filter of a K+ ion channel: structural heterogeneity, picosecond dynamics, and hydrogen-bonding

Water inside biological ion channels regulates the key properties of these proteins such as selectivity, ion conductance, and gating. In this Article we measure the picosecond spectral diffusion of amide I vibrations of an isotope labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100 - 2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D lineshapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent, or non-adjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations were observed on a picosecond timescale. These dynamics are in stark contrast with liquid water that remains highly dynamic even in the presence of ions at high concentrations.

Recent advances and research progress on microsystems and bioeffects of terahertz neuromodulation

Terahertz waves can interact with the nervous system of organisms under certain conditions. Compared to common optical modulation methods, terahertz waves have the advantages of low photon energy and low risk; therefore, the use of terahertz waves to regulate the nervous system is a promising new method of neuromodulation. However, most of the research has focused on the use of terahertz technology for biodetection, while relatively little research has been carried out on the biological effects of terahertz radiation on the nervous system, and there are almost no review papers on this topic. In the present article, we begin by reviewing principles and objects of research regarding the biological effects of terahertz radiation and summarizing the current state of related research from a variety of aspects, including the bioeffects of terahertz radiation on neurons in vivo and in vitro, novel regulation and detection methods with terahertz radiation devices and neural microelectrode arrays, and theoretical simulations of neural information encoding and decoding. In addition, we discuss the main problems and their possible causes and give some recommendations on possible future breakthroughs. This paper will provide insight and assistance to researchers in the fields of neuroscience, terahertz technology and biomedicine.

Two-Dimensional Infrared Spectroscopy of Isolated Molecular Ions

Two-dimensional infrared (2D IR) spectroscopy of mass-selected, cryogenically cooled molecular ions is presented. Nonlinear response pathways, encoded in the time-domain photodissociation action response of weakly bound N2 messenger tags, were isolated using pulse shaping techniques following excitation with four collinear ultrafast IR pulses. 2D IR spectra of Re(CO)3(CH3CN)3+ ions capture off-diagonal cross-peak bleach signals between the asymmetric and symmetric carbonyl stretching transitions. These cross peaks display intensity variations as a function of pump-probe delay time due to coherent coupling between the vibrational modes. Well-resolved 2D IR features in the congested fingerprint region of protonated caffeine (C8H10N4O2H+) are also reported. Importantly, intense cross-peak signals were observed at 3 ps waiting time, indicating that tag-loss dynamics are not competing with the measured nonlinear signals. These demonstrations pave the way for more precise studies of molecular interactions and dynamics that are not easily obtainable with current condensed-phase methodologies.

Structural Dependence of Extended Amide III Vibrations in Two-Dimensional Infrared Spectra

Two-dimensional infrared (2D-IR) spectroscopy is a powerful experimental method for probing the structure and dynamics of proteins in aqueous solution. So far, most experimental studies have focused on the amide I vibrations, for which empirical vibrational exciton models provide a means of interpreting such experiments. However, such models are largely lacking for other regions of the vibrational spectrum. To close this gap, we employ an efficient quantum-chemical methodology for the calculation of 2D-IR spectra, which is based on anharmonic theoretical vibrational spectroscopy with localized modes. We apply this approach to explore the potential of 2D-IR spectroscopy in the extended amide III region. Using calculations for a dipeptide model as well as alanine polypeptides, we show that distinct 2D-IR cross-peaks in the extended amide III region can potentially be used to distinguish α-helix and β-strand structures. We propose that the extended amide III region could be a promising target for future 2D-IR experiments.

Stepwise hydration of [CH3COOMg]+ studied by cold ion trap infrared spectroscopy: insights into interactions in the magnesium channel selection filters

The magnesium channel controls Mg2+ concentration in the cell and plays an indispensable role in biological functions. The crystal structure of the Magnesium Transport E channel suggested that Mg2+ hydrated by 6 water molecules is transported through a selection filter consisting of COO- groups on two Asp residues. This Mg2+ motion implies successive pairing with -OOC-R and dissociation mediated by water molecules. For another divalent ion, however, it is known that RCOO-⋯Ca2+ cannot be separated even with 12 water molecules. From this discrepancy, we probe the structure of Mg2+(CH3COO-)(H2O)4-17 clusters by measuring the infrared spectra and monitoring the vibrational frequencies of COO- with the help of quantum chemistry calculations. The hydration by (H2O)6 is not enough to induce ion separation, and partially-separated or separated pairs are formed from 10 water molecules at least. These results suggest that the ion separation between Mg2+ and carboxylate ions in the selection-filter of the MgtE channel not only results from water molecules in their first hydration shell, but also from additional factors including water molecules and protein groups in the second solvation shell of Mg2+.

Probing Ion Configurations in the KcsA Selectivity Filter with Single-Isotope Labels and 2D IR Spectroscopy

The potassium ion (K+) configurations of the selectivity filter of the KcsA ion channel protein are investigated with two-dimensional infrared (2D IR) spectroscopy of amide I vibrations. Single 13C-18O isotope labels are used, for the first time, to selectively probe the S1/S2 or S2/S3 binding sites in the selectivity filter. These binding sites have the largest differences in ion occupancy in two competing K+ transport mechanisms: soft-knock and hard-knock. According to the former, water molecules alternate between K+ ions in the selectivity filter while the latter assumes that K+ ions occupy the adjacent sites. Molecular dynamics simulations and computational spectroscopy are employed to interpret experimental 2D IR spectra. We find that in the closed conductive state of the KcsA channel, K+ ions do not occupy adjacent binding sites. The experimental data is consistent with simulated 2D IR spectra of soft-knock ion configurations. In contrast, the simulated spectra for the hard-knock ion configurations do not reproduce the experimental results. 2D IR spectra of the hard-knock mechanism have lower frequencies, homogeneous 2D lineshapes, and multiple peaks. In contrast, ion configurations of the soft-knock model produce 2D IR spectra with a single peak at a higher frequency and inhomogeneous lineshape. We conclude that under equilibrium conditions, in the absence of transmembrane voltage, both water and K+ ions occupy the selectivity filter of the KcsA channel in the closed conductive state. The ion configuration is central to the mechanism of ion transport through potassium channels.

Inhibition of vibrational energy flow within an aromatic scaffold via heavy atom effect

The regulation of intramolecular vibrational energy redistribution (IVR) to influence energy flow within molecular scaffolds provides a way to steer fundamental processes of chemistry, such as chemical reactivity in proteins and design of molecular diodes. Using two-dimensional infrared (2D IR) spectroscopy, changes in the intensity of vibrational cross-peaks are often used to evaluate different energy transfer pathways present in small molecules. Previous 2D IR studies of para-azidobenzonitrile (PAB) demonstrated that several possible energy pathways from the N3 to the cyano-vibrational reporters were modulated by Fermi resonance, followed by energy relaxation into the solvent [Schmitz et al., J. Phys. Chem. A 123, 10571 (2019)]. In this work, the mechanisms of IVR were hindered via the introduction of a heavy atom, selenium, into the molecular scaffold. This effectively eliminated the energy transfer pathway and resulted in the dissipation of the energy into the bath and direct dipole-dipole coupling between the two vibrational reporters. Several structural variations of the aforementioned molecular scaffold were employed to assess how each interrupted the energy transfer pathways, and the evolution of 2D IR cross-peaks was measured to assess the changes in the energy flow. By eliminating the energy transfer pathways through isolation of specific vibrational transitions, through-space vibrational coupling between an azido (N3) and a selenocyanato (SeCN) probe is facilitated and observed for the first time. Thus, the rectification of this molecular circuitry is accomplished through the inhibition of energy flow using heavy atoms to suppress the anharmonic coupling and, instead, favor a vibrational coupling pathway.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

Resolving the Impact of Hydrogen Bonding on the Phylloquinone Cofactor through Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy (2DIR) was applied to phylloquinone (PhQ), an important biological cofactor, to elucidate the impact of hydrogen bonding on the ultrafast dynamics and energetics of the carbonyl stretching modes. 2DIR measurements were performed on PhQ dissolved in hexanol, which served as the hydrogen bonding solvent, and hexane, which served as a non-hydrogen bonding control. Molecular dynamics simulations and quantum chemical calculations were performed to aid in spectral assignment and interpretation. From the position of the peaks in the 2DIR spectra, we extracted the transition frequencies for the fundamental, overtone, and combination bands of hydrogen bonded and non-hydrogen bonded carbonyl groups of PhQ in the 1635-1680 cm^-1 region. We find that hydrogen bonding to a single carbonyl group acts to decouple the two carbonyl units of PhQ. Through analysis of the time-resolved 2DIR data, we find that hydrogen bonding leads to faster vibrational relaxation as well as an increase in the inhomogeneous broadening of the carbonyl groups. Overall, this work demonstrates how hydrogen bonding to the carbonyl groups of PhQ presents in the 2DIR spectra, laying the groundwork to use PhQ as a 2DIR probe to characterize the ultrafast fluctuations in the local environment of natural photosynthetic complexes.

Ultrafast Two-Dimensional Infrared Spectroscopy Resolved a Structured Lysine 159 on the Cytoplasmic Surface of the Microbial Photoreceptor Bacteriorhodopsin

Bacteriorhodopsin (bR) is a light-driven microbial receptor, and lysine 159 (K159) is a charged residue on the cytoplasmic (CP) side of its E-F loop. However, its conformation and function remain unknown due to fast surface dynamics. By utilizing a ¹³C, ¹⁵N-labeled lysine (K) as an isotope probe, we created a network of site-specific amide-I vibrational signatures (backbone carbonyl stretch) to identify the frequency contribution of the labeled residues to the amide-I excitonic band structure. Thus, the red-shifted amide-I frequency in the ¹³C, ¹⁵N-lysine-labeled bR (uK-bR) to the unlabeled bR (WT-bR) could be differentiated and examined by ultrafast two-dimensional vibrational echo infrared (2D IR) spectroscopy. Our results showed that the backbone carbonyl of K159 is located at a high frequency of ca. 1693 cm^-1 and has a vibrational excited-state relaxation time shorter than the bulk helical amide-I mode at the same frequency, suggesting that K159 may possess a hydrogen-bonded γ-turn structure with E161, one of the carboxylate residues on the CP surface of bR. The 2D solid-state NMR study of uK-bR also revealed conformational dependent lysine residues, from which K159 was found to involve the turn motif. This γ-turn structure maintained by K159 may help to stabilize the E-F loop and support E161 in attracting protons from the bulk during the late stage of the bR photocycle. The combined spectroscopic approach illustrated in this work may be applied to map residue-specific local structures and dynamics of other receptors and large proteins.

A bottom-up approach to the ion recognition mechanism of K+ channels from laser spectroscopy of hydrated partial peptide-alkali metal ion complexes

K⁺ channels allow selective permeation of K⁺, but not physiologically abundant Na⁺, at almost diffusion limit rates. The conduction mechanism of K⁺ channels is still controversial, with experimental and computation studies supporting two distinct conduction mechanisms: either with or without water inside the channel. Here, we employ a bottom-up approach on hydrated alkali metal complexes of a model peptide of K⁺ channels, Ac-Tyr-NHMe, to characterize metal-peptide, metal-water, and water-peptide interactions that govern the selectivity of K⁺ channels at a molecular level. Both the extension to the series of alkali metal ions and to temperature-dependent studies (approaching physiological values) have revealed the clear difference between permeable and non-permeable ions in the spectral features of the ion complexes. Furthermore, the impact of hydration is discussed in relation to the K⁺ channels by comparisons of the non-hydrated and hydrated complexes.

Measuring Protein Conformation at Aqueous Interfaces with 2D Infrared Spectroscopy of Emulsions

Determining the secondary and tertiary structures of proteins at aqueous interfaces is crucial for understanding their function, but measuring these structures selectively at the interface is challenging. Here we demonstrate that two-dimensional infrared (2D-IR) spectroscopy of protein stabilized emulsions offers a new route to measuring interfacial protein structure with high levels of detail. We prepared hexadecane/water oil-in-water emulsions stabilized by model LK peptides that are known to fold into either α-helix or β-sheet conformations at hydrophobic interfaces and measured 2D-IR spectra in a transmission geometry. We saw clear spectral signatures of the peptides folding at the interface, with no detectable residue from remaining bulk peptides. Using 2D spectroscopy gives us access to correlation and dynamics data, which enables structural assignment in cases where linear spectroscopy fails. Using the emulsions allows one to study interfacial spectra with standard transmission geometry spectrometers, bringing the richness of 2D-IR to the interface with no additional optical complexity.

A Multi-Scale Approach to Model K+ Permeation Through the KcsA Channel

K+ channels allow a very efficient passage of K+ ions through the membrane while excluding Na+ ions, and these properties are essential for life. The 3D structure of the KcsA K+ channel, solved more than 20 years ago, allows to address many relevant aspects of K+ permeation and selectivity mechanisms at the molecular level. Recent crystallographic data and molecular dynamics (MD) studies suggest that no water is normally present inside the selectivity filter (SF), which can instead accommodate four adjacent K+ ions. Using a multi-scale approach, whereby information taken from a low-level simulation approach is used to feed a high-level model, we studied the mechanism of K+ permeation through KcsA channels. More specifically, we used MD to find stable ion configurations under physiological conditions. They were characterized by two adjacent K+ ions occupying the more central positions of the SF (sites S2 and S3), while the other two K+ ions could be found at the external and internal entrances to the SF. Sites S1 and S4 were instead not occupied by K+. A continuum Bikerman–Poisson–Boltzmann model that takes into account the volume of the ions and their dehydration when entering the SF fully confirmed the MD results, showing peaks of K+ occupancy at S2, S3, and the external and internal entrances, with S1 and S4 sites being virtually never occupied by K+. Inspired by the newly found ion configuration in the SF at equilibrium, we developed a simple kinetic permeation model which, fed with kinetic rate constants assessed from molecular meta-dynamics, reproduced the main permeation properties of the KcsA channel found experimentally, including sublinear current-voltage and saturating conductance-concentration relationships. This good agreement with the experimental data also implies that the ion configuration in the SF we identified at equilibrium would also be a key configuration during permeation.

Quantum coherence on selectivity and transport of ion channels

Recently, it has been suggested that ion channel selectivity filter may exhibit quantum coherence, which may be appropriate to explain ion selection and conduction processes. Potassium channels play a vital role in many physiological processes. One of their main physiological functions is the efficient and highly selective transfer of K⁺ ions through the membranes into the cells. To do this, ion channels must be highly selective, allowing only certain ions to pass through the membrane, while preventing the others. The present research is an attempt to investigate the relationship between hopping rate and maintaining coherence in ion channels. Using the Lindblad equation to describe a three-level system, the results in different quantum regimes are examined. We studied the distillable coherence and the second order coherence function of the system. The oscillation of distillable coherence from zero, after the decoherence time, and also the behavior of the coherence function clearly show the point that the system is coherent in ion channels with high throughput rates.

Epithelial Sodium Channel Inhibition by Amiloride Addressed with THz Spectroscopy and Molecular Modeling

THz spectroscopy is important for the study of ion channels because it directly addresses the low frequency collective motions relevant for their function. Here we used THz spectroscopy to investigate the inhibition of the epithelial sodium channel (ENaC) by its specific blocker, amiloride. Experiments were performed on A6 cells' suspensions, which are cells overexpressing ENaC derived from Xenopus laevis kidney. THz spectra were investigated with or without amiloride. When ENaC was inhibited by amiloride, a substantial increase in THz absorption was noticed. Molecular modeling methods were used to explain the observed spectroscopic differences. THz spectra were simulated using the structural models of ENaC and ENaC-amiloride complexes built here. The agreement between the experiment and the simulations allowed us to validate the structural models and to describe the amiloride dynamics inside the channel pore. The amiloride binding site validated using THz spectroscopy agrees with previous mutagenesis studies. Altogether, our results show that THz spectroscopy can be successfully used to discriminate between native and inhibited ENaC channels and to characterize the dynamics of channels in the presence of their specific antagonist.

Machine learning recognition of protein secondary structures based on two-dimensional spectroscopic descriptors

Protein secondary structure discrimination is crucial for understanding their biological function. It is not generally possible to invert spectroscopic data to yield the structure. We present a machine learning protocol which uses two-dimensional UV (2DUV) spectra as pattern recognition descriptors, aiming at automated protein secondary structure determination from spectroscopic features. Accurate secondary structure recognition is obtained for homologous (97%) and nonhomologous (91%) protein segments, randomly selected from simulated model datasets. The advantage of 2DUV descriptors over one-dimensional linear absorption and circular dichroism spectra lies in the cross-peak information that reflects interactions between local regions of the protein. Thanks to their ultrafast (∼200 fs) nature, 2DUV measurements can be used in the future to probe conformational variations in the course of protein dynamics.

Site-directed fluorescence approaches to monitor the structural dynamics of proteins using intrinsic Trp and labeled with extrinsic fluorophores

Comprehensive understanding of a protein’s function depends on having reliable, sophisticated tools to study protein structural dynamics in physiologically-relevant conditions. Here, we present an effective, robust step-by-step protocol to monitor the structural dynamics (including hydration dynamics) of a protein utilizing various site-directed fluorescence (SDFL) approaches. This protocol should be widely applicable for studying soluble proteins, intrinsically-disordered proteins, and membrane proteins.

For complete details on the use and execution of this protocol, please refer to Das et al. (2020), Das and Raghuraman (2021), and Chatterjee et al. (2021).

•

A step-by-step protocol to monitor the structural dynamics of proteins using SDFL

•

Applicable to proteins with intrinsic Trp and labeled with extrinsic fluorophores

•

This protocol should be widely applicable for soluble and membrane proteins

Direct Detection of Bound Ammonium Ions in the Selectivity Filter of Ion Channels by Solid-State NMR

The flow of ions across cell membranes facilitated by ion channels is an important function for all living cells. Despite the huge amount of structural data provided by crystallography, elucidating the exact interactions between the selectivity filter atoms and bound ions is challenging. Here, we detect bound ¹⁵N-labeled ammonium ions as a mimic for potassium ions in ion channels using solid-state NMR under near-native conditions. The non-selective ion channel NaK showed two ammonium peaks corresponding to its two ion binding sites, while its potassium-selective mutant NaK2K that has a signature potassium-selective selectivity filter with four ion binding sites gave rise to four ammonium peaks. Ions bound in specific ion binding sites were identified based on magnetization transfer between the ions and carbon atoms in the selectivity filters. Magnetization transfer between bound ions and water molecules revealed that only one out of four ions in the selectivity filter of NaK2K is in close contact with water, which is in agreement with the direct knock-on ion conduction mechanism where ions are conducted through the channel by means of direct interactions without water molecules in between. Interestingly, the potassium-selective ion channels investigated here (NaK2K and, additionally, KcsA-Kv1.3) showed remarkably different chemical shifts for their bound ions, despite having identical amino acid sequences and crystal structures of their selectivity filters. Molecular dynamics simulations show similar ion binding and conduction behavior between ammonium and potassium ions and identify the origin of the differences between the investigated potassium channels.

The bare necessities of plant K+ channel regulation

Potassium (K+) channels serve a wide range of functions in plants from mineral nutrition and osmotic balance to turgor generation for cell expansion and guard cell aperture control. Plant K+ channels are members of the superfamily of voltage-dependent K+ channels, or Kv channels, that include the Shaker channels first identified in fruit flies (Drosophila melanogaster). Kv channels have been studied in depth over the past half century and are the best-known of the voltage-dependent channels in plants. Like the Kv channels of animals, the plant Kv channels are regulated over timescales of milliseconds by conformational mechanisms that are commonly referred to as gating. Many aspects of gating are now well established, but these channels still hold some secrets, especially when it comes to the control of gating. How this control is achieved is especially important, as it holds substantial prospects for solutions to plant breeding with improved growth and water use efficiencies. Resolution of the structure for the KAT1 K+ channel, the first channel from plants to be crystallized, shows that many previous assumptions about how the channels function need now to be revisited. Here, I strip the plant Kv channels bare to understand how they work, how they are gated by voltage and, in some cases, by K+ itself, and how the gating of these channels can be regulated by the binding with other protein partners. Each of these features of plant Kv channels has important implications for plant physiology.

Computational spectroscopy of complex systems

Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.

Ultrafast Structure and Vibrational Dynamics of a Cyano-Containing Non-Fullerene Acceptor for Organic Solar Cells Revealed by Two-Dimensional Infrared Spectroscopy

Non-fullerene molecules, such as ITIC (3,9-bis(2-methylene-(3-(1,1-dicyanomethylene) indanone)-5,5,11,11-tetrakis(4-hexylphenyl)dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']-dithiophene), are among the most promising non-fullerene acceptors for organic solar cells (OSCs). In this work, using the cyano stretching mode as a vibrational marker, the structural and vibrational energy dynamics of ITIC were examined on an ultrafast time scale with two-dimensional infrared spectroscopy. Two IR-active modes studied here mainly correspond to two anti-symmetric combinations of symmetric and asymmetric stretching vibrations of two C≡N modes originating from two -C(CN)₂ chromophores that are located across the ITIC system, which were found to have significantly larger off-diagonal anharmonicity than their corresponding diagonal anharmonicities. This indicates strong anharmonic vibrational coupling between the two modes, which is supported by ab initio anharmonic frequency computations. Transient IR results indicate picosecond intramolecular vibrational energy transfer between the two C≡N modes upon excitation. The structural basis for these vibrational and energetic features is the conjugating molecular frame that is composed of a network of single/double bonds connecting the two -C(CN)₂ chromophores and may enable efficient vibration delocalization, in addition to its well-known electron delocalization capability. The importance of the results for the OSC applications is discussed.

Transport threshold in a quantum model for the KscA ion channel

The mechanism behind the high throughput rate in K⁺channels is still an open problem. However, recent simulations have shown that the passage of potassium through the K⁺channel core, the so-called selectivity filter (SF), is water-free against models where the strength of Coulomb repulsion freezes ions conduction. Thus, it has been suggested that coherent quantum hopping might be relevant in mediating ion conduction. Within the quantum approach and the hypothesis of desolvated ions along the pathway, we start with several particles in a source to see how they go across a SF, modeled by a linear chain of sites, to be collected in a drain. We show that the average SF occupancy is three ions, and the ion transfer rate is ∼10⁸ions s^-1, results which agree with the recent findings in the literature.

Ion-dependent structure, dynamics, and allosteric coupling in a non-selective cation channel

The selectivity filter (SF) determines which ions are efficiently conducted through ion channel pores. NaK is a non-selective cation channel that conducts Na⁺ and K⁺ with equal efficiency. Crystal structures of NaK suggested a rigid SF structure, but later solid-state NMR and MD simulations questioned this interpretation. Here, we use solution NMR to characterize how bound Na⁺ vs. K⁺ affects NaK SF structure and dynamics. We find that the extracellular end of the SF is flexible on the ps-ns timescale regardless of bound ion. On a slower timescale, we observe a structural change between the Na⁺ and K⁺-bound states, accompanied by increased structural heterogeneity in Na⁺. We also show direct evidence that the SF structure is communicated to the pore via I88 on the M2 helix. These results support a dynamic SF with multiple conformations involved in non-selective conduction. Our data also demonstrate allosteric coupling between the SF and pore-lining helices in a non-selective cation channel that is analogous to the allosteric coupling previously demonstrated for K⁺-selective channels, supporting the generality of this model.

NaK is a bacterial non-selective cation channel. Here, the authors use solution NMR to show that selectivity filter (SF) in NaK is dynamic, with structural differences between the Na+ and K + -bound states. The conformation of the SF is communicated to the pore-lining helices similarly as in the K + -selective channels.

Conductance selectivity of Na+ across the K+ channel via Na+ trapped in a tortuous trajectory

Ion selectivity of the potassium channel is crucial for regulating electrical activity in living cells; however, the mechanism underlying the potassium channel selectivity that favors large K⁺ over small Na⁺ remains unclear. Generally, Na⁺ is not completely excluded from permeation through potassium channels. Herein, the distinct nature of Na⁺ conduction through the prototypical KcsA potassium channel was examined. Single-channel current recordings revealed that, at a high Na⁺ concentration (200 mM), the channel was blocked by Na⁺, and this blocking was relieved at high membrane potentials, suggesting the passage of Na⁺ across the channel. At a 2,000 mM Na⁺ concentration, single-channel Na⁺ conductance was measured as one-eightieth of the K⁺ conductance, indicating that the selectivity filter allows substantial conduit of Na⁺ Molecular dynamics simulations revealed unprecedented atomic trajectories of Na⁺ permeation. In the selectivity filter having a series of carbonyl oxygen rings, a smaller Na⁺ was distributed off-center in eight carbonyl oxygen-coordinated sites as well as on-center in four carbonyl oxygen-coordinated sites. This amphipathic nature of Na⁺ coordination yielded a continuous but tortuous path along the filter. Trapping of Na⁺ in many deep free energy wells in the filter caused slow elution. Conversely, K⁺ is conducted via a straight path, and as the number of occupied K⁺ ions increased to three, the concerted conduction was accelerated dramatically, generating the conductance selectivity ratio of up to 80. The selectivity filter allows accommodation of different ion species, but the ion coordination and interactions between ions render contrast conduction rates, constituting the potassium channel conductance selectivity.

Protein Dynamics by Two-Dimensional Infrared Spectroscopy

Proteins function as ensembles of interconverting structures. The motions span from picosecond bond rotations to millisecond and longer subunit displacements. Characterization of functional dynamics on all spatial and temporal scales remains challenging experimentally. Two-dimensional infrared spectroscopy (2D IR) is maturing as a powerful approach for investigating proteins and their dynamics. We outline the advantages of IR spectroscopy, describe 2D IR and the information it provides, and introduce vibrational groups for protein analysis. We highlight example studies that illustrate the power and versatility of 2D IR for characterizing protein dynamics and conclude with a brief discussion of the outlook for biomolecular 2D IR.

Isotope-Edited Amide II Mode: A New Label for Site-Specific Vibrational Spectroscopy

Vibrational spectroscopy is a powerful tool used to analyze biological and chemical samples. However, in proteins, the most predominant peaks that arise from the backbone amide groups overlap one another, hampering site-specific analyses. Isotope editing has provided a robust, noninvasive approach to overcome this hurdle. In particular, the 1-¹³C═¹⁶O and 1-¹³C═¹⁸O labels that shift the amide I vibrational mode have enabled 1D- and 2D-IR spectroscopy to characterize proteins with excellent site-specific resolution. Herein, we expand the vibrational spectroscopy toolkit appreciably by introducing the 1-¹³C[Formula: see text]¹⁵N probe at specific locations along the protein backbone. A new, isotopically edited amide II peak is observed clearly in the spectra despite the presence of unlabeled modes arising from the rest of the protein. The experimentally determined shift of -30 cm^-1 is reproduced by DFT calculations providing further credence to the mode assignment. Since the amide II mode arises from different elements than the amide I mode, it affords molecular insights that are both distinct and complementary. Moreover, multiple labeling schemes may be used simultaneously, enhancing vibrational spectroscopy's ability to provide detailed molecular insights.

Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport

The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and provide an outlook for prioritizing strategic research directions.

Thermodynamics of ion binding and occupancy in potassium channels

Potassium channels modulate various cellular functions through efficient and selective conduction of K⁺ ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K⁺ ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations.

Polarization shapes the energy landscape of ion conduction in potassium channels.

Potassium and sodium ion complexes with a partial peptide of the selectivity filter in K+ channels studied by cold ion trap infrared spectroscopy: the effect of hydration

Potassium channels allow K+ to rapidly diffuse, while the selectivity filter (SF) actively blocks Na+. The presence of water in the SF during ion translocation remains under debate due to the experimental and computational challenges in characterizing the interactions between water, ions, and the SF. Our bottom-up approach has been applied to a system composed of a partial peptide of the SF (Ac-tyrosine-NHMe) with a metal ion and a single water molecule to probe these interactions. The IR photodissociation spectra of M+Ac-tyrosine-NHMe(H2O) (M = Na, K) combined with quantum chemical calculations revealed that the water molecule binding sites are ion-dependent. In addition, the ion-peptide distances are elongated significantly for the K+ complex in comparison to the Na+ complex by the addition of a single water molecule. This striking structural difference with the water molecule is discussed in relation to ion selectivity and translocation within the K+ channel.

Shaker-IR K+ channel gating in heavy water: Role of structural water molecules in inactivation

It has been reported earlier that the slow (C-type) inactivated conformation in K_v channels is stabilized by a multipoint hydrogen-bond network behind the selectivity filter. Furthermore, MD simulations revealed that structural water molecules are also involved in the formation of this network locking the selectivity filter in its inactive conformation. We found that the application of an extracellular, but not intracellular, solution based on heavy water (D₂O) dramatically slowed entry into the slow inactivated state in Shaker-IR mutants (T449A, T449A/I470A, and T449K/I470C, displaying a wide range of inactivation kinetics), consistent with the proposed effect of the dynamics of structural water molecules on the conformational stability of the selectivity filter. Alternative hypotheses capable of explaining the observed effects of D₂O were examined. Increased viscosity of the external solution mimicked by the addition of glycerol had a negligible effect on the rate of inactivation. In addition, the inactivation time constants of K⁺ currents in the outward and the inward directions in asymmetric solutions were not affected by a H₂O/D₂O exchange, negating an indirect effect of D₂O on the rate of K⁺ rehydration. The elimination of the nonspecific effects of D₂O on our macroscopic current measurements supports the hypothesis that the rate of structural water exchange at the region behind the selectivity filter determines the rate of slow inactivation, as proposed by molecular modeling.

Structural plasticity of the selectivity filter in a nonselective ion channel

The sodium potassium ion channel (NaK) is a nonselective ion channel that conducts both sodium and potassium across the cellular membrane. A new crystallographic structure of NaK reveals conformational differences in the residues that make up the selectivity filter between the four subunits that form the ion channel and the inner helix of the ion channel. The crystallographic structure also identifies a side-entry, ion-conduction pathway for Na+ permeation that is unique to NaK. NMR studies and molecular dynamics simulations confirmed the dynamical nature of the top part of the selectivity filter and the inner helix in NaK as also observed in the crystal structure. Taken together, these results indicate that the structural plasticity of the selectivity filter combined with the dynamics of the inner helix of NaK are vital for the efficient conduction of different ions through the non-selective ion channel of NaK.

Computational IR Spectroscopy of Insulin Dimer Structure and Conformational Heterogeneity

We have investigated the structure and conformational dynamics of insulin dimer using a Markov state model (MSM) built from extensive unbiased atomistic molecular dynamics simulations and performed infrared spectral simulations of the insulin MSM to describe how structural variation within the dimer can be experimentally resolved. Our model reveals two significant conformations to the dimer: a dominant native state consistent with other experimental structures of the dimer and a twisted state with a structure that appears to reflect a ∼55° clockwise rotation of the native dimer interface. The twisted state primarily influences the contacts involving the C-terminus of insulin's B chain, shifting the registry of its intermolecular hydrogen bonds and reorganizing its side-chain packing. The MSM kinetics predict that these configurations exchange on a 14 μs time scale, largely passing through two Markov states with a solvated dimer interface. Computational amide I spectroscopy of site-specifically ¹³C¹⁸O labeled amides indicates that the native and twisted conformation can be distinguished through a series of single and dual labels involving the B24F, B25F, and B26Y residues. Additional structural heterogeneity and disorder is observed within the native and twisted states, and amide I spectroscopy can also be used to gain insight into this variation. This study will provide important interpretive tools for IR spectroscopic investigations of insulin structure and transient IR kinetics experiments studying the conformational dynamics of insulin dimer.

The Persistent Question of Potassium Channel Permeation Mechanisms

Potassium channels play critical roles in many physiological processes, providing a selective permeation route for K⁺ ions in and out of a cell, by employing a carefully designed selectivity filter, evolutionarily conserved from viruses to mammals. The structure of the selectivity filter was determined at atomic resolution by x-ray crystallography, showing a tight coordination of desolvated K⁺ ions by the channel. However, the molecular mechanism of K⁺ ions permeation through potassium channels remains unclear, with structural, functional and computational studies often providing conflicting data and interpretations. In this review, we will present the proposed mechanisms, discuss their origins, and will critically assess them against all available data. General properties shared by all potassium channels are introduced first, followed by the introduction of two main mechanisms of ion permeation: soft and direct knock-on. Then, we will discuss critical computational and experimental studies that shaped the field. We will especially focus on molecular dynamics (MD) simulations, that provided mechanistic and energetic aspects of K⁺ permeation, but at the same time created long-standing controversies. Further challenges and possible solutions are presented as well.

The molecular pH-response mechanism of the plant light-stress sensor PsbS

Plants need to protect themselves from excess light, which causes photo-oxidative damage and lowers the efficiency of photosynthesis. Photosystem II subunit S (PsbS) is a pH sensor protein that plays a crucial role in plant photoprotection by detecting thylakoid lumen acidification in excess light conditions via two lumen-faced glutamates. However, how PsbS is activated under low-pH conditions is unknown. To reveal the molecular response of PsbS to low pH, here we perform an NMR, FTIR and 2DIR spectroscopic analysis of Physcomitrella patens PsbS and of the E176Q mutant in which an active glutamate has been replaced. The PsbS response mechanism at low pH involves the concerted action of repositioning of a short amphipathic helix containing E176 facing the lumen and folding of the luminal loop fragment adjacent to E71 to a 310-helix, providing clear evidence of a conformational pH switch. We propose that this concerted mechanism is a shared motif of proteins of the light-harvesting family that may control thylakoid inter-protein interactions driving photoregulatory responses.

A Genetically Encoded Two-Dimensional Infrared Probe for Enzyme Active-Site Dynamics

While two-dimensional infrared (2D-IR) spectroscopy is uniquely suitable for monitoring femtosecond (fs) to picosecond (ps) water dynamics around static protein structures, its utility for probing enzyme active-site dynamics is limited due to the lack of site-specific 2D-IR probes. We demonstrate the genetic incorporation of a novel 2D-IR probe, m-azido-L-tyrosine (N3Y) in the active-site of DddK, an iron-dependent enzyme that catalyzes the conversion of dimethylsulfoniopropionate to dimethylsulphide. Our results show that both the oxidation of active-site iron to Fe^III , and the addition of denaturation reagents, result in significant decrease in enzyme activity and active-site water motion confinement. As tyrosine residues play important roles, including as general acids and bases, and electron transfer agents in many key enzymes, the genetically encoded 2D-IR probe N3Y should be broadly applicable to investigate how the enzyme active-site motions at the fs-ps time scale direct reaction pathways to facilitating specific chemical reactions.

Lessons from combined experimental and theoretical examination of the FTIR and 2D-IR spectroelectrochemistry of the amide I region of cytochrome c

Amide I difference spectroscopy is widely used to investigate protein function and structure changes. In this article, we show that the common approach of assigning features in amide I difference signals to distinct secondary structure elements in many cases may not be justified. Evidence comes from Fourier transform infrared (FTIR) and 2D-IR spectroelectrochemistry of the protein cytochrome c in the amide I range, in combination with computational spectroscopy based on molecular dynamics (MD) simulations. This combination reveals that each secondary structure unit, such as an alpha-helix or a beta-sheet, exhibits broad overlapping contributions, usually spanning a large part of the amide I region, which in the case of difference absorption experiments (such as in FTIR spectroelectrochemistry) may lead to intensity-compensating and even sign-changing contributions. We use cytochrome c as the test case, as this small electron-transferring redox-active protein contains different kinds of secondary structure units. Upon switching its redox-state, the protein exhibits a different charge distribution while largely retaining its structural scaffold. Our theoretical analysis suggests that the change in charge distribution contributes to the spectral changes and that structural changes are small. However, in order to confidently interpret FTIR amide I difference signals in cytochrome c and proteins in general, MD simulations in combination with additional experimental approaches such as isotope labeling, the insertion of infrared labels to selectively probe local structural elements will be required. In case these data are not available, a critical assessment of previous interpretations of protein amide I 1D- and 2D-IR difference spectroscopy data is warranted.

Infrared spectroscopy probes ion binding geometries

Infrared (IR) spectroscopy is a well-established technique for probing the structure, behavior, and surroundings of molecules in their native environments. Its characteristics-most specifically high structural sensitivity, ready applicability to aqueous samples, and broad availability-make it a valuable enzymological technique, particularly for the interrogation of ion binding sites. While IR spectroscopy of the "garden variety" (steady state at room temperature with wild-type proteins) is versatile and powerful in its own right, the combination of IR spectroscopy with specialized experimental schemes for leveraging ultrafast time resolution, protein labeling, and other enhancements further extends this utility. This book chapter provides the fundamental physical background and literature context essential for harnessing IR spectroscopy in the general context of enzymology with specific focus on interrogation of ion binding. Studies of lanthanide ions binding to calmodulin are highlighted as illustrative examples of this process. Appropriate sample preparation, data collection, and spectral interpretation are discussed from a detail-oriented and practical perspective with the goal of facilitating the reader's rapid progression from reading words in a book to collecting and analyzing their own data in the lab.

Characterizing the Water Wire in the Gramicidin Channel Found by Monte Carlo Sampling Using Continuum Electrostatics and in Molecular Dynamics Trajectories with Conventional or Polarizable Force Fields

Water molecules play a key role in all biochemical processes. They help define the shape of proteins, and they are reactant or product in many reactions and are released as ligands are bound. They facilitate the transfer of protons through transmembrane proton channel, pump and transporter proteins. Continuum electrostatics (CE) force fields used by program Multiconformation CE (MCCE) capture electrostatic interactions in biomolecules with an implicit solvent, which captures the averaged solvent water equilibrium properties. Hybrid CE methods can use explicit water molecules within the protein surrounded by implicit solvent. These hybrid methods permit the study of explicit hydrogen bond networks within the protein and allow analysis of processes such as proton transfer reactions. Yet hybrid CE methods have not been rigorously tested. Here, we present an explicit treatment of water molecules in the Gramicidin A (gA) channel using MCCE and compare the resulting distributions of water molecules and key hydration features against those obtained with explicit solvent Molecular Dynamics (MD) simulations with the nonpolarizable CHARMM36 and polarizable Drude force fields. CHARMM36 leads to an aligned water wire in the channel characterized by a large absolute net water dipole moment; the MCCE and Drude analysis lead to a small net dipole moment as the water molecules change orientation within the channel. The correct orientation is not as yet known, so these calculations identify an open question.

Structure Changes of a Membrane Polypeptide under an Applied Voltage Observed with Surface-Enhanced 2D IR Spectroscopy

The structures of many membrane-bound proteins and polypeptides depend on the membrane potential. However, spectroscopically studying their structures under an applied field is challenging, because a potential is difficult to generate across more than a few bilayers. We study the voltage-dependent structures of the membrane-bound polypeptide, alamethicin, using a spectroelectrochemical cell coated with a rough, gold film to create surface plasmons. The plasmons sufficiently enhance the 2D IR signal to measure a single bilayer. The film is also thick enough to conduct current and thereby apply a potential. The 2D IR spectra resolve features from both 3₁₀- and α-helical structures and cross-peaks connecting the two. We observe changes in the peak intensity, not their frequencies, upon applying a voltage. A similar change occurs with pH, which is known to alter the angle of alamethicin relative to the surface normal. The spectra are modeled using a vibrational exciton Hamiltonian, and the voltage-dependent spectra are consistent with a change in angle of the 3₁₀- and α-helices in the membrane from 55 to 44°and from 31 to 60°, respectively. The 3₁₀- and α-helices are coupled by approximately 10 cm^-1. These experiments provide new structural information about alamethicin under a potential difference and demonstrate a technique that might be applied to voltage-gated membrane proteins and compared to molecular dynamics structures.

Polarization Effects in Water-Mediated Selective Cation Transport across a Narrow Transmembrane Channel

Despite the progress in modeling complex molecular systems of ever-increasing complexity, a quantitatively accurate computational treatment of ion permeation through narrow membrane channels remains challenging. An important factor to reach this goal is induced electronic polarization, which is likely to impact the permeation rate of small ions through narrow molecular pores. In this work, we extended the recently developed polarizable force field based on the classical Drude oscillators to assess the role of induced polarization effects on the energetics of sodium and potassium ion transport across the gramicidin A (gA) ion channel. The inclusion of induced polarization lowers barriers present in 1D potential of mean force (PMF) for cation permeation by ∼50% compared to those obtained with the additive force field. Conductance properties calculated with 1D PMFs from Drude simulations are in better agreement with experimental results. Polarization of single-file water molecules and protein atoms forming the narrow pore has a direct impact on the free-energy barriers and cation-specific solid-state NMR chemical shifts. Sensitivity analysis indicates that small changes to water-channel interactions can alter the free energy barrier for ion permeation. These results, illustrating polarization effects present in the complex electrostatic environment of the gA channel, have broad implications for revising proposed mechanisms of ion permeation and selectivity in a variety of ion channels.