Enthalpic Components Research Articles

Changing Your Martini Can Still Give You a Hangover.

The Martini 3.0 coarse-grained force field, which was parametrized to better capture transferability in top-down coarse-grained models, is analyzed to assess its accuracy in representing thermodynamic and structural properties with respect to the underlying atomistic representation of the system. These results are compared to those obtained following the principles of statistical mechanics that start from the same underlying atomistic system. To this end, the potentials of mean force for lateral association in Martini 3.0 binary lipid bilayers are decomposed into their entropic and enthalpic components and compared to those of corresponding atomistic bilayers that have been projected onto equivalent coarse-grained mappings but evolved under the fully atomistic forces. This is accomplished by applying the reversible work theorem to lateral pair correlation functions between coarse-grained lipid beads taken at a range of different temperatures. The entropy-enthalpy decompositions provide a metric by which the underlying statistical mechanical properties of Martini can be investigated. Overall, Martini 3.0 is found to fail to properly partition entropy and enthalpy for the PMFs compared to the mapped all-atom results, despite changes made to the force field from the Martini 2.0 version. This outcome points to the fact that the development of more accurate top-down coarse-grained models such as Martini will likely necessitate temperature-dependent terms in the corresponding CG force-field; although necessary, this may not be sufficient to improve Martini. In addition to the entropy-enthalpy decompositions, Martini 3.0 produces an incorrect undulation spectrum, in particular at intermediate length scales of biophysical pertinence.

Tracing the slow Arrhenius process deep in the glassy state-quantitative evaluation of the dielectric relaxation of bulk samples and thin polymer films in the temperature domain.

The slow Arrhenius process (SAP) is a dielectric mode connected to thermally activated equilibration mechanisms, allowing for a fast reduction in free energy in liquids and glasses. The SAP, however, is still poorly understood, and so far, this process has mainly been investigated at temperatures above the glass transition. By employing a combination of methods to analyze dielectric measurements under both isochronal and isothermal conditions, we were able to quantitatively reproduce the dielectric response of the SAP of different polymers and to expand the experimental regime over which this process can be observed down to lower temperatures, up to 70K below the glass transition. Employing thin films of thicknesses varying between 10 and 800nm, we further verified that the peak shape and activation energy of the SAP of poly(4-bromostyrene) are not sensitive to temperature, nor do they vary upon confinement at the nanoscale level. These observations confirm the preliminary trends reported for other polymers. We find that one single set of parameters-meaning the activation barrier and the pre-exponential factor, respectively, linked to the enthalpic and entropic components of the process-can describe the dynamics of the SAP in both the supercooled liquid and glassy states, in bulk and thin films. These results are discussed in terms of possible molecular origins of the slow Arrhenius process in polymers.

Enthalpies and entropies of hydration from Monte Carlo simulations.

The changes in free energy, enthalpy, and entropy for transfer of a solute from the gas phase into solution are the fundamental thermodynamic quantities that characterize the solvation process. Owing to the development of methods based on free-energy perturbation theory, computation of free energies of solvation has become routine in conjunction with Monte Carlo (MC) statistical mechanics and molecular dynamics (MD) simulations. Computation of the enthalpy change and by inference the entropy change is more challenging. Two methods are considered in this work corresponding to direct averaging for the solvent and solution and to computing the temperature derivative of the free energy in the van't Hoff approach. The application is for neutral organic solutes in TIP4P water using long MC simulations to improve precision. Definitive results are also provided for pure TIP4P water. While the uncertainty in computed free energies of hydration is ca. 0.05 kcal mol-1, it is ca. 0.4 kcal mol-1 for the enthalpy changes from either van't Hoff plots or the direct method with sampling for 5 billion MC configurations. Partial molar volumes of hydration are also computed by the direct method; they agree well with experimental data with an average deviation of 3 cm3 mol-1. In addition, the results permit breakdown of the errors in the free energy changes from the OPLS-AA force field into their enthalpic and entropic components. The excess hydrophobicity of organic solutes is enthalpic in origin.

The impact of chemical properties of the solid-liquid-adsorbate interfaces on the entropy-enthalpy compensation involved in adsorption.

