5,6-Dihydroxyindole-2-carboxylic acid: a first principles density functional study

We report first principles density functional calculations for 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and several oxidised forms. DHICA and 5,6-dihydroxyindole (DHI) are believed to be the basic building blocks of the eumelanins. Our results show that carboxylation has a significant effect on the physical properties of the molecules. In particular, the relative stabilities and the highest occupied molecular orbital-lowest unoccupied molecular orbital gaps (calculated with the DeltaSCF method) of the various redox forms are strongly affected. We predict that, in contrast to DHI, the density of unpaired electrons, and hence the ESR signal, in DHICA is negligibly small. (C) 2004 Elsevier B.V. All rights reserved.

Why H Atom Prefers the On-Top Site and Alkali Metals Favor the Middle Hollow Site on the Basal Plane of Graphite

In this work, the different adsorption properties of H and alkali metal atoms on the basal plane of graphite are studied and compared using a density functional method on the same model chemistry level. The results show that H prefers the on-top site while alkali metals favor the middle hollow site of graphite basal plane due to the unique electronic structures of H, alkali metals, and graphite. H has a higher electronegativity than carbon, preferring to form a covalent bond with C atoms, whereas alkaline metals have lower electronegativity, tending to adsorb on the highest electrostatic potential sites. During adsorption, there are more charges transferred from alkali metal to graphite than from H to graphite.

Chemical and structural disorder in eumelanins: A possible explanation for broadband absorbance

We report the results of an experimental and theoretical study of the electronic and structural properties of a key eumelanin precursor-5,6,-dihydroxyindole-2-carboxylic acid ( DHICA) - and its dimeric forms. We have used optical spectroscopy to follow the oxidative polymerization of DHICA to eumelanin and observe red shifting and broadening of the absorption spectrum as the reaction proceeds. First principles density functional theory calculations indicate that DHICA oligomers ( possible reaction products of oxidative polymerization) have the gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital red-shifted gaps with respect to the monomer. Furthermore, different bonding configurations ( leading to oligomers with different structures) produce a range of gaps. These experimental and theoretical results lend support to the chemical disorder model where the broadband monotonic absorption characteristic of all melanins is a consequence of the superposition of a large number of nonhomogeneously broadened Gaussian transitions associated with each of the components of a melanin ensemble. These results suggest that the traditional model of eumelanin as an amorphous organic semiconductor is not required to explain its optical properties and should be thoroughly reexamined. These results have significant implications for our understanding of the physics, chemistry, and biological function of these important biological macromolecules. Indeed, one may speculate that the robust functionality of melanins in vitro is a direct consequence of its heterogeneity, i.e., chemical disorder is a "low cost" natural resource in these systems

Determination of chromophore charge states in the low pH color transition of the fluorescent protein Rtms5(H146S) via time-dependent DFT

The red fluorescent protein Rtms5H146S displays a transition from blue (absorbance λmax 590 nm) to yellow (absorbance λmax not, vert, similar453 nm) upon titration to low pH. The pKa of the reaction depends on the concentration of halide, offering promise for new expressible halide sensors. The protonation state involved in the low pH form of the chromophore remains, however, ambiguous. We report calculated excitation energies of different protonation states of an RFP chromophore model. These suggest that the relevant titration site is the phenoxy moiety of the chromophore, and the relevant base and conjugate acid are anionic and neutral chromophore species, respectively.

Trans-cis Isomerism and acylimine formation in DsRed chromophore models: Intrinsic rotation barriers

The chromophore of the red fluorescent protein DsRed contains an acylimine substituent to a GFP-like chromophore structure. The acylimine is formed from the trans peptide linkage between residues F65 and Q66 in immature DsRed, but has a cis configuration in the mature protein. The relationship between acylimine formation and trans–cis isomerization is unresolved. We have calculated bond rotation profiles for models of mature and immature DsRed chromophores using B3LYP DFT. The isomerization barrier is substantially reduced in acylimine-substituted models, providing prima facie evidence that acylimine formation precedes trans–cis isomerization in DsRed chromophores.

