Modeling Enantiomeric Separations as an Interfacial Process Using Amylose Tris(3,5-dimethylphenyl carbamate) (ADMPC) Polymers Coated on Amorphous Silica.

Chiral high-performance liquid chromatography (HPLC) is commonly performed to isolate the biologically active enantiomer of a drug from the ineffective or even harmful ones. Understanding the molecular-level recognition that underlies this process is necessary for trimming down the very large number of possible combinations of chiral stationary phases, solvent systems, and other experimental HPLC conditions, a particularly important consideration when only small quantities of the racemate are available. Fully atomistic molecular dynamics (MD) simulation is a useful tool to provide this molecular-level understanding and predict experimental separation factors under a given set of conditions. To predict the chiral separation results for drug enantiomers by amylose tris(3,5-dimethylphenyl carbamate) (ADMPC) chiral stationary phase, we design a model of multiple ADMPC polymer strands coated on an amorphous silica slab. Using various MD techniques, we successfully coat ADMPCs onto the surface without losing the structural character of the backbone in the presence of the solvent system. Not only is this model more representative of the polymer surface on a solid support that is encountered by the enantiomers, but it also provides more opportunities for chiral molecules interacting with ADMPC, provides the possibility for large drug molecules to interact with two polymer strands at the same instant, and provides better agreement with experiment when we use the overall average quantities as the predictive metric. For a better understanding of why some metrics are better predictors than others, we use charts of the distribution of hydrogen-bonding lifetimes for various donor-acceptor pairs that contribute to the interaction events determining the relative retention times for the enantiomers. We also examine the contribution of ring-ring interactions to the molecular recognition process and ultimately to the differential retention of enantiomers. The results are more consistent than previous models and resolve the problematic case of two drugs, thalidomide and valsartan.

The Composition of the Mobile Phase Affects the Dynamic Chiral Recognition of Drug Molecules by the Chiral Stationary Phase.

More than half of all pharmaceuticals are chiral compounds. Although the enantiomers of chiral compounds have the same chemical structure, they can exhibit marked differences in physiological activity; therefore, it is important to remove the undesirable enantiomer. Chromatographic separation of chiral enantiomers is one of the best available methods to get enantio-pure substances, but the optimization of the experimental conditions can be very time-consuming. One of the most widely used chiral stationary phases, amylose tris(3,5-dimethylphenyl carbamate) (ADMPC), has been extensively investigated using both experimental and computational methods; however, the dynamic nature of the interaction between enantiomers and ADMPC, as well as the solvent effects on the ADMPC-enantiomer interaction, are currently absent from models of the chiral recognition mechanism. Here we use QM/MM and molecular dynamics (MD) simulations to model the enantiomers of flavanone on ADMPC in either methanol or heptane/2-propanol (IPA) (90/10) to elucidate the chiral recognition mechanism from a new dynamic perspective. In atomistic MD simulations, the 12-mer model of ADMPC is found to hold the 4/3 left-handed helical structure in both methanol and heptane/IPA (90/10); however, the ADMPC polymer is found to have a more extended average structure in heptane/IPA (90/10) than in methanol. This results from the differences in the distribution of solvent molecules close to the backbone of ADMPC leads to changes in the distribution of the (φ, ψ) dihedral angles of the glycoside bond (between adjacent monomers) that define the structure of the polymer. Our simulations have shown that the lifetime of hydrogen bonds formed between ADMPC and flavanone enantiomers in the MD simulations are able to reproduce the elution order observed in experiments for both the methanol and the heptane/IPA solvent systems. Furthermore, the ratios of hydrogen-bonding-lifetime-related properties also capture the solvent effects, in that heptane/IPA (90/10) is found to make the separation between the two enantiomers of flavanone less effective than methanol, which agrees with the experimental separation factors of 0.9 versus 0.4 for R/S, respectively.

Simulated Permeation and Characterization of PEGylated Gold Nanoparticles in a Lipid Bilayer System.

