Comparative studies of two transthyretin variants with protective effects on familial amyloidotic polyneuropathy: TTR R104H and TTR T119M

Recently, a new nonpathogenic transthyretin (TTR) variant-TTR R104H (TTR H104)-has been described in heterozygotic and compound heterozygotic individuals from a Japanese family with familial amyloidotic polyneuropathy (FAP). The compound heterozygotic individual, a carrier of TTR V30M (TTR M30) and TTR R104H (TTR M30/H104) presented a very mild form of FAP with slow progression of the disease. TTR and retinol binding protein (RBP) levels were found to be increased in serum from TTR H104 carriers. These characteristics are very similar to those found in compound heterozygotic carriers of TTR V30M-T119M (TTR M30/M119). To structurally compare these variants, we performed stability and thyroxine (T(4)) binding studies. TTR M30/H104 showed an increased resistance to dissociation into monomers similar to TTR M30/M119. This suggests that the His104 substitution has the same stabilizing effect on tetrameric TTR as the Met119 substitution. Concerning T(4) binding, TTR H104 presents a T(4) binding affinity lower than that of TTR M119, but still higher than normal TTR. However, TTR from the compound heterozygotic carrier of TTR M30/H104 presented a T(4) binding affinity lower than normal. The results indicate that the His 104 substitution induces structural alterations that increase the stability of the tetramer in compound heterozygotes for TTR M30 despite a lower affinity for T(4) binding. Thus, stability of TTR and binding affinity for T(4) may not be related. More detailed characterization of these variants is needed to clarify the structural alterations responsible for their increased stability.     

Structural and energetic insight into the isoform-selective inhibitors of tumour marker Hsp90 against Grp94

The heat shock protein 90 (Hsp90) represents a new avenue for cancer therapy. A novel benzolactam inhibitor, compound 31, showed a great selectivity for Hsp90α and Hsp90β against Grp94. However, the structural basis for the great selectivity of compound 31 for Hsp90α/β versus Grp94 remains poorly understood. In this study, we carried out molecular docking, molecular dynamics simulations and binding free energy calculations (MM-GBSA) to address the isoform selective property. Molecular docking studies indicated the different binding modes of the Hsp90 and Grp94 with compound 31. The MM-GBSA calculations revealed that the hydrophobic interactions between compound 31 and proteins contributed the most to the binding affinity and the Grp94/compound 31 complex could result in a less energy-favourable complex compared with the Hsp90α/compound 31 and the Hsp90β/compound 31 complexes. This may render the weak binding of compound 31 to the Grp94. This study may be helpful for the future design of novel and selective Hsp90 inhibitors.     

Intrinsic Thermodynamics and Structure Correlation of Benzenesulfonamides with a Pyrimidine Moiety Binding to Carbonic Anhydrases I,II, VII,XII, and XIII

The early stage of drug discovery is often based on selecting the highest affinity lead compound. To this end the structural and energetic characterization of the binding reaction is important. The binding energetics can be resolved into enthalpic and entropic contributions to the binding Gibbs free energy. Most compound binding reactions are coupled to the absorption or release of protons by the protein or the compound. A distinction between the observed and intrinsic parameters of the binding energetics requires the dissection of the protonation/deprotonation processes. Since only the intrinsic parameters can be correlated with molecular structural perturbations associated with complex formation, it is these parameters that are required for rational drug design. Carbonic anhydrase (CA) isoforms are important therapeutic targets to treat a range of disorders including glaucoma, obesity, epilepsy, and cancer. For effective treatment isoform-specific inhibitors are needed. In this work we investigated the binding and protonation energetics of sixteen [(2-pyrimidinylthio)acetyl]benzenesulfonamide CA inhibitors using isothermal titration calorimetry and fluorescent thermal shift assay. The compounds were built by combining four sulfonamide headgroups with four tailgroups yielding 16 compounds. Their intrinsic binding thermodynamics showed the limitations of the functional group energetic additivity approach used in fragment-based drug design, especially at the level of enthalpies and entropies of binding. Combined with high resolution crystal structural data correlations were drawn between the chemical functional groups on selected inhibitors and intrinsic thermodynamic parameters of CA-inhibitor complex formation.     