A novel and rational approach towards designing aluminium chelators

Aluminium (Al) is known to be neurotoxic and has been associated with the aetiology of Alzheimer's Disease. To date, only desferrioxamine (DFO), a trihydroxamic acid siderophore has been used in the clinical environment for the removal of Al from the body. However, this drug is expensive, orally inactive and is associated with many side effects. These studies employed a theoretical approach, with the use of quantum mechanics (QM) via semi-empirical molecular orbital (MO) calculations, and a practical approach using U87-MG glioblastoma cells as a model for evaluating the influence of potential chelators on the passage of aluminium into cells. Preliminary studies involving the Cambridge Structural Database (CSD) identified that Al prefers binding to bidentate ligands in a 3:1 manner, whereby oxygen was the exclusive donating atom. Statistically significant differences in M-O bond lengths when compared to other trivalent metal ions such as Fe3+ were established and used as an acceptance criterion for subsequent MO calculations. Of the semi-empirical methods parameterised for Al, the PM3 Hamiltonian was found to give the most reliable final optimised geometries of simple 3:1 Al complexes. Consequently the PM3 Hamiltonian was used for evaluating the Hf of 3:1 complexes with more complicated ligands. No correlation exists between published stability constants and individual parameters calculated via PM3 optimisations, although investigation of the dicarboxylates reveals a correlation of 0.961 showing promise for affinity prediction of closely related ligands. A simple and inexpensive morin spectrofluorescence assay has been developed and optimised producing results comparable to atomic absorption spectroscopy methods for the quantitative analysis of Al. This assay was used in subsequent in vitro models, initially on E. coli, which indicated that Al inhibits the antimicrobial action of ciprofloxacin, a potent quinolone antibiotic. Ensuing studies using the second model, U87-MG cells, investigated the influence of chelators on the transmembrane transport of Al, identifying 1,2-diethylhydroxypyridin-4-one as a ligand showing greatest potential for chelating Al in the clinical situation. In conclusion, these studies have explored semi-empirical MO Hamiltonians and an in-vitro U87-MG cell line, both as possible methods for predicting effective chelators of Al.

Non-classical inhibitors of dihydrofolate reductase

This thesis comprises two main objectives. The first objective involved the stereochemical studies of chiral 4,6-diamino-1-aryl-1,2-dihydro-s-triazines and an investigation on how the different conformations of these stereoisomers may affect their binding affinity to the enzyme dihydrofolate reductase (DHFR). The ortho-substituted 1-aryl-1,2-dihydro-s-triazines were synthesised by the three component method. An ortho-substitution at the C6' position was observed when meta-azidocycloguanil was decomposed in acid. The ortho-substituent restricts free rotation and this gives rise to atropisomerism. Ortho-substituted 4,6-diamino-1-aryl-2-ethyl-1,2-dihydro-2-methyl-s-triazine contains two elements of chirality and therefore exists as four stereoisomers: (S,aR), (R,aS), (R,aR) and (S,aS). The energy barriers to rotation of these compounds were calculated by a semi-empirical molecular orbital program called MOPAC and they were found to be in excess of 23 kcal/mol. The diastereoisomers were resolved and enriched by C18 reversed phase h.p.l.c. Nuclear overhauser effect experiments revealed that (S,aR) and (R,aS) were the more stable pair of stereoisomers and therefore existed as the major component. The minor diastereoisomers showed greater binding affinity for the rat liver DHFR in in vitro assay. The second objective entailed the investigation into the possibility of retaining DHFR inhibitory activity by replacing the classical diamino heterocyclic moiety with an amidinyl group. 4-Benzylamino-3-nitro-N,N-dimethyl-phenylamidine was synthesised in two steps. One of the two phenylamidines indicated weak inhibition against the rat liver DHFR. This weak activity may be due to the failure of the inhibitor molecule to form strong hydrogen bonds with residue Glu-30 at the active site of the enzyme.

