A first-principles density-functional calculation of the electronic and vibrational structure of the key melanin monomers

We report first-principles density-functional calculations for hydroquinone (HQ), indolequinone (IQ), and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of biomacromolecules with important biological functions (including photoprotection) and with the potential for certain bioengineering applications. We have used the difference of self-consistent fields method to study the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, HL. We show that HL is similar in IQ and SQ, but approximately twice as large in HQ. This may have important implications for our understanding of the observed broadband optical absorption of the eumelanins. The possibility of using this difference in HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to nondestructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behavior of the eumelanins.

Capillary coexistence and criticality in mesopores: Modification of the Kelvin theory

The classical model of capillary equilibrium in cylindrical pores is modified here by the introduction of molecular concepts and the solid fluid interaction potential. The new approach accurately predicts capillary coexistence and criticality, with results quantitatively matching those from density functional theory for nitrogen adsorption, while also predicting condensation pressures in agreement with reported experimental findings for MCM-41. The larger critical pore size for nitrogen adsorption in these materials, however, suggests a modification of the potential function parameters, evaluated here from data for hydroxylated silica.

Antimalarial Treatment May Have a Time-Dependent Effect on Lupus Survival

Objective. To evaluate the beneficial effect of antimalarial treatment on lupus survival in a large, multiethnic, international longitudinal inception cohort. Methods. Socioeconomic and demographic characteristics, clinical manifestations, classification criteria, laboratory findings, and treatment variables were examined in patients with systemic lupus erythematosus (SLE) from the Grupo Latino Americano de Estudio del Lupus Eritematoso (GLADEL) cohort. The diagnosis of SLE, according to the American College of Rheumatology criteria, was assessed within 2 years of cohort entry. Cause of death was classified as active disease, infection, cardiovascular complications, thrombosis, malignancy, or other cause. Patients were subdivided by antimalarial use, grouped according to those who had received antimalarial drugs for at least 6 consecutive months (user) and those who had received antimalarial drugs for <6 consecutive months or who had never received antimalarial drugs (nonuser). Results. Of the 1,480 patients included in the GLADEL cohort, 1,141 (77%) were considered antimalarial users, with a mean duration of drug exposure of 48.5 months (range 6-98 months). Death occurred in 89 patients (6.0%). A lower mortality rate was observed in antimalarial users compared with nonusers (4.4% versus 11.5%; P < 0.001). Seventy patients (6.1%) had received antimalarial drugs for 6-11 months, 146 (12.8%) for 1-2 years, and 925 (81.1%) for >2 years. Mortality rates among users by duration of antimalarial treatment (per 1,000 person-months of followup) were 3.85 (95% confidence interval [95% CI] 1.41-8.37), 2.7 (95% CI 1.41-4.76), and 0.54 (95% CI 0.37-0.77), respectively, while for nonusers, the mortality rate was 3.07 (95% CI 2.18-4.20) (P for trend < 0.001). After adjustment for potential confounders in a Cox regression model, antimalarial use was associated with a 38% reduction in the mortality rate (hazard ratio 0.62, 95% CI 0.39-0.99). Conclusion. Antimalarial drugs were shown to have a protective effect, possibly in a time-dependent manner, on SLE survival. These results suggest that the use of antimalarial treatment should be recommended for patients with lupus.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Analysis of adsorption data of graphitized thermal carbon black with a DFT-lattice gas theory

In this paper we analyzed the adsorption of gases and vapors on graphitised thermal carbon black by using a modified DFT-lattice theory, in which we assume that the behavior of the first layer in the adsorption film is different from those of second and higher layers. The effects of various parameters on the topology of the adsorption isotherm were first investigated, and the model was then applied in the analysis of adsorption data of numerous substances on carbon black. We have found that the first layer in the adsorption film behaves differently from the second and higher layers in such a way that the adsorbate-adsorbate interaction energy in the first layer is less than that of second and higher layers, and the same is observed for the partition function. Furthermore, the adsorbate-adsorbate and adsorbate-adsorbent interaction energies obtained from the fitting are consistently lower than the corresponding values obtained from the viscosity data and calculated from the Lorentz-Berthelot rule, respectively.

A new method for analysing the equilibrium and time-dependent behaviour of Markovian models

