On the magnitude of the electrostatic contribution to ligand-DNA interactions.

A model based on the nonlinear Poisson-Boltzmann equation is used to study the electrostatic contribution to the binding free energy of a simple intercalating ligand, 3,8-diamino-6-phenylphenanthridine, to DNA. We find that the nonlinear Poisson-Boltzmann model accurately describes both the absolute magnitude of the pKa shift of 3,8-diamino-6-phenylphenanthridine observed upon intercalation and its variation with bulk salt concentration. Since the pKa shift is directly related to the total electrostatic binding free energy of the charged and neutral forms of the ligand, the accuracy of the calculations implies that the electrostatic contributions to binding are accurately predicted as well. Based on our results, we have developed a general physical description of the electrostatic contribution to ligand-DNA binding in which the electrostatic binding free energy is described as a balance between the coulombic attraction of a ligand to DNA and the disruption of solvent upon binding. Long-range coulombic forces associated with highly charged nucleic acids provide a strong driving force for the interaction of cationic ligands with DNA. These favorable electrostatic interactions are, however, largely compensated for by unfavorable changes in the solvation of both the ligand and the DNA upon binding. The formation of a ligand-DNA complex removes both charged and polar groups at the binding interface from pure solvent while it displaces salt from around the nucleic acid. As a result, the total electrostatic binding free energy is quite small. Consequently, nonpolar interactions, such as tight packing and hydrophobic forces, must play a significant role in ligand-DNA stability.

Aggregate community properties and the strength of species’ interactions

Although humanity depends on the continued, aggregate functioning of natural ecosystems, few studies have explored the impact of community structure on the stability of aggregate community properties. Here we derive the stability of the aggregate property of community biomass as a function of species’ competition coefficients for a two-species model. The model predicts that the stability of community biomass is relatively independent of the magnitude of the interaction strengths. Instead, the degree of asymmetry of the interactions appears to be key to community stability.

A unified theory of carcinogenesis based on order–disorder transitions in DNA structure as studied in the human ovary and breast

Fourier transform-infrared/statistics models demonstrate that the malignant transformation of morphologically normal human ovarian and breast tissues involves the creation of a high degree of structural modification (disorder) in DNA, before restoration of order in distant metastases. Order–disorder transitions were revealed by methods including principal components analysis of infrared spectra in which DNA samples were represented by points in two-dimensional space. Differences between the geometric sizes of clusters of points and between their locations revealed the magnitude of the order–disorder transitions. Infrared spectra provided evidence for the types of structural changes involved. Normal ovarian DNAs formed a tight cluster comparable to that of normal human blood leukocytes. The DNAs of ovarian primary carcinomas, including those that had given rise to metastases, had a high degree of disorder, whereas the DNAs of distant metastases from ovarian carcinomas were relatively ordered. However, the spectra of the metastases were more diverse than those of normal ovarian DNAs in regions assigned to base vibrations, implying increased genetic changes. DNAs of normal female breasts were substantially disordered (e.g., compared with the human blood leukocytes) as were those of the primary carcinomas, whether or not they had metastasized. The DNAs of distant breast cancer metastases were relatively ordered. These findings evoke a unified theory of carcinogenesis in which the creation of disorder in the DNA structure is an obligatory process followed by the selection of ordered, mutated DNA forms that ultimately give rise to metastases.

Recognition principle of the TAP transporter disclosed by combinatorial peptide libraries

Transport of peptides across the membrane of the endoplasmic reticulum for assembly with MHC class I molecules is an essential step in antigen presentation to cytotoxic T cells. This task is performed by the major histocompatibility complex-encoded transporter associated with antigen processing (TAP). Using a combinatorial approach we have analyzed the substrate specificity of human TAP at high resolution and in the absence of any given sequence context, revealing the contribution of each peptide residue in stabilizing binding to TAP. Human TAP was found to be highly selective with peptide affinities covering at least three orders of magnitude. Interestingly, the selectivity is not equally distributed over the substrate. Only the N-terminal three positions and the C-terminal residue are critical, whereas effects from other peptide positions are negligible. A major influence from the peptide backbone was uncovered by peptide scans and libraries containing d amino acids. Again, independent of peptide length, critical positions were clustered near the peptide termini. These approaches demonstrate that human TAP is selective, with residues determining the affinity located in distinct regions, and point to the role of the peptide backbone in binding to TAP. This binding mode of TAP has implications in an optimized repertoire selection and in a coevolution with the major histocompatibility complex/T cell receptor complex.

The selectivity of the hair cell’s mechanoelectrical-transduction channel promotes Ca2+ flux at low Ca2+ concentrations

