Search PubMed⌕ Search

Biomedical subjects

W J Welsh

Publications and source records attributed to W J Welsh.

32 records · Page 2Linked to original sources

Pharmaceutical fingerprinting: evaluation of neural networks and chemometric techniques for distinguishing among same-product manufacturers.

The present study was undertaken to evaluate several computer-based classifiers as potential tools for pharmaceutical fingerprinting by utilizing normalized data obtained from HPLC trace organic impurity patterns. To assess the utility of this approach, samples of L-tryptophan (LT) drug substance were analyzed from commercial production lots of six different manufacturers. The performance of several artificial neural network (ANN) architectures was compared with that of two standard chemometric methods, K-nearest neighbors (KNN) and soft independent modeling of class analogy (SIMCA), as well as with a panel of human experts. The architecture of all three computer-based classifiers was varied with respect to the number of input variables. The ANNs were also optimized with respect to the number of nodes per hidden layer and to the number of hidden layers. A novel preprocessing scheme known as the Window method was devised for converting the output of 899 data entries extracted from each chromatogram into an appropriate input file for the classifiers. Analysis of the test set data revealed that an ANN with 46 inputs (i.e., ANN-46) was superior to all other classifiers evaluated, with 93% of the chromatograms correctly classified. Among the classifiers studied in detail, the order of performance was ANN-46 (93%) > SIMCA-46 (87%) > KNN-46 (85%) = ANN-899 (85%) > "human experts" (83%) > SIMCA-899 (78%) > or = ANN-22 (77%) = KNN-22 (77%) > or = KNN-899 (76%) > SIMCA-22 (73%). These results confirm that ANNs, particularly when used in conjunction with the Window preprocessing scheme, can provide a fast, accurate, and consistent methodology applicable to pharmaceutical fingerprinting. Particular attention was paid to variations in the HPLC patterns of same-manufacturer samples due to differences in LT production lots, HPLC columns, and even run-days to quantify how these factors might hinder correct classifications. The results from these classification studies indicate that the chromatograms evidenced variations across LT manufacturers, across the three HPLC columns and, for one manufacturer, across lots. The extent of column-to-column variations is particularly noteworthy in that all three columns had identical specifications with respect to their stationary-phase characteristics and two of the three columns were from the same vendor.

Chemistry, Pharmaceutical↗

A comparative molecular field analysis study of N-benzylpiperidines as acetylcholinesterase inhibitors.

A series of 1-benzyl-4-[2-(N-benzoylamino)ethyl]piperidine derivatives and of N-benzylpiperidine benzisoxazoles has been investigated using the comparative molecular field analysis (CoMFA) approach. These compounds have been found to inhibit the metabolic breakdown of the neurotransmitter acetylcholine (ACh) by the enzyme acetylcholinesterase (AChE) and hence alleviate memory deficits in patients with Alzheimer's Disease by potentiating cholinergic transmission. Development of the CoMFA model considered two separate alignments: (i) alignment I which emphasized the electrostatic fitting of the subject compounds and (ii) alignment II which emphasized their steric fitting. In addition, the inhibitor compounds were considered both as neutral species and as N-piperidine-protonated species. The resulting 3D-QSAR indicates a strong correlation between the inhibitory activity of these N-benzylpiperidines and the steric and electronic factors which modulate their biochemical activity. A CoMFA model with considerable predictive ability was obtained.

Acetylcholinesterase↗

Molecular modeling of substrate-enzyme reactions for the cysteine protease papain.

