Molecular mechanics calculations on some phosphine oxide metal complexes in aqueous and organic solvents

Molecular mechanics calculations were done on tetrahedral phosphine oxide zinc complexes in simulated water, benzene and hexane phases using the DREIDING II force field in the BIOGRAF molecular modeling program. The SUN workstation computer (SUN_ 4c, with SPARK station 1 processor) was used for the calculations. Experimental structural information used in the parameterization was obtained from the September 1989 version of the Cambridge Structural Database. 2 Steric and solvation energies were calculated for complexes of the type ZnCl2 (RlO)2' The calculations were done with and without inclusion of electrostatic interactions. More reliable simulation results were obtained without inclusion of charges. In the simulated gas phase, the steric energies increase regularly with number of carbons in the alkyl group, whereas they go through a maximum when solvent shells are included in the calculation. Simulated distribution ratios vary with chain length and type of chain branching and the complexes are found to be more favourable for extraction by benzene than by hexane, in accord with experimental data. Also, in line with what would be expected for a favorable extraction, calculations without electrostatics predict that the complexes are better solvated by the organic solvents than by water.

Assessing the Reproducibility of Clustering of Molecular Dynamics Conformations on Self-Organizing Maps

The goal of most clustering algorithms is to find the optimal number of clusters (i.e. fewest number of clusters). However, analysis of molecular conformations of biological macromolecules obtained from computer simulations may benefit from a larger array of clusters. The Self-Organizing Map (SOM) clustering method has the advantage of generating large numbers of clusters, but often gives ambiguous results. In this work, SOMs have been shown to be reproducible when the same conformational dataset is independently clustered multiple times (~100), with the help of the Cramérs V-index (C_v). The ability of C_v to determine which SOMs are reproduced is generalizable across different SOM source codes. The conformational ensembles produced from MD (molecular dynamics) and REMD (replica exchange molecular dynamics) simulations of the penta peptide Met-enkephalin (MET) and the 34 amino acid protein human Parathyroid Hormone (hPTH) were used to evaluate SOM reproducibility. The training length for the SOM has a huge impact on the reproducibility. Analysis of MET conformational data definitively determined that toroidal SOMs cluster data better than bordered maps due to the fact that toroidal maps do not have an edge effect. For the source code from MATLAB, it was determined that the learning rate function should be LINEAR with an initial learning rate factor of 0.05 and the SOM should be trained by a sequential algorithm. The trained SOMs can be used as a supervised classification for another dataset. The toroidal 10×10 hexagonal SOMs produced from the MATLAB program for hPTH conformational data produced three sets of reproducible clusters (27%, 15%, and 13% of 100 independent runs) which find similar partitionings to those of smaller 6×6 SOMs. The χ^2 values produced as part of the C_v calculation were used to locate clusters with identical conformational memberships on independently trained SOMs, even those with different dimensions. The χ^2 values could relate the different SOM partitionings to each other.

Bounds on edit metric codes with combinatorial DNA constraints

The design of a large and reliable DNA codeword library is a key problem in DNA based computing. DNA codes, namely sets of fixed length edit metric codewords over the alphabet {A, C, G, T}, satisfy certain combinatorial constraints with respect to biological and chemical restrictions of DNA strands. The primary constraints that we consider are the reverse--complement constraint and the fixed GC--content constraint, as well as the basic edit distance constraint between codewords. We focus on exploring the theory underlying DNA codes and discuss several approaches to searching for optimal DNA codes. We use Conway's lexicode algorithm and an exhaustive search algorithm to produce provably optimal DNA codes for codes with small parameter values. And a genetic algorithm is proposed to search for some sub--optimal DNA codes with relatively large parameter values, where we can consider their sizes as reasonable lower bounds of DNA codes. Furthermore, we provide tables of bounds on sizes of DNA codes with length from 1 to 9 and minimum distance from 1 to 9.

A comparative study of protoheme and heme d catalases : role of the heme and the heme pocket in catalysis and ligand binding

