Nanoconfinement of Carbon Dioxide within Interfacial Aqueous/Ionic Liquid Systems.

Nanoporous, gas-selective membranes have shown encouraging results for the removal of CO2 from flue gas, yet the optimal design for such membranes is often unknown. Therefore, we used molecular dynamics simulations to elucidate the behavior of CO2 within aqueous and ionic liquid (IL) systems ([EMIM][TFSI] and [OMIM][TFSI]), both confined individually and as an interfacial aqueous/IL system. We found that within aqueous systems the mobility of CO2 is reduced due to interactions between the CO2 oxygens and hydroxyl groups on the pore surface. Within the IL systems, we found that confinement has a greater effect on the [EMIM][TFSI] system as opposed to the [OMIM][TFSI] system. Paradoxically, the larger and more asymmetrical [OMIM]+ molecule undergoes less efficient packing, resulting in fewer confinement effects. Free energy surfaces of the nanoconfined aqueous/IL interface demonstrate that CO2 will transfer spontaneously from the aqueous to the IL phase.

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory.

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here, we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F- over Cl- in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory-experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

Energetics of high temperature degradation of fentanyl into primary and secondary products.

Fentanyl is a synthetic opioid used for managing chronic pain. Due to its higher potency (50-100×) than morphine, fentanyl is also an abused drug. A sensor that could detect illicit fentanyl by identifying its thermally degraded fragments would be helpful to law enforcement. While experimental studies have probed the thermal degradation of fentanyl, little theoretical work has been done to understand the mechanism. Here, we studied the thermal degradation pathways of fentanyl using extensive ab initio molecular dynamics simulations combined with enhanced sampling via multiple-walker metadynamics. We calculated the free energy profile for each bond suggested earlier as a potential degradation point to map the thermodynamic driving forces. We also estimated the forward attempt rate of each bond degradation reaction to gain information about degradation kinetics.

Molecular Cage Reports on Its Contents: Spectroscopic Signatures of Cryo-Cooled K+- and Ba2+-Benzocryptand Complexes.

UV photofragment spectroscopy and IR-UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K+ or Ba2+ imbedded inside it (labeled as K+-BzCrypt, Ba2+-BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH2-CH2-O-CH2-CH2-O-CH2-CH2-N chains provides a UV chromophore whose S0-S1 transition we probe. K+-BzCrypt and Ba2+-BzCrypt have their S0-S1 origin transitions at 35,925 and 36,446 cm-1, respectively, blue-shifted by 174 and 695 cm-1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR-UV dip spectroscopy. The alkyl CH stretch region (2800-3000 cm-1) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K+-BzCrypt and Ba2+-BzCrypt, identifying two conformations (A and B) that lie below all others in energy. We extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH2 groups adjacent to these atoms. The localized CH2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm-1. As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.

Channelrhodopsin C1C2: Photocycle kinetics and interactions near the central gate.

Channelrhodopsins (ChR) are light-sensitive cation channels used in optogenetics, a technique that applies light to control cells (e.g., neurons) that have been modified genetically to express those channels. Although mutations are known to affect pore kinetics, little is known about how mutations induce changes at the molecular scale. To address this issue, we first measured channel opening and closing rates of a ChR chimera (C1C2) and selected variants (N297D, N297V, and V125L). Then, we used atomistic simulations to correlate those rates with changes in pore structure, hydration, and chemical interactions among key gating residues of C1C2 in both closed and open states. Overall, the experimental results show that C1C2 and its mutants do not behave like ChR2 or its analogous variants, except V125L, making C1C2 a unique channel. Our atomistic simulations confirmed that opening of the channel and initial hydration of the gating regions between helices I, II, III, and VII of the channel occurs with 1) the presence of 13-cis retinal; 2) deprotonation of a glutamic acid gating residue, E129; and 3) subsequent weakening of the central gate hydrogen bond between the same glutamic acid E129 and asparagine N297 in the central region of the pore. Also, an aspartate (D292) is the unambiguous primary proton acceptor for the retinal Schiff base in the hydrated channel.

