Surface Enhanced Raman Spectroscopy Detection of Biomolecules Using EBL Fabricated Nanostructured Substrates

Fabrication and characterization of conjugate nano-biological systems interfacing metallic nanostructures on solid supports with immobilized biomolecules is reported. The entire sequence of relevant experimental steps is described, involving the fabrication of nanostructured substrates using electron beam lithography, immobilization of biomolecules on the substrates, and their characterization utilizing surface-enhanced Raman spectroscopy (SERS). Three different designs of nano-biological systems are employed, including protein A, glucose binding protein, and a dopamine binding DNA aptamer. In the latter two cases, the binding of respective ligands, D-glucose and dopamine, is also included. The three kinds of biomolecules are immobilized on nanostructured substrates by different methods, and the results of SERS imaging are reported. The capabilities of SERS to detect vibrational modes from surface-immobilized proteins, as well as to capture the protein-ligand and aptamer-ligand binding are demonstrated. The results also illustrate the influence of the surface nanostructure geometry, biomolecules immobilization strategy, Raman activity of the molecules and presence or absence of the ligand binding on the SERS spectra acquired.     

The exponential model for a regulatory enzyme: Its extension to describe the binding of two ligands

The exponential model for a regulatory enzyme (Ainsworth, 1977a) is extended to deal explicitly with the presence in solution of a second ligand. This is achieved by introducing exponential interaction coefficients which respectively describe how the affinity of the free and bound forms of the protein for the ligand depend on its fractional saturation by the second ligand. The basic equations, so derived, are applied to binding experiments where the ligands bind independently or competitively and to rate experiments where the ligands represent two substrates or one substrate and a modifier which may be either competitive or non-competitive in type. The conditions required to display linkage between the binding of the two ligands are established and it is also shown that rate data may display a maximum as one ligand concentration is varied at a fixed concentration of the other. The equations that are derived are tested by application to experimental data and the conditions that have to be met to justify such an application are discussed.     

Theoretical aspects of DNA-protein interactions: co-operative and non-co-operative binding of large ligands to a one-dimensional homogeneous lattice

The interaction of proteins binding non-specifically to DNA, as well as the properties of many other interacting ligand-lattice systems important in molecular biology, requires a fundamentally different type of theoretical analysis than that provided by the classical Scatchard independent-binding-site treatment. Exact and relatively simple equations describing the binding of both non-interacting and interacting (co-operative) ligands to a homogeneous one-dimensional lattice are derived in terms of ligand site size, intrinsic binding constant and ligand-ligand co-operativity (equations (10) and (15) in the text). The mathematical approach is based on simple conditional probabilities, and reveals some largely unrecognized characteristics of such lattice binding systems. The results indicate that the binding of any non-interacting ligand covering more than one lattice residue results in non-linear (convex downward) Scatchard plots. The introduction of positive ligand-ligand co-operativity antagonizes this non-linearity, and eventually leads to plots of the opposite curvature. The maxima, limiting slopes, and intercepts of such plots can be used to estimate the required binding parameters. The method can be extended to systems involving heterogeneous ligands, and some types of heterogeneous lattices. Procedures for applying the method to a variety of interacting systems are presented, and a preliminary analysis is carried out for some selected sets of data from the literature.     

Mono-, bi-, or tridentate ligands? The labeling of peptides with 99mTc-carbonyls

