Synthetic site-directed ligands

Complexes of nucleotides, peptides and aromatic hapten-like compounds with immunoglobulin fragments were studied by X-ray analysis. After tri- or hexanucleotides of deoxythymidylate were diffused into triclinic crystals of a Fab (BV04-01) with specificity for single-stranded DNA, extensive changes were detected throughout the structure of the protein. The Fab co-crystallized with a tri- or pentanucleotide in a different space group (monoclinic), an observation sometimes correlated with alterations in the structure of the 'native' protein. Structural analyses of the co-crystals are in progress for direct comparisons with the unliganded Fab. In crystals of a human (Mcg) Bence-Jones dimer, synthetic opioid peptides, chemotactic peptides or dinitrophenyl (DNP) derivatives could be diffused into a large conical binding cavity. The conformations of both the ligand and the protein were usually altered during the binding process. At the base of the cavity tyrosine residues could be displaced like trap-doors to permit entry of some opioid peptides and DNP compounds into a deep binding pocket. In co-crystals of the dimer and bis(DNP)lysine, two ligand molecules were bound in tandem, one in the main cavity and the second in the deep pocket. One ligand adopted an extended conformation, with the epsilon-DNP ring near the floor of the main cavity and the alpha-DNP group in solvent outside the binding site. There were no significant conformational changes in the protein. In contrast, the second ligand was very compact, with DNP rings immersed in the deep pocket, and the binding site was expanded to accommodate the oversized ligand. Peptides designed to be specific for the main cavity were incrementally constructed from minimal binding units by M. Geysen, G. Trippick, S. Rodda and their colleagues. A pentapeptide optimized for binding by this method was diffused into a crystal of the dimer and found by Fourier difference analysis to lodge exclusively in the main cavity as predicted. Binding regions in the BV04-01 Fab and the Mcg dimer were markedly different in size and shape. The Fab had a groove-type site, in which a layer of sidechains acted like a false floor over regions analogous to the cavity and deep pocket of the Bence-Jones dimer.     

Structural basis for lower lysine methylation state-specific readout by MBT repeats of L3MBTL1 and an engineered PHD finger

Human L3MBTL1, which contains three malignant brain tumor (MBT) repeats, binds monomethylated and dimethylated lysines, but not trimethylated lysines, in several histone sequence contexts. In crystal structures of L3MBTL1 complexes, the monomethyl- and dimethyllysines insert into a narrow and deep cavity of aromatic residue-lined pocket 2, while a proline ring inserts into shallower pocket 1. We have also engineered a single Y to E substitution within the aromatic cage of the BPTF PHD finger, resulting in a reversal of binding preference from trimethyl- to dimethyllysine in an H3K4 sequence context. In both the "cavity insertion" (L3MBTL1) and "surface groove" (PHD finger) modes of methyllysine recognition, a carboxylate group both hydrogen bonds and ion pairs to the methylammonium proton. Our structural and binding studies of these two modules provide insights into the molecular principles governing the decoding of lysine methylation states, thereby highlighting a methylation state-specific layer of histone mark readout impacting on epigenetic regulation.     

Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin

We report the results of an extended molecular dynamics simulation on the migration of photodissociated carbon monoxide in wild-type sperm whale myoglobin. Our results allow following one possible ligand migration dynamics from the distal pocket to the Xe1 cavity via a path involving the other xenon binding cavities and momentarily two additional packing defects along the pathway. Comparison with recent time resolved structural data obtained by Laue crystallography with subnanosecond to millisecond resolution shows a more than satisfactory agreement. In fact, according to time resolved crystallography, CO, after photolysis, can occupy the Xe1 and Xe4 cavities. However, no information on the trajectory of the ligand from the distal pocket to the Xe1 is available. Our results clearly show one possible path within the protein. In addition, although our data refer to a single trajectory, the local dynamics of the ligand in each cavity is sufficiently equilibrated to obtain local structural and thermodynamic information not accessible to crystallography. In particular, we show that the CO motion and the protein fluctuations are strictly correlated: free energy calculations of the migration between adjacent cavities show that the migration is not a simple diffusion but is kinetically or thermodynamically driven by the collective motions of the protein; conversely, the protein fluctuations are influenced by the ligand in such a way that the opening/closure of the passage between adjacent cavities is strictly correlated to the presence of CO in its proximity. The compatibility between time resolved crystallographic experiments and molecular dynamics simulations paves the way to a deeper understanding of the role of internal dynamics and packing defects in the control of ligand binding in heme proteins.     

