Conformational studies on pyrimidine 5'-monophosphates and 3',5'-diphosphates. Effect of the phosphate groups on the backbone conformation of polynucleotides

In continuation of our studies on the effect of the base and the phosphate groups on the glycosyl and the sugar-phosphate backbone conformation, we have carried out semi-empirical potential energy calculations on the common 5′- and 3′5′-ribopyrimidine mono- and diphosphates by considering simultaneous rotations about the glycosyl (χ) and the C(4′)–C(5′) (ψ) bonds. This calculation provides an assessment of the nature and orientation of the base on the sugar–phosphate backbone conformation of nucleotides and polynucleotides. It is found that the attractive inetractions between the 5′-phosphate group and the base mutually stabilize the antiand the gauche-gauche (gg) conformations about χ and ψ, respectively, in 5′-ribopyrimidine nucleotides. The introduction of the 3′-phosphate group as in 3′,5′-ribopyrimidine diphosphates, still leaves the anti-gg as the most favored conformation with the important difference that the probability of occurrence of the anti, gauche-trans (gt) is how substantially increased. This is dependent to a large extent on the sugar conformation and to a lesser extent on the base. Uracil and thymine show a greater probability for the anti-gt than cytosine. The syn conformation is considerably less likely and its occurrence is also dependent on the base type, cytosine showing a lesser tendency than uracil and thymine. For the syn base, the most favourec conformation for ψ is gt, since gg is sterically disallowed and tg is destabilized by electrostatic repulsive interactions between the 3′ and 5′-phosphate groups. Thus, there is a striking correlation between the glycoysl and the backbone C(4′)–C(5′) bond conformations. The rest of the bonds of the backbone are considerable less dependent on the glycosyl conformation. These studies reveal that in poly-ribopyrimidine nucletides the majority of the nucleotide residues are expected to occur in the anti-gg conformation.     

Configuration-dependent properties of randomly coiling polynucleotide chains. I. A comparison of theoretical energy estimates

Various theoretical estimates of the conformational energy associated with polynucleotides in solution have been compared with each other and also with the experimentally observed conformations found in X-ray crystallographic investigations of low-molecular-weight nucleic acid analogs. In view of the disparities between these data, certain configuration-dependent properties (i.e., the mean-square unperturbed end-to-end distance 〈r²〉₀ and the average vicinal nmr coupling constant 〈J〉) appropriate to randomly coiling polynucleotides described by either the energy estimates or by the crystallographically preferred conformations have also been calculated and compared with the known solution behavior of polynucleotide chains. Both the theoretical energy surfaces and the X-ray data show good agreement with the nmr coupling constant indications of the preferred rotations about the O-C and C-C bonds of the chain backbone. The principal discrepancies between the theoretical methods and X-ray data arise in their ability to predict successfully the preferred rotations about the two phosphodiester bonds of the chain backbone and the unperturbed dimensions of randomly coiling polynucleotide chains.     

Conformational studies on polynucleotide chains. III. Intramolecular energy maps and comparison with experiments

Conformational energy maps for the four combinations of two consecutive torsional angles of the backbone structure of polydeoxyribonucleotides are presented. Both the C(2′)-endo and the C(3′)-endo conformation of sugar rings were considered. The energies were evaluated with an analytical expression representing the best fit to ab initio energies computed in the Hartree-Fock approximation, and consisting of a contribution from nonbonded interactions of the Lennard-Jones 6-12 type and an intrinsic torsional potential. It is shown that the minima of these maps are in excellent agreement with the most stable conformations as obtained from x-ray crystallographic analysis of nucleic acids and polynucleotides.     

Configuration-dependent properties of randomly coiling polynucleotide chains. II. The role of the phosphodiester linkage

