Small viscosity asymptotics for the inertial range of local structure and for the wall region of wall-bounded turbulent shear flow.

The small viscosity asymptotics of the inertial range of local structure and of the wall region in wallbounded turbulent shear flow are compared. The comparison leads to a sharpening of the dichotomy between Reynolds number dependent scaling (power-type) laws and the universal Reynolds number independent logarithmic law in wall turbulence. It further leads to a quantitative prediction of an essential difference between them, which is confirmed by the results of a recent experimental investigation. These results lend support to recent work on the zero viscosity limit of the inertial range in turbulence.

Global properties of the mapping between local amino acid sequence and local structure in proteins.

Local protein structure prediction efforts have consistently failed to exceed approximately 70% accuracy. We characterize the degeneracy of the mapping from local sequence to local structure responsible for this failure by investigating the extent to which similar sequence segments found in different proteins adopt similar three-dimensional structures. Sequence segments 3-15 residues in length from 154 different protein families are partitioned into neighborhoods containing segments with similar sequences using cluster analysis. The consistency of the sequence-to-structure mapping is assessed by comparing the local structures adopted by sequence segments in the same neighborhood in proteins of known structure. In the 154 families, 45% and 28% of the positions occur in neighborhoods in which one and two local structures predominate, respectively. The sequence patterns that characterize the neighborhoods in the first class probably include virtually all of the short sequence motifs in proteins that consistently occur in a particular local structure. These patterns, many of which occur in transitions between secondary structural elements, are an interesting combination of previously studied and novel motifs. The identification of sequence patterns that consistently occur in one or a small number of local structures in proteins should contribute to the prediction of protein structure from sequence.

Assembly of τ protein into Alzheimer paired helical filaments depends on a local sequence motif (306VQIVYK311) forming β structure

We have searched for a minimal interaction motif in τ protein that supports the aggregation into Alzheimer-like paired helical filaments. Digestion of the repeat domain with different proteases yields a GluC-induced fragment comprising 43 residues (termed PHF43), which represents the third repeat of τ plus some flanking residues. This fragment self assembles readily into thin filaments without a paired helical appearance, but these filaments are highly competent to nucleate bona fide PHFs from full-length τ. Probing the interactions of PHF43 with overlapping peptides derived from the full τ sequence yields a minimal hexapeptide interaction motif of 306VQIVYK311 at the beginning of the third internal repeat. This motif coincides with the highest predicted β-structure potential in τ. CD and Fourier transform infrared spectroscopy shows that PHF43 acquires pronounced β structure in conditions of self assembly. Point mutations in the hexapeptide region by proline-scanning mutagenesis prevent the aggregation. The data indicate that PHF assembly is initiated by a short fragment containing the minimal interaction motif forming a local β structure embedded in a largely random-coil protein.

Structure and function in rhodopsin: Topology of the C-terminal polypeptide chain in relation to the cytoplasmic loops*

Cysteine mutagenesis and site-directed spin labeling in the C-terminal region of rhodopsin have been used to probe the local structure and proximity of that region to the cytoplasmic loops. Each of the native amino acids in the sequence T335–T340 was replaced with Cys, one at a time. The sulfhydryl groups of all mutants reacted rapidly with the sulfhydryl reagent 4,4′-dithiodipyridine, which indicated a high degree of solvent accessibility. Furthermore, to probe the proximity relationships, a series of double Cys mutants was constructed. One Cys in all sets was at position 338 and the other was at a position in the sequence S240–V250 in the EF interhelical loop, at position 65 in the AB interhelical loop, or at position 140 in the CD interhelical loop. In the dark state, no significant disulfide formation was observed between C338 and C65 or C140 under the conditions used, whereas a relatively rapid disulfide formation was observed between C338 and C242 or C245. Spin labels in the double Cys mutants showed the strongest magnetic interactions between the nitroxides attached to C338 and C245 or C246. Light activation of the double mutant T242C/S338C resulted in slower disulfide formation, whereas interactions between nitroxides at C338 and C245 or C246 decreased. These results suggest the proximity of the C-terminal residue C338 to residues located on the outer face of a cytoplasmic helical extension of the F helix with an apparent increase of distance upon photoactivation.

NMR solution structure of the human prion protein

