Viral dynamics in vivo: limitations on estimates of intracellular delay and virus decay.

Anti-viral drug treatment of human immunodeficiency virus type I (HIV-1) and hepatitis B virus (HBV) infections causes rapid reduction in plasma virus load. Viral decline occurs in several phases and provides information on important kinetic constants of virus replication in vivo and pharmacodynamical properties. We develop a mathematical model that takes into account the intracellular phase of the viral life-cycle, defined as the time between infection of a cell and production of new virus particles. We derive analytic solutions for the dynamics following treatment with reverse transcriptase inhibitors, protease inhibitors, or a combination of both. For HIV-1, our results show that the phase of rapid decay in plasma virus (days 2-7) allows precise estimates for the turnover rate of productively infected cells. The initial quasi-stationary phase (days 0-1) and the transition phase (days 1-2) are explained by the combined effects of pharmacological and intracellular delays, the clearance of free virus particles, and the decay of infected cells. Reliable estimates of the first three quantities are not possible from data on virus load only; such estimates require additional measurements. In contrast with HIV-1, for HBV our model predicts that frequent early sampling of plasma virus will lead to reliable estimates of the free virus half-life and the pharmacological properties of the administered drug. On the other hand, for HBV the half-life of infected cells cannot be estimated from plasma virus decay.

Simultaneous determination of helical unwinding angles and intrinsic association constants in ligand–DNA complexes: The interaction between DNA and calichearubicin B

We present a helical unwinding assay for reversibly binding DNA ligands that uses closed circular DNA, topoisomerase I (Topo I), and two-dimensional agarose gel electrophoresis. Serially diluted Topo I relaxation reactions at constant DNA/ligand ratio are performed, and the resulting apparent unwinding of the closed circular DNA is used to calculate both ligand unwinding angle (φ) and intrinsic association constant (Ka). Mathematical treatment of apparent unwinding is formally analogous to that of apparent extinction coefficient data for optical binding titrations. Extrapolation to infinite DNA concentration yields the true unwinding angle of a given ligand and its association constant under Topo I relaxation conditions. Thus this assay delivers simultaneous structural and thermodynamic information describing the ligand–DNA complex. The utility of this assay has been demonstrated by using calichearubicin B (CRB), a synthetic hybrid molecule containing the anthraquinone chromophore of (DA) and the carbohydrate domain of calicheamicin γ1I. The unwinding angle for CRB calculated by this method is −5.3 ± 0.5°. Its Ka value is 0.20 × 106 M−1. For comparison, the unwinding angles of ethidium bromide and DA have been independently calculated, and the results are in agreement with canonical values for these compounds. Although a stronger binder to selected sites, CRB is a less potent unwinder than its parent compound DA. The assay requires only small amounts of ligand and offers an attractive option for analysis of DNA binding by synthetic and natural compounds.

Development of an in vitro mRNA decay system for Escherichia coli: Poly(A) polymerase I is necessary to trigger degradation

Using a novel Escherichia coli in vitro decay system in which polysomes are the source of both enzymes and mRNA, we demonstrate a requirement for poly(A) polymerase I (PAP I) in mRNA turnover. The in vitro decay of two different mRNAs (trxA and lpp) is triggered by the addition of ATP only when polysomes are prepared from a strain carrying the wild-type gene for PAP I (pcnB+). The relative decay rates of these two messages are similar in vitro and in vivo. Poly(A) tails are formed on both mRNAs, but no poly(A) tails are detected on the 3′ end of mature 23S rRNA. The size distribution of poly(A) tails generated in vitro, averaging 50 nt in length, is comparable to that previously reported in vivo. PAP I activity is associated exclusively with the polysomes. Exogenously added PAP I does not restore mRNA decay to PAP I− polysomes, suggesting that, in vivo, PAP I may be part of a multiprotein complex. The potential of this in vitro system for analyzing mRNA decay in E. coli is discussed.

Effective intercellular communication distances are determined by the relative time constants for cyto/chemokine secretion and diffusion

