Measuring electrostatic fields in both hydrogen-bonding and non-hydrogen-bonding environments using carbonyl vibrational probes

The covalently bound diazo group as an infrared probe for hydrogen bonding environments

The bandwidth or position of diazo groups can probe the polarizability and HBA, or density of HBD groups.

Insights into acid dissociation of HCl and HBr with internal electric fields

A critical electric field exerted by the solvent on the ionizable group leads to acid dissociation.

Display of functional proteins on supramolecular peptide nanofibrils using a split-protein strategy

The display of functional proteins on self-assembled peptide nanofibrils is accomplished by noncovalent attachment using a split-protein strategy.

A Critical Test of the Electrostatic Contribution to Catalysis with Noncanonical Amino Acids in Ketosteroid Isomerase

The vibrational Stark effect (VSE) has been used to measure the electric field in the active site of ketosteroid isomerase (KSI). These measured fields correlate with ΔG(⧧) in a series of conventional mutants, yielding an estimate for the electrostatic contribution to catalysis (Fried et al. Science 2014, 346, 1510-1513). In this work we test this result with much more conservative variants in which individual Tyr residues in the active site are replaced by 3-chlorotyrosine via amber suppression. The electric fields sensed at the position of the carbonyl bond involved in charge displacement during catalysis were characterized using the VSE, where the field sensitivity has been calibrated by vibrational Stark spectroscopy, solvatochromism, and MD simulations. A linear relationship is observed between the electric field and ΔG(⧧) that interpolates between wild-type and more drastic conventional mutations, reinforcing the evaluation of the electrostatic contribution to catalysis in KSI. A simplified model and calculation are developed to estimate changes in the electric field accompanying changes in the extended hydrogen-bond network in the active site. The results are consistent with a model in which the O-H group of a key active site tyrosine functions by imposing a static electrostatic potential onto the carbonyl bond. The model suggests that the contribution to catalysis from the active site hydrogen bonds is of similar weight to the distal interactions from the rest of the protein. A similar linear correlation was also observed between the proton affinity of KSI's active site and the catalytic rate, suggesting a direct connection between the strength of the H-bond and the electric field it exerts.

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

Potassium channels are responsible for the selective permeation of K+ ions across cell membranes. K+ ions permeate in single file through the selectivity filter, a narrow pore lined by backbone carbonyls that compose four K+ binding sites. Here, we report on the two-dimensional infrared (2D IR) spectra of a semisynthetic KcsA channel with site-specific heavy (13C18O) isotope labels in the selectivity filter. The ultrafast time resolution of 2D IR spectroscopy provides an instantaneous snapshot of the multi-ion configurations and structural distributions that occur spontaneously in the filter. Two elongated features are resolved, revealing the statistical weighting of two structural conformations. The spectra are reproduced by molecular dynamics simulations of structures with water separating two K+ ions in the binding sites, ruling out configurations with ions occupying adjacent sites.

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.

Pick and Choose the Spectroscopic Method to Calibrate the Local Electric Field inside Proteins

Electrostatic interactions in proteins play a crucial role in determining the structure-function relation in biomolecules. In recent years, fluorescent probes have been extensively employed to interrogate the polarity in biological cavities through dielectric constants or semiempirical polarity scales. A choice of multiple spectroscopic methods, not limited by fluorophores, along with a molecular level description of electrostatics involving solute-solvent interactions, would allow more flexibility to pick and choose the experimental technique to determine the local electrostatics within protein interiors. In this work we report that ultraviolet/visible-absorption, infrared-absorption, or (13)C NMR can be used to calibrate the local electric field in both hydrogen bonded and non-hydrogen bonded protein environments. The local electric field at the binding site of a serum protein has been determined using the absorption wavelength as well as the carbonyl stretching frequency of its natural steroid substrate, testosterone. Excellent agreement is observed in the results obtained from two independent spectroscopic techniques.

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.

Conformational Heterogeneity in the Michaelis Complex of Lactate Dehydrogenase: An Analysis of Vibrational Spectroscopy Using Markov and Hidden Markov Models

Lactate dehydrogenase (LDH) catalyzes the interconversion of pyruvate and lactate. Recent isotope-edited IR spectroscopy suggests that conformational heterogeneity exists within the Michaelis complex of LDH, and this heterogeneity affects the propensity toward the on-enzyme chemical step for each Michaelis substate. By combining molecular dynamics simulations with Markov and hidden Markov models, we obtained a detailed kinetic network of the substates of the Michaelis complex of LDH. The ensemble-average electric fields exerted onto the vibrational probe were calculated to provide a direct comparison with the vibrational spectroscopy. Structural features of the Michaelis substates were also analyzed on atomistic scales. Our work not only clearly demonstrates the conformational heterogeneity in the Michaelis complex of LDH and its coupling to the reactivities of the substates, but it also suggests a methodology to simultaneously resolve kinetics and structures on atomistic scales, which can be directly compared with the vibrational spectroscopy.

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.

