Use of density functional calculations to predict the regioselectivity of drugs and molecules metabolized by aldehyde oxidase.

Aldehyde oxidase is a molybdenum hydroxylase that catalyzes the oxidation of aldehydes and nitrogen-containing heterocycles. The enzyme plays a dual role in the metabolism of physiologically important endogenous compounds and the biotransformation of xenobiotics. Using density functional theory methods, geometry optimization of tetrahedral intermediates of drugs and druglike compounds was examined to predict the likely metabolites of aldehyde oxidase. The calculations suggest that the lowest energy tetrahedral intermediate resulting from the initial substrate corresponds to the observed metabolite >or=90% of the time. Additional calculations were performed on a series of heterocyclic compounds where the products resulting from metabolism by xanthine oxidase and aldehyde oxidase differ in many instances. Again, the lowest energy tetrahedral intermediate corresponded to the observed product of aldehyde oxidase metabolism >or=90% for the compounds examined, while the observed products of xanthine oxidase were not well predicted.

Inhibition of recombinant cytochrome P450 isoforms 2D6 and 2C9 by diverse drug-like molecules.

The affinities of a diverse set of 500 drug-like molecules to cytochrome P450 isoforms 2C9 and 2D6 were measured using recombinant expressed enzyme. The dose-response curve of each compound was fitted with a series of equations representing typical or various types of atypical kinetics. Atypical kinetics was identified where the Akaike Information Criterion, plus other criteria, suggested the kinetics was more complex than expected for a Michaelis-Menten model. Approximately 20% of the compounds were excluded due to poor solubility, and approximately 15% were excluded due to fluorescence interference. Of the remaining compounds, roughly half were observed to bind with an affinity of 200 microM or lower for each of the two isoforms. Atypical kinetics was observed in 18% of the compounds that bind to cytochrome 2C9, but less than 2% for 2D6. The resulting collection of competitive inhibitors and inactive compounds were analyzed for trends in binding affinity. For CYP2D6, a clear relationship between polar surface area and charge was observed, with the most potent inhibitors having a formal positive charge and a low percent polar surface area. For CYP2C9, no clear trend between activity and physicochemical properties could be seen for the group as a whole; however, certain classes of compounds have altered frequencies of activity and atypical kinetics.

Empirical regioselectivity models for human cytochromes P450 3A4, 2D6, and 2C9.

Cytochromes P450 3A4, 2D6, and 2C9 metabolize a large fraction of drugs. Knowing where these enzymes will preferentially oxidize a molecule, the regioselectivity, allows medicinal chemists to plan how best to block its metabolism. We present QSAR-based regioselectivity models for these enzymes calibrated against compiled literature data of drugs and drug-like compounds. These models are purely empirical and use only the structures of the substrates, in contrast to those models that simulate a specific mechanism like hydrogen radical abstraction, and/or use explicit models of active sites. Our most predictive models use three substructure descriptors and two physical property descriptors. Descriptor importances from the random forest QSAR method show that other factors than the immediate chemical environment and the accessibility of the hydrogen affect regioselectivity in all three isoforms. The cross-validated predictions of the models are compared to predictions from our earlier mechanistic model (Singh et al. J. Med. Chem. 2003, 46, 1330-1336) and predictions from MetaSite (Cruciani et al. J. Med. Chem. 2005, 48, 6970-6979).

Theoretical examination of Mg(2+)-mediated hydrolysis of a phosphodiester linkage as proposed for the hammerhead ribozyme.

The hammerhead ribozyme is an RNA molecule capable of self-cleavage at a unique site within its sequence. Hydrolysis of this phosphodiester linkage has been proposed to occur via an in-line attack geometry for nucleophilic displacement by the 2'-hydroxyl on the adjoining phosphorus to generate a 2',3'-cyclic phosphate ester with elimination of the 5'-hydroxyl group, requiring a divalent metal ion under physiological conditions. The proposed S(N)2(P) reaction mechanism was investigated using density functional theory calculations incorporating the hybrid functional B3LYP to study this metal ion-dependent reaction with a tetraaquo magnesium (II)-bound hydroxide ion. For the Mg(2+)-catalyzed reaction, the gas-phase geometry optimized calculations predict two transition states with a kinetically insignificant, yet clearly defined, pentacoordinate intermediate. The first transition state located for the reaction is characterized by internal nucleophilic attack coupled to proton transfer. The second transition state, the rate-determining step, involves breaking of the exocyclic P-O bond where a metal-ligated water molecule assists in the departure of the leaving group. These calculations demonstrate that the reaction mechanism incorporating a single metal ion, serving as a Lewis acid, functions as a general base and can afford the necessary stabilization to the leaving group by orienting a water molecule for catalysis.

The circumsphere as a tool to assess distortion in [4Fe-4S] atom clusters.

