From snapshot to movie: φ analysis of protein folding transition states taken one step further

Kinetic anomalies in protein folding can result from changes of the kinetic ground states (D, I, and N), changes of the protein folding transition state, or both. The 102-residue protein U1A has a symmetrically curved chevron plot which seems to result mainly from changes of the transition state. At low concentrations of denaturant the transition state occurs early in the folding reaction, whereas at high denaturant concentration it moves close to the native structure. In this study we use this movement to follow continuously the formation and growth of U1A's folding nucleus by φ analysis. Although U1A's transition state structure is generally delocalized and displays a typical nucleation–condensation pattern, we can still resolve a sequence of folding events. However, these events are sufficiently coupled to start almost simultaneously throughout the transition state structure.

High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology

The worldwide threat of tuberculosis to human health emphasizes the need to develop novel approaches to a global epidemiological surveillance. The current standard for Mycobacterium tuberculosis typing based on IS6110 restriction fragment length polymorphism (RFLP) suffers from the difficulty of comparing data between independent laboratories. Here, we propose a high-resolution typing method based on variable number tandem repeats (VNTRs) of genetic elements named mycobacterial interspersed repetitive units (MIRUs) in 12 human minisatellite-like regions of the M. tuberculosis genome. MIRU-VNTR profiles of 72 different M. tuberculosis isolates were established by PCR analysis of all 12 loci. From 2 to 8 MIRU-VNTR alleles were identified in the 12 regions in these strains, which corresponds to a potential of over 16 million different combinations, yielding a resolution power close to that of IS6110-RFLP. All epidemiologically related isolates tested were perfectly clustered by MIRU-VNTR typing, indicating that the stability of these MIRU-VNTRs is adequate to track outbreak episodes. The correlation between genetic relationships inferred from MIRU-VNTR and IS6110-RFLP typing was highly significant. Compared with IS6110-RFLP, high-resolution MIRU-VNTR typing has the considerable advantages of being fast, appropriate for all M. tuberculosis isolates, including strains that have a few IS6110 copies, and permitting easy and rapid comparison of results from independent laboratories. This typing method opens the way to the construction of digital global databases for molecular epidemiology studies of M. tuberculosis.

NOMAD: a versatile strategy for in vitro DNA manipulation applied to promoter analysis and vector design.

Molecular analysis of complex modular structures, such as promoter regions or multi-domain proteins, often requires the creation of families of experimental DNA constructs having altered composition, order, or spacing of individual modules. Generally, creation of every individual construct of such a family uses a specific combination of restriction sites. However, convenient sites are not always available and the alternatives, such as chemical resynthesis of the experimental constructs or engineering of different restriction sites onto the ends of DNA fragments, are costly and time consuming. A general cloning strategy (nucleic acid ordered assembly with directionality, NOMAD; WWW resource locator http:@Lmb1.bios.uic.edu/NOMAD/NOMAD.htm l) is proposed that overcomes these limitations. Use of NOMAD ensures that the production of experimental constructs is no longer the rate-limiting step in applications that require combinatorial rearrangement of DNA fragments. NOMAD manipulates DNA fragments in the form of "modules" having a standardized cohesive end structure. Specially designed "assembly vectors" allow for sequential and directional insertion of any number of modules in an arbitrary predetermined order, using the ability of type IIS restriction enzymes to cut DNA outside of their recognition sequences. Studies of regulatory regions in DNA, such as promoters, replication origins, and RNA processing signals, construction of chimeric proteins, and creation of new cloning vehicles, are among the applications that will benefit from using NOMAD.

From local to global in quasiconformal structures.

We exhibit a large class of metric spaces whose infinitesimal quasiconformal structure is strong enough to capture the global quasiconformal structure. A sufficient condition for this to happen is described in terms of a Poincaré-type inequality.

Observation via one-dimensional 13Calpha NMR of local conformational substates in thermal unfolding equilibria of a synthetic analog of the GCN4 leucine zipper.

