Search PubMed⌕ Search

Biomedical subjects

Stefan Balaz

Publications and source records attributed to Stefan Balaz.

12 recordsLinked to original sources
PubMed2007

Drug-membrane interactions studied in phospholipid monolayers adsorbed on nonporous alkylated microspheres.

Characterization of interactions with phospholipids is an integral part of the in vitro profiling of drug candidates because of the roles the interactions play in tissue accumulation and passive diffusion. Currently used test systems may inadequately emulate the bilayer core solvation properties (immobilized artificial membranes [IAM]), suffer from potentially slow transport of some chemicals (liposomes in free or immobilized forms), and require a tedious separation (if used for free liposomes). Here the authors introduce a well-defined system overcoming these drawbacks: nonporous octadecylsilica particles coated with a self-assembled phospholipid monolayer. The coating mimics the structure of the headgroup region, as well as the thickness and properties of the hydrocarbon core, more closely than IAM. The monolayer has a similar transition temperature pattern as the corresponding bilayer. The particles can be separated by filtration or a mild centrifugation. The partitioning equilibria of 81 tested chemicals were dissected into the headgroup and core contributions, the latter using the alkane/water partition coefficients. The deconvolution allowed a successful prediction of the bilayer/water partition coefficients with the standard deviation of 0.26 log units. The plate-friendly assay is suitable for high-throughput profiling of drug candidates without sacrificing the quality of analysis or details of the drug-phospholipid interactions.

Adsorption↗
PubMed2007

Structural determinants of binding of aromates to extracellular matrix: a multi-species multi-mode CoMFA study.

For small molecules acting in tissues, including signaling peptides, effectors, inhibitors, and other drug candidates, nonspecific binding to the extracellular matrix (ECM) is a critical phenomenon affecting their disposition, toxicity, and other effects. A commercially available ECM mimic, forming a solidified layer at the bottom of the vials, was used to measure the association constants of 25 simple aromatic compounds to two forms of ECM proteins, solidified (s-ECM) and dissolved (d-ECM) in the buffer during the incubation. Except for small homologous series, the binding data did not correlate with the lipophilicity and acidity of the compounds, contrary to a common expectation for nonspecific binding. To elucidate the putative structures of averaged binding sites of s-ECM and d-ECM, comparative molecular field analysis (CoMFA) was applied in a modified version taking into consideration that multiple modes and multiple species may be involved. The method shapes a receptor site model from a set of grid points in which the interaction energies between a probe atom and superimposed ligands are calculated. Electrostatic and steric energies in the grid points are characterized by regression coefficients. The forward-selection nonlinear regression analysis was used to optimize the coefficients in the novel multi-species, multi-mode CoMFA models. These models showed satisfactory descriptive and predictive abilities for both s-ECM and d-ECM binding data, which were better than those obtained with the standard, one-mode CoMFA analysis. The calibrated models, defining the electrostatic and van der Waals regions of putative binding sites, are suitable for the prediction of ECM binding for untested chemicals.

Extracellular Matrix↗
PubMed2006

Quantitative characterization of binding of small molecules to extracellular matrix.

Extracellular matrix (ECM) is a major tissue component that, besides its cell support function, is implicated in cell-cell signaling, wound repair, cell adhesion, and other cell and tissue functions. For small molecules acting in tissues, including chemicals, signaling peptides, effectors, inhibitors, and other man-made and physiological compounds, non-specific binding to ECM is a critical phenomenon affecting their disposition. We describe here a method for a quantitative characterization of the ECM binding, using a solidified ECM layer incubated with medium containing studied small molecules. Working conditions of Matrigel, a commercial basement membrane preparation, were optimized in terms of the protein concentration, surface area, gel layer thickness, solidification time, and mixing speed. The release of proteins from the solidified layer into the buffer was monitored and taken into account. Two major proteins, laminin and collagen IV, dissolve at different rates. The Matrigel stability data, obtained under varying incubation conditions and gentle mixing, can also be useful in other ECM-related research. The experimental binding data, averaged over all binding sites, were analyzed assuming a fast linear binding. The binding constants were determined for 10 small organic molecules for both dissolved proteins and the solidified layer. The binding constants tend to increase with lipophilicity of the compounds, as characterized by the 1-octanol/water partition coefficients.

Animals↗
PubMed2006

Partitioning of organic compounds in phases imitating the headgroup and core regions of phospholipid bilayers.

