Methemoglobin Formation and Characterization of Hemoglobin Adducts of Carcinogenic Aromatic Amines and Heterocyclic Aromatic Amines.

Arylamines (AAs) and heterocyclic aromatic amines (HAAs) are structurally related carcinogens formed during the combustion of tobacco or cooking of meat. They undergo cytochrome P450 mediated N-hydroxylation to form metabolites which bind to DNA and lead to mutations. The N-hydroxylated metabolites of many AAs also can undergo a co-oxidation reaction with oxy-hemolgobin (HbO2) to form methemoglobin (met-Hb) and the arylnitroso intermediates, which react with the β-Cys(93) chain of Hb to form Hb-arylsulfinamide adducts.

Highly predictive support vector machine (SVM) models for anthrax toxin lethal factor (LF) inhibitors.

Anthrax is a highly lethal, acute infectious disease caused by the rod-shaped, Gram-positive bacterium Bacillus anthracis. The anthrax toxin lethal factor (LF), a zinc metalloprotease secreted by the bacilli, plays a key role in anthrax pathogenesis and is chiefly responsible for anthrax-related toxemia and host death, partly via inactivation of mitogen-activated protein kinase kinase (MAPKK) enzymes and consequent disruption of key cellular signaling pathways. Antibiotics such as fluoroquinolones are capable of clearing the bacilli but have no effect on LF-mediated toxemia; LF itself therefore remains the preferred target for toxin inactivation.

Ligand-induced expansion of the S1' site in the anthrax toxin lethal factor.

The Bacillus anthracis lethal factor (LF) is one component of a tripartite exotoxin partly responsible for persistent anthrax cytotoxicity after initial bacterial infection. Inhibitors of the zinc metalloproteinase have been investigated as potential therapeutic agents, but LF is a challenging target because inhibitors lack sufficient selectivity or possess poor pharmaceutical properties. These structural studies reveal an alternate conformation of the enzyme, induced upon binding of specific inhibitors, that opens a previously unobserved deep pocket termed S1'(∗) which might afford new opportunities to design selective inhibitors that target this subsite.

Probing the S2' Subsite of the Anthrax Toxin Lethal Factor Using Novel N-Alkylated Hydroxamates.

The lethal factor (LF) enzyme secreted by Bacillus anthracis is a zinc hydrolase that is chiefly responsible for anthrax-related cell death. Although many studies of the design of small molecule LF inhibitors have been conducted, no LF inhibitor is yet available as a therapeutic agent. Inhibitors with considerable chemical diversity have been developed and investigated; however, the LF S2' subsite has not yet been systematically explored as a potential target for lead optimization.

Tactical Approaches to Interconverting GPCR Agonists and Antagonists.

There are many reported examples of small structural modifications to GPCR-targeted ligands leading to major changes in their functional activity, converting agonists into antagonists or vice versa. These shifts in functional activity are often accompanied by negligible changes in binding affinity. The current perspective focuses on outlining and analyzing various approaches that have been used to interconvert GPCR agonists, partial agonists, and antagonists in order to achieve the intended functional activity at a GPCR of therapeutic interest.

Statistical analysis, optimization, and prioritization of virtual screening parameters for zinc enzymes including the anthrax toxin lethal factor.

The anthrax toxin lethal factor (LF) and matrix metalloproteinase-3 (MMP-3, stromelysin-1) are popular zinc metalloenzyme drug targets, with LF primarily responsible for anthrax-related toxicity and host death, while MMP-3 is involved in cancer- and rheumatic disease-related tissue remodeling. A number of in silico screening techniques, most notably docking and scoring, have proven useful for identifying new potential drug scaffolds targeting LF and MMP-3, as well as for optimizing lead compounds and investigating mechanisms of action. However, virtual screening outcomes can vary significantly depending on the specific docking parameters chosen, and systematic statistical significance analyses are needed to prioritize key parameters for screening small molecules against these zinc systems.

Anthrax toxin lethal factor domain 3 is highly mobile and responsive to ligand binding.

The secreted anthrax toxin consists of three components: the protective antigen (PA), edema factor (EF) and lethal factor (LF). LF, a zinc metalloproteinase, compromises the host immune system primarily by targeting mitogen-activated protein kinase kinases in macrophages. Peptide substrates and small-molecule inhibitors bind LF in the space between domains 3 and 4 of the hydrolase.

Analysis of the Errors in the Electrostatically Embedded Many-Body Expansion of the Energy and the Correlation Energy for Zn and Cd Coordination Complexes with Five and Six Ligands and Use of the Analysis to Develop a Generally Successful Fragmentation Strategy.

In the present paper, we apply the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE) to the calculation of zinc-ligand and cadmium-ligand bond dissociation energies, and we analyze the errors due to various fragmentation schemes in a variety of neutral, positively charged, and negatively charged Zn and Cd coordination complexes. As a result of the analysis, we are able to present a new, simple, and unambiguous fragmentation strategy. Following this strategy, we show that both methods perform well for zinc-ligand and cadmium-ligand bond dissociation energies for all systems studied in the paper, including a model of the catalytic site of the zinc-bearing anthrax toxin lethal factor (LF), which has garnered substantial attention as a target for drug development.