The dependence of the unperturbed dimensions of randomly coiling polynucleotides on the rotations about the phosphodiester linkages of the chain has been examined in order to understand the conformational discrepancies, set forth in paper I, regarding these angles (ω′ and ω). Large values of the characteristic ratio 〈r²〉₀/nl² , which agree with the experimental behavior of the chain, are obtained only if a sizeable proportion of the polymer residues have trans ω′ values. The asymmetric torsional potential that is believed to arise from gauche effects associated with the P-O bonds has been approximated using a hard core model. The calculated characteristic ratio exhibits a strong dependence upon the magnitude of this torsional barrier (separating trans and gauche conformations) and shows agreement with experimental values for polyribonucleotides only if this energy difference is 1 kcal/mol or less.     

Stereochemistry of nucleic acids and their constituents. IV. Allowed and preferred conformations of nucleosides,nucleoside mono-, di-, tri-, tetraphosphates,nucleic acids and polynucleotides

A uniform notation and convention is suggested to describe the torsional angles in nucleic acids and their derivatives. The torsional angle χ, relating the stereochemistry of the base with respect to the sugar, shows more variation for the β-purine glycosides than for the β-pyrimidine glycosides. This variation is attributed to the fact that the β-purine derivatives may form intramolecular O(5′)-H…N(3) hydrogen bonding. The χ values for the α-purine and α-pyrimidine glycosides show preference for the –syn-clinal (or anti) conformation. The mode of puckering of the sugar also influences the χ value. The various possible conformations for the furanose ring are described by the torsional angles τ₀ τ₁, τ₂, τ₃, τ₄, about the five ring bonds. From an analysis of the torsional angles (ω, ?, ψ, ψ′, ?′, ω′) about the sugar phosphate bonds in the x-ray structures of the known nucleosides, nucleotides, phosphodiesters, nucleic acids, and related compounds, and from a consideration of molecular models, it is found that the possible conformations for the backbone of helical nucleic acids is strikingly limited. Most importantly, the preferred conformation of the nucleotide unit in poly nucleotides and nucleic acids turns out to be the same as that found for the nucleotide in the crystal structure. It is observed that base “stacking” is a consequence of the restricted backbone conformation. The torsional angles are illustrated in the form of conformational “wheels”. Interrelation between the torsion angles about successive pairs of sugar-phosphate bonds are presented in the form of conformational maps: ω,?; ?,ψ; ψ.ψ′; ψ′,?′; ?′,ω′; ω′,ω. The ω′,ω map shows the perferred conformations about the inter-nucleotide bonds of right- and left-handed helices and the possible conformations of phosphodiesters. The preferred conformation of the pyrophosphate and triphosphate is that in which the phosphate oxygens display a staggered arrangement when viewed along the P–P axis. A plausible structure and conformation for the ATPM^2? backbound complex is presented. This structure differs from that proposed by SzentGyorgi in that the metal (only transition metals are considered here) is not bound to the NH₂ nitrogen of adenine, but rather is simultaneously bound to N(7) of the ring and three phosphates (α, β, γ), or N(7) of the ring and two phosphates (β, γ). The remaining metal coordination may be satisfied by solvent–metal or enzyme–metal bonds.     

The molecular mechanics program DUPLEX: Computing structures of carcinogen modified DNA by surveying the potential energy surface

The nucleic acids molecular mechanics program DUPLEX has been designed with useful features for surveying the potential energy surface of polynucleotides, especially ones that are modified by polycyclic aromatic carcinogens. The program features helpful strategies for addressing the multiple minimum problem: (1) the reduced variable domain of torsion angle space; (2) search strategies that emphasize large scale searches for smaller subunits, followed by building to larger units by a variety of strategies; (3) the use of penalty functions to aid the minimizer in locating selected structural types in first stage minimizations; penalty functions are released in terminal minimizations to yield final unrestrained minimum energy conformations. Predictive capability is illustrated by DNA modified by activated benzo[a]pyrenes.     

The occurrence of the syn-C3′ endo conformation and the distorted backbone conformations for C4′-C5′ and P-O5′ in oligo and polynucleotides

Abstract

