Nitrile and thiocyanate IR probes: Quantum chemistry calculation studies and multivariate least-square fitting analysis

Strategy for Modeling the Electrostatic Responses of the Spectroscopic Properties of Proteins

For better understanding and more efficient use of the spectroscopic probes (vibrational and NMR) of the local electrostatic situations inside proteins, appropriate modeling of the properties of those probes is essential. The present study is devoted to examining the strategy for constructing such models. A more well-founded derivation than the ones in previous studies is given in constructing the models. Theoretical analyses are conducted on two representative example cases related to proteins, i.e., the peptide group of the main chains and the CO and NO ligands to the Fe²⁺ ion of heme, with careful treatment of the behavior of electrons in the electrostatic responses and with verification of consistency with observable quantities. It is shown that, for the stretching frequencies and NMR chemical shifts, it is possible to construct reasonable electrostatic interaction models that encompass the situations of hydration and uniform electric field environment and thus are applicable also to the cases of nonuniform electrostatic situations, which are highly expected for inside of proteins.

Ultrafast dynamics of ionic liquids in colloidal dispersion

Ionic liquid (IL)-surfactant complexes have significance both in applications and fundamental research, but their underlying dynamics are not well understood. We apply polarization-controlled two-dimensional infrared spectroscopy (2D-IR) to study the dynamics of [BMIM][SCN]/surfactant/solvent model systems. We examine the effect of the choice of surfactants and solvent, and the IL-to-surfactant ratio (W-value), with a detailed analysis of the orientation and structural dynamics of each system. Different surfactants create very different environments for the entrapped ILs, ranging from a semi-static micro-environment to a fluxional environment that evolves even faster than the bulk IL. The oil-phase also clearly affects the microscopic dynamics. The anisotropy decay for entrapped ILs completes within 10 ps, which is similar to free thiocyanate ion in water, while a significant reorientation-induced spectral diffusion (RISD) effect is observed. The entrapped ionic liquid are highly dynamic for all W-values, and no core-shell structure is observed. We hypothesize that, instead of an ionic liquid-reverse micelle (IL-RM), the microscopic structure of this system is small colloidal dispersions or pairs of IL and surfactants. A detailed analysis of the polarization-controlled 2D-IR spectra of AOT system reveals a potential ion-exchange mechanism.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

Electrostatic frequency maps for amide-I mode of β-peptide: Comparison of molecular mechanics force field and DFT calculations

The spectroscopy of amide-I vibrations has been widely utilized for the understanding of dynamical structure of polypeptides. For the modeling of amide-I spectra, two frequency maps were built for β-peptide analogue (N-ethylpropionamide, NEPA) in a number of solvents within different schemes (molecular mechanics force field based, GM map; DFT calculation based, GD map), respectively. The electrostatic potentials on the amide unit that originated from solvents and peptide backbone were correlated to the amide-I frequency shift from gas phase to solution phase during map parameterization. GM map is easier to construct with negligible computational cost since the frequency calculations for the samples are purely based on force field, while GD map utilizes sophisticated DFT calculations on the representative solute-solvent clusters and brings insight into the electronic structures of solvated NEPA and its chemical environments. The results show that the maps' predicted amide-I frequencies present solvation environmental sensitivities and exhibit their specific characters with respect to the map protocols, and the obtained vibrational parameters are in satisfactory agreement with experimental amide-I spectra of NEPA in solution phase. Although different theoretical schemes based maps have their advantages and disadvantages, the present maps show their potentials in interpreting the amide-I spectra for β-peptides, respectively.

Vibrational Probes: From Small Molecule Solvatochromism Theory and Experiments to Applications in Complex Systems

