A fast new approach to pharmacophore mapping and its application to dopaminergic and benzodiazepine agonists

Summary In the absence of a 3D structure of the target biomolecule, to propose the 3D requirements for a small molecule to exhibit a particular bioactivity, one must supply both a bioactive conformation and a superposition rule for every active compound. Our strategy identifies both simultaneously. We first generate and optimize all low-energy conformations by any suitable method. For each conformation we then use ALAD-DIN to calculate the location of points to be considered as part of the superposition. These points include atoms in the molecule and projections from the molecule to hydrogen-bond donors and acceptors or charged groups in the binding site. These positions and the relative energy of each conformation are the input to our new program DISCO. It uses a clique-detection method to find superpositions that contain a least one conformation of each molecule and user-specified numbers of point types and chirality. DISCO is fast; for example, it takes about 1 min CPU to propose pharmacophores from 21 conformations of seven molecules. We typically run DISCO several times to compare alternative pharmacophore maps. For D₂ dopamine agonists DISCO shows that the newer 2-aminothiazoles fit the traditional pharmacophore. Using site points correctly identifies the bioactive enantiomers of indoles to compare with catechols whereas using only ligand points leads to selecting the inactive enantiomer for the pharmacophore map. In addition, DISCO reproduces pharmacophore maps of benzodiazepines in the literature and proposes subtle improvements. Our experience suggests that clique-detection methods will find many applications in computational chemistry and computer-assisted molecular design.     

Applications of simulated annealing to the conformational analysis of flexible molecules

We describe in this article our solution to the global minimum problem which uses the simulated annealing algorithm of Kirkpatrick. This method is a Metropolis (e-ΔE/kT) Monte Carlo sampling of conformation space with simultaneous constraint of the search by lowering the temperature T so that the search converges on the global minimum. The Anneal-Conformer program has been extensively tested with peptides and organic molecules using either the Amber or MM2 force fields. A history file of the simulated annealing process allows reconstruction of the random walk in conformation space for subsequent examination. Thus plots of distance and dihedral angle changes during the search for the global minimum can be examined to deduce molecular shape and flexibility. A separate program Conf-Gen reads the history file and extracts all low energy conformations visited during the run.     

An ant algorithm for the conformational analysis of flexible molecules

Originally, the ant system was developed for optimization in discrete search spaces such as the traveling salesman problem. We detail our adaptation of the algorithm to optimization in the continuous search space of conformational analysis. The parameters of the algorithm were tuned using a simple test molecule, undecane, and a drug molecule, imatinib. The algorithm is further tested on four more drug or drug-like molecules, on vitamin A and on alanine tetrapeptide.     

A new approach to variational CI calculations on systems containing non-interacting molecules

A new method to perform variational CI calculations on systems containing non-interacting molecules is presented.     

The application of Artificial Intelligence to the conformational analysis of strained molecules

An algorithm for predicting the conformations of strained molecules using Artificial Intelligence techniques is described and illustrated with some typical examples.     

A molecular mechanics approach to mapping the conformational space of diaryl and triarylphosphines

A molecular mechanics force field has been developed which accurately reproduces experimental solid state structures and conformer interconversion barriers for a series of sterically congested diaryl and triaryl phosphines and some of their chalcogenide and Cr(CO)5 derivatives.     

A new approach to the coordination of anions,novel polycobalticinium macrocyclic receptor molecules

The syntheses and electrochemistry of four novel polycobalticinium macrocyclic receptor molecules are reported. Preliminary anion coordination studies in one case reveal the receptor to be redox responsive to the binding of a Br⁻ guest anion, shifting the respective cobalticinium reduction wave to more negative potentials.     

RigFit: A new approach to superimposing ligand molecules

