Determination of key receptor–ligand interactions of dopaminergic arylpiperazines and the dopamine D2 receptor homology model

Interest in structure-based G-protein-coupled receptor (GPCR) ligand discovery is huge, given that almost 30 % of all approved drugs belong to this category of active compounds. The GPCR family includes the dopamine receptor subtype D2 (D2DR), but unfortunately—as is true of most GPCRs—no experimental structures are available for these receptors. In this publication, we present the molecular model of D2DR based on the previously published crystal structure of the dopamine D3 receptor (D3DR). A molecular modeling study using homology modeling and docking simulation provided a rational explanation for the behavior of the arylpiperazine ligand. The observed binding modes and receptor–ligand interactions provided us with fresh clues about how to optimize selectivity for D2DR receptors.

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

Charge sensitivity approach to mutual polarization of reactants: molecular mechanics perspective

Charge sensitivity analysis (CSA) in force-field atoms resolution was applied to describe the mutual polarization of reactants as well as charge-transfer (CT) effects. An inclusion complex of β-cyclodextrin with salicylic acid was used as a model system. Three CSA models were taken into account and verified on a Born–Oppenheimer molecular dynamics (BOMD) trajectory. The models differed in terms of the equilibrium conditions imposed on the system. It was demonstrated that mutual polarization is an important source of stabilization, in contrast to the results obtained from static charge calculations. The energy lowering induced by CT was small and comparable to the CT stabilization that occurs in hydrogen-bonded systems. All models correctly described the main topological features of the BOMD energy surface. CSA in force-field atoms resolution qualitatively reproduced the charge reorganization accompanying hydrogen-bond formation. It was shown that CSA parameters are very sensitive to the bond formation process, which suggests that they could be applied in reactive force fields as detectors of newly formed chemical bonds.

Simulations on the possibility of formation of complexes between fluorouracil drug and cucurbit[n]urils: ab initio van der Waals DFT study

The binding geometry of fluorouracil/cucurbit[n]urils (CB[n]s) complexes with n?=?5–8 is investigated using the first-principles van der Waals density functional (vdW-DF) method, involving full geometry optimization. Such host-guest complexes are typically calculated using conventional DFT method, which significantly underestimates non-local dispersion forces (or vdW contributions) and therefore affects interactions between respected entities. We address here the role of vdW forces for the fluorouracil and CB[n]s molecules which can form directional hydrogen bonds with each other. It was found that the inclusion of dispersion interactions significantly affects the host-guest binding properties and the hydrogen bonding between the molecules provides the main binding mechanism, while results in the same geometries for the considered complexes. The 0.84 eV binding energy, for the thermodynamically favorable state, reveals that the interaction of fluorouracil with CB[n]s is an exothermic interaction and typical for strong hydrogen bonding. Furthermore, we have investigated the binding nature of these host-guest systems in aqueous solution with ab initio MD simulations adopting vdW-DF method. These findings afford evidence for the applicability of the vdW-DF approach and provide a realistic benchmark for the investigation of the host-guest complexes.

Assessing the quantum mechanical level of theory for prediction of linear and nonlinear optical properties of push-pull organic molecules

In this paper, we assessed the quantum mechanical level of theory for prediction of linear and nonlinear optical (NLO) properties of push-pull organic molecules. The electric dipole moment (μ), mean polarizability (〈α〉) and total static first hyperpolarizability (β_t) were calculated for a set of benzene, styrene, biphenyl and stilbene derivatives using HF, MP2 and DFT (31 different functionals) levels and over 71 distinct basis sets. In addition, we propose two new basis sets, NLO-V and aNLO-V, for NLO properties calculations. As the main outcomes it is shown that long-range corrected DFT functionals such as M062X, ωB97, cam-B3LYP, LC-BLYP and LC-ωPBE work satisfactorily for NLO properties when appropriate basis sets such as those proposed here (NLO-V or aNLO-V) are used. For most molecules with β ranging from 0 to 190 esu, the average absolute deviation was 13.2 esu for NLO-V basis sets, compared to 27.2 esu for the standard 6-31 G(2d) basis set. Therefore, we conclude that the new basis sets proposed here (NLO-V and aNLO-V), together with the cam-B3LYP functional, make an affordable calculation scheme to predict NLO properties of large organic molecules.

Theoretical study of the solvatochromism of a donor-acceptor bithiophene

The solvation and the solvatochromic behavior of the 5-(methylthio)-5′-nitro-2,2′-bithiophene 1 in diethyl ether, dichloromethane, acetonitrile, methanol and formamide was theoretically investigated with an iterative molecular and quantum mechanics (QM/MM) approach. Calculated longest-wavelength solvatochromic absorption band of 1, obtained as averages of statistically uncorrelated configurations, including the solute and explicit solvent molecules of the first and second solvation layer, were in excellent agreement with the experimental results.