Many large-scale stochastic systems, such as telecommunications networks, can be modelled using a continuous-time Markov chain. However, it is frequently the case that a satisfactory analysis of their time-dependent, or even equilibrium, behaviour is impossible. In this paper, we propose a new method of analyzing Markovian models, whereby the existing transition structure is replaced by a more amenable one. Using rates of transition given by the equilibrium expected rates of the corresponding transitions of the original chain, we are able to approximate its behaviour. We present two formulations of the idea of expected rates. The first provides a method for analysing time-dependent behaviour, while the second provides a highly accurate means of analysing equilibrium behaviour. We shall illustrate our approach with reference to a variety of models, giving particular attention to queueing and loss networks. (C) 2003 Elsevier Ltd. All rights reserved.

Synthesis, characterization, electrochemical behavior and in vitro protein tyrosine kinase inhibitory activity of the cymene-halogenobenzohydroxamato [Ru(eta(6)-cymene)(bha)Cl] complexes

The ruthenium(II)-cymene complexes [Ru(eta(6)-cymene)(bha)Cl] with substituted halogenobenzohydroxamato (bha) ligands (substituents = 4-F, 4-Cl, 4-Br, 2,4-F-2, 3,4-F-2, 2,5-F-2, 2,6-F-2) have been synthesized and characterized by elemental analysis, IR, H-1 NMR, C-13 NMR, cyclic voltammetry and controlled-potential electrolysis, and density functional theory (DFT) studies. The compositions of their frontier molecular orbitals (MOs) were established by DFT calculations, and the oxidation and reduction potentials are shown to follow the orders of the estimated vertical ionization potential and electron affinity, respectively. The electrochemical E-L Lever parameter is estimated for the first time for the various bha ligands, which can thus be ordered according to their electron-donor character. All complexes exhibit very strong protein tyrosine kinase (PTK) inhibitory activity, even much higher than that of genistein, the clinically used PTK inhibitory drug. The complex containing the 2,4-difluorobenzohydroxamato ligand is the most active one, and the dependences of the PTK activity of the complexes and of their redox potentials on the ring substituents are discussed. (C) 2012 Elsevier B.V. All rights reserved.

Microscopic Theory of Anchoring Transitions at the Surfaces of Pure Liquid-Crystals and their Mixtures .1. The Fowler Approximation

We have generalized earlier work on anchoring of nematic liquid crystals by Sullivan, and Sluckin and Poniewierski, in order to study transitions which may occur in binary mixtures of nematic liquid crystals as a function of composition. Microscopic expressions have been obtained for the anchoring energy of (i) a liquid crystal in contact with a solid aligning surface; (ii) a liquid crystal in contact with an immiscible isotropic medium; (iii) a liquid crystal mixture in contact with a solid aligning surface. For (iii), possible phase diagrams of anchoring angle versus dopant concentration have been calculated using a simple liquid crystal model. These exhibit some interesting features including re-entrant conical anchoring, for what are believed to be realistic values of the molecular parameters. A way of relaxing the most drastic approximation implicit in the above approach is also briefly discussed.

Simulation and theory of hybrid aligned liquid crystal films

We present a study of the effects of nanoconfinement on a system of hard Gaussian overlap particles interacting with planar substrates through the hard-needle-wall potential, extending earlier work by two of us [D. J. Cleaver and P. I. C. Teixeira, Chem. Phys. Lett. 338, 1 (2001)]. Here, we consider the case of hybrid films, where one of the substrates induces strongly homeotropic anchoring, while the other favors either weakly homeotropic or planar anchoring. These systems are investigated using both Monte Carlo simulation and density-functional theory, the latter implemented at the level of Onsager's second-virial approximation with Parsons-Lee rescaling. The orientational structure is found to change either continuously or discontinuously depending on substrate separation, in agreement with earlier predictions by others. The theory is seen to perform well in spite of its simplicity, predicting the positional and orientational structure seen in simulations even for small particle elongations.

Automatic learning for the classification of chemical reactions and in statistical thermodynamics