The vibrational frequency of a chosen normal mode is one of the most accurately measurable spectroscopic properties of molecules in condensed phases. Accordingly, infrared absorption and Raman scattering spectroscopy have provided valuable information on both distributions and ensemble-average values of molecular vibrational frequencies, and these frequencies are now routinely used to investigate structure, conformation, and even absolute configuration of chemical and biological molecules of interest. Recent advancements in coherent time-domain nonlinear vibrational spectroscopy have allowed the study of heterogeneous distributions of local structures and thermally driven ultrafast fluctuations of vibrational frequencies. To fully utilize IR probe functional groups for quantitative bioassays, a variety of biological and chemical techniques have been developed to site-specifically introduce vibrational probe groups into proteins and nucleic acids. These IR-probe-labeled biomolecules and chemically reactive systems are subject to linear and nonlinear vibrational spectroscopic investigations and provide information on the local electric field, conformational changes, site-site protein contacts, and/or function-defining features of biomolecules. A rapidly expanding library of data from such experiments requires an interpretive method with atom-level chemical accuracy. However, despite prolonged efforts to develop an all-encompassing theory for describing vibrational solvatochromism and electrochromism as well as dynamic fluctuations of instantaneous vibrational frequencies, purely empirical and highly approximate theoretical models have often been used to interpret experimental results. They are, in many cases, based on the simple assumption that the vibrational frequency of an IR reporter is solely dictated by electric potential or field distribution around the vibrational chromophore. Such simplified description of vibrational solvatochromism generally referred to as vibrational Stark effect theory has been considered to be quite appealing and, even in some cases, e.g., carbonyl stretch modes in amide, ester, ketone, and carbonate compounds or proteins, it works quantitatively well, which makes it highly useful in determining the strength of local electric field around the IR chromophore. However, noting that the vibrational frequency shift results from changes of solute-solvent intermolecular interaction potential along its normal coordinate, Pauli exclusion repulsion, polarization, charge transfer, and dispersion interactions, in addition to the electrostatic interaction between distributed charges of both vibrational chromophore and solvent molecules, are to be properly included in the theoretical description of vibrational solvatochromism. Since the electrostatic and nonelectrostatic intermolecular interaction components have distinctively different distance and orientation dependences, they affect the solvatochromic vibrational properties in a completely different manner. Over the past few years, we have developed a systematic approach to simulating vibrational solvatochromic data based on the effective fragment potential approach, one of the most accurate and rigorous theories on intermolecular interactions. We have further elucidated the interplay of local electric field with the general vibrational solvatochromism of small IR probes in either solvents or complicated biological systems, with emphasis on contributions from non-Coulombic intermolecular interactions to vibrational frequency shifts and fluctuations. With its rigorous foundation and close relation to quantitative interpretation of experimental data, this and related theoretical approaches and experiments will be of use in studying and quantifying the structure and dynamics of biomolecules with unprecedented time and spatial resolution when combined with time-resolved vibrational spectroscopy and chemically sensitive vibrational imaging techniques.

Extended timescale 2D IR probes of proteins: p-cyanoselenophenylalanine

The importance of dynamics to the function of proteins is well appreciated, but the difficulty in their measurement impedes investigation into their precise role(s). 2D IR spectroscopy is a developing approach for the study of dynamics and has motivated efforts to develop spectrally resolved IR probe groups that enable its application for measuring the dynamics at specific sites in a protein. A challenge with this approach is that the timescales accessible are limited by the vibrational lifetimes of the probes. Toward development of better probes for 2D IR spectroscopy of protein dynamics, we report the characterization of p-cyano-seleno-phenylalanine (CNSePhe), a derivative of the well established IR probe p-cyano-phenylalanine (CNPhe), by FT IR, pump-probe, and 2D IR spectroscopy. The incorporation of the heavy Se atom decouples the CN vibration from the rest in the molecule. Although this leads to a reduction of the transition dipole strength, and thus a reduction in signal intensity, it also dramatically increases the vibrational lifetime, enabling collection of 2D IR spectra for analysis of molecular dynamics on much longer timescales. Interestingly, we also find that the lifetime for CNSePhe shows increased sensitivity to the presence of hydrogen bonding interactions with the CN, suggesting that the probe should be useful for interpretation of CN spectra and possibly for the study of solvation.

A water-mediated allosteric network governs activation of Aurora kinase A

