Importance of Nuclear Quantum Effects for Molecular Cocrystals with Short Hydrogen Bonds.

Many efforts have been recently devoted to the design and investigation of multicomponent pharmaceutical solids, such as salts and cocrystals. The experimental distinction between these solid forms is often challenging. Here, we show that the transformation of a salt into a cocrystal with a short hydrogen bond does not occur as a sharp phase transition but rather a smooth shift of the positional probability of the hydrogen atoms. A combination of solid-state NMR spectroscopy, X-ray diffraction, and diffuse reflectance measurements with density functional theory calculations that include nuclear quantum effects (NQEs) provides evidence of temperature-induced hydrogen atom shift in cocrystals with short hydrogen bonds. We demonstrate that for the predictions of the salt/cocrystal solid forms with short H-bonds, the computations have to include NQEs (particularly hydrogen nuclei delocalization) and temperature effects.

Mutations at hypothetical binding site 2 in insulin and insulin-like growth factors 1 and 2 result in receptor- and hormone-specific responses.

Information on how insulin and insulin-like growth factors 1 and 2 (IGF-1 and -2) activate insulin receptors (IR-A and -B) and the IGF-1 receptor (IGF-1R) is crucial for understanding the difference in the biological activities of these peptide hormones. Cryo-EM studies have revealed that insulin uses its binding sites 1 and 2 to interact with IR-A and have identified several critical residues in binding site 2. However, mutagenesis studies suggest that Ile-A10, Ser-A12, Leu-A13, and Glu-A17 also belong to insulin's site 2. Here, to resolve this discrepancy, we mutated these insulin residues and the equivalent residues in IGFs. Our findings revealed that equivalent mutations in the hormones can result in differential biological effects and that these effects can be receptor-specific. We noted that the insulin positions A10 and A17 are important for its binding to IR-A and IR-B and IGF-1R and that A13 is important only for IR-A and IR-B binding. The IGF-1/IGF-2 positions 51/50 and 54/53 did not appear to play critical roles in receptor binding, but mutations at IGF-1 position 58 and IGF-2 position 57 affected the binding. We propose that IGF-1 Glu-58 interacts with IGF-1R Arg-704 and belongs to IGF-1 site 1, a finding supported by the NMR structure of the less active Asp-58-IGF-1 variant. Computational analyses indicated that the aforementioned mutations can affect internal insulin dynamics and inhibit adoption of a receptor-bound conformation, important for binding to receptor site 1. We provide a molecular model and alternative hypotheses for how the mutated insulin residues affect activity.

Dimerization of Acetic Acid in the Gas Phase-NMR Experiments and Quantum-Chemical Calculations.

Due to the nature of the carboxylic group, acetic acid can serve as both a donor and acceptor of a hydrogen bond. Gaseous acetic acid is known to form cyclic dimers with two strong hydrogen bonds. However, trimeric and various oligomeric structures have also been hypothesized to exist in both the gas and liquid phases of acetic acid. In this work, a combination of gas-phase NMR experiments and advanced computational approaches were employed in order to validate the basic dimerization model of gaseous acetic acid. The gas-phase experiments performed in a glass tube revealed interactions of acetic acid with the glass surface. On the other hand, variable-temperature and variable-pressure NMR parameters obtained for acetic acid in a polymer insert provided thermodynamic parameters that were in excellent agreement with the MP2 (the second order Møller-Plesset perturbation theory) and CCSD(T) (coupled cluster with single, double and perturbative triple excitation) calculations based on the basic dimerization model. A slight disparity between the theoretical dimerization model and the experimental data was revealed only at low temperatures. This observation might indicate the presence of other, entropically disfavored, supramolecular structures at low temperatures.

A versatile insulin analog with high potency for both insulin and insulin-like growth factor 1 receptors: Structural implications for receptor binding.

