Characterization of the estradiol-binding site structure of human pancreas-specific protein disulfide isomerase: indispensable role of the hydrogen bond between His278 and the estradiol 3-hydroxyl group

Estradiol (E(2)), a female sex hormone, has important biological functions. Human pancreas-specific protein disulfide isomerase (PDIp), a protein folding catalyst, was recently found to be able to bind E(2). Here we report the characterization of its E(2)-binding site by using biochemical methods coupled with molecular modeling tools. Analysis of various truncated PDIp proteins showed that the b-b' fragment contains an intact E(2)-binding site that has the same binding affinity as the full-length PDIp protein, with apparent K(d) values of approximately 170 nM. Computational modeling and docking analyses revealed that the E(2)-binding site in the b-b' fragment is located in a hydrophobic pocket composed mainly of the b' domain and partially of the b domain. The hydrogen bond, formed between the 3-hydroxyl group of E(2) (donor) and PDIp's His278 (acceptor), is indispensable for its binding. By contrast, the 17β-hydroxyl group of E(2) is of negligible importance for E(2) binding. This binding model was jointly confirmed by a series of experiments, such as selective mutation of the binding site amino acid residues and selective modification of the ligand structures.     

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.     

Molecular modeling of cardiac glycoside binding by the human sequence monoclonal antibody 1B3

The amino acid sequences of the heavy- and light-chain variable regions of the high-affinity human sequence antidigoxin monoclonal antibody 1B3 (mAb 1B3) were determined, and a structural model for the mAb's variable region was developed by homology modeling techniques. The structural model provided the basis for computationally docking digoxin and eight related cardiac glycosides into the putative binding site of mAb 1B3. Analysis of the consensus binding mode obtained for digoxin showed that the cardenolide moiety of digoxin is deeply embedded in a predominantly hydrophobic, narrow cavity, whereas the terminal, gamma-carbohydrate group is solvent-exposed. The docking results indicated that the primary driving forces for digoxin binding by mAb 1B3 are hydrophobic interactions with the digoxin steroid ring system and hydrogen bonds with the digitoxose groups. The binding model accounts for the experimentally observed variations in mAb 1B3 binding affinity for various structural analogs of digoxin used previously to develop a 3D structure-activity relationship model of drug binding (Farr CD, Tabet MR, Ball WJ Jr, Fishwild DM, Wang X, Nair AC, Welsh WJ. Three-dimensional quantitative structure-activity relationship analysis of ligand binding to human sequence antidigoxin monoclonal antibodies using comparative molecular field analysis. J Med Chem 2002;45:3257-3270). In particular, the hydrogen bond pattern is consistent with the unique sensitivity of mAb 1B3's binding affinity to the number of sugar residues present in a cardiac glycoside. The hydrophobic environment about the steroid moiety of digoxin is compatible with the mAb's reduced affinity for ligands that possess hydrophilic hydroxyl and acetyl group modifications in this region. The model also indicated that most of the amino acid residues in contact with the ligand reside in or about the three complementarity determining regions (CDRs) of the heavy chain and the third CDR of the light chain. A comparison of the 1B3 binding model with the crystal structures of two murine antidigoxin mAbs revealed similar binding patterns used by the three mAbs, such as a high frequency of occurrence of aromatic, hydrophobic residues in the CDRs and a dominant role of the heavy chain CDR3 in antigen binding.     

Molecular mechanisms of estrogen recognition and 17-keto reduction by human 17beta-hydroxysteroid dehydrogenase 1

The reduction of inactive estrone (E1) to the active estrogen 17beta-estradiol (E2) is catalyzed by type 1 17beta-hydroxysteroid dehydrogenase (17HSD1). Crystallographic studies, modeling and activity measurement of mutants and chimeric enzymes have led to the understanding of its mechanism of action and the molecular basis for the estrogenic specificity. An electrophilic attack on the C17-keto oxygen by the Tyr 155 hydroxyl is proposed for initiation of the transition state. The active site is a hydrophobic pocket with catalytic residues at one end and the recognition machinery on the other. Tyr 155, Lys 159 and Ser 142 are essential for the activity. The presence of certain other amino acids near the substrate recognition end of the active site including His 152 and Pro 187 is critical to the shape complementarity of estrogenic ligands. His 221 and Glu 282 form hydrogen bonds with 3-hydroxyl of the aromatic A-ring of the ligand. This mechanism of recognition of E1 by 17HSD1 is similar to that of E2 by estrogen receptor alpha. In a ternary complex with NADP(+) and equilin, an equine estrogen with C7=C8 double bond, the orientation of C17=O of equilin relative to the C4-hydride is more acute than the near normal approach of the hydride for the substrate. In the apo-enzyme structure, a substrate-entry loop (residues 186-201) is in an open conformation. The loop is closed in this complex and Phe 192 and Met 193 make contacts with the ligand. Residues of the entry loop could be partially responsible for the estrogenic specificity.     