Despite extensive studies on the thermodynamic mechanism governing molecular adsorption at the solid-water interface, a comprehensive understanding of the crucial role of interface properties in mediating the entropy-enthalpy compensation during adsorption is lacking, particularly at a quantitative level. Herein, we employed two types of surface models (hydroxyapatite and graphene) along with a series of amino acids to successfully elucidate how distinct interfacial features dictate the delicate balance between entropy and enthalpy variations. The adsorption of all amino acids on the hydroxyapatite surface is an enthalpy-dominated process, where the water-induced enthalpic component of the free energy and the surface-adsorbate electrostatic interaction term alternatively act as the driving force for adsorption in different regions of the surface. Although favorable interactions are observed between amino acids and the graphene surface, the entropy-enthalpy compensation exhibits dependence on the molecular size of the adsorbates. For small amino acids, favorable enthalpy changes predominantly determine their adsorption behavior; however, larger amino acids tend to bind more tightly with the graphene surface, which is thermodynamically dominated by the entropy variations despite the structural characteristics of amino acids. This study reveals specific entropy-enthalpy mechanisms underlying amino acid adsorption at the solid-liquid interface, providing guidance for surface design and synthesis of new biomolecules.

Influence of temperature on non-exchange sorption of aromatic amino acids by anion exchanger AV-17-8

The equilibrium in thearomatic amino acid (tryptophan, phenylalanine, tyrosine and histidine) – anion exchanger AV-17-8 (Cl-) systems at different temperatures was investigated In order to identify the influence of both increased and decreased temperatures on the sorption equilibrium the temperature range was 277-343 K was studied. Sorption isotherms were obtained at 277, 293, and 343 K. The thermodynamic description of sorption systems was carried out using an approach based on the idea that the sorption process is stoichiometric. The integral and differential energy characteristics and the compensation effect between the enthalpy and entropy components for the studied sorption systems were calculated. A change in temperature led to a change in the physicochemical properties of the studied sorbates and sorbent. In the temperature range 277-243 K, the solvent content in the sorbent phase remained unchanged. A low solvent content at a temperature of 333 K and an increase in its amount at 343 K were experimentally established. It should be considered that at different temperatures the ionization constants of functional groups and the regions of existence of ionic forms of amino acids change, but the sorption mechanism does not change, since the recharging of amino acids did not occur. It was established that the process is exothermic for tryptophan and phenylalanine. With increasing temperature, the amount of sorbed substance decreased. The absorption energy was low due to the significant compensatory effect of enthalpy and entropy factors, the differential Gibbs energy was practically independent of temperature due to the stabilizing effect of hydrophobic forces in the sorbent phase. For the poorly soluble amino acid tyrosine, when the temperature increased from 277 to 293 K, a slight increase in sorption was observed, and with a further increase, a decrease to trace amounts was revealed. In the temperature range of 277-293 K, the sorption process was characterized by the formation of a supersaturated solution in the sorbent phase and was endothermic. The sorption of hydrophilic histidine increased with increasing temperature, which correlates with the solvent content in the sorbent phase, including low content at 333 K. The process of histidine sorption is endothermic, characterized by a strong dependence of the differential Gibbs energy on temperature and the absence of a compensation effect between the enthalpic and entropic components.

Enthalpy-entropy compensation in the slow Arrhenius process.

The Meyer-Neldel compensation law, observed in a wide variety of chemical reactions and other thermally activated processes, provides a proportionality between the entropic and the enthalpic components of an energy barrier. By analyzing 31 different polymer systems, we show that such an intriguing behavior is encountered also in the slow Arrhenius process, a recently discovered microscopic relaxation mode, responsible for several equilibration mechanisms both in the liquid and the glassy state. We interpret this behavior in terms of the multiexcitation entropy model, indicating that overcoming large energy barriers can require a high number of low-energy local excitations, providing a multiphonon relaxation process.

Efficient Empirical Valence Bond Simulations with GROMACS.

We describe a protocol to perform empirical valence bond (EVB) simulations using GROMACS software. EVB is a fast and reliable method that allows one to determine the reaction free-energy profiles in complex systems, such as enzymes, by employing classical force fields to represent a chemical reaction. Therefore, running EVB simulations is basically as fast as any classical molecular dynamics simulation, and the method uses standard free-energy calculations to map the free-energy change along a given reaction path. To exemplify and validate our EVB implementation, we replicated two cases of our earlier enzyme simulations. One of these addresses the decomposition of the activation free energy into its enthalpic and entropic components, and the other is focused on calculating the overall catalytic effect of the enzyme compared to the same reaction in water. These two examples give virtually identical results to those obtained with programs that were specifically designed for EVB simulations and show that the GROMACS implementation is robust and can be used for very large systems.