The catalytic activity of many protein kinases is controlled by conformational changes of a conserved Asp-Phe-Gly (DFG) motif. We used an infrared probe to track the DFG motif of the mitotic kinase Aurora A (AurA) and found that allosteric activation by the spindle-associated protein Tpx2 involves an equilibrium shift toward the active DFG-in state. Förster resonance energy transfer experiments show that the activation loop undergoes a nanometer-scale movement that is tightly coupled to the DFG equilibrium. Tpx2 further activates AurA by stabilizing a water-mediated allosteric network that links the C-helix to the active site through an unusual polar residue in the regulatory spine. The polar spine residue and water network of AurA are essential for phosphorylation-driven activation, but an alternative form of the water network found in related kinases can support Tpx2-driven activation, suggesting that variations in the water-mediated hydrogen bond network mediate regulatory diversification in protein kinases.

Vibrational Stark Effects of Carbonyl Probes Applied to Reinterpret IR and Raman Data for Enzyme Inhibitors in Terms of Electric Fields at the Active Site

IR and Raman frequency shifts have been reported for numerous probes of enzyme transition states, leading to diverse interpretations. In the case of the model enzyme ketosteroid isomerase (KSI), we have argued that IR spectral shifts for a carbonyl probe at the active site can provide a connection between the active site electric field and the activation free energy (Fried et al. Science 2014, 346, 1510-1514). Here we generalize this approach to a much broader set of carbonyl probes (e.g., oxoesters, thioesters, and amides), first establishing the sensitivity of each probe to an electric field using vibrational Stark spectroscopy, vibrational solvatochromism, and MD simulations, and then applying these results to reinterpret data already in the literature for enzymes such as 4-chlorobenzoyl-CoA dehalogenase and serine proteases. These results demonstrate that the vibrational Stark effect provides a general framework for estimating the electrostatic contribution to the catalytic rate and may provide a metric for the design or modification of enzymes. Opportunities and limitations of the approach are also described.

Quantum mechanical calculation of electric fields and vibrational Stark shifts at active site of human aldose reductase

Recent advance in biophysics has made it possible to directly measure site-specific electric field at internal sites of proteins using molecular probes with C = O or C≡N groups in the context of vibrational Stark effect. These measurements directly probe changes of electric field at specific protein sites due to, e.g., mutation and are very useful in protein design. Computational simulation of the Stark effect based on force fields such as AMBER and OPLS, while providing good insight, shows large errors in comparison to experimental measurement due to inherent difficulties associated with point charge based representation of force fields. In this study, quantum mechanical calculation of protein's internal electrostatic properties and vibrational Stark shifts was carried out by using electrostatically embedded generalized molecular fractionation with conjugate caps method. Quantum calculated change of mutation-induced electric field and vibrational Stark shift is reported at the internal probing site of enzyme human aldose reductase. The quantum result is in much better agreement with experimental data than those predicted by force fields, underscoring the deficiency of traditional point charge models describing intra-protein electrostatic properties.

Accurate Calculation of Electric Fields Inside Enzymes

The specific electric field generated by a protease at its active site is considered as an important source of the catalytic power. Accurate calculation of electric field at the active site of an enzyme has both fundamental and practical importance. Measuring site-specific changes of electric field at internal sites of proteins due to, eg, mutation, has been realized by using molecular probes with CO or CN groups in the context of vibrational Stark effect. However, theoretical prediction of change in electric field inside a protein based on a conventional force field, such as AMBER or OPLS, is often inadequate. For such calculation, quantum chemical approach or quantum-based polarizable or polarized force field is highly preferable. Compared with the result from conventional force field, significant improvement is found in predicting experimentally measured mutation-induced electric field change using quantum-based methods, indicating that quantum effect such as polarization plays an important role in accurate description of electric field inside proteins. In comparison, the best theoretical prediction comes from fully quantum mechanical calculation in which both polarization and inter-residue charge transfer effects are included for accurate prediction of electrostatics in proteins.

Synthesis and Evaluation of the Sensitivity and Vibrational Lifetimes of Thiocyanate and Selenocyanate Infrared Reporters

