Structural changes by sulfoxidation of phenothiazine drugs

Summary The side-chain conformations of psychoactive phenothiazine drugs in crystals are different from those of biologically inactive ring sulfoxide metabolites. This study examines the potential energies, molecular conformations and electrostatic potentials in chlorpromazine, levomepromazine (methotrimeprazine), their sulfoxide metabolites and methoxypromazine. The purpose of the study was to examine the significance of the different crystal conformations of active and inactive phenothiazine derivatives, and to determine why phenothiazine drugs lose most of their biological activity by sulfoxidation. Quantum mechanics and molecular mechanics calculations demonstrated that conformations with the side chain folded over the ring structure had lowest potential energy in vacuo, both in the drugs and in the sulfoxide metabolites. In the sulfoxides, side chain conformations corresponding to the crystal structure of chlorpromazine sulfoxide were characterized by stronger negative electrostatic potentials around the ring system than in the parent drugs. This may weaken the electrostatic interaction of sulfoxide metabolites with negatively charged domains in dopamine receptors, and cause the sulfoxides to be virtually inactive in dopamine receptor binding and related pharmacological tests.     

Ab initio MO-SCF calculation on a model of anomalous water

Ab initio calculations are carried out on the planar hexamer model of anomalous water using a variety of basis sets. The results are compared to results of similar calculations on free water and water in the ice I-like puckered ring configurations.

---

Zusammenfassung Es werden ab initio-Rechnungen für das planare Hexamer-Modell von anomalem Wasser unter Verwendung mehrerer BasissÄtze durchgeführt. Die Ergebnisse werden mit Resultaten Ähnlicher Rechnungen über freies Wasser und Wasser in der I-Eis Ähnlichen Konfiguration verzerrter Ringe verglichen.

---

Résumé Calculs ab initio sur le modèle hexamérique plan de l'eau anormale en utilisant différentes bases. Les résultats sont comparés à ceux de calculs similaires sur l'eau libre et sur l'eau dans une configuration cyclique compacte analogue à celle de la glace I.

     

THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method

The Weighted Histogram Analysis Method (WHAM), an extension of Ferrenberg and Swendsen's Multiple Histogram Technique, has been applied for the first time on complex biomolecular Hamiltonians. The method is presented here as an extension of the Umbrella Sampling method for free-energy and Potential of Mean Force calculations. This algorithm possesses the following advantages over methods that are currently employed: (1) It provides a built-in estimate of sampling errors thereby yielding objective estimates of the optimal location and length of additional simulations needed to achieve a desired level of precision; (2) it yields the “best” value of free energies by taking into account all the simulations so as to minimize the statistical errors; (3) in addition to optimizing the links between simulations, it also allows multiple overlaps of probability distributions for obtaining better estimates of the free-energy differences. By recasting the Ferrenberg–Swendsen Multiple Histogram equations in a form suitable for molecular mechanics type Hamiltonians, we have demonstrated the feasibility and robustness of this method by applying it to a test problem of the generation of the Potential of Mean Force profile of the pseudorotation phase angle of the sugar ring in deoxyadenosine. © 1992 by John Wiley & Sons, Inc.     

An ab initio molecular orbital study of intramolecular H-bonding: 1,3 propanediol

An ab initio molecular orbital calculation has been carried out for three different conformations of 1,3 propanediol, one of which permits intramolecular H-bond studied by ab initio quantum mechanical methods. The ΔE for H-bonding formation is compated to be 0.9 kcal/mole and the charge redistributions and molecular orbital energy changes are compared to those found in intermolecular H-bonds.     

Free energy perturbation calculations on models of active sites: Applications to adenosine deaminase inhibitors

Free energy perturbation calculations were performed to determine the free energy of binding associated with the presence of perhaps an unusual hydroxyl group in the transition state analog of nebularine, an inhibitor of the enzyme adenosine deaminase. The presence of a single hydroxyl group in this inhibitor has been found to contribute ?9.8 kcal/mol to the free energy of binding, with a 10⁸-fold increase in the binding affinity by the enzyme. In this work, we calculate the difference in solvation free energy for the 1,6-dihydropurine complex versus that of the 6-hydroxyl-1,6-dihydropurine complex to determine if this marked increase in binding affinity is attributed to an unusually hydrophobic hydroxyl group. The calculated ΔG associated for the solvation free energy is ?11.8 kcal/mol. This large change in the solvation free energy suggests that this hydroxyl is instead unusually hydrophilic and that the difference in free energy of interaction for the two inhibitors to the enzyme must be at least ca. 20 kcal/mol. Although the crystal structure for adenosine deaminase is currently not known, we attempt to mimic the nature of the active site by constructing models which simulate the enzyme-inhibitor complex. We present a first attempt at determining the change in free energy of binding for a system in which structural data for the enzyme is incomplete. To do this, we construct what we believe is a minimal model of the binding between adenosine deaminase and an inhibitor. The active site is simulated as a single charged carboxyl group which can form a hydrogen bond with the hydroxyl group of the analog. Two different carboxyl anion models are used. In the first model, the association is modeled between an acetic acid anion and the modified inhibitor. The second model consists of a hydrophobic amino acid pocket with an interior Glu residue in the active site. From these models we calculate the change in free energy of association and the overall change in free energy of binding. We calculate the free energies of interaction both in the absence and presence of water. We conclude from this that the presence of a single suitably placed-CO^?₂ group probably cannot explain the binding effect of the-OH group and that additional interactions will be found in the adenosine deaminase active site.     

The structure of the H3O+ (hydronium) ion

Near Hartree-Fock level ab initio molecular orbital calculations on H₃O⁺ and a minimum energy structure with θ(HOH) = 112.5° and r(OH) = 0.963 Å and an inversion barrier of 1.9 kcal/mole. By comparing these results to calculations on NH₃ and H₂O, where precise experimental geometries are known, we estimate the “true” geometry of isolated H₃O⁺ to have a structure with θ(HOH) = 110-112°, r(OH) = 0.97–0.98 Å and an inversion barrier of 2–3 kcal/mole. Our prediction for the proton affinity of water is ≈ 170 kcal/mole, which is somewhat smaller than the currently accepted value.     

Synthesis of periphery-functionalized dendritic molecules using polylithiated dendrimers as starting material

A general method for the functionalization of Si-Cl terminated carbosilane dendritic molecules via organolithium or organomagnesium reagents is described. Quantitative exchange of the bromine atoms of 4-bromophenyl-functionalized dendrimers affords polylithiated species that are valuable starting materials for further functionalization, e.g., into pyridyl alcohols. The latter were successfully applied as catalyst precursors in a ruthenium-mediated ring-closure metathesis reaction.     

Oscillatory melting temperature of the vortex smectic phase in layered superconductors

We report on transport measurements of YBa 2Cu 3O (7-delta) single crystals with different oxygen contents in the geometry B, J ||ab (J perpendicularB). Our data show that the vortices become confined between the Cu-O planes below a well-defined temperature at which the effective size 2xi of the vortex core is approximately equal to the period of the Cu-O layers. This confinement strongly increases the vortex liquid freezing temperature. A new melting line is found separating a vortex liquid and a smectic phase, which shows an oscillatory field dependence reflecting differences between commensurate and incommensurate smectic states.