Calculation of binding affinities of HIV-1 RT and β-secretase inhibitors using the linear interaction energy method with explicit and continuum solvation approaches

The linear interaction energy (LIE) approach has been applied to estimate the binding free energies of representative sets of HIV-1 RT and β-Secretase inhibitors, using both molecular dynamics (MD) and tethered energy minimization sampling protocols with the OPLS-AA potential, using a range of solvation methodologies. Generalized Born (GB), ‘shell’ and periodic boundary condition (PBC) solvation were used, the latter with reaction field (RF) electrostatics. Poisson-Boltzmann (PB) and GB continuum electrostatics schemes were applied to the simulation trajectories for each solvation type to estimate the electrostatic ligand-water interaction energy in both the free and bound states. Reasonable agreement of the LIE predictions was obtained with respect to experimental binding free energy estimates for both systems: for instance, ‘PB’ fits on MD trajectories carried out with PBC solvation and RF electrostatics led to models with standard errors of 1.11 and 1.03 kcal mol⁻¹ and coefficients of determination, r ² of 0.76 and 0.75 for the HIV-1 RT and β-Secretase sets. However, it was also found that results from MD sampling using PBC solvation provided only slightly better fits than from simulations using shell or Born solvation or tethered energy minimization sampling. Figure Evolution of the running averages for compound H11 (binding to HIV-1RT) of the bound state ligand-water and ligand-protein interaction energies. The ligand-water electrostatic terms are twice the corresponding GB and PB electrostatic solvation free energies. The ligand-receptor van der Waals and Coulombic interaction energies are also shown, in addition to the ligand-water van der Waals interaction term. The terms were calculated (without application of a cut-off) from a trajectory sampled under PBC solvation with reaction field electrostatics Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Intracellular trafficking of the β-secretase and processing of amyloid precursor protein

The main component of the amyloid plaques found in the brains of those with Alzheimer's disease (AD) is a polymerized form of the β-amyloid peptide (Aβ) and is considered to play a central role in the pathogenesis of this neurodegenerative disorder. Aβ is derived from the proteolytic processing of the amyloid precursor protein (APP). Beta site APP-cleaving enzyme, BACE1 (also known as β-secretase) is a membrane-bound aspartyl protease responsible for the initial step in the generation of Aβ peptide and is thus a prime target for therapeutic intervention. Substantive evidence now indicates that the processing of APP by BACE1 is regulated by the intracellular sorting of the enzyme and, moreover, perturbations in these intracellular trafficking pathways have been linked to late-onset AD. In this review, we highlight the recent advances in the understanding of the regulation of the intracellular sorting of BACE1 and APP and illustrate why the trafficking of these cargos represent a key issue for understanding the membrane-mediated events associated with the generation of the neurotoxic Aβ products in AD.     

What determines the van der Waals coefficient beta in the LIE (linear interaction energy) method to estimate binding free energies using molecular dynamics simulations?

Recently a semiempirical method has been proposed by Aqvist et al. to calculate absolute and relative binding free energies. In this method, the absolute binding free energy of a ligand is estimated as deltaGbind = alpha + beta, where Vel(bound) and Vvdw(bound) are the electrostatic and van der Waals interaction energies between the ligand and the solvated protein from an molecular dynamics (MD) trajectory with ligand bound to protein and Vel(free) and Vel(free) and Vvdw(free) are the electrostatic and van der Waals interaction energies between the ligand and the water from an MD trajectory with the ligand in water. A set of values, alpha = 0.5 and beta = 0.16, was found to give results in good agreement with experimental data. Later, however, different optimal values of beta were found in studies of compounds binding to P450cam and avidin. The present work investigates how the optimal value of beta depends on the nature of binding sites for different protein-ligand interactions. By examining seven ligands interacting with five proteins, we have discovered a linear correlation between the value of beta and the weighted non-polar desolvation ratio (WNDR), with a correlation coefficient of 0.96. We have also examined the ability of this correlation to predict optimal values of beta for different ligands binding to a single protein. We studied twelve neutral compounds bound to avidin. In this case, the WNDR approach gave a better estimate of the absolute binding free energies than results obtained using the fixed value of beta found for biotin-avidin. In terms of reproducing the relative binding free energy to biotin, the fixed-beta value gave better results for compounds similar to biotin, but for compounds less similar to biotin, the WNDR approach led to better relative binding free energies.     