AM1 quantum mechanical reaction coordinate (RC) calculations were run to simulate the rate-limiting deacylation (hydrolysis) reaction for a series of para-X-PhC(O)NHCH2-C(Y)-S-papain intermediates, where X = OCH3, CH3, H, Cl, NO2 and Y = O (thioester) or S (dithioester), for which a large body of structural, kinetic, and spectroscopic data is available. Several reaction zones, in particular the so-designated Large Zone and Small Zone, were extracted for these RC simulations from the fully solvated and energy-minimized X-ray crystal structure of papain (pdb9pap) bound to the appropriate substrate moiety. The major structural difference between these two zones was the absence of the oxyanion hole in the latter. For both the thioester and dithioester cases, the calculated Ea value associated with the parent (X = H) acyl-enzyme intermediate was lower by ca. 10 kcal/mol for the Large Zone than for the Small Zone. The magnitude of this difference suggests that the oxyanion hole plays a functional if not essential role in stabilizing the anionic tetrahedral intermediate with the cysteine proteases. The calculated Ea value was lower by ca. 10 kcal/mol for the thioester [-C(O)-S-] than for the corresponding dithioester [-C(S)-S-], in qualitative agreement with kinetic data for this series of substrates which reveal that the specific rate constant for deacylation k3 is ca. 60 times larger for the former. This difference is also consistent with both AM1 and 6-31G* calculations on model intermediates, which indicate that the weaker polarity of the dithioester compared with the thioester [i.e., -C(<--S)-S-versus-C(-->O)-S-] renders the former a much poorer site for nucleophilic attack. The anionic tetrahedral intermediate is energetically more stable for the dithioester than for the corresponding thioester, a finding that is discussed in terms of its kinetic and mechanistic implications. The mode of attack by the H2O nucleophile is "concerted" rather than "sequential" in terms of the order of proton abstraction by His-159 and nucleophilic attack on the acyl-enzyme intermediate. While the presumably key Sthiol . . . N nonbonded contact remained almost constant (ca. 2.90 A) up to formation of the [TS] structure, the substrate torsion angles phi and psi rotated significantly as the hybridization around the reaction site transforms from sp2 to sp3 during formation of the tetrahedral intermediate. The AM1-calculated frontier molecular orbitals for model thioester and dithioester acyl-enzyme intermediates generally associate the HOMOs with the reaction site and the LUMOs with the benzamide moiety. Computer graphics images corroborate our view that, in relation to the Sthiol . . . N interaction, the HOMOs and LUMOs should be identified, respectively, with Sthiol and N rather than the reverse, as suggested by other workers.

Acylation↗

Conformational analysis of the lipophilic antifolate trimetrexate.

The conformational properties of the lipophilic antifolate trimetrexate (TMQ) were calculated and compared to the structurally-analogous prototypical antifolate methotrexate (MTX) using both empirical force-field and AM1 quantum mechanical methods. The conformational preferences of TMQ and MTX are diametrically opposed with respect to the bridge-system set of torsion angles tau 1, tau 2: TMQ prefers gauche, trans while MTX prefers approximately trans, gauche. These predictions are consistent with the observed crystal structures of TMQ (i.e., tau 1 = 79 degrees, tau 2 = 178 degrees) and of DHFR-bound MTX (i.e., tau 1 = -157 degrees, tau 2 = 57 degrees in L. casei). The crystal structure of MTX.4H2O deviates from this pattern with tau 1 closer to cis (i.e., 39 degrees) than the predicted trans, yet this near-cis conformation is driven by intermolecular hydrogen-bonding and electrostatic forces operative in the MTX crystal. As a consequence of these strong intermolecular forces, MTX incurs 1.8 kcal/mole in conformational-strain energy in its crystalline form. In contrast, TMQ experiences virtually no conformational strain in its crystalline form. This disparity is attributed to two distinctions between TMQ and MTX: (i) MTX crystallizes as a zwitterion while TMQ crystallizes as the free base, and (ii) the hydrophilic glutamate tail in MTX is replaced by three lipophilic trimethoxy groups in TMQ. The corresponding conformational-strain energy of DHFR-bound MTX is 2.0 kcal/mole while that of DHFR-bound TMQ is only 0.65 kcal/mole based on the assumption that the latter adopts the same bridge conformation as the former. This cost in conformational-strain energy for TMQ and MTX is paid at the expense of their respective free energies of binding of DHFR. Consequently, the present study offers the possibility of designing a new class of antifolates which are conformationally strain-free when bound to DHFR and thereby more effective as chemotherapeutic agents.

Antineoplastic Agents↗

Experimental and theoretical evidence supporting the role of Gly363 in blood coagulation factor IXa (Gly193 in chymotrypsin) for proper activation of the proenzyme.

Factor IX is the zymogen of the serine protease factor IXa involved in blood coagulation. In addition to a catalytic domain homologous to the chymotrypsin family, it has Ca2+, phospholipid, and factor VIIIa binding regions needed for full biologic activity. We isolated a nonfunctional factor IX protein designated factor IXEagle Rock (IXER) from a patient with hemophilia B. The variant protein is indistinguishable from normal factor IX (IXN) in its migration on sodium dodecyl sulfate-gel electrophoresis, isoelectric point in urea, carbohydrate content and distribution, number of gamma-carboxyglutamic acid residues, and beta-OH aspartic acid content, and in its binding to an anti-IXN monoclonal antibody which has been shown previously to inhibit the interaction of factor VIIIa with factor IXaN. Further, IXER is cleaved to yield a factor IXa-like molecule by factor XIa/Ca2+ at a rate similar to that observed for IXN. However, in contrast to IXaN, IXaER does not bind to antithrombin-III (specific inhibitor of IXaN) and does not catalyze the activation of factor X (substrate) to factor Xa. To identify the mutation in IXER, all eight exons of IXN and IXER gene were amplified by the polymerase chain reaction technique and cloned. A single point mutation (G----T) which results in the replacement of Val for Gly363 in the catalytic domain of IXER was identified. Gly363 in factor IXa corresponds to the universally conserved Gly193 in the active site sequence of the chymotrypsin serine protease family. X-ray crystallographic data in the literature demonstrate a critical role of this Gly in stabilizing the active conformation of chymotrypsin/trypsin in two major ways: 1) in the formation of the substrate binding site; and 2) in the development of the oxyanion hole. Our computer structural data support a concept that the Gly363----Val change prevents the development of the active site conformation in factor IXa such that the substrate binding site and the oxyanion hole are not formed in the mutated enzyme.