Modeling the effect of H-bonding interactions and molecular packing on the molecular structure of [Ag(ethylnicotinate)2]NO3 complex

The gas phase molecular structure of a single isolated molecule of [Ag(Etnic)₂NO₃];1 where Etnic = Ethylnicotinate was calculated using B3LYP method. The H-bonding interaction between 1 with one (complex 2) and two (complex 3) water molecules together with the dimeric formula [Ag(Etnic)₂NO₃]₂;4 and the tetrameric formula [Ag(Etnic)₂NO₃]₄;5 were calculated using the same level of theory to model the effect of intermolecular interactions and molecular packing on the molecular structure of the titled complex. The H-bond dissociation energies of complexes 2 and 3 were calculated to be in the range of 12.220–14.253 and 30.106–31.055 kcal?mol^?1, respectively, indicating the formation of relatively strong H-bonds between 1 and water molecules. The calculations predict bidentate nitrate ligand in the case of 1 and 2, leading to distorted tetrahedral geometry around the silver ion with longer Ag–O distances in case of 2 compared to 1, while 3 has a unidentate nitrate ligand leading to a distorted trigonal planar geometry. The packing of two [Ag(Etnic)₂NO₃] complex units; 4 does not affect the molecular geometry around Ag(I) ion compared to 1. In the case of 5, the two asymmetric units of the formula [Ag(Etnic)₂NO₃] differ in the bonding mode of the nitrate group, where the geometry around the silver ion is distorted tetrahedral in one unit and trigonal planar in the other. The calculations predicted almost no change in the charge densities at the different atomic sites except at the sites involved in the C–H?O interactions as well as at the coordinated nitrogen of the pyridine ring.

Enhanced Dissolution and Stability of Lansoprazole by Cyclodextrin Inclusion Complexation: Preparation, Characterization, and Molecular Modeling

In this study, lansoprazole (LSP)/cyclodextrin (CD) inclusion complexes were prepared using a fluid bed coating technique, with β-cyclodextrin (β-CD) and 2-hydroxypropyl-β-cyclodextrin (HPCD) as the host molecules, respectively, to simultaneously improve the dissolution and stability of LSP. The dissolution rate and stability of LSP was dramatically enhanced by inclusion complexation regardless of CD type. LSP/HPCD inclusion complex was more stable under illumination than LSP/β-CD inclusion complex. Differential scanning calorimetry and powder X-ray diffractometry proved the absence of crystallinity in both LSP/CD inclusion complexes. Fourier transform infrared spectroscopy together with molecular modeling indicated that the benzimidazole of LSP was included in the cavity of both CDs, while LSP was more deeply included in HPCD than β-CD. The enhanced photostability was due to the inclusion of the sulfinyl moiety into the HPCD cavity. CD inclusion complexation could improve the dissolution and stability of LSP.KEY WORDS: cyclodextrin, dissolution, inclusion complex, lansoprazole, molecular modeling, stability     

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Rational design of the survivin/CDK4 complex by combining protein–protein docking and molecular dynamics simulations

Survivin, the smallest inhibitor of apoptosis protein (IAP), is a valid target for cancer research. It mediates both the apoptosis pathway and the cell cycle and has been proposed to form a complex with the cyclin-dependent kinase protein CDK4. The resulting complex transports CDK4 from the cytosol to the nucleus, where CDK4 participates in cell division. Survivin has been recognized as a node protein that interacts with several partners; disruption of the formed complexes can lead to new anticancer compounds. We propose a rational model of the survivin/CDK4 complex that fulfills the experimental evidence and that can be used for structure-based design of inhibitors modifying its interface recognition. In particular, the suggested complex involves the alpha helical domain of survivin and resembles the mode of binding of survivin in the survivin/borealin X-ray structure. The proposed model has been obtained by combining protein–protein docking, fractal-based shape complementarity, electrostatics studies and extensive molecular dynamics simulations.

Theoretical study on cooperative effects between X⋯N and X⋯Carbene halogen bonds (X = F,Cl,Br and I)

Quantum chemical calculations are performed to study the interplay between halogen?nitrogen and halogen?carbene interactions in NCX?NCX?CH₂ complexes, where X?=?F, Cl, Br and I. Molecular geometries and interaction energies of dyads and triads are investigated at the MP2/aug-cc-pVTZ level of theory. It is found that the X?N and X?C_carbene interaction energies in the triads are larger than those in the dyads, indicating that both the halogen bonding interactions are enhanced. The estimated values of cooperative energy E _coop are all negative with much larger E _coop in absolute value for the systems including iodine. The nature of halogen bond interactions of the complexes is analyzed using parameters derived from the quantum theory atoms in molecules methodology and energy decomposition analysis.