The NMR structures of the recombinant human prion protein, hPrP(23–230), and two C-terminal fragments, hPrP(90–230) and hPrP(121–230), include a globular domain extending from residues 125–228, for which a detailed structure was obtained, and an N-terminal flexibly disordered “tail.” The globular domain contains three α-helices comprising the residues 144–154, 173–194, and 200–228 and a short anti-parallel β-sheet comprising the residues 128–131 and 161–164. Within the globular domain, three polypeptide segments show increased structural disorder: i.e., a loop of residues 167–171, the residues 187–194 at the end of helix 2, and the residues 219–228 in the C-terminal part of helix 3. The local conformational state of the polypeptide segments 187–193 in helix 2 and 219–226 in helix 3 is measurably influenced by the length of the N-terminal tail, with the helical states being most highly populated in hPrP(23–230). When compared with the previously reported structures of the murine and Syrian hamster prion proteins, the length of helix 3 coincides more closely with that in the Syrian hamster protein whereas the disordered loop 167–171 is shared with murine PrP. These species variations of local structure are in a surface area of the cellular form of PrP that has previously been implicated in intermolecular interactions related both to the species barrier for infectious transmission of prion disease and to immune reactions.

The folding mechanism of larger model proteins: role of native structure.

The folding mechanism of a 125-bead heteropolymer model for proteins is investigated with Monte Carlo simulations on a cubic lattice. Sequences that do and do not fold in a reasonable time are compared. The overall folding behavior is found to be more complex than that of models for smaller proteins. Folding begins with a rapid collapse followed by a slow search through the semi-compact globule for a sequence-dependent stable core with about 30 out of 176 native contacts which serves as the transition state for folding to a near-native structure. Efficient search for the core is dependent on structural features of the native state. Sequences that fold have large amounts of stable, cooperative structure that is accessible through short-range initiation sites, such as those in anti-parallel sheets connected by turns. Before folding is completed, the system can encounter a second bottleneck, involving the condensation and rearrangement of surface residues. Overly stable local structure of the surface residues slows this stage of the folding process. The relation of the results from the 125-mer model studies to the folding of real proteins is discussed.

Chromophore-assisted light inactivation and self-organization of microtubules and motors

Chromophore-assisted light inactivation (CALI) offers the only method capable of modulating specific protein activities in localized regions and at particular times. Here, we generalize CALI so that it can be applied to a wider range of tasks. Specifically, we show that CALI can work with a genetically inserted epitope tag; we investigate the effectiveness of alternative dyes, especially fluorescein, comparing them with the standard CALI dye, malachite green; and we study the relative efficiencies of pulsed and continuous-wave illumination. We then use fluorescein-labeled hemagglutinin antibody fragments, together with relatively low-power continuous-wave illumination to examine the effectiveness of CALI targeted to kinesin. We show that CALI can destroy kinesin activity in at least two ways: it can either result in the apparent loss of motor activity, or it can cause irreversible attachment of the kinesin enzyme to its microtubule substrate. Finally, we apply this implementation of CALI to an in vitro system of motor proteins and microtubules that is capable of self-organized aster formation. In this system, CALI can effectively perturb local structure formation by blocking or reducing the degree of aster formation in chosen regions of the sample, without influencing structure formation elsewhere.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25–30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras–p120GAP–GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH2 group of the Gln-61 and the γ-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for P—O bond cleavage.

A single G-to-C change causes human centromere TGGAA repeats to fold back into hairpins.

Recently, we established that satellite III (TGGAA)n tandem repeats, which occur at the centromeres of human chromosomes, pair with themselves to form an unusual "self-complementary" antiparallel duplex containing (GGA)2 motifs in which two unpaired guanines from opposite strands intercalate between sheared G.A base pairs. In separate studies, we have also established that the GCA triplet does not form bimolecular (GCA)2 motifs but instead promotes the formation of hairpins containing a GCA-turn motif in which the loop contains a single cytidine closed by a sheared G.A pair. Since TGCAA is the most frequent variant of TGGAA found in satellite III repeats, we reasoned that the potential of this variant to form GCA-turn miniloop fold-back structures might be an important factor in modulating the local structure in natural (TGGAA)n repeats. We report here the NMR-derived solution structure of the heptadecadeoxynucleotide (G)TGGAATGCAATGGAA(C) in which a central TGCAA pentamer is flanked by two TGGAA pentamers. This 17-mer forms a rather unusual and very stable hairpin structure containing eight base pairs in the stem, only four of which are Watson-Crick pairs, and a loop consisting of a single cytidine residue. The stem contains a (GGA)2 motif with intercalative 14G/4G stacking between two sheared G.A base pairs; the loop end of the stem consists of a sheared 8G.10A closing pair with the cytosine base of the 9C loop stacked on 8G. The remarkable stability of this unusual hairpin structure (Tm = 63 degrees C) suggests that it probably plays an important role in modulating the folding of satellite III (TGGAA)n repeats at the centromere.