This Thesis describes the application of automatic learning methods for a) the classification of organic and metabolic reactions, and b) the mapping of Potential Energy Surfaces(PES). The classification of reactions was approached with two distinct methodologies: a representation of chemical reactions based on NMR data, and a representation of chemical reactions from the reaction equation based on the physico-chemical and topological features of chemical bonds. NMR-based classification of photochemical and enzymatic reactions. Photochemical and metabolic reactions were classified by Kohonen Self-Organizing Maps (Kohonen SOMs) and Random Forests (RFs) taking as input the difference between the 1H NMR spectra of the products and the reactants. The development of such a representation can be applied in automatic analysis of changes in the 1H NMR spectrum of a mixture and their interpretation in terms of the chemical reactions taking place. Examples of possible applications are the monitoring of reaction processes, evaluation of the stability of chemicals, or even the interpretation of metabonomic data. A Kohonen SOM trained with a data set of metabolic reactions catalysed by transferases was able to correctly classify 75% of an independent test set in terms of the EC number subclass. Random Forests improved the correct predictions to 79%. With photochemical reactions classified into 7 groups, an independent test set was classified with 86-93% accuracy. The data set of photochemical reactions was also used to simulate mixtures with two reactions occurring simultaneously. Kohonen SOMs and Feed-Forward Neural Networks (FFNNs) were trained to classify the reactions occurring in a mixture based on the 1H NMR spectra of the products and reactants. Kohonen SOMs allowed the correct assignment of 53-63% of the mixtures (in a test set). Counter-Propagation Neural Networks (CPNNs) gave origin to similar results. The use of supervised learning techniques allowed an improvement in the results. They were improved to 77% of correct assignments when an ensemble of ten FFNNs were used and to 80% when Random Forests were used. This study was performed with NMR data simulated from the molecular structure by the SPINUS program. In the design of one test set, simulated data was combined with experimental data. The results support the proposal of linking databases of chemical reactions to experimental or simulated NMR data for automatic classification of reactions and mixtures of reactions. Genome-scale classification of enzymatic reactions from their reaction equation. The MOLMAP descriptor relies on a Kohonen SOM that defines types of bonds on the basis of their physico-chemical and topological properties. The MOLMAP descriptor of a molecule represents the types of bonds available in that molecule. The MOLMAP descriptor of a reaction is defined as the difference between the MOLMAPs of the products and the reactants, and numerically encodes the pattern of bonds that are broken, changed, and made during a chemical reaction. The automatic perception of chemical similarities between metabolic reactions is required for a variety of applications ranging from the computer validation of classification systems, genome-scale reconstruction (or comparison) of metabolic pathways, to the classification of enzymatic mechanisms. Catalytic functions of proteins are generally described by the EC numbers that are simultaneously employed as identifiers of reactions, enzymes, and enzyme genes, thus linking metabolic and genomic information. Different methods should be available to automatically compare metabolic reactions and for the automatic assignment of EC numbers to reactions still not officially classified. In this study, the genome-scale data set of enzymatic reactions available in the KEGG database was encoded by the MOLMAP descriptors, and was submitted to Kohonen SOMs to compare the resulting map with the official EC number classification, to explore the possibility of predicting EC numbers from the reaction equation, and to assess the internal consistency of the EC classification at the class level. A general agreement with the EC classification was observed, i.e. a relationship between the similarity of MOLMAPs and the similarity of EC numbers. At the same time, MOLMAPs were able to discriminate between EC sub-subclasses. EC numbers could be assigned at the class, subclass, and sub-subclass levels with accuracies up to 92%, 80%, and 70% for independent test sets. The correspondence between chemical similarity of metabolic reactions and their MOLMAP descriptors was applied to the identification of a number of reactions mapped into the same neuron but belonging to different EC classes, which demonstrated the ability of the MOLMAP/SOM approach to verify the internal consistency of classifications in databases of metabolic reactions. RFs were also used to assign the four levels of the EC hierarchy from the reaction equation. EC numbers were correctly assigned in 95%, 90%, 85% and 86% of the cases (for independent test sets) at the class, subclass, sub-subclass and full EC number level,respectively. Experiments for the classification of reactions from the main reactants and products were performed with RFs - EC numbers were assigned at the class, subclass and sub-subclass level with accuracies of 78%, 74% and 63%, respectively. In the course of the experiments with metabolic reactions we suggested that the MOLMAP / SOM concept could be extended to the representation of other levels of metabolic information such as metabolic pathways. Following the MOLMAP idea, the pattern of neurons activated by the reactions of a metabolic pathway is a representation of the reactions involved in that pathway - a descriptor of the metabolic pathway. This reasoning enabled the comparison of different pathways, the automatic classification of pathways, and a classification of organisms based on their biochemical machinery. The three levels of classification (from bonds to metabolic pathways) allowed to map and perceive chemical similarities between metabolic pathways even for pathways of different types of metabolism and pathways that do not share similarities in terms of EC numbers. Mapping of PES by neural networks (NNs). In a first series of experiments, ensembles of Feed-Forward NNs (EnsFFNNs) and Associative Neural Networks (ASNNs) were trained to reproduce PES represented by the Lennard-Jones (LJ) analytical potential function. The accuracy of the method was assessed by comparing the results of molecular dynamics simulations (thermal, structural, and dynamic properties) obtained from the NNs-PES and from the LJ function. The results indicated that for LJ-type potentials, NNs can be trained to generate accurate PES to be used in molecular simulations. EnsFFNNs and ASNNs gave better results than single FFNNs. A remarkable ability of the NNs models to interpolate between distant curves and accurately reproduce potentials to be used in molecular simulations is shown. The purpose of the first study was to systematically analyse the accuracy of different NNs. Our main motivation, however, is reflected in the next study: the mapping of multidimensional PES by NNs to simulate, by Molecular Dynamics or Monte Carlo, the adsorption and self-assembly of solvated organic molecules on noble-metal electrodes. Indeed, for such complex and heterogeneous systems the development of suitable analytical functions that fit quantum mechanical interaction energies is a non-trivial or even impossible task. The data consisted of energy values, from Density Functional Theory (DFT) calculations, at different distances, for several molecular orientations and three electrode adsorption sites. The results indicate that NNs require a data set large enough to cover well the diversity of possible interaction sites, distances, and orientations. NNs trained with such data sets can perform equally well or even better than analytical functions. Therefore, they can be used in molecular simulations, particularly for the ethanol/Au (111) interface which is the case studied in the present Thesis. Once properly trained, the networks are able to produce, as output, any required number of energy points for accurate interpolations.

