Residual delay maps unveil global patterns of atmospheric nonlinearity and produce improved local forecasts

We use residual-delay maps of observational field data for barometric pressure to demonstrate the structure of latitudinal gradients in nonlinearity in the atmosphere. Nonlinearity is weak and largely lacking in tropical and subtropical sites and increases rapidly into the temperate regions where the time series also appear to be much noisier. The degree of nonlinearity closely follows the meridional variation of midlatitude storm track frequency. We extract the specific functional form of this nonlinearity, a V shape in the lagged residuals that appears to be a basic feature of midlatitude synoptic weather systems associated with frontal passages. We present evidence that this form arises from the relative time scales of high-pressure versus low-pressure events. Finally, we show that this nonlinear feature is weaker in a well regarded numerical forecast model (European Centre for Medium-Range Forecasts) because small-scale temporal and spatial variation is smoothed out in the grided inputs. This is significant, in that it allows us to demonstrate how application of statistical corrections based on the residual-delay map may provide marked increases in local forecast accuracy, especially for severe weather systems.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

Congruent Docking of Dimeric Kinesin and ncd into Three-dimensional Electron Cryomicroscopy Maps of Microtubule–Motor ADP Complexes

We present a new map showing dimeric kinesin bound to microtubules in the presence of ADP that was obtained by electron cryomicroscopy and image reconstruction. The directly bound monomer (first head) shows a different conformation from one in the more tightly bound empty state. This change in the first head is amplified as a movement of the second (tethered) head, which tilts upward. The atomic coordinates of kinesin·ADP dock into our map so that the tethered head associates with the bound head as in the kinesin dimer structure seen by x-ray crystallography. The new docking orientation avoids problems associated with previous predictions; it puts residues implicated by proteolysis-protection and mutagenesis studies near the microtubule but does not lead to steric interference between the coiled-coil tail and the microtubule surface. The observed conformational changes in the tightly bound states would probably bring some important residues closer to tubulin. As expected from the homology with kinesin, the atomic coordinates of nonclaret disjunctional protein (ncd)·ADP dock in the same orientation into the attached head in a map of microtubules decorated with dimeric ncd·ADP. Our results support the idea that the observed direct interaction between the two heads is important at some stages of the mechanism by which kinesin moves processively along microtubules.

Coexistence of linear zones and pinwheels within orientation maps in cat visual cortex

Revealing the layout of cortical maps is important both for understanding the processes involved in their development and for uncovering the mechanisms underlying neural computation. The typical organization of orientation maps in the cat visual cortex is radial; complete orientation cycles are mapped around orientation singularities. In contrast, long linear zones of orientation representation have been detected in the primary visual cortex of the tree shrew. In this study, we searched for the existence of long linear sequences and wide linear zones within orientation preference maps of the cat visual cortex. Optical imaging based on intrinsic signals was used. Long linear sequences and wide linear zones of preferred orientation were occasionally detected along the border between areas 17 and 18, as well as within area 18. Adjacent zones of distinct radial and linear organizations were observed across area 18 of a single hemisphere. However, radial and linear organizations were not necessarily segregated; long (7.5 mm) linear sequences of preferred orientation were found embedded within a typical pinwheel-like organization of orientation. We conclude that, although the radial organization is dominant, perfectly linear organization may develop and perform the processing related to orientation in the cat visual cortex.

A novel modifier gene for plasma von Willebrand factor level maps to distal mouse chromosome 11

Type 1 von Willebrand disease (VWD), characterized by reduced levels of plasma von Willebrand factor (VWF), is the most common inherited bleeding disorder in humans. Penetrance of VWD is incomplete, and expression of the bleeding phenotype is highly variable. In addition, plasma VWF levels vary widely among normal individuals. To identify genes that influence VWF level, we analyzed a genetic cross between RIIIS/J and CASA/Rk, two strains of mice that exhibit a 20-fold difference in plasma VWF level. DNA samples from F2 progeny demonstrating either extremely high or extremely low plasma VWF levels were pooled and genotyped for 41 markers spanning the autosomal genome. A novel locus accounting for 63% of the total variance in VWF level was mapped to distal mouse chromosome 11, which is distinct from the murine Vwf locus on chromosome 6. We designated this locus Mvwf for “modifier of VWF.” Additional genotyping of as many as 2407 meioses established a high resolution genetic map with gene order Cola1-Itg3a-Ngfr-Mvwf/Gip-Hoxb9-Hoxb1-Cbx·rs2-Cox5a-Gfap. The Mvwf candidate interval between Ngfr and Hoxb9 is ≈0.5 centimorgan (cM). These results demonstrate that a single dominant gene accounts for the low VWF phenotype of RIIIS/J mice in crosses with several other strains. The pattern of inheritance suggests a gain-of-function mutation in a unique component of VWF biosynthesis or processing. Characterization of the human homologue for Mvwf may have relevance for a subset of type 1 VWD cases and may define an important genetic factor modifying penetrance and expression of mutations at the VWF locus.

Unexpected crucial role of residue 225 in serine proteases

