The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (alpha and epsilon) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A-like base-pair stacking, adopt backbone conformations distinct from those in PNA.DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.     

Conformational analysis of protein and nucleic acid fragments with the new grid search algorithm FOUND

The new computer algorithm FOUND, which is implemented as an integrated module of the DYANA structure calculation program, is capable of performing systematic local conformation analyses by exhaustive grid searches for arbitrary contiguous fragments of proteins and nucleic acids. It uses torsion angles as the only degrees of freedom to identify all conformations that fulfill the steric and NMR-derived conformational restraints within a contiguous molecular fragment, as defined either by limits on the maximal restraint violations or by the fragment-based DYANA target function value. Sets of mutually dependent torsion angles, for example in ribose rings, are treated as a single degree of freedom. The results of the local conformation analysis include allowed torsion angle ranges and stereospecific assignments for diastereotopic substituents, which are then included in the input of a subsequent structure calculation. FOUND can be used for grid searches comprising up to 13 torsion angles, such as the backbone of a complete -helical turn or dinucleotide fragments in nucleic acids, and yields a significantly higher number of stereospecific assignments than the precursor grid search algorithm HABAS.     

Structure of a triple helical DNA with a triplex-duplex junction

Extended purine sequences on a DNA strand can lead to the formation of triplex DNA in which the third strand runs parallel to the purine strand. Triplex DNA structures have been proposed to play a role in gene expression and recombination and also have potential application as antisense inhibitors of gene expression. Triplex structures have been studied in solution by NMR, but have hitherto resisted attempts at crystallization. Here, we report a novel design of DNA sequences, which allows the first crystallographic study of DNA segment containing triplexes and its junction with a duplex. In the 1.8 A resolution structure, the sugar-phosphate backbone of the third strand is parallel to the purine-rich strand. The bases of the third strand associate with the Watson and Crick duplex via Hoogsteen-type interactions, resulting in three consecutive C(+).GC, BU.ABU (BU = 5-bromouracil), and C(+).GC triplets. The overall conformation of the DNA triplex has some similarity to the B-form, but is distinct from both A- and B-forms. There are large changes in the phosphate backbone torsion angles (particularly gamma) of the purine strand, probably due to the electrostatic interactions between the phosphate groups and the protonated cytosine. These changes narrow the minor groove width of the purine-Hoogsteen strands and may represent sequence-specific structural variations of the DNA triplex.     

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I' beta-turn forming tetrapeptide and an alpha-helical tetradecapeptide.

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I'beta-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt alpha-helical conformations stabilized by 11 successive 5 --> 1 hydrogen bonds. In addition, a single 4 --> 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (phi, psi) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle N-C(alpha)-C' (tau) and the observed backbone phi,psi values. For tau > 106 degrees, helices are observed, while fully extended structures are characterized by tau < 106 degrees. The mean tau values for extended and folded conformations for the Dpg residue are 103.6 degrees +/- 1.7 degrees and 109.9 degrees +/- 2.6 degrees, respectively.     

The crystal structure of the RNA/DNA hybrid r(GAAGAGAAGC). d(GCTTCTCTTC) shows significant differences to that found in solution.

The crystal structure of the RNA/DNA hybrid r(GAAGAGAAGC). d(GCTTCTCTTC) has been solved and refined at 2.5 A resolution. The refinement procedure converged at R = 0.181 for all reflections in the range 20.0-2.5 A. In the crystal, the RNA/DNA hybrid duplex has an A' conformation with all but one of the nucleotide sugar moieties adopting a C3'- endo (N) conformation. Both strands in the double helix adopt a global conformation close to the A-form and the width of the minor groove is typical of that found in the crystal structures of other A-form duplexes. However, differences are observed between the RNA and DNA strands that make up the hybrid at the local level. In the central portion of the duplex, the RNA strand has backbone alpha, beta and gamma torsion angles that alternate between the normal gauche -/ trans / gauche + conformation and an unusual trans / trans / trans conformation. Coupled with this so-called 'alpha/gamma flipping' of the backbone torsion angles, the distance between adjacent phosphorous atoms on the RNA strand systematically varies. Neither of these phenomena are observed on the DNA strand. The structure of the RNA/DNA hybrid presented here differs significantly from that found in solution for this and other sequences. Possible reasons for these differences and their implications for the current model of RNase H activity are discussed.     