The nucleoside constituents of nucleic acids prefer the anti conformation (1). When the sugar pucker is taken into account the nucleosides prefer the C2′endo-anti conformation. Of the nearly 300 nucleosides known, about 250 are in the anti conformation and 50 are in the syn-conformation, i.e., anti to syn conformation is 5:1. The nucleotide building blocks of nucleic acids show the same trend as nucleosides. Both the deoxy-guanosine and ribo- guanosine residues in nucleosides and nucleotides prefer the syn-C2′endo conformation with an intra-molecular hydrogen bond (for nucleosides) between the O5′- H and the N3 of the base and, a few syn-C3′endo conformations are also observed. Evidence is presented for the occurrence of the C3′endo-syn conformation for guanines in mis-paired double helical right-handed structures with the distorted sugar phosphate C4′-C5′ and P-O5′ bonds respectively, from g+ (gg) and g- to trans. Evidence is also provided for guanosine nucleotides in left-handed double-helical (Z-DNA) oligo and polynucleotides which has the same syn-C3′endo conformation and the distorted backbone sugar-phosphate bonds (C4′-C5′ and P- O5′) as in the earlier right-handed case.     

Conformations of cyclic 2,3-nucleotides and 2,3-Cyclic-5-diphosphates. Interrelation between the phase angle of pseudorotation of the sugar and the torsions about the glycosyl C(1)-N and the backbone C(4)-C(5) bonds

Semiempirical potential energy calculations have been carried out for cyclic 2′,3′-nucleotides and their 5′-phosphorylated derivatives, which are the intermediates in the hydrolysis of RNA. Calculations have been performed for both purine and pyrimidine bases for the observed O(1′)-endo, O(1′)-exo and the unpuckered planar sugar ring conformations. It is found that the mode of sugar pucker largely determines the preferred conformations of these molecules. For cyclic 2′,3′-nucleotides themselves, the O(1′)-endo sugars show a preference for the syn glycosyl conformation while the O(1′)-exo sugars exclusively favor the anti conformation regardless of whether the base is a purine or pyrimidine. For the unpuckered planar sugar, the syn conformation is favored for purines and anti for pyrimidines. Both the gauche (+) (60°) and trans (180°) conformations about the C(4′)? C(5′) bond are favored for O(1′)-endo sugars, while the gauche (?) (300°) and trans (180°) are favored for O(1′)-exo sugars. On the contrary, the 5′-phosphorylated cyclic 2′,3′-nucleotides of both purines and pyrimidines show a preference for the anti-gauche (+) conformational combination about the glycosyl and C(4′)? C(5′) bonds for the O(1′)-endo sugars and the anti-trans combination for the O(1′)-exo sugars. The correlation between the phase angle of the sugar ring and the favored torsions about the glycosyl and the backbone C(4′)? C(5′) bonds as one traverses along the pseudorotational pathway of the sugar ring is examined.     

The occurence of the syn-C3' endo conformation and the distorted backbone conformations for C4'-C5' and P-O5' in oligo and polynucleotides.

The nucleoside constituents of nucleic acids prefer the anti conformation (1). When the sugar pucker is taken into account the nucleosides prefer the C2'endo-anti conformation. Of the nearly 300 nucleosides known, about 250 are in the anti conformation and 50 are in the syn-conformation, i.e., anti to syn conformation is 5:1. The nucleotide building blocks of nucleic acids show the same trend as nucleosides. Both the deoxy-guanosine and riboguanosine residues in nucleosides and nucleotides prefer the syn-C2'endo conformation with an intra-molecular hydrogen bond (for nucleosides) between the O5'-H and the N3 of the base and, a few syn-C3'endo conformations are also observed. Evidence is presented for the occurrence of the C3'endo-syn conformation for guanines in mis-paired double helical right-handed structures with the distorted sugar phosphate C4'-C5' and P-O5' bonds respectively, from g+ (gg) and g- to trans. Evidence is also provided for guanosine nucleotides in left-handed double-helical (Z-DNA) oligo and polynucleotides which has the same syn-C3'endo conformation and the distorted backbone sugar-phosphate bonds (C4'-C5' and P-O5') as in the earlier right-handed case.     

