Following ligand binding and ligand reactions in proteins via Raman crystallography

Raman crystallography permits the monitoring of chemical events in single-protein crystals in real time. Using a Raman microscope, it is possible to obtain protein Raman spectroscopic data of unprecedented quality and stability. The latter features allow us to obtain the Raman spectrum for small molecules soaking into crystals under normal (nonresonance) Raman conditions. Thus, via an approach utilizing Raman difference spectroscopy, we can quantitate the amount of ligand in the crystal, determine the chemistry of inhibitor-protein interactions, and follow chemical reactions in the active site on the time scale of minutes. While providing unique chemical insights, these data also provide an invaluable guide for determining the conditions for flash-freezing crystals for X-ray crystallographic analysis. In addition, the Raman difference spectra often contain contributions from protein modes due to protein conformational changes occurring upon ligand binding. These features allow us to probe events ranging from small cooperative conformational changes to massive and unexpected secondary structure changes in the crystal. An experimental advantage of Raman crystallography is that the data can be collected from crystals in situ, in sitting or hanging drops, under the conditions used to grow the crystals.     

Quantitative measurement of water diffusion lifetimes at a protein/DNA interface by NMR

Hydration site lifetimes of slowly diffusing water molecules at the protein/DNA interface of the vnd/NK-2 homeodomain DNA complex were determined using novel three-dimensional NMR techniques. The lifetimes were calculated using the ratios of ROE and NOE cross-relaxation rates between the water and the protein backbone and side chain amides. This calculation of the lifetimes is based on a model of the spectral density function of the water-protein interaction consisting of three timescales of motion: fast vibrational/rotational motion, diffusion into/out of the hydration site, and overall macromolecular tumbling. The lifetimes measured ranged from approximately 400 ps to more than 5 ns, and nearly all the slowly diffusing water molecules detected lie at the protein/DNA interface. A quantitative analysis of relayed water cross-relaxation indicated that even at very short mixing times, 5 ms for ROESY and 12 ms for NOESY, relay of magnetization can make a small but detectable contribution to the measured rates. The temperature dependences of the NOE rates were measured to help discriminate direct dipolar cross-relaxation from chemical exchange. Comparison with several X-ray structures of homeodomain/DNA complexes reveals a strong correspondence between water molecules in conserved locations and the slowly diffusing water molecules detected by NMR. A homology model based on the X-ray structures was created to visualize the conserved water molecules detected at the vnd/NK-2 homeodomain DNA interface. Two chains of water molecules are seen at the right and left sides of the major groove, adjacent to the third helix of the homeodomain. Two water-mediated hydrogen bond bridges spanning the protein/DNA interface are present in the model, one between the backbone of Phe8 and a DNA phosphate, and one between the side chain of Asn51 and a DNA phosphate. The hydrogen bond bridge between Asn51 and the DNA might be especially important since the DNA contact made by the invariant Asn51 residue, seen in all known homeodomain/DNA structures, is critical for binding affinity and specificity.     

Dissecting the contribution of diffusion and interactions to the mobility of nuclear proteins

Quantitative characterization of protein interactions under physiological conditions is vital for systems biology. Fluorescence photobleaching/activation experiments of GFP-tagged proteins are frequently used for this purpose, but robust analysis methods to extract physicochemical parameters from such data are lacking. Here, we implemented a reaction-diffusion model to determine the contributions of protein interaction and diffusion on fluorescence redistribution. The model was validated and applied to five chromatin-interacting proteins probed by photoactivation in living cells. We found that very transient interactions are common for chromatin proteins. Their observed mobility was limited by the amount of free protein available for diffusion but not by the short residence time of the bound proteins. Individual proteins thus locally scan chromatin for binding sites, rather than diffusing globally before rebinding at random nuclear positions. By taking the real cellular geometry and the inhomogeneous distribution of binding sites into account, our model provides a general framework to analyze the mobility of fluorescently tagged factors. Furthermore, it defines the experimental limitations of fluorescence perturbation experiments and highlights the need for complementary methods to measure transient biochemical interactions in living cells.     