Two novel 2'-deoxyadenosine (dA) analogues, Si₂-dA-SCN and Si₂-dA-SeCN, and two novel phenylalanine (Phe) analogues, Boc-Me-PheCH₂SCN and Boc-Me-PheCH₂SeCN, have been synthesized and the thiocyanate (SCN) and selenocyanate (SeCN) functional groups evaluated as vibrational reporters. The syntheses of Si₂-dA-SCN and Si₂-dA-SeCN were accomplished in three steps in 16% and 32% overall yields, respectively, and the syntheses of Boc-Me-PheCH₂SCN and Boc-Me-PheCH₂SeCN were completed in four steps in 8.9% and 2.3% overall yields, respectively. The SCN and SeCN stretch vibrational modes were shown to be sensitive to the local environment by frequency shifts and full-width half-maximum (fwhm) changes in response to tetrahydrofuran (THF) and THF/water solvent mixtures. The vibrational lifetimes of the Si₂-dA-SeCN (237±12 ps) and Boc-Me-PheCH₂SeCN (295±31 ps) in THF solution were determined by ultrafast infrared pump-probe spectroscopy to be 1.5 to 3 times longer than those for Si₂-dA-SCN (140±6 ps) and Boc-Me-PheCH₂SCN (102±4 ps). The longer lifetimes for the SeCN analogues were attributed to the better insulating effects of the heavier selenium atom compared to the sulfur atom. The solvent sensitivity and longer vibrational lifetimes compared to other vibrational reporters suggest that SCN and SeCN vibrational reporters are well suited to studying several dynamic processes including protein and nucleic acid hydration and conformational changes, however stability issues may require post-synthetic modification methods to incorporate these reporters into biomacromolecules.

Roles of the scalar and vector components of the solvation effects on the vibrational properties of hydrogen- or halogen-bond accepting stretching modes

Solvation-induced vibrational frequency shifts and infrared (IR) intensity changes of the hydrogen- or halogen-bond accepting stretching modes, especially their dependence on the angular position of the hydrogen- or halogen-bond donating molecule, are examined theoretically. Calculations are carried out for some modes of hydrogen- or halogen-bonding molecular complexes, including the S[double bond, length as m-dash]O stretch of dimethyl sulfoxide-(13)C2H2O, the C[triple bond, length as m-dash]N stretch of acetonitrileH2O, and the amide I' mode of the N-methylacetamide-d1BrNC 1 : 1 complex. It is shown that, in all the example cases dealt with in this study, the frequency shift depends rather strongly on the hydrogen- or halogen-bond angle (e.g., S[double bond, length as m-dash]OH angle), with a larger low-frequency shift as the hydrogen or halogen bond becomes more bent, indicating the generality of the results obtained for the amide I' mode of the N-methylacetamide-d1(2)H2O 1 : 1 complex in a previous study. Contrary to our vague expectation, the frequency shift is not well correlated to the hydrogen- or halogen-bond distance or strength, but nevertheless, it is well reproduced by an electrostatic interaction model if it is carefully constructed by considering the scalar and vector components separately in a reasonable way. On the basis of this electrostatic interaction model, the reason why our vague expectation is not realized is clarified, and a unified understanding is achieved on the hydration-induced high-frequency shift of the C[triple bond, length as m-dash]N stretch and the low-frequency shifts of the S[double bond, length as m-dash]O stretch and amide I'. With regard to the IR intensity, it is shown that, in some of the example cases, it also has rather strong angular position dependence. The mechanism of the IR intensity changes is estimated by analyzing the dipole derivative vector, especially its angular relation with the hydrogen or halogen bond.

General Applicable Frequency Map for the Amide-I Mode in β-Peptides

In this work, a general applicable amide-I vibrational frequency map (GA map) for β-peptides in a number of common solvents was constructed, based on a peptide derivative, N-ethylpropionamide (NEPA). The map utilizes force fields at the ab initio computational level to accurately describe molecular structure and solute-solvent interactions, and also force fields at the molecular mechanics level to take into account long-range solute-solvent interactions. The results indicate that the GA map works reasonably for mapping the vibrational frequencies of the amide-I local-modes for β-peptides, holding promises for understanding the complicated infrared spectra of the amide-I mode in β-polypeptides.

Vibrational solvatochromism of nitrile infrared probes: beyond the vibrational Stark dipole approach

Systematic probing of local environments around biopolymers is important for understanding their functions.

Understanding and Manipulating Electrostatic Fields at the Protein-Protein Interface Using Vibrational Spectroscopy and Continuum Electrostatics Calculations

