Towards Quantum-Chemical Modeling of the Activity of Anesthetic Compounds.

The modeling of the activity of anesthetics is a real challenge because of their unique electronic and structural characteristics. Microscopic approaches relevant to the typical features of these systems have been developed based on the advancements in the theory of intermolecular interactions. By stressing the quantum chemical point of view, here, we review the advances in the field highlighting differences and similarities among the chemicals within this group. The binding of the anesthetics to their partners has been analyzed by Symmetry-Adapted Perturbation Theory to provide insight into the nature of the interaction and the modeling of the adducts/complexes allows us to rationalize their anesthetic properties. A new approach in the frame of microtubule concept and the importance of lipid rafts and channels in membranes is also discussed.

Affinity of fentanyl and its derivatives for the σ1-receptor.

Fentanyl and its 11 commercially available derivatives were investigated as to their affinity for the σ1 receptor. The parent compound is a rather poor binder (IC50 = 4973 nM), but its close derivatives (benzylfentanyl or p-fluorofentanyl) have submicromolar affinities. Modelling provides a structural basis for the observed trends in activity.

Molecular dynamics of fentanyl bound to µ-opioid receptor.

The molecular dynamics simulations of fentanyl complexed with the µ-opioid receptor (µOR) were studied using both inactive 4DKL and active 5C1M opioid receptor crystal structures. Analogous simulations in morphine with or without a ligand were done for comparison. Simulations of the inactive states were carried out in the absence and presence of the Na+ ion. The obtained fentanyl's binding mode agrees with some of the mutagenesis data, and it overlaps with that of morphine only to a minor extent. Notably, fentanyl stabilizes different rotameric states of Trp2936.48 than observed for morphine or unliganded receptor. Another difference is tighter arrangement of the interaction between Asp1473.32 and Tyr3267.43 (a link between helices TM3 and TM7) in the presence of fentanyl. Principal component analysis reveals differences in the trajectories dependent on the ligand bound. The differences found could be linked to ligand-dependent efficacy with respect to receptor intracellular signaling events.

Fentanyl Family at the Mu-Opioid Receptor: Uniform Assessment of Binding and Computational Analysis.

Interactions of 21 fentanyl derivatives with µ-opioid receptor (µOR) were studied using experimental and theoretical methods. Their binding to µOR was assessed with radioligand competitive binding assay. A uniform set of binding affinity data contains values for two novel and one previously uncharacterized derivative. The data confirms trends known so far and thanks to their uniformity, they facilitate further comparisons. In order to provide structural hypotheses explaining the experimental affinities, the complexes of the studied derivatives with µOR were modeled and subject to molecular dynamics simulations. Five common General Features (GFs) of fentanyls' binding modes stemmed from these simulations. They include: GF1) the ionic interaction between D147 and the ligands' piperidine NH⁺ moiety; GF2) the N-chain orientation towards the µOR interior; GF3) the other pole of ligands is directed towards the receptor outlet; GF4) the aromatic anilide ring penetrates the subpocket formed by TM3, TM4, ECL1 and ECL2; GF5) the 4-axial substituent (if present) is directed towards W318. Except for the ionic interaction with D147, the majority of fentanyl-µOR contacts is hydrophobic. Interestingly, it was possible to find nonlinear relationships between the binding affinity and the volume of the N-chain and/or anilide's aromatic ring. This kind of relationships is consistent with the apolar character of interactions involved in ligand⁻receptor binding. The affinity reaches the optimum for medium size while it decreases for both large and small substituents. Additionally, a linear correlation between the volumes and the average dihedral angles of W293 and W133 was revealed by the molecular dynamics study. This seems particularly important, as the W293 residue is involved in the activation processes. Further, the Y326 (OH) and D147 (Cγ) distance found in the simulations also depends on the ligands' size. In contrast, neither RMSF measures nor D114/Y336 hydrations show significant structure-based correlations. They also do not differentiate studied fentanyl derivatives. Eventually, none of 14 popular scoring functions yielded a significant correlation between the predicted and observed affinity data (R < 0.30, n = 28).

Molecular dynamics simulations suggest why the A2058G mutation in 23S RNA results in bacterial resistance against clindamycin.