GPCRs through the keyhole: the role of protein flexibility in ligand binding to β-adrenoceptors

G protein-coupled receptors (GPCRs) are proteins of pharmaceutical importance, with over 30% of all drugs in clinical use targeting them. Increasing numbers of X-ray crystal (XRC) structures of GPCRs offer a wealth of data relating to ligand binding. For the β-adrenoceptors (β-ARs), XRC structures are available for human β₂- and turkey β₁-subtypes, in complexes with a range of ligands. While these structures provide insight into the origins of ligand structure-activity relationships (SARs), questions remain. The ligands in all published complexed XRC structures lack extensive substitution, with no obvious way the ligand-binding site can accommodate β₁-AR-selective antagonists with extended side-chains para- to the common aryloxypropanolamine pharmacophore. Using standard computational docking tools with such ligands generally returns poses that fail to explain known SARs. Application of our Active Site Pressurisation modelling method to β-AR XRC structures and homology models, however, reveals a dynamic area in the ligand-binding pocket that, through minor changes in amino acid side chain orientations, opens a fissure between transmembrane helices H4 and H5, exposing intra-membrane space. This fissure, which we term the “keyhole”, is ideally located to accommodate extended moieties present in many high-affinity β₁-AR-selective ligands, allowing the rest of the ligand structure to adopt a canonical pose in the orthosteric binding site. We propose the keyhole may be a feature of both β₁- and β₂-ARs, but that subtle structural differences exist between the two, contributing to subtype-selectivity. This has consequences for the rational design of future generations of subtype-selective ligands for these therapeutically important targets.     

Structure-based discovery of a novel inhibitor targeting the β-catenin/Tcf4 interaction

Overactivation or overexpression of β-catenin in the Wnt (wingless) signaling pathway plays an important role in tumorigenesis. Interaction of β-catenin with T-cell factor (Tcf) DNA binding proteins is a key step in the activation of the proliferative genes in response to upstream signals of this Wnt/β-catenin pathway. Recently, we identified a new small molecule inhibitor, named BC21 (C(32)H(36)Cl(2)Cu(2)N(2)O(2)), which effectively inhibits the binding of β-catenin with Tcf4-derived peptide and suppresses β-catenin/Tcf4 driven reporter gene activity. This inhibitor decreases the viability of β-catenin overexpressing HCT116 colon cancer cells that harbor the β-catenin mutation, and more significantly, it inhibits the clonogenic activity of these cells. Down-regulation of c-Myc and cyclin D1 expression, the two important effectors of the Wnt/β-catenin signaling, is confirmed by treating HCT116 cells with BC21. This compound represents a new and modifiable potential anticancer candidate that targets β-catenin/Tcf-4 interaction.     

Inhibitors of the p53/hdm2 protein–protein interaction—path to the clinic

A growing number of the elements identified in intracellular signaling events that affect cell growth and transformation are proteins that physically interact with each other via domains or specifically recognized amino acid sequences. Some of these intracellular protein–protein interactions are attractive targets for anticancer targeted therapy, but progress in this field has been compromised by the paucity of compounds with suitable biological profiles and pharmacological properties. This Letter covers salient achievements in the identification and development of inhibitors of the p53–hdm2 protein–protein interaction, and highlights different screening techniques and structure-based design approaches that may be brought to bear on the discovery and development of inhibitors of other therapeutically relevant intracellular protein–protein interactions.     

Contribution of the γ-secretase subunits to the formation of catalytic pore of presenilin 1 protein

γ-Secretase is an intramembrane-cleaving protease related to the etiology of Alzheimer disease. γ-Secretase is a membrane protein complex composed of presenilin (PS) and three indispensable subunits: nicastrin, Aph-1, and Pen-2. PS functions as a protease subunit forming a hydrophilic catalytic pore structure within the lipid bilayer. However, it remains unclear how other subunits are involved in the pore formation. Here, we show that the hydrophilic pore adopted with an open conformation has already been formed by PS within the immature γ-secretase complex. The binding of the subunits induces the close proximity between transmembrane domains facing the catalytic pore. We propose a model in which the γ-secretase subunits restrict the arrangement of the transmembrane domains of PS during the formation of the functional structure of the catalytic pore.     

