The accuracy limit of chemical shift predictions for species in aqueous solution.

Interpreting NMR experiments benefits from first-principles predictions of chemical shifts. Reaching the accuracy limit of theory is relevant for unambiguous structural analysis and dissecting theoretical approximations. Since accurate chemical shift measurements are based on using internal reference compounds such as trimethylsilylpropanesulfonate (DSS), a detailed comparison of experimental with theoretical data requires simultaneous consideration of both target and reference species ensembles in the same solvent environment. Here we show that ab initio molecular dynamics simulations to generate liquid-state ensembles of target and reference compounds, including explicitly their short-range solvation environments and combined with quantum-mechanical solvation models, allows for predicting highly accurate 1H (∼0.1-0.5 ppm) and aliphatic 13C (∼1.5 ppm) chemical shifts for aqueous solutions of the model compounds trimethylamine N-oxide (TMAO) and N-methylacetamide (NMA), referenced to DSS without any system-specific adjustments. This encompasses the two peptide bond conformations of NMA identified by NMR. The results are used to derive a general-purpose guideline set for predictive NMR chemical shift calculations of NMA in the liquid state and to identify artifacts of force field models. Accurate predictions are only obtained if a sufficient number of explicit water molecules is included in the quantum-mechanical calculations, disproving a purely electrostatic model of the solvent effect on chemical shifts.

Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory.

Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.

Structural Sampling and Solvation Models for the Simulation of Electronic Spectra: Pyrazine as a Case Study.

The impact of sampling methods on spectral broadening in the gas phase and on the convergence of spectra in aqueous solution when using microsolvation, continuum solvation, and hybrid models is studied using pyrazine as a test case. For the sake of comparing classical Maxwell-Boltzmann and Wigner samplings in the gas phase, static and time-resolved X-ray absorption spectra after photoexcitation to the lowest 1B2u(ππ*) state, as well as the static UV-vis absorption spectrum, are considered. In addition, the UV-vis absorption spectrum of pyrazine in aqueous solution is also computed in order to systematically investigate its convergence with the number of explicitly included solvent shells with and without taking bulk solvation effects into account with the conductor-like screening model to represent implicit water beyond such explicit solute complexes. Concerning the static and time-resolved X-ray absorption spectra of pyrazine at the carbon K-edge as well as its UV-vis absorption spectrum in the gas phase, we find that these spectra obtained with Wigner and Maxwell-Boltzmann samplings substantially agree. For the UV-vis absorption spectrum in the aqueous solution, only the first two energetically low-lying bands converge quickly with the size of the explicitly included solvation shells, either with or without an additional continuum solvation taken into account. In stark contrast, calculations of the higher-lying excitations relying on finite microsolvated clusters without additional continuum solvation severely suffer from unphysical charge-transfer excitations into Rydberg-like orbitals at the cluster/vacuum interface. This finding indicates that computational UV-vis absorption spectra covering sufficiently high-lying states converge only if continuum solvation of the explicitly microsolvated solutes is included in the models.

A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water.

The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of "vertically desolvated" spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run.

Solvent Effects in the Ultraviolet and X-ray Absorption Spectra of Pyridazine in Aqueous Solution.

Electrostatic interaction of the solvent with the solute and fluctuations of the solvent configurations may make excitation energies of the solute different from those in the gas phase. These effects may dominate photoinduced or chemical reaction dynamics in solution systems and can be observed as shifts or broadening of peaks in absorption spectra. In this work, the nitrogen K-edge X-ray absorption spectra were measured for pyridazine in the gas phase and in aqueous solution. The ultraviolet and X-ray absorption spectra of pyridazine in aqueous solution, as well as those in the gas phase, were then calculated with models based on the algebraic-diagrammatic construction through second order [ADC(2)] with the resolution-of-identity (RI) approximation and compared with the spectra obtained in experiments. For aqueous solution, explicit local solvation structures were extracted from an ab initio molecular dynamics (AIMD) trajectory of pyridazine in bulk water, and RI-ADC(2) was combined with the conductor-like screening model (COSMO). The experimental absorption spectra of pyridazine in aqueous solution were reproduced with good accuracy by theoretical treatment of an ensemble containing the explicit local solvation structures of pyridazine with relevant water molecules combined with the COSMO solvation model of water for long-range solvation.

Dynamics of Water in the Solvation Shell of an Iodate Ion: A Born-Oppenheimer Molecular Dynamics Study.