Synthesis of a 33-residue, capped leucine zipper analogous to that in GCN4 is reported. Histidine and arginine residues are mutated to lysine to reduce the unfolding temperature. CD and ultracentrifugation studies indicate that the molecule is a two-stranded coiled coil under benign conditions. Versions of the same peptide are made with 99% 13Calpha at selected sites. One-dimensional 13C NMR spectra are assigned by inspection and used to study thermal unfolding equilibria over the entire transition from 8 to 73 degrees C. Spectra at the temperature extremes establish the approximate chemical shifts for folded and unfolded forms at each labeled site. Resonances for each amino acid appear at both locations at intermediate T, indicating that folded and unfolded forms interconvert slowly (> >2 ms) on the NMR time scale. Moreover, near room temperature, the structured form's resonance is double at several, but not all, sites, indicating at least two slowly interconverting, structured, local conformational substates. Analysis of the dynamics is possible. For example, near room temperature at the Val-9, Ala-24, and Gly-31 positions, the equilibrium constant for interconversion of the two structured forms is near unity and the time scale is > or= 10-20 ms.

Mathematical analysis of activation thresholds in enzyme-catalyzed positive feedbacks: application to the feedbacks of blood coagulation.

A hierarchy of enzyme-catalyzed positive feedback loops is examined by mathematical and numerical analysis. Four systems are described, from the simplest, in which an enzyme catalyzes its own formation from an inactive precursor, to the most complex, in which two sequential feedback loops act in a cascade. In the latter we also examine the function of a long-range feedback, in which the final enzyme produced in the second loop activates the initial step in the first loop. When the enzymes generated are subject to inhibition or inactivation, all four systems exhibit threshold properties akin to excitable systems like neuron firing. For those that are amenable to mathematical analysis, expressions are derived that relate the excitation threshold to the kinetics of enzyme generation and inhibition and the initial conditions. For the most complex system, it was expedient to employ numerical simulation to demonstrate threshold behavior, and in this case long-range feedback was seen to have two distinct effects. At sufficiently high catalytic rates, this feedback is capable of exciting an otherwise subthreshold system. At lower catalytic rates, where the long-range feedback does not significantly affect the threshold, it nonetheless has a major effect in potentiating the response above the threshold. In particular, oscillatory behavior observed in simulations of sequential feedback loops is abolished when a long-range feedback is present.

Calcium block of Na+ channels and its effect on closing rate

Calcium ion transiently blocks Na+ channels, and it shortens the time course for closing of their activation gates. We examined the relation between block and closing kinetics by using the Na+ channels natively expressed in GH3 cells, a clonal line of rat pituitary cells. To simplify analysis, inactivation of the Na+ channels was destroyed by including papain in the internal medium. All divalent cations tested, and trivalent La3+, blocked a progressively larger fraction of the channels as their concentration increased, and they accelerated the closing of the Na+ channel activation gate. For calcium, the most extensively studied cation, there is an approximately linear relation between the fraction of the channels that are calcium-blocked and the closing rate. Extrapolation of the data to very low calcium suggests that closing rate is near zero when there is no block. Analysis shows that, almost with certainty, the channels can close when occupied by calcium. The analysis further suggests that the channels close preferentially or exclusively from the calcium-blocked state.

Dissociation of HIV-1 from follicular dendritic cells during HAART: Mathematical analysis

Follicular dendritic cells (FDC) provide a reservoir for HIV type 1 (HIV-1) that may reignite infection if highly active antiretroviral therapy (HAART) is withdrawn before virus on FDC is cleared. To estimate the treatment time required to eliminate HIV-1 on FDC, we develop deterministic and stochastic models for the reversible binding of HIV-1 to FDC via ligand–receptor interactions and examine the consequences of reducing the virus available for binding to FDC. Analysis of these models shows that the rate at which HIV-1 dissociates from FDC during HAART is biphasic, with an initial period of rapid decay followed by a period of slower exponential decay. The speed of the slower second stage of dissociation and the treatment time required to eradicate the FDC reservoir of HIV-1 are insensitive to the number of virions bound and their degree of attachment to FDC before treatment. In contrast, the expected time required for dissociation of an individual virion from FDC varies sensitively with the number of ligands attached to the virion that are available to interact with receptors on FDC. Although most virions may dissociate from FDC on the time scale of days to weeks, virions coupled to a higher-than-average number of ligands may persist on FDC for years. This result suggests that HAART may not be able to clear all HIV-1 trapped on FDC and that, even if clearance is possible, years of treatment will be required.

Tight frames of k-plane ridgelets and the problem of representing objects that are smooth away from d-dimensional singularities in Rn