Residue 225 in serine proteases of the chymotrypsin family is Pro or Tyr in more than 95% of nearly 300 available sequences. Proteases with Y225 (like some blood coagulation and complement factors) are almost exclusively found in vertebrates, whereas proteases with P225 (like degradative enzymes) are present from bacteria to human. Saturation mutagenesis of Y225 in thrombin shows that residue 225 affects ligand recognition up to 60,000-fold. With the exception of Tyr and Phe, all residues are associated with comparable or greatly reduced catalytic activity relative to Pro. The crystal structures of three mutants that differ widely in catalytic activity (Y225F, Y225P, and Y225I) show that although residue 225 makes no contact with substrate, it drastically influences the shape of the water channel around the primary specificity site. The activity profiles obtained for thrombin also suggest that the conversion of Pro to Tyr or Phe documented in the vertebrates occurred through Ser and was driven by a significant gain (up to 50-fold) in catalytic activity. In fact, Ser and Phe are documented in 4% of serine proteases, which together with Pro and Tyr account for almost the entire distribution of residues at position 225. The unexpected crucial role of residue 225 in serine proteases explains the evolutionary selection of residues at this position and shows that the structural determinants of protease activity and specificity are more complex than currently believed. These findings have broad implications in the rational design of enzymes with enhanced catalytic properties.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

NMR structures of three single-residue variants of the human prion protein

The NMR structures of three single-amino acid variants of the C-terminal domain of the human prion protein, hPrP(121–230), are presented. In hPrP(M166V) and hPrP(R220K) the substitution is with the corresponding residue in murine PrP, and in hPrP(S170N) it is with the corresponding Syrian hamster residue. All three substitutions are in the surface region of the structure of the cellular form of PrP (PrPC) that is formed by the C-terminal part of helix 3, with residues 218–230, and a loop of residues 166–172. This molecular region shows high species variability and has been implicated in specific interactions with a so far not further characterized “protein X,” and it is related to the species barrier for transmission of prion diseases. As expected, the three variant hPrP(121–230) structures have the same global architecture as the previously determined wild-type bovine, human, murine, and Syrian hamster prion proteins, but with the present study two localized “conformational markers” could be related with single amino acid exchanges. These are the length and quality of definition of helix 3, and the NMR-observability of the residues in the loop 166–172. Poor definition of the C-terminal part of helix 3 is characteristic for murine PrP and has now been observed also for hPrP(R220K), and NMR observation of the complete loop 166–172 has so far been unique for Syrian hamster PrP and is now also documented for hPrP(S170N).

Roles of a conserved arginine residue of DsbB in linking protein disulfide-bond-formation pathway to the respiratory chain of Escherichia coli

The active-site cysteines of DsbA, the periplasmic disulfide-bond-forming enzyme of Escherichia coli, are kept oxidized by the cytoplasmic membrane protein DsbB. DsbB, in turn, is oxidized by two kinds of quinones (ubiquinone for aerobic and menaquinone for anaerobic growth) in the electron-transport chain. We describe the isolation of dsbB missense mutations that change a highly conserved arginine residue at position 48 to histidine or cysteine. In these mutants, DsbB functions reasonably well aerobically but poorly anaerobically. Consistent with this conditional phenotype, purified R48H exhibits very low activity with menaquinone and an apparent Michaelis constant (Km) for ubiquinone seven times greater than that of the wild-type DsbB, while keeping an apparent Km for DsbA similar to that of wild-type enzyme. From these results, we propose that this highly conserved arginine residue of DsbB plays an important role in the catalysis of disulfide bond formation through its role in the interaction of DsbB with quinones.

In vitro roles of invariant helix–turn–helix motif residue R383 in σ54 (σN)

In vitro DNA-binding and transcription properties of σ54 proteins with the invariant Arg383 in the putative helix–turn–helix motif of the DNA-binding domain substituted by lysine or alanine are described. We show that R383 contributes to maintaining stable holoenzyme–promoter complexes in which limited DNA opening downstream of the –12 GC element has occurred. Unlike wild-type σ54, holoenzymes assembled with the R383A or R383K mutants could not form activator-independent, heparin-stable complexes on heteroduplex Sinorhizobium meliloti nifH DNA mismatched next to the GC. Using longer sequences of heteroduplex DNA, heparin-stable complexes formed with the R383K and, to a lesser extent, R383A mutant holoenzymes, but only when the activator and a hydrolysable nucleotide was added and the DNA was opened to include the –1 site. Although R383 appears inessential for polymerase isomerisation, it makes a significant contribution to maintaining the holoenzyme in a stable complex when melting is initiating next to the GC element. Strikingly, Cys383-tethered FeBABE footprinting of promoter DNA strongly suggests that R383 is not proximal to promoter DNA in the closed complex. This indicates that R383 is not part of the regulatory centre in the σ54 holoenzyme, which includes the –12 promoter region elements. R383 contributes to several properties, including core RNA polymerase binding and to the in vivo stability of σ54.

Formation of temporal-feature maps by axonal propagation of synaptic learning