The labeling of targeting peptides with (99m)Tc is a useful concept for the diagnosis of various diseases such as cancer. Although in research for at least one decade, only a very few radiopharmaceuticals based on peptides are in clinical use. The difficulty of labeling, and the resulting authenticity of the new vector, is largely responsible for this observation. In this overview, we present an alternate strategy based on the organometallic fac-[(99m)Tc(CO)(3)](+) core for introducing (99m)Tc in biomolecules in general and in peptides in particular. The three coordination sites available in [(99m)Tc(OH(2))(3)(CO)(3)](+) can be occupied with many different ligand types, pendant to a biomolecule and serving as the anchor group for labeling. This makes the appropriate choice difficult. We intend to present some useful concepts for the practice. Monodentate chelators are robust but bear the risk of multiple binding of biomolecules. Coordinating a bidentate ligand of choice prior to labeling bypasses this problem and enables a systematic drug discovery by variation of the bidentate ligand. Bidentate ligands attached to the biomolecule are stronger but occasionally require protection of the remaining site by a monodentate ligand. Both approaches refer to a mixed-ligand [2+1] approach. Tridentate chelators are the most efficient but need some protecting group chemistry in order to achieve selectivity for the coupling process. Examples with cysteine and histidine are presented. This article aims to provide versatile and reproducible approaches for the labeling of biomolecules while not focusing on particular systems. It should be left to the readers to derive a strategy for their own peptide.     

亲和层析技术在生物科学中的应用及发展

近几十年来,亲和层析技术发展十分迅速,广泛应用于生物分子(如结合蛋白、酶、抑制剂、抗原、抗体、激素、激素受体、糖蛋白、核酸及多糖类等)及组织(如细胞、细胞器、病毒等)的分离和纯化,是蛋白质组学研究中重要的技术之一.介绍了亲和层析的基本类型及配体合成的研究进展,概述了亲和层析技术在蛋白质组学以及在其他方面的应用和发展动态.     

92 Novel interpretation of the isotherms of multimodal ligands binding with DNA

The binding of ligands with DNA is a key moment in a whole range of cellular processes that provide not only the normal cell vital activity but also the development of some pathological processes. Depending on ligand type, structure of DNA adsorption centers, and physical–chemical conditions of the surrounding, the ligand may bind to DNA by several modes [1]. Particularly, adsorption isotherm of multimodal ligands binding to DNA in Scatchard’s coordinates has a concave shape with two brightly expressed linear areas in the region of small fillings. The analysis of such type of adsorption isotherm for determining of important binding parameters such as binding constant and number of adsorption centers (the part of DNA polymer with which one ligand molecule binds) presents difficulties. Practically in all cases, the analysis of such adsorption isotherm is carried out by linear parts of curves. Such analysis mode of experimental points is approximate method, since all registered of experimental points are roughly divided into two groups and they are treated by linear binding isotherm and therefore the binding parameters are determined. In the present work, the non-linear adsorption isotherm in Scatchard‘s coordinates is obtained which allowed, provided, the more precise treatment of all experimental points by unique curve which includes linear regions as well. Such mode of treatment of experimental points makes more precise the determination of not only binding constant and number of adsorption centers that correspond to the one ligand molecule binding, but also additional binding parameter – a proportion of adsorption centers of each binding to DNA type of multimodal ligand.     

Steps toward high specific activity labeling of biomolecules for therapeutic application: preparation of precursor [(188)Re(H(2)O)(3)(CO)(3)](+) and synthesis of tailor-made bifunctional ligand systems