Stress and strain in staphylococcal nuclease.

Protein molecules generally adopt a tertiary structure in which all backbone and side chain conformations are arranged in local energy minima; however, in several well-refined protein structures examples of locally strained geometries, such as cis peptide bonds, have been observed. Staphylococcal nuclease A contains a single cis peptide bond between residues Lys 116 and Pro 117 within a type VIa beta-turn. Alternative native folded forms of nuclease A have been detected by NMR spectroscopy and attributed to a mixture of cis and trans isomers at the Lys 116-Pro 117 peptide bond. Analyses of nuclease variants K116G and K116A by NMR spectroscopy and X-ray crystallography are reported herein. The structure of K116A is indistinguishable from that of nuclease A, including a cis 116-117 peptide bond (92% populated in solution). The overall fold of K116G is also indistinguishable from nuclease A except in the region of the substitution (residues 112-117), which contains a predominantly trans Gly 116-Pro 117 peptide bond (80% populated in solution). Both Lys and Ala would be prohibited from adopting the backbone conformation of Gly 116 due to steric clashes between the beta-carbon and the surrounding residues. One explanation for these results is that the position of the ends of the residue 112-117 loop only allow trans conformations where the local backbone interactions associated with the phi and psi torsion angles are strained. When the 116-117 peptide bond is cis, less strained backbone conformations are available. Thus the relaxation of the backbone strain intrinsic to the trans conformation compensates for the energetically unfavorable cis X-Pro peptide bond. With the removal of the side chain from residue 116 (K116G), the backbone strain of the trans conformation is reduced to the point that the conformation associated with the cis peptide bond is no longer favorable.     

Protein backbone angle prediction with machine learning approaches

MOTIVATION: Protein backbone torsion angle prediction provides useful local structural information that goes beyond conventional three-state (alpha, beta and coil) secondary structure predictions. Accurate prediction of protein backbone torsion angles will substantially improve modeling procedures for local structures of protein sequence segments, especially in modeling loop conformations that do not form regular structures as in alpha-helices or beta-strands. RESULTS: We have devised two novel automated methods in protein backbone conformational state prediction: one method is based on support vector machines (SVMs); the other method combines a standard feed-forward back-propagation artificial neural network (NN) with a local structure-based sequence profile database (LSBSP1). Extensive benchmark experiments demonstrate that both methods have improved the prediction accuracy rate over the previously published methods for conformation state prediction when using an alphabet of three or four states. AVAILABILITY: LSBSP1 and the NN algorithm have been implemented in PrISM.1, which is available from www.columbia.edu/~ay1/. SUPPLEMENTARY INFORMATION: Supplementary data for the SVM method can be downloaded from the Website www.cs.columbia.edu/compbio/backbone.     

Lipid conformation in crystalline bilayers and in crystals of transmembrane proteins

Dihedral torsion angles evaluated for the phospholipid molecules resolved in the X-ray structures of transmembrane proteins in crystals are compared with those of phospholipids in bilayer crystals, and with the phospholipid conformations in fluid membranes. Conformations of the lipid glycerol backbone in protein crystals are not restricted to the gauche C1-C2 rotamers found invariably in phospholipid bilayer crystals. Lipid headgroup conformations in protein crystals also do not conform solely to the bent-down conformation, with gauche-gauche configuration of the phospho-diester, that is characteristic of phospholipid bilayer membranes. This suggests that the lipids that are resolved in crystals of membrane proteins are not representative of the entire lipid-protein interface. Much of the chain configurational disorder of the membrane-bound lipids in crystals arises from energetically disallowed skew conformations. This indicates a configurational heterogeneity in the lipids at a single binding site: eclipsed conformations occur also in some glycerol backbone torsion angles and C-C torsion angles in the lipid headgroups. Stereochemical violations in the protein-bound lipids are evidenced by one-third of the ester carboxyl groups in non-planar configurations, and certain of the carboxyls in the cis configuration. Some of the lipid structures in protein crystals have the incorrect enantiomeric configuration of the glycerol backbone, and many of the branched methyl groups in structures of the phytanyl chains associated with bacteriorhodopsin crystals are in the incorrect S-configuration.     