Equilenin: a specific fluorescent probe for steroid—protein interactions in sex steroid-binding protein

Equilenin, a naturally fluorescent steroid, has high binding affinity for human sex steroid-binding protein (SBP). At 4°C the equilibrium association constant is ～6 × 10⁷ M^?1. The fluorescence excitation and emission spectra of the steroid—protein complex indicate that both hydrophobic interactions and hydrogen bonding of the 3'-hydroxyl group of the estrogen are important in its binding to the protein. Equilenin has a substantially different 3-dimensional spatial configuration compared with the normally bound androgens, and yet exhibits very tight binding to SBP. This suggests that SBP undergoes a conformational change to accomodate equilenin.     

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.     

Molecular structure of the rat vitamin D receptor ligand binding domain complexed with 2-carbon-substituted vitamin D3 hormone analogues and a LXXLL-containing coactivator peptide

We have determined the crystal structures of the ligand binding domain (LBD) of the rat vitamin D receptor in ternary complexes with a synthetic LXXLL-containing peptide and the following four ligands: 1alpha,25-dihydroxyvitamin D(3); 2-methylene-19-nor-(20S)-1alpha,25-dihydroxyvitamin D(3) (2MD); 1alpha-hydroxy-2-methylene-19-nor-(20S)-bishomopregnacalciferol (2MbisP), and 2alpha-methyl-19-nor-1alpha,25-dihydroxyvitamin D(3) (2AM20R). The conformation of the LBD is identical in each complex. Binding of the 2-carbon-modified analogues does not change the positions of the amino acids in the ligand binding site and has no effect on the interactions in the coactivator binding pocket. The CD ring of the superpotent analogue, 2MD, is tilted within the binding site relative to the other ligands in this study and to (20S)-1alpha,25-dihydroxyvitamin D(3) [Tocchini-Valentini et al. (2001) Proc. Natl. Acad. Sci. U.S.A. 98, 5491-5496]. The aliphatic side chain of 2MD follows a different path within the binding site; nevertheless, the 25-hydroxyl group at the end of the chain occupies the same position as that of the natural ligand, and the hydrogen bonds with histidines 301 and 393 are maintained. 2MbisP binds to the receptor despite the absence of the 25-hydroxyl group. A water molecule is observed between His 301 and His 393 in this structure, and it preserves the orientation of the histidines in the binding site. Although the alpha-chair conformer is highly favored in solution for the A ring of 2AM20R, the crystal structures demonstrate that this ring assumes the beta-chair conformation in all cases, and the 1alpha-hydroxyl group is equatorial. The peptide folds as a helix and is anchored through hydrogen bonds to a surface groove formed by helices 3, 4, and 12. Electrostatic and hydrophobic interactions between the peptide and the LBD stabilize the active receptor conformation. This stablization appears necessary for crystal growth.     

Differential recognition of resveratrol isomers by the human estrogen receptor-alpha: molecular dynamics evidence for stereoselective ligand binding