Two kit preparations of the organometallic precursor [(188)Re(H(2)O)(3)(CO)(3)](+) in aqueous media are presented. Method A uses gaseous carbon monoxide and amine borane (BH(3).NH(3)) as the reducing agent. In method B CO(g) is replaced by K(2)[H(3)BCO(2)] that releases carbon monoxide during hydrolysis. Both procedures afford the desired precursor in yields >85% after 10 min at 60 degrees C. HPLC and TLC analyses revealed 7 +/- 3% of unreacted (188)ReO(4)(-) and <5% of colloidal (188)ReO(2). Solutions of up to 14 GBq/mL Re-188 have been successfully carbonylated with these two methods. The syntheses of two tailor-made bifunctional ligand systems for the precursor [(188)Re(H(2)O)(3)(CO)(3)](+) are presented. The tridentate chelates consist of a bis[imidazol-2-yl]methylamine or an iminodiacetic acid moiety, respectively. Both types of ligand systems have been prepared with alkyl spacers of different length and a pendent primary amino or carboxylic acid functionality, enabling the amidic linkage to biomolecules. The tridentate coordination of the ligands to the rhenium-tricarbonyl core could be elucidated on the macroscopic level by X-ray structure analyses and 1D and 2D NMR experiments of two representative model complexes. On the nca level, the ligands allow labeling yields >95% with [(188)Re(H(2)O)(3)(CO)(3)](+) under mild reaction conditions (PBS buffer, 60 degrees C, 60 min) at ligand concentrations between 5 x 10(-4) M and 5 x 10(-5) M. Thus, specific activities of 22-220 GBq pe micromol of ligand could be achieved. Incubation of the corresponding Re-188 complexes in human serum at 37 degrees C revealed stabilities between 80 +/- 4% and 45 +/- 10% at 24 h, respectively, and 63 +/- 3% and 34 +/- 3% at 48 h postincubation in human serum depending on the chelating system. Decomposition product was mainly (188)ReO(4)(-). The routine kit-preparation of the precursor [(188)Re(H(2)O)(3)(CO)(3)](+) in combination with tailor-made ligand systems enables the organometallic labeling of biomolecules with unprecedented high specific activities.     

Spatio-temporal dependence of the signaling response in immune-receptor trafficking networks regulated by cell density: a theoretical model

Cell signaling processes involve receptor trafficking through highly connected networks of interacting components. The binding of surface receptors to their specific ligands is a key factor for the control and triggering of signaling pathways. In most experimental systems, ligand concentration and cell density vary within a wide range of values. Dependence of the signal response on cell density is related with the extracellular volume available per cell. This dependence has previously been studied using non-spatial models which assume that signaling components are well mixed and uniformly distributed in a single compartment. In this paper, a mathematical model that shows the influence exerted by cell density on the spatio-temporal evolution of ligands, cell surface receptors, and intracellular signaling molecules is developed. To this end, partial differential equations were used to model ligand and receptor trafficking dynamics through the different domains of the whole system. This enabled us to analyze several interesting features involved with these systems, namely: a) how the perturbation caused by the signaling response propagates through the system; b) receptor internalization dynamics and how cell density affects the robustness of dose-response curves upon variation of the binding affinity; and c) that enhanced correlations between ligand input and system response are obtained under conditions that result in larger perturbations of the equilibrium ligand + surface receptor [Please see text] ligand - receptor complex. Finally, the results are compared with those obtained by considering that the above components are well mixed in a single compartment.     

Finite element solution of the steady-state Smoluchowski equation for rate constant calculations

This article describes the development and implementation of algorithms to study diffusion in biomolecular systems using continuum mechanics equations. Specifically, finite element methods have been developed to solve the steady-state Smoluchowski equation to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to mouse acetylcholinesterase. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different reaction criteria. The calculated rates were compared with experimental data and show very good agreement when the correct reaction criterion is used. Additionally, these finite element methods require significantly less computational resources than existing particle-based Brownian dynamics methods.     

Reconciliation of classical and reacted‐site probability approaches to allowance for ligand multivalence in binding studies

The objective of this investigation is to engender greater confidence in the validity of binding equations derived for multivalent ligands on the basis of reacted‐site probability theory. To that end, a demonstration of the theoretical interconnection between expressions derived by the classical stepwise equilibria and reacted‐site probability approaches for univalent ligands is followed by the use of the traditional stepwise procedure to derive binding equations for bivalent and trivalent ligands. As well as demonstrating the unwieldy nature of the classical binding equation for multivalent ligand systems, that exercise has allowed numerical simulation to be used to illustrate the equivalence of binding curves generated by the two approaches. The advantages of employing a redefined binding function for multivalent ligands are also confirmed by subjecting the simulated results to a published analytical procedure that has long been overlooked. Copyright © 2013 John Wiley & Sons, Ltd.     

EQUIL: Simulation and data analysis of binding reactions with arbitrary chemical models

