Protonation of a histidine copper ligand in fern plastocyanin

Plastocyanin is a small blue copper protein that shuttles electrons as part of the photosynthetic redox chain. Its redox behavior is changed at low pH as a result of protonation of the solvent-exposed copper-coordinating histidine. Protonation and subsequent redox inactivation could have a role in the down regulation of photosynthesis. As opposed to plastocyanin from other sources, in fern plastocyanin His90 protonation at low pH has been reported not to occur. Two possible reasons for that have been proposed: pi-pi stacking between Phe12 and His90 and lack of a hydrogen bond with the backbone oxygen of Gly36. We have produced this fern plastocyanin recombinantly and examined the properties of wild-type protein and mutants Phe12Leu, Gly36Pro, and the double mutant with NMR spectroscopy, X-ray crystallography, and cyclic voltammetry. The results demonstrate that, contrary to earlier reports, protonation of His90 in the wild-type protein does occur in solution with a pKa of 4.4 (+/-0.1). Neither the single mutants nor the double mutant exhibit a change in protonation behavior, indicating that the suggested interactions have no influence. The crystal structure at low pH of the Gly36Pro variant does not show His90 protonation, similar to what was found for the wild-type protein. The structure suggests that movement of the imidazole ring is hindered by crystal contacts. This study illustrates a significant difference between results obtained in solution by NMR and by crystallography.     

Investigation of the copper binding site on the human islet amyloid polypeptide hormone

The metal ion binding sites of human islet amyloid polypeptide (hIAPP) have been investigated to explain the biological activity difference in the fibril formation process. The structures of [hIAPP...Cu (or Al)](n+) and [hIAPP17-30...Cu]2+ complex were investigated by electrospray ionization-mass spectrometry (ESI-MS). The fragmentation patterns of [hIAPP...Cu [or Al)](n+) and [hIAPP17-30...Cu]2+ complex were analyzed by tandem mass spectrometry (MS/MS) and multi-stage mass spectrometry (MS3) spectra. The [hIAPP+Cu+H]3+, [hIAPP+Al+H]4+ and [hIAPP17-30+Cu]2+ complexes were observed in MS spectra. The Cu binding site of hIAPP is suggested to be the N22-F-G-A-I26 part for the [hIAPP+Cu+H]3+ gas-phase complex. The original hIAPP conformation was supposed to be changed by the interaction between the Cu ion and the N22-F-G-A-I26 part in the [hIAPP+Cu+H]3+ gas-phase complex.     

Probing the dynamic nature of water molecules and their influences on ligand binding in a model binding site

The model binding site of the cytochrome c peroxidase (CCP) W191G mutant is used to investigate the structural and dynamic properties of the water network at the buried cavity using computational methods supported by crystallographic analysis. In particular, the differences of the hydration pattern between the uncomplexed state and various complexed forms are analyzed as well as the differences between five complexes of CCP W191G with structurally closely related ligands. The ability of docking programs to correctly handle the water molecules in these systems is studied in detail. It is found that fully automated prediction of water replacement or retention upon docking works well if some additional preselection is carried out but not necessarily if the entire water network in the cavity is used as input. On the other hand, molecular interaction fields for water calculated from static crystal structures and hydration density maps obtained from molecular dynamics simulations agree very well with crystallographically observed water positions. For one complex, the docking and MD results sensitively depend on the quality of the starting structure, and agreement is obtained only after redetermination of the crystal structure and refinement at higher resolution.     

Evaluation of quail egg white riboflavin binding protein as a chiral selector in high-performance liquid chromatography and capillary electrophoresis

A new chiral stationary phase for high-performance liquid chromatography of quail egg white riboflavin binding protein is presented. Several chiral acidic, basic and uncharged drugs were analysed and the influence of the mobile phase's parameters on the retention times and enantioselectivity was evaluated. On the basis of the results obtained, the same protein was studied as a background electrolyte additive in free solution capillary electrophoresis, in order to evaluate if capillary electrophoresis (CE) could be used as a rapid scouting technique for screening the enantioselectivity of novel proteins without immobilisation on a solid support. To investigate if it is possible to directly compare the results obtained by each technique, the CE experiments were planned on the basis of both the findings and ideas originated in liquid chromatography.     

Molecular modelling prediction of ligand binding site flexibility