The geometry proposition that "four points not in a plane describe one and only one sphere" provides a novel tool for analyzing protein-induced distortions in [4Fe-4S] clusters. A geometrically perfect reference structure comprises interlaced, regular tetrahedra of Fe, S, and S gamma atoms having T(d) symmetry. Three circumspheres are defined by the three sets of four atoms, the circumcenters of which are unique points within the cluster. The structure is thus re-defined by the positions of the circumcenters in xyz space and the r, theta, phi of each atom on its respective sphere. Analysis of 12 high-resolution structures of protein-bound and small molecule [4Fe-4S](SR)(4) clusters revealed: (a) the circumcenters are generally non-coincident by approximately 0.01 to approximately 0.06 A; (b) the Fe radius, r(Fe), is nominally independent of core oxidation state, having values between 1.66 to 1.69 A, whereas r(S) and r(SG), which have ranges of 2.18-2.24 A and 3.87-3.94 A, respectively, both increase by as much as approximately 3% upon reduction from the 3+ to the 1+ core valence; (c) deviation of some atoms from the theta, phi of a perfect tetrahedron can be large, approximately 10 degrees, and sets of atoms can show patterns of motion on their spheres that result from changes in Fe-S bond lengths. Density functional theory calculations suggest that the [4Fe-4S] core itself requires rather little energy to distort (approximately 2 kcal/mol), whereas significantly more energy is required to distort the Sgamma shell (~4 kcal/mol) to that of cluster I in Clostridium acidurici ferredoxin.

Density functional and reduction potential calculations of Fe4S4 clusters.

Density functional theory geometry optimizations and reduction potential calculations are reported for all five known oxidation states of [Fe(4)S(4)(SCH(3))(4)](n)()(-) (n = 0, 1, 2, 3, 4) clusters that form the active sites of iron-sulfur proteins. The geometry-optimized structures tend to be slightly expanded relative to experiment, with the best comparison found in the [Fe(4)S(4)(SCH(3))(4)](2)(-) model cluster, having bond lengths 0.03 A longer on average than experimentally observed. Environmental effects are modeled with a continuum dielectric, allowing the solvent contribution to the reduction potential to be calculated. The calculated protein plus solvent effects on the reduction potentials of seven proteins (including high potential iron proteins, ferredoxins, the iron protein of nitrogenase, and the "X", "A", and "B" centers of photosystem I) are also examined. A good correlation between predicted and measured absolute reduction potentials for each oxidation state of the cluster is found, both for relative potentials within a given oxidation state and for the absolute potentials for all known couples. These calculations suggest that the number of amide dipole and hydrogen bonding interactions with the Fe(4)S(4) clusters play a key role in modulating the accessible redox couple. For the [Fe(4)S(4)](0) (all-ferrous) system, the experimentally observed S = 4 state is calculated to lie lowest in energy, and the predicted geometry and electronic properties for this state correlate well with the EXAFS and Mössbauer data. Cluster geometries are also predicted for the [Fe(4)S(4)](4+) (all-ferric) system, and the calculated reduction potential for the [Fe(4)S(4)(SCH(3))(4)](1)(-)(/0) redox couple is in good agreement with that estimated for experimental model clusters containing alkylthiolate ligands.

Metal substitution in the active site of nitrogenase MFe(7)S(9) (M = Mo(4+), V(3+), Fe(3+)).

The unifying view that molybdenum is the essential component in nitrogenase has changed over the past few years with the discovery of a vanadium-containing nitrogenase and an iron-only nitrogenase. The principal question that has arisen for the alternative nitrogenases concerns the structures of their corresponding cofactors and their metal-ion valence assignments and whether there are significant differences with that of the more widely known molybdenum-iron cofactor (FeMoco). Spin-polarized broken-symmetry (BS) density functional theory (DFT) calculations are used to assess which of the two possible metal-ion valence assignments (4Fe(2+)4Fe(3+) or 6Fe(2+)2Fe(3+)) for the iron-only cofactor (FeFeco) best represents the resting state. For the 6Fe(2+)2Fe(3+) oxidation state, the spin coupling pattern for several spin state alignments compatible with S = 0 were generated and assessed by energy criteria. The most likely BS spin state is composed of a 4Fe cluster with spin S(a) = (7)/(2) antiferromagnetically coupled to a 4Fe' cluster with spin S(b) = (7)/(2). This state has the lowest DFT energy for the isolated FeFeco cluster and displays calculated Mössbauer isomer shifts consistent with experiment. Although the S = 0 resting state of FeFeco has recently been proposed to have metal-ion valencies of 4Fe(2+)4Fe(3+) (derived from experimental Mössbauer isomer shifts), our isomer shift calculations for the 4Fe(2+)4Fe(3+) oxidation state are in poorer agreement with experiment. Using the Mo(4+)6Fe(2+)Fe(3+) oxidation level of the cofactor as a starting point, the structural consequences of replacement of molybdenum (Mo(4+)) with vanadium (V(3+)) or iron (Fe(3+)) in the cofactor have been investigated. The size of the cofactor cluster shows a dependency on the nature of the heterometal and increases in the order FeMoco < FeVco < FeFeco.

Insights into properties and energetics of iron-sulfur proteins from simple clusters to nitrogenase.

Some of the principal physical features of iron-sulfur clusters in proteins are analyzed, including metal-ligand covalency, spin polarization, spin coupling, valence delocalization, valence interchange and small reorganization energies, with emphasis on recent spectroscopic and theoretical work. The current state of structural, spectroscopic, and computational knowledge for the iron-sulfur clusters in the nitrogenase iron and iron-molybdenum proteins is examined by comparison and contrast to 'simpler' ironclusters. The differing interactions of the nitrogenase iron and iron-molybdenum clusters compared with those of other iron-sulfur clusters with the protein and solvent environment are also explored.