The iodate ion has an anisotropic structure and charge distribution. It has a pyramidal shape with the iodine atom located at the peak of the pyramid. The water molecules interact differently with the positively charged iodine and the negatively charged oxygen atoms of this anion, giving rise to two distinct solvation shells. In the present study, we have performed ab initio Born-Oppenheimer molecular dynamics simulations to investigate the dynamics of water molecules in the iodine and oxygen solvation shells of the iodate ion and compared the behavior with those of the bulk. The dynamics of water is calculated for both the BLYP and the dispersion-corrected BLYP-D3 functionals at room temperature. The dynamics of water in the solvation shells at higher temperatures of 353 and 330 K has also been investigated for the BLYP and BLYP-D3 functionals, respectively. The hydrogen bond dynamics, vibrational spectral diffusion, orientational and translational diffusion, and residence dynamics of water molecules in the two solvation shells are looked at in the current study. The ion-water hydrogen bond dynamics is found to be somewhat faster than that for water-water hydrogen bonds in the bulk, which can be attributed to a ring-like electron distribution on the iodate oxygens. The dynamical trends are connected to the water structure making/breaking properties of the positively charged iodine and negatively charged oxygen sites of the anion. Furthermore, orientational jumps of the iodate ion and also those of surrounding water molecules which are hydrogen bonded to the oxygen atoms of the iodate ion are also investigated. It is found that the nature of these orientational jumps can be different from those reported earlier for planar polyoxyanions such as the nitrate ion.

Caspase-3 mediated programmed cell death by a gold-stabilised peptide carbene.

The synthesis of novel N-heterocyclic carbene complexes derived from a tripeptide ligand (L), containing non-natural amino acid, thiazolylalanine is described here. The peptide ligand was reacted with suitable precursors to generate gold and mercury carbene complexes. The plausible structures of both complexes were predicted by spectroscopic data and DFT calculations. The binding energy data was also analyzed to predict their stability. The gold carbene complex (1A), showed activity against MCF7 breast cancer cell line due to mitochondrial triggered caspase-3 mediated programmed cell death. Its internalization inside cells could be observed due to autofluorescence. This study affords a methodology for successful generation of peptide carbene complexes for their therapeutic potential.

Conformation-Induced Dynamical Heterogeneity of Water in the Solvation Shell of Zwitterionic γ-Aminobutyric Acid.

The structure and dynamics of water molecules around the carboxylate and amino groups of γ-aminobutyric acid (GABA), a primary neurotransmitter in mammals, are investigated by means of ab initio molecular dynamics simulation. Zwitterionic GABA has two major conformations in water, namely, the open and the closed conformations. The angle-averaged one-dimensional structures of water in the solvation shells around the carboxylate and amino groups are found to be quite similar for the closed and open conformations of the solute. The two-dimensional structural correlations, which describe the solvation shell structure with better resolution, reveal some differences in the arrangement of water molecules around the solute for its open and closed conformations. It is found that the dynamics of solvation shells in the two conformations vary only slightly. However, the existence of trapped water between the oppositely charged carboxylate and amino groups of GABA in its closed form is found to give rise to a very different dynamical behavior as compared to the overall solvation shell in the same conformation as well as that in the open conformation. Thus, dynamical heterogeneity at a local level is induced by a change in the conformation of zwitterionic GABA. Such trapped water is not seen in the open form of the solute. A similar type of "connecting water" has also been observed for microsolvated ß-alanine in a recent experimental and theoretical study ( Ghassemizadeh J. Phys. Chem. B , 2019 , 123 , 4392 - 4399 ). Thus, the current study shows the variation of binding properties of water with the change in the conformation of zwitterionic GABA, which in turn changes the dynamics of water at a local level. The conformation-induced changes in the water dynamics constitute an example of dynamical heterogeneity of water, which is normally observed in the presence of large biomolecules like proteins, DNA, etc.

On the issue of closed versus open forms of gamma-aminobutyric acid (GABA) in water: Ab initio molecular dynamics and metadynamics studies.

Gamma-aminobutyric acid (GABA), a primary neurotransmitter, accomplishes its activities by binding to different receptor sites in different conformations. It is known to have two major conformers: the closed and open forms. Earlier studies on preferred conformation of GABA in water revealed differing results with some reporting the open form while others inferring the closed form to be more stable. We found the existence of many open forms and only one closed form of GABA in water through ab initio metadynamics simulation. Some of the open conformers are equally or more stable while others are less stable than the closed form. Free energy barriers reveal that different conformers are interconvertible at room temperature in typical experimental time scales. Ab initio molecular dynamics simulations are performed to further investigate the inter-conversion of various conformers of GABA in water and their dipole moments and also to make connections to experiments on the conformation of GABA in water.

Nature of hydration shells of a polyoxy-anion with a large cationic centre: The case of iodate ion in water.