Internal electric fields in small water clusters [(H2O)n; n = 2-6]

The electric field experienced by a water molecule within a water cluster depends on its position relative to the rest of the water molecules. The stabilization energies and the red-shifts in the donor O-H stretching vibrations in the water clusters increase with the cluster size concomitant with the increase in the electric field experienced by the donor O-H of a particular water molecule due to the hydrogen bonding network. The red-shifts in O-H stretching frequencies show a spread of about ±100 cm(-1) against the corresponding electric fields. Deviations from linearity were marked in the region of 100-160 MV cm(-1), which can be attributed to the strain in the hydrogen bonding network, especially for structures with DDAA and DDA motifs. The linear Stark effect holds up to 200 MV cm(-1) of internal electric field for the average red-shifts in the O-H stretching frequencies, with a Stark tuning rate of 2.4 cm(-1) (MV cm(-1))(-1), suggesting the validity of the classical model in small water clusters.

Electrostatic Interactions Are Key to C═O n-π* Shifts: An Experimental Proof

Carbonyl n-π* transitions are known to undergo blue shift in polar and hydrogen-bonding solvents. Using semiempirical expressions, previous studies hypothesized several factors like change in dipole moment and hydrogen-bond strength upon excitation to cause the blue shift. Theoretically, ground-state electrostatics has been predicted to be the key to the observed shifts, however, an experimental proof has been lacking. Our experimental results demonstrate a consistent linear correlation between IR (ground-state phenomenon) and n-π* frequency shifts (involves both ground and excited electronic-states) of carbonyls in hydrogen-bonded and non-hydrogen-bonded environments. The carbonyl hydrogen-bonding status is experimentally verified from deviation in n-π*/fluorescence correlation. The IR/n-π* correlation validates the key role of electrostatic stabilization of the ground state toward n-π* shifts and demonstrates the electrostatic nature of carbonyl hydrogen bonds. n-π* shifts show linear sensitivity to calculated electrostatic fields on carbonyls. Our results portray the potential for n-π* absorption to estimate local polarity in biomolecules and to probe chemical reactions involving carbonyl activation/stabilization.

Organized Molecular Interface-Induced Noncrystallizable Polymer Ultrathin Nanosheets with Ordered Chain Alignment

The orientation and state of organization of polymer chains play significant roles in determining the final properties and functions of these materials. Unlike most semicrystalline polymers, which have an inherent driving force toward crystallization, the means to control chain packing in noncrystallizable polymers is still restricted and remains a challenge. We report herein a 2D soft template-directed fabrication for ultrathin polyacrylamide nanosheets with a thickness as low as 3.5 nm and large dimensions (>20 μm). More importantly, the polymer chains in the nanosheets produced are well aligned with a clear interchain spacing. The formation of polymer nanosheets with ordered chain alignment was performed in a special solution containing a periodic sandwich structure of lamellar bilayer membranes and water layers that are hundreds of nanometers thick. It functions as a 2D orientation template to align the monomers in an orderly manner along the in-plane direction of the bimolecular membrane via hydrogen bonding.

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.

2D-IR spectroscopy of hydrogen-bond-mediated vibrational excitation transfer

Inter-molecular vibrational energy transfer in the hydrogen-bonded complexes of methyl acetate and 4-cyanophenol is studied by dual-frequency 2D-IR spectroscopy.

Excited states in large molecular systems through polarizable embedding

Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.

How Sensitive is the Amide I Vibration of the Polypeptide Backbone to Electric Fields?

Site-selective isotopic labeling of amide carbonyls offers a nonperturbative means to introduce a localized infrared probe into proteins. Although this strategy has been widely used to investigate various biological questions, the dependence of the underlying amide I vibrational frequency on electric fields (or Stark tuning rate) has not been fully determined, which prevents it from being used in a quantitative manner in certain applications. Herein, through the use of experiments and molecular dynamics simulations, the Stark tuning rate of the amide I vibration of an isotopically labeled backbone carbonyl in a transmembrane α-helix is determined to be approximately 1.4 cm(-1) /(MV/cm). This result provides a quantitative basis for using this vibrational model to assess local electric fields in proteins, among other applications. For instance, by using this value, we are able to show that the backbone region of a dipeptide has a surprisingly low dielectric constant.

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.

