Polarity of Hydrated Phosphatidylcholine Headgroups.

The headgroup (H) stratum (sometimes called the polar region) of membrane bilayers is a relevant yet poorly understood solvation phase for small molecules and macromolecules interacting with the membranes. Solvation of compounds in bilayer strata is characterized experimentally by wide- and small-angle X-ray scattering, neutron diffraction, and various NMR techniques. The quantification is tedious and only available for a limited set of small molecules. Our recently published model of liposome partitioning of small molecules shows that solvation of compounds in the H-stratum of fluid phosphatidylcholine (PC) bilayers correlates well with their solvation in hydrated diacetyl phosphatidylcholine (DAcPC), and solvation in the core (C) depends in a similar way on that in n-hexadecane. These two correlations became a basis for a model describing the location of compounds in the H- and C-strata and at the connecting interface as a nonlinear function of the fragment solvation characteristics of the compounds. In this study, refractivity of hydrated DAcPC phases with varying water contents was measured and polarity was determined using the steady-state fluorescence of indole and Nile Red. The results were compared with the published data obtained by other techniques for PC bilayers in liposomes or on solid supports. The demonstrated qualitative agreement, as well as the polarity and refractivity dependencies on the DAcPC concentration, supports the suitability of hydrated DAcPC as the H-stratum surrogate. Interestingly, depending on hydrations typical for the H-strata of fluid PC bilayers, the dielectric constant could decrease significantly from 31.0 to 7.3 for 16 and 8 water molecules per headgroup, respectively. Although additional experiments are needed for confirmation, this observation could help set proper dielectric constant magnitudes in continuum-based computational models of accumulation and crossing of the PC bilayers with varying hydration levels thanks to the temperature or the structure of fatty acid chains.

Mechanism of Histone H3K4me3 Recognition by the Plant Homeodomain of Inhibitor of Growth 3.

Aberrant access to genetic information disrupts cellular homeostasis and can lead to cancer development. One molecular mechanism that regulates access to genetic information includes recognition of histone modifications, which is carried out by protein modules that interact with chromatin and serve as landing pads for enzymatic activities that regulate gene expression. The ING3 tumor suppressor protein contains a plant homeodomain (PHD) that reads the epigenetic code via recognition of histone H3 tri-methylated at lysine 4 (H3K4me3), and this domain is lost or mutated in various human cancers. However, the molecular mechanisms targeting ING3 to histones and the role of this interaction in the cell remain elusive. Thus, we employed biochemical and structural biology approaches to investigate the interaction of the ING3 PHD finger (ING3PHD) with the active transcription mark H3K4me3. Our results demonstrate that association of the ING3PHD with H3K4me3 is in the sub-micromolar range (KD ranging between 0.63 and 0.93 µm) and is about 200-fold stronger than with the unmodified histone H3. NMR and computational studies revealed an aromatic cage composed of Tyr-362, Ser-369, and Trp-385 that accommodate the tri-methylated side chain of H3K4. Mutational analysis confirmed the critical importance of Tyr-362 and Trp-385 in mediating the ING3PHD-H3K4me3 interaction. Finally, the biological relevance of ING3PHD-H3K4me3 binding was demonstrated by the failure of ING3PHD mutant proteins to enhance ING3-mediated DNA damage-dependent cell death. Together, our results reveal the molecular mechanism of H3K4me3 selection by the ING3PHD and suggest that this interaction is important for mediating ING3 tumor suppressive activities.

Response to "comment on 'structural determinants of drug partitioning in surrogates of phosphatidylcholine bilayer strata'".

We used the solvatochromic correlation to explain the influence of characteristics of studied compounds on the partition coefficients (P) measured using n-hexadecane (C16) and the novel headgroup surrogate (diacetyl phosphatidylcholine, DAcPC), and compared them with those in other systems, including the C16/water (W) system. The comment analyzes why our correlation for the C16/W system has the standard deviation (SD) higher than that published previously. The main reason is that in our, much smaller, data set the measured P values are complemented by the P values predicted by a reliable, unrelated method. We believe that this approach is acceptable for the aforementioned comparison. We did not use just experimental values, as suggested in the comment, because the solvatochromic correlation, although exhibiting 35% reduction in the SD, was accompanied by a sign change of one of the regression coefficients. The recommended use of special solvatochromic solute characteristics for a few compounds and replacement of a predicted PC16/W value by the experimental value resulted in improved correlations. The observed differences between our correlation and those published in the comment and in a previous article do not affect our main conclusions regarding the solvation of solutes in the surrogates (DAcPC and C16) of intrabilayer strata.