Insulin and insulin-like growth factor 1 (IGF-1) are closely related hormones involved in the regulation of metabolism and growth. They elicit their functions through activation of tyrosine kinase-type receptors: insulin receptors (IR-A and IR-B) and IGF-1 receptor (IGF-1R). Despite similarity in primary and three-dimensional structures, insulin and IGF-1 bind the noncognate receptor with substantially reduced affinity. We prepared [d-HisB24, GlyB31, TyrB32]-insulin, which binds all three receptors with high affinity (251 or 338% binding affinity to IR-A respectively to IR-B relative to insulin and 12.4% binding affinity to IGF-1R relative to IGF-1). We prepared other modified insulins with the aim of explaining the versatility of [d-HisB24, GlyB31, TyrB32]-insulin. Through structural, activity, and kinetic studies of these insulin analogs, we concluded that the ability of [d-HisB24, GlyB31, TyrB32]-insulin to stimulate all three receptors is provided by structural changes caused by a reversed chirality at the B24 combined with the extension of the C terminus of the B chain by two extra residues. We assume that the structural changes allow the directing of the B chain C terminus to some extra interactions with the receptors. These unusual interactions lead to a decrease of dissociation rate from the IR and conversely enable easier association with IGF-1R. All of the structural changes were made at the hormones' Site 1, which is thought to interact with the Site 1 of the receptors. The results of the study suggest that merely modifications of Site 1 of the hormone are sufficient to change the receptor specificity of insulin.

Antifungal activity of analogues of antimicrobial peptides isolated from bee venoms against vulvovaginal Candida spp.

Candida albicans is the main causative agent of vulvovaginal candidiasis (VVC), a common mycosis in women, relapses of which are difficult to manage due to biofilm formation. This study aimed at developing novel non-toxic compounds active against Candida spp. biofilms. We synthesised analogues of natural antifungal peptides LL-III (LL-III/43) and HAL-2 (peptide VIII) originally isolated from bee venoms and elucidated their structures by nuclear magnetic resonance spectroscopy. The haemolytic, cytotoxic, antifungal and anti-biofilm activities of LL-III/43 and peptide VIII were then tested. LL-III/43 and VIII showed moderate cytotoxicity to HUVEC-2 cells and had comparable inhibitory activity against C. albicans and non-albicans spp. The lowest minimum inhibitory concentration (MIC90) of LL-III/43 was observed towards Candida tropicalis (0.8 µM). That was 8-fold lower than that of antimycotic amphotericin B. Both peptides can be used to inhibit Candida spp. bio film f ormation. Biofilm inhibitory concentrations (BIC50) ranged from 0.9 to 58.6 µM and biofilm eradication concentrations (BEC50) for almost all tested Candida spp. strains ranged from 12.8 to 200 µM. Als o pro ven were the peptides' abilities to reduce the area colonised by biofilms , inhibit hyphae formation and permeabilise cell membranes in biofil ms . LL-III/43 and VIII are promising candidates for further development as therapeutics against VVC.

Proton transfer in guanine-cytosine base pair analogues studied by NMR spectroscopy and PIMD simulations.

It has been hypothesised that proton tunnelling between paired nucleobases significantly enhances the formation of rare tautomeric forms and hence leads to errors in DNA replication. Here, we study nuclear quantum effects (NQEs) using deuterium isotope-induced changes of nitrogen NMR chemical shifts in a model base pair consisting of two tautomers of isocytosine, which form hydrogen-bonded dimers in the same way as the guanine-cytosine base pair. Isotope effects in NMR are consequences of NQEs, because ro-vibrational averaging of different isotopologues gives rise to different magnetic shielding of the nuclei. The experimental deuterium-induced chemical shift changes are compared with those calculated by a combination of path integral molecular dynamics (PIMD) simulations with DFT calculations of nuclear shielding. These calculations can directly link the observable isotope-induced shifts with NQEs. A comparison of the deuterium-induced changes of 15N chemical shifts with those predicted by PIMD simulations shows that inter-base proton transfer reactions do not take place in this system. We demonstrate, however, that NMR isotope shifts provide a unique possibility to study NQEs and to evaluate the accuracy of the computational methods used for modelling quantum effects in molecules. Calculations based on the PBE functional from the general-gradient-approximation family provided significantly worse predictions of deuterium isotope shifts than those with the hybrid B3LYP functional.

Influence of Intramolecular Charge Transfer and Nuclear Quantum Effects on Intramolecular Hydrogen Bonds in Azopyrimidines.