Substantial energetic improvement with minimal structural perturbation in a high affinity mutant antibody

Here, we compare an antibody with the highest known engineered affinity (K(d)=270 fM) to its high affinity wild-type (K(d)=700 pM) through thermodynamic, kinetic, structural, and theoretical analyses. The 4M5.3 anti-fluorescein single chain antibody fragment (scFv) contains 14 mutations from the wild-type 4-4-20 scFv and has a 1800-fold increase in fluorescein-binding affinity. The dissociation rate is approximately 16,000 times slower in the mutant; however, this substantial improvement is offset somewhat by the association rate, which is ninefold slower in the mutant. Enthalpic contributions to binding were found by calorimetry to predominate in the differential binding free energy. The crystal structure of the 4M5.3 mutant complexed with antigen was solved to 1.5A resolution and compared with a previously solved structure of an antigen-bound 4-4-20 Fab fragment. Strikingly, the structural comparison shows little difference between the two scFv molecules (backbone RMSD of 0.6A), despite the large difference in affinity. Shape complementarity exhibits a small improvement between the variable light chain and variable heavy chain domains within the antibody, but no significant improvement in shape complementarity of the antibody with the antigen is observed in the mutant over the wild-type. Theoretical modeling calculations show electrostatic contributions to binding account for -1.2 kcal/mol to -3.5 kcal/mol of the binding free energy change, of which -1.1 kcal/mol is directly associated with the mutated residue side-chains. The electrostatic analysis reveals several mechanistic explanations for a portion of the improvement. Collectively, these data provide an example where very high binding affinity is achieved through the cumulative effect of many small structural alterations.     

Docking study and binding free energy calculation of poly (ADP-ribose) polymerase inhibitors

Recently, the massively parallel computation of absolute binding free energy with a well-equilibrated system (MP-CAFEE) has been developed. The present study aimed to determine whether the MP-CAFEE method is useful for drug discovery research. In the drug discovery process, it is important for computational chemists to predict the binding affinity accurately without detailed structural information for protein / ligand complex. We investigated the absolute binding free energies for Poly (ADP-ribose) polymerase-1 (PARP-1) / inhibitor complexes, using the MP-CAFEE method. Although each docking model was used as an input structure, it was found that the absolute binding free energies calculated by MP-CAFEE are well consistent with the experimental ones. The accuracy of this method is much higher than that using molecular mechanics Poisson-Boltzmann / surface area (MM / PBSA). Although the simulation time is quite extensive, the reliable predictor of binding free energies would be a useful tool for drug discovery projects.     

Soluble HLA-DQ2 expressed in S2 cells copurifies with a high affinity insect cell derived protein

We here describe that soluble HLA-DQ2 (sDQ2) molecules, when expressed in Drosophila melanogaster S2 insect cells without a covalently tethered peptide, associate tightly with the D. melanogaster calcium binding protein DCB-45. The interaction between the proteins is stable in S2 cell culture and during affinity purification, which is done at high salt concentrations and pH 11.5. After affinity purification, the sDQ2/DCB-45 complex exists in substantial quantities next to a small amount of free heterodimeric sDQ2 and large amounts of aggregated sDQ2 free of DCB-45. Motivated by the stable complex formation and our interest in the development of reagents which inhibit HLA-DQ2 peptide binding, we have further characterized the sDQ2/DCB-45 interaction. Several lines of evidence indicate that an N-terminal fragment of DCB-45 is involved in the interaction with the peptide binding groove of sDQ2. Further mapping of this fragment of 54 residues identified a pentadecapeptide with high affinity for sDQ2 which may serve as a lead compound for the design of HLA-DQ2 blockers.     

