Sensitivity of flow simulations in a cerebral aneurysm.

The flow in a cerebral aneurysm on the Anterior Communicating Artery is simulated based on patient medical imaging. The sensitivity of the flow pattern is assessed with respect to variation of outflow split and geometric parameters. A comparison is made between the results from unsteady pulsatile computations and steady-state computations. The results show that for this case there is no significant sensitivity to pulsatile flow-rate variations. Steady-state and pulsatile simulations both predict a strong sensitivity to outflow split and geometric smoothing. It is therefore proposed to analyse the stability of aneurysmal flow with an extensive sensitivity study based on steady-state computations and use this sensitivity for rupture prediction.

Specificity of an immunochromatographic test for anthrax.

OBJECTIVE: To evaluate the specificity of an immunochromatographic test (ICT) for anthrax in cattle. DESIGN: A comparison of an ICT with blood smear and culture in uninfected cattle. PROCEDURE: Two hundred and forty blood samples were collected from dead cattle at two knackeries within Victoria and tested on-site with an ICT for the detection of protective antigen (PA) of Bacillus anthracis. Blood smears were prepared on-site and blood samples transported to the laboratory for aerobic and anaerobic culture. The results of the ICT were compared with blood smear and culture. Animals were regarded as not infected with B anthracis if the organism was not detected in a stained blood smear or on culture. Ten healthy yearling cattle were vaccinated with live spore anthrax vaccine and blood samples collected on days 0 to 7 and day 15 were tested in the ICT for the presence of PA. RESULTS: All blood samples from the 240 knackery cattle were ICT, smear and culture negative. All blood samples from the 10 vaccinated cattle were ICT negative. CONCLUSION: The ICT is a test with high specificity in cattle (98.5 to 100%; 95% CI) and recent vaccination of cattle does not give rise to positive reactions.

SpikeCell: a deterministic spiking neuron.

We present a model of spiking neuron that emulates the output of the usual static neurons with sigmoidal activation functions. It allows for hardware implementations of standard feedforward networks, trained off-line with any classical learning algorithm (i.e. back-propagation and its variants). The model is validated on hand-written digits recognition, and image classification tasks. A digital architecture is proposed and evaluated. The area needed for implementing the spiking neuron on a chip is 10 times smaller than that for the corresponding static neuron. The accuracy of the network's output increases with time, and reaches that of the emulated static neural network after an adequate integration period. Single errors in the spike trains, or interruption of the relaxation process, due for example to irradiation in harsh environments, are harmless.

4',6-diamidino-2-phenylindole (DAPI) interacts with rare structures of GC polymers.

The binding of 4',6-diamidino-2-phenylindole (DAPI) to double-stranded GC polymers either in the alternating or in homopolymer sequence was investigated using fluorescence techniques. We employed fluctuation correlation spectroscopy, which measures the diffusion coefficient of fluorescent particles, to demonstrate that the fluorescence was originating from relatively slowly diffusing entities. These entities display a very large heterogeneity of diffusing coefficients, indicating that molecular aggregation is extensive in our samples. We used frequency domain fluorometry to characterize the fluorescence lifetime of the species, while varying the GC polymer-dye coverage systematically. At very low coverage we observed a relatively bright fluorescent component with a lifetime value of approximately 4 ns. The stoichiometry of binding of this bright species was such that it can only arise from rare molecular structures, either unusual loops or large molecular aggregates. The amount and characteristics of this bright fluorescent component were different between the homo and the alternating polymer, indicating that the difference in sequence of the two polymers is responsible for the different aggregates which are then detected in the fluorescence experiment. At large GC polymer coverage we observed a relatively wide distribution of fluorescent species with short lifetime values, in the range between 0.12 and 0.2 ns. Given the stoichiometry of binding of this fluorescent component, we concluded that it could arise either from intercalative and/or non-specific binding to the DNA double-stranded molecules. We comment on the origin of the rare but brightly fluorescent binding sites and discuss the potential to detect such unusual DNA structures.

Solution and interface aggregation states of Crotalus atrox venom phospholipase A2 by two-photon excitation fluorescence correlation spectroscopy.