PEGylated gold nanoparticles are considered suitable nanocarriers for use in biomedical applications and targeted drug delivery systems. In our previous investigation with the alkanethiol-functionalized gold nanoparticle, we found that permeation across a protein-free phospholipid membrane resulted in damaging effects of lipid displacement and water and ion leakage. In the present study, we carry out a series of coarse-grained molecular simulations to explore permeation of lipid bilayer systems by a PEGylated gold nanoparticle, especially at the bulk-liquid-lipid interface as well as the interface between the two lipid leaflets. Initially, we examine molecular-level details of a PEGylated gold nanoparticle (constructed from cycled annealing) in water and find a distribution of ligand configurations (from mushroom to brush states) present in nanoparticles with medium to high surface coverage. We also find that the characteristic properties of the PEGylated gold nanoparticle do not change when it is placed in a salt solution. In our permeation studies, we investigate events of water and ion penetration as well as lipid translocation while varying the ligand length, nanoparticle surface coverage, and ion concentration gradient of our system. Results from our studies show the following: (1) The number of water molecules in the interior of the membrane during ligand-coated nanoparticle permeation increases with PEGn-SH surface coverage, ligand length, and permeation velocity but is not sensitive to the ion concentration gradient. (2) Lipid molecules do not leave the membrane; instead they complete trans-bilayer lipid flip-flop with longer ligands and higher surface coverages. (3) The lack of formation of stable water pores prevents ion translocation. (4) The PEGylated nanoparticle causes less damage to the membrane overall due to favorable interactions with the lipid headgroups which may explain why experimentalists observe endocytosis of PEGylated nanocarriers in vivo.

Surface-functionalized nanoparticle permeation triggers lipid displacement and water and ion leakage.

Functionalized nanoparticles (NPs) are considered suitable carriers for targeted drug delivery systems. However, the ion and water leakage induced by permeation of these nanoparticles is a challenge in these drug delivery methods because of cytotoxic effects of some ions. In this study, we have carried out a series of coarse-grained molecular dynamics simulations to investigate the effect of length of ligands on permeation of a nanoparticle across a protein-free phospholipid bilayer membrane. Water and ion penetration as well as incidence of lipid flip-flop events and loss of lipid molecules from the membrane are explored in this study while varying the nanoparticle size, length of ligand, ion concentration gradient, pressure differential across the membrane, and nanoparticle permeation velocity. Some results from our studies include (1) the number of water molecules in the interior of the membrane during ligand-coated nanoparticle permeation increases with nanoparticle size, ligand length, pressure differential, and permeation velocity but is not sensitive to the ion concentration gradient; (2) some lipid molecules leave the membrane by being entangled with ligands of the NP instead of completing the flip-flop that permits them to rejoin the membrane, thereby leading to fewer flip-flop events; and (3) the formation of water columns or water "fingers" provides a mechanism of ion transport across lipid bilayer membranes, but such ion penetration events are less likely for sodium ions than chloride ions and less likely for nanoparticles with longer-ligands.

Xe Recovery from Nuclear Power Plants Off-Gas Streams: Molecular Simulations of Gas Permeation through DD3R Zeolite Membrane.

Recent experimental work has shown zeolite membrane-based separation as a promising potential technology for Kr/Xe gas mixtures due to its much lower energy requirements in comparison to cryogenic distillation, the conventional separation method for such mixtures. Such a separation is also economically rewarding because Xe is in high demand, as a valuable product for many applications/processes. In this work, we have used Molecular Dynamics (MD) simulations to study the effects of different conditions, i.e., temperature, pressure, and gas feed composition, on Kr/Xe separation performance via DD3R zeolite membranes. We provide a comprehensive study of the permeation of the different gas species, density profiles, and diffusion coefficients. Molecular simulations show that if the feed is changed from pure Kr/Xe to an equimolar mixture, the Kr/Xe separation factor increases, which agrees with experiments. In addition, when Ar is introduced as a sweep gas, the adsorption of both Kr and Xe increases, while the permeation of pure Kr increases. A similar behavior is observed with equimolar mixtures of Kr/Xe with Ar as the sweep gas. High-separation Kr/Xe selectivity is observed at 50 atm and 425 K but with low total permeation rates. Changing pressure and temperature are found to have profound effects on optimizing the separation selectivity and the permeation throughput.

Nanoparticle permeation induces water penetration, ion transport, and lipid flip-flop.