Ligand interactions with lactose repressor protein and the repressor-operator complex: the effects of ionization and oligomerization on binding

Specific interactions between proteins and ligands that modify their functions are crucial in biology. Here, we examine sugars that bind the lactose repressor protein (LacI) and modify repressor affinity for operator DNA using isothermal titration calorimetry and equilibrium DNA binding experiments. High affinity binding of the commonly-used inducer isopropyl-beta,D-thiogalactoside is strongly driven by enthalpic forces, whereas inducer 2-phenylethyl-beta,D-galactoside has weaker affinity with low enthalpic contributions. Perturbing the dimer interface with either pH or oligomeric state shows that weak inducer binding is sensitive to changes in this distant region. Effects of the neutral compound o-nitrophenyl-beta,D-galactoside are sensitive to oligomerization, and at elevated pH this compound converts to an anti-inducer ligand with slightly enhanced enthalpic contributions to the binding energy. Anti-inducer o-nitrophenyl-beta,D-fucoside exhibits slightly enhanced affinity and increased enthalpic contributions at elevated pH. Collectively, these results both demonstrate the range of energetic consequences that occur with LacI binding to structurally-similar ligands and expand our insight into the link between effector binding and structural changes at the subunit interface.     

Novel 3-chloro-6-nitro-1H-indazole derivatives as promising antileishmanial candidates: synthesis,biological activity,and molecular modelling studies

An efficient pathway was disclosed for the synthesis of 3-chloro-6-nitro-1H-indazole derivatives by 1,3-dipolar cycloaddition on dipolarophile compounds 2 and 3. Faced the problem of separation of two regioisomers, a click chemistry method has allowed us to obtain regioisomers of triazole-1,4 with good yields from 82 to 90% were employed. Also, the antileishmanial biological potency of the compounds was achieved using an MTT assay that reported compound 13 as a promising growth inhibitor of Leishmania major. Molecular docking demonstrated highly stable binding with the Leishmania trypanothione reductase enzyme and produced a network of hydrophobic and hydrophilic interactions. Molecular dynamics simulations were performed for TryR-13 complex to understand its structural and intermolecular affinity stability in a biological environment. The studied complex remained in good equilibrium with a structure deviation of ∼1–3 Å. MM/GBSA binding free energies illustrated the high stability of TryR-13 complex. The studied compounds are promising leads for structural optimisation to enhance the antileishmanial activity.     

[3-cis-3,5-Dimethyl-(1-piperazinyl)alkyl]-bis-(4′-fluorophenyl)amine analogues as novel probes for the dopamine transporter

In a continuing effort to identify novel probes with which to study the dopamine transporter (DAT), we discovered that the σ receptor antagonist, rimcazole, binds with moderate affinity (K_i=224 nM) to the DAT. The results from previous SAR studies suggested that substitution of the carbazole ring system of rimcazole with bis-(4′-fluorophenyl)amine might improve binding affinity and selectivity for the DAT. Thus, a novel series of [3-cis-3,5-dimethyl-(1-piperazinyl)alkyl]bis-(4′-fluorophenyl)amines were synthesized. The most potent compound in this series (9b) displaced [³H]WIN 35,428 binding in rat caudate-putamen (K_i=17.6 nM) with comparable affinity to GBR 12909. Despite high-affinity binding at DAT, and structural similarity to GBR 12909, preliminary studies suggest 9b behaves more like rimcazole than GBR 12909 and does not demonstrate cocaine-like psychostimulant behavior in mice.     

Docking studies on kinesin spindle protein inhibitors: an important cooperative ‘minor binding pocket’ which increases the binding affinity significantly