Pharmacological evaluation of ocular β-adrenoceptors in rabbit by tissue segment binding method

AimsThis study evaluates ocular (iris, ciliary body and ciliary process) and nonocular (atria and lung) β-adrenoceptors in rabbit to characterize the plasma membrane β-adrenoceptors and binding affinities of β-adrenoceptor antagonists.Main methodsThe tissue segment binding method with a hydrophilic radioligand (?)-4-[3-t-butylamino-2-hydroxypropoxy]-[5,7-³H]benzimidazol-2-one ([³H]-CGP12177) was employed.Key findingsSpecific and saturable binding of [³H]-CGP12177 to intact tissue segments was detected by using (±)-propranolol to define nonspecific binding, showing a single population of plasma membrane binding sites with high affinity. Competition experiments with selective β₁- and β₂-adrenoceptor antagonists revealed a single population of β₂-adrenoceptors in ocular tissues and of β₁-adrenoceptors in atria, but mixed populations of β₁- and β₂-adrenoceptors in 70% and 30%, respectively, in lung. A competition curve for timolol was biphasic in lung and its binding affinity for β₂-adrenoceptors was approximately 158-fold higher than for β₁-adrenoceptors, indicating the β₂-selectivity of timolol. In contrast, competition curves for stereoisomers of befunolol, carteolol, and propranolol were monophasic in all tissues. The (?)-enantiomers of these antagonists were more potent than corresponding (+)-enantiomers in displacing from [³H]-CGP12177 binding, and the isomeric potency ratios of befunolol and carteolol were less than those of propranolol.SignificanceThis study with tissue segment binding method suggests that the binding affinity of (?)-enantiomers of β-adrenoceptor antagonists for plasma membrane β-adrenoceptors (β₁-adrenoceptors of atria, β₂-adrenoceptors of ocular tissues, and mixed β₁-/β₂-adrenoceptors of lung) is higher than that of corresponding (+)-enantiomers and their stereoselectivity is different between β-adrenoceptor antagonists.     

Conformational transition in the substrate binding domain of β-secretase exploited by NMA and its implication in inhibitor recognition: BACE1-myricetin a case study

BACE1 is a key protease involved in the proteolysis of amyloid precursor protein (APP) that generates a toxic peptide amyloid beta (Aβ), a pathological feature of Alzheimer's disease (AD). The enzyme is believed to possess an open and a closed conformation that corresponds to its free and inhibitor-bound form respectively. Here, we study the dynamic transition of BACE1 employing normal mode analysis (NMA) using a simplified elastic network model (ENM). Estimation of the catalytic cavity volume on the structures of BACE1 encoded by the lowest frequency normal mode reveals the dynamical transition of the enzyme from the open to the closed conformer. Detailed analysis reveals that concerted movement of different loop segments in the active site of the protein, namely flap regions, 10s loop, A loop and F loop, squeeze the catalytic cavity between the N-terminal and C-terminal lobe of the substrate binding domain of BACE1. We also propose that the NMA encoded multiple receptor conformations (MRC) of BACE1 elucidate the pharmacophoric feature necessary to inhibit the enzyme by a polyphenol, myricetin. van der Waals interaction is found to be the main driving force that guides the ligand induced conformational switching to the closed conformer. We suggest that NMA derived MRC of BACE1 is an efficient way to treat the receptor flexibility in docking and thus can be further applied in virtual screening and structure based drug design.     

Discovery of aminoheterocycles as a novel β-secretase inhibitor class: pH dependence on binding activity part 1

We have developed a novel series of heteroaromatic BACE-1 inhibitors. These inhibitors interact with the enzyme in a unique fashion that allows for potent binding in a non-traditional paradigm. In addition to the elucidation of their binding profile, we have discovered a pH dependent effect on the binding affinity as a result of the intrinsic pK_a of these inhibitors and the pH of the BACE-1 enzyme binding assay.     

Can one measure the free energy of binding of the histone octamer to different DNA sequences by salt-dependent reconstitution?

