Single amino acid substitutions determine mouse CD8 allotypes: epitope mapping of mouse CD8.
Explore the source record for details and available documents.
Biomedical subjects
Publications and source records attributed to A Sarai.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
In the DNA-binding domain of the c-myb protooncogene product (c-Myb) which consists of three repeats of 51-52 amino acids, there are 3 perfectly conserved tryptophans in each repeat. Site-directed mutagenesis of these tryptophans showed that any single or multiple mutations of tryptophan to hydrophilic residues or alanine abolished or greatly reduced the sequence-specific DNA-binding activity, but mutations to hydrophobic amino acids retained considerable activity. Raman spectroscopic study showed that these tryptophans were buried in the protein core. These 3 tryptophans are proposed to form a cluster in the hydrophobic core in each repeat. This hypothetical structure is referred to as the "tryptophan cluster," and it may represent a characteristic property of a group of DNA-binding proteins including the myb- and ets-related proteins.
The role of solvation on the sequence dependent conformational variabilities in DNA has been studied by calculating hydration free energies from solvent accessible surface areas for several base steps, as a function of various helical parameters, roll, twist and propeller twist. The results of roll calculations suggest opposite trends for AA and GG steps, with the former tending to have a compressed minor groove and the latter a compressed major groove. These trends are consistent with the experimental findings on sequence preferences and the nature of anisotropic bending of DNA observed in nucleosomes (Drew, H.R. and Travers, A.A., J. Mol. Biol. 186, 773-790 (1985); Satchwell, S.C., Drew, H.R. and Travers, A.A., J. Mol. Biol. 191, 659-675 (1986)) and CAP-DNA interactions (Gartenberg, M.R. and Crothers, D.M., Nature 333, 824-829, (1988)). Solvation energy profiles also indicate preferences for the base pairs in GG and AA steps to adopt low and high propeller twists, respectively. Such agreements may either reflect a coincidence of solvation effects with other energy terms or a dominance of solvent effects. The results are discussed in the context of the crystallographic observations of structural tendencies.
By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.
Molecular recognition is considered as an information transduction process, and information theory is applied to analyze the efficiency of the information transduction during molecular recognition. It is shown that a certain amount of information is gained in the recognition process, and it can be expressed in terms of thermodynamic quantities in the binding process of molecules. The information gain, which can be interpreted as the amount of information extracted from the sequence or structure by a molecule, represents the intrinsic ability of a molecule to recognize specific sequence or structure out of large ambient ensemble based on physical interactions. In contrast with thermodynamic quantities themselves, the information gain is a normalized quantity, and thus serves as a good measure of specificity in the molecular recognition. The information gain, which can be evaluated experimentally, enables us to identify the specific interactions and compare the specificity among different molecules quantitatively. In contrast with macroscopic recognition, specificity of molecular recognition may be severely restricted by the thermodynamic nature of molecular interactions. The thermodynamic meaning of the information gain and its biological implications to molecular recognition are discussed.
The distributions of the junction sequences of homooligomer tracts of various lengths have been examined in prokaryotic DNA sequences and compared with those of eukaryotes. The general trends in the nearest and next to nearest neighbors to the tracts are similar for both groups. In both prokaryotes and eukaryotes A/T runs are preferentially flanked on either the 5' or the 3' ends by A and/or T. G/C runs are preferentially flanked by G and/or C. There is discrimination against A/T runs flanked by G or C and G/C runs flanked by A or T. However, whereas the distribution of prokaryotic homooligomer tract junction sequences was quite homogeneous, large variations were observed in the 5-fold larger eukaryotic database, increasing in magnitude from tracts of length 2 to 3 to 4 base pairs long. Possible DNA conformational implications and in particular DNA curvature and packaging aspects of prokaryotes and eukaryotes are discussed.
We have studied, by conformational analysis, the sequence dependence of DNA conformational transition between B- and A-forms. We have considered intramolecular interactions between base pairs, without backbone, to examine their role in the conformational transition between B- and A-forms, and found that base pairs themselves usually have intrinsic conformational preferences for the B- or A-form. Calculation of all ten possible base steps shows that the base combinations, CC (or GG), GC, AT, and TA, have tendencies to assume the A-conformation. Results show that it is particularly easy to slide along the long axis of the base pair for these steps, with AT and CC showing especially flat energies. These calculations show that a preference for the B- or A-conformation depends on the electrostatic energy parameters, in particular, on dielectric and shielding constants; the A-conformation is preferred for low dielectric constant or low shielding. Both the A- and B-conformations are mainly stabilized by electrostatic interactions between favorably juxtaposed atomic charges on base pairs; however, the B-conformation generally has more favorable van der Waals interactions than the A-form. These sequence-dependent conformational preference and environmental effects agree roughly with experimental observations, suggesting that the origin of the conformational polymorphism is attributable to the intrinsic conformational preference of base pairs.
Results of systematic base-substitution experiments suggest that the lambda repressor dimer, made of identical subunits, recognizes the "pseudo(2-fold)symmetric" operator sequence asymmetrically. Base substitutions within the consensus half of the operator affect binding more than base substitutions within the nonconsensus half of the operator. Furthermore, changing the nonconsensus base pairs to the consensus base pairs does not increase, but decreases, binding. Evidently, the two subunits of the lambda repressor dimer bind to the two halves of the operator differently. This is consistent with the recently determined crystal structure of the complex, which shows that the relative positioning of the amino acids to the DNA bases are slightly different in the two halves of the operator. The sequence-specific interactions indicated by the systematic base-substitution experiments correlate well with the locations of the specific contacts found in the complex. Thus, the amino acids of lambda repressor, mainly of alpha 3-helix and the N-terminus arm, seem to directly read-out the DNA sequence by forming specific hydrogen bonds and hydrophobic contacts to the DNA bases. The observed asymmetric recognition suggests that no recognition code governs amino acids and DNA bases in protein-DNA interactions.
We measured quantitatively the binding affinities of purified Cro repressor to the chemically synthesized wild-type and mutant OR1 operators, consisting of all three possible base-pair substitutions and of thymine to uracil substitutions at each base-pair position of the 17-base-pair operator sequence. The sequence-specific interactions between Cro repressor and the operator DNA occur at the symmetrically disposed outer 7-base-pair positions of each half operator and at the central base-pair position. The binding of Cro is almost symmetrical with respect to the pseudo-twofold symmetry of the binding site. The binding free energy changes calculated from the affinity changes are mostly additive for specific Cro binding. Also the binding affinities of Cro to the operators or any other DNA sequences can be predicted by simple addition of free energy changes of single base substitutions. We isolated cro mutants by site-directed mutagenesis and studied their DNA binding to the wild-type and base-substituted mutant operators. The sequence-specific contacts derived from such studies are significantly different from the models proposed by Ohlendorf et al. [Ohlendorf, D. H., Anderson, W. F., Takeda, Y. & Matthews, B. W. (1982) Nature (London) 298, 719-723] and by Hochschild et al. [Hochschild, A., Douhan, J., III, & Ptashne, M. (1986) Cell 47, 807-816].
Previous studies of the dinucleotides flanking both the 5' and 3' ends of homooligomer tracts have shown that some flanks are consistently preferred over others (1,2). In the first preferred group, the homooligomer tracts are flanked by the same nucleotide and/or the complementary nucleotides, e.g.,ATAn,TTAn,CCGn, where n = 2-5. Runs flanked by nucleotides with which they cannot base pair are distinctly disfavored. (In this group An/Tn are flanked by C and/or G; Gn/Cn are flanked by A/T, e.g.,CGAn,TnGG,GnAT). The frequencies of runs flanked by A or T, and G or C ("mixed"group) are as expected. Here we seek the origin of this effect and its relevance to protein-DNA interactions. Surprisingly, within the first group, runs flanked by their complements with a pyrimidine-purine junction (e.g.,TTAn,CnGG) are greatly preferred. The frequencies of their purine-pyrimidine junction mirror-images is just as expected. This effect, as well as additional ones enumerated below, is seen universally in eukaryotes and in prokaryotes, although it is stronger in the former. Detailed analysis of regulatory regions shows these strong trends, particularly in GC sequences. The potential relationship to DNA conformation and DNA-protein interaction is discussed.
Conformational analysis of DNA shows that the origin of the B-form double helix can be attributed in large part to the atomic charge pattern in the base pairs. The charge patterns favor specific helical stacking of the base pairs. Base pairs alone--without backbones--have a strong tendency to form helix, indicating that the backbones play a rather passive role in determining the basic helical structure of DNA. It is mainly the electrostatic interactions determined by the charge pattern on base pairs that stabilize a particular helical conformation. The charge pattern in the base pairs appears to be responsible for much of the sequence dependence of DNA conformation, rather than steric clashes.
Recent studies of homooligomer tracts suggest different characteristics from random sequence DNA (dA).(dT) and (dG).(dC) tracts are frequent in upstream regions and in some cases have been shown to be essential for regulation. Here we examine homooligomer occurrences in non-coding and coding eukaryotic sequences, focusing on the context in which the homooligomers occur. This analysis of sequences in the junction areas yields distinct and consistent characteristics. In particular, the nucleotide interrupting a run is most frequently complementary to the run. The base next to it is most frequently identical to the one constituting the run. For A or T runs the least frequent nearest and next to nearest neighbors are G or C. For G or C tracts the least frequent are A or T. Complementary oligomers behave similarly. These and additional trends are strongest for run lengths greater than or equal to 3. The computations are carried out on the whole eukaryotic database of greater than 4 x 10(6) nucleotides, separately for coding and non-coding regions. These same trends are evident for both groups, but are somewhat stronger for the non-coding regions. The context in which the homooligomers occur may yield some clues to DNA conformation and its biological implications.
We propose some specific DNA conformations that explain, in terms of molecular conformations, the anomalous gel electrophoretic behavior of the sequences (VA4T4X), and (V2A3T3X2)i where V and X are either G or C. Previously (J. Biomole. Struct. Dyn. 4, 41, 1986) we considered hydrophobic interactions among aliphatic hydrocarbon groups in A/T sequences. In the sequences (T)n.(A)n, the T's are slightly bent to yield structures with tightly stacked methyl groups along one side of the major groove. By folding together the two pairs of stacked methyls on the opposite sides of the major groove. TTAA might yield a relatively sharp bend. On this basis, we show below that the sequences (VT4A4X)i might form a very tightly coiled super-helix whereas the sequences (VA4T4X)i form a broad super-helix of radius approximately 120 A for i = 25. The sequence (V2A3T3X2)i forms a slightly smaller radius super-helix. The time of passage through the gel has been taken to be inversely proportional to the smallest dimension of the molecule. Specifically we are taking the ratio of the apparent molecular weight to the actual molecular weight to be related to the moment of inertia I1 about the smallest principal axis of the molecular conformation. We find a good fit to the experimental gel mobility data of Hagerman (2) if we assume this ratio to be proportional to (I1)1/5.
We have considered hydrophobic interactions among aliphatic hydrocarbon groups in A/T sequences. The slightly overwound sequences (T)n.(A)n yield structures with tightly stacked methyl groups along one side of the major groove. The sequence TTAA may yield a sharp bend by folding together the two pairs of stacked methyls on the opposite sides of the major groove. Thus the sequence can affect the formation of either a smooth bend or a sharp kink. These sequence dependent local conformations may be related to a number of biological results.
The B-to-Z transition in supercoiled circular DNA is modeled as a strain-induced nonlinear excitation process. Using a model, in which DNA is regarded as a chain of units with a bistable energy function along the twisting coordinate together with a harmonic inter-unit interaction, we show that a Z region and the accompanying two B-Z junctions of finite width appear naturally as a solution of nonlinear equations, when the strain exceeds a critical value. We examine the B-Z transition behaviour as a function of twist under various situations. We also analyse available experimental results on B-Z transition in supercoiled plasmid with G-C insertions by this mechanistic model in order to estimate the magnitude of model parameters. The energy barrier of the B-Z transition is estimated to be of the order of 1 kcal/mole per base pair. The analysis shows that if the length of the insertion is less than a certain value, the entire insertion converts to Z form at a transition point, but if the insertion is much longer, the B-Z transition exhibits a different behavior, in which part of the insertion flips to Z form and the Z region expands linearly upon changing linking number.
Explore the source record for details and available documents.
The rate of the electron-transfer reaction between bacteriopheophytin and the first quinone in isolated reaction centers of Rhodopseudomonas sphaeroides has an unusual temperature dependence. The rate increases about threefold with decreasing temperature between 300 and 25 K, and decreases abruptly at temperatures below 25 K. Partial deuteration of the reaction centers alters the temperature dependence of the rate constant. Qualitative features of the temperature dependence can be understood in the context of a theory of nonadiabatic electron transfer (Sarai, 1980. Biochim. Biophys. Acta 589:71-83). We conclude that very low-energy (10-50 cm-1) processes, perhaps skeletal vibrations of the protein, are important to electron transfer. Higher-energy vibrations, possibly involving the pyrrolic N--H bonds of bacteriopheophytin, also are important in this process.