Vibrational spectra of an RDX film over an aluminum substrate from molecular dynamics simulations and density functional theory

We report calculated vibrational spectra in the range of 0–3,500 cm^?1 of RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) molecules adsorbed on a model aluminum surface. A molecular film was modeled using two approaches: (1) density functional theory (DFT) was used to optimize a single RDX molecule interacting with its periodic images, and (2) a group of nine molecules extracted from the crystal structure was deposited on the surface and interacted with its periodic images via molecular dynamics (MD) simulations. In both cases, the molecule was initialized in the AAA conformer geometry having the three nitro groups in axial positions, and kept that conformation in the DFT examination, but some molecules were found to change to the AAE conformer (two nitro groups in axial and one in equatorial position) in the MD analysis. The vibrational spectra obtained from both methods are similar to each other, except in the regions where collective RDX intermolecular interactions (captured by MD simulations) are important, and compare fairly well with experimental findings.

Characterization, phase-solubility, and molecular modeling of inclusion complex of 5-nitroindazole derivative with cyclodextrins

The slightly water-soluble 5-nitroindazole derivative (5-NI) and its inclusion with either β-cyclodextrin (βCD) or Heptakis (2,6-di-O-methyl)-β-cyclodextrin (DMβCD) were investigated. The stoichiometric ratios and stability constants describing the extent of formation of the complexes were determined by phase-solubility measurements obtaining type-A_L diagrams in both cases. According to the continuous variation method (Job’s plot) a 1:1 stoichiometry has been proposed for the complexes. Also electrochemical studies were carried out on both CDs complexes, where the observed change in the E_PC value for DMβCD indicated a lower feasibility of the nitro group reduction. The detailed spatial configuration is proposed based on two-dimensional NMR methods. These results are further interpreted using molecular modeling studies. The latter results are in good agreement with the experimental data.     

Microhydration of caesium compounds: Cs,CsOH, CsI and Cs2I2 complexes with one to three H2O molecules of nuclear safety interest

The structure and properties of natural gas hydrates containing hydrocarbons, CO₂, and N₂ molecules were studied by using computational quantum chemistry methods via the density functional theory approach. All host cages involved in I, II, and H types structures where filled with hydrocarbons up to pentanes, CO₂ and N₂ molecules, depending on their size, and the structures of these host–guest systems optimized. Structural properties, vibrational spectra, and density of states were analyzed together with results from atoms-in-a-molecule and natural bond orbitals methods. The inclusion of dispersion terms in the used functional plays a vital role for obtaining reliable information, and thus, B97D functional was shown to be useful for these systems. Results showed remarkable interaction energies, not strongly affected by the type of host cage, with molecules tending to be placed at the center of the cavities when host cages and guest molecules cavities are of similar size, but with molecules approaching hexagonal faces for larger cages. Vibrational properties show remarkable features in certain regions, with shiftings rising from host-guest interactions, and useful patterns in the terahertz region rising from water surface vibrations strongly coupled with guest molecules. Likewise, calculations on crystal systems for the I and H types were carried out using a pseudopotential approach combined with Grimme’s method to take account of dispersion.

Cytotoxic effect and molecular docking of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide—a novel topoisomerase II inhibitor

The preliminary cytotoxic effect of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide hydrochloride (1)—a potent topoisomerase II inhibitor—was measured using a MTT assay. It was found that the compound decreased the number of viable cells in both estrogen receptor-positive MCF-7 and estrogen receptor-negative MDA-MB-231breast cancer cells, with IC₅₀ values of 146?±?2 and 132?±?2 μM, respectively. To clarify the molecular basis of the inhibitory action of 1, molecular docking studies were carried out. The results suggest that 1 targets the ATP binding pocket.

Theoretical studies of the interaction between enflurane and water

Increase of the atmospheric concentration of halogenated organic compounds is partially responsible for a change of the global climate. In this work we have investigated the interaction between halogenated ether and water, which is one of the most important constituent of the atmosphere. The structures of the complexes formed by the two most stable conformers of enflurane (a volatile anaesthetic) with one and two water molecules were calculated by means of the counterpoise CP-corrected gradient optimization at the MP2/6–311++G(d,p) level. In these complexes the CH…O_w hydrogen bonds are formed, with the H…O_w distances varying between 2.23 and 2.32 Å. A small contraction of the CH bonds and the blue shifts of the ν(CH) stretching vibrations are predicted. There is also a weak interaction between one of the F atoms and the H atom of water, with the H_w…F distances between 2.41 and 2.87 Å. The CCSD(T)/CBS calculated stabilization energies in these complexes are between ?5.89 and ?4.66 kcal?mol^?1, while the enthalpies of formation are between ?4.35 and ?3.22 kcal?mol^?1. The Cl halogen bonding between enflurane and water has been found in two complexes. The intermolecular (Cl···O) distance is smaller than the sum of the corresponding van der Waals radii. The CCSD(T)/CBS stabilization energies for these complexes are about ?2 kcal?mol^?1.