Nanoparticles are generally considered excellent candidates for targeted drug delivery. However, ion leakage and cytotoxicity induced by nanoparticle permeation is a potential problem in such drug delivery schemes because of the toxic effect of many ions. In this study, we have carried out a series of coarse-grained molecular dynamics simulations to investigate the water penetration, ion transport, and lipid molecule flip-flop in a protein-free phospholipid bilayer membrane during nanoparticle permeation. The effect of ion concentration gradient, pressure differential across the membrane, nanoparticle size, and permeation velocity have been examined in this work. Some conclusions from our studies include (1) The number of water molecules in the interior of the membrane during the nanoparticle permeation increases with the nanoparticle size and the pressure differential across the membrane but is unaffected by the nanoparticle permeation velocity or the ion concentration gradient. (2) Ion transport is sensitive to the size of nanoparticle as well as the ion concentration gradient between two sides of the membrane; no anion/cation selectivity is observed for small nanoparticle permeation, while anions are preferentially translocated through the membrane when the size of nanoparticle is large enough. (3) Incidences of lipid molecule flip-flop increases with the size of nanoparticle and ion concentration gradient and decreases with the pressure differential and the nanoparticle permeation velocity.

Novel methods to characterise spatial distribution and enantiomeric composition of usnic acids in four Icelandic lichens.

Usnic acid is an antibiotic metabolite produced by a wide variety of lichenized fungal lineages. The enantiomers of usnic acid have been shown to display contrasting bioactivities, and hence it is important to determine their spatial distribution, amounts and enantiomeric ratios in lichens to understand their roles in nature and grasp their pharmaceutical potential. The overall aim of the study was to characterise the spatial distribution of the predominant usnic acid enantiomer in lichens by combining spatial imaging and chiral chromatography. Specifically, separation and quantification of usnic acid enantiomers in four common lichens in Iceland was performed using a validated chiral chromatographic method. Molecular dynamics simulation was carried out to rationalize the chiral separation mechanism. Spatial distribution of usnic acid in the lichen thallus cross-sections were analysed using Desorption Electrospray Ionization-Imaging Mass Spectrometry (DESI-IMS) and fluorescence microscopy. DESI-IMS confirmed usnic acid as a cortical compound, and revealed that usnic acid can be more concentrated around the algal vicinity. Fluorescence microscopy complemented DESI-IMS by providing more detailed distribution information. By combining results from spatial imaging and chiral separation, we were able to visualize the distribution of the predominant usnic acid enantiomer in lichen cross-sections: (+)-usnic acid in Cladonia arbuscula and Ramalina siliquosa, and (-)-usnic acid in Alectoria ochroleuca and Flavocetraria nivalis. This study provides an analytical foundation for future environmental and functional studies of usnic acid enantiomers in lichens.

Molecular-Level "Observations" of the Behavior of Gold Nanoparticles in Aqueous Solution and Interacting with a Lipid Bilayer Membrane.

We use coarse-grained molecular dynamics simulations to "observe" details of interactions between ligand-covered gold nanoparticles and a lipid bilayer model membrane. In molecular dynamics simulations, one puts the individual atoms and groups of atoms of the physical system to be "observed" into a simulation box, specifies the forms of the potential energies of interactions between them (ultimately quantum based), and lets them individually move classically according to Newton's equations of motion, based on the forces arising from the assumed potential energy forms. The atoms that are chemically bonded to each other stay chemically bonded, following known potentials (force fields) that permit internal degrees of freedom (internal rotation, torsion, vibrations), and the interactions between nonbonded atoms are simplified to Lennard-Jones forms (in our case) and coulombic (where electrical charges are present) in which the parameters are previously optimized to reproduce thermodynamic properties or are based on quantum electronic calculations. The system is started out at a reasonable set of coordinates for all atoms or groups of atoms, and then permitted to develop according to the equations of motion, one small step (usually 10 fs time step) at a time, for millions of steps until the system is at a quasi-equilibrium (usually reached after hundreds of nanoseconds). We then let the system play out its motions further for many nanoseconds to observe the behavior, periodically taking snapshots (saving all positions and energies), and post-processing the snapshots to obtain various average descriptions of the system. Alkanethiols of various lengths serve as examples of hydrophobic ligands and methyl-terminated PEG with various numbers of monomer units serve as examples of hydrophilic ligands. Spherical gold particles of various diameters as well as gold nanorods form the core to which ligands are attached. The nanoparticles are characterized at the molecular level, especially the distributions of ligand configurations and their dependence on ligand length, and surface coverage. Self-assembly of the bilayer from an isotropic solution and observation of membrane properties that correspond well to experimental values validate the simulations. The mechanism of permeation of a gold NP coated with either a hydrophobic or a hydrophilic ligand, and its dependence on surface coverage, ligand length, core diameter, and core shape, is investigated. Lipid response such as lipid flip-flops, lipid extraction, and changes in order parameter of the lipid tails are examined in detail. The mechanism of permeation of a PEGylated nanorod is shown to occur by tilting, lying down, rotating, and straightening up. The nature of the information provided by molecular dynamics simulations permits understanding of the detailed behavior of gold nanoparticles interacting with lipid membranes which in turn helps to understand why some known systems work better than others and aids the design of new particles and improvement of methods for preparing existing ones.