Hydration Mimicry by Membrane Ion Channels.

Ions transiting biomembranes might pass readily from water through ion-specific membrane proteins if these protein channels provide environments similar to the aqueous solution hydration environment. Indeed, bulk aqueous solution is an important reference condition for the ion permeation process. Assessment of this hydration mimicry concept depends on understanding the hydration structure and free energies of metal ions in water in order to provide a comparison for the membrane channel environment. To refine these considerations, we review local hydration structures of ions in bulk water and the molecular quasi-chemical theory that provides hydration free energies. In doing so, we note some current views of ion binding to membrane channels and suggest new physical chemical calculations and experiments that might further clarify the hydration mimicry concept.

Partitioning of Seven Different Classes of Antibiotics into LPS Monolayers Supports Three Different Permeation Mechanisms through the Outer Bacterial Membrane.

The outer membrane (OM) of Gram-negative (G-) bacteria presents a barrier for many classes of antibacterial agents. Lipopolysaccharide (LPS), present in the outer leaflet of the OM, is stabilized by divalent cations and is considered to be the major impediment for antibacterial agent permeation. However, the actual affinities of major antibiotic classes toward LPS have not yet been determined. In the present work, we use Langmuir monolayers formed from E. coli Re and Rd types of LPS to record pressure-area isotherms in the presence of antimicrobial agents. Our observations suggest three general types of interactions. First, some antimicrobials demonstrated no measurable interactions with LPS. This lack of interaction in the case of cefsulodin, a third-generation cephalosporin antibiotic, correlates with its low efficacy against G- bacteria. Ampicillin and ciprofloxacin also show no interactions with LPS, but in contrast to cefsulodin, both exhibit good efficacy against G- bacteria, indicating permeation through common porins. Second, we observe substantial intercalation of the more hydrophobic antibiotics, novobiocin, rifampicin, azithromycin, and telithromycin, into relaxed LPS monolayers. These largely repartition back to the subphase with monolayer compression. We find that the hydrophobic area, charge, and dipole all show correlations with both the mole fraction of antibiotic retained in the monolayer at the monolayer-bilayer equivalence pressure and the efficacies of these antibiotics against G- bacteria. Third, amine-rich gentamicin and the cationic antimicrobial peptides polymyxin B and colistin show no hydrophobic insertion but are instead strongly driven into the polar LPS layer by electrostatic interactions in a pressure-independent manner. Their intercalation stably increases the area per molecule (by up to 20%), which indicates massive formation of defects in the LPS layer. These defects support a self-promoted permeation mechanism of these antibiotics through the OM, which explains the high efficacy and specificity of these antimicrobials against G- bacteria.

Bio-inspired incorporation of phenylalanine enhances ionic selectivity in layer-by-layer deposited polyelectrolyte films.

The addition of a common amino acid, phenylalanine, to a Layer-by-Layer (LbL) deposited polyelectrolyte (PE) film on a nanoporous membrane can increase its ionic selectivity over a PE film without the added amino acid. The addition of phenylalanine is inspired by detailed knowledge of the structure of the channelrhodopsins family of protein ion channels, where phenylalanine plays an instrumental role in facilitating sodium ion transport. The normally deposited and crosslinked PE films increase the cationic selectivity of a support membrane in a controllable manner where higher selectivity is achieved with thicker PE coatings, which in turn also increases the ionic resistance of the membrane. The increased ionic selectivity is desired while the increased resistance is not. We show that through incorporation of phenylalanine during the LbL deposition process, in solutions of NaCl with concentrations ranging from 0.1 to 100 mM, the ionic selectivity can be increased independently of the membrane resistance. Specifically, the addition is shown to increase the cationic transference of the PE films from 81.4% to 86.4%, an increase on par with PE films that are nearly triple the thickness while exhibiting much lower resistance compared to the thicker coatings, where the phenylalanine incorporated PE films display an area specific resistance of 1.81 Ω cm2 in 100 mM NaCl while much thicker PE membranes show a higher resistance of 2.75 Ω cm2 in the same 100 mM NaCl solution.