We have investigated the efficacy of generating multiple sidechain conformations using a rotamer library in order to find the experimentally observed ligand binding site conformation of a protein in the presence of a bound ligand. We made use of a recently published algorithm that performs an exhaustive conformational search using a rotamer library to enumerate all possible sidechain conformations in a binding site. This approach was applied to a dataset of proteins whose structures were determined by X-ray and NMR methods. All chosen proteins had two or more structures, generally involving different bound ligands. By taking one of these structures as a reference, we were able in most cases to successfully reproduce the experimentally determined conformations of the other structures, as well as to suggest alternative low-energy conformations of the binding site. In those few cases where this procedure failed, we observed that the bound ligand had induced a high-energy conformation of the binding site. These results suggest that for most proteins that exhibit limited backbone motion, ligands tend to bind to low energy conformations of their binding sites. Our results also reveal that it is possible in most cases to use a rotamer search-based approach to predict alternative low-energy protein binding site conformations that can be used by different ligands. This opens the possibility of incorporating alternative binding site conformations to improve the efficacy of docking and structure-based drug design algorithms.     

A single hydroxyl group governs ligand site selectivity in human ileal bile acid binding protein

The recognition between proteins and their native ligands is fundamental to biological function. In vivo, human ileal bile acid binding protein (I-BABP) encounters a range of bile salts that vary in the number and position of steroidal hydroxyl groups and the presence and type of side-chain conjugation. Therefore, it is necessary to understand how chemical variability in the ligand affects the energetic and structural aspects of its recognition. Here we report studies of the binding site selectivity of I-BABP for glycocholic (GCA) and glycochenodeoxycholic (GCDA) acids using isotope-enriched bile salts along with two-dimensional heteronuclear NMR methods. When I-BABP is presented with either GCA or GCDA alone, the ligands bind to both sites. However, when presented with an equimolar mixture of the two bile salts, GCDA binds exclusively to site 1 and GCA to site 2. This remarkable selectivity is governed by the presence or absence of a single hydroxyl group at the C-12 position of the steroid tetracycle. The basis for this site selectivity appears to be energetic rather then steric.     

Geometrically feasible binding modes of a flexible ligand molecule at the receptor site

A very efficient algorithm for determining the geometrically feasible binding modes of a flexible ligand molecule at the receptor site is presented. It is based on distance geometry but maintains the requirements of three dimensions. The distance geometry manipulation can superimpose two bodies without explicitly calculating the necessary rigid rotation and translation. The whole conformation space of a flexible molecule can be efficiently examined by considering only a finite number of conformational points. The method is suitable only when the criterion for superposition is some minimum distance limit. It cannot, however, give the exact distance between two points in two different bodies.     

A molecular graphics study exploring a putative ligand binding site of theβ-adrenoceptor

Summary The recent elucidation of the primary structure of the cell membrane-bound-adrenoceptor has prompted us to explore putative ligand binding sites on this physiologically important receptor. By minimizing the energies of the prototype ligand propranolol, (part of) the receptor and the proposed ligand-receptor complex with the aid of force field and quantum chemical calculations, we identified amino acid residue Trp³¹³ as a highly probable candidate for interaction with the aromatic moiety of propranolol. The charge distribution on the indole nucleus of another-blocker, pindolol, with higher affinity for the-adrenoceptor, enables an even stronger interaction with the tryptophan residue. The carboxylic amino acid residue Glu³⁰⁶, located near the extracellular space of the cell membrane, interacts favorably with the positively charged nitrogen atom in the aliphatic side chain of the ligands. Finally, this putative model is discussed in the light of recent findings in mutagenesis studies, and compared to other ideas with respect to ligand-receptor interactions.     

Calculation of absolute ligand binding free energy to a ribosome-targeting protein as a function of solvent model

A comparative analysis is provided of the effect of different solvent models on the calculation of a potential of mean force (PMF) for determining the absolute binding affinity of the small molecule inhibitor pteroic acid bound to ricin toxin A-chain (RTA). Solvent models include the distance-dependent dielectric constant, several different generalized Born (GB) approximations, and a hybrid explicit/GB-based implicit solvent model. We found that the simpler approximation of dielectric screening and a GB model, with Born radii fitted to a switching-window dielectric-boundary surface Poisson solvent model, severely overpredicted the binding affinity as compared to the experimental value, estimated to range from -4.4 to -6.0 kcal/mol. In contrast, GB models that are parametrized to fit the Lee-Richards molecular surface performed much better, predicting binding free energy within 1-3 kcal/mol of experimental estimates. However, the predicted free-energy profiles of these GB models displayed alternative binding modes not observed in the crystal structure. Finally, the most rigorous and computationally costly approach in this work, which used a hybrid explicit/implicit solvent model, correctly determined a binding funnel in the PMF near the crystallographic bound state and predicted an absolute binding affinity that was 2 kcal/mol more favorable than the estimated experimental binding affinity.     