Measuring chirality in NMR in the presence of a static electric field.

The scalar Hamiltonian of nuclear spins in the presence of a static electric field supports chirality. However, the eigenvalues of the Hamiltonian are not chiral; hence, chirality is not manifested in the usual NMR experiment. In this work, we show that the magnetization response to certain radio frequency pulse sequences exhibits chirality as well as handedness.

Ion permeation dynamics in carbon nanotubes.

Molecular dynamics simulations are carried out to investigate the permeation of ions and water in a membrane consisting of single wall carbon nanotubes possessing no surface charges connecting two reservoirs. Our simulations reveal that there are changes in the first hydration shell of the ions upon confinement in tubes of 0.82 or 0.90 nm effective internal diameter. Although the first minimum in the g(r) is barely changed in the nanotube compared to in the bulk solution, the hydration number of Na(+) ion is reduced by 1.0 (from 4.5 in bulk to 3.5 in the 0.90 nm tube) and the hydration number is reduced further in the 0.82 nm tube. The changes in the hydration shell of Cl(-) ion are negligible, within statistical errors. The water molecules of the first hydration shell of both ions exchange less frequently inside the tube than in the bulk solution. We compare ion trajectories for ions in the same tube under identical reservoir conditions but with different numbers of ions in the tubes. This permits investigation of changes in structure and dynamics which arise from multiple ion occupancy in a carbon nanotube possessing no surface charges. We also investigated the effects of tube flexibility. Ions enter the tubes so as to form a train of ion pairs. We find that the radial distribution profiles of Na(+) ions broaden significantly systematically with increasing number of ion pairs in the tube. The radial distribution profiles of Cl(-) ions change only slightly with increasing number of ions in the tube. Trajectories reveal that Na(+) ions do not pass each other in 0.90 nm tubes, while Cl(-) ions pass each other, as do ions of opposite charge. An ion entering the tube causes the like-charged ions preceding it in the tube to be displaced along the tube axis and positive or negative ions will exit the tube only when one or two other ions of the same charge are present in the tube. Thus, the permeation mechanism involves multiple ions and Coulomb repulsion among the ions plays an essential role.

A note on chirality in NMR spectroscopy.

Using simple symmetry arguments we give proofs of the derivations of the manifestation of chirality in the chemical shift and spin-spin coupling constant in nuclear magnetic resonance and relate our proofs to earlier discussions.

Xe nuclear magnetic resonance line shapes in channels decorated with paramagnetic centers.

To make predictions of the Xe NMR line shapes for Xe in channels decorated with paramagnetic centers, we consider a model system using the O(2) molecule as the paramagnetic center. The previously calculated quantum mechanical Xe@O(2) hyperfine tensor for various configurations of Xe in the presence of O(2) provides a model for the hyperfine response of Xe atom to the presence of a paramagnetic center. The averaging is carried out using the same grand canonical Monte Carlo methodology as for calculating NMR line shapes for Xe in diamagnetic channels, modified to include the effects of the hyperfine tensor response. We explore the temperature dependence of the Xe line shapes, the dependence on the concentration, and the symmetry of distribution of embedded paramagnetic centers, on the orientation of the paramagnetic center axis with respect to the channel axis, and on the radial distance of the paramagnetic center from the axis of the channel. We predict Xe line shape signatures of the presence and orientation of paramagnetic centers and deduce which tensor elements provide measures of concentration and radial distance of paramagnetic centers from the channel axis.

Molecular reorientation of CD(4) in gas-phase mixtures.