Stacking-unstacking of the dinucleoside monophosphate guanylyl-3'',5''-uridine studied with molecular dynamics.

Molecular dynamics simulations were carried out on two conformations of the dinucleoside monophosphate guanylyl-3',5'-uridine (GpU) in aqueous solution with one sodium counterion. One stacked conformation and one with the C3'-O3'-P-O5' backbone torsion angle twisted 180 degrees to create an unstacked conformation. We observed a relatively stable behavior of the stacked conformation, which remained stacked throughout the simulation, whereas the unstacked conformation showed major changes in the backbone torsion and glycosidic angles. During the simulation the unstacked conformation transformed into a more stacked form and then back again to an unstacked one. The calculated correlation times for rotational diffusion from the molecular dynamics simulations are in agreement with fluorescence anisotropy and nuclear magnetic resonance data. As expected, the correlation times for rotational diffusion of the unstacked conformation were observed to be longer than for the stacked conformation. The 2'OH group may contribute in stabilizing the stacked conformation, where the O2'-H...O4' hydrogen bond occurred in 82.7% of the simulation.     

Crystal and solution structures of the B-DNA dodecamer d(CGCAAATTTGCG) probed by Raman spectroscopy: heterogeneity in the crystal structure does not persist in the solution structure

The self-complementary dodecamer d(CGCAAATTTGCG) crystallizes as a double helix of the B form and manifests a Raman spectrum with features not observed in Raman spectra of either DNA solutions or wet DNA fibers. A number of Raman bands are assigned to specific nucleoside sugar and phosphodiester conformations associated with this model B-DNA crystal structure. The Raman bands proposed as markers of the crystalline B-DNA structure are compared and contrasted with previously proposed markers of Z-DNA and A-DNA crystals. The results indicate that the three canonical forms of DNA can be readily distinguished by Raman spectroscopy. However, unlike Z-DNA and A-DNA, which retain their characteristic Raman fingerprints in aqueous solution, the B-DNA Raman spectrum is not completely conserved between crystal and solution states. The Raman spectra reveal greater heterogeneity of nucleoside conformations (sugar puckers) in the DNA molecules of the crystal structure than in those of the solution structure. The results are consistent with conversion of one-third of the dG residues from the C2'-endo/anti conformation in the solution structure to another conformation, deduced to be C1'-exo/anti, in the crystal. The dodecamer crystal also exhibits unusually broad Raman bands at 790 and 820 cm-1, associated with the geometry of the phosphodiester backbone and indicating a wider range of (alpha, zeta) backbone torsion angles in the crystal than in the solution structure. The results suggest that backbone torsion angles in the CGC and GCG sequences, which flank the central AAATTT sequence, are significantly different for crystal and solution structures, the former containing the greater diversity.(ABSTRACT TRUNCATED AT 250 WORDS)     

Conformational analysis of a dinucleotide photodimer with the aid of the genetic algorithm.

The solution structure of the photodimer cis,syn-dUp[]dT is derived with the aid of the genetic algorithm. The conformational space available for the molecule is sampled efficiently using the computer program DENISE and tested against a set of constraints available from nmr experiments. The dominant conformation in solution found with this approach can be described by the following combinations of sugar-phosphate backbone torsion angles: epsilon(t), zeta(t), alpha(+), beta(-ac), and gamma(t). The conformation of the sugars and glycosidic torsion angles are S type and syn, respectively. The cyclobutane ring and pyrimidines are puckered. In addition, other conformations that exist in equilibrium with the first are found. It is concluded that the cyclobutane-pyrimidine system is rigid, whereas the sugar-phosphate backbone is flexible. The solution structures are compared with the crystal structure of the strongly related cyano-ethyl ester of cis,syn-dTp[]dT.     