Spectroscopic and density functional studies of the red copper site in nitrosocyanin: role of the protein in determining active site geometric and electronic structure

The electronic structure of the red copper site in nitrosocyanin is defined relative to that of the well understood blue copper site of plastocyanin by using low-temperature absorption, circular dichroism, magnetic circular dichroism, resonance Raman, EPR and X-ray absorption spectroscopies, combined with DFT calculations. These studies indicate that the principal electronic structure change in the red copper site is the sigma rather than the pi donor interaction of the cysteine sulfur with the Cu 3d(x2-y2) redox active molecular orbital (RAMO). Further, MCD data show that there is an increase in ligand field strength due to an increase in coordination number, whereas resonance Raman spectra indicate a weaker Cu-S bond. The latter is supported by the S K-edge data, which demonstrate a less covalent thiolate interaction with the RAMO of nitrosocyanin at 20% relative to plastocyanin at 38%. EXAFS results give a longer Cu-S(Cys) bond distance in nitrosocyanin (2.28 A) compared to plastocyanin (2.08 A) and also show a large change in structure with reduction of the red copper site. The red copper site is the only presently known blue copper-related site with an exogenous water coordinated to the copper. Density functional calculations reproduce the experimental properties and are used to determine the specific protein structure contributions to exogenous ligand binding in red copper. The relative orientation of the CuNNS and the CuSC(beta) planes (determined by the protein sequence) is found to be key in generating an exchangeable coordination position at the red copper active site. The exogenous water ligation at the red copper active site greatly increases the reorganization energy (by approximately 1.0 eV) relative to that of the blue copper protein site, making the red site unfavorable for fast outer-sphere electron transfer, while providing an exchangeable coordination position for inner-sphere electron transfer.     

Investigation on the binding site in heparin by spectrophotometry

Heparin has a variety of biological activities, most of them due to heparin’s high sulfate groups. To gain insight into the mechanism of activation of the spectroscopic probe with sulfate groups of heparin in vitro, we have used a cationic dye by a spectrophotometric method. It is considered that the combination of heparin with methylene blue is due to noncovalent binding forces. Dye binding requires an organic chain structure form with sulfate groups. The solution equilibria of the reaction system are discussed. A new linear regression equation has been proposed, in which the maximum binding number N expresses the binding ability of methylene blue (MB) with sulfate groups of heparin. The linear regression equation can estimate this parameter.     

蛋白质分子与配体结合过程构象变化的研究

蛋白质分子与配体的作用模式主要有直接的环区结合及铰链式结合两种方式。针对这两种不同的作用方式,我们提出采用不同的策略进行结合过程的构象研究。对于直接的环区结合模式,通过建立环区主链构象库,来实现蛋白质环区与配体的准柔性对接,并以链霉抗生物素蛋白体系为例对构象库建立的可行性进行了验证计算。对铰链结合方式,采用分步对接的方法进行计算,并具体应用于HIV蛋白酶与其小分子配体的结合过程。计算结果表明,这两种处理方法分别能较好地模拟不同类型的蛋白质与配体结合的的构象变化。     

Effect of a tridentate ligand on the structure, electronic structure, and reactivity of the copper(I) nitrite complex: role of the conserved three-histidine ligand environment of the type-2 copper site in copper-containing nitrite reductases