Structure-based prediction of drug distribution across the headgroup and core strata of a phospholipid bilayer using surrogate phases.

Solvation of drugs in the core (C) and headgroup (H) strata of phospholipid bilayers affects their physiological transport rates and accumulation. These characteristics, especially a complete drug distribution profile across the bilayer strata, are tedious to obtain experimentally, to the point that even simplified preferred locations are only available for a few dozen compounds. Recently, we showed that the partition coefficient (P) values in the system of hydrated diacetyl phosphatidylcholine (DAcPC) and n-hexadecane (C16), as surrogates of the H- and C-strata of the bilayer composed of the most abundant mammalian phospholipid, PC, agree well with the preferred bilayer location of compounds. High P values are typical for lipophiles accumulating in the core, and low P values are characteristic of cephalophiles preferring the headgroups. This simple pattern does not hold for most compounds, which usually have more even distribution and may also accumulate at the H/C interface. To model complete distribution, the correlates of solvation energies are needed for each drug state in the bilayer: (1) for the H-stratum it is the DAcPC/W P value, calculated as the ratio of the C16/W and C16/DAcPC (W for water) P values; (2) for the C-stratum, the C16/W P value; (3) for the H/C interface, the P values for all plausible molecular poses are characterized using the fragment DAcPC/W and C16/W solvation parameters for the parts of the molecule embedded in the H- and C-strata, respectively. The correlates, each scaled by two Collander coefficients, were used in a nonlinear, mass-balance based model of intrabilayer distribution, which was applied to the easily measurable overall P values of compounds in the DMPC (M = myristoyl) bilayers and monolayers as the dependent variables. The calibrated model for 107 neutral compounds explains 94% of experimental variance, achieves similar cross-validation levels, and agrees well with the nontrivial, experimentally determined bilayer locations for 27 compounds. The resulting structure-based prediction system for intrabilayer distribution will facilitate more realistic modeling of passive transport and drug interactions with those integral membrane proteins, which have the binding sites located in the bilayer, such as some enzymes, influx and efflux transporters, and receptors. If only overall bilayer accumulation is of interest, the 1-octanol/W P values suffice to model the studied set.

Structural determinants of drug partitioning in surrogates of phosphatidylcholine bilayer strata.

The knowledge of drug concentrations in bilayer headgroups, core, and at the interface between them is a prerequisite for quantitative modeling of drug interactions with many membrane-bound transporters, metabolizing enzymes and receptors, which have the binding sites located in the bilayer. This knowledge also helps understand the rates of trans-bilayer transport because balanced interactions of drugs with the bilayer strata lead to high rates, while excessive affinities for any stratum cause a slowdown. Experimental determination of bilayer location is so tedious and costly that the data are only available for some fifty compounds. To extrapolate these valuable results to more compounds at a higher throughput, surrogate phases have been used to obtain correlates of the drug affinities for individual strata. We introduced a novel system, consisting of a diacetyl phosphatidylcholine (DAcPC) solution with the water content of the fluid bilayer as the headgroup surrogate and n-hexadecane (C16) representing the core. The C16/DAcPC partition coefficients were measured for 113 selected compounds, containing structural fragments that are frequently occurring in approved drugs. The data were deconvoluted into the ClogP-based fragment solvation characteristics and processed using a solvatochromic correlation. Increased H-bond donor ability and excess molar refractivity of compounds promote solvation in the DAcPC phase as compared to bulk water, contrary to H-bond acceptor ability, dipolarity/polarizability, and volume. The results show that aromates have more balanced distribution in bilayer strata, and thus faster trans-bilayer transport, than similar alkanes. This observation is in accordance with the frequent occurrence of aromatic rings in approved drugs and with the role of rigidity of drug molecules in promoting intestinal absorption. Bilayer locations, predicted using the C16/DAcPC system, are in excellent agreement with available experimental data, in contrast to other surrogate systems.