Resveratrol (RSVL) is a phytoestrogen that occurs naturally in two forms (trans- (E) and cis- (Z)). We have conducted molecular dynamics (MD) studies to differentially characterize the estrogen receptor-alpha (ER-alpha) binding profiles of RSVL stereoisomers. Favorable orientations for RSVL isomers at the ER-alpha pocket were first inferred from (1) alignment with pharmacophoric elements of the pure ER-alpha agonists estradiol (E2) and (2) assessment of ligand recognition by the ER-alpha binding domain. Subsequently, these orientations for RSVL isomers were subjected to MD analyses versus E2. A 100-picosecond MD simulation revealed that E2 contributed four stable hydrogen bonds with the key ER-alpha pocket residue: Arg394, Glu353, His524, and Leu525. Further, E2 displayed favorable binding energy, conformational energy change (DeltaE), and movement of the binding pocket residues (RMSd). Compared to E2, (E)-RSVL lacked a hydrogen bond (HB) with His524 but formed three additional bonds with Gly521, Phe404, and Met343 of the ER-alpha pocket. Further, (E)-RSVL conferred more favorable energy of interaction, less favorable DeltaE, but comparable RMSd values. In contrast, (Z)-RSVL orientations missed hydrogen bonding (HB) with His524 and Leu525, two essential ligand binding residues, and/or produced considerably less favorable-binding energy, -DeltaE, and -RMSd values than did (E)-RSVL. In conclusion, the present study demonstrates the utility of this MD model in distinguishing between RSVL stereoisomers. The weak binding of (Z)-RSVL by the human ER-alpha binding is congruent with its inferior ligand profiles in ER-endowed biological systems. Further, evidence is provided for a considerable variation in the mode of recognition of the mixed agonist/antagonist (E)-RSVL, and the pure agonist E2.     

A nonisotopic estrogen receptor-based assay to detect estrogenic compounds

We have used the ligand binding domain of the recombinant human estrogen receptor (hER) to develop a nonisotopic assay for detection of estrogenic compounds. The assay is based on competition of the estrogenic ligand with 17beta-estradiol for binding to the receptor, which leaves 17beta-estradiol free to bind to an anti-17beta-estradiol antibody. Unbound anti-17beta-estradiol antibody then binds to immobilized 17beta-estradiol-protein conjugate (to which hER is unable to bind for steric reasons), and is detected by an enzyme-labeled anti-rabbit IgG antibody. We used the assay to detect estrogenic compounds (mainly members of the flavonoid group of plant polyphenols) in a variety of commonly consumed plant foods.     

Non-additivity of Functional Group Contributions in Protein-Ligand Binding: A Comprehensive Study by Crystallography and Isothermal Titration Calorimetry

Additivity of functional group contributions to protein-ligand binding is a very popular concept in medicinal chemistry as the basis of rational design and optimized lead structures. Most of the currently applied scoring functions for docking build on such additivity models. Even though the limitation of this concept is well known, case studies examining in detail why additivity fails at the molecular level are still very scarce. The present study shows, by use of crystal structure analysis and isothermal titration calorimetry for a congeneric series of thrombin inhibitors, that extensive cooperative effects between hydrophobic contacts and hydrogen bond formation are intimately coupled via dynamic properties of the formed complexes. The formation of optimal lipophilic contacts with the surface of the thrombin S3 pocket and the full desolvation of this pocket can conflict with the formation of an optimal hydrogen bond between ligand and protein. The mutual contributions of the competing interactions depend on the size of the ligand hydrophobic substituent and influence the residual mobility of ligand portions at the binding site. Analysis of the individual crystal structures and factorizing the free energy into enthalpy and entropy demonstrates that binding affinity of the ligands results from a mixture of enthalpic contributions from hydrogen bonding and hydrophobic contacts, and entropic considerations involving an increasing loss of residual mobility of the bound ligands. This complex picture of mutually competing and partially compensating enthalpic and entropic effects determines the non-additivity of free energy contributions to ligand binding at the molecular level.     

Molecular dynamics simulations of the Apo-, Holo-, and acyl-forms of Escherichia coli acyl carrier protein

Acyl carrier protein (ACP) is an essential co-factor protein in fatty acid biosynthesis that shuttles covalently bound fatty acyl intermediates in its hydrophobic pocket to various enzyme partners. To characterize acyl chain-ACP interactions and their influence on enzyme interactions, we performed 19 molecular dynamics (MD) simulations of Escherichia coli apo-, holo-, and acyl-ACPs. The simulations were started with the acyl chain in either a solvent-exposed or a buried conformation. All four short-chain (< or = C10) and one long-chain (C16) unbiased acyl-ACP MD simulation show the transition of the solvent-exposed acyl chain into the hydrophobic pocket of ACP, revealing its pathway of acyl chain binding. Although the acyl chain resides inside the pocket, Thr-39 and Glu-60 at the entrance stabilize the phosphopantetheine linker through hydrogen bonding. Comparisons of the different ACP forms indicate that the loop region between helices II and III and the prosthetic linker may aid in substrate recognition by enzymes of fatty acid synthase systems. The MD simulations consistently show that the hydrophobic binding pocket of ACP is best suited to accommodate an octanoyl group and is capable of adjusting in size to accommodate chain lengths as long as decanoic acid. The simulations also reveal a second, novel binding mode of the acyl chains inside the hydrophobic binding pocket directed toward helix I. This study provides a detailed dynamic picture of acyl-ACPs that is in excellent agreement with available experimental data and, thereby, provides a new understanding of enzyme-ACP interactions.     