It is postulated that the copper(I) nitrite complex is a key reaction intermediate of copper containing nitrite reductases (Cu-NiRs), which catalyze the reduction of nitrite to nitric oxide (NO) gas in bacterial denitrification. To investigate the structure-function relationship of Cu-NiR, we prepared five new copper(I) nitrite complexes with sterically hindered tris(4-imidazolyl)carbinols [Et-TIC = tris(1-methyl-2-ethyl-4-imidazolyl)carbinol and iPr-TIC = tris(1-methyl-2-isopropyl-4-imidazolyl)carbinol] or tris(1-pyrazolyl)methanes [Me-TPM = tris(3,5-dimethyl-1-pyrazolyl)methane; Et-TPM = tris(3,5-diethyl-1-pyrazolyl)methane; and iPr-TPM = tris(3,5-diisopropyl-1-pyrazolyl)methane]. The X-ray crystal structures of all of these copper(I) nitrite complexes were mononuclear eta(1)-N-bound nitrite complexes with a distorted tetrahedral geometry. The electronic structures of the complexes were investigated by absorption, magnetic circular dichroism (MCD), NMR, and vibrational spectroscopy. All of these complexes are good functional models of Cu-NiR that form NO and copper(II) acetate complexes well from reactions with acetic acid under anaerobic conditions. A comparison of the reactivity of these complexes, including previously reported (iPr-TACN)Cu(NO2) [iPr-TACN = 1,4,7-triisopropyl-1,4,7-triazacyclononane], clearly shows the drastic effects of the tridentate ligand on Cu-NiR activity. The copper(I) nitrite complex with the Et-TIC ligand, which is similar to the highly conserved three-histidine ((His)3) ligand environment in the catalytic site of Cu-NiR, had the highest Cu-NiR activity. This result suggests that the (His)3 ligand environment is essential for acceleration of the Cu-NiR reaction. The highest Cu-NiR activity for the Et-TIC complex can be explained by the structural and spectroscopic characterizations and the molecular orbital calculations presented in this paper. Based on these results, the functional role of the (His)3 ligand environment in Cu-NiR is discussed.     

Identification of a novel high affinity copper binding site in the APP(145-155) fragment of amyloid precursor protein

The copper(II) binding features of the APP(145-155) and APP(145-157) fragments of the amyloid precursor protein, Ac-Glu-Thr-His-Leu-His-Trp-His-Thr-Val-Ala-Lys-NH2 and Ac-Glu-Thr-His-Leu-His-Trp-His-Thr-Val-Ala-Lys-Glu-Thr-NH2 were studied by NMR spectroscopy and NMR findings were supported by UV-vis, CD and EPR spectra. Potentiometric measurements were performed only for the more soluble Ac-Glu-Thr-His-Leu-His-Trp-His-Thr-Val-Ala-Lys-Glu-Thr-NH2 peptide fragment. The following was shown: (i) the imidazole rings of all the three His residues are involved in metal coordination; (ii) metal binding induces ionisation of Leu-148 and His-149 amide nitrogens that complete the donor set to copper(II) in the species dominant at neutral pH; (iii) the unusual coordination scheme of the His-Xxx-His-Xxx-His consensus sequence justifies the high specificity for Cu(II) when compared to SOD-like or albumin-like peptides or even in amyloid Abeta fragments. The present findings may represent the key for interpreting the observed requirement of His residues conservation for the redox cycling between Cu(II) and Cu(I) by soluble APP.     

Backbone dynamics in the DNA HhaI protein binding site

The dynamics of the phosphodiester backbone in the [5'-GCGC-3'] 2 moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs nonstereospecifically deuterated on the 5' methylene group of nucleotides within the [5'-GCGC-3'] 2 moiety indicated that all of these positions are structurally flexible. Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone and furanose ring of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through not only a loss of mobility but also a change of direction of motion. These NMR data indicate that the [5'-GCGC-3'] 2 moiety is dynamic, with the largest amplitude motions occurring nearest the methylation site. The change of orientation of this moiety in DNA upon methylation may make the molecule less amenable to binding to the HhaI endonuclease.     

Evaluation of quail egg white riboflavin binding protein as a chiral selector in capillary electrophoresis by applying a modified partial filling technique.

A preliminary evaluation of the enantioselective properties of quail egg yolk riboflavin binding protein (qRfBP) was carried out in capillary electrophoresis by using the complete filling technique. The most promising results obtained by this screening of nineteen chiral drugs were singled out with the aim of optimizing enantiomer separations by applying the partial filling technique, which allows operating at much higher protein concentrations without detection problems. The building of the separation zone in the partial filling technique has been modified in order to enable on-line monitoring, before each run, of the actual protein plug application velocity and, consequently, the building of a plug of the desired length. The electrophoretic conditions chosen gave opposite migration directions for the chiral selector and the analytes, with qRfBP migrating away from the detector. A polyvinyl alcohol-coated capillary was first totally filled with protein and the optimal plug length was obtained by further applying negative pressure together with positive voltage for the time needed. Separations of basic drugs were optimized by using protein concentrations ranging from 200 microM up to 900 microM and different plug lengths, while the running buffer pH (6.0), temperature (25 degrees C) and operating voltage (+20 kV) were kept constant. The enantioresolution of all solutes was affected by both the chiral selector concentration and protein plug length. Baseline separations were obtained for oxprenolol, prilocaine and bupivacaine.