Theoretical study of the hydrogen-bonded complexes serotonin–water/hydrogen peroxide

Eight H-bonded complexes between serotonin (5-hydroxy-tryptamine) and water/hydrogen peroxide were studied at the B3LYP and HF levels of theory, using the 6-31+G(d) basis set. A thermodynamic analysis was performed in order to find the most stable complex. The calculated bonding parameters showed that the most stable H-bonded complex is formed between serotonin and hydrogen peroxide by means of the intermolecular H-bond –H₂N...H–OOH. Fig. a Theoretical study of the hydrogen-bonded supersystems serotonin-water/hydrogen peroxide     

A theoretical study of some new analogues of the anti-cancer drug camptothecin

The enzyme topoisomerase I (topo I), which is essential for cell replication, transiently causes a DNA single strand break and makes a complex with it. The anti-cancer agent camptothecin (CPT) binds to the topo I–DNA complex and stabilizes it, preventing resealing of the broken DNA strand and cell growth. Considering the structural factors of CPT that are believed to be involved in stabilizing the topo I–DNA complex via hydrogen bonding and stacking interactions, designs of two new analogues of CPT (topo I inhibitors) have been suggested. The molecular geometries of CPT, two of its analogues and certain other related molecules included in the study were fully optimized in both gas phase and aqueous media at the B3LYP/6-311++G(d,p) level of density functional theory. Solvation effects of aqueous media were treated using the polarizable continuum model (PCM). Net CHelpG charges and surface molecular electrostatic potentials (MEP) near the atomic sites of the molecules were studied. Structural analogy and surface MEP values suggests that the two new CPT analogues studied here would be potent topoisomerase I inhibitors. Figure Optimized structures of CPT and two of its new analogues, 10 and 11.     

Possible dynamic anchor points in a benzoxazinone derivative–human oxytocin receptor system — a molecular docking and dynamics calculation

In this study, we performed a molecular docking and dynamics simulation for a benzoxazinone–human oxytocin receptor system to determine the possible hydrophobic and electrostatic interaction points in the dynamic complex. After the homology modeling, the ligand was docked into the putative active using AutoDock 3.05. After the application of energetic and structural filters, the complexes obtained were further refined with a simulated annealing protocol (AMBER8) to remove steric clashes. Three complexes were selected for subjection to the molecular dynamics simulation (5 ns), and the results on the occurrence of average anchor points showed a stable complex between the benzoxazinone derivative and the receptor. The complex could be used as a good starting point for further analysis with site-directed mutagenesis, or further computational research. Figure The location of the ligands (complex B – blue; complex E – red; and complex F – green) in the transmembrane regions (TM1 – red; TM2 – blue; TM3 – yellow; TM4 – purple; TM5 – orange; TM6 – cyan; TM7 – pink) of the hOTR. For clarity, the EC and IC loops are not shown Electronic Supplementary Material Supplementary material is available at     

Thrombin inhibitors with novel P1 binding pocket functionality: free energy of binding analysis

The high incidence of thrombembolic diseases justifies the development of new antithrombotics. The search for a direct inhibitor has resulted in the synthesis of a considerable number of low molecular weight molecules that inhibit human α-thrombin potently. However, efforts to develop an orally active drug remain in progress as the most active inhibitors with a highly basic P1 moiety exhibit an unsatisfactory bioavailability profile. In our previous work we solved several X-ray structures of human α-thrombin in complexes with (1) novel bicyclic arginine mimetics attached to the glycylproline amide and pyridinone acetamide scaffold and (2) inhibitors with a novel aza scaffold and with charged or neutral P1 moieties. In the present contribution, we correlate the structures of the complex between these inhibitors and the protein with the calculated free energy of binding. The energy of solvation was calculated using the Poisson–Boltzmann approach. In particular, the requirements for successful recognition of an inhibitor at the protein’s active site pocket S1 are discussed. Figure We report here on free energy of binding analysis of thrombin inhibitors with novel aza scaffold and novel bicyclic arginine mimetics in S1 pocket of thrombin     

Reverse-docking study of the TADDOL-catalyzed asymmetric hetero-Diels–Alder reaction

