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 step 3 and step 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 timeaveraged MD simulation results helped understand the differences in binding modes of related compounds.

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, and 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 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 H-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 on how histone recognition by the bromodomain directs the biological function of BRPF1, ultimately targeting the MOZ HAT complex to chromatin substrates.

No abstract

No abstract

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

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.

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.

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 transport and partitioning. Solvation in strata of the phospholipid bilayer is an important drug characteristics 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.