An approach to rational ligand-design based on a thermodynamic analysis

Thermodynamic analysis is an effective tool in screening of lead-compounds for development of potential drug candidates. In most cases, a ligand achieve high affinity and specificity to a target protein by means of both favorable enthalpy and entropy terms, which can be reflected in binding profiles of Isothermal Titration Calorimetry (ITC). A favorable enthalpy change suggests the contribution of noncovalent contacts such as hydrogen bonding and van der Waals interaction between a ligand and its target protein. In general, optimization of binding enthalpy is more difficult than that of entropies in ligand-design; therefore, it is desirable to choose firstly a lead-compound based on its binding enthalpic gain. In this paper, we demonstrate the utility of thermodynamic approach to ligand screening using anti-ciguatoxin antibody 10C9 as a model of a target protein which possesses a large hydrophobic pocket. As a result of this screening, we have identified three compounds that could bind to the antigen-binding pocket of 10C9 with a few kcal/mol of favorable binding enthalpy. Comparison of their structure with the proper antigen ciguatoxin CTX3C revealed that 10C9 rigorously identifies their cyclic structure and a characteristic hydroxyl group. ITC measurement might be useful and powerful for a rational ligand screening and the optimization of the ligand; the enthalpic gain is an effective index for ligand-design studies.     

Interaction of aminoadamantane derivatives with the influenza A virus M2 channel-Docking using a pore blocking model

Interaction of aminoadamantanes with the influenza A virus M2 proton channel was analyzed by docking simulations of a series of synthetic aminoadamantane derivatives, of differing binding affinity, into the crystal structure of the transmembrane (M2TM) pore. The pore blocking model tested in the ‘gas phase’ describes qualitatively the changes on the relative binding affinities of the compounds (although a series of highly hydrophobic ligands which seem to have little capacity for different specific interactions with their receptor). The docking calculations predicted poses in which the adamantane ring is surrounded mainly by the alkyl side chains of Val27 or Ala30 and the ligand’s amino group is generally hydrogen bonded with hydroxyls of Ser31 or carbonyls of Val27 or carbonyls of Ala30, the former (Ser31) being the most stable and most frequently observed. The binding of the ligand is a compromise between hydrogen bonding ability, which is elevated by a primary amino group, and apolar interactions, which are increased by the ability of the lipophilic moiety to adequately fill a hydrophobic pocket within the M2TM pore. A delicate balance of these hydrophobic contributions is required for optimal interaction.     

The effects of E7 and E11 mutations on the kinetics of ligand binding to R state human hemoglobin

Association and dissociation rate constants for O2, CO, and methyl isocyanide binding to native and distal pocket mutants of R state human hemoglobin were measured using ligand displacement and partial photolysis techniques. Individual rate constants for the alpha and beta subunits were resolved by comparisons between the kinetic behavior of the native and mutant proteins. His-E7 was replaced with Gly and Gln in both alpha and beta subunits and with Phe in beta subunits alone. In separate experiments Val-E11 was replaced with Ala, Leu, and Ile in each globin chain. The parameters describing ligand binding to R state alpha subunits are sensitive to the size and polarity of the amino acids at positions E7 and E11. The distal histidine in this subunit inhibits the bimolecular rate of binding of both O2 and CO, sterically hinders bound CO and methyl isocyanide, and stabilizes bound O2 by hydrogen bonding. The Val-E11 side chain in alpha chains also appears to be part of the kinetic barrier to O2 and CO binding since substitution with Ala causes approximately 10-fold increases in the association rate constants for the binding of these diatomic ligands. However, substitution of Val-E11 by Ile produces only small decreases in the rates of ligand binding to alpha subunits. For R state beta subunits, the bimolecular rates of O2 and CO binding are intrinsically large, approximately 2-5-fold greater than those for alpha subunits, and with the exception of Val-E11----Ile mutation, little affected by substitutions at either the E7 or E11 positions. For the beta Val-E11----Ile mutant the association rate and equilibrium constants for all three ligands decreased 10-50-fold. All of these results agree with Shaanan's conclusions that the distal pocket in liganded beta subunits is more open whereas in alpha subunits bound ligands are more sterically hindered by adjacent distal residues (Shaanan, B. (1983) J. Mol. Biol. 171, 31-59). In the case of O2 binding to alpha subunits, the unfavorable steric effects are compensated by the formation of a hydrogen bond between the nitrogen atom of His-E7 and bound dioxygen.     