The dimeric Crotalus atrox venom PLA2 is part of the secreted phospholipase A2 (PLA2) enzyme family that interacts at the lipid-solution interface to hydrolyze the sn-2 acyl ester bond of phospholipids. We have employed fluorescence correlation spectroscopy (FCS) to study the monomer-dimer equilibrium of the C. atrox venom PLA2 in solution, in the presence of urea, and in the presence of monomeric and micellar n-dodecylphosphocholine (C12-PN), a phosphatidylcholine analogue. Dilution experiments show that PLA2 is an extremely tight dimer, Kd < or = 0.01 nM, in solution. Urea was introduced to weaken the subunit's association, and an estimate for the PLA(2) dimer dissociation constant in buffer was obtained by linear extrapolation. The derived dissociation constant was at least several orders of magnitude greater than that suggested from the dilution experiments, indicating a complex interaction between urea and the PLA2 dimer. FCS data indicate that the PLA2 dimer begins to dissociate at 10 mM C12-PN in 10 mM Ca2+ and at 5 mM C12-PN in 1 mM EDTA. The PLA2 tryptophan fluorescence displayed spectral shifts and intensity changes upon interacting with C12-PN. On the basis of the FCS and tryptophan fluorescence results, we postulate an intermediate state where the two monomers are in loose interaction within a protein-lipid comicelle. As the concentration of C12-PN was increased, complete dissociation of the dimer was observed, inferred from the doubling of the particle number, and the average diffusion constant decreased to approximately 60 microm2/s, consistent with PLA2 associated with a C12-PN micelle. The presence of Ca2+ makes the comicelle intermediate more stable, retarding the separation of the monomers in the micellar suspension. Our data clearly indicate that PLA2, though a strong dimer in the absence of lipids, is dissociated by micellar C12-PN and supports the monomer hypothesis for PLA2 action.

Fluorescence fluctuation analysis for the study of interactions between oligonucleotides and polycationic polymers.

The interactions between a cationic polymer, poly(2-dimethylamino)ethyl methacrylate (pDMAEMA), and negatively charged rhodamine-labeled 25-mer phosphodiester oligonucleotides (Rh-ONs) were studied by fluorescence fluctuation spectroscopy and other techniques. The composition of the pDMAEMA/Rh-ON complexes was investigated as a function of the charge ratio (+/-) by increasing the pDMAEMA concentration and keeping the Rh-ON concentration constant. We applied two different methods for analyzing the fluorescence fluctuation profiles of the pDMAEMA/Rh-ON complexes, which depended on their composition. First, we analyzed the data with the photon counting histogram (PCH) technique, which determines the molecular brightness and the concentration of fluorophores (Chen et al, 1999). A particular challenge for the data analysis is the occurrence of sudden fluorescence bursts in the fluorescence fluctuation profiles, which are linked to the appearance of multimolecular complexes (i. e. when several Rh-ONs were present in one complex). A quantitative interpretation of the analysis for the complexes remains challenging and is connected to the rarity of the fluorescent bursts, which do not provide sufficient data statistics. To specifically address the problem of the fluorescent bursts we employed a method described by Van Craenenbroeck et al. (1999). This method, applicable only when data were integrated over much longer time bins, allowed us to estimate the number of fluorescence bursts which could be considered as a relative measure of the amount of multimolecular complexes present. When monomolecular complexes were formed, i. e. at high values of the charge ratio, highly intense fluorescence peaks were not present and the interpretation of the PCH analysis was more straightforward. The molecular brightness of the species (epsilon), as revealed from PCH analysis, was greater than epsilon for the free Rh-ONs, indicating that the Rh-ONs were attached to pDMAEMA chains.

Molecular heterogeneity of O-acetylserine sulfhydrylase by two-photon excited fluorescence fluctuation spectroscopy.