The mechanoelectrical-transduction channel of the hair cell is permeable to both monovalent and divalent cations. Because Ca2+ entering through the transduction channel serves as a feedback signal in the adaptation process that sets the channel’s open probability, an understanding of adaptation requires estimation of the magnitude of Ca2+ influx. To determine the Ca2+ current through the transduction channel, we measured extracellular receptor currents with transepithelial voltage-clamp recordings while the apical surface of a saccular macula was bathed with solutions containing various concentrations of K+, Na+, or Ca2+. For modest concentrations of a single permeant cation, Ca2+ carried much more receptor current than did either K+ or Na+. For higher cation concentrations, however, the flux of Na+ or K+ through the transduction channel exceeded that of Ca2+. For mixtures of Ca2+ and monovalent cations, the receptor current displayed an anomalous mole-fraction effect, which indicates that ions interact while traversing the channel’s pore. These results demonstrate not only that the hair cell’s transduction channel is selective for Ca2+ over monovalent cations but also that Ca2+ carries substantial current even at low Ca2+ concentrations. At physiological cation concentrations, Ca2+ flux through transduction channels can change the local Ca2+ concentration in stereocilia in a range relevant for the control of adaptation.

Kinesin’s processivity results from mechanical and chemical coordination between the ATP hydrolysis cycles of the two motor domains

Kinesin is a processive motor protein: A single molecule can walk continuously along a microtubule for several micrometers, taking hundreds of 8-nm steps without dissociating. To elucidate the biochemical and structural basis for processivity, we have engineered a heterodimeric one-headed kinesin and compared its biochemical properties to those of the wild-type two-headed molecule. Our construct retains the functionally important neck and tail domains and supports motility in high-density microtubule gliding assays, though it fails to move at the single-molecule level. We find that the ATPase rate of one-headed kinesin is 3–6 s−1 and that detachment from the microtubule occurs at a similar rate (3 s−1). This establishes that one-headed kinesin usually detaches once per ATP hydrolysis cycle. Furthermore, we identify the rate-limiting step in the one-headed hydrolysis cycle as detachment from the microtubule in the ADP⋅Pi state. Because the ATPase and detachment rates are roughly an order of magnitude lower than the corresponding rates for two-headed kinesin, the detachment of one head in the homodimer (in the ADP⋅Pi state) must be accelerated by the other head. We hypothesize that this results from internal strain generated when the second head binds. This idea accords with a hand-over-hand model for processivity in which the release of the trailing head is contingent on the binding of the forward head. These new results, together with previously published ones, allow us to propose a pathway that defines the chemical and mechanical cycle for two-headed kinesin.

Fractionation factors and activation energies for exchange of the low barrier hydrogen bonding proton in peptidyl trifluoromethyl ketone complexes of chymotrypsin

NMR investigations have been carried out of complexes between bovine chymotrypsin Aα and a series of four peptidyl trifluoromethyl ketones, listed here in order of increasing affinity for chymotrypsin: N-Acetyl-l-Phe-CF3, N-Acetyl-Gly-l-Phe-CF3, N-Acetyl-l-Val-l-Phe-CF3, and N-Acetyl-l-Leu-l-Phe-CF3. The D/H fractionation factors (φ) for the hydrogen in the H-bond between His 57 and Asp 102 (His 57-Hδ1) in these four complexes at 5°C were in the range φ = 0.32–0.43, expected for a low-barrier hydrogen bond. For this series of complexes, measurements also were made of the chemical shifts of His 57-Hɛ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 8.97–9.18), the exchange rate of the His 57-Hδ1 proton with bulk water protons (284–12.4 s−1), and the activation enthalpies for this hydrogen exchange (14.7–19.4 kcal⋅mol−1). It was found that the previously noted correlations between the inhibition constants (Ki 170–1.2 μM) and the chemical shifts of His 57-Hδ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 18.61–18.95) for this series of peptidyl trifluoromethyl ketones with chymotrypsin [Lin, J., Cassidy, C. S. & Frey, P. A. (1998) Biochemistry 37, 11940–11948] could be extended to include the fractionation factors, hydrogen exchange rates, and hydrogen exchange activation enthalpies. The results support the proposal of low barrier hydrogen bond-facilitated general base catalysis in the addition of Ser 195 to the peptidyl carbonyl group of substrates in the mechanism of chymotrypsin-catalyzed peptide hydrolysis. Trends in the enthalpies for hydrogen exchange and the fractionation factors are consistent with a strong, double-minimum or single-well potential hydrogen bond in the strongest complexes. The lifetimes of His 57-Hδ1, which is solvent shielded in these complexes, track the strength of the hydrogen bond. Because these lifetimes are orders of magnitude shorter than those of the complexes themselves, the enzyme must have a pathway for hydrogen exchange at this site that is independent of dissociation of the complexes.

Characteristics of virus-specific CD8+ T cells in the liver during the control and resolution phases of influenza pneumonia

