Computer-aided drug design and ADMET predictions for identification and evaluation of novel potential farnesyltransferase inhibitors in cancer therapy

We have used various computational methodologies including molecular dynamics, density functional theory, virtual screening, ADMET predictions and molecular interaction field studies to design and analyze four novel potential inhibitors of farnesyltransferase (FTase). Evaluation of two proposals regarding their drug potential as well as lead compounds have indicated them as novel promising FTase inhibitors, with theoretically interesting pharmacotherapeutic profiles, when Compared to the very active and most cited FTase inhibitors that have activity data reported, which are launched drugs or compounds in clinical tests. One of our two proposals appears to be a more promising drug candidate and FTase inhibitor, but both derivative molecules indicate potentially very good pharmacotherapeutic profiles in comparison with Tipifarnib and Lonafarnib, two reference pharmaceuticals. Two other proposals have been selected with virtual screening approaches and investigated by LIS, which suggest novel and alternatives scaffolds to design future potential FTase inhibitors. Such compounds can be explored as promising molecules to initiate a research protocol in order to discover novel anticancer drug candidates targeting farnesyltransferase, in the fight against cancer. (C) 2009 Elsevier Inc. All rights reserved.

Molecular dynamics, flexible docking, virtual screening, ADMET predictions, and molecular interaction field studies to design novel potential MAO-B inhibitors