Ab initio and force field molecular dynamics study of bulk organophosphorus and organochlorine liquid structures.

We performed ab initio molecular dynamics (AIMD) simulations to benchmark bulk liquid structures and to evaluate results from all-atom force field molecular dynamics (FFMD) simulations with the generalized Amber force field (GAFF) for organophosphorus (OP) and organochlorine (OC) compounds. Our work also addresses the current and important topic of force field validation, applied here to a set of nonaqueous organic liquids. Our approach differs from standard treatments, which validate force fields based on thermodynamic data. Utilizing radial distribution functions (RDFs), our results show that GAFF reproduces the AIMD-predicted asymmetric liquid structures moderately well for OP compounds that contain bulky alkyl groups. Among the OCs, RDFs obtained from FFMD overlap well with AIMD results, with some offsets in position and peak structuring. However, re-parameterization of GAFF for some OCs is needed to reproduce fully the liquid structures predicted by AIMD. The offsets between AIMD and FFMD peak positions suggest inconsistencies in the developed force fields, but, in general, GAFF is able to capture short-ranged and long-ranged interactions of OPs and OCs observed in AIMD. Along with the local coordination structure, we also compared enthalpies of vaporization. Overall, calculated bulk properties from FFMD compared reasonably well with experimental values, suggesting that small improvements within the FF should focus on parameters that adjust the bulk liquid structures of these compounds.

Electrostatic lock in the transport cycle of the multidrug resistance transporter EmrE.

EmrE is a small, homodimeric membrane transporter that exploits the established electrochemical proton gradient across the Escherichia coli inner membrane to export toxic polyaromatic cations, prototypical of the wider small-multidrug resistance transporter family. While prior studies have established many fundamental aspects of the specificity and rate of substrate transport in EmrE, low resolution of available structures has hampered identification of the transport coupling mechanism. Here we present a complete, refined atomic structure of EmrE optimized against available cryo-electron microscopy (cryo-EM) data to delineate the critical interactions by which EmrE regulates its conformation during the transport process. With the model, we conduct molecular dynamics simulations of the transporter in explicit membranes to probe EmrE dynamics under different substrate loading and conformational states, representing different intermediates in the transport cycle. The refined model is stable under extended simulation. The water dynamics in simulation indicate that the hydrogen-bonding networks around a pair of solvent-exposed glutamate residues (E14) depend on the loading state of EmrE. One specific hydrogen bond from a tyrosine (Y60) on one monomer to a glutamate (E14) on the opposite monomer is especially critical, as it locks the protein conformation when the glutamate is deprotonated. The hydrogen bond provided by Y60 lowers the [Formula: see text] of one glutamate relative to the other, suggesting both glutamates should be protonated for the hydrogen bond to break and a substrate-free transition to take place. These findings establish the molecular mechanism for the coupling between proton transfer reactions and protein conformation in this proton-coupled secondary transporter.

Enhancing Paraoxon Binding to Organophosphorus Hydrolase Active Site.

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.

Quantum Calculations of VX Ammonolysis and Hydrolysis Pathways via Hydrated Lithium Nitride.

Recently, lithium nitride (Li3N) has been proposed as a chemical warfare agent (CWA) neutralization reagent for its ability to produce nucleophilic ammonia molecules and hydroxide ions in aqueous solution. Quantum chemical calculations can provide insight into the Li3N neutralization process that has been studied experimentally. Here, we calculate reaction-free energies associated with the Li3N-based neutralization of the CWA VX using quantum chemical density functional theory and ab initio methods. We find that alkaline hydrolysis is more favorable to either ammonolysis or neutral hydrolysis for initial P-S and P-O bond cleavages. Reaction-free energies of subsequent reactions are calculated to determine the full reaction pathway. Notably, products predicted from favorable reactions have been identified in previous experiments.

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes.

