Characterization of the binding interaction between atrazine and human serum albumin: Fluorescence spectroscopy, molecular dynamics and quantum biochemistry.

Atrazine (ATR), one of the most used herbicides worldwide, causes persistent contamination of water and soil due to its high resistance to degradation. ATR is associated with low fertility and increased risk of prostate cancer in humans, as well as birth defects, low birth weight and premature delivery. Describing ATR binding to human serum albumin (HSA) is clinically relevant to future studies about pharmacokinetics, pharmacodynamics and toxicity of ATR, as albumin is the most abundant carrier protein in plasma and binds important small biological molecules. In this work we characterize, for the first time, the binding of ATR to HSA by using fluorescence spectroscopy and performing simulations using molecular docking, classical molecular dynamics and quantum biochemistry based on density functional theory (DFT). We determine the most likely binding sites of ATR to HSA, highlighting the fatty acid binding site FA8 (located between subdomains IA-IB-IIA and IIB-IIIA-IIIB) as the most important one, and evaluate each nearby amino acid residue contribution to the binding interactions explaining the fluorescence quenching due to ATR complexation with HSA. The stabilization of the ATR/FA8 complex was also aided by the interaction between the atrazine ring and SER454 (hydrogen bond) and LEU481(alkyl interaction).

The vibrational properties of the bee-killer imidacloprid insecticide: A molecular description.

The chemical imidacloprid belongs to the neonicotinoids insecticide class, widely used for insect pest control mainly for crop protection. However, imidacloprid is a non-selective agrochemical to the insects and it is able to kill the most important pollinators, the bees. The high toxicity of imidacloprid requires controlled release and continuous monitoring. For this purpose, high performance liquid chromatography (HPLC) is usually employed; infrared and Raman spectroscopy, however, are simple and viable techniques that can be adapted to portable devices for field application. In this communication, state-of-the-art quantum level simulations were used to predict the infrared and Raman spectra of the most stable conformer of imidacloprid. Four molecular geometries were investigated in vacuum and solvated within the Density Functional Theory (DFT) approach employing the hybrid meta functional M06-2X and the hybrid functional B3LYP. The M062X/PCM model proved to be the best to predict structural features, while the values of harmonic vibrational frequencies were predicted more accurately using the B3LYP functional.

Encapsulation of nor-ß-lapachone into poly(d,l)-lactide-co-glycolide (PLGA) microcapsules: full characterization, computational details and cytotoxic activity against human cancer cell lines.

In this work, we characterize nor-ß-lapachone-loaded (NßL-loaded) microcapsules prepared using an emulsification/solvent extraction technique. Features such as surface morphology, particle size distribution, zeta potential, optical absorption, Raman and Fourier transform infrared spectra, thermal analysis data, drug encapsulation efficiency, drug release kinetics and in vitro cytotoxicity were studied. Spherical microcapsules with a size of 1.03 ± 0.46 µm were produced with an encapsulation efficiency of approximately 19%. Quantum DFT calculations were also performed to estimate typical interaction energies between a single nor-ß-lapachone molecule and the surface of the microparticles. The NßL-loaded PLGA microcapsules exhibited a pronounced initial burst release. After the in vitro treatment with NßL-loaded microcapsules, a clear phagocytosis of the spheres was observed in a few minutes. The cytotoxic activity against a set of cancer cell lines was investigated.

DFT Calculations with van der Waals Interactions of Hydrated Calcium Carbonate Crystals CaCO3·(H2O, 6H2O): Structural, Electronic, Optical, and Vibrational Properties.

The role of hydration on the structural, electronic, optical, and vibrational properties of monohydrated (CaCO3·H2O, hexagonal, P31, Z = 9) and hexahydrated (CaCO3·6H2O, monoclinic, C2/c, Z = 4) calcite crystals is assessed with the help of published experimental and theoretical data applying density functional theory within the generalized gradient approximation and a dispersion correction scheme. We show that the presence of water increases the main band gap of monohydrocalcite by 0.4 eV relative to the anhydrous structure, although practically not changing the hexahydrocalcite band gap. The gap type, however, is modified from indirect to direct as one switches from the monohydrated to the hexahydrated crystal. A good agreement was obtained between the simulated vibrational infrared and Raman spectra and the experimental data, with an infrared signature of hexahydrocalcite relative to monohydrocalcite being observed at 837 cm(-1). Other important vibrational signatures of the lattice, water molecules, and CO3(2-) were identified as well. Analysis of the phonon dispersion curves shows that, as the hydration level of calcite increases, the longitudinal optical-transverse optical phonon splitting becomes smaller. The thermodynamics properties of hexahydrocalcite as a function of temperature resemble closely those of calcite, while monohydrocalcite exhibits a very distinct behavior.

