Cu-Catalyzed Three-Component Alkene Carboamination: Mechanistic Insights and Rational Design to Overcome Limitations.

Herein, we report mechanistic investigations into the Cu-catalyzed three-component carboamination of alkenes with α-halo carbonyls and aryl amines via an oxocarbenium intermediate. Monitoring the reaction reveals the formation of transient atom transfer radical addition (ATRA) intermediates with both electron-neutral and deficient vinyl arenes as well as unactivated alkenes. Based on our experimental studies and density functional theory calculations, the oxocarbenium is generated through atom transfer and subsequent intramolecular substitution. Further, mechanistic factors that dictate the regioselectivity of the nucleophilic attack onto the oxocarbenium to afford the γ-amino ester, γ-iminolactone, or γ-lactone are discussed. A strategy to overcome scope limitation with respect to unactivated alkenes is developed using the mechanistic insights gained herein. Finally, we demonstrate that under modified conditions, our Cu catalyst enables the ATRA reaction between a variety of alkyl halides and vinyl arenes/α-olefins, and we present a one-pot, two-step carbofunctionalization with an array of nucleophiles through ATRA/SN2.

Structural Insight on Supramolecular Polyion Salts: Inositol Hexaphosphate Enclosed in Cationic Macrocyclic Clusters.

Supramolecular macrocyclic forces have been used to trap phytate, myo-inositol-1,2,3,4,5,6-hexakisphosphate, a key bioanion with multiple roles in metabolic processes. Due to the complex chemistry of six multivalent phosphates surrounding the small, cyclic inositol framework, crystallographic information of simple phytate salts has been elusive. This report represents a combined crystallographic, theoretical, and solution binding investigation of a supramolecular macrocyclic complex of phytate. Together, the results provide significant insight to phytate's intramolecular and intermolecular interactions at the microenvironment level. The macrocycle-phytate aggregates consist of phytate anionic pairs, each partly sandwiched by two 24-membered, amide/amine-based cationic macrocycles. The phytate ion pairs hold the tetrameric macrocyclic array together by six strong intermolecular hydrogen bonds. Both phytates crystallize in 1a5e phosphate conformations (one axial (P2) and five equatorial phosphates). Solution NMR binding studies in 1 : 1 DMSO-d6 : D2 O indicate 2 : 1 macrocycle:phytate associations, suggesting that the sandwich-like nature of the complex holds together in solution. DFT studies indicate the likely occurrence of dynamic intramolecular interchange of phosphate protons, as well as important roles for the axial (P2) phosphate in both intramolecular and intermolecular hydrogen bonding interactions.

Newly identified C-H⋯O hydrogen bond in histidine.

New Cδ-H⋯O histidine hydrogen bonding interactions in various proteins are identified by neutron diffraction and computationally characterized. Neutron diffraction data shows several H-bond motifs with the Cδ-H moiety in histidine side chains, including interactions in ß-sheets and with coordinated waters, mostly with histidinium and τ-tautomers. In yellow protein, an active site histidine H-bonds via Cδ-H to a main chain carbonyl while the Cε-H bond coordinates a water molecule. Although the H-bonding ability of Cε-H bonds in histidine have been previously identified, analysis of neutron diffraction structures reveals Cε-H H-bonds in notable active site interactions: for the proximal histidine in myoglobin; a zinc-bound histidine in human carbonic anhydrase II; within the Ser-Asp-His catalytic triad of the trypsin active site; and a histidine in the proton shuttle mechanism of RNase A, in addition to more general roles of coordinating water and forming H-bonds with carbonyl groups in ß-sheets within a number of proteins. Properties of these H-bonds were computationally investigated using 5-methylimidazole and 5-methylimidazolium as models for histidine and histidinium. The π- and τ-tautomeric states of 5-methylimidazole were investigated, as both histidine tautomers are observed in the crystal structures. The newly characterized Cε-H⋯O and Cδ-H⋯O model complexes with water and acetone meet the overwhelming majority of IUPAC H-bonding criteria. 5-Methylimidazolium forms complexes that are nearly twice as strong as the respective neutral τ-5-methylimidazole and π-5-methylimidazole complexes. While the τ- and π-tautomers form Cε-H⋯O complexes of similar strength, the τ-Cδ-H⋯O interaction is approximately twice as strong as the π-Cδ-H⋯O interaction. Calculated charges on C-H (and N-H) hydrogens not participating in the H-bond are only slightly perturbed upon complex formation, implying that formation of one H-bond does not diminish the molecule's capacity for further H-bond formation at other sites in the imidazole ring. Overall, findings indicate that the Cδ-H⋯O interaction may be important for ß-sheet stability, conformation, interactions with solvent, and mechanisms in the active site. Recognition of C-H bond polarity and hydrogen bonding ability in histidine may improve molecular modeling and provide further insight into the diverse roles of histidine in protein structure-function-dynamics.

