UV Photoelectron Spectroscopy and Ab Initio Characterization of Valence Orbital Structures and Conformations of Neutral Phosphate Esters

Abstract

The Hel UV photoelectron spectrum of trimethyl phosphate (TMP) has been measured and interpreted with the aid of SCF molecular orbital calculations carried out with STO-3G, STO-3G* and 4–31G basis functions. The photoelectron spectrum of TMP is more accurately reproduced by results from 4–31G calculations than by results from STO-3G or STO-3G* calculations. However, all three basis sets yield results which predict the same assignment of the photoelectron spectrum. Results at the 4–31G level indicate that whether calculations are based on crystallographic bond angles and bond lengths or on STO-3G optimized geometries has little effect on the energetic ordering of the upper occupied orbitals. The energetic ordering of orbitals is also found to be only weakly dependent upon the torsional angle φ, describing rotation of ester groups about P-O bonds and upon the torsional angle ψ, describing rotation of methyl groups about C-O bonds. For trimethyl phosphate, with C₃ symmetry, the vertical ionization potentials of the upper occupied orbitals are 10.81 eV (8e), 11.4 eV (9a), 11.93 eV (7e), 12.6–12.9 eV (8a and 6e), 14.4 eV (7a) and 15.0–16.0 eV(5e and 6a). Calculations at the 4–31G level indicate that many of the highest occupied orbitals in neutral dimethyl phosphate and methyl phosphate have energies and electron distributions similar to orbitals in TMP.

For TMP, a search for optimized values of φ and ψ has been carried out at the STO-3G* level. In agreement with previous NMR studies and with classical potential calculations, the STO- 3G* results indicate that both the gauche φ= 53.1 °) and anticlinal (φ = 141.9°) conformations are thermally accessible. Also in agreement with the classical potential calculations, the STO-3G* results predict that in the all gauche conformation energy is minimized when the methyl groups assume a staggered geometry (ψ= 60° to 80°) and that an energy maximum occurs for an eclipsed geometry (ψ = 0° to 20°). A study of the dependence of optimized values of O-P-O ester bond angles on the torsional angles, φ, was carried out at the STO-3G, STO-3G* and 4–31G levels. The results demonstrate that for C₃ symmetry, the coupling of O-P-O angles to φ is influenced by repulsive steric interactions.     

Ab initio molecular orbital calculations on furanose sugars: a study with the 6-31G basis set

Ab initio molecular orbital calculations were performed on 2-deoxy-beta-D-glycero-tetrofuranose (1) using the 6-31G* basis set to evaluate the effect of ring conformation on the molecular parameters (bond lengths, angles, and torsions). Geometric optimizations were conducted on the planar and ten envelope conformers of 1, and these data were compared to those obtained from previous calculations using the STO-3G and 3-21G basis sets. Conformational energy profiles derived from 3-21G and 6-31G* data were found to be qualitatively comparable. The effect of furanose ring conformation on key bond lengths (e.g., C-H, C-O), bond angles (e.g., COC), and bond torsions (e.g., the exoanomeric C-1-O-1 torsion) was examined, and a qualitative agreement was observed between the 3-21G and 6-31G* analyses. The results indicate that, for semi-quantitative ab initio studies of intact carbohydrates, the 3-21G basis set is sufficient, and that the STO-3G basis set should not be employed unless crude structural approximations are desired. The observed concerted behavior of C-O bond lengths in the vicinity of the anomeric carbon of the aldofuranose ring has suggested a possible role of C-1-O-1 bond orientation in affecting the mechanism of glycoside bond hydrolysis.     

Role of azaamino acid residue in beta-turn formation and stability in designed peptide.