A cell’s ability to effectively communicate with a neighboring cell is essential for tissue function and ultimately for the organism to which it belongs. One important mode of intercellular communication is the release of soluble cyto- and chemokines. Once secreted, these signaling molecules diffuse through the surrounding medium and eventually bind to neighboring cell’s receptors whereby the signal is received. This mode of communication is governed both by physicochemical transport processes and cellular secretion rates, which in turn are determined by genetic and biochemical processes. The characteristics of transport processes have been known for some time, and information on the genetic and biochemical determinants of cellular function is rapidly growing. Simultaneous quantitative analysis of the two is required to systematically evaluate the nature and limitations of intercellular signaling. The present study uses a solitary cell model to estimate effective communication distances over which a single cell can meaningfully propagate a soluble signal. The analysis reveals that: (i) this process is governed by a single, key, dimensionless group that is a ratio of biological parameters and physicochemical determinants; (ii) this ratio has a maximal value; (iii) for realistic values of the parameters contained in this dimensionless group, it is estimated that the domain that a single cell can effectively communicate in is ≈250 μm in size; and (iv) the communication within this domain takes place in 10–30 minutes. These results have fundamental implications for interpretation of organ physiology and for engineering tissue function ex vivo.

Equilibrium constants and free energies in unfolding of proteins in urea solutions