Initiation sites of protein folding by NMR analysis.

Detailed characterization of denatured states of proteins is necessary to understand the interactions that funnel the large number of possible conformations along fast routes for folding. Nuclear magnetic resonance experiments based on the nuclear Overhauser effect (NOE) detect hydrogen atoms close in space and provide information about local structure. Here we present an NMR procedure that detects almost all sequential NOEs between amide hydrogen atoms (HN-HN NOE), including those in random coil regions in a protein, barnase, in urea solutions. A semi-quantitative analysis of these HN-HN NOEs identified partly structured regions that are in remarkable agreement with those found to form early on the reaction pathway. Our results strongly suggest that the folding of barnase initiates at the first helix and the beta-turn between the third and the fourth strands. This strategy of defining residual structure has also worked for cold-denatured barstar and guanidinium hydrochloride-denatured chymotrypsin inhibitor 2 and so should be generally applicable.

On the existence and implications of an inverse folding code in proteins.

The existence of a code relating the set of possible sequences at a given position in a protein backbone to the local structure at that location is investigated. It is shown that only 73% of 4-C alpha structure fragments in a sample of 114 protein structures exhibit a preference for a particular set of sequences. The remaining structures can accommodate essentially any sequence. The structures that encode specific sequence distributions include the classical "secondary" structures, with the notable exception of planar (beta) bends. It is suggested that this has implications as to the mechanism of folding in proteins with extensive sheet/barrel structure. The possible role of structures that do not encode specific sequences as mutation hot spots is noted.

Crystal structure of 2,5-diketo-d-gluconic acid reductase A complexed with NADPH at 2.1-Å resolution

The three-dimensional structure of Corynebacterium 2,5-diketo-d-gluconic acid reductase A (2,5-DKGR A; EC 1.1.1.-), in complex with cofactor NADPH, has been solved by using x-ray crystallographic data to 2.1-Å resolution. This enzyme catalyzes stereospecific reduction of 2,5-diketo-d-gluconate (2,5-DKG) to 2-keto-l-gulonate. Thus the three-dimensional structure has now been solved for a prokaryotic example of the aldo–keto reductase superfamily. The details of the binding of the NADPH cofactor help to explain why 2,5-DKGR exhibits lower binding affinity for cofactor than the related human aldose reductase does. Furthermore, changes in the local loop structure near the cofactor suggest that 2,5-DKGR will not exhibit the biphasic cofactor binding characteristics observed in aldose reductase. Although the crystal structure does not include substrate, the two ordered water molecules present within the substrate-binding pocket are postulated to provide positional landmarks for the substrate 5-keto and 4-hydroxyl groups. The structural basis for several previously described active-site mutants of 2,5-DKGR A is also proposed. Recent research efforts have described a novel approach to the synthesis of l-ascorbate (vitamin C) by using a genetically engineered microorganism that is capable of synthesizing 2,5-DKG from glucose and subsequently is transformed with the gene for 2,5-DKGR. These modifications create a microorganism capable of direct production of 2-keto-l-gulonate from d-glucose, and the gulonate can subsequently be converted into vitamin C. In economic terms, vitamin C is the single most important specialty chemical manufactured in the world. Understanding the structural determinants of specificity, catalysis, and stability for 2,5-DKGR A is of substantial commercial interest.

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.

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.

A physical basis for protein secondary structure

A physical theory of protein secondary structure is proposed and tested by performing exceedingly simple Monte Carlo simulations. In essence, secondary structure propensities are predominantly a consequence of two competing local effects, one favoring hydrogen bond formation in helices and turns, the other opposing the attendant reduction in sidechain conformational entropy on helix and turn formation. These sequence specific biases are densely dispersed throughout the unfolded polypeptide chain, where they serve to preorganize the folding process and largely, but imperfectly, anticipate the native secondary structure.