Estimating sequestered parasite population dynamics in cerebral malaria

Clinical investigation of malaria is hampered by the lack of a method for estimating the number of parasites that are sequestered in the tissues, for it is these parasites that are thought to be crucial to the pathogenesis of life-threatening complications such as cerebral malaria. We present a method of estimating this hidden population by using clinical observations of peripheral parasitemia combined with an age-structured mathematical model of the parasite erythrocyte cycle. Applying the model to data from 217 Gambian children undergoing treatment for cerebral malaria we conclude that although artemether clears parasitemia more rapidly than quinine, the clearance of sequestered parasites is similar for the two drugs. The estimated sequestered mass was found to be a more direct predictor of fatal outcome than clinically observed parasitemia. This method allows a sequential analysis of sequestered parasite population dynamics in children suffering from cerebral malaria, and the results offer a possible explanation for why artemether provides less advantage than might have been expected over quinine in reducing mortality despite its rapid effect on circulating parasites.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200°C

The rubredoxin protein from the hyperthermophilic archaebacterium Pyrococcus furiosus was examined by a hydrogen exchange method. Even though the protein does not exhibit reversible thermal unfolding, one can determine its stability parameters—free energy, enthalpy, entropy, and melting temperature—and also the distribution of stability throughout the protein, by using hydrogen exchange to measure the reversible cycling of the protein between native and unfolded states that occurs even under native conditions.

Regular and stochastic dynamics in the real quadratic family

We prove the Regulat or Stochastic Conjecture for the real quadratic family which asserts that almost every real quadratic map Pc, c ∈ [−2, 1/4], has either an attracting cycle or an absolutely continuous invariant measure.

Different dynamics in nuclear entry of subunits of the repair/transcription factor TFIIH

We report here the different ways in which four subunits of the basal transcription/repair factor TFIIH (XPB, XPD, p62 and p44) and the damage recognition XPC repair protein can enter the nucleus. We examined their nuclear localization by transiently expressing the gene products tagged with the enhanced green fluorescent protein (EGFP) in transfected 3T3 cells. In agreement with the identification of more than one putative nuclear localization signal (NLS) in their protein sequences, XPB, XPC, p62 and p44 chimeras were rapidly sorted to the nucleus. In contrast, the XPD–EGFP chimeras appeared mainly localized in the cytoplasm, with a minor fraction of transfectants showing the EGFP-based fluorescence also in the nucleus. The ability of the XPD chimeras to enter the nucleus was confirmed by western blotting on fractionated cell extracts and by functional complementation of the repair defect in the UV5 rodent cells, mutated in the XPD homologous gene. By deletion mutagenesis, we were unable to identify any sequence specific for nuclear localization. In particular, deletion of the putative NLS failed to affect subcellular localization and, conversely, the C-terminal part of XPD containing the putative NLS showed no specific nuclear accumulation. These findings suggest that the nuclear entry of XPD depends on its complexation with other proteins in the cytoplasm, possibly other components of the TFIIH complex.

Cell Cycle-Dependent Changes in Microtubule Dynamics in Living Cells Expressing Green Fluorescent Protein-α TubulinV⃞

LLCPK-1 cells were transfected with a green fluorescent protein (GFP)-α tubulin construct and a cell line permanently expressing GFP-α tubulin was established (LLCPK-1α). The mitotic index and doubling time for LLCPK-1α were not significantly different from parental cells. Quantitative immunoblotting showed that 17% of the tubulin in LLCPK-1α cells was GFP-tubulin; the level of unlabeled tubulin was reduced to 82% of that in parental cells. The parameters of microtubule dynamic instability were compared for interphase LLCPK-1α and parental cells injected with rhodamine-labeled tubulin. Dynamic instability was very similar in the two cases, demonstrating that LLCPK-1α cells are a useful tool for analysis of microtubule dynamics throughout the cell cycle. Comparison of astral microtubule behavior in mitosis with microtubule behavior in interphase demonstrated that the frequency of catastrophe increased twofold and that the frequency of rescue decreased nearly fourfold in mitotic compared with interphase cells. The percentage of time that microtubules spent in an attenuated state, or pause, was also dramatically reduced, from 73.5% in interphase to 11.4% in mitosis. The rates of microtubule elongation and rapid shortening were not changed; overall dynamicity increased 3.6-fold in mitosis. Microtubule release from the centrosome and a subset of differentially stable astral microtubules were also observed. The results provide the first quantitative measurements of mitotic microtubule dynamics in mammalian cells.

Quantitative ER ↔ Golgi Transport Kinetics and Protein Separation upon Golgi Exit Revealed by Vesicular Integral Membrane Protein 36 Dynamics in Live CellsV⃞