The association between a negatively charged ligand and the electronegative binding pocket of its receptor.

Many examples exist of charged amino acids that play a role in attracting or holding a charged ligand toward or inside an oppositely charged binding pocket of the protein. For example, the enzymes superoxide dismutase, triose-phosphate isomerase, and acetylcholinesterase can steer ligands toward their oppositely charged binding pockets or gorges. Interestingly, in our Brownian dynamics simulations of a phosphate-binding protein, we discovered that negatively charged phosphate (HPO(2-)(4)) could make its way into the negatively charged binding pocket. In fact, the phosphate-binding protein exhibits counterintuitive kinetics of association. That is, one would expect that the rate of association would increase on increases to the ionic strength since the interaction between the ligand, with a charge of -2, and the electronegative binding pocket would be repulsive and greater screening should reduce this repulsion and increase the rate of association. However, the opposite is seen-i.e., the rate of association decreases on increases in the ionic strength. We used Brownian dynamics techniques to compute the diffusion limited association rate constants between the negatively charged phosphate ligand and several open forms of PBP (wild-type and several mutants based on an x-ray structure of open-form PBP, mutant T141D). With the appropriate choices of reaction criteria and molecular parameters, the ligand was able to diffuse into the binding pocket. A number of residues influence binding of the ligand within the pocket via hydrogen bonds or salt bridges. Arg135 partially neutralizes the charges on the HPO(2-)(4) ligand in the binding pocket, allowing it to enter. It is also found that the positive electrostatic patches above and below the binding entrance of PBP contribute the major attractive forces that direct the ligand toward the surface of the protein near the binding site.     

The tumour suppressor L(3)mbt inhibits neuroepithelial proliferation and acts on insulator elements

In Drosophila, defects in asymmetric cell division often result in the formation of stem-cell-derived tumours. Here, we show that very similar terminal brain tumour phenotypes arise through a fundamentally different mechanism. We demonstrate that brain tumours in l(3)mbt mutants originate from overproliferation of neuroepithelial cells in the optic lobes caused by derepression of target genes in the Salvador-Warts-Hippo (SWH) pathway. We use ChIP-sequencing to identify L(3)mbt binding sites and show that L(3)mbt binds to chromatin insulator elements. Mutating l(3)mbt or inhibiting expression of the insulator protein gene mod(mdg4) results in upregulation of SWH pathway reporters. As l(3)mbt tumours are rescued by mutations in bantam or yorkie or by overexpression of Expanded, the deregulation of SWH pathway target genes is an essential step in brain tumour formation. Therefore, very different primary defects result in the formation of brain tumours, which behave quite similarly in their advanced stages.     

Exploring the molecular basis of action of ring D aromatic steroidal antiestrogens

Salpichrolides are natural plant steroids that contain an unusual six‐membered aromatic ring D. We recently reported that some of these compounds, and certain analogs with a simplified side chain, exhibited antagonist effects toward the human estrogen receptor (ER), a nuclear receptor whose endogenous ligand has an aromatic A ring (estradiol). Drugs acting through the inhibition or modulation of ERs are frequently used as a hormonal therapy for ER(+) breast cancer. Previous results suggested that the aromatic D ring was a key structural motif for the observed activity; thus, this modified steroid nucleus may provide a new scaffold for the design of novel antiestrogens. Using molecular dynamics (MD) simulation we have modeled the binding mode of the natural salpichrolide A and a synthetic analog with an aromatic D ring within the ERα. These results taken together with the calculated energetic contributions associated to the different ligand‐binding modes are consistent with a preferred inverted orientation of the steroids in the ligand‐binding pocket with the aromatic ring D occupying a position similar to that observed for the A ring of estradiol. Major changes in both dynamical behavior and global positioning of H11 caused by the loss of the ligand–His524 interaction might explain, at least in part, the molecular basis of the antagonism exhibited by these compounds. Using steered MD we also found a putative unbinding pathway for the steroidal ligands through a cavity formed by residues in H3, H7, and H11, which requires only minor changes in the overall receptor conformation. Proteins 2015; 83:1297–1306. © 2015 Wiley Periodicals, Inc.     