Redox-active cytotoxic diorganotin(IV) cycloalkylhydroxamate complexes with different ring sizes: Reduction behaviour and theoretical interpretation

Two series of new diorganotin(IV) cycloalkylhydroxamate complexes with different ring sizes (cyclopropyl, cyclobutyl, cyclopentyl and cyclohexyl), formulated as the mononuclear [R2Sn(HL)(2)] (1:2) (a, R=Bu-n and Ph) and the polymeric [R2SnL](n) (1:1) (b, R=Bu-n) compounds, were prepared and fully characterized. Single crystal X-ray diffraction for [(Bu2Sn)-Bu-n{C5H9C(O)NHO}(2)] (3a) discloses the cis geometry and strong intermolecular NH center dot center dot center dot O interactions. The in vitro cytotoxic activities of the complexes were evaluated against HL-60, Bel-7402, BGC-823 and KB human tumour cell lines, the greater activity concerning [(Bu2Sn)-Bu-n(HL)(2)] [HL=C3H5C(O)NHO (1a), C6H11C(O)NHO (4a)] towards BGC-823. The complexes undergo, by cyclic voltammetry and controlled-potential electrolysis, one irreversible overall two-electron cathodic process at a reduction potential that does not appear to correlate with the antitumour activity. The electrochemical behaviour of [R2Sn(C5H9C(O)NHO)(2)] [R=Bu-n (3a), Ph (7a)] was also investigated using density functional theory (DFT) methods, showing that the ultimate complex structure and the mechanism of its formation are R dependent: for the aromatic (R = Ph) complex, the initial reduction step is centred on the phenyl ligands and at the metal, being followed by a second reduction with Sn-O and Sn-C ruptures, whereas for the alkyl (R=Bu-n) complex the first reduction step is centred on one of the hydroxamate ligands and is followed by a second reduction with Sn-O bond cleavages and preservation of the alkyl ligands. In both cases, the final complexes are highly coordinative unsaturated Sn-II species with the cis geometry, features that can be of biological significance.