For each pair (n, k) with 1 ≤ k < n, we construct a tight frame (ρλ : λ ∈ Λ) for L2 (Rn), which we call a frame of k-plane ridgelets. The intent is to efficiently represent functions that are smooth away from singularities along k-planes in Rn. We also develop tools to help decide whether k-plane ridgelets provide the desired efficient representation. We first construct a wavelet-like tight frame on the X-ray bundle χn,k—the fiber bundle having the Grassman manifold Gn,k of k-planes in Rn for base space, and for fibers the orthocomplements of those planes. This wavelet-like tight frame is the pushout to χn,k, via the smooth local coordinates of Gn,k, of an orthonormal basis of tensor Meyer wavelets on Euclidean space Rk(n−k) × Rn−k. We then use the X-ray isometry [Solmon, D. C. (1976) J. Math. Anal. Appl. 56, 61–83] to map this tight frame isometrically to a tight frame for L2(Rn)—the k-plane ridgelets. This construction makes analysis of a function f ∈ L2(Rn) by k-plane ridgelets identical to the analysis of the k-plane X-ray transform of f by an appropriate wavelet-like system for χn,k. As wavelets are typically effective at representing point singularities, it may be expected that these new systems will be effective at representing objects whose k-plane X-ray transform has a point singularity. Objects with discontinuities across hyperplanes are of this form, for k = n − 1.

Identification of inflections in T-cell counts among HIV-1-infected individuals and relationship with progression to clinical AIDS

Studies of circulating T (CD3+) lymphocytes have shown that on a population basis T-cell numbers remain stable for many years after HIV-1 infection (blind T-cell homeostasis), but decline rapidly beginning approximately 1.5–2.5 years before the onset of clinical AIDS. We derived a general method for defining the loss of homeostasis on the individual level and for determining the prevalence of homeostasis loss according to HIV status and the occurrence of AIDS in more than 5,000 men enrolled in the Multicenter AIDS Cohort Study. We used a segmented regression model for log10 CD3+ cell counts that included separate T-cell trajectories before and after a time (the T-cell inflection point) where the loss of T-cell homeostasis was most likely to have occurred. The average slope of CD3+ lymphocyte counts before the inflection point was close to zero for HIV− and HIV+ men, consistent with blind T-cell homeostasis. After the inflection point, the HIV+ individuals who developed AIDS generally showed a dramatic decline in CD3+ cell counts relative to HIV− men and HIV+ men not developing AIDS. A CD3+ cell decline of greater than 10 percent per year was present in 77% of HIV+ men developing AIDS but in only 23% of HIV+ men with no onset of AIDS. Our findings at the individual level support the blind T-cell homeostasis hypothesis and provide strong evidence that the loss of homeostasis is an important mechanism in the pathogenesis of the severe immunodeficiency that characterizes the late stages of HIV infection.

Root Xylem Embolisms and Refilling. Relation to Water Potentials of Soil, Roots, and Leaves, and Osmotic Potentials of Root Xylem Sap1

Embolism and refilling of vessels was monitored directly by cryomicroscopy of field-grown corn (Zea mays L.) roots. To test the reliability of an earlier study showing embolism refilling in roots at negative leaf water potentials, embolisms were counted, and root water potentials (Ψroot) and osmotic potentials of exuded xylem sap from the same roots were measured by isopiestic psychrometry. All vessels were full at dawn (Ψroot −0.1 MPa). Embolisms were first seen in late metaxylem vessels at 8 am. Embolized late metaxylem vessels peaked at 50% at 10 am (Ψroot −0.1 MPa), fell to 44% by 12 pm (Ψroot −0.23 MPa), then dropped steadily to zero by early evening (Ψroot −0.28 MPa). Transpiration was highest (8.5 μg cm−2 s−1) between 12 and 2 pm when the percentage of vessels embolized was falling. Embolized vessels were refilled by liquid moving through their lateral walls. Xylem sap was very low in solutes. The mechanism of vessel refilling, when Ψroot is negative, requires further investigation. Daily embolism and refilling in roots of well-watered plants is a normal occurrence and may be a component of an important hydraulic signaling mechanism between roots and shoots.

Quantitative analysis of chromosomal CGH in human breast tumors associates copy number abnormalities with p53 status and patient survival

We present a general method for rigorously identifying correlations between variations in large-scale molecular profiles and outcomes and apply it to chromosomal comparative genomic hybridization data from a set of 52 breast tumors. We identify two loci where copy number abnormalities are correlated with poor survival outcome (gain at 8q24 and loss at 9q13). We also identify a relationship between abnormalities at two loci and the mutational status of p53. Gain at 8q24 and loss at 5q15-5q21 are linked with mutant p53. The 9q and 5q losses suggest the possibility of gene products involved in breast cancer progression. The analytical techniques are general and also are applicable to the analysis of array-based expression data.