Modeling Vibrational Spectra of Ester Carbonyl Stretch in Water and DMSO Based on Molecular Dynamics Simulation

On the basis of molecular dynamics simulation, we model the ester carbonyl stretch FTIR signals of methyl acetate in D2O and DMSO. An ab initio map is constructed at the B3LYP/6-311++G** level to relate the carbonyl stretch frequency to the external electric field. Using this map, fluctuating Hamiltonian of the carbonyl stretch is constructed from the MD simulation trajectory. The IR spectra calculated based on this Hamiltonian are found to be in good agreement with the experiment. For methyl acetate in D2O, hydrogen bonding on alkoxy oxygen causes a blue shift of frequency, while that on carbonyl oxygen causes a red shift. Two peaks observed in FTIR signals originate from the balance of these two effects. Furthermore, in both D2O and DMSO solutions, correlations are found between the instantaneous electric field on C═O and the frequencies. Broader line width of the signal in D2O suggests a more inhomogeneous electric field distribution due to the complicated hydrogen-bonding environment.

Live-cell vibrational imaging of choline metabolites by stimulated Raman scattering coupled with isotope-based metabolic labeling

Choline is a small molecule that occupies a key position in the biochemistry of all living organisms. Recent studies have strongly implicated choline metabolites in cancer, atherosclerosis and nervous system development. To detect choline and its metabolites, existing physical methods such as magnetic resonance spectroscopy and positron emission tomography are often limited by the poor spatial resolution and substantial radiation dose. Fluorescence imaging, although with submicrometer resolution, requires introduction of bulky fluorophores and thus is difficult in labeling the small choline molecule. By combining the emerging bond-selective stimulated Raman scattering microscopy with metabolic incorporation of deuterated choline, herein we have achieved high resolution imaging of choline-containing metabolites in living mammalian cell lines, primary hippocampal neurons and the multicellular organism C. elegans. Different subcellular distributions of choline metabolites are observed between cancer cells and non-cancer cells, which may reveal a functional difference in the choline metabolism and lipid-mediated signaling events. In neurons, choline incorporation is visualized within both soma and neurites, where choline metabolites are more evenly distributed compared to proteins. Furthermore, choline localization is also observed in the pharynx region of C. elegans larvae, consistent with its organogenesis mechanism. These applications demonstrate the potential of isotope-based stimulated Raman scattering microscopy for future choline-related disease detection and development monitoring in vivo.

Measuring electric fields and noncovalent interactions using the vibrational stark effect

Over the past decade, we have developed a spectroscopic approach to measure electric fields inside matter with high spatial (<1 Å) and field (<1 MV/cm) resolution. The approach hinges on exploiting a physical phenomenon known as the vibrational Stark effect (VSE), which ultimately provides a direct mapping between observed vibrational frequencies and electric fields. Therefore, the frequency of a vibrational probe encodes information about the local electric field in the vicinity around the probe. The VSE method has enabled us to understand in great detail the underlying physical nature of several important biomolecular phenomena, such as drug-receptor selectivity in tyrosine kinases, catalysis by the enzyme ketosteroid isomerase, and unidirectional electron transfer in the photosynthetic reaction center. Beyond these specific examples, the VSE has provided a conceptual foundation for how to model intermolecular (noncovalent) interactions in a quantitative, consistent, and general manner. The starting point for research in this area is to choose (or design) a vibrational probe to interrogate the particular system of interest. Vibrational probes are sometimes intrinsic to the system in question, but we have also devised ways to build them into the system (extrinsic probes), often with minimal perturbation. With modern instruments, vibrational frequencies can increasingly be recorded with very high spatial, temporal, and frequency resolution, affording electric field maps correspondingly resolved in space, time, and field magnitude. In this Account, we set out to explain the VSE in broad strokes to make its relevance accessible to chemists of all specialties. Our intention is not to provide an encyclopedic review of published work but rather to motivate the underlying framework of the methodology and to describe how we make and interpret the measurements. Using certain vibrational probes, benchmarked against computer models, it is possible to use the VSE to measure absolute electric fields in arbitrary environments. The VSE approach provides an organizing framework for thinking generally about intermolecular interactions in a quantitative way and may serve as a useful conceptual tool for molecular design.

Identification of Subvisible Particles in Biopharmaceutical Formulations Using Raman Spectroscopy Provides Insight into Polysorbate 20 Degradation Pathway

