Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Population dynamics, demographic stochasticity, and the evolution of cooperation

A basic evolutionary problem posed by the Iterated Prisoner’s Dilemma game is to understand when the paradigmatic cooperative strategy Tit-for-Tat can invade a population of pure defectors. Deterministically, this is impossible. We consider the role of demographic stochasticity by embedding the Iterated Prisoner’s Dilemma into a population dynamic framework. Tit-for-Tat can invade a population of defectors when their dynamics exhibit short episodes of high population densities with subsequent crashes and long low density periods with strong genetic drift. Such dynamics tend to have reddened power spectra and temporal distributions of population size that are asymmetric and skewed toward low densities. The results indicate that ecological dynamics are important for evolutionary shifts between adaptive peaks.

HIV-1 and T cell dynamics after interruption of highly active antiretroviral therapy (HAART) in patients with a history of sustained viral suppression

Identifying the immunologic and virologic consequences of discontinuing antiretroviral therapy in HIV-infected patients is of major importance in developing long-term treatment strategies for patients with HIV-1 infection. We designed a trial to characterize these parameters after interruption of highly active antiretroviral therapy (HAART) in patients who had maintained prolonged viral suppression on antiretroviral drugs. Eighteen patients with CD4+ T cell counts ≥ 350 cells/μl and viral load below the limits of detection for ≥1 year while on HAART were enrolled prospectively in a trial in which HAART was discontinued. Twelve of these patients had received prior IL-2 therapy and had low frequencies of resting, latently infected CD4 cells. Viral load relapse to >50 copies/ml occurred in all 18 patients independent of prior IL-2 treatment, beginning most commonly during weeks 2–3 after cessation of HAART. The mean relapse rate constant was 0.45 (0.20 log10 copies) day−1, which was very similar to the mean viral clearance rate constant after drug resumption of 0.35 (0.15 log10 copies) day−1 (P = 0.28). One patient experienced a relapse delay to week 7. All patients except one experienced a relapse burden to >5,000 RNA copies/ml. Ex vivo labeling with BrdUrd showed that CD4 and CD8 cell turnover increased after withdrawal of HAART and correlated with viral load whereas lymphocyte turnover decreased after reinitiation of drug treatment. Virologic relapse occurs rapidly in patients who discontinue suppressive drug therapy, even in patients with a markedly diminished pool of resting, latently infected CD4+ T cells.

Antigen-driven CD4+ T cell and HIV-1 dynamics: Residual viral replication under highly active antiretroviral therapy

Antigen-induced stimulation of the immune system can generate heterogeneity in CD4+ T cell division rates capable of explaining the temporal patterns seen in the decay of HIV-1 plasma RNA levels during highly active antiretroviral therapy. Posttreatment increases in peripheral CD4+ T cell counts are consistent with a mathematical model in which host cell redistribution between lymph nodes and peripheral blood is a function of viral burden. Model fits to patient data suggest that, although therapy reduces HIV replication below replacement levels, substantial residual replication continues. This residual replication has important consequences for long-term therapy and the evolution of drug resistance and represents a challenge for future treatment strategies.

Dynamics of fluorescence fluctuations in green fluorescent protein observed by fluorescence correlation spectroscopy

We have investigated the pH dependence of the dynamics of conformational fluctuations of green fluorescent protein mutants EGFP (F64L/S65T) and GFP-S65T in small ensembles of molecules in solution by using fluorescence correlation spectroscopy (FCS). FCS utilizes time-resolved measurements of fluctuations in the molecular fluorescence emission for determination of the intrinsic dynamics and thermodynamics of all processes that affect the fluorescence. Fluorescence excitation of a bulk solution of EGFP decreases to zero at low pH (pKa = 5.8) paralleled by a decrease of the absorption at 488 nm and an increase at 400 nm. Protonation of the hydroxyl group of Tyr-66, which is part of the chromophore, induces these changes. When FCS is used the fluctuations in the protonation state of the chromophore are time resolved. The autocorrelation function of fluorescence emission shows contributions from two chemical relaxation processes as well as diffusional concentration fluctuations. The time constant of the fast, pH-dependent chemical process decreases with pH from 300 μs at pH 7 to 45 μs at pH 5, while the time-average fraction of molecules in a nonfluorescent state increases to 80% in the same range. A second, pH-independent, process with a time constant of 340 μs and an associated fraction of 13% nonfluorescent molecules is observed between pH 8 and 11, possibly representing an internal proton transfer process and associated conformational rearrangements. The FCS data provide direct measures of the dynamics and the equilibrium properties of the protonation processes. Thus FCS is a convenient, intrinsically calibrated method for pH measurements in subfemtoliter volumes with nanomolar concentrations of EGFP.

Flickering fusion pores comparable with initial exocytotic pores occur in protein-free phospholipid bilayers

For the act of membrane fusion, there are two competing, mutually exclusive molecular models that differ in the structure of the initial pore, the pathway for ionic continuity between formerly separated volumes. Because biological “fusion pores” can be as small as ionic channels or gap junctions, one model posits a proteinaceous initial fusion pore. Because biological fusion pore conductance varies widely, another model proposes a lipidic initial pore. We have found pore opening and flickering during the fusion of protein-free phospholipid vesicles with planar phospholipid bilayers. Fusion pore formation appears to follow the coalescence of contacting monolayers to create a zone of hemifusion where continuity between the two adherent membranes is lipidic, but not aqueous. Hypotonic stress, causing tension in the vesicle membrane, promotes complete fusion. Pores closed soon after opening (flickering), and the distribution of fusion pore conductance appears similar to the distribution of initial fusion pores in biological fusion. Because small flickering pores can form in the absence of protein, the existence of small pores in biological fusion cannot be an argument in support of models based on proteinaceous pores. Rather, these results support the model of a lipidic fusion pore developing within a hemifused contact site.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].