Structural determinants of drug partitioning in n-hexadecane/water system.

Surrogate phases have been widely used as correlates for modeling transport and partitioning of drugs in biological systems, taking advantage of chemical similarity between the surrogate and the phospholipid bilayer as the elementary unit of biological phases, which is responsible for most of the transport and partitioning. Solvation in strata of the phospholipid bilayer is an important drug characteristic because it affects the rates of absorption and distribution, as well as the interactions with the membrane proteins having the binding sites located inside the bilayer. The bilayer core can be emulated by n-hexadecane (C16), and the headgroup stratum is often considered a hydrophilic phase because of the high water content. Therefore, we tested the hypothesis that the C16/water partition coefficients (P) can predict the bilayer locations of drugs and other small molecules better than other surrogate systems. Altogether 514 PC16/W values for nonionizable (458) and completely ionized (56) compounds were collected from the literature or measured, when necessary. With the intent to create a fragment-based prediction system, the PC16/W values were factorized into the fragment solvation parameters (f) and correction factors based on the ClogP fragmentation scheme. A script for the PC16/W prediction using the ClogP output is provided. To further expand the prediction system and reveal solvation differences, the fC16/W values were correlated with their more widely available counterparts for the 1-octanol/water system (O/W) using solvatochromic parameters. The analysis for 50 compounds with known bilayer location shows that the available and predicted PC16/W and PO/W values alone or the PC16/O values representing their ratio do not satisfactorily predict the preference for drug accumulation in bilayer strata. These observations indicate that the headgroups stratum, albeit well hydrated, does not have solvation characteristics similar to water and is also poorly described by the O/W partition characteristics.

Rigorous treatment of multispecies multimode ligand-receptor interactions in 3D-QSAR: CoMFA analysis of thyroxine analogs binding to transthyretin.

For a rigorous analysis of the receptor-ligand binding, speciation of the ligands caused by ionization, tautomerism, covalent hydration, and dynamic stereoisomerism needs to be considered. Each species may bind in several orientations or conformations (modes), especially for flexible ligands and receptors. A thermodynamic description of the multispecies (MS), multimode (MM) binding events shows that the overall association constant is equal to the weighted sum of the sums of microscopic association constants of individual modes for each species, with the weights given by the unbound fractions of individual species. This expression is a prerequisite for a precise quantitative characterization of the ligand-receptor interactions in both structure-based and ligand-based structure-activity analyses. We have implemented the MS-MM correlation expression into the comparative molecular field analysis (CoMFA), which deduces a map of the binding site from structures and binding affinities of a ligand set, in the absence of experimental structural information on the receptor. The MS-MM CoMFA approach was applied to published data for binding to transthyretin of 28 thyroxine analogs, each forming up to four ionization species under physiological conditions. The published X-ray structures of several analogs, exhibiting multiple binding modes, served as templates for the MS-MM superposition of thyroxine analogs. Additional modes were generated for compounds with flexible alkyl substituents, to identify bound conformations. The results demonstrate that the MS-MM modification improved predictive abilities of the CoMFA models, even for the standard procedure with MS-MM selected species and modes. The predicted prevalences of individual modes and the generated receptor site model are in reasonable agreement with the available X-ray data. The calibrated model can help in the design of inhibitors of transthyretin amyloid fibril formation.

Biodegradation mechanism of biphenyl by a strain of Pseudomonas stutzeri.

The mechanism of biphenyl biodegradation by Pseudomonas stutzeri was studied. Growth curves on biphenyl were measured, along with dissolution kinetics of biphenyl and production of biosurfactants by the bacterium. Biphenyl was supplied as pure crystals (the crystal biphenyl), adsorbed to Chromosorb G (the Chromosorb G biphenyl) or XAD-4 particles (the XAD-4 biphenyl). No lag phase was observed on the crystal biphenyl, while this period lasted 3.5-6.0 hours on the Chromosorb G biphenyl. The linear specific growth rates (LSGRs) ranged from 2.00 × 10(-4) to 293 × 10(-4) dry weight grams/L/h (d.w.g/L/h) on these two substrates. The LSGR values were directly proportional to the interfacial area between the solid substrate and the microbial suspension. The XAD-4 biphenyl was not bioavailable to the bacterium. The zero-order rates of dissolution ranged from 2.3 × 10(-5)to 8.0 × 10(-4) g/L/h on the crystal biphenyl and the Chromosorb G biphenyl. No biosurfactant production was observed on any biphenyl substrate. Growth curves, results of dissolution measurements and the scanning electron microscope (SEM) images indicate that Pseudomonas stutzeri takes up biphenyl from the aqueous phase as well as the solid phase of the crystal and the Chromosorb G biphenyl. The mechanism of uptake depends on the initial biphenyl concentration and the particle diameter of the biphenyl substrate.

