PubMed Health⌕ Search

Biomedical subjects

T Mietzner

Publications and source records attributed to T Mietzner.

5 recordsLinked to original sources

Methodological developments and strategies for a fast flexible superposition of drug-size molecules.

An alternative to experimental high through-put screening is the virtual screening of compound libraries on the computer. In absence of a detailed structure of the receptor protein, candidate molecules are compared with a known reference by mutually superimposing their skeletons and scoring their similarity. Since molecular shape highly depends on the adopted conformation, an efficient conformational screening is performed using a knowledge-based approach. A comprehensive torsion library has been compiled from crystal data stored in the Cambridge Structural Database. For molecular comparison a strategy is followed considering shape associated physicochemical properties in space such as steric occupancy, electrostatics, lipophilicity and potential hydrogen-bonding. Molecular shape is approximated by a set of Gaussian functions not necessarily located at the atomic positions. The superposition is performed in two steps: first by a global alignment search operating on multiple rigid conformations and then by conformationally relaxing the best scored hits of the global search. A normalized similarity scoring is used to allow for a comparison of molecules with rather different shape and size. The approach has been implemented on a cluster of parallel processors. As a case study, the search for ligands binding to the dopamine receptor is given.

Algorithms↗

Three-dimensional quantitative similarity-activity relationships (3D QSiAR) from SEAL similarity matrices.

The program SEAL is suited to describe the electrostatic, steric, hydrophobic, and hydrogen bond donor and acceptor similarity of different molecules in a quantitative manner. Similarity scores AF can be calculated for pairs of molecules, using either a certain molecular property or a sum of weighted properties. Alternatively, their mutual similarity can be derived from distances d or covariances c between SEAL-based property fields that are calculated in a regular grid. For a set of N chemically related molecules, such values form an N x N similarity matrix which can be correlated with biological activities, using either regression analysis and an appropriate variable selection procedure or partial least-squares (PLS) analysis. For the Cramer steroid data set, the test set predictivities (r2pred = 0.53-0.84) of different PLS models, based on a weighted sum of molecular properties, are superior to published results of CoMFA and CoMSIA studies (r2pred = 0.31-0.40), regardless of whether a common alignment or individual, pairwise alignments of all molecules are used in the calculation of the similarity matrices. Training and test set selections have a significant influence on the external predictivities of the models. Although the SEAL similarity score between two molecules is a single number, its value is based on the 3D properties of both molecules. The term 3D quantitative similarity-activity analyses (3D QSiAR) is proposed for approaches which correlate 3D structure-derived similarity matrices with biological activities.

Drug Design↗

Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity.

An alternative approach is reported to compute property fields based on similarity indices of drug molecules that have been brought into a common alignment. The fields of different physicochemical properties use a Gaussian-type distance dependence, and no singularities occur at the atomic positions. Accordingly, no arbitrary definitions of cutoff limits and deficiencies due to different slopes of the fields are encountered. The fields are evaluated by a PLS analysis similar to the CoMFA formalism. Two data sets of steroids binding to the corticosteroid-binding-globulin and thermolysin inhibitors were analyzed in terms of the conventional CoMFA method (Lennard-Jones and Coulomb potential fields) and the new comparative molecular similarity indices analysis (CoMSIA). Models of comparative statistical significance were obtained. Field contribution maps were produced for the different models. Due to cutoff settings in the CoMFA fields and the steepness of the potentials close to the molecular surface, the CoMFA maps are often rather fragmentary and not contiguously connected. This makes their interpretation difficult. The maps obtained by the new CoMSIA approach are superior and easier to interpret. Whereas the CoMFA maps denote regions apart from the molecules where interactions with a putative environment are to be expected, the CoMSIA maps highlight those regions within the area occupied by the ligand skeletons that require a particular physicochemical property important for activity. This is a more significant guide to trace the features that really matter especially with respect to the design of novel compounds.

Models, Molecular↗

A fast and efficient method to generate biologically relevant conformations.

Mutual binding between a ligand of low molecular weight and its macromolecular receptor demands structural complementarity of both species at the recognition site. To predict binding properties of new molecules before synthesis, information about possible conformations of drug molecules at the active site is required, especially if the 3D structure of the receptor is not known. The statistical analysis of small-molecule crystal data allows one to elucidate conformational preferences of molecular fragments and accordingly to compile libraries of putative ligand conformations. A comparison of geometries adopted by corresponding fragments in ligands bound to proteins shows similar distributions in conformations space. We have developed an automatic procedure that generates different conformers of a given ligand. The entire molecule is decomposed into its individual ring and open-chain torsional fragments, each used in a variety of favorable conformations. The latter ones are produced according to the library information about conformational preferences. During this building process, an extensive energy ranking is applied. Conformers ranked as energetically favorable are subjected to an optimization in torsion angle space. During minimization, unfavorable van der Waals interactions are removed while keeping the open-chain torsion angles as close as possible to the experimentally most frequently observed values. In order to assess how well the generated conformers map conformation space, a comparison with experimental data has been performed. This comparison gives some confidence in the efficiency and completeness of this approach. For some ligands that had been structurally characterized by protein crystallography the program was used to generate sets of some 10 to 100 conformers. Among these, geometries are found that fall convincingly close to the conformations actually adopted by these ligands at the binding site.

Adenosine Monophosphate↗

Different approaches toward an automatic structural alignment of drug molecules: applications to sterol mimics, thrombin and thermolysin inhibitors.

A relative comparison of the binding properties of different drug molecules requires their mutual superposition with respect to various alignment criteria. In order to validate the results of different alignment methods, the crystallographically observed binding geometries of ligands in the pocket of a common protein receptor have been used. The alignment function in the program SEAL that calculates the mutual superposition of molecules has been optimized with respect to these references. Across the reference data set, alignments could be produced that show mean rms deviations of approximately 1 A compared to the experimental situation. For structures with obvious skeletal similarities a multiple-flexible fit, linking common pharmacophoric groups by virtual springs, has been incorporated into the molecular mechanics program MOMO. In order to combine conformational searching with comparative alignments, the optimized SEAL approach has been applied to sets of conformers generated by MIMUMBA, a program for conformational analysis. Multiple-flexible fits have been calculated for inhibitors of ergosterol biosynthesis. Sets of different thrombin and thermolysin inhibitors have been conformationally analyzed and subsequently aligned by a combined MIMUMBA/SEAL approach. Since for these examples crystallographic data on their mutual alignment are available, an objective assessment of the computed results could be performed. Among the generated conformers, one geometry could be selected for the thrombin and thermolysin inhibitors that approached reasonably well the experimentally observed alignment.

Amino Acid Sequence↗