We describe the methodology, parametrization, and application of a conformational search method, called ConfGen, designed to efficiently generate bioactive conformers. We define efficiency as the ...ability to generate a bioactive conformation within a small total number of conformations using a reasonable amount of computer time. The method combines physics-based force field calculations with empirically derived heuristics designed to achieve efficient searching and prioritization of the ligand’s conformational space. While many parameter settings are supported, four modes spanning a range of speed and quality trades-offs are defined and characterized. The validation set used to test the method is composed of ligands from 667 crystal structures covering a broad array of target and ligand classes. With the fastest mode, ConfGen uses an average of 0.5 s per ligand and generates only 14.3 conformers per ligand, at least one of which lies within 2.0 Å root-mean-squared deviation of the crystal structure for 96% of the ligands. The most computationally intensive mode raises this recovery rate to 99%, while taking 8 s per ligand. Combining multiple search modes to “fill-in” holes in the conformation space or energy minimizing using an all-atom force field each lead to improvements in the recovery rates at higher resolutions. Overall, ConfGen is at least as good as competing programs at high resolution and demonstrates higher efficiency at resolutions sufficient for many downstream applications, such as pharmacophore modeling.
Sampling low energy conformations of macrocycles is challenging due to the large size of many of these molecules and the constraints imposed by the macrocycle. We present a new conformational search ...method (implemented in MacroModel) that uses brief MD simulations followed by minimization and normal-mode search steps. The method was parametrized using a set of 100 macrocycles from the PDB and CSD. It was then tested on a publicly available data set for which there are published results using alternative methods; we found that when the same force field is used (in this case MMFFs in vacuum), our method tended to identify conformations with lower energies than what the other methods identified. The performance on a new set of 50 macrocycles from the PDB and CSD was also quite good; the mean and median RMSD values for just the ring atoms were 0.60 and 0.33 Å, respectively. However, the RMSD values for macrocycles with more than 30 ring-atoms were quite a bit larger compared to the smaller macrocycles. Possible origins for this and ideas for improving the performance on very large macrocycles are discussed.
Chemical investigation of the marine red alga (Rhodophyta) Ceratodictyon spongiosum containing the symbiotic sponge Sigmadocia symbiotica collected from Biaro Island, Indonesia, yielded two isomers ...of a new and bioactive thiazole-containing cyclic heptapeptide, cis,cis-ceratospongamide (1) and trans, trans-ceratospongamide (2). Isolation of these peptides was assisted by bioassay-guided fractionation using a brine shrimp toxicity assay (Artemia salina). The structures of the ceratospongamides, which each consist of two L-phenylalanine residues, one (L-isoleucine)-L-methyloxazoline residue, one L-proline residue, and one (L-proline)thiazole residue, were established through extensive NMR spectroscopy, including (1)H-(13)C HMQC-TOCSY, and (1)H-(15)N HMBC experiments, as well as chemical degradation and chiral analysis. cis,cis- and trans,trans-ceratospongamide are stable conformational isomers of the two proline amide bonds. Molecular modeling of these two ceratospongamide isomers showed the trans, trans isomer to be quite planar, whereas the cis,cis isomer has a more puckered overall conformation. trans,trans-Ceratospongamide exhibits potent inhibition of sPLA(2) expression in a cell-based model for antiinflammation (ED(50) 32 nM), whereas the cis,cis isomer is inactive. trans,trans-Ceratospongamide was also shown to inhibit the expression of a human-sPLA(2) promoter-based reporter by 90%.
Hectochlorin (1) was isolated from marine isolates of Lyngbya majuscula collected from Hector Bay, Jamaica, and Boca del Drago Beach, Bocas del Toro, Panama. The planar structure was deduced by one- ...and two-dimensional NMR spectroscopy. X-ray crystallography was used to determine the absolute stereochemistry of hectochlorin as 2S,3S,14S,22S. Hectochlorin is equipotent to jasplakinolide (5) in its ability to promote actin polymerization, but unlike jasplakinolide, is unable to displace a fluorescent phalloidin analogue from polymerized actin. In addition, hectochlorin shows both a unique profile of cytotoxicity by the COMPARE algorithm and potent inhibitory activity toward the fungus Candida albicans. Structurally, hectochlorin resembles dolabellin and the recently reported lyngbyabellin class of compounds.
Epik version 7 is a software program that uses machine learning for predicting the pKa values and protonation state distribution of complex, druglike molecules. Using an ensemble of atomic graph ...convolutional neural networks (GCNNs) trained on over 42,000 pKa values across broad chemical space from both experimental and computed origins, the model predicts pKa values with 0.42 and 0.72 pKa unit median absolute and root mean square errors, respectively, across seven test sets. Epik version 7 also generates protonation states and recovers 95% of the most populated protonation states compared to previous versions. Requiring on average only 47 ms per ligand, Epik version 7 is rapid and accurate enough to evaluate protonation states for crucial molecules and prepare ultra-large libraries of compounds to explore vast regions of chemical space. The simplicity and time required for the training allow for the generation of highly accurate models customized to a program's specific chemistry.
Epik version 7 is a software program that uses machine learning for predicting the p
values and protonation state distribution of complex, druglike molecules. Using an ensemble of atomic graph ...convolutional neural networks (GCNNs) trained on over 42,000 p
values across broad chemical space from both experimental and computed origins, the model predicts p
values with 0.42 and 0.72 p
unit median absolute and root mean square errors, respectively, across seven test sets. Epik version 7 also generates protonation states and recovers 95% of the most populated protonation states compared to previous versions. Requiring on average only 47 ms per ligand, Epik version 7 is rapid and accurate enough to evaluate protonation states for crucial molecules and prepare ultra-large libraries of compounds to explore vast regions of chemical space. The simplicity and time required for the training allow for the generation of highly accurate models customized to a program's specific chemistry.
Epik version 7 is a software program that uses machine learning for predicting the pK a values and protonation state distribution of complex, druglike molecules. Using an ensemble of atomic graph ...convolutional neural networks (GCNNs) trained on over 42,000 pK a values across broad chemical space from both experimental and computed origins, the model predicts pK a values with 0.42 and 0.72 pK a unit median absolute and root mean square errors, respectively, across seven test sets. Epik version 7 also generates protonation states and recovers 95% of the most populated protonation states compared to previous versions. Requiring on average only 47 ms per ligand, Epik version 7 is rapid and accurate enough to evaluate protonation states for crucial molecules and prepare ultra-large libraries of compounds to explore vast regions of chemical space. The simplicity and time required for the training allow for the generation of highly accurate models customized to a program’s specific chemistry.
Dihydroorotate dehydrogenase (DHODH) has been clinically validated as a target for the development of new antimalarials. Experience with clinical candidate triazolopyrimidine DSM265 (1) suggested ...that DHODH inhibitors have great potential for use in prophylaxis, which represents an unmet need in the malaria drug discovery portfolio for endemic countries, particularly in areas of high transmission in Africa. We describe a structure-based computationally driven lead optimization program of a pyrrole-based series of DHODH inhibitors, leading to the discovery of two candidates for potential advancement to preclinical development. These compounds have improved physicochemical properties over prior series frontrunners and they show no time-dependent CYP inhibition, characteristic of earlier compounds. Frontrunners have potent antimalarial activity in vitro against blood and liver schizont stages and show good efficacy in Plasmodium falciparum SCID mouse models. They are equally active against P. falciparum and Plasmodium vivax field isolates and are selective for Plasmodium DHODHs versus mammalian enzymes.
Tyrosine sulfate-mediated interactions play an important role in HIV-1 entry. After engaging the CD4 receptor at the cell surface, the HIV-1 gp120 glycoprotein binds to the CCR5 co-receptor via an ...interaction that requires two tyrosine sulfates, at positions 10 and 14 in the CCR5-N terminus. Building on previous structure determinations of this interaction, here we report the targeting of these tyrosine sulfate binding sites for drug design through in silico screening of small molecule libraries, identification of lead compounds, and characterization of biological activity. A class of tyrosine sulfate-mimicking small molecules containing a “phenyl sulfonate-linker-aromatic” motif was identified that specifically inhibited binding of gp120 to the CCR5-N terminus as well as to sulfated antibodies that recognize the co-receptor binding region on gp120. The most potent of these compounds bound gp120 with low micromolar affinity and its CD4-induced conformation with K D's as tight as ∼50 nM. Neutralization experiments suggested the targeted site to be conformationally inaccessible prior to CD4 engagement. Primary HIV-1 isolates were weakly neutralized, preincubation with soluble CD4 enhanced neutralization, and engineered isolates with increased dependence on the N terminus of CCR5 or with reduced conformational barriers were neutralized with IC50 values as low as ∼1 μM. These results reveal the potential of targeting the tyrosine sulfate interactions of HIV-1 and provide insight into how mechanistic barriers, evolved by HIV-1 to evade antibody recognition, also restrict small-molecule-mediated neutralization.
Abstract
MDM2 and MDMX function as key regulators of p53 by binding to its N terminus, inhibiting its transcriptional activity, and promoting its degradation. In particular, MDM2 is overexpressed in ...some of human tumors, and with MDMX contributes directly to loss of p53 function during the development of nearly 50% of human cancers. Due to p53 inactivation, MDM2 in many tumors confers tumor survival; therefore it is an important molecular target for anticancer therapy. Several studies showed that reactivation of wild type p53 in tumor cells can be obtained by disrupting the MDM2/p53 interaction with peptidic, peptidomimetic, and small molecule p53-mimetics. Specific successful examples include the Nutlins and spirooxindole analogs (MI-219 and MI-63). Amongst the peptidic and peptidomimetic inhibitors examined to date, none is nearly as effective as Nutlins and MI-219 in tumor killing in vitro. Hence, new inhibitors against MDM2 and/or MDMX are needed: as cell permeable chemical probes of the p53 pathway in cancer biology, and as templates for structure-based rational design of p53 activators for future therapeutic use. As part of our drug discovery program to identify antagonists of the p53/MDM2 and p53/MDMx protein-protein interactions, a high-throughput in-silico screen of a 3.2 millions virtual library of compounds (from Schrödinger, Inc.). A physical restraint was applied during the screen, in order to mimic binding to the hydrophobic cleft of MDM2 normally occupied by three p53 side chains (F19, W23, and L26) that are critical for MDM2/p53 binding. The top highest ranked 160 compounds were then assessed for their ability to block p53 interaction with MDM2 and MDMx in an ELISA assay. This resulted in the identification of E12/DP3-117, a small molecule disruptor of the p53/MDM2 protein-protein interaction with an IC50 value of 47 ± 14 μM. We will report the synthesis and biological evaluation of focused libraries based on the initial hit and on compounds showing improved activity. Structure activity relationship studies around the hits will be disclosed as well as the outcomes of further rounds of chemical design and biological assessment. Binding of E12/DP3-117 to MDM2 is currently being assessed via co-crystallization and other biophysical techniques. We will describe the use of the crystal structure of p53-like mutant peptides in complex with the N-terminal domains of Mdm2, as the basis for rational design of more potent MDM2 small-molecule/peptide hybrid inhibitors.
Citation Format: {Authors}. {Abstract title} abstract. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 3242. doi:10.1158/1538-7445.AM2011-3242