I explain why many recently reported measurements for the "free energy" of positioning of the histone octamer on different DNA sequences are likely to be in error: i.e. because histone octamers do not equilibrate between different DNA molecules at low salt, but only at high salt. Thus, the reported "free energies" refer to an equilibrium at high salt, under nearly dissociating conditions between protein and DNA, and they are likely to be much too small on an absolute scale. There are many other lines of evidence to suggest that the preferences of the histone octamer for different DNA sequences are rather strong and of importance in biological systems.     

Sorting nexin 17 prevents lysosomal degradation of β1 integrins by binding to the β1-integrin tail

Integrin functions are controlled by regulating their affinity for ligand, and by the efficient recycling of intact integrins through endosomes. Here we demonstrate that the Kindlin-binding site in the β1-integrin cytoplasmic domain serves as a molecular switch enabling the sequential binding of two FERM-domain-containing proteins in different cellular compartments. When β1 integrins are at the plasma membrane, Kindlins control ligand-binding affinity. However, when they are internalized, Kindlins dissociate from integrins and sorting nexin 17 (SNX17) is recruited to free β1-integrin tails in early endosomes to prevent β1-integrin degradation, leading to their recycling back to the cell surface. Our results identify SNX17 as a β1-integrin-tail-binding protein that interacts with the free Kindlin-binding site in endosomes to stabilize β1 integrins, resulting in their recycling to the cell surface where they can be reused.     

KCTD11 inhibits progression of lung cancer by binding to β-catenin to regulate the activity of the Wnt and Hippo pathways

KCTD11 has been reported to be a potential tumour suppressor in several tumour types. However, the expression of KCTD11 and its role has not been reported in human non-small cell lung cancer (NSCLC). Whether its potential molecular mechanism is related to its BTB domain is also unknown. The expression of KCTD11 in 139 NSCLC tissue samples was detected by immunohistochemistry, and its correlation with clinicopathological factors was analysed. The effect of KCTD11 on the biological behaviour of lung cancer cells was verified in vitro and in vivo. Its effect on the epithelial-mesenchymal transition(EMT)process and the Wnt/β-catenin and Hippo/YAP pathways were observed by Western blot, dual-luciferase assay, RT-qPCR, immunofluorescence and immunoprecipitation. KCTD11 is under-expressed in lung cancer tissues and cells and was negatively correlated with the degree of differentiation, tumour-node-metastasis (TNM) stage and lymph node metastasis. Low KCTD11 expression was associated with poor prognosis. KCTD11 overexpression inhibited the proliferation and migration of lung cancer cells. Further studies indicated that KCTD11 inhibited the Wnt pathway, activated the Hippo pathway and inhibited EMT processes by inhibiting the nuclear translocation of β-catenin and YAP. KCTD11 lost its stimulatory effect on the Hippo pathway after knock down of β-catenin. These findings confirm that KCTD11 inhibits β-catenin and YAP nuclear translocation as well as the malignant phenotype of lung cancer cells by interacting with β-catenin. This provides an important experimental basis for the interaction between KCTD11, β-catenin and YAP, further revealing the link between the Wnt and Hippo pathways.     

The interrelationship between ligand binding and self-association of the folate binding protein. The role of detergent–tryptophan interaction

Background

The folate binding protein (FBP) regulates homeostasis and intracellular trafficking of folic acid, a vitamin of decisive importance in cell division and growth. We analyzed whether interrelationship between ligand binding and self-association of FBP plays a significant role in the physiology of folate binding.

Methods

Self-association behavior of apo- and holo-FBP was addressed through size exclusion chromatography, SDS-PAGE, mass spectrometry, surface plasmon resonance and fluorescence spectroscopy.

Results

Especially holo-FBP exhibits concentration-dependent self-association at pH 7.4 (pI), and is more prone to associate into stable complexes than apo-FBP. Even more pronounced was the tendency to complexation between apo-FBP and holo-FBP in accord with a model predicting association between apo and holo monomers [19]. This will lead to removal of apo monomers from the reaction scheme resulting in a weak incomplete ligand binding similar to that observed at FBP concentrations < 10 nM. The presence of synthetic and natural detergents normalized folate binding kinetics and resulted in appearance of monomeric holo-FBP. Fluorescence spectroscopy indicated molecular interactions between detergent and tryptophan residues located in hydrophobic structures of apo-FBP which may participate in protein associations.

