Search PubMed⌕ Search

Biomedical subjects

William J Welsh

Publications and source records attributed to William J Welsh.

45 records · Page 3Linked to original sources

Three-dimensional quantitative structure--property relationship (3D-QSPR) models for prediction of thermodynamic properties of polychlorinated biphenyls (PCBs): enthalpy of sublimation.

Three-dimensional quantitative structure--property relationship (3D-QSPR) models have been constructed using comparative molecular field analysis (CoMFA) to correlate the sublimation enthalpies at 298.15 K of a series of polychlorinated biphenyls (PCBs) with their CoMFA-calculated physicochemical properties. Various alignment schemes, such as atom fit, as is, and inertial were employed in this study. Separate CoMFA models were developed using different partial charge formalisms, namely, electrostatic potential (ESP) and Gasteiger-Marsili (GM) charges. Among the four different CoMFA models constructed for sublimation enthalpy (Delta(sub)H(m)(298.15 K)), the model that combined atom fit alignment and ESP charges yielded the greatest self-consistency (r(2) = 0.976) and internal predictive ability (r(cv)(2) = 0.750). This CoMFA model was used to predict Delta(sub)H(m)(298.15 K) of PCBs for which the corresponding experimental values are unavailable in the literature.

Journal Article↗

Three-dimensional quantitative structure-property relationship (3D-QSPR) models for prediction of thermodynamic properties of polychlorinated biphenyls (PCBs): enthalpy of vaporization.

Three-dimensional Quantitative Structure-Property Relationship (QSPR) models have been derived using Comparative Molecular Field Analysis (CoMFA) to correlate the vaporization enthalpies of a representative set of polychlorinated biphenyls (PCBs) at 298.15 K with their CoMFA-calculated physicochemical properties. Various alignment schemes, such as inertial, as is, and atom fit, were employed in this study. The CoMFA models were also developed using different partial charge formalisms, namely, electrostatic potential (ESP) charges and Gasteiger-Marsili (GM) charges. The most predictive model for vaporization enthalpy (Delta(vap)H(m)(298.15 K)), with atom fit alignment and Gasteiger-Marsili charges, yielded r2 values 0.852 (cross-validated) and 0.996 (conventional). The vaporization enthalpies of PCBs increased with the number of chlorine atoms and were found to be larger for the meta- and para-substituted isomers. This model was used to predict Delta(vap)H(m)(298.15 K) of the entire set of 209 PCB congeners.

Journal Article↗

Three-dimensional quantitative structure-property relationship (3D-QSPR) models for prediction of thermodynamic properties of polychlorinated biphenyls (PCBs): enthalpies of fusion and their application to estimates of enthalpies of sublimation and aqueous solubilities.

Comparative Molecular Field Analysis (CoMFA) has been used to develop three-dimensional quantitative structure-property relationship (3D-QSPR) models for the fusion enthalpy at the melting point (Delta(fus)H(m)(T(fus))) of a representative set of polychlorinated biphenyls (PCBs). Various alignment schemes, such as inertial, as is, atom fit, and field fit, were used in this study to evaluate the predictive capabilities of the models. The CoMFA models have also been derived using partial atomic charges calculated from the electrostatic potential (ESP) and Gasteiger-Marsili (GM) methods. The combination of atom fit alignment and GM charges yielded the greatest self-consistency (r(2) = 0.955) and internal predictive ability (r(cv)(2) = 0.783). This CoMFA model was used to predict Delta(fus)H(m)(T(fus)) of the entire set of 209 PCB congeners, including 193 PCB congeners for which experimental values are unavailable. The CoMFA-predicted values, combined with previous estimations of vaporization and sublimation enthalpies, were used to construct a thermodynamic cycle that validated the internal self-consistency of the predictions for these three thermodynamic properties. The CoMFA-predicted values of fusion enthalpy were also used to calculate aqueous solubilities of PCBs using Mobile Order and Disorder Theory. The agreement between calculated and experimental values of solubility at 298.15 K, characterized by a standard deviation of +/- 0.41 log units, demonstrates the utility of CoMFA-predicted values of fusion enthalpies to calculate aqueous solubilities of PCBs.