Flexibility of the furanose ring in nucleic acids: a compilation of crystal data

A compilation of crystal structure data on deoxyribo- and ribonucleosides and their higher derivatives is presented. The aim of this paper is to highlight the flexibility of deoxyribose and ribose rings. So far, the conformational parameters of nucleic acids constituents of ribose and deoxyribose have not been analysed separately. This paper aims to correlate the conformational parameters with the nature and puckering of the sugar. Deoxyribose puckering occurs in the C2′ endo region while ribose puckering is observed both in the C3′ endo and C2′ endo regions. A few endocyclic and exocyclic bond angles depend on the puckering and the nature of the sugar. The majority of structures have an anti conformation about the glycosyl bond. There appears to be a puckering dependence on the torsion angle about the C4′C5′ bonds. Such stereochemical information is useful in model building studies of polynucleotides and nucleic acids.     

Prediction of the native conformation of a polypeptide by a statistical-mechanical procedure. II. Average backbone structure of enkephalin

The average conformation of Met-enkephalin was determined by using an adaptive, importance-sampling Monte Carlo algorithm (SMAPPS—Statistical Mechanical Algorithm for Predicting Protein Structure). In the calculation, only the backbone dihedral angles (? and ψ) were allowed to vary; i.e., all side-chain (χ) and peptide-bond (ω) dihedral angles were kept fixed at the values corresponding to a low-energy structure of the pentapeptide. The total conformational energy for each randomly generated structure of the polypeptide was obtained by summing over the interaction energies of all pairs of nonbonded atoms of the whole molecule. The interaction energies were computed by the program ECEPP/2 (Empirical Conformational Energy Program for Peptides). Solvent effects were not included in the computation. The calculation was repeated until a total of 10 independent average conformations were established. The regions of conformational space occupied by the average structures were compared with the regions of low conditional free energy obtained by SMAPPS in the first paper of this series. Such a comparison provides an analysis of the capacity of SMAPPS to adjust the Monte Carlo search to regions of highest probability. The results demonstrate that the ability of SMAPPS to focus the Monte Carlo search is excellent. Finally, the 10 independent average conformations and the mean of the 10 average structures were utilized as the initial conformations for a direct energy minimization of the pentapeptide. Of the 11 final energy-minimized structures, three of the conformations were found to be equivalent to the conformation of lowest energy determined previously. In addition, all but two of the remaining energy-minimized structures were found to correspond to one of the two other conformations of high probability obtained in the first paper of this series. These results indicate that a set of independent average conformations can provide a rational, unbiased choice for the initial conformation, to be used in a direct energy minimization of a polypeptide. The final energy-minimized structures consequently constitute a set of low-energy conformations, which include the global energy minimum.     

Backbone conformations in deoxyribonucleic acids. Conformational role of the 2′-hydroxyl group on the internucleotide phosphodiester conformations