Purpose

To study composition and heterogeneity of insoluble subvisible particles in Mab formulations resulting from degradation of polysorbate 20 and to develop a better understanding of the mechanisms of polysorbate degradation leading to particle formation.

Methods

In this study, we exploit the potential of Raman microscopy for chemical identification of particles in monoclonal antibody formulations. Through a combination of experiments and density functional theory (DFT) calculations, we identified unique spectral marker bands for insoluble degradation products of polysorbate 20. We first applied our methodology to identify particles in model systems containing complex mixtures of fatty acids and then to subvisible particles in antibody formulations stored at 5°C for several years.

Results

Most of the subvisible particles identified were comprised of mixtures of fatty acids with no observable signal from fatty acid esters consistent with hydrolysis being the predominant degradation mechanism leading to particulate formation under these storage conditions.

Conclusions

Our methodology is generally applicable for identification of particles in antibody formulations and, in particular, has the potential to give detailed information about particle heterogeneity and insight into mechanistic aspects of particle formation.

Electronic supplementary material

The online version of this article (doi:10.1007/s11095-015-1670-x) contains supplementary material, which is available to authorized users.

Extreme electric fields power catalysis in the active site of ketosteroid isomerase

Enzymes use protein architecture to impose specific electrostatic fields onto their bound substrates, but the magnitude and catalytic effect of these electric fields have proven difficult to quantify with standard experimental approaches. Using vibrational Stark effect spectroscopy, we found that the active site of the enzyme ketosteroid isomerase (KSI) exerts an extremely large electric field onto the C=O chemical bond that undergoes a charge rearrangement in KSI's rate-determining step. Moreover, we found that the magnitude of the electric field exerted by the active site strongly correlates with the enzyme's catalytic rate enhancement, enabling us to quantify the fraction of the catalytic effect that is electrostatic in origin. The measurements described here may help explain the role of electrostatics in many other enzymes and biomolecular systems.

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.

Experimental Determination of the Electrostatic Nature of Carbonyl Hydrogen-Bonding Interactions Using IR-NMR Correlations

Hydrogen-bonding plays a fundamental role in the structure, function, and dynamics of various chemical and biological systems. Understanding the physical nature of interactions and the role of electrostatics in hydrogen-bonding has been the focus of several theoretical and computational research. We present an experimental approach involving IR-(13)C NMR correlations to determine the electrostatic nature of carbonyl hydrogen-bonding interactions. This report provides a direct experimental evidence of the classical nature of hydrogen-bonding interaction in carbonyls, independent of any theoretical approximation. These results have important implications in chemistry and biology and can be applied to probe the reaction mechanisms involving carbonyl activation/stabilization by hydrogen bonds using spectroscopic techniques.

Ester carbonyl vibration as a sensitive probe of protein local electric field

The ability to quantify the local electrostatic environment of proteins and protein/peptide assemblies is key to gaining a microscopic understanding of many biological interactions and processes. Herein, we show that the ester carbonyl stretching vibration of two non-natural amino acids, L-aspartic acid 4-methyl ester and L-glutamic acid 5-methyl ester, is a convenient and sensitive probe in this regard, since its frequency correlates linearly with the local electrostatic field for both hydrogen-bonding and non-hydrogen-bonding environments. We expect that the resultant frequency-electric-field map will find use in various applications. Furthermore, we show that, when situated in a non-hydrogen-bonding environment, this probe can also be used to measure the local dielectric constant (ε). For example, its application to amyloid fibrils formed by Aβ(16-22) revealed that the interior of such β-sheet assemblies has an ε value of approximately 5.6.

Excited state electronic structures of 5,10-methenyltetrahydrofolate and 5,10-methylenetetrahydrofolate determined by Stark spectroscopy

Folates are ubiquitous cofactors that participate in a wide variety of critical biological processes. 5,10-Methenyltetrahydrofolate and its photodegradation product 5,10-methylenetetrahydrofolate are both associated with the light-driven DNA repair protein DNA photolyase and its homologues (e.g., cryptochromes). The excited state electronic properties of these folate molecules have been studied here using Stark spectroscopy and complementary quantum calculations. The tetrahydrofolates have relatively large difference dipole moments (ca. 6-8 Debye) and difference polarizabilities (ca. 100 Å(3)). This extensive excited state charge redistribution appears to be due largely to the pendant p-aminobenzoic acid group, which helps shuttle charge over the entirety of the molecule. Simple calculations based on the experimental difference dipole moments suggest that tetrahydrofolates should have large two photon cross sections sufficient to enable two photon microscopy to selectively detect and follow folate-containing proteins both in vitro and in vivo.