Model of three-dimensional structure of vitamin D receptor and its binding mechanism with 1alpha,25-dihydroxyvitamin D(3).

Comparative modeling of the vitamin D receptor three-dimensional structure and computational docking of 1alpha,25-dihydroxyvitamin D(3) into the putative binding pocket of the two deletion mutant receptors: (207-423) and (120-422, Delta [164-207]) are reported and evaluated in the context of extensive mutagenic analysis and crystal structure of holo hVDR deletion protein published recently. The obtained molecular model agrees well with the experimentally determined structure. Six different conformers of 1alpha,25-dihydroxyvitamin D(3) were used to study flexible docking to the receptor. On the basis of values of conformational energy of various complexes and their consistency with functional activity, it appears that 1alpha,25-dihydroxyvitamin D(3) binds the receptor in its 6-s-trans form. The two lowest energy complexes obtained from docking the hormone into the deletion protein (207-423) differ in conformation of ring A and orientation of the ligand molecule in the VDR pocket. 1alpha,25-Dihydroxyvitamin D(3) possessing the A-ring conformation with axially oriented 1alpha-hydroxy group binds receptor with its 25-hydroxy substituent oriented toward the center of the receptor cavity, whereas ligand possessing equatorial conformation of 1alpha-hydroxy enters the pocket with A ring directed inward. The latter conformation and orientation of the ligand is consistent with the crystal structure of hVDR deletion mutant (118-425, Delta [165-215]). The lattice model of rVDR (120-422, Delta [164-207]) shows excellent agreement with the crystal structure of the hVDR mutant. The complex obtained from docking the hormone into the receptor has lower energy than complexes for which homology modeling was used. Thus, a simple model of vitamin D receptor with the first two helices deleted can be potentially useful for designing a general structure of ligand, whereas the advanced lattice model is suitable for examining binding sites in the pocket.     

Molecular dynamics simulations of the whey protein beta-lactoglobulin.

Molecular dynamics simulations have been used to model the motions and conformational behavior of the whey protein bovine beta-lactoglobulin. Simulations were performed for the protein by itself and complexed to a single retinol ligand located in a putative interior binding pocket. In the absence of the retinol ligand, the backbone loops around the opening of this interior pocket shifted inward to partially close off this cavity, similar to the shifts observed in previously reported molecular dynamics simulations of the uncomplexed form of the homologous retinol binding protein. The protein complexed with retinol does not exhibit the same conformational shifts. Conformational changes of this type could serve as a recognition signal allowing in vivo discrimination between the free and retinol complexed forms of the beta-lactoglobulin molecule. The unusual bending of the single alpha-helix observed in the simulations of retinol binding protein were not observed in the present calculations.     

VLDL-受体的配体结合结构域结构分析

极低密度脂蛋白受体(VLDL-R)的配体结合域具有8个富含半胱氨酸的配体结合重复序列(ligand-binding repeats,LBR),被认为是与配体结合的部位。该受体与含7个类似重复序列的低密度脂蛋白受体(LDL-R)的配体结合特性明显不同。为了明确VLDL-R中8个LBR在配体结合中的作用并探讨结合位点的结构,本研究采用计算机辅助蛋白质结构预测方法,在二级结构分析的基础上,通过同源建模方法预测受体N-端328个氨基酸的配体结合域空间结构,结果显示该区域呈现弧形口袋样结构,其中前3个LBR结构紧凑,呈棒状,负电荷相对集中,推测这一特征结构是配体结合的重要结构基础,结构分析同时表明在LBR5与LBR6之间连接区的弹性结构可以赋予结合位点一定的伸缩性,利于其与不同配体的结合。本研究预测结果首次提出了VLDL-R配体结构域结构及结合位点的结构特征,并与已有的实验结果一致。     

2.1 A crystal structure of human PXR in complex with the St. John's wort compound hyperforin