Determinants for macromolecular crowding-induced thermodynamic stabilization of acid-denatured cytochrome c to molten globules

The macromolecular crowding effect transforms the acid-denatured ferricytochrome c (cyt cIII) (UA-state) to molten-globule (MGMC-state) at pH 1.85. Crowding-induced stabilization free energy (ΔΔG) and preferential hydration ΔΓW were estimated for the UA → MGMC transition. The magnitudes of ΔΔG and ΔΓW were found to be decreased as dextran 70 (D70) > dextran 40 (D40) > ficoll 70 (F70), which demonstrates that ΔΔG and ΔΓW track the molecular size and shape of the crowder towards refolding and stabilization of UA-state to MGMC-state. Analysis of effects of crowders (D40, D70, F70) on thermal and chemical-denaturations of acid-denatured cyt cIII provided several important information, (i) macromolecular crowding increased the thermodynamic stability of acid-denatured cyt cIII, (ii) concentration, size and shape of crowder control the crowding-induced thermodynamic stabilization of MGMC-state, (iii) crowding effect increased the thermal-denaturation midpoint (Tm) with a slight change in enthalpy (ΔHm), suggesting that the steric-excluded volume effect contributes to the crowding-induced increased thermal stability of the acid-denatured protein. Analysis of entropy − enthalpy plots for D40, D70, and F70 reveals that in addition to the steric-excluded volume effect, the enthalpic contribution is also added to the macromolecular crowding-induced stabilization of acid-denatured cyt cIII. The dilute-medium, compound-crowder, purely entropic-crowder and purely enthalpic-crowder curves were obtained for acid-denatured cyt cIII for D70, D40 and F70. The crossover temperature, Tx was calculated from the dilute and compound-crowder curves. The Tx values measured for D40, D70, and F70 were found to be ∼ 250.15 K, 272.15 K, and 275.15 K, respectively, which suggests that the Tx value depends on the size and shape of the crowder. Furthermore, the observation of a lower value of Tx and a minor enthalpic component for D40, D70, and F70 is likely due to the formation of weaker soft interactions of acid-denatured cyt cIII with D40, D70, and F70.

Excluded Volume and Weak Interactions in Crowded Solutions Modulate Conformations and RNA Binding of an Intrinsically Disordered Tail.

The cellular milieu is a solution crowded with a significant concentration of different components (proteins, nucleic acids, metabolites, etc.). Such a crowded environment affects protein conformations, dynamics, and interactions. Intrinsically disordered proteins and regions are particularly sensitive to these effects. Here, we investigate the impact on an intrinsically disordered tail that flanks a folded domain, the N-terminal domain, and the RNA-binding domain of the SARS-CoV-2 nucleocapsid protein. We mimic the crowded environment of the cell using polyethylene glycol (PEG) and study its impact on protein conformations using single-molecule Förster resonance energy transfer. We found that high-molecular-weight PEG induces a collapse of the disordered N-terminal tail, whereas low-molecular-weight PEG induces a chain expansion. Our data can be explained by accounting for two opposing contributions: favorable interactions between the protein and crowder molecules and screening of excluded volume interactions. We further characterized the interaction between protein and RNA in the presence of crowding agents. While for all PEG molecules tested, we observed an increase in the binding affinity, the trend is not monotonic as a function of the degree of PEG polymerization. This points to the role of nonspecific protein-PEG interactions on binding in addition to the entropic effects due to crowding. To separate the enthalpic and entropic components of the effects, we investigated the temperature dependence of the association constants in the absence and presence of crowders. Finally, we compared the effects of crowding across mutations in the disordered region and found that the threefold difference in association constants for two naturally occurring variants of the SARS-CoV-2 nucleocapsid protein is reduced to almost identical affinities in the presence of crowders. Overall, our data provide new insights into understanding and modeling the contribution of crowding effects on disordered regions, including the impact of interactions between proteins and crowders and their interplay when binding a ligand.

SophosQM: Accurate Binding Affinity Prediction in Compound Optimization.