Environmental Pollutants↗

Fractal fingerprinting of chromatographic profiles based on wavelet analysis and its application to characterize the quality grade of medicinal herbs.

Extracting chemical fingerprints is an important step for representing and interpreting chromatographic data. In this paper, the chromatographic profile is decomposed into components at different resolution levels using wavelet analysis, then the fractal dimensions of these components are computed as the chemical fingerprints. The chromatographic fingerprint is characterized by the vector composed of these chemical fingerprints, which can represent the chemical patterns of different categories of complex samples. Computer simulations reveal that the fractal fingerprints are more stable than the original chromatographic profile data with respect to variations of peak retention time. To demonstrate the validity of this method, the evaluation of the quality of the medicinal herb Angelica sinensis (Oliv.) diels is investigated. Principal component analysis of the fractal fingerprints indicates that samples belonging to the same quality grade are clustered together, while those belonging to different quality grades are separated. Using these fractal fingerprints taken from the chromatographic scans as inputs for an artificial neural network (ANN). The quality grades of two sets of the herbs were verified by cross-validation, indicating that 96.7% of the herbs are correctly identified with respect to their quality grades evaluated by experienced experts, and 100.0% of the herbs are correctly identified with respect to their quality grades determined by pharmacodynamical evaluation.

Algorithms↗

Identification of a minimal subset of receptor conformations for improved multiple conformation docking and two-step scoring.

Docking and scoring are critical issues in virtual drug screening methods. Fast and reliable methods are required for the prediction of binding affinity especially when applied to a large library of compounds. The implementation of receptor flexibility and refinement of scoring functions for this purpose are extremely challenging in terms of computational speed. Here we propose a knowledge-based multiple-conformation docking method that efficiently accommodates receptor flexibility thus permitting reliable virtual screening of large compound libraries. Starting with a small number of active compounds, a preliminary docking operation is conducted on a large ensemble of receptor conformations to select the minimal subset of receptor conformations that provides a strong correlation between the experimental binding affinity (e.g., Ki, IC50) and the docking score. Only this subset is used for subsequent multiple-conformation docking of the entire data set of library (test) compounds. In conjunction with the multiple-conformation docking procedure, a two-step scoring scheme is employed by which the optimal scoring geometries obtained from the multiple-conformation docking are re-scored by a molecular mechanics energy function including desolvation terms. To demonstrate the feasibility of this approach, we applied this integrated approach to the estrogen receptor alpha (ERalpha) system for which published binding affinity data were available for a series of structurally diverse chemicals. The statistical correlation between docking scores and experimental values was significantly improved from those of single-conformation dockings. This approach led to substantial enrichment of the virtual screening conducted on mixtures of active and inactive ERalpha compounds.

Protein Binding↗

Enrichment of ligands for the serotonin receptor using the Shape Signatures approach.

Shape Signatures, a new 3-dimensional molecular comparison method, has been adapted to rank ligands of the serotonin receptors. A set of 825 agonists and 400 antagonists together with approximately 10,000 randomly chosen compounds from the NCI database were used in this study. Both 1D and 2D Shape Signature databases were created, and enrichment studies were carried out. Results from these studies reveal that the 1D Shape Signature approach is highly efficient in separating agonists from a mixture of molecules which includes compounds randomly selected from the NCI database taken as inactives. It is also equally effective at separating agonists and antagonists from a pool of active ligands for the serotonin receptor. Parallel enrichment studies using 2D shape signatures showed high selectivity with more restricted coverage due to the high specificity of 2D signatures. The influence of conformational variation of the shape signature on enrichment was explored by docking a subset of ligands into the crystal structure of serotonin N-acetyltransferase. Enrichment studies on the resulting "docked" conformations produced only slightly improved results compared with the CORINA-generated conformations.