Clindamycin, a lincosamide antibiotic, binds to 23S ribosomal RNA and inhibits protein synthesis. The A2058G mutation in 23S RNA results in bacterial resistance to clindamycin. To understand the influence of this mutation on short-range interactions of clindamycin with 23S RNA, we carried out full-atom molecular dynamics simulations of a ribosome fragment containing clindamycin binding site. We compared the dynamical behavior of this fragment simulated with and without the A2058G mutation. Molecular dynamics simulations suggest that clindamycin in the native ribosomal binding site is more internally flexible than in the A2058G mutant. Only in the native ribosome fragment did we observe intramolecular conformational change of clindamycin around its C7-N1-C10-C11 dihedral. In the mutant, G2058 makes more stable hydrogen bonds with clindamycin hindering its conformational freedom in the ribosome-bound state. Clindamycin binding site is located in the entrance to the tunnel through which the newly synthesized polypeptide leaves the ribosome. We observed that in the native ribosome fragment, clindamycin blocks the passage in the tunnel entrance, whereas in the mutated fragment the aperture is undisturbed due to a different mode of binding of clindamycin in the mutant. Restricted conformational freedom of clindamycin in a position not blocking the tunnel entrance in the A2058G mutant could explain the molecular mechanism of bacterial resistance against clindamycin occurring in this mutant.

Structural Insights into σ1 Receptor Interactions with Opioid Ligands by Molecular Dynamics Simulations.

Despite considerable advances over the past years in understanding the mechanisms of action and the role of the σ1 receptor, several questions regarding this receptor remain unanswered. This receptor has been identified as a useful target for the treatment of a diverse range of diseases, from various central nervous system disorders to cancer. The recently solved issue of the crystal structure of the σ1 receptor has made elucidating the structure-activity relationship feasible. The interaction of seven representative opioid ligands with the crystal structure of the σ1 receptor (PDB ID: 5HK1) was simulated for the first time using molecular dynamics (MD). Analysis of the MD trajectories has provided the receptor-ligand interaction fingerprints, combining information on the crucial receptor residues and frequency of the residue-ligand contacts. The contact frequencies and the contact maps suggest that for all studied ligands, the hydrophilic (hydrogen bonding) interactions with Glu172 are an important factor for the ligands' affinities toward the σ1 receptor. However, the hydrophobic interactions with Tyr120, Val162, Leu105, and Ile124 also significantly contribute to the ligand-receptor interplay and, in particular, differentiate the action of the agonistic morphine from the antagonistic haloperidol.

Tailoring the Schiff base photoswitching - a non-adiabatic molecular dynamics study of substituent effect on excited state proton transfer.

Small molecular systems exhibiting Excited State Intramolecular Proton Transfer (ESIPT) attract considerable attention due to their possible role as ultrafast, efficient, and photostable molecular photoswitches. Here, by means of static potential energy profile scan and on-the-fly non-adiabatic dynamics simulations we study the photodeactivation process of a minimal-chromophore aromatic Schiff base, salicylidene methylamine (SMA), and its two derivatives 6-cyano-salicylidene methylamine (6-CN-SMA) and 3-hydroxy-salicylidene methylamine (3-OH-SMA). We show that the dominant character of the lowest excited singlet state - ππ* vs. nπ* - plays a crucial role in the system's photophysics and controls the ESIPT efficiency. We also show that the relative alignment of the ππ* and nπ* states may be controlled through chemical substitutions made to the aromatic ring of the Schiff-base molecule. We believe that our findings will improve the rational-design strategies employed for the ESIPT systems, especially in the context of their possible photoswitching.

Conformational space of clindamycin studied by ab initio and full-atom molecular dynamics.

Molecular dynamics (MD) simulations allow determining internal flexibility of molecules at atomic level. Using ab initio Born-Oppenheimer molecular dynamics (BOMD), one can simulate in a reasonable time frame small systems with hundreds of atoms, usually in vacuum. With quantum mechanics/molecular mechanics (QM/MM) or full-atom molecular dynamics (FAMD), the influence of the environment can also be simulated. Here, we compare three types of MD calculations: ab initio BOMD, hybrid QM/MM, and classical FAMD. As a model system, we use a small antibiotic molecule, clindamycin, which is one of the lincosamide antibiotics. Clindamycin acquires two energetically stable forms and we investigated the transition between these two experimentally known conformers. We performed 60-ps BOMD simulations in vacuum, 50-ps QM/MM, and 100-ns FAMD in explicit water. The transition between two antibiotic conformers was observed using both BOMD and FAMD methods but was not noted in the QM/MM simulations.