The structural perturbation induced by C(alpha)-->N(alpha) exchange in azaamino acid-containing peptides was predicted by ab initio calculation of the 6-31G* and 3-21G* levels. The global energy-minimum conformations for model compounds, For-azaXaa-NH2 (Xaa=Gly, Ala, Leu) appeared to be the beta-turn motif with a dihedral angle of phi= +/- 90 degrees, psi=0 degrees. This suggests that incorporation of the azaXaa residue into the i+2 position of designed peptides could stabilize the beta-turn structure. The model azaLeu-containing peptide, Boc-Phe-azaLeu-Ala-OMe, which is predicted to adopt a beta-turn conformation was designed and synthesized in order to experimentally elucidate the role of the azaamino acid residue. Its structural preference in organic solvents was investigated using 1H NMR, molecular modelling and IR spectroscopy. The temperature coefficients of amide protons, the characteristic NOE patterns, the restrained molecular dynamics simulation and IR spectroscopy defined the dihedral angles [ (phi i+1, psi i+1) (phi i+2, psi i+2)] of the Phe-azaLeu fragment in the model peptide, Boc-Phe-azaLeu-Ala-OMe, as [(-59 degrees, 127 degrees) (107 degrees, -4 degrees)]. This solution conformation supports a betaII-turn structural preference in azaLeu-containing peptides as predicted by the quantum chemical calculation. Therefore, intercalation of the azaamino acid residue into the i+2 position in synthetic peptides is expected to provide a stable beta-turn formation, and this could be utilized in the design of new peptidomimetics adopting a beta-turn scaffold.     

Conformational studies of cyclic tetrapeptides. Evidence for a bis gamma-turn conformation for chlamydocin and Ala4-chlamydocin in nonpolar solvents

The conformations of chlamydocin and cyclo (Ala-Aib-Phe-D-Pro) (Ala4-chlamydocin) in chloroform have been investigated by nuclear magnetic resonance, infrared and circular dichroism spectroscopy. The data obtained from these experiments establish an all transoid, bis gamma-turn conformation for both compounds in chloroform with the following torsional angles (+/- 20 degrees): Ala4-chlamydocin: Aib, phi + 60 degrees, psi - 50 degrees; omega + 160 degrees; Phe phi - 120 degrees, psi + 120 degrees, omega - 160 degrees; D-Pro phi + 60 degrees, psi - 55 degrees, omega + 160 degrees; Ala phi - 110 degrees, psi + 110 degrees, omega - 160 degrees. Chlamydocin adopts a closely related conformation in neat chloroform. Nuclear Overhauser Effect (NOE) data are utilized to assign amide bond geometries in the cyclic tetrapeptide ring system.     

Empirical torsional potential functions from protein structure data. Phi- and psi-potentials for non-glycyl amino acid residues.

The torsional potential functions Vt(phi) and Vt(psi) around single bonds N--C alpha and C alpha--C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (phi, psi)-plane with the value of Vtot(phi, psi), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in L-configuration, are Vt(phi) = 1.0 cos (phi + 60 degrees); Vt(psi) = 0.5 cos (psi + 60 degrees) - 1.0 cos (2 psi + 30 degrees) - 0.5 cos (3 psi + 30 degrees). The dipeptide energy maps Vtot(phi, psi) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line psi = 0 degrees. These functions derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.     

Conformational preferences of cyclopropyl peptides. Crystal structure of (E)-DL-1-benzamido-1-methoxycarbonyl-2-chlorcyclopropane (BCP)

Crystal structure analysis of (E)-DL-1-benzamido-1-methoxycarbonyl-2-chlorocyclopropane (C12H12NO3Cl) is reported. The phi' (about N1-C1 bond) and psi' (about C1-C11 bond) torsional angles for this compound are -62.5 degrees and -33.0 degrees, respectively, and are close to the phi, psi values of the 3(10) helix and the alpha-helix. Semi-empirical potential energy calculations are performed on a cyclopropyl dipeptide which is a special case of alpha,alpha-disubstituted dipeptide where the alpha-carbon and the two substituent carbon atoms form a 3-membered ring. Our calculations show that different types of helics: alpha-, gamma-, pi-, omega-, 3(10-) and delta-helices, are energetically favorable. Another interesting possibility is the formation of a cyclic pentapeptide with five-fold symmetry. The effect of substitutions on C beta atom are also studied with the help of potential energy maps. Selective substitutions on C beta atom may be used effectively to restrict either phi or psi values into a very narrow range.     