Catalase dismutes H20 2 to O2 and H20. In successive twoelectron reactions H20 2 induces both oxidation and reduction at the heme group. In the first step the protoheme prosthetic group of beef liver catalase forms compound I, in which the heme has been oxidized from Fe3+ to Fe4+=0 and a porphyrin radical has been created. Compound II is formed by the oneelectron reduction of comp I. It retains Fe4+=0 but lacks the porphyrin radical and is catalytically inert. Molecular structures are available for Escherichia coli Hydroperoxidase II, Micrococcus Iysodeiktus, Penicillium vitale and beef liver enzymes, which contain different hemes and heme pockets. In the present work, the pockets and substrate access channels of protoheme (beef liver & Micrococcus) and heme d (HPII of E. coli and Penicillium) catalases have been analysed using Quanta™ and CharmMTM molecular modeling packages on the Silicon Graphics Iris Indigo 2 computer. Experimental studies have been carried out with two catalases, HPII (and its mutants) and beef liver. Fluoride and formate' are inhibitors of both enzymes, and their binding is modulated by the heme and by distal residues N201 & H128. Both HPII and beef liver enzymes form compound I with H202 or peracetate. The reduction of beef liver enzyme compound I to II and the decay of compound II are accelerated by fluoride. The decay of compound II is also accelerated by formate, and this reagent acts as a 2-electron donor towards compound I of both enzymes. It is concluded that heme d enzymes (Penicillium and HPII of E. coli) are formed by autocatalytic transformation of protoheme in a modified pocket which contains a characteristic serine residue as well as a partially occluded heme channel. They are less active than protoheme enzymes but also do not form the inactive compound II species. Binding of peroxide as well as fluoride and formate is prevented by mutation of H128 and modulated by mutation of N201.

Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Investigation of the composition of linear alkylbenzenes with emphasis on the identification and quantitation of some trace compounds using GS/MS system in both electron impact and chemical ionization modes

Linear alkylbenzenes, LAB, formed by the Alel3 or HF catalyzed alkylation of benzene are common raw materials for surfactant manufacture. Normally they are sulphonated using S03 or oleum to give the corresponding linear alkylbenzene sulphonates In >95 % yield. As concern has grown about the environmental impact of surfactants,' questions have been raised about the trace levels of unreacted raw materials, linear alkylbenzenes and minor impurities present in them. With the advent of modem analytical instruments and techniques, namely GCIMS, the opportunity has arisen to identify the exact nature of these impurities and to determine the actual levels of them present in the commercial linear ,alkylbenzenes. The object of the proposed study was to separate, identify and quantify major and minor components (1-10%) in commercial linear alkylbenzenes. The focus of this study was on the structure elucidation and determination of impurities and on the qualitative determination of them in all analyzed linear alkylbenzene samples. A gas chromatography/mass spectrometry, (GCIMS) study was performed o~ five samples from the same manufacturer (different production dates) and then it was followed by the analyses of ten commercial linear alkylbenzenes from four different suppliers. All the major components, namely linear alkylbenzene isomers, followed the same elution pattern with the 2-phenyl isomer eluting last. The individual isomers were identified by interpretation of their electron impact and chemical ionization mass spectra. The percent isomer distribution was found to be different from sample to sample. Average molecular weights were calculated using two methods, GC and GCIMS, and compared with the results reported on the Certificate of Analyses (C.O.A.) provided by the manufacturers of commercial linear alkylbenzenes. The GC results in most cases agreed with the reported values, whereas GC/MS results were significantly lower, between 0.41 and 3.29 amu. The minor components, impurities such as branched alkylbenzenes and dialkyltetralins eluted according to their molecular weights. Their fragmentation patterns were studied using electron impact ionization mode and their molecular weight ions confirmed by a 'soft ionization technique', chemical ionization. The level of impurities present i~ the analyzed commercial linear alkylbenzenes was expressed as the percent of the total sample weight, as well as, in mg/g. The percent of impurities was observed to vary between 4.5 % and 16.8 % with the highest being in sample "I". Quantitation (mg/g) of impurities such as branched alkylbenzenes and dialkyltetralins was done using cis/trans-l,4,6,7-tetramethyltetralin as an internal standard. Samples were analyzed using .GC/MS system operating under full scan and single ion monitoring data acquisition modes. The latter data acquisition mode, which offers higher sensitivity, was used to analyze all samples under investigation for presence of linear dialkyltetralins. Dialkyltetralins were reported quantitatively, whereas branched alkylbenzenes were reported semi-qualitatively. The GC/MS method that was developed during the course of this study allowed identification of some other trace impurities present in commercial LABs. Compounds such as non-linear dialkyltetralins, dialkylindanes, diphenylalkanes and alkylnaphthalenes were identified but their detailed structure elucidation and the quantitation was beyond the scope of this study. However, further investigation of these compounds will be the subject of a future study.