The nuclear xenobiotic receptor PXR is activated by a wide variety of clinically used drugs and serves as a master regulator of drug metabolism and excretion gene expression in mammals. St. John's wort is used widely in Europe and the United States to treat depression. This unregulated herbal remedy leads to dangerous drug-drug interactions, however, in patients taking oral contraceptives, antivirals, or immunosuppressants. Such interactions are caused by the activation of the human PXR by hyperforin, the psychoactive agent in St. John's wort. In this study, we show that hyperforin induces the expression of numerous drug metabolism and excretion genes in primary human hepatocytes. We present the 2.1 A crystal structure of hyperforin in complex with the ligand binding domain of human PXR. Hyperforin induces conformational changes in PXR's ligand binding pocket relative to structures of human PXR elucidated previously and increases the size of the pocket by 250 A(3). We find that the mutation of individual aromatic residues within the ligand binding cavity changes PXR's response to particular ligands. Taken together, these results demonstrate that PXR employs structural flexibility to expand the chemical space it samples and that the mutation of specific residues within the ligand binding pocket of PXR tunes the receptor's response to ligands.     

A consensus-binding structure for adenine at the atomic level permits searching for the ligand site in a wide spectrum of adenine-containing complexes

Attempts to derive structural features of ligand-binding sites have traditionally involved seeking commonalities at the residue level. Recently, structural studies have turned to atomic interactions of small molecular fragments to extract common binding-site properties. Here, we explore the use of larger ligand elements to derive a consensus binding structure for the ligand as a whole. We superimposed multiple molecular structures from a nonredundant set of adenosine-5'-triphosphate (ATP) protein complexes, using the adenine moiety as template. Clustered binding-site atoms of compatible atomic classes forming attractive contacts with the adenine probe were extracted. A set of atomic clusters characterizing the adenine binding pocket was then derived. Among the clusters are three vertices representing the interactions of adenine atom N6 with its protein-binding niche. These vertices, together with atom C6 of the purine ring system, complete the set of four vertices for the pyramid-like structure of the N6 anchor atom. Also, the sequence relationship for the adenine-binding loop interacting with the C2-N6 end of the conjugated ring system is expanded to include a third hydrophilic cluster interacting with atom N1. A search procedure involving interatomic distances between cluster centers was formulated and applied to seek putative binding sites in test cases. The results show that a consensus network of clusters, based on an adenine probe and an ATP-complexed training set of proteins, is sufficient to recognize the experimental cavity for adenine in a wide spectrum of ligand-protein complexes.     

Mbt,a Drosophila PAK protein,combines with Cdc42 to regulate photoreceptor cell morphogenesis

The Drosophila gene mushroom bodies tiny (mbt) encodes a putative p21-activated kinase (PAK), a family of proteins that has been implicated in a multitude of cellular processes including regulation of the cytoskeleton, cell polarisation, control of MAPK signalling cascades and apoptosis. The mutant phenotype of mbt is characterised by fewer neurones in the brain and the eye, indicating a role of the protein in cell proliferation, differentiation or survival. We show that mutations in mbt interfere with photoreceptor cell morphogenesis. Mbt specifically localises at adherens junctions of the developing photoreceptor cells. A structure-function analysis of the Mbt protein in vitro and in vivo revealed that the Mbt kinase domain and the GTPase binding domain, which specifically interacts with GTP-loaded Cdc42, are important for Mbt function. Besides regulation of kinase activity, another important function of Cdc42 is to recruit Mbt to adherens junctions. We propose a role for Mbt as a downstream effector of Cdc42 in photoreceptor cell morphogenesis.     

Sampling and energy evaluation challenges in ligand binding protein design