Comparative study of molecular dynamics,diffusion, and permeability for ligands in biomembranes of different lipid composition

A comparative study of several model lipid bilayers of different composition, which included analysis of kinetic parameters of model lipid bilayers and permeability of bilayer membranes for small molecules, has been carried out. The conformity of results of numeric experiments to experimental data (structure of membrane lipid bilayers, lateral diffusion coefficients, and relative permeability of biomembranes for ligands) is discussed in the framework of a standard molecular dynamics protocol.     

Analysis of intracellular receptor/ligand sorting. Calculation of mean surface and bulk diffusion times within a sphere.

Cell surface receptors bind extracellular ligand molecules and transport those ligands into the cell by a process termed receptor-mediated endocytosis. Receptor and ligand molecules are sorted from one another after endocytosis, apparently within a structure consisting of intracellular vesicles and connected thin tubules. The experimental observation is that most free (unbound) ligand molecules are found in the lumen of the vesicles and receptors are located primarily within the tubules. Because equilibrium and geometric considerations do not explain this segregation, a kinetic scheme involving the passive diffusion of molecules from a vesicle into a tubule is investigated. Two possible sorting mechanisms are considered: first, that receptors are able to move into tubules more rapidly than ligand molecules due to an advantage in dimensionality and, second, that receptors diffusing into tubules are trapped there while ligands are not. Mean diffusion times for receptor and ligand movement into a tubule are calculated by solving Poisson's equation in two and three dimensions, respectively, on the surface of and within a sphere. Using estimated parameter values, we found that only the second scheme is able to account for the experimentally observed sorting. An estimate is obtained for the length of time a tubule and vesicle must be connected in order to remove a significant number of receptors into a tubule. The fraction of free ligand that is "mis-sorted" with the recycling receptor population and thus exocytosed is also determined.     

The size of the unstirred layer as a function of the solute diffusion coefficient.

By monitoring the concentration distribution of several solutes that are diffusing at the same time under given mixing conditions, it was established that the unstirred layer (USL) has no clearly defined boundary. For the cases of solute permeation and water movement across planar bilayer lipid membranes, respectively, experiments carried out with double-barreled microelectrodes have shown that the thickness of the USL depends on which species is diffusing. Small molecules with a larger diffusion coefficient encounter an apparently thicker USL than larger molecules with a smaller diffusion coefficient. The ratio of the USL thicknesses of two different substances is equal to the third root of the ratio of the respective diffusion coefficients. This experimental finding is in good agreement with theoretical predictions from the theory of physicochemical hydrodynamics.     

A biological interpretation of transient anomalous subdiffusion. I. Qualitative model

Anomalous subdiffusion has been reported for two-dimensional diffusion in the plasma membrane and three-dimensional diffusion in the nucleus and cytoplasm. If a particle diffuses in a suitable infinite hierarchy of binding sites, diffusion is well known to be anomalous at all times. But if the hierarchy is finite, diffusion is anomalous at short times and normal at long times. For a prescribed set of binding sites, Monte Carlo calculations yield the anomalous diffusion exponent and the average time over which diffusion is anomalous. If even a single binding site is present, there is a very short, almost artifactual, period of anomalous subdiffusion, but a hierarchy of binding sites extends the anomalous regime considerably. As is well known, an essential requirement for anomalous subdiffusion due to binding is that the diffusing particle cannot be in thermal equilibrium with the binding sites; an equilibrated particle diffuses normally at all times. Anomalous subdiffusion due to barriers, however, still occurs at thermal equilibrium, and anomalous subdiffusion due to a combination of binding sites and barriers is reduced but not eliminated on equilibration. This physical model is translated directly into a plausible biological model testable by single-particle tracking.     

Anomalous diffusion of fluorescent probes inside living cell nuclei investigated by spatially-resolved fluorescence correlation spectroscopy