The reverse-docking of a TADDOL catalyst to rigid transition-state (TS) representations of an asymmetric hetero-Diels–Alder reaction is described. The resulting docking poses represent a simplified geometric model of the TS for the catalyzed reaction. The conformational space of the catalyst in proximity to the catalyst-free TS models is sampled stochastically and the energetically favored poses are subjected to a clustering procedure to highlight structural attributes compatible with organocatalysis. Each pose is scored and ranked based on its molecular-mechanics docking energy. The reverse-docking procedure reveals a clear energetic trend in favor of the experimentally preferred product enantiomers. Analysis of the best poses suggests a geometric model that is consistent with principles of molecular recognition, catalysis, and experimental data.      

Theoretical study on the ground state intramolecular proton transfer (IPT) and solvation effect in two Schiff bases formed by 2-aminopyridine with 2-hydroxy-1- naphthaldehyde and 2-hydroxy salicylaldehyde

The tautomerization mechanism the isolated and monohydrated forms of two Schiff bases 1 and 2, and the effect of solvation on the proton transfer from enol-imine form to the keto-enamine form have been investigated using the B3LYP hybrid density functional method at the 6-31G** basis set level. The barrier heights for H₂O-assisted reactions are significantly lower than that of unassisted tautomerization reaction in the gas phase. Nonspecific solvent effects have also been taken into account by using the continuum model (IPCM) of four different solvent. The tautomerization energies and the potential energy barriers are decreased by increasing solvent polarity. Figure The tautomerization mechanism the isolated and monohydrated forms of two Schiff bases 1 and 2, and the effect of solvation on the proton transfer from enol-imine form to the keto-enamine form have been investigated using the B3LYP hybrid density functional method at the 6-31G** basis set level     

Study of the influence of temperature on the dynamics of the catalytic cleft in 1,3-1,4-β-glucanase by molecular dynamics simulations

The dependence of some molecular motions in the enzyme 1,3-1,4-β-glucanase from Bacillus licheniformis on temperature changes and the role of the calcium ion in them were explored. For this purpose, two molecular dynamics simulated trajectories along 4 ns at low (300 K) and high (325 K) temperatures were generated by the GROMOS96 package. Several structural and thermodynamic parameters were calculated, including entropy values, solvation energies, and essential dynamics (ED). In addition, thermoinactivation experiments to study the influence of the calcium ion and some residues on the activity were conducted. The results showed the release of the calcium ion, which, in turn, significantly affected the movements of loops 1, 2, and 3, as shown by essential dynamics. These movements differ at low and high temperatures and affect dramatically the activity of the enzyme, as observed by thermoinactivation studies. The first two authors contributed equally to this work     

The decomposition of α-aminophosphine oxides to phosphonic acid derivatives (P<Superscript>III</Superscript>)

Aminophosphine oxides and aminophosphonates are, in general, very stable compounds. However, following phosphorus–carbon bond cleavage in aqueous acidic media these compounds sometimes decompose to phosphonic acids derivatives (P^III). Despite some controversy in the literature, careful analysis supported by theoretical studies leads to the conclusion that decomposition to P^III derivatives proceeds via an elimination reaction. Figure The decomposition of α-aminophosphine oxides to phosphonic acid derivatives (P^III)     

Structural model of an antistasin/notch-like fusion protein from the cocoon wall of the aquatic leech, Theromyzon tessulatum

The aquatic leech, Theromyzon tessulatum, secretes a proteinaceous cocoon with extraordinary physical properties (e.g., proteolytic, thermal resiliency). The deduced amino acid sequence of a major protein (Tcp—Theromyzon cocoon protein) from the T. tessulatum cocoon wall has been used to model the endogenous structure of the Tcp protein. The Tcp protein sequence comprises six internal repeats, each containing 12 ordered Cys residues. Amino acid alignments suggest that the region Cys1→6 is homologous to antistasin, a leech anticoagulant, and Cys7→12 is homologous to an epidermal growth factor-like domain found in notch-class proteins, which play critical roles in development, signaling, and adhesion throughout the Animalia. Modeling of individual domains (i.e., antistasin and notch) positions multiple hydrophobic and charged residues on the surface. When the antistasin and notch domains were fused, hydrophobic pockets appeared that may facilitate a polymerization mechanism. Collectively, the predicted features of our Tcp model are consistent with the physical properties of the leech cocoon wall.      