The steroid hormone 17α‐hydroxylprogesterone (17‐OHP) is a biomarker for congenital adrenal hyperplasia and hence there is considerable interest in development of sensors for this compound. We used computational protein design to generate protein models with binding sites for 17‐OHP containing an extended, nonpolar, shape‐complementary binding pocket for the four‐ring core of the compound, and hydrogen bonding residues at the base of the pocket to interact with carbonyl and hydroxyl groups at the more polar end of the ligand. Eight of 16 designed proteins experimentally tested bind 17‐OHP with micromolar affinity. A co‐crystal structure of one of the designs revealed that 17‐OHP is rotated 180° around a pseudo‐two‐fold axis in the compound and displays multiple binding modes within the pocket, while still interacting with all of the designed residues in the engineered site. Subsequent rounds of mutagenesis and binding selection improved the ligand affinity to nanomolar range, while appearing to constrain the ligand to a single bound conformation that maintains the same “flipped” orientation relative to the original design. We trace the discrepancy in the design calculations to two sources: first, a failure to model subtle backbone changes which alter the distribution of sidechain rotameric states and second, an underestimation of the energetic cost of desolvating the carbonyl and hydroxyl groups of the ligand. The difference between design model and crystal structure thus arises from both sampling limitations and energy function inaccuracies that are exacerbated by the near two‐fold symmetry of the molecule.     

Solution structure of ileal lipid binding protein in complex with glycocholate.

Ileal lipid binding protein (ILBP) is a cytosolic lipid-binding protein that binds both bile acids and fatty acids. We have determined the solution structure of porcine ILBP in complex with glycocholate by homonuclear and heteronuclear two-dimensional NMR spectroscopy. The conformation of the protein-ligand complex was determined by restrained energy minimization and simulated annealing calculations after docking the glycocholate ligand into the protein structure. The overall tertiary structure of ILBP is highly analogous to the three-dimensional structures of several other intracellular lipid binding proteins (LBPs). Like the apo-structure, the bile-acid complex of ILBP is composed of 10 anti-parallel beta-strands that form a water-filled clam-shell structure, and two short alpha-helices. Chemical shift data indicated that the bile acid ligand is bound inside the protein cavity. Furthermore, 13C-edited heteronuclear single-quantum correlation-NOESY experiments showed NOE contacts between several aromatic residues located in the proposed bile acid portal region and the 13C-labeled ligand. A single bile acid molecule is bound inside the protein, with the steroid moiety penetrating deep into the water-accessible internal cavity, such that ring A is located right above the plane of the Trp49 indole ring. The carboxylate tail of the ligand is protruding from the proposed bile acid portal into the surrounding aqueous solution. The body of the steroid moiety is oriented with the nonpolar face in contact with the mostly hydrophobic residues of beta-strands C, D and E, while the polar face shows contacts with the side-chains of Tyr97, His99, Glu110 and Arg121 in beta-strands H, I and J. Thus, the conformational arrangement of the ligand complex suggests that the binding affinity of ILBP for bile acid molecules is based mainly on strong hydrophobic interactions inside the protein cavity. Furthermore, this binding mode explains how ILBP can transport unconjugated and conjugated bile acids.     

Anatomy of protein pockets and cavities: measurement of binding site geometry and implications for ligand design.

Identification and size characterization of surface pockets and occluded cavities are initial steps in protein structure-based ligand design. A new program, CAST, for automatically locating and measuring protein pockets and cavities, is based on precise computational geometry methods, including alpha shape and discrete flow theory. CAST identifies and measures pockets and pocket mouth openings, as well as cavities. The program specifies the atoms lining pockets, pocket openings, and buried cavities; the volume and area of pockets and cavities; and the area and circumference of mouth openings. CAST analysis of over 100 proteins has been carried out; proteins examined include a set of 51 monomeric enzyme-ligand structures, several elastase-inhibitor complexes, the FK506 binding protein, 30 HIV-1 protease-inhibitor complexes, and a number of small and large protein inhibitors. Medium-sized globular proteins typically have 10-20 pockets/cavities. Most often, binding sites are pockets with 1-2 mouth openings; much less frequently they are cavities. Ligand binding pockets vary widely in size, most within the range 10(2)-10(3)A3. Statistical analysis reveals that the number of pockets and cavities is correlated with protein size, but there is no correlation between the size of the protein and the size of binding sites. Most frequently, the largest pocket/cavity is the active site, but there are a number of instructive exceptions. Ligand volume and binding site volume are somewhat correlated when binding site volume is < or =700 A3, but the ligand seldom occupies the entire site. Auxiliary pockets near the active site have been suggested as additional binding surface for designed ligands (Mattos C et al., 1994, Nat Struct Biol 1:55-58). Analysis of elastase-inhibitor complexes suggests that CAST can identify ancillary pockets suitable for recruitment in ligand design strategies. Analysis of the FK506 binding protein, and of compounds developed in SAR by NMR (Shuker SB et al., 1996, Science 274:1531-1534), indicates that CAST pocket computation may provide a priori identification of target proteins for linked-fragment design. CAST analysis of 30 HIV-1 protease-inhibitor complexes shows that the flexible active site pocket can vary over a range of 853-1,566 A3, and that there are two pockets near or adjoining the active site that may be recruited for ligand design.     