With a longer-term goal of addressing the comparative behavior of the aqueous halides F-, Cl-, Br-, and I- on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl-, Br-, and I-, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

First-principles modeling of chemistry in mixed solvents: Where to go from here?

Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.

Insertion of Dengue E into lipid bilayers studied by neutron reflectivity and molecular dynamics simulations.

The envelope (E) protein of Dengue virus rearranges to a trimeric hairpin to mediate fusion of the viral and target membranes, which is essential for infectivity. Insertion of E into the target membrane serves to anchor E and possibly also to disrupt local order within the membrane. Both aspects are likely to be affected by the depth of insertion, orientation of the trimer with respect to the membrane normal, and the interactions that form between trimer and membrane. In the present work, we resolved the depth of insertion, the tilt angle, and the fundamental interactions for the soluble portion of Dengue E trimers (sE) associated with planar lipid bilayer membranes of various combinations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-rac-glycerol (POPG), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and cholesterol (CHOL) by neutron reflectivity (NR) and by molecular dynamics (MD) simulations. The results show that the tip of E containing the fusion loop (FL) is located at the interface of the headgroups and acyl chains of the outer leaflet of the lipid bilayers, in good agreement with prior predictions. The results also indicate that E tilts with respect to the membrane normal upon insertion, promoted by either the anionic lipid POPG or CHOL. The simulations show that tilting of the protein correlates with hydrogen bond formation between lysines and arginines located on the sides of the trimer close to the tip (K246, K247, and R73) and nearby lipid headgroups. These hydrogen bonds provide a major contribution to the membrane anchoring and may help to destabilize the target membrane.

Strontium and barium in aqueous solution and a potassium channel binding site.

Ion hydration structure and free energy establish criteria for understanding selective ion binding in potassium (K+) ion channels and may be significant to understanding blocking mechanisms as well. Recently, we investigated the hydration properties of Ba2+, the most potent blocker of K+ channels among the simple metal ions. Here, we use a similar method of combining ab initio molecular dynamics simulations, statistical mechanical theory, and electronic structure calculations to probe the fundamental hydration properties of Sr2+, which does not block bacterial K+ channels. The radial distribution of water around Sr2+ suggests a stable 8-fold geometry in the local hydration environment, similar to Ba2+. While the predicted hydration free energy of -331.8 kcal/mol is comparable with the experimental result of -334 kcal/mol, the value is significantly more favorable than the -305 kcal/mol hydration free energy of Ba2+. When placed in the innermost K+ channel blocking site, the solvation free energies and lowest energy structures of both Sr2+ and Ba2+ are nearly unchanged compared with their respective hydration properties. This result suggests that the block is not attributable to ion trapping due to +2 charge, and differences in blocking behavior arise due to free energies associated with the exchange of water ligands for channel ligands instead of free energies of transfer from water to the binding site.

Atomistic Study of Intramolecular Interactions in the Closed-State Channelrhodopsin Chimera, C1C2.