A novel thermodynamic approach to the reversible unfolding of proteins in aqueous urea solutions has been developed based on the premise that urea ligands are bound cooperatively to the macromolecule. When successive stoichiometric binding constants have values larger than expected from statistical effects, an equation for moles of bound urea can be derived that contains imaginary terms. For a very steep unfolding curve, one can then show that the fraction of protein unfolded, B̄, depends on the square of the urea concentration, U, and is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\bar {B}=\frac{{\mathit{A}}^{{\mathit{2}}}_{{\mathit{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}n\bar {B}}{\mathit{U}}^{{\mathit{2}}}}{{\mathrm{1\hspace{.167em}+\hspace{.167em}}}{\mathit{A}}^{{\mathrm{2}}}_{{\mathrm{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}\bar {B}}{\mathit{U}}^{{\mathrm{2}}}}{\mathrm{.}}\end{equation*}\end{document} A12 is the binding constant as B̄→ 0, and λ is a parameter that reflects the augmentation in affinities of protein for urea as the moles bound increases to the saturation number, n. This equation provides an analytic expression that reproduces the unfolding curve with good precision, suggests a simple linear graphical procedure for evaluating A12 and λ, and leads to the appropriate standard free energy changes. The calculated ΔG° values reflect the coupling of urea binding with unfolding of the protein. Some possible implications of this analysis to protein folding in vivo are described.

Fluorescence correlation spectroscopy reveals fast optical excitation-driven intramolecular dynamics of yellow fluorescent proteins

Fast excitation-driven fluctuations in the fluorescence emission of yellow-shifted green fluorescent protein mutants T203Y and T203F, with S65G/S72A, are discovered in the 10−6–10−3-s time range, by using fluorescence correlation spectroscopy at 10−8 M. This intensity-dependent flickering is conspicuous at high pH, with rate constants independent of pH and viscosity with a minor temperature effect. The mean flicker rate increases linearly with excitation intensity for at least three decades, but the mean dark fraction of the molecules undergoing these dynamics is independent of illumination intensity over ≈6 × 102 to 5 × 106 W/cm2. These results suggest that optical excitation establishes an equilibration between two molecular states of different spectroscopic properties that are coupled only via the excited state as a gateway. This reversible excitation-driven transition has a quantum efficiency of ≈10−3. Dynamics of external protonation, reversibly quenching the fluorescence, are also observed at low pH in the 10- to 100-μs time range. The independence of these two bright–dark flicker processes implies the existence of at least two separate dark states of these green fluorescent protein mutants. Time-resolved fluorescence measurements reveal a single exponential decay of the excited state population with 3.8-ns lifetime, after 500-nm excitation, that is pH independent. Our fluorescence correlation spectroscopy results are discussed in terms of recent theoretical studies that invoke isomerization of the chromophore as a nonradiative channel of the excited state relaxation.

Femtosecond dynamics of DNA-mediated electron transfer

Diverse biophysical and biochemical studies have sought to understand electron transfer (ET) in DNA in part because of its importance to DNA damage and its repair. However, the dynamics and mechanisms of the elementary processes of ET in this medium are not fully understood and have been heavily debated. Two fundamental issues are the distance over which charge is transported and the time-scale on which the transport through the π-stack of the DNA base pairs may occur. With femtosecond resolution, we report direct observation in DNA of ultrafast ET, initiated by excitation of tethered ethidium (E), the intercalated electron acceptor (A); the electron donor (D) is 7-deazaguanine (Z), a modified base, placed at different, fixed distances from A. The ultrafast ET between these reactants in DNA has been observed with time constants of 5 ps and 75 ps and was found to be essentially independent of the D–A separation (10–17 Å). However, the ET efficiency does depend on the D–A distance. The 5-ps decay corresponds to direct ET observed from 7-deazaguanine but not guanine to E. From measurements of orientation anisotropies, we conclude that the slower 75-ps process requires the reorientation of E before ET, similar to E/nucleotide complexes in water. These results reveal the nature of ultrafast ET and its mechanism: in DNA, ET cannot be described as in proteins simply by a phenomenological parameter, β. Instead, the involvement of the base pairs controls the time scale and the degree of coherent transport.

Influence of follicular dendritic cells on decay of HIV during antiretroviral therapy

Drug treatment of HIV type 1 (HIV-1) infection leads to a rapid initial decay of plasma virus followed by a slower second phase of decay. To investigate the role of HIV-1 retained on follicular dendritic cells (FDCs) in this process, we have developed and analyzed a mathematical model for HIV-1 dynamics in lymphoid tissue (LT) that includes FDCs. Analysis of clinical data using this model indicates that decay of HIV-1 during therapy may be influenced by release of FDC-associated virus. The biphasic character of viral decay can be explained by reversible multivalent binding of HIV-1 to receptors on FDCs, indicating that the second phase of decay is not necessarily caused by long-lived or latently infected cells. Furthermore, viral clearance and death of short-lived productively infected cells may be faster than previously estimated. The model, with reasonable parameter values, is consistent with kinetic measurements of viral RNA in plasma, viral RNA on FDCs, productively infected cells in LT, and CD4+ T cells in LT during therapy.

Polysomal ribonuclease 1 exists in a latent form on polysomes prior to estrogen activation of mRNA decay

Estrogen induces a global change in the translation profile of Xenopus hepatocytes, replacing serum protein synthesis with production of the yolk protein precursor vitellogenin. This is accomplished by the coordinate destabilization of serum protein mRNAs and the transcriptional induction and subsequent stabilization of vitellogenin mRNA. Previous work identified an endonuclease activity whose appearance on polysomes correlated with the disappearance of serum protein mRNAs. This enzyme, polysomal ribonuclease 1 (PMR1), is a novel member of the peroxidase gene family. The current study examined the association of PMR1 with its mRNA targets on polysomes and mRNPs. The highest amount of polysome-bound PMR1 was observed prior to estrogen induction of mRNA decay. Its distribution on sucrose density gradients matched the absorbance profile of polysome-bound mRNA, suggesting that PMR1 forms a latent complex with mRNA. Following dissociation with EDTA the 62 kDa PMR1 sedimented with a larger complex of >670 kDa. Estrogen induces a 22-fold increase in unit enzymatic activity of polysome-bound PMR1, and a time-dependent loss of PMR1 from polysomes in a manner that mirrors the disappearance of albumin mRNA. These data suggest that the key step in the extensive estrogen-induced change in mRNA decay in Xenopus liver is activation of a latent mRNA endonuclease associated with its target mRNA.

Determination of the binding constants of the centromere protein Cbf1 to all 16 centromere DNAs of Saccharomyces cerevisiae

Cbf1p is a Saccharomyces cerevisiae chromatin protein belonging to the basic region helix–loop–helix leucine zipper (bHLHzip) family of DNA binding proteins. Cbf1p binds to a conserved element in the 5′-flanking region of methionine biosynthetic genes and to centromere DNA element I (CDEI) of S.cerevisiae centromeric DNA. We have determined the apparent equilibrium dissociation constants of Cbf1p binding to all 16 CDEI DNAs in gel retardation assays. Binding constants of full-length Cbf1p vary between 1.7 and 3.8 nM. However, the dissociation constants of a Cbf1p deletion variant that has been shown to be fully sufficient for Cbf1p function in vivo vary in a range between 3.2 and 12 nM. In addition, native polyacrylamide gel electrophoresis revealed distinct changes in the 3D structure of the Cbf1p/CEN complexes. We also show that the previously reported DNA binding stimulation activity of the centromere protein p64 functions on both the Cbf1 full-length protein and a deletion variant containing only the bHLHzip domain of Cbf1p. Our results suggest that centromeric DNA outside the consensus CDEI sequence and interaction of Cbf1p with adjacent centromere proteins contribute to the complex formation between Cbf1p and CEN DNA.

Nonsense-mediated Decay of mRNA for the Selenoprotein Phospholipid Hydroperoxide Glutathione Peroxidase Is Detectable in Cultured Cells but Masked or Inhibited in Rat Tissues

Previous studies of mRNA for classical glutathione peroxidase 1 (GPx1) demonstrated that hepatocytes of rats fed a selenium-deficient diet have less cytoplasmic GPx1 mRNA than hepatocytes of rats fed a selenium-adequate diet. This is because GPx1 mRNA is degraded by the surveillance pathway called nonsense-mediated mRNA decay (NMD) when the selenocysteine codon is recognized as nonsense. Here, we examine the mechanism by which the abundance of phospholipid hydroperoxide glutathione peroxidase (PHGPx) mRNA, another selenocysteine-encoding mRNA, fails to decrease in the hepatocytes and testicular cells of rats fed a selenium-deficient diet. We demonstrate with cultured NIH3T3 fibroblasts or H35 hepatocytes transiently transfected with PHGPx gene variants under selenium-supplemented or selenium-deficient conditions that PHGPx mRNA is, in fact, a substrate for NMD when the selenocysteine codon is recognized as nonsense. We also demonstrate that the endogenous PHGPx mRNA of untransfected H35 cells is subject to NMD. The failure of previous reports to detect the NMD of PHGPx mRNA in cultured cells is likely attributable to the expression of PHGPx cDNA rather than the PHGPx gene. We conclude that 1) the sequence of the PHGPx gene is adequate to support the NMD of product mRNA, and 2) there is a mechanism in liver and testis but not cultured fibroblasts and hepatocytes that precludes or masks the NMD of PHGPx mRNA.