SOMMER: self-organising maps for education and research

SOMMER is a publicly available, Java-based toolbox for training and visualizing two- and three-dimensional unsupervised self-organizing maps (SOMs). Various map topologies are implemented for planar rectangular, toroidal, cubic-surface and spherical projections. The software allows for visualization of the training process, which has been shown to be particularly valuable for teaching purposes. Spread of a spherical self-organizing map (SOM) in a three-dimensional data space     

Structure-based 3D-QSAR studies on thiazoles as 5-HT<Subscript>3</Subscript> receptor antagonists

Structure-based 3D-QSAR studies were performed on 20 thiazoles against their binding affinities to the 5-HT₃ receptor with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The thiazoles were initially docked into the binding pocket of a human 5-HT_3A receptor homology model, constructed on the basis of the crystal structure of the snail acetylcholine binding protein (AChBP), using the GOLD program. The docked conformations were then extracted and used to build the 3D-QSAR models, with cross-validated values 0.785 and 0.744 for CoMFA and CoMSIA, respectively. An additional five molecules were used to validate the models further, giving satisfactory predictive values of 0.582 and 0.804 for CoMFA and CoMSIA, respectively. The results would be helpful for the discovery of new potent and selective 5-HT₃ receptor antagonists.      

Activation energies of selenoxide elimination from Se-substituted selenocysteine

Transition states for selenoxide elimination have been determined for a series of Se-substituted selenocysteine (RSeCys) derivatives that have potential use in the prevention and treatment of cancer, either directly or in conjunction with cisplatin (to reduce its nephrotoxic effects). Reduced activation barriers vs R=Me and R=Ph are found when the alkyl chain length is increased or when activating groups are para to the selenide. Ortho substitution of Lewis bases stabilizes the transition state by directly donating electron density to the selenoxide. The results suggest that RSeCys derivatives incorporating the properties of glutathione peroxidase mimics will, upon oxidation, rapidly eliminate selenenic acid, a precursor to chemopreventative selenols. Scheme Mechanism of selenoxide elimination from Se-substituted SeCys.     

Topochemical models for the prediction of voltage-gated sodium channel binding activity of hydantoins and related non-hydantoins

The relationship of Wiener’s topochemical index—a distance based topochemical index, molecular connectivity topochemical index—an adjacency based topochemical index and eccentric connectivity topochemical index—an adjacency-cum-distance based topochemical index with sodium channel binding activity has been investigated. A dataset comprising 50 hydantoins and related non-hydantoins was selected. The dataset was divided equally into training and test sets. The values of the three topochemical indices for all the compounds present in both the training and test sets were computed using an in-house computer program. The resulting data was analyzed and suitable models were developed after identification of the active ranges in the training set. Subsequently, a biological activity was assigned to each compound involved in the training set using these models, which was then compared with the reported sodium channel binding activity. An accuracy of prediction of the order of >99% was observed using the proposed models. Cross-validation of these models using the test set revealed an exceptionally high accuracy of ∼95%. 3,5-disubstituted-5-phenylhydantoin     

Ion pair aggregates and reactions; experiment and theory

Two methods have been developed for determining the aggregation equilibrium constants of lithium enolates based on the change in UV-vis spectrum with concentration and the effect of aggregation on proton-transfer equilibria. Dimers and tetramers are common. Substitution α to the carbonyl group generally reduces aggregation. Kinetic studies show that S_N2 alkylations generally occur with the monomers, even in the presence of large amounts of aggregate. The qualitative chemistry is well modeled by ab initio computations at the HF level with modest basis sets; in these studies solvation is modeled by a combination of coordination of lithium cation with an ether oxygen and the electrostatic interaction of the resulting dipoles and quadrupoles with the solvent dielectric continuum. This review is based on the fifth annual Paul vR Schleyer Lecture presented at Professor Schleyer's 75th birthday symposium, Athens GA, Feb 26, 2005. Carbon Acidity 123. For no. 122 see Streitwieser A, Kaufman MJ, Bors DA, MacArthur CA, Murphy JT, Guibe F, Arkivoc (2005) (vi) 200–210Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid–liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. Figure Log–log plot of the self-diffusion coefficient D versus the chain length N. The error bars are the standard deviation measured in three repeated simulations     

