Studies on ligand binding of kidney bean leghemoglobin.

Absorption spectra of different ligand derivative;s of kidney bean leghemoglobin alpha have been recorded. The effect of pH on the absorption spectra of kidney bean leghemoglobin alpha has been studied. The pK of the acid-alkaline transition of the heme-linked water molecule is 8.25 and the pK for the acid dissociation of the heme group is 4.03. Affinities of kidney bean leghemoglobin for two different types of ligands have been studied in comparison with soybean leghemoglobins alpha and c and sperm whale myoglobin. All these leghemoglobins have similar affinities for the small anionic ligand fluoride ion, and they are only slightly more accessible to this ligand than is sperm whale myoglobin. Differences in the primary structure or in conformation of these proteins are reflected in the affinity for the bulky ligand imidazole. The accessibility to imidazole increases in the order sperm whale myoglobin less than soybean Lbalpha less than soybean Lbc less than kidney bean Lbalpha, and also low spin Lbalpha less than high spin Lbalpha. The results are discussed with respect to the amino acid sequences of the leghemoglobins.     

基于导电聚吡咯生物电化学传感器的研究进展

导电聚合物具有良好的导电性能,可以作为分子导线使电子在生物活性物质与电极间直接传递,是构建生物传感器的一种新型材料.聚吡咯(PPy)具有导电性、生物相容性、易固定等特点,在传感器中用于固定生物活性物质有着良好的应用前景.该文简要介绍了导电聚吡咯的合成方法及掺杂机理,重点评述了聚吡咯用于固定生物活性物质构建生物传感器的多...     

Effect of ligand binding on the intraminimum dynamics of proteins

Effects of ligand binding on protein dynamics are studied via molecular dynamics (MD) simulations on two different enzymes, dihydrofolate reductase (DHFR) and triosephosphate isomerase (TIM), in their unliganded (free) and liganded states. Domain motions in MD trajectories are analyzed by collectivities and rotation angles along the principal components (PCs). DHFR in the free state has well‐defined domain rotations, whereas rotations are slightly damped in the binary complex with nicotinamide adenine dinucleotide phosphate (NADPH), and remarkably distorted in the presence of NADP⁺, showing that NADP⁺ is solely responsible for the loss of correlation of the domains in DHFR. Although mean square fluctuations of MD simulations in the same PC subspaces are similar for different ligation states, linear stochastic time series models show that backbone flexibility along the first five PCs is decreased upon NADPH and NADP⁺ binding in subpicosecond scale. This shows that mobility of the protein along the PCs is closely related with intraminimum dynamics, and alterations in ligation states may change the intraminimum dynamics significantly. Low vibrational frequencies of the alpha‐carbon atoms of DHFR are determined from the time series models of a larger number of low indexed PCs, and it is found that number of modes in the lowest frequencies is reduced upon ligand binding. A similar result is obtained for TIM in the unliganded and dihydroxyacetone phosphate bound states. We suggest that stochastic time series modeling is a promising method to be used in determining subtle perturbations in protein dynamics. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Influence of ligand binding on the conformation of Torpedo californica acetylcholinesterase

With the aim of identifying structural changes in acetylcholinesterase, induced by ligand binding, we use a completely automatic procedure to analyse the differences between the backbone conformation of the free enzyme and those in eight complexes of Torpedo californica acetylcholinesterase, with various quaternary ammonium ligands, and with the protein inhibitor fasciculin. In order to discriminate between structural changes due to ligand binding and those arising from model imprecision, we also examine protein–ligand and protein–water contacts. Except for the peptide flip in the complex with huperzine A, the backbones of other complexes with quaternary ammonium ligands display negligible changes relative to the free enzyme. Another exception is the complex with the bisquaternary ammonium ligand decamethonium, where several loops display above average deformations, but only two, those spanning residues 334–348 and residues 277–304, seem to move as a result of ligand binding. Movement of the ω loop (residues 61–95) is detected only in the complex with the protein fasciculin. Received: 21 July 2000 / Accepted: 18 September 2000 / Published online: 28 February 2001     