O-acetylserine sulfhydrylase, a homo-dimeric enzyme from Salmonella typhimurium, covalently binds one pyridoxal 5'-phosphate molecule per subunit as a fluorescent coenzyme. Different tautomers of the Schiff base between the coenzyme and lysine 41 generate structured absorption and fluorescence spectra upon one-photon excitation. We investigated the protein population heterogeneity by fluorescence correlation spectroscopy and lifetime techniques upon two-photon excitation. We sampled the fluorescence intensity from a small number of molecules (approximately 10) and analyzed the distribution of photon counts to separately determine the number and the fluorescence brightness of excited protein molecules. The changes in the average number of molecules and in the fluorescence brightness with the excitation wavelength indicate the presence of at least two fluorescent species, with two-photon excitation maxima at 660 and 800 nm. These species have been identified as the enolimine and ketoenamine tautomers of the protein-coenzyme internal aldimine. Their relative abundance is estimated to be 4:1, whereas the ratio of their two-photon cross sections is reversed with respect to the single-photon excitation case. Consistent results are obtained from the measurement of the lifetime decays, which are sensitive to the excited-state heterogeneity. At least two components were detected, with lifetimes of approximately 2.5 and 0.5 ns. The lifetimes are very close to the values measured in bulk solutions upon one-photon excitation and attributed to the ketoenamine tautomer and to a dipolar species formed upon proton dissociation in the excited state.

The human estrogen receptor alpha dimer binds a single SRC-1 coactivator molecule with an affinity dictated by agonist structure.

Nuclear receptors act as ligand-inducible transcription factors. Agonist binding leads to interaction with coactivator proteins, and to the assembly of the general transcription machinery. In addition to structural information, a thorough understanding of transcriptional activation by the nuclear receptors requires the characterization of the thermodynamic parameters governing these protein/protein interactions. In this study we have quantitatively characterized the interactions of full-length baculovirus expressed human estrogen receptor alpha (ERalpha), as well as ERalpha hormone binding domain (ERHBD) with a fragment of the coactivator protein SRC-1 (amino acid residues 570 to 780). Fluorescence anisotropy and fluorescence correlation spectroscopy of fluorescently labeled SRC-1(570-780) demonstrate unambiguously that the stoichiometry of the SRC-1/ERalpha/estradiol complex is one coactivator molecule per ERalpha dimer. The affinity of the estradiol or estriol bound ERalpha/SRC-1 complexes was found to be significantly higher than that observed in the presence of estrone. No binding was observed in the absence of ligand or in the presence of antagonists. Distinct anisotropy values for the ERalpha-SRC-1 complexes with different agonists suggest distinct conformations of the complexes depending upon agonist structure.

Surveillance for aquabirnavirus in fish hatcheries in Victoria.

OBJECTIVE: To determine whether or not aquabirnavirus is present in Victorian fish hatcheries. PROCEDURE: Milt and ovarian fluids were collected from brood stock at 12 hatcheries and cultured in two sensitive cell lines for the presence of viruses. RESULTS: No cytopathic effect was detected indicating the absence of virus. CONCLUSION: There is no evidence of infection with aquabirnavirus or Epizootic Haematopoietic Necrosis virus in Victorian fish hatcheries.

Probing ligand protein binding equilibria with fluorescence fluctuation spectroscopy.

We examine the binding of fluorescent ligands to proteins by analyzing the fluctuation amplitude g(0) of fluorescence fluctuation experiments. The normalized variance g(0) depends on the molecular brightness and the concentration of each species in the sample. Thus a single g(0) measurement is not sufficient to resolve individual species. Titration of the ligand with protein establishes the link between molecular brightness and concentration by fitting g(0) to a binding model and allows the separation of species. We first apply g(0) analysis to binary dye mixtures with brightness ratios of 2 and 4 to demonstrate the feasibility of this technique. Next we consider the influence of binding on the fluctuation amplitude g(0). The dissociation coefficient, the molecular brightness ratio, and the stochiometry of binding strongly influence the fluctuation amplitude. We show that proteins with a single binding site can be clearly differentiated from proteins with two independent binding sites. The binding of fluorescein-labeled digoxigenin to a high-affinity anti-digoxin antibody was studied experimentally. A global analysis of the fluctuation amplitude and the fluorescence intensity not only recovered the dissociation coefficient and the number of binding sites, but also revealed the molecular heterogeneity of the hapten-antibody complex. Two species were used to model the molecular heterogeneity. We confirmed the molecular heterogeneity independently by fluorescence lifetime experiments, which gave fractional populations and molecular brightness values that were virtually identical to those of the g(0) analysis. The identification and characterization of molecular heterogeneity have far-reaching consequences for many biomolecular systems. We point out the important role fluctuation experiments may have in this area of research.