Conformational analysis of a segment in bacterioopsin by two-dimensional (1)H-NMR spectroscopy]

The spatial structure of a synthetic peptide, an analogue of the membrane spanning segment B (residues 34-65) of bacterioopsin from Halobacterium halobium, has been refined. Backbone torsion angles were derived from intensities of short-range interproton NOEs. These, together with a complete set of the NOEs integral intensities formed the basis for the three-dimensional structure refinement by the energy minimization with consideration of NOE penalty functions. Analysis indicates the right-handed alpha-helical conformation of segment B extending from Asp-38 to Tyr-64 with a kink of the helical axis (27 degrees) at Pro-50. The most stable region with an average root-mean-square deviation of 0.43 A between the backbone atoms includes residues 42-60 in six energy refined structures. The N-terminal part of segment B (residues 34-37) has no ordered conformation. The inferred structure is in close agreement with the electron cryomicroscopy structure of bacteriorhodopsin, differing from it in conformations of most of the side chains.     

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.     

Site-specific identification of non-beta-strand conformations in Alzheimer's beta-amyloid fibrils by solid-state NMR

The most well-established structural feature of amyloid fibrils is the cross-beta motif, an extended beta-sheet structure formed by beta-strands oriented perpendicular to the long fibril axis. Direct experimental identification of non-beta-strand conformations in amyloid fibrils has not been reported previously. Here we report the results of solid-state NMR measurements on amyloid fibrils formed by the 40-residue beta-amyloid peptide associated with Alzheimer's disease (Abeta(1-40)), prepared synthetically with pairs of (13)C labels at consecutive backbone carbonyl sites. The measurements probe the peptide backbone conformation in residues 24-30, a segment where a non-beta-strand conformation has been suggested by earlier sequence analysis, cross-linking experiments, and molecular modeling. Data obtained with the fpRFDR-CT, DQCSA, and 2D MAS exchange solid-state NMR techniques, which provide independent constraints on the phi and psi backbone torsion angles between the labeled carbonyl sites, indicate non-beta-strand conformations at G25, S26, and G29. These results represent the first site-specific identification and characterization of non-beta-strand peptide conformations in an amyloid fibril.     

Determination of nucleic acid backbone conformation by 1H NMR.

In DNA or RNA duplexes, the six-bond C3'-O3'-P-O5'-C5'-C4'-C3' backbone linkage connecting adjacent residues contains six torsion angles (epsilon, zeta, alpha, beta, gamma, delta) but only four protons. This seriously limits the ability to define the backbone conformation by NMR using purely 1H-1H distance geometry (DG) methods. The problem is further compounded by the inability to assign two of the four backbone protons, namely the poorly resolved H5' and H5' protons, and invariably leads to DG structures with poorly defined backbone conformations. We have developed and tested a reliable method to constrain the beta, gamma, and epsilon (and indirectly alpha and zeta) backbone torsion angles by lower-bound NOE distances to unassigned H5'/H5' resonances combined with either 1H line widths or the conservative use of sigma J measurements; the method relies only on 1H 2-D NMR data, does not involve any structural assumptions, and leads to much improved backbone convergence among DG structures. The C4'-C5' torsion angle gamma is constrained by lower-bound NOE distances from H2' and from H6/H8 to any H5'/H5', as well as by sigma JH4, coupling measurements in the 3.9-4.4 ppm region; delta is constrained by H1'-H4' NOE distances and by H3'-H4' and H3'-H2' J couplings in COSY data; epsilon is partially constrained by H3' line width and/or further constrained by subtracting the minimum possible sigma JH3'-H from the observed sigma JH3' (COSY) to arrive at the maximum possible JH3'-P, which is then converted to H3'-P distance bounds. The angle beta is partially constrained via H5'-P and H5'-P distance bounds consistent with the maximum H5'-P and H5'-P J couplings derived from the observed H5' and H5' line widths, while alpha and zeta are indirectly constrained by lower distance bounds on the observed (n)H1' to (n + 1)H5'/H5' NOEs combined with the prior partial constraints on beta, gamma, delta, and epsilon. The combined effects of these additional constraints in determining distance geometry structures have been demonstrated using a 12-base duplex, [d(GCCGTTAACGGC)]2. Coordinate RMSDs per atom between structures refined with these constraints from random-embedded DG structures, from ideal A-DNA, and from B-DNA starting structures were less than 0.4 A for the central 8 base pairs indicating good convergence. All backbone angles for the central 8 base pairs are very well constrained with less than 10 degrees variation in any of the 48 torsion angles.     