Structural studies of trans-N2S2 copper macrocycles

The X-ray crystal structures of two related trans-N2S2 copper macrocycles are reported. One was isolated with the copper in the divalent form and the other with copper in its univalent form affording a valuable insight into the changes of geometry and metrical parameters that occur during redox processes in macrocyclic copper complexes. A variable temperature NMR study of the copper(I) complex is reported, indicative of a chair-boat conformational change within the alkyl chain backbone of the macrocycle. It was possible to extract the relevant kinetic and thermodynamic parameters (?G‡, 57.8 kJ mol-1; ?H‡, 52.1 kJ mol-1; ?S‡, -19.2 J K-1 mol-1) for this process at 298 K. DFT molecular orbital calculations were used to confirm these observations and to calculate the energy difference (26.2 kJmol-1) between the copper(I) macrocycle in a planar and a distorted tetrahedral disposition.

Development, validation and uncertainty analysis of quantitative structure and activity relationship models for Log P of disinfection by -products

Hydrophobicity as measured by Log P is an important molecular property related to toxicity and carcinogenicity. With increasing public health concerns for the effects of Disinfection By-Products (DBPs), there are considerable benefits in developing Quantitative Structure and Activity Relationship (QSAR) models capable of accurately predicting Log P. In this research, Log P values of 173 DBP compounds in 6 functional classes were used to develop QSAR models, by applying 3 molecular descriptors, namely, Energy of the Lowest Unoccupied Molecular Orbital (ELUMO), Number of Chlorine (NCl) and Number of Carbon (NC) by Multiple Linear Regression (MLR) analysis. The QSAR models developed were validated based on the Organization for Economic Co-operation and Development (OECD) principles. The model Applicability Domain (AD) and mechanistic interpretation were explored. Considering the very complex nature of DBPs, the established QSAR models performed very well with respect to goodness-of-fit, robustness and predictability. The predicted values of Log P of DBPs by the QSAR models were found to be significant with a correlation coefficient R2 from 81% to 98%. The Leverage Approach by Williams Plot was applied to detect and remove outliers, consequently increasing R 2 by approximately 2% to 13% for different DBP classes. The developed QSAR models were statistically validated for their predictive power by the Leave-One-Out (LOO) and Leave-Many-Out (LMO) cross validation methods. Finally, Monte Carlo simulation was used to assess the variations and inherent uncertainties in the QSAR models of Log P and determine the most influential parameters in connection with Log P prediction. The developed QSAR models in this dissertation will have a broad applicability domain because the research data set covered six out of eight common DBP classes, including halogenated alkane, halogenated alkene, halogenated aromatic, halogenated aldehyde, halogenated ketone, and halogenated carboxylic acid, which have been brought to the attention of regulatory agencies in recent years. Furthermore, the QSAR models are suitable to be used for prediction of similar DBP compounds within the same applicability domain. The selection and integration of various methodologies developed in this research may also benefit future research in similar fields.

Chlorine contribution to quantitative structure and activity relationship models of disinfection by -products' quantum chemical descriptors and toxicities