Fifteen KSP inhibitors were docked into the receptor and the binding mode was analyzed for the first time. It was considered that in addition to the main binding pocket all the inhibitors merged in, there exists a cooperative minor binding pocket, which could be explored for significantly increased binding affinity. In addition, a good linear relationship between the biological activities and the lowest binding free energies has also been found. This may help in predicting the binding affinity of newly designed KSP inhibitors. Figure Two binding pockets considered after the analysis. Seven docked ligands (compound 1–7) were overlapped at the binding site. All inhibitors tested interacted with the main pocket, while CK0106023, interacted also with the cooperative minor pocket mainly surrounded by Arg221 and Ala218. Coloring of the binding site surface are different ends of each amino acid residue: blue represents amino group while red means carboxyl     

Exploring additivity effects of double mutations on the binding affinity of protein‐protein complexes

Additivity in binding affinity of protein‐protein complexes refers to the change in free energy of binding (ΔΔG_bind) for double (or multiple) mutations which is approximately equal to the sum of their corresponding single mutation ΔΔG_bind values. In this study, we have explored the additivity effect of double mutants, which shows a linear relationship between the binding affinity of double and sum of single mutants with a correlation of 0.90. However, the comparison of ΔΔG_bind values showed a mean absolute deviation of 0.86 kcal/mol, and 25.6% of the double mutants show a deviation of more than 1 kcal/mol, which are identified as non‐additive. The additivity effects have been analyzed based on the influence of structural features such as accessible surface area, long range order, binding propensity change, surrounding hydrophobicity, flexibility, atomic contacts between the mutations and distance between the 2 mutations. We found that non‐additive mutations tend to be closer to each other and have more contacts. We have also used machine learning methods to discriminate additive and non‐additive mutations using structure‐based features, which showed the accuracies in the range of 0.77–0.92 for protein‐protein complexes belonging to different functions. Further, we have compared the additivity effects of protein stability along with binding affinity and explored the similarities and differences between them. The results obtained in this study provide insights into the effects of various structural features on binding affinity of double mutants, and will aid the development of accurate methods to predict the binding affinity of double mutants.     

Binding thermodynamics of substituted diaminopyrimidine renin inhibitors

Renin is an aspartyl protease involved in the production of angiotensin II, a potent vasoconstrictor. Renin inhibitors can prevent blood vessel constriction and therefore could be useful for the treatment of hypertension. High-throughput screening efforts identified a small molecule renin inhibitor with a core substituted diaminopyrimidine ring. Parallel medicinal chemistry efforts based on this lead resulted in compound 1. A complex of 1 bound to renin was crystallized, and structural data were obtained by X-ray diffraction. The structure indicated that there were adjacent unoccupied binding pockets. Synthetic efforts were initiated to extend functionality into these pockets so as to improve affinity and adjust pharmacokinetic parameters. Thermodynamics data for inhibitor binding to renin were also collected using isothermal titration calorimetry. These data were used to help guide inhibitor optimization by suggesting molecular alterations to improve binding affinity from both thermodynamic and structural perspectives. The addition of a methoxypropyl group extending into the S3 subpocket improved inhibitor affinity and resulted in greater binding enthalpy. Initial additions to the pyrimidine ring template that extended into the large hydrophobic S2 pocket did not improve affinity and dramatically altered the thermodynamic driving force for the binding interaction. Binding of the core template was enthalpically driven, whereas binding of initial inhibitors with S2 extensions was both enthalpically and entropically driven but lost significant binding enthalpy. Additional electrostatic interactions were then incorporated into the S2 extension to improve binding enthalpy while taking advantage of the favorable entropy.     

Structural basis of substrate binding specificity revealed by the crystal structures of polyamine receptors SpuD and SpuE from Pseudomonas aeruginosa