The sex steroid binding protein (SBP or SHBG) of human plasma: identification of Tyr-57 and Met-107 in the steroid binding site

Tyrosine-57 (Y57) and methionine-107 (M107) have been identified in the binding site of the sex steroid binding protein (SBP) (or sex hormone binding globulin) of human plasma by replacing the two amino acids with a number of residues of varying structure. Replacement of Y57 with phenylalanine resulted in a fourfold increase in the K(d) of 5 alpha-dihydrotestosterone but left the K(d) of 17 beta-estradiol unchanged. Except in two cases, no further loss in binding took place when replacing Y57 with other residues, suggesting that the phenolic group of Y57 may form a hydrogen bond with the ligand. Replacement of M107 with isoleucine increased the 5 alpha-dihydrotestosterone K(d) fourfold to a value equal to that of rabbit SBP, which contains isoleucine at the corresponding position; however, the K(d) of 17 beta-estradiol remained unchanged. Replacement of M107 with threonine resulted in a tenfold decrease in 5 alpha-dihydrotestosterone binding affinity, whereas replacement with leucine left the K(d) unchanged. These data indicate that substitutions on the beta-carbon of the amino acid side-chain at position 107 causes significant loss of binding affinity but, as in the case of Y57, the activity was not totally eliminated. We conclude that Y57 and M107 form part of a structural motif within the steroid binding site and specifically contribute binding energy to ring A of 5 alpha-dihydrotestosterone but not to ring A of 17 beta-estradiol. We also propose that the integrated contribution of several side chains may be required to optimize the ligand affinity of the steroid binding site. This proposal may fit a 'lock and key' model where little movement of the side chains occurs during binding as might be expected for a rigid structure like the steroid nucleus.     

Structural requirements of para-alkylphenols to bind to estrogen receptor.

Octyl- and nonylphenols in the environment have been proposed to function as estrogens. To gain insight into their structural essentials in binding to the estrogen receptor, a series of phenols with saturated alkyl groups at the para position, HO-C6H4-CnH2n+1 (n = 0-12), were examined for their ability to displace [3H]17beta-estradiol in the recombinant human estrogen receptor, which was expressed in Sf9 cells using the vaculovirus expression system. All tested para-alkylphenols were found to bind fully to the estrogen receptors in a dose-dependent manner. The interaction of alkylphenols with the receptor became stronger with increase in the number of the alkyl carbons and the activity was maximized with n = 9 of nonylphenol. Phenol (n = 0) exhibited weak but full binding to the receptor, whereas anisole with a protected phenolic hydroxyl group was completely inactive. Also, alkanes such as n-octane, 2,2, 4-trimethylpentane corresponding to tert-octane, and n-nonane exhibited no binding. The results indicate that the binding of para-alkylphenols to the estrogen receptor is due to the effect of covalent bonding of two constituents of the phenol and alkyl groups, which correspond to the A-ring and hydrophobic moiety of the steroid structure, respectively. When alkylphenols were examined for their receptor binding conformation by 1H-NMR measurements and ab initio molecular orbital calculations, it was suggested that nonbranched alkyl groups are in an extended conformation, while branched alkyl groups are in a folded conformation. These results suggest that branched and nonbranched alkyl moieties of alkylphenols interact differently with the lipophilic ligand binding cavity of the estrogen receptor when compared to the binding of 17beta-estradiol.     

NADP-Dependent enzymes. I: Conserved stereochemistry of cofactor binding