The normal modes of the gramicidin-A dimer channel.

The dynamics of the gramicidin-A dimer channel is studied in the harmonic approximation by a vibrational analysis of the atomic motions relative to their equilibrium positions. The system is represented by an empirical potential energy function, and all degrees of freedom (bonds lengths, bond angles, and torsional angles) are allowed to vary. The thermal fluctuations in the backbone dihedral angles phi and psi, atomic root mean square displacements, and the correlations between the different amide planes are computed. It is found that only adjacent dihedral psi i and phi i+1 are strongly correlated, while different hydrogen-bonded amide planes are only weakly correlated. Modes with relatively low vibrational frequencies (75-175 cm-1) make the dominant contributions to the carbonyl librations. The general flexibility of the structure and the role of carbonyl librations in the ion transport mechanism are discussed.     

Effects of structure on alpha C-H bond enthalpies of amino acid residues: relevance to H transfers in enzyme mechanisms and in protein oxidation.

The bond dissociation enthalpies (BDE) of all of the amino acid residues, modeled by HC(O)NHCH(R)C(O)NH(2) (PH(res)), were determined at the B3LYP/6-31G//B3LYP/6-31G level, coupled with isodesmic reactions. The results for neutral side chains with phi, psi angles approximately 180 degrees, approximately 180 degrees in ascending order, to an expected accuracy of +/-10 kJ mol(-)(1), are Asn 326; cystine 330; Asp 332; Gln 334; Trp 337; Arg 340; Lys 340; Met 343; His 344; Phe 344; Tyr 344; Leu 344; Ala 345; Cys 346; Ser 349; Gly 350; Ile 351; Val 352; Glu 354; Thr 357; Pro-cis 358; Pro-trans 369. BDEs calculated at the ROMP2/6-31G//B3LYP/6-31G level exhibit the same trends but are approximately 7 kJ mol(-)(1) higher. All BDEs are smaller than those of typical secondary or tertiary C-H bonds due to the phenomenon of captodative stabilization. The stabilization is reduced by changes in the phi,psi angles. As a result the BDEs increase by about 10 kJ mol(-)(1) in beta-sheet and 40 kJ mol(-)(1) in alpha-helical environments, respectively. In effect the alpha C-H BDEs can be "tuned" from about 345 to 400 kJ mol(-)(1) by adjusting the local environment. Some very significant effects of this are seen in the current literature on H-transfer processes in enzyme mechanisms and in oxidative damage to proteins. These observations are discussed in terms of the findings of the present study.     

Facile transition between 3(10)- and alpha-helix: structures of 8-, 9-, and 10-residue peptides containing the -(Leu-Aib-Ala)2-Phe-Aib- fragment.