The conformations accessible to the internucleotide phosphodiester group in deoxydinucleoside monophosphates, deoxydinucleoside triphosphates, and deoxypolynucleotides have been explored in detail by potential energy calculations. The two most predominant conformations for the nucleotide moiety (³E and ²E) and their possible combinations (³E^?3E, ³E^?2E, ²E^?2E, ²E^?3E) have been employed, similar to our earlier studies on polyribonucleotides. The internucleotide P-O bond torsions are very sensitive to the sugar pucker (³E and ²E) and sugar type (ribose and 2′-deoxyribose) on the 3′-residue of dinucleoside phosphates. The preferred phosphodiester conformations found for the deoxydinucleoside monophosphates and triphosphates, in general, follow the same pattern as those obtained for ribose sugars when the sugar on the 3′-side of the molecule has the ³E sugar-ring conformation. The internucleotide P-O bonds show a greater degree of conformational freedom when the 3′-sugar has the ²E pucker. The double gauche g^?g^? conformation for the phosphodiester, which leads to the overlap of the adjacent bases, is shown to be one of the energetically most favored conformations for all the sequence of sugar puckers. It is found that the ²E^?2E sequence of sugar puckers shows a greater energetic preference for the stacked helical conformation (g^?g^?) than the (³E^?3E) and the mixed sugar-pucker combinations. This effect becomes more pronounced in going from a dinucleoside monophosphate to a dinucleoside triphosphate suggesting that the 2′-deoxy sugars favor the ²E sugar pucker in di-, oligo-, and polydeoxyribonucleotide structures. In addition to g^?g^?, the conformations g⁺g^?, tg^?, g^?t, tg⁺, and g⁺t are also found to be possible for the phosphodiester in a polydeoxyribonucleotide and their populations depend to some extent on the sugar-pucker sequence. It is shown that the short-range intramolecular interactions involving the sugar and the phosphate groups dictate to a large extent the backbone conformations of nucleic acids and polynucleotides.     

A protein folding potential that places the native states of a large number of proteins near a local minimum

Background

We present a simple method to train a potential function for the protein folding problem which, even though trained using a small number of proteins, is able to place a significantly large number of native conformations near a local minimum. The training relies on generating decoys by energy minimization of the native conformations using the current potential and using a physically meaningful objective function (derivative of energy with respect to torsion angles at the native conformation) during the quadratic programming to place the native conformation near a local minimum.

Results

We also compare the performance of three different types of energy functions and find that while the pairwise energy function is trainable, a solvation energy function by itself is untrainable if decoys are generated by minimizing the current potential starting at the native conformation. The best results are obtained when a pairwise interaction energy function is used with solvation energy function.

Conclusions

We are able to train a potential function using six proteins which places a total of 42 native conformations within ~4 Å rmsd and 71 native conformations within ~6 Å rmsd of a local minimum out of a total of 91 proteins. Furthermore, the threading test using the same 91 proteins ranks 89 native conformations to be first and the other two as second.     

Biophysical studies on molecular mechanism of abortificient action of prostaglandins—IV. Conformation energy calculations on PGA1, PGB1 and PGE1

This paper reports the conformation energy (CE) calculations on three forms of prostaglandins (PGs) PGA₁, PGB₁ and PGE₁ on the basis of the empirical potential energy functions, for the simultaneous rotations around C₇–C₈ (θ), C₁₂–C₁₃ (β) and C₁₄–C₁₅ (β) bonds [Fig. 1(a)]. The isoenergy contours plotted for θβ rotations for the different β values show the existence of two low energy regions for thg equal to about 90° and 240° in all the three cases. The absolute minimum was obtained for thg = 240° and almost coincided with the crystallographic conformation for PGE₁ and PGA₁. In the case of PGB₁ series of low energy conformations were obtained with the thg values equal to about 90° and 270°, but none of them coincided with the observed crystallographic conformation. The paper discusses the comparison of the different low energy conformations in these three molecules, their biological relevance and the cause of disagreement in the case of PGB₁ with the crystallographic data.     

Conformational analysis of chitobiose and chitosan

The relaxed potential energy surfaces of chitobiose were calculated based on the MM3-force field by optimizing dimer structures on a 10° grid spacing of the torsional angles about the glycosidic bonds (Φ,Ψ). The 36 conformations; the four combinations of the hydroxymethyl group orientations coupled with the nine of the secondary group ones— were assumed for each Φ,Ψ conformation. The four conformations, each differing in the hydroxymethyl group orientations, were considered for the whole Φ,Ψ space, and all the 36 conformations, for the restricted space of low energy. While the resulting energy map and the structures of the energy minima were similar to those proposed for cellobiose in many respects, more restricted energy profile was suggested for the relaxed map of chitobiose where differences in the energy level between the global minimum and the local minima were within 5.4 kcal/mol, compared with the equivalent value of 3.6 kcal/mol for cellobiose. Further depression of the global minimum occurred when the acidic residue was used. The Monte Carlo samples of the chitosan chain were generated based on the relaxed map to predict the unperturbed coil dimension in solution. The chitosan chains showed Gaussian behavior at x = 500 (x, degree of polymerization) and gave the characteristic ratio C_x, of about 70, which was much larger than the experimental values observed for the chitosan and cellulosic chains. © 1994 John Wiley & Sons, Inc.     