Biological function emerges in large part from the interactions of biomacromolecules in the complex and dynamic environment of the living cell. For this reason, macromolecular interactions in biological systems are now a major focus of interest throughout the biochemical and biophysical communities. The affinity and specificity of macromolecular interactions are the result of both structural and electrostatic factors. Significant advances have been made in characterizing structural features of stable protein-protein interfaces through the techniques of modern structural biology, but much less is understood about how electrostatic factors promote and stabilize specific functional macromolecular interactions over all possible choices presented to a given molecule in a crowded environment. In this Feature Article, we describe how vibrational Stark effect (VSE) spectroscopy is being applied to measure electrostatic fields at protein-protein interfaces, focusing on measurements of guanosine triphosphate (GTP)-binding proteins of the Ras superfamily binding with structurally related but functionally distinct downstream effector proteins. In VSE spectroscopy, spectral shifts of a probe oscillator's energy are related directly to that probe's local electrostatic environment. By performing this experiment repeatedly throughout a protein-protein interface, an experimental map of measured electrostatic fields generated at that interface is determined. These data can be used to rationalize selective binding of similarly structured proteins in both in vitro and in vivo environments. Furthermore, these data can be used to compare to computational predictions of electrostatic fields to explore the level of simulation detail that is necessary to accurately predict our experimental findings.

Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone

We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.

Ultrafast vibrational spectroscopy (2D-IR) of CO2 in ionic liquids: Carbon capture from carbon dioxide's point of view

The CO2ν3 asymmetric stretching mode is established as a vibrational chromophore for ultrafast two-dimensional infrared (2D-IR) spectroscopic studies of local structure and dynamics in ionic liquids, which are of interest for carbon capture applications. CO2 is dissolved in a series of 1-butyl-3-methylimidazolium-based ionic liquids ([C4C1im][X], where [X](-) is the anion from the series hexafluorophosphate (PF6 (-)), tetrafluoroborate (BF4 (-)), bis-(trifluoromethyl)sulfonylimide (Tf2N(-)), triflate (TfO(-)), trifluoroacetate (TFA(-)), dicyanamide (DCA(-)), and thiocyanate (SCN(-))). In the ionic liquids studied, the ν3 center frequency is sensitive to the local solvation environment and reports on the timescales for local structural relaxation. Density functional theory calculations predict charge transfer from the anion to the CO2 and from CO2 to the cation. The charge transfer drives geometrical distortion of CO2, which in turn changes the ν3 frequency. The observed structural relaxation timescales vary by up to an order of magnitude between ionic liquids. Shoulders in the 2D-IR spectra arise from anharmonic coupling of the ν2 and ν3 normal modes of CO2. Thermal fluctuations in the ν2 population stochastically modulate the ν3 frequency and generate dynamic cross-peaks. These timescales are attributed to the breakup of ion cages that create a well-defined local environment for CO2. The results suggest that the picosecond dynamics of CO2 are gated by local diffusion of anions and cations.

Vibrational solvatochromism. II. A first-principle theory of solvation-induced vibrational frequency shift based on effective fragment potential method

Vibrational solvatochromism is a solvation-induced effect on fundamental vibrational frequencies of molecules in solutions. Here we present a detailed first-principle coarse-grained theory of vibrational solvatochromism, which is an extension of our previous work [B. Błasiak, H. Lee, and M. Cho, J. Chem. Phys. 139(4), 044111 (2013)] by taking into account electrostatic, exchange-repulsion, polarization, and charge-transfer interactions. By applying our theory to the model N-methylacetamide-water clusters, solute-solvent interaction-induced effects on amide I vibrational frequency are fully elucidated at Hartree-Fock level. Although the electrostatic interaction between distributed multipole moments of solute and solvent molecules plays the dominant role, the contributions from exchange repulsion and induced dipole-electric field interactions are found to be of comparable importance in short distance range, whereas the charge-transfer effect is negligible. The overall frequency shifts calculated by taking into account the contributions of electrostatics, exchange-repulsion, and polarization terms are in quantitative agreement with ab initio results obtained at the Hartree-Fock level of theory.

Site-specific infrared probes of proteins