The optimization of compounds' binding affinity for a biological target is a crucial aspect of the drug development process. Being able to accurately predict binding energies in advance of synthesizing compounds would have a massive impact on the speed of the drug discovery process. The ideal binding affinity prediction method should combine accuracy, reliability, and speed. In this paper, we present SophosQM, a quantum mechanics (QM)-based approach, which can accurately predict the binding affinities of compounds to proteins. The binding affinity predictive models generated by SophosQM are based on the fragment molecular orbital (FMO) method to estimate the enthalpic component of the binding free energy, and a macroscopic descriptor, clog P, is used as an approximation of the entropic component. The affinity prediction is performed using multilinear regression, fitting the experimental values against the FMO-computed enthalpic term and clog P. The quality of the prediction can be assessed in terms of the correlation coefficient between experimental and predicted values. In this work, the method's reliability and accuracy are exemplified by applying SophosQM to 70 compounds binding to six different targets of pharmaceutical relevance. Overall, the results show a very satisfactory performance with a global correlation coefficient in the order of 0.9. Our predictions also show a satisfactory performance compared to data based on free energy perturbation. Finally, SophosQM can also be applied in high-throughput mode by using semiempirical QM methods to evaluate large portions of chemical space, while retaining a good level of accuracy, but decreasing the computing time to just a few seconds per compound.

Miscibility of Phosphatidylcholines in Bilayers: Effect of Acyl Chain Unsaturation

The miscibility of phospholipids in a hydrated bilayer is an issue of fundamental importance for understanding the organization of biological membranes. Despite research on lipid miscibility, its molecular basis remains poorly understood. In this study, all-atom MD simulations complemented by Langmuir monolayer and DSC experiments have been performed to investigate the molecular organization and properties of lipid bilayers composed of phosphatidylcholines with saturated (palmitoyl, DPPC) and unsaturated (oleoyl, DOPC) acyl chains. The experimental results showed that the DOPC/DPPC bilayers are systems exhibiting a very limited miscibility (strongly positive values of excess free energy of mixing) at temperatures below the DPPC phase transition. The excess free energy of mixing is divided into an entropic component, related to the ordering of the acyl chains, and an enthalpic component, resulting from the mainly electrostatic interactions between the headgroups of lipids. MD simulations showed that the electrostatic interactions for lipid like-pairs are much stronger than that for mixed pairs and temperature has only a slight influence on these interactions. On the contrary, the entropic component increases strongly with increasing temperature, due to the freeing of rotation of acyl chains. Therefore, the miscibility of phospholipids with different saturations of acyl chains is an entropy-driven process.

Associated Water Dynamics Might Be a Key Factor Affecting Protein Stability in the Crowded Milieu.

Over the past 20 years, the most studied and debated aspect of macromolecular crowding is how it affects protein stability. Traditionally, it is explained by a delicate balance between the stabilizing entropic effect and the stabilizing or destabilizing enthalpic effect. However, this traditional crowding theory cannot explain experimental observations like (i) negative entropic effect and (ii) entropy-enthalpy compensation. Herein, we provide experimental evidence that associated water dynamics plays a crucial role in controlling protein stability in the crowded milieu for the first time. We have correlated the modulation of associated water dynamics with the overall stability and its individual components. We showed that rigid associated water would stabilize the protein through entropy but destabilize it through enthalpy. In contrast, flexible associated water destabilizes the protein through entropy but stabilizes through enthalpy. Consideration of entropic and enthalpic modulation through crowder-induced distortion of associated water successfully explains the negative entropic part and entropy-enthalpy compensation. Furthermore, we argued that the relationship between the associated water structure and protein stability should be better understood by individual entropic and enthalpic components instead of the overall stability. Although a huge effort is necessary to generalize the mechanism, this report provides a unique way of understanding the relationship between protein stability and associated water dynamics, which might be a generic phenomenon and should trigger much research in this area.

Unraveling the enthalpy of macromolecular crowding.

Concentrated solutions of synthetic polymers are often used to stabilize proteins, but there are few informed models that provide a detailed explanation of the effect. Classic crowding models are based on spherical, colloid-based solution theories, but such models only account for steric effects. To aid in the establishment of new, more informed, model, we acquired an extensive dataset that quantifies the effect of polyethylene glycol molecular weight and concentration on the temperature dependence of protein stability using 19F NMR. More specifically, we investigated the stability of the 7-kDa metastable N-terminal SH3 domain of Drosophila signal transduction protein drk (SH3) in polyethylene glycol solutions and quantified the modified standard-state free energy of unfolding and its enthalpic and entropic components. We used these data to build a new model which shows that interactions of solution components with both the folded and unfolded states must be considered. With this more nuanced understanding, we can better predict how polymers will affect the stability of protein-based drugs, how different biological molecules impact protein stability in cells and perhaps inform the design of even more useful stabilizing agents.