The ubiquitous redox cofactors nicotinamide adenine dinucleotides [NAD and NADP] are very similar molecules, despite their participation in substantially different biochemical processes. NADP differs from NAD in only the presence of an additional phosphate group esterified to the 2′-hydroxyl group of the ribose at the adenine end and yet NADP is confined with few exceptions to the reactions of reductive biosynthesis, whereas NAD is used almost exclusively in oxidative degradations. The discrimination between NAD and NADP is therefore an impressive example of the power of molecular recognition by proteins. The many known tertiary structures of NADP complexes affords the possibility for an analysis of their discrimination. A systematic analysis of several crystal structures of NAD(P)-protein complexes show that: 1) the NADP coenzymes are more flexible in conformation than those of NAD; 2) although the protein-cofactor interactions are largely conserved in the NAD complexes, they are quite variable in those of NADP; and 3) in both cases the pocket around the nicotinamide moiety is substrate dependent. The conserved and variable interactions between protein and cofactors in the respective binding pockets are reported in detail. Discrimination between NAD and NADP is essentially a consequence of the overall pocket and not of a few residues. A clear fingerprint in NAD complexes is a carboxylate side chain that chelates the diol group at the ribose near the adenine, whereas in NADP complexes an arginine side chain faces the adenine plane and interacts with the phosphomonoester. The latter type of interaction might be a general feature of recognition of nucleotides by proteins. Other features such as strand-like hydrogen bonding between the NADP diphosphate moeties and the protein are also significant. The NADP binding pocket properties should prove useful in protein engineering and design. © 1997 Wiley-Liss Inc.     

Electrostatic interactions of androgens and progesterone derivatives with rainbow trout estrogen receptor

In primary cultures of immature male rainbow trout (rt) hepatocytes, vitellogenin (Vg) gene expression is regulated by E(2) via the estrogen receptor (ER). However, steroids other than estrogens can also stimulate Vg gene expression. These steroids are hardly converted into E(2) during incubation and their stimulatory activity is completely inhibited by tamoxifen implying rtER involvement. These steroids have no or a slightly positive charge on the Connolly surface. In contrast, steroids that failed to stimulate Vg gene expression had a strong positive or negative charge around rings C and D due to polarization. The amino acid sequences of the ligand binding domains (LBD) of rtER and human ER alpha have 57.7% homology; only one amino acid differs in the presumed steroid binding site. We modeled the three-dimensional structure of the LBD of rtER using X-ray crystallographic data for hER alpha in order to investigate the fit (structural and electrostatic) between steroid and rtER. Two factors are essential for binding to rtER: (i) hydroxyl or carbonyl groups near C3 and C17 of the steroids (hydrophilic regions) that can form hydrogen bonds with His(489), Arg(359), and Glu(318), (ii) a hydrophobic steroid nucleus that interacts with a hydrophobic region of the rtER LBD through van der Waals forces. If polar functional groups are present, the hydrophobic interaction between steroid and the rtER LBD is considerably weakened.     

Optimized Hydrophobic Interactions and Hydrogen Bonding at the Target-Ligand Interface Leads the Pathways of Drug-Designing

Background

Weak intermolecular interactions such as hydrogen bonding and hydrophobic interactions are key players in stabilizing energetically-favored ligands, in an open conformational environment of protein structures. However, it is still poorly understood how the binding parameters associated with these interactions facilitate a drug-lead to recognize a specific target and improve drugs efficacy. To understand this, comprehensive analysis of hydrophobic interactions, hydrogen bonding and binding affinity have been analyzed at the interface of c-Src and c-Abl kinases and 4-amino substituted 1H-pyrazolo [3, 4-d] pyrimidine compounds.

Methodology

In-silico docking studies were performed, using Discovery Studio software modules LigandFit, CDOCKER and ZDOCK, to investigate the role of ligand binding affinity at the hydrophobic pocket of c-Src and c-Abl kinase. Hydrophobic and hydrogen bonding interactions of docked molecules were compared using LigPlot program. Furthermore, 3D-QSAR and MFA calculations were scrutinized to quantify the role of weak interactions in binding affinity and drug efficacy.

Conclusions

The in-silico method has enabled us to reveal that a multi-targeted small molecule binds with low affinity to its respective targets. But its binding affinity can be altered by integrating the conformationally favored functional groups at the active site of the ligand-target interface. Docking studies of 4-amino-substituted molecules at the bioactive cascade of the c-Src and c-Abl have concluded that 3D structural folding at the protein-ligand groove is also a hallmark for molecular recognition of multi-targeted compounds and for predicting their biological activity. The results presented here demonstrate that hydrogen bonding and optimized hydrophobic interactions both stabilize the ligands at the target site, and help alter binding affinity and drug efficacy.