Spin-lattice relaxation times were measured for the deuterons in CD(4) in pure gas and in mixtures with the following buffer gases: Ar, Kr, Xe, HCl, N(2), CO, CO(2), CF(4), and SF(6). Effective collision cross sections sigma(theta, 2) for the molecular reorientation of CD(4) in collisions with these ten molecules are obtained as a function of temperature. These cross sections are compared with the corresponding cross sections sigma(J) obtained from (1)H spin-rotation relaxation in mixtures of CH(4) with the same set of buffer gases. Various classical reorientation models typically applied in liquids predict different ratios of the reduced correlation times for the reorientation of spherical tops. The Langevin model comes closest to predicting the magnitude of the sigma(theta, 2)/sigma(J) ratio that we obtain for CD(4).

Diastereomeric Xe chemical shifts in tethered cryptophane cages.

Cryptophane cages serve as host molecules to a Xe atom. Functionalization of cryptophane-A has permitted the development of Xe as a biosensor. Synthetic routes used to prepare cryptophanes result in racemic mixtures of the chiral cages. In the preparation of a tethered cryptophane-A cage for biosensor applications, some achiral and chiral substituents such as left-handed amino acids have been used. When the substituent is achiral, the NMR signal of the Xe atom in the functionalized cage in solution is a single isotropic peak, since the Xe shielding tensor components in the R and L cages differ by no more than the signs of the off-diagonal elements. Chiral substituents can split the cage-encapsulated Xe NMR signal into one or more sets of doublets, depending on the number of asymmetric centers in the substituent. We carry out quantum mechanical calculations of Xe nuclear magnetic shielding for the Xe atom at the same strategic position within an L cryptophane-A cage, under the influence of chiral potentials that represent r or l substituents outside the cage. Calculations of the Xe shielding response in the Lr and Ll diastereomeric pairs permit the prediction of the relative order of the Xe chemical shifts in solutions containing the Rl and Ll diastereomers. Where the substituent itself possesses two chiral centers, comparison of the calculated isotropic shielding responses in the Llr, Lrl, Rll, and Lrr systems, respectively, permits the prediction of the Xe spectrum of diastereomeric systems in solutions containing Llr, Rlr, Lll, and Rll systems. Assignment of the peaks observed in the experimental Xe NMR spectra is therefore possible, without having to undertake the difficult synthetic route that produces a single optically pure enantiomer.

Xe chemical shift tensor in silicalite and SSZ-24.

We report, for the first time, a theoretical prediction of the (129)Xe nuclear magnetic resonance chemical shift tensor of xenon atom in a single crystal of silicalite at near-zero occupancy and the temperature dependence of the Xe NMR chemical shift tensor for the polycrystalline silicalite at maximum occupancy. The former is a measure of the sensitivity of the Xe tensor components to the local structure of the channels without Xe-Xe contributions. The latter is a measure of the sensitivity of the Xe-Xe tensor components to the Xe-Xe distributions, as determined by the Xe-Xe potential function in competition with the Xe-silicalite potential function. Both theoretical predictions can be compared against Xe NMR experiments: the first against the Xe spectra collected as a function of rotation of the single crystal about the three crystalline axes in a magnetic field, and the second against variable temperature Xe NMR studies (below room temperature) of polycrystalline silicalite at maximum Xe occupancy. With the same parameter set (Xe-O potential and shielding functions), we predict the line shapes of Xe in SSZ-24 zeolite under various conditions of occupancy and temperature.

The nuclear magnetic resonance line shapes of Xe in the cages of clathrate hydrates.

We report, for the first time, a prediction of the line shapes that would be observed in the (129)Xe nuclear magnetic resonance (NMR) spectrum of xenon in the cages of clathrate hydrates. We use the dimer tensor model to represent pairwise contributions to the intermolecular magnetic shielding tensor for Xe at a specific location in a clathrate cage. The individual tensor components from quantum mechanical calculations in clathrate hydrate structure I are represented by contributions from parallel and perpendicular tensor components of Xe-O and Xe-H dimers. Subsequently these dimer tensor components are used to reconstruct the full magnetic shielding tensor for Xe at an arbitrary location in a clathrate cage. The reconstructed tensors are employed in canonical Monte Carlo simulations to find the Xe shielding tensor component along a particular magnetic field direction. The shielding tensor component weighted according to the probability of finding a crystal fragment oriented along this direction in a polycrystalline sample leads to a predicted line shape. Using the same set of Xe-O and Xe-H shielding functions and the same Xe-O and Xe-H potential functions we calculate the Xe NMR spectra of Xe atom in 12 distinct cage types in clathrate hydrates structures I, II, H, and bromine hydrate. Agreement with experimental spectra in terms of the number of unique tensor components and their relative magnitudes is excellent. Agreement with absolute magnitudes of chemical shifts relative to free Xe atom is very good. We predict the Xe line shapes in two cages in which Xe has not yet been observed.