The conformational analysis of adenosine triphosphate by classical potential energy calculations

The conformational analysis of adenosine triphosphate was conducted by using classical potential energy calculations. All rotatable bonds were examined, i.e., no dihedral angles were fixed at predetermined conformations except for the ribofuranose ring, which was held in the C(3′)-endo conformation—the conformation observed for adenosine in the crystal state. The energy terms included in the total energy expression consist of nonbonded pairwise interaction, electrostatic pairwise interaction, free energy of solvation, and torsional bond potentials. Two separate approaches were used in the conformational analyses. The first consisted of a sequential fragment approach were four bonds were rotated simultaneously at 30° increments. Each fragment overlapped the preceding one by at least one bond. All rotors were then simultaneously examined at their minima and at ±15°. The second approach consisted of a coarse grid search where all rotors were examined simultaneously, but only at staggered positions. The low-energy conformations thus obtained were then used as starting conformations for a minimization routine based on the method of conjugate directions. The first approach required about 40 hr of central processing unit (CPU) computer time, while the coarse grid/minimization approach required about 4 hr of CPU time. Both the sequential fragment approach and the minimization approach yielded lowest-energy conformations which are remarkably similar to the solid-state conformation of C(3′)-endo ATP.     

Conformational stability in dinucleoside phosphate crystals. Semiempirical potential energy calculations for uridylyl-3'-5'-adenosine monophosphate (UpA) and guanylyl-3',5'-cytidine monophosphate (GpC)

Classical potential energy calculations were performed for the dinucleoside phosphates UpA and GpC. Two widely accessible low-energy regions of conformation space were found for the ω′, ω pair. That of lowest energy contains conformations similar to helical RNA, with ω′ and ω in the vicinity of 300° and 280°, respectively. All five experimental observations of crystalline GpC, two of ApU, and the helical fragment of ApApA fall in this range. The second lowest region has ω′ and ω at about 20° and 80°, respectively, which is in the general region of one experimentally observed crystalline conformer of UpA, and the nonhelical region of ApApA. It is concluded that GpC and ApU, which were crystallized as either sodium or calcium salts, are shielded from each other in the crystal by the water of hydration and are therefore free to adopt their predicted in vacuo minimum energy helical conformations. By contrast, crystalline UpA had only 1/2 water per molecule, and was forced into higher energy conformations in order to maximize intermolecular hydrogen bonding.     

Crystallographic studies of drug-nucleic acid interactions: Proflavine intercalation between the non-complementary base-pairs of cytidilyl-3′,5′-adenosine