The type III secretion system (T3SS) of Pseudomonas aeruginosa is a key virulence determinant whose expression is induced by polyamine signals from mammalian host. SpuD and SpuE were postulated to be spermidine-preferential binding proteins, which regulate the polyamine content in this bacterial pathogen. In this study, we found that SpuD is a putrescine-preferential binding protein, while SpuE binds to spermidine exclusively. We have determined the crystal structures of SpuD in free form and in complex with putrescine and SpuE in free form and in complex with spermidine. Upon ligand binding, SpuD and SpuE undergo an "open-to-closed" conformational switch with the resultant closed ligand-bound forms, SpuD-putrescine and SpuE-spermidine, similar to their Escherichia coli counterparts PotF-putrescine and PotD-spermidine, respectively. Structural comparison suggested that two aromatic residues, Trp271 of SpuE and Phe273 of SpuD in segment II region, are the key structural determinants for putrescine/spermidine recognition specificity. Mutagenesis combined with isothermal titration calorimetry showed that substitution of Trp271 by Phe enabled SpuE to gain substantial binding affinity for putrescine, while replacement of Phe273 by Trp reduced the binding affinity of SpuD toward putrescine by 250-fold. Altogether, these results revealed the molecular mechanism governing polyamine recognition specificity by SpuD and SpuE and provide the basis for further structural and functional studies of polyamine signal importation system in P. aeruginosa.     

Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations

Molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have been performed to explore the dynamic behaviors of cytochrome P450 2A6 (CYP2A6) binding with nicotine analogs (that are typical inhibitors) and to calculate their binding free energies in combination with Poisson–Boltzmann surface area (PBSA) calculations. The combined MD simulations and QM/MM-PBSA calculations reveal that the most important structural parameters affecting the CYP2A6-inhibitor binding affinity are two crucial internuclear distances, that is, the distance between the heme iron atom of CYP2A6 and the coordinating atom of the inhibitor, and the hydrogen-bonding distance between the N297 side chain of CYP2A6 and the pyridine nitrogen of the inhibitor. The combined MD simulations and QM/MM-PBSA calculations have led to dynamic CYP2A6-inhibitor binding structures that are consistent with the observed dynamic behaviors and structural features of CYP2A6-inhibitor binding, and led to the binding free energies that are in good agreement with the experimentally-derived binding free energies. The agreement between the calculated binding free energies and the experimentally-derived binding free energies suggests that the combined MD and QM/MM-PBSA approach may be used as a valuable tool to accurately predict the CYP2A6-inhibitor binding affinities in future computational design of new, potent and selective CYP2A6 inhibitors.     

SCORE: A New Empirical Method for Estimating the Binding Affinity of a Protein-Ligand Complex

A new method is presented to estimate the binding affinity of a protein-ligand complex with known three-dimensional structure. The method, SCORE, uses an empirical scoring function to describe the binding free energy, which includes terms to account for van der Waals contact, metal-ligand bonding, hydrogen bonding, desolvation effect, and deformation penalty upon the binding process. The coefficients of each term are obtained by multivariate regressional analysis of a diverse training set of 170 protein-ligand complexes. The final scoring function reproduces the binding free energies of the whole training set with a cross-validated deviation of 6.3 kJ/mol. The predictive ability of the function is further tested by a set of 11 endothiapepsin complexes and the internal consistency of the function is demonstrated in a stepwise procedure named Evolutionary Test. A major innovation of this method is the introduction of an atomic binding score which allows the researcher to inspect and optimize the lead compound rationally in a structure-based drug design scheme.     

AT4 receptor structure-binding relationship: N-terminal-modified angiotensin IV analogues

The effect of structural changes in the N-terminal amino acid of AIV, with respect to AT₄ receptor binding, was examined by competition with [¹²⁵I]AIV in bovine adrenal membranes. Analogues with modifications of the first residue -amino group possessed lower affinities than the primary amine-containing parent compound. Peptides with a residue 1 -carbon in the conformation exhibited poor affinity for the AT₄ receptor. Modifications of the residue 1 R-group demonstrate that a straight chain aliphatic moiety containing four carbons is optimal for receptor-ligand binding, as evidenced by the extremely high affinity of [Nle¹]AIV (K_i = 3.59±0.51 pM). Replacement of the 1–2 peptide bond of AIV with the methylene bond isostere Ψ (CH₂-NH), increased the K_i approximately fivefold, indicating that the peptide bond may be replaced wihle maintaining relatively high-affinity receptor binding.     