Electric field control of proton-transfer molecular switching: molecular dynamics study on salicylidene aniline.

In this letter, we propose a novel, ultrafast, efficient molecular switch whose switching mechanism involves the electric field-driven intramolecular proton transfer. By means of ab initio quantum chemical calculations and on-the-fly dynamics simulations, we examine the switching performance of an isolated salicylidene aniline molecule and analyze the perspectives of its possible use as an electric field-controlled molecular electronics unit.

Calculations of NMR properties for sI and sII clathrate hydrates of methane, ethane and propane.

Calculations of NMR parameters (the absolute shielding constants and the spin-spin coupling constants) for 5(12), 5(12)6(2) and 5(12)6(4) cages enclathrating CH4, C2H6 and C3H8 molecules are presented. The DFT/B3LYP/HuzIII-su3 level of theory was employed. The (13)C shielding constants of guest molecules are close to available experimental data. In two cases (the ethane in 5(12) and the propane in 5(12)6(2) cages) the (13)C shielding constants are reported for the first time. Inversion of the methyl/methylene (13)C and (1)H shielding constants order is found for propane in the 5(12)6(2) cage. Topological criteria are used to interpret the changes of values of NMR parameters of water molecules and they establish a connection between single cages and bulk crystal.

Excited-state intramolecular proton transfer: photoswitching in salicylidene methylamine derivatives.

The effect of chemical substitutions on the photophysical properties of the salicylidene methylamine molecule (SMA) (J. Jankowska, M. F. Rode, J. Sadlej, A. L. Sobolewski, ChemPhysChem, 2012, 13, 4287-4294) is studied with the aid of ab initio electronic structure methods. It is shown that combining π-electron-donating and π-electron-withdrawing substituents results in an electron-density push-and-pull effect on the energetic landscape of the ground and the lowest excited ππ* and nπ* singlet states of the system. The presented search for the most appropriate SMA derivatives with respect to their photoswitching functionality offers an efficient prescreening tool for finding chemical structures before real synthetic realization.

Prediction of (L)-methionine VCD spectra in the gas phase and water solution.

In this paper we provide a computational study of the l-methionine conformational landscape and VCD spectra in the gas phase and a water environment simulated by implicit PCM and the hybrid model, i.e., a combination of explicit "microsolvation" and implicit models. In the gas phase, two groups of conformers differing in H-bonding, i.e., OH···NH2 and NH···O═C, could be distinguished based solely on the IR ν(OH) and ν(NH) stretching vibrations range. On the other hand, VCD better reflected chain differences. The most stable OH···NH2 conformer was predicted to be easily detected, and the presence of two out of four NH···O═C conformers could be confirmed. Three zwitterionic methionine conformers were shown to dominate in water. Their VCD spectra, simulated within the hybrid model at the B3LYP-IEF-PCM/aug-cc-pVDZ level of theory, indicated that they could be recognized in the mixture. Use of the hybrid model is crucial for good reproduction of the hydrogen bonding pattern in the VCD spectra of methionine in water solution. However, the 1300-800 cm(-1) region of the skeleton vibrations of methionine appeared to be relatively insensitive to the model of the solvent.

Photophysics of Schiff bases: theoretical study of salicylidene methylamine.

The proton-transfer reaction in a model aromatic Schiff base, salicylidene methylamine (SMA), in the ground and in the lowest electronically-excited singlet states, is theoretically analyzed with the aid of second-order approximate coupled-cluster model CC2, time-dependent density functional theory (TD-DFT) using the Becke, three-parameter Lee-Yang-Parr (B3LYP) functional, and complete active space perturbation theory CASPT2 electronic structure methods. Computed vertical-absorption spectra for the stable ground-state isomers of SMA fully confirm the photochromism of SMA. The potential-energy profiles of the ground and the lowest excited singlet state are calculated and four photophysically relevant isomeric forms of SMA; α, ß, γ, and δ are discussed. The calculations indicate two S(1)/S(0) conical intersections which provide non-adiabatic gates for a radiationless decay to the ground state. The photophysical scheme which emerges from the theoretical study is related to recent experimental results obtained for SMA and its derivatives in the low-temperature argon matrices (J. Grzegorzek, A. Filarowski, Z. Mielke, Phys. Chem. Chem. Phys. 2011, 13, 16596-16605). Our results suggest that aromatic Schiff bases are potential candidates for optically driven molecular switches.