A structural transition from a 3(10)-helix to an alpha-helix has been characterized at high resolution for an octapeptide segment located in 3 different sequences. Three synthetic peptides, decapeptide (A) Boc-Aib-Trp-(Leu-Aib-Ala)2-Phe-Aib-OMe, nonapeptide (B) Boc-Trp-(Leu-Aib-Ala)2-Phe-Aib-OMe, and octapeptide (C) Boc-(Leu-Aib-Ala)2-Phe-Aib-OMe, are completely helical in their respective crystals. At 0.9 A resolution, R factors for A, B, and C are 8.3%, 5.4%, and 7.3%, respectively. The octapeptide and nonapeptide form ideal 3(10)-helices with average torsional angles phi(N-C alpha) and psi(C alpha-C') of -57 degrees, -26 degrees C and -60 degrees, -27 degrees for B. The 10-residue peptide (A) begins as a 3(10)-helix and abruptly changes to an alpha-helix at carbonyl O(3), which is the acceptor for both a 4-->1 hydrogen bond with N(6)H and a 5-->1 hydrogen with N(7)H, even though the last 8 residues have the same sequence in all 3 peptides. The average phi, psi angles in the decapeptide are -58 degrees, -28 degrees for residues 1-3 and -63 degrees, -41 degrees for residues 4-10. The packing of helices in the crystals does not provide any obvious reason for the transition in helix type. Fourier transform infrared studies in the solid state also provide evidence for a 3(10)- to alpha-helix transition with the amide I band appearing at 1,656-1,657 cm-1 in the 9- and 10-residue peptides, whereas in shorter sequences the band is observed at 1,667 cm-1.     

Determination of local conformation of proteins from 1H-NMR spectroscopy data

A method is proposed to determine conformations of amino acid residues of the protein and effective correlation time tau c from cross-peak intensities in two-dimensional nuclear Overhauser enhancement (NOESY) spectra. The method consists in fitting complete relaxation matrix of dipeptide unit protons to experimental cross-peak intensities by varying phi, psi, chi torsional angles and tau c. To verify the method, NOESY spectra of basic pancreatic trypsin inhibitor (BPTI) were theoretically generated at mixing times tau m = 25-300 ms and tau c = 4 ns and used for local structure determination. The method works well with optimum for measurement of NOE intensities tau m 100-200 ms. As a result, the backbone phi, psi torsion angles were unambiguously determined at tau m = 100 ms for all but Gly residues of BPTI, and chi 1 angles were determined for the majority of side chains. The obtained dipeptide unit conformations are very close to the BPTI crystallographic structure: root mean square deviation (RMSD) of interproton distances within dipeptide units, on the average, is 0.08 A (maximal deviation 0.44 A), and RMSD of phi and psi angles are 18 and 9 degrees, respectively (maximal deviations are 44 and 22 degrees).     

Structure and conformation of linear peptides. XII. Structure of tryptophanyl-glycyl-glycine dihydrate

The crystal structure of a tripeptide, tryptophanyl-glycyl-glycine dihydrate (C15H18N4O4.2H2O, molecular weight = 354) has been determined. The crystals are orthorhombic, space group P2(1)2(1)2(1), with a = 7.875 (1) A, b = 9.009(1), c = 24.307(1) and Z = 4. The final R-index is 0.058 for 1488 reflections [sin theta)/lambda less than or equal to 0.6 A-1) with I greater than 2 sigma (I). The molecule exists as a zwitterion, with terminal NH3+ and COO- groups. The peptide units are trans and nearly perpendicular to the plane of the carboxyl group. The backbone torsion angles are: psi 1 = 132.7 degrees, omega 1 = 174.2 degrees, phi 2 = 88.2 degrees, psi 2 = 8.6 degrees, omega 2 = -179.8 degrees, phi 3 = -85.2 degrees, psi 31 = -178.1 degrees, psi 32 = 5.0 degrees. For the sidechain of tryptophan, chi 1 = -171.6 degrees, chi 2 = 101.0 degrees.     

Predominant torsional forms adopted by dipeptide conformers in solution: parameters for molecular recognition.