Tuning Cationic Micelle Properties with an Antioxidant Additive: A Molecular Perspective.

In this work, we characterize the micellization and morphology transition induced in aqueous cetyltrimethylammonium bromide (CTAB) solution by the addition of the antioxidant propyl gallate (PG) using tensiometry, rheology, and small-angle neutron scattering (SANS) techniques combined with the molecular dynamics (MD) simulation approach. The adsorption of CTAB at the air-water interface in the presence of varying [PG] revealed a progressive decrease in the critical micelle concentration (CMC), while the changes in different interfacial parameters indicated enhancement of the hydrophobicity induced by PG in the CTAB micellar system. The dynamic rheology behavior indicated an increase in the flow viscosity (η) as a function of [PG]. Moreover, the rheological components (storage modulus, G', and loss modulus, G″) depicted the viscoelastic features. SANS measurements depicted the existence of ellipsoidal micelles with varying sizes and aggregation number (Nagg) as a function of [PG] and temperature. Computational simulation performed using density functional theory (DFT) calculations and molecular dynamics (MD) provided an insight into the atomic composition of the examined system. The molecular electrostatic potential (MEP) analysis depicted a close proximity of CTAB, i.e., emphasized favorable interactions between the quaternary nitrogen of CTAB and the hydroxyl group of the PG monomer, further validated by the two-dimensional nuclear Overhauser enhancement spectroscopy (2D-NOESY), which showed the penetration of PG inside the CTAB micelles. In addition, various dynamic properties, viz., the radial distribution function (RDF), the radius of gyration (Rg), and solvent-accessible surface area (SASA), showed a significant microstructural evolution of the ellipsoidal micelles in the examined CTAB-PG system, where the changes in the micellar morphology with a more elongated hydrophobic chain and the increased Rg and SASA values indicated the notable intercalation of PG in the CTAB micelles.

Contrasting effect of 1-butanol and 1,4-butanediol on the triggered micellar self-assemblies of C16-type cationic surfactants.

The self-assembly in aqueous solutions of three quaternary salt-based C16-type cationic surfactants with different polar head groups and identical carbon alkyl chain viz., cetylpyridinium bromide (CPB), cetyltrimethylammonium tosylate (CTAT), and cetyltriphenylphosphonium bromide (CTPPB) in the presence of 1-butanol (BuOH) and 1,4-butanediol (BTD) was investigated using tensiometry, 2D-nuclear Overhauser enhancement spectroscopy (2D-NOESY) and small angle neutron scattering (SANS) techniques. The adsorption parameters and micellar characteristics evaluated at 303.15 K distinctly showed that BuOH promotes the mixed micelle formation while BTD interfered with the micellization phenomenon. The SANS data fitted using an ellipsoid (as derived by Hayter and Penfold using the Ornstein-Zernike equation and the mean spherical approximation) and wormlike micellar models offered an insight into the micelle size/shape and aggregation number (Nagg) in the examined systems. The evaluated descriptors presented a clear indication of the morphology transition in cationic micelles as induced by the addition of the two alcohols. We also offer an investigation into the acceptable molecular interactions governing the differences in micelle morphologies, using the non-invasive 2D-NOESY technique and molecular modeling. The experimental observations elucidated from computational simulation add novelty to this work. Giving an account to the structural complexity in the three cationic surfactants, the molecular dynamics (MD) simulation was performed for CPB micelles in an aqueous solution of alcohols that highlighted the micelle solvation and structural transition, which is further complemented in terms of critical packing parameter (PP) for the examined systems.

Dynamic Preference for NADP/H Cofactor Binding/Release in E. coli YqhD Oxidoreductase.

YqhD, an E. coli alcohol/aldehyde oxidoreductase, is an enzyme able to produce valuable bio-renewable fuels and fine chemicals from a broad range of starting materials. Herein, we report the first computational solution-phase structure-dynamics analysis of YqhD, shedding light on the effect of oxidized and reduced NADP/H cofactor binding on the conformational dynamics of the biocatalyst using molecular dynamics (MD) simulations. The cofactor oxidation states mainly influence the interdomain cleft region conformations of the YqhD monomers, involved in intricate cofactor binding and release. The ensemble of NADPH-bound monomers has a narrower average interdomain space resulting in more hydrogen bonds and rigid cofactor binding. NADP-bound YqhD fluctuates between open and closed conformations, while it was observed that NADPH-bound YqhD had slower opening/closing dynamics of the cofactor-binding cleft. In the light of enzyme kinetics and structural data, simulation findings have led us to postulate that the frequently sampled open conformation of the cofactor binding cleft with NADP leads to the more facile release of NADP while increased closed conformation sampling during NADPH binding enhances cofactor binding affinity and the aldehyde reductase activity of the enzyme.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics.

The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

Entropy connects water structure and dynamics in protein hydration layer.

The enzyme Candida Antarctica lipase B (CALB) serves here as a model for understanding connections among hydration layer dynamics, solvation shell structure, and protein surface structure. The structure and dynamics of water molecules in the hydration layer were characterized for regions of the CALB surface, divided around each α-helix, ß-sheet, and loop structure. Heterogeneous hydration dynamics were observed around the surface of the enzyme, in line with spectroscopic observations of other proteins. Regional differences in the structure of the biomolecular hydration layer were found to be concomitant with variations in dynamics. In particular, it was seen that regions of higher density exhibit faster water dynamics. This is analogous to the behavior of bulk water, where dynamics (diffusion coefficients) are connected to water structure (density and tetrahedrality) by excess (or pair) entropy, detailed in the Rosenfeld scaling relationship. Additionally, effects of protein surface topology and hydrophobicity on water structure and dynamics were evaluated using multiregression analysis, showing that topology has a somewhat larger effect on hydration layer structure-dynamics. Concave and hydrophobic protein surfaces favor a less dense and more tetrahedral solvation layer, akin to a more ice-like structure, with slower dynamics. Results show that pairwise entropies of local hydration layers, calculated from regional radial distribution functions, scale logarithmically with local hydration dynamics. Thus, the Rosenfeld relationship describes the heterogeneous structure-dynamics of the hydration layer around the enzyme CALB. These findings raise the question of whether this may be a general principle for understanding the structure-dynamics of biomolecular solvation.

Evaluating electronic structure methods for accurate calculation of 19 F chemical shifts in fluorinated amino acids.

The ability of electronic structure methods (11 density functionals, HF, and MP2 calculations; two basis sets and two solvation models) to accurately calculate the 19 F chemical shifts of 31 structures of fluorinated amino acids and analogues with known experimental 19 F NMR spectra has been evaluated. For this task, BHandHLYP, ωB97X, and Hartree-Fock with scaling factors (provided within) are most accurate. Additionally, the accuracy of methods to calculate relative changes in fluorine shielding across 23 sets of structural variants, such as zwitterionic amino acids versus side chains only, was also determined. This latter criterion may be a better indicator of reliable methods for the ultimate goal of assigning and interpreting chemical shifts of fluorinated amino acids in proteins. It was found that MP2 and M062X calculations most accurately assess changes in shielding among analogues. These results serve as a guide for computational developments to calculate 19 F chemical shifts in biomolecular environments. © 2017 Wiley Periodicals, Inc.

Tautomeric stabilities of 4-fluorohistidine shed new light on mechanistic experiments with labeled ribonuclease A.

Ribonuclease A is the oldest model for studying enzymatic mechanisms, yet questions remain about proton transfer within the active site. Seminal work by Jackson et al. (Science, 1994) labeled Ribonuclease A with 4-fluorohistidine, concluding that active-site histidines act as general acids and bases. Calculations of 4-fluorohistidine indicate that the π-tautomer is predominant in all simulated environments (by ~17 kJ/mol), strongly suggesting that fluoro-labeled ribonuclease A functions with His119 in π-tautomer. The tautomeric form of His119 during proton transfer and tautomerism as a putative mechanistic step in wild-type RNase A remain open questions and should be considered in future mechanistic studies.

Demystifying fluorine chemical shifts: electronic structure calculations address origins of seemingly anomalous (19)F-NMR spectra of fluorohistidine isomers and analogues.

Fluorine NMR spectroscopy is a powerful tool for studying biomolecular structure, dynamics, and ligand binding, yet the origins of (19)F chemical shifts are not well understood. Herein, we use electronic structure calculations to describe the changes in (19)F chemical shifts of 2F- and 4F-histidine/(5-methyl)-imidazole upon acid titration. While the protonation of the 2F species results in a deshielded chemical shift, protonation of the 4F isomer results in an opposite, shielded chemical shift. The deshielding of 2F-histidine/(5-methyl)-imidazole upon protonation can be rationalized by concomitant decreases in charge density on fluorine and a reduced dipole moment. These correlations do not hold for 4F-histidine/(5-methyl)-imidazole, however. Molecular orbital calculations reveal that for the 4F species, there are no lone pair electrons on the fluorine until protonation. Analysis of a series of 4F-imidazole analogues, all with delocalized fluorine electron density, indicates that the deshielding of (19)F chemical shifts through substituent effects correlates with increased C-F bond polarity. In summary, the delocalization of fluorine electrons in the neutral 4F species, with gain of a lone pair upon protonation may help explain the difficulty in developing a predictive framework for fluorine chemical shifts. Ideas debated by chemists over 40 years ago, regarding fluorine's complex electronic effects, are shown to have relevance for understanding and predicting fluorine NMR spectra.

Probing Adaptation of Hydration and Protein Dynamics to Temperature.

Protein dynamics is strongly influenced by the surrounding environment and physiological conditions. Here we employ broadband megahertz-to-terahertz spectroscopy to explore the dynamics of water and myoglobin protein on an extended time scale from femto- to nanosecond. The dielectric spectra reveal several relaxations corresponding to the orientational polarization mechanism, including the dynamics of loosely bound, tightly bound, and bulk water, as well as collective vibrational modes of protein in an aqueous environment. The dynamics of loosely bound and bulk water follow non-Arrhenius behavior; however, the dynamics of water molecules in the tightly bound layer obeys the Arrhenius-type relation. Combining molecular simulations and effective-medium approximation, we have determined the number of water molecules in the tightly bound hydration layer and studied the dynamics of protein as a function of temperature. The results provide the important impact of water on the biochemical functions of proteins.

Self-Assembly and Micellar Transition in CTAB Solutions Triggered by 1-Octanol.

This study exploits higher-order micellar transition ranging from ellipsoidal to rodlike to wormlike induced by 1-octanol (C8OH) in an aqueous solution of cetyltrimethylammonium bromide (CTAB), characterizing phase behavior, rheology, and small-angle neutron scattering (SANS). The phase diagram for the ternary system CTAB-C8OH-water was constructed, which depicted the varied solution behavior. Such performance was further inferred from the rheology study (oscillatory-shear frequency sweep (ω) and viscosity (η)) that displayed an interesting solution behavior of CTAB solutions as a function of C8OH. It was observed that at low C8OH concentrations, the solutions appeared viscous/viscoelastic fluids that changed to an elastic gel with an infinite relaxation time at higher concentrations of C8OH, thereby confirming the existence of distinct micelle morphologies. Small-angle neutron scattering (SANS) provided various micellar parameters such as aggregation numbers (Nagg) and micellar size/shape. The experimental results were further validated with a computational simulation approach. The molecular dynamic (MD) study offered an insight into the molecular interactions and aggregation behavior through different analyses, including radial distribution function (RDF), radius of gyration (Rg), and solvent-accessible surface area (SASA).

The study and application of biomolecules in deep eutectic solvents.

Biomolecules have been thoroughly investigated in a multitude of solvents historically in order to accentuate or modulate their superlative properties in an array of applications. Ionic liquids have been extensively explored over the last two decades as potential replacements for traditional organic solvents, however, they are sometimes associated with a number of limitations primarily related to cost, convenience, accessibility, and/or sustainability. One potential solvent which is gaining considerable traction in recent years is the so-called deep eutectic solvent which holds a number of striking advantages, including biodegradability, inherently low toxicity, and a facile, low-cost, and solventless preparation from widely available natural feedstocks. In this review, we highlight recent progress and insights into biomolecular behavior within deep eutectic solvent-containing systems, including discussions of their demonstrated utility and prospects for the biostabilization of proteins and nucleic acids, free enzyme and whole-cell biocatalysis, various extraction processes (e.g., aqueous biphasic systems, nanosupported separations), drug solubilization, lignocellulose biomass treatment, and targeted therapeutic drug delivery. All indications point to the likelihood that these emerging solvents have the capacity to satisfy the requirements of environmental responsibility while unlocking biomolecular proficiency in established biomedical and biotechnological pursuits as well as a number of academic and industrial ventures not yet explored.

Infrared spectra of a model phenol-amine proton transfer complex in nanoconfined CH3Cl.

The vibrational spectra of a model phenol-amine proton transfer complex dissolved in CH3Cl solvent confined in a 12 A radius spherical hydrophobic cavity were calculated using mixed quantum-classical molecular dynamics simulations. The reaction free energy of the proton transfer complex was varied in order to explore the contributions to the vibrational absorption band from product and reactant species. The vibrational spectra of the model proton transfer complex resulted in motionally narrowed spectral linewidths with two distinct peaks for products and reactants in cases where the system undergoes chemical exchange. It was found that the n=1 and n=2 vibrational excited states combine to form diabatic states such that the spectra have contributions from both n=0 --> n=1 and n=0 --> n=2 transitions. A strong relationship between the instantaneous vibrational frequency and a collective solvent coordinate was found that assists in understanding the origin of the spectral features.

Evaluating computational predictions of the relative stabilities of polymorphic pharmaceuticals.

The ability of computational methods to describe the relative energies of polymorphic pharmaceuticals is investigated for a diverse array of compounds. The initial molecular geometries were taken from crystal structures, and energy differences between polymorphic pairs were calculated with various geometry optimization methods. Results using molecular mechanics were compared to experimental calorimetric data and periodic density functional theory (DFT) calculations. The best agreement with experimental heats of transition was shown with energies calculated from geometry optimizations using the Compass force field. Calculations that optimized atomic positions with the Compass force field gave correct energy rankings for all 11 polymorphic pairs studied, with an average deviation of 0.61 kcal/mol from experimental results. These findings suggest that computational methods are poised to predict enthalpy differences between polymorphic forms with levels of accuracy that are quite acceptable when proper approaches are employed.

How Does Solvation Layer Mobility Affect Protein Structural Dynamics?

Solvation is critical for protein structural dynamics. Spectroscopic studies have indicated relationships between protein and solvent dynamics, and rates of gas binding to heme proteins in aqueous solution were previously observed to depend inversely on solution viscosity. In this work, the solvent-compatible enzyme Candida antarctica lipase B, which functions in aqueous and organic solvents, was modeled using molecular dynamics simulations. Data was obtained for the enzyme in acetonitrile, cyclohexane, n-butanol, and tert-butanol, in addition to water. Protein dynamics and solvation shell dynamics are characterized regionally: for each α-helix, ß-sheet, and loop or connector region. Correlations are seen between solvent mobility and protein flexibility. So, does local viscosity explain the relationship between protein structural dynamics and solvation layer dynamics? Halle and Davidovic presented a cogent analysis of data describing the global hydrodynamics of a protein (tumbling in solution) that fits a model in which the protein's interfacial viscosity is higher than that of bulk water's, due to retarded water dynamics in the hydration layer (measured in NMR τ2 reorientation times). Numerous experiments have shown coupling between protein and solvation layer dynamics in site-specific measurements. Our data provides spatially-resolved characterization of solvent shell dynamics, showing correlations between regional solvation layer dynamics and protein dynamics in both aqueous and organic solvents. Correlations between protein flexibility and inverse solvent viscosity (1/η) are considered across several protein regions and for a rather disparate collection of solvents. It is seen that the correlation is consistently higher when local solvent shell dynamics are considered, rather than bulk viscosity. Protein flexibility is seen to correlate best with either the local interfacial viscosity or the ratio of the mobility of an organic solvent in a regional solvation layer relative to hydration dynamics around the same region. Results provide insight into the function of aqueous proteins, while also suggesting a framework for interpreting and predicting enzyme structural dynamics in non-aqueous solvents, based on the mobility of solvents within the solvation layer. We suggest that Kramers' theory may be used in future work to model protein conformational transitions in different solvents by incorporating local viscosity effects.

High-Precision Megahertz-to-Terahertz Dielectric Spectroscopy of Protein Collective Motions and Hydration Dynamics.

The low-frequency collective vibrational modes in proteins as well as the protein-water interface have been suggested as dominant factors controlling the efficiency of biochemical reactions and biological energy transport. It is thus crucial to uncover the mystery of the hydration structure and dynamics as well as their coupling to collective motions of proteins in aqueous solutions. Here, we report dielectric properties of aqueous bovine serum albumin protein solutions as a model system using an extremely sensitive dielectric spectrometer with frequencies spanning from megahertz to terahertz. The dielectric relaxation spectra reveal several polarization mechanisms at the molecular level with different time constants and dielectric strengths, reflecting the complexity of protein-water interactions. Combining the effective-medium approximation and molecular dynamics simulations, we have determined collective vibrational modes at terahertz frequencies and the number of water molecules in the tightly bound and loosely bound hydration layers. High-precision measurements of the number of hydration water molecules indicate that the dynamical influence of proteins extends beyond the first solvation layer, to around 7 Å distance from the protein surface, with the largest slowdown arising from water molecules directly hydrogen-bonded to the protein. Our results reveal critical information of protein dynamics and protein-water interfaces, which determine biochemical functions and reactivity of proteins.

The Biophysical Probes 2-fluorohistidine and 4-fluorohistidine: Spectroscopic Signatures and Molecular Properties.

Fluorinated amino acids serve as valuable biological probes, by reporting on local protein structure and dynamics through 19F NMR chemical shifts. 2-fluorohistidine and 4-fluorohistidine, studied here with DFT methods, have even more capabilities for biophysical studies, as their altered pKa values, relative to histidine, allow for studies of the role of proton transfer and tautomeric state in enzymatic mechanisms. Considering the two tautomeric forms of histidine, it was found that 2-fluorohistidine primarily forms the common (for histidine) τ-tautomer at neutral pH, while 4-fluorohistidine exclusively forms the less common π-tautomer. This suggests the two isomers of fluorohistidine can also serve as probes of tautomeric form within biomolecules, both by monitoring NMR chemical shifts and by potential perturbation of the tautomeric equilibrium within biomolecules. Fluorine also enables assignment of tautomeric states in crystal structures. The differences in experimental pKa values between the isomers was found to arise from solvation effects, providing insight into the polarization and molecular properties of each isomer. Results also encompass 13C and 19F NMR chemical shifts, from both tautomers of 2-fluorohistidine and 4-fluorohistidine in a number of different environments. This work can serve as a guide for interpretation of spectroscopic results in biophysical studies employing 2-fluorohistidine and 4-fluorohistidine.

To Keep or Not to Keep? The Question of Crystallographic Waters for Enzyme Simulations in Organic Solvent.

The use of enzymes in non-aqueous solvents expands the use of biocatalysts to hydrophobic substrates, with the ability to tune selectivity of reactions through solvent selection. Non-aqueous enzymology also allows for fundamental studies on the role of water and other solvents in enzyme structure, dynamics, and function. Molecular dynamics simulations serve as a powerful tool in this area, providing detailed atomic information about the effect of solvents on enzyme properties. However, a common protocol for non-aqueous enzyme simulations does not exist. If you want to simulate enzymes in non-aqueous solutions, how many and which crystallographic waters do you keep? In the present work, this question is addressed by determining which crystallographic water molecules lead most quickly to an equilibrated protein structure. Five different methods of selecting and keeping crystallographic waters are used in order to discover which crystallographic waters lead the protein structure to reach an equilibrated structure more rapidly in organic solutions. It is found that buried waters contribute most to rapid equilibration in organic solvent, with slow-diffusing waters giving similar results.