On the path integral formulation and the evaluation of quantum statistical averages

Four problems of physical interest have been solved in this thesis using the path integral formalism. Using the trigonometric expansion method of Burton and de Borde (1955), we found the kernel for two interacting one dimensional oscillators• The result is the same as one would obtain using a normal coordinate transformation, We next introduced the method of Papadopolous (1969), which is a systematic perturbation type method specifically geared to finding the partition function Z, or equivalently, the Helmholtz free energy F, of a system of interacting oscillators. We applied this method to the next three problems considered• First, by summing the perturbation expansion, we found F for a system of N interacting Einstein oscillators^ The result obtained is the same as the usual result obtained by Shukla and Muller (1972) • Next, we found F to 0(Xi)f where A is the usual Tan Hove ordering parameter* The results obtained are the same as those of Shukla and Oowley (1971), who have used a diagrammatic procedure, and did the necessary sums in Fourier space* We performed the work in temperature space• Finally, slightly modifying the method of Papadopolous, we found the finite temperature expressions for the Debyecaller factor in Bravais lattices, to 0(AZ) and u(/K/ j,where K is the scattering vector* The high temperature limit of the expressions obtained here, are in complete agreement with the classical results of Maradudin and Flinn (1963) .

Single-Molecule Magnets and Multifunctional Molecular Magnetic Materials Based on Polynuclear Metal Complexes

Our work on single molecule magnets and multifunctional magnetic materials is presented in four projects. In the first project we show for first time that heteroatomic-type pseudohalides, such as OCN-, can be employed as structure-directing ligands and ferromagnetic couplers in higher oxidation state metal cluster chemistry. The initial use of cyanato groups in Mn cluster chemistry has afforded structurally interesting MnII/III14 (1) and MnII/III/IV16 (2) clusters in which the end-on bridging cyanates show a preference in binding through their O-atom. The Mn14 compound shows entirely visible out-of-phase alternating currect signals below 5 K and large hysteresis loops below 2 K. Furthermore, the amalgamation of azido groups with the triethanolamine tripodal ligand in manganese carboxylate cluster chemistry has led to the isolation of a new ferromagnetic, high-nuclearity and mixed-valence MnII/III15Na2 (3) cluster with a large ground-state spin value of S = 14. In the second project we demonstrate a new synthetic route to purely inorganic-bridged, transition metal-azido clusters [CoII7 (4) and NiII7 (5)] and coordination polymers [{FeII/III2}n (6)] which exhibit strong ferromagnetic, SMM and long-range magnetic ordering behaviors. We also show that access to such a unique ferromagnetic class of inorganic, N-rich and O-free materials is feasible through the use of Me3SiN3 as the azido-ligand precursor without requiring the addition of any organic chelating/bridging ligand. In the last projects we have tried to bring together molecular magnetism and optics via the synthesis of multifunctional magnetic materials based on 3d- or 4f-metal ions. We decided to approach such challenge from two different directions: firstly, in our third project, by the deliberate replacement of non-emissive carboxylato ligands in known 3d-SMMs with their fluorescent analogues, without perturbing the metal-core structure and SMM properties (complexes 7, 8, and 9). The second route (last project) involves the use of naphthalene or pyridine-based polyalcohol bridging ligands for the synthesis of new polynuclear LnIII metal clusters (Ln = lanthanide) with novel topologies, SMM behaviors and luminescent properties arising from the increased efficiency of the “antenna” organic group. This approach has led us to the isolation of two new families of LnIII8 (complexes 10-13) and LnIII4 (complexes 14-20) clusters.