Molecular mechanics modelling of Pt(II) complexes with antitumor activity. Influence of the type and the positions of the ring substituents on the conformational energies and thermodynamic stabilities

Summary The influence of the type and the positions of the ring substituents on the conformational energies and thermodynamic stabilities of a series of Pt(II) complexes of the general formula [1,2-bis(hydroxyphenyl)ethylenediamine]PtL ₂(L ₂=2Cl^–, 2I^–, SO ₄ ^2– ) has been studied by molecular mechanics. The calculations were carried out for the ligand conformations (R,S/S,R)-, (R,S/S,R)-, (R,R)-, and (S,S)-. The obtained energies and thermodynamic stabilities are in agreement with experimental data on the reactivity and antitumor activity of the compounds.

---

Molecular Modelling von Pt(II)-Komplexen mit Antitumoraktivität. Einfluß von Art und Position von Substituenten auf Konformationsenergien und thermodynamische Stabilitäten

Zusammenfassung Die konformationellen Energien und thermodynamischen Stabilitäten einer Reihe von Pt(II)-Komplexen mit der allgemeinen Formel [1,2-bis(Hydroxyphenyl) ethylendiamin]PtL ₂(L ₂=2Cl^–, 2I^–, SO ₄ ^2– ) wurden mittels molekularmechanischer Methoden in Abhängigkeit von Art und Stellung der Substituenten an den Phenylringen untersucht. Die Berechnungen wurden für die Ligandenkonformationen (R,S/S,R)-, (R,S/S,R)-, (R,R)- und (S,S)- durchgeführt. Die erhaltenen Energien und Stabilitäten stimmen mit experimentellen Daten über Reaktivität und Antitumoraktivität der Verbindungen überein.

     

Cation π interaction between acetylcholine and the benzene ring

 Ab initio self-consistent-field second-order M?ller–Plesset perturbation theory computations including basis set superposition error and zero-point vibrational energy corrections have been performed on the complexation of benzene with the polar head of acetylcholine (ACh). The ACh–benzene complex is about 0.5 kcal/mol less stable than the corresponding tetramethylammonium (TMA)–benzene complex, with a structure a little distorted with respect to the latter. The electronic structure of ACh is little modified by the ligand. Overall, the replacement of one methyl group of TMA by the acetyl tail of ACh does not affect strongly the complexation to benzene, as far as the main interaction is concerned. Received: 1 April 1999 / Accepted: 19 October 1999 / Published online: 14 March 2000     

De novo prediction of the ground state structure of transition metal complexes using semiempirical and ab initio quantum mechanics. Coordination isomerism