Computer Simulation↗

Using surrogate modeling in the prediction of fibrinogen adsorption onto polymer surfaces.

We present a Surrogate (semiempirical) Model for prediction of protein adsorption onto the surfaces of biodegradable polymers that have been designed for tissue engineering applications. The protein used in these studies, fibrinogen, is known to play a key role in blood clotting. Therefore, fibrinogen adsorption dictates the performance of implants exposed to blood. The Surrogate Model combines molecular modeling, machine learning and an Artificial Neural Network. This novel approach includes an accounting for experimental error using a Monte Carlo analysis. Briefly, measurements of human fibrinogen adsorption were obtained for 45 polymers. A total of 106 molecular descriptors were generated for each polymer. Of these, 102 descriptors were computed using the Molecular Operating Environment (MOE) software based upon the polymer chemical structures, two represented different monomer types, and two were measured experimentally. The Surrogate Model was developed in two stages. In the first stage, the three descriptors with the highest correlation to adsorption were determined by calculating the information gain of each descriptor. Here a Monte Carlo approach enabled a direct assessment of the effect of the experimental uncertainty on the results. The three highest-ranking descriptors, defined as those with the highest information gain for the sample set, were then selected as the input variables for the second stage, an Artificial Neural Network (ANN) to predict fibrinogen adsorption. The ANN was trained using one-half of the experimental data set (the training set) selected at random. The effect of experimental error on predictive capability was again explored using a Monte Carlo analysis. The accuracy of the ANN was assessed by comparison of the predicted values for fibrinogen adsorption with the experimental data for the remaining polymers (the validation set). The mean value of the Pearson correlation coefficient for the validation data sets was 0.54 +/- 0.12. The average root-mean-square (relative) error in prediction for the validation data sets is 38%. This is an order of magnitude less than the range of experimental values (i.e., 366%) and compares favorably with the average percent relative standard deviation of the experimental measurements (i.e., 17.9%). The effects of each of the user-defined parameters in the ANN were explored. None were observed to have a significant effect on the results. Thus, the Surrogate Model can be used to accurately and unambiguously identify polymers whose fibrinogen absorption is at the limits of the range (i.e., low or high) which is an essential requirement for assessing polymers for regenerative tissue applications.

Adsorption↗

Highly predictive CoMFA and CoMSIA models for two series of stromelysin-1 (MMP-3) inhibitors elucidate S1' and S1-S2' binding modes.

Three-dimensional quantitative structure-activity relationship models have been derived using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) for two training sets of arylsulfonyl isoquinoline-based and thazine/thiazepine-based matrix metalloproteinase inhibitors (MMPIs). The crystal structure of stromelysin-1 (MMP-3) was used to pinpoint areas on the ligands and receptors where steric and electrostatic effects (for CoMFA) and steric, electrostatic, hydrogen-bond donor, hydrogen-bond acceptor, and hydrophobic effects (for CoMSIA) correlate with an increase or decrease in experimental biological activity. The most predictive CoMFA and CoMSIA models were obtained using training-series subsets that sampled a wide range of activities, together with docking and scoring, inertial alignment, investigation of various partial charge formalisms, and manual adjustment of each compound within the active site. The models developed in this study are in agreement with experimentally observed MMP-3 structure-activity relationship data and offer new insights into binding modes involving the partly solvent-exposed S1-S2' subpocket and certain zinc-chelating groups.

Hydrogen Bonding↗

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. Most of the new compounds demonstrate medium to high predicted biological activity and good bioavailability as estimated by the octanol-water partition coefficient ClogP. Compound 102 in particular exhibits extremely favorable predicted activity against MMP-3; is moderately bioavailable; satisfies Lipinski's Rule of Five; and shows promise for further optimization, synthesis, and experimental evaluation as a potential adjunct anticancer or antirheumatic therapeutic.

Binding Sites↗