In order to get insights into the binding of dyes and mutagens with denatured and single-stranded nucleic acids and the possible implications in frameshift mutagenesis, a 1:1 complex between the non-self-complementary dinucleoside monophosphate cytidilyl-3′,5′-adenosine (CpA) and proflavine was crystallized. The crystals belong to the tetragonal space group P4₂2₁2 with cell constants $\text{a = b = 19.38(1)}\text{A}\text{?}\text{and}\text{c = 27.10(1)}\text{A}\text{?}$. The asymmetric unit contains one CpA, one proflavine and nine water molecules by weight. The structure was determined using Patterson and direct methods and refined to an R-value of 11% using 2454 diffractometer intensities.The non-self-complementary dinucleoside monophosphate CpA forms a selfpaired parallel chain dimer with a proflavine molecule intercalated between the protonated cytosine-cytosine (C · C) pair and the neutral adenine-adenine (A · A) pair. The dimer complex exhibits a right-handed helical twist and an irregular girth. The neutral A · A pair is doubly hydrogen-bonded through the N(6) and N(7) sites (C(1′)C(1′) distance: 10.97(2) Å) and the protonated C · C pair is triply hydrogen-bonded with a proton shared between the N(3) sites (C(1′)C(1′) distance: 9.59(2) Å). To accommodate the intercalating dye, the sugars of successive nucleotide residues adopt the two fundamental conformations (5′ end: 3′-endo, 3′ end: 2′-endo), the backbone adopts torsion angle values that fluctuate within their preferred conformational domains: the PO bonds (ω, ω′) adopt the characteristic helical (gauche^?-gauche^?) conformation, the CO bonds (φ, φ′) are both in the trans domain and the C(4′)C(5′) bonds (ψ) are in the gauche⁺ region. The bases of both residues are disposed in the preferred anti domain with the glycosyl torsion angles (χ) correlated to the puckering mode of the sugar so that the cytidine residue is C(3′)-endo, low χ (12 dg), and the adenosine residue is C(2′)-endo, high χ (84 °). The intercalated proflavine stacks more extensively with the C · C pair than the A · A pair. Between 4₂-related CpA proflavine units there is a second proflavine which stacks well with both the A · A and the C · C pairs sandwiching it. Both proflavine molecules are positionally disordered. In each of its two disordered sites, the intercalated proflavine forms hydrogen-bonded interactions with only one sugar-phosphate backbone. A total of 26 water sites has been characterized of which only two are fully occupied. These hydration sites are involved in an intricate network of hydrogen bonds with both the dye and CpA and provide insights on the various modes of interactions between water molecules and between water molecules and nucleic acids.The structure of the proflavine-CpA complex shows that intercalation of planar drugs can occur between non-complementary base-pairs. This result can be relevant for understanding the strong binding of acridine dyes to denatured DNA, single-stranded RNA, and single-stranded polynucleotides. Also, the ability of proflayine to promote self-pairs of adenine and cytosine bases could provide a chemical basis for an alternative mechanism of frameshift mutagenesis.     

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.     

Conformational geometry and vibrational frequencies of nucleic acid chains.

Normal coordinate analysis of diethyl phosphate has been made, which predicts all observed Raman frequencies in the range 170–1300 cm^?1. The force constants from this calculation have been transferred to a vibrational calculation for a simplified model of the backbone of nucleic acids, which also involves the ? O? PO₂^?? O phosphate group and the ? C₅′? C₄′? C₃′? linkage of the ribose. The coordinates of these atoms are those recently given by Arnott and Hukins, which place the ribose ring of B-DNA in a C₃′-exo conformation. This simple polymer model appears to be able to describe adequately the frequency-dependent changes observed in the Raman spectra arising from the backbone vibrations of nucleic acid in going from the B- to A-form. The symmetric ? O? P? O? diester stretch increases in frequency from about 787 cm^?1 in the B-form to 807 cm^?1 in the A-form. The increased frequency characteristic of the A-form is due to the combining of the diester stretch with vibrations involving the C₅′, C₄′, and C₃′ nuclei. The frequency of the symmetric ? O? P? O? diester stretch is shown to be very dependent on the conformation of the ribose ring, indicating that in polynucleotides the ribose ring takes on one of two rigid conformations: C₃′-endo for A-form or C₃′-exo for B-form and “disordered” polynucleotides. The calculation lends confirmation to the atomic coordinates of Arnott and Hukins since the use of other geometries with the same force constants failed to give results in agreement with experimental evidence. The calculations also demonstrate the lowering effect of hydration on the anionic PO stretching frequencies. Experimental results show that the 814-cm^?1 band observed in the spectra of 5′GMP gel arises from a different vibrational mode than that of the 814-cm^?1 band of A-DNA.