Persistence analysis of the static and dynamical helix deformations of DNA oligonucleotides: Application to the crystal structure and molecular dynamics simulation of d(CGCGAATTCGCG)2

A theory and graphical presentation for the analysis of helix structure and deformations in oligonucleotides is presented. The parameters “persistence” and “flexibility” as defined in the configurational statistics of polymers of infinite length are reformulated at the oligonucleotide level in an extension of J. A. Schellman's method [(1974) Biopolymers, Vol. 17, pp. 217–226], and used as a basis for a systematic “Persistence Analysis” of the helix deformation properties for all possible subsequences in the structure. The basis for the analysis is a set of link vectors referenced to individual base pairs, and is limited to sequences exhibiting only perturbed rod-like behavior, i.e., below the threshold for supercoiling. The present application of the method is concerned with a physical model for the angular component of bending, so the link vectors are defined as the unit components of a global helix axis obtained by the procedure “Curves” of R. Lavery and H. Sklenar [(1988) J. Biomol. Struct. Dynam., Vol. 6, pp. 63–91; (1989) J. Biomol. Struct. Dynam., Vol. 6, pp. 655–667]. A discussion, of the relationship between global bending and relative orientation of base pairs is provided. Our approach is illustrated by analysis of some model oligonucleotide structures with intrinsic kinks, the crystal structure of the dodecamer d (CGCGAATTCGCG)₂, and the results of two molecular dynamics simulations on this dodecamer using two variations of the GROMOS force field. The results indicate that essentially all aspects of curvature in short oligonucleotides can be determined, such as the position and orientation of each bend, the sharpness or smoothness, and the location and linearity of subsequences. In the case of molecular dynamics simulations, where a Boltzmann ensemble of structures is analyzed, the spatial extent of the deformations (flexibility) is also considered. © 1993 John Wiley & Sons, Inc.     

A database of 32 DNA triplets to study triple helices by molecular mechanics and dynamics

We present here a database of 32 deoxyribonucleotide triplets, that can be used as building blocks of triple helix forming deoxyribonucleotides on a computer. This database is made of all the pairing schemes of the triplets ATT, GCC⁺, ATA and GCG where the third base forms two hydrogen bonds with the purine of the first two Watson-Crick strands. The essential features of the known triple helices were preserved in the resulting structures. A triple helix can be easily built from any combination of these basic triplets. Four homogeneous and alternate triple helices thus obtained were studied by molecular mechanics and dynamics in vacuo. The results are in agreement with known experimental observations for ATT and suggest a possible structure for the GCG triple helix. In order to characterize the geometry of the structures obtained, the definitions of nucleic acid structure parameters (R.E. Dickerson et al., EMBO J. 8 (1989) 1–4) have been extended to triple helical polynucleotides.     

Interactions of molecules with nucleic acids. I. An algorithm to generate nucleic acid structures with an application to the B-DNA structure and a counterclockwise helix.

An algorithm is developed that enables the routine determination of backbone conformations of nucleic acids. All atomic positions including hydrogen are specified in accord with experimental bond lengths and angles but with theoretically determined conformational angles. For two Watson-Crick base pairs at a separation of 3.38 Å, and perpendicular to a common helical axis, minimum energy configurations are found for all 10 combinations at helical angles of α ～ 36°–38°, corresponding to the B-DNA structure with C(2′)-endo sugar puckers. Backbone configurations exist only within the range 35.5° ? α ? 42°, which suggests the origin of the 10-fold helix. Calculated stacking energies for the B-DNA structure increases for each of the clustered groups of base pairs: G·C with G·C, G·C with A·T, and A·T with A·T, and they are in approximate agreement with experimental observations. The counter-clockwise helix is examined, and physically meaningful structures are found only when the helical axes of successive base pairs are disjointed.     

Comparison of analyses of DNA curvature