Controlled Release of Nor-ß-lapachone by PLGA Microparticles: A Strategy for Improving Cytotoxicity against Prostate Cancer Cells.

Prostate cancer is one of the most common malignant tumors in males and it has become a major worldwide public health problem. This study characterizes the encapsulation of Nor-ß-lapachone (NßL) in poly(d,l-lactide-co-glycolide) (PLGA) microcapsules and evaluates the cytotoxicity of the resulting drug-loaded system against metastatic prostate cancer cells. The microcapsules presented appropriate morphological features and the presence of drug molecules in the microcapsules was confirmed by different methods. Spherical microcapsules with a size range of 1.03 ± 0.46 µm were produced with an encapsulation efficiency of approximately 19%. Classical molecular dynamics calculations provided an estimate of the typical adsorption energies of NßL on PLGA. Finally, the cytotoxic activity of NßL against PC3M human prostate cancer cells was demonstrated to be significantly enhanced when delivered by PLGA microcapsules in comparison with the free drug.

Two Binding Geometries for Risperidone in Dopamine D3 Receptors: Insights on the Fast-Off Mechanism through Docking, Quantum Biochemistry, and Molecular Dynamics Simulations.

Risperidone is an atypical antipsychotic used in the treatment of schizophrenia and of symptoms of irritability associated with autism spectrum disorder (ASD). Its main action mechanism is the blockade of D2-like receptors acting over positive and negative symptoms of schizophrenia with small risk of extrapyramidal symptoms (EPS) at doses corresponding to low/moderate D2 occupancy. Such a decrease in the side effect incidence can be associated with its fast unbinding from D2 receptors in the nigrostriatal region allowing the recovery of dopamine signaling pathways. We performed docking essays using risperidone and the D3 receptor crystallographic data and results suggested two possible distinct orientations for risperidone at the binding pocket. Orientation 1 is more close to the opening of the binding site and has the 6-fluoro-1,2 benzoxazole fragment toward the bottom of the D3 receptor cleft, while orientation 2 is deeper inside the binding pocket with the same fragment toward to the receptor surface. In order to unveil the implications of these two binding orientations, classical molecular dynamics and quantum biochemistry computations within the density functional theory formalism and the molecular fractionation with conjugate caps framework were performed. Quantum mechanics/molecular mechanics suggests that orientation 2 (considering the contribution of Glu90) is slightly more energetically stable than orientation 1 with the main contribution coming from residue Asp110. The residue Glu90, positioned at the opening of the binding site, is closer to orientation 1 than 2, suggesting that it may have a key role in stability through attractive interaction with risperidone. Therefore, although orientations 1 and 2 are both likely to occur, we suggest that the occurrence of the first may contribute to the reduction of side effects in patients taking risperidone due to the reduction of dopamine receptor occupancy in the nigrostriatal region through a mechanism of fast dissociation. The atypical effect may be obtained simply by either delaying D3R full blockage by spatial hindrance of orientation 1 at the binding site or through an effective blockade followed by orientation 1 fast dissociation. While the molecular interpretation suggested in this work shed some light on the potential molecular mechanisms accounting for the reduced extrapyramidal symptoms observed during risperidone treatment, further studies are necessary in order to evaluate the implications of both orientations during the receptor activation/inhibition. Altogether these data highlight important hot spots in the dopamine receptor binding site bringing relevant information for the development of novel/derivative agents with atypical profile.

A quantum biochemistry investigation of willardiine partial agonism in AMPA receptors.

We employ quantum biochemistry methods based on the Density Functional Theory (DFT) approach to unveil the detailed binding energy features of willardiines co-crystallized with the AMPA receptor. Our computational results demonstrate that the total binding energies of fluorine-willardiine (FW), hydrogen-willardiine (HW), bromine-willardiine (BrW) and iodine-willardiine (IW) to the iGluR2 ligand-pocket correlate with the agonist binding energies, whose experimental sequential data match our computational counterpart, excluding the HW case. We find that the main contributions to the total willardiine-iGluR2 binding energy are due to the amino acid residues in decreasing order Glu705 > Arg485 > Ser654 > Tyr450 > T655. Furthermore, Met708, which is positioned close to the 5-substituent, attracts HW and FW, but repels BrW and IW. Our results contribute significantly to an improved understanding of the willardiine-iGluR2 binding mechanisms.

Antipsychotic haloperidol binding to the human dopamine D3 receptor: beyond docking through QM/MM refinement toward the design of improved schizophrenia medicines.