Influence of ligand density on antibody binding capacity of cation-exchange adsorbents

Adsorption properties of a set of polymethacrylate-based cation exchangers designed for purification of monoclonal antibodies were investigated. The materials differed significantly in the density of sulphoisobutyl ligand groups. The ligand density had a pronounced effect on the static adsorption capacity of a polyclonal human immunoglobulin G. An optimal ligand density was observed at any pH and NaCl concentration tested when sharp optima were observed at low pH and ionic strength values. This was caused by effective clogging of pore mouth at high ligand densities. An anomalous effect of ionic strength was observed for the adsorbents with the high ligand density when the adsorption capacity increased with the addition of NaCl at low pH.     

Identification of ligand binding sites on proteins using a multi-scale approach

Identification of a ligand binding site on a protein is pivotal to drug discovery. To date, no reliable and computationally feasible general approach to this problem has been published. Here we present an automated efficient method for determining binding sites on proteins for potential ligands without any a priori knowledge. Our method is based upon the multiscale concept where we deal with a hierarchy of models generated using a k-means clustering algorithm for the potential ligand. This is done in a simple approach whereby a potential ligand is represented by a growing number of feature points. At each increasing level of detail, a pruning of potential binding site is performed. A nonbonding energy function is used to score the interactions between molecules at each step. The technique was successfully employed to seven protein-ligand complexes. In the current paper we show that the algorithm considerably reduces the computational effort required to solve this problem. This approach offers real opportunities for exploiting the large number of structures that will evolve from structural genomics.     

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.     

Assessment of programs for ligand binding affinity prediction

The prediction of the binding free energy between a ligand and a protein is an important component in the virtual screening and lead optimization of ligands for drug discovery. To determine the quality of current binding free energy estimation programs, we examined FlexX, X-Score, AutoDock, and BLEEP for their performance in binding free energy prediction in various situations including cocrystallized complex structures, cross docking of ligands to their non-cocrystallized receptors, docking of thermally unfolded receptor decoys to their ligands, and complex structures with "randomized" ligand decoys. In no case was there a satisfactory correlation between the experimental and estimated binding free energies over all the datasets tested. Meanwhile, a strong correlation between ligand molecular weight-binding affinity correlation and experimental predicted binding affinity correlation was found. Sometimes the programs also correctly ranked ligands' binding affinities even though native interactions between the ligands and their receptors were essentially lost because of receptor deformation or ligand randomization, and the programs could not decisively discriminate randomized ligand decoys from their native ligands; this suggested that the tested programs miss important components for the accurate capture of specific ligand binding interactions.     

A review on the ligand binding studies by isothermal titration calorimetry

Thermodynamics of biomacromolecule ligand interaction is very important to understand the structure function relationship in proteins. One of the most powerful techniques useful to obtain additional information about the structure of proteins in biophysical chemistry field is isothermal titration calorimetry (ITC). An ITC experiment is a titration of a biomacromolecule solution by a solution containing a reactant (ligand) at constant temperature to obtain the exchanged heat of the reaction. The total concentration of ligand is the independent variable under experimental control. There are many reports on data analysis for ITC to find the number of binding sites (g), the equilibrium constant (K), the Gibbs free energy of binding process (ΔG), the enthalpy of binding (ΔH) and the entropy of binding (ΔS). Moreover, ITC gives information about the type of reaction, electrostatic and hydrophobic interactions, including determination of cooperativity characterization in binding process by calculating the Hill coefficient (n). A double reciprocal plot and a graphical fitting method are two simple methods used in the enzyme inhibition and metal binding to a protein. Determination of a binding isotherm needs more ITC experiments and more complex data analysis. Protein denaturation by ligand includes two processes of binding and denaturation so that ITC data analysis are more complex. However, the enthalpy of denaturation process obtained by ITC help to understand the fine structure of a protein.     