Computational maps are of central importance to a neuronal representation of the outside world. In a map, neighboring neurons respond to similar sensory features. A well studied example is the computational map of interaural time differences (ITDs), which is essential to sound localization in a variety of species and allows resolution of ITDs of the order of 10 μs. Nevertheless, it is unclear how such an orderly representation of temporal features arises. We address this problem by modeling the ontogenetic development of an ITD map in the laminar nucleus of the barn owl. We show how the owl's ITD map can emerge from a combined action of homosynaptic spike-based Hebbian learning and its propagation along the presynaptic axon. In spike-based Hebbian learning, synaptic strengths are modified according to the timing of pre- and postsynaptic action potentials. In unspecific axonal learning, a synapse's modification gives rise to a factor that propagates along the presynaptic axon and affects the properties of synapses at neighboring neurons. Our results indicate that both Hebbian learning and its presynaptic propagation are necessary for map formation in the laminar nucleus, but the latter can be orders of magnitude weaker than the former. We argue that the algorithm is important for the formation of computational maps, when, in particular, time plays a key role.

Control of gating mode by a single amino acid residue in transmembrane segment IS3 of the N-type Ca2+ channel

N-type Ca2+ channels can be inhibited by neurotransmitter-induced release of G protein βγ subunits. Two isoforms of Cav2.2 α1 subunits of N-type calcium channels from rat brain (Cav2.2a and Cav2.2b; initially termed rbB-I and rbB-II) have different functional properties. Unmodulated Cav2.2b channels are in an easily activated “willing” (W) state with fast activation kinetics and no prepulse facilitation. Activating G proteins shifts Cav2.2b channels to a difficult to activate “reluctant” (R) state with slow activation kinetics; they can be returned to the W state by strong depolarization resulting in prepulse facilitation. This contrasts with Cav2.2a channels, which are tonically in the R state and exhibit strong prepulse facilitation. Activating or inhibiting G proteins has no effect. Thus, the R state of Cav2.2a and its reversal by prepulse facilitation are intrinsic to the channel and independent of G protein modulation. Mutating G177 in segment IS3 of Cav2.2b to E as in Cav2.2a converts Cav2.2b tonically to the R state, insensitive to further G protein modulation. The converse substitution in Cav2.2a, E177G, converts it to the W state and restores G protein modulation. We propose that negatively charged E177 in IS3 interacts with a positive charge in the IS4 voltage sensor when the channel is closed and produces the R state of Cav2.2a by a voltage sensor-trapping mechanism. G protein βγ subunits may produce reluctant channels by a similar molecular mechanism.

An essential amino acid residue in the protein translocation channel revealed by targeted random mutagenesis of SecY

The SecY/Sec61α family of membrane proteins are the central subunits of the putative protein translocation channel. We introduced random mutations into a segment of Escherichia coli SecY within its cytoplasmic domain 5, which was shown previously to be important for the SecA-dependent translocation activity. Mutations were classified into those retaining function and those gaining a dominant-interfering ability caused by a loss of function. These analyses showed that Arg-357, Pro-358, Gly-359, and Thr-362 are functionally important; Arg-357, conserved in almost all organisms, was identified as an indispensable residue.

Growth factors regulate phototransduction in retinal rods by modulating cyclic nucleotide-gated channels through dephosphorylation of a specific tyrosine residue

Illumination of vertebrate rod photoreceptors leads to a decrease in the cytoplasmic cGMP concentration and closure of cyclic nucleotide-gated (CNG) channels. Except for Ca2+, which plays a negative feedback role in adaptation, and 11-cis-retinal, supplied by the retinal pigment epithelium, all of the biochemical machinery of phototransduction is thought to be contained within rod outer segments without involvement of extrinsic regulatory molecules. Here we show that insulin-like growth factor-I (IGF-I), a paracrine factor released from the retinal pigment epithelium, alters phototransduction by rapidly increasing the cGMP sensitivity of CNG channels. The IGF-I-signaling pathway ultimately involves a protein tyrosine phosphatase that catalyzes dephosphorylation of a specific residue in the α-subunit of the rod CNG channel protein. IGF-I conjointly accelerates the kinetics and increases the amplitude of the light response, distinct from events that accompany adaptation. These effects of IGF-I could result from the enhancement of the cGMP sensitivity of CNG channels. Hence, in addition to long-term control of development and survival of rods, growth factors regulate phototransduction in the short term by modulating CNG channels.

Imaging neuroscience: Principles or maps?

This article reviews some recent trends in imaging neuroscience. A distinction is made between making maps of functional responses in the brain and discerning the rules or principles that underlie their organization. After considering developments in the characterization of brain imaging data, several examples are presented that highlight the context-sensitive nature of neuronal responses that we measure. These contexts can be endogenous and physiological, reflecting the fact that each cortical area, or neuronal population, expresses its dynamics in the context of interactions with other areas. Conversely, these contexts can be experimental or psychological and can have a profound effect on the regional effects elicited. In this review we consider experimental designs and analytic strategies that go beyond cognitive subtraction and speculate on how functional imaging can be used to address both the details and principles underlying functional integration and specialization in the brain.