Solvation free energies of drugs, peptides, and other small molecules in the core and headgroup regions of phospholipid bilayers determine their conformations, accumulation, and transport properties. The transfer free energy includes the energy terms for the formation of a cavity for the solute, the interactions of the solute with phospholipids, electrostatic interactions of the solute with the membrane, and dipole potentials and entropy terms. The interaction energies with phospholipids can be estimated by correlating the partitioning in surrogate solvent systems and in the bilayer. As the headgroup surrogate, we use diacetylphosphatidylcholine (DAcPC), the acetylated headgroup of the most abundant mammalian phospholipid, phosphatidylcholine, which forms a homogeneous solution with acceptable viscosity when mixed with water in ratios similar to those in the fully hydrated bilayer. The two-phase system of n-hexadecane (C16) as the core surrogate and hydrated DAcPC was used to monitor partitioning of 16 nonionizable compounds. On the bilogarithmic scale, the C16/DAcPC partition coefficients correlate neither with those in the C16/water and 1-octanol/water systems nor with their difference, which is frequently used as a parameter of hydrogen bonding for prediction of the bilayer location of the solutes. The C16/DAcPC system provides a satisfactory emulation of the solvation properties of the bilayer regions, as reflected in correct predictions of the bilayer location for those of the studied chemicals, for which this information is available.

1-Octanol↗
PubMed2005

A combination of docking, QM/MM methods, and MD simulation for binding affinity estimation of metalloprotein ligands.

To alleviate the problems in the receptor-based design of metalloprotein ligands due to inadequacies in the force-field description of coordination bonds, a four-tier approach was devised. Representative ligand-metalloprotein interaction energies are obtained by subsequent application of (1) docking with metal-binding-guided selection of modes, (2) optimization of the ligand-metalloprotein complex geometry by combined quantum mechanics and molecular mechanics (QM/MM) methods, (3) conformational sampling of the complex with constrained metal bonds by force-field-based molecular dynamics (MD), and (4) a single point QM/MM energy calculation for the time-averaged structures. The QM/MM interaction energies are, in a linear combination with the desolvation-characterizing changes in the solvent-accessible surface areas, correlated with experimental data. The approach was applied to structural correlation of published binding free energies of a diverse set of 28 hydroxamate inhibitors to zinc-dependent matrix metalloproteinase 9 (MMP-9). Inclusion of steps 3 and 4 significantly improved both correlation and prediction. The two descriptors explained 90% of variance in inhibition constants of all 28 inhibitors, ranging from 0.08 to 349 nM, with the average unassigned error of 0.318 log units. The structural and energetic information obtained from the time-averaged MD simulation results helped understand the differences in binding modes of related compounds.

Binding Sites↗
PubMed2005

Processing multimode binding situations in simulation-based prediction of ligand-macromolecule affinities.

The linear response (LR) approximation and similar approaches belong to practical methods for estimation of ligand-receptor binding affinities. The approaches correlate experimental binding affinities with the changes upon binding of the ligand electrostatic and van der Waals energies and of solvation characteristics. These attributes are expressed as ensemble averages that are obtained by conformational sampling of the protein-ligand complex and of the free ligand by molecular dynamics or Monte Carlo simulations. We observed that outliers in the LR correlations occasionally exhibit major conformational changes of the complex during sampling. We treated the situation as a multimode binding case, for which the observed association constant is the sum of the partial association constants of individual states/modes. The resulting nonlinear expression for the binding affinities contains all the LR variables for individual modes that are scaled by the same two to four adjustable parameters as in the one-mode LR equation. The multimode method was applied to inhibitors of a matrix metalloproteinase, where this treatment improved the explained variance in experimental activity from 75% for the unimode case to about 85%. The predictive ability scaled accordingly, as verified by extensive cross-validations.

Computer Simulation↗
PubMed2005

A comparison of the binding sites of matrix metalloproteinases and tumor necrosis factor-alpha converting enzyme: implications for selectivity.

MMPs and TACE (ADAM-17) assume independent, parallel, or opposite pathological roles in cancer, arthritis, and several other diseases. For therapeutic purposes, selective inhibition of individual MMPs and TACE is required in most cases due to distinct roles in diseases and the need to preserve activities in normal states. Toward this goal, we compared force-field interaction energies of five ubiquitous inhibitor atoms with flexible binding sites of 24 known human MMPs and TACE. The results indicate that MMPs 1-3, 10, 11, 13, 16, and 17 have at least one subsite very similar to TACE. S3 subsite is the best target for development of specific TACE inhibitors. Specific binding to TACE compared to most MMPs is promoted by placing a negatively charged ligand part at the bottom of S2 subsite, at the entrance of S1' subsite, or the part of S3' subsite that is close to catalytic zinc. Numerous other clues, consistent with available experimental data, are provided for design of selective inhibitors.

ADAM Proteins↗
PubMed2004

Similarity of binding sites of human matrix metalloproteinases.