Effect of strong electric field on the conformational integrity of insulin

A series of molecular dynamics (MD) simulations up to 1 μs for bovine insulin monomer in different external electric fields were carried out to study the effect of external electric field on conformational integrity of insulin. Our results show that the secondary structure of insulin is kept intact under the external electric field strength below 0.15 V/nm, but disruption of secondary structure is observed at 0.25 V/nm or higher electric field strength. Although the starting time of secondary structure disruption of insulin is not clearly correlated with the strength of the external electric field ranging between 0.15 and 0.60 V/nm, long time MD simulations demonstrate that the cumulative effect of exposure time under the electric field is a major cause for the damage of insulin's secondary structure. In addition, the strength of the external electric field has a significant impact on the lifetime of hydrogen bonds when it is higher than 0.60 V/nm. The fast evolution of some hydrogen bonds of bovine insulin in the presence of the 1.0 V/nm electric field shows that different microwaves could either speed up protein folding or destroy the secondary structure of globular proteins deponding on the intensity of the external electric field.

Modeling organochlorine compounds and the σ-hole effect using a polarizable multipole force field

The charge distribution of halogen atoms on organochlorine compounds can be highly anisotropic and even display a so-called σ-hole, which leads to strong halogen bonds with electron donors. In this paper, we have systematically investigated a series of chloromethanes with one to four chloro substituents using a polarizable multipole-based molecular mechanics model. The atomic multipoles accurately reproduced the ab initio electrostatic potential around chloromethanes, including CCl₄, which has a prominent σ-hole on the Cl atom. The van der Waals parameters for Cl were fitted to the experimental density and heat of vaporization. The calculated hydration free energy, solvent reaction fields, and interaction energies of several homo- and heterodimer of chloromethanes are in good agreement with experimental and ab initio data. This study suggests that sophisticated electrostatic models, such as polarizable atomic multipoles, are needed for accurate description of electrostatics in organochlorine compounds and halogen bonds, although further improvement is necessary for better transferability.

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.

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.

Calculations of the electric fields in liquid solutions

The electric field created by a condensed-phase environment is a powerful and convenient descriptor for intermolecular interactions. Not only does it provide a unifying language to compare many different types of interactions, but it also possesses clear connections to experimental observables, such as vibrational Stark effects. We calculate here the electric fields experienced by a vibrational chromophore (the carbonyl group of acetophenone) in an array of solvents of diverse polarities using molecular dynamics simulations with the AMOEBA polarizable force field. The mean and variance of the calculated electric fields correlate well with solvent-induced frequency shifts and band broadening, suggesting Stark effects as the underlying mechanism of these key solution-phase spectral effects. Compared to fixed-charge and continuum models, AMOEBA was the only model examined that could describe nonpolar, polar, and hydrogen bonding environments in a consistent fashion. Nevertheless, we found that fixed-charge force fields and continuum models were able to replicate some results of the polarizable simulations accurately, allowing us to clearly identify which properties and situations require explicit polarization and/or atomistic representations to be modeled properly, and to identify for which properties and situations simpler models are sufficient. We also discuss the ramifications of these results for modeling electrostatics in complex environments, such as proteins.

Solvent-induced infrared frequency shifts in aromatic nitriles are quantitatively described by the vibrational Stark effect

The physical properties of solvents strongly affect the spectra of dissolved solutes, and this phenomenon can be exploited to gain insight into the solvent-solute interaction. The large solvatochromic shifts observed for many dye molecules in polar solvents are due to variations in the solvent reaction field, and these shifts are widely used to estimate the change in the dye's dipole moment upon photoexcitation, which is typically on the order of ∼1-10 D. In contrast, the change in dipole moment for vibrational transitions is approximately 2 orders of magnitude smaller. Nonetheless, vibrational chromophores display significant solvatochromism, and the relative contributions of specific chemical interactions and electrostatic interactions are debated, complicating the interpretation of vibrational frequency shifts in complex systems such as proteins. Here we present a series of substituted benzonitriles that display widely varying degrees of vibrational solvatochromism. In most cases, this variation can be quantitatively described by the experimentally determined Stark tuning rate, coupled with a simple Onsager-like model of solvation, reinforcing the view that vibrational frequency shifts are largely caused by electrostatic interactions. In addition, we discuss specific cases where continuum solvation models fail to predict solvatochromic shifts, revealing the necessity for more advanced theoretical models that capture local aspects of solute-solvent interactions.