Sequence dependence of DNA conformational flexibility

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.     

Distinct patterns in homooligomer tract sequence context in prokaryotic and eukaryotic DNA

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.     

Equilibrium folding and unfolding pathways for a model protein

The folding–unfolding process of reduced bovine pancreatic trypsin inhibitor was investigated with an idealized model employing approximate free energies. The protein is regarded to consist of only C^α and C^β atoms. The backbone dihedral angles are the only conformational variables and are permitted to take discrete values at every 10°. Intraresidue energies consist of two terms: an empirical part taken from the observed frequency distributions of (?,ψ) and an additional favorable energy assigned to the native conformation of each residue. Interresidue interactions are simplified by assuming that there is an attractive energy operative only between residue pairs in close contact in the native structure. A total of 230,000 molecular conformations, with no atomic overlaps, ranging from the native state to the denatured state, are randomly generated by changing the sampling bias. Each conformation is classified according to its conformational energy, F; a conformational entropy, S(F) is estimated for each value of F from the number of samples. The dependence of S(F) on energy reveals that the folding–unfolding transition for this idealized model is an “all-or-none” type; this is attributable to the specific long-range interactions. Interresidue contact probabilities, averaged over samples representing various stages of folding, serve to characterize folding intermediates. Most probable equilibrium pathways for the folding–unfolding transition are constructed by connecting conformationally similar intermediates. The specific details obtained for bovine pancreatic trypsin inhibitor are as follows: (1) Folding begins with the appearance of nativelike medium-range contacts at a β-turn and at the α-helix. (2) These grow to include the native pair of interacting β-strands. This state includes intact regular secondary conformations, as well as the interstrand sheet contacts, and corresponds to an activated state with the highest free energy on the pathway. (3) Additional native long-range contacts are completely formed either toward the amino terminus or toward the carboxyl terminus. (4) In a final step, the missing contacts appear. Although these folding pathways for this model are not consistent with experimental reports, it does indicate multiple folding pathways. The method is general and can be applied to any set of calculated conformational energies and furthermore permits investigation of gross folding features.     

Conformations of folded proteins in restricted spaces

A new method is presented to examine the complete range of folded topologies accessible in the compact state of globular proteins. The procedure is to generate all conformations, with volume exclusion, upon a lattice in a space restricted to the individual protein's known compact conformational space. Using one lattice point per residue, we find 10(2)-10(4) possible compact conformations for the five small globular proteins studied. Subsequently, these conformations are evaluated in terms of residue-specific, pairwise contact energies that favor nonbonded, hydrophobic interactions. Native structures for the five proteins are always found within the best 2% of all conformers generated. This novel method is simple and general and can be used to determine a small group of most favorable overall arrangements for the folding of specific amino acid sequences within a restricted space.     

Solvent effect on binding thermodynamics of biopolymers

The indirect solvent-induced effect on the free energy of binding of biopolymers is examined within the framework of classical statistical mechanics. We focus specifically on the role of the solute-solvent hydrogen bonding. In particular, we have estimated the first order solvent effect on the indirect interaction between two biopolymers. We find that the solvent-induced interactions between two hydrophilic groups through water-bridged hydrogen bonds could significantly enhance the binding free energy. Some preliminary estimates indicate that this effect is significant and perhaps could be crucial in molecular recognition processes. Furthermore, we have calculated, from crystal structure data, the distance distribution between all the oxygens and nitrogens on the surface of some proteins that do not belong to the binding domain. In most cases we found an enhanced peak in the range of 4-5 A, which is where we expect to find strong solvent-induced interactions.     

Strong Patterns in Homooligomer Tracts Occurrences in Non-Coding and in Potential Regulatory Sites in Eukaryotic Genomes

Abstract

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., ATA_n, TTA_n, CCG_n, where n=2–5. Runs flanked by nucleotides with which they cannot base pair are distinctly disfavored. (In this group A/T_n are flanked by C and/or G; G_n/C_n are flanked by A/T, e.g., CGA_n, T_nGG, G., AT). The frequencies of runs flanked by AorT, 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., TTA_n, C_nGG) 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.     

164 Extracting dynamics information from multiple structures

Meaningful dynamics information can be extracted from multiple experimental structures of the same, or closely related, proteins or RNAs. The covariance matrix of atom positions is decomposable into its principal components, and in this way, it is possible to rank-order the changes in the set of structures, and to determine what the most significant changes are. Usually, only a few principal components dominate the motions of the structures, and these usually relate to the functional dynamics. This dynamics information provides strong evidence for the plasticity of protein and RNA structures, and also suggests that these structures almost always have a highly limited repertoire of motions. In some cases, such as HIV protease, the dominant motions are opening and closing over the active site. For myoglobin, the changes are much smaller, reflecting in part the small changes in sequence, but nonetheless they show characteristic details that depend on the species. Sets of structures can also be used to derive the effective microscopic forces that are forcing a given conformational transition.     

Relationship Between Curved DNA Conformations and Slow Gel Migration

Abstract

We propose some specific DNA conformations that explain, in terms of molecular conformations, the anomalous gel electrophoretic behavior of the sequences (VA₄T₄X)₁, and (V₂A₃X₂)₁ where V and X are either G or C. Previously (J. Biomole. Struct. Dyn. 4, 41, 1986) we considered hydrophobic interactions a mong 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 (VT₄A₄X)₁ might form a very tightly coiled super-helix whereas the sequences (VA₄T₄X)₁ form a broad super-helix of radius ~ 120 A for i = 25. The sequence (V₂A₃T₃X₂)₁ forms a slightly smaller radius super-helix. The time of passage through the gel has been taken to be inversely proportional to the smallesuiimension 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 I₁ 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 (I₁)^1/5.     

Dynamic Allostery Mediated by a Conserved Tryptophan in the Tec Family Kinases

Bruton’s tyrosine kinase (Btk) is a Tec family non-receptor tyrosine kinase that plays a critical role in immune signaling and is associated with the immunological disorder X-linked agammaglobulinemia (XLA). Our previous findings showed that the Tec kinases are allosterically activated by the adjacent N-terminal linker. A single tryptophan residue in the N-terminal 17-residue linker mediates allosteric activation, and its mutation to alanine leads to the complete loss of activity. Guided by hydrogen/deuterium exchange mass spectrometry results, we have employed Molecular Dynamics simulations, Principal Component Analysis, Community Analysis and measures of node centrality to understand the details of how a single tryptophan mediates allostery in Btk. A specific tryptophan side chain rotamer promotes the functional dynamic allostery by inducing coordinated motions that spread across the kinase domain. Either a shift in the rotamer population, or a loss of the tryptophan side chain by mutation, drastically changes the coordinated motions and dynamically isolates catalytically important regions of the kinase domain. This work also identifies a new set of residues in the Btk kinase domain with high node centrality values indicating their importance in transmission of dynamics essential for kinase activation. Structurally, these node residues appear in both lobes of the kinase domain. In the N-lobe, high centrality residues wrap around the ATP binding pocket connecting previously described Catalytic-spine residues. In the C-lobe, two high centrality node residues connect the base of the R- and C-spines on the αF-helix. We suggest that the bridging residues that connect the catalytic and regulatory architecture within the kinase domain may be a crucial element in transmitting information about regulatory spine assembly to the catalytic machinery of the catalytic spine and active site.