General significance

Self-association into multimers may protect binding sites, and in case of holo-FBP even folate from biological degradation. High-affinity folate binding in body secretions, typically containing 1–10 nM FBP, requires the presence of natural detergents, i.e. cholesterol and phospholipids, to avoid complexation between apo- and holo-FBP.     

Adaptor protein 2-mediated endocytosis of the β-secretase BACE1 is dispensable for amyloid precursor protein processing

The β-site amyloid precursor protein (APP)-cleaving enzyme 1 (BACE1) is a transmembrane aspartyl protease that catalyzes the proteolytic processing of APP and other plasma membrane protein precursors. BACE1 cycles between the trans-Golgi network (TGN), the plasma membrane, and endosomes by virtue of signals contained within its cytosolic C-terminal domain. One of these signals is the DXXLL-motif sequence DISLL, which controls transport between the TGN and endosomes via interaction with GGA proteins. Here we show that the DISLL sequence is embedded within a longer [DE]XXXL[LI]-motif sequence, DDISLL, which mediates internalization from the plasma membrane by interaction with the clathrin-associated, heterotetrameric adaptor protein 2 (AP-2) complex. Mutation of this signal or knockdown of either AP-2 or clathrin decreases endosomal localization and increases plasma membrane localization of BACE1. Remarkably, internalization-defective BACE1 is able to cleave an APP mutant that itself cannot be delivered to endosomes. The drug brefeldin A reversibly prevents BACE1-catalyzed APP cleavage, ruling out that this reaction occurs in the endoplasmic reticulum (ER) or ER-Golgi intermediate compartment. Taken together, these observations support the notion that BACE1 is capable of cleaving APP in late compartments of the secretory pathway.     

Study of the residues involved in the binding of β1 to β3 subunits in the sodium channel

The voltage-gated sodium channel (VGSC) is a complex, which is composed of one pore-forming α subunit and at least one β subunit. Up to now, five β subunits are known: β1/β1A, β1B, β2, β3, and β4, encoded by four genes (SCN1B∼SCN4B). It is critical to have a deep understanding of the interaction between β1 and β3 subunits, two subunits which frequently appear in many diseases concurrently. In this study, we had screened out the new template of β1 subunit for homology modelling, which shares higher similarity to β3. Docking studies of the β1 and β3 homology model were conducted, and likely β1 and β3 binding loci were investigated. The results revealed that β1–β3 is more likely to form a di-polymer than β1–β1 based on molecular interaction analysis, including potential energy analysis, Van der Waals (VDW) energy analysis and electrostatic energy analysis, and in addition, consideration of the hydrogen bonds and hydrophobic contacts that are involved. Based on these analyses, the residues His122 and Lys140 of β1 and Glu 66, Asn 131, Asp 118, Glu 120, Glu133, Asn135, Ser 137 of β3 were predicted to play a functional role.     

Kinetic advantage of the interaction between the fatty acid β-oxidation enzymes and the complexes of the respiratory chain

Respiration-linked oxidation of 3-hydroxybutyryl-CoA, crotonyl-CoA and saturated fatty acyl (C₄, C₈ and C₁₄)-CoA esters was studied in different mitochondrial preparations. Oxidation of acyl-CoA esters was poor in intact mitochondria; however, it was significant, as well as, NAD⁺ and CoA-dependent in gently and in vigorously sonicated mitochondria. The respiration-linked oxidation of crotonyl-CoA and 3-hydroxybutyryl-CoA proceeded at much higher rates (over 700%) in gently disrupted mitochondria than in completely disrupted mitochondria. The redox dye-linked oxidation of crotonyl-CoA (with inhibited respiratory chain) was also higher in gently disrupted mitochondria (149%) than in disrupted ones. During the respiration-linked oxidation of 3-hydroxybutyryl-CoA the steady-state NADH concentrations in the reaction chamber were determined, and found to be 8 μM in gently sonicated and 15 μM in completely sonicated mitochondria in spite of the observation that the gently sonicated mitochondria oxidized the 3-hydroxybutyryl-CoA much faster than the completely sonicated mitochondria. The NAD⁺-dependence of 3-hydroxybutyryl-CoA oxidation showed that a much smaller NAD⁺ concentration was enough to half-saturate the reaction in gently disrupted mitochondria than in completely disrupted ones. Thus, these observations indicate the positive kinetic consequence of organization of β-oxidation enzyme in situ. Respiration-linked oxidation of bytyryl-, oxtanoyl- and palmitoyl-CoA was also studied and these CoA intermediates were oxidized at approx. 50% of the rate of crotonyl- and 3-hydroxybutyryl-CoA in the gently disrupted mitochondria. In vigorously disrupted mitochondria the oxidation rate of these saturated acyl-CoA intermediates was hardly detectable indicating that the connection between the acyl-CoA dehydrogenase and the respiratory chain had been disrupted.     