Resolving heterogeneity on the single molecular level with the photon-counting histogram.

The diffusion of fluorescent particles through a small, illuminated observation volume gives rise to intensity fluctuations caused by particle number fluctuations in the open observation volume and the inhomogeneous excitation-beam profile. The intensity distribution of these fluorescence fluctuations is experimentally captured by the photon-counting histogram (PCH). We recently introduced the theory of the PCH for diffusing particles (Chen et al., Biophys. J., 77:553-567), where we showed that we can uniquely describe the distribution of photon counts with only two parameters for each species: the molecular brightness of the particle and the average number of particles within the observation volume. The PCH is sensitive to the molecular brightness and thus offers the possibility to separate a mixture of fluorescent species into its constituents, based on a difference in their molecular brightness alone. This analysis is complementary to the autocorrelation function, traditionally used in fluorescence fluctuation spectroscopy, which separates a mixture of species by a difference in their diffusion coefficient. The PCH of each individual species is convoluted successively to yield the PCH of the mixture. Successful resolution of the histogram into its components is largely a matter of the signal statistics. Here, we discuss the case of two species in detail and show that a concentration for each species exists, where the signal statistics is optimal. We also discuss the influence of the absolute molecular brightness and the brightness contrast between two species on the resolvability of two species. A binary dye mixture serves as a model system to demonstrate that the molecular brightness and the concentration of each species can be resolved experimentally from a single or from several histograms. We extend our study to biomolecules, where we label proteins with a fluorescent dye and show that a brightness ratio of two can be resolved. The ability to resolve a brightness ratio of two is very important for biological applications.

Fluorescence fluctuation spectroscopy.

The analysis of the intensity fluctuation of a fluorescence signal from a relatively small volume and from a few molecules contains information about the distribution of different species present in the solution and about kinetic parameters of the system. The same information is generally averaged out when the fluorescence experiment is performed in a much larger volume, typically a cuvette experiment. The fundamental reason for this difference is that the fluctuations of the fluorescence signal from a few molecules directly reflect the molecular nature of the matter. Only recently, with the advent of confocal microscopy and two-photon excitation, it has become practical to achieve small excitation volumes in which only a few fluorescent molecules are present. We introduce the concept of fluctuation spectroscopy and highlight some of the technical aspects. We discuss different analysis methods used in fluctuation spectroscopy and evaluate their use for studying protein-protein interactions.

The photon counting histogram in fluorescence fluctuation spectroscopy.

Fluorescence correlation spectroscopy (FCS) is generally used to obtain information about the number of fluorescent particles in a small volume and the diffusion coefficient from the autocorrelation function of the fluorescence signal. Here we demonstrate that photon counting histogram (PCH) analysis constitutes a novel tool for extracting quantities from fluorescence fluctuation data, i.e., the measured photon counts per molecule and the average number of molecules within the observation volume. The photon counting histogram of fluorescence fluctuation experiments, in which few molecules are present in the excitation volume, exhibits a super-Poissonian behavior. The additional broadening of the PCH compared to a Poisson distribution is due to fluorescence intensity fluctuations. For diffusing particles these intensity fluctuations are caused by an inhomogeneous excitation profile and the fluctuations in the number of particles in the observation volume. The quantitative relationship between the detected photon counts and the fluorescence intensity reaching the detector is given by Mandel's formula. Based on this equation and considering the fluorescence intensity distribution in the two-photon excitation volume, a theoretical expression for the PCH as a function of the number of molecules in the excitation volume is derived. For a single molecular species two parameters are sufficient to characterize the histogram completely, namely the average number of molecules within the observation volume and the detected photon counts per molecule per sampling time epsilon. The PCH for multiple molecular species, on the other hand, is generated by successively convoluting the photon counting distribution of each species with the others. The influence of the excitation profile upon the photon counting statistics for two relevant point spread functions (PSFs), the three-dimensional Gaussian PSF conventionally employed in confocal detection and the square of the Gaussian-Lorentzian PSF for two photon excitation, is explicitly treated. Measured photon counting distributions obtained with a two-photon excitation source agree, within experimental error with the theoretical PCHs calculated for the square of a Gaussian-Lorentzian beam profile. We demonstrate and discuss the influence of the average number of particles within the observation volume and the detected photon counts per molecule per sampling interval upon the super-Poissonian character of the photon counting distribution.