Infrared spectroscopy has played an instrumental role in the study of a wide variety of biological questions. However, in many cases, it is impossible or difficult to rely on the intrinsic vibrational modes of biological molecules of interest, such as proteins, to reveal structural and environmental information in a site-specific manner. To overcome this limitation, investigators have dedicated many recent efforts to the development and application of various extrinsic vibrational probes that can be incorporated into biological molecules and used to site-specifically interrogate their structural or environmental properties. In this review, we highlight recent advancements in this rapidly growing research area.

Site-selective Characterization of Src Homology 3 Domain Molecular Recognition with Cyanophenylalanine Infrared Probes

Local heterogeneity of microenvironments in proteins is important in biological function, but difficult to characterize experimentally. One approach is the combination of infrared (IR) spectroscopy and site-selective incorporation of probe moieties with spectrally resolved IR absorptions that enable characterization within inherently congested protein IR spectra. We employed this method to study molecular recognition of a Src homology 3 (SH3) domain from the yeast protein Sho1 for a peptide containing the proline-rich recognition sequence of its physiological binding partner Pbs2. Nitrile IR probes were introduced at four distinct sites in the protein by selective incorporation of p-cyanophenylalanine via the amber codon suppressor method and characterized by IR spectroscopy. Variation among the IR absorption bands reports on heterogeneity in local residue environments dictated by the protein structure, as well as on residue-dependent changes upon peptide binding. The study informs on the molecular recognition of SH3 ^Sho1 and illustrates the speed and simplicity of this approach for characterization of select microenvironments within proteins.

Solvent effect on the vibrational spectrum of Michler's ketone. Experimental and theoretical investigations

We examined solvent effect on the IR and Raman spectra of MK in several solvents of different polarity and proticity, for understanding of intermolecular interactions, focusing on solvent effect in detail. It has been found that change of solvent polarity has an ambiguous influence on solvatochromism of MK. We have observed that not only vibrations of carbonyl group are affected by the solvent polarity, but also mode ν(CN) and ν(CC) in IR and Raman spectra of MK. Experimental investigations have been supported by the quantum-mechanical computations to gain more insight into the solvatochromic behavior of Michler's ketone. Calculations have been carried using Kohn-Sham formulation of Density Functional Theory (DFT) and the Polarizable Continuum Model (PCM) was employed to account for solute solvent interactions.

Modeling biophysical and biological properties from the characteristics of the molecular electron density, electron localization and delocalization matrices, and the electrostatic potential

The electron density and the electrostatic potential are fundamentally related to the molecular hamiltonian, and hence are the ultimate source of all properties in the ground- and excited-states. The advantages of using molecular descriptors derived from these fundamental scalar fields, both accessible from theory and from experiment, in the formulation of quantitative structure-to-activity and structure-to-property relationships, collectively abbreviated as QSAR, are discussed. A few such descriptors encode for a wide variety of properties including, for example, electronic transition energies, pK(a)'s, rates of ester hydrolysis, NMR chemical shifts, DNA dimers binding energies, π-stacking energies, toxicological indices, cytotoxicities, hepatotoxicities, carcinogenicities, partial molar volumes, partition coefficients (log P), hydrogen bond donor capacities, enzyme-substrate complementarities, bioisosterism, and regularities in the genetic code. Electronic fingerprinting from the topological analysis of the electron density is shown to be comparable and possibly superior to Hammett constants and can be used in conjunction with traditional bulk and liposolubility descriptors to accurately predict biological activities. A new class of descriptors obtained from the quantum theory of atoms in molecules' (QTAIM) localization and delocalization indices and bond properties, cast in matrix format, is shown to quantify transferability and molecular similarity meaningfully. Properties such as "interacting quantum atoms (IQA)" energies which are expressible into an interaction matrix of two body terms (and diagonal one body "self" terms, as IQA energies) can be used in the same manner. The proposed QSAR-type studies based on similarity distances derived from such matrix representatives of molecular structure necessitate extensive investigation before their utility is unequivocally established.

Vibrational solvatochromism: towards systematic approach to modeling solvation phenomena