Thermal fluctuations and infrared spectra of the formamide-formamidine complex.

Combined effects of hydrogen bonding and thermal fluctuations on the structure and infrared spectra of the formamide-formamidine dimer, FM···FI, are studied using ab initio molecular dynamics simulations. The equilibrium structure of the dimer is stabilized by two hydrogen bonds that form a pattern reminiscent of that found in the adenine-thymine base pair. The structure of the hydrogen bonds at 300 K is subject to large fluctuations, with the hydrogen atom being tightly bound to the donor in the covalent bonding scenario. The hydrogen bond acceptor has a tendency to detach farther away from the D-H pair, approaching the dimer dissociation limit. Moreover, the N-H···O hydrogen bond breaks occasionally, thus giving rise to an "open" structure of the dimer, while the N-H···N bond stays largely intact at this temperature. Thermal fluctuations result in the minor red shifts of the monomer vibrational frequencies indicative of the anharmonicity of the potential energy surface. In contrast, the IR frequencies of the two symmetric NH(2) vibrational modes of the FM and FI monomers are shifted substantially toward the red upon hydrogen bond formation in the FM···FI dimer. Dynamical effects studied here are relevant, in particular, to the hydrogen bonding of nucleic acids at finite temperatures.

On vibrational circular dichroism chirality transfer in electron donor-acceptor complexes: a prediction for the quinine···BF3 system.

Vibrational circular dichroism (VCD) chirality transfer occurs when an achiral molecule interacts with a chiral one and becomes VCD-active. Unlike for H-bonds, for organic electron donor-acceptor (EDA) complexes this phenomenon remains almost unknown. Here, the VCD chirality transfer from chiral quinine to achiral BF3 is studied at the B3LYP/aug-cc-pVDZ level. Accessibility of four quinine electron donor sites changes with conformation. Therefore, the quinine conformational landscape was explored and a considerable agreement between X-ray and the most stable conformer geometries was achieved. The BF3 complex through the aliphatic quinuclidine N atom is definitely dominating and is predicted to be easily recognizable in the VCD spectrum. Out of several VCD chirality transfer modes, the ν(s)(BF3) mode, the most intense in the entire VCD spectrum, satisfies the VCD mode robustness criterion and can be used for monitoring the chirality transfer phenomenon in quinine···BF3 system.

Structure, energetics, and infrared spectra of weakly bound HC(2n+1)N···HCl complexes. A theoretical study.

Three model systems, HCN···HCl, HC(3)N···HCl, and HC(5)N···HCl, have been investigated computationally with the use of the second-order Möller-Plesset (MP2) and the coupled cluster (with single and double excitations and noniterative inclusion of triples) methods. The global minima are linear hydrogen-bonded structures with HCl as a proton donor. Bent structures are proton-side complexes with HCl as an electron donor, while the bifurcated hydrogen bonds are predicted for t-shape complexes. One of the most important findings in this paper is that, according to symmetry-adapted perturbation analysis, the induction-to-dispersion ratios are the biggest for linear complexes, and it is the most noticeable difference between linear, bent, and t-shape structures.

Quantum mechanical studies of lincosamides.

Lincosamides are a class of antibiotics used both in clinical and veterinary practice for a wide range of pathogens. This group of drugs inhibits the activity of the bacterial ribosome by binding to the 23S RNA of the large ribosomal subunit and blocking protein synthesis. Currently, three X-ray structures of the ribosome in complex with clindamycin are available in the Protein Data Bank, which reveal that there are two distinct conformations of the pyrrolidinyl propyl group of the bound clindamycin. In this work, we used quantum mechanical methods to investigate the probable conformations of clindamycin in order to explain the two binding modes in the ribosomal 23S RNA. We studied three lincosamide antibiotics: clindamycin, lincomycin, and pirlimycin at the B3LYP level with the 6-31G** basis set. The focus of our work was to connect the conformational landscape and electron densities of the two clindamycin conformers found experimentally with their physicochemical properties. For both functional conformers, we applied natural bond orbital (NBO) analysis and the atoms in molecules (AIM) theory, and calculated the NMR parameters. Based on the results obtained, we were able to show that the structure with the intramolecular hydrogen bond C=O…H-O is the most stable conformer of clindamycin. The charge transfer between the pyrrolidine-derivative ring and the six-atom sugar (methylthiolincosamide), which are linked via an amide bond, was found to be the dominant factor influencing the high stability of this conformer.