Connection between the taxonomic substates and protonation of histidines 64 and 97 in carbonmonoxy myoglobin.

Infrared spectra of heme-bound CO in sperm whale carbonmonoxy myoglobin and two mutants (H64L and H97F) were studied in the pH range from 4.2 to 9.5. Comparison of the native protein with the mutants shows that the observed pH effects can be traced to protonations of two histidine residues, H64 and H97, near the active site. Their imidazole sidechains experience simple, uncoupled Henderson-Hasselbalch type protonations, giving rise to four different protonation states. Because two of the protonation states are linked by a pH-independent equilibrium, the overall pH dependence of the spectra is described by a linear combination of three independent components. Global analysis, based on singular value decomposition and matrix least-squares algorithms enabled us to extract the pK values of the two histidines and the three basis spectra of the protonating species. The basis spectra were decomposed into the taxonomic substates A(0), A(1), and A(3), previously introduced in a heuristic way to analyze CO stretch spectra in heme proteins at fixed pH (see for instance, Biophys. J. 71:1563-1573). Moreover, an additional, weakly populated substate, called A(x), was identified. Protonation of H97 gives rise to a blue shift of the individual infrared lines by about 2 cm(-1), so that the A substates actually appear in pairs, such as A(0) and A(0)(+). The blue shift can be explained by reduced backbonding from the heme iron to the CO. Protonation of the distal histidine, H64, leads to a change of the infrared absorption from the A(1) or A(3) substate lines to A(0). This behavior can be explained by a conformational change upon protonation that moves the imidazole sidechain of H64 away from the CO into the high-dielectric solvent environment, which avoids the energetically unfavorable situation of an uncompensated electric charge in the apolar, low-dielectric protein interior. Our results suggest that protonation reactions serve as an important mechanism to create taxonomic substates in proteins.

Electron transfer and protein dynamics in the photosynthetic reaction center.