Monoamine oxidase is a flavoenzyme bound to the mitochondrial outer membranes of the cells, which is responsible for the oxidative deamination of neurotransmitter and dietary amines. It has two distinct isozymic forms, designated MAO-A and MAO-B, each displaying different substrate and inhibitor specificities. They are the well-known targets for antidepressant, Parkinson`s disease, and neuroprotective drugs. Elucidation of the x-ray crystallographic structure of MAO-B has opened the way for the molecular modeling studies. In this work we have used molecular modeling, density functional theory with correlation, virtual screening, flexible docking, molecular dynamics, ADMET predictions, and molecular interaction field studies in order to design new molecules with potential higher selectivity and enzymatic inhibitory activity over MAO-B.

Adsorption in mesopores: A molecular-continuum model with application to MCM-41

The classical model of surface layering followed by capillary condensation during adsorption in mesopores, is modified here by consideration of the adsorbate solid interaction potential. The new theory accurately predicts the capillary coexistence curve as well as pore criticality, matching that predicted by density functional theory. The model also satisfactorily predicts the isotherm for nitrogen adsorption at 77.4 K on MCM-41 material of various pore sizes, synthesized and characterized in our laboratory, including the multilayer region, using only data on the variation of condensation pressures with pore diameter. The results indicate a minimum mesopore diameter for the surface layering model to hold as 14.1 Å, below which size micropore filling must occur, and a minimum pore diameter for mechanical stability of the hemispherical meniscus during desorption as 34.2 Å. For pores in-between these two sizes reversible condensation is predicted to occur, in accord with the experimental data for nitrogen adsorption on MCM-41 at 77.4 K.

Characterization of pore size distributions of mesoporous materials from adsorption isotherms

In this article, a new hybrid model for estimating the pore size distribution of micro- and mesoporous materials is developed, and tested with the adsorption data of nitrogen, oxygen, and argon on ordered mesoporous materials reported in the literature. For the micropore region, the model uses the Dubinin-Rudushkevich (DR) isotherm with the Chen-Yang modification. A recent isotherm model of the authors for nonporous materials, which uses a continuum-mechanical model for the multilayer region and the Unilan model for the submonolayer region, has been extended for adsorption in mesopores. The experimental data is inverted using regularization to obtain the pore size distribution. The present model was found to be successful in predicting the pore size distribution of pure as well as binary physical mixtures of MCM-41 synthesized with different templates, with results in agreement with those from the XRD method and nonlocal density functional theory. It was found that various other recent methods, as well as the classical Broekhoff and de Beer method, underpredict the pore diameter of MCM-41. The present model has been successfully applied to MCM-48, SBA's, CMK, KIT, HMS, FSM, MTS, mesoporous fly ash, and a large number of other regular mesoporous materials.

Characterization of pore wall heterogeneity in nanoporous carbons using adsorption: the slit pore model revisited

In this paper, we propose a new nonlocal density functional theory characterization procedure, the finite wall thickness model, for nanoporous carbons, whereby heterogeneity of pore size and pore walls in the carbon is probed simultaneously. We determine the pore size distributions and pore wall thickness distributions of several commercial activated carbons and coal chars, with good correspondence with X-ray diffraction. It is shown that the conventional infinite wall thickness approach overestimates the pore size slightly. Pore-pore correlation has been shown to have a negligible effect on prediction of pore size and pore wall thickness distributions for small molecules such as argon used in characterization. By utilizing the structural parameters (pore size and pore wall thickness distribution) in the generalized adsorption isotherm (GAI) we are able to predict adsorption uptake of supercritical gases in BPL and Norit RI Extra carbons, in excellent agreement with experimental adsorption uptake data up to 60 MPa. The method offers a useful technique for probing features of the solid skeleton, hitherto studied by crystallographic methods.

Photochemical Reactions of fac-[Mn(CO)(3)(phen)imidazole](+): Evidence for Long-Lived Radical Species Intermediates

The electronic absorption spectrum of fac[Mn(CO)(3)(phen)imH](+), fac-1 in CH(2)Cl(2) is characterized by a strong absorption band at 378 nm (epsilon(max) = 3200 mol(-1) L cm(-1)). On the basis of quantum mechanical calculations, the visible absorption band has been assigned to ligand-to-ligand charge-transfer (LLCT, im -> phen) and metal-to-ligand charge-transfer (MLCT, Mn -> phen) charge transfer transition. When fac-1 in CH(2)Cl(2) is irradiated with 350 nm continuous light, the absorption features are gradually shifted to represent those of the meridional complex mer-[Mn(CO)(3)(phen)imH](+), mer-1 (lambda(max) = 556 nm). The net photoreaction under these conditions is a photoisomerization, although, the presence of the long-lived radical species was also detected by (1)H NMR and FTIR spectroscopy. 355 nm continuous photolysis of fac-1 in CH(3)CN solution also gives the long-lived intermediate which is readly trapped by metylviologen (MV(2+)) giving rise to the formation of the one-electron reduced methyl viologen (MV(center dot+)). The UV-vis spectra monitored during the slow (45 min) thermal back reaction exhibited isosbestic conversion at 426 nm. On the basis of spectroscopic techniques and quantum mechanical calculations, the role of the radicals produced is analyzed.

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters,. but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.

Synthesis, X-ray analysis and DFT study of 2-amino-3-(N-oxipiridin-4-ilsulfanil)-propionic acid dihydrate

The title compound, C(8)H(14)N(2)O(5)S 2(H(2)O), 2-amino-3-(N-oxipiridin-4-ilsulfanil)-propionic acid dihydrate, is obtained by the reaction of cysteine and 4-nitropyridine N-oxide in dimethylformamide, removing the NO(2) group from the benzene ring and releasing nitrous acid into the solution. The molecule exists as a Zwitterion. Hydrogen bond interactions involving the title molecule and water molecules allow the formation of R(5)(5)(23) edge fused rings parallel to (010). Water molecules are connected independently, forming infinite chains (wires), in square wave form, along the b-axis. The chirality of the cysteine molecule used in the synthesis is retained in the title molecule. A density functional theory (DFT) optimized structure at the B3LYP/6-311G(3df,2p) level allows comparison of calculated and experimental IR spectra.

A modified pore-filling isotherm for liquid-phase adsorption in activated carbon

This article modifies the usual form of the Dubinin-Radushkevich pore-filling model for application to liquid-phase adsorption data, where large molecules are often involved. In such cases it is necessary to include the repulsive part of the energy in the micropores, which is accomplished here by relating the pore potential to the fluid-solid interaction potential. The model also considers the nonideality of the bulk liquid phase through the UNIFAC activity coefficient model, as well as structural heterogeneity of the carbon. For the latter the generalized adsorption integral is used while incorporating the pore-size distribution obtained by density functional theory analysis of argon adsorption data. The model is applied here to the interpretation of aqueous phase adsorption isotherms of three different esters on three commercial activated carbons. Excellent agreement between the model and experimental data is observed, and the fitted Lennard-Jones size parameter for the adsorbate-adsorbate interactions compares well with that estimated from known critical properties, supporting the modified approach. On the other hand, the model without consideration of bulk nonideality, or when using classical models of the characteristic energy, gives much poorer bts of the data and unrealistic parameter values.

The use of liquid phase adsorption isotherms for characterization of activated carbons

The characterization of three commercial activated carbons was carried out using the adsorption of various compounds in the aqueous phase. For this purpose the generalized adsorption isotherm was employed, and a modification of the Dubinin-Radushkevich pore filling model, incorporating repulsive contributions to the pore potential as well as bulk liquid phase nonideality, was used as the local isotherm. Eight different flavor compounds were used as adsorbates, and the isotherms were jointly fitted to yield a common pore size distribution for each carbon. The bulk liquid phase nonideality was incorporated through the UNIFAC activity coefficient model, and the repulsive contribution to the pore potential was incorporated through the Steele 10-4-3 potential model. The mean micropore network coordination number for each carbon was also determined from the fitted saturation capacity based on percolation theory. Good agreement between the model and the experimental data was observed. In addition, excellent agreement between the bimodal gamma pore size distribution and density functional theory-cum-regularization-based pore size distribution obtained by argon adsorption was also observed, supporting the validity of the model. The results show that liquid phase adsorption, using adsorptive molecules of different sizes, can be an effective means of characterizing the pore size distribution as well as connectivity. Alternately, if the carbon pore size distribution is independently known, the method can be used to measure critical molecular sizes. (C) 2001 Elsevier Science.

Carbene and nitrene rearrangements: A theoretical study of cyclic allenes and carbenes, carbodiimides, and azirines

B3LYP/6-31G(d) calculations of structures, energies, and infrared spectra of several rearrangement products of (hetero)aromatic nitrenes and carbenes are reported. 3-Isoquinolylnitrene 36 ring closes to the azirine 37 prior to ring expansion to the potentially stable but unobserved seven-membered-ring carbodiimide 38 and diazacycloheptatrienylidene C-s-39S. A new, stable cycloheptatrienylidene, C-s-19S, is located on the naphthylcarbene energy surface. 4-Quinolylnitrene undergoes reaction via the azirine 50 in solution, but ring expansion to the stable seven-membered-ring ketenimine 47 under Ar matrix photolysis conditions. There is excellent agreement between calculated infrared spectra of 1,5-diazacyclohepta-1,2,4,6-tetraene 54 (obtained by photolysis of 4-pyridyl azide), 1-azacyclohepta-1,2,4,6-tetraene 5, 1-azacyclohepta-1,3,5,6-tetraene 55, and 1-azacyclohepta-1,3,4,6-tetraene 56 and the available experimental data.

2-pyrazinylnitrene and 4-pyrimidylnitrene. Ring expansion to 1,3,5-triazacyclohepta-1,2,4,6-tetraene and ring opening to (2-isocyanovinyl)carbodiimide

Tetrazolo[1,5-a]pyrazine/2-azidopyrazine 9T/9A undergo photolysis in Ar matrix at cryogenic temperatures to yield 1,3,5-triazacyclohepta-1,2,4,6-tetraene 21 as the first observable intermediate, and 1-cyanoimidazole 11 and (2-isocyanovinyl)carbodiimide 22 as the final products. The latter tautomerizes to 2-(isocyanovinyl)cyanamide 23 on warming to 40 K. The same intermediate 21 and the same final products are obtained on matrix photolysis of the isomeric tetrazolo[1,5-c]pyrimidine/4-azidopyrimidine 24T/24A. These photolysis results as well as those of the previously reported thermal ring contraction of N-15-labeled 2-pyrazinyl- and 4-pyrimidylnitrenes to 1-cyanoimidazoles can all be rationalized in terms of selective ring opening of 21 or nitrine 10 to a nitrile ylide zwitterion 28 prior to formation of the final products, 11 and 22. The results are supported by high-level ab initio and DFT calculations (CASPT2-CASSCF(6,6), G3(MP2), and B3LYP/6-31+G*) of the energies and IR spectra of the intermediates and products.

From molecular complexes to zwitterions and final products. Reactions between C3O2 and amines

The formation of molecular complexes (prereactive intermediates) between C3O2 and amines (ammonia, dimethylamine, trimethylamine, and 4-(dimethylamino)pyridine) as well as the subsequent transformation of the complexes into C3O2-amine zwitterions in cryogenic matrixes (ca. 40 K) has been observed. In the case of dimethylamine, the formation of tetramethylmalonamide has also been documented. Calculations using density functional theory (B3LYP/6-31G(2d, p)) are used to assign all above species and are in excellent agreement with the IR spectra.

Ab initio modelling of basal plane oxidation of graphenes and implications for modelling char combustion

Ab initio calculations have been performed to determine the energetics of oxygen atoms adsorbed onto graphene planes and the possible reaction path extracting carbon atorns in the form of carbon monoxide. Front the energetics it is confirmed that this reaction path will not significantly contribute to the gasification of well ordered carbonaceous chars. Modelling results which explore this limit Lire presented. (C) 2002 Elsevier Science Ltd, All rights reserved.

Use of liquid phase adsorption for characterizing pore network connectivity in activated carbon

A simple percolation theory-based method for determination of the pore network connectivity using liquid phase adsorption isotherm data combined with a density functional theory (DFT)-based pore size distribution is presented in this article. The liquid phase adsorption experiments have been performed using eight different esters as adsorbates and microporous-mesoporous activated carbons Filtrasorb-400, Norit ROW 0.8 and Norit ROX 0.8 as adsorbents. The density functional theory (DFT)-based pore size distributions of the carbons were obtained using DFT analysis of argon adsorption data. The mean micropore network coordination numbers, Z, of the carbons were determined based on DR characteristic plots and fitted saturation capacities using percolation theory. Based on this method, the critical molecular sizes of the model compounds used in this study were also obtained. The incorporation of percolation concepts in the prediction of multicomponent adsorption equilibria is also investigated, and found to improve the performance of the ideal adsorbed solution theory (IAST) model for the large molecules utilized in this study. (C) 2002 Elsevier Science B.V. All rights reserved.