Structural and energetic insights into sequence-specific interaction in DNA–drug recognition: development of affinity predictor and analysis of binding selectivity

Although the molecular mechanism and thermodynamic profile of a wide variety of chemical agents have been examined intensively in the past decades in terms of specific recognition of their protein receptors, to date the physicochemical nature of DNA–drug recognition and association still remains largely unexplored. The present study focused on understanding the structural basis, energetic landscape, and biological implications underlying the binding of small-molecule ligands to their cognate or non-cognate DNA receptors. First, a new method to capture the structural features of DNA–drug complex architecture was proposed and then used to correlate the extracted features with binding affinity of the complexes. By employing this method, a statistical regression-based predictor was developed to quantitatively evaluate the interaction potency of drug compounds with DNA in a fast and reliable manner. Subsequently, we use the predictor to examine the binding behavior of a number of structure-available, affinity-known DNA–drug complexes as well as a large pool of randomly generated DNA decoys in complex with the same drugs. It was found that (1) as compared with protein–DNA recognition, small-molecule agents have relatively low specificity in selecting their cognate DNA targets from the background of numerous random decoys; (2) the abundance of A–T base pairs in the DNA core motif exhibits a significant positive correlation with the affinity of drug ligand binding to the DNA receptor; and (3) high affinity seems not to be closely related to high selectivity for a DNA-targeting drug, although high-affinity drug entities have a greater possibility of being ranked computationally as top binders. We hope that this work will provide a preliminary insight into the molecular origin of sequence-specific interactions in DNA–drug recognition.

Structure and spectral characteristics of diquat-cucurbituril complexes from density functional theory

Electronic structure, ¹H NMR and infrared spectra of diquat (6,7-dihydrodipyrido[1,2-b:1′,2′-e] pyrazine-5,8-diium or DQ²⁺) encapsulated by cucurbit[n]uril (n?=?7,8) hosts are obtained using the density functional theory. Theoretical calculations have shown that both CB[7] or CB[8] host possesses strong affinity toward DQ²⁺ compared to its reduced cation or neutral species. Calculated ¹H NMR spectra reveal that H_α protons on bi-pyridinium rings of DQ²⁺@CB[8] complex are de-shielded owing to C=O?H interactions. On the other hand aromatic (H_β and H_δ) of DQ²⁺ within the CB[8] cavity exhibit significant shielding. The complexation of CB[8] with DQ²⁺ splits the carbonyl stretching vibration (1788 cm^?1) into two distinct vibrations which correspond to 1765 cm^?1 arising from hydrogen bonded carbonyls and the 1792 cm^?1 band from non-interacting ones. Further, the CN stretching vibration in DQ²⁺ exhibits a frequency blue-shift of 6 cm^?1 on its encapsulation within the CB[8] cavity. The direction of frequency shift has been explained on the basis of natural bond orbital analyses.

Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies

We describe a procedure for performing quantitative analyses of fields f(r) on molecular surfaces, including statistical quantities and locating and evaluating their local extrema. Our approach avoids the need for explicit mathematical representation of the surface and can be implemented easily in existing graphical software, as it is based on the very popular representation of a surface as collection of polygons. We discuss applications involving the volumes, surface areas and molecular surface electrostatic potentials, and local ionization energies of a group of 11 molecules.

Discovery of potent inhibitors for interleukin-2-inducible T-cell kinase: structure-based virtual screening and molecular dynamics simulation approaches

In our study, a structure-based virtual screening study was conducted to identify potent ITK inhibitors, as ITK is considered to play an important role in the treatment of inflammatory diseases. We developed a structure-based pharmacophore model using the crystal structure (PDB ID: 3MJ2) of ITK complexed with BMS-50944. The most predictive model, SB-Hypo1, consisted of six features: three hydrogen-bond acceptors (HBA), one hydrogen-bond donor (HBD), one ring aromatic (RA), and one hydrophobic (HY). The statistical significance of SB-Hypo1 was validated using wide range of test set molecules and a decoy set. The resulting well-validated model could then be confidently used as a 3D query to screen for drug-like molecules in a database, in order to retrieve new chemical scaffolds that may be potent ITK inhibitors. The hits retrieved from this search were filtered based on the maximum fit value, drug-likeness, and ADMET properties, and the hits that were retained were used in a molecular docking study to find the binding mode and molecular interactions with crucial residues at the active site of the protein. These hits were then fed into a molecular dynamics simulation to study the flexibility of the activation loop of ITK upon ligand binding. This combination of methodologies is a valuable tool for identifying structurally diverse molecules with desired biological activities, and for designing new classes of selective ITK inhibitors.