Protein ligand docking based on empirical method for binding affinity estimation

An empirical protein-ligand binding affinity estimation method, SCORE, was incorporated into a popular docking program, DOCK4. The combined program, ScoreDock, was used to reconstruct the 200 protein-ligand complex structures and found to give good results for the complexes with high binding affinities. A quality assessment method for docking results from ScoreDock was developed based on the whole test set and tested by additionally selected complexes. The method significantly improves the docking accuracy and was shown to be reliable in docking quality assessment. As a docking tool in structural based drug design, ScoreDock can screen out final hits directly based on the predicted negative logarithms of dissociation equilibrium constants of protein-ligand complexes, and can explicitly deal with structure water molecules, as well as metal atoms.     

A note on the analysis of ligand binding by the 'double-logarithmic' plot

The double-logarithmic plot has been employed for the fluorometric analysis of ligand binding. This analysis is valid only for 1:1 binding of a non-aggregating ligand and if the plot is based on the free, and not the total, ligand concentration.     

Gaussian mapping of chemical fragments in ligand binding sites

We present a new approach to automatically define a quasi-optimal minimal set of pharmacophoric points mapping the interaction properties of a user-defined ligand binding site. The method is based on a fitting algorithm where a grid of sampled interaction energies of the target protein with small chemical fragments in the binding site is approximated by a linear expansion of Gaussian functions. A heuristic approximation selects from this expansion the smallest possible set of Gaussians required to describe the interaction properties of the binding site within a prespecified accuracy. We have evaluated the performance of the approach by comparing the computed Gaussians with the positions of aromatic sites found in experimental protein-ligand complexes. For a set of 53 complexes, good correspondence is found in general. At a 95% significance level, approximately 65% of the predicted interaction points have an aromatic binding site within 1.5 A. We then studied the utility of these points in docking using the program DOCK. Short docking times, with an average of approximately 0.18 s per conformer, are obtained, while retaining, both for rigid and flexible docking, the ability to sample native-like binding modes for the ligand. An average 4-5-fold speed-up in docking times and a similar success rate is estimated with respect to the standard DOCK protocol.     

Synthesis and ligand binding properties of triptycene-linked porphyrin arrays

Multiporphyrin arrays are a complex class of molecules with numerous potential applications in energy transfer, photomedicine, and light harvesting. We have developed a facile/versatile route to a class of triptycene-linked porphyrin arrays via both Suzuki and Sonogashira cross-coupling methods, which makes use of the rigid three-pronged orientation of triptycene to construct trimeric porphyrin arrays linked either in the meso or β-position with various linker groups. In order to understand the properties of these potential antenna systems and probe their potential applications, the coordination behavior of zinc(II) derivatives with mono- and bidentate N-donor ligands was investigated. Depending on ligand concentration, both one- and two-point binding was observed with a bidentate ligand. Also/in addition, different cavity sizes, obtained by the use of different linker groups, resulted in differences in the binding properties of each trimeric system.     

Energetics of ligand binding to the DNA minor groove

In the present work the decomposition of the total Gibbs free energy of ligand-DNA binding onto various physical terms was accomplished for the group of nine DNA minor groove binders (MGB ligands) differing in both structure and charge state. The decomposition protocol includes the analysis of the most complete set of physical factors known to contribute to the complexation process, viz. the net change in the number of degrees of freedom (translational, rotational, vibrations of the chemical bonds and vibrations of the ligand as a whole within the binding site), the conformational entropy, van der Waals, electrostatic and hydrophobic interactions, the polyelectrolyte contribution and the net effect of changes in the number of hydrogen bonds. All of these processes are further decomposed into the interaction with the solvent and the interaction of the ligand with DNA. The principal outcome of the decomposition is the possibility of performing a comparative analysis of the energetic contribution of various physical terms and provide an answer to the question concerning what physical factors stabilize or destabilize the complexes of MGB ligands with DNA.     

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.