If structural knowledge of a receptor under consideration is lacking, drug design approaches focus on similarity or dissimilarity analysis of putative ligands. In this context the mutual ligand superposition is of utmost importance. Methods that are rapid enough to facilitate interactive usage, that allow to process sets of conformers and that enable database screening are of special interest here. The ability to superpose molecular fragments instead of entire molecules has proven to be helpful too. The RigFit approach meets these requirements and has several additional advantages. In three distinct test applications, we evaluated how closely we can approximate the observed relative orientation for a set of known crystal structures, we employed RigFit as a fragment placement procedure, and we performed a fragment-based database screening. The run time of RigFit can be traded off against its accuracy. To be competitive in accuracy with another state-of-the-art alignment tool, with which we compare our method explicitly, computing times of about 6s per superposition on a common day workstation are required. If longer run times can be afforded the accuracy increases significantly. RigFit is part of the flexible superposition software FlexS which can be accessed on the WWW [http://cartan.gmd.de/FlexS].     

Automated clustering of probe molecules from solvent mapping of protein surfaces: new algorithms applied to hot-spot mapping and structure-based drug design

Use of solvent mapping, based on multiple-copy minimization (MCM) techniques, is common in structure-based drug discovery. The minima of small-molecule probes define locations for complementary interactions within a binding pocket. Here, we present improved methods for MCM. In particular, a Jarvis-Patrick (JP) method is outlined for grouping the final locations of minimized probes into physical clusters. This algorithm has been tested through a study of protein-protein interfaces, showing the process to be robust, deterministic, and fast in the mapping of protein "hot spots." Improvements in the initial placement of probe molecules are also described. A final application to HIV-1 protease shows how our automated technique can be used to partition data too complicated to analyze by hand. These new automated methods may be easily and quickly extended to other protein systems, and our clustering methodology may be readily incorporated into other clustering packages.     

Systematic approach to design organic magnetic molecules: strongly coupled diradicals with ethylene coupler

The intramolecular magnetic coupling constant (J) values of diradical systems linked with two monoradicals through a coupler (para-substituted phenyl acetylene (Model I), meta-substituted phenyl acetylene (Model II), ethylene (Model III)) were investigated by unrestricted density functional theory calculations. We divided eight monoradicals into α-group and β-group according to Mulliken spin density values of the connected atoms. The overall trends in the strength of magnetic interactions of diradicals were found to be identical in three different model systems. The diradicals with para-substituted phenyl acetylene coupler resulted in almost twice stronger intramolecular magnetic coupling interactions of the corresponding diradicals as compared to the meta-substituted one with opposite magnetism. NN-Ethylene-PO (nitronyl nitroxide radical coupled to phenoxyl radical via ethylene coupler) was calculated to have the strongest magnetic coupling constant with ferromagnetism, and to be even stronger (more than twice) than NN-ethylene-NN (nitronyl nitroxide diradical with ethylene coupler), which was reported to have strong antiferromagnetic interactions in a previous experiment. It was found that the spin density values of the connected atoms are closely related to the determination of magnetic interactions and J values. The spin states of the ground state in diradical systems were explained by means of the spin alternation rule.     

A new approach to calculations of intensities of vibrational overtone transitions

A new approach to the theory of intensities of vibrational overtone transitions is formulated in terms of vibrational coupling of the molecular ground state to excited electronic states. Model calculations indicate an important role of nuclear geometry of the excited electronic states in determining overtone spectra. It is shown that the observed overtone spectrum of the CH stretching mode in benzene can be reproduced theoretically with the assumption that CH bond lengths in the e_lu electronic state are shortened relative to the ground configuration. A simple rule for qualitative prediction of the overtone spectra for diatomic molecules is proposed.     

A CoMFA analysis with conformational propensity: An attempt to analyze the SAR of a set of molecules with different conformational flexibility using a 3D-QSAR method

CoMFA analysis, a widely used 3D-QSAR method, has limitations to handle a set of SAR data containing diverse conformational flexibility since it does not explicitly include the conformational entropic effects into the analysis. Here, we present an attempt to incorporate the conformational entropy effects of a molecule into a 3D-QSAR analysis. Our attempt is based on the assumption that the conformational entropic loss of a ligand upon making a ligand-receptor complex is small if the ligand in an unbound state has a conformational propensity to adopt an active conformation in a complex state. For a QSAR analysis, this assumption was interpreted as follows: a potent ligand should have a higher conformational propensity to adopt an `active-conformation'-like structure in an unbound state than an inactive one. The conformational propensity value was defined as the populational ratio, N<sub>active</sub>/N<sub>stable</sub>, of the number of energetically stable conformers, N<sub>stable</sub>, to the number of `active-conformation'-like structures, N_active. The latter number was calculated by counting the number of conformers that satisfied the structural parameters deduced from the active conformation. A set of SAR data of imidazoleglycerol phosphate dehydratase inhibitors containing 20 molecules with different conformational flexibility was used as a training set for developing a 3D structure-activity relationship by a CoMFA analysis with the conformational propensity value. This resulted in a cross-validated squared correlation coefficient of the CoMFA model with the conformational propensity value (R ² _cross = 0.640) higher than that of the standard CoMFA model (R ² _cross = 0.431). Then we evaluated the quality of the CoMFA models by predicting the inhibitory activity for a new molecule.     

A topological approach to the interaction of lanthanide shift reagents with the carbonyl group : Its relevance to conformational analysis

A topological approach to the LIS simulation process of compounds bearing the carbonyl functional group is reported. The Lanthanide was always found to complex in the less hindered position, along the sp² lone-pair orbital of the O atom. The relevance to chemistry of such an approach is that only chemical significant positions are considered as feasible sites of complexation of the Lanthanide with the substrate, disregarding all other positions. More important, the results obtained have thrown some light on the factors that may influence the equilibria of conformationally mobile systems, in the presence of LSR.The reliability of the method employed has been tested in several disparate carbonyl derivatives.     

A new approach to the relationship between bond energy and electronegativity

The familiar equation whereby Pauling related heteronuclear bond energies D_A–B to the electronegativity difference Δχ (=∣χ_A − χ_B∣) and the homonuclear bond energies D_A–A and D_B–B has been the subject of critical scrutiny for at least half a century. A modification to this equation that combines the concepts of electronegativity and hardness/softness can be rewritten in terms of two quantities x and y, both having absolute significance. Both homo- and heteronuclear bond energies can be rationalised from these new equations. The quantities x are linearly related to Pauling electronegativities, while y appears to be a measure of an atom’s intrinsic bonding potential, related to size and availability of valence orbitals.     

Binding free energy predictions of farnesoid X receptor (FXR) agonists using a linear interaction energy (LIE) approach with reliability estimation: application to the D3R Grand Challenge 2

Computational protein binding affinity prediction can play an important role in drug research but performing efficient and accurate binding free energy calculations is still challenging. In the context of phase 2 of the Drug Design Data Resource (D3R) Grand Challenge 2 we used our automated eTOX ALLIES approach to apply the (iterative) linear interaction energy (LIE) method and we evaluated its performance in predicting binding affinities for farnesoid X receptor (FXR) agonists. Efficiency was obtained by our pre-calibrated LIE models and molecular dynamics (MD) simulations at the nanosecond scale, while predictive accuracy was obtained for a small subset of compounds. Using our recently introduced reliability estimation metrics, we could classify predictions with higher confidence by featuring an applicability domain (AD) analysis in combination with protein–ligand interaction profiling. The outcomes of and agreement between our AD and interaction-profile analyses to distinguish and rationalize the performance of our predictions highlighted the relevance of sufficiently exploring protein–ligand interactions during training and it demonstrated the possibility to quantitatively and efficiently evaluate if this is achieved by using simulation data only.     

Waste-recycling Monte Carlo with optimal estimates: application to free energy calculations in alloys

The estimator proposed recently by Delmas and Jourdain for waste-recycling Monte Carlo achieves variance reduction optimally with respect to a control variate that is evaluated directly using the simulation data. Here, the performance of this estimator is assessed numerically for free energy calculations in generic binary alloys and is compared to those of other estimators taken from the literature. A systematic investigation with varying simulation parameters of a simplified system, the anti-ferromagnetic Ising model, is first carried out in the transmutation ensemble using path-sampling. We observe numerically that (i) the variance of the Delmas-Jourdain estimator is indeed reduced compared to that of other estimators; and that (ii) the resulting reduction is close to the maximal possible one, despite the inaccuracy in the estimated control variate. More extensive path-sampling simulations involving an FeCr alloy system described by a many-body potential additionally show that (iii) gradual transmutations accommodate the atomic frustrations; thus, alleviating the numerical ergodicity issue present in numerous alloy systems and eventually enabling the determination of phase coexistence conditions.     

A neoceptor approach to unraveling microscopic interactions between the human A2A adenosine receptor and its agonists

Strategically mutated neoceptors, e.g., with anionic residues in TMs 3 and 7 intended for pairing with positively charged amine-modified nucleosides, were derived from the antiinflammatory A(2A) adenosine receptor (AR). Adenosine derivatives functionalized at the 5', 2, and N(6) positions were synthesized. The T88D mutation selectively enhanced the binding of the chain-length-optimized 5'-(2-aminoethyl)uronamide but not 5'-(2-hydroxyethyl)uronamide, suggesting a critical role of the positively charged amine. Combination of this modification with the N(6)-(2-methylbenzyl) group enhanced affinity at the Q89D- and N181D- but not the T88D-A(2A)AR. Amino groups placed near the 2- or N(6)-position only slightly affected the binding to mutant receptors. The 5'-hydrazide MRS3412 was 670- and 161-fold enhanced, in binding and functionally, respectively, at the Q89D-A(2A)AR compared to the wild-type. Thus, we identified and modeled pairs of A(2A)AR-derived neoceptor-neoligand, which are pharmacologically orthogonal with respect to the native species.     

Free energy derivatives: A new method for probing the convergence problem in free energy calculations

The convergence behavior of free energy calculations has been explored in more detail than in any previously reported work, using a model system of two neon atoms in a periodic box of water. We find that for thermodynamic integration-type free energy calculations as much as a nanosecond or more molecular dynamics sampling is required to obtain a fully converged value for a single λ point of the integrand. The concept of “free energy derivatives” with respect to the individual parameters of the force field is introduced. This formalism allows the total convergence of the simulation to be deconvoluted into components. A determination of the statistical “sampling ratio” from these simulations indicates that for window-type free energy calculations carried out in a periodic waterbox of typical size at least 0.6 ps of sampling should be performed at each window (0.7 ps if constraint contributions to the free energy are being determined). General methods to estimate and reduce the error in thermodynamic integration and free energy perturbation calculations are discussed. We show that the difficulty in applying such methods is determining a reliable estimate of the correlation length from a short series of data. © 1994 by John Wiley & Sons, Inc.