Quantitative Structure-Activity Relationship (QSAR) has been applied extensively in predicting toxicity of Disinfection By-Products (DBPs) in drinking water. Among many toxicological properties, acute and chronic toxicities of DBPs have been widely used in health risk assessment of DBPs. These toxicities are correlated with molecular properties, which are usually correlated with molecular descriptors. The primary goals of this thesis are: (1) to investigate the effects of molecular descriptors (e.g., chlorine number) on molecular properties such as energy of the lowest unoccupied molecular orbital (E LUMO) via QSAR modelling and analysis; (2) to validate the models by using internal and external cross-validation techniques; (3) to quantify the model uncertainties through Taylor and Monte Carlo Simulation. One of the very important ways to predict molecular properties such as ELUMO is using QSAR analysis. In this study, number of chlorine (NCl ) and number of carbon (NC) as well as energy of the highest occupied molecular orbital (EHOMO) are used as molecular descriptors. There are typically three approaches used in QSAR model development: (1) Linear or Multi-linear Regression (MLR); (2) Partial Least Squares (PLS); and (3) Principle Component Regression (PCR). In QSAR analysis, a very critical step is model validation after QSAR models are established and before applying them to toxicity prediction. The DBPs to be studied include five chemical classes: chlorinated alkanes, alkenes, and aromatics. In addition, validated QSARs are developed to describe the toxicity of selected groups (i.e., chloro-alkane and aromatic compounds with a nitro- or cyano group) of DBP chemicals to three types of organisms (e.g., Fish, T. pyriformis, and P.pyosphoreum) based on experimental toxicity data from the literature. The results show that: (1) QSAR models to predict molecular property built by MLR, PLS or PCR can be used either to select valid data points or to eliminate outliers; (2) The Leave-One-Out Cross-Validation procedure by itself is not enough to give a reliable representation of the predictive ability of the QSAR models, however, Leave-Many-Out/K-fold cross-validation and external validation can be applied together to achieve more reliable results; (3) E LUMO are shown to correlate highly with the NCl for several classes of DBPs; and (4) According to uncertainty analysis using Taylor method, the uncertainty of QSAR models is contributed mostly from NCl for all DBP classes.

Ground and Electronic Excited States from Pairing Matrix Fluctuation and Particle-Particle Random Phase Approximation

The accurate description of ground and electronic excited states is an important and challenging topic in quantum chemistry. The pairing matrix fluctuation, as a counterpart of the density fluctuation, is applied to this topic. From the pairing matrix fluctuation, the exact electron correlation energy as well as two electron addition/removal energies can be extracted. Therefore, both ground state and excited states energies can be obtained and they are in principle exact with a complete knowledge of the pairing matrix fluctuation. In practice, considering the exact pairing matrix fluctuation is unknown, we adopt its simple approximation --- the particle-particle random phase approximation (pp-RPA) --- for ground and excited states calculations. The algorithms for accelerating the pp-RPA calculation, including spin separation, spin adaptation, as well as an iterative Davidson method, are developed. For ground states correlation descriptions, the results obtained from pp-RPA are usually comparable to and can be more accurate than those from traditional particle-hole random phase approximation (ph-RPA). For excited states, the pp-RPA is able to describe double, Rydberg, and charge transfer excitations, which are challenging for conventional time-dependent density functional theory (TDDFT). Although the pp-RPA intrinsically cannot describe those excitations excited from the orbitals below the highest occupied molecular orbital (HOMO), its performances on those single excitations that can be captured are comparable to TDDFT. The pp-RPA for excitation calculation is further applied to challenging diradical problems and is used to unveil the nature of the ground and electronic excited states of higher acenes. The pp-RPA and the corresponding Tamm-Dancoff approximation (pp-TDA) are also applied to conical intersections, an important concept in nonadiabatic dynamics. Their good description of the double-cone feature of conical intersections is in sharp contrast to the failure of TDDFT. All in all, the pairing matrix fluctuation opens up new channel of thinking for quantum chemistry, and the pp-RPA is a promising method in describing ground and electronic excited states.

Chlorine Contribution to Quantitative Structure and Activity Relationship Models of Disinfection By-Products' Quantum Chemical Descriptors and Toxicities