Theoretical predictions of the spectroscopic parameters in noble-gas molecules: HXeOH and its complex with water.

We employ state-of-the-art methods and basis sets to study the effect of inserting the Xe atom into the water molecule and the water dimer on their NMR parameters. Our aim is to obtain predictions for the future experimental investigation of novel xenon complexes by NMR spectroscopy. Properties such as molecular structure and energetics have been studied by supermolecular approaches using HF, MP2, CCSD, CCSD(T) and MP4 methods. The bonding in HXeOH···H(2)O complexes has been analyzed by Symmetry-Adapted Perturbation Theory to provide the intricate insight into the nature of the interaction. We focus on vibrational spectra, NMR shielding and spin-spin coupling constants-experimental signals that reflect the electronic structures of the compounds. The parameters have been calculated at electron-correlated and Dirac-Hartree-Fock relativistic levels. This study has elucidated that the insertion of the Xe atom greatly modifies the NMR properties, including both the electron correlation and relativistic effects, the (129)Xe shielding constants decrease in HXeOH and HXeOH···H(2)O in comparison to Xe atom; the (17)O, as a neighbour of Xe, is deshielded too. The HXeOH···H(2)O complex in its most stable form is stabilized mainly by induction and dispersion energies.

Nuclear magnetic resonance parameters for methane molecule trapped in clathrate hydrates.

The calculations of the nuclear shielding and spin-spin coupling constants were carried out for two models of clathrate hydrates, 5(12) and 5(12)6(8), using the density functional theory three-parameter Becke-Lee-Yang-Parr method with the basis set aug-cc-pVDZ (optimization) and HuzIII-su3 (NMR parameters). Particular attention has been devoted to evaluate the influence of a geometrical arrangement, the effect of long-range interactions on the NMR shielding of methane molecule, and to predict whether (13)C and (1)H chemical shifts can distinguish between guests in two clathrate hydrates cages. The correlation of the changes in the (17)O shielding constants depend strongly on the hydrogen-bonding topology. The intermolecular hydrogen-bond transmitted (1h)J(OH) spin-spin coupling constants are substantial. The increase of their values is connected with the elongation of the intramolecular O-H bond and the shortening of the intermolecular O···H distance. These data suggests that hydrogen bonds between double donor-single acceptor (DDA)-type water molecules acting as a proton acceptor from single donor-double acceptor (DAA)-type water molecules are stronger than ones formed by DAA-type water molecules acting as an acceptor for a DDA water proton. These state-of-the-art calculations confirmed the earlier experimental findings of the cage-dependency of (13)C chemical shift of methane.

Calculated nuclear magnetic resonance parameters for multiproton-exchange and nonbonded-hydrogen rotation processes in cyclic water clusters.

In this work we report, for the first time, calculations of nuclear magnetic resonance parameters for the processes of multiproton-exchange and nonbonded-proton rotations in small, cyclic water clusters. The simultaneous proton exchange induces a large decrease in the oxygen shielding constants in both clusters, with a mean value of -52.6 ppm for the water trimer and -50.1 ppm for the water tetramer. The (1(h))J(OH) coupling constant between an oxygen nucleus and exchanging proton decreases (in absolute value) along the path, changes sign, finally reaching a value of 5-7 Hz. The changes in the NMR parameters induced by the nonbonded proton rotations are smaller. The calculated dependencies of the intermolecular spin-spin coupling constants upon rotation reveal that the largest changes are expected for the couplings transmitted through the hydrogen bond between the rotating and neighboring molecule which acts as a proton donor. The symmetry-adapted perturbation theory (SAPT) interaction energy calculations for each dimer forming the water trimer have allowed us to relate a strength of interactions within pairs of water molecules with coupling constant values. The predicted maximal values of the interaction-energy terms (energetically unfavorable orientations of the constituent dimers) along paths correlate with the extremal values of the spin-spin coupling constants, which mostly correspond to the maximal couplings along pathways.