As the dopamine D3R receptor is a promising target for schizophrenia treatment, an improved understanding of the binding of existing antipsychotics to this receptor is crucial for the development of new potent and more selective therapeutic agents. In this work, we have used X-ray cocrystallization data of the antagonist eticlopride bound to D3R as a template to predict, through docking essays, the placement of the typical antipsychotic drug haloperidol at the D3R receptor binding site. Afterward, classical and quantum mechanics/molecular mechanics (QM/MM) computations were employed to improve the quality of the docking calculations, with the QM part of the simulations being accomplished by using the density functional theory (DFT) formalism. After docking, the calculated QM improved total interaction energy EQMDI = -170.1 kcal/mol was larger (in absolute value) than that obtained with classical molecular mechanics improved (ECLDI = -156.3 kcal/mol) and crude docking (ECRDI = -137.6 kcal/mol) procedures. The QM/MM computations reveal the pivotal role of the Asp110 amino acid residue in the D3R haloperidol binding, followed by Tyr365, Phe345, Ile183, Phe346, Tyr373, and Cys114. Besides, it highlights the relevance of the haloperidol hydroxyl group axial orientation, which interacts with the Tyr365 and Thr369 residues, enhancing its binding to dopamine receptors. Finally, our computations indicate that functional substitutions in the 4-clorophenyl and in the 4-hydroxypiperidin-1-yl fragments (such as C3H and C12H hydrogen replacement by OH or COOH) can lead to haloperidol derivatives with distinct dopamine antagonism profiles. The results of our work are a first step using in silico quantum biochemical design as means to impact the discovery of new medicines to treat schizophrenia.

Optical absorption of the antitrypanocidal drug benznidazole in water.

UV-vis optical absorption spectra of the antitrypanocidal drug benznidazole solvated in water were measured for various concentrations. The spectra show a prominent peak around 3.80 eV, while deconvolution of the UV-vis optical absorption spectra revealed six bands centered at 3.60, 3.83, 4.15, 4.99, 5.60, and 5.76 eV. Benznidazole electronic transitions were obtained after density functional theory (DFT) calculations within the polarized continuum (PCM) model for water solvation. Molecular geometry optimizations were carried out, and the measured absorption peaks were related to specific molecular orbital transitions obtained within the time dependent DFT (TD-DFT) with excellent agreement between theory and experiment.

Inactivation of ovine cyclooxygenase-1 by bromoaspirin and aspirin: a quantum chemistry description.

The resulting noncovalent bonding of the salicylic acid to ovine COX-1 after bromoaspirin and aspirin acetylation by Ser530 is investigated within the scope of density functional theory considering a 6.5 Å radius binding pocket. We have not only took full advantage of published X-ray structural data for the ovine COX-1 cocrystallized with bromoaspirin, but we also have improved that data through computation, finding good estimates for the hydrogen atom positions at the residues of the binding pocket, and repositioning the Ser530Ac[Br;H] lateral chain and salicylic acid by total energy minimization procedures employing LDA and GGA+D exchange-correlation functionals. Using bromoaspirin as a template, we have simulated the positioning of aspirin in the binding pocket, estimating its interaction energy with each of its neighbor COX-1 residues. We demonstrate that the binding energies of bromoaspirin and aspirin to COX-1 are very close when second-order quantum refinements of the structural data are performed, which points to an explanation on why the IC(50) values for the 126 µM COX-1 activity of both bromoaspirin and aspirin are practically the same. Attracting and repelling residues were identified, being shown that Arg120 is the most effective residue attracting the salicylic acid, followed by Ala527, Leu531, Leu359, and Ser353. On the other hand, Glu524 was found the most effective repulsive residue (strength interaction comparable to Arg120).

Explaining statin inhibition effectiveness of HMG-CoA reductase by quantum biochemistry computations.

By taking advantage of the crystallographic data of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) complexed with statins, a quantum biochemistry study based on the density functional theory is performed to estimate the interaction energy for each statin when one considers binding pockets of different sizes. Assuming a correlation between statin potency and the strength of the total HMGR-statin binding energy, clinical data as well as IC(50) values of these cholesterol-lowering drugs are successfully explained only after stabilization of the calculated total binding energy for a larger size of the ligand-interacting HGMR region, one with a radius of at least 12.0 Å. Actually, the binding pocket radius suggested by classic works, which was based solely on the interpretation of crystallographic data of the HMGR-statin complex, is smaller than that necessary to achieve total binding energy convergence in our simulations. Atorvastatin and rosuvastatin are shown to be the most strongly bound HMGR inhibitors, while simvastatin and fluvastatin are the weakest ones. A binding site, interaction energy between residues and statin atoms, and residues domain (BIRD) panel is constructed, indicating clear quantum biochemistry-based routes for the development of new statin derivatives.