Dissection of the primary and secondary response to an influenza A virus established that the liver contains a substantial population of CD8+ T cells specific for the immunodominant epitope formed by H-2Db and the influenza virus nucleoprotein peptide fragment NP366–374 (DbNP366). The numbers of CD8+ DbNP366+ cells in the liver reflected the magnitude of the inflammatory process in the pneumonic lung, though replication of this influenza virus is limited to the respiratory tract. Analysis of surface phenotypes indicated that the liver CD8+ DbNP366+ cells tended to be more “activated” than the set recovered from lymphoid tissue but generally less so than those from the lung. The distinguishing characteristic of the lymphocytes from the liver was that the prevalence of the CD8+ DbNP366+ set was always much higher than the percentage of CD8+ T cells that could be induced to synthesize interferon γ after short-term, in vitro stimulation with the NP366–374 peptide, whereas these values were generally comparable for virus-specific CD8+ T cells recovered from other tissue sites. Also, the numbers of apoptotic CD8+ T cells were higher in the liver. The results overall are consistent with the idea that antigen-specific CD8+ T cells are destroyed in the liver during the control and resolution phases of this viral infection, though this destruction is not necessarily an immediate process.

The degree of nonlinearity and anisotropy of blood vessel elasticity

Blood vessel elasticity is important to physiology and clinical problems involving surgery, angioplasty, tissue remodeling, and tissue engineering. Nonlinearity in blood vessel elasticity in vivo is important to the formation of solitons in arterial pulse waves. It is well known that the stress–strain relationship of the blood vessel is nonlinear in general, but a controversy exists on how nonlinear it is in the physiological range. Another controversy is whether the vessel wall is biaxially isotropic. New data on canine aorta were obtained from a biaxial testing machine over a large range of finite strains referred to the zero-stress state. A new pseudo strain energy function is used to examine these questions critically. The stress–strain relationship derived from this function represents the sum of a linear stress–strain relationship and a definitely nonlinear relationship. This relationship fits the experimental data very well. With this strain energy function, we can define a parameter called the degree of nonlinearity, which represents the fraction of the nonlinear strain energy in the total strain energy per unit volume. We found that for the canine aorta, the degree of nonlinearity varies from 5% to 30%, depending on the magnitude of the strains in the physiological range. In the case of canine pulmonary artery in the arch region, Debes and Fung [Debes, J. C. & Fung, Y. C.(1995) Am. J. Physiol. 269, H433–H442] have shown that the linear regime of the stress–strain relationship extends from the zero-stress state to the homeostatic state and beyond. Both vessels, however, are anisotropic in both the linear and nonlinear regimes.

Cell surface expression of the epithelial Na channel and a mutant causing Liddle syndrome: A quantitative approach

The epithelial amiloride-sensitive sodium channel (ENaC) controls transepithelial Na+ movement in Na+-transporting epithelia and is associated with Liddle syndrome, an autosomal dominant form of salt-sensitive hypertension. Detailed analysis of ENaC channel properties and the functional consequences of mutations causing Liddle syndrome has been, so far, limited by lack of a method allowing specific and quantitative detection of cell-surface-expressed ENaC. We have developed a quantitative assay based on the binding of 125I-labeled M2 anti-FLAG monoclonal antibody (M2Ab*) directed against a FLAG reporter epitope introduced in the extracellular loop of each of the α, β, and γ ENaC subunits. Insertion of the FLAG epitope into ENaC sequences did not change its functional and pharmacological properties. The binding specificity and affinity (Kd = 3 nM) allowed us to correlate in individual Xenopus oocytes the macroscopic amiloride-sensitive sodium current (INa) with the number of ENaC wild-type and mutant subunits expressed at the cell surface. These experiments demonstrate that: (i) only heteromultimeric channels made of α, β, and γ ENaC subunits are maximally and efficiently expressed at the cell surface; (ii) the overall ENaC open probability is one order of magnitude lower than previously observed in single-channel recordings; (iii) the mutation causing Liddle syndrome (β R564stop) enhances channel activity by two mechanisms, i.e., by increasing ENaC cell surface expression and by changing channel open probability. This quantitative approach provides new insights on the molecular mechanisms underlying one form of salt-sensitive hypertension.

Autoinhibition of serotonin cells: An intrinsic regulatory mechanism sensitive to the pattern of usage of the cells