The role of loop dynamics in the prediction of ligand-protein binding enthalpy.

The enthalpic and entropic components of ligand-protein binding free energy reflect the interactions and dynamics between ligand and protein. Despite decades of study, our understanding and hence our ability to predict these individual components remains poor. In recent years, there has been substantial effort and success in the prediction of relative and absolute binding free energies, but the prediction of the enthalpic (and entropic) contributions in biomolecular systems remains challenging. Indeed, it is not even clear what kind of performance in terms of accuracy could currently be obtained for such systems. It is, however, relatively straight-forward to compute the enthalpy of binding. We thus evaluated the performance of absolute enthalpy of binding calculations using molecular dynamics simulation for ten inhibitors against a member of the bromodomain family, BRD4-1, against isothermal titration calorimetry data. Initial calculations, with the AMBER force-field showed good agreement with experiment (R2 = 0.60) and surprisingly good accuracy with an average of root-mean-square error (RMSE) = 2.49 kcal mol-1. Of the ten predictions, three were obvious outliers that were all over-predicted compared to experiment. Analysis of various simulation factors, including parameterization, buffer concentration and conformational dynamics, revealed that the behaviour of a loop (the ZA loop on the periphery of the binding site) strongly dictates the enthalpic prediction. Consistent with previous observations, the loop exists in two distinct conformational states and by considering one or the other or both states, the prediction for the three outliers can be improved dramatically to the point where the R2 = 0.95 and the accuracy in terms of RMSE improves to 0.90 kcal mol-1. However, performance across force-fields is not consistent: if OPLS and CHARMM are used, different outliers are observed and the correlation with the ZA loop behaviour is not recapitulated, likely reflecting parameterization as a confounding problem. The results provide a benchmark standard for future study and comparison.

Effects of Salinity and Temperature on the Flexibility and Function of a Polyextremophilic Enzyme.

The polyextremophilic β-galactosidase enzyme of the haloarchaeon Halorubrum lacusprofundi functions in extremely cold and hypersaline conditions. To better understand the basis of polyextremophilic activity, the enzyme was studied using steady-state kinetics and molecular dynamics at temperatures ranging from 10 °C to 50 °C and salt concentrations from 1 M to 4 M KCl. Kinetic analysis showed that while catalytic efficiency (kcat/Km) improves with increasing temperature and salinity, Km is reduced with decreasing temperatures and increasing salinity, consistent with improved substrate binding at low temperatures. In contrast, kcat was similar from 2-4 M KCl across the temperature range, with the calculated enthalpic and entropic components indicating a threshold of 2 M KCl to lower the activation barrier for catalysis. With molecular dynamics simulations, the increase in per-residue root-mean-square fluctuation (RMSF) was observed with higher temperature and salinity, with trends like those seen with the catalytic efficiency, consistent with the enzyme's function being related to its flexibility. Domain A had the smallest change in flexibility across the conditions tested, suggesting the adaptation to extreme conditions occurs via regions distant to the active site and surface accessible residues. Increased flexibility was most apparent in the distal active sites, indicating their importance in conferring salinity and temperature-dependent effects.

Enthalpic and entropic contributions to the activation free energy of single noncovalent bonds in molecular systems: A computational methodology

Insight into enthalpic and entropic contributions to the activation free energy of noncovalent bonds is crucial for understanding various phenomena in molecular systems, from formation of different molecular conformers to complex biophysical and biochemical processes. Standard experimental and computational methods provide only the average strength of noncovalent bonds. In this paper, we present a simple, formal, computational framework to determine both energy contributions for a single bond. Our approach is based on steered molecular dynamics simulations, the dynamic force spectroscopy method, and a modified theoretical model of force-induced bond rupture. To demonstrate the methodology, the enthalpic and entropic components for a single intramolecular hydrogen bond of a peptide helix were determined. The sum of the contributions is in agreement with the previously reported Gibbs energy for intramolecular hydrogen bonds. The proposed methodological framework is universal and can be applied to other noncovalent interactions in various molecular systems.

Physical insight into the entropy-driven ion association.