Identification of novel anthrax toxin countermeasures using in silico methods.

Anthrax is an acute infectious disease caused by the spore-forming, gram-positive, rod-shaped bacterium Bacillus anthracis. The anthrax toxin lethal factor (LF) is the primary anthrax toxin component responsible for cytotoxicity and host death and has been a heavily researched target for design of postexposure therapeutics in the event of a bioterror attack. Various computer-aided drug design methodologies have proven useful for pinpointing new antianthrax drug scaffolds, optimizing existing leads and probes, and elucidating key mechanisms of action.

Development of a comprehensive, validated pharmacophore hypothesis for anthrax toxin lethal factor (LF) inhibitors using genetic algorithms, Pareto scoring, and structural biology.

Anthrax is an acute infectious disease caused by the spore-forming bacterium Bacillus anthracis. The anthrax toxin lethal factor (LF), an 89-kDa zinc hydrolase secreted by the bacilli, is the toxin component chiefly responsible for pathogenesis and has been a popular target for rational and structure-based drug design. Although hundreds of small-molecule compounds have been designed to target the LF active site, relatively few reported inhibitors have exhibited activity in cell-based assays, and no LF inhibitor is currently available to treat or prevent anthrax.

Electrostatically Embedded Many-Body Expansion for Neutral and Charged Metalloenzyme Model Systems.

The electrostatically embedded many-body (EE-MB) method has proven accurate for calculating cohesive and conformational energies in clusters, and it has recently been extended to obtain bond dissociation energies for metal-ligand bonds in positively charged inorganic coordination complexes. In the present paper, we present four key guidelines that maximize the accuracy and efficiency of EE-MB calculations for metal centers. Then, following these guidelines, we show that the EE-MB method can also perform well for bond dissociation energies in a variety of neutral and negatively charged inorganic coordination systems representing metalloenzyme active sites, including a model of the catalytic site of the zinc-bearing anthrax toxin lethal factor, a popular target for drug development.

Assessment and Validation of the Electrostatically Embedded Many-Body Expansion for Metal-Ligand Bonding.

The electrostatically embedded many-body method has been very successful for calculating cohesive energies and relative conformational energies of clusters, and here we extend it to calculate bond breaking energies for metal-ligand bonds in inorganic coordination chemistry. We find that, on average, the electrostatically embedded pairwise additive method is able to predict bond energies yielded by conventional full-system calculations done at the same level of theory to within 2.5 kcal/mol and that the electrostatically embedded three-body method consistently yields energies within 1.

Evaluation of Density Functionals, SCC-DFTB, Neglect of Diatomic Differential Overlap (NDDO) Models and Molecular Mechanics Methods for Prolyl-Leucyl-Glycinamide (PLG) and Structural Derivatives.

Prolyl-leucyl-glycinamide (PLG) is a unique endogenous peptide that modulates dopamine receptor subtypes of the D(2) receptor family within the CNS. We seek to elucidate the structural basis and molecular mechanism by which PLG and its analogues modulate dopamine receptors, toward the development of new therapeutics to treat Parkinson's disease, tardive dyskinesia and schizophrenia. As a first step toward establishing a validated protocol for accurate computational modeling of PLG and associated peptidomimetic analogues, we evaluated the accuracy of density functional theory (DFT), wavefunction theory (WFT), and molecular mechanics (MM) calculations for PLG and for a library of structurally related small molecules.

Design, synthesis and evaluation of analogs of initiation factor 4E (eIF4E) cap-binding antagonist Bn7-GMP.

Aberrant regulation of cap-dependent translation has been frequently observed in the development of cancer. Association of the cap-binding protein eIF4E with N(7)-methylated guanosine capped mRNA is the rate limiting step governing translation initiation; and therefore represents an attractive process for cancer drug discovery. Previously, replacement of the 7-Me group of the Me(7)-guanosine monophosphate with a benzyl group has been found to increase binding affinity to eIF4E.

Total synthesis and evaluation of C26-hydroxyepothilone D derivatives for photoaffinity labeling of beta-tubulin.

Three photoaffinity labeled derivatives of epothilone D were prepared by total synthesis, using efficient novel asymmetric synthesis methods for the preparation of two important synthetic building blocks. The key step for the asymmetric synthesis of (S,E)-3-(tert-butyldimethylsilyloxy)-4-methyl-5-(2-methylthiazol-4-yl)pent-4-enal involved a ketone reduction with (R)-Me-CBS-oxazaborolidine. For the synthesis of (5S)-5,7-di[(tert-butyldimethylsilyl)oxy]-4,4-dimethylheptan-3-one an asymmetric Noyori reduction of a beta-ketoester was employed.

Identification of novel non-hydroxamate anthrax toxin lethal factor inhibitors by topomeric searching, docking and scoring, and in vitro screening.