The present paper describes the predominant conformational forms adopted by dipeptides in aqueous solution. More than 50 dipeptides were subjected to conformational analysis using SYBYL Random Search. The resultant collections of conformers for individual dipeptides, for small groups with related side chain residues and for large groups of about 50 dipeptides were visualized graphically and analysed using a novel three-dimensional pseudo-Ramachandran plot. The distribution of conformers, weighted according to the percentage of each in the total conformer pool, was found to be restricted to nine main combinations of backbone psi (psi) and phi (phi) torsion angles. The preferred psi values were in sectors A7 (+150 degrees to +/-180 degrees), A10 (+60 degrees to +90 degrees) and A4 (-60 degrees to -90 degrees), and these were combined with preferred phi values in sectors B12 (-150 degrees to +/-180 degrees), B9 (-60 degrees to -90 degrees) and B2 (+30 degrees to +60 degrees). These combinations of psi and phi values are distinct from those found in common secondary structures of proteins. These results show that although dipeptides can each adopt many conformations in solution, each possesses a profile of common conformers that is quantifiable. A similarly weighted distribution of dipeptide conformers according to distance between amino-terminal nitrogen and carboxyl-terminal carbon shows how the preferred combinations of backbone torsional angles result in particular N-C geometries for the conformers. This approach gives insight into the important conformational parameters of dipeptides that provide the basis for their molecular recognition as substrates by widely distributed peptide transporters. It offers a basis for the rational design of peptide-based bioactive compounds able to exploit these transporters for targeting and delivery.     

A sequence preference for nucleation of alpha-helix--crystal structure of Gly-L-Ala-L-Val and Gly-L-Ala-L-Leu: some comments on the geometry of leucine zippers

The synthetic peptide Gly-L-Ala-L-Val (C10H19N3O4.3H2O; GAV) crystallizes in the monoclinic space group P21, with a = 8.052(2), b = 6.032(2), c = 15.779(7) A, beta = 98.520(1) degree, V = 757.8 A3, Dx = 1.312 g cm-3, and Z = 2. The peptide Gly-L-Ala-L-Leu (C11H21N3O4.3H2O; GAL) crystallizes in the orthorhombic space group P212121, with a = 6.024(1), b = 8.171(1), c = 32.791(1) A, V = 1614 A3, Dx = 1.289 g cm-3, and Z = 4. Their crystal structures were solved by direct methods using the program SHELXS-86, and refined to an R index of 0.05 for 1489 reflections for GAV and to an R index of 0.05 for 1563 reflections for GAL. The tripeptides exist as a zwitterion in the crystal and assume a near alpha-helical backbone conformation with the following torsion angles: psi 1 = -150.7 degrees; phi 2, psi 2 = -68.7 degrees, -38.1 degrees; phi 3, psi 32 = -74.8 degrees, -44.9 degrees, 135.9 degrees for GAV; psi 1 = -150.3 degrees; phi 2, psi 2 = -67.7 degrees, -38.9 degrees; phi 3, psi 31, psi 32 = -72.2 degrees, -45.3 degrees, 137.5 degrees for GAL. Both the peptide units in both of the tripeptides show significant deviation from planarity [omega 1 = -171.3(6) degrees and omega 2 = -172.0(6) degrees for GAV; omega 1 = -171.9(5) degrees and omega 2 = -173.2(6) degrees for GAL]. The side-chain conformational angles chi 21 and chi 22 are -61.7(5) degrees and 175.7(5) degrees, respectively, for valine, and the side-chain conformations chi 12 and chi 23's are -68.5(5) degrees and (-78.4(6) degrees, 159.10(5) degrees) respectively, for leucine. Each of the tripeptide molecule is held in a near helical conformation by a water molecule that bridges the NH3+ and COO- groups, and acts as the fourth residue needed to complete the turn by forming two hydrogen bonds. Two other water molecules form intermolecular hydrogen bonds in stabilizing the helical structure so that the end result is a column of molecules that looks like an alpha-helix.     

Dihedral angles of tripeptides in solution directly determined by polarized Raman and FTIR spectroscopy