Quantitative Structure-Activity Relationship (QSAR) has been applied extensively in predicting toxicity of Disinfection By-Products (DBPs) in drinking water. Among many toxicological properties, acute and chronic toxicities of DBPs have been widely used in health risk assessment of DBPs. These toxicities are correlated with molecular properties, which are usually correlated with molecular descriptors. The primary goals of this thesis are: 1) to investigate the effects of molecular descriptors (e.g., chlorine number) on molecular properties such as energy of the lowest unoccupied molecular orbital (ELUMO) via QSAR modelling and analysis; 2) to validate the models by using internal and external cross-validation techniques; 3) to quantify the model uncertainties through Taylor and Monte Carlo Simulation. One of the very important ways to predict molecular properties such as ELUMO is using QSAR analysis. In this study, number of chlorine (NCl) and number of carbon (NC) as well as energy of the highest occupied molecular orbital (EHOMO) are used as molecular descriptors. There are typically three approaches used in QSAR model development: 1) Linear or Multi-linear Regression (MLR); 2) Partial Least Squares (PLS); and 3) Principle Component Regression (PCR). In QSAR analysis, a very critical step is model validation after QSAR models are established and before applying them to toxicity prediction. The DBPs to be studied include five chemical classes: chlorinated alkanes, alkenes, and aromatics. In addition, validated QSARs are developed to describe the toxicity of selected groups (i.e., chloro-alkane and aromatic compounds with a nitro- or cyano group) of DBP chemicals to three types of organisms (e.g., Fish, T. pyriformis, and P.pyosphoreum) based on experimental toxicity data from the literature. The results show that: 1) QSAR models to predict molecular property built by MLR, PLS or PCR can be used either to select valid data points or to eliminate outliers; 2) The Leave-One-Out Cross-Validation procedure by itself is not enough to give a reliable representation of the predictive ability of the QSAR models, however, Leave-Many-Out/K-fold cross-validation and external validation can be applied together to achieve more reliable results; 3) ELUMO are shown to correlate highly with the NCl for several classes of DBPs; and 4) According to uncertainty analysis using Taylor method, the uncertainty of QSAR models is contributed mostly from NCl for all DBP classes.

Characterization of High-Performance Organic Dyes for Dye-Sensitized Solar Cell: A DFT/TDDFT Study

Dye-sensitized solar cell (DSSC) is currently a promising technology that makes solar energy efficient and cost-effective to harness. In DSSC, metal free dyes, such indoline-containing D149 and D205, are proved to be potential alternatives for traditional metal organic dyes. In this work, a DFT/TDDFT characterization for D149 and D205 were carried out using different functionals, including B3LYP, MPW1K, CAM-B3LYP and PBE0. Three different conformers for D149 and four different conformers for D205 were identified and calculated in vacuum. The performance of different functionals on calculating the maximum absorbance of the dyes in vacuum and five common solvents (acetonitrile, chloroform, ethanol, methanol, and THF) were examined and compared to determine the suitable computational setting for predicting properties of these two dyes. Furthermore, deprotonated D149 and D205 in solvents were also considered, and the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) were calculated, which elucidates the substitution effect on the rhodanine ring of D149 and D205 dyes on their efficiency. Finally, D149 and D205 molecules were confirmed to be firmly anchored on ZnO surface by periodic DFT calculations. These results would shed light on the design of new highly efficiency metal-free dyes.

Calculated ionization potentials of the linear alkanes

The equivalent orbital (EO) method is investigated and used for predicting outer and inner ionization potentials of the linear alkanes. The calculated ionization potentials are in good agreement with those observed in photoelectron spectra provided that a set of 12 parameters is used in the theory. An optimization technique is used to find the best values for thle parameters and a single transferable parameter set can be found which is applicable to all the n-alkanes. A good fit to the experimental results can only be obtained if the uppermost molecular orbital of the n-alkanes is an antisymmetrical orbital built up from CH equivalent orbitals.

The OPIT system part I.

The OPIT program is briefly described. OPIT is a basis-set-optimising, self-consistent field, molecular orbital program for calculating properties of closed-shell ground states of atoms and molecules. A file handling technique is then put forward which enables core storage to be used efficiently in large FORTRAN scientific applications programs. Hashing and list processing techniques, of the type frequently used in writing system software and computer operating systems, are here applied to the creation of data files (integral label and value lists etc.). Files consist of a chained series of blocks which may exist in core or on backing store or both. Efficient use of core store is achieved and the processes of file deletion, file re-writing and garbage collection of unused blocks can be easily arranged. The scheme is exemplified with reference to the OPIT program. A subsequent paper will describe a job scheduling scheme for large programs of this sort.