Understanding sterol-membrane interactions part I: Hartree-Fock versus DFT calculations of 13C and 1H NMR isotropic chemical shifts of sterols in solution and analysis of hydrogen-bonding effects

1H and 13C NMR chemical shifts are exquisitely sensitive probes of the local environment of the corresponding nuclei. Ultimately, direct determination of the chemical shifts of sterols in their membrane environment has the potential to reveal their molecular interactions and dynamics, in particular concerning the hydrogen-bonding partners of their OH groups. However, this strategy requires an accurate and efficient means to quantify the influence of the various interactions on chemical shielding. Herein the validity of Hartree-Fock and DFT calculations of the 13C and 1H NMR chemical shifts of cholesterol and ergosterol are compared with one another and with experimental chemical shifts measured in solution at 500 MHz. A computational strategy (definition of basis set, simpler molecular models for the sterols themselves and their molecular complexes) is proposed and compared with experimental data in solution. It is shown in particular that the effects of hydrogen bonding with various functional groups (water as a hydrogen-bond donor and acceptor, acetone) on NMR chemical shifts in CDCl3 solution can be accurately reproduced with this computational approach.     

Characterization of the monovalent ion position and hydrogen-bond network in guanine quartets by DFT calculations of NMR parameters

Conformational stability of G-quartets found in telomeric DNA quadruplex structures requires the coordination of monovalent ions. Here, an extensive Hartree-Fock and density functional theory analysis of the energetically favored position of Li+, Na+, and K+ ions is presented. The calculations show that at quartet-quartet distances observed in DNA quadruplex structures (3.3 A), the Li+ and Na+ ions favor positions of 0.55 and 0.95 A outside the plane of the G-quartet, respectively. The larger K+ ion prefers a central position between successive G-quartets. The energy barrier separating the minima in the quartet-ion-quartet model are much smaller for the Li+ and Na+ ions compared with the K+ ion; this suggests that K+ ions will not move as freely through the central channel of the DNA quadruplex. Spin-spin coupling constants and isotropic chemical shifts in G-quartets extracted from crystal structures of K+- and Na+-coordinated DNA quadruplexes were calculated with B3LYP/6-311G(d). The results show that the sizes of the trans-hydrogen-bond couplings are influenced primarily by the hydrogen bond geometry and only slightly by the presence of the ion. The calculations show that the R(N2N7) distance of the N2-H2...N7 hydrogen bond is characterized by strong correlations to both the chemical shifts of the donor group atoms and the (h2)J(N2N7) couplings. In contrast, weaker correlations between the (h3)J(N1C6') couplings and single geometric factors related to the N1-H1...O6=C6 hydrogen bond are observed. As such, deriving geometric information on the hydrogen bond through the use of trans-hydrogen-bond couplings and chemical shifts is more complex for the N1-H1...O6=C6 hydrogen bond than for the N2-H2...N7 moiety. The computed trans-hydrogen-bond couplings are shown to correlate with the experimentally determined couplings. However, the experimental values do not show such strong geometric dependencies.     

Uncovering Intramolecular π‐Type Hydrogen Bonds in Solution by NMR Spectroscopy and DFT Calculations

Reaction between the phosphinito bridged diplatinum species [(PHCy₂)Pt(μ‐PCy₂){κ²P,O‐μ‐P(O)Cy₂}Pt(PHCy₂)](Pt–Pt) ( 1 ), and (trimethylsilyl)acetylene at 273 K affords the σ‐acetylide complex [(PHCy₂)(η¹‐Me₃SiC≡C)Pt(μ‐PCy₂)Pt(PHCy₂){κP‐P(OH)Cy₂}](Pt–Pt) ( 2 ) featuring an intramolecular π‐type hydrogen bond. Scalar and dipolar couplings involving the POH proton were detected by 2D NMR experiments. Relativistic DFT calculations of the geometry, relative energy, and NMR properties of model systems of 2 confirmed the structural assignment and allowed the energy of the π‐type hydrogen bond to be estimated (ca. 22 kJ mol^?1).     

Self-Assembly of DNA and RNA Building Blocks Explored by Nitrogen-14 NMR Crystallography: Structure and Dynamics

The isotopic enrichment of nucleic acids with nitrogen-15 is often carried out by solid-phase synthesis of oligonucleotides using phosphoramidite precursors that are synthetically demanding and expensive. These synthetic challenges, combined with the overlap of chemical shifts, explain the lag of nitrogen-15 NMR studies of nucleic acids behind those of proteins. For the structural characterization of DNA and RNA-related systems, new NMR methods that exploit the naturally occurring 99.9 % abundant nitrogen-14 isotope are therefore highly desirable. In this study, we have investigated nitrogen-14 spectra of self-assembled quartets based on the nucleobase guanine in the solid state by means of magic-angle spinning NMR spectroscopy. The network of dipolar proton–nitrogen couplings between neighboring stacked purine units is probed by 2D spectra based on ¹H→¹⁴N→¹H double cross-polarization. Interplane dipolar contacts are identified between the stacked G quartets. The assignment is supported by density functional theory (DFT) calculations of the anisotropic chemical shifts and quadrupolar parameters. The experimental spectra are fully consistent with internuclear distances obtained in silico. Averaging of chemical shifts due to internal motions can be interpreted by semiempirical calculations. This method can easily be extended to synthetic G quartets based on nucleobase or nucleoside analogs and potentially to oligonucleotides.     

Torsional and Electronic Factors Control the C−H⋅⋅⋅O Interaction

The precise role of non‐conventional hydrogen bonds such as the C?H???O interaction in influencing the conformation of small molecules remains unresolved. Here we survey a series of β‐turn mimetics using X‐ray crystallography and NMR spectroscopy in conjunction with quantum calculation, and conclude that favourable torsional and electronic effects are important for the population of states with conformationally influential C?H???O interactions. Our results also highlight the challenge in attempting to deconvolute a myriad of interdependent noncovalent interactions in order to focus on the contribution of a single one. Within a small molecule that is designed to resemble the complexity of the environment within peptides and proteins, the interplay of different steric burdens, hydrogen‐acceptor/‐donor properties and rotational profiles illustrate why unambiguous conclusions based solely on NMR chemical shift data are extremely challenging to rationalize.     

Determining the geometry of hydrogen bonds in solids with picometer accuracy by quantum-chemical calculations and NMR spectroscopy.

The structure of multiply hydrogen-bonded systems is determined with picometer accuracy by a combined solid-state NMR and quantum-chemical approach. On the experimental side, advanced 1H-15N dipolar recoupling NMR techniques are capable of providing proton-nitrogen distances of up to about 250 pm with an accuracy level of +/-1 pm for short distances (i.e., around 100 pm) and +/-5 pm for longer ones (i.e., 180 to 250 pm). The experiments were performed under fast magic-angle spinning, which ensures sufficient dipolar decoupling and spectral resolution of the 1H resonance lines. On the quantum-chemical side, the structures of the hydrogen-bonded systems were computationally optimised, yielding complete sets of nitrogen-proton and proton-proton distances, which are essential for correctly interpreting the experimental NMR data. In this way, nitrogen-proton distances were determined with picometer accuracy, so that vibrational averaging effects on dipole-dipole couplings need to be considered. The obtained structures were finally confirmed by the complete agreement of computed and experimental 'H and '5N chemical shifts. This demonstrates that solid-state NMR and quantum-chemical methods ideally complement each other and, in a combined manner, represent a powerful approach for reliable, high-precision structure determination whenever scattering techniques are inapplicable.     

Highly Enantioselective Extraction of Underivatized Amino Acids by the Uryl‐Pendant Hydroxyphenyl‐Binol Ketone

The hydroxyphenyl chiral ketone, (S)‐ 3 , reacts with D ‐amino acids bearing hydrophobic side chains exclusively over the L ‐amino acids in a two‐phase liquid–liquid extraction, and thus acts as a highly stereoselective extractant. Calculations for the energy‐minimized structures for the imine diastereomers and the comparison of the selectivities with other phenyl ketones, (S)‐ 4 and (S)‐ 5 , demonstrate that the hydrogen bond between the carboxylate group and the phenolic hydroxyl group contributes to the remarkable enantioselectivities. The multiple hydrogen bonds present in the imine of (S)‐ 3 reinforce the rigidity, and results in the difference between the stabilities of the imine diastereomers. The imine could be hydrolyzed in methanolic HCl solution, and the extraction of the evaporated residues revived the organic layer of (S)‐ 3 , which could enter into a new extractive cycle and leaves the D ‐amino acid with enantiomeric excess (ee) values of over 97 % in the aqueous layer.     

Spin–Spin Coupling Constants Transmitted through Ir?H⋅⋅⋅H?N Dihydrogen Bonds

Relativity matters: Calculations of NMR shielding tensors and spin–spin coupling constants transmitted through Ir? H???H? N dihydrogen bonds are presented. The picture shows one of the simplified models employed. It is shown that the spin–orbit relativistic effects influence the NMR shielding constants far more than the spin–spin coupling constants.

     

Asymmetry induction by cooperative intermolecular hydrogen bonds in surface-anchored layers of achiral molecules.

The mesoscale induction of two-dimensional supramolecular chirality (formation of 2D organic domains with a single handedness) was achieved by self-assembly of 1,2,4-benzenetricarboxylic (trimellitic) acid on a Cu(100) surface at elevated temperatures. The combination of spectroscopic [X-ray photoelectron spectroscopy (XPS) and near-edge X-ray absorption fine structure (NEXAFS)], real-space-probe [scanning tunneling microscopy (STM)], and computational [density functional theory (DFT)] methods allows a comprehensive characterization of the obtained organic adlayers, where details of molecular adsorption geometry, intermolecular coupling, and surface chemical bonding are elucidated. The trimellitic acid species, comprising three functional carboxylic groups, form distinct stable mirror-symmetric hydrogen-bonded domains. The chiral ordering is associated with conformational restriction in the domains: molecules anchor to the substrate with an ortho carboxylate group, providing two para carboxylic acid moieties for collective lateral interweaving through H bonding, which induces a specific tilt of the molecular plane. The ease of molecular symmetry switching in domain formation makes homochiral-signature propagation solely limited by the terrace width. The molecular layer modifies the morphology of the underlying copper substrate and induces mum-sized strictly homochiral terraces.     

Modeling of heavy-atom-ligand NMR spin-spin coupling in solution: molecular dynamics study and natural bond orbital analysis of Hg-C coupling constants

Ab initio molecular dynamics (MD) and relativistic density functional NMR methods were applied to calculate the one‐bond Hg? C NMR indirect nuclear spin–spin coupling constants (J) of [Hg(CN)₂] and [CH₃HgCl] in solution. The MD averages were obtained as J(¹⁹⁹Hg? ¹³C)=3200 and 1575 Hz, respectively. The experimental Hg? C spin–spin coupling constants of [Hg(CN)₂] in methanol and [CH₃HgCl] in DMSO are 3143 and 1674 Hz, respectively. To deal with solvent effects in the calculations, finite “droplet” models of the two systems were set up. Solvent effects in both systems lead to a strong increase of the Hg? C coupling constant. From a relativistic natural localized molecular orbital (NLMO) analysis, it was found that the degree of delocalization of the Hg 5d_σ nonbonding orbital and of the Hg? C bonding orbital between the two coupled atoms, the nature of the trans Hg? C/Cl bonding orbital, and the s character of these orbitals, exhibit trends upon solvation of the complexes that, when combined, lead to the strong increase of J(Hg? C).     

Discriminating Halogen‐Bonding from Other Noncovalent Interactions by a Combined NOE NMR/DFT Approach

Herein a combined NOE NMR/DFT methodology to discriminate between adducts held together by halogen bonding (XB) and other noncovalent interactions (non‐XB, such as lone pair/π), based on the determination of the XB donors′ and acceptors′ relative orientation, is proposed. In particular, ¹⁹F,¹H HOESY NMR spectroscopy experiments and DFT calculations on different XB donors, such as perfluorohexyl iodide ( I1 ), iodopentafluorobenzene ( I2 ) and bromopentafluorobenzene ( Br ), combined with different Lewis bases, such as 1,4‐diazabicyclo[2.2.2]octane ( DABCO ) and 2,4,6‐trimethylpyridine ( Me₃Py ), were performed. The results clearly show that in the case DABCO / I1 the XB adduct is practically the only one present in solution, whereas for the other pairs a certain amount of non‐XB adduct is present. Combining DFT and HOESY results, the amount of non‐XB adducts can be roughly quantified under our experimental conditions as 4 % for DABCO / I2 , between 10 and 20 % for Me₃Py / I1 and Me₃Py / I2 , and 44 % for DABCO / Br.     

First investigation of non-classical dihydrogen bonding between an early transition-metal hydride and alcohols: IR, NMR, and DFT approach

The interaction of [NbCp(2)H(3)] with fluorinated alcohols to give dihydrogen-bonded complexes was studied by a combination of IR, NMR and DFT methods. IR spectra were examined in the range from 200-295 K, affording a clear picture of dihydrogen-bond formation when [NbCp(2)H(3)]/HOR(f) mixtures (HOR(f) = hexafluoroisopropanol (HFIP) or perfluoro-tert-butanol (PFTB)) were quickly cooled to 200 K. Through examination of the OH region, the dihydrogen-bond energetics were determined to be 4.5+/-0.3 kcal mol(-1) for TFE (TFE = trifluoroethanol) and 5.7+/-0.3 kcal mol(-1) for HFIP. (1)H NMR studies of solutions of [NbCp(2)H(2)(B)H(A)] and HFIP in [D(8)]toluene revealed high-field shifts of the hydrides H(A) and H(B), characteristic of dihydrogen-bond formation, upon addition of alcohol. The magnitude of signal shifts and T(1) relaxation time measurements show preferential coordination of the alcohol to the central hydride H(A), but are also consistent with a bifurcated character of the dihydrogen bonding. Estimations of hydride-proton distances based on T(1) data are in good accord with the results of DFT calculations. DFT calculations for the interaction of [NbCp(2)H(3)] with a series of non-fluorinated (MeOH, CH(3)COOH) and fluorinated (CF(3)OH, TFE, HFIP, PFTB and CF(3)COOH) proton donors of different strengths showed dihydrogen-bond formation, with binding energies ranging from -5.7 to -12.3 kcal mol(-1), depending on the proton donor strength. Coordination of proton donors occurs both to the central and to the lateral hydrides of [NbCp(2)H(3)], the former interaction being of bifurcated type and energetically slightly more favourable. In the case of the strong acid H(3)O(+), the proton transfer occurs without any barrier, and no dihydrogen-bonded intermediates are found. Proton transfer to [NbCp(2)H(3)] gives bis(dihydrogen) [NbCp(2)(eta(2)-H(2))(2)](+) and dihydride(dihydrogen) complexes [NbCp(2)(H)(2)(eta(2)-H(2))](+) (with lateral hydrides and central dihydrogen), the former product being slightly more stable. When two molecules of TFA were included in the calculations, in addition to the dihydrogen-bonded adduct, an ionic pair formed by the cationic bis(dihydrogen) complex [NbCp(2)(eta(2)-H(2))(2)](+) and the homoconjugated anion pair (CF(3)COO...H...OOCCF(3))(-) was found as a minimum. It is very likely that these ionic pairs may be intermediates in the H/D exchange between the hydride ligands and the OD group observed with the more acidic alcohols in the NMR studies.     

Theoretical analysis of the terahertz spectrum of the high explosive PETN.

The experimental solid-state terahertz (THz) spectrum (3 to 120 cm(-1)) of the high explosive pentaerythritol tetranitrate (PETN, C(5)H(6)N(4)O(12)) has been modeled using solid-state density functional theory (DFT) calculations. Solid-state DFT, employing the BP density functional, is in best qualitative agreement with the features in the previously reported THz spectrum. The crystal environment of PETN includes numerous intermolecular hydrogen-bonding interactions that contribute to large (up to 80 cm(-1)) calculated shifts in molecular normal-mode positions in the solid state. Comparison of the isolated-molecule and solid-state normal-mode calculations for a series of density functionals reveals the extent to which the inclusion of crystal-packing interactions and the relative motions between molecules are required for correctly reproducing the vibrational structure of solid-state THz spectra. The THz structure below 120 cm(-1) is a combination of both intermolecular (relative rotations and translations) and intramolecular (torsions, large amplitude motions) vibrational motions. Vibrational-mode analyses indicate that the first major feature (67.2 cm(-1)) in the PETN THz spectrum contains all of the optical rotational and translational cell modes and no internal (molecular) vibrational modes.     

Metallobiliverdin Radicals—DFT Studies

Several aspects of the molecular and electronic structure of biliverdin derivatives have been studied using density functional theory (DFT). The calculations have been performed for complexes of trianion (BvO₂)^3? and dianion [BvO(OH)]^2?, derived from two tautomeric forms of biliverdin, BvO₂H₃ and [BvO(OH)]H₂, with redox innocent metal ions: lithium(I ), zinc(II ), and gallium(III ). One‐electron‐oxidized and ‐reduced forms of each complex (cation and anion radicals) have been also considered. The molecular structures of all species investigated are characterized by a helical arrangement of tetrapyrrolic ligands with the metal ion lying in the plane formed by the two central pyrrole rings. The spin density distribution in four types of metallobiliverdin radicals—[(BvO₂^.)Mⁿ⁺]^n‐2, [{BvO(OH)^.}Mⁿ⁺]^n‐1 (cation radicals), [(BvO2^.)Mⁿ⁺]^n‐4, [{BvO(OH)^.}Mⁿ⁺]^n‐3 (anion radicals)—has been investigated. In general, the absolute values of spin density on meso carbon atoms were larger than for the β‐carbon atoms. Sign alteration of spin density has been found for meso positions, and also for the β‐carbon atoms of at least two pyrrole rings. The calculated spin density maps accounted for the essential NMR spectroscopic features of iron biliverdin derivatives, including the considerable isotropic shifts detected for the meso resonances and shift alteration at the meso and β‐positions.     

The calculation of indirect nuclear spin-spin coupling constants in large molecules

We present calculations of indirect nuclear spin-spin coupling constants in large molecular systems, performed using density functional theory. Such calculations, which have become possible because of the use of linear-scaling techniques in the evaluation of the Coulomb and exchange-correlation contributions to the electronic energy, allow us to study indirect spin-spin couplings in molecules of biological interest, without having to construct artificial model systems. In addition to presenting a statistical analysis of the large number of short-range coupling constants in large molecular systems, we analyse the asymptotic dependence of the indirect nuclear spin-spin coupling constants on the internuclear separation. In particular, we demonstrate that, in a sufficiently large one-electron basis set, the indirect spin-spin coupling constants become proportional to the inverse cube of the internuclear separation, even though the diamagnetic and paramagnetic spin-orbit contributions to the spin-spin coupling constants separately decay as the inverse square of this separation. By contrast, the triplet Fermi contact and spin-dipole contributions to the indirect spin-spin coupling constants decay exponentially and as the inverse cube of the internuclear separation, respectively. Thus, whereas short-range indirect spin-spin coupling constants are usually dominated by the Fermi contact contribution, long-range coupling constants are always dominated by the negative diamagnetic spin-orbit contribution and by the positive paramagnetic spin-orbit contribution, with small spin-dipole and negligible Fermi contact contributions.