To quantitatively investigate the trafficking of the transmembrane lectin VIP36 and its relation to cargo-containing transport carriers (TCs), we analyzed a C-terminal fluorescent-protein (FP) fusion, VIP36-SP-FP. When expressed at moderate levels, VIP36-SP-FP localized to the endoplasmic reticulum, Golgi apparatus, and intermediate transport structures, and colocalized with epitope-tagged VIP36. Temperature shift and pharmacological experiments indicated VIP36-SP-FP recycled in the early secretory pathway, exhibiting trafficking representative of a class of transmembrane cargo receptors, including the closely related lectin ERGIC53. VIP36-SP-FP trafficking structures comprised tubules and globular elements, which translocated in a saltatory manner. Simultaneous visualization of anterograde secretory cargo and VIP36-SP-FP indicated that the globular structures were pre-Golgi carriers, and that VIP36-SP-FP segregated from cargo within the Golgi and was not included in post-Golgi TCs. Organelle-specific bleach experiments directly measured the exchange of VIP36-SP-FP between the Golgi and endoplasmic reticulum (ER). Fitting a two-compartment model to the recovery data predicted first order rate constants of 1.22 ± 0.44%/min for ER → Golgi, and 7.68 ± 1.94%/min for Golgi → ER transport, revealing a half-time of 113 ± 70 min for leaving the ER and 1.67 ± 0.45 min for leaving the Golgi, and accounting for the measured steady-state distribution of VIP36-SP-FP (13% Golgi/87% ER). Perturbing transport with AlF4− treatment altered VIP36-SP-GFP distribution and changed the rate constants. The parameters of the model suggest that relatively small differences in the first order rate constants, perhaps manifested in subtle differences in the tendency to enter distinct TCs, result in large differences in the steady-state localization of secretory components.

Functional dynamics in the active site of the ribonuclease binase

Binase, a member of a family of microbial guanyl-specific ribonucleases, catalyzes the endonucleotic cleavage of single-stranded RNA. It shares 82% amino acid identity with the well-studied protein barnase. We used NMR spectroscopy to study the millisecond dynamics of this small enzyme, using several methods including the measurement of residual dipolar couplings in solution. Our data show that the active site of binase is flanked by loops that are flexible at the 300-μs time scale. One of the catalytic residues, His-101, is located on such a flexible loop. In contrast, the other catalytic residue, Glu-72, is located on a β-sheet, and is static. The residues Phe-55, part of the guanine base recognition site, and Tyr-102, stabilizing the base, are the most dynamic. Our findings suggest that binase possesses an active site that has a well-defined bottom, but which has sides that are flexible to facilitate substrate access/egress, and to deliver one of the catalytic residues. The motion in these loops does not change on complexation with the inhibitor d(CGAG) and compares well with the maximum kcat (1,500 s−1) of these ribonucleases. This observation indicates that the NMR-measured loop motions reflect the opening necessary for product release, which is apparently rate limiting for the overall turnover.

Domain structure and dynamics in the helical filaments formed by RecA and Rad51 on DNA

Both the bacterial RecA protein and the eukaryotic Rad51 protein form helical nucleoprotein filaments on DNA that catalyze strand transfer between two homologous DNA molecules. However, only the ATP-binding cores of these proteins have been conserved, and this same core is also found within helicases and the F1-ATPase. The C-terminal domain of the RecA protein forms lobes within the helical RecA filament. However, the Rad51 proteins do not have the C-terminal domain found in RecA, but have an N-terminal extension that is absent in the RecA protein. Both the RecA C-terminal domain and the Rad51 N-terminal domain bind DNA. We have used electron microscopy to show that the lobes of the yeast and human Rad51 filaments appear to be formed by N-terminal domains. These lobes are conformationally flexible in both RecA and Rad51. Within RecA filaments, the change between the “active” and “inactive” states appears to mainly involve a large movement of the C-terminal lobe. The N-terminal domain of Rad51 and the C-terminal domain of RecA may have arisen from convergent evolution to play similar roles in the filaments.

Excited-state dynamics in photosystem II: Insights from the x-ray crystal structure

The heart of oxygenic photosynthesis is photosystem II (PSII), a multisubunit protein complex that uses solar energy to drive the splitting of water and production of molecular oxygen. The effectiveness of the photochemical reaction center of PSII depends on the efficient transfer of excitation energy from the surrounding antenna chlorophylls. A kinetic model for PSII, based on the x-ray crystal structure coordinates of 37 antenna and reaction center pigment molecules, allows us to map the major energy transfer routes from the antenna chlorophylls to the reaction center chromophores. The model shows that energy transfer to the reaction center is slow compared with the rate of primary electron transport and depends on a few bridging chlorophyll molecules. This unexpected energetic isolation of the reaction center in PSII is similar to that found in the bacterial photosystem, conflicts with the established view of the photophysics of PSII, and may be a functional requirement for primary photochemistry in photosynthesis. In addition, the model predicts a value for the intrinsic photochemical rate constant that is 4 times that found in bacterial reaction centers.

Nonlinear dynamics in ventricular fibrillation.

Electrogram recordings of ventricular fibrillation appear complex and possibly chaotic. However, sequences of beat-to-beat intervals obtained from these recordings are generally short, making it difficult to explicitly demonstrate nonlinear dynamics. Motivated by the work of Sugihara on atmospheric dynamics and the Durbin-Watson test for nonlinearity, we introduce a new statistical test that recovers significant dynamical patterns from smoothed lag plots. This test is used to show highly significant nonlinear dynamics in a stable canine model of ventricular fibrillation.

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.