The amide I mode of the peptide linkage is highly delocalized in peptides and protein segments due to through-bond and through-space vibrationally coupling between adjacent peptide groups. J. Phys. Chem. B. 104:11316-11320) used coherent femtosecond infrared (IR) spectroscopy to determine the excitonic coupling energy and the orientational angle between the transition dipole moments of the interacting amide I modes of cationic tri-alanine in D(2)O. Recently, the same parameters were determined for all protonation states of tri-alanine by analyzing the amide I bands in the respective IR and isotropic Raman spectra (. J. Am. Chem. Soc. 119:1720-1726.). In both studies, the dihedral angles phi and psi were then obtained by utilizing the orientational dependence of the coupling energy obtained from ab initio calculations on tri-glycine in vacuo (. J. Raman Spectrosc. 29:81-86) to obtain an extended 3(1) helix-like structure for the tripeptide. In the present paper, a novel algorithm for the analysis of excitonic coupling between amide I modes is presented, which is based on the approach by Schweitzer-Stenner et al. but avoids the problematic use of results from ab initio calculations. Instead, the dihedral angles are directly determined from infrared and visible polarized Raman spectra. First, the interaction energy and the corresponding degree of wave-function mixing were obtained from the amide I profile in the isotropic Raman spectrum. Second, the depolarization ratios and the amide I profiles in the anisotropic Raman and IR-absorption spectra were used to determine the orientational angle between the peptide planes and the transition dipole moments, respectively. Finally, these two geometric parameters were utilized to determine the dihedral angles phi and psi between the interacting peptide groups. Stable extended conformations with dihedral angles in the beta-sheet region were obtained for all protonation states of tri-alanine, namely phi(+) = -126 degrees, psi(+) = 178 degrees; phi(+/-) = -110 degrees, psi(+/-) = 155 degrees; and phi(-) = -127 degrees, psi(-) = 165 degrees for the cationic, zwitterionic, and anionic state, respectively. These values reflect an extended beta-helix structure. Tri-glycine was found to be much more heterogeneous in that different extended conformers coexist in the cationic and zwitterionic state, which yield a noncoincidence between isotropic and anisotropic Raman scattering. Our study introduces vibrational spectroscopy as a suitable tool for the structure analysis of peptides in solution and tripeptides as suitable model systems for investigating the role of local interactions in determining the propensity of peptide segments for distinct secondary structure motifs.     

The beta-turn scaffold of tripeptide containing an azaphenylalanine residue

The conformational preferences of azaphenylalanine-containing peptide were investigated using a model compound, Ac-azaPhe-NHMe with ab initio method at the HF/3-21G and HF/6-31G(*) levels, and the seven minimum energy conformations with trans orientation of acetyl group and the 4 minimum energy conformations with cis orientation of acetyl group were found at the HF/6-31G(*) level if their mirror images were not considered. An average backbone dihedral angle of the 11 minimum energy conformations is phi=+/-91 degrees +/-24 degrees , psi =+/-18 degrees +/-10 degrees (or +/-169 degrees +/-8 degrees ), corresponding to the i+2 position of beta-turn (delta(R)) or polyproline II (beta(P)) structure, respectively. The chi(1) angle in the aromatic side chain of azaPhe residue adopts preferentially between +/-60 degrees and +/-130 degrees, which reflect a steric hindrance between the N-terminal carbonyl group or the C-terminal amide group and the aromatic side chain with respect to the configuration of the acetyl group. These conformational preferences of Ac-azaPhe-NHMe predicted theoretically were compared with those of For-Phe-NHMe to characterize the structural role of azaPhe residue. Four tripeptides containing azaPhe residue, Boc-Xaa-azaPhe-Ala-OMe [Xaa=Gly(1), Ala(2), Phe(3), Asn(4)] were designed and synthesized to verify whether the backbone torsion angles of azaPhe reside are still the same as compared with theoretical conformations and how the preceding amino acids of azaPhe residue perturb the beta-turn skeleton in solution. The solution conformations of these tripeptide models containing azaPhe residue were determined in CDCl(3) and DMSO solvents using NMR and molecular modeling techniques. The characteristic NOE patterns, the temperature coefficients of amide protons and small solvent accessibility for the azapeptides 1-4 reveal to adopt the beta-turn structure. The structures of azapeptides containing azaPhe residue from a restrained molecular dynamics simulation indicated that average dihedral angles [(phi(1), psi(1)), (phi(2), psi(2))] of Xaa-azaPhe fragment in azapeptide, Boc-Xaa-azaPhe-Ala-OMe were [(-68 degrees, 135 degrees ), (116 degrees, -1 degrees )], and this implies that the intercalation of an azaPhe residue in tripeptide induces the betaII-turn conformation, and the volume change of a preceding amino acid of azaPhe residue in tripeptides would not perturb seriously the backbone dihedral angle of beta-turn conformation. We believe such information could be critical in designing useful molecules containing azaPhe residue for drug discovery and peptide engineering.     