A stereochemical model for phospholipids. II: Conformational analysis and molecular packing of phosphatidylethanolamine (PE)

A partition energy method procedure was applied to select the energetically favoured conformations of phosphatidylethanolamine (PE) as polar constituents of phospholipid molecules. The result indicated a large degree of freedom for the two torsion angles of the ester bond of the phosphate and a gauche, gauche star conformation for the ethane bond.A packing process of the molecule was carried out through a potential energy calculation by considering the conformers selected above, using previously published procedure and conventions. All the arrangements which possess the best packing energy values were characterised by an orientation of the PN dipolar segment parallel to the lattice plain. Rotation of the internal torsion angles and rotation in the eulerian space of the molecule produced differences in the charged groups that interact. An additional minimum was present in the energy packing process of those conformers which have the first torsion angle of the phosphate in a trans conformation. This minimum, which corresponds to an orientation of the molecule orthogonal to the lattice plane, requires a complete neutralisation of the point charges on the system.The results of the calculation underline the importance of changes in the behaviour of the polar group of the phospholipids in the packing process.     

TANGLE: two-level support vector regression approach for protein backbone torsion angle prediction from primary sequences

Protein backbone torsion angles (Phi) and (Psi) involve two rotation angles rotating around the C(α)-N bond (Phi) and the C(α)-C bond (Psi). Due to the planarity of the linked rigid peptide bonds, these two angles can essentially determine the backbone geometry of proteins. Accordingly, the accurate prediction of protein backbone torsion angle from sequence information can assist the prediction of protein structures. In this study, we develop a new approach called TANGLE (Torsion ANGLE predictor) to predict the protein backbone torsion angles from amino acid sequences. TANGLE uses a two-level support vector regression approach to perform real-value torsion angle prediction using a variety of features derived from amino acid sequences, including the evolutionary profiles in the form of position-specific scoring matrices, predicted secondary structure, solvent accessibility and natively disordered region as well as other global sequence features. When evaluated based on a large benchmark dataset of 1,526 non-homologous proteins, the mean absolute errors (MAEs) of the Phi and Psi angle prediction are 27.8° and 44.6°, respectively, which are 1% and 3% respectively lower than that using one of the state-of-the-art prediction tools ANGLOR. Moreover, the prediction of TANGLE is significantly better than a random predictor that was built on the amino acid-specific basis, with the p-value<1.46e-147 and 7.97e-150, respectively by the Wilcoxon signed rank test. As a complementary approach to the current torsion angle prediction algorithms, TANGLE should prove useful in predicting protein structural properties and assisting protein fold recognition by applying the predicted torsion angles as useful restraints. TANGLE is freely accessible at http://sunflower.kuicr.kyoto-u.ac.jp/~sjn/TANGLE/.     

Theoretical conformational analysis of methylamide of N-acetyl-L-lysine]

The spatial structure of the methylamide of N-acetyl-L-lysine has been analysed taking into account non-bonded and electrostatic interactions, torsional energy, bond angles distortion and hydrogen bonding. Conformational capacities of the backbone and mutual dependence of spatial structures of the backbone and the side chain was described by conformational maps obtained by energy minimisation, the dihedral angles and the bond angles of the side chain being varied for every phi, psi point. Every possible combination for phi, psi, x1-x5-angles was used corresponding to the stable form of the backbone and to torsion potential minima of the initial approximations in the calculation of preferred conformations of the molecule. Comparisons are made between stable forms of the methylamide of N-acetyl-L-lysine and Lys residues in proteins with known structure.