We have investigated spatial variations of the diffusion behavior of the green fluorescent protein mutant EGFP (F64L/S65T) and of the EGFP-beta-galactosidase fusion protein in living cells with fluorescence correlation spectroscopy. Our fluorescence correlation spectroscopy device, in connection with a precision x-y translation stage, provides submicron spatial resolution and a detection volume smaller than a femtoliter. The fluorescence fluctuations in cell lines expressing EGFP are caused by molecular diffusion as well as a possible internal and a pH-dependent external protonation process of the EGFP chromophore. The latter processes result in two apparent nonfluorescent states that have to be taken into account when evaluating the fluorescence correlation spectroscopy data. The diffusional contribution deviates from ideal behavior and depends on the position in the cell. The fluorescence correlation spectroscopy data can either be evaluated as a two component model with one fraction of the molecules undergoing free Brownian motion with a diffusion coefficient approximately five times smaller than in aqueous solution, and another fraction diffusing one or two orders of magnitude slower. This latter component is especially noticeable in the nuclei. Alternatively, we can fit the data to an anomalous diffusion model where the time dependence of the diffusion serves as a measure for the degree of obstruction, which is large especially in nuclei. Possible mechanisms for this long tail behavior include corralling, immobile obstacles, and binding with a broad distribution of binding affinities. The results are consistent with recent numerical models of the chromosome territory structure in the cell nucleus.     

Transport effects on the kinetics of protein-surface binding.

A detailed model is presented for protein binding to active surfaces, with application to the binding of avidin molecules to a biotin-functionalized fiber optic sensor in experiments reported by S. Zhao and W. M. Reichert (American Chemical Society Symposium Series 493, 1992). Kinetic data for binding in solution are used to assign an intrinsic catalytic rate coefficient k to the biotin-avidin pair, deconvoluted from transport and electrostatic factors via application of coagulation theory. This intrinsic chemical constant is built into a reaction-diffusion analysis of surface binding where activity is restricted to localized sites (representing immobilized biotin molecules). The analysis leads to an effective catalytic rate coefficient keff characterizing the active surface. Thereafter, solution of the transport problem describing absorption of avidin molecules by the macroscopic sensor surface leads to predictions of the avidin flux, which are found to be in good agreement with the experimental data. The analysis suggests the following conclusions. 1) Translational diffusion limitations are negligible for avidin-biotin binding in solution owing to the small (kinetically limiting) value k = 0.00045 m/s. 2) The sparse distribution of biotin molecules and the presence of a repulsive hydration force produce an effective surface-average catalytic rate coefficient keff of order 10(-7) m/s, much smaller than k. 3) Avidin binding to the fiber optic sensor occurs in an intermediate regime where the rate is influenced by both kinetics and diffusion.     

Resolving Cytosolic Diffusive States in Bacteria by Single-Molecule Tracking

The trajectory of a single protein in the cytosol of a living cell contains information about its molecular interactions in its native environment. However, it has remained challenging to accurately resolve and characterize the diffusive states that can manifest in the cytosol using analytical approaches based on simplifying assumptions. Here, we show that multiple intracellular diffusive states can be successfully resolved if sufficient single-molecule trajectory information is available to generate well-sampled distributions of experimental measurements and if experimental biases are taken into account during data analysis. To address the inherent experimental biases in camera-based and MINFLUX-based single-molecule tracking, we use an empirical data analysis framework based on Monte Carlo simulations of confined Brownian motion. This framework is general and adaptable to arbitrary cell geometries and data acquisition parameters employed in two-dimensional or three-dimensional single-molecule tracking. We show that, in addition to determining the diffusion coefficients and populations of prevalent diffusive states, the timescales of diffusive state switching can be determined by stepwise increasing the time window of averaging over subsequent single-molecule displacements. Time-averaged diffusion analysis of single-molecule tracking data may thus provide quantitative insights into binding and unbinding reactions among rapidly diffusing molecules that are integral for cellular functions.     

FCS cell surface measurements--photophysical limitations and consequences on molecular ensembles with heterogenic mobilities.