Ab initio computational study of beta-cellobiose conformers using B3LYP/6-311++G**

The molecular structure of 27 conformers of beta-cellobiose were studied in vacuo through gradient geometry optimization using B3LYP density functionals and the 6-311++G** basis set. The conformationally dependent geometry changes and energies were explored as well as the hydrogen-bonding network. The lowest electronic energy structures found were not those suggested from available crystallographic and NMR solution data, where the glycosidic dihedral angles fall in the region (phi, psi) approximately (40 degrees, -20 degrees ). Rather, 'flipped' conformations in which the dihedral angles are in the range (phi, psi) approximately (180 degrees, 0 degrees ) are energetically more stable by approximately 2.5 kcal/mol over the 'experimentally accepted' structure. Further, when the vibrational free energy, deltaG, obtained from the calculated frequencies, is compared throughout the series, structures with (phi, psi) in the experimentally observed range still have higher free energy ( approximately 2.0 kcal/mol) than 'flipped' forms. The range of bridging dihedral angles of the 'normal' conformers, resulting from the variance in the phi dihedral is larger than that found in the 'flipped' forms. Due to this large flat energy surface for the normal conformations, we surmise that the summation of populations of these conformations will favor the 'normal' conformations, although evidence suggests that polar solvent effects may play the dominant role in providing stability for the 'normal' forms. Even though some empirical studies previously found the 'flipped' conformations to be lowest in energy, these studies have been generally discredited because they were in disagreement with experimental results. Most of the DFT/ab initio conformations reported here have not been reported previously in the ab initio literature, in part because the use of less rigorous theoretical methods, i.e. smaller basis sets, have given results in general agreement with experimental data, that is, they energetically favored the 'normal' forms. These are the first DFT/ab initio calculations at this level of theory, apparently because of the length and difficulty of carrying out optimizations at these high levels.     

Ab initio Investigation of the Molecular Structure of Methyl Methoxymethyl Phosphonate,a Promising Nuclease-resistant Alternative of the Phosphodiester Linkage

Abstract

Conformational flexibility of the methyl methoxymethyl phosphonate anion (CH₃-O-PO₂- CH₂-O-CH₃)^?, a nuclease resistant alternative to the phosphodiester linkage in DNA, have been investigated by ab initio quantum mechanical calculations. The potential of backbone torsional degrees of freedom of methyl methoxymethyl phosphonate anion (MMP) was determined at the Hartree-Fock (HF) 3–21G* level using the adiabatic mapping technique. Energies, geometries, and effective atomic charges of different conformers were calculated at HF/6–31G* and MP2/6–31G* levels of theory. These were compared to the results obtained for dimethyl phosphate calculated at the same level. The impact on DNA structure from inserting a methylene group between phosphorus and oxygen of the nucleoside sugar moiety was examined via distance- and angle-constrained geometry optimizations. Due to its high flexibility, MMP has been shown to be compatible with both A and B forms of DNA.     

Flexibility of bovine pancreatic trypsin inhibitor