Popular programs for characterizing DNA structure include Curves 5.1 (Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 63-91, 1988; Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 655-67, 1989) and Freehelix98 (Dickerson, R. E., Nucleic Acids Res. 26, 1906-1926, 1998), along with the more recent 3DNA (X. J. Lu, Z. Shakked and W. K. Olson., J. Mol. Biol. 300, 819-840 (2000). Given input of structural coordinates, all of these programs return values of the local helical parameters, such as roll, tilt, twist, etc. The first two programs also provide characterization of global curvature. Madbend (Strahs, D. and Schlick, T., J. Mol. Biol. 301, 643-663, 2000), a program that computes global curvature from local roll, tilt, and twist parameters, can be applied to the output of all three structural programs. We have compared the curvature predicted by the three programs with and without the use of Madbend. Global bend magnitudes and directions as well as values of helical kinks were calculated for four high-resolution DNA structures and four model DNA helices. Global curvature determined by Curves 5.1 without Madbend was found to differ from values obtained using Freehelix98 with or without Madbend or 3DNA and Curves 5.1 with Madbend. Using model helices, this difference was attributed the fact that Curves 5.1 is the only program sensitive to changes in axial displacement, such as shift and slide. Madbend produced robust values of bend magnitude and direction, and displayed little sensitivity to axis displacement or the source of local helical parameters. Madbend also appears to be the method of choice for bending comparisons of high-resolution structures with results from cyclization kinetics, a method that measures DNA curvature as a vectorial sum of local roll and tilt angles.     

Length does matter for cGAS

Recognition of foreign nucleic acids by the immune system is essential to host protection against many viral and bacterial infections. It relies on the capacity of innate immune sensors to selectively distinguish self‐ and non‐self‐nucleic acids, on the basis of a variety of parameters including base modifications, sequence composition, length or subcellular localisation. In this issue of EMBO Reports, Luecke et al 1 describe that the sensing of cytoplasmic double‐stranded DNA by the cyclic GMP–AMP (cGAMP) synthase (cGAS) is much more sensitive for longer fragments, when low doses of cytoplasmic DNA are used. This finding repositions length as the predominant factor governing the discrimination between self‐ and non‐self‐cytoplasmic DNA.     

Flexibility of 3′,5′ Deoxyribonucleooside Diphosphates

Abstract

In 3′,5′ deoxyribonucleoside diphosphates, in addition to the nature of the base and the sugar puckering, there are six single bond rotations. However, from the analysis of crystal structure data on the constituents of nucleic acids, only three rotational angles, that are about glycosyl bond, about C4′-C5′ and about C3′-O3′ bonds, are flexible. For a given sugar puckering and a base, potential energy calculations using non-bonded, electrostatic and torsional functions were carried out by varying the three torsion angles. The energies are represented as isopotential energy surfaces. Since the availability of the real-time color graphics, it is possible to analyse these isopotential energy surfaces. The calculations were carried out for C3′ exo and C3′ endo puckerings for deoxyribose and also for four bases. These calculations throw more light not only on the allowed regions for the three rotational angles but also on the relationships among them. The dependence of base and the puckering of the sugar on these rotational angles and thereby the flexibility of the 3′,5′ deoxyribonucleoside diphosphates is discussed. From our calculations, it is now possible to follow minimum energy path for interconversion among various conformers.     

Mutagen–nucleic acid complexes at the polynucleotide duplex level in solution: Intercalation of proflavine into poly(dA-dT) and the melting transition of the complex

The nmr chemical shifts and line widths of the nucleic acid base and sugar proton resonances and the proflavine ring protons can be monitored through the melting transition of the proflavine + poly(dA-dT) complex, phosphate/dye (P/D) ratio = 24 and 8 in 1M salt solution. The nucleic acid and mutagen protons in the complex are in fast exchange between duplex and strand states with the midpoint of the melting transition monitored at the nucleic acid resonances increasing from 72.6°C for poly(dA-dT) to 78.1°C for the P/D = 24 complex and 83.4°C for the P/D = 8 complex in 1M salt solution. The melting transition monitored by the proflavine resonances were 80.0°C for the P/D = 24 complex and 84.3°C for the P/D = 8 complex in 1M salt solution. Since the nucleic acid is in excess at high P/D ratios, the nucleic acid transitions are an average for the opening of mutagen-free and mutagen-bound base-pair regions, while the proflavine transitions monitor the melting of mutagen-bound base-pair regions. The observed 0.75 to 0.95 ppm unfield shift at all four proflavine protons on formation of the complex with poly(dA-dT) provides direct evidence for intercalation of the mutagen between base pairs of the nucleic acid duplex. We have deduced the approximate overlap geometry between the proflavine ring and nearest-neighbor base pairs at the intercalation site from a comparison between experimental proflavine complexation shifts and those calculated for various stacking orientations. The experimental chemical shift of the poly(dA-dT) adenine H-2 resonance in the duplex state in the absence and presence of proflavine suggests that intercalation occurs preferentially at dT-dA sites. The selective chemical shift changes at the sugar H-2′,2″ and H-3′ resonances of the poly(dA-dT) duplex on complex formation demonstrates changes in the sugar pucker and/or torsion angles of the sugar phosphate backbone at the intercalation site.     

Early spring warming as one of the factors responsible for expansion of aquatic fern Salvinia natans (L.) All. in the Vistula delta (south Baltic Sea coast)

From the mid-19th century to the end of the 20th century, Salvinia natans (L.) All. occurred very rarely in the Vistula Delta (northern Poland), but from the beginning of the 21st century it was present in almost every watercourse and had formed very abundant populations. We examined the influence of temperature on the abundance of this plant and the efficiency of macrospore germination. Field work was carried out in 10 permanent plots every 14 days for 5 years. Macrospores germinate at water temperature of 12.4 ± 0.2°C or higher; at 20°C they develop more effectively than at 15°C. Usually, ice cover on the rivers melts in the second half of March. At this time, macro- and microspores emerge on the water surface and germinate in April. They occur in the water surface film at 15.1 ± 2.4°C and massively die during spring frost. After 1989, March and April mean temperature in the Vistula Delta rose by 1.6°C versus 1901–1988, and by 1.9°C versus 1851–1988. In 1951–1988, the mean temperature for March and April was +4.6°C and was characterized by considerable interannual variation (SD = 1.64), whereas in 1989–2009, it rose to +5.7°C and the variation range narrowed (SD = 1.24). We found that macrospores are active earlier during the warm and mild summers, germination is more effective, survival of young stages is higher, the growing season is longer, and the number of vegetative offsprings in a year is larger.     

The Opening of a Single Base Without Perturbations of Neighboring Nucleotides: A Study on Crystal B-DNA duplex d(CGCGAATTCGCG)2

Abstract

In this work we explore the possibility of the opening of a single base without perturbation of its neighboring nucleotides. Low energy base opening into the grooves can be accomplished by rotation of the relevant backbone and glycosidic bond torsion angles. The pathway has been determined by identifying ζ torsion angle as the reaction coordinate together with the accompanying geometric requirement that guides the displacement of other torsion angles. Our study on Dickerson dodecamer duplex d(CGCGAATTCGCG)₂ showed that all bases with normal equilibrium ζ can be rotated by ~ 30°, corresponding to ~ 3.5Å base displacement, towards the major groove. Such an opening extent is comparable with estimated amplitudes of local angular motions in DNA bases from NMR experiments, which might facilitate proton exchange. The computed base opening energy barrier is also comparable with measured base pair opening enthalpy. These results indicate possible relevance of the pathway studied in this work with experimentally observed base pair opening process. Our analysis also showed a preference for base opening along the major groove and an abnormal behavior for bases with unusual equilibrium ζ torsion angle.     

TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts

NMR chemical shifts in proteins depend strongly on local structure. The program TALOS establishes an empirical relation between ¹³C, ¹⁵N and ¹H chemical shifts and backbone torsion angles ϕ and ψ (Cornilescu et al. J Biomol NMR 13 289–302, 1999). Extension of the original 20-protein database to 200 proteins increased the fraction of residues for which backbone angles could be predicted from 65 to 74%, while reducing the error rate from 3 to 2.5%. Addition of a two-layer neural network filter to the database fragment selection process forms the basis for a new program, TALOS+, which further enhances the prediction rate to 88.5%, without increasing the error rate. Excluding the 2.5% of residues for which TALOS+ makes predictions that strongly differ from those observed in the crystalline state, the accuracy of predicted ϕ and ψ angles, equals ±13°. Large discrepancies between predictions and crystal structures are primarily limited to loop regions, and for the few cases where multiple X-ray structures are available such residues are often found in different states in the different structures. The TALOS+ output includes predictions for individual residues with missing chemical shifts, and the neural network component of the program also predicts secondary structure with good accuracy. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Biphasic helix–coil transition of the ethidium bromide·poly(dA-dT) and the propidium diiodide·poly(dA-dT) Complexes. Stabilization of base-pair regions centered about the intercalation site

The nonexchangeable base and sugar proton nmr resonances and the 260 and 278-nm uv-absorbance bands of the nucleic acid were utilized to monitor the temperature-dependent duplex-to-strand transition of the alternating purine–pyrimidine deoxyribopolynucleotide poly(dA-dT) in the absence and presence of ethidium bromide (EB) at phosphate/drug = 50, 28, and 15 and propidium diiodide (PI) at P/D = 50, 25, 15, 10, and 5 in 0.1 M salt between 50° and 100°C. The nmr and optical methods monitor a biphasic duplex-to strand transition for the drug–poly(dA-dT) complexes. We have monitored the dissociation of the drug from the complex at the ethidium bromide phenanthridine ring and side-chain proton nmr resonances and the propidium diiodide 494 and 535-nm uv-absorbance bands and demonstrate that dissociation of the drug corresponds to the higher temperature transition in the biphasic nucleic acid melting curves. The lower temperature cooperative transition is assigned to the opening of drug-free AT base-pair regions in the drug–poly(dA-dT) complex and exhibits an increase in transition midpoint and a decrease in cooperativity with increasing drug concentration. The higher temperature cooperative transition is assigned to the opening of AT base-pair regions centered about the bound drug in the complex and exhibits an increase in the transition midpoint on raising the drug concentration. The large upfield shifts of the phenanthridine ring (but not side chain) protons of ethidium bromide on complex formation demonstrate intercalation of the drug between base pairs of the poly(dA-dT) duplex. The nucleic acid base and sugar resonances of poly(dA-dT) in 0.1 M phosphate undergo chemical shift changes between 0° and 50°C indicative of premelting conformational transition(s).     

Conformational Studies of Nucleic Acids: IV. The Conformational Energetics of Oligonucleotides: d(ApApApA) and ApApApA

Abstract

Utilizing a new method for modeling furanose pseudorotation (D. A Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: δ (C5′-C4′-C3′-03′)χ (04′-C1′-N9- C4), where the energetic basis for an approximately linear δ-χ correlation can be seen; ζ (C3′- 03′-P-05′)-δ, in which the experimentally observed linear correlation between ζ and δ in DNA (220° < ζ <280°) is clearly predicted; ζ-ε (C4′-C3′-03′-P), which shows that e increases with decreasing ζ <260°; α (03′-P-05′-C5′)-γ (05′-C5′-C4′-C3′) where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion Ô is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in a and y. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.     

Cell surface hydrophobicity of colistin-susceptible vs resistant Acinetobacter baumannii determined by contact angles: methodological considerations and implications

Contact angle analysis of cell surface hydrophobicity (CSH) describes the tendency of a water droplet to spread across a lawn of filtered bacterial cells. Colistin‐induced disruption of the Gram‐negative outer membrane necessitates hydrophobic contacts with lipopolysaccharide (LPS). We aimed to characterize the CSH of Acinetobacter baumannii using contact angles, to provide insight into the mechanism of colistin resistance. Contact angles were analysed for five paired colistin‐susceptible and resistant Ac. baumannii strains. Drainage of the water droplet through bacterial layers was demonstrated to influence results. Consequently, measurements were performed 0·66 s after droplet deposition. Colistin‐resistant cells exhibited lower contact angles (38·8±2·8–46·8±1·3°) compared with their paired colistin‐susceptible strains (40·7±3·0–48·0±1·4°; anova ; P < 0·05). Contact angles increased at stationary phase (50·3±2·9–61·5±2·5° and 47·4±2·0–50·8±3·2°, susceptible and resistant, respectively, anova ; P < 0·05) and in response to colistin 32 mg l^?1 exposure (44·5±1·5–50·6±2·8° and 43·5±2·2–48·0±2·2°, susceptible and resistant, respectively; anova ; P < 0·05). Analysis of complemented strains constructed with an intact lpxA gene, or empty vector, highlighted the contribution of LPS to CSH. Compositional outer‐membrane variations likely account for CSH differences between Ac. baumannii phenotypes, which influence the hydrophobic colistin–bacterium interaction. Important insight into the mechanism of colistin resistance has been provided. Greater consideration of contact angle methodology is necessary to ensure accurate analyses are performed.     

Geometry of experimentally observed RNA residues in tRNA and dinucleoside monophosphates: the effect of small variations in the backbone angles.

The polymerization of various experimentally observed conformers of RNA from tRNA and some dinucleoside monophosphates have been examined with a program that computes the basic helix parameters directly from the six backbone torsion angles ω′, ?′, ψ′, ψ, ?, ω to give n (= 360/θ), the number of residues per turn; h, the rise per residue; and r, the radius of the phosphate atoms from the helix axis. The single-stranded regions of tRNA that have A-form residues have a notably lower value of n than the double-stranded regions. The G-U “wobble” base pair is shown to be an energetically strained left-handed form. The A-form dinucleoside monophosphates also have a low value of n. A model of UpAl polymerized as a fourfold left-handed helix with the bases on the outside and phosphates on the inside is investigated for its sharp 90° turn angle characteristics. UpA2 cannot be polymerized due to a low values of h (1.31 Å) and r (2.72 Å), which cause steric hindering. An eightfold model of poly(rA) is discussed as are the nonhelical residues of tRNA. Finally, the effects of small changes in dihedral angles and bond lengths and angles on the helical parameters are investigated and discussed by way of explaining this behavior.     

Back to Units of Protein Folding

Abstract

In response to the criticism by A. Finkelstein (J. Biomol. Struct. Dyn. 20, 311–314, 2002) of our Communication (J. Biomol. Struct. Dyn. 20, 5–6, 2002) several issues are dealt with. Importance of the notion of elementary folding unit, its size and structure, and the necessity of further characterization of the units for the elucidation of the protein folding in vivo are discussed. The criticism (J. Biomol. Struct. Dyn. 20, 311–314, 2002) on the hierarchical protein folding is also briefly addressed.     

An Automatic Search for Similar Spatial Arrangements of α-Helices and β-Strands in Globular Proteins

Abstract

A fast search algorithm to reveal similar polypeptide backbone structural motifs in proteins is proposed. It is based on the vector representation of a polypeptide chain fold in which the elements of regular secondary structures are approximated by linear segments (Abagyan and Maiorov, J. Biomol. Struct. Dyn. 5, 1267–1279 (1988)). The algorithm permits insertions and deletions in the polypeptide chain fragments to be compared. The fast search algorithm implemented in FASEAR program is used for collecting βαβ supersecondary structure units in a number of α/β proteins of Brookhaven Data Bank. Variation of geometrical parameters specifying backbone chain fold is estimated. It appears that the conformation of the majority of the fragments, although almost all of them are right-handed, is quite different from that of standard βαβ units. Apart from searching for specific type of secondary structure motif, the algorithm allows automatically to identify new recurrent folding patterns in proteins. It may be of particular interest for the development of tertiary template approach for prediction of protein three-dimensional structure as well for constructing artificial polypeptides with goal-oriented conformation.     

Experimental support for a right-handed vertical double helix

The experimentally observed geometry of a miniature double helix of high anti (χCN fixed ? 120°) nucleic acid structure indicates substantial interstrand base stacking with little intrastrand base stacking. This geometry is consistent with the right-handed vertical double helix in which the base planes are parallel to helical axis proposed by Olson [(1977) Proc. Natl. Acad. Sci. USA 74 , 1775] from theoretical calculations. The experimental data do not agree with the left-handed model in which the base planes are perpendicular to helical axis that has been proposed by Yathindra and Sundaralingam [(1976) Nucl. Acids Res. 3 , 729]. The theoretically computed chemical shift changes for the various double-helical configurations show that for the system examined, only the vertical model can explain the experimentally observed shift trends. The melting curve for the helix–coil transition for high anti dinucleoside monophosphate has been observed for the first time. Even though the experimental data support vertical double helices when χCN is fixed at 120°, data on naturally occurring nucleic acid structures indicate that they have no proclivity to enter into vertically stabilized double-helical arrays.     

Optimization of interatomic angle geometry in torsion angle space

Summary Interatomic angle constraints are usually specified as distance constraints in torsion angle distance geometry. Such an approach is inaccurate and often inadequate. We provide a direct definition of the inter-atomic angle constraint term, which can be incorporated in the target function. The first derivative of this term with respect to the torsion angle has been described for all possible cases. This feature has been implemented in the nucleic acid distance geometry program TANDY [Ajay Kumar et al. (1991)J. Biomol NMR,1, 363–378], and has been tested on base pairing in the DNA fragment, d(AT)₂. The results clearly indicate the need and adequacy of such angle constraints. Other applications that would also benefit from this technique have been identified.