NMR studies of 4-19F-phenylalanine-labeled intestinal fatty acid binding protein: evidence for conformational heterogeneity in the native state

(19)F-Nuclear magnetic resonance (NMR) studies have been carried out after incorporation of 4-(19)F-phenylalanine into the intestinal fatty acid binding protein (IFABP), a protein composed of two beta-sheets containing a large hydrophobic cavity into which ligands bind. NMR spectra have been obtained with both the ligand-free and ligand-bound (oleate) forms. There are 29 residues involved in van der Waals or hydrophobic interactions or both to form a U-shaped ligand binding pocket (Sacchettni J. C., Scapin G., Gopaul D., and Gordon J. I. (1992) J. Biol. Chem. 267, 23534-23545). The protein contains eight phenylalanines, and all are included in those residues that line the pocket. Peak assignments were made using site-specific incorporation of 4-(19)F-phenylalanine. Fluorine is a highly sensitive probe to monitor the conformation and dynamics of the side chains in native state. We find that chemical exchange in the binding pocket exists in the native apo- and holo-state. Of the eight phenylalanine residues, Phe2, Phe47, Phe62, Phe68, and Phe93 are arranged on one side of the binding pocket, and all exist in two conformations with Phe2, Phe47, and Phe62 showing exchange cross-peaks with minor conformation in (19)F-(19)F nuclear Overhauser effect (NOESY) spectra. The line widths of Phe68 and Phe93 are broader than those of other phenylalanine residues and can be deconvoluted into two peaks. Phe47, Phe62, Phe68, Phe93, and Trp82 have been proposed to be involved in the early stage of collapse (Ropson, I. J., and Frieden, C. (1992) Proc. Natl. Acad. Sci U.S.A. 89, 7222-7226), but a temperature study suggests that Phe47 behaves differently than other residues and may be more involved in a later stage of folding, for example, side chain stabilization. In the holo-form, Phe17 shows an extra exchange cross-peak in addition to those exchange cross-peaks observed in apo-form. Holo-IFABP exhibits broader line width than the apo-form, suggesting more flexibility of the binding cavity upon ligand binding.     

Mutational analysis of the ligand binding domain of the low density lipoprotein receptor

The ligand binding domain of the low density lipoprotein (LDL) receptor contains seven imperfect repeats of a 40-amino acid cysteine-rich sequence. Each repeat contains clustered negative charges that have been postulated as ligand-binding sites. The adjacent region of the protein, the growth factor homology region, contains three cysteine-rich repeats (A-C) whose sequence differs from those in the ligand binding domain. To dissect the contribution of these different cysteine-rich repeats to ligand binding, we used oligonucleotide-directed mutagenesis to alter expressible cDNAs for the human LDL receptor which were then introduced into monkey COS cells by transfection. We measured the ability of the mutant receptors to bind LDL, which contains a single protein ligand for the receptor (apoB-100), and beta-migrating very low density lipoprotein (beta-VLDL), which contains apoB-100 plus multiple copies of another ligand (apoE). The results show that repeat 1 is not required for binding of either ligand. Repeats 2 plus 3 and repeats 6 plus 7 are required for maximal binding of LDL, but not beta-VLDL. Repeat 5 is required for binding of both ligands. Repeat A in the growth factor homology region is required for binding of LDL, but not beta-VLDL. Repeat B is not required for ligand binding. These results support a model for the LDL receptor in which various repeats play additive roles in ligand binding, each repeat making a separate contribution to the binding event.