QM/MM linear response method distinguishes ligand affinities for closely related metalloproteins.

Design of selective ligands for closely related targets is becoming one of the most important tasks in the drug development. New tools, more precise than fast scoring functions and less demanding than sophisticated Free Energy Perturbation methods, are necessary to help accomplish this goal. The methods of intermediate complexity, characterizing individual contributions to the binding energy, have been an area of intense research in the past few years. Our recently developed quantum mechanical/molecular mechanical (QM/MM) modification of the Linear Response (LR) method describes the binding free energies as the sum of empirically weighted contributions of the QM/MM interaction energies and solvent-accessible surface areas for the time-averaged structures of hydrated complexes, obtained by molecular dynamics (MD) simulations. The method was applied to published data on 27 inhibitors of matrix metalloproteinase-3 (MMP-3). The two descriptors explained 90% of variance in the inhibition constants with RMSE of 0.245 log units. The QM/MM treatment is indispensable for characterization of the systems lacking suitable force-field expressions. In this case, it provided characteristics of H-bonds of the inhibitors to Glu202, charges of binding site atoms, and accurate coordination geometries of the ligands to catalytic zinc. The geometries were constrained during the MD simulations, which characterized conformational flexibility of the complexes and helped in the elucidation of the binding differences for related compounds. A comparison of the presented QM/MM LR results with those previously published for inhibition of MMP-9 by the same set of ligands showed that the QM/MM LR approach was able to distinguish subtle differences in binding affinities for MMP-3 and MMP-9, which did not exceed one order of magnitude. This precision level makes the approach a useful tool for design of selective ligands to similar targets, because the results can be safely extrapolated to maximize selectivity.

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.

Scalable Synthesis and Purification of Acetylated Phosphatidyl Choline Headgroup.

The acetylated headgroup of the most abundant mammalian phospholipid, 1,2-diacetyl-3-sn-phosphatidyl choline (DAcPC), has several important applications in research. For instance, it can be dissolved in the same amount of water as in the fluid PC bilayer, to create a surrogate of a PC headgroup stratum, for studying solvation of small molecules and the influence of their structure on the process. In contrast to PC derivatives with longer acyl chains, DAcPC does not self-aggregate, rendering the aqueous solution homogeneous and suitable for simplified analyses of interactions of molecules with the headgroups. Several studies have been published where DAcPC was used in a crudely purified form. Here we describe a one-step preparation of DAcPC from commercially available bulk chemicals and purification of the product by crystallization and washing. The process gives a good yield and is easily scalable. The availability of enantiopure, crystalline DAcPC could open the door to more extensive biochemical, pharmacological, and nutritional studies of this interesting chemical.

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.

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.

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.

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.

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%.

Sorption of hydrophobic organic compounds on natural sorbents and organoclays from aqueous and non-aqueous solutions: a mini-review.