Channelrhodopsins (ChR1 and ChR2) are light-activated ion channels that enable photomobility of microalgae from the genus Chlamydomonas. Despite common use of ChR2 in optogenetics for selective control and monitoring of individual neurons in living tissue, the protein structures remain unresolved. Instead, a crystal structure of the ChR chimera (C1C2), an engineered combination of helices I-V from ChR1, without its C-terminus, and helices VI-VII from ChR2, is used as a template for ChR2 structure prediction. Surprisingly few studies have focused in detail on the chimera. Here, we present atomistic molecular dynamics studies of the closed-state, non-conducting C1C2 structure and protonation states. A new and comprehensive characterization of interactions in the vicinity of the gating region of the pore, namely between residues E90, E123, D253, N258, and the protonated Schiff base (SBH), as well as nearby residues K93, T127, and C128, indicates that the equilibrated C1C2 structure with both E123 and D253 deprotonated closely resembles the available crystal structure. In agreement with experimental studies on C1C2, no direct or water-mediated hydrogen bonding between an aspartate and a cysteine (D156-O…S-C128) that would define a direct-current gate in C1C2 was observed in our simulations. Finally, we show that a single hydrogen bond between a glutamic acid (E90) and an asparagine (N258) residue suffices to keep the gate of C1C2 closed and to disable free water and ion passage through the putative pore, in contrast to the double bond proposed earlier for ChR2. We anticipate that this work will provide context for studies of both the gating process and water and ion transport in C1C2, and will spark interest in further experimental studies on the chimera.

Method for measuring the unbinding energy of strongly-bound membrane-associated proteins.

We describe a new method to measure the activation energy for unbinding (enthalpy ΔH*u and free energy ΔG*u) of a strongly-bound membrane-associated protein from a lipid membrane. It is based on measuring the rate of release of a liposome-bound protein during centrifugation on a sucrose gradient as a function of time and temperature. The method is used to determine ΔH*u and ΔG*u for the soluble dengue virus envelope protein (sE) strongly bound to 80:20 POPC:POPG liposomes at pH5.5. ΔH*u is determined from the Arrhenius equation whereas ΔG*u is determined by fitting the data to a model based on mean first passage time for escape from a potential well. The binding free energy ΔGb of sE was also measured at the same pH for the initial, predominantly reversible, phase of binding to a 70:30 PC:PG lipid bilayer. The unbinding free energy (20±3kcal/mol, 20% PG) was found to be roughly three times the binding energy per monomer, (7.8±0.3kcal/mol for 30% PG, or est. 7.0kcal/mol for 20% PG). This is consistent with data showing that free sE is a monomer at pH5.5, but assembles into trimers after associating with membranes. This new method to determine unbinding energies should be useful to understand better the complex interactions of integral monotopic proteins and strongly-bound peripheral membrane proteins with lipid membranes.

The Impact of Protonation on Early Translocation of Anthrax Lethal Factor: Kinetics from Molecular Dynamics Simulations and Milestoning Theory.

We report atomically detailed molecular dynamics simulations of the permeation of the lethal factor (LF) N-terminal segment through the anthrax channel. The N-terminal chain is unstructured and leads the permeation process for the LF protein. The simulations were conducted in explicit solvent with milestoning theory, making it possible to extract kinetic information from nanosecond to millisecond time scales. We illustrate that the initial event is strongly influenced by the protonation states of the permeating amino acids. While the N-terminal segment passes easily at high protonation state through the anthrax channel (and the ϕ clamp), the initial permeation represents a critical step, which can be irreversible and establishes a hook in the channel mouth.

Quasi-chemical theory of F-(aq): The "no split occupancies rule" revisited.

We use ab initio molecular dynamics (AIMD) calculations and quasi-chemical theory (QCT) to study the inner-shell structure of F-(aq) and to evaluate that single-ion free energy under standard conditions. Following the "no split occupancies" rule, QCT calculations yield a free energy value of -101 kcal/mol under these conditions, in encouraging agreement with tabulated values (-111 kcal/mol). The AIMD calculations served only to guide the definition of an effective inner-shell constraint. QCT naturally includes quantum mechanical effects that can be concerning in more primitive calculations, including electronic polarizability and induction, electron density transfer, electron correlation, molecular/atomic cooperative interactions generally, molecular flexibility, and zero-point motion. No direct assessment of the contribution of dispersion contributions to the internal energies has been attempted here, however. We anticipate that other aqueous halide ions might be treated successfully with QCT, provided that the structure of the underlying statistical mechanical theory is absorbed, i.e., that the "no split occupancies" rule is recognized.