Tissue components hydrolyzing matrix metalloproteinases (MMPs) exhibit a high sequence similarity (56-64% in catalytic domains) and yet a significant degree of functional specificity. The hexapeptide-binding sites of 24 known human MMPs were compared in terms of their force field interaction energies with five probes that are most frequently encountered in substrates and inhibitors. The probes moved along a grid enclosing partially flexible binding sites in rigid catalytic domains that were represented by published experimental structures and comparative models and new comparative models for nine most recently characterized MMPs. For individual MMPs, representative interaction energies were obtained as averages for all suitable experimental structures. Correlations of the representative energies for all MMP pairs were succinctly catalogued for individual probes, subsites, and correlation levels. Among the probes (neutral sp(3) carbon and sp(3) oxygen, positive sp(3) nitrogen and hydrogen, and negative carbonyl oxygen), the last probe is least distinctive. Similarities of subsites are decreasing as S1 ' > S2 > S3 ' > S1 approximately S3 > S2 '. Most interesting, occupancies of subsites in published structures of MMP-inhibitor complexes follow an almost parallel trend, alluding to overall low selectivity of known MMP inhibitors. Flexible subsite S1 ' that appears as the specificity pocket in rigid x-ray structures is actually very similar among individual MMPs. Several correlations indicated that MMPs 3, 8, and 12 have similar binding sites. Modeling results are corroborated with published experimental data on MMP inhibition and substrate specificities. The results provide numerous clues for development of specific inhibitors and substrates, as well as for selection of MMPs for testing that provides maximum information without redundant experiments.

Binding Sites↗
PubMed2004

A practical approach to docking of zinc metalloproteinase inhibitors.

Forty zinc-dependent metalloproteinase/ligand complexes with known crystal structures were re-docked using five docking/scoring approaches (DOCK, FlexX, DrugScore, GOLD, and AutoDock). Correct geometry of the coordination bonds between the ligand's zinc binding group (ZBG) and the catalytic zinc is important for docking accuracy and scoring reliability. More than 75% of docked poses with RMSD less than 2A were found to have appropriate ZBG binding, but for poor ZBG binding, about 95% of poses failed to dock correctly. Elimination of poses with inappropriate zinc binding resulted in better binding energy predictions that were further improved by dividing the ligands into subsets according to the ZBG (carboxylates, hydroxamates, and phosphorus containing groups). After a subset re-scoring using the regression functions obtained for individual subsets, DrugScore was able to explain 77% and the consensus scoring scheme X-CSCORE even 88% of variance in binding energies. The approach combining ZBG-based pose selection and subset re-scoring improved the hit rate in virtual screening for metalloproteinase inhibitors for all tested methods by 4-16%.

Algorithms↗
PubMed2002

Subcellular pharmacokinetics and its potential for library focusing.

Subcellular pharmacokinetics (SP) optimizes biology-related factors in the design of libraries for high throughput screening by defining comparatively narrow ranges of properties (lipophilicity, amphiphilicity, acidity, reactivity, 3D-structural features) of the included compounds. The focusing ensures appropriate absorption, distribution, metabolism, excretion, and toxicity (ADMET) in those test biosystems, which are more complex than isolated receptors, and in humans. The SP deploys conceptual models that include transport and accumulation in a series of membranes, protein binding, hydrolysis, and other reactions with cell constituents. The kinetics of drug disposition is described as a non-linear disposition function of drug structure and properties. The SP capabilities are illustrated here using a model-based quantitative structure-activity relationship of toxicity of phenolic compounds against Tetrahymena pyriformis as dependent on lipophilicity and acidity. The resulting SP models clearly outperform empirical models in predictive ability outside the parameter space, as revealed by the leave-extremes-out cross-validation technique with omission of compounds beyond pre-defined lipophilicity and acidity ranges. The SP models do not change substantially if the parameters space is shrunk within some limits. In contrast, the shapes of empirical models vary widely depending upon the fraction of the data set used for their optimization. Once calibrated for a given biosystem, the SP models provide a detailed recipe for tailoring the drug properties to ensure optimum ADMET. The focusing is more accurate than with traditional empirical QSAR studies, assessment of drug-likeness, or the rules for identification of compounds with permeability problems.

Biological Transport↗
PubMedDate not supplied

Multimode ligand binding in receptor site modeling: implementation in CoMFA.

Receptor site modeling methods usually use one binding mode (conformation and/or orientation) for each ligand in a 1:1 complex with receptor. Multiple modes should be considered instead because (1). they have frequently been observed experimentally; (2). in a series, ligands can bind in single yet different modes; and (3). a series may only exhibit one but unknown mode and a few plausible modes must be examined. For multimode binding, the observed ligand/receptor association constant is the sum of the association constants that characterize individual binding modes. This relation, when applied to Comparative Molecular Field Analysis (CoMFA), results in a dependence of the observed binding energy on the probe energies that is nonlinear in optimized parameters. The dependence was linearized to allow parameter optimization by the partial least-squares method that was used iteratively until self-consistency. In addition to the standard CoMFA output, the procedure objectively selects one or a few optimal binding modes out of a dozen or more modes that are considered for each ligand. The approach was applied to published data for binding of 34 polychlorinated dibenzofurans to the aryl hydrocarbon receptor. Descriptive and predictive abilities of the 16-mode model were significantly better than for the one-, two-, and four-mode models. Predominantly, edge-aligned modes were selected that are seldom used in CoMFA. Since inclusion of multimode binding only changes the form of the correlation equation and does not affect the number of optimized parameters, the improvement is believed to be due to a more realistic description.

Algorithms↗