Katie R. Mitchell-Koch

On the reorientation and hydrogen-bond dynamics of alcohols.

The mechanism of the OH bond reorientation in liquid methanol and ethanol is examined. It is found that the extended jump model, recently developed for water, describes the OH reorientation in these liquids. The slower reorientational dynamics in these alcohols compared to water can be explained by two key factors. The alkyl groups on the alcohol molecules exclude potential partners for hydrogen bonding exchanges, an effect that grows with the size of the alkyl chain. This increases the importance of the reorientation of intact hydrogen bonds, which also slows with increasing size of the alcohol and becomes the dominant reorientation pathway.

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.

Demystifying fluorine chemical shifts: electronic structure calculations address origins of seemingly anomalous 19F-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 fluorines complex electronic effects, are shown to have relevance for understanding and predicting fluorine NMR spectra.

Molecular Modeling of Cetylpyridinium Bromide, a Cationic Surfactant, in Solutions and Micelle.

Cationic surfactants are widely used in biological and industrial processes. Notably, surfactants with pyridinium salts, such as cetylpyridinium bromide (CPB), have diverse applications. The cetylpyridium cation has a quaternary nitrogen in the aromatic heterocyclic ring of the headgroup and 16 carbons in the hydrocarbon tail. At present and in the past, it has been widely used in germicides. Recently, several interesting applications of CPB have been explored, including its use in protein folding, polymerization, enzyme studies, and gene delivery as well as in pharmaceuticals as a drug delivery tool. A molecular-level understanding of CPB and its micelle in solution can enhance its development in such applications. Herein, we have proposed the first united-atom force field for CPB that yields stable micellar aggregates in molecular dynamics (MD) simulations. The force field is validated through classical MD simulations of the CPB monomer in pure water and 1-octanol as well as in an aqueous CPB micelle. We have performed principal component analysis (PCA) and calculated the translational and rotational diffusion coefficients, spatial distribution of solvent, counterion distribution, and rotational correlation time of CPB molecule in solutions and in micelle, comparing these data to previous experimental and theoretical results for a strong validation of the force field. PCA confirms that the pyridinium ring remains planar, whereas the movement of the hydrophobic tail region leads to conformational changes during the simulations. The collective modes of the pyridinum ring were identical for CPB molecule in solution and micelle, but conformational dynamics of the CPB tail were restricted in the micelle relative to motions in water and 1-octanol. Using this force field, a spherical CPB micelle was shown to be stable throughout the course of simulation, and its solvation and structural properties are characterized.

New insights into the dynamics of zwitterionic micelles and their hydration waters by gigahertz-to-terahertz dielectric spectroscopy

Gigahertz-to-terahertz spectroscopy of macromolecules in aqueous environments provides an important approach for identifying their global and transient molecular structures, as well as directly assessing hydrogen-bonding. We report dielectric properties of zwitterionic dodecylphosphocholine (DPC) micelles in aqueous solutions over a wide frequency range, from 50 MHz to 1.12 THz. The dielectric relaxation spectra reveal different polarization mechanisms at the molecular level, reflecting the complexity of DPC micelle-water interactions. We have made a deconvolution of the spectra into different components and combined them with the effective-medium approximation to separate delicate processes of micelles in water. Our measurements demonstrate reorientational motion of the DPC surfactant head groups within the micelles, and two levels of hydration water shells, including tightly and loosely bound hydration water layers. From the dielectric strength of bulk water in DPC solutions, we found that the number of waters in hydration shells is approximately constant at 950 ± 45 water molecules per micelle in DPC concentrations up to 400 mM, and it decreases after that. At terahertz frequencies, employing the effective-medium approximation, we estimate that each DPC micelle is surrounded by a tightly bound layer of 310 ± 45 water molecules that behave as if they are an integral part of the micelle. Combined with molecular dynamics simulations, we determine that tightly bound waters are directly hydrogen-bonded to oxygens of DPC, while loosely bound waters reside within 4 Å of micellar atoms. The dielectric response of DPC micelles at terahertz frequencies yields, for the first time, experimental information regarding the largest scale, lowest frequency collective motions in micelles. DPC micelles are a relatively simple biologically relevant system, and this work paves the way for more insight into future studies of hydration and dynamics of biomolecular systems with gigahertz-to-terahertz spectroscopy.

To keep or not to keep? the question of crystallographic waters for enzyme simulations in organic solvent

Abstract 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.

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.

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.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme–substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

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.