The native conformation of a protein may be expressed in terms of the dihedral angles, phi's and psi's for the backbone, and kappa's for the side chains, for a given geometry (bond lengths and bond angles). We have developed a method to obtain the dihedral angles for a low-energy structure of a protein, starting with the X-ray structure; it is applied here to examine the degree of flexibility of bovine pancreatic trypsin inhibitor. Minimization of the total energy of the inhibitor (including nonbonded, electrostatic, torsional, hydrogen bonding, and disulfide loop energies) yields a conformation having a total energy of -221 kcal/mol and a root mean square deviation between all atoms of the computed and experimental structures of 0.63 A. The optimal conformation is not unique, however, there being at least two other conformations of low-energy (-222 and -220 kcal/mol), which resemble the experimental one (root mean square deviations of 0.66 and 0.64 A, respectively). These three conformations are located in different positions in phi, psi space, i.e., with a total deviation of 81 degrees, 100 degrees and 55 degrees from each other (with a root mean square deviation of several degrees per dihedral angle from each other). The nonbonded energies of the backbones, calculated along lines in phi, psi space connecting these three conformations, are all negative, without any intervening energy barriers (on an energy contour map in the phi, psi plane). Side chains were attached at several representative positions in this plane, and the total energy was minimized by varying the kappa's. The energies were of approximately the same magnitude as the previous ones, indicating that the conformation of low energy is flexible to some extent in a restricted region of phi, psi space. Interestingly, the difference delta phi i+1 in phi i+1 for the (i + 1)th residue from one conformation to another is approximately the same as -delta psi i for the ith residue; i.e., the plane of the peptide group between the ith and (i + 1)th residues re-orient without significant changes in the positions of the other atoms. The flexibility of the orientations of the planes of the peptide groups is probably coupled in a cooperative manner to the flexibility of the positions of the backbone and side-chain atoms.     

Serine and threonine residues bend alpha-helices in the chi(1) = g(-) conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.     

An observation of non-superimposable stereogeometrical features in a non-chiral one-component beta-Ala model peptide

This paper describes the chemical synthesis and crystal molecular conformation of a non-chiral beta-Ala containing model peptide Boc-beta-Ala-Acc5-OCH3. The analysis revealed the existence of two crystallographically independent molecules A and B, in the asymmetric unit. Unexpectedly, while the magnitudes of the backbone torsion angles in both molecules are remarkably similar, the signs of the corresponding torsion angles are reverse therefore, inclining us to suggest the existence of non-superimposable stereogeometrical features in a non-chiral one-component beta-Ala model system. The critical mu torsion angle around CbetaH2-CalphaH2 bond of the beta-Ala residue represents a typical gauche orientation i.e., mu = 67.7 degrees in A and mu = -61.2 degrees in B, providing the molecule an overall crescent shaped topology. The observed conformation contrasts markedly to those determined for the correlated non-chiral model peptides: Boc-beta-Ala-Acc6-OCH3 and Boc-beta-Ala-Aib-OCH3 signifying the role of stereocontrolling elements since the stereochemically constrained Calpha, alpha-disubstituted glycyl residues (e.g., Acc5, Acc6, and the prototype Aib) are known to strongly restrict the peptide backbone conformations in the 3(10)/alpha-helical-regions ( phi approximately +/-60+/-20 degrees, psi approximately +/-30+/-20 degrees) of the Ramachandran map. Unpredictably, the preferred, phi, psi torsion angles of the Acc5 residue fall outside the helical regions of the Ramachandran map and exhibit opposite-handed twists for A and B. The implications of the semi-extended conformation of the Acc5 residue in the construction of backbone-modified novel scaffolds and peptides of biological relevance are highlighted. Taken together, the results indicate that in short linear beta-Ala containing peptides specific structural changes can be induced by selective substitution of non-coded linear- or cyclic symmetrically Calpha,alpha-disubstituted glycines, reinstating the hypothesis that in addition to conformational restrictions, the chemical nature of the neighboring side-chain substituents and local environments collectively influences the stabilization of folding-unfolding behavior of the two methylene units of a beta-Ala residue.