Computer Graphics↗

CNDO/2 molecular orbital calculations on the antifolate DAMP and some related species: structural geometries, ring distortions, change distributions and conformational characteristics.

Geometry-optimized CNDO/2 molecular orbital calculations were carried out on 2, 4-diamino-5-(1-adamantyl 1)-6-methyl pyrimidine (DAMP), a potent inhibitor of mammalian dihydrofolate reductase which is now in clinical trials, and on its inactive 5-(1-naphthyl) analogue (DNMP-1). Crystallographic data show that DAMP (as the ethylsulfonate salt) has a severely distorted, N1 protonated, pyrimidine ring and has steric crowding of the 6-methyl and adamantyl hydrogens whereas DNMP-2 (as a methanol complex) has a planar, nonprotonated pyrimidine ring that is nearly perpendicular to the naphthalene ring. The CNDO/2 results largely reproduce the crystal structure geometry and show that the ring distortions in DAMP are initiated by steric conflicts between the adamantyl group and the 4- and 6-substituents on the ring. In DNMP-1, the non-interfering naphthyl ring induces little strain within the pyrimidine ring and the effect of protonation is negligible. Rotation about the bond joining the two ring groups is restricted in DAMP by a broad barrier of ca. 8.0 kcal mol-1, and no conformation was successful in relieving steric conflicts and hence reducing the ring distortions. In DNMP-1, rotation is less hindered overall with a broad region of accessible conformational space and a maximum barrier of ca. 7.2 kcal mol-1 for the coplanar conformation. The electronic charge distributions of DAMP and DNMP-1 are almost identical and protonation is preferred at N1 rather than at N3 by ca. 3.7 kcal mol-1 for both DAMP and DNMP-1. The calculations establish that the present methodology can be useful as a predictive tool with regard to the structure and conformational characteristics of these and related species.

Adamantane↗

Application of a pruning algorithm to optimize artificial neural networks for pharmaceutical fingerprinting.

The present study investigates an application of artificial neural networks (ANNs) for use in pharmaceutical fingerprinting. Several pruning algorithms were applied to decrease the dimension of the input parameter data set. A localized fingerprint region was identified within the original input parameter space from which a subset of input parameters was extracted leading to enhanced ANN performance. The present results confirm that ANNs can provide a fast, accurate, and consistent methodology applicable to pharmaceutical fingerprinting.

Algorithms↗

Evaluation of quantitative structure-activity relationship methods for large-scale prediction of chemicals binding to the estrogen receptor.

Three different QSAR methods, Comparative Molecular Field Analysis (CoMFA), classical QSAR (utilizing the CODESSA program), and Hologram QSAR (HQSAR), are compared in terms of their potential for screening large data sets of chemicals as endocrine disrupting compounds (EDCs). While CoMFA and CODESSA (Comprehensive Descriptors for Structural and Statistical Analysis) have been commercially available for some time, HQSAR is a novel QSAR technique. HQSAR attempts to correlate molecular structure with biological activity for a series of compounds using molecular holograms constructed from counts of sub-structural molecular fragments. In addition to using r2 and q2 (cross-validated r2) in assessing the statistical quality of QSAR models, another statistical parameter was defined to be the ratio of the standard error to the activity range. The statistical quality of the QSAR models constructed using CoMFA and HQSAR techniques were comparable and were generally better than those produced with CODESSA. It is notable that only 2D-connectivity, bond and elemental atom-type information were considered in building HQSAR models. Since HQSAR requires no conformational analysis or structural alignment, it is straightforward to use and lends itself readily to the rapid screening of large numbers of compounds. Among the QSAR methods considered, HQSAR appears to offer many attractive features, such as speed, reproducibility and ease of use, which portend its utility for prioritizing large numbers of potential EDCs for subsequent toxicological testing and risk assessment.

Animals↗