Structure-based virtual screening for insect ecdysone receptor ligands using MM/PBSA

The ecdysone receptor (EcR) is an insect nuclear receptor that is activated by the molting hormone, 20-hydroxyecdysone. Because synthetic EcR ligands disrupt the normal growth of insects, they are attractive candidates for new insecticides. In this study, the Molecular Mechanics/Poisson–Boltzmann Surface Area (MM/PBSA) method was used to predict the binding activity of EcR ligands. Validity analyses using 40 known EcR ligands showed that the binding activity was satisfactorily predicted when the ligand conformational free energy term was introduced. Subsequently, this MM/PBSA method was applied to structure-based hierarchical virtual screening, and 12 candidate compounds were selected from a database of 3.8 million compounds. Five of these compounds were active in a cell-based competitive binding assay. The most potent compound is a simple proline derivative with low micromolar binding activity, representing a valuable lead compound for further structural optimization.     

Structure of the D142N mutant of the family 18 chitinase ChiB from Serratia marcescens and its complex with allosamidin

Catalysis by ChiB, a family 18 chitinase from Serratia marcescens, involves a conformational change of Asp142 which is part of a characteristic D(140)XD(142)XE(144) sequence motif. In the free enzyme Asp142 points towards Asp140, whereas it rotates towards the catalytic acid, Glu144, upon ligand binding. Mutation of Asp142 to Asn reduced k(cat) and affinity for allosamidin, a competitive inhibitor. The X-ray structure of the D142N mutant showed that Asn142 points towards Glu144 in the absence of a ligand. The active site also showed other structural adjustments (Tyr10, Ser93) that had previously been observed in the wild-type enzyme upon substrate binding. The X-ray structure of a complex of D142N with allosamidin, a pseudotrisaccharide competitive inhibitor, was essentially identical to that of the wild-type enzyme in complex with the same compound. Thus, the reduced allosamidin affinity in the mutant is not caused by structural changes but solely by the loss of electrostatic interactions with Asp142. The importance of electrostatics was further confirmed by the pH dependence of catalysis and allosamidin inhibition. The pH-dependent apparent affinities for allosamidin were not correlated with k(cat), indicating that it is probably better to view the inhibitor as a mimic of the oxazolinium ion reaction intermediate than as a transition state analogue.     

Increase of (Ca2+ +Mg2+)-ATPase activity of renal basolateral membranes by platelet-derived growth factor through a specific receptor

Studies were made on the direct effect of platelet-derived growth factor (PDGF) on the high-affinity (Ca2+ +Mg2+)-ATPase, a membrane bound Ca2+-extrusion pump enzyme of the basolateral membranes (BLM) of canine kidney (Km for free Ca2+ = 1.0 x 10(-7) M, Vmax = 180 nmol Pi/mg/min). At 1 x 10(-7) M free Ca2+, PDGF (10(-10)-10(-8) M) stimulated the enzyme activity significantly. Addition of 5 - 200 microM suramin, a compound that blocks binding of PDGF to its receptors on cell membranes, inhibited the stimulatory effect of PDGF dose-dependently (IC50 = 40 microM). A high affinity specific receptor for PDGF (Kd = 4.4 x 10(-10) M, Bmax = 460 fmol/mg protein) was detected on BLM preparations by radioreceptor assay with 125I-PDGF and unlabelled PDGF. Suramin (10-1000 microM) also inhibited the binding of PDGF to BLM preparations dose-dependently. From these results, it is proposed that PDGF stimulates (Ca2+ +Mg2+)-ATPase activity of kidney BLM preparations by enhancing its affinity for free Ca2+ through a specific receptor.