Vibrational solvatochromic frequency shift of IR probe is an effect of interaction between local electric field and IR probe in condensed phases. Despite prolonged efforts to develop empirical maps for vibrational frequency shifts and transition dipoles of IR probes, a systematic approach to ab initio calculation of vibrational solvatochromic charges and multipoles has not been developed. Here, we report on density functional theory (DFT) calculations of N-methylacetamide (NMA) frequency shifts using implicit and coarse-grained models. The solvatochromic infrared spectral shifts are estimated based on the distributed multipole analysis of electronic densities calculated for gas-phase equilibrium structure of NMA. Thus obtained distributed solvatochromic multipole parameters are used to calculate the amide I vibrational frequency shifts of NMA in water clusters that mimic the instantaneous configurations of the liquid water. Our results indicate that the spectral shifts are primarily electrostatic in nature and can be quantitatively reproduced using the proposed model with semi-quantitative accuracy when compared to the corresponding DFT results.

A conserved water-mediated hydrogen bond network defines bosutinib's kinase selectivity

Kinase inhibitors are important cancer drugs, but they tend to display limited target specificity, and their target profiles are often challenging to rationalize in terms of molecular mechanism. Here we report that the clinical kinase inhibitor bosutinib recognizes its kinase targets by engaging a pair of conserved structured water molecules in the active site and that many other kinase inhibitors share a similar recognition mechanism. Using the nitrile group of bosutinib as an infrared probe, we show that the gatekeeper residue and one other position in the ATP-binding site control access of the drug to the structured water molecules and that the amino acids found at these positions account for the kinome-wide target spectrum of the drug. Our work highlights the importance of structured water molecules for inhibitor recognition, reveals a new role for the kinase gatekeeper and showcases an effective approach for elucidating the molecular origins of selectivity patterns.

Molecular mechanics force field-based general map for the solvation effect on amide I probe of peptide in different micro-environments

A general electrostatic potential map based on molecular mechanics force field for modeling the amide I frequency is presented. This map is applied to N-methylacetamide (NMA) and designed to be transferable in different micro-environments. The electrostatic potentials from solvent and peptide side chain are projected on the amide unit of NMA to induce the frequency shift of amide I mode. It is shown that the predicted amide I frequency reproduces the experimental data satisfactorily, especially when NMA in polar solvents. The amide I frequency shift is largely determined by the solvents in aqueous solution while it is dominated by the local structure of peptide in other solvent environments. The map parameters are further applied on NMA-MeOH system and the obtained IR spectra show doublet peak profile with negligible deviation from the experimental data, suggesting the usefulness of this general map for providing information about vibrational parameters of amide motions of peptide in different environments.

Conserved electrostatic fields at the Ras–effector interface measured through vibrational Stark effect spectroscopy explain the difference in tilt angle in the Ras binding domains of Raf and RalGDS

Vibrational Stark effect (VSE) spectroscopy was used to measure the electrostatic fields present at the interface of the human guanosine triphosphatase (GTPase) Ras docked with the Ras binding domain (RBD) of the protein kinase Raf.

Vibrational dynamics and solvatochromism of the label SCN in various solvents and hemoglobin by time dependent IR and 2D-IR spectroscopy

The vibrational label SCN is used to report on local structural dynamics in a protein revealing spectral diffusion on a picosecond scale. The SCN spectra are compared to the response of methylthiocyanate in solvents with different polarity and hydrogen-bonding capabilities.

Aqueous Mg2+ and Ca2+ ligand exchange mechanisms identified with 2DIR spectroscopy

Biological systems must discriminate between calcium and magnesium for these ions to perform their distinct biological functions, but the mechanism for distinguishing aqueous ions has yet to be determined. Ionic recognition depends upon the rate and mechanism by which ligands enter and leave the first solvation shell surrounding these cations. We present a time-resolved vibrational spectroscopy study of these ligand exchange dynamics in aqueous solution. The sensitivity of the CN-stretch frequency of NCS(-) to ion pair formation has been utilized to investigate the mechanism and dynamics of ligand exchange into and out of the first solvation shell of aqueous magnesium and calcium ions with multidimensional vibrational (2DIR) spectroscopy. We have determined that anion exchange follows a dissociative mechanism for Mg(2+) and an associative mechanism for Ca(2+).