Anthrax is an infectious disease caused by Bacillus anthracis, a Gram-positive, rod-shaped, anaerobic bacterium. The lethal factor (LF) enzyme is secreted by B. anthracis as part of a tripartite exotoxin and is chiefly responsible for anthrax-related cytotoxicity.

Energies, Geometries, and Charge Distributions of Zn Molecules, Clusters, and Biocenters from Coupled Cluster, Density Functional, and Neglect of Diatomic Differential Overlap Models.

We present benchmark databases of Zn-ligand bond distances, bond angles, dipole moments, and bond dissociation energies for Zn-containing small molecules and Zn coordination compounds with H, CH3, C2H5, NH3, O, OH, H2O, F, Cl, S, and SCH3 ligands. The test set also includes clusters with Zn-Zn bonds. In addition, we calculated dipole moments and binding energies for Zn centers in coordination environments taken from zinc metalloenzyme X-ray structures, representing both structural and catalytic zinc centers.

Total synthesis and evaluation of 22-(3-azidobenzoyloxy)methyl epothilone C for photoaffinity labeling of beta-tubulin.

The total synthesis of 22-(3-azidobenzoyloxy)methyl epothilone C is described as a potential photoaffinity probe to elucidate the beta-tubulin binding site. A sequential Suzuki-aldol-Yamaguchi macrolactonization strategy was utilized employing a novel derivatized C1-C6 fragment. The C22-functionalized analog exhibited good activity in microtubule assembly assays, but cytotoxicity was significantly reduced.

Assessment of Density Functionals, Semiempirical Methods, and SCC-DFTB for Protonated Creatinine Geometries.

Creatinine concentrations in blood and urine can be used to detect renal insufficiencies and muscle diseases, but current chemical sensors cannot measure creatinine with sufficient selectivity and robustness because they lack a receptor that binds protonated creatinine (creatininium) selectively enough. As a first step toward identifying potential receptors for creatininium, we examine the accuracy of density functional theory (DFT) and wave function theory (WFT) calculations for creatininium cation geometries, evaluated against reference parameters from experiment. We tested twenty-one local and nonlocal density functionals, Hartree-Fock theory, four semiempirical molecular orbital (SEMO) methods of the neglect of differential overlap (NDO) type, and one self-consistent-field nonorthogonal tight-binding method (SCC-DFTB) as implemented in two closely related software packages.

Molecular Modeling of Geometries, Charge Distributions, and Binding Energies of Small, Drug-Like Molecules Containing Nitrogen Heterocycles and Exocyclic Amino Groups in the Gas Phase and Aqueous Solution.

We have tested a variety of approximate methods for modeling 30 systems containing mixtures of nitrogen heterocycles and exocyclic amines, each of which is studied with up to 31 methods in one or two phases (gaseous and aqueous). Fifteen of the systems are protonated, and 15 are not. We consider a data set consisting of geometric parameters, partial atomic charges, and water binding energies for the methotrexate fragments 2-(aminomethyl)pyrazine and 2,4-diaminopyrimidine, as well as their cationic forms 1-2-(aminomethyl)pyrazine and 1-2,4-diaminopyrimidine.

Total synthesis and evaluation of C25-benzyloxyepothilone C for tubulin assembly and cytotoxicity against MCF-7 breast cancer cells.

The total synthesis of C25-benzyloxy epothilone C is described. A sequential Suzuki-Aldol-Yamaguchi macrolactonization strategy was utilized employing a novel derivatized C8-C12 fragment. The C25-benzyloxy analog exhibited significantly reduced biological activity in microtubule assembly and cytotoxicity assays.

Zn Coordination Chemistry: Development of Benchmark Suites for Geometries, Dipole Moments, and Bond Dissociation Energies and Their Use To Test and Validate Density Functionals and Molecular Orbital Theory.

We present nonrelativistic and relativistic benchmark databases (obtained by coupled cluster calculations) of 10 Zn-ligand bond distances, 8 dipole moments, and 12 bond dissociation energies in Zn coordination compounds with O, S, NH3, H2O, OH, SCH3, and H ligands. These are used to test the predictions of 39 density functionals, Hartree-Fock theory, and seven more approximate molecular orbital theories. In the nonrelativisitic case, the M05-2X, B97-2, and mPW1PW functionals emerge as the most accurate ones for this test data, with unitless balanced mean unsigned errors (BMUEs) of 0.

A preliminary in silico lead series of 2-phthalimidinoglutaric acid analogues designed as MMP-3 inhibitors.

Matrix metalloproteinases (MMPs) have been the subject of intense research because of their roles in tumor metastasis and in the rise and spread of degenerative diseases such as osteo- and rheumatoid arthritis. A preliminary class of 140 druglike, small-molecule matrix metalloproteinase-3 inhibitors, intended as starting scaffolds for optimization and synthesis, has been designed in silico using a series of highly predictive three-dimensional quantitative structure-activity relationship models, including comparative molecular field analysis and comparative molecular similarity indices analysis, with docking and scoring. Thalidomide was chosen as the skeleton on which to base the new lead series, as it moderately inhibits MMP-3, is antiangiogenic, and lends itself easily to structural modifications.