We have developed an algorithm for simulation and analysis of arbitrary chemical systems in equilibrium, with emphasis on ligand binding reactions. The program EQUIL can treat reactions involving multiple ligands, multiple binding sites, ternary complex models, allosteric effectors, competitive and noncompetitive binding, conformational changes, cooperativity, and generally any scheme that can be represented as a set of chemical equations. EQUIL is based on a general thermodynamic model of chemical equilibria; it does not involve nonlinear transformation of experimental data, but it does require the user to define the model of interaction between ligands and receptors by writing down the appropriate chemical reactions. EQUIL contains features of particular importance to ligand binding experiments: variable binding capacities, nonspecific binding, and the ability to simultaneously analyze data from different types of experiments. Furthermore, the simulation feature of EQUIL allows the user to investigate the feasibility of experiments that could possibly distinguish between different reaction models. We illustrate the use of this program on personal computers to analyze and simulate simple and complicated interactions between ligands and receptors.     

New plot showing the existence of cooperativity, anticooperativity, and an arbitrary value of the neighbor-exclusion parameter in systems of ligand binding to linear biopolymers

A new way of plotting the isotherm for ligand binding to linear biopolymers is presented. In this plot the isotherm for noncooperative binding of ligands of length n (n-mers) becomes a straight line and the existence of cooperativity and anticooperativity between bound ligands is detected by appearance of opposite convexity of the curved isotherm. It is also usable for cases of fractional n values. Usefulness of the new plot in determining precise mechanisms of binding is shown using experimental data.     

Kinetics of binding the mRNA cap analogues to the translation initiation factor eIF4E under second-order reaction conditions

The kinetics of binding of five analogues of the 5'-mRNA cap, differing in size and electric charge, to the eukaryotic initiation factor eIF4E, at 20 degrees C, pH 7.2, and ionic strength of 150 mM, were measured, after mixing solutions of comparable concentrations of the reagents, in a stopped-flow spectrofluorimeter. The registered stopped-flow signals were fitted using an efficient software package, called Dyna Fit, based on a numerical solution of the kinetic rate equations for assumed reaction mechanisms. One-, two-, and three-step binding models were considered. The quality of fits for these models were compared using two statistical criteria: Akaike's Information Criterion and Bayesian Information Criterion. Based on resulting probabilities of the models, it was concluded that for all investigated ligands a one-step binding model has essentially no support in the experimental observations. Our conclusions were also analysed from the perspective of kinetic transients obtained for cap-eIF4E systems under the so called pseudo-first order reaction condition, which result in the linear correlation of the observed association rate constant with ligand concentration. The existence of such a linear correlation is usually considered as proof of a one-step binding mechanism. The kinetic and optical parameters, derived from fitting a two-step cap-binding model with the DynaFit, were used to simulate kinetic transients under pseudo-first order reaction conditions. It appeared that the observed association rate constants derived from these simulated transients are also linearly correlated with the ligand concentration. This indicated that these linear dependencies are not sufficient to conclude a one-step binding.     

Analysis of intracellular receptor/ligand sorting in endosomes

After binding to specific cell surface receptors, many extracellular ligand molecules are internalized via the process termed receptor-mediated endocytosis. Within the cell, in endosomes, a sorting process occurs: receptors and ligands are directed along various intracellular pathways. The extent of this intracellular separation of receptors from ligands has been shown experimentally to vary with receptor and ligand properties such as binding affinity and valency. In this paper, we propose and analyze a simple model mechanism for the sorting process based on binding and dissociation kinetics along with diffusive molecular transport. We show that the outcome of the sorting process can be directly linked to measurable parameters such as the intrinsic rate constants for the binding to, dissociation from, and crosslinking of receptors by ligands. We further show that this mechanism is able to account for the wide range of reported experimental observations. Manipulation of ligand and receptor properties guided by the results presented here may enable the outcome of the sorting process to be controlled.     