BACKGROUND: Fluorescence Correlation Spectroscopy is a powerful method to analyze densities and diffusive behavior of molecules in membranes, but effects of photodegradation can easily be overlooked. METHOD: Based on experimental photophysical parameters, calculations were performed to analyze the consequences of photobleaching in fluorescence correlation spectroscopy (FCS) cell surface experiments, covering a range of standard measurement conditions. RESULTS: Cumulative effects of photobleaching can be prominent, although an absolute majority of the fluorescent molecules would pass the laser excitation beam without being photo-bleached. Given a distribution of molecules on a cell surface with different diffusive properties, the fraction of molecules that is actually analyzed depends strongly on the excitation intensities and measurement times, as well as on the size of the reservoir of freely diffusing molecules. Both the slower and the faster diffusing molecules can be disfavored. CONCLUSIONS: Apart from quantifying photobleaching effects, the calculations suggest that the effects can be used to extract additional information, for instance about the size of the reservoirs of free diffusion. By certain choices of measurement conditions, it may be possible to more specifically analyze certain species within a population, based on their different diffusive properties, different areas of free diffusion, or different kinetics of possible transient binding.     

An automated continuous-flow dynamic dialysis technique for investigating protein-ligand binding

An automated, continuous-flow dynamic dialysis technique has been developed to investigate protein-ligand binding. The method depends on a comparison of the diffusion of the low molecular mass ligand, in the presence and absence of protein, through a semipermeable membrane. The ligand passes from the sample compartment of a dialysis cell into the sink compartment through which a constant flow of eluting buffer is maintained. Digitized spectrophotometric determinations of the ligand concentration in the eluting buffer at successive, equally spaced time intervals, punched onto paper tape, provide the primary data (normally about 1000 data points). A mathematical treatment of the data based on a model of the diffusion system, whereby the protein-ligand binding isotherm may be evaluated, is discussed. The validity of the method is demonstrated from studies of the binding of phenol red to bovine serum albumin (BSA) at 15, 20, and 25°C. The method yields a large number of points on the binding isotherm (usually several hundred) which, in terms of a Scatchard model, provide values for the number of binding sites on the BSA molecule and binding constants for the phenol red-BSA interaction. The results obtained are consistent with values reported in the chemical literature but which are based on much scantier data.     

Complexes of porcine odorant binding protein with odorant molecules belonging to different chemical classes

Porcine odorant binding protein (pOBP) is a monomer of 157 amino acid residues, purified in abundance from pig nasal mucosa. In contrast to the observation on lipocalins as retinol binding protein (RBP), major urinary protein (MUP) or bovine odorant binding protein (bOBP), no naturally occurring ligand was found in the beta-barrel cavity of pOBP. Porcine OBP was therefore chosen as a simple model for structure/function studies with odorant molecules. In competition experiments with tritiated pyrazine, the affinity of pOBP towards several odorant molecules belonging to different chemical classes has been found to be of the micromolar order, with a 1:1 stoichiometry. The X-ray structures of pOBP complexed to these molecules were determined at resolution between 2.15 and 1.4 A. As expected, the electron density of the odorant molecules was observed into the hydrophobic beta-barrel of the lipocalin. Inside this cavity, very few specific interactions were established between the odorant molecule and the amino acid side-chains, which did not undergo significant conformational change. The high B-factors observed for the odorant molecules as well as the existence of alternative conformations reveal a non-specific mode of binding of the odorant molecules in the cavity.     

Calculation of diffusion-limited kinetics for the reactions in collision coupling and receptor cross-linking

Both enzyme (e.g., G-protein) activation via a collision coupling model and the formation of cross-linked receptors by a multivalent ligand involve reactions between two molecules diffusing in the plasma membrane. The diffusion of these molecules is thought to play a critical role in these two early signal transduction events. In reduced dimensions, however, diffusion is not an effective mixing mechanism; consequently, zones in which the concentration of particular molecules (e.g., enzymes, receptors) becomes depleted or enriched may form. To examine the formation of these depletion/ accumulation zones and their effect on reaction rates and ultimately the cellular response, Monte Carlo techniques are used to simulate the reaction and diffusion of molecules in the plasma membrane. The effective reaction rate at steady state is determined in terms of the physical properties of the tissue and ligand for both enzyme activation via collision coupling and the generation of cross-linked receptors. The diffusion-limited reaction rate constant is shown to scale with the mean square displacement of a receptor-ligand complex. The rate constants determined in the simulation are compared with other theoretical predictions as well as experimental data.     