The structural nature of the solvation shells of an iodate ion, which is known to be a polyoxy-anion with a large cationic centre, is investigated by means of Born-Oppenheimer molecular dynamics (BOMD) simulations using BLYP and the dispersion corrected BLYP-D3 functionals. The iodate ion is found to have two distinct solvation regions around the positively charged iodine (iodine solvation shell or ISS) and the negatively charged oxygens (oxygen solvation shell or OSS). We have looked at the spatial, orientational, and hydrogen bond distributions of water in the two solvation regions. It is found that the water orientational profile in the ISS is typical of a cation hydration shell. The hydrogen bonded structure of water in the OSS is found to be very similar to that of the bulk water structure. Thus, the iodate ion essentially behaves like a positively charged iodine ion in water as if there is no anionic part. This explains why the cationic character of the iodate ion was prominently seen in earlier studies. The arrangement of water molecules in the two solvation shells and in the intervening regions around the iodate ion is further resolved by looking at structural cross-correlations. The electronic properties of the solvation shells are also looked at by calculating the solute-solvent orbital overlap and dipole moments of the solute and solvation shell water. We have also performed BOMD simulations of iodate ion-water clusters at experimentally relevant conditions. The simulation results are found to be in agreement with experimental results. © 2018 Wiley Periodicals, Inc.

Born-Oppenheimer Molecular Dynamics Simulations of a Bromate Ion in Water Reveal Its Dual Kosmotropic and Chaotropic Behavior.

The solvation structure and dynamics of a bromate (BrO3-) ion in water are studied by means of Born-Oppenheimer molecular dynamics simulations at two different temperatures using the Becke-Lee-Yang-Parr functional with Grimme D3 dispersion corrections. The bromate ion possesses a pyramidal structure, and it has two types of solvation sites, namely, the bromine and oxygen atoms. We have looked at different radial and orientational distributions of water molecules around the bromate ion and also investigated their hydrogen bonding properties. The solvation structure of the bromate ion is also compared with that of the iodate (IO3-) ion, which is structurally rather similar to the bromate ion and was found to have some unusual solvation properties in water. It is found that the bromate ion follows a similar trend as that followed by the iodate ion as far as the solvation structure is concerned. However, the effect of the former on surrounding water is found to be much weaker than that of the latter. On the dynamical side, we have looked at diffusion, residence dynamics, and also the orientational and hydrogen bond relaxation of water molecules around the BrO3- ion and compared them with those of the bulk. Dynamical results are presented for both H2O and D2O around the BrO3- ion. Interpretation of the dynamical results in terms of structure-making (kosmotropic)/-breaking (chaotropic) properties of the BrO3- ion reveals that the bromine atom of this ion acts as a water structure breaker, whereas the three oxygens act as water structure makers. Thus, in spite of being a single ion, the bromate ion has dual characteristics and the experimentally observed kosmotropic ability of this ion is actually a trade-off between a chaotropic site (the bromine atom) and three kosmotropic sites (three oxygen atoms) that are present in the ion.

Ab Initio Molecular Dynamics Simulation of the Phosphate Ion in Water: Insights into Solvation Shell Structure, Dynamics, and Kosmotropic Activity.

The structure and dynamics of solvation shells of the phosphate ion in deuterated water are studied by means of Born-Oppenheimer molecular dynamics simulation. The total number of molecules in the first and second solvation shells is found to be close to the effective hydration number reported experimentally. The OD bonds that are hydrogen bonded to the phosphate ion are found to be red shifted as compared to bulk water, which is consistent with experimental results. However, the two OD bonds of the same water molecule in the first hydration shell are found to be vibrationally distinct, which can be attributed to different strengths of the ion-water and water-water hydrogen bonds near the ion. Also, the hydrogen bonds formed by the second solvation shell OD bonds are somewhat stronger than the bulk. This finding shows a long ranged effect of the phosphate ion on water and also gives insights into the water structure making property of this anion. The dynamics of water in the first solvation shell is found to be significantly slower than that of the bulk. We have investigated the origin of the orientational slowing down of the first solvation shell water molecules and made connections to similar results observed experimentally.

Fluorescent Biotin Analogues for Microstructure Patterning and Selective Protein Immobilization.

Benzyl substitution on ureido nitrogens of biotin led to manifestation of aggregation-induced emission, which was studied by steady-state fluorescence, microscopy, and TD-DFT, providing a rationale into the observed photophysical behavior. Besides exhibiting solvatochromism, the biotin derivatives revealed emission peaks centered at ∼430 and 545 nm, which has been attributed to the π-π stacking interactions. Our TD-DFT results also correlate the spectroscopic data and quantify the nature of transitions involved. The isothermal titration calorimetry data substantiates that the binding of the biotin derivatives with avidin are pretty strong. These derivatives on lithographic patterning present a platform for site specific strept(avidin) immobilization, thus opening avenues for potential applications exploiting these interactions. The fluorescent biotin derivatives can thus find applications in cellular biology and imaging.