The ion association is widely believed to be dominated by the favorable entropy change arising from the release of water molecules from ion hydration shells. However, no direct thermodynamic evidence exists to validate the reliability and suitability of this view. Herein, we employ complicated free energy calculations to rigorously split the free energy including its entropic and enthalpic components into the water-induced contributions and ion-ion interaction terms for several ion pairs from monatomic to polyatomic ions, spanning the size range from small kosmotropes to large chaotropes (Na+ , Cs+ , Ca2+ , F- , I- , CO3 2- , and HPO4 2- ). Our results successfully reveal that though ion associations are indeed determined by a delicate balance between the favorable entropy variation and the repulsive enthalpy change, the entropy gain dominated by the solvent occurs only for the monatomic ion pairing. The water-induced entropic contribution significantly goes against the ion pairing between polyatomic anion and cation, which is, alternatively, dominated by the favorable entropy from the ion-ion interaction term, due to the configurational arrangement of polyatomic anions involved in ion association. The structural and dynamic analysis demonstrates that the entropy penalty from the water phase is primarily ascribed to the enhanced stability of water molecules around the cation imposed by the incoming anion. Our study successfully provides a fundamental understanding of water-mediated ion associations and highlights disparate lengthscale dependencies of the dehydration thermodynamics on the specific types of ions.

Macromolecular Crowding Is More than Hard-Core Repulsions.

Cells are crowded, but proteins are almost always studied in dilute aqueous buffer. We review the experimental evidence that crowding affects the equilibrium thermodynamics of protein stability and protein association and discuss the theories employed to explain these observations. In doing so, we highlight differences between synthetic polymers and biologically relevant crowders. Theories based on hard-core interactions predict only crowding-induced entropic stabilization. However, experiment-based efforts conducted under physiologically relevant conditions show that crowding can destabilize proteins and their complexes. Furthermore, quantification of the temperature dependence of crowding effects produced by both large and small cosolutes, including osmolytes, sugars, synthetic polymers, and proteins, reveals enthalpic effects that stabilize or destabilize proteins.Crowding-induced destabilization and the enthalpic component point to the role of chemical interactions between and among the macromolecules, cosolutes, and water. We conclude with suggestions for future studies.

Effect of Surface Nanostructures and Speciation on Undercooling for Low-Temperature Solder Alloys

Role of surface structure/composition in altering solidification behavior is demonstrated through undercooling of core–shell metal particles. Autonomous surface speciation in the passivating oxide plays a critical role in the relaxation of a molten metal due to divergence in the concentration and composition of oxidizing species across the thickness of the oxide layer. In an unreactive environment, surface speciation is dictated by flux, cohesive energy density, and surface energy minimization. Under oxidizing conditions (e.g., ambient), however, reduction potential, curvature, and surface plasticity dictate the spatial order and concentration(s) across thin passivating oxide layers. It is therefore important to redefine solubility beyond the limitations of Hume–Rothery rules by substituting electronegativity for redox potential and cohesive energy density. Increasing number of components in an alloy does not necessarily lead to increased undercooling, but a maximum is observed around two to three components. These results lead to an empirical observation that ΔGLS (Gibbs’ free energy for liquid–solid transition) around the freezing point can be understood from enthalpy and surface tension balance, with the degree of undercooling being a proportionality term for the enthalpic component. This work extends our understanding of classical nucleation theory by illustrating the importance of asymmetry in surface work in phase change, largely due to changes in structure of nanoscale passivating oxides. This phenomena is critical in development of low-temperature solders.

Exploring Unconventional SAM Analogues To Build Cell-Potent Bisubstrate Inhibitors for Nicotinamide N-Methyltransferase.

Nicotinamide N-methyltransferase (NNMT) methylates nicotinamide and has been associated with various diseases. Herein, we report the first cell-potent NNMT bisubstrate inhibitor II399, demonstrating a Ki of 5.9 nM in a biochemical assay and a cellular IC50 value of 1.9 μM. The inhibition mechanism and cocrystal structure confirmed II399 engages both the substrate and cofactor binding pockets. Computational modeling and binding data reveal a balancing act between enthalpic and entropic components that lead to II399's low nM binding affinity. Notably, II399 is 1 000-fold more selective for NNMT than closely related methyltransferases. We expect that II399 would serve as a valuable probe to elucidate NNMT biology. Furthermore, this strategy provides the first case of introducing unconventional SAM mimics, which can be adopted to develop cell-potent inhibitors for other SAM-dependent methyltransferases.