The chemical shifts of Xe in the cages of clathrate hydrate Structures I and II.

We report, for the first time, a calculation of the isotropic NMR chemical shift of 129Xe in the cages of clathrate hydrates Structures I and II. We generate a shielding surface for Xe in the clathrate cages by quantum mechanical calculations. Subsequently this shielding surface is employed in canonical Monte Carlo simulations to find the average isotropic Xe shielding values in the various cages. For the two types of cages in clathrate hydrate Structure I, we find the intermolecular shielding values [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-214.0 ppm, and [sigma(Xe@5(12)6(2) cage)-sigma(Xe atom)]=-146.9 ppm, in reasonable agreement with the values -242 and -152 ppm, respectively, observed experimentally by Ripmeester and co-workers between 263 and 293 K. For the 5(12) and 5(12)6(4) cages of Structure II we find [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-206.7 ppm, and [sigma(Xe@5(12)6(4) cage)-sigma(Xe atom)]=-104.7 ppm, also in reasonable agreement with the values -225 and -80 ppm, respectively, measured in a Xe-propane type II mixed clathrate hydrate at 77 and 220-240 K by Ripmeester et al.

The Xe shielding surfaces for Xe interacting with linear molecules and spherical tops.

The 129Xe nuclear magnetic resonance spectrum of xenon in gas mixtures of Xe with other molecules provides a test of the ab initio surfaces for the intermolecular shielding of Xe in the presence of the other molecule. We examine the electron correlation contributions to the Xe-CO2, Xe-N2, Xe-CO, Xe-CH4, and Xe-CF4 shielding surfaces and test the calculations against the experimental temperature dependence of the density coefficients of the Xe chemical shift in the gas mixtures at infinite dilution in Xe. Comparisons with the gas phase data permit the refinement of site-site potential functions for Xe-N2, Xe-CO, and Xe-CF4 especially for atom-Xe distances in the range 3.5-6 A. With the atom-atom shielding surfaces and potential parameters obtained in the present work, construction of shielding surfaces and potentials for applications such as molecular dynamics averaging of Xe chemical shifts in liquid solvents containing CH3, CH2, CF3, and CF2 groups is possible.

Molecular dynamics averaging of Xe chemical shifts in liquids.

The Xe nuclear magnetic resonance chemical shift differences that afford the discrimination between various biological environments are of current interest for biosensor applications and medical diagnostic purposes. In many such environments the Xe signal appears close to that in water. We calculate average Xe chemical shifts (relative to the free Xe atom) in solution in eleven liquids: water, isobutane, perfluoro-isobutane, n-butane, n-pentane, neopentane, perfluoroneopentane, n-hexane, n-octane, n-perfluorooctane, and perfluorooctyl bromide. The latter is a liquid used for intravenous Xe delivery. We calculate quantum mechanically the Xe shielding response in Xe-molecule van der Waals complexes, from which calculations we develop Xe (atomic site) interpolating functions that reproduce the ab initio Xe shielding response in the complex. By assuming additivity, these Xe-site shielding functions can be used to calculate the shielding for any configuration of such molecules around Xe. The averaging over configurations is done via molecular dynamics (MD). The simulations were carried out using a MD technique that one of us had developed previously for the simulation of Henry's constants of gases dissolved in liquids. It is based on separating a gaseous compartment in the MD system from the solvent using a semipermeable membrane that is permeable only to the gas molecules. We reproduce the experimental trends in the Xe chemical shifts in n-alkanes with increasing number of carbons and the large chemical shift difference between Xe in water and in perfluorooctyl bromide. We also reproduce the trend for a given solvent of decreasing Xe chemical shift with increasing temperature. We predict chemical shift differences between Xe in alkanes vs their perfluoro counterparts.

Nuclear magnetic shielding and chirality IV. The odd and even character of the shielding response to a chiral potential.

We investigate the odd and even character of the shielding response in a chiral molecule (modeled by a Ne8 helix) when subjected to a chiral potential. We establish that the diastereomeric splittings are a measure of odd powers of Vodd. Implications for diastereomeric, splittings of Xe in handed cages with handed tethers are discussed.