On the nucleation and growth of amyloid beta-protein fibrils: detection of nuclei and quantitation of rate constants.

We have studied the fibrillogenesis of synthetic amyloid beta-protein-(1-40) fragment (A beta) in 0.1 M HCl. At low pH, A beta formed fibrils at a rate amenable to detailed monitoring by quasi-elastic light-scattering spectroscopy. Examination of the fibrils with circular dichroism spectroscopy and electron microscopy showed them to be highly similar to those found in amyloid plaques. We determined the hydrodynamic radii of A beta aggregates during the entire process of fibril nucleation and growth. Above an A beta concentration of approximately 0.1 mM, the initial rate of elongation and the final size of fibrils were independent of A beta concentration. Below an A beta concentration of 0.1 mM, the initial elongation rate was proportional to the peptide concentration, and the resulting fibrils were significantly longer than those formed at higher concentration. We also found that the surfactant n-dodecylhexaoxyethylene glycol monoether (C12E6) slowed nucleation and elongation of fibrils in a concentration-dependent manner. Our observations are consistent with a model of A beta fibrillogenesis that includes the following key steps: (i) peptide micelles form above a certain critical A beta concentration, (ii) fibrils nucleate within these micelles or on heterogeneous nuclei (seeds), and (iii) fibrils grow by irreversible binding of monomers to fibril ends. Interpretation of our data enabled us to determine the sizes of fibril nuclei and A beta micelles and the rates of fibril nucleation (from micelles) and fibril elongation. Our approach provides a powerful means for the quantitative assay of A beta fibrillogenesis.

Electrogenic light reactions in photosystem I: resolution of electron-transfer rates between the iron-sulfur centers.

Flash-induced voltage changes (electrogenic events) in photosystem I particles from spinach, oriented in a phospholipid layer, have been studied at room temperature on a time scale ranging from 1 micros to several seconds. A phospholipid layer containing photosystem I particles was adsorbed to a Teflon film separating two aqueous compartments. Voltage changes were measured across electrodes immersed in the compartments. In the absence of added electron donors and acceptors, a multiphasic voltage increase, associated with charge separation, was followed by a decrease, associated with charge recombination. Several kinetic phases were resolved: a rapid (<1 micros) increase, ascribed to electron transfer from the primary electron donor P700 to the iron-sulfur electron acceptor FB, was followed by a slower, biphasic increase with time constants of 30 and 200 micros. The 30-micros phase is assigned to electron transfer from FB to the iron-sulfur center FA. The voltage decrease had a time constant of 90 ms, ascribed to charge recombination from FA to P700. Upon chemical prereduction of FA and FB the 30- and 200-micros phases disappeared and the decay time constant was accelerated to 330 micros, assigned to charge recombination from the phylloquinone electron acceptor (A1) or the iron-sulfur center FX to P700.