Intramolecular hydrogen bonds (IMHBs) in 5-azopyrimidines are investigated by NMR spectroscopy and DFT computations that involve nuclear quantum effects. A series of substituted 5-phenylazopyrimidines with one or two hydrogen bond donors able to form IMHBs with the azo group is prepared by azo coupling. The barrier of interconversion between two rotamers of the compounds with two possible IMHBs is determined by variable temperature NMR spectroscopy and it is demonstrated that the barrier is significantly affected by intramolecular charge transfer. Through-hydrogen-bond scalar coupling is investigated in 15N labeled compounds and the stability of the IMHBs is correlated with experimental NMR parameters and rationalized by path integral molecular dynamics simulations that involve nuclear quantum effects. Detailed information on the hydrogen bond geometry upon hydrogen-to-deuterium isotope exchange is obtained from a comparison of experimental and calculated NMR data.

Exploring Systematic Discrepancies in DFT Calculations of Chlorine Nuclear Quadrupole Couplings.

Previous studies have revealed significant discrepancies between density functional theory (DFT)-calculated and experimental nuclear quadrupolar coupling constants (CQ) for chlorine atoms, particularly in ionic solids. Various aspects of the computations are systematically investigated here, including the choice of the DFT functional, basis set convergence, and geometry optimization protocol. The effects of fast (fs) time-scale dynamics are probed using molecular dynamics (MD) and nuclear quantum effects (NQEs) are considered using path-integral MD calculations. It is shown that the functional choice is the most important factor related to improving the accuracy of the quadrupolar coupling calculations, and that functionals beyond the generalized gradient approximation (GGA) level, such as hybrid and meta-GGA functionals, are required for good correlations with experiment. The influence of molecular dynamics and NQEs is less important than the functional choice in the studied systems. A method which involves scaling the calculated quadrupolar coupling constant is proposed here; its application leads to good agreement with experimental data.

Paramagnetic encoding of molecules.

Contactless digital tags are increasingly penetrating into many areas of human activities. Digitalization of our environment requires an ever growing number of objects to be identified and tracked with machine-readable labels. Molecules offer immense potential to serve for this purpose, but our ability to write, read, and communicate molecular code with current technology remains limited. Here we show that magnetic patterns can be synthetically encoded into stable molecular scaffolds with paramagnetic lanthanide ions to write digital code into molecules and their mixtures. Owing to the directional character of magnetic susceptibility tensors, each sequence of lanthanides built into one molecule produces a unique magnetic outcome. Multiplexing of the encoded molecules provides a high number of codes that grows double-exponentially with the number of available paramagnetic ions. The codes are readable by nuclear magnetic resonance in the radiofrequency (RF) spectrum, analogously to the macroscopic technology of RF identification. A prototype molecular system capable of 16-bit (65,535 codes) encoding is presented. Future optimized systems can conceivably provide 64-bit (~10^19 codes) or higher encoding to cover the labelling needs in drug discovery, anti-counterfeiting and other areas.

Double Hydrogen Bonding Dimerization Propensity of Aqueous Hydroxy Acids Investigated Using Vibrational Optical Activity.

Lactic and malic acids are key substances in a number of biochemical processes in living cells and are also utilized in industry. Vibrational spectroscopy represents an efficient and sensitive way to study their structure and interactions. Since water is the natural environment, proper understanding of their molecular dynamics in aqueous solutions is of critical importance. To this end, we employed Raman spectroscopy and Raman optical activity (ROA) to study the conformation of l-lactic and l-malic acids in water (while varying pH, temperature, and concentration), with special emphasis on their double hydrogen bonding dimerization propensity. Raman and ROA experimental data were supported by extensive theoretical calculations of the vibrational properties and by additional experiments (IR absorption, vibrational circular dichroism, and NMR). Conformational behavior of the acids in water was described by molecular dynamics simulations. Reliability of the results was verified by calculating the vibrational properties of populated conformers and by comparing thus obtained spectral features with the experimental data. Calculations estimated the incidence of H-bonded dimers in water to be low in lactic acid and comparable to monomers in malic acid. The "hybrid" approach presented here reveals limitations of relying on the experimental spectra alone to study dimer formation.

Determination of nucleobase-pairing free energies from rotamer equilibria of 2-(methylamino)pyrimidines.

A fast straightforward method for the determination of free energies of modified nucleobase pairing is proposed. The method is based on the nuclear magnetic resonance (NMR) spectroscopy monitoring of conformational changes of 2-(methylamino)-pyrimidines upon intermolecular binding.