Kinetics of reaction zone formation with radial diffusion of ligands over a receptor-coated surface

Theoretical calculations are presented, describing the kinetics of reaction zone formation with radial diffusion of ligands over a receptor coated surface. Calculated concentration distributions of ligands diffusing radially over a receptor-coated surface are combined with different types of receptor-ligand reactions, taking place at the surface, in order to obtain theoretical relations between the initial concentration of ligand in the source, the diameter of the receptor-ligand reaction zone and reaction time. These relations are compared to experimental data, using bovine serum albumin (BSA) as immobilized receptor and anti-BSA antibodies as diffusing ligand. The theory predicts how the diffusion constant of the ligand and the detection level of the visualization method may be determined and how to discriminate between different kinetics of the receptor-ligand reaction. The practical use of the theory in experimental studies of receptor-ligand interaction is discussed.     

Mathematical model for the effects of adhesion and mechanics on cell migration speed.

Migration of mammalian blood and tissue cells over adhesive surfaces is apparently mediated by specific reversible reactions between cell membrane adhesion receptors and complementary ligands attached to the substratum. Although in a number of systems these receptors and ligand molecules have been isolated and identified, a theory capable of predicting the effects of their properties on cell migration behavior currently does not exist. We present a simple mathematical model for elucidating the dependence of cell speed on adhesion-receptor/ligand binding and cell mechanical properties. Our model can be applied to propose answers to questions such as: does an optimal adhesiveness exist for cell movement? How might changes in receptor and ligand density and/or affinity affect the rate of migration? Can cell rheological properties influence movement speed? This model incorporates cytoskeletal force generation, cell polarization, and dynamic adhesion as requirements for persistent cell movement. A critical feature is the proposed existence of an asymmetry in some cell adhesion-receptor property, correlated with cell polarity. We consider two major alternative mechanisms underlying this asymmetry: (a) a spatial distribution of adhesion-receptor number due to polarized endocytic trafficking and (b) a spatial variation in adhesion-receptor/ligand bond strength. Applying a viscoelastic-solid model for cell mechanics allows us to represent one-dimensional locomotion with a system of differential equations describing cell deformation and displacement along with adhesion-receptor dynamics. In this paper, we solve these equations under the simplifying assumption that receptor dynamics are at a quasi-steady state relative to cell locomotion. Thus, our results are strictly valid for sufficiently slow cell movement, as typically observed for tissue cells such as fibroblasts. Numerical examples relevant to experimental systems are provided. Our results predict how cell speed might vary with intracellular contractile force, cell rheology, receptor/ligand kinetics, and receptor/ligand number densities. A biphasic dependence is shown to be possible with respect to some of the system parameters, with position of the maxima essentially governed by a balance between transmitted contractile force and adhesiveness. We demonstrate that predictions for the two alternative asymmetry mechanisms can be distinguished and could be experimentally tested using cell populations possessing different adhesion-receptor numbers.     

Biological activities of polypyridyl-type ligands: implications for bioinorganic chemistry and light-activated metal complexes

Polypyridyl coordinating ligands are common in metal complexes used in medicinal inorganic chemistry. These ligands possess intrinsic cytotoxicity, but detailed data on this phenomenon are sparse, and cytotoxicity values vary widely and are often irreproducible. To provide new insights into the biological effects of bipyridyl-type ligands and structurally related metal-binding systems, reports of free ligand cytotoxicity were reviewed. The cytotoxicity of 25 derivatives of 2,2′-bipyridine and 1,10-phenanthroline demonstrates that there is no correlation between IC₅₀ values and ligand properties such as pKa, log D, polarizability volume, and electron density, as indicated by NMR shifts. As a result of these observations, as well as the various reported mechanisms of action of polypyridyl ligands, we offer the hypothesis that biological effects are governed by the availability of and affinity for specific metal ions within the experimental model.     

SpeedScreen: label-free liquid chromatography-mass spectrometry-based high-throughput screening for the discovery of orphan protein ligands