Estimation of binding parameters for the protein–protein interaction using a site-directed spin labeling and EPR spectroscopy

Sensitivity of the electron paramagnetic resonance (CW EPR) to molecular tumbling provides potential means for studying processes of molecular association. It uses spin-labeled macromolecules, whose CW EPR spectra may change upon binding to other macromolecules. When a spin-labeled molecule is mixed with its liganding partner, the EPR spectrum constitutes a linear combination of spectra of the bound and unbound ligand (as seen in our example of spin-labeled cytochrome c ₂ interacting with cytochrome bc ₁ complex). In principle, the fraction of each state can be extracted by the numerical decomposition of the spectrum; however, the accuracy of such decomposition may often be compromised by the lack of the spectrum of the fully bound ligand, imposed by the equilibrium nature of molecular association. To understand how this may affect the final estimation of the binding parameters, such as stoichiometry and affinity of the binding, a series of virtual titration experiments was conducted. Our non-linear regression analysis considered a case in which only a single class of binding sites exists, and a case in which classes of both specific and non-specific binding sites co-exist. The results indicate that in both models, the error due to the unknown admixture of the unbound ligand component in the EPR spectrum causes an overestimation of the bound fraction leading to the bias in the dissociation constant. At the same time, the stoichiometry of the binding remains relatively unaffected, which overall makes the decomposition of the EPR spectrum an attractive method for studying protein–protein interactions in equilibrium. Our theoretical treatment appears to be valid for any spectroscopic techniques dealing with overlapping spectra of free and bound component.     

Functional consequences of retro-inverso isomerization of a miniature protein inhibitor of the p53–MDM2 interaction

Peptide retro-inverso isomerization is thought to be functionally neutral and has been widely used as a tool for designing proteolytically stable d-isomers to recapitulate biological activities of their parent l-peptides. Despite success in a wide range of applications, exceptions amply exist that clearly defy this rule of thumb when parent l-peptides adopt an α-helical conformation in their bound state. The detrimental energetic effect of retro-inverso isomerization of an α-helical l-peptide on its target protein binding has been estimated to be 3.0–3.4 kcal/mol. To better understand how the retro-inverso isomer of a structured protein works at the molecular level, we chemically synthesized and functionally characterized the retro-inverso isomer of a rationally designed miniature protein termed stingin of 18 amino acid residues, which adopts an N-terminal loop and a C-terminal α-helix stabilized by two intra-molecular disulfide bridges. Stingin emulated the transactivation peptide of the p53 tumor suppressor protein and bound with high affinity and via its C-terminal α-helix to MDM2 and MDMX—the two negative regulators of p53. We also prepared the retro isomer and d-enantiomer of stingin for comparative functional studies using fluorescence polarization and surface plasmon resonance techniques. We found that retro-inverso isomerization of l-stingin weakened its MDM2 binding by 720 fold (3.9 kcal/mol); while enantiomerization of l-stingin drastically reduced its binding to MDM2 by three orders of magnitude, sequence reversal completely abolished it. Our findings demonstrate the limitation of peptide retro-inverso isomerization in molecular mimicry and reinforce the notion that the strategy works poorly with biologically active α-helical peptides due to inherent differences at the secondary and tertiary structural levels between an l-peptide and its retro-inverso isomer despite their similar side chain topologies at the primary structural level.¹