Construction and validation of an atomic model for bacterial TSPO from electron microscopy density,evolutionary constraints,and biochemical and biophysical data

The 18 kDa protein TSPO is a highly conserved transmembrane protein found in bacteria, yeast, animals and plants. TSPO is involved in a wide range of physiological functions, among which the transport of several molecules. The atomic structure of monomeric ligand-bound mouse TSPO in detergent has been published recently. A previously published low-resolution structure of Rhodobacter sphaeroides TSPO, obtained from tubular crystals with lipids and observed in cryo-electron microscopy, revealed an oligomeric structure without any ligand. We analyze this electron microscopy density in view of available biochemical and biophysical data, building a matching atomic model for the monomer and then the entire crystal. We compare its intra- and inter-molecular contacts with those predicted by amino acid covariation in TSPO proteins from evolutionary sequence analysis. The arrangement of the five transmembrane helices in a monomer of our model is different from that observed for the mouse TSPO. We analyze possible ligand binding sites for protoporphyrin, for the high-affinity ligand PK 11195, and for cholesterol in TSPO monomers and/or oligomers, and we discuss possible functional implications.     

Time-resolved protein crystallography.

Advances in synchrotron radiation technology have allowed exposure times from protein crystals of the order of milliseconds to be used routinely, and in exceptional circumstances exposure times of 100 ps have been obtained. However, many data sets take seconds to record because of the slow time scale of film change or crystal reorientation or translation when more than one exposure is required. This problem has been addressed by Amemiya et al. (1989). There has been considerable progress in methods to initiate reactions in protein crystals, especially the development of photolabile caged compounds but also temperature jump, pH jump, and diffusion. Although flash lamps deliver pulses of 100 mJ/ms, often several pulses are required to release sufficient product, and reaction initiation can take several seconds. Laser illumination can provide more powerful input, but the laser must be accommodated within the restricted space at the synchrotron station. The requirement to maintain synchrony among the molecules in the crystal lattice as the reaction proceeds and to ensure that the lifetime of intermediates is longer than data collection rates emphasizes the need for chemical characterization of the reaction under study. As Ringe advocated in the studies with chymotrypsin, it may be more profitable to devise conditions under which certain intermediates along the reaction pathway accumulate in the crystal and to record these in a series of discrete steps rather than continuous monitoring of the reaction. The Laue method is limited to those proteins that give well-ordered crystals and problems of transient disorder on initiation of reaction and problems of radiation damage need to be overcome or avoided by suitable experimental protocols.(ABSTRACT TRUNCATED AT 250 WORDS)     

Quantifying anisotropic solute transport in protein crystals using 3-D laser scanning confocal microscopy visualization

The diffusion of a solute, fluorescein into lysozyme protein crystals has been studied by confocal laser scanning microscopy (CLSM). Confocal laser scanning microscopy makes it possible to non-invasively obtain high-resolution three-dimensional (3-D) images of spatial distribution of fluorescein in lysozyme crystals at various time steps. Confocal laser scanning microscopy gives the fluorescence intensity profiles across horizontal planes at several depths of the crystal representing the concentration profiles during diffusion into the crystal. These intensity profiles were fitted with an anisotropic model to determine the diffusivity tensor. Effective diffusion coefficients obtained range from 6.2 x 10(-15) to 120 x 10(-15) m2/s depending on the lysozyme crystal morphology. The diffusion process is found to be anisotropic, and the level of anisotropy depends on the crystal morphology. The packing of the protein molecules in the crystal seems to be the major factor that determines the anisotropy.     

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may “select” protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational‐selection and induced‐change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational‐selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced‐change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on‐ and off‐rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single‐molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational‐selection and induced‐change processes in both binding and unbinding direction.