After periods of high-frequency firing, the normal rhythmically active serotonin (5HT)-containing neurosecretory neurons of the lobster ventral nerve cord display a period of suppressed spike generation and reduced synaptic input that we refer to as “autoinhibition.” The duration of this autoinhibition is directly related to the magnitude and duration of the current injection triggering the high-frequency firing. More interesting, however, is that the autoinhibition is inversely related to the initial firing frequency of these cells within their normal range of firing (0.5–3 Hz). This allows more active 5HT neurons to resume firing after shorter durations of inhibition than cells that initially fired at slower rates. Although superfused 5HT inhibits the spontaneous firing of these cells, the persistence of autoinhibition in saline with no added calcium, in cadmium-containing saline, and in lobsters depleted of serotonin suggests that intrinsic membrane properties account for the autoinhibition. A similar autoinhibition is seen in spontaneously active octopamine neurons but is absent from spontaneously active γ-aminobutyric acid cells. Thus, this might be a characteristic feature of amine-containing neurosecretory neurons. The 5HT cells of vertebrate brain nuclei share similarities in firing frequencies, spike shapes, and inhibition by 5HT with the lobster cells that were the focus of this study. However, the mechanism suggested to underlie autoinhibition in vertebrate neurons is that 5HT released from activated or neighboring cells acts back on inhibitory autoreceptors that are found on the dendrites and cell bodies of these neurons.

Identification of the proton pathway in bacterial reaction centers: Both protons associated with reduction of QB to QBH2 share a common entry point

The reaction center from Rhodobacter sphaeroides uses light energy for the reduction and protonation of a quinone molecule, QB. This process involves the transfer of two protons from the aqueous solution to the protein-bound QB molecule. The second proton, H+(2), is supplied to QB by Glu-L212, an internal residue protonated in response to formation of QA− and QB−. In this work, the pathway for H+(2) to Glu-L212 was studied by measuring the effects of divalent metal ion binding on the protonation of Glu-L212, which was assayed by two types of processes. One was proton uptake from solution after the one-electron reduction of QA (DQA→D+QA−) and QB (DQB→D+QB−), studied by using pH-sensitive dyes. The other was the electron transfer kAB(1) (QA−QB→QAQB−). At pH 8.5, binding of Zn2+, Cd2+, or Ni2+ reduced the rates of proton uptake upon QA− and QB− formation as well as kAB(1) by ≈an order of magnitude, resulting in similar final values, indicating that there is a common rate-limiting step. Because D+QA− is formed 105-fold faster than the induced proton uptake, the observed rate decrease must be caused by an inhibition of the proton transfer. The Glu-L212→Gln mutant reaction centers displayed greatly reduced amplitudes of proton uptake and exhibited no changes in rates of proton uptake or electron transfer upon Zn2+ binding. Therefore, metal binding specifically decreased the rate of proton transfer to Glu-L212, because the observed rates were decreased only when proton uptake by Glu-L212 was required. The entry point for the second proton H+(2) was thus identified to be the same as for the first proton H+(1), close to the metal binding region Asp-H124, His-H126, and His-H128.

An evaluation of the performance of cDNA microarrays for detecting changes in global mRNA expression

The cDNA microarray is one technological approach that has the potential to accurately measure changes in global mRNA expression levels. We report an assessment of an optimized cDNA microarray platform to generate accurate, precise and reliable data consistent with the objective of using microarrays as an acquisition platform to populate gene expression databases. The study design consisted of two independent evaluations with 70 arrays from two different manufactured lots and used three human tissue sources as samples: placenta, brain and heart. Overall signal response was linear over three orders of magnitude and the sensitivity for any element was estimated to be 2 pg mRNA. The calculated coefficient of variation for differential expression for all non-differentiated elements was 12–14% across the entire signal range and did not vary with array batch or tissue source. The minimum detectable fold change for differential expression was 1.4. Accuracy, in terms of bias (observed minus expected differential expression ratio), was less than 1 part in 10 000 for all non-differentiated elements. The results presented in this report demonstrate the reproducible performance of the cDNA microarray technology platform and the methods provide a useful framework for evaluating other technologies that monitor changes in global mRNA expression.

Strong DNA binding by covalently linked dimeric Lac headpiece: Evidence for the crucial role of the hinge helices

The combined structural and biochemical studies on Lac repressor bound to operator DNA have demonstrated the central role of the hinge helices in operator bending and the induction mechanism. We have constructed a covalently linked dimeric Lac-headpiece that binds DNA with four orders of magnitude higher affinity as compared with the monomeric form. This enabled a detailed biochemical and structural study of Lac binding to its cognate wild-type and selected DNA operators. The results indicate a profound contribution of hinge helices to the stability of the protein–DNA complex and highlight their central role in operator recognition. Furthermore, protein–DNA interactions in the minor groove appear to modulate hinge helix stability, thus accounting for affinity differences and protein-induced DNA bending among the various operator sites. Interestingly, the in vitro DNA-binding affinity of the reported dimeric Lac construct can de readily modulated by simple adjustment of redox conditions, thus rendering it a potential artificial gene regulator.

Very large radio surveys of the sky

Recent advances in electronics and computing have made possible a new generation of large radio surveys of the sky that yield an order-of-magnitude higher sensitivity and positional accuracy. Combined with the unique properties of the radio universe, these quantitative improvements open up qualitatively different and exciting new scientific applications of radio surveys.