A theoretical study on unusual intermolecular T-shaped X–H...π interactions between the singlet state HB=BH and HF, HCl, HCN or H2C2

The unusual T-shaped X–H...π hydrogen bonds are found between the B=B double bond of the singlet state HB=BH and the acid hydrogen of HF, HCl, HCN and H₂C₂ using MP2 and B3LYP methods at 6-311++G(2df,2p) and aug-cc-pVTZ levels. The binding energies follow the order of HB=BH...HF>HB=BH...HCl>HB=BH...HCN>HB=BH...H₂C₂. The hydrogen-bonded interactions in HB=BH...HX are found to be stronger than those in H₂C=CH₂...HX and OCB≡BCO...HX. The analyses of natural bond orbital (NBO) and the electron density shifts reveal that the nature of the T-shaped X–H...π hydrogen-bonded interaction is that much of the lost density from the π-orbital of B=B bond is shifted toward the hydrogen atom of the proton donor, leading to the electron density accumulation and the formation of the hydrogen bond. The atoms in molecules (AIM) theory have also been applied to characterize bond critical points and confirm that the B=B double bond can be a potential proton acceptor. The unusual T-shaped X–H...π hydrogen bonds are found between the B=B double bond of the singlet state HB=BH and the acid hydrogen of HF, HCl, HCN and H₂C₂     

An ab initio quantum mechanical drug designing procedure: application to the design of balanced dual ACE/NEP inhibitors

This article describes in a sequential fashion how ab initio quantum mechanical methods can be applied to study the pharmacophoric features of drugs. It also describes how accurate drug–receptor interaction calculations can guide the careful design of balanced dual inhibitors, which form an important class of second generation drugs. As an example, the authors have chosen balanced inhibitors of angiotensin converting enzyme/neutral endopeptidase (ACE/NEP) as modern antihypertensive drugs. A unified, accurate, in silico design approach is presented, encompassing all steps from pharmacophore derivation to complete understanding of mechanistic aspects leading to drug design.      

Theoretical 49Ti NMR chemical shifts

⁴⁹Ti chemical shifts for a total of 20 titanium complexes are reported, and several levels of theory are evaluated in order to identify a reliable approach for the calculation of titanium NMR data. The popular B3LYP/6–31G(d)//B3LYP/6–31G(d) proves to give very good agreement with experimental data over a range from 1,400 to −1,300 ppm. The MP2/6–31G(d)//MP2/6–31G(d) level computes even smaller average deviations but fails for TiI₄. This behavior together with its huge demand for computational resources requires careful handling of this theoretical level. In addition, NMR data for five titanium fulvene (or related) complexes are given. Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Conformational stability and three-dimensional model of the δ-opioid pharmacophore for the extended antiparallel dimer structure of Met-enkephalin in water

The conformational stability of the extended antiparallel dimer structure of Met-enkephalin in water was analyzed by examining the hydration structure of enkephalin using molecular dynamics simulations. The result shows that, despite of the hydrophicility of the terminal atoms in the pentapeptide, the main contributor for the stability of the dimer in water is the four intermolecular hydrogen bonds between the Gly² and Phe⁴ groups. The three-dimensional model of the δ-opioid pharmacophore for this dimer structure was also established. Such a model was demonstrated to match the δ-opioid pharmacophore query derived from the non-peptides SIOM, TAN-67, and OMI perfectly. This result thus strongly supports the assumption that the dimer structure of Met-enkephalin is a possible δ-receptor binding conformation. Figure Schematic model of the extended antiparallel dimer structure of Met-enkephalin     

DFT study of the gas phase proton transfer in guanine assisted by water, methanol, and hydrogen peroxide

A computational study of hydrogen-bonded complexes between the oxo-/hydroxo-amino N7/9H tautomers of guanine and water, methanol, and hydrogen peroxide has been performed at the B3LYP/6-31+G(d) level of theory. The mechanisms of the water-, methanol-, and hydrogen peroxide-assisted proton transfers in guanine were studied and compared with the intramolecular proton transfer in guanine in the gas phase. It was found that the assisted proton transfers pass through about three times lower energy barriers than those found for isolated guanine tautomers. Figure DFT study of the gas phase proton transfer in guanine assisted by water, methanol and hydrogen peroxide