We have measured the kinetics of electron transfer (ET) from the primary quinone (Q(A)) to the special pair (P) of the reaction center (RC) complex from Rhodobacter sphaeroides as a function of temperature (5-300 K), illumination protocol (cooled in the dark and under illumination from 110, 160, 180, and 280 K), and warming rate (1.3 and 13 mK/s). The nonexponential kinetics are interpreted with a quantum-mechanical ET model (Fermi's golden rule and the spin-boson model), in which heterogeneity of the protein ensemble, relaxations, and fluctuations are cast into a single coordinate that relaxes monotonically and is sensitive to all types of relaxations caused by ET. Our analysis shows that the structural changes that occur in response to ET decrease the free energy gap between donor and acceptor states by 120 meV and decrease the electronic coupling between donor and acceptor states from 2.7 x 10(-4) cm(-1) to 1.8 x 10(-4) cm(-1). At cryogenic temperatures, conformational changes can be slowed or completely arrested, allowing us to monitor relaxations on the annealing time scale (approximately 10(3)-10(4) s) as well as the time scale of ET (approximately 100 ms). The relaxations occur within four broad tiers of conformational substates with average apparent Arrhenius activation enthalpies of 17, 50, 78, and 110 kJ/mol and preexponential factors of 10(13), 10(15), 10(21), and 10(25) s(-1), respectively. The parameterization provides a prediction of the time course of relaxations at all temperatures. At 300 K, relaxations are expected to occur from 1 ps to 1 ms, whereas at lower temperatures, even broader distributions of relaxation times are expected. The weak dependence of the ET rate on both temperature and protein conformation, together with the possibility of modeling heterogeneity and dynamics with a single conformational coordinate, make RC a useful model system for probing the dynamics of conformational changes in proteins.

FTIR study of conformational substates in the CO adduct of cytochrome c oxidase from Rhodobacter sphaeroides.

Fourier transform infrared (FTIR) spectroscopy of cytochrome c oxidase from Rhodobacter sphaeroides reveals multiple CO stretch bands that are associated with different conformational substates of the enzyme. Here we report the temperature dependence of the infrared bands for the CO bound to the Fea3 heme iron and to CuB. We have also studied the kinetics of ligand return from Fea3 to CuB using temperature derivative spectroscopy (TDS). Two classes of substates (alpha/beta) can be distinguished from their different properties with regard to the width of the IR band, the temperature dependence of the peak position, and the peak of the enthalpy distribution. The pronounced temperature dependence of the stretch frequencies in the beta conformation and the lack thereof in the alpha conformation implies very different dynamic behavior in the active site and reflects structural differences between the two conformations, most likely a shift of the position of CuB in response to a change in its stereochemical environment. Similar conformational changes will be necessary during the catalytic cycle of the enzyme when dioxygen is bound in the active site.

Ligand binding to heme proteins. VI. Interconversion of taxonomic substates in carbonmonoxymyoglobin.

The kinetic properties of the three taxonomic A substates of sperm whale carbonmonoxy myoglobin in 75% glycerol/buffer are studied by flash photolysis with monitoring in the infrared stretch bands of bound CO at nu(A0) approximately 1967 cm-1, nu(A1) approximately 1947 cm-1, and nu(A3) approximately 1929 cm-1 between 60 and 300 K. Below 160 K the photodissociated CO rebinds from the heme pocket, no interconversion among the A substates is observed, and rebinding in each A substate is nonexponential in time and described by a different temperature-independent distribution of enthalpy barriers with a different preexponential. Measurements in the electronic bands, e.g., the Soret, contain contributions of all three A substates and can, therefore, be only approximately modeled with a single enthalpy distribution and a single preexponential. The bond formation step at the heme is fastest for the A0 substate, intermediate for the A1 substate, and slowest for A3. Rebinding between 200 and 300 K displays several processes, including geminate rebinding, rebinding after ligand escape to the solvent, and interconversion among the A substates. Different kinetics are measured in each of the A bands for times shorter than the characteristic time of fluctuations among the A substates. At longer times, fluctuational averaging yields the same kinetics in all three A substates. The interconversion rates between A1 and A3 are determined from the time when the scaled kinetic traces of the two substates merge. Fluctuations between A1 and A3 are much faster than those between A0 and either A1 or A3, so A1 and A3 appear as one kinetic species in the exchange with A0. The maximum-entropy method is used to extract the distribution of rate coefficients for the interconversion process A0 <--> A1 + A3 from the flash photolysis data. The temperature dependencies of the A substate interconversion processes are fitted with a non-Arrhenius expression similar to that used to describe relaxation processes in glasses. At 300 K the interconversion time for A0 <--> A1 + A3 is 10 microseconds, and extrapolation yields approximately 1 ns for A1 <--> A3. The pronounced kinetic differences imply different structural rearrangements. Crystallographic data support this conclusion: They show that formation of the A0 substate involves a major change of the protein structure; the distal histidine rotates about the C(alpha)-C(beta) bond, and its imidazole sidechain swings out of the heme pocket into the solvent, whereas it remains in the heme pocket in the A1 <--> A3 interconversion. The fast A1 <--> A3 exchange is inconsistent with structural models that involve differences in the protonation between A1 and A3.

Ligand binding to anti-fluorescyl antibodies: stability of the antigen binding site.

The problem of protein stability is addressed with spectroscopic studies of equilibrium and kinetic properties of the binding of fluorescein to high-affinity monoclonal anti-fluorescyl antibodies (Mab 4-4-20), Fab fragments, and single-chain antibodies (SCA). SCA molecules contain only the variable domains of the antibody and an amino acid linker. The influence of glycerol on the antigen binding reaction is studied by circular dichroism, fluorescence, and absorption spectroscopy. The presence of glycerol in the solvent lowers the affinity of SCA for the ligand drastically, and the affinity even decreases toward lower temperatures. This effect is not observed in Fab and Mab. Analysis of the temperature jump kinetics shows that the dissociation reaction can be modeled as a two-state transition. The CD spectra indicate that the domain structure of the SCA remains unaltered in the presence of glycerol. Therefore, it is concluded that glycerol promotes the dissociation of the two variable domains of SCA. In Fab and Mab, the constant domains provide additional stabilization of the molecular structure at the antigen binding site.