A high-throughput screening methodology tailored to the discovery of ligands for known and orphan proteins is presented. With this method, labeling of neither target protein nor screened compounds is required, as the ligands are affinity selected by incubation of the protein with mixtures of compounds in aqueous binding buffer. Unbound small-molecular-weight compounds are removed from the target protein:ligand complex by rapid size-exclusion chromatography in the 96-well format. The protein fraction is analyzed subsequently by liquid chromatography-mass spectrometry for detection and identification of the bound ligand. This screening method was validated with known protein:ligand model systems and optimized for selection of high-affinity binders in an industrial screening environment. All sample handling steps and the analytics are rapid, robust, and largely automated, adopting this technology to the needs of present high-throughput screening processes. This affinity-selection technology, termed SpeedScreen, is currently an integral part of our lead discovery process.     

Restricted diffusion of molecules in porous affinity chromatography adsorbents.

A restricted diffusion model is constructed and solved in order to study the permeability of large adsorbate molecules in the pores of affinity chromatography media, when the adsorbate molecules are adsorbed onto immobilized ligands. The combined effects of steric hindrance at the entrance to the pores and frictional resistance within the pores, as well as the effects of pore size distribution, pore connectivity of the adsorbent, molecular size of adsorbate and ligand, and the fractional saturation of adsorption sites (ligands), are considered. Affinity adsorbents with dilute and high ligand concentrations are examined, and the permeability of the adsorbate in porous networks of connectivity nT is studied by means of effective medium approximation (EMA) numerical solutions. As expected, the permeability of the adsorbate decreases as the size of the adsorbate and/or ligand molecule increases. The permeability also decreases when the fractional saturation of the ligands increases, as well as when the pore connectivity of the network decreases. The dependence of the permeability on the pore connectivity tends to be less marked in adsorbents with concentrated ligand than in porous media with dilute ligand concentration. The conditions are also presented for which the percolation threshold is attained in a number of different systems. The restricted diffusion model and results of this work may be of importance in studies involving the modeling, prediction of the dynamic behavior, design, and control of affinity chromatography (biospecific adsorption) systems employing porous adsorbents. The theoretical results may also have important implications in the selection of a ligand as well as in the selection and construction of an affinity porous matrix, so that the adsorbate of interest can be efficiently separated from a given solution. Furthermore, with appropriate modifications this restricted diffusion model may be used in studies involving the immobilization of ligands or enzymes in porous solids.     

Solvation energy density occlusion approximation for evaluation of desolvation penalties in biomolecular interactions

In calculations involving many displacements of an interacting pair of biomolecules, such as brownian dynamics, the docking of a substrate/ligand to an enzyme/receptor, or the screening of a large number of ligands as prospective inhibitors for a particular receptor site, there is a need for rapid evaluation of the desolvation penalties of the interacting pair. Although continuum electrostatic treatments with distinct dielectric constants for solute and solvent provide an account of the electrostatics of solvation and desolvation, it is necessary to re-solve the Poisson equation, at considerable computational cost, for each displacement of the interacting pair. We present a new method that uses a formulation of continuum electrostatic solvation in terms of the solvation energy density and approximates desolvation in terms of the occlusion of this density. We call it the SEDO approximation. It avoids the need to re-solve the Poisson equation, as desolvation is now estimated by an integral over the occluded volume. Test calculations are presented for some simple model systems and for some real systems that have previously been studied using the Poisson equation approach: MHC class I protein-peptide complexes and a congeneric series of human immunodeficiency virus type 1 (HIV-1) protease--ligand complexes. For most of the systems considered, the trends and magnitudes of the desolvation component of interaction energies obtained using the SEDO approximation are in reasonable correlation with those obtained by re-solving the Poisson equation. In most cases, the error introduced by the SEDO approximation is much less than that of the often-used test-charge approximation for the charge-charge components of intermolecular interactions. Proteins 2001;43:12-27.