The ground state coordination isomers for 30 different trigonal bipyramidal transition metal complexes have been predicted using different levels of quantum mechanics: semiempirical (PM3(tm)), ab initio (MP2//HF), pure (BPW91) and hybrid (B3PW91) density functional theory (DFT) methods. For species where these methods failed to reproduce crystallographic data, hybrid quantum mechanics/molecular mechanics (QM/MM) methods were used to study more exact experimental models. Literature deficiencies regarding ground state multiplicity of these species were supplemented by spin predictions using previously tested PM3(tm) methods. Geometry optimization calculations were performed for each possible coordination isomer. The predicted ground state minima provided by the different methods are compared to each other and with crystallographic data. Pure DFT functionals outperformed hybrid functionals and MP2//HF. The very rapid PM3(tm) parameterization method provided accurate predictions in comparison to other levels of theory. An integrated MM/PM3(tm)/DFT de novo scheme accurately reproduced crystallographic data for species where the individual methods failed.     

Study of amino- and pyridyl-containing zinc(II) and cadmium(II) complexes by molecular mechanics using a force field based on the Gillespie-Kepert model

A new approach combining the molecular mechanics (MM) method and the Gillespie-Kepert model was applied to calculate the geometry and strain energy of zinc(II) and cadmium(II) complexes with amino- and pyridyl-containing ligands. High accuracy of calculations of the geometry was demonstrated for more than 20 complexes of these metals. Typical r.m.s. deviations between the calculated and experimental values (X-ray diffraction analysis) were 0.02 Å for bond lengths, 2° for bond angles, and 4° for torsion angles. The size-match selectivity of several macrocycles and polydentate open-chain ligands was studied. Correlations between the calculated strain energies of metal complexes and the experimental values of their stability constants and enthalpies of formation are discussed.     

Theoretical study of the reaction of CN with C2H2+

 A theoretical study of the reaction of CN with C₂H₂ ⁺ has been carried out at three levels of theory, namely G2, B3LYP and CCSD(T). The main conclusion is that this is a feasible process under interstellar conditions, but only linear species may be produced. The most favourable product is HCCCN⁺, followed by CCCNH⁺. Production of HCCNC⁺ is predicted to be slightly endothermic; therefore, the reaction of CN + C₂H₂ ⁺ may produce precursors of HC₃N and C₃N in space. Furthermore, the B3LYP level is found to perform rather well compared with G2 and even better than CCSD(T). Received: 14 September 1999 / Accepted: 3 February 2000 / Published online: 12 May 2000     

Theoretical study of the structural character of weakly bonding silicon carbonyl complexes

The structures, properties and the bonding character for sub-carbonyl Si, SiCO and Si(CO)₂, in singlet and triplet states have been investigated using complete-active-space self-consistent field (CASSCF), density functional theory and second-order M?ller–Plesset methods with a 6-311+G* basis set. The results indicate that the SiCO species possesses a ³∑⁻ ground state, and the singlet ¹Δ excited state is higher in energy than the ³∑⁻ state by 17.3 kcalmol⁻¹ at the CASSCF–MP2/6-311+G* level and by 16.4 kcalmol⁻¹ at the CCSD(T)/6-311+G* level. The SiCO ground state may be classified as silene (carbonylsilene), and its CO^δ− moiety possesses CO⁻ property. The formation of SiCO causes the weakening of CO bonds. The Si–C bond consists of a weak σ bond and two weak π bonds. Although the Si–C bond length is similar to that of typical Si–C bonds, the bond strength is weaker than the Si–C bonds in Si-containing alkanes; the calculated dissociation energy is 26.2 kcalmol⁻¹ at the CCSD(T)/6-311+G* level. The corresponding bending potential-energy surface is flat; therefore, the SiCO molecule is facile. For the bicarbonyl Si systems, Si(CO)₂, there exist two V-type structures for both states. The stablest state is the singlet state (¹A₁), and may be referred to the ground state. The triplet state (³B₁) is energetically higher in energy than the ¹A₁ state by about 40 kcalmol⁻¹ at the CCSD(T)/6-311 + G* level. The bond lengths in the ¹A₁ state are very close to those of the SiCO species, but the SiCO moieties are bent by about 10°, and the CSiC angles are only about 78°. The corresponding ³B₁ state has a CSiC angle of about 54° and a SiCO angle of about 165°, but its Si–C and C–O bonds are longer than those in the ¹A₁ state by about 0.07 and 0.03 ?, respectively. This Si(CO)₂ (¹A₁) has essentially silene character and should be referred to as a bicarbonyl silene. Comparison of the CO dissociation energies of SiCO and Si(CO)₂ in their ground states indicates that the first CO dissociation energy of Si(CO)₂ is smaller by about 7 kcalmol⁻¹ than that of SiCO; the average one over both CO groups is also smaller than that of SiCO. A detailed bonding analysis shows that the possibility is small for the existence of polycarbonyl Si with more than three CO. This prediction may also be true for similar carbonyl complexes containing other nonmetal and non-transition-metal atoms or clusters. Received: 17 April 2002 / Accepted: 11 August 2002 / Published online: 4 November 2002 Acknowledgements. This work was supported by the National Natural Science Foundation of China (29973022) and the Foundation for Key Teachers in University of the State Ministry of Education of China. Correspondence to: Y. Bu e-mail: byx@sdu.edu.ch     

Syntheses, crystal structure and spectroscopic characterization of new platinum(II) dithiocarbimato complexes

The reaction of 4-iodobenzenesulfonamide or 4-fluorobenzenesulfonamide with CS₂ and KOH in dimethylformamide yielded the potassium N-R-sulfonyldithiocarbimates, K₂(RSO₂NCS₂) [R = 4-IC₆H₄ (1) and 4-FC₆H₄ (2)]. These salts reacted with K₂[PtCl₄] in water/methanol to yield complex anions bis(N-R-sulfonyldithiocarbimato)platinate(II), which were isolated as their tetrabutylammonium salts, (Bu₄N)₂[Pt(RSO₂NCS₂)₂] [R = 4-IC₆H₄ (3) and 4-FC₆H₄ (4)]. The structures of 2–4 were determined by X-ray crystallography. The Pt²⁺ in both complexes 3 and 4 lies at the inversion centre and the PtS₄ moiety has a distorted square-planar configuration. The compounds were also characterized by IR, ¹H NMR and ¹³C NMR spectroscopies, and elemental analyses. The molar conductance data are consistent with the fact that 3 and 4 are dianionic complexes.     

Genistein complexes with amines. Part II: ab initio study of the complexes with piperazine and triethylamine

The piperazine and triethylamine complexes of genistein, exhibiting high immunosuppressant activity, were ab initio modeled at RHF/6-31G^** level and results were compared with those obtained for genistein–morpholine complexes by X-ray, NMR, and theoretical methods. The most stable genistein–piperazine complex is formed due to hydrogen bonding of genistein's OH group at position C7 to piperazine's nitrogen atom. In the most stable genistein–triethylamine complex genistein's OH group at position C4′ (position para to phenyl substituent) and trimethylamine nitrogen atom are engaged in hydrogen bond formation. The calculations confirmed our previous NMR conclusion that piperazine is more strongly complexed by genistein than is morpholine or triethylamine. The theoretical ¹³C NMR spectra correlate fairly well with the experimental spectra.     

Synthesis and structural characterization of novel bridged platinum(II) biscarbene complexes

Novel bridged platinum(II) biscarbene complexes are reported: 1,1′-dimethyl-3,3′-methylene-4-diimidazolin-2,2′-diylidene platinum(II) (3) and 1,1′-dimethyl-3,3′-ethylene-4-diimidazolin-2,2′-diylidene platinum(II) complexes 4 are directly accessible in high yields starting from platinum halides. The one-pot synthesis obviates the need for multi-step reactions via metal precursors or free carbenes. An X-ray crystal structure of 1,1′-dimethyl-3,3′-methylene-4-diimidazolin-2,2′-diylidene platinum(II) dibromide (3b) confirmed the structural similarity to the known corresponding palladium complexes. Since free 1,1′-di-R-3,3′-methylene-4-diimidazolin-2,2′-diylidenes are only available in low yields this synthetic route provides an easy access to the corresponding carbene complexes.     

Synthesis and coordination properties of palladium(II) and platinum(II) complexes with phosphonated triphenylphosphine derivatives

Two triphenylphosphine derivatives, diethyl [4-(diphenylphosphanyl)benzyl]phosphonate (3a) and tetraethyl {[5-(diphenylphosphanyl)-1,3-phenylene]dimethylene}bis(phosphonate) (3b), and also the corresponding free acids 4a and 4b were prepared. These ligands were characterized by ¹H, ¹³C and ³¹P NMR spectroscopy and mass spectrometry. A full set of their Pd(II) and Pt(II) complexes of the general formula [MCl₂L₂] and one dinuclear complex trans-[Pd₂Cl₄(3a)₂] were synthesized and their isomerization behaviour in solution was studied. The complexes were characterized by ¹H, ¹³C, ³¹P and ¹⁹⁵Pt NMR spectroscopy, mass spectrometry and far-IR spectroscopy. The X-ray structures of all complexes with 3a or 3b have usual slightly distorted square-planar geometry on the metal ion. Salts of phosphonic acids 4a and 4b and their complexes are freely soluble in aqueous solution; therefore, they can be potentially useful in aqueous or biphasic catalysis.     

Rhodium(III), palladium(II) and platinum(II) complexes of Bis(o-aminobenzenesulfonyl)ethylenediamine

New complexes of the formulae K₃[RhL ₃]·2 H₂O, [PdL]·H₂O and [M(LH₂)Cl₂] [whereM = Pd, Pt andLH₂ = bis(o-aminobenzenesulfonyl)ethylenediamine] have been prepared and characterized by conductivity measurements, thermogravimetric analysis, X-ray powder patterns and IR, Ligand Field and¹H-NMR spectroscopy.

---

Rhodium(III), Palladium(II)- und Platin(II)-Komplexe mit Bis(o-aminobenzolosulfonyl)ethylendiamin (Kurze Mitteilung)

Zusammenfassung Neue Komplexe der allgemeinen Formeln K₃[RhL ₃]·2H₂O, [PdL]·H₂O und [M(LH₂)Cl₂] mitM = Pd, Pt undLH₂ = Bis(o-aminobenzolosulfonyl)ethylendiamin wurden dargestellt und mit Konduktionsmessungen, thermogravimetrischen Analysen, Röntgenstrukturanalysen, IR, Ligandfeld- und¹H-NMR-Spektroskopie charakterisiert.

     

X-ray crystal structures of novel platinum(II) and palladium(II) complexes of dialkyl phosphonated phosphines

Two diethyl phosphonated phosphine ligands of formula Ph₂P(CH₂)₃PO₃Et₂ (ligand L) and Ph₂P(4-C₆H₄PO₃Et₂) (ligand L′) were used to prepare different complexes of platinum(II) (1, cis-PtCl₂L₂; 2, trans-PtCl₂L₂·H₂O; 3A and 3B, cis- and trans-PtCl₂L′₂) and palladium(II) (4, [PdCl₂L]₂; 5, trans-PdCl₂L₂·H₂O; 6, trans-PdCl₂L′₂·CH₂Cl₂). The single-crystal X-ray structure analyses of complexes 1, 2, 4-6 indicate that complexation involved only the phosphine end, whereas the strong polarization of the PO bond was highlighted by the formation of hydrogen bonds with a water molecule in 2 and 5, and with a dichloromethane molecule in 6, with an exceptionally short CH?O hydrogen bond length (C?O separation 3.094(3) Å).     

Molecular recognition by β-cyclodextrin derivatives: molecular dynamics, free-energy perturbation and molecular mechanics/ Poisson–Boltzmann surface area goals and problems

 The complexation of p-tert-butylphenyl p-tert-butylbenzoate and N-(p-tert-butylphenyl)-p-tert-butylbenzamide with a β-cyclodextrin derivative formed by two cyclodextrin units linked by a disulfide bridge on one of the C6 atoms has been studied by computational methods. The better amide solubility and the better internal interactions of the ester complex explain the experimentally observed better association constant for the ester. The free-energy perturbation methodology and molecular mechanics/Poisson–Boltzmann surface area analysis have been used to explain the problem and to compare the results. Received: 14 April 2002 / Accepted: 11 August 2002 / Published online: 4 November 2002 Acknowledgements. The Kollman group at the University of California San Francisco is gratefully acknowledged for support and encouragement throughout all this study. The authors thank UAB for inland and outland fellowships to I.B.. Financial support was obtained from grant no. PPQ2000-0369 from the “Ministerio de Ciencia y Tecnologia” (Spain). Intensive computations were performed either with the computers of the Kollman group or with those of CESCA-C⁴ (Catalonia, Spain). Correspondence to: C. Jaime e-mail: carlos.jaime@uab.es     

Partial Hessian vibrational analysis: the localization of the molecular vibrational energy and entropy

The partial Hessian vibrational analysis (PHVA), in which only a subblock of the Hesssian matrix is diagonalized to yield vibrational frequencies for partially optimized systems, is extended to the calculation of vibrational enthalpy and entropy changes for chemical reactions. The utility of this method is demonstrated for various deprotonation reactions by reproducing full HVA values to within 0.1–0.4 kcal/mol, depending on the number atoms included in the PHVA. When combined with the hybrid effective fragment potential method [Gordon MS, et al. (2001) J Phys Chem A 105:293–307], the PHVA method can provide (harmonic) free-energy changes for localized chemical reactions in very large systems. Received: 21 September 2001 / Accepted: 30 October 2001 / Published online: 22 March 2002     

Model building and molecular mechanics calculations of mitoxantrone-deoxytetranucleotide complexes: Molecular foundations of DNA intercalation as cytostatic active principle

Summary Several intercalation complexes of the antitumor-active drug mitoxantrone with base paired tetranucleotides were constructed by molecular modeling using computer graphics and molecular mechanics calculations. The mitoxantrone molecule favours DNA binding into CG intercalation site. The two side chains of the drug are orientated into the major groove and fixed by hydrogen bonds with the nucleotide bases. This molecular study can be helpful for understanding the mode of action of cytostatically active compounds and to design new structurally related compounds of the anthraquinone drug type.Dedicated to Prof. Dr. Dr. h.c. mult.H. Oelschläger on the occasion of his 75^th birthday     

Synthesis and characterization of new platinum(II) phosphinate complexes

The syntheses of platinum(II) complexes of bis(dimethylphosphinylmethylene)amine and bis(aminomethyl)phosphinic acid were investigated. In the case of bis(dimethyl-phosphinylmethylene)amine the reaction with K₂[PtCl₄] yields the potassium amino-trichloroplatinate K[PtCl₃L] (L?=?bis(dimethylphosphinylmethylene)amine), which was characterized by multinuclear (¹H, ¹³C, ³¹P, and ¹⁹⁵Pt) NMR spectroscopy in solution. Bis(aminomethyl)phosphinic acid reacts with K₂[PtCl₄] under strictly controlled pH conditions to give colorless crystals of the cisplatin analog K[PtCl₂L′] (L′?=?bis(aminomethyl)phosphinate). This complex was characterized by multinuclear NMR spectroscopy in solution as well as by single-crystal X-ray diffraction in the solid state. The bis(aminomethyl)phosphinate coordinates to platinum via both amino functions, thus acting as a chelating ligand.     

Acetylene-insertion reactions into Pt(II)–H and Pt(II)–SiH3 bonds. An ab initio MO study and analysis based on the vibronic coupling model

Acetylene insertion into Pt(II)–H and Pt(II)SiH₃ bonds of PtH(SiH₃)(PH₃) was investigated using ab initio molecular orbital and M?ller-Plesset perturbation theory methods. The insertion into PtH was predicted to proceed with a smaller activation energy (E _a =12.8 kcal/mol) than that into PtSiH₃ (E _a =20.9 kcal/mol). The reaction energy (ΔE) of the insertion into PtH is 10 kcal/mol smaller than that for the insertion into PtSiH₃, which reflects differences in bond energies between CH and CSi and between PtH and PtSiH₃. A comparison with ethylene insertion revealed that the acetylene insertion occurs more easily, and the latter reaction is more exothermic. A simple vibronic coupling model combined with Toyozawa's interaction mode analysis was used to examine interesting differences in E _a and ΔE between insertions into PtH and PtSiH₃, and between acetylene and ethylene insertions. This analysis suggests that the factors determining E _a are the stiffness of the PtH and PtSiH₃ bonds and the vibronic coupling strength of acetylene and ethylene. Received: 13 August 1998 / Accepted: 2 September 1998 / Published online: 15 February 1999     

Modeling the H5+ potential-energy surface: a first attempt

 Ab initio molecular electronic structure calculations are performed for H₅ ⁺ at the QCISD(T) level of theory, using a correlation-consistent quadruple-zeta basis set. Structures, vibrational frequencies and thermochemical properties are evaluated for ten stationary points of the H₅ ⁺ hypersurface and are compared with previous calculations. The features of the H₃ ⁺…H₂ interaction at intermediate and large intermolecular distances are also investigated. Furthermore, an analytical functional form for the potential-energy surface of H₅ ⁺ is derived using a first-order diatomics-in-molecule perturbation theory approach. Its topology is found to be qualitatively correct for the short-range interaction region. Received: 15 March 2001 / Accepted: 5 July 2001 / Published online: 11 October 2001