Renewed focus on the sorption of hydrophobic organic chemicals (HOCs) onto mineral surfaces and soil components is required due to the increased and wider range of organic pollutants being released into the environment. This mini-review examines the possibility of the contribution and mechanism of HOC sorption onto clay mineral sorbents such as kaolinite, and soil organic matter and the possible role of both in the prevention of environmental contamination by HOCs. Literature data indicates that certain siloxane surfaces can be hydrophobic. Therefore soils can retain HOCs even at low soil organic levels and the extent will depend on the structure of the pollutant and the type and concentration of clay minerals in the sorbent. Clay minerals are wettable by nonpolar solvents and so sorption of HOCs onto them from aqueous and non-aqueous solutions is possible. This is important for two reasons: firstly, the movement and remediation of soil environments will be a function of the concentration and type of clay minerals in the soil. Secondly, low-cost sorbents such as kaolinite and expandable clays can be added to soils or contaminated environments as temporary retention barriers for HOCs. Inorganic cations sorbed onto the kaolinite have a strong influence on the rate and extent of sorption of hydrophobic organic pollutants onto kaolinite. Structural sorbate classes that can be retained by the kaolinite matrix are limited by hydrogen bonding between hydroxyl groups of the octahedral alumosilicate sheet and the tetrahedral sheet with silicon. Soil organic carbon plays a key role in the sorption of HOCs onto soils, but the extent will be strongly affected by the structure of the organic soil matter and the presence of soot. Structural characterisation of soil organic matter in a particular soil should be conducted during a particular contamination event. Contamination by mining extractants and antibiotics will require renewed focus on the use of the QSAR approaches in the context of the sorption of HOCs onto clay minerals from aqueous and non-aqueous solutions.

Molecular insights into the recognition of N-terminal histone modifications by the BRPF1 bromodomain.

The monocytic leukemic zinc finger (MOZ) histone acetyltransferase (HAT) acetylates free histones H3, H4, H2A, and H2B in vitro and is associated with up-regulation of gene transcription. The MOZ HAT functions as a quaternary complex with the bromodomain-PHD finger protein 1 (BRPF1), inhibitor of growth 5 (ING5), and hEaf6 subunits. BRPF1 links the MOZ catalytic subunit to the ING5 and hEaf6 subunits, thereby promoting MOZ HAT activity. Human BRPF1 contains multiple effector domains with known roles in gene transcription, as well as chromatin binding and remodeling. However, the biological function of the BRPF1 bromodomain remains unknown. Our findings reveal novel interactions of the BRPF1 bromodomain with multiple acetyllysine residues on the N-terminus of histones and show that it preferentially selects for H2AK5ac, H4K12ac, and H3K14ac. We used chemical shift perturbation data from NMR titration experiments to map the BRPF1 bromodomain ligand binding pocket and identified key residues responsible for coordination of the post-translationally modified histones. Extensive molecular dynamics simulations were used to generate structural models of bromodomain-histone ligand complexes, to analyze hydrogen bonding and other interactions, and to calculate the binding free energies. Our results outline the molecular mechanism driving binding specificity of the BRPF1 bromodomain for discrete acetyllysine residues on the N-terminal histone tails. Together, these data provide insights into how histone recognition by the bromodomain directs the biological function of BRPF1, ultimately targeting the MOZ HAT complex to chromatin substrates.

Rigorous incorporation of tautomers, ionization species, and different binding modes into ligand-based and receptor-based 3D-QSAR methods.

Speciation of drug candidates and receptors caused by ionization, tautomerism, and/or covalent hydration complicates ligandand receptor-based predictions of binding affinities by 3-dimensional structure-activity relationships (3D-QSAR). The speciation problem is exacerbated by tendency of tautomers to bind in multiple conformations or orientations (modes) in the same binding site. New forms of the 3D-QSAR correlation equations, capable of capturing this complexity, can be developed using the time hierarchy of all steps that lie behind the monitored biological process - binding, enzyme inhibition or receptor activity. In most cases, reversible interconversions of individual ligand and receptor species can be treated as quickly established equilibria because they are finished in a small fraction of the exposure time that is used to determine biological effects. The speciation equilibria are satisfactorily approximated by invariant fractions of individual ligand and receptor species for buffered experimental or in vivo conditions. For such situations, the observed drug-receptor association constant of a ligand is expressed as the sum of products, for each ligand and receptor species pair, of the association microconstant and the fractions of involved species. For multiple binding modes, each microconstant is expressed as the sum of microconstants of individual modes. This master equation leads to new 3D-QSAR correlation equations integrating the results of all molecular simulations or calculations, which are run for each ligand-receptor species pair separately. The multispecies, multimode 3D-QSAR approach is illustrated by a ligand-based correlation of transthyretin binding of thyroxine analogs and by a receptor-based correlation of inhibition of MK2 by benzothiophenes and pyrrolopyrimidines.